基于有限元的某多管火箭炮模態(tài)分析＊

徐立黃，張貴林，張波濤，謝家松

（1南京理工大學機械工程學院，南京 210094；2湖北江山重工有限責任公司，湖北襄樊 441003）

0 引言

火箭炮是現(xiàn)代戰(zhàn)爭中重要的壓制性武器，但受射擊精度影響，限制了其應用范圍。影響射擊精度的主要因素之一是發(fā)射系統(tǒng)的振動，設法減少火箭炮發(fā)射瞬間的振動，將有效的減小火箭彈起始段的擾動，提高火箭炮的射擊密集度。要減小發(fā)射系統(tǒng)振動對射擊精度的影響，甚至利用發(fā)射系統(tǒng)的振動補償其他擾動因素造成的不利影響，首先必須了解和掌握其固有振動特性，以利于火箭炮結(jié)構(gòu)的優(yōu)化設計和改進。通過ANSYS對火箭炮進行模態(tài)分析，得到結(jié)構(gòu)的固有頻率和振型，能為某多管火箭炮的振動特性分析和結(jié)構(gòu)動力學特性的優(yōu)化設計提供依據(jù)。

1 模態(tài)分析的基本理論

模態(tài)分析理論吸取了振動理論、信號分析、數(shù)據(jù)處理、數(shù)理統(tǒng)計及自動控制理論中的有關(guān)“營養(yǎng)”，結(jié)合自身內(nèi)容的發(fā)展，形成了一套獨特的理論。自動控制中的傳遞函數(shù)概念的引入，對模態(tài)分析理論的發(fā)展起著很大的推動作用。

模態(tài)分析的核心內(nèi)容是確定描述系統(tǒng)動態(tài)特性的參數(shù)。對于一個有N個自由度的線性系統(tǒng)，其運動微分方程為：

對于自由振動（F（t）＝0），上式可寫為：

若忽略阻尼影響，式（2）可進一步簡化為：

式中：M為質(zhì)量矩陣；C為阻尼矩陣；K為剛度矩陣；X為位移向量；F（t）為作用力向量；t為時間。

自由振動時，結(jié)構(gòu)上各點作簡諧振動，各節(jié)點位移為：

由式（3）、式（4）得：

2 建立有限元模型

由于火箭炮結(jié)構(gòu)復雜，在ANSYS中直接建模難度較大，文中選用Pr o／E建立實體模型，然后直接導入ANSYS。合理簡化后的實體模型見圖1，有限元模型見圖2。

圖1 物理模型

圖2 有限元模型

建立物理模型時，需對火箭炮進行合理簡化，其原則是簡化不能影響結(jié)構(gòu)的實際動力學特性。遵循這一原則，忽略結(jié)構(gòu)中的一些圓角、倒角、螺紋孔等。考慮射擊時回轉(zhuǎn)與俯仰方向均被鎖定，因此將箱體與回轉(zhuǎn)體之間用于調(diào)整俯仰角的扇形齒輪簡化為剛性連接，耳軸與回轉(zhuǎn)體之間的軸承連接也簡化為剛性連接。

整個模型采用由20個節(jié)點定義的SOLID95單元，它能夠吸收不規(guī)則形狀的單元而精度沒有損失。SOLID95單元有可并立的位移形狀并且對于曲線邊界的模型能很好的適應，因此非常適合于復雜模型和由CAD導入的模型。模態(tài)分析對網(wǎng)格質(zhì)量的要求與靜力學分析并無大的區(qū)別，但可以忽略應力集中處對網(wǎng)格的細化處理，更加注重網(wǎng)格的均勻性，這是因為固有頻率和振型主要取決于結(jié)構(gòu)質(zhì)量分布和剛度分布，不存在類似應力集中的現(xiàn)象，采用均勻網(wǎng)格可使結(jié)構(gòu)剛度矩陣和質(zhì)量矩陣的元素不致相差太大，可減小數(shù)值計算誤差。

定向管為薄壁長圓管，若采用自由網(wǎng)格劃分則產(chǎn)生的單元數(shù)太多，計算耗時，甚至導致一般計算機無法進行計算。為減少單元和節(jié)點，定向管采用掃掠法劃分網(wǎng)格，為保證與支撐板部分網(wǎng)格正確連接，將定向管沿支撐板面切割后再劃分網(wǎng)格。其余部分結(jié)構(gòu)不規(guī)則，均采用自由網(wǎng)格劃分。對20節(jié)點單元退化為四面體單元的，將其轉(zhuǎn)化為10節(jié)點的SOLID92單元。通過劃分網(wǎng)格，共產(chǎn)生45460個單元，166075個節(jié)點。

