直列四缸发动机曲轴ansys分析报告

有限元分析课程报告直列四缸发动机曲轴有限元分析学号:分数：年 月 日目录1. 引言12. 曲轴有限元模型的建立23. 曲轴网格划分43.1确定物理场和网格划分法43.2确定全局网格参数设置43.3确定局部网格参数设置53.4网格质量检查53.5调整网格划分64. 曲轴强度分析94.1载荷工况94.2载荷与约束边界条件94.3有限元求解结果分析114.4加密收敛分析135. 曲轴自由模态分析165.1网格划分165.2约束条件165.3参数设置165.4结果分析17参考文献193/221. 引言曲轴（如图1所示）是发动机中最重要、承载最复杂的零件之一。其造价约占一台发动 机总价的20%-30%，它的使用寿命决定发动机的使用寿命。在工作过程中，曲轴承受着气缸 的气体压力与往复和旋转质量惯性力引起的周期性变化载荷。这些扭转力矩和大小、方向不 断变化的弯曲力，在曲轴产生交变的弯曲与切应力，使曲轴发生疲劳破坏。同时，在周期性 变化的载荷作用下，曲轴可能在发动机转速围发生共振，加速曲轴的疲劳破坏。由于曲轴是发动机中典型的易损件之一，其强度和振动特性都会影响到整机的工作性能， 因此对曲轴进行有限元分析，研究曲轴的应力和应变状态，了解其固有频率和模态振型，对 于设计和优化曲轴结构有重要意义。图1 曲轴实物图曲轴的基本结构如图2所示，主要由前端轴、主轴颈、连杆轴颈、曲柄、平衡重和后凸 缘盘等部分构成，前端轴主要用来安装止推垫圈、启动抓和皮带轮等部件；曲柄用来连接主 轴颈和连杆轴颈；平衡重的作用是平衡曲轴工作时产生的离心惯性力与其力矩；后凸缘盘用 来安装飞轮等部件。图2曲轴基本结构1-前端轴2-主轴颈3-平衡重4-曲柄5-连杆轴颈6-后凸缘盘2. 曲轴有限元模型的建立本文所研究曲轴为某直列四缸发动机整体式曲轴，含有5个主轴颈，4个连杆轴颈、曲 柄和平衡重。由于曲轴的具体尺寸参数设计较为复杂，本文不详细叙述。同时，由于曲轴结 构复杂，利用有限元软件进行建模时很难保证与图纸上的曲轴结构完全一致，因此建模时必 须简化模型。为了减少应力集中，曲轴上不同截面的结合处都有半径不同的倒角，如果在建 模时考虑这些倒角和油孔，则会使有限元的网格非常密集，大大增加了模型的单元数量且生 成的网格形状也不理想，降低了求解精度。因此，建模时忽略了对分析结果影响较小的特征， 如轴颈上的油孔和后凸缘盘上的环槽。现给出曲轴结构的主要设计尺寸如表1。表1曲轴结构设计尺寸设计参数主轴颈长度主轴颈直径连杆轴颈长度连杆轴颈直径曲柄厚度平衡重厚度尺寸/mm358040662522.5平衡重的二维设计尺寸如图3所示。图3 平衡重设计尺寸考虑到在ANSYS中直接进行曲轴建模的复杂性，本文通过SolidWorks建立三维几何模 型，并与ANSYS Workbench无缝连接，进行有限元分析。曲轴三维建模后如图4所示。图4曲轴三维实体图本文中曲轴材料为QT-800，其基体组织为珠光体或索氏体，主要用于所受载荷较大、 受力复杂的汽车、拖拉机的曲轴、连杆、凸轮轴等等场合。表2列出了 QT-800的材料性能 参数。表2QT-800材料参数材料抗拉强度/MPa抗拉屈服强度/MPa硬度/HBS弹性模量/密度/泊松比QT-800N800N480245-35578500.3在workbench中设置曲轴材料参数如图5所示。i-ropenjei-a1 Oudne fcw 2：&C1PlCiptrtiUrit7.65E-O9toine3a 函fwopc 强Ey4DerwerrmiVDings IModLius and PotsscnE Rath二5Voungs Modulus2E*05MPa6Poissons Rato0.37Bilk Mockikis.6667E05MRaeShear ModulusMPa9S Tensie Ytetl Strengtti湍MPaIDTensfe UtiTiate b-angthMPa图5曲轴材料参数1.3. 曲轴网格划分有限元网格划分是将整体模型结构离散化，是数值分析的前提，也是至关重要的一步， 它直接影响着后续数值计算分析结果的精确性。一个良好的网格划分，不但可以保证分析结 果的精度，同时也能够大大降低后期处理的计算量。本文中采用ANSYS Workbench15.0中自带的Mesh模块对曲轴进行网格划分。ANSYS Workbench15.0是一个单独的工具平台，为ANSYS的不同求解器提供相应的网格文件。自14.0 起，Workbench中可以直接划分网格(Direct Meshing)。直接划分网格最大的优点之一就 是能单独地划分几何体的网格，而之前划分网格时只能整个模型同时一起划分，显然，直接 控制网格划分具有更大的柔性。3.1确定物理场和网格划分法在网格划分之前必须确定物理场的类型，Workbench中提供Mechanical CFD、Electron-magnetics和Explicit四种物理场类型。由于本文对曲轴结构进行静力学和模态 分析，因此选择Mechanical网格划分法求解静力学结构场。3.2确定全局网格参数设置全局网格设置通常用于整体网格划分的部署，包括网格尺寸函数 size，边界层函数 Inflation，平滑度函数Smooth等等。现设置全局网格具体参数见表3,如图6所示。表3全局网格参数RelevanceElementSmoothingTransitionSpan AngleCenterSizeCenterCoarse Default MediumFastCoarseDetails of Mesh*Physics PreferenceMechankaCRelevance0Eglhe Mvaiced Size- FuncfcnOffRelevance CenterCoarseElement SizeDefaultSize SeedAdcve A繇刖ibiySmcxjchgMediumTranrsitionFastSpan Angle CenterCoar&eMinimum Edge Length4383如 mm图6全局网格参数设置界面3.3确定局部网格参数设置整个曲轴由5个主轴颈，4个连杆轴颈，4个曲柄和平衡重，1个前端轴和后凸缘盘。 为了整体较好的网格划分结果，需要对整个模型进行分割，不同部分进行不同方法划分网格。曲轴轴颈部分为一系列规则的圆柱体，采用Sweep Method划分法；曲柄、平衡重、前 端轴和后凸缘盘形状较为不规则，采用Automatic Method划分法，即在四面体和扫掠型划 分之间自动切换，当几何体不规则时，程序自动产生四面体。由于曲柄、平衡重、后凸缘盘（含螺纹孔）和前端轴（含键槽）形状较为复杂，对这几 处采用Body Sizing进行局部加密。整体曲轴网格划分后如图7所示，共包含45121个单元，82631个单元节点。图7曲轴网格划分3.4网格质量检查本文网格质量检查，采用Skewness （偏度）进行度量Skewness基于归一化的角误差， 适用于所有的面和单元形状，也适用于所有棱柱和棱锥，是最基本的网格质量检查项，其值 位于0与1之间，0最好，1最差。度量等级详见表4。5/22表4 Skewness网格质量度量等级0-0.250.25-0.500.50-0.800.80-0.950.95-0.980.98-1.00极好很好好可接受坏不可接受网格划分后，网格质量如表5所示。网格划分质量分布如图8所示，整体网格质量较好， 小于0.95的单元个数为42775，不可接受单元个数为3，占总体0.07%。，但可以发现网格质 量较差处多位于连杆轴颈处，因此后续分析中可能会出现应力奇异现象或者网格加密不收敛。 因此对于此网格划分，需要进一步调整，使其连接处网格质量提高。表5网格质量度量方法最小值最大值平均值标准差Skewness0.00160.9980.2900.175BL 1 k L li l. L iMtU*心心M代写FtaIm图8网格质量分布图3.5调整网格划分重新确定全局网格参数设置，采用高级尺寸函数(Advanced Sizing Function),此函数 主要控制曲面在曲率较大的地方的网格。其中Proximity and Curvature函数采用曲率法确 定细化边和曲面处的网格大小，同时又控制模型邻近区网格生成，对于窄薄处网格也有较好 的划分效果。更改后全局网格具体参数如图9所示。Details of1 Meshthe Advanced Size FuflOn: Projclmrly and CurxatureRelevance CenterCoarseInitial Size SeedActive As-semblySnnoo thingMfdiiLimTransitionRistSpan Angiie cenierCoarseCurvTure Normal An.,。河aim 口。.豹 5。)Munn Ceils Acros GapDefault (MMln SizeDefault 0.341060 mm)Proximity Min SlzsDefault 0.341060 mm)Face &lzeDefault (S4.1IK0 mmiMaocSiEeDefault (68.2120 mmjGrowth RateDefault (1.891)Minimum Edge Lenglh4弗削nrn图9 更改后全局网格参数由于先前网格划分在轴颈处网格质量较差，于是推测为sweep方法下六面体网格与四面 体网格过渡过程中，出现较差质量的网格。因此，此次采用全四面体网格划分方法，如图 10所示，最大单元尺寸为10mm。ctdih of Patdh Independent - MethodStopeScaping MettKidGeametry SelectionGeometry19 BodiesDefinrtlonSu|：ipres5edNoMethodTeUatedransAigorllhmPatdi ndependentElemEnt Mldside NlndeUse Global S-FttingAHvgWDefined ByMax Element 史eMaxE ement Siie10a mmFeature Angle30.0&Mesh Based DeTeaturlnganRefinementCLnmatureMin Sizt Umii1. mm图10四面体划分方式整体曲轴网格划分后如图11所示，共包含82408个单元，125538个单元节点。图11更改后曲轴网格划分网格划分质量分布如图12所示，整体网格质量较之前明显改善，且网格质量较差处位 于后凸缘盘边缘位置，此处不承受载荷，对后处理结果影响基本可以忽略。网格划分后，网 格质量与前一次划分做比较，如表6所示。图12 更改后网格质量分布图表6 两次网格划分网格质量比较划分方法度量方法最小值最大值平均值标准差1Skewness0.00160.9980.2900.1752Skewness0.00290.6910.2710.159因此采用调整后的方法得到整个曲轴的网格划分。4. 曲轴强度分析4.1载荷工况曲轴在工作过程中承受的载荷与边界条件有：主轴承支撑作用力、输出端反扭矩、平衡 块离心力、曲柄和连杆轴颈的离心力、连杆对连杆轴颈的载荷等。本文中直列四缸发动机点火顺序为1-3-4-2。当1、2、3、4缸分别点火时活塞和曲轴 位置示意如图9、10所示。图9 1、4缸点火示意图图10 2、3缸点火示意图对于直列四缸发动机，当活塞处于上止点位置时连杆轴颈载荷达到最大值，因此只需考 虑各个气缸分别处于压缩终了，活塞在上止点位置时的受力状况即可。本文中发动机的最大扭矩转速，最大功率。并以此计算工况。由发动机动力学计算得， 在点火上止点附近，由于爆发压力与活塞和连杆惯性的作用使得连杆轴颈受最大载荷。按动力学法计算载荷，并假设作用在轴颈上的载荷为分布载荷，且根据有限宽度轴颈油 膜压力应力分布规律并忽略油孔处压力峰值突变的影响，采用以沿轴线均匀分布，沿圆周方 向120围呈余弦分布的载荷边界条件，如图11所示。图11 连杆轴颈载荷边界条件4.2载荷与约束边界条件4.2.1转速条件在 Workbench 中 Static Structural 中插入项 Rotational Velocity,以主轴颈轴线为 转轴，转速设为188rad/s (ramped)，即，如图12所示。图12 转动条件4.2.2连杆轴颈载荷Workbench中自带了 Bearing Load项，即轴承载荷。其径向分量将根据投影面积来分 布压力载荷，轴向载荷分量沿着圆周均匀分布。由动力学计算可知，在点火上止点位置时， 由于气体爆炸压力与惯性作用使得作用于连杆轴颈表面，方向沿径向指向转轴的最大载荷为; 其他连杆轴颈主要受到连杆惯性力的作用，此时连杆轴颈所受最小载荷为。对于不同气缸点 火，对应有不同的连杆轴颈载荷条件，图13为1缸点火连杆轴颈示意图。图13 1缸点火连杆轴颈示意图4.2.3曲轴扭荷Workbench中可通过力矩/偶载荷(Moment)施加力矩/偶在任意实体表面，如果选择多个 表面则力矩/偶将分摊在这些面上。发动机工作中，由曲轴后端传出扭矩，曲轴必然受到一 定的扭转力，作用在曲轴前端轴和后凸缘盘上。根据公式，因此作用在曲轴上的最大扭矩可 以近似为，如图14所示。图14 扭矩载荷4.2.4支撑约束Workbench中常用的支撑约束为Cylindrical，即圆柱面约束，可以施加在圆柱表面， 并可以指定轴向、径向或者切向自由度的约束状况，但此约束通常仅适用于小变形(线性) 分析。由于曲轴在实际工作中，变形量极小，因此可采用此种约束。本文所研究曲轴为某直 列四缸发动机整体式曲轴，含有5个主轴颈，并通过轴瓦、轴承安装在发动机箱体，因此5 个主轴颈处可施加轴向和径向固定，切向自由的圆柱面约束，如图15所示。Pl Cylhdrfcal OtrnmDetails of Cylindrical Support-曰 SC&peScoping Method5旦如祀AxialTangentialSuppressedI Geometry曰 DefinitionTypeRadial5 FacesCylindrical SupportFixed图15主轴颈支撑约束4.3有限元求解结果分析对于主应力和主应变，根据力学基础理论，应力量可以转成只显示法向应力的形式。这 就是三个主应力，即，其值可以被指定，且三个主应力值带有方向。常用于塑性材料的理论 是最大等效应力和最大剪切应力。本文采用最大等效应力分析。等效应力Equivalent(von-Mises)，也称von-Mises应力，其定义为：总变形量Total Deformation，其定义为:4.3.1 1缸点火结果分析1缸点火时，曲轴的等效应力图和总变形量图分别如图16、17所示。其中，最大等效 应力为130.35MPa，位于1缸连杆轴颈与曲柄、主轴颈与平衡重的连接处。最大的变形量为 0.05471mm，位于连杆轴颈与平衡重的边缘处。图16等效应力图图17总变形量图4.3.2 2缸点火结果分析2缸点火时，曲轴的等效应力图和总变形量图分别如图18、19所示。其中，最大等效 应力为133.2MPa，位于2缸连杆轴颈与曲柄的连接处。最大的变形量为0.05224mm，位于连 杆轴颈与平衡重的边缘处。图18等效应力图图19总变形量图4.3.3 3缸点火结果分析3缸点火时，曲轴的等效应力图和总变形量图分别如图20、21所示。其中，最大等效 应力为127.5MPa，位于2缸连杆轴颈与曲柄的连接处。最大的变形量为0.05215mm，位于连 杆轴颈与平衡重的边缘处。图20等效应力图图21总变形量图4.3.4 4缸点火结果分析4缸点火时，曲轴的等效应力图和总变形量图分别如图22、23所示。其中，最大等效 应力为124.8MPa，位于2缸连杆轴颈与曲柄的连接处。最大的变形量为0.05425mm，位于连 杆轴颈与平衡重的边缘处。图22等效应力图图23总变形量图4.4加密收敛分析分析上述四种情况得，2缸点火时，最大等效应力为133.2MPa，最大变形量为0.05224mm。 因此为了，减少计算量，对2缸点火时，对曲轴加密网格划分，检查应力值与变形量是否稳 定于一个值。对于2缸连杆轴颈与曲柄进行局部加密，单元尺寸调整至最大6mm，加密后如图24所 示。共包含177312个单元，260480个单元节点。单元质量分布图如图25所示。图24局部加密后网格划分图图25单元质量分布图此时曲轴的等效应力图和总变形量图分别如图26、2 7所示。最大等效应力为149.3MPa, 最大的变形量为0.05272mm。图26等效应力图图27总变形量图加密后最大应力值从133.2MPa增大至149.3MPa，应力增幅约12%。最大的变形量从 0.05224增大至0.05272,变形增幅约0.92%。加密后，应力值仍处于未加密前的应力水平， 最大变形量几乎未变，因此可以认为采用Workbench进行的曲轴静力学强度分析数据基本可 靠。根据材料属性可得，QT-800的屈服强度为480MPa，塑性较好，对于塑性材料安全系数 一般选取1.2-2.5，本文取安全系数为2,因此，该曲轴强度校核合格。5. 曲轴自由模态分析振动模态是弹性结构固有的、整体的特性。通过模态分析方法可以了解结构在某一易受 影响的频率围的各阶主要模态的特性，从而进一步分析结构在此频段在外部或部各种振源作 用下产生的实际振动响应。因此，模态分析是结构动态设计与设备故障诊断的重要方法。对曲轴进行模态分析，可以确定曲轴的固有频率和振型。计算获取各阶固频率和观察振 型变形图，为曲轴避免共振设计提供参考。模态分析一般分为自由模态分析和约束模态分析，自由模态分析与对模型不施加任何约 束，而约束模态分析则需考虑结构所受的实际载荷，从而对模型添加合适的载荷约束，但如 果约束施加不当，则可能造成固有频率的巨大偏差或者遗漏。因此，本文采用自由模态分析。5.1网格划分网格划分仍采用先前调整后的网格划分法，即采用高级尺寸函数中的Proximity and Curvature函数进行全局设置，并采用全四面体的网格划分方法，最大网格尺寸为10mm。总 体网格质量良好，满足模态分析要求。5.2约束条件由于进行自由模态分析，对曲轴不施加任何约束。5.3参数设置Workbench模态分析(Modal)中，程序可设定的模态数为1-200，默认值为6。频率围为 0HzTe+08Hz。此设定从工程应用角度看完全足够。对于机械结构的共振问题，一般关心较低阶次的频率，且由于曲轴自由模态分析的前6 阶频率仅反映刚体振动，因此共取12阶模态进行分析，设置如图28所示。Details of Analysis Settings日OptionsMax Modes to Find12Limirt Search to RangeNo曰Solver ControlsDampedNoSalver Type-Program ControlledElRotordynamics Controls;iOurtpurt Controls*Anatysis; Data MsnagEment图28模态分析设置参数5.4结果分析曲轴的12阶自由模态频率图如图29所示。图29曲轴12阶自由模态频率图由于曲轴的前6阶自由模态频率接近于零，为刚体模态，对于本次分析意义不大，因此 将后6阶模态频率求解结果列入表7。表7曲轴非零模态频率与最大位移1阶2阶3阶4阶5阶6阶频率/Hz300.93418.08725.60804.81822.901104.60位移/mm12.7913.5914.2017.9813.5215.21非零模态频率振型图如图30图35所示。需要注意的是，振型显示模式变形量的放大 比例为2.3:1。QMocklT OTd Cefcnr-Jdiifl 7T j/pt: TDsbn WoflFrjua JCO.yJ Hz Unit ee 也IMS 1皿IJ.T9*2 陋箪 11,174 企炒7.12H 5i.TO8图301阶振型间成 D?iormrfri TiZT心 Dwlfr图312阶振型态，Jg顼Tc-tJ 匚yfonn 占#9 Typ Totd Cefcrmarion fwcqu?ry ?256-mt Llhiu mm 旬 i#i arm*:b4.1S7Mj)illil63)Q Medal如 10 Tjps: Told Crtfmafiian FfAquA:驻日 El H? UWe mm Ml 迥15.-9M1-4X1Q412JQ1510jQ2$flB556ZM594。眼掘胃加院加裾MSTotwijlTyif.fi! Tg k口血仙 誓 rguifl 毕 S225 H】 Jnic twi19:10*.临 MX52 3J1S4 T.Sm MJ1S3 gw 3-OW LS1ZOliK&fr MinQMqd*l 3 35二miTlkiH 12T炸吧：TW DiFcEddcfl FrHqufni H Mb Hr UCP 13.52& :1LHU 皿1公图323阶振型图334阶振型图345阶振型图356阶振型曲轴1、3、4、6阶非零模态从振型图中可以看出，以弯曲变形为主，而2、5阶非零模 对应的振型图可以看出，以弯曲扭转为主。曲轴最低固有频率，1阶非零模态频率为300.93Hz。该发动机在转速为1800 r/min时， 基频为60Hz，远低于最低固有频率，避开了共振频率。参考文献1 古成中，吴新跃.有限元网格划分与发展趋势J.计算机科学与探索，2008(03):248-259.2 仇亚萍，黄俐军，立飞.基于ANSYS的有限元网格划分方法J.机械管理开发，2007(06):76-77.3 娜.直列四缸柴油机减振分析与平衡机构优化设计D.工程技术大学，2011.4 俊峰.8L250系列柴油机曲轴强度计算与仿真研究D.理工大学，2010.5 必荣.基于ANSYS的小型柴油机曲轴应力分析J.拖拉机与农用运输车，2004(03)：30-32.6 邵康.直列四缸柴油机轴系动力学仿真分析D.大学，2009.7 学民，志琴.基于ANSYS Workbench的四缸发动机曲轴模态分析(英文)J.Journal of Measurement Science and Instrumentation，2015(03)：282-285.8 昌领，建义，清平，仇晨，罗晓兰.基于ANSYS的六缸压缩机曲轴模态分析与谐响应分析J.流体机械，2 012(08)：17-21；26.