熱工?；c控制第三次作業(yè)熱過(guò)程模化與控制程序

1、熱工?；c控制第三次作業(yè)1、求解方程：爐氣溫度數(shù)學(xué)模型：Tg g cg Tg g cg vgBQdWn cy t yh2ATgTwTy0.7Tg0.3Tw爐墻傳熱過(guò)程控制方程及定解條件：鋼板分為80個(gè)點(diǎn)，粘土磚分為232個(gè)點(diǎn):k 1 k 1Tw2 T w1k 1八、1 , kgh2 Tw1Ta ?X點(diǎn) 280:Twk1點(diǎn) 280:Twk1TkwiXkgk 1 k 1Twi 1TwikgXx。1 T：點(diǎn)81 （交界點(diǎn)）k 1 k 1Twi 1TwikgxT k 1 t kwi 1 wi1一 ktcg2T k 1 T kwi wiX；點(diǎn)82312:T kwiTx1-ktk 1 k 1Twi 1Tw

2、iktxxTwk1Twk ；k 1 k 1點(diǎn) 313:kt -wnh1 Tg Tjn1x定解條件：Ta 20,h1 100, h225.2、程序框圖:.程序清單：#include#include using namespace std; void main() #includefloat Bf,B1,Qdw,vy,cy,my,Twn,T0,ct,mt,kt,cg,mg,kg,v,h1,h2,Af,gx0y0,gxy,Ty,s1,s2,dt,dx;float Tgas,Tgas1,mgas,cgas;int n1,n2,n;float Tw400,A400,B400,C400,F400,U400

3、,V400;Qdw=41800000;Tgas=20; mgas=1.35; cgas=1300; cy=1250;my=1.25;Twn=20;T0=20;ct=960;mt=2000;kt=1.05;cg=485;mg=7800;kg=60;v=0.125;h1=100;h2=25;Af=1.5;Ty=20;dt=1;dx=0.001;s1=232;s2=80;/熱物性參數(shù)賦值coutB1;Bf=B1/3600;vy=Bf*10.975;n1=(int)s1/(dx*1000);n2=(int)s2/(dx*1000);n=n1+n2+1;for(int i=1;i=n;i+)Twi=T0

4、;for(int j=1;j+)Twn=Twn;gx0y0=(Bf*Qdw-vy*cy*(0.7*Tgas+0.3*Twn)-h1*Af*(Tgas-Twn)/(mgas*cgas*v);Tgas1=Tgas+gx0y0*dt;gxy=(Bf*Qdw-vy*cy*(0.7*Tgas1+0.3*Twn)-h1*Af*(Tgas1-Twn)/(mgas*cgas*v);Tgas=Tgas+0.5*dt*(gxy+gx0y0);B1=-kg-h2*dx;C1=kg;F1=-h2*dx*T0;U1=-C1/B1;V1=F1/B1;for(i=2;i=n2;i+)Ai=kg;Bi=-2*kg-mg*cg

5、*dx*dx/dt;Ci=kg;Fi=-mg*cg*dx*dx*Twi/dt;)An2+1=kg;Bn2+1=-kg-kt-(mt+mg)*(ct+cg)*dx*dx/(4*dt);Cn2+1=kt;Fn2+1=-(mt+mg)*(ct+cg)*dx*dx*Twn2+1/(4*dt); for(i=n2+2;i=n-1;i+)Ai=kt;Bi=-2*kt-mt*ct*dx*dx/dt;Ci=kt;Fi=-mt*ct*dx*dx*Twi/dt;)An=kt;Bn=-kt-h1*dx;Fn=-h1*dx*Tgas;for(i=2;i=1;i-)Twi=Ui*Twi+1+Vi;if (Tgas=80

6、0&Tgas=800.5)cout爐溫達(dá)到800c時(shí)間為j秒, 爐子外壁溫度為Tw1 C =850&Tgas=800&Tgas=800.5)cout爐溫達(dá)到800c時(shí)間為j秒, 爐子外壁溫度為Tw1 C =850&Tgas=850.4)cout爐溫達(dá)到850c時(shí)間為j秒, 爐子外壁溫度為Tw1 C =900&Tgas=900.4)cout爐溫達(dá)到900c時(shí)間為j秒, 爐子外壁溫度為Tw1 C =950&Tgas=950.3)cout爐溫達(dá)到950c時(shí)間為j秒, 爐子外壁溫度為Tw1 C endl;此時(shí)爐子內(nèi)壁溫度為此時(shí)爐子內(nèi)壁溫度為此時(shí)爐子內(nèi)壁溫度為此時(shí)爐子內(nèi)壁溫度為Twn CTwn CTwn

