李超 陳觀慈 侯光輝 王存珠 溫戈

摘要：本文采用數(shù)學(xué)建模的方法，引入四點(diǎn)接觸球軸承套圈剛度矩陣，以此考慮結(jié)構(gòu)剛度對四點(diǎn)接觸球軸承內(nèi)部載荷分布的影響。建立了四點(diǎn)接觸球軸承的非線性位移和力學(xué)平衡方程，該方程可采用牛頓拉普遜方法迭代求解。以風(fēng)力發(fā)電機(jī)組偏航變漿用四點(diǎn)接觸球軸承為例進(jìn)行了計(jì)算分析，驗(yàn)證了本文建立數(shù)學(xué)模型的計(jì)算快速性和精確性。

Abstract： In this paper， the stiffness matrix of four point contact ball bearing ring is introduced by mathematical modeling method， and the influence of structural stiffness on the load distribution of slewing bearing is considered. The nonlinear displacement and mechanical equilibrium equations of four point contact ball bearings are established， and the equations can be calculated by Newton-Raphson iteration procedure. Taking the four point contact ball bearing as an example， the calculation and analysis of the yaw and pitch bearing of the wind turbine are carried out， and the calculation speed and accuracy of the mathematical model established in this paper are verified.

關(guān)鍵詞：四點(diǎn)接觸球軸承；剛度矩陣；有限元；力學(xué)模型

Key words： Four point contact ball bearing；Stiffness matrix；Finite element；Mechanical model

中圖分類號：TH133.3 文獻(xiàn)標(biāo)識碼：A 文章編號：1006-4311（2018）19-0255-03

0 引言

四點(diǎn)接觸球軸承套圈尺寸較大，剛性較低，軸承套圈剛性對其內(nèi)部接觸載荷分布影響較大，因此進(jìn)行力學(xué)分析時其套圈剛性成為一個不可忽略的因素。目前國內(nèi)外已經(jīng)研究了關(guān)于四點(diǎn)接觸球軸承力學(xué)建模方法和接觸載荷分析的相關(guān)概念，但主要是針對不考慮結(jié)構(gòu)剛度的剛性模型[1-3]或有限元模型[4-5]，而針對考慮結(jié)構(gòu)剛度的四點(diǎn)接觸球軸承力學(xué)模型研究較少。本文引入四點(diǎn)接觸球軸承套圈剛度矩陣，以此考慮結(jié)構(gòu)剛度對其內(nèi)部載荷分布的影響，建立其非線性位移和力學(xué)平衡方程組，采用牛頓拉普遜方法求解此方程組。論文建立的數(shù)學(xué)模型兼顧了計(jì)算精確性和快速性，為四點(diǎn)接觸球軸承相關(guān)參數(shù)的優(yōu)化奠定了基礎(chǔ)。

1 考慮套圈結(jié)構(gòu)剛度的四點(diǎn)接觸球軸承力學(xué)模型

1.1 四點(diǎn)接觸球軸承力學(xué)系統(tǒng)

圖1為風(fēng)力發(fā)電機(jī)組用雙列四點(diǎn)接觸球軸承的截面圖。圖中：α0為軸承接觸角；Dw為滾動體直徑；g0為軸承游隙；Ri、Ro為內(nèi)外圈滾道溝曲率半徑。

前述雙列四點(diǎn)接觸球軸承工作時的受力和變形可用三個彈簧串聯(lián)表達(dá)，軸承外圈下端面固定，內(nèi)圈承受外部載荷。變量Kinn、Kout代表內(nèi)外圈的剛度矩陣，Kcontact是滾動體和滾道間的接觸剛度，對內(nèi)圈施加一個外載荷F，采用迭代程序獲得平衡位置，內(nèi)外圈彈性變形δinn和δout由接觸剛度和套圈結(jié)構(gòu)剛度共同決定。

圖2為四點(diǎn)接觸球軸承受載后的幾何模型，內(nèi)圈滾動體曲率中心沿x、y、z三個方向的偏移量分別為：u、v、w，繞x、y軸偏轉(zhuǎn)的角度分別為：φx、φy，這些運(yùn)動可由坐標(biāo)變換矩陣T給出。

其中Δi為軸承受載后的滾道溝心距。

如果δcontact<0，則滾動體與滾道沒有接觸，即該位置處的滾動體不承受載荷，反之承受載荷。每個滾動體對角線的接觸力使用赫茲接觸理論計(jì)算。

1.2 平衡方程與求解

