| 2 | 1/1 | 返回列表 |
| 查看: 657 | 回復: 1 | ||
bingshan2012木蟲 (正式寫手)
|
[求助]
大家?guī)兔聪逻@段粘度系數(shù)的UDF有沒有問題,謝了。
|
這是UDF中定義粘度系數(shù)的一段程序,大家?guī)兔纯从袥]有語法問題,謝謝,x數(shù)組為壓力值(單位pa),y數(shù)組為溫度值(單位T),z數(shù)組為不同壓力溫度下所對應的粘度系數(shù)值,為實驗數(shù)據(jù),因為只有1atm~16atm,300K~40000K的數(shù)據(jù),所以想在壓力小于1atm時按1atm處理,大于16atm按16atm處理,溫度小于300K按300K處理,大于40000K按40000K處理,slq3()是前面定義的一個插值函數(shù),主要是這幾個if,這樣寫沒有問題吧,幫忙看看, 。DEFINE_PROPERTY(SF6viscosity,cell,thread) { double PRE,TEM,VIS; double slq3(x,y,z,n,m,u,v); double x[5]={101325,202650,405300,810600,1621200}; double y[30]={300,500,1000,1400,1700,2000,2300,2600,3000,3500,4000,5000,6000,7000,8000,9000,10000,11000,12000,13000,14000,15000,16000,18000,20000,22000,26000,30000,35000,40000}; double z[5][30]={{1.22E-5,1.97E-5,3.62E-5,4.57E-5,5.63E-5,6.92E-5,7.94E-5,9.46E-5,1.10E-4,1.23E-4,1.36E-4,1.74E-4,2.09E-4,2.56E-4,2.94E-4,3.32E-4,3.60E-4,3.83E-4,4.00E-4,3.95E-4,3.42E-4,2.54E-4,1.69E-4,6.78E-5,3.24E-5,2.17E-5,1.78E-5,1.62E-5,0.88E-5,0.60E-5}, {1.22E-5,1.97E-5,3.62E-5,4.57E-5,5.55E-5,6.90E-5,7.77E-5,9.18E-5,1.09E-4,1.22E-4,1.36E-4,1.74E-4,2.09E-4,2.56E-4,2.95E-4,3.33E-4,3.60E-4,3.86E-4,4.07E-4,4.19E-4,3.95E-4,3.25E-4,2.39E-4,1.08E-4,5.03E-5,3.05E-5,2.20E-5,2.10E-5,1.32E-5,0.81E-5}, {1.22E-5,1.97E-5,3.62E-5,4.57E-5,5.37E-5,6.79E-5,7.66E-5,8.87E-5,1.07E-4,1.22E-4,1.36E-4,1.74E-4,2.09E-4,2.56E-4,2.94E-4,3.34E-4,3.64E-4,3.89E-4,4.12E-4,4.33E-4,4.32E-4,3.89E-4,3.15E-4,1.65E-4,7.98E-5,4.56E-5,2.80E-5,2.62E-5,1.96E-5,1.15E-5}, {1.22E-5,1.97E-5,3.62E-5,4.56E-5,5.30E-5,6.61E-5,7.60E-5,8.60E-5,1.05E-4,1.22E-4,1.36E-4,1.74E-4,2.09E-4,2.56E-4,2.95E-4,3.34E-4,3.65E-4,3.90E-4,4.15E-4,4.41E-4,4.55E-4,4.37E-4,3.86E-4,2.38E-4,1.24E-4,6.96E-5,3.78E-5,3.27E-5,2.76E-5,1.73E-5}, {1.22E-5,1.97E-5,3.62E-5,4.56E-5,5.25E-5,6.38E-5,7.54E-5,8.41E-5,1.02E-4,1.21E-4,1.35E-4,1.74E-4,2.09E-4,2.56E-4,2.95E-4,3.34E-4,3.66E-4,3.92E-4,4.18E-4,4.46E-4,4.69E-4,4.70E-4,4.44E-4,3.23E-4,1.88E-4,1.08E-4,5.23E-5,4.19E-5,3.73E-5,2.59E-5}}; PRE=C_P(cell,thread); if (PRE<101325) PRE=101325; else if (PRE>1621200) PRE=1621200; TEM=C_T(cell,thread); if (TEM<300) TEM=300; else if (TEM>40000) TEM=40000; VIS=slq3(x,y,z,5,30,PRE,TEM); return VIS; } |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 083200學碩321分一志愿暨南大學求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 南昌大學材料專碩311分求調(diào)劑 +5 | 77chaselx 2026-03-20 | 5/250 |
|
|
[考研] 考研調(diào)劑求學校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 一志愿蘇州大學材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 收復試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|