S. N. Atluri¹, H. T. Liu², Z. D. Han²

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.1, pp. 1-16, 2006, DOI:10.3970/cmes.2006.015.001

Abstract The Finite Difference Method (FDM), within the framework of the Meshless Local Petrov-Galerkin (MLPG) approach, is proposed in this paper for solving solid mechanics problems. A "mixed'' interpolation scheme is adopted in the present implementation: the displacements, displacement gradients, and stresses are interpolated independently using identical MLS shape functions. The system of algebraic equations for the problem is obtained by enforcing the momentum balance laws at the nodal points. The divergence of the stress tensor is established through the generalized finite difference method, using the scattered nodal values and a truncated Taylor expansion. The traction boundary conditions are imposed in… More >

R. Vignjevic¹, T. De Vuyst², J. C. Campbell¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.1, pp. 35-48, 2006, DOI:10.3970/cmes.2006.013.035

Abstract An approach to the treatment of contact problems involving frictionless sliding and separation under large deformations in meshless methods is proposed. The method is specially suited for non-structured spatial discretisation. The contact conditions are imposed using a contact potential for particles in contact. Inter-penetration is checked as a part of the neighbourhood search. In the case of conventional SPH contact conditions are enforced on the boundary layer 2h thick while in the case of the normalized SPH contact conditions are enforced for the particles lying on the contact surface. The implementation of the penalty based contact algorithm for the explicit… More >

R. C. Batra¹, H.-K. Ching¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.6, pp. 717-730, 2002, DOI:10.3970/cmes.2002.003.717

Abstract The Meshless Local Petrov-Galerkin (MLPG) method is used to analyze transient deformations near either a crack or a notch tip in a linear elastic plate. The local weak formulation of equations governing elastodynamic deformations is derived. It results in a system of coupled ordinary differential equations which are integrated with respect to time by a Newmark family of methods. Essential boundary conditions are imposed by the penalty method. The accuracy of the MLPG solution is established by comparing computed results for one-dimensional wave propagation in a rod with the analytical solution of the problem. Results are then computed for the… More >

Soon Wan Chung¹, Yoo Jin Choi¹, Seung Jo Kim¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 707-718, 2003, DOI:10.3970/cmes.2003.004.707

Abstract In this paper, an automatic finite element remeshing algorithm based on the bubble packing method is utilized for the purpose of numerical simulations of localized damage, because fine meshes are needed to represent the gradually concentrated damage. The bubble packing method introduces two parameters that easily control the remeshing criterion and the new mesh size. The refined area is determined by \textit {a posteriori} error estimation utilizing the value obtained from Superconvergent Patch Recovery. The isotropic ductile damage theory, founded on continuum damage mechanics, is used for this damage analysis. It was successfully shown in the numerical examples (upsetting and… More >

N. Sukumar¹, J. E. Bolander¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 691-706, 2003, DOI:10.3970/cmes.2003.004.691

Abstract In this paper, we explore the numerical approximation of discrete differential operators on non-uniform grids. The Voronoi cell and the notion of natural neighbors are used to approximate the Laplacian and the gradient operator on irregular grids. The underlying weight measure used in the numerical computations is the {\em Laplace weight function}, which has been previously adopted in meshless Galerkin methods. We develop a difference approximation for the diffusion operator on irregular grids, and present numerical solutions for the Poisson equation. On regular grids, the discrete Laplacian is shown to reduce to the classical finite difference scheme. Two techniques to… More >

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

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 665-678, 2003, DOI:10.3970/cmes.2003.004.665

Abstract The numerical implementation of the truly Meshless Local Petrov-Galerkin (MLPG) type weak-forms of the displacement and traction boundary integral equations is presented, for solids undergoing small deformations. In the accompanying part I of this paper, the general MLPG/BIE weak-forms were presented [Atluri, Han and Shen (2003)]. The 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,… More >

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 each one is surrounded by… 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 used. The first utilizes the… More >

Yan Zhou¹, Qingwei Ma^*

CMES-Computer Modeling in Engineering & Sciences, Vol.115, No.2, pp. 131-147, 2018, DOI: 10.3970/cmes.2018.00249

Abstract An identification technique for sharp interface and penetrated isolated particles is developed for simulating two-dimensional, incompressible and immiscible two-phase flows using meshless particle methods in this paper. This technique is based on the numerically computed density gradient of fluid particles and is suitable for capturing large interface deformation and even topological changes such as merging and breaking up of phases. A number of assumed particle configurations will be examined using the technique, including these with different level of randomness of particle distribution. The tests will show that the new technique can correctly identify almost all the interface and isolated particles,… More >

M. Dehghan¹, M. Abbaszadeh², A. Mohebbi³

CMES-Computer Modeling in Engineering & Sciences, Vol.107, No.6, pp. 481-516, 2015, DOI:10.3970/cmes.2015.107.481

Abstract In the current paper the two-dimensional time fractional Klein-Kramers equation which describes the subdiffusion in the presence of an external force field in phase space has been considered. The numerical solution of fractional Klein-Kramers equation is investigated. The proposed method is based on using finite difference scheme in time variable for obtaining a semi-discrete scheme. Also, to achieve a full discretization scheme, the Kansa's approach and meshless local Petrov-Galerkin technique are used to approximate the spatial derivatives. The meshless method has already proved successful in solving classic and fractional differential equations as well as for several other engineering and physical… More >