袁 黎，胡忠舉，梁紅強(qiáng)，彭 威，朱萍玉

(1.湖南科技大學(xué) 機(jī)電工程學(xué)院，湖南 湘潭 411201； 2.廣州大學(xué) 機(jī)械與電氣工程學(xué)院，廣東 廣州 510006)

基于工作變形分析的鏜桿設(shè)計(jì)研究*

袁 黎1，胡忠舉1，梁紅強(qiáng)1，彭 威1，朱萍玉2

(1.湖南科技大學(xué) 機(jī)電工程學(xué)院，湖南 湘潭 411201； 2.廣州大學(xué) 機(jī)械與電氣工程學(xué)院，廣東 廣州 510006)

基于振動(dòng)方程，采用傳遞矩陣法和模態(tài)疊加法，推導(dǎo)了鏜桿的工作變形公式；應(yīng)用算例對(duì)該工作變形公式進(jìn)行了驗(yàn)證。選取普通鏜桿和阻尼鏜桿作為研究對(duì)象，采用該工作變形公式快速分析兩種鏜桿的工作振型、模態(tài)貢獻(xiàn)量和頻響曲線，進(jìn)行對(duì)比分析，指導(dǎo)鏜桿的結(jié)構(gòu)設(shè)計(jì)。結(jié)果表明該工作變形分析作為快速分析鏜桿工作振動(dòng)特性的新方法，為鏜桿結(jié)構(gòu)參數(shù)設(shè)計(jì)計(jì)算提供了新的參考依據(jù)。

工作變形；鏜桿；模態(tài)貢獻(xiàn)量；頻響曲線

0 引 言

工作變形(Operational Deflection Shape， ODS)分析是一種新的結(jié)構(gòu)振動(dòng)分析方法，它區(qū)別于模態(tài)振型分析。模態(tài)是結(jié)構(gòu)的固有屬性，而ODS指結(jié)構(gòu)在某一工況或特定頻率下的受迫響應(yīng)，隨結(jié)構(gòu)載荷的變化而變化[1-4]。工程實(shí)際中，ODS可用于對(duì)某階模態(tài)振型作用溯源，識(shí)別在工作變形中的貢獻(xiàn)大小，進(jìn)而為結(jié)構(gòu)設(shè)計(jì)做指導(dǎo)。筆者針對(duì)有限元法仿真分析計(jì)算量較大，分析時(shí)間長，對(duì)計(jì)算機(jī)配置要求較高的問題，擬采用一種新方法對(duì)設(shè)計(jì)前后的鏜桿進(jìn)行振動(dòng)特性研究，可快速得到鏜桿的工作振型、模態(tài)貢獻(xiàn)量和頻響曲線等振動(dòng)特性參數(shù)，通過這些參數(shù)對(duì)比，即可指導(dǎo)鏜桿的結(jié)構(gòu)設(shè)計(jì)。

1 鏜桿工作變形函數(shù)推導(dǎo)

鏜桿可視為有限自由度系統(tǒng)。將其離散成N個(gè)質(zhì)點(diǎn)系，質(zhì)點(diǎn)系間通過剛架結(jié)構(gòu)元連接。根據(jù)傳遞矩陣法，質(zhì)點(diǎn)矢量位移在這些質(zhì)點(diǎn)之間傳遞[5-7]。連續(xù)鏜桿系統(tǒng)便轉(zhuǎn)化為由多個(gè)連接元組成的N自由度系統(tǒng)，如圖1所示。

圖1 鏜桿等效形式

忽略阻尼的影響,此N自由度系統(tǒng)振動(dòng)方程為：

(1)

質(zhì)量陣和剛度陣為：

(2)

式中：mij為質(zhì)量陣中的元素；kij為剛度陣中的元素。

其模態(tài)頻率和模態(tài)振型通過式(3)計(jì)算:

(3)

設(shè)鏜桿可視為線性時(shí)不變系統(tǒng)，則離散鏜桿系統(tǒng)任一點(diǎn)l處某時(shí)刻的響應(yīng)可表示為該剛架結(jié)構(gòu)元各階模態(tài)的線性疊加[8]：

(4)

式中：l為被測(cè)點(diǎn)；r為模態(tài)階次；φlr為被測(cè)點(diǎn)l處的振型系數(shù)；qr(ω)為比例因子；Mr為模態(tài)質(zhì)量；Kr為模態(tài)剛度。

