Z. D. Han¹, A. M. Rajendran², S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 1-12, 2005, DOI:10.3970/cmes.2005.010.001

Abstract A nonlinear formulation of the Meshless Local Petrov-Galerkin (MLPG) finite-volume mixed method is developed for the large deformation analysis of static and dynamic problems. In the present MLPG large deformation formulation, the velocity gradients are interpolated independently, to avoid the time consuming differentiations of the shape functions at all integration points. The nodal values of velocity gradients are expressed in terms of the independently interpolated nodal values of displacements (or velocities), by enforcing the compatibility conditions directly at the nodal points. For validating the present large deformation MLPG formulation, two example problems are considered: 1) large deformations and rotations of… More >

J. Sorić¹, Q. Li², T. Jarak¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.4, pp. 349-358, 2004, DOI:10.3970/cmes.2004.006.349

Abstract An efficient meshless formulation based on the Local Petrov-Galerkin approach for the analysis of shear deformable thick plates is presented. Using the kinematics of a three-dimensional continuum, the local symmetric weak form of the equilibrium equations over the cylindrical shaped local sub-domain is derived. The linear test function in the plate thickness direction is assumed. Discretization in the in-plane directions is performed by means of the moving least squares approximation. The linear interpolation over the thickness is used for the in-plane displacements, while the hierarchical quadratic interpolation is adopted for the transversal displacement in order to avoid the thickness locking… More >

J. Sladek¹, V. Sladek¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 309-318, 2004, DOI:10.3970/cmes.2004.006.309

Abstract Meshless methods based on the local Petrov-Galerkin approach are proposed for solution of steady and transient heat conduction problem in a continuously nonhomogeneous anisotropic medium. Fundamental solution of the governing partial differential equations and the Heaviside step function are used as the test functions in the local weak form. It is leading to derive local boundary integral equations which are given in the Laplace transform domain. The analyzed domain is covered by small subdomains with a simple geometry. To eliminate the number of unknowns on artificial boundaries of subdomains the modified fundamental solution and/or the parametrix with a convenient cut-off… More >

B. Jin^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 253-262, 2004, DOI:10.3970/cmes.2004.006.253

Abstract In this paper, we propose a new numerical scheme for the solution of the Laplace and biharmonic equations subjected to noisy boundary data. The equations are discretized by the method of fundamental solutions. Since the resulting matrix equation is highly ill-conditioned, a regularized solution is obtained using the truncated singular value decomposition, with the regularization parameter given by the L-curve method. Numerical experiments show that the method is stable with respect to the noise in the data, highly accurate and computationally very efficient. More >

E. Ferretti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 227-244, 2004, DOI:10.3970/cmes.2004.006.227

Abstract Defining the crack path in brittle and non-brittle crack is not easy, due to several unknowns. If the direction of crack propagation can be computed by means of one of the existing criteria, it is not known whether this direction will remain constant during crack propagation. A crack initiation leads to an enhanced stress field at crack tip. During propagation, the enhanced tip stress field propagates into the solid, locally interacting with the pre-existing stress field. This interaction can lead to modifications of the propagation direction, depending on the domain and crack geometry. Moreover, trajectory deviation affects the length of… 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 discontinuities, so that the boundary… 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 MLPG/BIEs are formulated, while also… 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 theories), which inevitably involves very… 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 >