王桂華，井海龍，史妍妍，陳云霞（.中航工業(yè)沈陽發(fā)動機設(shè)計研究所，沈陽005；.北京航空航天大學(xué)可靠性與系統(tǒng)工程學(xué)院，北京009）

?

基于響應(yīng)面法的齒輪容差設(shè)計方法研究

王桂華1，井海龍2，史妍妍1，陳云霞2
（1.中航工業(yè)沈陽發(fā)動機設(shè)計研究所，沈陽110015；2.北京航空航天大學(xué)可靠性與系統(tǒng)工程學(xué)院，北京100191）

摘要：隨著附件傳動系統(tǒng)不斷向高速、重載方向發(fā)展，齒面接觸可靠性問題越來越受關(guān)注。為提高齒輪可靠性，分析了齒輪傳動系統(tǒng)常見誤差類型及影響，開展了齒輪容差設(shè)計初步探討。根據(jù)齒輪制造、加工及系統(tǒng)動態(tài)誤差的產(chǎn)生原理和作用方式，基于齒輪接觸應(yīng)力仿真結(jié)果，確定了影響齒輪接觸應(yīng)力隨機性的5種主要誤差形式及其特征量分布規(guī)律，提出了基于響應(yīng)面法的齒輪容差設(shè)計方法，給出了齒輪系統(tǒng)公差設(shè)計的改進方向：在進行齒輪公差設(shè)計時，重點控制H、V方向的軸向平行度公差。

關(guān)鍵詞：齒輪；容差設(shè)計；可靠性；接觸應(yīng)力；隨機性；響應(yīng)面；傳動系統(tǒng)；敏感性分析

引用格式：王桂華，井海龍，史妍妍，等.基于響應(yīng)面法的齒輪容差設(shè)計方法研究[J].航空發(fā)動機，2016,42（3）1-6.WANG Guihua，JING Hailong，SHI Yanyan,et al.Methods of tolerance design for gears based on response surface methodology[J].Aeroengine，2016,42（3）：1-6.

0　引言

隨著航空發(fā)動機推重比不斷提高，附件傳動系統(tǒng)逐漸向高速、重載方向發(fā)展，越發(fā)嚴苛的工作條件也對傳動齒輪提出了更高的可靠性要求[1]。在航空發(fā)動機維修過程中發(fā)現(xiàn)，80%以上的齒輪失效是齒面接觸疲勞造成的[2]。因此，研究齒輪的齒面接觸應(yīng)力可靠性，對提高附件傳動系統(tǒng)及發(fā)動機整體工作的可靠性意義重大。

齒輪的工作可靠性受到多種隨機因素的綜合作用，其中制造、加工及系統(tǒng)動態(tài)誤差的存在是隨機性產(chǎn)生的主要原因。文獻[3]針對不同軸線平行度偏差對某附件機匣傳動系統(tǒng)軸系結(jié)構(gòu)的影響進行對比仿真分析；文獻[4]建立了帶有安裝與制造誤差的齒輪參數(shù)化模型，并對其進行動態(tài)仿真，得到齒輪在安裝與制造誤差影響下的動態(tài)接觸應(yīng)力；文獻[5]針對斜齒輪進行建模及仿真，分析了齒輪在齒側(cè)間隙誤差和軸平行度誤差下力學(xué)性能的變化；文獻[6]得到了齒形誤差及公差對齒輪振動的影響。

以上研究僅針對單一誤差對齒輪傳動系統(tǒng)齒輪性能的影響做了分析，而對于多種誤差協(xié)同作用對齒輪性能的影響，及每種誤差因素的影響程度，并沒有給出明確的結(jié)論。

本文根據(jù)齒輪制造、加工及系統(tǒng)動態(tài)誤差的產(chǎn)生原理和作用方式，確定了影響齒輪接觸應(yīng)力隨機性的5種主要誤差形式及其特征量分布規(guī)律，提出了基于響應(yīng)面方法的齒輪公差設(shè)計方法，并給出了齒輪系統(tǒng)公差設(shè)計的改進方向。

1　齒面接觸應(yīng)力隨機性模型的構(gòu)建

1.1齒輪加工誤差的影響分析

