王昆， 劉進福， 韓晶， 莫賢

（常州先進制造技術(shù)研究所，江蘇常州213164）

0 引言

目前，在企業(yè)組織產(chǎn)品焊接中，大都采用手工焊接，其不足之處在于：手工焊接質(zhì)量不穩(wěn)定，工人勞動強度大；焊接時產(chǎn)生的煙霧對工人的呼吸系統(tǒng)產(chǎn)生損害，焊接時發(fā)出的亮光會對工人的眼角膜產(chǎn)生損害；手工焊接的生產(chǎn)效率低。使用焊接機器人不僅能降低焊接對工人身體的危害，還能夠有效地提高工作效率，降低勞動成本。

焊接機器人的大臂是整個機器人受力最大的機構(gòu)，作為機器人的重要部件，其剛度和強度直接影響到機器人的使用精度和壽命。在設(shè)計過程中必須保證它的強度和剛度及在運動過程中振動頻率的影響。本文在有限元受力分析軟件ANSYS下對焊接機器人的大臂進行靜力分析，為機器人整體設(shè)計和優(yōu)化提供有效的理論支持。

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

有限元的模態(tài)分析就是建立模態(tài)模型進行數(shù)值分析的過程。對于一般的多自由度結(jié)構(gòu)系統(tǒng)而言，運動都可以由其自由振動的模態(tài)來合成。由于系統(tǒng)結(jié)構(gòu)阻尼主要發(fā)生在固定的或可動的連接部位，因此零件材料的內(nèi)阻尼只是吸收和損耗能量的一部分［5］。系統(tǒng)自由振動系統(tǒng)方程為：

其中：［M］和［K］分別為油缸坯料搬運系統(tǒng)的質(zhì)量矩陣和剛度矩陣，均為N×N階方陣；｛X｝表示油缸坯料搬運系統(tǒng)的位移向量，同時還可以將式（1）寫成位移向量的形式：

經(jīng)線性變換式｛x｝=［u］｛y｝對集合位置坐標｛x｝表示的耦合系統(tǒng)微分方程組進行解耦。因此，振型在坐標變換和解耦系統(tǒng)中發(fā)揮著重要的作用。為了得到振型的矩陣［u］，必須通過系統(tǒng)的特征值和特征向量，因此,假設(shè)系統(tǒng)的振動是一系列諧波振動頻率。系統(tǒng)運動方程為：

通過微分方程得：

將式（4）代入式（1），可得特征方程：

將括號中的矩陣行列式等于零，可得特征方程式即式（6）：

上式可以寫成以下形式：

2 有限元模型的建立

焊接機器人模型比較復雜，運用ANSYS Workbench軟件很難實現(xiàn)模型的三維建模，本文采用專業(yè)的三維建模軟件SolidWorks建立了機器人的模型，然后將模型另存為x_t文件，導入到ANSYS軟件中，模型的某些特征在劃分網(wǎng)格時，會造成網(wǎng)格劃分失敗，因此，在導入ANSYS軟件之前，應(yīng)刪除與分析意圖無關(guān)的特征，簡化三維模型，依據(jù)圣維南原理，采用等效質(zhì)量體代替。

1）定義材料屬性，劃分網(wǎng)格。大臂材料為鋁合金，彈性模量為10GPa，密度2700kg/m3，屈服極限275MPa，采用自由網(wǎng)格劃分，共有單元數(shù)40199，節(jié)點數(shù)68418，如圖1所示。

圖1 劃分網(wǎng)格

2）定義約束。大臂在實際工作過程中，只有繞Z軸的轉(zhuǎn)動，將全局坐標系轉(zhuǎn)化成柱坐標系，便于限制Z方向和R方向的位移。

3）定義載荷。最左側(cè)空腔用于安裝電機、減速機，為了簡化模型，將其刪除。添加50N的預(yù)緊力，最右側(cè)的安裝面用于安裝小臂等部件，在此添加120N的預(yù)緊力。

3 分析求解

1）靜力分析。將上述所有負載施加到建立好的大臂有限元模型上，并添加好限制性約束，進行求解計算，得到大臂的內(nèi)部應(yīng)力分布云圖和形變位移分布云圖，如圖2和圖3所示。

圖2 應(yīng)力圖解

圖3 位移圖解

從圖2的大臂應(yīng)力分布云圖中可以看出，大臂的最大應(yīng)力僅為0.85 MPa，位于機器人電機安裝位置，遠小于大臂的屈服極限，滿足設(shè)計要求。在圖3所示的大臂位移分布云圖中，大臂在最大靜載荷作用下的最大變形為1.18×10-5mm，該形變量非常微小，對運動機構(gòu)整體的運動精度影響很小。綜上所述，所設(shè)計的大臂結(jié)構(gòu)具有足夠的機械強度，完全符合使用要求。

2）模態(tài)分析。預(yù)應(yīng)力模態(tài)分析各階振型如圖4所示。根據(jù)模態(tài)分析的結(jié)果，可以得出：（1）大臂的固有頻率較高，一階固有頻率已達到136.89 Hz，運動機構(gòu)所采用的伺服電機額定轉(zhuǎn)速為3000r/min，其回轉(zhuǎn)頻率為50 Hz，遠小于大臂的一階固有頻率，因此大臂本身的剛度已滿足設(shè)計需要，不需再提高其固有頻率。（2）大臂的固有頻率雖已達到設(shè)計要求，但就大臂本身而言，由于其跨度大，左側(cè)電機安裝部分仍是其相對薄弱的地方。

圖4 模態(tài)振型

表1 X軸大臂的固有頻率和振型特征

4 結(jié)論

1）運用ANSYS建立了焊接機器人大臂的有限元模型，對整個結(jié)構(gòu)進行了靜力學分析，得到了大臂的應(yīng)力分布和位移分布：其中最大應(yīng)力位于大臂的中間部位，遠小于大臂的屈服極限，最大變形位于大臂前端，變形量為1.18×10-5mm，對整個系統(tǒng)的精度影響很?。?/p>

2）通過對大臂的模態(tài)分析發(fā)現(xiàn)：大臂的一階固有頻率遠遠大于電機的回轉(zhuǎn)頻率，大臂的剛度滿足設(shè)計要求；

3）通過對大臂的計算分析，大臂能夠滿足焊接機器人的使用要求，保證控制精度。

［參考文獻］

［1］張功學，田楊.基于ANSYS Workbench的變速自行車車架的有限元分析［J］.微型機與應(yīng)用，2009（6）：63-65.

［2］呂廷，石秀東，張秋菊，等.基于ANSYS的破碎機機架模態(tài)分析［J］.機械設(shè)計與制造，2008（11）：99-101.

［3］章華，孔曉玲，趙文龍，等.基于ANSYS Workbench的連續(xù)障礙圓餅的有限元分析［J］.機械工程師，2011（2）：81-83.

［4］陳金玉，楊來俠.基于ANSYS Workbench手機外殼有限元網(wǎng)格劃分研究［J］.現(xiàn)代制造技術(shù)與裝備，2008（1）：58－60.