找回密码
 注册
查看: 1935|回复: 0

UDF导入出错

[复制链接]
发表于 2012-3-21 15:36:29 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

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;
       
}
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表