| 查看: 2067 | 回复: 0 | ||
[求助]
udf读取数据
|
|
#include "udf.h" #include "math.h" real U_inf_variable[100]; real x[100]; real time; int vel; real cg_omega[2]; #define Nbpts 100 int i; int counter=0; FILE *fid; FILE *fid2; FILE *fid3; DEFINE_EXECUTE_ON_LOADING(readfile,libname) { fid=fopen("shuju.txt","r"  ;for (i=0;i<Nbpts;i++) { fscanf(fid,"%le ",&U_inf_variable); } fclose(fid); fid2=fopen("check_load_shuju.txt","w" ;for (i=0;i<Nbpts;i++) { fprintf(fid2,"%le %d\n",&U_inf_variable,counter); } fclose(fid2); } DEFINE_CG_MOTION(rotation,dt,cg_vel,cg_omega,time,dtime) { if(time<3) cg_omega[2]=-0.024532; else if(time>3&&time<=10) cg_omega[2]=-0.01214; } DEFINE_PROFILE(power_3,thread,position) { face_t f; time=CURRENT_TIME; begin_f_loop(f,thread) { F_PROFILE(f,thread,position)=U_inf_variable[counter]; } end_f_loop(f,thread) } DEFINE_EXECUTE_AT_END(execute) { face_t f; Domain*domain=Get_Domain(1); Thread *tf1=Lookup_Thread(domain,126); begin_f_loop(f,tf1) { F_CENTROID(x,f,tf1 ); vel= F_V(f,tf1); fid3=fopen("vel.dat", "a+" ;fprintf(fid3," %.6f\n",vel); } end_f_loop(f,thread) fclose(fid3); } 我想实现进口速度边界条件导入，经过一个时间步的迭代计算后，再从这个界面导出速度。。。请问大神我这个UDF为什么不能实现 |
回复此楼
