航空發(fā)動機套齒結(jié)構(gòu)配合穩(wěn)健性優(yōu)化設(shè)計

劉宏蕾，陳志英，周平

（北京航空航天大學(xué)能源與動力工程學(xué)院，北京100191）

航空發(fā)動機套齒結(jié)構(gòu)配合穩(wěn)健性優(yōu)化設(shè)計

劉宏蕾，陳志英，周平

（北京航空航天大學(xué)能源與動力工程學(xué)院，北京100191）

為了更合理地設(shè)計套齒結(jié)構(gòu)的配合關(guān)系，提高套齒結(jié)構(gòu)配合精度，適應(yīng)不對中運行工況，降低對加工誤差的敏感度，研究了航空發(fā)動機套齒結(jié)構(gòu)在傾角不對中的情況下的運動特點以及配合關(guān)系，并推導(dǎo)了含有傾角不對中的套齒結(jié)構(gòu)齒側(cè)間隙以及與之相關(guān)的套齒特征量的計算公式，對某型航空發(fā)動機的壓氣機軸徑套齒結(jié)構(gòu)進行了計算與分析。在此基礎(chǔ)上，考慮到參數(shù)的隨機性，應(yīng)用果蠅優(yōu)化算法對套齒結(jié)構(gòu)齒側(cè)間隙進行了穩(wěn)健性優(yōu)化設(shè)計。結(jié)果表明：計算公式中考慮傾角的影響能夠有效地避免運動干涉和裝配預(yù)應(yīng)力，穩(wěn)健性優(yōu)化設(shè)計可以有效地避免由于傾角不對中、加工誤差等參數(shù)隨機因素導(dǎo)致齒側(cè)間隙的設(shè)計值偏大。證明了所提出的計算和優(yōu)化方法可以使套齒結(jié)構(gòu)配合關(guān)系的設(shè)計更加合理，可為套齒結(jié)構(gòu)設(shè)計、裝配提供參考。

套齒結(jié)構(gòu)；傾角不對中；齒側(cè)間隙；果蠅優(yōu)化算法；穩(wěn)健性優(yōu)化；航空發(fā)動機

0 引言

套齒連接結(jié)構(gòu)在航空發(fā)動機中的應(yīng)用非常廣泛[1-2]。由于套齒結(jié)構(gòu)加工安裝誤差、軸承不同心以及機座高低偏差等因素會使套齒存在不對中[3]的情況。在旋轉(zhuǎn)機械轉(zhuǎn)子故障中，不對中故障占60%以上[4]。故障使得齒間發(fā)生運動干涉以及配合狀態(tài)的不穩(wěn)定，從而使系統(tǒng)可靠性降低。不共軸時的安全系數(shù)一般降低到共軸時的1/5～1/10[5]。轉(zhuǎn)子系統(tǒng)中套齒結(jié)構(gòu)不對中普遍存在，所以在設(shè)計階段就應(yīng)對配合關(guān)系進行合理分析和設(shè)計，對于提高轉(zhuǎn)子系統(tǒng)的穩(wěn)定性具有重要意義。

套齒等機械結(jié)構(gòu)裝配尺寸鏈及配合公差的分析與設(shè)計，屬于機械確定性裝配關(guān)系設(shè)計范疇，在參數(shù)確定狀態(tài)下裝配尺寸鏈及各零件的公差配合關(guān)系研究，目前已經(jīng)比較成熟，但近幾年國內(nèi)外正在發(fā)展研究一些新的裝配尺寸鏈及配合公差的不確定分析與設(shè)計[6-9]。由于加工誤差、工作環(huán)境的不確定性，使得必須要考慮參數(shù)的不確定性對套齒結(jié)構(gòu)配合的影響才能更加真實地反映配合狀態(tài)。

穩(wěn)健性設(shè)計是由日本質(zhì)量專家田口博士[10]提出的1種質(zhì)量改進技術(shù)，能夠有效減小產(chǎn)品質(zhì)量波動。近年來，工程領(lǐng)域引入穩(wěn)健性設(shè)計思想收到了很好的效果[11-13]。該設(shè)計的2項基本原則是使目標(biāo)值盡可能接近最優(yōu)均值且波動范圍盡可能小。

