Meshless Integral Method for Analysis of Elastoplastic Geotechnical Materials
The meshless integral method based on regularized boundary equation [1][2] is extended to analyze elastoplastic geotechnical materials. In this formulation, the problem domain is clouded with a node set using automatic node generation. The sub-domain and the support domain related to each node are also generated automatically using algorithms developed for this purpose. The governing integral equation is obtained from the weak form of elastoplasticity over a local sub-domain and the moving least-squares approximation is employed for meshless function approximation. The geotechnical materials are described by pressure-sensitive multi-surface Drucker-Prager/Cap plasticity constitutive law with hardening. A generalized collocation method is used to impose the essential boundary conditions and natural boundary conditions are incorporated in the system governing equations. A comparison of the meshless results with the FEM results shows that the meshless integral method is accurate and robust enough to solve geotechnical materials.