用XFEM建模移动不连续

各位机械师好，万博体育平台

我试图用XFEM建模不连续(晶界)的运动，因为它节省了网格适应所需的计算费用。不连续面两侧的物质行为是不同的。在任何一步之后，不连续都可能穿过一个元素。然后，为了计算元素的刚度，我必须对分割的区域分别积分，然后将它们相加。这在一维情况下很简单但在二维和三维情况下似乎很难。由于不连续面可以是任意形状，我必须找出元素中位于它两侧的点，对于这些点我可以使用符号距离函数。但是我如何在ABAQUS中实现这个。我有一个模糊的想法，我可能不得不使用UEL或UMAT。如果有人遇到类似的情况，请帮助我。

最好的问候,

爱