本文為了更好地確定套齒配合的設(shè)計參數(shù)，研究了傾角不對中的套齒結(jié)構(gòu)的運動特點、配合關(guān)系，推導(dǎo)了齒側(cè)間隙等套齒特征量的計算公式，進行了穩(wěn)健性優(yōu)化設(shè)計，以期降低對加工誤差參數(shù)隨機性的敏感度。

1 套齒結(jié)構(gòu)傾角不對中運動配合關(guān)系

在3點支承的方案中，由于軸和轉(zhuǎn)子支承零件加工、裝配誤差，以及安裝座位置的誤差，造成壓氣機轉(zhuǎn)子和渦輪轉(zhuǎn)子2軸線彼此不對中。由于壓氣機后軸承靠近套齒結(jié)構(gòu)，使得在這種不對中情況下，內(nèi)外套齒軸線有一定傾角，且內(nèi)、外套齒中心在徑向重合，稱為傾角不對中，如圖1所示。壓氣機軸與渦輪軸存在傾角φ（一般φ值很小）。

圖1 壓氣機轉(zhuǎn)子與渦輪轉(zhuǎn)子傾角不對中

航空漸開線套齒齒側(cè)既能起到驅(qū)動作用，又有自動定心的作用，在航空發(fā)動機中應(yīng)用最為廣泛。當(dāng)內(nèi)、外套齒存在不對中時，會影響套齒齒側(cè)的最小作用間隙，應(yīng)當(dāng)調(diào)整配合予以補償。為了合理地設(shè)計傾角不對中的套齒結(jié)構(gòu)，需要對這種狀況下套齒的運動狀態(tài)和配合關(guān)系進行分析研究。

1.1 運動關(guān)系

從運動狀態(tài)上來說，由于渦輪軸與壓氣機軸存在傾角φ，每旋轉(zhuǎn)1周，外套齒都會相對于內(nèi)套齒在齒長方向上作往復(fù)滑動，在圓周切線方向上繞齒中心做往復(fù)擺動。這就要求沿著齒長方向內(nèi)套齒的齒長要略大于外套齒的齒長；在周向，由于存在偏角引起外套齒在內(nèi)套齒往復(fù)擺動，這就需要在原來按照標(biāo)準(zhǔn)選定齒側(cè)間隙的基礎(chǔ)上增大設(shè)計側(cè)隙，以免內(nèi)、外套齒發(fā)生運動干涉。

1.2 配合關(guān)系

當(dāng)內(nèi)、外套齒完全對中時，其軸截面和A、B 2個端面如圖2所示。從圖中可見，內(nèi)、外套齒中心線重合，并且每對齒的齒側(cè)間隙相同。

圖2 內(nèi)、外套齒完全對中

如果內(nèi)、外套齒存在傾角（如圖3所示），并且假設(shè)此時傾角的方向是鉛垂方向，即內(nèi)套齒的軸線相對于外套齒在鉛垂面內(nèi)有個傾角φ，那么對鉛垂方向上、下2部分嚙合齒的齒側(cè)間隙影響較小，但使水平方向左、右2部分的嚙合齒的齒側(cè)間隙減小較大。具有最小齒側(cè)間隙的嚙合齒會承受最大的負(fù)荷。此時處于左右象限點的嚙合齒最先達到危險狀態(tài)。

圖3 內(nèi)、外套齒傾角不對中

由此可知，從配合狀態(tài)上來說，在垂直于傾角方向上齒側(cè)間隙減小最為迅速。在設(shè)計側(cè)隙時應(yīng)選處于該位置的鍵齒進行分析和計算。為減少裝配后由于該位置齒側(cè)間隙過小以至于卡死形成預(yù)應(yīng)力，一般需要加大設(shè)計間隙。

1.3 仿真驗證

為了更直觀地觀察以上規(guī)律，應(yīng)用數(shù)學(xué)建模仿真方法驗證。設(shè)定內(nèi)、外套齒的側(cè)隙為0.17 mm，如圖4所示。固定一端，在另外一端加1個小的垂直向下的壓力，在該壓力作用下，內(nèi)、外套齒軸線會產(chǎn)生鉛垂方向上的傾角。當(dāng)傾斜量超過現(xiàn)有側(cè)隙值時，最先發(fā)生接觸的地方就會產(chǎn)生應(yīng)力，如圖5、6所示。

