梁瑜洋,王益軒,高 丹,陳榮榮

(西安工程大學(xué) 機電工程學(xué)院,陜西 西安 710048)

龍門式數(shù)控銑床屬于重型、高精度、高自動化機床,在批量加工大型工件、空間曲面和特型零件等方面有廣泛應(yīng)用.自20世紀70年代以來,我國龍門銑床的加工功能和效率都有大幅提升[1].然而與發(fā)達國家相比,國內(nèi)龍門銑床在高速高效化、精密度和多軸聯(lián)動等方面仍有較大差距.目前,我國銑床設(shè)計多采用經(jīng)驗與類比方法[2],周期長、成本高、材耗大,而且很難滿足精度要求.因此必須借助現(xiàn)代先進設(shè)計方法,優(yōu)化龍門銑床結(jié)構(gòu),使其在滿足強度、剛度及振動特性要求的同時實現(xiàn)輕量化設(shè)計.龍門結(jié)構(gòu)是銑床的主要承力部件[3],其靜態(tài)和動態(tài)特性直接影響銑床的性能,因此必要對其動、靜態(tài)特性進行分析.

有限元仿真是現(xiàn)代機械結(jié)構(gòu)設(shè)計的重要技術(shù)手段,為現(xiàn)代機床設(shè)計提供了有效途徑[4].很多學(xué)者利用有限元技術(shù)對龍門銑床提出了優(yōu)化方案.例如文獻[2]通過對某龍門銑床橫梁部件的靜力分析,改進其靜態(tài)剛度,提高銑床的定位與加工精度.但未對龍門整體的靜動態(tài)特性進行綜合分析評價.本文利用虛擬樣機技術(shù),對龍門整體結(jié)構(gòu)進行靜動態(tài)分析,并預(yù)測了其疲勞特性,為龍門銑床的結(jié)構(gòu)設(shè)計提供詳細參考.

1 構(gòu)建數(shù)控銑床三維模型

該銑床模型采用工作臺固定,龍門架和橫梁移動式,包括床身、工作臺及龍門結(jié)構(gòu).龍門結(jié)構(gòu)又包括龍門框架、動梁、滾珠絲杠、托板、滑枕和數(shù)控銑頭等部件.銑床由滑臺支撐整個龍門結(jié)構(gòu),滑臺運動帶動龍門結(jié)構(gòu)實現(xiàn)縱向(X軸)進給;托板沿動梁導(dǎo)軌和絲杠運動,帶動滑枕、銑刀實現(xiàn)橫向(Y軸)進給;動梁通過立柱上的滑動導(dǎo)軌及滾珠絲杠傳動,實現(xiàn)銑床垂直方向(Z軸)運動.三個坐標軸由伺服電機控制聯(lián)動,完成工件的銑削加工.龍門銑床總體裝配模型如圖1所示.

圖1 龍門銑床總體裝配模型

2 ANSYS 仿真分析

ANSYS Workbench 14.5是ANSYS公司提供的協(xié)同仿真環(huán)境,它可與Pro/E、CATIA等三維軟件無縫連接,將三維模型直接導(dǎo)入該平臺,方便地進行各種分析.因此,可以充分利用該功能,對龍門結(jié)構(gòu)進行靜、動態(tài)特性分析,校核其強度、剛度,并獲取其振動特性.

2.1 龍門結(jié)構(gòu)有限元模型

龍門結(jié)構(gòu)高4 800mm,橫向尺寸為7 924mm,立柱截面長2 000mm,寬1 200mm.由于銑床的立柱和橫梁的縱向尺寸遠大于截面尺寸,為減輕質(zhì)量,提高固有頻率,改善動態(tài)特性,并保證足夠的靜剛度,在其表面合適位置布置通孔,內(nèi)部設(shè)置有水平方向和垂直方向的筋板[5].立柱、定梁和滑枕要求有一定的抗振性能,故選用阻尼大、吸振性好的灰鑄鐵材料HT300(密度7 350kg/m3,彈性模量130GPa,泊松比0.27),其余部分選用結(jié)構(gòu)鋼(密度7 850kg/m3,彈性模量200GPa,泊松比0.3).

