商保利，戰(zhàn)仁軍，商 鵬，張洪彪

(武警工程大學(xué)裝備工程學(xué)院，陜西西安710086)

1 引 言

九管發(fā)射器(圖1)是武警常用的車載非致命防暴武器，可用來發(fā)射38 mm催淚彈、防暴彈、煙霧彈、布袋彈等不同彈體，發(fā)射距離在100 m以上。該武器可以使人員短時間內(nèi)失去抵抗能力，其強(qiáng)大的威力不僅適用于快速制止、驅(qū)散非法集會，打擊流氓團(tuán)伙斗毆等聚眾鬧事活動，更能很好地適用于武裝巡邏、緝毒緝私、防暴搜捕、平息武裝叛亂等任務(wù)。其和單兵使用的防暴槍相比，具有發(fā)射距離遠(yuǎn)、發(fā)射精度高、發(fā)射速度快的特點。其在處置大規(guī)模群體性事件中發(fā)揮著重要作用，是武警部隊維護(hù)社會穩(wěn)定的一種大威力的防暴非致命武器。

發(fā)射器在發(fā)射時受到固有頻率、發(fā)動機(jī)轉(zhuǎn)速、發(fā)射頻率以及不同路況等因素的影響，其發(fā)射管往往會發(fā)生彎曲振動甚至產(chǎn)生共振。而這些振動會直接影響射擊的精度和系統(tǒng)的可靠性。因此，有必要對發(fā)射器進(jìn)行振動模態(tài)分析，獲取發(fā)射器的固有頻率，盡量避免共振的發(fā)生，為發(fā)射管的優(yōu)化設(shè)計和改進(jìn)提供合理的依據(jù)。

有限元方法是求解復(fù)雜物理數(shù)學(xué)問題的重要方法，能夠?qū)?fù)雜的結(jié)構(gòu)進(jìn)行數(shù)值模擬。有限元方法(FEM)的基礎(chǔ)是變分原理或加權(quán)余量法，其基本求解思想是把計算域劃分為有限個互不重疊的單元，在每個單元內(nèi)，選擇一些合適的節(jié)點作為求解函數(shù)的插值點，將微分方程中的變量改寫成由各變量或其導(dǎo)數(shù)的節(jié)點值與所選用的插值函數(shù)組成的線性表達(dá)式，借助于變分原理或加權(quán)余量法，將微分方程離散求解。有限元方法最早應(yīng)用于結(jié)構(gòu)力學(xué)，后來隨著計算機(jī)的發(fā)展慢慢用于流體力學(xué)的數(shù)值模擬。有限元方法已經(jīng)廣泛應(yīng)用于軍事裝備的建模仿真。其中吳東亞對某型坦克炮身管進(jìn)行了模態(tài)分析［1］。張海航進(jìn)行了多管火炮的模態(tài)分析［2］。但由于之前部隊對非致命武器射擊精確度要求較低，大家忽略了對非致命武器動力學(xué)的研究。這就在一定程度上造成了彈藥的浪費。近幾年，隨著暴騷亂等群體性事件的增多，非致命武器的適用范圍也越來越廣。這就要求必須提高非致命武器的射擊精度。既能實現(xiàn)有效打擊又可以避免不必要的浪費。

所以本文以武警常用的車載九管發(fā)射器為研究對象，以有限元方法為理論支撐，以CAD軟件Pro/E和大型常用的有限元軟件ANSYS 13．0為工具，對發(fā)射器的振動模態(tài)進(jìn)行分析，獲取發(fā)射管的固有頻率，為進(jìn)一步進(jìn)行動力學(xué)研究奠定基礎(chǔ)。

2 振動對發(fā)射穩(wěn)定性的影響

在武器發(fā)射的過程中，發(fā)射管的振動情況會增大彈著點的散布。在發(fā)射瞬間由于受到后坐力的影響發(fā)射管會產(chǎn)生俯仰運動。這樣就會影響下一發(fā)彈的起始擾動，進(jìn)而影響射擊精度。同時，如果車輛在較為復(fù)雜路面行駛，會引起車體的振動。若車體的振動頻率和發(fā)射管的固有頻率相同或相近，就極易引起共振。發(fā)射管的振動階數(shù)越低，對發(fā)射器射擊精度的影響就會越大。

2.1 靜止發(fā)射時的振動激勵

靜止發(fā)射時，振動的激勵來自于彈體發(fā)射的后坐力。在各方向分力的作用下，發(fā)射管產(chǎn)生不同方向的振動。若采用連發(fā)模式，第一發(fā)彈發(fā)射帶來的振動會對第二發(fā)彈的初始擾動帶來影響。第二發(fā)的發(fā)射時間應(yīng)避開第一發(fā)彈引起振動的最大振幅或者發(fā)射時間相差半個振動周期。

