李鵬松, 盛桂全, 孟永永

(東北電力大學(xué) 理學(xué)院, 吉林 吉林 132012)

(2iωc+2ζ+ωcτe-iωcτ)A′-2iδζA+HA2=0,(17)

再生型顫振系統(tǒng)的Hopf分岔分析與控制

李鵬松, 盛桂全, 孟永永

(東北電力大學(xué) 理學(xué)院, 吉林 吉林 132012)

考慮一類單自由度的非線性再生型切削顫振系統(tǒng), 利用多尺度法分析系統(tǒng)時滯參數(shù)對解穩(wěn)定性的影響, 并在此基礎(chǔ)上采用非線性位移反饋控制抑制Hopf分岔引起的顫振.理論分析和數(shù)值模擬結(jié)果驗證了所給控制方法在切削顫振控制模型中的有效性.

亞臨界Hopf分岔； 超臨界Hopf分岔； 非線性控制； 多尺度法

圖1 刀具切削過程中的再生效應(yīng)Fig.1 Regeneration effect of cutting process

機床、刀具和工件形成一個復(fù)雜的動力學(xué)結(jié)構(gòu), 在某種條件下, 系統(tǒng)會發(fā)生劇烈的振動.可將這種振動分成三類： 自由振動、強迫振動、自激振動.在機床動力學(xué)領(lǐng)域自激振動也稱為“顫振”.常見的顫振主要有由速度反饋引起的切削顫振、模態(tài)藕合引起的切削顫振以及再生型機理因素引起的顫振.在工件材料切削過程中發(fā)生顫振, 常會使加工過程變得不穩(wěn)定, 導(dǎo)致加工表面質(zhì)量和金屬切削率下降, 嚴(yán)重時甚至?xí)茐牡毒吆蜋C床.研究表明, 切削過程中的再生效應(yīng)是引起顫振的主要原因[1-4].一些外部擾動或工件表面的一個斑點即可引起切削力的微小變化, 從而導(dǎo)致動力系統(tǒng)的振動, 振動可在工件表面留下波紋狀的切削痕跡， 如圖1所示.波紋狀的表面將影響后續(xù)的切削, 當(dāng)?shù)毒邚墓ぜ砻媲谐∑瑫r, 振動導(dǎo)致薄片不平.如果振動振幅不斷衰退, 切削過程即為穩(wěn)定的.但在某種情況下, 振幅會變大, 同時產(chǎn)生顫振, 這種現(xiàn)象稱為再生型顫振.發(fā)生再生型顫振的系統(tǒng), 通常因時滯參數(shù)變化導(dǎo)致平衡點失穩(wěn)而引起Hopf分岔.因此, 研究再生型顫振系統(tǒng)離不開對Hopf分岔的研究[5-10].本文主要分析再生顫振系統(tǒng)在不同主軸旋轉(zhuǎn)周期情況下的運動穩(wěn)定性, 對系統(tǒng)失穩(wěn)的過程進行分岔分析及類型判定, 對發(fā)生Hopf分岔的系統(tǒng)設(shè)計合理的控制器, 從而將系統(tǒng)控制到穩(wěn)定, 并通過實際算例和仿真分析驗證所設(shè)計控制器的實用性和有效性.

1 再生型顫振系統(tǒng)

Hanna等[11]針對刀具切削系統(tǒng), 提出了非線性再生型顫振理論, 并建立了單自由度無量綱形式的振動方程

其中:xτ=x(t-τ),x(t)表示切削刀具尖端動態(tài)位移,τ表示主軸的旋轉(zhuǎn)周期;ζ是阻尼比;w是切割寬度;α2和α3表示非線性切削力的系數(shù).系統(tǒng)(1)的線性形式為

若式(2)的解形式設(shè)為x(t)=x0eλt, 則式(2)的特征方程為

在復(fù)數(shù)域內(nèi), 令λ=σ+iω, 代入式(3)并分離實部、虛部得

1+w+2ζσ+σ2-ω2-we-στcos(ωτ)=0, 2(ζ+σ)ω+we-στsin(ωτ)=0.

當(dāng)σ>0時, 隨t的增大x=x0e(σ±iω)t→∞, 系統(tǒng)(1)是不穩(wěn)定的； 當(dāng)σ<0時, 隨t的增大x=x0e(σ±iω)t→0, 系統(tǒng)(1)是穩(wěn)定的； 當(dāng)σ=0時,x=x0e(σ±iω)t被定義為系統(tǒng)(1)的穩(wěn)定臨界點, 即可能發(fā)生Hopf分岔的位置.為了找到該臨界點, 可令式(3)中的σ=0, 得到復(fù)特征方程為

