| 查看: 6801 | 回复: 12 | |||
water217银虫 (小有名气)
|
[求助]
如何用Compute_Force_And_Moment函数求流体对运动翼型的俯仰力矩 已有1人参与
|
||
|
使用fluent的动网格功能模拟了流体绕俯仰-升沉振荡二维翼型的流场后，想得到作用于翼型俯仰参考轴上的俯仰力矩随时间的变化关系，采用了函数Compute_Force_And_Moment(domain, thread, body_centroid, force, moment, FALSE);来求解，可是算下来的俯仰力矩明显有问题（图片中的俯仰力矩都是正值，可是俯仰振荡以后攻角有正有负，俯仰力矩肯定有负值的情况），更换模型俯仰轴的参考位置后，结果也没有什么变化（一个参考轴取距前缘0.3倍的弦长，一个取0.6倍的弦长，结果峰值差不多，变化规律基本一样），请问Compute_Force_And_Moment默认的力矩参考轴是什么位置？我的程序为什么输出的力矩系数有问题？ 这个程序里求得的力矩是相对于，大地坐标系，还是随体坐标系，谢谢 附程序： #include "udf.h" #include "mem.h" #include "dynamesh_tools.h" #include "math.h" /************************************************************/ FILE*fp; real current_time = 0.0000 ; static real vel[3]={0,0,0};/* 初始速度*/ static real omega[3]={0,0,0};/* 初始角速度*/ static real body_centroid[3]={0,0,0};/* 初始形心位置,需要修改*/ static real Chord=0.1,F_density=998.2,U_flow=2; static real pi=3.1416,frequency=2.0000,Y_amplitude=0.1,w_amplitude=0.8727; real NV_VEC( origin ),NV_VEC( force ),NV_VEC( moment ); real C_momentr; /************************************************************/ DEFINE_CG_MOTION(hydrofoil,dt,vel,omega,time,dtime) { Domain * domain ; Thread * thread ; face_t f; current_time = CURRENT_TIME ; vel[0] = 0.0 ; vel[1] =2*pi*frequency*Y_amplitude*cos(2*pi*frequency*current_time+pi/2); vel[2] = 0.0 ; omega[0] = 0.0 ; omega[1] = 0.0 ; omega[2] =2*pi*frequency*w_amplitude*cos(2*pi*frequency*current_time); thread=DT_THREAD(dt); domain = THREAD_DOMAIN (thread); /************************************************************/ body_centroid[0]=DT_CG(dt)[0];/* 获取形心坐标*/ body_centroid[1]=DT_CG(dt)[1];/* 获取形心坐标*/ Compute_Force_And_Moment(domain, thread, body_centroid, force, moment, FALSE);/* 计算水翼总的受力情况*/ C_moment=moment[2]/(0.5*F_density*U_flow*U_flow*Chord*Chord);/* 计算俯仰力矩系数*/ fp=fopen("D:\\moment_coefficient.dat", "a+"  ;fprintf(fp, "%12f",C_moment); fprintf(fp, "\n" ; fclose(fp); }  未命名.JPG [ Last edited by water217 on 2013-1-24 at 17:55 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 CFD_study |
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
国自然面上复盘~欢迎讨论
已经有5人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- compute force and moment宏
已经有7人回复
- 固体在流体中的运动，可以用什么软件来模拟和分析？
已经有4人回复
- PROE运动仿真问题
已经有3人回复
dhonghui
银虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1130.1
- 红花: 1
- 帖子: 78
- 在线: 23.2小时
- 虫号: 2015753
- 注册: 2012-09-20
- 性别: GG
- 专业: 流体力学
|
"This function needs size of 3 for force and moment. It works in parallel and the common arguments are the of the same type as your arguments. In addition you have to pass the domain as first arg and a boolean as last argument. This boolean has to be TRUE if you also call the function on the host else it has to be FALSE." Compute_Force_And_Moment (domain, tf, CG, force, moment, FALSE); domain = THREAD_DOMAIN (DT_THREAD ((Dynamic_Thread *)dt)); - pointer on Domain tf = Lookup_Thread(domain, ID); - pointer on the Thread in Domain I don't know if "force" includes viscose and pressure force or only pressure force!!! ************************************************** ************** compute viscose force in UDF you have to define in mem.h header or in local c-file this macro: #define F_VISCOSE(f,tf) F_STORAGE_R_N3V(f,tf,SV_WALL_SHEAR) Thread *tf = DT_THREAD(dt); /*pointer on thread*/ face_t f; /*number of faces*/ begin_f_loop(f, tf) { NV_VS(F_VISC, += , F_VISCOSE(f, tf), *, -1.0 * k); /* compute viscose force vector. k is a factor for depth in 2D computation and -1.0 corrects the sign of viscose force*/ } end_f_loop(f, tf); from http://www.cfd-online.com/Forums ... rce_and_moment.html |
9楼2014-09-17 14:11:33
dragonship
铁杆木虫 (正式写手)
- 应助: 159 (高中生)
- 金币: 5105.8
- 散金: 5
- 红花: 50
- 帖子: 844
- 在线: 196.3小时
- 虫号: 2193451
- 注册: 2012-12-18
- 性别: GG
- 专业: 海洋工程
2楼2013-01-30 09:32:36
3楼2013-02-02 11:27:56
4楼2013-02-05 22:31:35
water217
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 222
- 散金: 481
- 帖子: 140
- 在线: 82.8小时
- 虫号: 1973972
- 注册: 2012-09-04
- 性别: GG
- 专业: 机械结构强度学
5楼2013-02-20 17:54:58
water217
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 222
- 散金: 481
- 帖子: 140
- 在线: 82.8小时
- 虫号: 1973972
- 注册: 2012-09-04
- 性别: GG
- 专业: 机械结构强度学
6楼2013-02-21 09:30:32
universewang
金虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 3685.7
- 散金: 12
- 红花: 2
- 帖子: 773
- 在线: 216.1小时
- 虫号: 2453924
- 注册: 2013-05-08
- 性别: GG
- 专业: 制造系统与自动化
7楼2014-03-06 16:35:29
8楼2014-07-03 15:25:06
water217
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 222
- 散金: 481
- 帖子: 140
- 在线: 82.8小时
- 虫号: 1973972
- 注册: 2012-09-04
- 性别: GG
- 专业: 机械结构强度学
10楼2014-09-20 16:27:00
投票:
