□ 胡世軍 □ 王恩廣 □ 劉 學(xué) □ 侯劍波 □ 陳玉榮

1.蘭州理工大學(xué)數(shù)字制造技術(shù)與應(yīng)用省部共建教育部重點實驗室 蘭州 730050

2.蘭州理工大學(xué)機(jī)電工程學(xué)院 蘭州 730050

立柱作為數(shù)控滾齒機(jī)的關(guān)鍵部件，其結(jié)構(gòu)是否合理對機(jī)床的工作穩(wěn)定性、加工精度、使用壽命、加工效率等衡量機(jī)床的各個性能指標(biāo)都有很大程度上的影響［1］。滾齒機(jī)工作過程中，尤其是間斷非連續(xù)切削時，切削力是立柱產(chǎn)生振動的主要振源，當(dāng)激振頻率等于或者接近立柱固有頻率時將會引起立柱共振，使機(jī)床結(jié)構(gòu)剛度降低、動態(tài)特性變差，從而影響機(jī)床加工精度，因此有必要對其進(jìn)行詳細(xì)的動態(tài)分析［2］。由于立柱的結(jié)構(gòu)較為復(fù)雜，采用一般方法計算或是儀器進(jìn)行測量，難度較大，因此采用Pro/E、ANSYS軟件為平臺，充分利用兩者強(qiáng)大的三維設(shè)計、有限元分析功能及無縫連接的優(yōu)點［3］，對滾齒機(jī)立柱進(jìn)行有限元建模、模態(tài)分析，實現(xiàn)提高滾齒機(jī)立柱動態(tài)特性的目的。

1 模態(tài)分析的基本理論

在對數(shù)控滾齒機(jī)立柱動態(tài)特性分析中，物體的模態(tài)參數(shù)包括固有頻率和相對應(yīng)的振型，這兩個參數(shù)主要由物體的質(zhì)量、剛度和結(jié)構(gòu)決定，因此,可根據(jù)自由振動基本微分方程［4］對其分析。

式中：M為質(zhì)量矩陣；X為速度矢量；X為加速度矢量；C為阻尼矩陣；K為剛度矩陣；X為位移向量；F（t）為作用力向量；t為時間。

當(dāng)結(jié)構(gòu)阻尼較小時，對固有頻率與振型影響較小，可忽略不計。因此當(dāng)F（t）=0時，方程變?yōu)椋?/p>

結(jié)構(gòu)自由振動時，結(jié)構(gòu)上的各個節(jié)點作簡諧振動。它的解可以假設(shè)為以下形式:

式中：Φ為n階向量；ω為振動頻率；t0為由初始條件確定的時間常數(shù)。

將式（3）代入式（2）,得到一個廣義特征值問題,即：

結(jié)構(gòu)固有頻率高，說明單位質(zhì)量的剛度高，可作為結(jié)構(gòu)動態(tài)設(shè)計的優(yōu)化目標(biāo)，借鑒機(jī)械結(jié)構(gòu)學(xué)理論及設(shè)計經(jīng)驗，利用ANSYS Workbench軟件分析計算，改進(jìn)立柱結(jié)構(gòu)，提高剛度ki，盡可能減輕 mi，則固有頻率提高，從而使立柱的動態(tài)特性變優(yōu)。

2 立柱的有限元分析

2.1 立柱實體模型的建立

建立有效而正確的模型是進(jìn)行模態(tài)分析的關(guān)鍵環(huán)節(jié)，模型簡化正確與否直接關(guān)系到有限元計算結(jié)果的準(zhǔn)確性。因此，建立的有限元模型必須具有足夠的準(zhǔn)確性，不僅能反映機(jī)床整體的實際結(jié)構(gòu)，而且邊界條件、結(jié)構(gòu)簡化的設(shè)置要與實際情況一致。

利用Pro/E軟件建立滾齒機(jī)立柱的三維數(shù)字模型，來更真實地模擬立柱受力情況。由于立柱的實際結(jié)構(gòu)較復(fù)雜，因此對于明顯不會影響立柱強(qiáng)度、整體剛度的部位，如倒角、圓角、螺釘孔等予以簡化［5］，有利于網(wǎng)格的劃分,節(jié)省大量分析時間，如圖1所示。

