内燃机工作过程计算一、内燃机实际工作过程的数值计算1、基本原理与公式在推导气缸工作过程计算的基本微分方程式时,作如下的简化假定:①不考虑气缸内各点的压力、温度和浓度的差异,并认为在进气期间,流入气缸内的空气与气缸内的残余废气实现瞬时的完全混合,缸内状态是均匀的,亦即为单区计算模型。②工质为理想气体,其比热、内能仅与气体的温度和气体的组成有关。③气体流入或流出的气流为准稳定流动。④不计进气系统内压力和温度波动的影响。⑴能量平衡方程式根据热力学第一定律,得出能量守恒方程的一般形式如下:m:工质的质量u:工质比内能Q:燃烧放出的热量Q:通过系统边界传入或传出的热量W:作用在活塞上的机械功m:通过进气阀流入气缸的质量m:通过排气阀流出气缸的质量h和h:分别为进气阀前和气缸内气体的比焓⑵质量平衡方程式通过系统边界的质量有:喷入气缸的瞬时燃油质量、流入气缸的气体质量及流出气缸的气体质量,质量平衡的微分方程可写为:⑶气体状态方程式⑷气缸工作容积气缸工作容积()气缸瞬时容积()曲柄连杆比()气缸容积变化率(/度)1⑸传热计算公式气缸中的气体通过活塞顶面、气缸盖底面及气缸套的瞬时传热面进行热量传递。为简单起见通常将整个传热面分为活塞、气缸盖及气缸套三部分,用三部分传热面的平均壁温及按位置平均瞬时传热系数进行传热量的计算,计算公式为:(J/度):传热表面积i=1和i=2活塞顶部和气缸盖()i=3气缸套():传热表面的平均温度i=1活塞顶部(k)i=2气缸盖i=3气缸套注:为平均有效压力(bar)其中:为冲程数,额定功率(kw),为工作容积(),n发动机转速(转/分),i气缸数。传热系数经验公式:1、G.Sitkei公式alf=0.205*(1+b)*pow((P/1000000),0.7)*pow(Cm,0.7)/(pow(T,0.2)*pow(de,0.3));注:b:为系数。对直喷式的为0~0.15,本设计选用0.09。2无量纲。(m)h为活塞顶到气缸盖底面距离。(m/s)活塞平均速度p气缸内工质绝对压力(MPa)T气缸内工质绝对温度(K)2、佛劳姆(W.Pflanm)公式注:符号和单位同上⑹燃烧规律对于一般性能预测和增压柴油机的匹配计算,采用代用燃烧规律可满足要求。经常使用的代用燃烧函数是韦柏的半经验公式:x:燃料燃烧百分比y:无因次时间函数:燃烧起始角:瞬时曲轴转角:燃烧持续角m:燃烧品质指数由此可算出气缸中燃料燃烧时的瞬时放热率(J/度)3:每循环喷油量(kg):燃烧效率,对于柴油机稳定运算,可取⑺工质内能变化的计算柴油机工质为空气和燃烧废气的混合物,一般常用瞬时过量空气系数来表示,这样,工质的内能可表示为两个变量的函数可用于工质内能变化的经验公式很多,下面用的是Justi公式:(J/kg)⑻机械功的计算公式⑼气缸温度3、气缸内各过程的计算计算一般从压缩始点开始,压缩、燃烧和膨胀三部分为高压阶段,从排气门打开到进气门打开为止为排气期,以及延续到气门迭开期,接着是延续到进气门关闭的吸气期,这三个时期为工质更换阶段。⑴压缩期⑵燃烧期由于燃烧气缸内工质质量发生变化:4任一时刻气缸内工质质量为:工质成分也因燃烧发生变化,瞬时过量空气系数可得:⑶膨胀期膨胀期气缸工质的质量不变,方程式和压缩期的一样,只是各项数值不同。⑷换气过程下标u表示上游,d表示下游;s表示进气,e表示排气。54、常微分方程的数值解法常用的数值解法有三种:欧拉法、改进欧拉法和龙格—库塔法。⑴欧拉法:⑵改进欧拉法:⑶龙格—库塔法(略)5、工作过程计算程序编程步骤(1)定义初始和全局变量(2)编写Qb,Qw,dV,C,Cv各量的子函数(3)从进气门关闭角开始计算,首先给定该角度的气缸压力、温度、残余废气系数(初始条件),计算出该角度下的体积、工质质量和(2)步骤的各子函数,以及该角度下的温度变化率。(4)由式计算角度的温度,通常取=1度。(5)计算角度的体积、工质质量,再由计算此时的压力。返回到步骤(3),可根据角度的值继续下一个角度的计算直到排气门打开,结束计算。(6)将计算的结果保存到“outdata.txt”文件中。二、应用Lotus软件计算内燃机工作过程61、点图标进入LotusEngineSimulation软件。72、选LotusEngineSimulation和Newblank.Simfi...