受到機床加工精度等客觀因素的限制，在齒輪加工過程中不可避免會存在一定的誤差。根據(jù)對齒輪傳動性能的影響，歸納為主要影響傳遞運動準確性、傳動平穩(wěn)性和載荷分布均勻性的3組誤差項。其中，影響齒面接觸應(yīng)力可靠性的主要是第2、3組，見表1。

1.2齒輪裝配誤差的影響分析

發(fā)動機附件機匣通常采用多平行軸齒輪傳動系統(tǒng)，裝配中主要存在以下2類誤差：

（1）中心距偏差（Δfa）。是指在齒輪副的齒寬中間平面內(nèi)，實際中心距與公稱中心距之差。

（2）軸線平行度誤差。是指1對齒輪的軸線分別在2個基準平面上投影的平行度誤差（在等于齒寬的長度上測量），具體分為軸線平面內(nèi)的偏差（fΣδ）和垂直平面上的偏差（fΣβ），簡稱H、V方向的軸線平行度誤差。具體形式如圖1所示。

表1 影響齒面接觸應(yīng)力的加工誤差

圖1 齒輪系統(tǒng)裝配誤差[7]

1.3齒輪系統(tǒng)動態(tài)誤差的影響分析

齒輪系統(tǒng)動態(tài)誤差是指齒輪系統(tǒng)受到溫度、振動和工作載荷等綜合作用時，由其裝配環(huán)境的變化引起的誤差類型。該類誤差一般與齒輪裝配誤差相互耦合，可以累加入裝配誤差之中。

本文主要針對影響齒輪接觸應(yīng)力的誤差類型進行研究，根據(jù)每類誤差的形成原理和作用方法，結(jié)合工程經(jīng)驗，篩選出影響齒輪接觸應(yīng)力的2類加工誤差和3類裝配誤差，如圖2所示。

根據(jù)齒輪產(chǎn)品在機械加工裝配過程中各設(shè)計參數(shù)分布的特性，結(jié)合齒輪傳動系統(tǒng)的實際設(shè)計資料，以附件機匣1對齒輪為例，得到了其5類誤差參數(shù)，見表2。

圖2 影響齒輪接觸應(yīng)力的誤差類型

表2 齒輪公差信息

1.4構(gòu)建誤差齒輪模型

本文在標準齒輪副應(yīng)力仿真的基礎(chǔ)上，考慮齒輪存在的誤差類型，根據(jù)各類誤差的定義或常見形式，將其直接或間接轉(zhuǎn)化后加入到常規(guī)標準齒輪副中形成齒輪誤差模型，進行基于有限元的誤差齒輪應(yīng)力仿真。

根據(jù)定義，構(gòu)造齒輪副中心距偏差和軸線平行度誤差，如圖3、4所示；將分度圓齒厚誤差轉(zhuǎn)換成齒輪某1齒面的偏置，如圖5所示；在齒輪實際生產(chǎn)制造中，多數(shù)齒形誤差為中凹型，且齒形誤差一般處于齒面中部，所占齒面面積約為1/2，因此，針對齒輪模型構(gòu)造中凹齒形誤差輪齒，如圖6所示。

圖3 中心距偏差建模

圖4 軸線平行度誤差

1.5誤差齒輪的接觸應(yīng)力仿真分析

采用有限元方法進行上述5種誤差因素協(xié)同作用下的齒輪接觸應(yīng)力仿真分析。根據(jù)表2確定的誤差分布規(guī)律，隨機抽樣確定有限元仿真的誤差參數(shù)組合，同時利用有限元瞬態(tài)動力學(xué)仿真確定最大應(yīng)力出現(xiàn)的相位，從而建立處于最大應(yīng)力相位的齒輪誤差模型，進而采用有限元靜力學(xué)仿真進行基于齒輪誤差模型的應(yīng)力仿真分析。具體流程如圖7所示。

圖5 分度圓齒厚誤差

圖6 齒形誤差

圖7 基于誤差的齒輪應(yīng)力仿真流程

2　基于響應(yīng)面的齒輪系統(tǒng)公差設(shè)計

2.1基本原理

