MESHLESS METHOD FOR MODELING LARGE DEFORMATION WITH ELASTOPLASTICITY

Over the past two decades meshless methods have attracted much attention owing to their advantages in adaptivity, higher degree of solution field continuity, and capability to handle moving boundary and changing geometry. In this work, a meshless integral method based on the regularized boundary integral equation has been developed and applied to two-dimensional linear elasticity and elastoplasticity with small or large deformation. The development of the meshless integral method and its application to two-dimensional linear elasticity is described first. The governing integral equation is obtained from the weak form of elasticity over a local sub-domain, and the moving least-squares approximation is employed for meshless function approximation. This formulation incorporates: a subtraction method for singularity removal in the boundary integral equation, a special numerical integration for the calculation of integrals with weak singularity which further improves accuracy, a collocation method for the imposition of essential boundary conditions, and a method for incorporation of natural boundary conditions in the system governing equation. Next, elastoplastic material behavior with small deformation is introduced into the meshless integra