在仿真分析中,有限元模型的優(yōu)劣直接影響仿真結(jié)果的準確度.建立合適的有限元模型,要兼顧準確性和計算經(jīng)濟性.由于龍門結(jié)構(gòu)模型過于復(fù)雜,為節(jié)約求解時間,必須對模型進行適當簡化.因此,計算前,先將倒角、圓角、凸臺等小尺寸特征刪除,保留立柱、橫梁、托板和滑枕等零件,并用光杠代替絲杠.為了求解精確,本文采用映射面網(wǎng)格劃分法,利用sizing控制網(wǎng)格尺寸,為立柱、定梁和滑枕添加材料HT300,為動梁等其他結(jié)構(gòu)添加材料結(jié)構(gòu)鋼,設(shè)置單元類型,劃分網(wǎng)格,最終得到龍門結(jié)構(gòu)的有限元模型如圖2所示.

圖2 龍門結(jié)構(gòu)有限元模型

2.2 結(jié)構(gòu)靜力學(xué)分析

靜力分析用于計算結(jié)構(gòu)在固定載荷作用下的響應(yīng),通過對龍門結(jié)構(gòu)進行靜力分析,獲得其位移和應(yīng)力分布云圖,便于確定其剛度和強度.銑床工作過程中,切削力依次由銑刀、主軸傳動系統(tǒng)、滑枕、托板、橫梁傳至立柱,方向主要沿X方向和Y方向,而Z方向受力相對較小.根據(jù)實際大型銑床切削力范圍,取FX=FY=10 000N,FZ=2 000N.托板在橫梁中間位置時,橫梁靜剛度最差,故只取該工況進行分析.將銑刀所受切削力等效到滑枕上,同時考慮龍門自重對結(jié)構(gòu)的影響,劃分網(wǎng)格,求解,得到相應(yīng)的應(yīng)力云圖(如圖3)和位移云圖(如圖4).

為防止龍門在復(fù)雜應(yīng)力狀態(tài)下發(fā)生疲勞破壞,需要對其進行疲勞強度分析.結(jié)構(gòu)在對稱循環(huán)應(yīng)力狀態(tài)下最為危險,因此將疲勞分析應(yīng)力比定為-1.查閱文獻[6],得到結(jié)構(gòu)鋼(本文選用45鋼)的S-N曲線見圖5.分析完成后,求得龍門壽命和安全系數(shù)如圖6(a)～(b)所示.

圖3 靜力分析應(yīng)力云圖 圖4 靜力分析位移云圖

圖5 結(jié)構(gòu)鋼S-N曲線

(a) 龍門結(jié)構(gòu)壽命云圖 (b) 龍門結(jié)構(gòu)安全系數(shù)圖6 疲勞強度分析結(jié)果

由應(yīng)力云圖可知,龍門結(jié)構(gòu)整體均處于較低的應(yīng)力狀態(tài),在立柱與橫梁連接導(dǎo)軌處,應(yīng)力達到峰值,最大等效應(yīng)力為33.765MPa,遠遠小于HT300和結(jié)構(gòu)鋼的極限應(yīng)力,故該方案完全滿足強度要求,并且可以適當減小壁厚,或減少立柱和橫梁中肋板的數(shù)量,從而減輕結(jié)構(gòu)重量,達到優(yōu)化設(shè)計的目的.

龍門的切削變形很大程度上決定了銑床的加工精度.由位移云圖可知,立柱、定梁及動梁的變形很小,均在0.06mm以內(nèi),滑枕頂端變形最大,也僅為0.111 83mm.因此龍門結(jié)構(gòu)整體具有很好的切削剛度,從而保證了銑床的加工精度要求.在實際龍門銑床中,滑枕是安裝主軸、傳動軸等零件,限制去除率切削的關(guān)鍵部件,懸臂長,撓度大,抗彎性差,必須有足夠的剛度和吸振性[7],因此可適當增加壁厚或肋板,或采用其他優(yōu)質(zhì)材料.進一步查看重力與切削力耦合變形結(jié)果可知,滑枕Y方向變形較大,需要改善剛度.