基于響應(yīng)面的齒輪公差設(shè)計，首先根據(jù)齒輪誤差參數(shù)的分布類型進行抽樣，確定誤差參數(shù)仿真集，利用誤差齒輪的接觸應(yīng)力仿真分析方法確定誤差參數(shù)及對應(yīng)的齒輪接觸應(yīng)力大小，根據(jù)抽樣樣本及對應(yīng)的應(yīng)力仿真結(jié)果建立齒輪系統(tǒng)應(yīng)力與誤差參數(shù)的響應(yīng)面模型；其次根據(jù)響應(yīng)面模型進行誤差參數(shù)敏感性分析，確定其每種誤差參數(shù)的敏感性及單因素條件下的最壞情況，再進行綜合誤差因素最壞情況分析，確定綜合誤差因素下的最大應(yīng)力，若2種最壞情況的最大應(yīng)力任何1個不滿足要求，則計算其可靠度，若可靠度不滿足要求，則需找出最壞參數(shù)組合，進而根據(jù)單因素敏感性分析結(jié)果進行參數(shù)的優(yōu)化設(shè)計，具體流程如圖8所示。

2.2齒輪系統(tǒng)接觸應(yīng)力響應(yīng)面模型構(gòu)建

響應(yīng)面法[8]的基本原理是：首先假設(shè)1個包含未知系數(shù)的、由狀態(tài)變量與基本變量構(gòu)成的解析表達式，然后用擬合的方法確定未知系數(shù)以表達隱式函數(shù)或高度非線性函數(shù)。多項式系數(shù)的確定一般以試驗設(shè)計為基礎(chǔ)，應(yīng)用正交設(shè)計或均勻設(shè)計回歸得到特定因子的最小二乘估計。采用該方法時，隨機變量的個數(shù)越多，試驗次數(shù)便越多。

圖8 基于響應(yīng)面的齒輪系統(tǒng)公差設(shè)計

在結(jié)構(gòu)力學(xué)分析領(lǐng)域，響應(yīng)面函數(shù)模型常采用2階多項式形式[9]

式中：a、bi、ci、di為待定系數(shù)；xi和xj（i=1,2,…,n;j=1,2,…，n-1）為基本變量。

由于齒輪誤差種類較多且具有截尾分布的特點，在保證抽樣結(jié)果合理性的同時，考慮仿真分析的效率。本文采用超拉丁立方[10]抽樣得到80次仿真的樣本值，建立齒輪系統(tǒng)試驗誤差組合集，并進行誤差齒輪的接觸應(yīng)力仿真分析。

利用前70組數(shù)據(jù)進行響應(yīng)面建模，后10組數(shù)據(jù)進行響應(yīng)面模型的驗證。利用完全2階多項式對前70組數(shù)據(jù)擬合得到接觸應(yīng)力的響應(yīng)面模型

表3 響應(yīng)面模型驗證結(jié)果

式中：σ為接觸應(yīng)力；X11為分度圓齒厚誤差，單位為mm；X1為有限元仿真輸入?yún)?shù)，反映了分度圓齒厚誤差，單位為mm；X2為齒形誤差，單位為mm；X3為中心距偏差，單位為mm；X4為H方向軸向平行度誤差，單位為°；X5為V方向軸向平行度誤差，單位為°。

利用后10組數(shù)據(jù)進行接觸應(yīng)力模型的驗證，驗證結(jié)果見表3。

從表中可見，齒輪系統(tǒng)接觸應(yīng)力響應(yīng)面模型預(yù)測結(jié)果與仿真分析結(jié)果誤差在10%以內(nèi)，認為模型可以作為后續(xù)分析的依據(jù)。

利用協(xié)方差矩陣進行貢獻率分析，得到分度圓齒厚誤差、齒形誤差、中心距偏差和H、V方向軸向平行度誤差這5種誤差參數(shù)對應(yīng)力分析結(jié)果的貢獻率[11]分別為26.5874%、22.8871%、20.8440%、18.2140%和11.4675%。

圖9 齒厚誤差敏感性分析結(jié)果

2.3誤差參數(shù)敏感性分析

