坦克火炮身管模態(tài)分析及試驗(yàn)驗(yàn)證

郝丙飛， 王紅巖， 芮 強(qiáng)

(裝甲兵工程學(xué)院 機(jī)械工程系, 北京 100072)

【裝備理論與裝備技術(shù)】

坦克火炮身管模態(tài)分析及試驗(yàn)驗(yàn)證

郝丙飛， 王紅巖， 芮 強(qiáng)

(裝甲兵工程學(xué)院 機(jī)械工程系, 北京 100072)

針對(duì)火炮身管振動(dòng)特性對(duì)火炮高低向穩(wěn)定精度和射擊密集度造成的影響，以及身管彈性變形與模態(tài)的內(nèi)在關(guān)系，以某型主戰(zhàn)坦克火炮身管為例，建立身管結(jié)構(gòu)有限元模型；進(jìn)行模態(tài)分析，得到約束狀態(tài)下的固有頻率和振型；基于LMS測(cè)試系統(tǒng)對(duì)火炮身管進(jìn)行模態(tài)試驗(yàn)分析；試驗(yàn)與仿真計(jì)算結(jié)果對(duì)比，結(jié)果表明：所建身管有限元結(jié)構(gòu)模型能夠較好的符合身管實(shí)際的振動(dòng)特性，為火炮振動(dòng)響應(yīng)特性分析和模態(tài)優(yōu)化提供參考。

兵器科學(xué)與技術(shù)；火炮身管；有限元模態(tài)分析；模態(tài)試驗(yàn)

火炮身管是坦克炮控系統(tǒng)重要組成部件，其動(dòng)力學(xué)特性對(duì)炮控系統(tǒng)穩(wěn)定精度、射擊密集度、身管壽命等都有重要影響，模態(tài)分析和模態(tài)試驗(yàn)是研究身管振動(dòng)特性的重要方法之一。火炮振動(dòng)是其結(jié)構(gòu)動(dòng)態(tài)特性在外界激勵(lì)下的響應(yīng)，近年來隨著對(duì)火炮威力要求的不斷提高，往往通過改變身管的長(zhǎng)徑比提高火炮的初速[1-2]，在這種情況下身管的動(dòng)力學(xué)特性必然發(fā)生變化。為保證火炮具有良好的動(dòng)力學(xué)性能，在設(shè)計(jì)之初就應(yīng)考慮身管的模態(tài)。

針對(duì)火炮身管結(jié)構(gòu)及其振動(dòng)特性對(duì)火炮射擊精度的影響，吳東亞等[2-3]對(duì)加長(zhǎng)身管的火炮進(jìn)行了有限元仿真分析，為身管的優(yōu)化和改進(jìn)提出了參考方案；鄧劍等[4]建立了身管振動(dòng)動(dòng)力學(xué)模型，并提出通過改變模型中搖架前后襯套的距離改善身管的模態(tài)，為減小火炮的振動(dòng)提供依據(jù)；徐亞棟等[5]建立了復(fù)合材料的身管有限元模型，對(duì)復(fù)合材料身管進(jìn)行結(jié)構(gòu)優(yōu)化，通過提高材料強(qiáng)度和減輕重量的方式改變身管的模態(tài)。豆征等[6]采用有限元方法對(duì)身管模態(tài)進(jìn)行研究，對(duì)各階模態(tài)振型對(duì)火炮射擊精度和密集度的影響進(jìn)行分析，提出通過改變彈性模量與頻率的關(guān)系，提高火炮的射擊穩(wěn)定性；Dr.Stenphen Wilkerson等[7-8]分別通過試驗(yàn)和有限元仿真的方法對(duì)水平方向和垂直方向的M1A1坦克火炮身管進(jìn)行約束模態(tài)分析，對(duì)身管的非線性模型進(jìn)行了研究。

綜上所述，不少文獻(xiàn)都對(duì)火炮身管進(jìn)行了有限元模態(tài)分析，本文通過實(shí)車試驗(yàn)對(duì)有限元計(jì)算的模態(tài)結(jié)果進(jìn)行了驗(yàn)證，并進(jìn)行振動(dòng)特性分析。

1 彈性火炮身管的運(yùn)動(dòng)描述

坦克在越野路面上行駛時(shí)，底盤受地面輸入的激擾作用產(chǎn)生劇烈振動(dòng)，對(duì)安裝在底盤上的炮塔和火炮產(chǎn)生激勵(lì)作用，同時(shí)火炮身管作為安裝在炮塔上的懸臂梁結(jié)構(gòu)，在激勵(lì)的作用下產(chǎn)生彎曲振動(dòng)?；鹋谂诳诘恼駝?dòng)情況可以表達(dá)如下：

