□ 張世赟 □ 芮執(zhí)元 □ 李向龍

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

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

電解鋁預(yù)焙陽(yáng)極導(dǎo)電裝置專(zhuān)用摩擦焊機(jī)能夠?qū)崿F(xiàn)鋁導(dǎo)桿和鑄鋼爪的直接焊接，焊接接頭強(qiáng)度好，對(duì)接精度高，焊接過(guò)程高效、環(huán)保。作為焊機(jī)核心部件的主軸箱，承載著整個(gè)主軸系統(tǒng)，結(jié)構(gòu)復(fù)雜，動(dòng)、靜載荷分布不均，特別是要承受很大的頂鍛力作用。其剛度、強(qiáng)度及動(dòng)態(tài)特性直接影響到整機(jī)的加工精度、可靠性及使用壽命。

筆者以預(yù)焙陽(yáng)極導(dǎo)電裝置專(zhuān)用摩擦焊機(jī)為研究對(duì)象，應(yīng)用Pro/E軟件建立其有限元模型，根據(jù)實(shí)際情況添加約束，并以無(wú)預(yù)應(yīng)力狀態(tài)為邊界條件，進(jìn)行動(dòng)態(tài)特性分析，得到主軸箱前六階固有頻率和振型，并進(jìn)行模態(tài)和隨機(jī)振動(dòng)分析，分析結(jié)果對(duì)專(zhuān)用摩擦焊機(jī)主軸箱的優(yōu)化設(shè)計(jì)提供了重要依據(jù)。

1 主軸箱三維模型建立

為了利用Pro/E完成主軸箱三維實(shí)體建模，要在保證分析結(jié)果精確度的情況下，對(duì)主軸箱局部特征進(jìn)行適當(dāng)簡(jiǎn)化。筆者對(duì)主軸箱部分小孔和倒角進(jìn)行了適當(dāng)?shù)暮?jiǎn)化，完成后的主軸箱三維模型如圖1所示。

▲圖1 主軸箱三維模型

▲圖2 主軸箱網(wǎng)格劃分模型

2 主軸箱有限元模型建立

2.1 有限元法的基本原理

有限元法是把連續(xù)介質(zhì)離散成一組單元，使無(wú)限自由度問(wèn)題轉(zhuǎn)化成有限自由度問(wèn)題，利用每個(gè)離散單元內(nèi)的近似函數(shù)來(lái)分別表示求解域上需要求解的未知場(chǎng)函數(shù)，所有離散單元場(chǎng)函數(shù)的集合能近似表示整個(gè)連續(xù)體的場(chǎng)函數(shù)，從而求解此離散方程組。

求解的基本思路是將求解給定的泊松方程問(wèn)題轉(zhuǎn)化為求解泛函的極值問(wèn)題。

2.2 定義材料類(lèi)型及單元屬性

主軸箱材料采用了低合金高強(qiáng)度結(jié)構(gòu)鋼16Mn，材料的主要特點(diǎn)是冷沖壓、低溫、焊接及可切削性能較好，是目前我國(guó)應(yīng)用最為廣泛的低合金鋼，密度ρ=7 870 kg/m3；彈性模量 E=212 GPa，泊松比 μ=0.310。

2.3 劃分有限元網(wǎng)格

網(wǎng)格劃分是建立有限元模型的一個(gè)關(guān)鍵環(huán)節(jié)，它需要考慮的因素較多，工作量很大，所劃分的網(wǎng)格形式對(duì)進(jìn)行有限元計(jì)算的規(guī)模和計(jì)算精確都會(huì)產(chǎn)生直接影響，利用ANSYS Workbench對(duì)主軸箱網(wǎng)格劃分的模型如圖2所示。

2.4 設(shè)定邊界條件

根據(jù)摩擦焊機(jī)實(shí)際工作狀況，將主軸箱的邊界條件設(shè)定為約束主軸箱底板的X、Y、Z向位移。

3 主軸箱模態(tài)分析

3.1 模態(tài)分析理論基礎(chǔ)

根據(jù)模態(tài)分析后得到的結(jié)果參數(shù)可以對(duì)機(jī)械結(jié)構(gòu)動(dòng)態(tài)性能進(jìn)行評(píng)估。在一個(gè)機(jī)械振動(dòng)系統(tǒng)中，一般結(jié)構(gòu)必須滿(mǎn)足工作頻率不落在某階模態(tài)的半功率帶寬內(nèi)，即各階模態(tài)頻率必須遠(yuǎn)離工作頻率。而對(duì)結(jié)構(gòu)振動(dòng)貢獻(xiàn)大的振型，必須使其不能影響系統(tǒng)的正常工作。

有限元方法和模態(tài)分析相結(jié)合為機(jī)械結(jié)構(gòu)動(dòng)態(tài)設(shè)計(jì)提供了很好的途徑。主軸箱可看作一個(gè)多自由度振動(dòng)系統(tǒng)，對(duì)于多自由度系統(tǒng)，物理坐標(biāo)系中的自由振動(dòng)響應(yīng)為n（n為自由度的數(shù)目）個(gè)主振動(dòng)的疊加，其運(yùn)動(dòng)微分方程為：

式中：［M］為質(zhì)量矩陣；［C］為阻尼矩陣；［K］為剛度矩陣；｛｝為加速度矢量；｛˙｝為速度矢量；｛x｝為位移矢量；｛F（t）｝為激振力矩陣。

進(jìn)行主軸箱系統(tǒng)動(dòng)態(tài)分析時(shí)，可將其視為不受外激振力影響，且因?yàn)樽枘釋?duì)系統(tǒng)的影響很小，可以忽略，主軸箱無(wú)阻尼自由振動(dòng)運(yùn)動(dòng)微分方程為：