敏感性分析是指從定量分析的角度研究有關(guān)因素發(fā)生某種變化對某1個或1組關(guān)鍵指標影響程度的1種不確定分析技術(shù)。其實質(zhì)是通過逐一改變相關(guān)變量數(shù)值來解釋關(guān)鍵指標受這些因素變動影響大小的規(guī)律。

圖10 齒形誤差敏感性分析結(jié)果

分別將不同的齒輪誤差因素設(shè)為單一變量，通過固定其它4種誤差因素，對每種誤差因素進行齒輪接觸應(yīng)力敏感性分析。

首先以齒厚誤差為例，分別固定齒形誤差、中心距偏差和H、V方向軸向平行度誤差為均值，分析齒厚誤差對齒輪系統(tǒng)接觸應(yīng)力的影響，如圖9所示。

圖11 中心距偏差敏感性分析結(jié)果

圖12 軸向平行度偏差敏感性分析結(jié)果

從圖中可見，在固定齒形誤差、中心距偏差和H、V方向軸向平行度誤差這4個誤差因素的條件下，分度圓齒厚誤差在[0.07，0.14]范圍內(nèi)變化時，齒輪的接觸應(yīng)力范圍為[744.7，762.8]，應(yīng)力幅值相對于應(yīng)力均值的變化率為0.045%，變化不明顯。

以此類推，得到：

（1）齒形誤差在[0.003，0.006]范圍內(nèi)變化時，齒輪的接觸應(yīng)力范圍為[704.7，761.9]，應(yīng)力幅值相對于應(yīng)力均值的變化率為4.04%，變化不明顯，如圖10所示。

（2）中心距偏差在[-0.015，0.015]范圍內(nèi)變化時，齒輪的接觸應(yīng)力范圍為[761.1，814.5]，應(yīng)力幅值相對于應(yīng)力均值的變化率為3.39%，變化不明顯，如圖11所示。

（3）H方向軸向平行度誤差在[-0.0572，0.0572]范圍內(nèi)變化時，齒輪的接觸應(yīng)力范圍為[761.9，1045.4]，應(yīng)力幅值相對于應(yīng)力均值的變化率為15.7%，變化明顯；V方向軸向平行度誤差在[-0.0572，0.0572]范圍內(nèi)變化時，齒輪的接觸應(yīng)力范圍為[758.5，834.2]，最大應(yīng)力相對最小應(yīng)力變化4.75%，變化不明顯，如圖12所示。

根據(jù)單誤差因素敏感性分析結(jié)果，確定各誤差參數(shù)的敏感性由高到低為：H方向軸向平行度誤差、V方向的軸向平行度誤差與齒形誤差、中心距偏差和分度圓齒厚誤差。因此，在進行公差設(shè)計時，應(yīng)按上述敏感性從高到低考慮。

2.4綜合誤差因素最壞情況分析

最壞情況分析法是1種按照不常發(fā)生的最壞使用條件的組合為基礎(chǔ)，進行確保系統(tǒng)可靠性要求的設(shè)計方法[12]，是1種非概率統(tǒng)計方法，用來分析系統(tǒng)的設(shè)計參數(shù)和各種內(nèi)、外影響因素的參數(shù)處于最壞情況或最壞組合情況下的系統(tǒng)性能參數(shù)偏差。對于大規(guī)模復(fù)雜系統(tǒng)來說，用其分析各類參數(shù)處于最壞情況下系統(tǒng)的性能是否還滿足指標要求，并找出導(dǎo)致系統(tǒng)性能參數(shù)處于最壞情況下的設(shè)計參數(shù)和影響因素的組合規(guī)律。

應(yīng)用該方法的前提是假設(shè)各類影響因素對系統(tǒng)的影響是呈單調(diào)性的，即影響因素的參數(shù)取標稱值時系統(tǒng)處于最理想狀態(tài)，其參數(shù)偏離標稱值越大對系統(tǒng)的影響也越大[13]。由敏感性分析可知，本文研究的5個誤差因素只有H、V方向的軸向平行度誤差對系統(tǒng)的影響呈現(xiàn)單調(diào)性。結(jié)合接觸應(yīng)力的響應(yīng)面模型，采取仿真與數(shù)值分析相結(jié)合的方法進行最壞情況分析，H、V方向的軸向平行度誤差2個參數(shù)采取上、下極限值全面組合的方法，其余3個參數(shù)采取數(shù)值分析的方法進行3參數(shù)下的最大應(yīng)力分析。具體流程如圖13所示。