(1)

式中，y為高低向火炮的垂直振動(dòng)位移量(在絕對(duì)坐標(biāo)系下)，x為響應(yīng)測(cè)量點(diǎn)距離炮尾端面的水平距離，t為響應(yīng)時(shí)間，yt為耳軸中心高低向振動(dòng)位移，θ為火炮繞耳軸中心的俯仰角位移，xt為耳軸中心點(diǎn)距離炮尾端面的水平距離。Y為身管高低向由于彈性變形而產(chǎn)生的振動(dòng)位移。Y可以進(jìn)一步表達(dá)如下：

Y(x,t)=q1(t)Y1(x)+q2(t)Y2(x)+…+qn(t)Yn(x)

(2)

式中Yi(i=1,2,…,n)為前n階振型，qi為第n階振型的振動(dòng)參與因子。

火炮身管通過耳軸安裝在炮塔上近似為外伸的懸臂梁結(jié)構(gòu)，如圖1所示，主要產(chǎn)生橫向、縱向和徑向振動(dòng)。試驗(yàn)結(jié)果表明，縱向和徑向振動(dòng)對(duì)火炮的穩(wěn)定精度影響很小，可忽略不計(jì)。橫向振動(dòng)對(duì)火炮的穩(wěn)定精度影響較大。為列出橫向振動(dòng)的微分方程，在身管上選取微元dx受力分析，從梁的靜撓曲變形，得知：

(3)

式中：E為身管的彈性模量，J為微元繞中性層的慣性矩，M表示距O點(diǎn)為x處截面的彎矩。對(duì)式(3)進(jìn)行微分得到：

(4)

(5)

這里，Q為剪力，γ為材料密度，g為重力加速度，將式(4)和式(5)整理得[9]：

(6)

式中c2為EJg/(Aγ)。

圖1 身管懸臂梁簡(jiǎn)化模型

身管振動(dòng)的固有頻率與身管的材料、剛度以及質(zhì)量分布有關(guān)，振幅與初始條件、固有頻率有關(guān)。整個(gè)橫向振動(dòng)的形成是由基階振動(dòng)和幾個(gè)高階振動(dòng)組合而成?；A振動(dòng)頻率低、振幅大，對(duì)火炮的射擊精度影響較大，高階振動(dòng)頻率高、振幅小，對(duì)武器的性能產(chǎn)生的影響較小。

2 身管有限元結(jié)構(gòu)模態(tài)分析

首先，基于三維建模平臺(tái)Solidworks建立某型坦克火炮身管的三維數(shù)字模型，如圖2所示。

圖2 身管三維數(shù)字模型

模型的幾何參數(shù)、質(zhì)量屬性、材料屬性均按照工程設(shè)計(jì)圖紙和相關(guān)技術(shù)文件獲得，其主要參數(shù)如表1所示。

表1 身管模型結(jié)構(gòu)主要參數(shù)

2.1 建立火炮身管有限元模型

為了減小有限元網(wǎng)格的劃分難度和提高計(jì)算效率，在保證身管結(jié)構(gòu)主要?jiǎng)恿W(xué)特性的原則下，根據(jù)身管的結(jié)構(gòu)特點(diǎn)，參考有關(guān)模態(tài)計(jì)算經(jīng)驗(yàn)[10]，對(duì)身管的倒角、螺紋、局部細(xì)節(jié)特征等進(jìn)行簡(jiǎn)化。繪制有限元網(wǎng)格是有限元計(jì)算的基礎(chǔ)，網(wǎng)格質(zhì)量的好壞直接影響計(jì)算結(jié)果的準(zhǔn)確性和計(jì)算時(shí)間。由于結(jié)構(gòu)的固有頻率和振型主要取決于質(zhì)量和剛度的分布，所以模態(tài)分析時(shí)取較為均勻的劃分方式，本文以全六面體對(duì)火炮身管進(jìn)行網(wǎng)格劃分，網(wǎng)格大小為10 mm，共有節(jié)點(diǎn)162 552個(gè),單元總數(shù)為129 896，火炮身管有限元模型如圖3所示，在身管炮尾末端和耳軸處添加相應(yīng)約束。

