王海玲

（天津機電職業(yè)技術(shù)學(xué)院，天津 300350）

0 引言

液壓挖掘機在建筑、采礦及道路施工等工程領(lǐng)域有著廣泛的應(yīng)用，由于工作環(huán)境的惡劣性，挖掘機結(jié)構(gòu)件的強度和可靠性成為衡量挖掘機性能的重要指標(biāo)，對此眾多的研究者采用多體動力學(xué)、有限元等工具對挖掘機動臂、斗桿等重要構(gòu)件的強度展開仿真分析，但是缺乏仿真與工作裝置實際工況下的應(yīng)力對比，難以確保有限元模型的準(zhǔn)確性［1－3］。本文以某具體型號挖掘機為分析對象，采用貼電阻應(yīng)變片的方法，對動臂和斗桿關(guān)鍵位置進(jìn)行了應(yīng)力測試，并通過ANSYS有限元軟件對其進(jìn)行仿真計算，將仿真計算結(jié)果和測試值進(jìn)行比對，以驗證有限元模型的正確性。

1 靜應(yīng)力測試

1.1 測試方案

挖掘機應(yīng)力測量裝置如圖1所示。挖掘機工作裝置的靜態(tài)應(yīng)力測試采用在關(guān)心的位置貼電阻應(yīng)變片的方法，通過數(shù)采設(shè)備收集測試點位置的應(yīng)變，然后再根據(jù)力學(xué)原理換算為需要的應(yīng)力值。挖掘機工作載荷利用拉壓傳感器進(jìn)行加載，拉壓傳感器通過吊帶連接到挖掘機的鏟齒上，施加6×104N的拉力，力的位置和方向如圖1（a）中箭頭所示。整個應(yīng)力測試系統(tǒng)還包括LMS的數(shù)采前端和Test.Xpress測試分析軟件，如圖1（b）所示。

采用貼應(yīng)變片的方法進(jìn)行結(jié)構(gòu)靜應(yīng)力的測試，測量結(jié)果受貼片技術(shù)、測試環(huán)境、測試系統(tǒng)誤差等諸多因素的影響。為了保證測量過程的合理性和測量設(shè)備的有效性，先做簡單的驗證實驗。首先在一個矩形試塊上貼片、做靜力測試，然后建立對應(yīng)的有限元模型進(jìn)行分析，當(dāng)測試和有限元計算結(jié)果完全吻合后再實施本次測試工作。

圖1 挖掘機應(yīng)力測量裝置

1.2 測試點布置

應(yīng)力測試觀測點的布置要考慮挖掘機實際應(yīng)用的損傷點，盡量將應(yīng)變片布置在用戶關(guān)心的位置，比如動臂和斗桿上靠近鉸鏈的位置容易損壞，而斗桿底面靠近筋板的位置容易開裂。但是由于挖掘機結(jié)構(gòu)很大部分是板材焊接件，而焊接位置材料的屬性較焊接前會發(fā)生很大的變化，目前有限元分析軟件對此類問題的數(shù)值求解精確性很差，所以貼片一定要避開板材焊接處。此外貼片應(yīng)考慮將測試點布置于靜強度分析應(yīng)力較大的位置，且盡量靠近板材的中心線，因為該處的應(yīng)力變化比較均勻，能夠有效降低應(yīng)力突變引起的測量誤差［4］。綜合考慮作業(yè)工況、理論計算、受力分析及實際貼片條件，對工作裝置靜力測試點做出如下布置：在挖掘機動臂和斗桿上各設(shè)置3個測試點，即在動臂上邊、下邊、側(cè)面布置測試點1、2、3，在斗桿上邊、下邊、側(cè)面布置測試點4、5、6，各測試點位置如圖2所示。此外，由于本實驗裝置和測試系統(tǒng)要同時進(jìn)行挖掘機載荷譜的測試，實驗周期相對較長，要隨時關(guān)注由于應(yīng)變片膠水老化、貼合面開裂等問題引起的測試誤差，一旦誤差超過允許值需要重新打磨貼片。

圖2 挖掘機測試點布置

2 有限元建模與分析

2.1 三維幾何模型建立

工作裝置主要由動臂、斗桿、鏟斗組件以及驅(qū)動油缸等構(gòu)件組成，其中動臂和斗桿是最關(guān)鍵的承力部件，是有限元分析的重點，但是動臂和斗桿屬于焊接件，存在多處焊縫，因此要對其三維模型進(jìn)行修復(fù)，消除小的間隙、尖角以及實體穿透等。動臂和斗桿三維模型及內(nèi)部結(jié)構(gòu)如圖3所示。

圖3 動臂斗桿三維模型

2.2 有限元模型及材料屬性

為保證有限元分析結(jié)果的正確性，有限元網(wǎng)格質(zhì)量要滿足一定的要求，本次分析采用專業(yè)的有限元網(wǎng)格生成軟件HyperMesh進(jìn)行網(wǎng)格劃分，采用ANSYS提供的實體單元Solid185對幾何模型進(jìn)行網(wǎng)格劃分形成以六面體網(wǎng)格為主的有限元模型。ANSYS中網(wǎng)格質(zhì)量主要利用skewness（偏斜率）來度量，動臂和斗桿網(wǎng)格的偏斜率數(shù)值如圖4所示。

