直齒錐齒輪分岔脫嚙特性參數(shù)解域界結(jié)構(gòu)

田亞平, 楊江輝, 王瑞邦

(蘭州交通大學(xué) 機(jī)電工程學(xué)院, 蘭州 730070)

0 引 言

錐齒輪傳動(dòng)廣泛應(yīng)用于機(jī)車、航空、機(jī)床等行業(yè)傳動(dòng)裝置中．含間隙和時(shí)變因素的錐齒輪傳動(dòng)系統(tǒng)的非線性振動(dòng)對(duì)系統(tǒng)運(yùn)動(dòng)的穩(wěn)定性提出了嚴(yán)峻挑戰(zhàn),成為學(xué)者們研究的熱點(diǎn)．王三民等[1]建立了正交錐齒輪系統(tǒng)動(dòng)力學(xué)模型并分析了系統(tǒng)分岔和幅值跳躍特性．黃康等[2]考慮隨機(jī)因素探索了間隙的合理范圍．Wang等[3]考慮含齒面混合潤(rùn)滑研究了表面粗糙度對(duì)齒面沖擊的影響．Hua等[4]建立了軸承-錐齒輪轉(zhuǎn)子有限元模型,分析了軸承支承剛度對(duì)動(dòng)態(tài)嚙合力的影響．Cao等[5]分析了不同潤(rùn)滑狀態(tài)下錐齒輪傳動(dòng)系統(tǒng)的非線性動(dòng)態(tài)特性．李飛等[6]研究了摩擦和重合度對(duì)系統(tǒng)振動(dòng)幅值和嚙合力的影響．

小波法[7]、數(shù)值法[8]是求解非線性系統(tǒng)分岔特性的常用方法,但在全局分岔特性轉(zhuǎn)遷特性分析上適用性較差．胞映射理論是研究非線性動(dòng)態(tài)特性轉(zhuǎn)遷的常見(jiàn)方法．在參數(shù)域平面內(nèi),Gou等[9]采用胞映射研究了齒輪副模型的周期吸引子分形結(jié)構(gòu),Liu等[10]通過(guò)對(duì)齒輪解域界分析實(shí)現(xiàn)了混沌控制,林何等[11]分析了弧齒錐齒輪傳動(dòng)系統(tǒng)中周期運(yùn)動(dòng)的解域界結(jié)構(gòu)．因偽不動(dòng)點(diǎn)周期求解及追蹤C(jī)PNF法求解效率高,在含間隙的行星輪系[12]和單級(jí)齒輪系統(tǒng)[13]的周期求解中獲得了成功應(yīng)用．田亞平等[14]用CPNF法研究了單級(jí)齒輪系統(tǒng)的動(dòng)態(tài)特性解域界結(jié)構(gòu)．多參耦合下的解域界結(jié)構(gòu)能有效揭示其動(dòng)態(tài)特性間的耦合及轉(zhuǎn)遷規(guī)律．目前,針對(duì)含間隙直齒錐齒輪動(dòng)態(tài)特性參數(shù)解域界結(jié)構(gòu)的研究還鮮有報(bào)道．

據(jù)統(tǒng)計(jì),間隙和時(shí)變參數(shù)是導(dǎo)致動(dòng)車齒輪傳遞裝置(圖1)故障和車輛振動(dòng)超標(biāo)的主要因素．本文以該錐齒輪傳動(dòng)系統(tǒng)為研究對(duì)象,建立了7自由度錐齒輪系統(tǒng)動(dòng)力學(xué)模型,采用CPNF法求解周期、脫嚙、齒面沖擊、動(dòng)載特性的轉(zhuǎn)遷規(guī)律,探索參數(shù)組合下系統(tǒng)的周期分岔與齒面沖擊及脫嚙、動(dòng)載間的耦合關(guān)系,為錐齒輪結(jié)構(gòu)設(shè)計(jì)參數(shù)選擇提供理論依據(jù)．

1 動(dòng)力學(xué)模型

