| 查看: 3332 | 回復: 10 | ||
yangjg1991鐵蟲 (初入文壇)
|
[求助]
讀入UDF,一初始化FLUENT就閃退 已有2人參與
|
|
|
我現(xiàn)在做的是FLUENT動網格模擬,想讓流體材料的粘度隨著外部變量I不斷變化,而I是隨著時間變化的,就讀入I文件,編好了UDF,在FLUENT里面能interpreted,但只要一初始化initialize,F(xiàn)LUENT就閃退,不知道是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); } |
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
金蟲 (正式寫手)

金蟲 (正式寫手)
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調劑 +4 | 梁梁校校 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 296求調劑 +3 | www_q 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|
[考研] 0703化學調劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 085600材料與化工調劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 085600材料與化工求調劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 304求調劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 330求調劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 354求調劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
[考研] 328求調劑,英語六級551,有科研經歷 +3 | 生物工程調劑 2026-03-17 | 7/350 |
|
|
[考研] 材料專業(yè)求調劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 0703化學336分求調劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 278求調劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 334求調劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研] 考研化學學碩調劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 283求調劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 333求調劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 290求調劑 +3 | ADT 2026-03-13 | 3/150 |
|