圖4 側(cè)隙為0.17 mm的套齒

圖5 外套齒等效應(yīng)力分布

圖6 內(nèi)套齒等效應(yīng)力分布

從圖中可見，當(dāng)產(chǎn)生沿著鉛垂方向上的傾角時，在水平方向的左右兩端及其附近的齒最先發(fā)生卡死的現(xiàn)象，產(chǎn)生應(yīng)力，與以上分析相符。

2 套齒結(jié)構(gòu)傾角不對中齒側(cè)間隙設(shè)計

2.1 含有傾角不對中的內(nèi)外套齒特征量的推導(dǎo)

當(dāng)壓氣機軸與渦輪軸有偏角φ時，壓氣機與渦輪軸上的內(nèi)、外套齒之間會產(chǎn)生偏斜。而在航空減速器中，軸線偏斜的情況也很常見。一般認(rèn)為，當(dāng)傾角小于10'時，采用直齒形聯(lián)軸器；大于10'時，采用鼓形齒式聯(lián)軸器。

所以在齒側(cè)間隙設(shè)計時，除了依照截面參數(shù)和齒長查設(shè)計手冊選取正常的側(cè)隙△j之外，還應(yīng)該考慮由于傾角φ導(dǎo)致的外套齒相對于內(nèi)套齒的切向偏移量ft，進而引起側(cè)隙的增大。

以上分析可知，在垂直于傾角方向上的鍵齒齒側(cè)間隙減小最為迅速，亦即最容易達到危險狀態(tài)。假設(shè)傾角是沿垂直方向的，取水平方向上的1個外齒作為研究對象，沿著分度圓柱面展開這個齒，如圖7所示。其中b為齒長，S為齒厚。

圖7 嚙合齒展開

將不必要的輪廓線略去，留下關(guān)鍵的尺寸線以及幾何關(guān)系，如圖8所示。

圖8 關(guān)鍵尺寸幾何關(guān)系

則由于存在傾角而需要加大的間隙值ft為

式中：ft1、ft2為相互配合的齒槽與齒厚兩側(cè)的間隙值。

最終的套齒結(jié)構(gòu)齒側(cè)間隙為

在套齒結(jié)構(gòu)的設(shè)計、加工、裝配過程中會涉及到一些典型的特征量：內(nèi)套齒最小和最大作用齒槽寬EVmin和EVmax、外套齒最大和最小作用齒厚SVmin和SVmax、內(nèi)套齒最小和最大實際齒槽寬Emin和Emax、外套齒最大和最小實際齒厚Smin和Smax。其中，內(nèi)套齒最小作用齒槽寬和外套齒最大作用齒厚決定了套齒結(jié)構(gòu)作用側(cè)隙的最小值；內(nèi)套齒最大作用齒槽寬和外套齒最小作用齒厚決定了套齒結(jié)構(gòu)作用側(cè)隙的最大值；內(nèi)套齒最小實際齒槽寬和外套齒最大實際齒厚一般在單項檢驗法中使用；內(nèi)套齒最大實際齒槽寬和外套齒最小實際齒厚是零件合格與否的依據(jù)。

如果不考慮傾角不對中，套齒各特征量計算公式如下：

內(nèi)套齒最小作用齒槽寬

內(nèi)套齒最大作用齒槽寬

外套齒最大作用齒厚

外套齒最小作用齒厚

內(nèi)套齒最大實際齒槽寬

內(nèi)套齒最小實際齒槽寬

外套齒最大實際齒厚

外套齒最小實際齒厚

式中：m為模數(shù)；esv為齒厚極限上偏差；λ為綜合公差；T為加工公差；T+λ為總公差。

在套齒結(jié)構(gòu)傾角不對中的情況下，需要加大齒側(cè)間隙可加大齒槽寬和減小齒厚。一般情況下，考慮到航空套齒結(jié)構(gòu)中內(nèi)套齒的加工難度要大于外套齒，所以推導(dǎo)傾角不對中公式時采用移動外套齒的公差帶來實現(xiàn)不同齒側(cè)間隙配合。