含彈性支承的直齒錐齒輪傳動(dòng)系統(tǒng)動(dòng)力學(xué)模型如圖2所示[1]．兩齒輪軸線正交于原點(diǎn)O,建立坐標(biāo)系∑:{Oxyz}．兩輪支承等效在齒寬中點(diǎn)O1,O2,沿坐標(biāo)軸的支承剛度和阻尼為kij,cij(i=x,y,z;j=1,2)．ch,kh(t)和bh分別為齒輪副嚙合的阻尼、時(shí)變剛度和齒側(cè)間隙．設(shè)主動(dòng)輪1受到的驅(qū)動(dòng)力矩Tp由常量Tpm和變量Tpv構(gòu)成;從動(dòng)輪2受到穩(wěn)定的阻抗力矩Tgm．忽略齒輪的擺動(dòng)與彎曲,兩個(gè)剛性錐盤的振動(dòng)是分別以O(shè)1,O2為中心,沿坐標(biāo)軸的移動(dòng)和軸線的轉(zhuǎn)動(dòng),即8個(gè)自由度{X1,Y1,Z1,X2,Y2,Z2,θ1,θ2}．

考慮振動(dòng)和誤差,兩齒輪沿嚙合線的動(dòng)態(tài)綜合傳遞誤差Λn為

Λn=(X1-X2)a1-(Y1-Y2)a2-(Z1+r1θ1-Z2-r2θ2)a3-en(t),

(1)

式中,a1=cosδ1sinαn,a2=cosδ1cosαn,a3=cosαn,δ1為主動(dòng)輪節(jié)錐角,αn為法面壓力角,r1,r2為兩輪的節(jié)圓半徑．齒輪副的靜態(tài)綜合誤差en(t)展開成級(jí)數(shù)形式為

(2)

式中,Ωh,Al和Φl為嚙合頻率、l階諧波幅值和初相位．

齒輪副的嚙合力Fn及其沿坐標(biāo)軸的分力Fx,Fy和Fz分別為

(3)

用文獻(xiàn)[15]計(jì)算kh(t)并用級(jí)數(shù)展開為

(4)

式中,km,kkl和Φkl為嚙合剛度的均值、l階諧波幅值和初相位．結(jié)合動(dòng)車齒輪箱工程計(jì)算精度要求,本文選取一階級(jí)數(shù)進(jìn)行數(shù)值仿真計(jì)算．

間隙函數(shù)f(Λn)為

(5)

圖2 動(dòng)力學(xué)模型Fig. 2 The nonlinear dynamic model for a spiral bevel gear set

根據(jù)Newton第二定律,圖2所示的錐齒輪傳動(dòng)系統(tǒng)振動(dòng)方程為

(6)

式中,m1,m2,J1,J2為兩輪的質(zhì)量和轉(zhuǎn)動(dòng)慣量．

用沿嚙合線的動(dòng)態(tài)綜合傳遞誤差Λn作為新自由度,將式(6)中兩齒輪扭振方程合并得

(7)

將上述方程整理并量綱一化處理得

(8)

式中

xj=Xj/bn,yj=Yj/bn,zj=Zj/bn,λ=Λn/bn,Ωn=(km/me)1/2,

Ωij=(kij/mj)1/2,ξij=cij/(2mjΩn),kij= (Ωj/Ωn)2,i=x,y,z,j=1,2,

(9)

bn為標(biāo)稱間隙．

2 系統(tǒng)動(dòng)態(tài)特性仿真指標(biāo)

2.1 脫嚙占空比和齒背嚙合比

輪齒脫嚙和齒背接觸行為是引起齒面拍擊和附加動(dòng)載荷的主要因素．由文獻(xiàn)[16]知,齒輪嚙合的沖擊和脫嚙(I=0,1,2)狀態(tài)可由運(yùn)動(dòng)周期中的最小嚙合綜合誤差λmin與齒側(cè)間隙b判斷,其判斷式為

(10)

