| 5 | 1/1 | 返回列表 |
| 查看: 2996 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
血未冷鐵桿木蟲(chóng) (正式寫(xiě)手)
|
[求助]
方程求解出錯(cuò),報(bào)錯(cuò)“未定義函數(shù)或變量 'hi'”求助! 已有1人參與
|
|
|
function htp=calhtpbyGW1(Pcr,Pevap,Molmass,fi,Gwf,uf,ufg,xdry,di,T,deltaT,densf,densnp,densg,conductf,conductnp,Cpnp,Cpf,rp,ifg) htpxmin=1.0e+1; htpxmax=3.0e+4; opt=optimset( 'fzero'); opt=optimset( opt,'TolX', 1.0e-3); htp=fzero(@(htpx)calerror(htpx,Pcr,Pevap,Molmass,fi,Gwf,uf,ufg,xdry,di,T,deltaT,densf,densnp,densg,conductf,conductnp,Cpnp,Cpf,rp,ifg),[htpxmin htpxmax],opt); end function error=calerror(htpx,Pcr,Pevap,Molmass,fi,Gwf,uf,ufg,xdry,di,T,deltaT,densf,densnp,densg,conductf,conductnp,Cpnp,Cpf,rp,ifg) unf=(1+2.5*fi+6.25*fi)*uf; %運(yùn)用Batchelor公式計(jì)算納米流體動(dòng)力黏度 fiw=1/((1+(1-fi))*(densf/densnp));% 體積分?jǐn)?shù)換算成質(zhì)量分?jǐn)?shù) kB=1.3806505*10^-23; %玻爾茲曼常數(shù),單位J/K densnf=fiw*densnp+(1-fiw)*densf; %計(jì)算納米流體密度 Cpnf=(fi*densnp*Cpnp+(1-fi)*densf*Cpf)/densnf;%計(jì)算納米流體的比熱 conductnf=(conductnp+2*conductf-2*fi*(conductf-conductnp))/(conductnp+2*conductf+fi*(conductf-conductnp))+((fi*(densnp*Cpnp)/2*conductf))*(kB*T/3*3.14*rp*unf)^0.5;%運(yùn)用宣益民模型計(jì)算導(dǎo)熱系數(shù) Prnf=unf*Cpnf/conductnf; %計(jì)算納米流體的普朗特?cái)?shù) Rel=(Gwf*(1-xdry)*di/unf); %計(jì)算液相Re Frl=Gwf^2/(densf^2*9.81*di); %計(jì)算Fr數(shù) Xtt=((1-xdry)/xdry)^0.9*(densg/densnf)^0.5*(unf/ufg)^0.1; %計(jì)算Xtt數(shù) F=1+24000*(hi*deltaT/Gwf*ifg)^1.16+1.37*(1/Xtt)^0.86; S=1/(1+1.15*10^(-6)*F^2*Rel^1.17); %%計(jì)算核態(tài)沸騰抑制因子,無(wú)量綱 hl=0.023*Rel^0.8*Prnf^0.4*conductnf/di; %%計(jì)算液相對(duì)流換熱系數(shù),W/(m2.K) Prevapwf=Pevap/Pcr; q=htpx*deltaT; hnb=55*Prevapwf^0.12*(-log10(Prevapwf))^(-0.55)*Molmass^(-0.5)*q^0.67; %% cooper公式計(jì)算核態(tài)沸騰換熱系數(shù),W/(m2.K)(COOPER公式單位必須用W/(m2.K) hi0=Frl^(0.1-2*Frl)*F*hl+Frl^0.5*S*hnb; %%計(jì)算對(duì)流沸騰換熱系數(shù)hi,W/(m2.K) hi1=hl*F+hnb*S; %%計(jì)算對(duì)流沸騰換熱系數(shù)hi,W/(m2.K) if Frl<=0.05 hi=hi0; else hi=hi1; end error=hi-htpx; end |

鐵桿木蟲(chóng) (正式寫(xiě)手)


鐵桿木蟲(chóng) (正式寫(xiě)手)
|
那個(gè)地方打錯(cuò)了。。。不過(guò)還是沒(méi)能運(yùn)行出來(lái)。。初值htp=calhtpbyGW1(3661800,100000,152.93,0.0001,1600,0.00040588,0.000010826,0.3,0.00545,300.61,22.54,1457.6,0.00397,6.3918,0.075734,29,800,1026,10^(-8),2000)輸入后,出現(xiàn)了“錯(cuò)誤使用 fzero (line 274)區(qū)間端點(diǎn)處的函數(shù)值必須具有不同的符號(hào)!鼻蠼!! |

鐵桿木蟲(chóng) (正式寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|