P. C. Wallstedt¹, J. E. Guilkey¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 223-232, 2007, DOI:10.3970/cmes.2007.019.223

Abstract The standard velocity projection scheme for the Material Point Method (MPM) and a typical form of the GIMP Method are examined. It is demonstrated that the fidelity of information transfer from a particle representation to the computational grid is strongly dependent on particle density and location. In addition, use of non-uniform grids and even non-uniform particle sizes are shown to introduce error. An enhancement to the projection operation is developed which makes use of already available velocity gradient information. This enhancement facilitates exact projection of linear functions and reduces the dependence of projection accuracy on particle location and density for… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 145-162, 2007, DOI:10.3970/cmes.2007.019.145

Abstract A new method is developed to solve the Dirichlet problems for the two-dimensional Laplace equation in the doubly-connected domains, namely the meshless regularized integral equations method (MRIEM), which consists of three portions: Fourier series expansion, the Fredholm integral equations, and linear equations to determine the unknown boundary conditions onartificial circles. The boundary integral equations on artificial circles are singular-free and the kernels are degenerate. When boundary-type methods are inefficient to treat the problems with complicated domains, the new method can be applicable for such problems. The new method by using the Fourier series and the Fourier coefficients can be adopted… More >

Chein-Shan Liu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 99-110, 2007, DOI:10.3970/cmes.2007.019.099

Abstract A new meshless regularized integral equation method (MRIEM) is developed to solve the interior and exterior Dirichlet problems for the two-dimensional Laplace equation, which consists of three parts: Fourier series expansion, the second kind Fredholm integral equation and an analytically regularized solution of the unknown boundary condition on an artificial circle. We find that the new method is powerful even for the problem with complex boundary shape and with random noise disturbing the boundary data. More >

J. Perko¹, B. Šarler²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 55-68, 2007, DOI:10.3970/cmes.2007.019.055

Abstract This work introduces a procedure for automated determination of weight function free parameters in moving least squares (MLS) based meshless methods for non-uniform grids. The meshless method used in present work is Diffuse Approximate Method (DAM). The DAM is structured in 2D with the one or two parameter Gaussian weigh function, 6 polynomial basis and 9 noded domain of influence. The procedure consists of three main elements. The first is definition of the reference quality function which measures the difference between the MLS approximation on non-uniform and hypothetic uniform node arrangements. The second is the construction of the object function… More >

Weiran Yuan¹, Pu Chen^1,2, Kaishin Liu^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.2, pp. 115-134, 2007, DOI:10.3970/cmes.2007.017.115

Abstract In this paper we propose a direct solution method for the quasi-unsymmetric sparse matrix (QUSM) arising in the Meshless Local Petrov-Galerkin method (MLPG). QUSM, which is conventionally treated as a general unsymmetric matrix, is unsymmetric in its numerical values, but nearly symmetric in its nonzero distribution of upper and lower triangular portions. MLPG employs trial and test functions in different functional spaces in the local domain weak form of governing equations. Consequently the stiffness matrix of the resultant linear system is a QUSM. The new solver for QUSM conducts a two-level unrolling technique for LDU factorization method and can be… More >

U. Andreaus^1,3, R.C. Batra², M. Porfiri^{2, 3}

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.2, pp. 111-132, 2005, DOI:10.3970/cmes.2005.009.111

Abstract Structural health monitoring techniques based on vibration data have received increasing attention in recent years. Since the measured modal characteristics and the transient motion of a beam exhibit low sensitivity to damage, numerical techniques for accurately computing vibration characteristics are needed. Here we use a Meshless Local Petrov-Galerkin (MLPG) method to analyze vibrations of a beam with multiple cracks. The trial and the test functions are constructed using the Generalized Moving Least Squares (GMLS) approximation. The smoothness of the GMLS basis functions requires the use of special techniques to account for the slope discontinuities at the crack locations. Therefore, a… 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 >

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 >

Y. T. Gu, G. R. Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.4, pp. 463-476, 2001, DOI:10.3970/cmes.2001.002.463

Abstract A meshless method for the analysis of Kirchhoff plates based on the Meshless Local Petrov-Galerkin (MLPG) concept is presented. A MLPG formulation is developed for static and free vibration analyses of thin plates. Local weak form is derived using the weighted residual method in local supported domains from the 4th order partial differential equation of Kirchhoff plates. The integration of the local weak form is performed in a regular-shaped local domain. The Moving Least Squares (MLS) approximation is used to constructed shape functions. The satisfaction of the high continuity requirements is easily met by MLS interpolant, which is based on… More >