將式(4)分離實部與虛部得

設(shè)臨界點為(σ,τ,ω)=(0,τc,ωc), 此時系統(tǒng)(1)發(fā)生Hopf分岔, 因為式(2)的兩個復(fù)共軛特征根λ=σ±iω穿過虛軸.為了說明這種情況, 對式(3)做全微分, 可得λ關(guān)于τ的導(dǎo)數(shù)形式如下:

進而有

此時系統(tǒng)(1)發(fā)生Hopf分岔.由式(5)可求得

其中:w和ζ均為小于1的正的常數(shù);k=1,2,3,….

2 再生型顫振系統(tǒng)的Hopf分岔類型判定

引入新變量Tn=εnt(n=0,1,…), 即T0=t,T1=εt,T2=ε2t, …,ε表示一個無量綱的小參數(shù).其關(guān)于t的導(dǎo)數(shù)定義為

（9）

如果系統(tǒng)(1)的解含有慢尺度變量T1=εt, 則將在O(ε3)處出現(xiàn)永年項, 因此其解不應(yīng)含有慢尺度變量, 所以其解可近似表示為

令

其中δ為調(diào)諧參數(shù).將式(10),(11)代入式(1), 并比較方程兩邊ε,ε2和ε3的系數(shù), 可得下列攝動方程:

其中xi=xi(T0,T2);xiτ=xi(T0-τ,T2).式(12)的解可寫為

將式(15)代入式(13), 解得

將式(15),(16)代入式(14), 為消除永年項, 令e±iωcτ的系數(shù)為零, 可得

在復(fù)數(shù)域內(nèi), 令A(yù)(T2)=u(T2)+iv(T2), 則式(17)化為

分離式(18)的實部、虛部得

其中h(u,v)和g(u,v)為不低于2次的非線性項, 表示如下:

求得特征值為λ1,2=-βe±iυe.

當(dāng)βe=0,υe≠0時, 式(22)存在一對純虛根λ1,2=±iυe, 根據(jù)Hopf分岔定理, 此時系統(tǒng)發(fā)生Hopf分岔.由于βe和υe是時滯參數(shù)τ的函數(shù), 因此調(diào)整時滯參數(shù)τ可控制Hopf分岔發(fā)生, 也將影響Hopf分岔方向和周期解的穩(wěn)定性.應(yīng)用Poincare-Birkhoff范式定理對系統(tǒng)(1)發(fā)生的Hopf分岔類型進行判斷, 令

將式(20),(21)代入式(23)得

當(dāng)βe=0,υe≠0時,δ≠0,Q≠0, 所以根據(jù)βe=-2Qδζq2,υe=2Qδζq1有q2=0,q1≠0, 則式(24)變?yōu)?/p>

當(dāng)η<0時, 系統(tǒng)(1)發(fā)生超臨界Hopf分岔, 產(chǎn)生穩(wěn)定的極限環(huán)； 當(dāng)η>0時, 系統(tǒng)(1)發(fā)生亞臨界Hopf分岔, 產(chǎn)生不穩(wěn)定的極限環(huán).

3 數(shù)值仿真

針對系統(tǒng)(1), 取參數(shù)ζ=0.003 2,w=0.012,α2=1,α3=1,δ=1,ωc=2.011 6.當(dāng)參數(shù)τ∈(0,3.908 8)時, 系統(tǒng)(1)穩(wěn)定, 時間歷程與相軌跡如圖2所示, 此時τ=3.780 5; 當(dāng)參數(shù)τ∈(3.908 9,3.910 5)時, 發(fā)生超臨界Hopf分岔, 時間歷程與相軌跡如圖3所示, 此時τ=3.908 9,η=-0.036 6; 當(dāng)τ∈(3.910 6,5.468 4)時, 系統(tǒng)發(fā)生亞臨界Hopf分岔, 時間歷程與相軌跡如圖4所示, 此時τ=4.778 7,η=1.117 3; 當(dāng)τ∈(5.468 5,5.469 4)時, 系統(tǒng)再次發(fā)生超臨界Hopf分岔(圖略); 當(dāng)參數(shù)τ>5.469 4時系統(tǒng)穩(wěn)定.

圖2 當(dāng)τ=3.780 5時穩(wěn)定系統(tǒng)的運動情況Fig.2 Motion state of stable system for τ=3.780 5

圖3 當(dāng)τ=3.908 9時發(fā)生超臨界Hopf分岔系統(tǒng)的運動情況Fig.3 Movement of system with supercritical Hopf bifurcation for τ=3.908 9

圖4 當(dāng)τ=4.778 7時發(fā)生亞臨界Hopf分岔系統(tǒng)的運動情況Fig.4 Movement of system with subcritical Hopf bifurcation for τ=4.778 7

