| 查看: 2127 | 回复: 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为什么不能实现 |
回复此楼» 猜你喜欢
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有9人回复
-
闲聊
已经有3人回复
-
江苏省自然基金 什么时候出结果
已经有4人回复
-
如何查询2026年获批项目名称？
已经有3人回复
-
大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。
已经有11人回复
-
两块石头
已经有10人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
初秋的晨风
已经有6人回复