基于ABAQUS的飞轮强度计算与优化分析.docx

基于ABAQUS的飞轮强度计算与优化分析

廖礼平，陆金华

（柳州五菱柳机动力有限公司，广西柳州545005）

摘要：应客户要求，对飞轮的结构进行了重新设计。利用非线性有限元软件ABAQUS计算新状态下飞轮在热负荷、离心力作用

下的变形、温度和应力分布，讨论飞轮的疲劳强度，提供结构优化方案，为飞轮的可靠性设计提供依据。

关键词：非线性；ABAQUS；飞轮；疲劳强度；可靠性

前言

飞轮是发动机曲轴系中重要组成部分，起到的作用非常大。当发动机合成转矩大于总阻力矩的瞬间，飞轮

把剩余转矩（盈功）储存起来使发动机转速不会过大地升高，反之，当发动机合成转矩小于总阻力矩的瞬间，

欠缺的亏损转矩（亏功）就由飞轮储存的盈功释放出来补充，使发动机转速不会过大地下降。总之，飞轮主要

是保证发动机启动和最低最低稳定转速能正常运行，并作为扭矩输出与负载间的连接、盘车机构以及刻示曲轴

转角标志

【1】。

本文将考虑飞轮与离合器总成间的热传导，将所计算得到的各点温度载荷加载到模型上，并与离心力耦合

计算，从而得到飞轮的强度。

1飞轮强度分析

1.1分析模型

本文研究对象为飞轮总成。发动机工作过程中，飞轮与曲轴之间通过螺栓进行连接，因此采用CAD软件UG

删除飞轮螺栓孔的简化处理。将处理好的3D实体数模导入到SimLab进行网格划分，飞轮总成采用C3D10M单元，

网格单元大小为2mm，整个模型的节点总数为1075552，单元总数为709985。网格模型见图1。

图1网格模型

1.2边界条件

1.2.1热传导计算边界条件

热传导计算的边界条件主要由试验和经验值所得，数据见表1；边界设置见图2。

表1热传导边界条件

接触区域传热系数(w/mm^2.s.℃)环境温度（℃）

①与摩擦片接触区域2e-4300

②与曲轴配合面--100

③其它区域6.322e-005100

图2热传导边界设置

热传导计算到的温度分布结果见下图3：

图3热传导温度分布结果

1.2.2耦合边界条件

曲轴与飞轮接触端面旋转轴方向位移约束；将热传导计算中得到的各点温度作为温度载荷加载到模型上；

将旋转速度6000rpm产生的离心力加载到整个模型上，设置见图4。

图4耦合边界条件

1.3计算工况

整个分析过程分两个工况：热应力工况、耦合工况。

（1）热应力工况：仅考虑热传导产生的温度载荷作用；

（2）耦合工况：考虑温度载荷及离心力共同作用。

1.4材料属性

飞轮的材料为QT500-7，其属性见表2所列。

表2飞轮材料属性

名称数值单位

弹性模量154000MPa

泊松比0.275--

密度7.15e-009t/mm^3

导热系数3.52e-002w/mm.℃

比热4.61e+005J/ton.K

膨胀系数1.25e-0051/K

2计算结果

2.1热应力工况结果

图5为热应力工况下的应力分布，在与曲轴连接端面凸台上的应力最大，值为129MPa。图6为应力最大

区域的放大图。

图5热应力分布图6应力区域放大图

2.2耦合工况

根据疲劳理论，选取飞轮相关位置（图7--图9中的①②③处）的最大主应力、最小主应力、S,Mises应力

来计算应力幅值、平均应力，据此来判断疲劳强度是否满足使用要求。

（1）最大主应力结果

图7最大主应力结果

（2）最小主应力结果

图8最小主应力结果

（3）S,Mises应力结果

图9S，Mises应力

统计上述三个位置的最大主应力和最小主应力、S,Mises应力结果，计算出上述3个位置的应力幅值和平

均应力值，见表3。

表3应力结果统计（单位：MPa）

位置最大主应力最小主应力应力幅值平均应力S，Mises应力

①222.45.5108.45113.95191.7

②152.10.17676.1149.2

③4.6-228.8116.7-112.1227.6

将①②③位置的应力幅值、平均应力置于材料QT500的疲劳曲线内，如图10所示，根据goodman疲劳理

论可知，位置①处在红线上方，超出疲劳限值范围，会出现疲劳损伤。

3、结构优化及强度验证

3.1结构优化

由图10中海格图的疲劳强度结果可知：飞轮需要进行结构优化，以便达到强度要求。根据图9中所示的强

度结果，对飞轮做如下优化设计：第一，两个阶梯端面间倒0.5mm的圆角，见图11圈中部位所示；第二，飞轮

最小壁厚由2mm变为3mm，见图11直线厚度。

③

①

②

图10疲劳强度分布

图11结构优化方案

3.2强度验证

保持热传导、热应力工况和耦合工况的边界条件、载荷不变，网格单