與之對(duì)應(yīng)的特征方程為：

它是關(guān)于ω2的n次代數(shù)方程，其中ω2為振動(dòng)系統(tǒng)的主頻率，設(shè)無(wú)重根，解此方程得到ω的n個(gè)互異正根，即為振動(dòng)系統(tǒng)的n階主頻率。

3.2 有限元模態(tài)分析

經(jīng)ANSYS Workbench模態(tài)分析，可得主軸箱前六階固有頻率分別為：103.12 Hz、131.87 Hz、135.25 Hz、207.62 Hz、224.92 Hz、240.35 Hz，其六階模態(tài)振型分別如圖3～圖8所示。

由圖分析可知，主軸箱前六階固有頻率差值變化不大，主軸箱的低階模態(tài)頻率較高，激振頻率難以達(dá)到，可以很好地避免共振發(fā)生。第二、三階固有頻率分布較密集。通過(guò)前六階振型可觀察到，隨著模態(tài)陣型階數(shù)的升高，第一階振型主軸箱變形最大，主要是后立板沿X軸前后擺動(dòng)，隨后振型變化程度依次遞減。

4 主軸箱隨機(jī)振動(dòng)分析

隨機(jī)振動(dòng)是指振動(dòng)規(guī)律具有隨機(jī)性而無(wú)法用確定的函數(shù)表示的一類(lèi)振動(dòng)運(yùn)動(dòng)。為了對(duì)其進(jìn)行分析，要考慮運(yùn)用隨機(jī)分析手段，它是一種基于概率統(tǒng)計(jì)學(xué)的譜分析技術(shù)，在隨機(jī)激勵(lì)作用下對(duì)包括應(yīng)力、位移等物理量的概率分布進(jìn)行求解的分析手段。

本文在對(duì)專(zhuān)用摩擦焊機(jī)主軸箱模態(tài)分析的基礎(chǔ)上，利用其有限元模型作了進(jìn)一步的隨機(jī)振動(dòng)分析，分別得到1σ位移和1σ應(yīng)力響應(yīng)圖，如圖9、圖10所示。

根據(jù)所得結(jié)果分析可知，在主軸箱隨機(jī)振動(dòng)過(guò)程中，應(yīng)力分布不均勻，局部應(yīng)力較大，后立板與主軸后嵌套連接處及主軸后嵌套處發(fā)生變形的可能性較大，在優(yōu)化過(guò)程中應(yīng)予以修改。

▲圖3 一階模態(tài)振型

▲圖4 二階模態(tài)振型

▲圖5 三階模態(tài)振型

▲圖6 四階模態(tài)振型

▲圖7 五階模態(tài)振型

▲圖8 六階模態(tài)振型

▲圖9 1σ位移響應(yīng)圖

▲圖1 0 1σ應(yīng)力響應(yīng)圖

5 結(jié)論

利用Pro/E軟件建立電解鋁專(zhuān)用摩擦焊機(jī)主軸箱三維實(shí)體模型，導(dǎo)入ANSYS Workbench軟件中生成有限元模型，并進(jìn)行模態(tài)分析，得到前六階固有頻率和振型。之后進(jìn)行隨機(jī)振動(dòng)分析，得到主軸箱位移和應(yīng)力響應(yīng)圖。通過(guò)分析結(jié)果，識(shí)別主軸箱薄弱環(huán)節(jié)，對(duì)主軸箱的結(jié)構(gòu)優(yōu)化設(shè)計(jì)工作提供了重要參考，對(duì)提高其整體動(dòng)態(tài)性能具有重要意義。

［1］ 韋遼,李健.CA6140機(jī)床主軸箱模態(tài)和隨機(jī)振動(dòng)分析［J］.農(nóng)業(yè)裝備與車(chē)輛工程,2013,51（10）:49－51.

［2］ 趙建峰,魏鋒濤,宋俐.基于ANSYS的機(jī)床主軸箱結(jié)構(gòu)優(yōu)化設(shè)計(jì)［J］.機(jī)電產(chǎn)品開(kāi)發(fā)與創(chuàng)新,2010,23（5）:140－142.

［3］ D R Salgado，F(xiàn) J Alonso.Optimal Machine Tool Spindle Drive Gearbox Design［J］.International Journal of Manufacturing Technology,2008，37:851－860.

［4］ 張力,林建龍,項(xiàng)輝宇.模態(tài)分析與實(shí)驗(yàn) ［M］.北京:清華大學(xué)出版社,2011.

［5］ 凌桂龍,丁金濱,溫正.ANSYS Workbench 13.0從入門(mén)到精通［M］.北京:清華大學(xué)出版社,2012.

［6］ 張宏博,薄瑞峰,游小紅.基于拓?fù)鋬?yōu)化的深孔鉆床主軸箱的結(jié)構(gòu)重構(gòu)［J］.機(jī)床 與 液 壓 ,2013,41 （7）:150－153.

［7］ 杜勇,辛舟.5000kN連續(xù)驅(qū)動(dòng)摩擦焊機(jī)主軸箱的有限元分析［J］.機(jī)械設(shè)計(jì),2014,31（7）:51－53.

［8］ 韓江,孟超,姚銀鴿，等.大型數(shù)控落地鏜銑床主軸箱的有限元分析［J］.工藝與裝備,2009（10）:82－84.

［9］ 張曉龍,李功宇,吳俊華.TH65100主軸箱結(jié)構(gòu)動(dòng)力學(xué)分析及改進(jìn)設(shè)計(jì)研究［J］ .機(jī)床與液壓,2008,36（4） :354－356.