遲建華，王承剛，田志濤，李培宏

(1.大連開發(fā)區(qū)大開產(chǎn)業(yè)發(fā)展研究中心，遼寧 大連　116600；2.大連運明自動化技術(shù)有限公司，遼寧 大連　116600；3.大連機床集團有限責任公司，遼寧 大連　116620)

?

擺線針輪行星傳動多目標穩(wěn)健優(yōu)化設(shè)計

遲建華1，王承剛2，田志濤2，李培宏3

(1.大連開發(fā)區(qū)大開產(chǎn)業(yè)發(fā)展研究中心，遼寧 大連116600；2.大連運明自動化技術(shù)有限公司，遼寧 大連116600；3.大連機床集團有限責任公司，遼寧 大連116620)

摘要：采用粒子群優(yōu)化算法對擺線針輪行星傳動進行約束多目標優(yōu)化設(shè)計，求出滿足約束條件的Pareto解集?？紤]工程實際中誤差的客觀存在，進行了7因素9水平的均勻試驗設(shè)計及望目特性的信噪比計算，在此基礎(chǔ)上，依據(jù)多目標穩(wěn)健優(yōu)化設(shè)計原理在Pareto解集中選擇最穩(wěn)健的設(shè)計做為最終方案。分析了結(jié)果方案中目標函數(shù)(體積函數(shù)和修形誤差函數(shù))對各個設(shè)計變量的靈敏度，得出了目標函數(shù)最為敏感的因素及各因素對目標函數(shù)的影響程度用于指導公差的制定和實際加工制造。

關(guān)鍵詞：擺線針輪；行星傳動；優(yōu)化設(shè)計

0引言

目前，擺線針輪行星傳動在國內(nèi)外都被積極研究發(fā)展。文獻[1-2]利用普通平面機構(gòu)的速度瞬心的三心定理，通過針齒上嚙合點的位置和坐標變換得到擺線輪的齒廓曲線方程。文獻[3]根據(jù)曲面單參數(shù)的包絡(luò)方法建立了擺線針輪傳動的嚙合方程。文獻[4]根據(jù)微分幾何和齒輪嚙合原理，由針齒及給定的運動，建立行星輪共軛嚙合齒廓的通用方程，并系統(tǒng)分析了擺線針輪嚙合傳動的特性。文獻[5]在對擺線針輪行星傳動中擺線輪的修形方式進行分析的基礎(chǔ)上，提出了最佳修形齒廓的概念，對采用“負移距+正等距”修形方法獲得最佳修形齒廓原理進行了較深入地探討。有關(guān)擺線針輪行星傳動單目標優(yōu)化設(shè)計的論文現(xiàn)已多見,雖偶見有雙目標優(yōu)化設(shè)計的論文，仍未擺脫傳統(tǒng)多目標優(yōu)化方法的框架。他們的研究成果在不同程度上進一步揭示了擺線針輪行星傳動的內(nèi)在規(guī)律，對于提高擺針傳動設(shè)計質(zhì)量和加工質(zhì)量都具有重要的意義。但迄今為止，利用現(xiàn)代優(yōu)化技術(shù)對擺線針輪行星傳動進行多目標優(yōu)化設(shè)計的研究尚未見有論文發(fā)表。

本文將多目標粒子群優(yōu)化算法、穩(wěn)健設(shè)計、均勻試驗設(shè)計及靈敏度分析相結(jié)合，通過擺線針輪行星傳動算例試圖找到一條解決復雜約束、高維且有誤差情況下機械傳動多目標優(yōu)化設(shè)計的有效途徑。

1擺線針輪行星傳動的多目標優(yōu)化設(shè)計

1.1多目標粒子群優(yōu)化技術(shù)

擺線針輪行星傳動優(yōu)化設(shè)計是同時追求體積最小和修形誤差最小的兩目標約束優(yōu)化問題。兩目標函數(shù)的量綱不同，不能用非現(xiàn)代的加權(quán)系數(shù)的方法整合成一個目標函數(shù)來優(yōu)化。粒子群優(yōu)化技術(shù)(Particle Swarm Optimization,PSO)在這里得到了成功的應用。文獻[6-8]給出了無約束條件下兩目標粒子群尋優(yōu)原則。本文在此基礎(chǔ)上，給出了有約束條件下的尋優(yōu)程序框圖如圖1所示。該方法是在可行域內(nèi)先初始化一個粒子群，然后通過優(yōu)化問題中的各個目標函數(shù)來共同指導粒子在可行域內(nèi)向著各個目標函數(shù)不同時增大的方向運行，最終獲得非劣最優(yōu)解集。

圖1　尋優(yōu)程序框圖

