本次汇报的主要内容为:根据一篇论文中换热器的物理简化模型,和已有的参数数据,对模型进行了重现(应用了欧拉迭代法)。换热器模型简化1)图1换热器謬简化的物理模型dTm^-=Q-Qdt12当管内、外介质温度不同时,管壁金属具有显著的蓄热性能。式中m—单位长度管壁的金属质量;c—金属比热;Q—管外放热介质对单位管长的管壁金属1在单位时间内的放热量;Q—管壁金属向管内介质在单位管长和单位时间内的2放热量。管内介质在单位管长和单位时间内的放热量可
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
示:Q=k兀d(T-T)二ks(T-T)(2)2222222式中d一管子内径;s—单位长度管子的内表面积;T和T一分别表示金属222和管内介质的温度;k—放热系数,可表示为k沁KDn=KD0.8(3),式中K是22222常数,D是管内介质的流量。管外介质在单位管长和单位时间内的放热量可表示:TOC\o"1-5"\h\zQ=Q+Q⑷11d1f式中Q—介质向单位管长管壁的对流放热量;Q—介质向单位管长管壁的1d1f辐射放热量。且有Q=k兀d(T-T)二ks(T-T)(5)1d1d1y1d1yQ=k兀d(T-T)二ks(T-T)(6)1f1f1y1f1y式中k、k—分别为对流和辐射放热系数;T—管外介质温度;d—管子1d1fy1外径;s—单位长度管子外表面积。把方程(2)-(6)代入方程(1)中得:1mc--=ks(T-T)+ks(T-T)-ks(T-T)(7)dt1d1y1f1y222令A=mc,T=T=T,B=ks,B=ks,B=ks(7)式可以简化为:13y11d122231f18)AdT=£B(T-T)dtiii=1£BT令T=丄=丄(时间常数),q=»二,(8)式最后简化为关于管壁温度TcpfBAii=1的微分方程(T0为初值):9)dT+pT二qdt9)式存在解析解:T=—-—(q-Tp)e-ptpp0(10)解析解中有指数的形式,占用机时较多。仿真计算中大多不用解析算法,而是从(9)式出发,构造差分格式,用一定精度的近似解代替解析解,达到实时仿真的目的。差分格式算法(前进欧拉法)2.1基本原理方法前进欧拉法是基于向前微分近似(Tn+1-Tn)沁T'变形为:hnT=T+hf(T,t)T(t)=T(i=0,1,2,,n-1)i+1iii00其中T'=f(T,t),以n循环得到…iiiT=T+hT'100T=T+hf(T,t)1000T=T+hf(T,t)2111T=T+hf(T,t)•・nn-1•••n-1n-12.2具体构造的欧拉法迭代
公式
小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载
:£B(T-T)iinqT=T+h—二T+hp(-T)n+1nAnpn%显示欧拉法程序代码:(基本参数为p=0.02;q=10;h=1;T0=20;)%T=dsolve('DT+p*T-q=0','T(0)=T0','t')%T=(q-(q-T0*p)/exp(p*t))/p%微分方程的函数文件:functionz=mTf4(x,T)p=0.02;q=10;z=q-p*T;%准确解的函数文件(解析解):functionfz=fz(x)p=0.02;q=10;T0=20;fz=q/p-(1/p)*(q-T0*p)./exp(p*x);%欧拉格式算法:%初始值和步长假定p=0.02;q=10;h=1;T0=20;x=zeros(1,301);T=zeros(1,301);x(1)=0;T(1)=20;forn=1:300x(n+1)=x(n)+h;T(n+1)=T(n)+h*feval(@mTf4,x(n),T(n));endE=[x,T];T1=fz(x);plot(x,T,':og',x,T1,'--b');title('Euler格式与准确解比较图');Euler格式与准确解比较图500450400350300250200150100500050100150200250300从图中可以看出欧拉迭代法的近似解和微分方程的精确解的误差很小,从而证明了欧拉迭代算法满足了题目要求。存在问题1)迭代的差分格式算法替代微分方程的精确解,仅仅是因为计算时间问题?,还是有其他原因?2)原论文中采取了其他迭代算法(改进的欧拉法、梯形法和四阶龙格库塔法),这里只是选取了一种算法进行重现,也没有进行算法稳定性的对比
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
。汇报人:杨鹏志