| 查看: 3362 | 回复: 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); } |
回复此楼» 猜你喜欢
085600，材料与化工321分求调剂
已经有5人回复
-
求化学调剂
已经有3人回复
-
调剂求院校招收
已经有3人回复
-
材料与化工（0856）304求B区调剂
已经有8人回复
-
320分，材料与化工专业，求调剂
已经有11人回复
-
0856求调剂
已经有10人回复
-
266分，求材料冶金能源化工等调剂
已经有9人回复
-
085600 286分 材料求调剂
已经有8人回复
-
085701环境工程，267求调剂
已经有16人回复
-
286求调剂
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
9楼2016-07-05 12:58:11
2楼2016-04-22 19:12:21
xwx1990
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1155.3
- 散金: 313
- 红花: 7
- 帖子: 558
- 在线: 110.2小时
- 虫号: 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