隨著參數(shù)的變化, 系統(tǒng)(1)會先由穩(wěn)定到發(fā)生超臨界Hopf分岔, 再發(fā)生亞臨界Hopf分岔, 然后再在某個參數(shù)的鄰域內(nèi)發(fā)生超臨界Hopf分岔, 最后當(dāng)參數(shù)大于該參數(shù)的鄰域時系統(tǒng)漸近穩(wěn)定.實際上, 如果系統(tǒng)(1)發(fā)生Hopf分岔, 則會使系統(tǒng)失去穩(wěn)定性, 即刀具器械在實際工作過程中將遭受嚴(yán)重損壞.下面主要針對τ∈(0,10)區(qū)域, 為系統(tǒng)(1)設(shè)計合理、有效的控制器, 避免系統(tǒng)發(fā)生Hopf分岔.

4 再生型顫振系統(tǒng)的分岔控制

圖5 刀具振動模型Fig.5 Tool vibration models

刀具及其刀具固定器可近似簡化為一個單自由度的彈簧系統(tǒng), 其阻尼系數(shù)為C, 彈簧剛性系數(shù)為k, 刀具質(zhì)量為m.實際的工件切削厚度可由下式給出：

其中:μ表示切削寬度;M的大小與其切削參數(shù)有關(guān), 如切削速度、刀具的幾何形狀和f0大小等.假設(shè)M和μ均為常數(shù), 對式(28)關(guān)于f在f0處進行三階Taylor展開并代入方程(27), 對其進行無量綱化得到具有位移反饋系統(tǒng)的運動方程為

其中:xτ0=x(t-τ0),τ0<τ;α2,α3,a和k0為常數(shù).

取參數(shù)ζ=0.003 2,w=0.012,α2=1,α3=1,k0=2.540 3×10-5,τ0=3.0,a=9.當(dāng)參數(shù)τ∈(0,3.908 8)和τ>5.469 4時, 控制系統(tǒng)(31)穩(wěn)定; 當(dāng)參數(shù)τ∈(3.908 9,3.910 5)時, 原系統(tǒng)發(fā)生超臨界Hopf分岔, 而通過數(shù)值仿真可驗證控制系統(tǒng)(31)同樣穩(wěn)定, 原系統(tǒng)及控制系統(tǒng)的時間歷程如圖6所示, 此時τ=3.908 9; 當(dāng)τ∈(3.910 6,5.468 4)時, 原系統(tǒng)發(fā)生亞臨界Hopf分岔, 而通過數(shù)值仿真可驗證控制系統(tǒng)(31)仍然穩(wěn)定, 原系統(tǒng)及控制系統(tǒng)的時間歷程如圖7所示, 此時τ=4.778 7; 當(dāng)τ∈(5.468 5,5.469 4)時, 原系統(tǒng)發(fā)生超臨界Hopf分岔, 而通過數(shù)值仿真可驗證控制系統(tǒng)(31)仍然穩(wěn)定(圖略).

圖6 τ=3.908 9時臨界Hopf分岔的控制Fig.6 Supercritical Hopf bifurcation control for τ=3.908 9

圖7 τ=4.778 7時亞臨界Hopf分岔的控制Fig.7 Subcritical Hopf bifurcation control for τ=4.778 7

綜上所述, 本文針對再生切削顫振系統(tǒng), 利用多尺度法結(jié)合Poincare-Birkhoff范式定理, 分析了系統(tǒng)的分岔特性, 并針對系統(tǒng)發(fā)生的Hopf分岔, 設(shè)計了非線性延遲控制器, 將原系統(tǒng)發(fā)生的亞臨界Hopf分岔、超臨界Hopf分岔均控制為穩(wěn)定, 并通過數(shù)值模擬驗證了所設(shè)計控制器的有效性.

[1]趙愿玲, 左曙光.考慮驅(qū)動力影響的輪胎側(cè)向自激振動分析 [J].振動與沖擊, 2012, 31(22): 101-105.(ZHAO Yuanling, ZUO Shuguang.Lateral Self-excited Vibration Analysis for a Type Considering Driving Force [J].Journal of Vibration and Shock, 2012, 31(22): 101-105.)

[2]賈尚帥, 丁千, 劉煒.超音速彈翼非線性顫振分析與控制 [J].振動與沖擊, 2012, 31(13): 108-112.(JIA Shangshuai, DING Qian, LIU Wei.Nonlinear Flutter Analysis and Control of Supersonic Missile Wings [J].Journal of Vibration and Shock, 2012, 31(13): 108-112.)

