張麗娜，華順剛

（大連理工大學(xué) 機(jī)械工程學(xué)院，遼寧 大連 116024）

目前計(jì)算機(jī)工程分析技術(shù)，在機(jī)械、造船、航空航天、汽車等各個領(lǐng)域得到了廣泛應(yīng)用。在機(jī)械工程領(lǐng)域中，為了保證機(jī)械系統(tǒng)能夠高效安全地運(yùn)行，結(jié)構(gòu)零部件必須滿足一定的剛度強(qiáng)度要求，因此對機(jī)械系統(tǒng)主要零部件進(jìn)行動力學(xué)分析計(jì)算，可以保證和提高構(gòu)件的剛度強(qiáng)度及穩(wěn)定性。

目前應(yīng)用CAE技術(shù)進(jìn)行機(jī)械系統(tǒng)剛度強(qiáng)度分析的常用方法，包括有限元法和多體動力學(xué)仿真分析方法等。有限元方法，將結(jié)構(gòu)離散化生成有限元網(wǎng)格，通過設(shè)置單元屬性及邊界條件，建立合理的結(jié)構(gòu)數(shù)學(xué)模型進(jìn)行結(jié)構(gòu)動力學(xué)求解，得到結(jié)構(gòu)的固有特性、動態(tài)應(yīng)力與應(yīng)變結(jié)果等，可為結(jié)構(gòu)的設(shè)計(jì)優(yōu)化提供依據(jù)[1～2]；多體動力學(xué)仿真分析，考慮各部件之間的相互作用，以及零件的變形與彈性[3]，建立整機(jī)剛?cè)峄旌隙囿w動力學(xué)模型進(jìn)行系統(tǒng)動力學(xué)分析，可得到各構(gòu)件的作用載荷，全面地描述系統(tǒng)工作的全過程，預(yù)測整個系統(tǒng)的運(yùn)動與響應(yīng)情況[4～5]。

本文以火炮發(fā)射工況為例，對車體進(jìn)行了有限元模態(tài)疊加法瞬態(tài)動力學(xué)分析以及剛?cè)峄旌隙囿w動力學(xué)仿真分析。

1 機(jī)械零部件有限元瞬態(tài)動力學(xué)分析和多體動力學(xué)仿真

1.1 有限元瞬態(tài)動力學(xué)分析

有限元瞬態(tài)動力學(xué)分析，用于確定在隨時間變化載荷作用下結(jié)構(gòu)的動力學(xué)響應(yīng)，利用該方法可以計(jì)算獲得在穩(wěn)態(tài)載荷、瞬態(tài)載荷和簡諧載荷隨意組合作用下結(jié)構(gòu)的應(yīng)力、應(yīng)變等隨時間的變化。

瞬態(tài)動力學(xué)的基本運(yùn)動方程是：

其中，

[M]為質(zhì)量矩陣；

[C]為阻尼矩陣；

[K]為剛度矩陣；

有限元瞬態(tài)動力學(xué)分析有3種求解方法：完全法、縮減法以及模態(tài)疊加法[6]。

完全法采用完整的系數(shù)矩陣計(jì)算瞬態(tài)響應(yīng)，是3種方法中功能最強(qiáng)的，允許包括各類非線性特性，但其也是3種方法中開銷最大的；

縮減法采用主自由度及縮減矩陣來壓縮問題規(guī)模，該方法比完全法快且開銷小，但在整個瞬態(tài)分析過程中，時間步長必須保持恒定，不允許用自動時間步長，且唯一允許的非線性是簡單的點(diǎn)-點(diǎn)接觸；

模態(tài)疊加法利用由模態(tài)分析得到的振型（特征值）乘上因子并求和來計(jì)算結(jié)構(gòu)的響應(yīng)，對于很多問題模態(tài)疊加法比縮減法或完全法更快、開銷更小，但是不能施加非零強(qiáng)制位移，且不能采用自動時間步長。其中模態(tài)疊加法瞬態(tài)動力學(xué)分析的流程如圖1所示。

