自然单元法是一种新兴的无网格数值计算方法，其实质是基于自然相邻插值（C∞）的伽辽金法．文中推导了基于Lasserre凸多面体体积公式的三维自然邻结点坐标及其导数的算法，给出了三维自然单元法算法的流程图．该算法实际上可以用于任意维数的自然单元法计算．对于Lasserre算法带来的多余约束问题，提出了2种可行的解决算法．经验证算例，三维自然单元法的计算结果精度同六面体单元有限元法相当．
