基于接觸變形量修正模型的圓錐滾子軸承剛度計算

楊 娜，劉國政

基于接觸變形量修正模型的圓錐滾子軸承剛度計算

楊 娜1，劉國政2

（1.河南工學院 車輛與交通工程學院，河南 新鄉(xiāng) 453003；2.上汽大眾汽車有限公司，上海 201805）

軸承接觸受力以及剛度計算是一個非線性力學問題，文章基于修正的Palmgren模型對圓錐滾子軸承的滾子接觸力和接觸剛度進行了計算。通過對所有滾子的接觸力、力矩以及其對軸承內(nèi)外圈造成的接觸變形的計算和疊加，推導了圓錐滾子軸承的剛度模型。為了驗證模型，建立網(wǎng)格數(shù)量多、精度較高的圓錐滾子軸承有限元模型，計算了滾子以及外圈的應力分布情況，通過對有限元計算結(jié)果接觸力和接觸變形量的提取，識別了軸承剛度并與文章所建立模型的計算結(jié)果進行了對比，結(jié)果表明，所建立模型是準確的。

圓錐滾子軸承；剛度計算；修正模型；接觸變形

軸承可以說是現(xiàn)代工業(yè)的基礎，在機械系統(tǒng)的應用非常廣泛。軸承在工作過程中需要承受較大的接觸力，其接觸剛度的計算是否準確直接影響整個軸承系統(tǒng)動力學建模的準確度。目前對于軸承接觸剛度計算主要采用兩種方式：有限元法和基于經(jīng)驗公式的數(shù)值計算法[1]。有限元法的計算精度往往隨著網(wǎng)格數(shù)量等計算規(guī)模的擴大而增加[2]，但其計算量也迅速增加，計算收斂度急劇下降[3]。數(shù)值計算方法收斂性好，計算量小，但是因為其依賴經(jīng)驗公式，存在適應性差的問題[4]。本文擬基于一種修正的接觸變形量計算模型，對圓錐滾子軸承的接觸剛度進行計算，并探究計算方法的適應性問題。

1 軸承接觸受力建模

圓錐滾子軸承中的接觸屬于線接觸，對于線接觸中的接觸變形量存在多種的經(jīng)驗公式，其中比較著名的是PALMGREN在1959年提出的Palmgren經(jīng)驗公式[5]，如式（1）所示。

式中，為接觸變形量；為接觸區(qū)的法向載荷；為接觸線長度。

Palmgren經(jīng)驗公式應用于平面線接觸計算的精度還是比較高的。但對于軸承而言，其內(nèi)外支承圈存在曲率，不同尺寸的軸承，其曲率也是不同的，這就造成了該經(jīng)驗公式在計算圓錐滾子軸承的接觸問題時可能存在巨大的誤差。因此，洛陽軸承研究所的羅繼偉對Palmgren公式進行了改進，加入了軸承滾子直徑和軸承內(nèi)、外圈直徑等參數(shù)[6]，如式（2）所示。

式中，為接觸變形量；w為滾子直徑；為接觸區(qū)的法向載荷；為接觸線長度（即滾子長度）；=w/r，r為滾道直徑。式（2）中括號內(nèi)的正負號代表接觸面曲率情況，正號代表滾子與內(nèi)圈接觸，負號代表滾子與外圈接觸。

軸承結(jié)構(gòu)中存在保持架等結(jié)構(gòu)，其與軸承之間受力較小，因此，在本文的計算中忽略保持架結(jié)構(gòu)；另外認為滾子與內(nèi)外圈沒有間隙；實際使用過程中滾子與內(nèi)外圈之間極少存在滑動，因此，認為滾子與內(nèi)外圈是純滾動的；在動力學分析中，由于滾子的質(zhì)量引起的慣性力相對于接觸力來說要小得多，因此，本文在計算和建模過程中不考慮慣性力。

圖1 圓錐滾子軸承結(jié)構(gòu)與受力圖

