24小时热门版块排行榜

返回列表

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

water217

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 222
散金: 481
帖子: 140
在线: 82.8小时
虫号: 1973972
注册: 2012-09-04
性别: GG
专业: 机械结构强度学

[求助] 如何用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

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有11人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有16人回复
初秋的晨风已经有3人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有4人回复
面上没中，邀请各位路过的虫友分析一下分数已经有11人回复
Ei源刊怎么投已经有4人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有6人回复
两块石头已经有5人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

率性而为……

1楼 2013-01-24 16:54:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

water217

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 222
散金: 481
帖子: 140
在线: 82.8小时
虫号: 1973972
注册: 2012-09-04
性别: GG
专业: 机械结构强度学

引用回帖:

9楼: Originally posted by dhonghui at 2014-09-17 14:11:33
"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 ...

谢谢，您的回复太有价值了，将有有人看到本贴必将大有收获