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

1、熱工模化與控制第三次作業(yè)1、求解方程：爐氣溫度數(shù)學(xué)模型：爐墻傳熱過(guò)程控制方程及定解條件：鋼板分為80個(gè)點(diǎn)，粘土磚分為232個(gè)點(diǎn)：點(diǎn)1：；點(diǎn)280: ;點(diǎn)81（交界點(diǎn)）：；點(diǎn)82312： ；點(diǎn)313： .定解條件：2、程序框圖：yn結(jié)束輸出加熱時(shí)間和爐墻內(nèi)外溫度根據(jù)爐氣微分方程，用歐拉預(yù)報(bào)修正法求爐溫確定各物性參數(shù)及初始條件輸入爐溫和燃料流量b開始根據(jù)爐墻差分方程用tdma求爐墻各個(gè)點(diǎn)溫度3.程序清單：#include#includeusing namespace std;void main()float bf,b1,qdw,vy,cy,my,twn,t0,ct,mt,kt,cg,mg,kg,v

2、,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,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ù)賦

3、值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;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=t

4、gas+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*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/d

5、t;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=800&tgas=800.5)cout爐溫達(dá)到800時(shí)間為j秒，此時(shí)爐子內(nèi)壁溫度為twn爐子外壁溫度為tw1=850&tgas=850.4)cout爐溫達(dá)到850時(shí)間為j秒，此時(shí)爐子內(nèi)壁溫度為twn爐子外壁溫度為tw1=900&tgas=900.4)cout爐溫達(dá)到900時(shí)間為j秒，此時(shí)爐子內(nèi)壁溫度為twn爐子外壁溫度為tw1=950&tgas=950.3)cout爐溫達(dá)到950時(shí)間為j秒，此時(shí)爐子內(nèi)壁溫度為twn爐子外壁溫度為tw1=1000&tgas=1000.3)cout爐溫達(dá)到1000時(shí)間為j秒，此時(shí)爐子內(nèi)壁溫度為twn爐子外壁溫度為tw1=1050&tgas=1050.3)cout爐溫達(dá)到1050時(shí)間為j秒，此時(shí)爐子內(nèi)壁溫度為twn爐子外壁溫度為tw1=1100&tgas=1100.3)cout爐溫達(dá)到1100時(shí)間為j秒，此時(shí)