基于以上考慮，推導(dǎo)出套齒各特征量計算公式：

內(nèi)套齒最小作用齒槽寬

內(nèi)套齒最大作用齒槽寬

外套齒最大作用齒厚

外套齒最小作用齒厚

內(nèi)套齒最大實際齒槽寬

內(nèi)套齒最小實際齒槽寬

外套齒最小實際齒厚

外套齒最大實際齒厚

式中：Cmin為標(biāo)準(zhǔn)側(cè)隙，根據(jù)分度圓直徑以及所選取的基本偏差等級確定。

將含有傾角不對中的各裝配特征量用公差帶的形式表示，如圖9所示。

圖9 傾角不對中裝配特征量公差帶

在航空減速器中，由于尺寸空間的限制，也會出現(xiàn)一些強度裕度不大的情況，考慮到內(nèi)套齒的承載能力比外套齒的強，也可以用增大內(nèi)套齒齒槽寬的方法來確定套齒結(jié)構(gòu)傾角不對中時的特征量。

式中：Fp為齒距累積公差；Ff為齒形公差；其它特征量可按基本公式推出。

上述各式對于齒形可用于鼓形齒及2軸角偏斜大于10'情況，需要注意的是，對于鼓形套齒來說，外套齒相對于內(nèi)套齒的周向擺動引起的切向位移量ftkb與直齒的不同，其計算公式為

式中：Rt為鼓形齒沿齒長方向在分度圓截面上齒牙側(cè)面輪廓曲率半徑。

2.2 實例計算

某型航空發(fā)動機的壓氣機軸徑內(nèi)套齒基本參數(shù)見表1。

表1 內(nèi)套齒參數(shù)

計算所需公差值按照GB/T3478.1—1995中側(cè)隙d級和5級精度選取。具體參數(shù)見表2。

有無傾角的不對中情況的計算公式對比見表3。

假設(shè)壓氣機軸與渦輪軸的偏斜角φ=10'，將表1、2中的數(shù)據(jù)依次帶入表3中的公式，得到套齒結(jié)構(gòu)有無傾角不對中時各特征量，見表4。

表2 按照GB/T3478.1-1995選取的公差值

表3 有無傾角不對中情況的計算公式對比

表4 無傾角與傾角φ=10'不對中時的特征量

利用UG中干涉分析功能對表中的數(shù)據(jù)進行檢驗。取內(nèi)套齒最小作用齒槽寬EVmin與外套齒最大作用齒厚SVmax兩個極限尺寸進行分析。如果2個極限尺寸所對應(yīng)的套齒結(jié)構(gòu)不會發(fā)生運動干涉，那么2個套齒結(jié)構(gòu)取公差帶內(nèi)的尺寸就都不會發(fā)生運動干涉。

設(shè)置內(nèi)、外套齒為連桿，以外套齒的軸線中心點為中心設(shè)置轉(zhuǎn)動副，如圖10所示。設(shè)定10'的角度，檢查內(nèi)、外套齒是否干涉，結(jié)果如圖11所示。

圖10 檢查干涉參數(shù)設(shè)置

圖11 檢查干涉結(jié)果

將整個過程的時間步全部顯示出來，列出干涉的文件顯示沒有干涉，采用以上方法計算出含有傾角不對中的套齒結(jié)構(gòu)能很好地避免由于傾角引起的套齒間運動干涉，有效避免裝配應(yīng)力的產(chǎn)生。

3 考慮隨機因素齒側(cè)間隙穩(wěn)健性優(yōu)化設(shè)計

在成批生產(chǎn)的發(fā)動機中齒側(cè)間隙的設(shè)計按照最大傾角設(shè)計，使得套齒結(jié)構(gòu)產(chǎn)生過大的齒側(cè)間隙。而最大的傾角出現(xiàn)的概率可能很小，按照確定性設(shè)計的套齒間隙會使設(shè)計值偏大，在套齒運行工作時，過大的齒側(cè)間隙將會使套齒發(fā)生沖擊載荷，齒側(cè)間隙越大，沖擊載荷也越大。

