基于ANSYS_LS_DYNA的滚动轴承仿真与分析good- 摩擦系数
时间:2025-04-20
时间:2025-04-20
ANSYS LS_DYNA 滚动轴承 仿真
第24卷第9期2007年9月
机 械 设 计
JOURNALOFMACHINEDESIGNVol.24 No.9
Sep. 2007
基于ANSYS/LS-DYNA的滚动轴承仿真与分析
张乐乐,高祥,谭南林,樊莉
(北京交通大学轨道交通控制与安全国家重点实验室,北京 100044)
3
摘要:基于ANSYS/LS-DYNA建立了滚动轴承的有限元模型,有效处理了考虑摩擦条件的接触问题,实现了滚动轴承显式动力学的运动过程仿真。并以6003深沟球轴承为例,进行了试算和相关参数的仿真与分析,结果表明,与滚动轴承运动的实际情况基本吻合。不同工况下的参数仿真分析是进一步研究故障轴承运动仿真与分析的基础。
关键词:滚动轴承;显式动力学;仿真中图分类号:TH133;TB122 文献标识码:A 文章编号:1001-2354(2007)09-0062-03 滚动轴承是应用最为广泛的一类机械部件,也是最易受损
的部件之一。在轨道交通方面,滚动轴承是铁路机车车辆安全运行的关键部件之一,全构成重大影响,裂损等共16种,碰到。起的,,或是热处理工艺不过关,这种原始的隐患是对行车安全的最大威胁,其次,滚动轴承在运行中难免会出现不同程度的损伤,而在高速、重载、长交路下运行的列车,轴承的故障往往会迅速扩大,在短时间内造成热轴、燃
[5]
轴、
切轴,而最终导致列车颠覆等重大行车事故。
因此,通过仿真方法掌握轴承的故障特征,及时地、有针对性地监测、检查与维护是避免事故的有效手段[6]。以滚动轴承为研究对象,在物理模型的基础上,采用贴近实际的材料模型,考虑摩擦条件下的接触模型和约束条件,实现基于显式动力学的有限元仿真,并进行实例计算与参数仿真分析,进一步探讨正常轴承运动过程状态特性,是对有损轴承运动故障特征研究的基础。
at=
t-t
-2t
Δt)t+(2)(3)
t=
t-Δt
Δt)t+
将式(2)、式(3)代入式(1)可得求解各个离散时间点解的递推公式如下:
Δt2
Δt2
Δt
2
MΔtt+
=Q-K-
Δt
2
M
t
-
Δt
2
M-
t-Δt
(4)
在给定初始条件和一定的起步计算方法后就可以利用式
(4)求解各个离散时间点的位移值。
2 有限元模型的建立
滚动轴承内部的工作情况很复杂。一般情况下,与轴承配套的主机主轴作匀速转动,轴承的内圈随之作匀速转动,在力的作用下,滚动体在轴承内部一方面产生绕轴承轴线的公转运动,同时产生绕自身轴线的自转运动。由于滚动体的受载随位置的不同而不同,使得原本应该随套圈作匀速转动的滚动体公转运动变得随位置的变化而产生了变化。滚动轴承工作时,滚动体分为受载区滚动体和非受载区滚动体。受载区内的滚动体运动受控于滚道,同步于套圈的运动;非受载区内,由于存在着阻尼,滚动体的运动滞后于套圈。这种运动的同步与滞后的交替,使得滚动体与保持架兜孔之间产生带有微量撞击现象的不连续作用,从而导致保持架的运动复杂化。
针对滚动轴承运动的特点,在仿真过程中的基本假设如下:
(1)由于轴承的倒角和边棱对轴承内部应力分布的影响很
1 基本控制方程
ANSYS/LS-DYNA是以非线性动力分析为主,应用显式时
间积分的大型有限元分析程序。在动力学有限元分析中,系统
的求解方程为:
(t)+CMa¨
(t)+K(t)=Q(t)
(1)
(t)—式中:a¨——系统的节点加速度向量;
(t)———系统的节点速度向量;
M———质量矩阵;C———阻尼矩阵;K———刚度矩阵;
Q(t)———节点载荷向量。
小,可以忽略不计,因此实体模型不包括倒角与边棱;
(2)不考虑径向游隙和轴向游隙及油膜的影响;
(3)由于轴承的塑性变形很小,所以不考虑材料非线性,假
设轴承滚动体、内外圈均为线性弹性材料。
2.1 实体模型的离散化
在动力分析中,由于惯性力和阻尼力出现在平衡方程中,
最后得到的求解方程是常微分方程组。在ANSYS/LS-DYNA下,采用直接积分法中的中心差分格式对运动方程进行积分。在中心差分法中,加速度和速度可以用位移表示为:
3收稿日期:2007-01-15;修订日期:2007-06-11
以6003深沟球轴承为例进行分析研究,其结构尺寸见表1。采用扫掠网格划分,离散后的模型如图1所示,共有11450个单元,13170个节点。
作者简介:张乐乐(1973-),女,河北秦皇岛人,博士,副教授,研究方向:安全与故障机理仿真,非线性有限元理论与应用,CAD/CAM等。
ANSYS LS_DYNA 滚动轴承 仿真
2007年9月张乐乐,等:基于ANSYS/LS2DYNA的滚动轴承仿真与分析
63
表1 6003深沟球轴承的基本结构尺寸
内径
d/mm外径D/mm厚度B/mm球径r/mm
17
35
10
5.5
滚珠数Z/个
10
看出内圈的转速比较大,由速度矢量方向也可明显地看到内圈
为逆时针转动,滚动体的转速较小,在外力的作用下,滚动体和内外圈之间存在摩擦力,带动滚动体同向转动。在图3~图7中运用公式(5)计算出的μC也是不可见的,且不是唯一值,根据接触段的不同而不同,由程序在每一步计算得到。
图1 网格划分后的滚动轴承实体模型
2.2 材料模型参数的确定
6003深沟球轴承的内外圈及滚动体的材料均为GCr15钢图2=.6s保持架为冷轧钢板,材料模型。相关参数选择如表2。