24小时热门版块排行榜    

查看: 3883  |  回复: 1

william985

银虫 (小有名气)

[求助] 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]等进行运动方程计算。不知道理解的是否正确?求助
回复此楼
别老是眺望远方,脚下的路其实更重要
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuwenhua625

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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;
}

2楼2016-04-07 13:25:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 william985 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 5/250 2026-08-16 18:38 by lxmichacl
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 72/3600 2026-08-16 13:15 by 医学老男孩
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见