[发明专利]一种基于距离势函数可变形三维凸多面体块体离散单元法

摘要

本发明公开了一种基于距离势函数可变形三维凸多面体块体离散单元法，步骤一，选取研究对象，建立可变形离散单元系统；步骤二，确定可变形离散单元系统的时间步长Δt；步骤三，计算作用于接触单元与目标单元之间总的接触力；步骤四，将步骤三计算得出作用在接触单元以及目标单元上的总的接触力用形函数转化成载荷的等效节点力矢量R；步骤五，由步骤四计算得出的载荷的等效节点力矢量R，计算得出下一时刻t+Δt每个网格单元的位移，速度，以及加速度；步骤六，根据步骤五中网格单元的位移，更新下一时刻t+Δt每个网格单元节点坐标。本发明使数值模拟更加符合工程实际，可以精确地捕捉离散体系的运动过程，精确反应离散单元内部的真实应力和变形状态。

主权项

1.一种基于距离势函数三维可变形凸多面体块体离散单元法，其特征在于：包括以下步骤：步骤一，选取研究对象，建立可变形离散单元系统；步骤二，确定可变形离散单元系统的时间步长Δt；步骤三，在当前时刻t，采用No Binary Method接触检测方法对离散单元外围一层的网格单元进行接触检测，将其中一个离散单元外围一层产生接触的网格单元定义为目标单元，与之接触的离散单元外围一层的网格单元则定义为接触单元，并且根据距离势函数的定义，计算作用于接触单元与目标单元之间总的接触力；步骤四，将步骤三计算得出作用在接触单元以及目标单元上的总的接触力用形函数转化成载荷的等效节点力矢量步骤五，由步骤四计算得出的载荷的等效节点力矢量计算得出下一时刻t+Δt每个网格单元的位移，速度，以及加速度；步骤六，根据步骤五中网格单元的位移，更新下一时刻t+Δt每个网格单元节点坐标。