摘要：势函数离散单元法在接触力计算中存在势函数物理意义不明确、接触力计算与客观实际不符等重大缺陷。在已有研究的基础上,提出一种新型任意形状块体距离势离散单元法。该方法将势函数定义为接触体间归一化的距离函数,可以客观衡量单元间的嵌入程度,具有明确的物理意义,真实反映单元的接触刚度,使法向接触力计算更为合理。此外,新方法重新定义了切向方向矢量,提出切向接触力计算方法。新方法保留原势函数法所有优点,同时完善接触力计算受单元形式影响的不足,并突破三角形单元的限制,能够应用于任意凸多边形单元。通过若干算例说明新方法的正确性和有效性,可用于精确描述复杂非连续介质的运动过程。

关键词：数值分析 离散单元法 任意多边形块体单元 距离势函数 接触力 

单位：河海大学水利水电学院; 江苏南京210098; 天津大学水利工程仿真与安全国家重点实验室; 天津300072