整炮均為鋼材，其密度為7830kg／m3，彈性模量為209GPa，泊松比為0.3。模擬發(fā)射角為45°。這里需要強調(diào)的是ANSYS中單位制（Units）問題，ANSYS中有一個只能從命令行輸入的命令：“／UNIT，SI”，它的作用僅是標記作用，讓用戶有個地方做標記，它沒有任何單位轉(zhuǎn)換的功能。這要求單位由用戶去統(tǒng)一，這類似于振動實驗中的標定問題。如果在Pr o／E中建模的長度單位為 mm，這與密度單位kg／m3不符合。若在材料屬性中直接輸入7830，必將導致模型質(zhì)量太大，模態(tài)固有頻率太小。有人提出對密度單位進行換算。由牛頓第二定律F＝Ma可知，1 N＝1kg＊m／s2，那么質(zhì)量單位應該用噸，相應的密度單位就應該是噸每立方毫米了。但問題并未因此得到完美的解決，因為這時的彈性模量也需要調(diào)整，否則結(jié)果依舊會失常。文中采用在Pr o／E中建立模型的時候把單位改成國際單位制，同時，在裝配時調(diào)入模型之前先把裝配的模板的單位改成國際單位制，這樣，調(diào)入ANSYS中的單位就沒有問題了。

3 模態(tài)計算及分析

有限元模型建好后，進入ANSYS求解器，選擇分析類型為 Modal，在求解選項（GUI－Analysis Options）中指定模態(tài)提取方法為Block Lanzcos。這種模態(tài)提取方法采用Lanzcos算法，用一組向量來實現(xiàn)Lanzcos遞歸計算，同Subspace法相比，計算精度相當，但速度要快，存儲要求較低，只是所需內(nèi)存比Subspace法多50%。無論用EQSLV命令指定過哪一種求解器進行求解，Block Lanzcos法都將自動采用稀疏矩陣方程求解器。當計算某系統(tǒng)特征值譜所包含的一定范圍的固有頻率時，采用這種方法特別有效。計算時，求解從頻率譜中間位置到高頻范圍內(nèi)的固有頻率時的求解收斂速度和求解低階頻率時基本相同，因此當采用頻移頻率FREQB來提取從FREQB起始的n階模態(tài)時，該法提取大于FREQB的n階模態(tài)和提取n階低階模態(tài)的速度基本相同。同Subspace法一樣，Block Lanzcos法特別適用于大型對稱特征值求解問題。點擊確定按鈕，在隨后彈出的對話框中設置提取模態(tài)數(shù)No.of modes to extract為12。為了在后處理器中查看振型，必須進行模態(tài)擴展（也就是將振型寫入結(jié)果文件）。將Expansion Pass On／Off設置為ON。在number of Modes to Expand中指定要擴展的模態(tài)數(shù)。這里將其設置為12（與提取模態(tài)數(shù)相同），設置完成后，通過 Main Menu－Sol ution－Current LS對所建立的有限元模型進行計算。

通過計算，得到了模型的前12階模態(tài)，由于從第6階模態(tài)起出現(xiàn)局部模態(tài)，且為集中在定向器上的高階模態(tài)，因此這里只對前6階主要模態(tài)進行分析，詳見表1和圖3～圖8，圖中位移僅表示振型，非實際位移尺寸。

表1 固有頻率和振型

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

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

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

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

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

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

4 結(jié)論

1）該模態(tài)分析結(jié)果與同類武器（參考文獻［4］）相符合，證明計算結(jié)果的正確性。因此可作為諧譜響應分析、頻譜分析等其他動力學分析的基礎。可以為該火箭炮結(jié)構(gòu)動力特性的優(yōu)化設計提供理論依據(jù)，從而縮短研發(fā)周期，降低研發(fā)費用。

2）火箭炮發(fā)射時間間隔一般約為0.5s，即頻率為2 Hz左右，該發(fā)射系統(tǒng)主振頻率較好的避開了2 Hz及其鄰近倍頻。但考慮偏航和俯仰振動對射擊精度的敏感性，建議加強搖架的扭轉(zhuǎn)剛度和箱體俯仰連接剛度。如可將箱體底部加厚鋼板改為鋼架結(jié)構(gòu)，從而減輕質(zhì)量，提高固有頻率。

［1］ 傅志方.振動模態(tài)分析與參數(shù)辨識［M］.北京：機械工業(yè)出版社，1990.

［2］ 劉相新 孟憲頤.ANSYS基礎與應用教程［M］.北京：科學出版社，2005.

［3］ 張琪，劉莉.導彈固有特性的有限元分析［J］.彈箭與制導學報，2008，28（2）：6－7.

［4］ 潘宏俠，王福明.某火箭炮模態(tài)試驗分析及動態(tài)修改［J］.兵工學報武器分冊.1993（3）：58－64.