关于球面轴承转子运动中动网格问题

红色石头

我在做一个球面轴承，要模拟轴承转子的运动轨迹，轴承的转子的转速可设置在1万左右，想利用6DOF来实现，编的程序为
#include<stdio.h>
#include "udf.h"

DEFINE_SDOF_PROPERTIES(sdof_properties, prop, dt, time, dtime)
{   prop[SDOF_MASS]=0.08;         
        prop[SDOF_IXX]        = 2.019e-3;       
    prop[SDOF_IXY]        =  0;
        prop[SDOF_IXZ]        =  0;
        prop[SDOF_IYY]        =  2.021e-3;
        prop[SDOF_IYZ]        =  7.2766e-6;
        prop[SDOF_IZZ]        =  3.525e-3;
        prop[SDOF_ZERO_TRANS_Z] = TRUE;
       
}
但在运算的过程中，转子的轴心轨迹在半圆的时候网格就出现负体积，而导致无法计算，哪位大神能帮忙解决一下，或者有什么更好的方法来模拟这个问题，非常感谢！！