脫嚙占空比(non-meshing duty coefficient,δNMDC)和齒背嚙合比(back-meshing duty coefficient,δBMDC)是量化系統(tǒng)脫嚙沖擊嚴(yán)重性的主要指標(biāo),設(shè)在一個(gè)運(yùn)動(dòng)周期T內(nèi)脫嚙時(shí)間為t1,齒背嚙合時(shí)間為t2,則δNMDC和δBMDC為

(11)

2.2 動(dòng)載系數(shù)

齒輪強(qiáng)度設(shè)計(jì)中考慮齒輪振動(dòng)導(dǎo)致動(dòng)態(tài)嚙合力瞬態(tài)增大的因數(shù)稱為動(dòng)載系數(shù)(dynamic load coefficient,δDLC)．受支承剛度、箱體結(jié)構(gòu)和重合度等因素影響動(dòng)載系數(shù)計(jì)算十分復(fù)雜,為突出間隙和時(shí)變剛度對(duì)動(dòng)載系數(shù)的影響,定義運(yùn)動(dòng)周期內(nèi)最大嚙合力|fn(τ)|max與平均值fpm之比的平方根為動(dòng)載系數(shù):

(12)

3 參數(shù)解域結(jié)構(gòu)求解

參數(shù)解域結(jié)構(gòu)求解是基于胞映射原理將參數(shù)平面進(jìn)行胞元離散,在胞元內(nèi)對(duì)周期運(yùn)動(dòng)采用CPNF法延續(xù)追蹤判穩(wěn)獲得其周期運(yùn)動(dòng)的參數(shù)域界結(jié)構(gòu)．在追蹤過(guò)程中伴隨有脫嚙、沖擊和動(dòng)載特性的參數(shù)解域結(jié)構(gòu)求解．其基本思想及實(shí)施過(guò)程如下．

3.1 參數(shù)域胞元化

對(duì)多參γi控制的非自治動(dòng)力學(xué)系統(tǒng)可描述為

(13)

式中,X為n維狀態(tài)向量,τ為量綱一時(shí)間,γi∈RK為系統(tǒng)的k個(gè)控制參數(shù),f為映射法則．

為分析激勵(lì)參數(shù)γi與響應(yīng)間的關(guān)聯(lián)機(jī)制,在[xl,xu]×[yl,yu]二維參數(shù)域Θ內(nèi)對(duì)參數(shù)離散化為n×m個(gè)胞元,參數(shù)胞為cij(i=1,2,…,n;j=1,2,…,m)．其中,xl,yl和xu,yu分別表示x,y方向的起點(diǎn)和終點(diǎn),hi,hj為離散胞在x,y方向的胞元尺度．x,y方向的胞尺度hi,hj為

(14)

在參數(shù)域Θ內(nèi),從xl,yl開始得到x,y方向上參數(shù)的遞推關(guān)系xi,yj為

(15)

則域Θ內(nèi)離散后的各胞元(xi,yj)賦值為

(16)

3.2 動(dòng)力學(xué)行為延拓追蹤步驟

參數(shù)域Θ內(nèi)各胞元的動(dòng)力學(xué)行為采用CPNF法追蹤求解．首先沿x方向?qū)Τ跏及?x0,yj)用偽不動(dòng)點(diǎn)PNF法[10-11]進(jìn)行系統(tǒng)運(yùn)動(dòng)周期求解并判穩(wěn); 然后沿x方向按步長(zhǎng)hi進(jìn)行延拓追蹤、用Floquet乘子|λ|max判斷周期運(yùn)動(dòng)的穩(wěn)定性．當(dāng)|λ|max>1時(shí)周期運(yùn)動(dòng)失穩(wěn),在該胞元(xi,yj)處終止周期追蹤判斷分岔類型,并用PNF法重新求解新周期,若求解周期數(shù)超過(guò)設(shè)定值時(shí)認(rèn)定該胞元為混沌或擬周期狀態(tài)(可用Floquet乘子和Lyapunov指數(shù)判定),在(xi+1,yj)胞元處重新求解周期數(shù)并追蹤直至xu結(jié)束;當(dāng)完成x方向追蹤后,沿y方向遞進(jìn)一層繼續(xù)追蹤求解,直至遍歷全域結(jié)束．

