|
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
想请教高手为什么下面第一个UDF既可以编译也可以导入,但为什么第二个却只能编译,导入时提示系统找不到文件?
1)
#include "udf.h"
#include "dynamesh_tools.h"
#include "f_wall.h"
#define ZONE_ID1 6
#define ZONE_ID2 5
#define ZONE_ID3 4
double MASS=6.716635;
double IZZ=1.684491;
#define IXZ 0.0
#define IXY 0.0
#define IYZ 0.0
double omega_globle=0.1; /*公转速度 rad/s*/
double omega_local=1; /*自转速度 rad/s*/
double RR=5.0;
double PI=3.1415926535898;
static double cg1_omega_previous=0.0;
static double alfa1=-30.0/180*3.1415926535898;
static double alfa2=-150.0/180*3.1415926535898;
static double alfa3=90.0/180*3.1415926535898;
DEFINE_CG_MOTION(foil_1, dt, cg_vel, cg_omega, time, dtime)
{
alfa1=alfa1+omega_globle*dtime;
cg_vel[0]=-RR*omega_globle*sin(alfa1);
cg_vel[1]=RR*omega_globle*cos(alfa1);
cg_omega[2]=omega_globle+omega_local;
/*Message("alfa1=%f,x_velocity=%f,y_velocity=%f \n",alfa1,cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_1_out, dt, cg_vel, cg_omega, time, dtime)
{
cg_vel[0]=-RR*omega_globle*sin(alfa1);
cg_vel[1]=RR*omega_globle*cos(alfa1);
cg_omega[2]=omega_globle+omega_local;
/*Message("x_velocity=%f,y_velocity=%f \n",cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_2, dt, cg_vel, cg_omega, time, dtime)
{
alfa2=alfa2+omega_globle*dtime;
cg_vel[0]=-RR*omega_globle*sin(alfa2);
cg_vel[1]=RR*omega_globle*cos(alfa2);
cg_omega[2]=omega_globle+omega_local;
/*Message("alfa2=%f,x_velocity=%f,y_velocity=%f \n",alfa2,cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_2_out, dt, cg_vel, cg_omega, time, dtime)
{
cg_vel[0]=-RR*omega_globle*sin(alfa2);
cg_vel[1]=RR*omega_globle*cos(alfa2);
cg_omega[2]=omega_globle+omega_local;
/*Message("x_velocity=%f,y_velocity=%f \n",cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_3, dt, cg_vel, cg_omega, time, dtime)
{
alfa3=alfa3+omega_globle*dtime;
cg_vel[0]=-RR*omega_globle*sin(alfa3);
cg_vel[1]=RR*omega_globle*cos(alfa3);
cg_omega[2]=omega_globle+omega_local;
/*Message("alfa3=%f,x_velocity=%f,y_velocity=%f \n",alfa3,cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_3_out, dt, cg_vel, cg_omega, time, dtime)
{
cg_vel[0]=-RR*omega_globle*sin(alfa3);
cg_vel[1]=RR*omega_globle*cos(alfa3);
cg_omega[2]=omega_globle+omega_local;
/*Message("x_velocity=%f,y_velocity=%f \n",cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(cent_fluid, dt, cg_vel, cg_omega, time, dtime)
{
cg_omega[2] = omega_globle;
}
2)
#include "udf.h"
#include "dynamesh_tools.h"
#include "f_wall.h"
#define ZONE_ID1 6
#define ZONE_ID2 5
#define ZONE_ID3 4
double MASS=6.716635;
double IZZ=1.684491;
#define IXZ 0.0
#define IXY 0.0
#define IYZ 0.0
double omega_globle=0.1; /*公转速度 rad/s*/
double e=2;
double RR=10.0;
double PI=3.1415926535898;
static double cg1_omega_previous=0.0;
static double alfa1=90.0/180*3.1415926535898;
static double alfa2=-30.0/180*3.1415926535898;
static double alfa3=-150.0/180*3.1415926535898;
DEFINE_CG_MOTION(foil_1, dt, cg_vel, cg_omega, time, dtime)
{
alfa1=alfa1+omega_globle*dtime;
cg_vel[0]=-RR*omega_globle*sin(alfa1);
cg_vel[1]=RR*omega_globle*cos(alfa1);
cg_omega[2]=omega_globle-((e*omega_globle*(sin(alft1)+e))/(1+e*e+2*e*sin(alft1)));
/*Message("alfa1=%f,x_velocity=%f,y_velocity=%f \n",alfa1,cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_1_out, dt, cg_vel, cg_omega, time, dtime)
{
cg_vel[0]=-RR*omega_globle*sin(alfa1);
cg_vel[1]=RR*omega_globle*cos(alfa1);
cg_omega[2]=omega_globle-((e*omega_globle*(sin(alft1)+e))/(1+e*e+2*e*sin(alft1)));
/*Message("x_velocity=%f,y_velocity=%f \n",cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_2, dt, cg_vel, cg_omega, time, dtime)
{
alfa2=alfa2+omega_globle*dtime;
cg_vel[0]=-RR*omega_globle*sin(alfa2);
cg_vel[1]=RR*omega_globle*cos(alfa2);
cg_omega[2]=omega_globle-((e*omega_globle*(sin(alft2)+e))/(1+e*e+2*e*sin(alft2)));
/*Message("alfa2=%f,x_velocity=%f,y_velocity=%f \n",alfa2,cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_2_out, dt, cg_vel, cg_omega, time, dtime)
{
cg_vel[0]=-RR*omega_globle*sin(alfa2);
cg_vel[1]=RR*omega_globle*cos(alfa2);
cg_omega[2]=omega_globle-((e*omega_globle*(sin(alft2)+e))/(1+e*e+2*e*sin(alft2)));
/*Message("x_velocity=%f,y_velocity=%f \n",cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_3, dt, cg_vel, cg_omega, time, dtime)
{
alfa3=alfa3+omega_globle*dtime;
cg_vel[0]=-RR*omega_globle*sin(alfa3);
cg_vel[1]=RR*omega_globle*cos(alfa3);
cg_omega[2]=omega_globle-((e*omega_globle*(sin(alft3)+e))/(1+e*e+2*e*sin(alft3)));
/*Message("alfa3=%f,x_velocity=%f,y_velocity=%f \n",alfa3,cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(foil_3_out, dt, cg_vel, cg_omega, time, dtime)
{
cg_vel[0]=-RR*omega_globle*sin(alfa3);
cg_vel[1]=RR*omega_globle*cos(alfa3);
cg_omega[2]=omega_globle-((e*omega_globle*(sin(alft3)+e))/(1+e*e+2*e*sin(alft3)));
/*Message("x_velocity=%f,y_velocity=%f \n",cg_vel[0],cg_vel[1]);*/
}
DEFINE_CG_MOTION(cent_fluid, dt, cg_vel, cg_omega, time, dtime)
{
cg_omega[2] = omega_globle;
} |
|