由圖6結(jié)果可知,龍門結(jié)構(gòu)各部分壽命較為均勻,在107應(yīng)力循環(huán)以上,基本滿足等壽命設(shè)計要求.各部分疲勞安全系數(shù)最小值為4,符合疲勞強度設(shè)計要求,但是存在較大的剩余強度,需進行結(jié)構(gòu)優(yōu)化.

2.3 模態(tài)分析

模態(tài)分析是解決復(fù)雜結(jié)構(gòu)振動問題的主要方法,通過模態(tài)分析,獲取系統(tǒng)各階固有頻率及振型,使外界激勵頻率避開系統(tǒng)的固有頻率,可以有效減小振動幅值,避免共振破壞.文中通過模態(tài)分析,得到龍門結(jié)構(gòu)前20階固有頻率見表1,選取結(jié)果中典型振型見表2.

表1 龍門結(jié)構(gòu)的前20階固有頻率

表2 龍門振型描述

由分析結(jié)果可得出以下結(jié)論:(1) 滑枕和定梁的動態(tài)特性較差,尤其需要注意防止變形損壞.(2) 由于第1階振型固有頻率較低,為24.828Hz,當主軸轉(zhuǎn)速在1 490r/min附近時將發(fā)生共振,導(dǎo)致銑床工作異常,嚴重影響加工精度,甚至造成破壞.在不改變銑床結(jié)構(gòu)情況下,應(yīng)合理選擇主軸電機和主軸箱速比,使主軸轉(zhuǎn)速低于1 350r/min,或運行在各階固有頻率之間,從而使激振頻率避開結(jié)構(gòu)各階固有頻率,確保銑床在工作頻率范圍內(nèi)正常運轉(zhuǎn),并減小切削振動幅值,提高零件加工精度和表面加工質(zhì)量.

(a) 第一階振型 (b) 第三階振型圖7 龍門結(jié)構(gòu)振型圖

2.4 諧響應(yīng)分析

龍門銑床在工作過程中,將受到復(fù)雜的空間載荷,對其進行諧響應(yīng)分析[8],可以得到在不同頻率下龍門框架的響應(yīng)幅值,預(yù)測結(jié)構(gòu)的持續(xù)動力學(xué)特性,避免因共振、疲勞或受迫振動造成破壞.

在ANSYS Workbench 14.5中,為模型添加約束與載荷后,將掃頻范圍設(shè)置為0～120Hz,步數(shù)設(shè)為30,求得在外激振載荷作用下,結(jié)構(gòu)的位移響應(yīng)曲線及應(yīng)力響應(yīng)曲線如圖8(a)～(b)所示.

(a) 位移響應(yīng) (b) 應(yīng)力響應(yīng)圖8 諧響應(yīng)分析頻譜圖

由圖8(a)可知,外界激勵載荷在32Hz,44Hz,72Hz和112Hz附近時,位移響應(yīng)曲線達到峰值,特別在72Hz和112Hz時,振幅最大,應(yīng)使激勵載荷遠離此頻率,避免龍門結(jié)構(gòu)出現(xiàn)共振破壞或影響加工質(zhì)量.

由圖8(b)可知,外界激勵載荷在24Hz,44Hz,72Hz和112Hz附近時,應(yīng)力響應(yīng)曲線出現(xiàn)峰值,特別在44Hz和112Hz時,響應(yīng)幅值最大,應(yīng)使激勵載荷避開此頻率,避免應(yīng)力超過屈服極限造成結(jié)構(gòu)破壞.

2.5 響應(yīng)譜分析