3.3 胞元周期解求解過(guò)程

CPNF法周期求解判穩(wěn)的關(guān)鍵是求方程組(17)的穩(wěn)定解

(17)

式中,Φ=DP(X)為系統(tǒng)的轉(zhuǎn)遷矩陣[12], 其狀態(tài)方程的Jacobi矩陣?f(X,τ)/?X在λ≠±b的分段光滑區(qū)內(nèi)為

(18)

式中

A1=[-a1[kx12ξx1],a2[ky12ξy1],a3[kz12ξz1],

a4[kx22ξx2],-a5[ky22ξy2],-a3[kz22ξz2]],

(19)

因在λ=±b處Jacobi矩陣不存在,故結(jié)合計(jì)算精度選取合理步長(zhǎng)的有限差分法替代為

(20)

方程(17)以(Xk,I)為初值積分一個(gè)周期T獲得系統(tǒng)迭代點(diǎn)Xk+1和轉(zhuǎn)遷矩陣DP(Xk)．令Xk=Xk+1作為下次迭代的初始值,迭代到滿足精度要求的不動(dòng)點(diǎn)XF,就是式(13)的m倍周期運(yùn)動(dòng)X(τ)的不動(dòng)點(diǎn)．找到穩(wěn)定周期解后,引入分岔參數(shù)Ω,采用延續(xù)追蹤法進(jìn)行周期追蹤．設(shè)已求得Ω=Ωk(k=0, 1, 2, …)時(shí),周期T的不動(dòng)點(diǎn)為Xk,用Euler積分法得Ω=Ωk+1點(diǎn)對(duì)應(yīng)不動(dòng)點(diǎn)Xk+1的初始不動(dòng)點(diǎn)Xk+1,0的預(yù)測(cè)公式[14]:

(21)

式中,ΔΩ為追蹤步長(zhǎng),GX(Ωk,Xk)=I-DP(Ωk,Xk),GΩ(Ωk,Xk)=PΩ(Ωk,Xk)．PΩ(Ωk,Xk)由打靶法建立的微分方程組(22)求得

(22)

式中

fΩ(X,Ω,τ)=feΩ{2cos(Ωτ)+Ωsin(Ωτ)}[01×131]T+

αsin(Ωτ)[0 -1 0 1 0 1 0 1 0 -1 0 -1 0 -4]T/4．

(23)

以(Xk,Ωk)為初值積分一個(gè)映射周期T獲得以(Xk,PΩ(Ωk,Xk))為追蹤參數(shù)下的新不動(dòng)點(diǎn)．

4 系統(tǒng)動(dòng)態(tài)特性綜合分析與解域界

不失一般性,取量綱一仿真參數(shù):

ξi1=ξi2=0.01 (i=x,y,z),ξh1=ξh2=0.012 5,ξh=0.05,

ki1=ki2=1.0,kh1=kh2= 0.5,kh=1+αcos(Ωτ),

fpm=0.5,fav=0.0,fe= 0.2,b=1.0．

在參數(shù)解域界平面內(nèi),齒面沖擊/運(yùn)動(dòng)周期(I/P)類型用特定的顏色表示,其分岔曲線為解域界．齒輪系統(tǒng)常常穩(wěn)定運(yùn)行于短周期狀態(tài)下,為突出低周期(P<64)分岔特性,在I/P分岔圖中對(duì)周期數(shù)P≥65的長(zhǎng)周期、擬周期(quasi-P)和混沌(chaos)運(yùn)動(dòng)均用“I/N”表示．系統(tǒng)的倍化(PD bif)、鞍結(jié)(NS bif)、檫切[17](G bif)分岔定義參見(jiàn)文獻(xiàn)[14],不再贅述．當(dāng)系統(tǒng)從I/P向I/(3P)狀態(tài)轉(zhuǎn)遷時(shí),系統(tǒng)發(fā)生了齒輪系統(tǒng)典型的周期3分岔(3T bif),依次類推．

