发动机缸内工作过程计算资料.docVIP

发动机工作过程数值计算作业 缸内封闭过程数值计算 学 院：汽车学院 专 业：动力机械及工程 姓 名：滑海宁 学 号：2012122056 教 师：张春化 二〇一三年五月 气缸封闭过程的数值计算 发动机缸内工作过程的数值计算，是以内燃机缸内各工作阶段的物理模型为基础， 通过微分方程来对其各阶段工作过程进行数学描述， 然后通过程序编写求解微分方程， 得到缸内温度、压力等参数随曲轴转角的变化曲线。 基本热力学模型 图1 发动机缸内热力系统 在简化假设的基础上，取气缸为一个热力学系统，如图1所示。这个热力系统包括了质量交换项， 如排气 dmA，进气 dmE，喷入气缸内瞬时燃料质量 dmB ； 与能量交换项，如焓变， 功， 燃烧放热等。 图中 T 、P 、V、 m 及 u 分别为缸内瞬时气体温度、压力、体积、燃料质量及比内能。 简化假设 1.假定工质为理想混合气体； 2.假定缸内各处温度、压力及混合气浓度均匀； 3.用纯空气的气体常数代替混合气气体常数； 4.假定扫气完全，即不考虑残余废气； 5.不计漏气损失，并假定只有在燃烧始点才有燃油喷入气缸； 6.按代用燃烧规律进行喷油，并认为着火延迟等于零； 7.假定放热率为100%。 基本方程 1、气体状态方程 PV=mRT 2、压缩期() 能量方程：＝（－p） 质量方程：＝0 3、燃烧期() 能量方程：＝（＋－p－u－m） 质量方程：＝ 4、膨胀期() 能量方程：＝（－p） 质量方程：＝0 其他计算公式 1、气缸瞬时工作容积： = 式中：-气缸行程容积 VL-发动机排量 -压缩比 -曲柄连杆比 -活塞行程 -气缸直径 2、气缸容积随曲轴转角变化率： 3、单位曲轴转角的传热量： 式中：-活塞 -气缸盖 -气缸套 -传热表面积 ；； -传热表面平均温度 ；； -换热面平均的瞬时传热系数 利用Woschni公式： 式中：-活塞行程容积 -气体速度系数 压缩、膨胀阶段 换气阶段 -燃烧室形状系数，直喷式燃烧室=3.24X10-3 -活塞平均速度，计算得6.9m/s 、、分别为压缩始点的气缸压力、容积、温度。分别取值为102kPa, 0.89L, 375K -发动机倒拖时的气缸压力,取为100kPa 4、代用放热规律： 累积放热量： 瞬时放热率： 式中：-燃烧效率,取为100% -燃料低热值，取为42500kJ/kg -燃烧品质指数,取为1.0 -燃烧持续角，取为80°CA -循环喷油量,为3.5421X10-6kg -燃烧始点的曲轴转角，为338°CA 5、工质内能： u=u(,T ） 式中：-瞬时过量空气系数 1）压缩阶段，取。 2)燃烧阶段： 3)膨胀阶段： =const 式中：-理论空气量，取14.3kg/kg -喷入燃料量随曲轴转角的变化函数 -每循环喷油量 6、等容比热： 7、比内能对过量空气系数的偏导 8、燃烧期混合气的质量： mH=mL+mB 龙格—库塔法的计算步骤 已知 ， 步骤一： 步骤二： 步骤三： 步骤四： 步骤五： 步骤六： 步骤七： 步骤八： 步骤九： 步骤十： 步骤十一： 步骤十二： 步骤十三： 发动机基本参数 发动机形式：四冲程、直喷、水冷、单缸柴油机 发动机型号：ZH105W 压缩比：＝16.5 排量：=0.96L 活塞行程：s＝115mm 缸径：D＝105mm 循环喷油量：3.5421X10-5kg 燃烧持续角：80°CA 曲柄半径和连杆长度比 ＝0.3 喷油提前角： 配气相位： 进气提前（上止点前）12°CA 进气迟后（下止点后）38°CA 排气提前（下止点前）55°CA 排气迟后（上止点后）12°CA 程序框图 源程序（matlab） 压缩过程 自定义ysgc函数并生成m文件，源程序如下： function f=ysgc(fai,T) P1=102000;%压缩始点气缸内气体压力，单位pa V1=0.00089;%压缩始点气缸内气体体积(计算得)，单位m3 T1=375;%压缩始点气缸内气体的温度（计算得），单位K R=8.314;%气体常数，单位J/mol.k ML=29;%空气摩尔质量，单位g/mol nL=P1*V1/R/T1;%空气摩尔量，单位mol mL=nL*ML;%压缩始点气缸内气体质量，单位g lmd=10000;%过量空气系数 D=0.105;%气缸直径，单位m c1=2.29;%气体速度系数 cm=6.9;%活塞平均速度（计算得），单位m/s c2=