圖13 綜合因素最壞情況分析流程

分別對H、V方向軸向平行度誤差的極限值（± 0.0572°）進行全面組合，確定分度圓齒厚誤差、齒形誤差和中心距偏差3參數(shù)的響應(yīng)面，利用數(shù)值的方法確定其最大應(yīng)力，分析結(jié)果見表4。

表4 綜合因素最壞情況分析結(jié)果

齒輪的接觸疲勞極限為1500 MPa，其設(shè)計安全系數(shù)為1.0。從表中可見，在極限1、4狀態(tài)下，不滿足齒輪接觸應(yīng)力的要求。

2.5基于響應(yīng)面的齒輪系統(tǒng)可靠度分析

將齒輪的齒面安全系數(shù)作為可靠性分析的指標，由齒輪設(shè)計要求可知齒輪齒面的計算安全系數(shù)不能小于設(shè)計安全系數(shù)。因此，齒輪副的可靠度定義為在隨機誤差綜合影響下齒輪副的計算安全系數(shù)不大于設(shè)計安全系數(shù)的概率，表達式如下

式中：R為齒輪副可靠度；n為齒輪齒面的計算安全系數(shù)；ns為齒面設(shè)計安全系數(shù)。齒根的安全系數(shù)定義為[14]

式中：σHlim為材料的接觸疲勞極限；σH為齒輪齒面應(yīng)力。

本文基于蒙特卡洛思想[15]，根據(jù)齒輪誤差因素的分布特點，對其公差構(gòu)成的隨機變量進行大量重復(fù)抽樣，將樣本值帶入響應(yīng)面模型，確定每組誤差組合的接觸應(yīng)力響應(yīng)，得到每組誤差組合下的安全系數(shù)，利用式（4）計算其可靠度為0.9923。

3　結(jié)論

本文根據(jù)齒輪制造、加工及系統(tǒng)動態(tài)誤差的產(chǎn)生原理和作用方式，確定了影響齒輪接觸應(yīng)力隨機性的5種主要誤差形式及其特征量分布規(guī)律，提出了基于響應(yīng)面方法的齒輪公差設(shè)計方法，并給出了齒輪系統(tǒng)公差設(shè)計的改進方向。

分析結(jié)果表明：影響齒輪接觸應(yīng)力隨機性的誤差主要有5種，即分度圓齒厚誤差、齒形誤差、中心距偏差和H、V方向軸向平行度誤差。在進行齒輪公差設(shè)計時，應(yīng)重點控制H、V方向的軸向平行度誤差。

參考文獻：

[1]林基恕.航空燃氣渦輪發(fā)動機機械系統(tǒng)設(shè)計[M].北京：航空工業(yè)出版社，2005：157-160. LIN Jishu.Mechanical system design for aero gas turbine engines[M]. Beiijng：Aviation Industry Press，2005：157-160（.in Chinese）

[2]史妍妍.基于熱分析的附件機匣若干可靠性問題研究[D].沈陽：東北大學(xué)，2009. SHI Yanyan.Study on accessory gearbox reliability based on thermal analysis[D].Shenyang：Northeastern University，2009（.in Chinese）

[3]王宇寧，孫志禮，楊麗，等.考慮裝配誤差的軸系結(jié)構(gòu)靜力學(xué)對比分析[J].中國工程機械學(xué)報，2014，12（1）：13- 17.WANG Yuning，SUN Zhili，YANG Li，et al.Statics comparative analy-sis on shafting structures with consideration of assembly error [J].Chinese Journal of Construction，2014，12（1）：13-17.（in Chinese）

[4]佟操，孫志禮，馬小英，等.考慮安裝與制造誤差的齒輪動態(tài)接觸仿真[J].東北大學(xué)學(xué)報，2014，35（7）：996-1000. TONG Cao，SUN Zhili，MA Xiaoying，et al.Dynamic simulation of spur gears with assembly errors and machining errors [J].Journal of North-eastern University（Natural Science），2014，35（7）：996-1000.（in Ch-inese）