時(shí)變嚙合剛度系數(shù)α表征齒輪嚙合重合度,其值表征單雙對(duì)齒嚙合不均勻性程度．頻率比Ω表征系統(tǒng)轉(zhuǎn)速,其值表征系統(tǒng)轉(zhuǎn)速高低．故構(gòu)建Ω×α∈[1.2, 1.7]×[0.05, 0.55]參數(shù)平面(圖3),在平面內(nèi)均勻地選取501×501個(gè)胞元構(gòu)建I/P及相應(yīng)的δNMDC,δBMDC,δDLC解域界結(jié)構(gòu)．

I/P解域界結(jié)構(gòu)(圖3(a))表明:在參數(shù)平面內(nèi),以周期1(P=1,下同)運(yùn)動(dòng)為主,存在P=2,3,9,10,18,32和quasi-P、chaos等運(yùn)動(dòng)狀態(tài);3種齒面沖擊共存,以單邊沖擊為主．在α<0.15區(qū)域內(nèi),隨Ω遞增系統(tǒng)發(fā)生了SN bif和3T bif,使系統(tǒng)處于0/1,1/1,1/3運(yùn)動(dòng)狀態(tài)．在0.15<α<0.3區(qū)域內(nèi),系統(tǒng)通過(guò)PD bif、3T bif、HP bif等分岔方式進(jìn)入混沌運(yùn)動(dòng)．在0.3<α<0.55區(qū)域內(nèi),隨Ω遞增系統(tǒng)發(fā)生了G bif、3T bif、3T chaos和CIC bif 等分岔．在Ω=1.6附近表現(xiàn)出復(fù)雜的動(dòng)力學(xué)行為轉(zhuǎn)遷,即在該胞元內(nèi)動(dòng)力學(xué)行為極其不穩(wěn)定．沿α方向在Ω<1.25區(qū)間內(nèi),系統(tǒng)以P=1運(yùn)動(dòng)的3種沖擊轉(zhuǎn)遷為主．在1.25<Ω<1.5區(qū)間,系統(tǒng)的動(dòng)力學(xué)行為較為穩(wěn)定．當(dāng)Ω>1.5時(shí),系統(tǒng)通過(guò)HP bif、3T bif和CIC bif的形式使系處于周期-混沌-周期運(yùn)動(dòng)狀態(tài)．參數(shù)解域界結(jié)構(gòu)表明齒輪嚙合重合度是影響齒輪沖擊的主要因素,隨α遞增齒面脫嚙、齒背接觸的概率增大．在系統(tǒng)的共振頻率附近多周期、混沌現(xiàn)象加劇,并出現(xiàn)了齒輪傳動(dòng)系統(tǒng)典型的周期3分岔行為．

系統(tǒng)I/P解域界結(jié)構(gòu)對(duì)應(yīng)的齒面脫嚙、齒背接觸和動(dòng)載特性如圖3(b)—3(d)所示．圖3(b)表明,齒輪脫嚙程度隨周期分岔和齒面沖擊的轉(zhuǎn)遷而階躍突變,在1/1解域內(nèi)其脫嚙程度并不一致,在G bif分岔線附近δNMDC達(dá)到極值0.55左右．隨α增大其脫嚙程度加劇,而在1/N和2/N解域內(nèi)δNMDC值并未達(dá)到極值,即多周期或混沌運(yùn)動(dòng)下齒面脫嚙未必最嚴(yán)重．圖3(c)表明,在雙邊沖擊的2/1域界內(nèi)其齒背接觸最為嚴(yán)重達(dá)到極值0.38左右,在2/3、2/N域界內(nèi)其齒背接觸狀態(tài)隨α遞增而加劇但沒(méi)有達(dá)到極值,即齒輪嚙合的雙邊沖擊程度與系統(tǒng)混沌現(xiàn)象并無(wú)必然聯(lián)系．圖3(d)表明,系統(tǒng)的δDLC隨齒面沖擊和周期運(yùn)動(dòng)分岔出現(xiàn)了突變,在同一解域內(nèi)其動(dòng)載系數(shù)隨α遞增而增大,在無(wú)沖擊的0/1解域內(nèi)其值較小、在雙邊沖擊的2/1、2/N域界內(nèi)達(dá)到極值2.8．