2.2 劃分網(wǎng)格

▲圖1 立柱三維數(shù)字模型

如圖2所示，用Pro/E軟件建立滾齒機(jī)立柱的三維數(shù)字模型，存為STP格式導(dǎo)入ANSYS Workbench軟件中。由于立柱的結(jié)構(gòu)較為復(fù)雜，故采用智能網(wǎng)格劃分,智能網(wǎng)格劃分分為1～10級，默認(rèn)精度為6級，分析中選取精度為 6級。導(dǎo)入軟件中的立柱模型定義為10節(jié)點四面體的Solid187單元類型，對模型進(jìn)行自動劃分網(wǎng)格，網(wǎng)格劃分后，利用程序的網(wǎng)格檢查功能檢查網(wǎng)格的質(zhì)量，看是否存在畸形網(wǎng)格，進(jìn)行必要的修改。最后共劃分22 907個單元和42 093個節(jié)點。

立柱是由HT250整體澆鑄而成，HT250的金相組織主要由片狀石墨與珠光體組成，其在強(qiáng)度、耐磨性、熱穩(wěn)定方面均有較好的表現(xiàn)，其減振性能也較為優(yōu)異，同時還具有鑄造成本低、生產(chǎn)工藝簡單、生產(chǎn)周期短的優(yōu)勢［6］，其具體屬性見表 1。

表1 材料屬性

2.3 加載及求解

（1）約束條件。在典型的模態(tài)分析中唯一地有效載荷是零位移約束。如果指定了一個非零位移約束，程序?qū)⒁粤阄灰萍s束替代非零位移約束，同時也可以施加除位移約束之外的其它載荷，但它們將被忽略。根據(jù)機(jī)床實際的安裝與使用條件，可以把機(jī)床立柱的底部作為完全約束。

（2）模態(tài)提取。模態(tài)分析就是確定結(jié)構(gòu)的固有頻率和振型，從而有效預(yù)估結(jié)構(gòu)的振動特性。優(yōu)化結(jié)構(gòu)設(shè)計，利用ANSYS Workbench分析中的 Block Lanczos（蘭索斯）法對立柱進(jìn)行分析，采用Block Lanczos法提取模態(tài)，該方法適用于大型對稱特征求解問題,不僅可以很好地處理剛體振型,而且對于中、大型模型（50 000～100 000個自由度）的求解也非常有效［7］。

▲圖2 立柱有限元模型

（3）頻率選擇。一般由機(jī)床內(nèi)部振源引起的受迫振動頻率范圍約在0～500 Hz之間，而機(jī)床外部振源引起的受迫振動頻率,一般都是低頻振動,在幾Hz至幾十Hz范圍間，同時考慮到消除剛體振型,故將頻率取為1～500 Hz。

（4）階數(shù)選擇。由振動實驗表明,多自由度系統(tǒng)在初始干擾下自由振動時,總是低階振型占優(yōu)勢,故在保證一定計算精度的前提下,為了簡化計算,某些高階振型就可以忽略,主要研究低階振型。同樣對于具有無限多自由度的機(jī)床立柱振動系統(tǒng),主要對其低階振型進(jìn)行分析。本文在計算過程中，提取前五階振型進(jìn)行分析。

3 模態(tài)分析結(jié)果

用ANSYS Workbench對機(jī)床立柱劃分網(wǎng)格建立有限元模型,并進(jìn)行了必要的分析設(shè)置后,使用Solution求解器對其求解。經(jīng)計算，立柱的前五階固有頻率和振型見表2，相應(yīng)各階模態(tài)振型圖如圖3所示。

表2 立柱固有頻率與振型

▲圖3 立柱前五階模態(tài)振型云圖

激振頻率計算式［8］：

式中：n為主軸轉(zhuǎn)速；z為刀具齒數(shù)；c為采集通道數(shù)。

已知機(jī)床主軸轉(zhuǎn)速范圍為 50～500 r/min，z=4，c=3，可得作用在機(jī)床上的激振頻率范圍是10～100 Hz。

