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

» 猜你喜欢

寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

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

率性而为……

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

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

water217

银虫 (小有名气)

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

引用回帖:

3楼: Originally posted by watcher007 at 2013-02-02 11:27:56
楼主你Compute_Force_And_Moment(domain, thread, body_centroid, force, moment, FALSE)这个里面的为什么是用FALSE，我一直用的是TRUE，我刚才用一个简单的算例算了一下这个里面的force好像没错

我记得我开始用过true 好像是出错了然后用false就没错了现在也不知道当时错误到底是不是true or false的问题。