Q. Ma¹, C. Levy², M. Perl³

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.2, pp. 95-110, 2008, DOI:10.3970/cmes.2008.029.095

Abstract Networks of radial and longitudinally-coplanar, internal, surface cracks are typical in rifled, autofrettaged, gun barrels. In two previous papers, the separate effects of large arrays of either radial or longitudinally-coplanar semi-elliptical, internal, surface cracks in a thick-walled, cylindrical, pressure vessel under both ideal and realistic autofrettage were studied. When pressure is considered solely, radial crack density and longitudinal crack spacing were found to have opposing effects on the prevailing stress intensity factor, K_IP. Furthermore, the addition of the negative stress intensity factor (SIF), K_IA, resulting from the residual stress field due to autofrettage, whether ideal or… More >

M. Haji Mohammadi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.2, pp. 75-94, 2008, DOI:10.3970/cmes.2008.029.075

Abstract In this paper, the truly Meshless Local Petrov-Galerkin (MLPG) method is extended for computation of steady incompressible flows, governed by the Navier--Stokes equations (NSE), in vorticity-stream function formulation. The present method is a truly meshless method based on only a number of randomly located nodes. The formulation is based on two equations including stream function Poisson equation and vorticity advection-dispersion-reaction equation (ADRE). The meshless method is based on a local weighted residual method with the Heaviside step function and quartic spline as the test functions respectively over a local subdomain. Radial basis functions (RBF) interpolation… More >

E. L. Albuquerque¹, M. H. Aliabadi²

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.2, pp. 63-74, 2008, DOI:10.3970/cmes.2008.029.063

Abstract This paper presents a boundary element formulation for the analysis of thin shallow shells. Classical plate bending and plane elasticity formulations are coupled and effects of curvature are treated as body forces. The body forces are written as a sum of approximation functions multiplied by coefficients. Domain integrals that arise in the formulation are transformed into boundary integrals by the radial integration method. Two different approximation functions are employed, that is 1 + r and r² log r. The method is applied to several problems and the accuracy of each approximation function is assessed by comparison with More >

Yongjie Zhang¹, Boyle C. Cheng², Changho Oh¹, Jessica L. Spehar², James Burgess³

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.2, pp. 55-62, 2008, DOI:10.3970/cmes.2008.029.055

Abstract The spinal column plays a vital biomechanical role in the human body by providing structural support and facilitating motion. As degenerative changes occur in the spine, surgical treatment may be necessary in certain instances. Such treatments seek to address pain, frequently through the restriction of spinal motion. Traditional spinal implant devices are designed to restrict the motion of a functional spinal unit (FSU) but newer device designs allow for semi-constrained motion such as spinal arthroplasty devices. In this study, a sequence of fluoroscopic imaging data was recorded during the flexibility protocol with an interspinous process… More >

Leevan Ling¹, Tomoya Takeuchi²

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 45-54, 2008, DOI:10.3970/cmes.2008.029.045

Abstract The method of fundamental solutions is coupled with the boundary control technique to solve the Cauchy problems of the Laplace Equations. The main idea of the proposed method is to solve a sequence of direct problems instead of solving the inverse problem directly. In particular, we use a boundary control technique to obtain an approximation of the missing Dirichlet boundary data; the Tikhonov regularization technique and the L-curve method are employed to achieve such goal stably. Once the boundary data on the whole boundary are known, the numerical solution to the Cauchy problem can be More >

H.M. Lee¹, H.S. Park^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 29-44, 2008, DOI:10.3970/cmes.2008.029.029

Abstract This paper presents a computational model to estimate deformed shapes of beam structures using 3D coordinate information from terrestrial laser scanning (TLS). The model is composed of five components: 1) formulation of polynomial shape function, 2) application of boundary condition, 3) inducement of compatibility condition, 4) application of the least square method and 5) evaluation of error vector and determination of reasonable polynomial shape function. In the proposed model, the optimal degree of polynomial function is selected based on the complexity of beam structures, instead of using a specific degree of polynomial function. The chosen More >

V.C. Mariani¹, E.E.M. Alonso², S. Peters³

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 15-28, 2008, DOI:10.3970/cmes.2008.029.015

Abstract An application of Newton's method for linearization of advective terms given by the discretization on unstructured Voronoi meshes for the incompressible Navier-Stokes equations is proposed and evaluated in this article. One of the major advantages of the unstructured approach is its application to very complex geometrical domains and the mesh is adaptable to features of the flow. Moreover, in this work comparisons with the literature results in bi-dimensional lid-driven cavities for different Reynolds numbers allow us to assess the numerical properties of the new proposed finite-volume scheme. Results for the components of the velocity, and More >

H. T. Wang¹, Z. H. Yao²

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 1-14, 2008, DOI:10.3970/cmes.2008.029.001

Abstract Carbon nanotubes (CNTs) may provide ultimate enhancement in stiffness and strength for composite materials. This paper presents a rigid-fiber-based boundary integral equation formulation for the numerical simulation of debonding process and the corresponding strength of CNT reinforced composites. The CNT/matrix interfaces are assumed to fail when the interfacial shear force reaches a prescribed threshold, and the CNTs and matrix are considered to be detached in the failed areas. The matrix with one or several tens of originally well-bonded CNTs is subjected to an incremental tensile load and the effective stress-strain relations are readily obtained by… More >

S.Y. Long^1,2,3, K.Y. Liu^1,2,4, G.Y. Li¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.3, pp. 203-216, 2008, DOI:10.3970/cmes.2008.028.203

Abstract A meshless local Petrov-Galerkin method (MLPG) for the analysis of the elasto-plastic fracture problem is presented in this paper. The meshless method uses the moving least squares (MLS) to approximate the field functions. The shape function has not Kronecker Delta properties for the trial-function interpolation, and a direct interpolation method is adopted to impose essential boundary conditions. The MLPG method does not involve any domain and singular integrals if body force is ignored. It only involves a regular boundary integral. Two numerical examples show that results obtained by the present method have a good agreement More >

D. Kourounis¹, L.N. Gergidis¹, A. Charalambopoulos¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.3, pp. 185-202, 2008, DOI:10.3970/cmes.2008.028.185

Abstract The sensitivity of analytical solutions of the direct acoustic scattering problem in prolate spheroidal geometry on the wavenumber and shape, is extensively investigated in this work. Using the well known Vekua transformation and the complete set of radiating "outwards'' eigensolutions of the Helmholtz equation, introduced in our previous work ([Charalambopoulos and Dassios(2002)], [Gergidis, Kourounis, Mavratzas, and Charalambopoulos (2007)]), the scattered field is expanded in terms of it, detouring so the standard spheroidal wave functions along with their inherent numerical deficiencies. An approach is employed for the determination of the expansion coefficients, which is optimal in… More >