也就是說，基于確定性的設(shè)計是按照φ值最大來設(shè)計的，由于加工制造的原因，φ值、作用齒槽寬EV、作用齒厚SV、齒長b都是隨機分布的，如果φ值取最大的情況概率極其微小，該設(shè)計就偏于保守，在套齒這種以齒側(cè)間隙為關(guān)鍵技術(shù)指標(biāo)的配合中是不合理的，需要對其進行優(yōu)化。

果蠅優(yōu)化算法（Fruit Fly Optimization Algorithm，F(xiàn)FOA）由中國臺灣學(xué)者潘文超[14]提出，是1種基于果蠅覓食行為而推演出的尋求全局優(yōu)化的新方法。將該優(yōu)化算法用于工程中，可以避免優(yōu)化結(jié)果陷入局部最優(yōu)，而且收斂速度快[15]，所以采用果蠅優(yōu)化算法對傾角不對中的套齒結(jié)構(gòu)齒側(cè)間隙進行優(yōu)化。FFOA的算法流程如圖12所示。

圖12 FFOA的算法流程

其具體步驟為：

（1）果蠅群體位置隨機初始化。

（2）給利用嗅覺搜尋食物的果蠅個體賦予隨機方向與距離。

（3）由于無法知道食物位置，因此先估計與原點之距離(Dist)，再計算味道濃度判定值S，此值為距離的倒數(shù)。

（4）將味道濃度判定值S代入味道濃度判定函數(shù)，以求出該果蠅個體位置的味道濃度(Smelli)。

（5）找出此果蠅群體的中味道濃度最高的果蠅

（6）保留最佳味道濃度值與x、y坐標(biāo)，此時果蠅群體利用視覺往該位置飛去。

（7）進入迭代尋優(yōu)，重復(fù)執(zhí)行步驟（2）～（5），并判斷味道濃度是否優(yōu)于之前迭代味道濃度，若是則執(zhí)行步驟（6）。

3.1 套齒結(jié)構(gòu)配合齒側(cè)間隙穩(wěn)健性優(yōu)化設(shè)計模型

加工制造誤差的隨機性以及傾斜角的隨機性都會使側(cè)隙設(shè)計偏大，需要優(yōu)化諸如齒槽寬、外套齒的齒厚以及齒長等設(shè)計參數(shù)，使得齒側(cè)間隙的值盡可能小。

間隙的設(shè)計需要考慮按照標(biāo)準(zhǔn)應(yīng)留的齒側(cè)間隙值，還要考慮由于傾角不對中需要預(yù)留的齒側(cè)間隙值，除此之外再有的間隙就是由于加工制造裝配誤差所造成的間隙，稱作額外齒側(cè)間隙。額外齒側(cè)間隙會嚴(yán)重影響套齒結(jié)構(gòu)運行狀態(tài)的穩(wěn)定性，引起嚴(yán)重的沖擊載荷，所以取額外齒側(cè)間隙值為優(yōu)化目標(biāo)函數(shù)，希望其越小越好，其具體表達式如下

優(yōu)化目標(biāo)函數(shù)及約束條件如下

3.2 穩(wěn)健性優(yōu)化變量的選取

隨機變量可以分為2類：1類是可控隨機變量，在設(shè)計與制造過程中可以控制，因此能夠作為待優(yōu)化的參數(shù)，可以通過改變其數(shù)值來進行尋優(yōu)。作用齒槽寬EV，作用齒厚SV以及齒長b的尺寸及其公差帶的大小和位置都可以在設(shè)計和制造中加以控制，因此可以作為可控變量。需要說明的是，齒長b的確定和結(jié)構(gòu)尺寸限制以及承載能力的設(shè)計有關(guān)，但在如公差范圍這樣小的尺寸范圍內(nèi)，可以對其優(yōu)化。另外1類是不可控隨機變量，該類隨機變量不可以進行控制，只能作為隨機輸入?yún)?shù)而不能視為待優(yōu)化參數(shù)。傾角φ受安裝位置、工作環(huán)境的影響，不能對其進行控制，因此作為不可控變量。假設(shè)以上隨機變量均為正態(tài)分布?？紤]到航空套齒結(jié)構(gòu)的精度要求以及加工工藝所能達到的精度水平，各變量具體情況見表5。

