Jan Sladek¹, Vladimir Sladek¹, Chuanzeng Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 637-648, 2003, DOI:10.3970/cmes.2003.004.637

Abstract A new computational method for solving transient elastodynamic initial-boundary value problems in continuously non-homogeneous solids, based on the meshless local Petrov-Galerkin (MLPG) method, is proposed in the present paper. The moving least squares (MLS) is used for interpolation and the modified fundamental solution as the test function. The local Petrov-Galerkin method for unsymmetric weak form in such a way is transformed to the local boundary integral equations (LBIE). The analyzed domain is divided into small subdomains, in which a weak solution is assumed to exist. Nodal points are randomly spread in the analyzed domain and More >

E. J. Sellountos¹, D. Polyzos²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 619-636, 2003, DOI:10.3970/cmes.2003.004.619

Abstract A new meshless local Petrov-Galerkin (MLPG) method for solving two dimensional frequency domain elastodynamic problems is proposed. Since the method utilizes, in its weak formulation, either the elastostatic or the frequency domain elastodynamic fundamental solution as test function, it is equivalent to the local boundary integral equation (LBIE) method. Nodal points spread over the analyzed domain are considered and the moving least squares (MLS) interpolation scheme for the approximation of the interior and boundary variables is employed. Two integral equations suitable for the integral representation of the displacement fields in the local sub- domains are… More >

Q. Li¹, S. Shen¹, Z. D. Han¹, S. N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 571-586, 2003, DOI:10.3970/cmes.2003.004.571

Abstract In this paper, a truly meshless method, the Meshless Local Petrov-Galerkin (MLPG) Method, is developed for three-dimensional elasto-statics. The two simplest members of MLPG family of methods, the MLPG type 5 and MLPG type 2, are combined, in order to reduce the computational requirements and to obtain high efficiency. The MLPG5 method is applied at the nodes inside the 3-D domain, so that any domain integration is eliminated altogether, if no body forces are involved. The MLPG 2 method is applied at the nodes that are on the boundaries, and on the interfaces of material More >

L. F. Qian^1,2, R. C. Batra³, L. M. Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 519-534, 2003, DOI:10.3970/cmes.2003.004.519

Abstract We use a meshless local Petrov-Galerkin (MLPG) method to analyze three-dimensional infinitesimal elastodynamic deformations of a homogeneous rectangular plate subjected to different edge conditions. We employ a higher-order plate theory in which both transverse shear and transverse normal deformations are considered. Natural frequencies and the transient response to external loads have been computed for isotropic and orthotropic plates. Computed results are found to agree with those obtained from the analysis of the 3-dimensional problem either analytically or by the finite element method. More >

S. N. Atluri¹, Z. D. Han¹, S. Shen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 507-518, 2003, DOI:10.3970/cmes.2003.004.507

Abstract The general Meshless Local Petrov-Galerkin (MLPG) type weak-forms of the displacement & traction boundary integral equations are presented, for solids undergoing small deformations. These MLPG weak forms provide the most general basis for the numerical solution of the non-hyper-singular displacement and traction BIEs [given in Han, and Atluri (2003)], which are simply derived by using the gradients of the displacements of the fundamental solutions [Okada, Rajiyah, and Atluri (1989a,b)]. By employing the various types of test functions, in the MLPG-type weak-forms of the non-hyper-singular dBIE and tBIE over the local sub-boundary surfaces, several types of… More >

Z.Tang, S. Shen, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 177-196, 2003, DOI:10.3970/cmes.2003.004.177

Abstract A meshless numerical implementation is reported of the 2-D Fleck-Hutchinson phenomenological strain-gradient theory, which fits within the framework of the Toupin-Mindlin theories and deals with first-order strain gradients and the associated work-conjugate higher-order stresses. From a mathematical point of view, the two-dimensional Toupin-Mindlin strain gradient theory is a generalization of the Poisson-Kirchhoff plate theories, involving, in addition to the fourth-order derivatives of the displacements, also a second-order derivative. In the conventional displacement-based approaches in FEM, the interpolation of displacement requires C$^{1}$ --continuity (in order to ensure convergence of the finite element procedure for 4$^{th}$ order… More >

L. F. Qian^1,3, R. C. Batra², L. M. Chen³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 161-176, 2003, DOI:10.3970/cmes.2003.004.161

Abstract We use two meshless local Petrov-Galerkin formulations, namely, the MLPG1 and the MLPG5, to analyze infinitesimal deformations of a homogeneous and isotropic thick elastic plate with a higher-order shear and normal deformable plate theory. It is found that the two MLPG formulations give results very close to those obtained by other researchers and also by the three-dimensional analysis of the problem by the finite element method. More >

N. Mai-Duy, T. Tran-Cong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 85-102, 2003, DOI:10.3970/cmes.2003.004.085

Abstract This paper reports a mesh-free Indirect Radial Basis Function Network method (IRBFN) using Thin Plate Splines (TPSs) for numerical solution of Differential Equations (DEs) in rectangular and curvilinear coordinates. The adjustable parameters required by the method are the number of centres, their positions and possibly the order of the TPS. The first and second order TPSs which are widely applied in numerical schemes for numerical solution of DEs are employed in this study. The advantage of the TPS over the multiquadric basis function is that the former, with a given order, does not contain the… More >

E. Ferretti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 51-72, 2003, DOI:10.3970/cmes.2003.004.051

Abstract A numerical code for modeling crack propagation using the cell method is proposed. The Mohr-Coulomb criterion is used to compute the direction of crack propagation, and the new crack geometry is realized by an intra-element propagation technique. Automatic remeshing is then activated. Applications in Mode I and Mixed Mode are presented to illustrate the robustness of the implementation. More >

Dongwoo Sheen¹, Sangwon Seo², Jinwoo Cho³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 21-30, 2003, DOI:10.3970/cmes.2003.004.021

Abstract The reconstructing optimal microstructures of given homogenized coefficients of steady diffusion equation is studied. In the reconstruction, the governing equation of level set function is approximated by adding viscosity term and the numerical procedure for the evolution of the level set function for the solution is examined. The numerical experiments of reconstruction are obtained by applying a finite element method with locally fitted mesh. More >