圖1 模態(tài)疊加法瞬態(tài)動力學(xué)分析流程圖

1.2 多體動力學(xué)仿真分析

通常機(jī)械系統(tǒng)是一個復(fù)雜的多構(gòu)件耦合的機(jī)構(gòu)，各部件之間聯(lián)系緊密。為了全面分析系統(tǒng)的總體性能，得到更準(zhǔn)確的分析結(jié)果，可以將各零部件組裝，并考慮各個部件之間的相互作用，以及零部件的彈性建立剛?cè)峄旌系奶摂M樣機(jī)模型，針對其在實(shí)際應(yīng)用中的各種工況，進(jìn)行動力學(xué)仿真分析，預(yù)測產(chǎn)品的整體性能，進(jìn)而改進(jìn)產(chǎn)品設(shè)計(jì)、提高產(chǎn)品性能。

對機(jī)械系統(tǒng)進(jìn)行剛?cè)峄旌咸摂M樣機(jī)多體動力學(xué)仿真分析的流程如圖2所示。

圖2 多體動力學(xué)仿真分析流程圖

2 火炮發(fā)射工況車體有限元瞬態(tài)動力學(xué)分析

車載火炮武器發(fā)射是一個短暫而又復(fù)雜的過程，瞬間產(chǎn)生很大的沖擊載荷[7]。車體是車輛的重要組成部分，其結(jié)構(gòu)性能直接影響到整車的安全性和機(jī)動性，具有防護(hù)、安裝機(jī)件、承受沖擊載荷等重要功能[8]。

在火炮發(fā)射過程中，車身作為一個受力部件，在強(qiáng)度和剛度方面必須保證達(dá)到一定的要求，在保證其工作壽命和使用要求的同時，還需要保證乘員的安全，因此有必要對火炮發(fā)射過程中車體的剛度強(qiáng)度進(jìn)行分析。

本文采用有限元模態(tài)疊加法，對火炮在不同方向射角α 及高低射角β 發(fā)射工況下的車體，進(jìn)行了瞬態(tài)動力學(xué)分析，步驟如下：

（1）在三維CAD 軟件中創(chuàng)建車輛零部件幾何模型；

（2）抽取車體中面模型導(dǎo)入有限元軟件中，在相應(yīng)位置處添加支撐梁單元得到車體有限元模型；

（3）對車體進(jìn)行模態(tài)求解，求解過程中我們采用默認(rèn)的Block Lanczos 方法，并在懸掛系統(tǒng)處添加位移約束，重力加速度及所需的單元載荷；

（4）利用模態(tài)求解結(jié)果，并施加后坐力載荷進(jìn)行模態(tài)疊加法瞬態(tài)分析，對于隨時間變化的后坐力，采用分步加載方式；

（5）擴(kuò)展模態(tài)疊加解，得到結(jié)構(gòu)分析結(jié)果文件；

（6）在后處理器中觀察結(jié)果，得到各載荷步下車體的應(yīng)力、應(yīng)變結(jié)果。

為了接近實(shí)際，提高有限元分析計(jì)算的準(zhǔn)確性，在邊界條件處理中，我們采用彈簧單元combin14 并賦予相應(yīng)的彈性和阻尼系數(shù)來模擬懸掛系統(tǒng)。

對于車體所承受的配重力，我們采用了兩種處理方式：對于配重比較大的載荷（如發(fā)動機(jī)、減速器配重等）將其簡化為面壓力；對于配重較小的載荷（如隨車工具、前蓋等），我們建立剛性區(qū)域，將配重分散在幾個節(jié)點(diǎn)上加載。

對于瞬態(tài)分析中所需的面載荷及重力加速度，我們在模態(tài)分析中施加。程序會計(jì)算出一個載荷向量并將其寫入振型文件，這些載荷會在模態(tài)分析中忽略，在瞬態(tài)分析中調(diào)用載荷向量即可。