某圓錐滾子軸承的幾何結(jié)構(gòu)如圖1所示，圖中i和e分別是圓錐滾子與內(nèi)外圈的接觸角，是圓錐滾子的錐角，而且三者存在如式（3）所示關(guān)系。圓錐滾子受到三個外力：外圈的法向接觸力e、內(nèi)圈的法向接觸力i以及內(nèi)圈擋邊的接觸力f。其中f的方向平行于滾子的軸線，i和e的方向分別垂直于軸承的內(nèi)外圈。將上述三個力分解到垂直于滾子軸線的方向上，根據(jù)滾子受力平衡可得這個三個接觸力之間的關(guān)系如式（4）所示，可以看出滾子與內(nèi)外圈間的接觸力大小相等[7]。

我們再來分析單個滾子的受力情況，如圖2所示，滾子與內(nèi)外圈之間的剛度分別為i和e，滾子與內(nèi)外圈之間的彈性變形量分別為Ei和Ee，因此，滾子與內(nèi)外圈之間類似于兩個彈簧串聯(lián)。以軸承外圈的接觸位移為參照，將內(nèi)圈與滾子的接觸位移投影到外圈接觸力方向上，Ei在e方向上的投影為Eicos，因此，軸承內(nèi)外圈在e方向上總的相對位移如式（5）所示。

圖2 單個滾子受力分析圖

因為圓錐滾子的錐角較小，在計算過程中我們將其當作圓柱滾子處理，W為圓錐滾子的平均直徑。根據(jù)修正的接觸變形量計算公式（2），滾子與內(nèi)外圈的接觸位移分別如式（6）和式（7）所示。

將式（6）和（7）代入到式（5）中，可得軸承內(nèi)外圈總的相對位移，如式（8）所示。

融合DPI采集平臺，按照“單次分光、統(tǒng)一分流、標準化采集”統(tǒng)一標準部署融合DPI采集設備，通過并接方式，采集運營商的4張數(shù)據(jù)源網(wǎng)絡全部數(shù)據(jù)。運營商就能更快速準確的為各個政府部門提供源數(shù)據(jù)，保障我國信息和網(wǎng)絡安全，維護社會穩(wěn)定和保障經(jīng)濟發(fā)展，同時也能更好的為運營商各類分析平臺提供大數(shù)據(jù)，提升運營能力和降低成本。

利用式（8）可以求得軸承的外圈法向接觸力，如式（9）所示。

我們可以得到式（10），其中ne稱為外圈接觸的剛度系數(shù)。

2 圓錐滾子軸承剛度模型

2.1 軸向與徑向受力分析

圓錐滾子軸承的受力情況如圖3所示，在外力的作用下，軸承內(nèi)圈與外圈之間的軸向和徑向位移分別為δ和δ。取軸承內(nèi)圈的中心點為坐標系原點，內(nèi)圈的軸線為軸，正方向為從滾子大端指向小端。將軸承內(nèi)圈-滾子作為一個整體，僅考慮外圈與滾子之間的法向位移，則最頂端的滾子處，外圈的軸向和徑向位移分別為δ和δ，由幾何關(guān)系可以得到外圈的法向位移量，如式（11）所示。

式中，αe是軸承外圈的壓力角。

圖3中，是滾子的方位角，第個滾子的方位角φ=2π(-1)，根據(jù)幾何關(guān)系，第個滾子處，軸承外圈的徑向位移δcos，軸向位移，所以第個滾子處軸承外圈之間的法向位移量可以由式（12）求得

將式（12）代入到式（10）中，可以求出第個滾子與外圈的法向載荷e。如圖4所示，將軸承內(nèi)圈-滾子作為力學分析的隔離體，分別計算每個滾子法向接觸力e的軸向分力Q和徑向分力Q。再將各個滾子的軸向力和徑向力疊加，就能得到軸承整體所受到軸向力F和徑向力F，分別如式（13）和式（14）所示。

上述分析過程只考慮了軸承內(nèi)外圈之間的軸向和徑向位移。從圖4可以看出，各個滾子受到的軸向力，會對軸承中心產(chǎn)生一個附加力矩，迫使軸承內(nèi)圈產(chǎn)生轉(zhuǎn)角。根據(jù)力矩平衡，要想保持軸承受力平衡，必須給軸承內(nèi)圈提供一個和附加力矩大小相同，方向相反的力矩，如式（15）所示。下一小節(jié)中我們將詳細分析軸承所受力矩情況。

