三维弹塑性自然单元法算法实现

自然单元法是一种新兴的无网格数值计算方法，其实质是基于自然相邻插值（C∞）的伽辽金法。该方法计算精度与四边形或六面体单元有限元法相当，自然相邻插值函数比其他无网格法插值函数的计算速度快。由于自然相邻插值在凸域的边界上的相邻点之间是严格线性的．所以自然单元法在边界面的处理也相当简单。本文研究了在自然单元法中采用Von．Mises，Mohr-Coulomb和Drucker-Prager屈服准则解决三维弹塑性问题，并编制了相应计算程序，最后通过算例验证算法的正确性