我們對α=0°、45°、90°，β=0°、60°的各組合情況進(jìn)行了分析。由有限元分析結(jié)果我們可以得到：位于發(fā)動機(jī)、減速器下方的底面區(qū)域，在火炮發(fā)射的多種射角下，經(jīng)常是整個車體的最大應(yīng)力區(qū)域。在后坐力方向改變時，位于炮塔下方的車體底面位置，成為了某一時刻的車體最大應(yīng)力區(qū)域，但應(yīng)力值均在許用范圍內(nèi)。當(dāng)α=90°，β=60°，火炮后坐力達(dá)到最大時，車體底面某局部，成為了此時的車體最大應(yīng)力區(qū)域。

3 火炮發(fā)射工況車輛多體動力學(xué)仿真分析

多體動力學(xué)通過建立并分析機(jī)械系統(tǒng)的數(shù)學(xué)模型，得到系統(tǒng)構(gòu)件的運(yùn)動學(xué)規(guī)律和動力學(xué)響應(yīng)。考慮到車輛的高度耦合性以及車體彈性對整車行駛性能的影響，將車體考慮為柔性體，建立車輛的剛?cè)峄旌咸摂M樣機(jī)模型，進(jìn)行多體動力學(xué)仿真- 分析。

RecurDyn 采用相對坐標(biāo)系統(tǒng)和完全遞歸算法，非常適合于求解大規(guī)模及復(fù)雜接觸的多體系統(tǒng)動力學(xué)問題。該軟件包含多個專用模塊，我們利用其中的針對高速履帶車輛開發(fā)的Track-HM模塊，構(gòu)建履帶子系統(tǒng)模型，實(shí)現(xiàn)高速履帶車輛的運(yùn)動學(xué)、動力學(xué)仿真，并通過基于模態(tài)縮減技術(shù)的R-Flex 剛?cè)峄旌辖＜夹g(shù)，建立了車輛系統(tǒng)中的柔性車體[9]。

本文車輛剛?cè)峄旌咸摂M樣機(jī)模型建立過程如下：

（1）利用三維CAD 軟件建立車輛零部件幾何模型；

（2）對車體三維模型抽取中面，將中面模型導(dǎo)入有限元分析軟件中建立有限元模型；

（3）在相應(yīng)位置處添加支撐梁單元，對車體劃分網(wǎng)格后得到車體的有限元模型；

（4）對車體有限元模型進(jìn)行模態(tài)分析以得到柔性車體RFI 文件；

（5）將創(chuàng)建的RFI 文件導(dǎo)入RecurDyn 中；

（6）添加炮塔部分、懸掛系統(tǒng)、行走系統(tǒng)、配重及其他零部件；

（7）建立各零部件之間的連接關(guān)系，得到車輛的剛?cè)峄旌咸摂M樣機(jī)模型。

此處的模態(tài)求解雖然也在Ansys 中進(jìn)行，但與之前的模態(tài)求解設(shè)置不同。此處設(shè)置的模態(tài)求解數(shù)NMODES為無約束模態(tài)階數(shù)，約束模態(tài)階數(shù)由定義的接口節(jié)點(diǎn)數(shù)決定。求解過程中的約束條件通過指定接口節(jié)點(diǎn)并運(yùn)行Ansys 與RecurDyn的接口MAC文件來實(shí)現(xiàn)。

由于火炮發(fā)射過程中的后坐力變化比較復(fù)雜，用RecurDyn 中的函數(shù)表達(dá)式不容易實(shí)現(xiàn)，因此使用用戶子程序，利用高級匯編語言C 來定義。首先編寫用戶自定義子程序代碼，之后生成動態(tài)鏈接庫，在Recur-Dyn模型中對已創(chuàng)建的動態(tài)鏈接庫DLL 進(jìn)行調(diào)用。

