鄭暉,馮金,覃琨,劉曉凡,李嘉東
(中國航發商用航空發動機有限責任公司 機械系統部,上海202007)
齒輪傳動系統是機械傳動系統中最重要而且應用最廣泛的運動裝置,星型齒輪傳動是航空發動機上常用的一種功率分流齒輪系傳動,其結構如圖1所示。由于星型齒輪傳動結構采用定軸輪系,其系統剛度、強度和工作可靠性都較高,一般適用于航空減速器。星型輪系的傳動路徑是從太陽輪經行星輪到內齒圈,依靠多個行星輪將輪系的動力分流傳動,行星架則作為固定機構,將齒輪嚙合產生的轉矩傳遞到傳扭支架上。
當星型輪系由于加工誤差和裝配誤差產生不均載現象時,會導致系統應力分布不均,從而產生振動、噪聲、降低齒輪壽命等影響,降低了系統的性能和效率。對星型齒輪傳動系統動態不均載特性進行計算分析,有助于進一步地了解系統的傳動動態特性和均載特性,提高齒輪箱穩定性、壽命和效率。
目前國內很少有對行星架應變片的粘貼和標定方法方面的研究,大多數研究都集中在使用有限元軟件對風電齒輪箱和大型機械設備使用的減速器行星架進行分析。通過理論的計算與試驗的數據,使用軟件模擬行星架在齒輪傳動系統在實際運轉過程中的受力情況和邊界條件,得出行星架的應力和應變分布,從而對行星架進行優化處理[1-3],但是國內并沒有對行星架均載試驗過程中均載系數進行計算和監測的方法。

圖1 星型齒輪傳動結構圖
本文通過有限元分析與應變片粘貼的方式對行星架每個軸孔的受力情況進行分析與驗證。首先通過ANSYS軟件對行星架進行有限元分析,確定應變片粘貼的位置,將應變片貼至軸孔受力最大處,將補償片貼至軸孔受力最小處,并使用工裝對應變片進行標定,最后得出拉力與應變的關系進行分析。標定完的應變片可以在試驗過程中通過數采設備獲得應變數據,從而可以算出行星架每個軸孔在特定時刻(使用靜態數采設備)或者特定時段(使用動態數采設備)的均載情況,對星型輪系的合理設計與成功使用具有重要意義。
使用有限元分析軟件對行星架實際受力情況進行模擬,5組軸孔分別受到來自行星輪軸的周向力,所有孔的內外表面受力情況如圖2所示。行星架5個軸孔上下表面處,沿受力方向靠近軸孔邊緣處應力都較大,逆著受力方向靠近孔邊緣處和行星架邊緣處應力都較小。
行星輪軸與行星架為小間隙配合,在轉動方向不變的情況下,齒輪嚙合對行星架產生的周向力會使行星輪軸只有一側與行星架接觸,另一邊不接觸。結合有限元靜力計算分析,行星架上單個軸孔受到應力最大處在沿著受力方向經過軸孔圓心的連線上。在沿著行星架軸孔的受力方向不同距離的上下端面處粘貼多組應變片,大致位置如圖3所示,使用工裝對行星輪軸孔施加拉力。試驗測得,2號、9號和11號應變片位置處應力變化的線性度較好。所以最終將沿著受力方向且過圓心的線上選取最靠近軸孔的一點作為貼片點。

圖2 行星架受力有限元模擬結果

圖3 應變片位置選取
在行星齒輪箱實際運轉過程中,系統由于摩擦會產生大量的熱量,潤滑油和油霧會使系統溫度上升。材料會隨溫度上升而膨脹,導致應變片測量結果偏大。因此,在行星架應力最小且距離應變片距離較近處需要粘貼補償片,補償片會補償應變片由于溫度而產生的誤差。由于行星架受到的轉矩是沿著切線方向的,所以補償片的位置選定在垂直于受力方向且盡量靠近應變片的位置。最終,補償片粘貼在軸孔與內孔圓心的連線上且在軸孔和外圓之間。
行星架上每個孔的內表面在選取的粘貼點水平粘貼了3個應變片作為工作片,再貼了3個應變片作為補償片與3個工作片一一對應。行星架孔1應變片和補償片的粘貼位置如圖4所示,應變片和補償片編號從外向內依次為1號、2號、3號。行星架孔1應變片和補償片的粘貼位置如圖5所示,應變片和補償片編號從外向內依次為1號、2號、3號。2號與4號工作片與補償片距離行星架軸孔中心的距離都約為60 mm。
此試驗需要設計一套工裝,要求工裝有一根工裝軸可以施加沿軸孔實際受力方向的拉力并將所有受到的力傳遞到行星架軸孔上,工裝還需要有4根工裝軸將行星架定位并將受到的力全部傳遞到鑄鐵平臺上。

圖4 行星架孔1應變片實際貼片位置與編號方法