從立柱的前五階振型云圖分析可知，立柱的固有頻率避開了激振頻率，能有效防止結(jié)構(gòu)發(fā)生共振，但是立柱的上半部分變形比較大，從上往下逐漸減小，變形的最大位移基本處于立柱頂部靠前的區(qū)域。

立柱的固有頻率都在 100 Hz以上，其中前兩階屬于整體振型，立柱只是有稍微的擺動，說明立柱結(jié)構(gòu)的整體剛度較好。三階、四階和五階模態(tài)顯示的是局部振型，有著明顯的扭曲或彎曲，說明立柱局部剛度不均，這主要是立柱局部材料分布不合理，造成立柱的局部出現(xiàn)薄弱模態(tài)?？紤]到立柱上端懸掛主軸箱，將受到一個較大的彎曲力矩，因此必須通過安裝加強(qiáng)筋來加強(qiáng)立柱剛度，減少振動和變形。由于立柱內(nèi)部要留出空間來安裝配重，可以選用井字形加強(qiáng)筋，確保滾齒機(jī)立柱具有穩(wěn)定、強(qiáng)韌的基礎(chǔ)。

4 結(jié)論

利用ANSYS Workbench軟件實現(xiàn)了滾齒機(jī)立柱的參數(shù)化建模及有限元模型生成，并能夠根據(jù)結(jié)構(gòu)受力和約束情況較真實的進(jìn)行動態(tài)特性分析，獲得了立柱的固有頻率和振型，從中發(fā)現(xiàn)問題,及時消除隱患，縮短了產(chǎn)品的開發(fā)周期和投資成本。在滿足強(qiáng)度和剛度的條件下，可以減輕立柱結(jié)構(gòu)的質(zhì)量,進(jìn)行優(yōu)化設(shè)計。而ANSYS Workbench軟件可以幫助設(shè)計人員有效解決這些問題，并可以進(jìn)行多種設(shè)計方案比較，最終得出優(yōu)化設(shè)計方案。

［1］ 諸乃雄.機(jī)床動態(tài)設(shè)計原理與應(yīng)用［M］.上海:同濟(jì)大學(xué)出版社，1987.

［2］ 楊勇，王時龍，田志峰，等.大型數(shù)控滾齒機(jī)立柱動力學(xué)仿真分析［J］.中國機(jī)械工程，2013,24（11）:1473-1479.

［3］ Moaveni S.Finite Element Analysis-theory and Application with ANSYS［M］.NJ：Prentice Hall,2003.

［4］ 劉延柱，陳文良，陳立群，等.振動力學(xué)［M ］.北京：高等教育出版社，1998.

［5］ 謝黎明，王春彬，靳嵐，等.基于ANSYS的直驅(qū)擺頭動力刀架模態(tài)分析［J］.機(jī)械設(shè)計與制造，2013（4）:131-133.

［6］ 孫廣敏，安豐柱.基于ANSYS Workbench數(shù)控立式車床立柱模態(tài)分析［J］.現(xiàn)代制造技術(shù)與裝備，2013（4）:11-12.

［7］ 陳遠(yuǎn)龍，張亮，萬勝美，等.基于ANSYS的立式電解加工機(jī)床關(guān)鍵部件模態(tài)分析 ［J］.組合機(jī)床與自動化加工技術(shù)，2009（5）:42-45.

［8］ 劉杰，王海軍，王可，等.數(shù)控螺桿銑床床身結(jié)構(gòu)的動態(tài)優(yōu)化設(shè)計［J］.組合機(jī)床與自動化加工技術(shù)，2011 （3）:83-85.

［9］ 唐鵬飛，郭旭紅，張明利，等.基于ANSYS的重型機(jī)床立柱的靜態(tài)和模態(tài)分析［J］.機(jī)械制造與自動化，2012,41（6）:111-114.

［10］ 楊明亞，楊濤，湯本金，等.機(jī)床立柱動態(tài)特性的分析［J］.機(jī)械制造與自動化，2007,36 （1）:29-31.

［11］申保俊.數(shù)控機(jī)床立柱結(jié)構(gòu)有限元分析與優(yōu)化設(shè)計研究［J］.現(xiàn)代機(jī)械，2013（6）:6-9.