基于虚拟样机技术的摩托车发动机配气机构仿真分析

1、虚拟样机技术及其在摩托车发动机配气机构仿真分析中的应用邹志民前吞业内专家断言：我国摩托车发动机制造水平至少落后日本10年，设计开发能力落后15 年。因此，当前及今后的一段时间还需更仿制和引进国外摩托车发动机产品技术，而虚拟样 机技术为我们消化和吸收国外技术，追赶先进水平提供了一条捷径。以前我们仿制和引进别 人的东酋，一般是原封不动的照抄别人的零件，只知其然而不知其所以然，缺乏系统层面上 的理解，所以仿制的产品总赶不上原机的水平，这也就是我们的设计质量始终赶不上别人的 根木原因。采用虚拟样机技术，技术人员便可以对引进的样机进行深入研究，追踪样机的设 计思想和设计原理，甚至还可以发现样机的缺陷进行

2、优化改进，“青出于蓝而胜于兰”。国外 产品的高水平的设计源于多年的经验积累，而借助于在虚拟样机上的模拟，我们也会很快取 得一些经验成杲，及时用到新产品自主开发当中。c100发动机是国内弯梁车普遍使用的发动机，本文尝试采用虚拟样机技术对该型发动机 的配气机构进行了仿真分析，各种分析结果不单可以为优化指明方向，同时也有助于从广义 的系统层面上理解该机的设计理念，积累技术素材。一、虚拟样机技术(virtual prototyping technology)1、概述虚拟样机技术是指在计算机上建立一个虚似化的完整的产品模型，这个虚拟的模型具有 物理样机的功能，可以代替物理样机进行试验分析。对于传统产品开

3、发来说，产品的性能只 有在物理样机制造出來并进行试验后才能得到，如果性能达不到要求，只有修改设计，再做 样机进行试验，这样反复多次，周期很长，费用高昂，是制约产品开发的瓶颈。而虚拟样机 技术的出现，在产品的设计阶段就可以对性能仿真优化。从而大大节约了物理样机的制造和 试验的时间和费用，虚拟样机技术使产品开发真止产生了质的飞跃。运用虚拟样机技术，产品设计人员可以在齐种虚拟环境屮真实地模拟产品整体的运动及 受力情况，快速分析多种设计方莱，进行对物理样机而言难以进行或根本无法进行的试验， 直到获得系统级的优化设计方案。虚拟样机技术的应用贯串在整个设计过程当屮，它可以用 在概念设计和方案论证中，设计师

4、可以把自c的经验与想像结合在计算机内的虚拟样机里， 止想象力和创造力充分发挥。当虚拟样机用来代替物理样机验证设计时，不但可以缩短开发 周期，而且设计质量和效率得到了提高。2.虚拟样机技术的实现虚拟样机技术以三维实体建模（3d-cad： 3d-computer aided design） >多体系统动力学 分析（msda： multi-body system dynamic analysis）和有限元分析（fea： finite element analysis）为主线，以3d-cad、msda和fea软件为支撑。3d-cad软件提供系统各零部件的儿何信息，通过定义零部件的材料属性，cad

5、软件可以 自动计算出零部件的物理信息（质量、转动惯量等），进一步可以在计算机上定义零部件间的 连接关系并对机械系统进行虚拟装配，从而获得机械系统的虚拟样机（实体模型）。其代表性 软件有：pro/e> ug、solidworks mdt> inventor 等。msda软件在“实体模型”基础上进一步定义运动副、接触、速度、加速度、力、阻尼等 动力学参数形成“动力学分析模型”，仿真机械系统在各种工况下的运动和受力情况，通过分 析仿真结果，可以在计算机上很方便地发现设计缺陷加以改进，并立即验证改进效果。通过 不同设计方案的仿真分析比对，可以最终获得最优设计方案。其代表性软件有：msc.

6、adams、 dads、simpack、cosmosmotion 等。fea软件可以在零部件“实体模型”的基础上通过添加边界条件后进行静力学分析（分析 结构在静力作用下的应力与应变）和动力学分析（分析结构的固有频率、主振型（模态）、以及 在动载荷作用下的动挠度、动应力等）。其代表性软件有：ansys、msc/nastran. cosmosworksmsda软件的分析出的零部件的受力结果，为fea分析提供了准确的边界条件，而通过fea 分析，又可以提高和优化零件设计，从而实现产品从局部到整体、从单个零件到系统都得到 提高，达到了真正意义的虚拟样机的功能。-x虚拟样机技术在配气机构仿真分析中的应用