若采用多管齊發(fā)模式，管數(shù)不同則后坐力的大小不同，總的后坐力是各單管發(fā)射后坐力的矢量疊加。由于在不同大小的激勵作用下產(chǎn)生的振動頻率不同，如果某個振動頻率ωj(j表示發(fā)射管的數(shù)目)和發(fā)射器的固有頻率ω0相同或者相近，極易引起共振。

2.2 車輛行進(jìn)時的振動激勵

車輛行進(jìn)過程中，發(fā)射器的振動主要來自于路面不平所帶來的激勵。激勵通過輪胎、車架、懸架、車體的傳遞，作用于發(fā)射器，引起發(fā)射器的振動。以平均速度u通過路程為l的時間t=l/u，以平均速度通過波長為λ的路面時間T=λ/u為周期。T的倒數(shù)f=1/T=u/λ為頻率［3］。

路面不平度的時間頻率單邊功率密度函數(shù)為

由于T=l/u，故有

國際標(biāo)準(zhǔn)建議路面的空間頻率譜用下式擬合

故得到路面不平度時間頻率的單邊譜的擬合式為

3 建立發(fā)射器的有限元模型

由于九管發(fā)射器(圖1)的結(jié)構(gòu)比較復(fù)雜，很難實現(xiàn)在ANSYS中直接建模，所以首先在Pro/E中建立物理模型，然后將其導(dǎo)入ANSYS中進(jìn)行分析。

圖1 九管發(fā)射器模型Fig．1 Model of nine pipe launcher

由于結(jié)構(gòu)的固有頻率和模態(tài)振型主要取決于質(zhì)量分布和剛度，細(xì)節(jié)的影響不大。因此在建立三維數(shù)字化模型時，在不影響計算精度的前提下，適當(dāng)對模型進(jìn)行簡化，忽略倒角、圓角、螺紋孔等不必要的細(xì)節(jié)。

發(fā)射器的整體結(jié)構(gòu)過于復(fù)雜，而發(fā)射管對整個發(fā)射器的固有頻率起著主導(dǎo)作用，因此只對發(fā)射管進(jìn)行模態(tài)分析。根據(jù)發(fā)射管的實際尺寸在Pro/E中建立發(fā)射管的三維數(shù)字化模型(圖2)。然后利用接口將模型導(dǎo)入ANSYS的仿真平臺。整個模型采用8節(jié)點的solid185單元，它可模擬幾乎不可壓縮的彈塑和完全不可壓縮的超彈。采用smart size的方法對模型進(jìn)行智能網(wǎng)格化分。劃分后整個模型共有單元45 134個，節(jié)點23 446個。得到有限元模型如圖3所示。

圖2 三維數(shù)字化模型Fig．2 Three-dimensional digital model

圖3 有限元模型Fig．3 Finite element model

發(fā)射管的本質(zhì)模型是一個懸臂梁結(jié)構(gòu)，所以在發(fā)射管的后端施加全約束，根據(jù)物理屬性對模型進(jìn)行參數(shù)設(shè)置。發(fā)射管材料為鋼，其幾何參數(shù)、材料屬性由有關(guān)技術(shù)文件獲得。主要參數(shù)如表1所示。

表1 發(fā)射管模型的結(jié)構(gòu)尺寸和材料參數(shù)［4］Tab．1 Structure dimensions and material parameters of launch pipe model

根據(jù)材料參數(shù)對模型進(jìn)行材料參數(shù)相關(guān)設(shè)置。

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

模態(tài)分析理論吸取了振動理論、信號分析、數(shù)據(jù)處理、數(shù)理統(tǒng)計及自動控制中的有關(guān)“營養(yǎng)”，結(jié)合自身發(fā)展形成了一套獨特的理論。模態(tài)分析技術(shù)已經(jīng)廣泛地應(yīng)用在航天、航空、機(jī)械、造船、建筑、交通和兵器等工程領(lǐng)域。

