多管火箭發(fā)射動(dòng)力學(xué)研究

第25卷第5期兵工學(xué)報Vol 25 No 52004年9月ACTA ARMAMENTARIIsep.2004多管火箭發(fā)射動(dòng)力學(xué)研究芮筱亭*贠來(lái)峰**陸毓琪*陳衛東*王國平*趙克升**何斌陸文廣〔南京理工大學(xué)動(dòng)力工程學(xué)院江蘇南京210094·)(南京軍代局駐揚州地區軍代室)中國白城兵器試驗中心**)(哈爾濱工程大學(xué)”摘要從武器系統彈、炮、藥、環(huán)境一體化的角度研究了多管火箭武器發(fā)射動(dòng)力學(xué)建立了多管火箭武器剛體、彈性體耦合的多體系統發(fā)射動(dòng)力學(xué)模型和發(fā)射動(dòng)力學(xué)方程并給岀了增廣特征矢量的正交性條件力圖實(shí)現對多管火箭武器系統動(dòng)力響應的精確分析。用本文方法對某多管火箭武器動(dòng)力響應的仿真結果進(jìn)行了試驗驗證關(guān)鍵詞振動(dòng)與波迻管火箭炮火箭彈發(fā)射動(dòng)力學(xué)彈道正交性中圖分類(lèi)號TJ713;T760.13發(fā)射動(dòng)力學(xué)是研究武器系統在發(fā)射過(guò)程中的受 Marting模型中的一個(gè)致命弱點(diǎn)是未能考慮彈丸的力及運動(dòng)規律進(jìn)而研究控制受力與運動(dòng)規律旳—動(dòng)不平衡對彈丸靜不平衡的處理也不夠合理而這門(mén)新興學(xué)科是當代兵器領(lǐng)域科技工作者關(guān)注的熱兩個(gè)因素對彈丸起始擾動(dòng)研究恰恰是至關(guān)重要的點(diǎn)之_1國內外發(fā)射動(dòng)力學(xué)的研究文獻中相繼 Marting工作的另一個(gè)不足之處是描述彈丸運動(dòng)的提出了單自由度、多自由度等模型2實(shí)踐證明,動(dòng)參考系選取不佳從而使所建立的微分方程的形定向管的剛度對起始擾動(dòng)和密集度的影響極其明式非常復雜難以寫(xiě)出彈丸運動(dòng)微分方程的最終形顯必須作為彈性體來(lái)處理3式。本文針對 Marting方程的不足對 Marting模型為減少多管武器試驗用彈量實(shí)現其動(dòng)態(tài)設計,作了實(shí)質(zhì)性改進(jìn)建立了適用于各種類(lèi)型火箭彈的首先需解決含有剛體和彈性體(或剛性梯度很大發(fā)射動(dòng)力學(xué)方程。復雜耦合系統多管火箭武器)有振動(dòng)特性的計算冋題以及多管火箭剛彈耦合系統對任意激勵動(dòng)力響1多管火箭發(fā)射動(dòng)力學(xué)模型及受力分析應的精確分析。但是剛彈耦合系統由于剛體與彈1.1發(fā)射動(dòng)力學(xué)模型性體之間的耦合作用使得耦合系統的特征值問(wèn)題以某輪式多管火箭炮為研究對象發(fā)射環(huán)境下非自共軛系統的振型函數不具有通常意義下的正武器系統包括火箭炮、火箭彈和陣地祭件?；鸺诮恍詮亩y以用經(jīng)典的模態(tài)方法精確分析耦合系主要由定向管、搖架、平衡機、高低傳動(dòng)裝置、方向傳統的動(dòng)力響應問(wèn)題。為此本文引入了多體系統傳動(dòng)裝置、回轉盤(pán)、座圈、底架、車(chē)體、板簧固定器、輪胎遞矩陣法4。這種方法便于實(shí)現發(fā)射動(dòng)力學(xué)分析以及電傳動(dòng)系統等組成。的高度程式化易于工程問(wèn)題的解決得到了國內外根據多管火箭武器各個(gè)部件的自然屬性將其的重視劃分?？煞譃轶w”和鉸”兩1984年美國動(dòng)力學(xué)專(zhuān)家 Marting拋棄了彈帶大類(lèi)中國煤化彈性梁等而鉸"泛或后定心部沖心沿炮膛軸線(xiàn)運動(dòng)”的傳統假設建指犀<問(wèn)聯(lián)接關(guān)系,包括扭簧、彈立了彈丸在膛內作一般運動(dòng)的微分方程但簧、阻尼器等。如圖1.1所示根據體”和鉸"統2003年7月收稿2004年7月定稿κ國防科技預研重點(diǎn)項目(98sY05)第5期多管火箭發(fā)射動(dòng)力學(xué)研究557編號的原則將地面邊界編號為θ將六個(gè)車(chē)輪與地彈炮接觸力作為強迫力”考慮其對系統響應的影面及每個(gè)車(chē)輪與車(chē)體的聯(lián)接用反映三個(gè)方向相對線(xiàn)響而不計其對系統固有振動(dòng)特性的影響。這樣處理運動(dòng)的彈簧及與之并聯(lián)的阻尼器來(lái)等效編號為元的原因一是理論上不存在一般線(xiàn)性時(shí)變系統固有振件12灬613,14灬,8將六個(gè)車(chē)輪視為各具有動(dòng)特性嚴格的計算方法;二是用此法處理結果足以個(gè)自由度的集中質(zhì)量編號為元件78灬,,12將滿(mǎn)足工程精度要求(單發(fā)火箭彈質(zhì)量?jì)H僅是起落部除去車(chē)輪的車(chē)體(簡(jiǎn)稱(chēng)車(chē)體)除去俯仰部分的回轉分的1.5%~4.0%)理論與試驗結果表明發(fā)射過(guò)部分、除去已發(fā)射火箭彈和最新?lián)舭l(fā)火箭彈的起落程中相鄰兩發(fā)彈之間系統質(zhì)量及其分布變化引起的部分分別視為各具有六個(gè)自由度的剛體依次編號系統固有振動(dòng)頻率的變化很小故無(wú)需考慮運動(dòng)火為1921②3將方向機的作用以及元件19與元件箭彈在定向管中位置的變化而引起的火箭炮系統模21間的聯(lián)接、高低機和平衡機的作用以及元件21態(tài)的變化。非線(xiàn)性時(shí)變系統的動(dòng)力響應也可用離散與元件23的聯(lián)接用反映三個(gè)方向相對角運動(dòng)的無(wú)時(shí)間傳遞矩陣湛6求解。質(zhì)量的扭簧和反映三個(gè)方向相對線(xiàn)運動(dòng)的彈簧及與1.2受力分析之并聯(lián)的阻尼器來(lái)等效編號為20、22將將搖架與定多管火箭系統受到的力有重力、各部件之間以向管間的兩處聯(lián)接用反映三個(gè)方向相對角運動(dòng)的扭及火箭炮與地面之間的機械作用力、彈炮作用力包簧和反映三個(gè)方向相對線(xiàn)運動(dòng)的彈簧及與之并聯(lián)的括閉鎖力、火箭彈定心部與定向管間以及定向鈕與阻尼器來(lái)等效編號分別為17+7、18+7將第ⅰ導向槽之間的作用力等)發(fā)動(dòng)機推力、燃氣流沖擊個(gè)定向管的后端面自由邊界編號為19+7i將第i力等?；鸺诟鞑考g的作用力、火箭炮與地面之定向管尾部視為具有六個(gè)自由度的剛體編號為間的作用力均可等價(jià)為彈性鉸對應的作用力和相應20+7i將第i個(gè)定向管視為空間運動(dòng)彈性梁前支的阻尼力它們都是系統內力。本文用多體系統傳遞撐框與后支撐框之間的部分編號為21+7前支撐矩陣法求解多管火箭武器發(fā)射動(dòng)力學(xué),所以地面支框前面的部分編號為22+7i將第讠個(gè)定向管的前撐已被包含在系統動(dòng)力學(xué)模型中這些系統內力均端面自由邊界編號為23+7i?；鸺龔椀闹黧w為剛已在狀態(tài)矢量中予以考慮,無(wú)需進(jìn)行復雜而繁瑣的體其彈性效應等效為其定向鈕以及3個(gè)定心部與受力表達定向管的彈性接觸作用?？紤]火箭彈燃氣流對發(fā)射將除去已發(fā)射和已擊發(fā)火箭彈的多管火箭系統系統的作用等所以該多管火箭武器的發(fā)射動(dòng)力學(xué)的靜平衡位置作為系統位移的基準這樣系統重力模型為在地面支撐條件和燃氣流及彈炮耦合作用下(包括各部件的重力)在動(dòng)力學(xué)方程中的作用就已的由各種鉸相聯(lián)接的n+4個(gè)剛體、n個(gè)彈性體和6被系繚包括各部件)的靜變形所抵消故在這樣選個(gè)集中質(zhì)量組成的剛彈耦合多體系統。取的坐標系中將無(wú)需再考慮火箭炮的重力以及靜變22+7i形對應項的作用。使受力分析和多管火箭系統運動(dòng)21+7i求解大為簡(jiǎn)捷。0+T火箭彈與定向管間的作用力是一對作用力與反l7+7i19+7i作用力?；鸺龔椷\動(dòng)方程與火箭炮振動(dòng)方程正是通過(guò)定心部與定向管壁及定向鈕與導向槽之間的作用力而相互耦合的這些作用力包括定向鈕與定向槽間以及定心部與定向管間的法向碰撞力和切向摩擦力。被擊發(fā)火箭彈的重力通過(guò)彈炮相互作用傳遞給圖1.1某多管火箭武器發(fā)射動(dòng)力學(xué)模型火箭炮。ig. 1. 1 A model for the launch dynamics of mlrs發(fā)動(dòng)機推力可理論計算也可實(shí)驗測量;火箭燃氣流對發(fā)射裝置的沖擊力可由文執78]中的方法發(fā)射過(guò)程中由于搖架上火前彈個(gè)數的減少和計算火箭彈在定向管內的運動(dòng),發(fā)射系統為非線(xiàn)性時(shí)變中國煤化工了相關(guān)的試驗結果。系統。本文將非線(xiàn)性時(shí)變系統處理為分段線(xiàn)性時(shí)不2CNMHG變系統即考慮系統模態(tài)隨著(zhù)搖架上火箭彈個(gè)數的2.1體元件的動(dòng)力學(xué)方程減少而變化但不計發(fā)射過(guò)程中運動(dòng)火箭彈在定向集中質(zhì)量、剛體、彈性梁等體元件動(dòng)力學(xué)方程的管中位置變化引起的系統模態(tài)的變化將火箭彈在矩陣形式為定向管內的運動(dòng)、彈炮間隙、彈炮碰撞產(chǎn)生的非線(xiàn)性M,v; at Cva+K; v;=f558兵工學(xué)報第25卷式中,元件編號;體元件位移(包括線(xiàn)位移和角位移)的物理坐標列陣,下標t表示關(guān)于時(shí)間求u ,i導M、C、K,體元件的參數矩陣沂體元件受到D'I121i的外力(2.3)2.2集中質(zhì)量元件的動(dòng)力學(xué)方程xbj六個(gè)車(chē)輪元件78灬…,12)是一端輸入一端輸出的集中質(zhì)量將其空間運動(dòng)動(dòng)力學(xué)方程寫(xiě)成(2.1q2 o式的形式得到式中n,q,q2分別為體元件聯(lián)接點(diǎn)上受到的彈性?xún)攘Α?.3剛體元件的動(dòng)力學(xué)方程(2.2)元件19212320+7i均為空間一般運動(dòng)剛體其中元件19為六端輸入一端輸出剛體元件21∫=[fffT20+7為一端輸入一端輸出剛體元件23為一端輸式中=78灬12涅表示輸入端O表示輸出端;入兩端輸出剛體。將這些元件的動(dòng)力學(xué)方程寫(xiě)成D3為線(xiàn)性微分算子其表達式為(2.1)式的形式得到D1+乙D3|oF(j=19212320+7i)式中,m;剛體元件質(zhì)量江1第一輸入點(diǎn)匯質(zhì)心:的動(dòng)力學(xué)方程寫(xiě)成2.1)式的形式則有D外力作用點(diǎn),剛體元件關(guān)于I1點(diǎn)的慣量矩000陣m;=[mxm,m2于作用于剛體的外力矩坐標列陣洱=[f/作用于剛體的外力坐標M=0,0000m0列陣氵坐標叉乘矩陣;輸入端個(gè)數;L輸岀端000個(gè)數;D1線(xiàn)性微分算子其表達式為D11,0,D11.0y dDIog, d j(2.5)EAi(a/xt)0 EI(a/a xt式中m2I Ily,I i2n于m2Dmom=o1分中國煤化工別為體元件在輸入點(diǎn)和輸出點(diǎn)處受到的彈性?xún)攘NMHG矩f(x14)-(x1m(x1d)2.4梁元件振動(dòng)的動(dòng)力學(xué)方程f(x1t)+(/0x1)m元件21+7和22+7i作為彈性梁研究橫向0aMV+V劃t)+應的精確分析q(t)=f(5.2)(5)用本文方法對某多管火箭武器動(dòng)力響應的仿真結果得到了試驗驗證在(5.2)式兩邊對W(p=123)取內積由正(6)對某多管火箭武器的發(fā)射過(guò)程的仿真結果交性條件可得表明本文的處理方法足以滿(mǎn)足工程精度要求。(t)+(a+A2)(t)+or(t)=(f》M(p=12灬…)(5.3)參考文獻將3.1)式與(5.3)式聯(lián)立(3.1)式中的有關(guān)參量用模態(tài)坐標展開(kāi)可得關(guān)于彈丸運動(dòng)參量和廣1芮筱亭等.多體系統發(fā)射動(dòng)力學(xué).北京:國防工業(yè)出版社,義坐標《(tk=123灬…,)的常微分方程組。對中國煤化工北京北京理工大學(xué)出版社其進(jìn)行數值積分就可得到q(t)和彈丸的運動(dòng)參CNMHG數進(jìn)而用4.1)式得到多管火箭的動(dòng)力響應3王耀龍.對提高戰術(shù)火箭射擊密集度的分析.兵工學(xué)報火炮分冊1979(1):51-666計算和試驗結果驗證及結論4芮筱亭等.多體系統傳遞矩陣法在火炮動(dòng)力學(xué)中的應用.力學(xué)與實(shí)踐,1995,7(4):42-44利用本文的理論和方法對某多管火箭武器的發(fā)5 Marting T S, Robert s b. Projectile Motion in a Flexible Gt第5期多管火箭發(fā)射動(dòng)力學(xué)研究561Tube. ADA140737, 19842306 Rui xiaoting, Lu Yuqi,etc. Discrete time transfer matrix method8賀北斗林永明曹聽(tīng)榮.火箭發(fā)射裝置設計.北京國防工業(yè)出for multibody system dynamics. Advances in Computational版社1988.18Multibody Dynamics,DMEC/IST, Lisbon, portugal,1999。939徐明友火箭外彈道學(xué),北京浜兵器工業(yè)出版社A989.133-15410贠來(lái)峰芮筱亭陸毓琪等.多管火箭武器振動(dòng)特性研究兵工學(xué)7張福祥.火箭燃氣流動(dòng)力學(xué).北京國防工業(yè)出版社,1988.185報2004,2(1):9~13A STUDY ON THE LAUNCH DYNAMICS OF MULTIPLELAUNCH ROCKET SYSTEMRui Xiaoting Yun Laifeng Lu yugi Cheng weidongWang guoping Zhao Kesheng He Bin Lu wenguangNanjing University of Science and Technology Nanjing, 210094)( Military Representative Office in Yangzhou Area)Baicheng Ordnance Test Center of China***)( Harbin Engineering University**')Abstract Launch dynamics of Multiple Launch Rocket Systen( MlrS was studied in theory from thepoint of view of the whole system including the rockets launcher, propellant and environment. The launchdynamics model and launch dynamics equation of the multibody coupled rigid bodies and elastic bodies in theMLRS were established. Orthogonality conditions of augmented eigenvector were presented. Response ofthe mlrS was simulated, and the simulation of response of the MlrS was verified The basis to analysisof dispersion and reduction of rocket in the test of mlrS was established in this workKey words vibration and wave multiple launch rocket system rocket launch dynamics ballistics ,orthogonalit中國煤化工CNMHG