| 查看: 1194 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助帮忙改正一个FLUENT的UDF！急急急 已有1人参与
|
||
|
写了一个UDF，描述阀芯在弹簧力和流体压力作用下的运动，可是执行两次后压力总是变成很大的负值，进而速度、位移都变得很大，网格更新就不成功了。能否有人帮忙看一下？ 附UDF源码： #include <stdio.h> #include "udf.h" #define UDF_FILENAME "udf_loc_velo.txt" # define K_SPRING 31400 /* read current location and velocity from file */ DEFINE_CG_MOTION(spool, dt, vel, omega, time, dtime) { Thread *t = DT_THREAD (dt); face_t f; static real force; static real loc; static real velo; static real init_disp; static real s_force ; static real dv; static real total_force; static FILE *fp; /* Open data files */ fp=fopen(UDF_FILENAME,"r"  ;if(fp !=NULL) { float read_loc,read_velo; fscanf (fp,"%e %e",&read_loc,&read_velo); fclose(fp); loc=(real) read_loc; velo=(real) read_velo; } else { loc=0.0; velo=0.0; } /* compute force on piston wall */ force = 0.0; begin_f_loop (f, t) { real *AA; AA = F_AREA_CACHE (f, t); force += F_P (f, t) * AA[2]; } end_f_loop (f, t) /* add in spring force */ init_disp = 0.01; s_force = K_SPRING * (loc + init_disp); total_force = force - s_force-346-10*velo; /* compute change in velocity */ dv = dtime * total_force / 0.554; velo += dv; loc += velo * dtime; vel[2] = velo; write_loc_velo_file(loc,velo); Message ("\nUDF spool: time = %f, z_vel = %f, force = %f,total_f = %f, loc(m)= %f\n", time, velo, force,total_force, loc); } /* write current location and velocity in file*/ void write_loc_velo_file(real loc,real velo) { FILE*fp=fopen(UDF_FILENAME,"w" ;if(fp!=NULL) { fprintf(fp,"%e %e",loc,velo); fclose(fp); } else Message("\n Warning:cannot write %s file",UDF_FILENAME); } |
回复此楼» 猜你喜欢
时间戳又变了8-15
已经有23人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有61人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教udf无法初始化问题，各位大神帮帮忙
已经有11人回复
- 请大家帮忙指导一下这个计算题！
已经有3人回复
- Fluent中加载UDF文件出问题，求帮助！！！！！！！！！！
已经有6人回复
- 这轨道怎么看呀，各位大神帮帮忙，急急急急急.......
已经有7人回复
- 求大神帮检查下一个简单的UDF
已经有13人回复
- 动网格UDF 求帮忙
已经有5人回复
- 求助 EC50 怎么算！！！！ 求大神帮忙
已经有3人回复
- 求助大神！！1关于fluent模拟污染物散发并排除的问题
已经有4人回复
- 出口压力的UDF编写问题
已经有10人回复
- 全部身家求大神帮忙解决fluent压力脉冲怎么编写UDF的问题
已经有8人回复
- 500金币，帮忙找错，UDF并行
已经有15人回复
- 哪位大虾能帮编一个速度随时间变化的udf啊？
已经有6人回复
- 利用DEFINE_EXCHANGE_PROPERTY编的蒸发和冷凝的udf，大侠们帮我看看对不对
已经有9人回复
- 好心人帮帮忙吧，为什么这个UDF导入fluent以后，在边界条件处没有显示~~~
已经有17人回复
- udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看
已经有4人回复
- 【求助】求UDF高手帮忙看看我的问题
已经有18人回复
3楼2015-02-15 14:44:07
2楼2015-02-11 11:07:59
4楼2018-07-31 20:55:12
5楼2018-07-31 21:51:23