綜合式(2)～(4)，得鏜桿的工作變形函數(shù)：

(5)

2 算 例

設(shè)簡易鏜桿的橫截面形狀及尺寸如圖2所示。彈性模量為200 GPa，泊松比為0.3，密度為7 850 kg/m3，懸伸長度為225 mm。將鏜桿離散成9段，如圖3所示。

圖2 鏜桿橫截面 圖3 鏜桿離散示意圖

通過有限元仿真分析，可得到鏜桿在任意頻率下的工作振型[9]，分別取頻率為500 Hz、1 000 Hz、2 000 Hz、3 000 Hz時(shí)鏜桿的工作振型，如圖4所示。通過鏜桿的工作變形函數(shù)得到對(duì)應(yīng)的4個(gè)頻率下的工作振型如圖5所示。

圖4 有限元仿真結(jié)果

圖5 工作變形公式計(jì)算結(jié)果

由圖4、5的對(duì)比驗(yàn)證可知，通過工作變形函數(shù)得到的工作振型與仿真分析結(jié)果基本一致，即采用工作變形函數(shù)可正確地描繪出鏜桿的工作振型，且計(jì)算量小，對(duì)計(jì)算機(jī)配置要求很低，使用方便快捷。

3 應(yīng)用研究

選取普通鏜桿和阻尼鏜桿作為研究對(duì)象，如圖6所示，兩種鏜桿的長徑比相同，刀桿部分的形狀和長度也相同，阻尼鏜桿則是對(duì)普通鏜桿的頭部結(jié)構(gòu)形狀進(jìn)行了優(yōu)化。

圖6 研究對(duì)象

采用工作變形函數(shù)可對(duì)這兩種鏜桿進(jìn)行工作模態(tài)分析，任意選取三個(gè)工作頻率：800 Hz、1 700 Hz、2 600 Hz，即可得到該阻尼鏜桿在這三個(gè)工作頻率下的工作振型，如圖7所示。

圖7 工作振型對(duì)比

實(shí)線是普通鏜桿的工作振型，虛線是阻尼鏜桿的工作振型。由圖7可知，這兩種鏜桿在任意頻率下的工作振型相似，但是普通鏜桿的刀頭部分的振動(dòng)幅度明顯較大，即阻尼鏜桿的抗振能力更強(qiáng)。

通過模態(tài)判定準(zhǔn)則可分析各階模態(tài)對(duì)工作變形的貢獻(xiàn)量[10]。假定{x}和{y}為兩個(gè)長度相等的向量，模態(tài)判定準(zhǔn)則MACxy為：

(6)

經(jīng)計(jì)算，兩種鏜桿各階模態(tài)在800 Hz、1 700 Hz和2 600 Hz處ODS的貢獻(xiàn)量如表1、2所列。

表1 普通鏜桿模態(tài)貢獻(xiàn)量

表2 阻尼鏜桿模態(tài)貢獻(xiàn)量

由表1、2可知，工作頻率在800 Hz時(shí)，第1階固有模態(tài)對(duì)兩種鏜桿的振動(dòng)貢獻(xiàn)均較大；工作頻率在1 700 Hz時(shí)，第2階固有模態(tài)對(duì)兩種鏜桿的振動(dòng)貢獻(xiàn)均較大；工作頻率在2 600 Hz時(shí)，第1、2、3階固有模態(tài)對(duì)普通鏜桿的振動(dòng)貢獻(xiàn)較大，第1、2階固有模態(tài)對(duì)阻尼鏜桿的振動(dòng)貢獻(xiàn)較大。通過識(shí)別鏜桿固有模態(tài)在工作變形中的貢獻(xiàn)量，得到極易引起鏜桿顫振的模態(tài)固有頻率。

采用工作變形函數(shù)獲得兩種鏜桿刀頭端的頻率響應(yīng)曲線進(jìn)行對(duì)比，如圖8所示，0～2100 Hz范圍內(nèi)，阻尼鏜桿與普通鏜桿的頻率響應(yīng)曲線大部分吻合，即在該頻率范圍內(nèi)，阻尼鏜桿未體現(xiàn)出較好的動(dòng)態(tài)減振性能；但在2 100～4 700 Hz范圍內(nèi)，阻尼鏜桿的頻率響應(yīng)曲線明顯高于普通鏜桿的，即在該頻率范圍內(nèi)，阻尼鏜桿具備較好的動(dòng)態(tài)減振性能。

圖8 頻響曲線對(duì)比