模態(tài)分析的方法主要有兩種，一種是用有限元方法計算獲得，這種方法為數(shù)值模態(tài)分析。另外一種是通過實驗獲得采集信號和模態(tài)參數(shù)，這種方法為試驗?zāi)B(tài)分析。數(shù)值模態(tài)分析是將彈性結(jié)構(gòu)離散為有限個單元，根據(jù)單元的質(zhì)量和彈性特征在計算機(jī)上完成相關(guān)數(shù)學(xué)理論運算。通過這種方法，可以在產(chǎn)品生產(chǎn)出來之前對其性能進(jìn)行模擬，獲得產(chǎn)品的相關(guān)模態(tài)參數(shù)，根據(jù)獲得的數(shù)據(jù)對產(chǎn)品進(jìn)行優(yōu)化，可以縮短生產(chǎn)周期，降低生產(chǎn)成本［5］。本文采用此種方法進(jìn)行分析。

4.1 模態(tài)計算數(shù)學(xué)模型

根據(jù)達(dá)朗貝爾原理，對于一個有N個自由度的線性系統(tǒng)，其運動微分方程為:

根據(jù)達(dá)朗貝爾原理，對于一個有N個自由度的線性系統(tǒng)，其運動微分方程為:

(1)模態(tài)分析:F(t)=0。

(2)諧響應(yīng)分析:F(t)為周期載荷。

(3)瞬態(tài)動力學(xué)分析:F(t)為沖擊載荷。

(4)譜分析:F(t)為隨機(jī)載荷。

振動有F(t)=0，上式可以寫為:

若忽略阻尼影響，即C=0則可進(jìn)一步簡化為:

式中:［M］為質(zhì)量矩陣;［C］為阻尼矩陣;［K］為剛度矩陣;［X］為位移向量;{F(t)}為作用力向量;t為時間。

自由振動時各節(jié)點得的位移為:

其中:［Xmax］為位移振幅向量。

由式(3)、式(4)得:

式(5)是齊次方程，要想得到［Xmax］的非零解，只能使系數(shù)矩陣的行列式為0，即

由此式，可求出n個特征值和他們分別對應(yīng)的特征向量。其中特征值的平方根 ωi(i=1，2，…，n)就是結(jié)構(gòu)的第i階固有頻率，特征向量［］就是第i階的模態(tài)振型，它反映了結(jié)構(gòu)按頻率ωi振動時各自由度方向振幅間相對比例關(guān)系。［6－7］

4.2 模態(tài)仿真計算

在有限元模型建好后，進(jìn)入solution求解器，選擇分析類型為Modal。在ANSYS 13．0中，軟件提供了 Block Lanczos，Subspace，Power dynamics，縮減法，不對稱法，阻尼法6種模態(tài)提取方法。在選用提取方法時要根據(jù)模型的大小、計算機(jī)的運算能力等具體情況來定。Block Lanczos和其他的幾種方法相比計算速度更快，適合于大型對稱矩陣特征向量的求解，能使大部分問題得以解決，而且具有很高的精度，故本文采用此方法進(jìn)行提取。在這里對發(fā)射管的前12階模態(tài)進(jìn)行提取，在“No．of modes to extract”中輸入模態(tài)提取數(shù)目為12，在“No．of modes to expand”中輸入模態(tài)擴(kuò)展數(shù)目為12。

運用solution求解器進(jìn)行求解，通過計算得到了模型的前12階振動模態(tài)。由于低階模態(tài)對發(fā)射精度有較大影響，高階模態(tài)影響較小，所以在這里只分析1階到6階模態(tài)，可以得到模型的固有頻率和振型如表2所示。

表2 各階頻率和振型Tab．2 Frequencies and vibration

同時，各階的振型圖如圖4到圖9，圖中的位移僅表示振型，并非實際的位移大小。

圖4 一階模態(tài)振型Fig．4 1st modal shape

圖5 二階模態(tài)振型Fig．5 2nd modal shape

圖6 三階模態(tài)振型Fig．6 3rd modal shape

圖7 四階模態(tài)振型Fig．7 4th modal shape

4.3 結(jié)果分析

通過結(jié)果可以看出發(fā)射管的自由振動頻率分布在1 000Hz至1 500Hz之間，在使用過程中為避免共振，確保發(fā)射精度，必須盡量避開這一頻率。而對其振動產(chǎn)生影響的因素主要有發(fā)動機(jī)的頻率、彈體的發(fā)射頻率和路面的激勵。

圖8 五階模態(tài)振型Fig．8 5th modal shape

圖9 六階模態(tài)振型Fig．9 6th modal shape

(1)發(fā)動機(jī)的轉(zhuǎn)速通常在1 000～2 000 r/min，也就是振動頻率在 16．7～33．3 Hz［8］，與發(fā)射器固有頻率相差較大，不會產(chǎn)生共振，對發(fā)射精度的影響較小。

