| 查看: 2117 | 回复: 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为什么不能实现 |
回复此楼» 猜你喜欢
时间戳又变了8-15
已经有23人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有61人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