式中，m為滾子節(jié)圓半徑，即軸承內(nèi)外圈半徑的平均值。

圖4 圓錐滾子軸承受力與力矩

2.2 力矩分析

當軸承內(nèi)外圈存在相對轉(zhuǎn)角時，如圖5所示，第個滾子相對于外圈存在如式（16）所示的軸向位移。

與引起的軸向位移量δ相比，由引起的滾子徑向位移δ可以忽略[6]，因此，軸承滾子總的軸向位移量如式（17）所示。

則依據(jù)式（12）可以計算第個滾子處，軸承外圈相對于滾子的法向位移量。

考慮到徑向位移δ包括沿著和兩個方向的位移δ和δ，轉(zhuǎn)角也包括軸和繞軸兩個方向的轉(zhuǎn)角θ和θ，所以軸承外圈相對于滾子的法向位移如式（19）所示。

當δ＞0時，說明滾子發(fā)生變形，將上式代入到式子（10）中，就可以求出各個滾子與外圈的法向載荷e；但是當δ≤0時，表示軸承滾子的運動遠離軸承外圈，滾子不會發(fā)生接觸變形，此時e=0。

將每個滾子的徑向力、軸向力以及由軸向力產(chǎn)生的附加力矩進行疊加，可以得到軸承整體所受到的力和力矩，如式（21）所示。

假設軸承外圈固定不動，當軸承內(nèi)圈受到力F()和轉(zhuǎn)矩M()時，軸承內(nèi)圈會產(chǎn)生位移δ()和轉(zhuǎn)角θ()。依據(jù)上述力與位移的關(guān)系，我們可以對軸承的剛度進行計算和描述。

3 圓錐滾子軸承剛度模型驗證

為了驗證式（21）所示圓錐滾子軸承剛度模型的正確性，以某圓錐滾子軸承為例，建立軸承有限元模型，仿真分析其軸向和徑向剛度，并與軸承理論模型的計算結(jié)果作對比。

對圓錐滾子軸承幾何模型進行適當?shù)暮喕蛶缀吻謇韀8]，建立的有限元模型如圖6所示。

圖6 圓錐滾子軸承有限元模型

圖7和圖8分別是軸承受到軸向和徑向載荷時，滾子和外圈的應力云圖，可以看出滾子和外圈基本是線接觸，對于軸向載荷，每個滾子承受的載荷基本一致；軸承受到徑向力時，在力的作用線上的滾子受到的應力最大，其它滾子的受力情況關(guān)于力的作用線對稱，共有7個滾子承受載荷。

圖7 圓錐滾子軸承軸向受力100 kN

圖8 圓錐滾子軸承徑向受力100 kN

為了分析應力沿著軸承滾子母線方向的分布，沿著滾子母線方向提取滾子表面與外圈接觸區(qū)的應力，如圖9所示，由于滾子和內(nèi)外圈是有限長接觸，在滾子大端和小端邊緣處均會出現(xiàn)應力集中，滾子小端的應力大于大端。但是除了滾子的兩端，滾子母線上的應力變化不是很大，接觸應力基本平均分布。

圖9 滾子接觸應力提取路徑與分布情況

圖10和圖11為軸承剛度的有限元仿真結(jié)果與計算結(jié)果的對比圖，可以看出，整體上軸承的軸向剛度小于徑向剛度；在相同的位移下，無論是載荷還是剛度值，有限元仿真結(jié)果均稍大于理論模型，這是由于有限元建模時將滾子與軸承內(nèi)圈綁定在一起，導致兩者之間接觸面的剛度變大；軸承剛度隨著位移的增加呈現(xiàn)出明顯的非線性，尤其當位移在0.02 mm以內(nèi)時，當剛度位移量超過0.02 mm后，剛度基本呈現(xiàn)線性增加；軸承的徑向剛度大于軸向剛度。雖然非線性軸承模型與仿真結(jié)果存在一定的差異，但是兩者的變化趨勢基本一致，說明本文推導的非線性軸承理論模型是正確的，其綜合考慮了軸承內(nèi)外圈直徑對彈性變形量的影響、滾子錐角對滾子與內(nèi)外圈剛度系數(shù)的影響，以及各個滾子承受載荷的差異。

