- 积分
- 0
- 注册时间
- 2018-11-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 寻找tamia 于 2019-1-14 15:27 编辑
#include "udf.h"
DEFINE_SDOF_PROPERTIES(liheliju,prop,dt,time,dtime)
{
real x_cg[ND_ND];/*储存运动物体重心的位置变化量*/
real f_glob[ND_ND];/*储存运动物体在流场中的流体力*/
real m_glob[ND_ND];/*储存流体力对重心位置的力矩*/
int i;
FILE *fp;
real ntime;
Thread *t;/*定义边界条件*/
Domain *domain=Get_Domain(1);
t=Lookup_Thread(domain,14);
ntime=RP_Get_Real("flow-time");/*这里调用了获取流动模拟的物理时间函数*/
prop[SDOF_MASS]=81.67994;
prop[SDOF_IXX]=0.80599;
prop[SDOF_IYY]=31.97144;
prop[SDOF_IZZ]=32.27510;
if (ntime>0 && ntime<=1)
{
prop[SDOF_ZERO_ROT_Y]=TRUE;/*限制沿着Y的转动*/
prop[SDOF_ZERO_TRANS_Z]=TRUE;/*限制沿着Z的平动*/
}
if (ntime>1)
x_cg[0]=DT_CG(dt)[0];
x_cg[1]=DT_CG(dt)[1];
x_cg[2]=DT_CG(dt)[2];
Compute_Force_And_Moment(domain,t,x_cg,f_glob,m_glob,TRUE);
Message ("%e %e %e %e %e %e %e %e %e %e", ntime,x_cg[0],x_cg[1],x_cg[2],f_glob[0],f_glob[1],f_glob[2],m_glob[0],m_glob[1],m_glob[2]);
fp=fopen("shuju.txt","a");
fprintf(fp,"%e %e %e %e %e %e %e %e %e %e\n", ntime,x_cg[0],x_cg[1],x_cg[2],f_glob[0],f_glob[1],f_glob[2],m_glob[0],m_glob[1],m_glob[2]);
fclose(fp);
}
x_cg[0]=DT_CG(dt)[0];
x_cg[1]=DT_CG(dt)[1];
x_cg[2]=DT_CG(dt)[2];
报错 structure not implement
大神求助
|
|