圖5 行星架孔2應變片實際貼片位置與編號方法
行星架應變片標定工裝如圖6所示,底座由螺栓安裝在鑄鐵平臺上,墊板由螺栓安裝在底座上,支撐座和蓋板由螺栓安裝在墊板上。行星架安裝在支撐座內并由定位軸安裝在墊板上,施力軸安裝在行星架上且不與墊板接觸,牽引桿與拉桿通過鍵槽定位安裝在施力軸上,施力軸截面圖如圖7所示。
試驗時在拉桿上施加拉力,拉力通過施力軸全部作用在行星架上,拉力方向與星型齒輪系統受力方向相同。此試驗采用手拉葫蘆施加拉力,使用標定過的電子吊秤來讀取拉力。現場安裝照片如圖8所示,由于安裝空間不足和精度問題,定位軸只安裝了2根。
標定開始前,首先施加適當的拉力,去除行星架與工裝的殘留應力和安裝間隙。試驗開始時,先將鐵鏈保持在松弛狀態,待電子吊秤讀數為0時,采集第一組數據確保應變片為不受力狀態。隨后拉動手拉葫蘆,依次增大拉力至1000、2000、3000、4000、5000、6000、7000 kg,待每個測量點讀數穩定后停留至少2 min再采集應變片數據。然后從7000 kg減少拉力至1000 kg,并每隔1000 kg待拉力穩定后采集應變片數據。為了使數據有對比性,在拉力減少至1000 kg后,繼續加載到7000 kg最后卸載到1000 kg,并每隔1000 kg待拉力穩定后采集應變片數據。做完行星架上一組孔后,旋轉行星架,使用相同的方法采集其余4組軸孔上應變片拉力與應變片讀數。

圖6 行星架應變片標定工裝示意圖

圖7 行星架應變片標定工裝施力軸截面圖

圖8 行星架應變片標定工裝整體照片
由于行星架每組軸孔都有6組應變片,且需要完成2次加載和2次卸載,總計840個數據點。本文挑取了一組軸孔應變片數據進行分析,2次加載與卸載的數據如表1~表4所示。

表1 行星架標定試驗第一次加載應變片數據

表2 行星架標定試驗第一次卸載應變片數據

表3 行星架標定試驗第二次加載應變片數據

表4 行星架標定試驗第二次卸載應變片數據
每次加載與卸載的折線圖如圖9所示,每組應變片在單次加載與卸載過程與拉力中呈現出良好的線性關系。
然后對比各組應變片在每次加載與卸載時的線性趨勢,每次加載與卸載的擬合線性關系如圖10所示。
行星架受力后產生彈性形變,行星架軸孔受到的應力與應變呈線性關系;應變片是由敏感柵構成,其工作原理是基于應變效應制作的,即導體或半導體材料在外界力的作用下產生機械形變時,其電阻值也會相應地發生變化。在理想狀況下,行星架軸孔上施加的拉力應與應變片的讀數呈線性關系。在統計學中,對變量進行線性回歸分析,采用最小二乘法進行參數估計,R平方為回歸平方和與總離差平方和的比值,表示總離差平方和中可以由回歸平方和解釋的比。R平方介于0~1之間,越接近1,模型越精確且回歸擬合效果越好。在分析數據時,可以通過比較R平方的大小來確定拉力與應變片讀數的線性度情況。

圖9 行星架應變片加載與卸載數據線形圖
首先分析6個位置的應變片的線性度,從圖10可知,2號位置應變片的線性較差,R平方值約為0.9703~0.9842之間,1號與3號應變片的線性較好,R平方值可達到0.995左右,4號、5號、6號應變片的線性最好,平均能達到0.997以上。通過圖9可以比較直觀地看到,4號、5號、6號應變片在每次加載與卸載的過程中能保持較好的重疊性,3組應變片拉力與應變讀數的線性關系一致性較好;1號與3號應變片雖然本身線性度較好,但是1號、2號、3號應變片的重疊性較差,3組應變片拉力與應變讀數的線性關系一致性較差。

圖10 各組應變片加載與卸載線性圖
對比其余4組行星輪軸孔,每組軸孔的數據經過分析后得出一樣的結論,4號、5號、6號應變片顯示出非常好的線性度和重合度,1號、2號、3號應變片的線性度和重合度相比較4號、5號、6號應變片來說呈現出較差的趨勢。
從應變片粘貼位置來看,1號、2號、3號應變片粘貼在孔1處,而4號、5號、6號應變片粘貼在孔2處,雖然兩組應變片距離軸孔中心距離都為大約60 mm,但是由于孔1孔徑大于孔2,導致了1號、2號、3號應變片距離受力點相對比較近,而4號、5號、6號應變片則較遠。
根據圣維南原理,分布于彈性體上一小塊面積(或體積)內的荷載所引起的物體中的應力,在離荷載作用區稍遠的地方,基本上只同荷載的合力和合力矩有關;荷載的具體分布只影響荷載作用區附近的應力分布。1號、2號、3號應變片可能距離受力點比較近(荷載作用區),在加載和卸載的過程中,施力軸與軸孔之間的接觸面積在發生變化,導致1號、2號、3號應變片區域的應力分布發生變化,3組應變片的線性與重合性較差。相比較4號、5號、6號應變片,其距離與受力點較遠,應力變化只同荷載的合力與合力矩有關,所以線性和重合性較好。
本文通過有限元分析的方法模擬行星架在系統運轉過程中受到的轉矩,并發現行星架軸孔沿受力方向應力比較大,隨后在一個軸孔附近粘貼多組應變片,使用工裝對該軸孔施加拉力,通過拉力與應變片讀數的線性關系最終確定工作片和補償片的粘貼位置。
通過對行星架5組軸孔的試驗數據進行分析,6組應變片在加載與卸載過程中,拉力與應變片讀數呈現出較好的線性關系。在對比每組應變片的線性關系時發現,4號、5號、6號應變片拉力與應變片讀數的線性關系比較一致,而1號、2號、3號應變片線性關系一致性較差。
在分析問題原因時,發現此現象可能符合圣維南原理,1號、2號、3號應變片處在荷載作用區,受到載荷分布影響大,而4號、5號、6號應變片距離受力點較遠,基本上只同荷載的合力和合力矩有關。
最終,在星型齒輪系統正式運轉時,會優先考慮4號、5號、6號應變片數據進行均載特性的分析。