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

此文档收集于网络，如有侵权，请联系网站删除热工模化与控制第三次作业1、求解方程：炉气温度数学模型：炉墙传热过程控制方程及定解条件：钢板分为80个点，粘土砖分为232个点：点1：；点280: ;点81（交界点）：；点82312： ；点313： .定解条件：2、程序框图：YN结束输出加热时间和炉墙内外温度根据炉气微分方程，用欧拉预报修正法求炉温确定各物性参数及初始条件输入炉温和燃料流量B开始根据炉墙差分方程用TDMA求炉墙各个点温度3.程序清单：#include#includeusing namespace std;void main()float 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,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;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*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=800&Tgas=800.5)cout炉温达到800时间为j秒，此时炉子内壁温度为Twn炉子外壁温度为Tw1=850&Tgas=850.4)cout炉温达到850时间为j秒，此时炉子内壁温度为Twn炉子外壁温度为Tw1=900&Tgas=900.4)cout炉温达到900时间为j秒，此时炉子内壁温度为Twn炉子外壁温度为Tw1=950&Tgas=950.3)cout炉温达到950时间为j秒，此时炉子内壁温度为Twn炉子外壁温度为Tw1=1000&Tgas=1000.3)cout炉温达到1000时间为j秒，此时炉子内壁温度为Twn炉子外壁温度为Tw1=1050&Tgas=1050.3)cout炉温达到1050时间为j秒，此时炉子内壁温度为Twn炉子外壁温度为Tw1=1100&Tgas=1100.3)cout炉温达到1100时间为j秒，此时炉子内壁温度为Tw