版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

血未冷

鐵桿木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 7512.9
散金: 76
紅花: 1
帖子: 326
在線: 38.7小時
蟲號: 2637318
注冊: 2013-09-05
性別: GG
專業(yè): 工程熱物理與能源利用

[求助] 方程求解出錯，報錯“未定義函數(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公式計算納米流體動力黏度
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; %計算納米流體密度
Cpnf=(fi*densnp*Cpnp+(1-fi)*densf*Cpf)/densnf;%計算納米流體的比熱
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)用宣益民模型計算導(dǎo)熱系數(shù)
Prnf=unf*Cpnf/conductnf;             %計算納米流體的普朗特數(shù)
Rel=(Gwf*(1-xdry)*di/unf);             %計算液相Re
Frl=Gwf^2/(densf^2*9.81*di);             %計算Fr數(shù)
Xtt=((1-xdry)/xdry)^0.9*(densg/densnf)^0.5*(unf/ufg)^0.1;             %計算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);       %%計算核態(tài)沸騰抑制因子，無量綱
hl=0.023*Rel^0.8*Prnf^0.4*conductnf/di;                   %%計算液相對流換熱系數(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公式計算核態(tài)沸騰換熱系數(shù)，W/(m2.K)（COOPER公式單位必須用W/(m2.K)
hi0=Frl^(0.1-2*Frl)*F*hl+Frl^0.5*S*hnb;  %%計算對流沸騰換熱系數(shù)hi,W/(m2.K)
hi1=hl*F+hnb*S;    %%計算對流沸騰換熱系數(shù)hi,W/(m2.K)
if Frl<=0.05
hi=hi0;
else
hi=hi1;
  end


error=hi-htpx;
end

回復(fù)此樓

» 猜你喜歡

一志愿武理材料305分求調(diào)劑已經(jīng)有5人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有7人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有5人回復(fù)
354求調(diào)劑已經(jīng)有5人回復(fù)
材料專業(yè)求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有7人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有3人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

MATLAB 基礎(chǔ)問題求教已經(jīng)有6人回復(fù)

專心志，憂天下！

1樓 2015-05-22 10:21:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiegangmai

版主 (職業(yè)作家)

我沒頭銜

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
散金: 2394
紅花: 63
沙發(fā): 4
帖子: 4046
在線: 1330.2小時
蟲號: 653607
注冊: 2008-11-13
性別: GG
專業(yè): 機(jī)械動力學(xué)
管轄: 仿真模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
血未冷: 金幣+30, ★有幫助, 那個地方打錯了。。。 2015-05-22 11:25:52

calerror函數(shù)第13行：F=1+24000*(hi*deltaT/Gwf*ifg)^1.16+1.37*(1/Xtt)^0.86;
中hi是沒有定義。

贊一下(1人)

回復(fù)此樓

明德厚學(xué)、求是創(chuàng)新

2樓2015-05-22 11:14:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

血未冷

鐵桿木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 7512.9
散金: 76
紅花: 1
帖子: 326
在線: 38.7小時
蟲號: 2637318
注冊: 2013-09-05
性別: GG
專業(yè): 工程熱物理與能源利用

引用回帖:

2樓: Originally posted by xiegangmai at 2015-05-22 11:14:04
calerror函數(shù)第13行：F=1+24000*(hi*deltaT/Gwf*ifg)^1.16+1.37*(1/Xtt)^0.86;
中hi是沒有定義。

那個地方打錯了。。。不過還是沒能運(yùn)行出來。。初值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)了“錯誤使用 fzero (line 274)區(qū)間端點(diǎn)處的函數(shù)值必須具有不同的符號�！鼻蠼�!!

贊一下

回復(fù)此樓

專心志，憂天下！

3樓2015-05-22 11:25:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

血未冷

鐵桿木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 7512.9
散金: 76
紅花: 1
帖子: 326
在線: 38.7小時
蟲號: 2637318
注冊: 2013-09-05
性別: GG
專業(yè): 工程熱物理與能源利用

引用回帖:

2樓: Originally posted by xiegangmai at 2015-05-22 11:14:04
calerror函數(shù)第13行：F=1+24000*(hi*deltaT/Gwf*ifg)^1.16+1.37*(1/Xtt)^0.86;
中hi是沒有定義。

F=1+24000*(htpx*deltaT/Gwf*Hfg)^1.16+1.37*(1/Xtt)^0.86;這行應(yīng)該是這個

贊一下

回復(fù)此樓

專心志，憂天下！

4樓2015-05-22 11:28:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

血未冷

鐵桿木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 7512.9
散金: 76
紅花: 1
帖子: 326
在線: 38.7小時
蟲號: 2637318
注冊: 2013-09-05
性別: GG
專業(yè): 工程熱物理與能源利用

引用回帖:

2樓: Originally posted by xiegangmai at 2015-05-22 11:14:04
calerror函數(shù)第13行：F=1+24000*(hi*deltaT/Gwf*ifg)^1.16+1.37*(1/Xtt)^0.86;
中hi是沒有定義。

F=1+24000*(htpx*deltaT/Gwf*ifg)^1.16+1.37*(1/Xtt)^0.86;只改hi就行。。別的不用動。。剛才發(fā)錯了

贊一下

回復(fù)此樓

專心志，憂天下！

5樓2015-05-22 11:30:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiegangmai

版主 (職業(yè)作家)

我沒頭銜

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
散金: 2394
紅花: 63
沙發(fā): 4
帖子: 4046
在線: 1330.2小時
蟲號: 653607
注冊: 2008-11-13
性別: GG
專業(yè): 機(jī)械動力學(xué)
管轄: 仿真模擬

引用回帖:

3樓: Originally posted by 血未冷 at 2015-05-22 11:25:33
那個地方打錯了。。。不過還是沒能運(yùn)行出來。。初值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),2 ...

不了解你模型的專業(yè)背景，不知道你要計算什么。
錯誤提示已經(jīng)很明顯。

贊一下(1人)

回復(fù)此樓

明德厚學(xué)、求是創(chuàng)新

6樓2015-05-22 11:45:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

信彼南山

引用回帖:

6樓: Originally posted by xiegangmai at 2015-05-22 11:45:45
不了解你模型的專業(yè)背景，不知道你要計算什么。
錯誤提示已經(jīng)很明顯。...

7樓2015-05-22 21:40:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主血未冷的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿武理材料305分求調(diào)劑 +4	想上岸的鯉魚 2026-03-18	5/250	2026-03-18 17:21 by 216的故事
[考研] 344求調(diào)劑 +5	knight344 2026-03-16	6/300	2026-03-18 15:17 by knight344
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 281求調(diào)劑（0805） +5	煙汐憶海 2026-03-16	13/650	2026-03-18 14:30 by stone_128
[考研] 085601材料工程專碩求調(diào)劑 +6	慕寒mio 2026-03-16	6/300	2026-03-18 14:26 by 007_lilei
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +9	yangfz 2026-03-17	9/450	2026-03-18 12:38 by 盡舜堯1
[考研] 0703化學(xué)調(diào)劑 +4	pupcoco 2026-03-17	7/350	2026-03-18 12:14 by djl2006
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實驗室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321

24小時熱門版塊排行榜

[求助] 方程求解出錯，報錯“未定義函數(shù)或變量 'hi'”求助！ 已有1人參與

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

[求助] 方程求解出錯，報錯“未定義函數(shù)或變量 'hi'”求助！已有1人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: