《XD479Q发动机模拟平台搭建.doc》由会员分享，可在线阅读，更多相关《XD479Q发动机模拟平台搭建.doc（77页珍藏版）》请在三一办公上搜索。
1、目 录第一章 绪论21.1引言21.2数值模拟31.3模拟计算在内燃机发展中的作用41.4内燃机模拟计算的现状及发展51.5 GT-POWER软件介绍61.6本文主要内容及研究意义8第二章 流体流动模型基本原理92.1软件的基本原理92.1.1 管路（PIPE）模块基本原理112.1.2 接头（FLOWSPLIT）模块基本原理122.2缸内流动仿真原理132.3燃油喷射模块基本原理142.4燃烧模型基本原理142.5发动机摩擦及驱动附件损失模块基本原理15第三章 XD479Q发动机模拟平台搭建163.1 XD479Q发动机参数163.2整机模型的搭建173.3定义外部环境参数183.4定义空气
2、滤清器193.5 定义节气门203.6定义连接管件及孔实体213.7进气系统的创建及定义223.8定义进、排气门253.6气缸参数设定313.7燃油喷射模块设定353.9排气系统的创建及定义373.9发动机曲轴箱设定403.10原机外特性仿真结果及模型验证42第四章 配气相位对发动机性能影响464.1进气门处压力分析计算464.2配气相位对性能的影响474.3 XD479Q发动机配气相位的优化484.3.1 进气凸轮正时494.3.2 排气凸轮正时524.3.3配气相位角的正交试验554.4 优化后配气相位的验证58第五章 新凸轮型线设计605.1缓冲段曲线的设计605.2 工作段运动曲线设计
3、62第六章 总结65参考文献66附录I 气门升程随曲轴转角变化关系表68附录 进气凸轮升程曲线MATLAB源程序70附录 排气凸轮升程曲线MATLAB源程序73致 谢76第一章 绪论1.1引言随着国家对汽车油耗、排放等法规要求越来越严格，汽车企业生产的产品也不得不采取各种办法对产品进行技术升级和改造，以适应国家法规的要求。09年1月1日，燃油税开始实施。燃油税是指政府对燃油在零售环节征收的专项性质的税收。通过征税的办法从油价中提取一定比例作为养路等费用。其基本原理是：车辆类型及行驶里程长短，载货量大小是与耗油量的多少紧密相连的，耗油越多证明其享有使用公路的权力越多，因此，包含在油价中上交的燃油
4、税就随之增多，对公路养护所尽的义务也就越多。燃油税简而言之，就是将现有的养路费转换成燃油税，实行捆绑收费。由于现行的养路费一般是按吨位和运营收入两种计费方法收取，实际上形成了一种定额费。而对于用油大户尤其是汽车来说，道路使用率存在较大差距。“因为无法测算每台车的道路使用率，”长沙理工大学经济学教授朱锡庆表示，“燃油税通过将养路费捆绑进油价，将每辆汽车要交的养路费转换成税费，在道路等公共设施日益成为一种稀缺资源的大背景下，更多地体现了多用多缴，少用少缴的公平原则。”2009年6月，国家出台了鼓励新能源车消费的政策汽车产业振兴规划细则。国家将拿出200亿元补贴投向新能源汽车，支持大中城市推广混合动
5、力汽车、纯电动汽车、燃料电池汽车等新能源汽车。优先在城市公交、出租、公务、环卫、邮政、机场等领域推广使用新能源汽车。2009年至2011年，提升新能源汽车销量占乘用车销售总量至5%左右。 汽车在我国能源消耗和在大气污染排放中的分担率上占有很大比例，目前汽车每年消耗掉国产85%左右的汽油和23%左右的柴油，而汽车尾气排放带来的环境污染更是日趋严重。按照近年来的增长速度，有专家预计到2020年，我国道路上行驶的汽车数量将突破1.5亿大关。随着我国汽车保有量的迅速增加，汽车增长带来的能源和环保问题日益突出。一方面汽车进入普通家庭提高了老百姓的生活水平，但汽车大量增加带来了能源大量消耗和环境污染。在节
6、能减排大背景下，如何促进汽车产业与社会经济协调发展，促进人与大自然的和谐发展，成为我国汽车工业可持续发展的重大课题。XD479Q型汽油机结构先进、新颖、合理,升功率、升扭矩和比质量指标先进,耐久可靠,适用范围广,既适用于轻型客车,也适合于微型货车和家用轿车。本论文根据提供的数据在gt-power平台上搭建XD479Q的仿真模型，改变配气相位，模拟各条件下的动力性、经济性、排放质量等性能指标。对不同配气定时模拟出来的结果进行对比分析，选择最佳方案并设计出对应的凸轮。1.2数值模拟数值模拟也叫计算机模拟。它以电子计算机为手段，通过数值计算和图像显示的方法，达到对工程问题和物理问题乃至自然界各类问题
7、研究的目的1。在计算机上实现一个特定的计算，非常类似于履行一个物理实验。这时分析人员已跳出了数学方程的圈子来对待物理现象的发生，就像做一次物理实验。数值模拟实际上应该理解为用计算机来做实验。比如某一特定机翼的绕流，通过计算并将其计算结果在荧光屏上显示，可以看到流场的各种细节：如激波是否存在，它的位置、强度、流动的分离、表面的压力分布、受力大小及其随时间的变化等。通过上述方法，人们可以清楚地看到激波的运动、涡的生成与传播。总之数值模拟可以形象地再现流动情景，与做实验没有什么区别。从上面的例子可以看到，数值模拟包含以下几个步骤：首先要建立反映问题（工程问题、物理问题等）本质的数学模型。具体说就是要
8、建立反映问题各量之间的微分方程及相应的定解条件。这是数值模拟的出发点。没有正确完善的数学模型，数值模拟就无从谈起。牛顿型流体流动的数学模型就是著名的纳维斯托克斯方程及其相应的定解条件。数学模型建立之后，需要解决的问题是寻求高效率、高准确度的计算方法。由于人们的努力，目前巳发展了许多数值计算方法。计算方法不仅包括微分方程的离散化方法及求解方法，还包括贴体坐标的建立，边界条件的处理等。这些过去被人们忽略或回避的问题，现在受到越来越多的重视和研究。在确定了计算方法和坐标系后，就可以开始编制程序和进行计算。实践表明这一部分工作是整个工作的主体，占绝大部分时间。由于求解的问题比较复杂，比如 方程就是一个
9、非线性的十分复杂的方程，它的数值求解方法在理论上不够完善，所以需要通过实验来加以验证。正是在这个意义上讲，数值模拟又叫数值试验。在计算工作完成后，大量数据只能通过图像形象地显示出来。因此数值的图像显示也是一项十分重要的工作。目前人们已能把图作得像相片一样逼真。利用录像机或电影放映机可以显示动态过程，模拟的水平越来越高，越来越逼真。1.3模拟计算在内燃机发展中的作用大约在上世纪60年代中期以前，内燃机工作过程的计算基本上还是建立在理论循环的基础之上，对实际的工作过程作了很大的简化，例如用闭口循环来模拟开口循环；用等熵或多变指数来计算压缩和膨胀过程，用定容放热代替排气过程等。由于不考虑实际的换气过
10、程、燃烧过程、传热过程，所以只能作粗略的估算，计算结果的精确性在很大程度上依赖于经验数据的选取。因此长期以来，内燃机的设计、制造和调试主要依靠经验。直到60年代中期，由于计算机的应用才有可能对内燃机的实际工作过程进行接近实际的模拟计算，从而使内燃机的理论建立在一个新的基础之上。目前，内燃机实际工作过程的数值模拟已成为研究和开发内燃机的一种有效方法。所谓内燃机实际工作过程的数值计算（或称循环模拟计算），就是从内燃机各系统的物理模型出发，用微分方程对各系统的实际工作过程进行数学描述，然后编制计算程序，用电子计算机数值求解微分方程，以求得气缸内的热力过程、传热过程、燃气性质、气体流动、充量更换过程、
11、废气涡轮增压的特性及其与发动机的配合等问题，所以计算结果比较符合实际情况。目前，应用工作过程数值模拟大致可进行以下工作。(1) 预测发动机的性能，性能预测的内容很多，如预测发动机的标定工况性能，部分工况的运行特性，增压发动机的瞬态特性（变负荷时的响应特性），大气状况改变时的发动机性能，以及有害排放物（NOx,HC,CO）的预测等。（2）设计参数的优化 通过模拟计算可以研究各种结构参数（如增压压力、中冷器出口的空气、温度、排气背压等）以及运转参数（如发动机转速、负荷、冷却水温、活塞冷却油温等）对发动机性能的影响，并寻求最佳值。（3）在发动机调试阶段，利用实际工作过程的数值计算，可以揭示试验中没有
12、测量或不可测量的参数大小，指明参数的调整方向，从而减少试验工作量，缩短调试周期。(4)多维数值模拟能够提供空间任意位置和任意时刻的速度、压力、温度和组分浓度分布，在此基础上进行理论分析，可以为改进发动机结构提供预见性的指导。汽车工业的飞速发展，对汽车的环保、安全、经济性能等方面都提出了越来越高的要求，促使我们对发动机的研究越来越深入，涉及的范围越来越广泛。借助于数值模拟技术对发动机的工作过程进行模拟计算，对于缩短产品研发周期、提高产品设计质量、降低产品的研发成本等起到事半功倍的效果2。1.4内燃机模拟计算的现状及发展内燃机模拟计算对内燃机的研究和发展有三大贡献:一是更全面深入地理解构成模型的物
13、理过程;二是找出关键的控制变量以拟定更节省、更合理的实验方案;三是可以很方便的进行内燃机优化设计。在电子计算机用于内燃机的研究之前，为了能够从定性上对内燃机的性能做出估计，并进行有限的定量计算，不得不对内燃机的实际工作过程作出大量的简化与假设。而在内燃机设计中，则需要根据经验和类比，在大量选取参数的基础上，对工作循环中的几个特征点进行热力计算，以便做出数量上的估计。这种简单的热力计算虽然能够得到比较直观和可做定性分析的数据，但是却十分粗略，不能够全面地反映内燃机燃烧过程、缸内工质的流动及传热过程、进排气系统中的热力学和气体动力学过程以及内燃机与涡轮增压器的匹配性能等，更不能对变工况性能进行计算
14、与分析。随着内燃机性能的不断提高，产品更新的周期不断缩短，采用常规热力计算进行这种经验设计，己经远远满足不了现代高性能内燃机研制工作的需要。随着快速、大容量电子计算机和数据处理系统在内燃机研究中的应用，加上实验技术和测量仪器以及测试设备的改进，使得内燃机在实验和理论研究上有了较大的发展，同时也使得内燃机工作过程的计算机模拟不仅成为可能，而且取得了不少令人耳目一新的进展3。具体表现在:描述内燃机内部过程和与涡轮增压器匹配的一系列不受简化条件限制的微分方程组，现在己有可能进行数值求解，从而得出定量描述的结果;而那些由内燃机运行状态所决定的内部过程:如燃烧过程、气体流动过程、热交换过程以及进排气系统
15、与涡轮增压器匹配的气动过程等，也都有可能在计算中予以考虑;能够建立比较符合实际的物理模型，通过数学模拟予以表达，以求得各热力参数随时间变化的规律，用以分析内燃机的性能及其影响因素;能够系统地模拟结构参数、燃烧规律、配气相位、进排气系统中的流动阻力、中冷器特性以及涡轮增压器特性等与内燃机性能间的相互关系，并进行参数优化，寻求最佳的组合方案，为内燃机的设计、试验和性能改进提供了理论依据。这样，就使得现代内燃机的理论研究建立在了一个全新的基础上4。内燃机工作过程的模拟，不仅可以计算设计工况点，而且也可以计算变工况点，并可估计环境参数变化对内燃机性能的影响;不仅可以计算内燃机的稳态过程，而且也可以计算
16、瞬态过程，以研究各结构参数及性能参数与内燃机瞬态特性的关系，探求改善内燃机瞬态特性的技术措施;不仅可以在内燃机设计阶段通过计算模拟进行多方案的比较，且可以在内燃机调试阶段与测试技术相结合，指明调整的参数及其取值的大小;同时，在工作过程模拟的基础上，还可以进行工作过程的模拟优化，使内燃机工作过程模拟研究向优化设计推进一步。如此种种，是内燃机常规热力计算所无法比拟的5。1.5 GT-POWER软件介绍GT-POWER是包含在GT-suit中的一个CAE应用程序，GammaTechnofogies，Inc开发，具有强大的分析计算和设计功能，可以用于发动机稳态和瞬态的仿真计算，也可以用于发动机/动力系
17、统的控制分析。它适用于几乎所有形式的发动机，同时用它的部件可以搭建出用户设计的各种新型发动机，被发动机/汽车制造商和开发者广泛应用于发动机模拟和分析。GT-POWER以一维流体动力学为基础来描述发动机模型中的流动和热传递特性，它基于有限体积法，流体求解的基本观点是通过网格将系统离散为许多小的体积。在这个系统中，管子(pipe) 根据其离散长度被分割成许多小体积，而流体分叉(flow-split)只能用一个体积来描述，它们通过边界彼此相连。连续方程、动量方程和能量方程同时进行一维求解，因此流体流动方向上所有截面的值都是平均值。标量(如密度、压力等)在单元中心计算:矢量(如速度等)在单元表面(边界
18、)计算。在计算过程中充分考虑了因可然混合气的组份不同导致其热力学性质的差异，而且能用于研究一些进、排气系统结构因素(如分流、合流和弯曲等部分)对流动的影响。此外，它应用数学优化方法进行参数寻优，使得对发动机进行不需要人为经验控制的优化设计成为可能。GT-POWER包括用于建立发动机整机工作模型的前处理平台、GT-POWER解算器以及浏览和显示计算结果的后处理平台GT-POST。由于GT-POWER用的是面向对象的可视化建模方法，在系统建模的前期并不涉及到发动机各部件的几何结构尺寸，只是在利用GT-P0WER的模型库(TemPlatcLibrary)中各部件的模板(template)创建相应的对
19、象(object)和部件(Part)，井将各部件连接(Link)起来组装(assembly)成系统的图形化模型(MaP)之后，再通过输入和编辑各部件的属性(Property)确定他们的几何结构和尺寸，因而能显著减少前处理的工作量。模型库中每个模型都有自己用来描述各种输入参数的属性，许多属性值可以是随时间变化的变量，模型的大多数属性值要由用户自己定义。GT- POWER采用模块结构建立发动机工作过程的计算模型。发动机的元件(如气缸、空滤器、催化器、管接头和管道等)模块用方形图框表示，而元件之间必须用圆形图框的连接件连接。发动机的所有结构参数和特性参数在相应的元件模块和连接件模块中定义，连接件可以
20、有具体的物理定义(如气阀连接件和喷嘴连接件等)，也可以只具有象征意义(如发动机各部件及管道之间的连接件等)。建立了发动机系统模型之后，需要设置运行时间、解的收敛标准等解算器运行的控制条件。在指定了运算最长持续时间的前提下，可以将解算器运行特性设置成连续运行或者周期性运行，如果在达到最长持续时间之前数值模拟的结果已经满足了收敛标准，则模拟过程将会自行提前结束，用户定义一个最大时间步长，解算器也会根据数值稳定性标准计算出一个时间步长，实际采用的时间步长取这两个步长中的较小值。GT-POWER支持多算例，用户可以设置数值模拟的收敛标准和其他多算例中的变参数，可以选择自己感兴趣的运行结果，用图表显示出
21、来。GT-POWER还有许多有重要实践意义的特点:变量参数化、多算例模拟、局部装配模型(Subassemby)、系统优化、用户子程序、计算结果动画演示、程序中断的断点保存和恢复运行、在线帮助、多种单位支持、PC/Unix平台、模拟计算批处理。此外GT-POWER还可以和一些三维CFD软件(STAR-CD、FLUENT)或控制模拟软件包(如SIMULINK)结合起来进行祸合模拟，以使模拟的结果更准确.1.6本文主要内容及研究意义配气相位对发动机的动力性能有很大影响。如果配气相位不准确，会导致进气不充分，排气不顺畅，影响混合气的形成品质，造成燃烧不完全，使发动机的动力下降。由于废气和新鲜空气在运动
22、中的惯性，使进排气有一定的迟滞期。为了使气缸内的废气排除得更干净，新鲜空气进入更充分，需要进排气门有一个同时开启的时刻，即气门重叠角6。发动机低速时使用较小的气门重叠角，反之则使用较大的气门重叠角，这样就可以提高进排气效率，提高发动机的动力性。配气机构是内燃机的重要组成部分，它直接关系到内燃机运转的可靠性、振动和噪声，并影响内燃机的动力性和经济性等基本性能。它必须根据发动机工作循环及工作顺序的要求，按规律开启、关闭进、排气门，保证正确的配气相位。合理的配气相位是提高功率、降低油耗的一种有效方法，对改善发动机状况、节约能源、提高经济效益，降低排放有十分重要的现实意义。GT-POWER是一款包括高
23、级发动机设计概念在内的发动机仿真分析软件，可用于发动机性能模拟，以及气门升程曲线和气门正时的优化。本文利用GT-POWER，搭建了XD479Q汽油机仿真平台，然后，基于该平台计算研究配气定时对XD479Q汽油机充气效率和性能参数的影响规律，最后，重新设计了凸轮型线，研究配气定时变化对汽油机的充气效率、功率、扭矩和油耗，排放等参数的影响，从而找出最优的凸轮型线和配气定时，为实验研究提供参考。具体任务内容： 1、GT-power上搭建XD479Q发动机原机模拟平台。2、优化配气定时，对比分析配气定时对发动机充量系数和性能的影响。3、设计新的凸轮形线，将数据导入模型，模拟计算。4、分析总结，对比结果
24、，选择最佳方案。第二章 流体流动模型基本原理2.1软件的基本原理GT-POWER软件的基本思想是将发动机各系统进行模块化，将整体问题分散处理，有针对性的对具体问题进行集中建模，这样既能考虑模型的多样性，又能提高模型的精度。在对流体的处理上采用容积法将进排气管和气缸等串联起来的容积离散成小的单元，单元之间由边界相连，这样不仅能考虑流体的非稳态性，而且计算收敛快，精度高。其燃烧模型包括直喷柴油机的多区模型、基于韦伯函数放热率的汽油机模型和柴油机模型等多种型式。除了对热流体进行计算外，它还能充分地考虑曲轴、连杆、飞轮乃至于整车等机械部件对发动机性能的影响。最后集成的模型成为一个有机的整体，因而它可以
25、对发动机的非稳定工况进行动态模拟。GT-POWER软件所应用的是一维流体假设的动力学模型，它综合了发动机性能的分析代码，并几乎包含了发动机所有关键工况的细节模型，可以较完整地模拟发动机不同工况的性能变化。整个系统被划分为许多小的控制单元体，单元体上又划分许多相互交错的网格，网格是系统进行运算估值的基本单元。也就是说所有的标量（压力、温度等）都以本单元体的边缘为计算标准，系统内部划分的单元体也是如此。形状各异的气体通道都可转化为功能相当的标准管件，最后形成发动机的管网化模型，用于对发动机进行仿真。GT-POWER软件采用一维交错网络，基于一维气体动力学原理，在管路和其相关模块中以流动和传热理论为
26、基础，在特殊模块中应用液体流动和传热能力分析理论。流动模型同时包括连续方程、动量方程和能量方程。作为一维算法，设定所有的流体质量是均匀流动的。主要的变量参数是质量流量，密度和内能总和。这个系统被离散为多段网格，其中接头处被离散为一段网格，管路可被离散为一段或多段网格，这些网格通中心标量在计算中视为定值，主要的标量是密度和内能总和；辅助标量是压力，温度，焓值，组分浓度，质量分数流量等，如图2-2相关原理方程如下：连续性方程：能量方程： 焓方程： 动量方程：mflx：边界处的质量流量，是边界处有效面积，密度和速度的乘积m：单位体积的质量V：体积p：压力：密度A：流通面积e：总内能。等于内能加动能H
27、：焓值，：传热系数u：边界处的速度：表面摩擦系数：压力损失系数D：当量直径dx：边界处沿流动方向流体流过的厚度距离dp：经过dx距离的压力差离散化是将大部件分割成小总分以提高计算精度，有两种方法进行流动系统离散化。一种是将系统离散成多个管理（PIPE）和接头（FLOWSPLIT）,另一种是将长管路离散成多个短管路（离散长度应小于管路的总长度），后者相当于用多个短管代替相同长度的长管的作用。离散长度的选择规律：发动机性能分析时，进行系统离散长度大约是缸径的0.4倍，排气系统离散长度大约是缸径的0.55倍。进、排气离散长度的不同主要是由于瀑度导致声速变化引起的。噪声分析或转速超过6000rpm时，
28、离散长度应选择较小，为上述推荐值的一半左右。在一定的范围内，网格越粗，那么时间步长越大，计算速度越快，精度有一定的降低;同理网格越细，那么时间步长越细，计算速度越慢，精度有一定的提高。所以，在仿真中应该尽量使用网格的精细趋于均匀，这有利于提高计算精度和减少计算时间。2.1.1 管路（PIPE）模块基本原理管路（PIPE）模块设定为截面为圆形或方形的直管路，可通过设定不同的摩擦系数，传热系数，压力损失系数实现模拟不同几何形状的管路。对于锥度较小的圆锥管可采用改变压力损失系数的方法用直管代替，对于锥度大于30度或圆锥管长度小于离散长度的锥管需用直径突变的两段管路实现孔（ORIFICE）的结构，代替
29、圆锥管。管路的流动损失主要源自沿壁的摩擦损失，计算时考虑雷诺数和表面粗糙度。当管壁光滑时： （层流区，4000）当管壁粗糙且不是层流时:：雷诺数D：管路直径h：粗糙度管路的压力损失取决于其锥度，弯曲度，不规则接头等。而管径突变带来的压力损失不能单纯靠设定压力损失系数，需要考虑孔（ORIFICE）和膨胀直径（EXPANSION DIAMETER）7。2.1.2 接头（FLOWSPLIT）模块基本原理虽然GT-POWER是一维计算软件，但是接头模块（FLOWSPLIT）的设计主要考虑三维因素对流动的影响。在管路计算中，标量（如质量，能量等）位于有限容积的中心，而矢量（如质量流量等）位于边界处。若管
30、路（PIPE）中存在开口（PORT），就变成了接头（FLOWSPLIT）。接头（FLOWSPLIT）的计算原理与管路（PIPE）类似，标量仍位于有限容积的中心，而矢量方程在每个开口（PORT）处独立计算。接头几何形状由膨胀直径（EXPANSION DIAMETER流体进入管接头后能够达到的直径，主要用于计算流量系数和能量损失），特征长度（CHARACTERISTIC LENGTH流体进入管接头达到管壁或边界的长度，主要用于计算压力波的传递和反射）和管接头各开口的空间角度这三个因素决定。单纯的一维算法无法精确计算流体流出任意方向开口处的动量，接头（FLOWSPLIT）模块则是专门计算此种情况的，
31、流出出口的流体动量由按其出口方向特征速度矢量计算得出。图2-1为简单“三通”接头属性编辑界面，“1、2、3”代表三通的三个开口（PORT），膨胀直径（EXPANSION DIAMETER）如图中紫色标注所示，特征长度（CHARACTERISTIC LENGTH）如图中黑色标注所示，管接头各开口的空间角度为开口中心线与X、Y、Z轴的夹角，接头的容积可手动输入也可由GT-POWER软件的Discretizer功能自动生成。 图2-1摩擦，传热，壁温等计算原理与管路模块（PIPE）相同，不同处在于摩擦系数和传热系数计算中的速度值为接头开口处的速度。当接头（FLOWSPLIT）的几何形状不规则时，要精
32、确确定其容积，通过适当调节摩擦系数、传热系数和压力损失系数来实现不规则形状对流体流动带来的影响。2.2缸内流动仿真原理进、排气道建模原理与管路（PIPE）基本相同。气门处流量系数计算来自于与定压差质量流量的测量，这种测量通常是在气门安装在发动机以前进行的。气道的不规则形状（弯曲、直径变化），表面粗糙等均造成压力损失。气门流量系数的计算已经包含了气道压力损失的影响。为了防止计算过程中气道的损失被计算两次，因此在气道设定中摩擦系数和压力下降系数应设为“0”。气道的温度设定主要是由于发动机冷却系统决定的，流体流动对壁面温度的影响较小。此时的传热系数设定值高于正常值“1”，这主要用来体现气道弯度，喷油
33、嘴和气门杆的多余面积，气门和气门杆形成的紊流的影响。ENGCYLWALL模块推荐设定进气道温度450K，排气道温度550K，高温主要是来自于气门的传热。气缸内被分割成很多部分，中间区、挤流区、头部区和活塞凹槽区，在每个时间步长内，每个区域的轴向速度、径向速度和涡流速度都被应用到燃烧室几何形状，活塞位置，进排气门处流量/涡流/滚流等的计算。2.3燃油喷射模块基本原理在燃油喷射模块中，用“燃油雾化率”来定义燃油经喷嘴喷射瞬时雾化的比例。由于蒸发吸收热量改变了燃油的焓值，导致空气、雾化燃油和液体燃油三者的混合物冷却。此外雾化的燃油在气道内要另外多占一部分体积。进气道喷射模型设定燃油喷入气道后，没有雾
34、化的部分将一直以液体的形态存在，直到燃烧开始。燃烧计算过程中，雾化和液体的燃油始终按瞬态固定的比例燃烧。对于典型的进气道喷射模块，燃油雾化率初设为0.3（即燃油喷射瞬间30%的燃油气化）。燃油雾化蒸发体积膨胀对充气效率的影响不显著，主要因为雾化吸热对混合气的冷却抵消了液体蒸发膨胀的体积。这种特性对于气化潜热很大的醇类燃料更为重要，喷口布置的离进气门越远，雾化的时间越长。对于顺序喷射的汽油机，需要知道燃油喷射率，空燃比和喷油脉宽。2.4燃烧模型基本原理虽然点燃式发动机和压燃式发动机燃烧模型建模原理均利用缸压计算放热率预测发动机性能8，但由于点火方式不同给燃烧带来影响，GT-POWER软针对点燃式
35、发动机和压燃式发动机分别设定缸内燃烧分析模块。本论文仅应用“EngCylCombSIWiebe”单维韦伯燃烧模块对点燃式发动机进行整机性能模拟计算。主韦伯常数：燃烧始点： BMC =ln(1-BM) 燃烧中点常数，燃烧率设定为50%BSC =ln(1-BS) 燃烧始点常数，燃烧率设定为10%BEC =ln(1-BE) 燃烧终点常数，燃烧率设定为90%AA：上止点到韦伯曲线燃烧50%之间的曲轴转角D：韦伯曲线上燃烧10%90%之间的曲轴转角燃烧率计算：瞬时曲轴转角92.5发动机摩擦及驱动附件损失模块基本原理GT-POWER软件用“EngFrictionCF”模块设定发动机的摩擦损失功，驱动附件（
36、如水泵、散热器风扇等）的损失功也可连接到此模块里。计算公式如下：C：摩擦损失功的常数部分(Recommended values=0.3-0.5bar)PF：缸内最大压力系数(Recommended values=0.004-0.0.006)：缸内最大压力MPSF：活塞平均速度系数(Recommendedvalues=0.08-0.10bar/(m/sec)MPSSF：活塞平均速度平方系数(Recommended values=0.0006-0.0012bar/(m/sec)2)：活塞平均速度此公式是在未知摩擦损失功情况下预测整个发动机摩擦损失功的经验公式，摩擦损失功与缸内最大压力和活塞平均速度
37、有关。第三章 XD479Q发动机模拟平台搭建3.1 XD479Q发动机参数XD479Q型汽油机主要用作轻型客车、轻型载货车的前驱动力，也可以作为同类型汽车和其它机械配套动力使用。浙江吉利汽车公司的JL479Q具体相同的设计要求及研发目的，因此建模过程中XD479Q主要性能指标和结构参数可参照JL479Q设置，具体数据如表3-1所示：表3-1XD479Q型汽油机主要参数10冷却方式水冷缸数和布置直列4缸冲程数4燃烧室型式对向斜顶面气门机构4气门、顶置、双凸轮轴缸套型式无缸径行程78.969总排量，L1.349压缩比9.35最大输出功率,kW/rmin-163/6000最大扭矩,Nm/rmin-1
38、110/5200全负荷最低油耗率,g/(kWh)279怠速稳定转速,r/min700润滑油容量,L3.5进气门开,CA16(上止点前)进气门关,CA56(下止点后)排气门开,CA48(下止点前)排气门关,CA17(上止点后)气门间隙(冷态),进气门,mm0.20气门间隙(冷态),排气门,mm0.30点火提前角,CA上止点前10(700-800r/min)外形尺寸,长宽高,mm668592648干质量(不带离合器),kg1033.2整机模型的搭建 几何结构的精确描述是准确实现数值模拟的重要前提，但汽油机模型涉及到的实体很多，而且结构极其复杂，如按真实实体结构建立计算几何模型需要花费大量时间且难以
39、实现，也没有必要。因此根据计算经验，在保证对数值计算精度不产生很大影响的前提条件下，对实体结构进行一些等效简化处理。由于对整机进行分析涉及到的工作过程太多，同时也受到计算速度的影响，因此将进排气管路简化成一维的模型，将缸内工作过程简化成零维和准维模型，这样既能保证计算的精度，又可以节省大量的时间。流程图如图3.1所示，模型图如图3.2所示。3.2 XD479Q GT-POWER模型示意图3.3定义外部环境参数这个模块主要确定压强、温度、气体成份等描述外部边界环境状况的参数，参数如图3.3所示: 图 3.3指针变量air的设定如图3.4所示：图3.43.4定义空气滤清器 创建空气滤清器管件，为了
40、计算方便，将空气滤清器模拟成一维管道，通过设置压力损失系数实现它原有得功能。参数设置如下： 图3.53.5 定义节气门模板“EndEnvironment”用来描述上游边界条件（upstream boundary condition）。将由一个进气管来连接environment 到空气滤清器，一个管件（pipe）连接空气滤清器到进气歧管，歧管与空气滤清器之间连接了一个节气门。参数设置如下：图3.6图3.73.6定义连接管件及孔实体 在模型的离散化过程中，经常需要一些管件及孔实体作为连接件。只要溢流口和管件连接到一起，就会自动创建一个“def”孔。缺省孔（defaultorifice）定义模块实体
41、之间的连接时，GT-power 将计算前向和返回的分散率（forward and reverse dischargecoefficients）。尽管如此，在多数歧管中，从歧管到管道的转变都是非常平滑的（smooth）。使用缺省孔将为系统增加一个并不真实的额外的损失。因此，需要使用未缺省（non-default）孔，这里，分散率将是1。 从环境未端到进气歧管未端之间的连接管件均按如下参数进行设置：图3.8孔实体设置参数如下：图3.93.7进气系统的创建及定义 建立进气系统，首先要建立进气歧管。绘制进气歧管以及离散歧管的方法（the method used to discretizethe int
42、ake manifold）如下图所示。图3.10图3.11传统的分析中，离散化长度（discretization length）规定：进气系统中，一般是孔径的4 倍；排气系统中，是孔径的5.5倍。离散歧管的一般步骤如下：1、确定目标离散长度；2、从歧管到管道的溢流口的尺寸和位置。溢流口的尺寸和位置必须首先定义，因为它们通常包含一个容积，并且不能再次做离散。3、利用总体尺寸（overall dimensions）来计算管到溢流口的长度。利用进气系统管件的目标离散化长度。如果某些管件太短（小于5倍离散化长度），它们应该与相邻的合成一体离散。上面提到的歧管几何形状已经简化了，因此只需一个溢流口实体（
43、flowsplitTRight object）和一个管件实体（pipe object）来完成歧管模型。下面列出进气管道中溢流口实体和管件实体的参数设置：图3.12图3.133.8定义进、排气门XD479Q汽油机采用四气门结构，2 个进气门和进气道,2 个排气门与排气道均成错流分别布置在气缸盖两侧，排气凸轮轴通过齿型皮带由曲轴正时齿轮驱动，进气凸轮轴借助一对齿轮与排气凸轮轴同步反向旋转。进气或排气系统结构相同,其运动规律各自独立地受制于凸轮型线,同一缸两个进气门或两个排气门的关闭规律自然保持同步。进气门气门升程曲线影响到发动机曲轴每个角度处的气门开启程度，进而影响到进气压力特性的仿真结果。GT-
44、POWER 与之相关的边界条件包括了ValveCamPRConn模块的Lift Arrays、Angle Multiplier 和Cam Timing Angle。Lift Arrays 是气门升程lift array 与曲轴转角theta array 的对应关系，其输入采取theta array 为0 时，lift array 为最大值的形式。Angle Multiplier 是对气门从开启到关闭的角度差进行调节的参数，通过对theta array 乘上Angle Multiplier可以将Lift Arrays 中输入的气门从开启到关闭的角度差调节至与真实情况相同的值。Cam Timing
45、 Angle 是气门升程达到最大值点（即theta array 取0 点）与压缩上止点间的曲轴角度差。由于无法获知测试车辆的真实气门升程曲线，而只清楚其气门升程最大值为8.163mm，为了简化问题，Lift Arrays 中theta array 使用GT-POWER 设置的默认值，但对每点的lift array 乘上8.163mm 与默认lift array 最大值10.204mm 的比值0.8，得到的theta array 与lift array 对应关系如附录表1所示。Lift Arrays 中填入的气门从开启到关闭的角度差（即theta array 值的跨度）为156 度曲轴转角，而测
46、试车辆的角度差为252 度。按照比例关系，调整后的气门开启点与升程最大点间的曲轴角度差为=117.9本机的气门开启提前角为16 度曲轴转角，与压缩上止点成344 度曲轴转角轴角度差，因此Cam Timing Angle 值取344 度与117.9度之和，即461.9度，约为462度。取值时取凸轮转角，因此为231。调整Lift Arrays、Angle Multiplier 和Cam Timing Angle 值后得到的GT-POWER 输出进气门升程曲线如图3.14所示图3.14排气门在计算cam timing Angle的时候用的是凸轮转角，即为曲轴转角的一半，通过上述方法计算结果为257
47、度曲轴转角，取一半值即为124度。调整参数后排气门的升程曲线如图3.15所示。图3.15进气门的具体参数设定如下：图3.16图3.17图3.18排气门的具体参数设定如下：图3.19图3.20图 3.21图3.223.6气缸参数设定气缸模型建立得过程中，需要确定的是燃烧室的具体结构尺寸；确定气缸壁温度；热量在气缸内传递的方式及燃烧模型等。根据原机具体参数和GT-power推荐参数值建立燃烧室模型如下：图3.23图3.24图3.25 指针cylg1的设定本论文应用“EngCylCombSIWiebe”单维韦伯燃烧模型对点燃式发动机进行整机性能模拟计算。主要需要输入得数据是由燃烧0%50%，以及10%90%所对应的曲轴转角，软件自动调用计算程序，计算所需要的点火提前角，再调入燃烧模型计算公
