 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1121 | 回復(fù): 2 | ||
892226227銅蟲 (初入文壇)
|
[求助]
求助修改模擬列管換熱器的程序。感激不盡
|
|
求助：幫我修改一下這個模擬列管換熱器的程序。感激不盡 %已知 wh=68250/3600; wc=175000/3600; Th1=382; Tc1=275; A=130;%傳熱面積 Di=0.0225;%管內(nèi)徑 Do=0.025;%管外徑 D=0.7;%殼體內(nèi)徑 L=0.6;%管長 AT=0.0438;%管程通道截面積 AS=0.0525;%殼程通道截面積 NB=19;%殼程擋板數(shù) B=0.3;%擋板間距 Ri=0.0005;%管程污垢熱阻 Ro=0.0001;%殼程污垢熱阻 Ki=12.5;%渣油品特性參數(shù) Ko=12.5;%原油品特性參數(shù) D1=0.919;%渣相對密度 D0=0.850;%原相對密度 Tc20=285;%冷流體的估算出口溫度 %計算管中心距 t=Di/(1-AS/(B*D)); %計算殼體當(dāng)量直徑 De=4*(t^2-pi*Do^2/4)/(pi*Do); while 1 %最外層循環(huán)求冷熱物流出口溫度 %計算管內(nèi)熱物流性質(zhì) Thm=(Th1+Th20)/2; cpi=((0.7072+0.000551*D1)*Thm-0.318D1)*(0.055*Ki+0.35)*4.18; %計算管外冷物流性質(zhì) Tc2=fzero(@wd,295,[],Tc1,Th1,Th20,Thm,D0,cpi,Ko,wh,wc);%求只含Tc2的一元函數(shù)的零點(diǎn) %--------函數(shù)定義------- function s =wd(Tc2,Tc1,Th1,Th20,Thm,D0,cpi,Ko,wh,wc); s=(wh*cpi)*(Th1-Th20)-(wc*cpo)*(Tc2-Tc1); Tcm=(Tc1+Tc2)/2; cpo=((0.7072+0.000551*D0)*Tcm-0.318D0)*(0.055*Ko+0.35)*4.18; %計算傳熱量 Q=wh*(Th1-Th20); %計算管內(nèi)外流體的相對密度 roui=0.942+0.248*x+0.174*D1^2+0.0841/(x*D1)-0.312*x/D1-0.556*exp(-x); xi=1+Thm/100; rouo=0.942+0.248*x+0.174*D0^2+0.0841/(x*Do0)-0.312*x/D0-0.556*exp(-x); xo=1+Tcm/100; %管內(nèi)體積流速 ui=wh/(AT*roui); %管外體積流速 uo=wc/(AS*rouo); i=1; fprintf('no=%d\tuo\tui\n',i,uo.ui); %計算管內(nèi)外流體熱導(dǎo)率 lumdi=0.4213*(1-0.00054*Thm)/D1; lumdo=0.4213*(1-0.00054*Tcm)/D0; %計算管內(nèi)外流體運(yùn)動粘度 a=log(log(1500+1.22))-blog(50+273); b=log((log(1500+1.22))/log(120+1.22))/log((50+273)/(100+273)); miui=exp(exp(a+b*log(Thm+273)))-1.22; a=log(log(90+1.22))-blog(50+273); b=log((log(90+1.22))/log(13+1.22))/log((50+273)/(100+273)); miuo=exp(exp(a+b*log(Tcm+273)))-1.22; %粘度 ati=miui*roui; ato=miuo*rouo; %計算管內(nèi)外流體特征數(shù) Reo=Do*Wc/(ato*AS); Rei=Di*Wh/(ati*AT); Pro=cpo*ato/lumdo; Pri=cpi*ati/lumdi; %估算內(nèi)外壁溫 Hi1=0.023*Rei^0.8*Pri^0.33*lumdi/Di; %初值為不受運(yùn)動粘度影響的表面儲熱系數(shù) while 1 %迭代求內(nèi)壁溫 Twi=Tmc+Q/(Hi1*A*225/250); a=log(log(1500+1.22))-blog(50+273); b=log((log(1500+1.22))/log(120+1.22))/log((50+273)/(100+273)); miuiw=exp(exp(a+b*log(Twi+273)))-1.22; Hi=Hi1(miui/miuiw)^0.14; %受溫度影響的表面?zhèn)鳠嵯禂?shù) Tiw=Tmh+Q/(Hi*A*225/250); %替代 if abs(Tiw-Twi)<=0.5 %收斂 break; else Twi=Tiw; end %if語句結(jié)束 end %內(nèi)壁溫循環(huán)結(jié)束 %求管外的表面?zhèn)鳠嵯禂?shù)的初值 if Reo>2000 Ho1=0.36*Reo^0.55*Pro^0.33*lumdo/Do; else Ho1=0.527*Reo^0.5*Pro^0.33*lumdo/Do; end %if語句結(jié)束 while 1 %迭代求外壁溫 Two=Tmc+Q/(Ho1*A); a=log(log(1500+1.22))-blog(50+273); b=log((log(1500+1.22))/log(120+1.22))/log((50+273)/(100+273)); miuow=exp(exp(a+b*log(Two+273)))-1.22; Ho=Ho1(miuo/miuow)^0.14; Tow=Tmc+Q/(Ho*A); if abs(Tow-Two)<=0.5 %收斂 break; else Two=Tow;%替代 end %if語句結(jié)束 end %內(nèi)壁溫循環(huán)結(jié)束 %表面?zhèn)鳠嵯禂?shù) K=1/((1/Hi+Ri)*250/225+1/Ho+Ro); %計算傳熱效率 ibcn=2/((1+wh*cpi/wc*cpo)+sqrt(1+wh*cpi/wc*cpo)*(1+exp(-F))/(1-exp(-F))); F=ntu*sqrt(1+(wh*cpi/wc*cpo)^2); ntu=K*A/(wh*cpi); %核算Th2 Th2=Th1-ibcn(Th1-Tc1); Tc2=Tc1+(wh*cpi/wc*cpo)*(Th1-Th2); if abs(Tc20-Tc2)<=0.1 %收斂 break; else T20=Th2;%替代 end %if語句結(jié)束 fprintf('no=%d\tuii=%f\tReii=%f\tPrii=%f\tHii=%f\tKii=%f\n',i,ui,Rei,Pri,Hi,Ki); fprintf('no=%d\tuoi=%f\tReoi=%f\tProi=%f\tHoi=%f\tKoi=%f\n',i,uo,Reo,Pro,Ho,Ko); end % 最外層求冷熱物流出口溫度循環(huán)結(jié)束 fprintf('計算阻力…'); %計算換熱器阻力 %管程阻力 if Rei<100000 fi=0.4513*Rei^(-0.2663); else fi=0.2864*Rei^(-0.2258); end gi=ui*roui; dpt=(fi*L/Di+4)*gi^2*2*1.5/(2000*D1); %殼程阻力 if Reo<150 fo=120*Reo^(-0.993); elseif Reo>1500 fo=0.7664*Reo^(-0.0854); else fo=10*(15.312/(log(Reo))^4.735-0.44); end go=uo*rouo; dps=(D*(NB+1)*fo*0.15/De+10)*go^2/(2000*D0); |
回復(fù)此樓» 猜你喜歡
274求調(diào)劑
已經(jīng)有4人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有7人回復(fù)
-
材料專業(yè)求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有10人回復(fù)
-
344求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑
已經(jīng)有7人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有4人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助��！求一程序，用matlab程序做，用蒙特卡洛方法模擬
已經(jīng)有10人回復(fù)
- 【求助】對程序設(shè)計語言非常喜歡，但是我不會，想?yún)⒓觽�培訓(xùn)機(jī)構(gòu)，給點(diǎn)建議，感激不盡
已經(jīng)有13人回復(fù)
- 【求助】翅片式列管換熱器換熱面積問題
已經(jīng)有5人回復(fù)
- 【求助】關(guān)于GAMBIT中質(zhì)量不好的網(wǎng)格手動修改
已經(jīng)有5人回復(fù)
- 【求助】盤管換熱器的K值
已經(jīng)有4人回復(fù)
- 【求助】合成氨過程模擬 結(jié)果錯誤，怎么辦？并問下如何修改收斂次數(shù)！【已解決】
已經(jīng)有11人回復(fù)
- 【求助】物料屬性可以用UDF修改嗎
已經(jīng)有8人回復(fù)
- 【求助】換熱器計算
已經(jīng)有13人回復(fù)
2樓2011-05-23 23:17:50
只為找工作
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 417.2
- 散金: 20
- 帖子: 77
- 在線: 40.1小時
- 蟲號: 3518244
- 注冊: 2014-11-04
- 性別: GG
- 專業(yè): 內(nèi)流流體力學(xué)
3樓2015-12-15 22:10:54
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 892226227 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +3 | 守候夕陽CF 2026-03-18 | 3/150 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 281求調(diào)劑（0805） +5 | 煙汐憶海 2026-03-16 | 13/650 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|