響應(yīng)譜分析是一種基于模態(tài)分析的頻域分析,它可替代瞬態(tài)分析,用來獲取結(jié)構(gòu)在瞬態(tài)載荷作用下的最大響應(yīng)區(qū)域和響應(yīng)幅值[9].由模態(tài)分析結(jié)果,主軸轉(zhuǎn)速應(yīng)低于1 350r/min,為模擬數(shù)控銑床實際工作中所受瞬態(tài)載荷,查資料[10],為銑床添加Z方向的功率譜加速度見表3.在ANSYS Workbench 14.5中,選擇SRSS模態(tài)組合方法.對銑床進行響應(yīng)譜分析,結(jié)果如圖9(a)～(b)所示.

(a) Y方向位移云圖 (b) 龍門結(jié)構(gòu)應(yīng)力云圖圖9 響應(yīng)譜分析結(jié)果

表3 切削力加速度頻譜

在切削力載荷作用下,銑床最大變形區(qū)主要集中在滑枕和定梁部分,變形幅值為0.061 8mm.最大應(yīng)力區(qū)域出現(xiàn)在動梁及動梁和立柱的連接處,最大等效應(yīng)力為4.025 6MPa.該結(jié)果可為系統(tǒng)進一步優(yōu)化,減小振動和變形、降低噪聲和疲勞等提供依據(jù).

3 結(jié)束語

本文利用虛擬樣機技術(shù),對龍門銑床的關(guān)鍵部件龍門進行了靜態(tài)和動態(tài)特性分析,并針對仿真結(jié)果提出了合理的優(yōu)化建議,為數(shù)控龍門銑床設(shè)計優(yōu)化提供了參考.通過有限元分析計算,還可以定性定量地對橫梁、立柱、滑枕等重要部件的靜力學(xué)、振動特性和受迫激振響應(yīng)進行改進.實踐證明,該方法對大型機床設(shè)計中提高精度、降低成本具有重要價值.

參考文獻：

[1] 于亮.大型數(shù)控龍門銑床主軸柔性傳動系統(tǒng)優(yōu)化設(shè)計方法研究[D].沈陽:沈陽工業(yè)大學(xué),2013:1-3.

[2] 楊曼云.XK2130大(重)型數(shù)控龍門銑床橫梁性能有限元分析與研究[J].機電產(chǎn)品開發(fā)開發(fā)與創(chuàng)新,2010,23(1):170-172.

[3] 王一江,張曉燕,龐學(xué)慧.基于ANSYS的龍門銑床龍門結(jié)構(gòu)模態(tài)分析[J].計算機應(yīng)用技術(shù),2011,37(4):48-50.

[4] 陳旭,王益軒,李端,等.大型龍門數(shù)控纖維纏繞/鋪放機動態(tài)分析[J]. 西安工程大學(xué)學(xué)報,2013,27(2):215-217.

[5] WANG Yixuan,GUO Yanli,WANG Ying.Research on dynamics simulation of CNC milling machine based on virtual prototyping [J].International Review of Mechanical Engineering,2013，7(5):821-822.

[6] 秦大同,謝里陽.現(xiàn)代機械設(shè)計手冊(單行本)[M].北京:化學(xué)工業(yè)出版社,2013:85-95.

[7] 黃煜.從滑枕主傳動結(jié)構(gòu)設(shè)計上提高龍門銑床剛性和精度[J].裝備制造技術(shù),2008(10):120-121.

[8] 凌桂龍,丁金濱,溫正,等.ANSYS Workbench 13.0從入門到精通[M].北京:清華大學(xué)出版社,2012：57-65.

[9] 高華生,龍芋宏,杜懷輝.汽車駕駛室響應(yīng)譜分析[J].科學(xué)技術(shù)與工程,2012,12(35):9769-9772.

[10] 吳衛(wèi)國.高效精密切削及其振動特性的研究[D].鎮(zhèn)江:江蘇大學(xué),2007:65-68.