大行程機床動態(tài)特性靈敏度分析及實驗研究

劉海濤，王建華

(西安工業(yè)大學(xué) 機電工程學(xué)院，西安 710021)

?

大行程機床動態(tài)特性靈敏度分析及實驗研究

劉海濤，王建華

(西安工業(yè)大學(xué) 機電工程學(xué)院，西安 710021)

為大行程機床動態(tài)特性隨加工區(qū)間、姿態(tài)的不同將發(fā)生很大的改變，基于單一狀態(tài)下的動力學(xué)模型/方程不能描述其全加工空間內(nèi)動態(tài)特性的變化規(guī)律。文章以動梁式龍門機床為例，采用集總參數(shù)法建立了機床的動力學(xué)模型，并推導(dǎo)了其變系數(shù)線性動力學(xué)方程。采用靈敏度分析法研究了機床第1階固有頻率對于位置參數(shù)的敏感程度及變化規(guī)律。結(jié)果表明：系統(tǒng)第1階特征值對于位置參數(shù)z的靈敏度大于對于位置參數(shù)x的靈敏度，系統(tǒng)第1階固有頻率主要由位置參數(shù)z決定，且變化幅度可達(dá)40%以上。最后通過實驗測試驗證了模型的準(zhǔn)確性。

大行程機床；動態(tài)特性；靈敏度分析

0 引言

機床的動態(tài)特性(模態(tài)、頻響等)對于機床的加工精度和加工零件的表面質(zhì)量有著很大的影響[1-2]。大行程機床(行程大于1500mm)在加工空間內(nèi)其狀態(tài)是在不斷變化的，因此，其各種特性，例如靜態(tài)剛度、固有頻率、振型和頻響等，也是不斷變化的。要想在設(shè)計階段進(jìn)行優(yōu)化設(shè)計，在加工階段進(jìn)行合理的刀具路徑選擇，就必須獲得機床在整個加工空間內(nèi)的動態(tài)特性。

國內(nèi)外許多學(xué)者采用理論建模[3-5]和實驗研究[6]等方法，建立了各自的剛度場模型來描述這一問題，然而當(dāng)機床在加工過程中狀態(tài)發(fā)生改變時，各軸相對位置的變化將改變系統(tǒng)質(zhì)量、阻尼和剛度的分布。因此，高速機床在加工空間內(nèi)狀態(tài)改變的本質(zhì)是，反映系統(tǒng)質(zhì)量、阻尼和剛度大小和分布的質(zhì)量、阻尼和剛度矩陣將不再是常系數(shù)矩陣，所以，要想弄清高速機床動態(tài)特性的變化，就必須研究其質(zhì)量、阻尼和剛度矩陣在不同狀態(tài)下的變化規(guī)律，基于單一狀態(tài)下建立的機床動力學(xué)模型無法描述和分析這一變化規(guī)律，因此有必要針對機床狀態(tài)的變化建立一個動態(tài)模型來反映這一變化。

1 機床動態(tài)動力學(xué)建模

1.1 龍門機床集總參數(shù)模型

本文以典型的動梁式龍門機床為例來研究不同狀態(tài)下機床的變動特性問題，如圖1所示。由于立柱和床身剛性較大，橫梁沿Y軸運動時機床動態(tài)特性變化很小，因此，本文忽略機床沿Y方向的振動，只考慮橫梁-滑板-主軸箱-主軸系統(tǒng)在XOZ平面內(nèi)的平動及轉(zhuǎn)動。采用一致質(zhì)量梁單元來建立橫梁和主軸箱結(jié)構(gòu)，采用集中質(zhì)量單元來等效滑板和主軸單元，采彈簧-阻尼單元來建立各進(jìn)給軸導(dǎo)軌-滑塊結(jié)合部、絲杠傳動系統(tǒng)模型，如圖2所示。

圖1 動梁式龍門機床結(jié)構(gòu)示意圖

圖2 龍門機床集總參數(shù)模型

圖2所示的模型中，節(jié)點1和節(jié)點2之間的距離x以及節(jié)點5和節(jié)點6之間的距離z為變量，由于機床X軸的行程0.85m，Z軸行程為0.55m，所以x∈[0.735,1.585]，z∈[0.505,1.055]。這樣，系統(tǒng)在不同位置下反映其質(zhì)量、剛度大小及分布的質(zhì)量矩陣[M]、剛度矩陣[K]形式一致且可以表示為位置變量x、z的函數(shù)，從而有利于采用靈敏度分析的方法進(jìn)行分析。

1.2 單元質(zhì)量、剛度、阻尼矩陣的構(gòu)造

