基于MATLAB的中速磨煤機減速箱的優化設計*

郭 慶，龔九洲，吳運國，李 濤

（1.揚州電力設備修造廠，江蘇揚州 225003；2.中國科學院沈陽自動化研究所機器人學國家重點實驗室，遼寧沈陽 110016）

基于MATLAB的中速磨煤機減速箱的優化設計*

郭 慶1，龔九洲1，吳運國1，李 濤2

（1.揚州電力設備修造廠，江蘇揚州 225003；2.中國科學院沈陽自動化研究所機器人學國家重點實驗室，遼寧沈陽 110016）

以中速磨煤機減速箱為物理模型，建立了以體積最小為目標，以齒輪齒數、齒輪模數、齒輪齒寬為設計變量，以齒輪的強度要求等為約束條件的優化設計數學模型。將優化結果與原始方案對比，結果表明體積減少了原始方案的29%，達到了節省材料，結構緊湊的目的。

行星齒輪減速箱；優化設計；傳動比

0 前言

火電廠的HP1003型中速磨煤機采用了KMP300型行星齒輪減速箱。此減速箱主要是由一對傘齒輪和一組行星輪系組成。中速磨煤機通常情況下是立式布置，要求其減速箱在保證承載力和使用壽命的前提下結構緊湊。因此，按最小體積為目標對減速箱進行最優化設計，不僅可縮小體積，而且對減小質量、節約材料及降低成本等都是很有實效的，所以中速磨煤機減速箱的優化設計對其整機設計具有重要意義。在優化過程中，首先要對系統進行分析，建立分析模型，然后進行優化[1-3]。本文以最小體積作為目標函數，建立了優化設計數學模型。圖1所示為建立優化模型的具體過程。

1 優化設計數學模型

HP1003型減速箱的運動簡圖如圖2所示。其中，Z1、Z2為傘齒輪的齒數；Z3、Z4、Z5分別為太陽輪、內齒圈、行星輪的齒數，行星輪個數為3。已知，輸入功率P=520 kW，輸入轉速980 r/min，總傳動比i=29.75，傘齒輪、太陽輪和行星輪的材料均為17CrNiMo6，滲碳淬火處理，HRC58～62；內齒圈的材料為34CrNiMo6，調質處理，HB300～340。

1.1 目標函數

減速箱主要由傘齒輪、太陽輪、行星輪、內齒圈和行星支架等構成。其中傘齒輪、太陽輪、行星輪及內齒圈的體積影響著整個行星減速箱的體積和尺寸。本設計以齒輪體積最小作為設計追求的目標，其模型為[4]：

圖1 建模過程示意圖

式（1）中，F(X)為傘齒輪對、太陽輪、行星輪和內齒圈體積之和，n=3為行星輪個數；m1和b1分別為傘齒輪的模數和齒寬；m2和b2分別為行星輪系的模數和齒寬；V1、D1、Z1分別為傘齒輪1的體積、齒寬中點分度圓直徑和齒數；V2、D2、Z2分別為傘齒輪2的體積、齒寬中點分度圓直徑和齒數；V3、D3、Z3分別為太陽輪的體積、分度圓直徑和齒數；V4、D4、Z4分別為內齒圈的體積、分度圓直徑和齒數；V5、D5、Z5分別為行星輪的體積、分度圓直徑和齒數。S為內齒圈分度圓到外圓的距離，S=hf+δ，其中hf為齒根高，hf=1.25 m，δ為內齒圈壁厚，取δ=m。Dw為內齒圈外圓直徑，Dw= D4+2S。根據行星輪系同心條件，目標函數為：

圖2 減速箱運動簡圖

1.2 設計變量

根據目標函數式，取齒輪齒數Z1-Z5、傘齒輪中點法向模數m1和齒寬b1、行星輪系模數m2和齒寬b2這九個獨立參數為設計變量，即：

1.3 約束條件

（1）傳動比條件。通常傳動比的誤差范圍在6%以內，即|i/i0-1|＜6%。式（3）中，i0、i分別為優化前和優化后的減速箱總傳動比。

（2）同心條件。太陽輪和內齒圈的軸線是重合的，因此行星輪與內齒圈的中心距應等于行星輪與太陽輪的中心距，應滿足Z3+2Z5=Z4。

（3）裝配條件。行星輪均勻分布時，相鄰兩個行星輪之間的夾角為120°，為以防行星輪因齒輪干涉而裝不進的情況，太陽輪齒數Z3和內齒圈齒數Z4應滿足C=(Z3+Z4)/3。式中，C為正整數。

（4）鄰接條件。行星輪系是多個行星輪同時傳動，為使相鄰的兩個行星輪不互相碰撞，要求其齒頂圓間有一定的間隙，應滿足dac＜2α'arcsin（π/3）。