[3]陽輝, 黃筱調(diào), 方成剛.數(shù)控再生型切削顫振系統(tǒng)的分析及仿真 [J].機床與液壓, 2009, 37(1): 35-37.(YANG Hui, HUANG Xiaodiao, FANG Chenggang.The Analysis and Simulation of NC Regenerative Cutting Chatter System [J].Machine Tool and Hydraulics, 2009, 37(1): 35-37.)

[4]劉習(xí)軍, 王立剛, 賈啟芬.一種由干摩擦引起的車床切削顫振 [J].工程力學(xué), 2005, 22(1): 107-112.(LIU Xijun, WANG Ligang, JIA Qifen.Cutting Chatter of Machine Tools Induced by Dry Friction [J].Engineering Mechanics, 2005, 22(1): 107-112.)

[5]Pilkee Kim, Sanghyun Bae, Jongwon Seok.Bifucation Analysis on a Turning System with Large and State-Dependent Time Delay [J].Journal of Sound and Vibration, 2012, 331: 5562-5580.

[6]李鵬松, 陳書吉, 呂雪, 等.單參數(shù)電力系統(tǒng)亞臨界Hopf分岔控制 [J].吉林大學(xué)學(xué)報: 理學(xué)版, 2013, 51(4): 618-622.(LI Pengsong, CHEN Shuji, Lü Xue, et al.Subcritical Hopf Bifurcation Control of Power System with Single Parameter [J].Journal of Jilin University: Science Edition, 2013, 51(4): 618-622.)

[7]劉爽, 劉彬, 張業(yè)寬, 等.一類時滯非線性相對轉(zhuǎn)動系統(tǒng)的Hopf分岔與周期解的穩(wěn)定性 [J].物理學(xué)報, 2010, 59(1): 38-43.(LIU Shuang, LIU Bin, ZHANG Yekuan, et al.Hopf Bifurcation and Stability of Periodic Solutions in a Nonlinear Relative Rotation Dynamical System with Time Delay [J].Acta Physica Sinica, 2010, 59(1): 38-43.)

[8]陳衍茂, 劉濟科.非線性顫振系統(tǒng)中既是超臨界又是亞臨界的Hopf分岔點研究 [J].應(yīng)用數(shù)學(xué)和力學(xué), 2008, 29(2): 181-187.(CHEN Yanmao, LIU Jike.Supercritical as Well as Subcritical Hopf Bifurcation in Nonlinear Flutter Systems [J].Applied Mathematics and Mechanics, 2008, 29(2): 181-187.)

[9]XIE Qizhi, ZHANG Qichang, HAN Jianxin.Hopf Bifurcation for Delay Differential Equation with Application to Machine Tool Chatter [J].Applied Mathematical Modelling, 2012, 36(8): 3803-3812.

[10]Faria T.Stability and Bifurcation for a Delayed Predator-Prey Model and the Effect of Diffusion [J].J Math Anal Appl, 2001, 254(2): 433-463.

[11]Hanna N H, Tobias S A.A Theory of Nonlinear Regenerative Chatter [J].J Manuf Sci Eng, 1974, 96(1): 247-255.

[12]陸啟韶, 彭臨平, 楊卓琴.常微分方程與動力系統(tǒng) [M].北京: 北京航空航天大學(xué)出版社, 2010.(LU Qishao, PENG Linping, YANG Zhuoqin.Ordinary Differential Equations and Dynamical Systems [M].Beijing: Beihang University Press, 2010.)

AnalysisandControlofHopfBifurcationforRegenerativeChatterSystem

LI Pengsong, SHENG Guiquan, MENG Yongyong
(CollegeofSciences,NortheastDianliUniversity,Jilin132012,JilinProvince,China)

We considered a class of retarded one-degree-of-freedom nonlinear system, namely, the regenerative model of chatter in a lathe machine tool.The time-delay influence on the stability of the system solutions was analyzed with the aid of multiple scales.Nonlinear displacement feedback control can be used to suppress the self-sustained oscillation caused by Hopf bifurcation.The results of the theory analysis and numerical simulation show that the designed controller is valid for the cutting chatter system.

subcritical Hopf bifurcation; supercritical Hopf bifurcation; nonlinear control; method of multiple scales

2014-03-03.

李鵬松(1970—), 男, 漢族, 博士, 教授, 從事非線性動力系統(tǒng)的研究, E-mail： li-pengsong@163.com.

國家自然科學(xué)基金(批準(zhǔn)號： 11072085)和吉林省科技發(fā)展計劃項目(批準(zhǔn)號： 20130101065JC).

O193

A

1671-5489(2014)06-1155-07

10.13413/j.cnki.jdxblxb.2014.06.09

趙立芹)