一种连续非连续模型任意形状的光滑接触势计算方法

本发明公开了一种连续非连续模型任意形状的光滑接触势计算方法，包括如下步骤：S1、构建连续‑非连续耦合仿真模型，初始化模型参数，进行接触对判别，识别并存储存在可能性接触关系的单元对，建立接触单元对的数据列表；S2、对于每一接触单元对，进行主、从单元分配，并判断接触单元对是否发生重叠；S3、对于发生重叠的接触单元对，记录主单元嵌入从单元中的顶点坐标，该顶点称为嵌入顶点，建立嵌入顶点的权重函数模型；S4、基于步骤S3中建立的权重函数模型，计算嵌入顶点于从单元中的光滑势函数；S5、根据步骤S4计算得到的光滑势函数，计算从单元对主单元产生的法向接触力，并进一步将计算的接触力按照形函数分配为主、从单元的等效节点力。