在理想情況下，通過模態(tài)計(jì)算可以得到一個(gè)結(jié)構(gòu)完整的系統(tǒng)模態(tài)集[11-12]，但是在實(shí)際應(yīng)用中根本不必要。由于低階頻率對(duì)系統(tǒng)響應(yīng)的貢獻(xiàn)較大，所以實(shí)際應(yīng)用中采用模態(tài)截?cái)嗵幚?，將高階模態(tài)舍棄。

圖3 身管有限元網(wǎng)格劃分模型

2.2 身管有限元結(jié)構(gòu)模態(tài)分析

在HyperWorks環(huán)境中應(yīng)用Radioss求解器對(duì)約束狀態(tài)下的火炮身管模型進(jìn)行有限元模態(tài)求解，求解得到各階模態(tài)的頻率如表2所示，前三階振型如圖4所示。

表2 身管前5階固有頻率和振動(dòng)周期

圖4 身管高低向前三階模態(tài)振型云圖

通過圖4觀察，前三階模態(tài)振動(dòng)位移最大處均在炮口，這種振動(dòng)將嚴(yán)重影響彈丸射出炮口的初設(shè)角，影響坦克的射擊精度。

為了對(duì)比各階模態(tài)對(duì)炮口振動(dòng)的貢獻(xiàn)量，對(duì)身管進(jìn)行模態(tài)頻率響應(yīng)分析。在身管上選取典型節(jié)點(diǎn)位置，沿高低向施加200 N的掃頻集中力，頻率變化范圍0～1 000 Hz，在該力的作用下，炮口振動(dòng)的頻率響應(yīng)結(jié)果如圖5所示。

圖5 炮口振動(dòng)頻率響應(yīng)分析

由圖5可見，各階模態(tài)固有頻率的模態(tài)參與因子由低頻至高頻依次減小，前三階幅值比值近似為25∶5∶1，三階以上振型幅值可忽略不計(jì)。

各階模態(tài)的振動(dòng)參與因子隨系統(tǒng)輸入和時(shí)間的變化而變化。為描述其一般性的變化規(guī)律，對(duì)約束狀態(tài)下的身管進(jìn)行瞬態(tài)動(dòng)力學(xué)分析。在身管上選取典型節(jié)點(diǎn)，施加單位脈沖激勵(lì)，仿真炮口振動(dòng)響應(yīng)，通過濾波分離前三階振型幅值變化，如圖6所示。圖中幅值最大的是一階模態(tài)，二階模態(tài)幅值較小，三階模態(tài)幅值更小，接近一條直線。前兩階振型幅值即第一階模態(tài)和第二階模態(tài)之比值q1/q2隨時(shí)間的變化規(guī)律如圖7所示。

圖6 前三階模態(tài)振型幅值隨時(shí)間變化

圖7 前兩階模態(tài)振型幅值比隨時(shí)間變化

由圖6和圖7可見，前三階模態(tài)中，一階模態(tài)對(duì)炮口的振動(dòng)貢獻(xiàn)量最大，超過70%，在誤差允許的范圍內(nèi)。而二階模態(tài)和三階模態(tài)的振動(dòng)貢獻(xiàn)量更小，可以忽略不計(jì)。

3 身管模態(tài)試驗(yàn)驗(yàn)證

為了驗(yàn)證有限元結(jié)構(gòu)模態(tài)計(jì)算結(jié)果的準(zhǔn)確性，對(duì)火炮身管進(jìn)行試驗(yàn)驗(yàn)證，驗(yàn)證約束狀態(tài)下火炮身管各階固有頻率和主振型。身管模態(tài)試驗(yàn)裝置和試驗(yàn)過程傳感器布置如圖8所示，包括LMS測(cè)試系統(tǒng)、振動(dòng)加速度傳感器、試驗(yàn)力錘等。

圖8 試驗(yàn)裝置及傳感器布置

試驗(yàn)工況為：坦克炮控系統(tǒng)關(guān)閉，行軍閉鎖器閉鎖，火炮身管處于完全約束狀態(tài)，身管軸線俯仰角設(shè)置在很小的角度(5°左右)，熱護(hù)套去除，身管排氣裝置拆除。分別沿水平方向和高低向兩個(gè)方向?qū)ι砉苁┘渝N擊脈沖力，采集身管不同點(diǎn)、不同方向的振動(dòng)加速度響應(yīng)。LMS測(cè)試系統(tǒng)記錄力錘激振力信號(hào)和多通道加速度信號(hào)。