圖10 軸承軸向剛度計算與有限元仿真對比

圖11 軸承徑向剛度計算與有限元仿真對比

4 結(jié)論

本文基于修正的Palmgren經(jīng)驗公式，對軸承中單個滾子的受力情況進行了計算和分析，推導并建立了圓錐滾子軸承法向接觸力模型，定義了軸承滾子與外圈接觸剛度系數(shù)。

通過對圓錐滾子軸承每個滾子的徑向力、軸向力以及由軸向力產(chǎn)生的附加力矩進行計算和疊加，推導了圓錐滾子軸承整體軸向和徑向力，以及相應力矩的表達式，建立了圓錐滾子軸承剛度模型。

本文通過有限元方法對所建立的軸承剛度模型進行了驗證，兩者的計算結(jié)果無論從數(shù)值差還是變化趨勢相似度上，都表明本文所建立的圓錐滾子軸承剛度理論模型是正確的，探明了基于修正的Palmgren公式進行的軸承受力與剛度計算公式是可以適用于圓錐滾子軸承的，而利用本文所建立模型計算軸承受力與剛度等數(shù)據(jù)所耗費的時間要遠遠小于有限元法，證明本文研究方法可以有效改善圓錐滾子軸承的計算效率。

[1] 陳潤霖,趙少東,張夢豪,等.角接觸球軸承的時變剛度模型與承載特性分析[J].機械設計與研究,2022,38 (2):114-119.

[2] BOZCA M.Effects of Design Parameters on Static Equivalent Stress of Radial Rolling Bearings [J].Acta Polytechnica.2021,61(1):163-173.

[3] 張永奇.汽車驅(qū)動橋主減速器圓錐滾動軸承工作性能研究[D].長春:吉林大學,2012.

[4] 丁長安,張雷,周福章,等.線接觸彈性接觸變形的解析算法[J].摩擦學學報, 2001,21(2):135-138.

[5] SHUTING L.A Mathematical Model and Numeric Method for Contact Analysis of Rolling Bearings[J]. Mechanism and Machine Theory,2018,19(2):61-73.

[6] 羅繼偉,羅天宇.滾動軸承分析計算與應用[M].北京:機械工業(yè)出版社,2009.

[7] 王彥偉,羅繼偉,陳立平.圓錐滾子軸承接觸分析[J].華中科技大學學報(自然科學版),2007,35(9):19-21.

[8] 周馳,田程,丁煒琦,等.基于有限元法的準雙曲面齒輪時變嚙合特性研究[J].機械工程學報,2016,52(15): 36-43.

Stiffness Calculation of Tapered Roller Bearing Based on Modified Contact Deformation Model

YANG Na1, LIU Guozheng2

( 1.School of Vehicle and Traffic Engineering, Henan Institute of Technology, Xinxiang 453003, China;2.SAIC Volkswagen Company Limited, Shanghai 201805, China )

The calculation of bearing contact force and stiffness is a nonlinear mechanical problem. This paper calculates the roller contact force and contact stiffness of tapered roller bearing based on the modified palmgren model. Through the calculation and superposition of the contact forces and moments of all rollers and the contact deformation caused by them on the inner and outer rings of the bearing, the stiffness model of tapered roller bearings is derived in this paper. In order to verify the model, this paper establishes a finite element model of tapered roller bearing with a large number of grids and high accuracy, calculates the stress distribution of the roller and outer ring, identifies the bearing stiffness by extracting the contact force and contact deformation from the finite element calculation results, and compares the calculation results with the model established in this paper. The results show that the model established in this paper is accurate.

Tapered roller bearing; Stiffness calculation; Modified model; Contact deformation

U463.5

A

1671-7988(2023)10-119-06

10.16638/j.cnki.1671-7988.2023.010.025

楊娜（1991—），女，工程師，研究方向為汽車系統(tǒng)動力學，E-mail:yangna1216@163.com。