式（6）中，dac為行星輪齒輪頂圓直徑；α'ac為太陽輪和行星輪嚙合副的中心距。

（5）根切限制。17≤Z3。

（6）重合度要求。重合度的大小影響著傳動的承載能力和平穩性，而重合度越大，承載能力越高、平穩性越好，因此行星輪系傳動一般保證重合度大于或等于1.2，約束條件如下：式（8）、（9）中，αa、αb、αc分別是太陽輪、內齒圈和行星輪的齒頂圓壓力角，α'1為太陽輪和行星輪的嚙合角，α'2為內齒圈和行星輪的嚙合角。

（7）齒寬系數要求。5≤b2/m2≤17。

（8）最小模數要求。mmin=2。

（9）接觸強度條件[5]。σH≤[σ]H。

式（14）-（18）中，σH1、σH2、σH3、σH4、σH5分別為傘齒輪對、太陽輪、內齒圈、行星輪的計算接觸應力。

（10）齒根彎曲強度條件。σF≤[σ]F。

式（19）-（23）中，σF1、σF2、σF3、σF4、σF5分別為傘齒輪對、太陽輪、內齒圈、行星輪的計算齒根應力。

2 問題求解與結果分析

由以上分析，可以建立21個約束條件，其數學模型為：

滿足

式（24）、（25）中，gm(x)和hn(x)分別為不等式約束和等式約束。該數學模型屬于具有一個目標函數，9個設計變量，21個約束的非線性優化設計問題[6]。應用MATLAB工具箱，采用序列二次規劃法（SQP法）進行優化[7]。在模型的9個設計變量中，齒數為整型變量；模數應取標準值，是離散實型變量；齒寬考慮到工藝性應取整，是離散整型變量。在整個優化求解過程中，首先將離散型變量視為連續型變量，應用連續型變量的優化方法求解得到連續型優化值之后，然后圓整到最接近的有效離散值，最后把圓整后的方案進行比較確定最優結果。原始設計參數及優化后的設計變量最優值見表1。從表1中可以看出，以減速箱齒輪的體積為目標函數F(X)時，F(X)值較原設計減少29%，這表明優化設計的理論和方法用于生產實際，必將產生較大的經濟效益。

表1 設計參數優化計算結果

3 結論

本文以中速磨煤機減速箱為物理模型，建立了包括內齒圈體積在內的減速箱體積最小的目標函數，使優化數學模型更全面且符合實際，通過討論優化設計中約束條件及優化方法等問題，對其進行了優化設計。優化后的減速箱結構緊湊，體積體積小，節省了材料，降低了成本，且能滿足各項要求。從表1中可以看出，優化后體積減小了29%。本文的工作為減速箱的改進設計提供了理論依據，同時這種方法可以推廣到各種多級減速箱全局優化情況，更符合工程設計的需要，在工程界有著廣闊的應用前景。

［1］劉惟信.機械最優化設計［M］.北京：清華大學出版社，2002.

［2］張凱燕，莫云輝，鄧召義，等.改進遺傳算法的行星齒輪傳動多目標模糊物元可靠性優化［J］.上海大學學報，2007（2）：22-27.

［3］將宏春，孫巖.ZSJ-2800減速器傳動比最優分配［J］.中國重型裝備，2008（1）：21-22.

［4］濮良貴，紀名剛.機械設計［M］.北京：高等教育出版社，2006.

［5］成大先.機械設計手冊［M］.北京：化學工業出版社，2002.

［6］薛履中.工程最優化技術［M］.天津：天津大學出版社，1988.

［7］陳杰.MATLAB寶典［M］.北京：電子工業出版社，2007.

Optimal Design for Reduction Gearbox of Medium Speed Mill Based on MATLAB

GUO Qing1，GONG Jiu-zhou1，WU Yun-guo1，LI Tao2
（1.Yangzhou Electric Power Equipment Manufacture Factory，Yangzhou 225003，China；2.State Key Laboratory of Robotics，Shenyang Institute of Automation，Chinese Academy of Sciences，Shenyang110016，China）

The reduction gearbox of medium speed mill is used as physical model.Mathematical model is established with the minimum volume as the objective function，the number of teeth，the gear modulus and the tooth width as the design variables，the strength of the gear and so on as the constraints.The optimized results are compared with the original data and the results show that the volume is reduced by 29%，which can economize material and realize compact configuration.

planetary gear reducer；optimal design；transmission ratio

TH122

A

1009－9492(2014)07－0055－03

10.3969/j.issn.1009-9492.2014.07.016

郭 慶，男，1984年生，江蘇揚州人，碩士，工程師。研究領域：CAD/CAM技術，數字化設計。

(編輯：向 飛)

*國家863計劃資助項目（編號：2009AA04Z203）

2014－01－08