四點(diǎn)接觸球軸承單個滾動體與內(nèi)外圈滾道的接觸載荷為一對平衡力，且滾動體接觸載荷與套圈外載荷為平衡力，根據(jù)此原理即可建立套圈滾道接觸載荷與外載荷的平衡方程，如式（6）。整個系統(tǒng)的力學(xué)模型可用圖3表示。

其中Fx，F(xiàn)y，F(xiàn)z，Mx，My表示軸承承受的外載荷；φi表示第i個滾動體的位置角；αi表示軸承接觸角；Q1i、Q2i表示第一排滾動體在位置角φi處兩個接觸面法線方向的接觸載荷，Q3i、Q4i表示第二排滾動體在位置角φi處兩個接觸面法線方向的接觸載荷。

滾動體和滾道接觸載荷與軸承套圈的彈性變形相互影響，是非線性函數(shù)，因此要用圖4中牛頓拉普遜迭代程序計(jì)算套圈的變形和接觸載荷。

2 算例與結(jié)果分析

以某公司提供的風(fēng)力發(fā)電機(jī)組偏航變漿用四點(diǎn)接觸球軸承為例進(jìn)行了計(jì)算分析，該軸承的基本結(jié)構(gòu)參數(shù)和極限載荷見表2。

在相同外載荷下分別繪出四點(diǎn)接觸球軸承在剛性模型、柔性模型、有限元模型下接觸載荷和接觸角分布特點(diǎn)，如圖5和圖6所示。

從圖5可以看出，兩溝道滾珠接觸載荷分布趨勢基本一致，且每列溝道滾動體所受載荷分布不均，在圓周上呈現(xiàn)近似正弦分布。兩溝道與滾動體一半因分離而處于不受載狀態(tài)，另一半因壓緊而處于接觸狀態(tài)，大多數(shù)是兩點(diǎn)受力，少部分是四點(diǎn)受力，且四點(diǎn)接觸的滾動體承受的載荷比較小。軸承第一列滾動體對稱位置承受的載荷相差較大，第二列基本一致。

從圖6可以看出，兩溝道滾珠接觸角分布趨勢基本一致，且每列溝道滾動體接觸角不同，與滾動體的接觸載荷分布規(guī)律一樣，也成正弦分布。處于四點(diǎn)接觸的滾動體對應(yīng)的接觸角相比原始接觸角變小了，接觸角在0°到150°承載區(qū)域內(nèi)的變化幅度相對較小，在150°到360°承載區(qū)域內(nèi)變化幅度相對較大。

3 結(jié)論

本文以雙排四點(diǎn)接觸球軸承為例，建立了考慮套圈結(jié)構(gòu)剛度的大型轉(zhuǎn)盤軸承靜力學(xué)模型，建立了簡化滾動體的有限元模型。且柔性模型所得結(jié)果與有限元模型結(jié)果基本吻合，其計(jì)算時間成本優(yōu)于有限元模型，準(zhǔn)確度優(yōu)于剛性模型。該模型可擴(kuò)展到考慮支撐結(jié)構(gòu)剛度的四點(diǎn)接觸球軸承，分析其結(jié)構(gòu)剛度對其接觸載荷和接觸角分布的影響，并進(jìn)行相關(guān)結(jié)構(gòu)參數(shù)的優(yōu)化設(shè)計(jì)。

參考文獻(xiàn)：

[1]高學(xué)海，黃筱調(diào)，王華，陳捷.雙排四點(diǎn)接觸球轉(zhuǎn)盤軸承滾道接觸壓力分布[J].南京工業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2011（01）：80-83.

[2]徐紹仁，楊友中，白宇光，王燕霜，鄧四二.特大型負(fù)游隙四點(diǎn)接觸球軸承擬靜力學(xué)分析[J].軸承，2009（09）：1-3，7.

[3]陳國楨，徐紹仁，楊友中，王燕霜，鄧四二.特大型負(fù)游隙四點(diǎn)接觸球軸承接觸應(yīng)力分析[J].機(jī)械傳動，2009（05）：83-85，91，112.

[4]李云峰，劉彩霞，陳紅濤.雙排四點(diǎn)接觸球轉(zhuǎn)盤軸承的有限元分析[J].哈爾濱軸承，2015（02）：3-6，9.

[5]王存珠，陳觀慈，李肖杰，溫戈.螺栓及其預(yù)緊力對大型轉(zhuǎn)盤軸承接觸載荷分布的影響[J].軸承，2016（01）：24-27，33.