CN104537241B - 车轮钢圈疲劳分析方法 - Google Patents
车轮钢圈疲劳分析方法 Download PDFInfo
- Publication number
- CN104537241B CN104537241B CN201410839673.6A CN201410839673A CN104537241B CN 104537241 B CN104537241 B CN 104537241B CN 201410839673 A CN201410839673 A CN 201410839673A CN 104537241 B CN104537241 B CN 104537241B
- Authority
- CN
- China
- Prior art keywords
- mrow
- stress
- msub
- fatigue
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明旨在提供一种车轮钢圈疲劳分析方法,包括以下步骤:A、预设边界约束条件、材料参数与工况条件,用有限元法建立被测车轮钢圈模型;B、对模型施加6个方向的载荷,建立应力云图,确定各方向载荷下应力集中点作为危险点,在6个方向中确定危险载荷方向;C、测量各云图确定的危险点和经验已知的危险点的Mises应力、X、Y以及Z方向应力,确定最大方向应力为危险点的应力方向;D、结合预设的材料参数及步骤C得到的应力数据计算各个测试区域的疲劳寿命N,计算出的疲劳寿命中的最小值即为被测钢圈的最小疲劳寿命。该方法克服现有技术车轮钢圈设计开发周期长,成本高的缺陷,具有计算精确、运算效率高的特点。
Description
技术领域
本发明涉及汽车零部件涉及领域,具体涉及一种车轮钢圈疲劳分析方法
背景技术
车轮是衡量当今汽车产业发展水平的重要因素,汽车使用中一系列的安全性,包括轮胎磨损、燃油消耗及人员本身的操作方便程度都系于车轮在为汽车提供动力时的准确度,因此车轮为整车提供质量上的优越程度就变得很重要。鉴于在行驶过程中随机载荷中主要为疲劳载荷时常作用于车轮,由于车轮出现损毁的主要来源是工作行驶中产生的弯曲载荷,因此疲劳寿命成为车轮制造中主要研究的问题,即确保破坏与断裂尽量不发生在车轮的工作年限之内。现有车轮钢圈设计大多是采用经验公式计算或者参照样机进行,然后通过试样不断对结构进行修改优化,直到满足工作要求为止,因此存在开发周期长,成本高的缺陷。
发明内容
本发明旨在提供一种车轮钢圈疲劳分析方法,该方法克服现有技术车轮钢圈设计开发周期长,成本高的缺陷,具有计算精确、运算效率高的特点。
本发明技术方案如下:一种车轮钢圈疲劳分析方法,包括以下步骤:
A、预设边界约束条件、材料参数与工况条件,用有限元法建立被测车轮钢圈模型;
B、对模型施加6个方向的载荷,各方向彼此之间成60度角,建立应力云图,确定各方向载荷下应力集中点作为危险点,将6个方向中的载荷条件下的最大应力值最高的方向作为危险载荷方向;
C、测量各云图确定的危险点和经验已知的危险点的Mises应力、X方向应力、Y方向应力以及Z方向应力,确定最大方向应力为危险点的应力方向;
D、结合预设的材料参数及步骤C得到的应力数据计算各个测试区域的的疲劳寿命N,计算出的疲劳寿命中的最小值即为被测钢圈的最小疲劳寿命。
所述的步骤D中,疲劳寿命N的计算式如下:
其中,Δε为局部应变幅值,σf'为疲劳强度系数,εf'为疲劳延性系数,E为弹性模量,b为疲劳强度指数,c为疲劳延性指数,σmax为最大应力值;b’为结合了零件的表面加工与尺寸影响后的修正疲劳强度指数,其值为:
其中,σ-1为对称弯曲疲劳极限,ε为尺寸系数,β为表面系数,σm为平均应力。
作为其优选方案,所述经验已知的危险点包括:钢圈螺栓孔处、钢圈与轮辋焊缝处、轮辐与法兰轴连接大孔处。
本发明疲劳寿命N计算式化简过程如下:
改进史密斯公式为:
本发明车轮钢圈疲劳分析方法通过有限元法建立钢圈模型,并对其施加6个方向的载荷进行应力分析,将应力集中的区域作为危险点,后续的应力测试则针对危险点(薄弱环节)进行,保证能够获取到被测钢圈最小的疲劳寿命,减小测试误差;6个方向载荷的施加,可以确保本发明方案能够分析出在不同位置不同情况下钢圈应力的变化情况,保证测试分析的准确与全面;后续对于危险点测量得到各个方向的载荷情况下的Mises应力、X方向应力、Y方向应力以及Z方向应力,根据得到的应力数据以及被测钢圈的材料参数,通过史密斯公式计算得到被测钢圈的疲劳寿命;上述化简后的史密斯公式变量更少,计算更为简单,计算效率更高。
附图说明
图1为本发明车轮钢圈疲劳分析方法的流程图
图2为实施例1的6个载荷方向示意图
具体实施方式
下面结合附图和实施例具体说明本发明。
本实施例车轮钢圈疲劳分析方法包括以下步骤:
本实施例用有限元法建立被测车轮钢圈模型,预设工况为弯矩为20Kn.m,边界约束条件为在被测钢圈同侧挡圈圆环面上施加约束,对模型施加6个方向的载荷,各方向彼此之间成60度角,建立应力云图,确定各方向载荷下应力集中点作为危险点,本实施例应力云图上得到的危险点为通风孔小孔处以及通风孔缺口方向两侧,经验已知的危险点包括:钢圈螺栓孔处、钢圈与轮辋焊缝处、轮辐与法兰轴连接大孔处,将上述6个危险点分别定义为测试点1-6,利用Workbench中的Stress Tool工具获取测试点1-6的应力数据,如表1-6所示:
表1测试点1的应力计算结果
方向 | X轴方向应力(Mpa) | Y轴方向应力(Mpa) | Z轴方向应力(Mpa) | Mises应力(Mpa) |
1 | 8.5048 | 19.983 | 23.905 | 49.449 |
2 | 5.8275 | 11.439 | 14.124 | 28.417 |
3 | 2.7342 | 9.0199 | 9.9385 | 21.707 |
4 | -8.459 | -19.791 | -23.748 | 49.06 |
5 | -5.7169 | -11.334 | -14.008 | 28.106 |
6 | 2.7597 | 8.6929 | 9.8325 | 21.345 |
表2测试点2的应力计算结果
方向 | X轴方向应力(Mpa) | Y轴方向应力(Mpa) | Z轴方向应力(Mpa) | Mises应力(Mpa) |
1 | -4.5539 | -1.959 | -20.771 | 26.648 |
2 | -12.202 | -10.315 | -74.317 | 92.861 |
3 | -7.7859 | -8.6783 | -52.935 | 66.027 |
4 | 4.6138 | 1.8875 | 20.791 | 26.664 |
5 | 12.005 | 10.429 | 72.59 | 90.806 |
6 | 7.2712 | 8.3371 | 51.865 | 64.411 |
表3测试点3的应力计算结果
方向 | X轴方向应力(Mpa) | Y轴方向应力(Mpa) | Z轴方向应力(Mpa) | Mises应力(Mpa) |
1 | -4.5728 | -2.071 | -21.555 | 26.743 |
2 | 7.9253 | 8.3465 | 56.519 | 70.258 |
3 | 12.346 | 10.189 | 78.001 | 96.845 |
4 | 4.9438 | 2.2953 | 21.307 | 27.118 |
5 | -7.6496 | -8.5871 | -54.546 | 67.899 |
6 | -12.386 | -10.678 | -76.67 | 95.595 |
表4测试点4的应力计算结果
方向 | X轴方向应力(Mpa) | Y轴方向应力(Mpa) | Z轴方向应力(Mpa) | Mises应力(Mpa) |
1 | -0.78167 | -19.21 | -23.441 | 21.309 |
2 | -2.618 | -11.98 | -6.3281 | 9.2744 |
3 | -1.7888 | 7.5567 | 17.156 | 16.425 |
4 | 0.79081 | 19.411 | 23.482 | 21.385 |
5 | 2.5572 | 11.69 | 6.4606 | 9.0768 |
6 | 1.7654 | -7.3199 | -16.66 | 15.972 |
表5测试点5的应力计算结果
方向 | X轴方向应力(Mpa) | Y轴方向应力(Mpa) | Z轴方向应力(Mpa) | Mises应力(Mpa) |
1 | -0.55158 | 4.1786 | -3.7229 | 10.806 |
2 | -3.1035 | -1.1759 | -4.2292 | 5.3022 |
3 | -2.2228 | -5.4052 | -0.00526 | 11.134 |
4 | -0.17186 | -2.1964 | 2.0819 | 8.0038 |
5 | 0.63999 | -0.37207 | 0.99632 | 1.4141 |
6 | 0.83417 | 1.8488 | -0.93678 | 6.8572 |
表6测试点6的应力计算结果
方向 | X轴方向应力(Mpa) | Y轴方向应力(Mpa) | Z轴方向应力(Mpa) | Mises应力(Mpa) |
1 | -0.46127 | -31.1 | -25.805 | 28.397 |
2 | -0.21304 | -15.787 | -12.033 | 14.092 |
3 | 0.24903 | 15.32 | 13.874 | 14.477 |
4 | 0.45331 | 30.764 | 25.645 | 28.133 |
5 | 0.16579 | 14.804 | 11.606 | 13.346 |
6 | -0.22313 | -13.893 | -13.051 | 13.335 |
由上述数据,可以计算得到各个测试点的最大应力值σmax以及平均应力σm,并结合预设的材料参数代入式(1),即可求得各测试点疲劳寿命,上述疲劳寿命的最小值即为被测钢圈的疲劳寿命;
本实施例材料参数如表7、表8所示:
表7实施例1材料参数
疲劳强度系数σf'(Mpa) | 658.8 |
疲劳延性系数εf'(mm) | 0.2747 |
弹性模量E(Gpa) | 198.753 |
疲劳强度指数b | -0.0709 |
疲劳延性指数c | -0.4907 |
对称弯曲疲劳极限σ-1(Mpa) | 210 |
尺寸系数ε | 0.73 |
表面系数β | 1 |
表8实施例1材料参数中的局部应变幅值Δε
由上述参数结合应力数据计算得到本实施例被测钢圈的最小寿命为3051210次。
本实施例仅是对弯矩为20Kn.m的工况进行疲劳寿命分析举例,对于不同工况均可参照上述方法进行疲劳寿命分析。
Claims (1)
1.一种车轮钢圈疲劳分析方法,其特征在于包括以下步骤:
A、预设边界约束条件、材料参数与工况条件,用有限元法建立被测车轮钢圈模型;
B、对模型施加6个方向的载荷,各方向彼此之间成60度角,建立应力云图,确定各方向载荷下应力集中点作为危险点,将6个方向中的载荷条件下的最大应力值最高的方向作为危险载荷方向;
C、测量各云图确定的危险点和经验已知的危险点的Mises应力、X方向应力、Y方向应力以及Z方向应力,确定最大方向应力为危险点的应力方向;
D、结合预设的材料参数及步骤C得到的应力数据计算各个测试区域的的疲劳寿命N,计算出的疲劳寿命中的最小值即为被测钢圈的最小疲劳寿命;
所述的步骤D中,疲劳寿命N的计算式如下:
<mrow>
<mi>N</mi>
<mo>=</mo>
<mfrac>
<mroot>
<mfrac>
<mrow>
<msub>
<mi>E&sigma;</mi>
<mi>max</mi>
</msub>
<mi>&Delta;</mi>
<mi>&epsiv;</mi>
</mrow>
<mrow>
<mn>2</mn>
<mo>&CenterDot;</mo>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>3</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mrow>
<mn>3</mn>
<msup>
<mi>b</mi>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<mi>c</mi>
</mrow>
</mroot>
<mn>2</mn>
</mfrac>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,Δε为局部应变幅值,σf'为疲劳强度系数,εf'为疲劳延性系数,E为弹性模量,c为疲劳延性指数,σmax为最大应力值;b’为结合了零件的表面加工与尺寸影响后的修正疲劳强度指数,其值为:
<mrow>
<msup>
<mi>b</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>lg</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mi>&epsiv;</mi>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>lg</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>f</mi>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>lg&sigma;</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>lg</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>f</mi>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mi>b</mi>
</mrow>
其中,σ-1为对称弯曲疲劳极限,ε为尺寸系数,β为表面系数,σm为平均应力,b为疲劳强度指数;
所述经验已知的危险点包括:钢圈螺栓孔处、钢圈与轮辋焊缝处和轮辐与法兰轴连接大孔处;
其中σf'=658.8,εf'=0.2747,E=198.753,b=-0.0709,c=-0.4907,σ-1=210,ε=0.73,β=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410839673.6A CN104537241B (zh) | 2014-12-30 | 2014-12-30 | 车轮钢圈疲劳分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410839673.6A CN104537241B (zh) | 2014-12-30 | 2014-12-30 | 车轮钢圈疲劳分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104537241A CN104537241A (zh) | 2015-04-22 |
CN104537241B true CN104537241B (zh) | 2018-05-11 |
Family
ID=52852767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410839673.6A Expired - Fee Related CN104537241B (zh) | 2014-12-30 | 2014-12-30 | 车轮钢圈疲劳分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104537241B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104931240B (zh) * | 2015-05-25 | 2017-09-01 | 广西科技大学 | 一种钢圈故障自识别方法 |
CN106021688B (zh) * | 2016-05-13 | 2019-05-21 | 中车青岛四方机车车辆股份有限公司 | 一种轨道交通机车的车轮强度的评估方法以及装置 |
CN109308393B (zh) * | 2018-09-30 | 2021-04-23 | 中车齐齐哈尔车辆有限公司 | 一种车体疲劳寿命的评估方法和系统 |
CN109829185A (zh) * | 2018-12-13 | 2019-05-31 | 铠龙东方汽车有限公司 | 一种基于预紧力的螺母松动分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101393079A (zh) * | 2008-11-06 | 2009-03-25 | 上海交通大学 | 汽车车身结构疲劳寿命预测系统 |
CN101592552A (zh) * | 2009-07-01 | 2009-12-02 | 同济大学 | 预测汽车后悬架多轴疲劳的方法 |
CN102261988A (zh) * | 2011-04-25 | 2011-11-30 | 吴峰崎 | 一种起重机剩余寿命的测试方法 |
CN102854010A (zh) * | 2012-10-10 | 2013-01-02 | 湖南奔腾动力科技有限公司 | 一种基于道路循环工况的发动机零件疲劳寿命计算方法 |
-
2014
- 2014-12-30 CN CN201410839673.6A patent/CN104537241B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101393079A (zh) * | 2008-11-06 | 2009-03-25 | 上海交通大学 | 汽车车身结构疲劳寿命预测系统 |
CN101592552A (zh) * | 2009-07-01 | 2009-12-02 | 同济大学 | 预测汽车后悬架多轴疲劳的方法 |
CN102261988A (zh) * | 2011-04-25 | 2011-11-30 | 吴峰崎 | 一种起重机剩余寿命的测试方法 |
CN102854010A (zh) * | 2012-10-10 | 2013-01-02 | 湖南奔腾动力科技有限公司 | 一种基于道路循环工况的发动机零件疲劳寿命计算方法 |
Non-Patent Citations (6)
Title |
---|
"evaluation of fatigue life aluminum alloy wheels under bending loads";P.R. Raju etal;《fatigue &fracture of engineering materials & structures》;20090212;第32卷(第2期);第119-126页 * |
"fatigue life prediction for wheels by simulation of the rotating bending test";H.M.karandikar etal;《international congress & exposition》;19900201;第180-190页 * |
"汽车钢圈疲劳分析与寿命估算方法研究";刘新宇;《中国优秀硕士论文全文数据库 工程科技Ⅰ辑》;20071115(第5期);第20-25页第3.4-3.6节,第39-41页第4.3.2节 * |
"汽车钢圈的疲劳寿命分析及优化设计";朱红建 等;《中国优秀硕士论文全文数据库 工程科技Ⅱ辑》;20110315(第3期);第C35-43页 * |
"用改进的史密斯公式预测汽车车轮疲劳寿命";杨占春 等;《哈尔滨工业大学学报》;20001231;第32卷(第6期);第101页第2节 * |
"铝合金车轮弯曲疲劳寿命的仿真分析与试验研究";王良模 等;《南京理工大学学报(自然科学版)》;20091030;第33卷(第5期);第571-575页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104537241A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104537241B (zh) | 车轮钢圈疲劳分析方法 | |
CN101826125B (zh) | 一种设计麦弗逊悬架的方法 | |
CN201607377U (zh) | 用于测试轮胎/车轮组件的测试机 | |
CN102628769B (zh) | 一种含表面裂纹缺陷承压设备的定量风险分析方法 | |
CN103077320A (zh) | 一种基于粗糙集理论的焊接疲劳分析方法 | |
CN105699220A (zh) | 一种正交异性钢桥面板疲劳寿命的计算方法 | |
CN104956196B (zh) | 表征车辆的行动性的方法以及在车辆轮胎的选择中的使用 | |
CN104515685A (zh) | 基于道路载荷的扭转梁后桥耐久性评价方法 | |
CN101984333B (zh) | 载重车辆翻新轮胎胎体剩余寿命预报方法 | |
WO2009101841A1 (ja) | 空気圧低下検出方法におけるパラメータの設定方法 | |
Farrahi et al. | A comparative study on the fatigue life of the vehicle body spot welds using different numerical techniques: Inertia relief and Modal dynamic analyses | |
JP2017062205A (ja) | 溶接部の変形抵抗曲線の算出方法、溶接部を備えた部材の製造方法、プログラム、および、プログラムを記録したコンピュータ読み取り可能な記録媒体 | |
CN104749200A (zh) | 基于转向架焊接缺陷估算转向架使用寿命的方法 | |
CN104897497A (zh) | 一种轮胎接地压力分布分析提高耐磨性的方法 | |
CN107121236A (zh) | 基于矢量计算的轮胎静不平衡改善方法 | |
CN104933309A (zh) | 一种计算两端具有弹性约束的待测杆的屈曲载荷的方法 | |
CN107764493B (zh) | 一种振动耐久试验psd的获取方法及获取装置 | |
Gonera et al. | Analysis of car body and wheel geometry depending on mileage on the wear of tyres in passenger cars | |
Zhan et al. | Fatigue analysis of weld region in torsion beam rear suspension system | |
CN103592066A (zh) | 用于转向桥轴头力的测量方法及标定装置 | |
CN110245431B (zh) | 飞行器液压管路流固耦合振动断裂故障的溯源方法 | |
Mars et al. | Breaking the computational barrier to simulating full road load signals in fatigue | |
CN106326592A (zh) | 一种轻量化钢制车轮的设计方法 | |
Gobbi et al. | Design of a Lightweight Truck Wheel | |
Solf et al. | Analysis of the grinding force in plunging bevel gear grinding |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180511 Termination date: 20181230 |