| 5 | 1/1 | 返回列表 |
| 查看: 3365 | 回復(fù): 10 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
yangjg1991鐵蟲 (初入文壇)
|
[求助]
讀入U(xiǎn)DF,一初始化FLUENT就閃退 已有2人參與
|
|
|
我現(xiàn)在做的是FLUENT動網(wǎng)格模擬,想讓流體材料的粘度隨著外部變量I不斷變化,而I是隨著時(shí)間變化的,就讀入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ā)表 | |
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +6 | 溫柔的晚安 2026-03-25 | 7/350 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +6 | 芋泥小鈴鐺 2026-03-28 | 6/300 |
|
|
[考研] 0703化學(xué) +10 | 妮妮ninicgb 2026-03-27 | 10/500 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|