7、1、c100发动机配气机构的零部件构成1、挡油罩组件2、凸轮轴组件3、摇臂4、摇臂轴5、进气门6、排气门7、气门外弹簧8、气门内弹簧9、弹簧座10、弹簧圈11、锁夹12、调整螺钉13、调整螺母图1 c100发动机配气机构的零部件构成2、实体建模采用3d-cad软件严格按照零部件图纸精确建立了配气机构中的：汽缸头本体（仅完成了 分析相关部分）、凸轮轴、进气门（含气门锁夹和气门弹簧座）、排气门（含气门锁夹和气门弹 簧座）、气门摇僭、摇臂轴、气门调整螺钉、气门调整螺母等零部件的实体模型。由于进行仿 真分析时气门内外弹簧的弹力可由相应参数直接给出，因此不需要建立实体模型。建模时“凸 轮型线”采用“样条

8、曲线”的来描述。零部件实体建模完成后，再在建模软件屮选择零部件材料或输入材料的性质参数，软件 就可以口动计算出零部件的质量、转动惯量和质心位置等物理特征。由于是按图纸建立的模 型，因此，考察软件计算岀的零部件物理参数与图纸标注的物理参数或实际测得的物理参数 的一致性，是对建模正确性、准确性的一种检验方式。准确的实体模型是保证后续仿真分析 结果准确可靠的基础和前提。部分零部件的“模型质量”和“实际质量”统计表见表1：表1零部件质量参数对照表序号零部件名称模型质量实际质量备注1排气门0. 024kg0. 023kg均含：气门锁夹和气2进气门0.026kg0. 025kg门弹簧座3气门调整螺钉0.

9、0025kg0. 002kg4气门调整螺母0. 0016kg0. 002kg5气门摇臂0.041kg0. 045kg图2 c100发动机配气机构的实体装配(气门弹簧未建模)零部件实体建模完成并确认无误后，就可进行实体装配建模。图2所示为严格按照产品 图纸完成的实体装配。考虑到摩托车发动机的气门间隙只有0. 05mm,因此建立虚拟装配体时 忽略不计按无间隙处理了，即气门调整螺钉的球面与气门杆的端面直接采用了 “相切”约朿。装配体完成后，可以利用3d-cad软件提供的对测量工具对一些关键部位进行测量，用以 验证模型的正确性，同时也可以获得一些综合信息。例如：现代内燃机气门弹簧都采用内外 双弹簧，每

10、根弹簧都采用变节距结构，以充分利用空间，在有效保证弹力的同时减少共振。 一般的设计准则是：弹簧“力值一压缩量”曲线图上的拐点(p1)就在配气机构“自由”状态 时气门弹簧的长度处，p2点设计在气门最大升程处(取进排气门最大升程中哪一个大按哪一 个的数据)。c100发动机配气机构气门弹簧参数见表2。实体建模后，从模型上测量得到的 p1点处装配体上相关数据值和图纸要求值对比见表3。表2气门弹赞参数序号项目气门外弹簧气门内弹簧1弹簧自由长度(mm)35. 2532.412p1点长度（mm）25.4522.453p2点长度（mm）19. 6516. 654pl和p2点的长度差（mm）5.85.8表3 p

11、l点气门弹簧长度pl点要求气门外弹簧气门内弹簧装配体实测图纸耍求装配体实测图纸要求进气门25. 4325. 4522.4322.45排气门25. 4325. 4522.4322. 454. 动力学分析在“实体模型”的基础上定义后续仿真分析过程所需要的各种参数：部件（如：固定件、 运动件）、约束（如:平动副、转动副、凸轮副等运动副、移动、转动的相对运动）、力（如： 弹簧力、阻尼力、冲击力等）、以及必要的参考坐标系等就建立了 “分析模型”。'分析模型” 创建完成后就可以利用分析软件屮的仿真分析命令进行各种仿真分析，从而获得各种结果。 而且可以很方便改变各种初始条件，进行多参数或多方案比较，

12、这也是“虚拟样机技术”的 主要优点。配气机构的“凸轮轴”与“摇臂”间约朿为：标准凸轮副中的曲线一曲线碰撞约束，并 且依据凸轮轴与摇臂在发动机运转中实际可能的配合情况还选取了允许“间歇碰撞”的选项 （可以仿真凸轮和摇僭间的“飞脱”情况）。由于“气门调整螺钉”与“进气门”（或“排气门”）间的实物接触状态为“球面对平面” 接触，因此定义为：“3d碰撞接触”约束。同时将上面的“曲线-曲线碰撞约朿”和“3d碰撞接触约朿”间两实体的丁作状态由软 件默认的“干磨擦”状态，调整为“润滑”状态。“气门调整螺钉”以及“气门调整螺母”实物状态为固定在“摇臂” ±,随摇臂一起运 动，因此两件都用“附着在”约

13、束将他们固定在摇臂上。气门内弹簧和外弹簧可以通过“力库”中的“线性弹簧”模板生成，定义了作用位置之 后，再输入相应的“刚度”、“长度”、“预紧力”等参数就可以在配气机构的“分析模型”中 得到具有特定物理参数的虚拟弹簧。“凸轮轴”作为原动件，在其上还定义了一个恒定转速的运动。该转速为曲轴转速的一半。改变凸轮轴转速就可以仿真不同发动机转速下配气机构的运转情况。分析持续时间按保证凸轮轴正好转动1周（360°）确定。为保证有足够的仿真精度，同时 方便后续分析计算，360°圆周按1000帧处理，即不到半度就计算一次。再将其它的参数全部定义完成z后就可以进行仿真分析，仿真分析过程由软件

