| 查看: 2614 | 回复: 5 | ||
[求助]
根据阀两侧压力定义阀门的运动
|
|
求助：根据阀两侧压力的大小定义阀的运动 根据如下的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  运动终止网格.png  7.png  8.png  9.png  位置监视.png  压力1.png  压力2.png@wuming524 |
回复此楼» 猜你喜欢
面上本子正文33页，违规吗？会被低分嘛？
已经有10人回复
-
国自然上会要求
已经有9人回复
-
博士申请
已经有5人回复
-
评审有感
已经有6人回复
-
今年审到国自然15份，谈谈感受
已经有16人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
2楼2017-11-15 08:25:57
3楼2018-01-25 09:48:41
|
据我所知，如果使用动网格的话至少是要在阀座和运动的阀之间有一层网格的，否则会出现负体积，当然你可以把这层网格划的很小，近似看成闭合的。目前并没有查到UDF可以使interior在计算过程中改成wall的命令。 发自小木虫IOS客户端 |
4楼2018-03-29 12:48:21
5楼2018-04-01 22:24:37
6楼2018-04-12 14:02:06
