| 查看: 3699 | 回复: 10 | |||
[求助]
读入UDF，一初始化FLUENT就闪退 已有2人参与
|
|||
|
我现在做的是FLUENT动网格模拟，想让流体材料的粘度随着外部变量I不断变化，而I是随着时间变化的，就读入I文件，编好了UDF，在FLUENT里面能interpreted，但只要一初始化initialize，FLUENT就闪退，不知道是UDF程序有问题，还是操作步骤有问题，求大神帮助。 UDF程序如下： /* UDF for Herschel-Bulkley viscosity */ #include "udf.h" #include "math.h" #include"stdio.h" DEFINE_PROPERTY(cell_viscosity,f,t) { real vis, s_mag, s_mag_c, sigma_y,n,k,b; int i; double a[1001]; FILE* fp=fopen("I.txt","r" ) ; for(i=0;i<1001;i++) { fscanf(fp,"%lf",&b); a=b; }; i=0; begin_f_loop(f, t) { real time=CURRENT_TIME; int nDiv=(int)(time*100); { if(time==0) { s_mag = C_STRAIN_RATE_MAG(f,t); k= -1.574*pow(a,3)+7.102*pow(a,2)-10.24*a+4.728; n= 0.1573*pow(a,2)-0.607*a+0.6781; sigma_y = -2.145e4*pow(a,2)+7.505*a ; s_mag_c = 376.1*pow(a,3)-1332*pow(a,2)+1561*a+196.8 ; if (s_mag < s_mag_c) { vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c); } else { vis = sigma_y / s_mag + k*pow(s_mag/s_mag_c, (n-1)); } } else if(time!=0&&nDiv%2==0) { i++; s_mag = C_STRAIN_RATE_MAG(f,t); k= -1.574*pow(a,3)+7.102*pow(a,2)-10.24*a+4.728; n= 0.1573*pow(a,2)-0.607*a+0.6781; sigma_y = -2.145e4*pow(a,2)+7.505*a ; s_mag_c = 376.1*pow(a,3)-1332*pow(a,2)+1561*a+196.8 ; if (s_mag < s_mag_c) { vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c); } else { vis = sigma_y / s_mag + k*pow(s_mag/s_mag_c, (n-1)); } } else { s_mag = C_STRAIN_RATE_MAG(f,t); k= -1.574*pow(a,3)+7.102*pow(a,2)-10.24*a+4.728; n= 0.1573*pow(a,2)-0.607*a+0.6781; sigma_y = -2.145e4*pow(a,2)+7.505*a ; s_mag_c = 376.1*pow(a,3)-1332*pow(a,2)+1561*a+196.8 ; if (s_mag < s_mag_c) { vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c); } else { vis = sigma_y / s_mag + k*pow(s_mag/s_mag_c, (n-1)); } } } return vis; } end_f_loop(f,t); } |
回复此楼» 猜你喜欢
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
5楼2016-04-27 16:46:00
7楼2016-04-27 16:50:12
2楼2016-04-22 19:12:21
xwx1990
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1135.3
- 散金: 313
- 红花: 7
- 帖子: 559
- 在线: 110.4小时
- 虫号: 3092365
- 注册: 2014-03-27
- 性别: GG
- 专业: 流体力学
3楼2016-04-22 20:13:22
moons2
金虫 (正式写手)
- 应助: 144 (高中生)
- 金币: 1440.5
- 红花: 8
- 帖子: 431
- 在线: 239.4小时
- 虫号: 3784607
- 注册: 2015-04-03
- 性别: GG
- 专业: 流体力学
4楼2016-04-22 20:38:20
6楼2016-04-27 16:47:45
8楼2016-05-04 09:45:17
9楼2016-07-05 12:58:11
10楼2019-05-23 21:10:17
100