14、智能完成, 软件第一步首先计算整个系统的自由度，如果自由度为零，进行运动学分析。如果自由度不 为零，则软件通过分析初始条件，自动判定是进行动力学分析还是进行静力学分析。然后进 行仿真分析，分析完成后按要求输出参数的数据曲线、表格或生成实时动画。进行发动机各种不同转速下的仿真，可得到各种转速下的气门升程、速度、加速度、摇 臂和气门间碰撞力、摇臂与凸轮的接触力、气门落座碰撞力等。以及一些不常用的参数：凸 轮压力角，凸轮相对轮心滑动位移，凸轮和摇臂相对滑动速度、加速度，摇臂相对接触曲线 圆心滑动位移，摇臂、气门相对滑动位移、速度、加速度等。发动机转速8000rad/min（凸轮轴运转一周时间为：0.

15、015s）时的分析结果：4. 1、气门升程和配气相位o6543210o12o24003angle(deg)图3 c100发动机一个工作循环内进排气门升程曲线（粗线为排气门升程、细线为进气门升程）表4气门升程分析计算表序 号项目进气门排气门仿真结果设计值仿真结果设计值1最大气门升程（imn）5. 7795.85. 5855.621mm气门升程时刻207.497207211.655213配气相位(°)3丰满系数w0. 72324.2、接触力和冲击力图4进气门惯性力(图中粗线)和气门弹簧力(图中细线)(kopio 工00.0050.01time (see)o o o o o o o o o

16、 o o 0 8 6 4 20-015图5排气门惯性力和气门弹簧力的合力表5气门惯性力与气门弹簧力的合力序号项目进气门排气门备注1气门惯性力与气门弹簧力的合力最人值(n)967.81871.53仿真得到2气门惯性力与气门弹簧力的合力最小值(n)93.4842.87仿真得到由表5矢fh气门惯性力与气门弹簧力的合力在整个运动过程中均大于0,因此可以保证气门在运动过程屮“气门”与“气门调整螺钉”之间不发生“飞脱”。表6气门摇臂与凸轮的接触情况序号项日进气门排气门备注1作用力开始为零(飞脱)的时刻0.0103800.007200仿貞.得到2飞脱后又开始接触的时刻0.0108000.0072453又开始

17、接触瞬间的冲击力(n)34542039. 54飞脱持续时间(s)0. 004200.000045由仿真 数据计算 得到5e脱持续时间占整周运动时间的比 例2.8%0. 3%仿真分析表明改型发动机在8000rad/min时开始出现“气门摇臂”与“凸轮”飞脱情况,但持续的时间非常短。仅仅从计算机屏幕上显示的仿真动画上无法看到“飞脱”。当发动机到 达10000rad/min时从计算机屏幕就可以看到非常明显的“飞脱”现象。5、有限元分析各种有限元分析软件的使用流程基本相同，简单來说可分成三个阶段，前处理、处理和 后处理。前处理是建立有限元分析模型，完成单元网格划分；后处理则是采集处理分析结果, 使用户

18、能简便提取信息，了解计算结果。顶置凸轮配气机构中“气门摇臂”的工作条件、运动情况和受力状况都较为复朵。运动 过程中的角速度和角加速度都很高。它是一个双臂杠杆，两臂都承受较人弯曲应力，摇臂的 结构设计要求在最轻的重量下有最好的强度、刚度和固有频率，因此大多数采iht字形或工 字形的断面。5. 1固有频率分析节点(nodes)数为：81480、单元(elements)数为:54342、自由度(dof)为:240648表7气门摇臂模态分析结果模态号频率(赫兹)周期(秒)16538. 20.0001529526931. 10. 000144283112760.000088694124470.00008

19、0345210650.00004747由频率分析数据可以看出：气门摇臂的固有频率为6538.2hz,远远高于凸轮轴的转速范 围 05500rad/min(092hz)。5.2、静力分析结果静力分析的t况取为发动机转速10000rad/min和8000rad/min两种情况。由动力学分析可以确定10000rad/min和8000rad/min两种情况下气门摇臂受力的极限情 况，数据汇总如表8：表8气门摇臂受力发动机转速(rad/min)气门摇臂在气门侧受力最大值(n)气门摇臂在凸轮侧受力最大 值(n)100008719148888000967.813454将上述载荷加载在气门摇臂相应受力面上，并在气门摇臂上与摇臂轴配合表面上添加“不可 移动”的约束条件，同时确认材料认为“合金钢”、完成网格划分，就可以进行“静力分析”, 分析后可以得到“应力”、“应变”等结果：表9气门摇臂静力分析结果发动机 转速 (rad/min)气门摇臂上的 最大应力(n/m2)气门摇臂上的 最