24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jeework

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 14
在线: 10.9小时
虫号: 4887102
注册: 2016-08-01
专业: 工程热物理相关交叉领域

[求助] 根据阀两侧压力定义阀门的运动

求助：根据阀两侧压力的大小定义阀的运动
根据如下的UDF和图片上的条件定义阀门的运动，并利用pand变量保存阀体左右两侧的压差。但是计算过程中出现了压差在某一个时间步左侧大于右侧，下一个时间步右侧大于左侧的情况如图中使用message语句显示的红线部分的内容，在图片上的压力云图也可看出。导致最终阀门只向右走了不到一个网格，然后循环地向左运动又向又运动的情况，从图片上监视阀门左侧坐标位置的surf-mon-1也可看出阀门在第N个时间步长内向左动，N+1个时间步长内又向右运动，如此循环往复，而实际情况应该是左侧的气缸活塞加压，压力逐渐变大，阀门应该一直向右侧运动才是。请大家帮忙看下问题出在哪了，谢谢！
附：UDF如下：
#include "udf.h"
#define valvel_id 112
#define valver_id 19
#define valve_stop_l 0.2958 /*阀门运动的左死点*/
#define valve_stop_r 0.3087 /*阀门运动的右死点*/
real pand,val_left,val_right;
DEFINE_ADJUST(valveposition, domain)
{ cell_t c;
  face_t f;
  real ad[ND_ND];
  real ae[ND_ND];
  double force_l,force_r,z01,z02;
  Thread *t5 = Lookup_Thread(domain, valvel_id);
  Thread *t6 = Lookup_Thread(domain, valver_id);
   pand=0.0;
   z01=0.0;
   z02=0.0;
   force_l=0.0;
       begin_f_loop (f, t5)
   {    F_CENTROID(ad,f,t5);
      val_left=ad[0];
      force_l=force_l+F_P(f,t5)*2*M_PI*0.001*ad[1]; /*0.001为阀侧单个网格高度*/
   }
   end_f_loop (f, t5)
   z01=force_l+0.04909*9.8;
   force_r=0.0;
       begin_f_loop (f, t6)
   {    F_CENTROID(ae,f,t6);
      val_right=ae[0];
      force_r=force_r+F_P(f,t6)*2*M_PI*0.001*ae[1]; /*0.001为阀侧单个网格高度*/
   }
   end_f_loop (f, t6)
   z02=force_r;
   pand=(z02-z01)/z01;
   Message( "pand is %g", pand);
}

DEFINE_CG_MOTION(banmove,dt,vel,omega,time,dtime)
{  /* reset velocities */
  NV_S(vel, =, 0.0);
  NV_S(omega, =, 0.0);
if(val_left>=valve_stop_l && pand>0.02)
   {vel[0]=-10.0;
   Message( "move left pand is %g", pand);
   }
else if (val_right<=valve_stop_r && pand<-0.02)
   {vel[0]=10.0;
   Message( "move right pand is %g", pand);
   }
else
   {vel[0]=0;
   Message( "stationary pand is %g", pand);
   }
}

起始网格.png