Two-dimensional large deformation analysis of hyperelastic and elasto-plastic solids based on the Meshless Local Petrov-Galerkin method (MLPG) is presented. A material configuration based the nonlinear MLPG formulation is introduced for the large deformation analysis of both path-dependent and path-independent materials. The supports of the MLS approximation functions cover the same sets of nodes during material deformation, thus the shape function needs to be computed only in the initial stage. The multiplicative hyperelasto-plastic constitutive model is adopted to avoid objective time integration for stress update in large rota- tion. With this constitutive model, the computational formulations for path-dependent and path-independent materials become identical. Computational efficiency of the nonlinear MLPG method is discussed and optimized in several aspects to make the MLPG an O(N) algorithm. The numerical examples indicate that the MLPG method can solve large deformation problems accurately. Moreover, the MLPG computations enjoy better convergence rate than the FEM under very large particle distortion.
In this paper, an adaptive boundary element method (BEM) is presented for solving 3-D elasticity problems. The numerical scheme is accelerated by the new version of fast multipole method (FMM) and parallelized on distributed memory architectures. The resulting solver is applied to the study of representative volume element (RVE) for short fiberreinforced composites with complex inclusion geometry. Numerical examples performed on a 32-processor cluster show that the proposed method is both accurate and efficient, and can solve problems of large size that are challenging to existing state-of-the-art domain methods.
A boundary collocation method based on the least-square technique and a corresponding adaptive computation process have been developed for the plate bending problem. The trial functions are constructed using a series of the biharmonic polynomials, and the local error indicators are given by the residuals of the energy density on the boundary. In comparison with the conventional collocation methods, the solution accuracy in the present method can be improved in an economical and efficient way. In order to demonstrate the efficiency and advantages of the adaptive boundary collocation method proposed in this paper, two numerical examples are presented for circular plates subjected to uniform loads and restrained by mixed boundary conditions. The numerical results for the examples show good agreement with ones presented in the literature.