輪轂變形對星輪減速器嚙合特性的影響研究*

2020-11-04 01:29:40夏建芳劉董洋何啟成高放軍

機電工程 2020年10期

夏建芳,劉董洋,何啟成,高放軍

(1.中南大學機電工程學院,湖南長沙 410012;2.湖南機械科學研究院有限公司,湖南長沙 410007)

0 引言

星輪傳動是一種少齒內差、內平動的傳動系統。相比其他少齒差、外平動傳動系統，如三環減速器等，星輪傳動具有體積小、重量輕、傳動比大、承載能力大等優點，尤其適用于低速重載工況，因此在冶金、礦山、電力等行業中得到了廣泛的應用。

當前，國內外學者對星輪傳動這一類少齒差傳動系統均有了一定的研究。張俊等[1-2]根據子結構法及變形協調方程，對少齒差星輪型減速器進行了彈性靜力學及動力學分析；劉斌彬等[3]提出了多齒嚙合下，少齒差行星齒輪的齒根彎曲應力計算方法，建立了單齒有限元模型，并進行了修正；馮曉寧等[4]利用SolidWorks Simulation，計算了少齒差行星齒輪的齒根彎曲應力，并以此為依據，進一步研究了其超載能力；朱才朝等[5]研究了三環減速器的實際接觸齒對數；楊江兵等[6]對少齒差行星齒輪減速器進行了動態接觸仿真分析，得到了其時變嚙合剛度曲線及重合度，但其在建模過程中未考慮輪轂情況；LI[7]研究了少齒差傳動接觸問題的數值計算方法，建立了理論模型及有限元模型。另外，國內外其他一些學者[8-10]也同樣對各種少齒差傳動系統都進行了研究。

但以上研究均只考慮了少齒差傳動中的輪齒部分，且這類研究對輪轂部分大多是作剛性體來處理的。而在實際的星輪減速器中，行星齒輪輪轂通過傳動曲軸傳遞扭矩，且兩孔間存在薄壁，因此，在剛性不足的情況下，輪轂容易發生變形，必然會對齒輪嚙合造成不良影響。

在SolidWorks中，筆者對星輪傳動系統進行三維建模及虛擬裝配，在此基礎上運用ABAQUS，建立星輪傳動系統多組孔徑組合方案的有限元模型，對其進行靜力學動態仿真。

1 行星輪輸出型星輪減速器簡介

對于行星齒輪輸出型星輪減速器，在其運行過程中，其行星齒輪既繞輸入中心軸線作平動，又繞自身幾何中心作轉動。

星輪減速器機構簡圖如圖1所示。

圖1 星輪減速器機構簡圖1—輸入曲軸;2—右支撐盤;3—直連軸;4—左支撐盤;5—左星輪;6—傳動曲軸;7—機架;8—右星輪

圖1中，當輸入曲軸1順時針轉動時，左星輪5及右星輪8作平動；同時，因與固定內齒輪8的嚙合限制繞自身幾何中心線作逆時針自轉，各傳動曲軸6受星輪姿態變化的影響，繞減速器中心作逆時針轉動，從而帶動支撐盤4及輸出軸作逆時針轉動。在這一過程中，傳動曲軸的偏心方向始終保持與輸入軸一致。

行星輪的這種運動特性使其繞輸入軸線做平面復合周向運動，結果使得行星輪周期性地通過油池，而保證軸承、輪齒可以得到良好潤滑；但行星輪輪轂的形變，對行星輪和內齒輪之間的嚙合傳動將產生不良影響。其影響程度目前還沒有理論計算方法，大多還是基于有限元的數值仿真方法。

以某型星輪減速器為例，其齒輪設計參數如表1所示。

表1 齒輪設計參數

經簡化強度計算，及齒頂干涉檢查，該參數均滿足設計要求。

2 仿真模型

行星輪輸出型星輪傳動系統的行星齒輪輪轂含有6個周向均布的傳動曲軸孔及6個周向均布的直連軸孔，傳動曲軸孔的分布直徑為φ660 mm，直連軸孔的分布直徑為φ720 mm。在不改變φ660 mm、φ720 mm分布直徑的情況下，對傳動曲軸孔與直連軸孔的尺寸進行適度縮小。

傳動曲軸的直徑為φ110 mm，在原始設計方案中選取的是型號為23222CC/W33型調心滾子軸承，軸承外徑為φ200 mm。

在不改變傳動曲軸尺寸的情況下，可用的軸承型號如表2所示。

表2 傳動曲軸的可選軸承型號

其中，24022CC/W33型軸承寬度與原始方案軸承差異過大，24122CC/W33型軸承的寬度適宜。

此星輪減速器的額定負載為7.585×107N·mm，傳動曲軸的轉速為157.72 r/min，傳動曲軸孔的分布直徑為660 mm.

則單個軸承的載荷為：

(1)

式中：F—單個軸承的載荷；T—額定負載；D—傳動曲軸孔的分布圓直徑。

通過對式(1)計算，結果可得F=38 308.1 N。

傳動曲軸軸承僅承受徑向載荷，軸向載荷接近為0。根據其主要的應用場合，選取工況系數fd=1.8。因此，其當量動載荷為：

P=fd(XFr+YFa)

(2)

