发动机工作过程作业陈斌.doc

发动机工作过程作业 缸内封闭过程数值计算 学 院：汽车学院 专 业：流体机械及工程 姓 名：陈斌 学 号：2014122059 教 师：张春化 二〇一五年七月 气缸封闭过程数值计算 发动机缸内工作过程的数值计算，是以内燃机缸内各工作阶段的物理模型为基础， 通过微分方程来对其各阶段工作过程进行数学描述， 然后通过程序编写求解微分方程， 得到缸内温度、压力等参数随曲轴转角的变化曲线。 基本热力学模型为在简化假设的基础上，取气缸为一个热力学系统。这个热力系统包括了质量交换项， 如排气 dmA，进气 dmE，喷入气缸内瞬时燃料质量 dmB ； 与能量交换项，如焓变， 功， 燃烧放热等。 简化假设 1.假定工质为理想混合气体； 2.假定缸内各处温度、压力及混合气浓度均匀； 3.用纯空气的气体常数代替混合气气体常数； 4.假定扫气完全，即不考虑残余废气； 5.不计漏气损失，并假定只有在燃烧始点才有燃油喷入气缸； 6.按代用燃烧规律进行喷油，并认为着火延迟等于零； 7.假定放热率为100%。 二、基本方程 1、气体状态方程 PV=mRT 2、压缩阶段总微分方程 能量方程：＝（－p） 质量方程：＝0 3、燃烧阶段总微分方程 能量方程：＝（＋－p－u－m ） 质量方程：＝ 4、膨胀阶段总微分方程 能量方程：＝（－p） 质量方程：＝0 其他计算公式 1、气缸瞬时工作容积： V（φ）= －气缸工作容积 ε－压缩比 λs－曲柄半径连杆比 2、气缸容积随曲轴转角变化率： 3、传热计算公式： Ai－传热表面积，Twi－传热表面平均温度; α－按传热面平均的瞬时传热系数 A1=л*D2/4; A2=1.5*A1; A3=л*D*h TW1=560K; TW2=548K; TW3=548K 利用Woschni公式： －活塞行程容积 —气体速度系数,取为2.99 —燃烧室形状系数 —活塞平均系数，计算得8.4m/s P1 V1 T1——分别为压缩始点的气缸压力、容积、温度 分别取值为90kPa, 1.8249*10-3m3, 340K ——发动机倒拖时的气缸压力,取为100kPa 4、放热规律： 累积放热量： 瞬时放热率： －燃烧效率,取为100% HU－燃料低热值，取为42500kJ/kg m－燃烧品质指数,取为1.0 △φ －燃烧持续角，取为80°CA mB0—循环喷油量,为1.18474*10-4kg —燃烧始点的曲轴转角（350—359°CA，即上止点前1—10°CA,360-喷油提角） 5、工质内能表达式 u=f(,T ） -瞬时过量空气系数 压缩阶段可取，由此计算的比内能与时的比内能值差别很小。 燃烧阶段： 膨胀阶段： —理论空气量 取14.3kg/kg —喷入燃料量随曲轴转角的变化函数 6、瞬时比热 7、比内能对过量空气系数的偏导 四、发动机基本参数 发动机形式: 四冲程、直列、水冷、六缸柴油机，W形燃烧室 发动机型号：6135Q 压缩比 ＝16.5 活塞行程 s＝140mm 缸径 D＝135mm 曲柄半径和连杆长度比 ＝0.3 喷油提前角 8°CA 配气相位： 进气提前（上止点前）24°CA 进气迟后（下止点后）53.5°CA 排气提前（下止点前）52°CA 排气迟后（上止点后）26°CA 五、程序框图 六、源程序（C++） #includeiostream.h #includemath.h #define M 1.0 #define FA 80 #define P1 90 #define T1 340 #define FAVB 8 #define TW1 560 #define TW2 548 #define TW3 548 class P {public: double m; P(){m=0;} P(double a):m(a) {} double fv(double fai) {double FAI=fai*3.1416/180; return 3.1416*0.135*0.135*0.14/8*(2/15.5+1-cos(FAI)+(1-sqrt(1-0.09*sin(FAI)*sin(FAI)))/0.3); } double dv(double fai) {double FAI=fai*3.1416/180; return 3.1416*0.135*0.135*0.14/8*(sin(FAI)+0.3*sin(FAI)*cos(FAI)/sqrt(1-0.09*sin(FAI)*sin(FAI))); } do