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

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

血未冷

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

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

[求助] 方程求解出錯(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

回復(fù)此樓

» 猜你喜歡

求材料調(diào)劑已經(jīng)有6人回復(fù)
085600材料與化工已經(jīng)有5人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有6人回復(fù)
286求調(diào)劑已經(jīng)有9人回復(fù)
焦慮已經(jīng)有12人回復(fù)
344求調(diào)劑已經(jīng)有6人回復(fù)
266求調(diào)劑已經(jīng)有9人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有18人回復(fù)
314求調(diào)劑已經(jīng)有8人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有5人回復(fù)

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

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

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

專心志，憂天下！

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

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

xiegangmai

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

我沒(méi)頭銜

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
血未冷: 金幣+30, ★有幫助, 那個(gè)地方打錯(cuò)了。。。 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是沒(méi)有定義。

贊一下(1人)

回復(fù)此樓

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

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

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

血未冷

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

應(yīng)助: 1 (幼兒園)
金幣: 7512.9
散金: 76
紅花: 1
帖子: 326
在線: 38.7小時(shí)
蟲(chóng)號(hào): 2637318
注冊(cè): 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是沒(méi)有定義。

那個(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)。”求解!!

贊一下

回復(fù)此樓

專心志，憂天下！

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

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

血未冷

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

應(yīng)助: 1 (幼兒園)
金幣: 7512.9
散金: 76
紅花: 1
帖子: 326
在線: 38.7小時(shí)
蟲(chóng)號(hào): 2637318
注冊(cè): 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是沒(méi)有定義。

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

贊一下

回復(fù)此樓

專心志，憂天下！

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

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

血未冷

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

應(yīng)助: 1 (幼兒園)
金幣: 7512.9
散金: 76
紅花: 1
帖子: 326
在線: 38.7小時(shí)
蟲(chóng)號(hào): 2637318
注冊(cè): 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是沒(méi)有定義。

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

贊一下

回復(fù)此樓

專心志，憂天下！

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

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

xiegangmai

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

我沒(méi)頭銜

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

引用回帖:

3樓: Originally posted by 血未冷 at 2015-05-22 11:25:33
那個(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),2 ...

不了解你模型的專業(yè)背景，不知道你要計(jì)算什么。
錯(cuò)誤提示已經(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ì)算什么。
錯(cuò)誤提示已經(jīng)很明顯。...

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

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

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

返回列表

普通表情龍兔虎貓

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

[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	9/450	2026-03-18 15:31 by vgtyfty
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 268求調(diào)劑 +6	簡(jiǎn)單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無(wú)際的草原
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院�；瘜W(xué)工程與工藝，誠(chéng)求調(diào)劑 +7	小劉快快上岸 2026-03-11	8/400	2026-03-17 16:57 by ruiyingmiao
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant

24小時(shí)熱門(mén)版塊排行榜

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

» 猜你喜歡

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

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

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

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

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