楔形体自由落体进入粘性较大流体，流体被吸起来！！

soundyrao

最近在做楔形体自由落体进入油脂的问题，油脂的粘度很大，粘度在800左右，楔形体的重量为150g，初始时锥体尖端刚好与油脂表面接触，然后模拟锥体在自重作用下，刺入油脂，看锥体在5s时间内刺入锥体的深度。按照想法，做了数值模拟。发现很奇怪的是，把油脂粘度设为800，才刚计算，锥体不但没有刺入油脂，反而油脂被吸起来了（图1），之后把油脂粘度减到200，初始计算都比较合理，锥体逐步下降，锥体速度开始增大，之后减小，但在计算到4s时，锥体速度发生震荡，实在不明白，求各位大神指教~~~
油脂粘度800时，锥体图像见图1；粘度200，锥体下落4s，锥体图像间图2，锥体速度随流动时间变化间图3.模拟case文件请见附件

soundyrao

欢迎大家提意见，觉得有可能的原因，都请告诉我，谢谢

soundyrao

大神们，求指导~~~~~

baronsooo

Mark

598459670

你的UDF呢？

soundyrao

598459670 发表于 2017-5-25 08:18
你的UDF呢？

DEFINE_SDOF_PROPERTIES(falling_pyramid, prop, dt, time, dtime)
{
   prop[SDOF_MASS] = 0.15;
   prop[SDOF_ZERO_ROT_Z] = TRUE;
   prop[SDOF_ZERO_TRANS_X] = TRUE;
   printf ("\n2d_falling_pyramid: Updated 6DOF properties");
}

UDF是这样的，之所以限定Z的转动和X的位移是因为，如果不限定，锥体会发生转动

598459670

soundyrao 发表于 2017-5-25 09:38
DEFINE_SDOF_PROPERTIES(falling_pyramid, prop, dt, time, dtime)
{
   prop[SDOF_MASS] = 0.15;

不好意思，回晚了，你没有设置转动惯量吧