(2)在行車過程中，假設(shè)最高車速為80 km/h，路面不平度的最小波長為0．32(表3)。由于路面不平所帶來的激勵頻率為:在公路上行駛的激勵一般在15 Hz以下，在越野狀態(tài)下，路況較差，其產(chǎn)生的激勵也一般在50 Hz以下。與發(fā)射器固有頻率相差較大。

表3 路面的不平度波長［9］Tab．3 Pavement roughness wavelength

(3)彈體的發(fā)射周期為10 s左右，頻率在0．1 Hz，對發(fā)射影響較小。彈體的初速度為70 m/s，在發(fā)射管內(nèi)的飛行時間為 0．004 3 s［10］。而發(fā)射管的振動周期為 0．000 67～0．001 00 s。彈體的飛行時間避開了振動周期的1/2，即振幅的最大值，避免了共振的產(chǎn)生。

5 結(jié) 語

(1)本文通過模態(tài)分析獲得了發(fā)射管的固有頻率，分析了其可能產(chǎn)生共振的環(huán)境以及避免共振的措施，為下一步發(fā)射管的動力響應(yīng)分析提供了數(shù)據(jù)和技術(shù)基礎(chǔ)，具有可借鑒意義。

(2)以武警常用九管發(fā)射器為切入點，進(jìn)行動力學(xué)分析，在一定程度上彌補(bǔ)對非致命武器動力學(xué)研究不足的現(xiàn)狀，為其他車載非致命武器的系統(tǒng)動力學(xué)分析提供了思路和方法。

(3)本文只是獲得了發(fā)射器的固有頻率。不能很好地模擬其在整個發(fā)射過程的瞬態(tài)響應(yīng)情況。影響發(fā)射器發(fā)射精度的因素較多，需要在日后進(jìn)行進(jìn)一步研究。

［1］ 吳東亞，刑宏光．基于有限元的某型坦克炮身管模態(tài)分析［J］．科技導(dǎo)報，2008，26(23):44－47．WU Dongya，XING Hongguang．Modal analysis of a tank gun barrel based on finite element method［J］．Review of Science and Technology，2008，26(23):44－47．

［2］ 張海航，狄長安．某型多管火炮的振動模態(tài)分析［J］．兵工學(xué)報，2008，12(29):1514－1516．ZHANG Haihang，DI Changan．Analysis of vibration modal of a multibarrel cannon［J］．Acta Armamentaria，2008，12(29):1514－1516．

［3］ 毛保全．車載武器發(fā)射動力學(xué)［M］．北京:國防工業(yè)出版社，2010．MAO Baoquan．Launching dynamics of vehicle weapon［M］．Beijing:National Defense Industry Press，2010．

［4］ 曾正明．機(jī)械工程材料手冊［M］．北京:機(jī)械工業(yè)出版社，2010．ZENG Zhengming．Handbook of mechanical engineering materials［M］．Beijing:National Defense Industry Press，2010．

［5］ 袁安富，陳俊．ANSYS在模態(tài)分析中的應(yīng)用［J］．設(shè)計與研究，2007(8):79－81．YUAN Anfu，CHEN Jun．Application of ANSYSin model analysis［J］．Design and Research，2007(8):79－81．

［6］ 傅志方，華宏星．模態(tài)分析理論與應(yīng)用［M］．上海:上海交通大學(xué)出版社，2000．FU Zhifang，HUA Hongxing．Theory and application of modal analysis［M］．Shanghai:Shanghai Jiao Tong University Press，2000．

［7］ 徐立黃，張貴林．基于有限元的某多管火箭炮模態(tài)分析［J］．彈箭與制導(dǎo)學(xué)報，2010，6(30):110－115．XU Lihuang，ZHANG Guilin．The modal analysis of a MLRS based on FEM ［J］．Journal of Projectiles，Rockets，Missiles and Guidance，2010，6(30):110－115．

［8］ 喻凡，林逸．汽車系統(tǒng)動力學(xué)［M］．北京:機(jī)械工業(yè)出版社，2005．YU Fan，LIN Yi．Dynamics of vehicle system［M］．Beijing:Machinery Industry Press，2005．

［9］ 梁君，趙登峰．模態(tài)分析方法綜述［J］．現(xiàn)代制造工程，2006(8):139－141．LIANG Jun，ZHAO Dengfeng．Summary of method modal analysis［J］．Modern Manufacturing Engineering，2006(8):139－141．

［10］ 趙陜冬，馬永忠．非致命武器與警用器材［M］．北京:兵器工業(yè)出版社，2005．ZHAO Shandong，MA Yongzhong．Nonlethal weapons and police equipment［M］．Beijing:Weapon Industry Press，2005．