| 查看: 2075 | 回复: 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为什么不能实现 |
回复此楼» 猜你喜欢
求调剂
已经有3人回复
-
一志愿北京化工大学，初试成绩350求调剂
已经有14人回复
-
复试调剂
已经有4人回复
-
323分（计算机视觉和大模型项目）能直接上手
已经有3人回复
-
311分 22408 求调剂
已经有3人回复
-
320分人工智能调剂
已经有7人回复
-
一志愿郑大0705求调剂
已经有4人回复
-
0703化学
已经有10人回复
-
301求调剂
已经有10人回复
-
306分材料与化工求调剂
已经有6人回复
