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

时间：2020-09-30 15:07:27 来源：蒲公英阅读网 本文已影响人

相关热词搜索：曲轴 发动机 直列

　 有限元分析课程报告

　直列四缸发动机曲轴有限元分析

　 姓名：

　学号：

　 分数：

　 年

　月

　日

　目 录 1.引言 ....................................................................................................................................... 1 2.曲轴有限元模型的建立 ........................................................................................................ 2 3.曲轴网格划分 ........................................................................................................................ 4 3.1 确定物理场和网格划分法 ........................................................................................ 4 3.2 确定全局网格参数设置 ............................................................................................ 4 3.3 确定局部网格参数设置 ............................................................................................ 4 3.4 网格质量检查 ............................................................................................................ 5 3.5 调整网格划分 ............................................................................................................ 6 4.曲轴强度分析 ........................................................................................................................ 8 4.1 载荷工况 ..................................................................................................................... 8 4.2 载荷及约束边界条件 ................................................................................................. 8 4.3 有限元求解结果分析 .............................................................................................. 10 4.4 加密收敛分析 .......................................................................................................... 12 5.曲轴自由模态分析 .............................................................................................................. 14 5.1 网格划分 ................................................................................................................... 14 5.2 约束条件 .................................................................................................................. 14 5.3 参数设置 .................................................................................................................. 14 5.4 结果分析 .................................................................................................................. 14 参考文献 ................................................................................................................................. 17

　欢迎下载 1.引言 曲轴（如图 1 所示）是发动机中最重要、承载最复杂的零件之一。其造价约占一台发动机总价的 20%-30%，它的使用寿命决定发动机的使用寿命。在工作过程中，曲轴承受着气缸内的气体压力及往复和旋转质量惯性力引起的周期性变化载荷。这些扭转力矩和大小、方向不断变化的弯曲力，在曲轴内产生交变的弯曲及切应力，使曲轴发生疲劳破坏。同时，在周期性变化的载荷作用下，曲轴可能在发动机转速范围内发生共振，加速曲轴的疲劳破坏。

　由于曲轴是发动机中典型的易损件之一，其强度和振动特性都会影响到整机的工作性能，因此对曲轴进行有限元分析，研究曲轴的应力和应变状态，了解其固有频率和模态振型，对于设计和优化曲轴结构有重要意义。

　图 1

　 曲轴实物图 曲轴的基本结构如图 2 所示，主要由前端轴、主轴颈、连杆轴颈、曲柄、平衡重和后凸缘盘等部分构成，前端轴主要用来安装止推垫圈、启动抓和皮带轮等部件；曲柄用来连接主轴颈和连杆轴颈；平衡重的作用是平衡曲轴工作时产生的离心惯性力及其力矩；后凸缘盘用来安装飞轮等部件。

　1

　 2

　3

　4

　5

　 6

　 图 2

　曲轴基本结构

　1-前端轴

　2-主轴颈

　3-平衡重

　4-曲柄

　5-连杆轴颈

　6-后凸缘盘

　欢迎下载

　2.曲轴有限元模型的建立 本文所研究曲轴为某直列四缸发动机整体式曲轴，含有 5 个主轴颈，4 个连杆轴颈、曲柄和平衡重。由于曲轴的具体尺寸参数设计较为复杂，本文不详细叙述。同时，由于曲轴结构复杂，利用有限元软件进行建模时很难保证与图纸上的曲轴结构完全一致，因此建模时必须简化模型。为了减少应力集中，曲轴上不同截面的结合处都有半径不同的倒角，如果在建模时考虑这些倒角和油孔，则会使有限元的网格非常密集，大大增加了模型的单元数量且生成的网格形状也不理想，降低了求解精度。因此，建模时忽略了对分析结果影响较小的特征，如轴颈上的油孔和后凸缘盘上的环槽。现给出曲轴结构的主要设计尺寸如表 1。

　表 1

　曲轴结构设计尺寸 设计参数 主轴颈 长度 主轴颈 直径 连杆轴颈长度 连杆轴颈直径 曲柄厚度 平衡重 厚度 尺寸/mm 35 80 40 66 25 22.5

　平衡重的二维设计尺寸如图 3 所示。

　图 3

　 平衡重设计尺寸 考虑到在 ANSYS 中直接进行曲轴建模的复杂性，本文通过 SolidWorks 建立三维几何模型，并与 ANSYS Workbench 无缝连接，进行有限元分析。曲轴三维建模后如图 4 所示。

　欢迎下载 图 4

