热工模化与控制第三次作业热过程模化与控制程序

1、热工模化与控制第三次作业1、求解方程：炉气温度数学模型：Tg g cg Tg g cg vgBQdWn cy t yh2ATgTwTy0.7Tg0.3Tw炉墙传热过程控制方程及定解条件：钢板分为80个点，粘土砖分为232个点:k 1 k 1Tw2 T w1k 1八、1 , kgh2 Tw1Ta ?X点 280:Twk1点 280:Twk1TkwiXkgk 1 k 1Twi 1TwikgXx。1 T：点81 （交界点）k 1 k 1Twi 1TwikgxT k 1 t kwi 1 wi1一 ktcg2T k 1 T kwi wiX；点82312:T kwiTx1-ktk 1 k 1Twi 1Tw

2、iktxxTwk1Twk ；k 1 k 1点 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;/热物性参数赋值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炉温达到800c时间为j秒, 炉子外壁温度为Tw1 C =850&Tgas=800&Tgas=800.5)cout炉温达到800c时间为j秒, 炉子外壁温度为Tw1 C =850&Tgas=850.4)cout炉温达到850c时间为j秒, 炉子外壁温度为Tw1 C =900&Tgas=900.4)cout炉温达到900c时间为j秒, 炉子外壁温度为Tw1 C =950&Tgas=950.3)cout炉温达到950c时间为j秒, 炉子外壁温度为Tw1 C endl;此时炉子内壁温度为此时炉子内壁温度为此时炉子内壁温度为此时炉子内壁温度为Twn CTwn CTwn

7、 CTwn=1000&Tgas=1000.3)cout炉温达到1000 c时间为j秒，此时炉 Twn C炉子外壁温度为 Tw1 C =1050&Tgas=1050.3)cout炉温达到1050 c时间为j秒，此时炉 Twn C炉子外壁温度为 Tw1 C =1100&Tgas=1100.3)cout炉温达到1100 C时间为j秒，此时炉子内壁温度为 Twn C炉子外壁温度为 Tw1 C =1100)break;.运行结果：当炉子热负荷 B=2.5m3/hr时:T金型 I型rsl门达ww达达 八温温温温温温温热所何加28001c时间为268G秒,“古时间为3274机 904时怛为混47沙, 一乖的

8、471包机 由旬为T金型 I型rsl门达ww达达 八温温温温温温温热所何加28001c时间为268G秒,“古时间为3274机 904时怛为混47沙, 一乖的471包机 由旬为557州9500 1000VQ _ 105 Hle时间为656堂;11腼七时间为7666秒.此时炉子内壁遑度为砥4上抽T炉子外辘*磊翳震理软此向第壬内鹭噢K之;就!l-b壁温度据博随嚣耶瑞招奉壁温度为2瓦眄87c 壁温度为壁温度为20.前42七 隼温度F 20皿2y ,瘗楠蔗为2必3147T壁温度为靖w”1C 、壁温度为如.缩g2七当炉子热负荷 B=5m 3/hr时:歹夕外子子子4 17 3 37度度Em EmJ10 6

9、175 5 6 6为为为度 一度矗温温温温温温台壁壁9 16 8 鲁 6 * 01 00度温温目nn黄 温壁壁谿产.子炉炉 r-6 5179 T 9 5 5 3 B-T-到，达达差达达达 A.l温温温温温温 户户户产p B- B- 一司 环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壁壁壁壁帽内内 内内内勒子子 干子子子吃炉炉 SSLB此此 hELLLIK 让让让up，为度霸-2 _28ropvvr r -r - 6X 9 2 2 Tx 2 6 12g 7 11 5 0 0 当炉子热负荷 B=7.5m3/hr时:炉达达达达、A.温温旧皿温温温温la0-0050 0 5 0 0 19 9 1 1 1,5为为口7m T 日.fml fmL- - gDderer 负 Im,xmr-rr?14 8 7 7 0 fi 3 1 4。402 3 4 7 112此时猖子内壁遇魅口舫辆工内壁温度为此时而券内量温度为蚯128式旅一忙若曜温盟15以868务力 守炉?