[5]孫之星.考慮裝配誤差的齒輪接觸應(yīng)力仿真研究[D].太原：太原理工大學(xué)，2010. SUN Zhixing.Contact stress simulation study of gears with assembly error[D].Taiyuan：Taiyuan University of Technology，2010.（in Chinese）

[6]Bonori G，Pellicano F.Non-smooth dynamics of spur gears with manu-facturing errors [J]. Journal of Sound and Vibration，2007，30（6）：271-283.

[7]方昆凡.公差與配合使用手冊[M].北京：機械工業(yè)出版社，2012：35-40. FANG Kunfan.Manual of tolerances and fits [M].Beijing：Machinery Industry Press，2012：35-40.（in Chinese）

[8]Chamoli S.ANN and RSM approach for modeling and optimization of designing parameters for a V down perforated baffle roughened rectan-gular channel [J]. Alexandria Engineering Journal，2015，54（3）：429-446.

[9] Myers R H，Montgomery D C，Vining G G，et al. Response surface methodology:aretrospective and literature survey[J].Journal of Quality Technology，2004，36（1）：53-77.

[10]Gangel M，Seiler M J，Collins A. Latin hypercube sampling and the identification of the foreclosure contagion threshold [J].Journal of Behavioral Finance，2013，14（2）：149-159.

[11]PianosiF，Sarrazin F，W agener T.A matlab toolbox for global sensi-tivity analysis [J]. Environmental Modeling & Software，2015，70：80-85.

[12]陳云霞，段朝陽.飛控系統(tǒng)最壞情況分析方法研究[J].航空學(xué)報，2005，26（5）：647-651. CHEN Yunxia，DUAN Chaoyang.Study on worst case analysis method for flight control system [J].Acta Aeronautica et Astronautica Sinica，2005，26（5）：647-651.（in Chinese）

[13]陳云霞.性能與可靠性一體化建模和分析方法研究[D].北京：北京航空航天大學(xué)，2004. CHEN Yunxia.Study on the method of integrated modeling and anal-ysis for performance and reliability [D].Beijing：Beihang University，2004.（in Chinese）

[14]孫自龍，閆鳳英，劉麗.齒輪加工中剃齒的誤差分析及處理[J].金屬加工（冷加工），2014（5）：63-65. SUN Zilong，YAN Fengying，LIU Li.Error analysis and treatment for gear shaving during gearmachining [J].Metal W orking（Cold），2014 （5）：63-65.（in Chinese）

[15]W alter J C，Barkema G T.An introduction to Monte Carlo methods[J]. Physica A: Statistical Mechanics and Its Applications，2015，418：78-87.

（編輯：栗樞）

Methods of Tolerance Design for Gears Based on Response Surface Methodology

WANG Gui-hua1，JING Hai-long2，SHIYan-yan1，CHEN Yun-xia2
（1.AVIC Shenyang EngineDesign and Research Institute，Shenyang 10015；2.School of Reliability and Systems Engineering，Beihang University，Beijing 100191）

Abstract:With the accessory driving system takes the way towards the high speed and heavy load,the contact reliability of aviation gears make more sense. In order to improve the reliability,the common error types and effects of gear driving system were analyzed,and the study of tolerance design was proposed.According to gear rtypes and their characteristics of the random error of the contact stress of the gear were determined on the basis of the simulation results.Moreover,a method of gear tolerance design based on response surface was proposed and the improvement direction of tolerance design of gear system was given:the tolerance for the direction of H and V should be mainly controlled during the tolerance design.

Key words:gear；tolerance design；reliability；contact stress；randomness；response surface；driving system;sensitivity analysis

中圖分類號：V215.1

文獻標識碼：A

doi：10.13477/j.cnki.aeroengine.2016.03.001

收稿日期：2015-12-30基金項目：航空動力基礎(chǔ)研究項目資助

作者簡介：王桂華（1965），女，博士，自然科學(xué)研究員，從事航空發(fā)動機強度與可靠性設(shè)計工作；E-mail：wangguih@sina.com。