我們對α=0°、45°、90°，β=0°、60°的不同情況，進(jìn)行了動力學(xué)仿真分析，得到了不同射角情況下的部件受力以及柔性車體的應(yīng)力及應(yīng)變結(jié)果。對結(jié)果進(jìn)行分析，我們得到車體底面經(jīng)常是車體最大應(yīng)力區(qū)域，后坐力方向改變時，車體最大應(yīng)力區(qū)域也會相應(yīng)的改變；當(dāng)高低射角β =60°時，位于炮塔下方的底面區(qū)域的應(yīng)力明顯增大，但均在許用應(yīng)力范圍內(nèi)；當(dāng)方向射角α 改變時，左右兩側(cè)扭桿所受扭矩隨方向射角的改變而改變，車體恢復(fù)穩(wěn)定后，左右側(cè)扭桿扭矩恢復(fù)相同。

比較車體有限元瞬態(tài)動力學(xué)和車輛剛?cè)峄旌咸摂M樣機(jī)多體動力學(xué)兩種方法的結(jié)果，可以得到兩種方法下，火炮發(fā)射時車體的應(yīng)力分布相似，但有限元瞬態(tài)動力學(xué)得到的應(yīng)力值比多體動力學(xué)仿真分析的結(jié)果要略大。由于虛擬樣機(jī)技術(shù)考慮了部件之間的相互聯(lián)系，懸掛系統(tǒng)的緩沖作用等，對車體受力分析的邊界條件更接近實(shí)際，因此結(jié)果更合理。

4 結(jié)束語

本文對受沖擊載荷的車體，分別進(jìn)行了有限元模態(tài)疊加法瞬態(tài)動力學(xué)分析和剛?cè)峄旌隙囿w動力學(xué)分析，并對其分析的一般方法和步驟進(jìn)行了總結(jié)。并以火炮發(fā)射過程為例，利用兩種方法對車體進(jìn)行了動力學(xué)分析，兩種分析方法得到的結(jié)果相似。這兩種方法在產(chǎn)品設(shè)計(jì)過程中，均被廣泛應(yīng)用，在分析系統(tǒng)結(jié)構(gòu)的固有特性及主要部件的力學(xué)性能，深入了解系統(tǒng)內(nèi)部結(jié)構(gòu)的動力學(xué)特性時，可以采用有限元方法；在需要研究系統(tǒng)的總體性能，預(yù)測整個工況系統(tǒng)的動態(tài)響應(yīng)情況時，采用多體動力學(xué)分析方法。本文的研究為其他機(jī)械零部件剛強(qiáng)度分析提供一定的參考。

[1]張衛(wèi)國，權(quán) 龍，程 珩，楊 敬：基于真實(shí)載荷的挖掘機(jī)工作裝置瞬態(tài)動力學(xué)分析[J].機(jī)械工程學(xué)報，2011，47（12）：144-149.

[2]呂彭民，丁 智.架橋機(jī)主梁結(jié)構(gòu)瞬態(tài)動力學(xué)分析[J].長安大學(xué)學(xué)報(自然科學(xué)版)，2005，25（4）：80-83.

[3]林建生，王 珊，張寶歡，張 霖.內(nèi)燃機(jī)多連桿機(jī)構(gòu)的多體動力學(xué)分析[J].天津大學(xué)學(xué)報，2007，40（6）：640-643.

[4]王 玉，羊 玢：基于RecurDyn的液壓挖掘機(jī)的建模與動態(tài)分析[J].重慶理工大學(xué)學(xué)報(自然科學(xué))，2011，25（8）：10-13.

[5]聶旭濤，郭隆德，劉伯林.基于ADAMS 風(fēng)洞柔壁噴管動力學(xué)仿真分析[J].實(shí)驗(yàn)流體力學(xué)，2011，25（2）：73-76.

[6]龔曙光.ANSYS 工程應(yīng)用實(shí)例解析[M].北京：機(jī)械工業(yè)出版社，2003.

[7]馬吉勝.自行火炮發(fā)射動力學(xué)仿真[J].計(jì)算機(jī)仿真，2003，20（6）：8-10.

[8]閆清東,張連第,趙毓芹.坦克構(gòu)造與設(shè)計(jì)[M].北京：北京理工大學(xué)出版社，2007.

[9]焦曉娟，張湝渭，彭斌彬.RecurDyn 多體系統(tǒng)優(yōu)化仿真技術(shù)[M].北京：清華大學(xué)出版社，2010.