　曲轴三维实体图 本文中曲轴材料为 QT-800，其基体组织为珠光体或索氏体，主要用于所受载荷较大、受力复杂的汽车、拖拉机的曲轴、连杆、凸轮轴等等场合。表 2 列出了 QT-800 的材料性能参数。

　表 2

　 QT-800 材料参数 材料 抗拉强度/MPa 抗拉屈服强度/MPa 硬度/HBS 弹性模量/ 密度/ 泊松比 QT-800 ≥800 ≥480 245-355

　7850 0.3

　在 workbench 中设置曲轴材料参数如图 5 所示。

　图 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

　全局网格参数 Relevance Center Element Size Smoothing Transition Span Angle Center Coarse Default Medium Fast Coarse

　 图 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。

　表 4 Skewness 网格质量度量等级 0-0.25 0.25-0.50 0.50-0.80 0.80-0.95 0.95-0.98 0.98-1.00 极好 很好 好 可接受 坏 不可接受

　网格划分后，网格质量如表 5 所示。网格划分质量分布如图 8 所示，整体网格质量较好，小于 0.95 的单元个数为 42775，不可接受单元个数为 3，占总体 0.07‰,但可以发现网格质量较差处多位于连杆轴颈处，因此后续分析中可能会出现应力奇异现象或者网格加密不收敛。因此对于此网格划分，需要进一步调整，使其连接处网格质量提高。

　表 5

　网格质量 度量方法 最小值 最大值 平均值 标准差 Skewness 0.0016 0.998 0.290 0.175

　欢迎下载

　图 8

　网格质量分布图

　3.5 调整网格划分 重新确定全局网格参数设置，采用高级尺寸函数(Advanced Sizing Function),此函数主要控制曲面在曲率较大的地方的网格。其中 Proximity and Curvature 函数采用曲率法确定细化边和曲面处的网格大小，同时又控制模型邻近区网格生成，对于窄薄处网格也有较好的划分效果。更改后全局网格具体参数如图 9 所示。

　 图 9

　 更改后全局网格参数 由于先前网格划分在轴颈处网格质量较差，于是推测为 sweep 方法下六面体网格与四面体网格过渡过程中，出现较差质量的网格。因此，此次采用全四面体网格划分方法，如图10 所示，最大单元尺寸为 10mm。

　图 10

　四面体划分方式

　欢迎下载 整体曲轴网格划分后如图 11 所示，共包含 82408 个单元，125538 个单元节点。

　图 11

　更改后曲轴网格划分 网格划分质量分布如图 12 所示，整体网格质量较之前明显改善，且网格质量较差处位于后凸缘盘边缘位置，此处不承受载荷，对后处理结果影响基本可以忽略。网格划分后，网格质量与前一次划分做比较，如表 6 所示。

　图 12

　 更改后网格质量分布图

　表 6

　 两次网格划分网格质量比较 划分方法 度量方法 最小值 最大值 平均值 标准差 1 Skewness 0.0016 0.998 0.290 0.175 2 Skewness 0.0029 0.691 0.271 0.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 所示。

　图 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、27 所示。最大等效应力为 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。频率范围为 0Hz-1e+08Hz。此设定从工程应用角度看完全足够。

　对于机械结构的共振问题，一般关心较低阶次的频率，且由于曲轴自由模态分析的前 6阶频率仅反映刚体振动，因此共取 12 阶模态进行分析，设置如图 28 所示。

　图 28

　模态分析设置参数

　5.4 结果分析 曲轴的 12 阶自由模态频率图如图 29 所示。

　欢迎下载

　图 29

　曲轴 12 阶自由模态频率图 由于曲轴的前 6 阶自由模态频率接近于零，为刚体模态，对于本次分析意义不大，因此将后 6 阶模态频率求解结果列入表 7。

　表 7

　曲轴非零模态频率及最大位移

　1 阶 2 阶 3 阶 4 阶 5 阶 6 阶 频率/Hz 300.93 418.08 725.60 804.81 822.90 1104.60 位移/mm 12.79 13.59 14.20 17.98 13.52 15.21

　非零模态频率振型图如图 30—图 35 所示。需要注意的是，振型显示模式变形量的放大比例为 2.3:1。

　图 30

　 1 阶振型

　图 31

　 2 阶振型

　图 32

　 3 阶振型

　欢迎下载

　图 33

　 4 阶振型

　图 34

　 5 阶振型

　图 35

　 6 阶振型

　曲轴 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].流体机械,2012(08):17-21;26.