采用一致質(zhì)量梁單元來描述模型中主要結(jié)構(gòu)件橫梁和主軸箱的特性，即單元1、2、3、4。當(dāng)主軸箱沿X方向運動到位置x，沿Z方向運動到位置z時，單元1和單元2的質(zhì)量矩陣為：

式中，ρb為橫梁等效密度；Ab為橫梁截面面積；L為橫梁的長度。

單元1和單元2的剛度矩陣為：

式中，Izb為橫梁等效截面慣性矩；E為橫梁材料的彈性模量(鑄鐵，120GPa)。

同理可以得到單元3和4的質(zhì)量矩陣和剛度矩陣。

結(jié)合面單元的剛度矩陣為：

式中，kx、kz、kθ分別為結(jié)合面x方向剛度、z方向剛度以及繞Y軸的扭轉(zhuǎn)剛度。本模型只考慮結(jié)合面的阻尼而忽略材料自身的阻尼，其阻尼矩陣為：

式中，cx、cz和cθ分別為結(jié)合面x方向阻尼、z方向阻尼以及繞Y軸的扭轉(zhuǎn)阻尼。

滑板和主軸被簡化為集中質(zhì)量單元，其質(zhì)量矩陣為：

將各單元矩陣組裝起來即可得到整個機床系統(tǒng)動力學(xué)方程矩陣。

2 機床系統(tǒng)動態(tài)特性靈敏度分析

2.1 特征值一階靈敏度分析

對于第1節(jié)所描述的動力學(xué)系統(tǒng)，系統(tǒng)的對于參數(shù)u的一階特征值靈敏度[7-8]為：

其中,λr和ψr分別為系統(tǒng)第r階特征值和特征向量，ar為常數(shù)。因此，系統(tǒng)對于位置參數(shù)x和位置參數(shù)z的靈敏度分別為:

圖3為位置參數(shù)x和z變化時，機床系統(tǒng)第1階固有頻率的靈敏度值?？梢钥闯觯到y(tǒng)特征值對于位置參數(shù)z的靈敏度S(λ/z)都大于特征值對于位置參數(shù)x的靈敏度S(λ/x)，因此，系統(tǒng)的第1階特征值，即固有頻率主要由位置參數(shù)z決定，但位置參數(shù)x也會對靈敏度S(λ/z)產(chǎn)生影響。

(a)S(λ1/x)

(b)S(λ1/z)

2.2 不同加工區(qū)間固有頻率(特征值)變化

基于本文中的動態(tài)模型，系統(tǒng)的質(zhì)量、剛度矩陣為位置參數(shù)x和z的函數(shù)，前后兩個狀態(tài)下系統(tǒng)質(zhì)量和剛度矩陣的變化量可以表示為位置參數(shù)變化量Δx和Δz的函數(shù)[M]=[M0]+[ΔM(Δx,Δz)]和[K]=[K0]+[ΔK(Δx,Δz)]，其中Δx=xi+1-xi，Δz=zi+1-zi。因此，固定位置參數(shù)變化量Δx和Δz，采用一階矩陣攝動理論可以很方便的來求解系統(tǒng)在全加工空間內(nèi)的特征值[9-11]。

式中，λ0i為變化前系原始系統(tǒng)第i階特征值，λ1i為特征值一階攝動量，Ψ0i為原始系統(tǒng)特征值矩陣第i列所組成的向量。

由圖4可以看出，系統(tǒng)的第1階固有頻率主要隨位置參數(shù)z變化而變化，這一點與靈敏度得到的結(jié)論相同，且最大值與最小值之間變化的幅度可達(dá)40%以上，變化相當(dāng)大，因此，對于大行程機床而言，這一變化在設(shè)計階段和實際的加工過程中不能忽略。

圖4 系統(tǒng)第1階固有頻率隨位置參數(shù)x和z的變化

3 驗證實驗

為了驗證計模型的準(zhǔn)確性，采用16通道LMS模態(tài)測試系統(tǒng)，在動梁式龍門機床上布置6個三向加速度傳感器，采用錘擊法對機床幾個極端位置的頻率相應(yīng)進(jìn)行了測試，如圖5所示。

在LMS.Test.Lab中對模態(tài)測試得到的頻響函數(shù)進(jìn)行處理，獲得機床的1階固有頻率和振型[12]。計算的結(jié)果與實驗測試的結(jié)果比較如圖6所示?？梢钥闯觯m然由于模型的簡化、結(jié)合面參數(shù)不精確等因素，計算值與測試值有一定的誤差，最大誤差為9%左右，但是總體趨勢符合的較好。