1.2應用實例

(1)設(shè)計原始數(shù)據(jù)：功率P=30kW，輸入轉(zhuǎn)速nH=1500r/min，傳動比i=25，擺線輪、針齒銷、針齒套、柱銷及柱銷套均采用軸承鋼GCr15，表面硬度60HRC，許用接觸應力σHp=1200MPa，針齒銷的許用彎曲應力σFp=1200MPa。

(2)選取8個設(shè)計變量：

X=[X1,X2,X3,X4,X5,X6,X7,X8]T=

[rp,bc,K1,K2,Δrrp,Δrp,Zp,δ]T

式中：rp—針齒中心圓半徑；

bc—擺線輪齒寬；

K1—短幅系數(shù)；

K2—針徑系數(shù)；

Δrrp—正等距修形量；

Δrp—正移距修形量；

Zp—針輪齒數(shù)；

δ—轉(zhuǎn)角修形量。

xci—轉(zhuǎn)角修形擺線輪齒形第i個點的坐標

Zc—擺線輪齒數(shù)

φ—轉(zhuǎn)臂相對于某一針齒中心矢徑的轉(zhuǎn)角

a—中心距(偏心距)

(4)確定約束條件：依據(jù)文獻[9]和工程經(jīng)驗確定齒面接觸疲勞強度條件、針齒銷抗彎強度條件、剛度條件、限定設(shè)計變量的上下界、控制同時嚙合有效傳力齒數(shù)、修形量約束及不產(chǎn)生頂切約束條件。

(5)程序設(shè)計：用MATLAB進行程序設(shè)計計算，循環(huán)100次，結(jié)果得到100組可行解及相應的目標函數(shù)值。

2穩(wěn)健設(shè)計

采用多目標粒子群優(yōu)化算法得到非劣解集以后，還必須考慮誤差的客觀存在，對眾多的非劣解采用均勻試驗設(shè)計計算信噪比，再依據(jù)信噪比決策出一組對誤差最不敏感的穩(wěn)健解作為最終方案[10]。

2.1均勻試驗設(shè)計

均勻試驗設(shè)計是考慮如何將設(shè)計點均勻地散布在試驗范圍內(nèi)，使得能用較少的試驗點獲得最多的信息。

表1　7因素9水平均勻設(shè)計表

2.2望目特性的信噪比計算

在實際應用中，望目特性的信噪比按下式計算[10]：

(1)

n—樣本點的數(shù)量

對于多屬性的設(shè)計方案，如方案具有[f1(X),f2(X),…fm(X)]m個屬性。對m個屬性分別求其望目特征的信噪比，得[SN1,SN2,…SNm]，則該方案的信噪比為：

SNA=Min[SN1,SN2,…SNm]

(2)

2.3穩(wěn)健設(shè)計最終方案

信噪比SN*=53.1566

設(shè)計變量

X*=[rp,bc,K1,K2,Δrrp,Δrp,Zp,δ]T=

[205.0000,21.4627,0.8049,1.7891,0.4086,0.2083,30.0000,0.0056]T

體積目標函數(shù)值：2.8336e+006mm3

結(jié)果修形誤差目標函數(shù)值：0.0105°

3靈敏度分析

式中：m—函數(shù)的個數(shù)；n—變量的個數(shù)。

分別求出目標函數(shù)f1(X)和f2(X)在X*處對各設(shè)計變量的偏導數(shù)：

由此可見，總體積目標函數(shù)f1(X)隨針齒中心圓半徑rp和擺線輪齒寬bc單調(diào)增加，f1(X)對擺線輪齒寬bc比針齒中心圓半徑rp更為敏感。修形誤差函數(shù)f2(X)對轉(zhuǎn)角修形量δ極為敏感，f2(X)隨δ的微小增加而急劇增加。因此，必須嚴格控制δ的公差，且δ應為負偏差，以使正移距優(yōu)化組合的齒形曲線與轉(zhuǎn)角修形齒形曲線在工作段上盡量吻合，保證傳動的精確性和平穩(wěn)性。對f2(X)影響次之的因素是短幅系數(shù)K1，f2(X)隨K1的增加而下降，因此K1應制定正偏差。對f2(X)影響再次之的因素是正移距修形量Δrp，f2(X)隨Δrp的增加而單調(diào)增加，因此Δrp應制定負偏差。對f2(X)影響最小的因素是針齒中心圓半徑rp，其次是針徑系數(shù)K2。另外還可以看到，f2(X)隨正等距修形量Δrrp的增加而單調(diào)下降，因此Δrrp應制定正偏差。

4結(jié)論