表5 穩(wěn)健性優(yōu)化隨機變量

3.3 優(yōu)化過程和結(jié)果分析

利用FFOA對考慮隨機因素的套齒結(jié)構(gòu)配合齒側(cè)間隙進行穩(wěn)健性優(yōu)化設(shè)計，優(yōu)化前后的裝配關(guān)鍵特征量的對比見表6。

表6 優(yōu)化前后數(shù)據(jù)對比mm

從表中可見，相對于優(yōu)化之前，優(yōu)化后的公差帶長度沒有變化，均為0.045 mm，表明加工制造費用和之前一樣。而額外齒側(cè)間隙值由之前的0.0446 mm減小到0.0201 mm，減小了55%，所以進行穩(wěn)健性優(yōu)化的效果非常明顯，可以在現(xiàn)有加工水平的基礎(chǔ)上提高套齒結(jié)構(gòu)齒側(cè)間隙對抗參數(shù)波動干擾能力。

4 結(jié)論

在套齒結(jié)構(gòu)配合關(guān)系設(shè)計之初就應(yīng)該充分考慮轉(zhuǎn)子系統(tǒng)不對中所引起的額外齒側(cè)間隙，計算公式中考慮傾角的影響能夠有效地避免運動干涉和裝配預(yù)應(yīng)力。

采用果蠅優(yōu)化算法對齒側(cè)間隙進行穩(wěn)健性優(yōu)化設(shè)計，可以有效地避免由于傾角不對中以及工作環(huán)境、加工誤差等參數(shù)隨機因素的存在導(dǎo)致齒側(cè)間隙的設(shè)計值偏大，使設(shè)計更加合理。

以上對于轉(zhuǎn)子系統(tǒng)中存在傾角不對中的套齒結(jié)構(gòu)特征量的計算方法，以及考慮到參數(shù)的隨機性特點，對齒側(cè)間隙穩(wěn)健性優(yōu)化設(shè)計的方法可以為套齒結(jié)構(gòu)設(shè)計、裝配提供參考，在實際的設(shè)計、裝配、使用過程中需要重視。

[1] 陳光.航空發(fā)動機結(jié)構(gòu)設(shè)計分析[M].北京：北京航空航天大學(xué)出版社，2014：349-393.CHEN Guang.Structural design analysis of aeroengine[M].Beijing：Beihang University Press,2014：349-393.（in Chinese）

[2] 航空發(fā)動機設(shè)計手冊總編委會.航空發(fā)動機設(shè)計手冊（第13分冊）：減速器[M].北京：航空工業(yè)出版社，2001：251-275.Aeroengine Design Manual General Editorial Board.Aeroengine design manual(album 13)：deceleration device[M].Beijing：Aviation Industry Press，2001：251-275.（in Chiense）

[3] 王聰梅.航空發(fā)動機典型零件機械加工[M].北京：航空工業(yè)出版社，2014：128-185.WANG Congmei.Mechanical processing of typical parts of aeroengine[M].Beijing：Aviation Industry Press,2014：128-185.（in Chinese）

[4] 劉占生.轉(zhuǎn)子系統(tǒng)聯(lián)軸器不對中研究綜述[J].汽輪機技術(shù),2007，49（5）：321-325.LIU Zhansheng.Research on the coupling misalignment of rotor system[J].Steam Turbine Technology,2007,49（5）：321-325.（in Chinese）

[5] 劉永福.航空發(fā)動機構(gòu)造[M].北京：國防工業(yè)出版社，1995：194-201.LIU Yongfu.Aviation engine construction[M].Beijing：National Defense Industry Press，1995：194-201.（in Chinese）

[6] 王瑜.公差穩(wěn)健優(yōu)化設(shè)計的研究[J].計算機集成制造系統(tǒng)，2007，13(11)：2081-2085.WANG Yu.Research on robust optimal design of tolerance[J].Computer Integrated Manufacturing System,2007，13(11)：2081-2085.（in Chinese）