LMS.test.lab軟件中集成了專門用于模態(tài)試驗(yàn)分析的modal analysis模塊，提供了兩種模態(tài)分析方法：Time MDOF和PolyMAX。Time MDOF是從時(shí)域的角度進(jìn)行參數(shù)識(shí)別，多用于小阻尼結(jié)構(gòu)分析；PolyMAX是從頻域的角度進(jìn)行參數(shù)識(shí)別，適用于大阻尼結(jié)構(gòu)分析。本試驗(yàn)選用PolyMAX法進(jìn)行模態(tài)識(shí)別。

PolyMAX是基于模型穩(wěn)定圖法進(jìn)行模態(tài)參數(shù)識(shí)別。模態(tài)穩(wěn)定圖是通過計(jì)算得到的一組結(jié)構(gòu)模態(tài)中發(fā)現(xiàn)有物理意義的結(jié)構(gòu)模態(tài)。模態(tài)穩(wěn)定圖的橫坐標(biāo)為頻率，縱坐標(biāo)為結(jié)構(gòu)模態(tài)階數(shù)(擬合模態(tài)的個(gè)數(shù))。通過對(duì)不同模態(tài)階數(shù)的計(jì)算顯示在穩(wěn)定圖中?；鹋谏砉苁艿搅﹀N激振產(chǎn)生多點(diǎn)振動(dòng)加速度，通過LMS系統(tǒng)可建立輸入力錘激振與多點(diǎn)振動(dòng)加速度響應(yīng)的對(duì)應(yīng)關(guān)系，從而得到多組頻響函數(shù)，進(jìn)行多次疊加擬合獲得系統(tǒng)頻響函數(shù)。身管約束狀態(tài)的前五階固有頻率和阻尼比如表3所示。前三階振型如圖9所示。

得到身管的模態(tài)參數(shù)后，為了評(píng)價(jià)模態(tài)參數(shù)估計(jì)的質(zhì)量，對(duì)結(jié)構(gòu)的前5階模態(tài)振型進(jìn)行相關(guān)性分析。通過MAC值(模態(tài)置信度值)來評(píng)價(jià)模態(tài)參數(shù)估計(jì)的質(zhì)量，MAC矩陣用于計(jì)算振型的獨(dú)立性，MAC=(macij)，MAC的元素macij定義如下：

(7)

φi和φj為分別為第i階和第j階的計(jì)算模態(tài)振型。其中，對(duì)角線的MAC值越接近100，非對(duì)角線MAC值越接近0，說明各階模態(tài)的相關(guān)性越好，最后以相關(guān)性分析結(jié)果較好的模態(tài)參數(shù)作為模態(tài)試驗(yàn)的分析結(jié)果。前5階振型的MAC矩陣相關(guān)性分析如表4所示。

表3 身管高低向前5階模態(tài)頻率和阻尼比

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

Mode1Mode2Mode3Mode4Mode5Mode110011．8386．6375．5610．537Mode211．83810013．5341．1800．133Mode36．63713．53410056．1153．150Mode45．5611．18056．1151002．660Mode50．5370．1343．1502．660100

通過表4中身管前5階試驗(yàn)?zāi)B(tài)相關(guān)性可知，表的對(duì)角線值均為100，非對(duì)角線的值接近于0，但是，還存在兩個(gè)比較大的值，為56.115，分析原因第3階和第4階振型分別為橫向振動(dòng)和徑向振動(dòng)，二者在數(shù)值上大小相差不大，所以模態(tài)試驗(yàn)中會(huì)產(chǎn)生一定的相干。

通過模態(tài)試驗(yàn)計(jì)算得到火炮身管約束狀態(tài)下的固有頻率和主振型，取前三階有限元仿真計(jì)算模態(tài)頻率與試驗(yàn)結(jié)果對(duì)比如表5所示。

表5 模態(tài)頻率計(jì)算結(jié)果對(duì)比

通過表5分析可得，試驗(yàn)得到的模態(tài)結(jié)果與有限元仿真計(jì)算得到的模態(tài)結(jié)果基本吻合，誤差均在10%以內(nèi)，充分表明所建有限元模型的可信性。

4 結(jié)論

本文分析了坦克火炮身管的運(yùn)動(dòng)特性，建立了身管的有限元模型，仿真計(jì)算火炮身管約束狀態(tài)下的模態(tài)頻率和模態(tài)振型，通過模態(tài)試驗(yàn)對(duì)有限元模態(tài)計(jì)算結(jié)果進(jìn)行了驗(yàn)證，結(jié)論如下：

1) 坦克在越野環(huán)境下機(jī)動(dòng)，火炮身管產(chǎn)生剛性運(yùn)動(dòng)和彈性振動(dòng)，二者相疊加，使火炮炮口的振動(dòng)加劇。