(1)將多目標粒子群優(yōu)化算法的最優(yōu)解評估選取原則應用于復雜約束、高維多目標尋優(yōu)求解中，使得所設(shè)計的擺線針輪行星傳動即能夠總體積最小，又能夠修形誤差最小，使設(shè)計的產(chǎn)品即有較低的經(jīng)濟成本又有較高的加工質(zhì)量。

(2)從多目標優(yōu)化結(jié)果的非劣解集中選擇最終方案時，考慮在生產(chǎn)過程中的誤差，首先采用均勻試驗設(shè)計計算各目標的信噪比，再選擇信噪比最大(最穩(wěn)健)的設(shè)計方案作為最終方案。算例驗證了采用基于信噪比的多目標決策方法的必要性和有效性。研究還表明，容差越大，非劣解集(備選方案)中真正的可行解越少。因此，為了使設(shè)計方案更穩(wěn)健可靠，可適度取稍大點兒的容差。

(3)通過對擺線針輪行星傳動進行靈敏度分析，得到了總體體積和修形誤差對各設(shè)計變量的敏感程度。對利用靈敏度分析制定公差、指導機械工程實踐進行了嘗試性工作。

[參考文獻]

[1] Shin J H,Kwon S M.On the lobc profile design in a cycloid reducer using instant velocity center[J]. Journal of Mech Theory,2006,41:596-616.

[2] Shin J H,Chang S, Kwon S M,et al.New shape design method of an epicycloidal gear for an epicycloid drive[J].Proceeding of SPIE,2005,6040:60401G1-G6.

[3] Lai T S.Geometric design of roller drive with cylindrical meshing elements[J]. Journal of Mech Theory,2005,40:55-67.

[4] 陳兵奎，房婷婷，李朝陽，等.擺線針輪行星傳動共軛嚙合理論[J].中國科學E輯：科學技術(shù)，2008，38(1)：148-160.

[5] 關(guān)天民.擺線針輪行星傳動中擺線輪最佳修形量的確定方法[J].中國機械工程，2002，13(10)：811-814.

[6] 魏靜萱.解決單目標和多目標優(yōu)化問題的進化算法[D].西安:西安電子科技大學，2009.

[7] 張敏.約束優(yōu)化和多目標優(yōu)化的進化算法研究[D].合肥:中國科學技術(shù)大學，2008.

[8] 張利彪，周春光，馬銘，等.基于粒子群算法求解多目標優(yōu)化問題[J].計算機研究與發(fā)展，2004，41(7)：1286-1291.

[9] 齒輪手冊編委會.齒輪手冊[M].北京：機械工業(yè)出版社，2005.

[10] 安偉剛.多目標優(yōu)化方法研究及其工程應用[D].西安:西北工業(yè)大學，2005.

[11] 孫志禮，陳良玉.實用機械可靠性設(shè)計理論與方法[M].北京：科學出版社，2003.

[12] 張靜，張柏林，劉永均.基于靈敏度分析的多學科設(shè)計優(yōu)化解耦方法[J].西南交通大學學報，2007,42(5):563-566.

(編輯趙蓉)

Multi-objective Moderate Optimization in a Cycloid Pin Gear Planetary Drive

CHI Jian-hua1,WANG Cheng-gang2,TIAN Zhi-tao2,LI Pei-hong3

(1.Dalian Development Area Dakai Industrial Development Research Center,Dalian Liaoning 116600,China;Dalian Yunming Auto-Technology Co.，Ltd., Dalian Liaoning 116600,China)

Abstract：We could derive to meet the constraints of the Pareto solution set for cycloid pin gear planetary drive by the particle swarm optimization algorithm for constrained multi-objective optimization. As there is objective error in engineering, on the basis of seven factors, nine levels of uniform design and the look mesh characteristics of SNR calculation,We take the most robust design as the final stage from the Pareto which the multi-objective robust design principle is in. By the analysis of the results of the program in the objective function (volume function and the modification error function) to the sensitivity of the various design variables, it concluded that the most sensitive factors of the objective function and their impact on the objective function which can be used to guide the development of tolerance and the actual processing.

Key words：cycloid pin gear;planetary drive;optimization

中圖分類號：TH132.41；TG506

文獻標識碼：A

作者簡介：遲建華(1957—)，男，山東昌邑人，大連開發(fā)區(qū)大開產(chǎn)業(yè)發(fā)展研究中心工程師，研究方向為機械設(shè)計和工藝，智能制造裝備產(chǎn)業(yè)區(qū)域性發(fā)展戰(zhàn)略，(E-mail)chijianhua7722@163.com。

收稿日期：2015-12-10

文章編號：1001-2265(2016)02-0033-03

DOI:10.13462/j.cnki.mmtamt.2016.02.010