7、 CTwn=1000&Tgas=1000.3)cout爐溫達(dá)到1000 c時(shí)間為j秒，此時(shí)爐 Twn C爐子外壁溫度為 Tw1 C =1050&Tgas=1050.3)cout爐溫達(dá)到1050 c時(shí)間為j秒，此時(shí)爐 Twn C爐子外壁溫度為 Tw1 C =1100&Tgas=1100.3)cout爐溫達(dá)到1100 C時(shí)間為j秒，此時(shí)爐子內(nèi)壁溫度為 Twn C爐子外壁溫度為 Tw1 C =1100)break;.運(yùn)行結(jié)果：當(dāng)爐子熱負(fù)荷 B=2.5m3/hr時(shí):T金型 I型rsl門達(dá)ww達(dá)達(dá) 八溫溫溫溫溫溫溫?zé)崴渭?8001c時(shí)間為268G秒,“古時(shí)間為3274機(jī) 904時(shí)怛?yàn)榛?7沙, 一乖的

8、471包機(jī) 由旬為T金型 I型rsl門達(dá)ww達(dá)達(dá) 八溫溫溫溫溫溫溫?zé)崴渭?8001c時(shí)間為268G秒,“古時(shí)間為3274機(jī) 904時(shí)怛?yàn)榛?7沙, 一乖的471包機(jī) 由旬為557州9500 1000VQ _ 105 Hle時(shí)間為656堂;11靦七時(shí)間為7666秒.此時(shí)爐子內(nèi)壁遑度為砥4上抽T爐子外轆*磊翳震理軟此向第壬內(nèi)鷺噢K之;就!l-b壁溫度據(jù)博隨囂耶瑞招奉壁溫度為2瓦眄87c 壁溫度為壁溫度為20.前42七 隼溫度F 20皿2y ,瘞楠蔗為2必3147T壁溫度為靖w”1C 、壁溫度為如.縮g2七當(dāng)爐子熱負(fù)荷 B=5m 3/hr時(shí):歹夕外子子子4 17 3 37度度Em EmJ10 6

9、175 5 6 6為為為度 一度矗溫溫溫溫溫溫臺(tái)壁壁9 16 8 魯 6 * 01 00度溫溫目nn黃 溫壁壁谿產(chǎn).子爐爐 r-6 5179 T 9 5 5 3 B-T-到，達(dá)達(dá)差達(dá)達(dá)達(dá) A.l溫溫溫溫溫溫 戶戶戶產(chǎn)p B- B- 一司 環(huán)5792 7 9 14 5 13 7 8 2 7 9-11 3 4 6 71一吊同 s fm fm 1 日 1 -B 丁二二、-1、上日.0 0 曜口。翼眈吹 n .- 0 0 0 0 0 0 0505001 J28 8 9 9 1 1 1壁壁壁壁帽內(nèi)內(nèi) 內(nèi)內(nèi)內(nèi)勒子子 干子子子吃爐爐 SSLB此此 hELLLIK 讓讓讓up，為度霸-2 _28ropvvr r -r - 6X 9 2 2 Tx 2 6 12g 7 11 5 0 0 當(dāng)爐子熱負(fù)荷 B=7.5m3/hr時(shí):爐達(dá)達(dá)達(dá)達(dá)、A.溫溫舊皿溫溫溫溫la0-0050 0 5 0 0 19 9 1 1 1,5為為口7m T 日.fml fmL- - gDderer 負(fù) Im,xmr-rr?14 8 7 7 0 fi 3 1 4。402 3 4 7 112此時(shí)猖子內(nèi)壁遇魅口舫輛工內(nèi)壁溫度為此時(shí)而券內(nèi)量溫度為蚯128式旅一忙若曜溫盟15以868務(wù)力 守爐?