圖5 錘擊法對機床幾個極端位置進(jìn)行模態(tài)測試

圖6 4個極端位置機床1階固有頻率值

4 結(jié)論

本文采用集總參數(shù)法建立了大行程機床的動態(tài)動力學(xué)模型，推導(dǎo)了其變系數(shù)線性動力學(xué)方程，采用靈敏度分析法來描述系統(tǒng)特征值隨位置參數(shù)的變化規(guī)律，并進(jìn)行了實驗驗證，得到的主要結(jié)論如下：

(1)系統(tǒng)第1階特征值對于位置參數(shù)z的靈敏度S(λ/z)大于特征值對于位置參數(shù)x的靈敏度S(λ/x)；

(2)系統(tǒng)第1階固有頻率主要由位置參數(shù)z決定，這與靈敏度分析的結(jié)果相同，且對于本機床而言，變化的幅度可達(dá)40%以上。

[1] Y Altintas, M Weck. Chatter Stability of Metal Cutting and Grinding [J]. CIRP Annals Manufacturing Technology, 2004 (53):619-642.

[2] Y Altintas,G Stepan,D Merdol. Chatter stability of milling in frequency and discrete time domain[J]. CIRP Journal of Manufacturing Science and Technology, 2008(1):35-44.

[3] Jun Wu, Jinsong Wang, Liping Wang,et al. Study on the stiffness of a 5-DOF hybrid machine tool with actuation redundancy [J]. Mechanism and Machine Theory,2009(44):289-305.

[4] Sanger DJ, Chen JQ, Zhang SJ, et al. A general method for the stiffness analysis of manipulator mechanisms [J]. Journal of Mechanical Engineering Science, 2000, 214(5): 673-685 .

[5] M Law, S Ihlenfeldt, M Wabner,et al. Position-dependent dynamics and stability of serial-parallel kinematic machines [J]. CIRP Annals-Manufacturing Technology, article in press,2013(62)：375-378.

[6] B Paijmans, W Symens, H V Brussel, et al. Experimental identification of affine LPV models for mechatronic systems with one varying parameter [J]. European Journal of Control, 2008, 14 (1): 16-29.

[7] 李德葆．結(jié)構(gòu)修改—靈敏度分析和修改結(jié)構(gòu)重分析[J]．機械強度，1990，12(2) ：44-50．

[8] 王書娟．結(jié)構(gòu)參數(shù)小幅變化后橋梁固有模態(tài)修正的矩陣攝動法[D]．長春:吉林大學(xué),2005．

[9] 杜思義．基于頻率變化與單元矩陣攝動理論的結(jié)構(gòu)損傷識別方法研究[D]．重慶:重慶大學(xué)，2005．

[10] 陳塑寰．結(jié)構(gòu)振動分析的矩陣攝動理論[M]．重慶:重慶出版社，1991．

[11] 陳塑寰．結(jié)構(gòu)動態(tài)設(shè)計的矩陣攝動理論[M]．北京:科學(xué)出版社，2011．

[12] 伍偉，黃美發(fā)，韋禮凱,等．?dāng)?shù)控龍門銑床整機動態(tài)特性分析[J] ．組合機床與自動化加工技術(shù)，2015(11):19-21.

(編輯 李秀敏)

Dynamic Characteristics Sensitivity Analysis and Experimental Research of Large Stroke Machine Tools

LIU Hai-tao，WANG Jian-hua

(School of Mechatronic Engineering, Xi'an Technological University, Xi'an 710021, China)

The dynamic characteristics of large stroke machine tools will highly change when its machining space and cutting pose vary. So that the dynamic model/equation based on single state cannot describe the variation of dynamic characteristics in its whole machining space. This paper take a gantry beam-moving machine tool as an example to demonstrate the method of lumped parameter dynamic model of machine tool, and the coefficient-varying linear dynamic equation is derived. The sensitivity analyzing method is used to studied the sensitivity and variation rules of the first natural frequencies with respect to position parameters. The results show that the sensitivities of the first natural frequencies with respect to position parameter z is greater than that with respect to position parameter x. The first natural frequencies are dominated by position parameter z, and the variation magnitude can be up to above 40%. At last the model accuracy is validated by experimental measurement.

large stroke machine tool; dynamic characteristics; sensitivity analysis

1001-2265(2016)11-0035-03

10.13462/j.cnki.mmtamt.2016.11.010

2015-12-08;

2016-01-19

劉海濤(1981—)，男，山東青島人，西安工業(yè)大學(xué)講師，博士，研究方向為機床性能仿真、精密測量、產(chǎn)品流程自動化，(E-mail)429059044@qq.com。

TH166;TG659

A