石 林,馬 威,楊龍軍
(江西洪都航空工業(yè)集團有限責(zé)任公司,南昌 330024)
基于ADAMS的無人機零長發(fā)射動力學(xué)仿真及分析*
石 林,馬 威,楊龍軍
(江西洪都航空工業(yè)集團有限責(zé)任公司,南昌 330024)
針對無人機零長發(fā)射系統(tǒng)研制需要,以無人機零長助推火箭發(fā)射段為主要研究對象,利用ADAMS建立了無人機發(fā)射動力學(xué)仿真模型,分別對質(zhì)心偏心和助推火箭推力偏心工況進行動力學(xué)仿真,獲取了發(fā)射過程中無人機質(zhì)心位移和角度的變化規(guī)律,分析了無人機與發(fā)射架的安全距離。結(jié)果表明:無人機質(zhì)心偏心和推力偏心影響其運動姿態(tài),且推力偏心的影響更明顯,無人機在發(fā)射過程中能夠安全離開發(fā)射架。
無人機;零長發(fā)射;仿真;模型
無人機的快速發(fā)展引發(fā)了層出不窮的軍事學(xué)術(shù)、裝備技術(shù)等相關(guān)問題的研究,并在未來會有越來越廣泛的應(yīng)用[1]?;鸺屏汩L發(fā)射方式是地面發(fā)射無人機時用得較多的發(fā)射方式,該發(fā)射方式具有便捷、快速,對發(fā)射場地?zé)o嚴(yán)格要求,適合于野戰(zhàn)環(huán)境等優(yōu)點,與發(fā)射方式配套的發(fā)射系統(tǒng)同時由于體積小、機動性強、適應(yīng)性強等特點,使其在無人機系統(tǒng)中具有重要而特殊的地位[2-3]。
對于研制階段的無人機系統(tǒng),開展零長發(fā)射試驗來確定各種設(shè)計參數(shù)還比較困難,采用ADAMS仿真軟件,以多體系統(tǒng)動力學(xué)理論為基礎(chǔ),把無人機和發(fā)射系統(tǒng)簡化為多剛體系統(tǒng)[4],對無人機質(zhì)心偏心和推力偏心工況進行發(fā)射動力學(xué)仿真,為驗證無人機零長發(fā)射可行性和安全性提供了有效可行的方法,對無人機總體作戰(zhàn)性能的提高具有重要的意義。
無人機零長發(fā)射系統(tǒng)按照結(jié)構(gòu)特點分為發(fā)射架和無人機兩大部分。無人機機身前方和后部共有4個支撐軸約束在發(fā)射架上,無人機下方后端裝有助推火箭,無人機零長發(fā)射系統(tǒng)如圖1所示。發(fā)射架主要由發(fā)射架體、翻轉(zhuǎn)機構(gòu)、后支撐、剪切銷機構(gòu)等組成,其中發(fā)射架體是發(fā)射架的主承力件,為發(fā)射架各機構(gòu)的安裝平臺,翻轉(zhuǎn)機構(gòu)是無人機前支撐軸約束點,后支撐是無人機后支撐軸約束點,當(dāng)助推火箭達到一定的推力情況下,剪切銷剪斷,無人機脫離翻轉(zhuǎn)機構(gòu)的約束而離開發(fā)射架。
圖1 無人機與發(fā)射架
2.1 模型假設(shè)
為了研究問題的需要,計算中模型作了如下假設(shè):
1)將無人機和發(fā)射架作為剛體,在運動過程中,形狀、質(zhì)量、質(zhì)心都不發(fā)生改變[5];
2)仿真忽略助推火箭與無人機脫離現(xiàn)象;
3)不考慮振動問題,忽略發(fā)射架和無人機的振動;
4)不考慮風(fēng)的影響。
2.2 模型導(dǎo)入及約束添加
將無人機與發(fā)射架的三維數(shù)模簡化,然后導(dǎo)入ADAMS中,將無人機和發(fā)射架均作為剛體,添加約束條件,發(fā)射架體與大地固定連接,后支撐與發(fā)射架體固連,助推火箭和無人機作為一個整體,用固定副連接;剪切銷機構(gòu)轉(zhuǎn)桿與發(fā)射架體之間為轉(zhuǎn)動關(guān)系,翻轉(zhuǎn)機構(gòu)與發(fā)射架體之間為轉(zhuǎn)動關(guān)系,用旋轉(zhuǎn)副連接;由于無人機前、后支撐軸與翻轉(zhuǎn)機構(gòu)、后支撐之間相對運動時會發(fā)生碰撞,因此需要創(chuàng)建接觸(contact)進行約束。
其中,以無人機質(zhì)心為原點,以平行于地面的無人機行進方向為X軸,以垂直向上方向為Y軸,根據(jù)笛卡爾坐標(biāo)系建立Z軸。在ADAMS中建立無人機的發(fā)射動力學(xué)仿真模型。
2.3 設(shè)置仿真參數(shù)
無人機發(fā)射仰角為10°,無人機軸線和助推火箭軸線夾角為13°,助推火箭推力按常溫工作狀態(tài)下的推力曲線進行加載,在ADAMS中推力加載調(diào)用格式為:AKISPL(time,0,SPLINE_CW,0)。
取翻轉(zhuǎn)機構(gòu)扭簧扭矩為17.5 N·m,取剪切銷機構(gòu)扭簧扭矩為1.8 N·m;閉鎖力為4 000 N,根據(jù)推力曲線可知,靶機解鎖時間約為0.005 s,將剪切銷設(shè)置成固定副,利用仿真控制腳本語句來描述,控制描述語句為:SIMULATE/DYNAMIC,END=0.005,STEPS=10;DEACTIVATE/JOINT,ID=28。設(shè)置仿真時間End Time為0.5 s,仿真步長Steps為1000。
2.4 零長發(fā)射動力學(xué)方程的建立
ADAMS建模方法采用第二類拉格朗日方程,選取系統(tǒng)內(nèi)每個剛體質(zhì)心的笛卡爾坐標(biāo)和確定剛體方位的歐拉角作為廣義坐標(biāo),即:
則對每個剛體,建立對應(yīng)于6個廣義坐標(biāo)帶乘子的拉格朗日方程為:
完整約束方程時:f(q,t)=0
3.1 無人機發(fā)射仿真結(jié)果
1)理想工況
在理想工況下,此時質(zhì)心偏心0 mm,助推火箭推力偏心0′時,無人機質(zhì)心的位移和角度隨時間的變化規(guī)律,如圖2和圖3所示。
圖2 無人機質(zhì)心位移曲線
圖3 無人機角度曲線
2)無人機質(zhì)心偏心計算
分別取無人機質(zhì)心側(cè)向偏心0.2 mm、0.4 mm、0.6 mm工況,對其在0.5 s發(fā)射階段進行動力學(xué)仿真,得出無人機發(fā)射過程中的位移和角度隨時間的變化規(guī)律。仿真結(jié)果表明:隨著無人機質(zhì)心偏移距離的增大,無人機質(zhì)心X、Y向位移的變化不明顯,Z向位移偏移量隨之增加;俯仰角的變化不明顯,滾轉(zhuǎn)角和偏航角的偏移量也隨之增加。由于滾轉(zhuǎn)方向轉(zhuǎn)動慣量遠小于偏航方向,使?jié)L轉(zhuǎn)角的偏移量大于偏航角,因此,質(zhì)心偏心對滾轉(zhuǎn)角影響更大。
Z向位移、滾轉(zhuǎn)角和偏航角的變化規(guī)律如圖4~圖6所示。
圖4 無人機質(zhì)心Z向位移曲線
圖5 無人機質(zhì)心X向滾轉(zhuǎn)角曲線
圖6 無人機質(zhì)心Y向偏航角曲線
3)助推火箭推力偏心計算
分別取無人機助推火箭推力偏心3′,6′,9′工況,對其在0.5 s發(fā)射階段進行動力學(xué)仿真,得出無人機發(fā)射過程中的位移和角度隨時間的變化規(guī)律。仿真結(jié)果表明:隨著無人機助推火箭推力偏心角度的增大,無人機質(zhì)心X、Y向位移的變化不明顯,Z向的位移偏移量隨之增加;俯仰角的變化不明顯,滾轉(zhuǎn)角和偏航角的偏移量也隨之增加。若推力偏心角度較大,導(dǎo)致滾轉(zhuǎn)角和偏航角偏移量過大,不利于無人機的穩(wěn)定發(fā)射,因此,在研制設(shè)計時需要控制助推火箭推力的偏心角度。
Z向位移、滾轉(zhuǎn)角和偏航角的變化規(guī)律如圖7~圖9所示。
圖7 無人機質(zhì)心Z向位移曲線
圖8 無人機質(zhì)心X向滾轉(zhuǎn)角度曲線
圖9 無人機質(zhì)心Y向偏航角度曲線
3.2 發(fā)射安全距離分析
無人機發(fā)射初始階段,無人機助推器與發(fā)射架距離存在一定風(fēng)險,助推火箭與發(fā)射架可能發(fā)生干涉的位置如圖10所示。
圖10 無人機與發(fā)射架距離較短位置
由仿真結(jié)果可得:各工況下,無人機在0.203s時刻,助推火箭后端到后支撐側(cè)向的距離L1和L2達到最小值;無人機在0.129 2s時刻,尾翼到導(dǎo)軌的距離L3達到最小值,間距如表1所示。
表1 助推火箭與發(fā)射架間距
表1表明,在質(zhì)心和火箭推力偏心工況下,隨側(cè)向偏移量增大,助推火箭到側(cè)支架的距離L1和L2之間的偏差越來越大;在質(zhì)心偏心0.6mm工況下偏移量為3.05mm,在推力偏心9′工況下,偏移量達到了10.35mm;尾翼到導(dǎo)軌的距離L3變化不明顯,偏移值在設(shè)計值范圍內(nèi),滿足發(fā)射安全性要求。無人機與發(fā)射架之間的距離L1、L2和L3均大于最小安全設(shè)計值,因此無人機在發(fā)射過程中能夠安全離開發(fā)射架,不會發(fā)生干涉。
仿真結(jié)果表明:無人機質(zhì)心偏心和助推火箭推力偏心對無人機位移和角度有一定的影響,且推力偏心的影響更明顯,因此在無人機系統(tǒng)研制過程中,質(zhì)心和推力偏心均不能過大,控制助推火箭與無人機的安裝精度尤為重要;根據(jù)發(fā)射安全距離分析可知:無人機與發(fā)射架不會發(fā)射干涉,最小安全距離也滿足零長發(fā)射要求。
[1] 盛懷潔.無人機發(fā)射回收方式面面觀 [J].無人機,2004(3):42-44.
[2] 余勇.無人機是如何起飛和回收的 [J].兵工科技,2006(10):60-61.
[3] 馬威,過海峰,姜寧.海軍通用無人機及其起降方式分析 [J].飛航導(dǎo)彈,2006(12):37-40.
[4] 張勝三,顧銀芳.傾斜零長發(fā)射技術(shù)研究 [J].導(dǎo)彈與航天運載技術(shù),2009(6):41-45.
[5] 廖莎莎,吳成.機載導(dǎo)彈發(fā)射動力學(xué)建模與虛擬樣機仿真 [J].北京理工大學(xué)學(xué)報,2011,31(9):1013-1017.
SimulationandAnalysisofUAVZero-lengthLaunchingDynamicsBasedonADAMS
SHI Lin,MA Wei,YANG Longjun
(Jiangxi Hongdu Aviation Industry Group Corporation Limited,Nanchang 330024,China)
According to the requirement of UAV zero-length launching system development,the launch section of UAV zero-length lauching was taken as the main research object,and Adams was used to establish UAV launching dynamics simulation model.The dynamic simulation of centroid eccentricity and the thrust eccentric condition of booster rocket was carried out,and the centroid displacement and the angle change law of UAV in launching process were obtained.And the safety distance between the UAV and launcher was analyzed.The results showed that the centroid eccentricity and thrust eccentricity affected the motion attitude of UAV,and the influence of thrust eccentricity was more obvious.The UAV could safely leave the launcher during launching.
UAV; zero-length launching; simulation; model
10.15892/j.cnki.djzdxb.2017.02.009
2016-05-05
石林(1980-),男,江西鄱陽人,高級工程師,研究方向:發(fā)射、機械、液壓、仿真。
V553
A