CN114580093A - 一种计算滚道局部故障球轴承动力学响应的方法 - Google Patents
一种计算滚道局部故障球轴承动力学响应的方法 Download PDFInfo
- Publication number
- CN114580093A CN114580093A CN202111676302.7A CN202111676302A CN114580093A CN 114580093 A CN114580093 A CN 114580093A CN 202111676302 A CN202111676302 A CN 202111676302A CN 114580093 A CN114580093 A CN 114580093A
- Authority
- CN
- China
- Prior art keywords
- fault
- raceway
- ball bearing
- local
- calculating
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Rolling Contact Bearings (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提出了一种计算滚道局部故障球轴承动力学响应的方法,获取滚动轴承的基本参数和故障参数;计算故障跨度角和局部故障引起的额外间隙;计算滚动体中心与局部故障中心的角度差,结合滚动体是否位于故障区,确定故障发生后滚动体与滚道的变形量;建立滚道局部故障球轴承动力学模型;采用4阶龙格‑库塔法求解滚道局部故障球轴承动力学方程,得到球轴承在垂直方向和水平方向的径向位移,完成滚道局部故障球轴承动力学响应的计算。本发明能够模拟多种局部故障状况,系统、高效地计算滚道局部故障球轴承的动力学响应,对球轴承的故障诊断具有重要价值。
Description
技术领域
本发明涉及故障球轴承动力学响应求解领域,具体涉及一种计算滚道局部故障球轴承动力学响应的方法。
背景技术
球轴承包含轴承变刚度、赫兹接触和轴承间隙等多个非线因素,会导致球轴承出现复杂的动力学行为(张智勇,2015)。球轴承的点蚀、剥落、裂纹和磨损类故障是最为常见的滚道类故障形式,这类故障可简化为局部类故障(尉询楷,2014)。球轴承发生局部故障后,非线性因素增加,导致球轴承的动力学响应更加复杂。Tandon和Choudhury认为轴承局部故障产生的冲击有持续时间,他们选用矩形、三角形和半正弦形三种常见脉冲表征故障产生的冲击,最终冲击力等于冲击函数乘以故障处接触力(TandonN et al.1997)。Tadina和Boltezar对故障轴承系统建模时,将多体动力学方法与有限元方法相结合得到了系统的控制方程,并考虑了不同边界性质对接触刚度的影响(Tadina M,Boltezar M.2011)。以上对球轴承局部故障的研究中,前者通过使用不同的冲击力序列对局部故障进行模拟,其缺点是无法直接得到冲击力的大小,不能对局部故障进行定量描述;后者使用有限元方法模拟局部故障,其缺点为无法模拟故障球轴承的动态行为。本发明克服了以上缺点,利用局部故障引起的额外间隙描述局部故障产生的冲击,并能够模拟多种局部故障状况(包括故障类型、故障尺寸和故障位置),并且系统、高效地计算滚道局部故障球轴承动力学响应。
发明内容
本发明的目的在于提供一种计算滚道局部故障球轴承动力学响应的方法。
完成本发明目的的技术解决方案为:一种计算滚道局部故障球轴承动力学响应的方法,包括以下步骤:
步骤1、获取滚动轴承的基本参数和故障参数;
步骤2、计算故障跨度角和局部故障引起的额外间隙;
步骤3、计算滚动体中心与局部故障中心的角度差,结合滚动体是否位于故障区,确定故障发生后滚动体与滚道的变形量;
步骤4、代入步骤1滚动轴承的基本参数和故障参数以及步骤2故障跨度角和局部故障引起的额外间隙,并结合滚动体是否位于故障区,建立滚道局部故障球轴承动力学模型;
步骤5、采用4阶龙格-库塔法求解滚道局部故障球轴承动力学方程,得到球轴承在垂直方向和水平方向的径向位移,完成滚道局部故障球轴承动力学响应的计算。
进一步的,步骤1,获取滚动轴承的基本参数和故障参数,其中:
所述基本参数包括几何尺寸、滚珠个数、接触刚度、轴承间隙和等效阻尼;
所述故障参数包括故障类型、故障尺寸和故障位置;
所述故障类型包括内圈滚道故障和外圈滚道故障。
进一步的,步骤2,计算故障跨度角和局部故障引起的额外间隙,具体方法为:
步骤2.1、计算故障跨度角θe:
步骤2.3、计算局部故障引起的额外间隙Δh:
式(1)至(3)中,d为局部故障尺寸;di、do分别为轴承内圈、外圈滚道直径;D为滚动体直径。
进一步的,步骤3,计算滚动体中心与局部故障中心的角度差,结合滚动体是否位于故障区,确定故障发生后滚动体与滚道的变形量,具体方法为:
步骤3.1、计算每一个滚动体中心与局部故障中心的角度差θbd:
步骤3.2、如果滚动体位于故障区,则滚动体与正常滚道的变形量需要减去局部故障引起的额外间隙,否则不变:
进一步的,代入步骤1滚动轴承的基本参数和故障参数以及步骤2故障跨度角和局部故障引起的额外间隙,并结合滚动体是否位于故障区,建立滚道局部故障球轴承动力学模型,具体方法为:
且
θi=2π(i-1)/Nb+ωct (9)
δi=xcosθi+ysinθi-δ0 (10)
式(6)至(11)中,x、y为球轴承在垂直方向和水平方向的径向位移;Fx、Fy为球轴承在x、y方向轴承反力;m为球轴承等效质量;c为等效阻尼;W为球轴承重力;θi为第i个滚动体的瞬时角位置;2δ0为轴承径向工作游隙;Fi *为滚动体与滚道间的赫兹接触力,K为赫兹接触刚度;G(·)是表示接触状况的Heaviside函数,描述滚动体与滚道接触情形,发生接触G(·)值为1,失去接触其值为0;Nb为滚动体个数。
一种计算滚道局部故障球轴承动力学响应的系统,基于所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时,基于所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时,基于所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
本发明与现有技术相比,其显著优点为:1)利用滚动体中心在故障区下降的高度表示局部故障引起的额外间隙,克服了人为主观规定额外间隙给模型带来的不利影响;2)包含了局部故障的(故障类别、故障尺寸和故障位置)等多种故障参数,对开展不同故障情况下球轴承动力学响应的研究带来了很大的便利;3)采用4阶龙格-库塔法求解滚道局部故障球轴承的动力学方程,系统、高效地计算滚道局部故障球轴承的动力学响应,不仅可以计算稳态响应,还可以计算瞬态响应。
附图说明
图1是本发明实施例轴承内圈、外圈滚道局部故障示意图。
图2是本发明实施例局部故障引起的额外间隙几何变化示意图,(a)是外圈故障示意图,(b)是内圈故障示意图。
图3是本发明实施例球轴承模型图,(a)是球轴承剖面图,(b)是动力学原理图。
图4是本发明实施例外圈局部故障系统的加速度响应图,(a)是水平方向加速度响应图,(b)是垂直方向加速度响应图。
图5是本发明实施例内圈局部故障系统的加速度响应图,(a)是水平方向加速度响应图,(b)是垂直方向加速度响应图。
图6是本发明实施例改变外圈故障位置角后系统的加速度响应图,(a)是水平方向加速度响应图,(b)是垂直方向加速度响应图。
图7是本发明计算滚道局部故障球轴承动力学响应方法的流程图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
如图7所示,一种计算滚道局部故障球轴承动力学响应的方法,包括以下步骤:
步骤1,输入滚动轴承的基本参数(几何尺寸、滚珠个数、接触刚度、轴承间隙、等效阻尼)和故障参数(内圈或外圈故障、故障尺寸、故障位置);
步骤2,根据步骤1提供的滚动轴承的基本参数和故障参数,计算故障跨度角和局部故障引起的额外间隙,具体方法为:
步骤2.1,计算故障跨度角θe:
步骤2.3,计算局部故障引起的额外间隙Δh:
式(1)至(3)中,d为局部故障尺寸;di、do分别为轴承内圈、外圈滚道直径;D为滚动体直径。
步骤3中,计算每一个滚动体中心与局部故障中心的角度差,并判断滚动体是否位于故障区,如果位于故障区则滚动体与滚道的变形量减去额外间隙,否则不变,具体方法为:
步骤3.1,计算每一个滚动体中心与局部故障中心的角度差θbd:
步骤3.2,判断滚动体是否位于故障区,如果滚动体位于故障区,则滚动体与正常滚道的变形量需要减去局部故障引起的额外间隙,否则不变:
步骤4,代入步骤1滚动轴承的基本参数和故障参数以及步骤2故障跨度角和局部故障引起的额外间隙,并结合步骤3中判断滚动体是否位于故障区,建立滚道局部故障球轴承动力学模型:
且
θi=2π(i-1)/Nb+ωct (9)
δi=xcosθi+ysinθi-δ0 (10)
式(6)至(11)中,x、y为球轴承在垂直方向和水平方向的径向位移;Fx、Fy为球轴承在x、y方向轴承反力;m为球轴承等效质量;c为等效阻尼;W为球轴承重力;θi为第i个滚动体的瞬时角位置;2δ0为轴承径向工作游隙;Fi *为滚动体与滚道间的赫兹接触力,K为赫兹接触刚度;G(·)是表示接触状况的Heaviside函数,描述滚动体与滚道接触情形(发生接触G(·)值为1,失去接触其值为0);Nb为滚动体个数。
步骤5,采用4阶龙格-库塔法求解步骤4中滚道局部故障球轴承动力学方程,得到球轴承在垂直方向和水平方向的径向位移,完成滚道局部故障球轴承动力学响应的计算。
本发明还提出一种计算滚道局部故障球轴承动力学响应的系统,基于所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时,基于所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时,基于所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
实施例
为了验证本发明方法的有效性,进行如下仿真。给定深沟球轴承6308几何参数如表1所示,局部故障参数如表2所示。
表1深沟球轴承6308几何参数
表2局部故障参数
对于如图1所示的局部故障分布示意图,计算外圈、内圈故障跨度角分别为0.0175和0.02,根据图2滚动体中心与故障区的几何位置关系,计算外圈、内圈滚道局部故障的故障角分别为0.0933和0.0666,计算外圈、内圈滚道局部故障引起的额外间隙分别为7.9μm和4.1μm。计算滚动体中心与局部故障中心的角度差,判断滚动体是否位于故障区。图3为球轴承剖面图和动力学原理图。代入滚动轴承的基本参数和故障参数以及故障跨度角和局部故障引起的额外间隙,采用4阶龙格-库塔法求解滚道局部故障球轴承动力学方程,得到滚道局部故障球轴承的动力学响应,其中图4为外圈故障动力学响应,加速度呈现规则的冲击;图5为内圈故障动力学响应,加速度幅值随时间变化。改变外圈故障位置角为π/60,可以计算得到此故障情况下系统的动力学响应,如图6所示,垂直方向加速度增强,水平方向加速度减弱。
综上所述,本发明克服了人为主观规定局部故障引起的额外间隙给模型带来的不利影响,并且包含了局部故障的(故障类别、故障尺寸和故障位置)等多种故障参数,能够系统、高效地计算滚道局部故障球轴承的稳态响应和瞬态响应,对球轴承的故障诊断具有重要价值。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (8)
1.一种计算滚道局部故障球轴承动力学响应的方法,其特征在于,包括以下步骤:
步骤1、获取滚动轴承的基本参数和故障参数;
步骤2、计算故障跨度角和局部故障引起的额外间隙;
步骤3、计算滚动体中心与局部故障中心的角度差,结合滚动体是否位于故障区,确定故障发生后滚动体与滚道的变形量;
步骤4、代入步骤1滚动轴承的基本参数和故障参数以及步骤2故障跨度角和局部故障引起的额外间隙,并结合滚动体是否位于故障区,建立滚道局部故障球轴承动力学模型;
步骤5、采用4阶龙格-库塔法求解滚道局部故障球轴承动力学方程,得到球轴承在垂直方向和水平方向的径向位移,完成滚道局部故障球轴承动力学响应的计算。
2.根据权利要求1所述的计算滚道局部故障球轴承动力学响应的方法,其特征在于,步骤1,获取滚动轴承的基本参数和故障参数,其中:
所述基本参数包括几何尺寸、滚珠个数、接触刚度、轴承间隙和等效阻尼;
所述故障参数包括故障类型、故障尺寸和故障位置;
所述故障类型包括内圈滚道故障和外圈滚道故障。
5.根据权利要求1所述的计算滚道局部故障球轴承动力学响应的方法,其特征在于,步骤4,代入步骤1滚动轴承的基本参数和故障参数以及步骤2故障跨度角和局部故障引起的额外间隙,并结合滚动体是否位于故障区,建立滚道局部故障球轴承动力学模型,具体方法为:
且
θi=2π(i-1)/Nb+ωct (9)
δi=xcosθi+ysinθi-δ0 (10)
式(6)至(11)中,x、y为球轴承在垂直方向和水平方向的径向位移;Fx、Fy为球轴承在x、y方向轴承反力;m为球轴承等效质量;c为等效阻尼;W为球轴承重力;θi为第i个滚动体的瞬时角位置;2δ0为轴承径向工作游隙;Fi *为滚动体与滚道间的赫兹接触力,K为赫兹接触刚度;G(·)是表示接触状况的Heaviside函数,描述滚动体与滚道接触情形,发生接触G(·)值为1,失去接触其值为0;Nb为滚动体个数。
6.一种计算滚道局部故障球轴承动力学响应的系统,其特征在于,基于权利要求1-5任一项所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
7.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时,基于权利要求1-5任一项所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
8.一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时,基于权利要求1-5任一项所述的计算滚道局部故障球轴承动力学响应的方法,完成滚道局部故障球轴承动力学响应的计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111676302.7A CN114580093A (zh) | 2021-12-31 | 2021-12-31 | 一种计算滚道局部故障球轴承动力学响应的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111676302.7A CN114580093A (zh) | 2021-12-31 | 2021-12-31 | 一种计算滚道局部故障球轴承动力学响应的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114580093A true CN114580093A (zh) | 2022-06-03 |
Family
ID=81769880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111676302.7A Pending CN114580093A (zh) | 2021-12-31 | 2021-12-31 | 一种计算滚道局部故障球轴承动力学响应的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114580093A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115270342A (zh) * | 2022-08-01 | 2022-11-01 | 安徽大学 | 一种剥落故障深沟球轴承多自由度动力学建模方法 |
-
2021
- 2021-12-31 CN CN202111676302.7A patent/CN114580093A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115270342A (zh) * | 2022-08-01 | 2022-11-01 | 安徽大学 | 一种剥落故障深沟球轴承多自由度动力学建模方法 |
CN115270342B (zh) * | 2022-08-01 | 2023-08-15 | 安徽大学 | 一种剥落故障深沟球轴承多自由度动力学建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cui et al. | HVSRMS localization formula and localization law: Localization diagnosis of a ball bearing outer ring fault | |
Liu | A dynamic modelling method of a rotor-roller bearing-housing system with a localized fault including the additional excitation zone | |
Liu et al. | A combined acoustic and dynamic model of a defective ball bearing | |
Li et al. | Analysis of varying contact angles and load distributions in defective angular contact ball bearing | |
Qin et al. | A fault dynamic model of high-speed angular contact ball bearings | |
Rafsanjani et al. | Nonlinear dynamic modeling of surface defects in rolling element bearing systems | |
Niu et al. | A systematic study of ball passing frequencies based on dynamic modeling of rolling ball bearings with localized surface defects | |
Behzad et al. | A new model for estimating vibrations generated in the defective rolling element bearings | |
Kong et al. | A vibration model of ball bearings with a localized defect based on the hertzian contact stress distribution | |
Qin et al. | Multiple-degree-of-freedom dynamic model of rolling bearing with a localized surface defect | |
Qin et al. | Dynamics modelling for deep groove ball bearings with local faults based on coupled and segmented displacement excitation | |
Kogan et al. | Toward a 3D dynamic model of a faulty duplex ball bearing | |
Bastami et al. | A comprehensive evaluation of the effect of defect size in rolling element bearings on the statistical features of the vibration signal | |
CN114580093A (zh) | 一种计算滚道局部故障球轴承动力学响应的方法 | |
Parmar et al. | Effect of dynamic misalignment on the vibration response, trajectory followed and defect-depth achieved by the rolling-elements in a double-row spherical rolling-element bearing | |
Zhang et al. | Dynamic modeling and analysis of rolling bearing with compound fault on raceway and rolling element | |
Kumbhar et al. | An overview of dynamic modeling of rolling-element bearings | |
Ghalamchi et al. | Modeling and dynamic analysis of spherical roller bearing with localized defects: analytical formulation to calculate defect depth and stiffness | |
Kerst et al. | A semi-analytical bearing model considering outer race flexibility for model based bearing load monitoring | |
LetUriondo et al. | Multi-body modelling of rolling element bearings and performance evaluation with localised damage | |
Houpert et al. | Bearing life calculations in rotating and oscillating applications | |
Govardhan et al. | Fault diagnosis of dynamically loaded bearing with localized defect based on defect-induced excitation | |
Karacay et al. | Vibrations of a grinding spindle supported by angular contact ball bearings | |
Jamadar et al. | An in-situ synthesized model for detection of defective roller in rolling bearings | |
CN115270342A (zh) | 一种剥落故障深沟球轴承多自由度动力学建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |