| 查看: 3330 | 回復(fù): 10 | ||
yangjg1991鐵蟲 (初入文壇)
|
[求助]
讀入U(xiǎn)DF,一初始化FLUENT就閃退 已有2人參與
|
|
|
我現(xiàn)在做的是FLUENT動(dòng)網(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ā)表 | |
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|---|---|---|---|---|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 288求調(diào)劑 +11 | 于海海海海 2026-03-19 | 11/550 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|