Jin Ma, Yang Liu, Hongbing Lu, Ranga Komanduri¹

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

Abstract A multiscale simulation technique coupling three scales, namely, the molecular dynamics (MD) at the atomistic scale, the discrete dislocations at the meso scale and the generalized interpolation material point (GIMP) method at the continuum scale is presented. Discrete dislocations are first coupled with GIMP using the principle of superposition (van der Giessen and Needleman (1995)). A detection band seeded in the MD region is used to pass the dislocations to and from the MD simulations (Shilkrot, Miller and Curtin (2004)). A common domain decomposition scheme for each of the three scales was implemented for parallel More >

J. Ma², H. Lu², B. Wang², R. Hornung³, A. Wissink³, R. Komanduri^2,*

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 101-118, 2006, DOI:10.3970/cmes.2006.014.101

Abstract A new method for multiscale simulation bridging two scales, namely, the continuum scale using the generalized interpolation material point (GIMP) method and the atomistic scale using the molecular dynamics (MD), is presented and verified in 2D. The atomistic strain from the molecular dynamics simulation is determined through interpolation of the displacement field into an Eulerian background grid using the same generalized interpolation functions as that in the GIMP method. The atomistic strain is consistent with that determined from the virial theorem for interior points but provides more accurate values at the boundary of the MD… More >

Jin Ma, Hongbing Lu, Ranga Komanduri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.3, pp. 213-228, 2006, DOI:10.3970/cmes.2006.012.213

Abstract The generalized interpolation material point (GIMP) method, recently developed using a C¹ continuous weighting function, has solved the numerical noise problem associated with material points just crossing the cell borders, so that it is suitable for simulation of relatively large deformation problems. However, this method typically uses a uniform mesh in computation when one level of material points is used, thus limiting its effectiveness in dealing with structures involving areas of high stress gradients. In this paper, a spatial refinement scheme of the structured grid for GIMP is presented for simulations with highly localized stress gradients.… More >

Vladimir Sladek¹, Jan Sladek¹, Chuanzeng Zhang²

CMC-Computers, Materials & Continua, Vol.4, No.3, pp. 177-188, 2006, DOI:10.3970/cmc.2006.004.177

Abstract This paper concerns the stability, convergence of accuracy and cost efficiency of four various formulations for solution of boundary value problems in non-homogeneous elastic solids with functionally graded Young's modulus. The meshless point interpolation method is employed with using various basis functions. The interaction among the elastic continuum constituents is considered in the discretized formulation either by collocation of the governing equations or by integral satisfaction of the force equilibrium on local sub-domains. The exact benchmark solutions are used in numerical tests. More >

J. Sladek¹, V. Sladek¹, Ch. Zhang², F. Garcia-Sanche³, M. Wünsche²

CMC-Computers, Materials & Continua, Vol.4, No.2, pp. 109-118, 2006, DOI:10.3970/cmc.2006.004.109

Abstract Piezoelectric materials have wide range engineering applications in smart structures and devices. They have usually anisotropic properties. Except this complication electric and mechanical fields are coupled each other and the governing equations are much more complex than that in the classical elasticity. Thus, efficient computational methods to solve the boundary or the initial-boundary value problems for piezoelectric solids are required. In this paper, the Meshless local Petrov-Galerkin (MLPG) method with a Heaviside step function as the test functions is applied to solve two-dimensional (2-D) piezoelectric problems. The mechanical fields are described by the equations of… More >

J. Sladek¹, V. Sladek, Ch.Zhang²

Structural Durability & Health Monitoring, Vol.1, No.2, pp. 131-144, 2005, DOI:10.3970/sdhm.2005.001.131

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for crack analysis in two-dimensional (2-d), anisotropic and linear elastic solids with continuously varying material properties. Both quasi-static and transient elastodynamic problems are considered. For time-dependent problems, the Laplace-transform technique is utilized. A unit step function is used as the test function in the local weak-form. It is leading to local boundary integral equations (LBIEs) involving only a domain-integral in the case of transient dynamic problems. The analyzed domain is divided into small subdomains with a circular shape. The moving least-squares (MLS) method is More >

V. Sladek¹, J. Sladek¹, M. Tanaka²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.1, pp. 69-84, 2005, DOI:10.3970/cmes.2005.007.069

Abstract An efficient numerical method is proposed for 2-d potential problems in anisotropic media with continuously variable material coefficients. The method is based on the local integral equations (utilizing a fundamental solution) and meshfree approximation of field variable. A lot of numerical experiments are carried out in order to study the numerical stability, accuracy, convergence and efficiency of several approaches utilizing various interpolations. 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… More >

S. G. Bardenhagen^1,2, E. M. Kober³

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.6, pp. 477-496, 2004, DOI:10.3970/cmes.2004.005.477

Abstract The Material Point Method (MPM) discrete solution procedure for computational solid mechanics is generalized using a variational form and a Petrov–Galerkin discretization scheme, resulting in a family of methods named the Generalized Interpolation Material Point(GIMP) methods. The generalizationpermits identification with aspects of other point or node based discrete solution techniques which do not use a body–fixed grid, i.e. the “meshless methods”. Similarities are noted and some practical advantages relative to some of these methods are identified. Examples are used to demonstrate and explain numerical artifact noise which can be expected inMPM calculations. Thisnoiseresultsin non-physical local More >

Yoshihiro Ochiai¹, Vladimir Sladek²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.6, pp. 525-536, 2004, DOI:10.3970/cmes.2004.006.525

Abstract The conventional boundary element method (BEM) uses internal cells for the domain integralsCwhen solving nonlinear problems or problems with domain effects. This paper is concerned with conversion of the domain integral into boundary ones and some non-integral terms in a three-dimensional BIEM, which does not require the use of internal cells. This method uses arbitrary internal points instead of internal cells. The method is based on a three-dimensional interpolation method in this paper by using a polyharmonic function with volume distribution. In view of this interpolation method, the three-dimensional numerical integration is replaced by boundary More >