内燃机燃烧放热程序及结果.pdf

内燃机燃烧放热率计算分析 姓名：吴彬；学号：J200864681 ；考生类别：在职硕士；所在院系：能源学院 经测量获得缸内压力变化曲线见下图左，采用七点五次光顺后的缸内压力变化 曲线见下图右。 图1 缸内气体压力变化曲线 采用放热计算程序（见附录），计算后的到燃烧过程放热率曲线、燃油消耗率 曲线以及缸内压力变化曲线，分别见图2 - 4 。 图2 缸内燃烧过程放热率曲线 图3 随曲轴转角变化燃油消耗率变化曲线 图4 随曲轴转角变化缸内温度变化曲线 附录 燃烧计算程序 close all; clear;clc; %从文件读取数据 [filename,pathname]=uigetfile({*.txt;*.dat},选择示功图数据文件); str_filename=[pathname,filename]; [Deg,A] = textread(str_filename,%s %s,headerlines,2);%读取数据 Deg= str2num(char(Deg)); %读取角度数据 A = str2num(char(A)); %读取未光顺压力值 Mpa %七点五次光顺 N=numel(A); Y(1)=(39*A(1)+8*A(2)-4*(A(3)+A(4)-A(6))+A(5)-2*A(7))/42; Y(2)=(8*A(1)+19*A(2)+16*A(3)+6*A(4)-4*A(5)-7*A(6)+4*A(7))/42; Y(3)=(-4*(A(1)+A(6))+16*A(2)+19*A(3)+12*A(4)+2*A(5)+7*A(7))/42; if N7 return; end for k=4:(N-3) Y(k)=(-2*(A(k-3)+A(k+3))+3*(A(k-2)+A(k+2))+6*(A(k-1)+A(k+1))+7*A(k))/21; end Y(N-2)=(A(N-6)-4*(A(N-5)+A(N))+A(N-4)+12*A(N-3)+19*A(N-2)+16*A(N-1))/42; Y(N-1)=(4*A(N-6)-7*A(N-5)+4*A(N-4)+6*A(N-3)+16*A(N-2)+19*A(N-1)+8*A(N))/42; Y(N)=(-2*A(N-6)+4*(A(N-5)-A(N-3)-A(N-2))+A(N-4)+8*A(N-1)+39*A(N))/42; %预设初始计算参数 gc=0.870; gh=0.126; go=0.004; %柴油平均质量成分（碳、氢、氧） tao=4; fai=1.06; %冲程、扫气系数、 n=2500; i=6; %转速，气缸数 bc=235.6; %燃油消耗率 Pe=161.51; %有效功率 Gma=0.02; %残余废气系数 d=0.102; %缸径 s=0.12; %冲程 skm=17.3; %压缩比 cor=0.192; %连杆长