Y. H. Liu^1,2, S. S. Chen¹, J. Li¹, Z. Z. Cen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.3, pp. 145-156, 2008, DOI:10.3970/cmes.2008.031.145

Abstract A novel meshless method for structural dynamic analysis is presented and discussed in this paper. It is called meshless local natural neighbour interpolation (MLNNI) method, which uses a meshless spatial approximation based only on nodes. The MLNNI is derived from the generalized meshless local Petrov-Galerkin (MLPG) method as a special case. Local weak forms are developed using weighted residual method locally from the dynamic partial differential equation. In the construction of trial functions, the natural neighbour interpolation (NNI) is employed to simplify the treatment of the essential boundary conditions. The domain integration is evaluated over included Delaunay triangles in each… 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 with that by FEM software-ANSYS.… More >

S. S. Chen¹, Y. H. Liu^1,2, Z. Z. Cen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.1, pp. 39-56, 2008, DOI:10.3970/cmes.2008.028.039

Abstract In most engineering applications, solutions derived from the lower-bound theorem of plastic limit analysis are particularly valuable because they provide a safe estimate of the load that will cause plastic collapse. A solution procedure based on the meshless local Petrov-Galerkin (MLPG) method is proposed for lower-bound limit analysis. This is the first work for lower-bound limit analysis by this meshless local weak form method. In the construction of trial functions, the natural neighbour interpolation (NNI) is employed to simplify the treatment of the essential boundary conditions. The discretized limit analysis problem is solved numerically with the reduced-basis technique. The self-equilibrium… More >

Kai Wang¹, Shenjie Zhou^1,2, Zhifeng Nie¹, Shengli Kong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.2, pp. 107-122, 2008, DOI:10.3970/cmes.2008.026.107

Abstract The natural neighbour Petrov-Galerkin method (NNPG) is one of the special cases of the generalized meshless local Petrov-Galerkin method (MLPG). This paper demonstrates the NNPG can be successfully used in design sensitivity analysis in 2D elasticity. The design sensitivity analysis method based on the local weak form (DSA-LWF) in the NNPG context is proposed. In the DSA-LWF, the local weak form of governing equation is directly differentiated with respect to design variables and discretized with NNPG to obtain the sensitivities of structural responds. The calculation of derivatives of shape functions with respect to design variables is avoided. No background meshes… More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 75-90, 2008, DOI:10.3970/cmes.2008.023.075

Abstract In the MLPG_R (Meshless Local Petrove-Galerkin based on Rankine source solution) method, one needs a meshless interpolation scheme for an unknown function to discretise the governing equation. The MLS (moving least square) method has been used for this purpose so far. The MLS method requires inverse of matrix or solution of a linear algebraic system and so is quite time-consuming. In this paper, a new scheme, called simplified finite difference interpolation (SFDI), is devised. This scheme is generally as accurate as the MLS method but does not need matrix inverse and consume less CPU time to evaluate. Although this scheme… More >

J. Sladek¹, V. Sladek¹, Ch. Zhang², P. Solek³

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 217-234, 2007, DOI:10.3970/cmes.2007.022.217

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for the solution of boundary value problems for coupled thermo-electro-mechanical fields. Transient dynamic governing equations are considered here. To eliminate the time-dependence in these equations, the Laplace-transform technique is applied. Material properties of piezoelectric materials are influenced by a thermal field. It is leading to an induced nonhomogeneity and the governing equations are more complicated than in a homogeneous counterpart. Two-dimensional analyzed domain is subdivided into small circular subdomains surrounding nodes randomly spread over the whole domain. A unit step function is used as the test functions in the… More >

J. Sladek¹, V. Sladek¹, Ch. Zhang², P. Solek³, L. Starek³

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 247-262, 2007, DOI:10.3970/cmes.2007.019.247

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for crack analysis in two-dimensional (2-D) and three-dimensional (3-D) axisymmetric piezoelectric solids with continuously varying material properties. Axial symmetry of geometry and boundary conditions reduces the original 3-d boundary value problem into a 2-d problem. Stationary problems are considered in this paper. The axial cross section is discretized into small circular subdomains surrounding nodes randomly spread over the analyzed domain. A unit step function is used as the test functions in the local weak-form. Then, the derived local integral equations (LBIEs) involve only contour-integrals on the surfaces of subdomains.… More >

Amit Shaw¹, D Roy²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 69-98, 2007, DOI:10.3970/cmes.2007.019.069

Abstract A novel discretization strategy and derivative reproduction based on reproducing kernel (RK) particle approximations of functions are proposed. The proposed scheme is in the form of an RK interpolation that offers significant numerical advantages over a recent version of the strategy by Chen et al. (2003), wherein the authors added a set of primitive functions to the reproducing kernel (enrichment) functions. It was also required that the support size of the primitive function be less than the smallest distance between two successive grid points. Since the primitive function was required to vary from 0 to 1 within half of this… More >

T. Jarak¹, J. Sorić¹, J. Hoster¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 235-246, 2007, DOI:10.3970/cmes.2007.018.235

Abstract A meshless computational method based on the local Petrov-Galerkin approach for the analysis of shell structures is presented. A concept of a three dimensional solid, allowing the use of completely 3-D constitutive models, is applied. Discretization is carried out by using both a moving least square approximation and polynomial functions. The exact shell geometry can be described. Thickness locking is eliminated by using a hierarchical quadratic approximation over the thickness. The shear locking phenomena in case of thin structures and the sensitivity to rigid body motions are minimized by applying interpolation functions of sufficiently high order. The numerical efficiency of… More >

J. Ma¹, H. Lu¹, B. Wang¹, S. Roy¹, R. Hornung², A. Wissink², R. Komanduri^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 135-152, 2005, DOI:10.3970/cmes.2005.008.135

Abstract In the simulation of a wide range of mechanics problems including impact/contact/penetration and fracture, the material point method (MPM), Sulsky, Zhou and Shreyer (1995), demonstrated its computational capabilities. To resolve alternating stress sign and instability problems associated with conventional MPM, Bardenhagen and Kober (2004) introduced recently the generalized interpolation material point (GIMP) method and implemented for one-dimensional simulations. In this paper we have extended GIMP to 2D and applied to simulate simple tension and indentation problems. For simulations spanning multiple length scales, based on the continuum mechanics approach, we present a parallel GIMP computational method using the Structured Adaptive Mesh… More >