基于遺傳算法的擺線針輪行星減速器優化設計
摘要:為了解決擺線針輪行星減速器的傳動效率低,結構不夠緊湊的問題,通過推導其傳動效率的公式,并對設計變量和約束條件進行研究,建立了以體積最小、效率最高為目標函數的優化設計數學模型。采用遺傳算法,運用MATLAB軟件編制遺傳優化程序,調用遺傳工具箱進行優化。優化設計結果與傳統的設計相比,體積減小16.6,效率增加2.06%,優化效果顯著,從而為擺線針輪行星減速器優化設計提供一定的理論和方法。
減速器是機械傳動中的核心部件,漸開線減速器具有廣泛的應用,但有些場合卻不能滿足使用要求。擺線針輪行星減速器具有體積小、質量輕、傳動比范圍大、壽命長、精度保持穩定、傳動平穩等一系列的優點而得到了廣泛的應用。但擺線針輪行星減速器的傳動效率卻比漸開線行星減速器的傳動效率低,因此研究擺線針輪行星減速器的參數設計,優化體積和傳動效率具有重要意義。
由于擺線針輪行星減速器涉及的參數比較多,關系比較復雜,影響的因素也較多,這就給擺線針輪行星減速器的優化設計帶來了一定的困難。從現在所發表的擺線針輪行星減速器的優化設計論文[1][2]142[3]來看尚存一定的不足,比如目標函數單一,不能有效提高擺線針輪行星減速器的傳動效率;各約束條件單位不統一,給優化設計的實現造成了一定的難度。因此全面和深入地研究擺線針輪行星減速器的優化設計問題,對擺線針輪行星減速器技術發展有一定的實用意義。本文中我們推導出擺線針行星減速器傳動效率的公式,以體積最小和效率最高作為目標函數,利用遺傳算法對擺線針輪行星減速器進行多目標優化設計。
1 擺線針輪減速器優化設計數學模型
1.1 設計變量
擺線針輪傳動結構如圖1所示。現研究一級擺線針輪行星減速器優化問題,取擺線針輪行星減速器的針輪分布圓直徑DZ;針齒銷直徑dz′;擺線輪寬度B;短幅系數K1;擺線輪內孔直徑D1;柱銷直徑ds′作為設計變量,即
1.2 目標函數
(1)本文以在功率和傳動比一定條件下,并選定輸出機構銷軸數Zs,選擇減速器體積最小為目標函 數[2]142,即
minf1(x)=π4(DZ+dZ′+2Δ1)2(2B+δ) (2)
式中,Δ1為針齒套壁厚,一般取Δ1=2~10mm;δ為擺線輪之間的間隔,δ=b-B,b為轉臂軸承寬度。
(2)以擺線針輪行星減速器的傳動效率為另一個目標函數。
當針輪固定時,轉臂輸入時(針輪輸出時,轉臂輸入時類似),擺線針輪行星減速器的傳動效率η 可近似用下式計算[4]330。
η=ηxcηBηWηB′2 (3)
式中,ηxc為當針輪固定,減速部分的嚙合效率,ηxc=1-0.04πZcK1(1-dZ′DZ)1+0.04πK1(1-dZ′DZ);ηB為轉臂軸承的效率,其概略值為,ηB=0.99~0.995;ηW為輸出機構的效率,ηW =1- 0.08DZK1ds′π[(DZ-dz′-2Δ1-K1DZ/Zb)+D1](d′s+2Δ2);ηB′為輸入或輸出軸上的滾動軸承的嚙合效率,其概略值為ηB′=0.99~0.999。
由文獻[4]324-330經整理后,其效率計算式為
考慮到η∈(0~1),將效率最大值轉換成最小值,則效率最大的目標函數是
式中,Δ2為柱銷套壁厚,一般取Δ2=2.5~5mm;Zb為針輪齒數,Zb=i+1;Zc為擺線輪齒數Zc=i。
1.3 約束條件的建立
(1)擺線輪齒廓不根切條件的約束條件
式中,Ds為柱銷孔中心分布圓直徑,Ds=dfc+D12;dfc為擺線輪齒根圓直徑,dfc=DZ-(dz′+2Δ1)-2 K1RZ/Zb,所以Ds=[DZ-(dz′+2Δ1)-K1DZ/Zb+D1]/2;dsk為銷軸孔直徑,dsk=ds+2A,ds為銷軸套外徑,ds=ds′+2Δ2,A 為偏心距A =K1DZ/2Zb,所以dsk=ds′+2Δ2+K1DZ/Zb;Δ′為銷軸孔與擺線輪內孔之間所留的一定厚度,Δ′=0.03DZ。
所以推得
(4)擺線輪與針齒的接觸強度條件
由文獻[5]74-82推導整理得
式中,MV為輸出軸的阻力矩,MV =97 360ZcNHnH,NH為輸入功率(kW);nH為主軸轉速(r/min);Y2max為最大接觸應力處的位置系數(由文[5]77表10查取);[σj]為許用接觸應力(MPa)。得
(5)針齒銷的彎曲強度條件
對于雙支點針齒銷的彎曲應力,由文獻[5]82-85推導整理得
故約束條件為
式中,[σW ]為針齒銷許用彎曲應力(MPa);δ′為擺線輪于針齒殼側面見的間隙(mm),一般取δ′=4mm;δ為擺線輪間隔環的厚度(mm),δ=b-B,b為軸承的寬度;Δ為針齒殼側面的壁厚(mm),一般為B≥Δ≥dZ′。
(6)柱銷的彎曲強度限制條件為
W 機構中的柱銷所受的最大彎曲應力由文獻[5]85-87推導得
(8)短幅系數K 應在0.45~0.8范圍選取,即
g10(X)=0.45-K≤0 (21)
g11(X)=K-0.8≤0 (22)
(9)擺線輪壁厚的約束條件[4]318
擺線輪的厚度一般為B=(0.05~0.1)DZ,故
g12(X)=0.05DZ-B≤0 (23)
g13(X)=B-0.1DZ≤0 (24)
綜上所述,擺線針輪減速器的多目標優化設計數學模型可表示為
2 遺傳優化
由于遺傳優化具有魯棒性好、全局尋優的特點[6],為了避免陷入局部最優解,本文采用遺傳算法對擺線針輪行星減速器進行多目標優化設計。
2.1 求解多目標優化問題的遺傳算法
遺傳優化中為了綜合考慮各分目標函數在相對重要程度方面的差異和量綱上的差異,使各個分目標函數能均勻一致地趨向于各自的最優解,采用加權法[7]26。其目標函數表示為
F(x)=w1f1(x)+w2f2(x) (26)
綜上所述,擺線針輪行星傳動是多目標函數的優化問題。由于體積和效率兩目標函數的量綱不同,值的大小相差很大,構造統一目標函數時難以充分反映各子目標函數的重要程度,因此先將體積和效率的目標函數作無量綱處理[8]。由于在擺線針輪行星減速器中體積和效率具有同等的地位,所以我們將各單目標優化值的倒數取作權重系數。
ωi=1/fi*(x)fi*(x)=minfi(x) (27)
式中,ωi是作了無量綱處理后的第i個子目標函數的 權重。
經過這樣的處理后,使各子目標函數的量綱得到了統一,在同一量綱的基礎上再對各子目標函數賦予相應的權數,最后就將多目標優化轉化為了單目標優化。按照此方法,評價函數不僅不會受到各分目標函數值相對大小的影響,還能充分反映出各單目標函數值偏離各自最優值的程度。
2.2 遺傳優化算法的實現
本文采用MATLAB 語言編制遺傳算法程序[7]26[9],利用MATLAB7.0中的遺傳算法工具箱中遺傳優化函數ga對擺線針輪行星減速器進行遺傳優化。先求得f1(x)、f2(x)分別為單目標函數時的最小值,然后再求多目標的最小值F(x)。
3 優化實例及結果分析
3.1 實例
以某公司生產的B12型減速機為設計實例,減速機的輸入額定功率P1=11kW,輸入轉數n1=1440r/min,傳動比i=29,輸出機構環形槽數Zs=10,擺線輪和針齒的材料均采用GCr15軸承鋼,HRC58~62,許用接觸應力[σj]=1 100Ma,許用彎曲應力[σW ]=200MPa。
3.2 優化求解
采用MATLAB 遺傳工具箱進行遺傳優化,具體步驟如下:
①編寫以體積最小為單目標函數的M文件,通過調用遺傳優化函數ga,求得最小體積為2 044 181;②編寫以效率最高為單目標函數的M文件,通過調用遺傳優化函數ga,求得最高效率為95.36%;③編寫以體積最小和效率最高為多目標函數的M 文件,如圖2所示,其中體積的權重為1/2 044 181,效率的權重為1/ (1-0.953 6);④調用遺傳算法主函數ga.m。
計算結果表明:采用遺傳算法優化設計要優于常規設計方法。在保證主要傳動元件設計要求的同時,體積減小16.6%,效率增加2.06% 。所以,利用遺傳算對擺線針輪行星減速器進行優化設計效果是顯著的。
4 結論
分析影響擺線針輪行星減速器的主要參數,建立以體積最小和效率最高為優化設計的目標函數,并系統研究了優化的約束函數。用MATLAB語言編制了遺傳算法程序,運用MATLAB遺傳工具箱對其進行遺傳優化設計,在滿足強度條件的前提下,獲得比傳統方法更好的結果,使其結構更緊湊,傳動效率更高。可為擺線針輪行星減速器設計提供理論支撐。