| 查看: 2035 | 回复: 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为什么不能实现 |
回复此楼» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
