函數(shù)Compute_Force_And_Moment

2017-01-16  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

FLUENT中可以調(diào)用函數(shù) Compute_Force_And_Moment()進(jìn)行力和力矩的計(jì)算。該函數(shù)的聲明位于f_wall.h頭文件內(nèi)。沒有任何對(duì)該函數(shù)的官方什么。

該函數(shù)的描述:

Compute_Force_And_Moment (d, t, x_cg, f_body, m_body, TRUE)

fluent軟件中一個(gè)內(nèi)部函數(shù),主要是用以計(jì)算三維物體的受力。其中t指向所要求解的物體邊界,x_cg是物體邊界的中心位置,f_body是物體所受水作用力,m_body是物體所受外力對(duì)于中心位置的力矩。

函數(shù)說明:

1. 這個(gè)宏不需要執(zhí)行循環(huán)操作;

2. domain是計(jì)算區(qū)域,tf1是面索,也就是你要計(jì)算受力的面,x_cg是面的形心位置,f_glob與m_glob分別指在慣性系下的力與力矩分量;

3. 該宏既可以計(jì)算3D受力,也可用于2D受力計(jì)算;

4. TRUE與FALSE是邏輯類型的量,TRUE就是非零的整數(shù)即可,當(dāng)然也可以直接寫1,FALSE當(dāng)然就是0了;如果是1則調(diào)用該宏,0的話則不調(diào)用;

5. 可以在并行版中應(yīng)用,不過得注意將UDF改寫為并行版本的。

以上資料來源:http://www.efluid.com.cn/bbs/showtopic-4306.aspx

范例:

#include "udf.h"
#include "mem.h"
#include "dynamesh_tools.h"

static real v_x = 0.0,v_y=0.0,omega_z=0.0;
DEFINE_CG_MOTION(piston,dt,vel,omega,time,dtime)
{
face_t f;
cell_t c;
real f_glob[ND_ND],m_glob[ND_ND],x_cg[ND_ND],dv_x,dv_y,domega_z;
Domain *domain= Get_Domain (1);
Thread *tf1 = Lookup_Thread (domain, 3);
int i;

NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);
if (!Data_Valid_P())
return;

for(i=0;i<=ND_ND;i++)
{
f_glob[i]=0;
m_glob[i]=0;
}

for(i=0;i<ND_ND;i++)
x_cg[i]=DT_CG(dt)[i];

if(time<=0.001)
{
vel[1]=34.5;
Message ("ntime=%f ,x=%.10lf ,y=%.10lf , force_x=%.1f ,force_y=%.1f ,moment_x=%f,moment_y=%f,moment_z=%f ,Vx=%f ,Vy=%f ,omega_z=%fn", time,x_cg[0],x_cg[1],f_glob[0],f_glob[1],m_glob[0],m_glob[1],m_glob[2],v_x,v_y,omega_z);
return;
}


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

dv_x = dtime * f_glob[0] / 39400;
v_x += dv_x;

dv_y = dtime * f_glob[1] / 39400;
v_y += dv_y;

domega_z=dtime*m_glob[2]/390000;
omega_z+=domega_z;

Message ("ntime=%f ,x=%.10lf ,y=%.10lf , force_x=%.1f ,force_y=%.1f ,moment_x=%f,moment_y=%f,moment_z=%f ,Vx=%f ,Vy=%f ,omega_z=%fn", time,x_cg[0],x_cg[1],f_glob[0],f_glob[1],m_glob[0],m_glob[1],m_glob[2],v_x,v_y,omega_z);


vel[0] = v_x,vel[1]=v_y,omega[2]=omega_z;
}

范例資料來源:http://blog.sina.com.cn/s/blog_58216001010007ur.html

上述范例是根據(jù)力和力矩求解計(jì)算物體的運(yùn)動(dòng),來源于網(wǎng)絡(luò),未經(jīng)過本人測(cè)試。


開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:函數(shù)Compute_Force_And_Moment Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報(bào)名:
  • 客服在線請(qǐng)直接聯(lián)系我們的客服,您也可以通過下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號(hào)金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com