气动仿真作业2008-5-14介绍.doc

气动仿真作业 2008-5-14 ——张业明 BY0703130 1 题目描述 条件：给定一个体积为30L的容器，对其充气， 声速流导C=0.2dm3/(s·bar)； 有效截面积Se=1mm3； 临界压力比b=0.4； 入口压力Ps=0.5MPa； 室温为20℃ 大气压力为1bar。 问题： 多变过程（n=1，n=1.4），求P（t）和θ（t）； 传热率h=20W/m2·K。 2 问题分析 （1）多变过程 多变过程中满足的方程式如下： (1) 这里，n是多变指数。等温过程时，n=1；绝热过程时，n=k=1.4。 对（1）进行微分，可得 (2) 质量流量G的公式: (3) 质量流量G的单位为[kg/s]。 式中：标准状态下的体积流量[dm3/s(ANR)] C ：声速流导[dm3/(s·bar)] P1 ：上流绝对压力[bar] P2 ：下流绝对压力[bar] ：上流空气的绝对温度[K] b ：临界压力比 ：标准状态下的密度为1.185[kg/m3] 声速流导C和有效截面积Se尽管单位不同，但都是表示流通能力的参数。两者之间存在如下的换算关系： Se[mm2]=5.03×C[dm3/(s·bar)]≈ 5×C[dm3/(s·bar)] (4) 由方程（4）可以求得声速流导C,再由（3）可以得到质量流量G，将G代入（2）得到关于t的微分方程。求解得到。 （2）传热过程（热传率一定） 根据能量守恒定律，容腔内空气的内能的变化量dU取决于流进或流出的空气的焓、外界的热交换量 （5） 容腔壁通常由金属制成，金属的热容量比空气的热容量大很多，金属内壁的温度可假设为不受空气充放气影响，与室温相同。这样，容腔内空气与容腔内壁的热交换可用下式表示 （6） 这里，是空气与壁面的热传率，是容腔内壁的表面积。将（6）代入（5），可得 （7） 这里，是上流的温度，充气的时候是气源温度，放气的时候是容腔内的空气温度。将（7）左侧展开后整理得到 （8） 由压力微分方程（状态方程） 对时间微分 (9) 由（8）和（9）联立后整理，得到如下的微分方程组 （10） 是上流的温度，充气的时候是气源温度，放气的时候是容腔内的空气温度；是室温。 将G代入（10）求得。 3 问题的Matlab仿真求解 （1）多变过程的仿真程序及仿真曲线 主求解程序solve1SYS.m function solve1SYS() p=1;t0=293; options = odeset(RelTol,1e-4,AbsTol,[1e-4 1e-4]); [T,Y] = ode45(@solve1,[0 5],[p t0],options); subplot(2,1,1); plot(T,Y(:,1),r-); xlabel(t/s),ylabel(P/bar);title(P(t));legend(P(t)); subplot(2,1,2); plot(T,Y(:,2),b-); xlabel(t/s),ylabel(θ/K);title(θ(t));legend(θ(t)); 问题描述文件solve1.m function dy=solve1(t,y) %空气常数：R=287[J/(kg.K)]%容腔体积：v=30[L]%声速流导：C=0.2[dm^3/(s.bar)]%有效截面积：Se=5.03*C=1.006 k=1.4; %绝热 %k=1; %等温 C=0.2; %[dm^3/(s.bar)] R=287; %[J/(kg.K)] V=30; %[L] p1=6; %[bar] t1=293; %[K] b=0.4; % density=1.185e-3 ; %[kg/dm*3] dy=zeros(2,1); if(y(1)/p1=b) G=density*C*p1*sqrt(293/t1); elseif(y(1)/p1=1) G=density*C*p1*sqrt(293/t1)*sqrt(1-((y(1)/p1-b)/(1-b))^2); else G=0; end dy(1)=k*G*R*y(2)/V; dy(2)=(k-1)*G*R*y(2)^2/y(1)/V; 绝热条件下的仿真结果： 等温条件下的仿真结果： (2)传热过程的仿真程序及仿真曲线 主求解程序solve2SYS.m function solve2SYS() p=