由圖4可知：動臂的偏斜率最大為0.78，平均值為7.35×10－2；斗桿的偏斜率最大為0.95，平均值為6.18×10－2。根據(jù)skewness網(wǎng)格質(zhì)量度量等級標(biāo)準(zhǔn)，0.95以下的為可接受，0.8以下的是較好的，0.5以下是優(yōu)秀的，總體來說網(wǎng)格質(zhì)量滿足要求。動臂和斗桿材料參數(shù)如表1所示。

圖4 動臂和斗桿網(wǎng)格質(zhì)量

表1 動臂和斗桿材料參數(shù)

2.3 有限元模型約束及載荷

通常情況下，做靜力分析要在有限元模型中對結(jié)構(gòu)添加與實際工況相一致的約束，如果結(jié)構(gòu)約束非常復(fù)雜，模型中施加的約束不合理或者偏離實際的話，會導(dǎo)致有限元分析結(jié)果出現(xiàn)較大誤差。為了規(guī)避約束方式引入的分析偏差，采用ANSYS慣性釋放功能。慣性釋放是ANSYS中的一個高級應(yīng)用，允許對完全無約束的結(jié)構(gòu)進(jìn)行靜力分析，實質(zhì)上就是用結(jié)構(gòu)的慣性力來平衡外力。盡管結(jié)構(gòu)沒有約束，分析時仍假設(shè)其處于一種“靜態(tài)”的平衡狀態(tài)。采用慣性釋放功能進(jìn)行靜力分析時，只需要對一個節(jié)點進(jìn)行6個自由度的約束（虛支座），針對虛支座，程序首先計算在外力作用下每個節(jié)點在每個方向上的加速度，然后將加速度轉(zhuǎn)化為慣性力反向施加到每個節(jié)點上，由此構(gòu)造一個平衡的力系（支座反力等于零），求解得到的位移描述所有節(jié)點相對于該支座的相對運動。

為了得到動臂和斗桿各鉸接點的載荷，在多體動力學(xué)軟件ADAMS中施加與試驗現(xiàn)場大小及方向均一致的載荷［5］，即按照圖1（a）在鏟斗齒間施加6×104N的力。經(jīng)ADAMS計算得到的各鉸接點的載荷如圖5所示。

圖5 動臂和斗桿各鉸接點載荷

采用ANSYS慣性釋放功能計算得到動臂和斗桿的拉應(yīng)力云圖，如圖6所示。由圖6可知：在當(dāng)前工況下，動臂和斗桿整體應(yīng)力分布比較均勻，動臂應(yīng)力沿縱向由中間向兩端逐漸較小，斗桿應(yīng)力較大區(qū)域出現(xiàn)在與動臂的鉸接點兩側(cè)；動臂的最大應(yīng)力值出現(xiàn)在動臂上蓋板的耳板處，該處是動臂的危險截面，而斗桿的最大應(yīng)力值出現(xiàn)在與動臂油缸的鉸接處，此處也是實際應(yīng)用中的易損點。

圖6 動臂和斗桿應(yīng)力云圖

6個測試點的有限元應(yīng)力計算值如圖7所示，各測試點應(yīng)力有限元仿真值與測試值對比如表2所示。

圖7 各測試點的應(yīng)力值

由表2可知：各測點應(yīng)力有限元仿真值與測試值最大誤差為測試點1處，差值為6.9%；最小誤差為測試點2處，差值為1.1%。這說明有限元靜應(yīng)力分析與實測值不可避免地存在誤差，誤差來源于仿真模型與實際裝置的差異、貼片測量誤差、載荷計算誤差等因素。但是仿真與測試的誤差整體在7%以內(nèi)，誤差比較小，故可以認(rèn)為有限元模型能夠正確反映實際工況中結(jié)構(gòu)的應(yīng)力分布情況。

表2 各測試點應(yīng)力有限元仿真值與測試值對比

3 結(jié)論

以某具體型號挖掘機為研究主體，采用貼電阻應(yīng)變片的方法對動臂和斗桿關(guān)鍵點進(jìn)行應(yīng)力應(yīng)變測試，并將有限元模型仿真計算結(jié)果與測試值做比對。分析結(jié)果表明：采用ANSYS慣性釋放功能可以很好地解決約束不當(dāng)造成的應(yīng)力集中和計算誤差的問題；有限元應(yīng)力計算結(jié)果和測試結(jié)果吻合較好，可以作為挖掘機工作裝置改進(jìn)和優(yōu)化的參考依據(jù)。進(jìn)一步工作是增加應(yīng)力測試的貼片范圍和貼片數(shù)量，通過應(yīng)力分布實驗數(shù)據(jù)進(jìn)一步校準(zhǔn)挖掘機有限元仿真模型，提高仿真計算的精度。