(a) I/P(b) δNMDC

(c) δBMDC(d) δDLC圖3 Ω×α平面動(dòng)態(tài)特性解域界結(jié)構(gòu)Fig. 3 Dynamic properties solution domain boundary structures in the Ω×α 2-parameter plane

為驗(yàn)證CPNF法求解的正確性,沿α方向以0.05等距選取截面,用Runge-Kutta法對(duì)Ω參數(shù)分岔進(jìn)行數(shù)值計(jì)算獲得三維I/P分岔圖(圖4)．

圖4 Ω×α雙參分岔圖Fig. 4 The bifurcation diagram in the Ω×α 2-parameter plane

圖4中,黑、紅、藍(lán)色映射點(diǎn)分別對(duì)應(yīng)無(wú)沖擊、單邊沖擊和雙邊沖擊,其I/P分岔過(guò)程與圖3(a)相匹配驗(yàn)證了解界域結(jié)構(gòu)的正確性．因CPNF法對(duì)周期運(yùn)動(dòng)求解追蹤過(guò)程不用計(jì)算瞬態(tài),故較Runge-Kutta法求解效率高．

為更清晰地展示周期分岔與齒面脫嚙的關(guān)聯(lián)關(guān)系,取α=0.20,0.40的分岔圖如圖5所示．圖5(a)中的系統(tǒng)通過(guò)SN bif、3T bif、HP bif和激變(CIC bif)分岔方式使系統(tǒng)處于0/1, 1/1, 2/3, 1/3, quasi-P,2/chaos,1/chaos,1/1運(yùn)動(dòng)．圖5(b)中系統(tǒng)通過(guò)SN bif、3T bif、PD bif和C bif分岔使系統(tǒng)處于2/1,1/1,2/3,2/chaos運(yùn)動(dòng)狀態(tài)．HP bif和SN bif可通過(guò)CPNF追蹤中Floquet乘子|λ|max來(lái)判定,而G bif與SN bif是通過(guò)齒面沖擊類型和Floquet乘子共同判別,若穿越λ=±b截面時(shí)Floquet乘子從實(shí)數(shù)軸方向穿越單位圓則發(fā)生了鞍結(jié)分岔,否則為擦切分岔．

取α=0.20截面(圖5(a))的部分Poincaré映射和相圖(圖6)驗(yàn)證CPNF求解的分岔過(guò)程．Ω=1.531時(shí)系統(tǒng)為1/1運(yùn)動(dòng)(圖6(a))經(jīng)3T分岔為2/3運(yùn)動(dòng)(圖6(b)),P=3轉(zhuǎn)遷為3個(gè)焦點(diǎn)(圖6(c))、Hopf圈(圖6(d))、扭曲鎖相(圖6(e))直至形成混沌運(yùn)動(dòng)(圖6(f))．即圖5驗(yàn)證了通過(guò)周期3分岔和Hopf分岔使系統(tǒng)處于1/1,2/3,1/1,擬周期和混沌運(yùn)動(dòng)．

(a) α=0.20(b) α=0.40圖5 分岔圖Fig. 5 The bifurcation diagram via α

(a) Ω=1.531(b) Ω=1.535(c) Ω=1.563

(d) Ω=1.575(e) Ω=1.600(f) Ω=1.650圖6 Poincaré映射圖和相圖Fig. 6 Poincaré and phase maps

圖7為圖5(b)中α=0.40時(shí)I/P部分的Poincaré映射和相圖．Ω=1.523時(shí)系統(tǒng)處于0/1運(yùn)動(dòng)(圖7(a))經(jīng)周期3分岔為2/3運(yùn)動(dòng)(圖7(b)),P=3再經(jīng)周期3分岔為P=9(圖7(c))、P=18(圖7(d))、……、混沌運(yùn)動(dòng)(圖7(e))．圖7驗(yàn)證了系統(tǒng)經(jīng)周期3分岔進(jìn)入混沌的過(guò)程．

