声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2171|回复: 0

[Fluent应用] 请教UDF问题!

[复制链接]
发表于 2006-6-5 19:13 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
请问在UDF的CG_MOTION宏中作用在移动的面上的力是这个面上的力还是压强,如果说是压力是作用在每个网格上的压力还是整个面上的力。<BR>下面这个程序是已经毕业的师兄写的,它上面的600师兄说是移动面上的网格数,可是我看帮助上并没有这个说法。<BR>这个程序是个压缩弹簧的的运动的程序,0.88是弹簧的刚度,4.88是弹簧的预压量!<BR>请问各位高手!<BR>#include "udf.h"<BR>#include "dynamesh_tools.h"<BR><BR>static real v_prev=0;<BR>static real S=0.00001;<BR><BR>DEFINE_CG_MOTION(coopermove,dt,vel,omega,time,dtime)<BR>{<BR>  Thread *t;<BR>  float mass=0.4e-3;<BR>  real force;<BR>  real spring=0.88;<BR>  face_t f;<BR>  real dv;<BR>  real ds;<BR><BR>  ds = 0.2*dtime;<BR>  S+=ds;<BR>  force = spring * (4.88-S*1000)/600;<BR><BR>  NV_S(vel, =, 0.0);<BR>  <BR>  t=DT_THREAD((Dynamic_Thread*)dt);<BR>  <BR>  dv = dtime * (force/mass);<BR>  v_prev += dv;<BR>  <BR>  vel[2] = v_prev;<BR>}<BR><BR>DEFINE_CG_MOTION(tuiganmove,dt,vel,omega,time,dtime)<BR>{<BR>  NV_S(vel, =, 0.0);<BR>  vel[2] = v_prev;<BR>}
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 21:53 , Processed in 0.090468 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表