24小时热门版块排行榜

>论坛更新日志 (3600)
>虫友互识 (870)
>导师招生 (219)
>文献求助 (218)
>硕博家园 (105)
>休闲灌水 (94)
>考博 (76)
>博后之家 (51)
>基金申请 (45)
>找工作 (38)
>论文投稿 (38)
>招聘信息布告栏 (37)
>公派出国 (36)
>考研 (32)
>教师之家 (30)
>外文书籍求助 (26)

返回列表

william985

银虫 (小有名气)

应助: 4 (幼儿园)
金币: 253.8
散金: 46
红花: 1
帖子: 62
在线: 23.7小时
虫号: 3888664
注册: 2015-05-24
性别: GG
专业: 海洋工程

[求助] UDF函数中的Compute_Force_And_Moment的用法已有1人参与

想要利用Compute_Force_And_Moment调用每一个迭代步计算的力和力矩信息（刚体运动，下一步要调用整个刚体所受的力与力矩），不太了解这个函数的用法。
利用它求出的是不是整个定义的动区域物体所受合力或力矩，还是整个物体所受的力和力矩需要利用cell_loop或者face_loop来循环迭加求出？
之前只找到一个小程序，个人觉得求出来的是整个动区域的物体所受合力，合力矩，直接以Compute_Force_And_Moment(d, t, CG, force, moment, FALSE)的形式调用该函数，计算结果存储到force和moment向量中，下一步直接用force[0]等进行运动方程计算。不知道理解的是否正确？求助

回复此楼

» 猜你喜欢

AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有3人回复
有院领导为了换新车，用横向课题经费买了俩车已经有9人回复
酰胺脱乙酰基已经有13人回复
博士延得我，科研能力直往上蹿已经有8人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有4人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有10人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有5人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有6人回复
遇见不省心的家人很难过已经有22人回复

别老是眺望远方,脚下的路其实更重要

1楼 2016-04-06 22:42:16

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

xuwenhua625

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 203.7
散金: 20
红花: 2
帖子: 216
在线: 118.2小时
虫号: 4122526
注册: 2015-10-06
性别: GG
专业: 流体力学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
william985: 金币+10, ★★★很有帮助, 昨天搜了一下,发现这个函数的帖子还是有几个讲的很不错的,感谢帮助 2016-04-07 22:38:26

这是我做过的一个case，可以参考一下

CODE:

DEFINE_CG_MOTION(foilright,dt,vel,omega,time,dtime)

{

   real x_A=0.0,y_A=0.0,x_N=0.0,y_N=0.0;

   real theta_A,theta_N;

   real omegap = 0.0;

   real presure;

   real x[ND_ND];

/////////////////////////////////////////////////////////////////////////////

    face_t f ;

    cell_t c ;

        real f_glob[ND_ND],m_glob[ND_ND],x_cg[ND_ND];

    real f1,f3,f5;

    Domain *domain= Get_Domain (1);

    Thread *tf1 = Lookup_Thread (domain, 33);

/////////////////////////////////////////////////////////////////////////////

   theta_A = -theta0*cos(w*CURRENT_TIME);

   theta_N = -theta0*cos(w*(CURRENT_TIME+CURRENT_TIMESTEP));

   omegap = (theta_N-theta_A)/CURRENT_TIMESTEP;

   x_A= Radius*cos(w*CURRENT_TIME);

   x_N= Radius*cos(w*(CURRENT_TIME+CURRENT_TIMESTEP));

   y_A= Radius*sin(w*CURRENT_TIME);

   y_N= Radius*sin(w*(CURRENT_TIME+CURRENT_TIMESTEP));

   vel[0] =(x_N - x_A) / CURRENT_TIMESTEP;

   vel[1] =(y_N - y_A) / CURRENT_TIMESTEP;

   omega[0] =omega[1] =0.0; 

   omega[2] = omegap;

   // f1  f3  f5 output /////////////////////////////////////////////////////

   if(time>=dtime){

   x_cg[0]=DT_CG(dt)[0];

    x_cg[1]=DT_CG(dt)[1];

    Compute_Force_And_Moment (domain, tf1, x_cg, f_glob, m_glob,1);

    f1 = f_glob[0]/(0.5*rou*cord*U*U);

    f3 = f_glob[1]/(0.5*rou*cord*U*U);

    f5 = m_glob[2]/(0.25*rou*cord*cord*U*U);

    fpright=fopen("foil_right.dat", "a+");

    fprintf(fpright,"%.6f %.5f %.5f %.5f\n",CURRENT_TIME*w/(2*PI),f1,f3,f5);

        fclose(fpright);}

   //////////////////////////////////////////////////////////////////////////

  return;

}

赞一下