2) 火炮身管模態(tài)試驗(yàn)結(jié)果與有限元模型計(jì)算火炮身管的模態(tài)頻率和模態(tài)振型結(jié)果基本一致，表明身管有限元模型能夠很好地表達(dá)炮口的振動(dòng)特性。

3) 約束狀態(tài)下，身管的高低向低階模態(tài)振型主要是橫向振動(dòng)，橫向彎曲是身管的主要振動(dòng)變形形式。

[1] 譚樂斌,張相炎,管紅根,等.火炮概論[M].北京：北京理工大學(xué)出版社，2003：40-42.

[2] 吳東亞,邢宏光,崔軍，等.基于有限元的某型坦克炮身管模態(tài)分析[J].科技導(dǎo)報(bào),2008,26(23):44-47.

[3] 吳東亞,邢宏光,崔軍.基于虛擬樣機(jī)技術(shù)的模型坦克炮炮身模態(tài)分析[J].火力與指揮控制,2011,36(7):65-71.

[4] 鄧劍,郭保全,韓海濤.基于搖架約束的炮口振動(dòng)仿真分析[J].計(jì)算機(jī)仿真，2012,29(9)：14-16.

[5] 徐亞棟,錢林方,陳龍淼.復(fù)合材料身管三維有限元分析與結(jié)構(gòu)優(yōu)化[J].彈箭與制導(dǎo)學(xué)報(bào),2005,25(4):209-212.

[6] 豆征,陳永剛,李勇.基于FEM法的身管振動(dòng)特性研究[J].機(jī)械設(shè)計(jì)與研究,2013,29(3):126-130.

[7] DR.STEPHEN WILKERSON,MR.MORRIS BORMAN,MR.TING LI.A Model Survey of the M1A1 Main Weapon System[R].U.S Army Research Laboratory.AMSRL_WT_PD,1993.

[8] BIRD J S.Measurement of Tank Gun Dynamics in Support of Dynamic Muzzle Referencing System[R].Ottawa.Defence Research Establishment OTTAWA.1990.

[9] 聞邦椿,劉樹英,陳照波.機(jī)械振動(dòng)理論及應(yīng)用[M].北京：高等教育出版社,2010:72-79.

[10] 陳可,周孔亢,徐凌,等.低速汽車整車模態(tài)試驗(yàn)研究[J].機(jī)械設(shè)計(jì)與制造,2007(12)：31.

[11] 蘇忠亭,徐達(dá),楊明華,等.基于模態(tài)試驗(yàn)的某火炮身管有限元模型修正[J].振動(dòng)與沖擊，2012,31(24):54-59.

[12] 呂佼珂,程廣偉,胡勇.基于workbench的車載高炮炮身模態(tài)分析[J].火炮發(fā)射與控制學(xué)報(bào),2016,37(1):36-39.

ModalAnalysisandExperimentalVerificationofTankGunBasedonFEMMethod

HAO Bingfei, WANG Hongyan, RUI Qiang

(Department of Mechanical Engineering, Academy of Armored Force Engineering, Beijing 100072, China)

Aiming at the vibration characteristics of gun barrel that makes a big influence on the stability precision and shooting intensity in elevation, and as well as the inherent relationship between the elastic deformation of the barrel and modal, taking a certain type of main battle tank as a typical example, finite element model of body structure about modal analysis is built by finite element method with performing of modal analysis. The natural frequencies and modes shape in constrained were calculated. Modal experiment analysis of gun barrel was based on LMS test system. The FEM simulation results were compared with the experiment result. The result show that the finite element model of the barrel can better reflect the vibration characteristics of the barrel, which provides a model reference for the vibration response characteristic analysis and modal optimization of the gun.

ordnance science and technology; gun barrel; finite element modal analysis; modal experiment

2017-06-20；

2017-07-10

郝丙飛(1992—)，男，碩士研究生,主要從事裝甲車輛論證與仿真研究。

10.11809/scbgxb2017.10.016

本文引用格式：郝丙飛,王紅巖,芮強(qiáng).坦克火炮身管模態(tài)分析及試驗(yàn)驗(yàn)證[J].兵器裝備工程學(xué)報(bào)，2017(10):74-78.

formatHAO Bingfei, WANG Hongyan, RUI Qiang.Modal Analysis and Experimental Verification of Tank Gun Based on FEM Method[J].Journal of Ordnance Equipment Engineering,2017(10):74-78.

TJ38

A

2096-2304(2017)10-0074-05

(責(zé)任編輯周江川)