4 結(jié) 語

推導(dǎo)了一種新的工作變形函數(shù)，并采用算例對(duì)比驗(yàn)證了該工作變形函數(shù)的正確性；以普通鏜桿和阻尼鏜桿作為研究對(duì)象，采用該工作變形函數(shù)對(duì)比研究了兩種鏜桿的工作振型、模態(tài)貢獻(xiàn)量和頻率響應(yīng)曲線等動(dòng)態(tài)特性，得到了阻尼鏜桿適用的頻率范圍，驗(yàn)證了阻尼鏜桿的減振性能。

[1] Devriendt C, Steenackers G, DE S G, et al. From Operating Deflection Shapes Towards Mode Shapes Using Transmissibility Measurements[J]. Mechanical Systems and Signal Processing, 2010,24(3):665-677.

[2] Wonki B, Yongsoo K, Jedol D, et al. Scaling the Operating Deflection Shapes Obtained from Scanning Laser Doppler Vibrometer[J]. Journal of Nondestructive Evaluation,2011,30(2):91-98.

[3] 陳永會(huì)，李志譚，李海虹.基于ODS測(cè)試方法的鉆架模型試驗(yàn)?zāi)B(tài)分析研究[J].太原科技大學(xué)學(xué)報(bào)，2005，26(4)：268-272.

[4] 程 珩，趙 遠(yuǎn).工作模態(tài)分析在減速器監(jiān)測(cè)和診斷中的應(yīng)用 [J]. 振動(dòng)、測(cè)試與診斷，2010，30(2)：197-200.

[5] 蘭賢輝，李育錫，蘇武會(huì).輪轉(zhuǎn)子系統(tǒng)穩(wěn)態(tài)不平衡響應(yīng)分析的整體傳遞矩陣法[J]. 機(jī)械科學(xué)與技術(shù)，2009，28(12)：1613-1615，1620.

[6] 彭澤軍，楊兆建.多支承轉(zhuǎn)子軸承系統(tǒng)的負(fù)荷敏感度研究[J].工程設(shè)計(jì)學(xué)報(bào)，2002，9(1)：12-16.

[7] Rui X T, Edwin K, Bao R, et al. Discrete Time Transfer Matrix Method for Dynamics of Multibody System with Flexible Beams Moving in Space[J].Acta Mechanica Sinica, 2012, 28(2):490-504.

[8] 傅志方，華宏星.模態(tài)分析理論與應(yīng)用[M].上海:上海交通大學(xué)出版社，2000.

[9] 孫廣敏，安豐柱.基于ANSYS Workbench對(duì)立式車床橫梁的模態(tài)分析[J].機(jī)械研究與應(yīng)用，2013(4)：40-41，46.

[10] 郭 榮，周 鋐.某型轎車白車身實(shí)驗(yàn)?zāi)B(tài)分析及動(dòng)態(tài)特性評(píng)價(jià)[J].機(jī)械設(shè)計(jì)，2010(8)：18-22.

Design Research of Boring Bar Based on ODS Analysis

YUAN Li1, HU Zhong-ju1, LIANG Hong-qiang1, PENG Wei1, ZHU Ping-yu2

(1.CollegeofElectromechanicalEngineering,HunanUniversityofScienceandTechnology,XiangtanHunan411201,China;(2.SchoolofMechanicalandElectricEngineering,GuangzhouUniversity,GuangzhouGuangdong510006,China)

Based on vibration equation, a ODS formula of boring bar was deduced with transfer matrix method and mode superposition method. An example was used to verify the accuracy of this method. A general boring bar and a damping boring bar were selected as objects of study. The ODS, modal contribution and frequency response curve of this two kinds of boring bar were researched rapidly with this method. Comparisons of the results were used to direct structural design of boring bar. The results show that ODS analysis based on this ODS formula is a new method to analyze operational vibration characteristic of boring bar rapidly. It provides a new reference basis for structure parameter design calculation of boring bar.

ODS; boring bar; modal contribution; frequency response curven

2014-01-24

國家自然科學(xué)基金資助項(xiàng)目(編號(hào)：51105140)；湖南科技大學(xué)研究生創(chuàng)新基金項(xiàng)目(編號(hào)：S120016)

袁 黎(1987-)，男，湖南郴州人，在讀碩士，主要從事先進(jìn)制造工藝與裝備研究方面的工作.

TH113.1

A

1007-4414(2014)02-0017-03