式中：P—當量動載荷；fd—工況系數；X—徑向動載荷系數；Fr—軸承的徑向載荷；Y—軸向動載荷系數；Fa—軸承的軸向載荷。

求解式(2)可得，P=68 954.4 N。

查《機械設計手冊》可得，24122CC/W33型調心滾子軸承的基本額定動載荷為4.58×105N，壽命指數ε=10/3，則軸承的基本額定壽命為：

(3)

式中：L10—軸承的額定壽命；n—轉速；Cr—軸承的基本額定動載荷；P—軸承的當量動載荷；ε—壽命指數。

求解式(3)可得：L10=58 206 h。

經尺寸檢查和軸承壽命初步校核可知，24122CC/W33型調心滾子軸承達到要求。因此，在孔徑方案中，筆者選取傳動曲軸孔徑為φ180 mm，變動幅度為10%.

直連軸孔的作用是使連接左、右支撐盤的直連軸穿過行星輪，且不與行星輪接觸，無配合關系；初始尺寸為φ140 mm，在孔徑方案中，選取直連軸孔尺寸逐步縮小10%。

綜上所述，僅縮小直連軸孔的孔徑方案如表3所示。

表3 僅縮小直連軸孔方案的孔徑尺寸

僅縮小傳動曲軸軸承孔的方案如表4所示。

表4 僅縮小傳動曲軸軸承孔方案的孔徑尺寸

2.1 模型簡化

實際情況下，行星輪輸出型星輪減速器傳動輪系零件多，行星輪與內齒輪齒數差小，且在載荷作用下多齒嚙合，6根傳動曲軸與行星輪裝配，輸入主軸(曲軸)與兩個行星輪裝配，這些轉配關系與多齒接觸導致邊界條件極為復雜，屬于高度非線性問題。

若直接按照真實結構及其裝配關系建立有限元模型，會導致計算規模巨大，而使計算無法進行，且過度的非線性將導致計算過程難以收斂。

因此，在滿足工程精度要求的前提下，筆者對模型進行如下簡化：

(1)行星輪輸出型星輪傳動系統的兩相行星齒輪結構相同，傳遞的轉矩相等，且與機座固聯的內齒輪具有足夠的剛度。為降低計算規模，筆者建立一對齒輪嚙合模型，負載取額定負載的一半；

(2)將傳動曲軸視為剛體，忽略傳動曲軸與行星輪之間軸承變形，將軸承結構與傳動曲軸視為一個幾何體，建立解析剛體模型；

(3)忽略嚙合齒對的齒面摩擦力；

(4)將輸入曲軸轉化為Beam-Hinge連接器，以模擬主軸剛性和轉動。

綜上所述，該模型的零件包括內齒輪齒圈、行星齒輪(含輪轂)、傳動曲軸(解析剛體)、輸入曲軸(連接器)和支撐盤(連接器)。

2.2 有限元模型

根據真實結構及上述模型簡化方法，在三維建模軟件SolidWorks中，筆者建立內齒輪齒圈、行星齒輪齒圈、行星輪輪轂等的三維實體模型，并將之導入ABAQUS中；在ABAQUS中，再建立傳動曲軸解析剛體模型(實體殼)，并在ABAQUS中，將這些實體按原有裝配關系予以裝配。

輪齒部分劃分六面體結構化網格，行星齒輪及內齒輪嚙合區域的40個輪齒網格作細化處理。以行星齒輪的輪齒部分為例，在劃分網格時通過切分、布種將齒面接觸區域和齒根彎曲區域的網格進行加密處理。加密區域網格尺寸為1 mm，平均長寬比為2.12，形狀合格；細化的40個輪齒共劃分328 000個單元。經試算，將網格進一步細化后所得結果與此網格模型結果相差在3%以內，因此，該網格可以滿足計算精度，同時控制計算時間。

輪轂部分劃分六面體掃掠網格，全局網格尺寸為3 mm。整體網格模型和局部放大展示如圖2所示(模型單元總數為968 756個)。

圖2 整體網格模型

內齒輪和行星齒輪的材料均為42CrMo，在ABAQUS中定義彈性模量為206 000 MPa，泊松比為0.29。

此處選用靜力通用分析步。因為存在大變形和接觸，勾選分析步設置中的幾何非線性。在初始分析步中，定義全局通用接觸，在解析剛體軸外表面與輪轂軸孔圓柱面之間，定義表面-表面接觸。

輸入主軸驅動行星輪以L型Beam-Hinge組合連接器來定義，其中，Beam連接器的輸入端位于輸入曲軸中心線上，Hinge連接器一端與Beam輸出端相連，另一端與行星輪中心孔耦合點相連。傳動曲軸驅動支撐盤的定義方法類似，其Beam連接器的輸入端位于剛性曲軸的中心線上，Hinge連接器位于傳動曲軸的回轉中心線上，其一端與Beam輸出端相連，另一端與支撐盤的回轉中心點(即施加載荷的參考點)建立Beam連接器。

分析過程總共定義5個分析步(含初始步)，依次實現嚙合接觸及傳動曲軸與配合孔接觸(內齒轉動0.001 rad)、預加初始載荷105N·mm、施加滿轉矩載荷7.58 5107N·mm、滿載荷轉動預定角度1 rad。