取α=0.40截面(圖5(b))對(duì)應(yīng)的齒面脫嚙比、齒背接觸比和動(dòng)載系數(shù)分布圖(圖8)．齒面脫嚙在齒面沖擊轉(zhuǎn)遷截面發(fā)生了突變,在1/1周期區(qū)間其脫嚙程度隨Ω遞增而減弱(圖8(a))．齒背接觸在雙邊沖擊區(qū)域隨Ω遞增而加劇,在陣發(fā)性混沌區(qū)域存在突變現(xiàn)象(圖8(b))．動(dòng)載系數(shù)隨齒面沖擊和周期分岔產(chǎn)生了突變,隨齒面沖擊加劇其動(dòng)載系數(shù)俱增(圖8(c)),在單邊沖擊區(qū)其δDLC隨Ω遞增而遞減,在雙邊沖擊區(qū)隨Ω遞增而遞增．即齒面沖擊和周期分岔是影響齒面脫嚙、齒背接觸及動(dòng)載特性的主要因素．

圖8 脫嚙比、動(dòng)載系數(shù)分布圖Fig. 8 The δNMDC, δBMDC and δDLC distribution diagrams

時(shí)變嚙合剛度和頻率比平面內(nèi)的I/P,δDLC,δBMDC和δNMDC解域界結(jié)構(gòu)揭示了系統(tǒng)隨參數(shù)分岔、脫嚙、沖擊、動(dòng)載荷的轉(zhuǎn)遷規(guī)律及其關(guān)聯(lián)關(guān)系．針對(duì)不同速度車型要求,依據(jù)解域界結(jié)構(gòu)對(duì)轉(zhuǎn)速和重合度參數(shù)進(jìn)行匹配優(yōu)化．選取參數(shù)時(shí)盡量避開混沌、雙邊沖擊和動(dòng)載較大的參數(shù)區(qū)域,盡可能在0/1區(qū)域選擇合理的重合度或轉(zhuǎn)速參數(shù)能有效地降低嚙合沖擊、脫嚙和動(dòng)載對(duì)齒輪箱的影響,可延長(zhǎng)齒輪箱的服役周期、降低設(shè)備維護(hù)費(fèi)用．

5 結(jié) 論

本文針對(duì)7自由度直齒錐齒輪傳動(dòng)系統(tǒng)的周期分岔、齒面沖擊、脫嚙和動(dòng)載特性問(wèn)題,采用胞映射理論構(gòu)建了重合度和轉(zhuǎn)速激勵(lì)參數(shù)平面,用CPNF法進(jìn)行了數(shù)值求解獲得了I/P,δNMDC,δBMDC和δDLC解界域結(jié)構(gòu),得到了如下結(jié)論:

1) 參數(shù)平面內(nèi)基于胞映射理論的CPNF法解域界結(jié)構(gòu)求解方法是一種高效的數(shù)值求解方法．

2) 含間隙的錐齒輪系統(tǒng)存在鞍結(jié)、倍化、Hopf、激變、擦切和周期3等分岔現(xiàn)象,隨時(shí)變嚙合剛度系數(shù)增加其齒面沖擊加劇,在嚙頻附近系統(tǒng)的周期運(yùn)動(dòng)和齒面沖擊最為嚴(yán)重．

3) 系統(tǒng)分岔和齒面沖擊轉(zhuǎn)遷導(dǎo)致其齒面脫嚙和齒背接觸及動(dòng)載行為產(chǎn)生了突變,在各I/P區(qū)域內(nèi)其值也隨參數(shù)波動(dòng)．隨時(shí)變嚙合剛度系數(shù)增大其動(dòng)態(tài)特性加劇．在1/1和2/N區(qū)域其脫嚙最為嚴(yán)重,在2/1區(qū)域內(nèi)齒背接觸最為嚴(yán)重,其動(dòng)載系數(shù)最大．