[7] 閻艷,余美瓊,王國新,等.平面尺寸鏈公差分析算法研究[J].北京理大學(xué)學(xué)報，2011,31（7）：799-802.YAN Yan,YU Meiqiong,WANG Guoxin,et al.Study on the algorithm of plane dimension chain tolerance analysis[J].Journal of Beijing University of Science and Technology，2011，31（7）：799-802.（in Chinese）

[8] Siva K M，Islam M N，Lenin N,et al.Optimum tolerance synthesis for complex assembly with alternative process selection using Lagrange multiplier method[J].International Journal of Engineering，2010，3（4）：380-402.

[9] Jinn T T.Robust optimal-parameter design approach for tolerance design problems[J].Engineering Optimization,2010，42（12）：1079-1093.

[10] Taguchi G.Taguchi on robust technology development：bringing quality engineering upstream[M].New York：ASME Press，1993：1-15.

[11] 陳立周.穩(wěn)健設(shè)計[M].北京：機械工業(yè)出版社，1999：2-67.CHEN Lizhou.Robust design[M].Beijing：Mechanical Industry Press，1999：2-67.（in Chinese）

[12] Ghanmi S，Bouazizi M L，Bouhaddi N.Robustness of mechanical systems against uncertainties[J].Finite Elements in Analysis and Design，2007，43（9）：715-731.

[13] 程遠勝，鐘玉湘，游建軍.概率及非概率不確定性條件下結(jié)構(gòu)魯棒設(shè)計方法[J].工程力學(xué)，2005，8（4）：10-14.CHENG Yuansheng，ZHONG Yuxiang，YOU Jianjun.Probabilistic and non-probabilistic uncertainties under the condition of robust design of structure[J].Engineering Mechanics，2005，8（4）：10-14.（in Chinese）

[14] Pan W T.A new fruit fly optimization algorithm：taking the financial distress model as an example[J].Knowledge Based Systems，2012，26(2)：69-74.

[15] 周平，白廣忱.基于神經(jīng)網(wǎng)絡(luò)與果蠅優(yōu)化算法的渦輪葉片低循環(huán)疲勞壽命健壯性設(shè)計[J].航空動力學(xué)報，2013，28（5）：1013-1018.ZHOU Ping，BAI Guangchen.Based on neural network in drosophila and optimization algorithm of the turbine blade low cycle fatigue life of robust design[J].Journal of Aerospace Power，2013，28（5）：1013-1018.（in Chinese）

Robust Optimization Design of Aeroengine Mating Spline Structure

LIU Hong-lei，CHEN Zhi-ying，ZHOU Ping
（School of Energy and Power Engineering，Beihang University，Beijing 100191 China）

In order to design the matching relationship of spline structure,improve the matching precision,adapt to the misalignment of the working conditions and reduce the sensitivity of the mismachining tolerance,the characteristics of spline structure of the aircraft engine casing in the dip angle and the matching relationship were studied.The calculation formula of the tooth flank clearance and the related characteristics of the spline structure were derived.On this basis,the optimal design of the algorithm for the optimization of the tooth flank clearance of the sleeve was designed.Results show that the calculation formula considered angle influence can effectively avoid the movement interference and assembly prestressing force,and robust design and optimization can effectively reduce tooth side gap due to the influence of random parameters such as angle misalignment and mismachining tolerance.It is proved that the calculation and optimization method can make the design more reasonable,which can provide reference for the design and assembly of spline structure.

spline structure;angle misalignment;tooth side gap;fruit fly optimization algorithm;robust optimization;aeroengine

V 232.9

A

10.13477/j.cnki.aeroengine.2016.06.007

2016-05-21基金項目：國家自然科學(xué)基金（51275024）資助

劉宏蕾（1990），男，碩士，研究方向為航空發(fā)動機結(jié)構(gòu)強度可靠性和動態(tài)裝配；E-mail：15933906068@163.com。

劉宏蕾，陳志英，周平.航空發(fā)動機套齒結(jié)構(gòu)配合穩(wěn)健性優(yōu)化設(shè)計[J].航空發(fā)動機，2016,42（6）：44-50.LIUHonglei，CHENZhiying，ZHOUPing.Robustoptimizationdesignofaeroenginematingsplinestructure[J].Aeroengine，2016,42（6）：44-50.

（編輯：栗樞）