• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    助推-滑翔導(dǎo)彈彈道優(yōu)化與總體參數(shù)分析

    2012-12-25 08:47:24張青斌
    彈道學(xué)報(bào) 2012年3期
    關(guān)鍵詞:彈頭滑翔射程

    劉 欣,楊 濤,張青斌

    (國(guó)防科學(xué)技術(shù)大學(xué) 航天與材料工程學(xué)院,長(zhǎng)沙410073)

    助推-滑翔導(dǎo)彈通過火箭助推加速,依靠滑翔彈頭的高升阻比外形,利用空氣動(dòng)力控制飛行軌跡,實(shí)現(xiàn)遠(yuǎn)距離的非彈道式再入機(jī)動(dòng)飛行.與彈道導(dǎo)彈相比,由于其在增大射程、突防、再入段機(jī)動(dòng)等方面具有優(yōu)勢(shì),助推-滑翔導(dǎo)彈逐漸成為近來的研究熱點(diǎn)[1].

    彈道優(yōu)化設(shè)計(jì)是助推-滑翔導(dǎo)彈總體設(shè)計(jì)的關(guān)鍵技術(shù)之一.助推-滑翔導(dǎo)彈的彈道可初步分為主動(dòng)段、自由飛行段、再入段,目前國(guó)內(nèi)外相關(guān)研究集中在滑翔彈頭再入彈道優(yōu)化設(shè)計(jì)方面,針對(duì)導(dǎo)彈全彈道優(yōu)化的研究較少.文獻(xiàn)[2]利用近似方法對(duì)全球到達(dá)高超聲速飛行器的彈道特性進(jìn)行了分析,文獻(xiàn)[3]利用序列二次規(guī)劃法對(duì)助推-滑翔導(dǎo)彈縱向全彈道進(jìn)行最大射程優(yōu)化.由于助推-滑翔導(dǎo)彈彈道的復(fù)雜性,利用解析法進(jìn)行總體參數(shù)分析十分困難,且誤差很大,結(jié)合各種優(yōu)化技術(shù)對(duì)助推-滑翔導(dǎo)彈的彈道特性與總體參數(shù)進(jìn)行分析是一種趨勢(shì).

    本文利用Gauss偽譜法對(duì)全彈道進(jìn)行優(yōu)化設(shè)計(jì),并分析了起飛質(zhì)量、主動(dòng)段終端傾角等參數(shù)對(duì)導(dǎo)彈射程的影響.

    1 助推器參數(shù)估算

    為了完成主動(dòng)段彈道的計(jì)算,建立起飛質(zhì)量與主動(dòng)段終點(diǎn)參數(shù)乃至導(dǎo)彈射程的聯(lián)系,需要對(duì)助推器的基本參數(shù)進(jìn)行估算.

    假設(shè)助推-滑翔導(dǎo)彈采用單級(jí)固體火箭助推.若初始速度為零,根據(jù)齊奧爾柯夫斯基公式,單級(jí)火箭的理想速度增量為[4]

    式中:u′e為有效排氣速度,數(shù)值上等于真空比沖Ip;mf為載荷質(zhì)量,m0為火箭起飛質(zhì)量,m0-mf為助推器質(zhì)量;結(jié)構(gòu)系數(shù)ε為助推器結(jié)構(gòu)質(zhì)量與助推器質(zhì)量之比.考慮飛行過程中由引力、氣動(dòng)力和大氣靜壓引起的飛行器速度損失Δv,則要使火箭主動(dòng)段終點(diǎn)速度達(dá)到vf,火箭的理想速度增量應(yīng)為vid=vf+Δv,可求得起飛質(zhì)量與目標(biāo)速度、載荷質(zhì)量的關(guān)系:

    依據(jù)文獻(xiàn)[5]中的相關(guān)論述,可取u′e=Ip=2 800m/s,ε=0.1,Δv=0.2vf.對(duì)于給定的彈頭,式(2)建立了起飛質(zhì)量與主動(dòng)段終點(diǎn)速度的關(guān)系.

    假設(shè)助推系統(tǒng)的任務(wù)是將質(zhì)量為907.8kg的滑翔彈頭[6]助推至速度4 500m/s.根據(jù)式(2)可求得起飛質(zhì)量m0=18 013kg,進(jìn)一步可計(jì)算出助推器裝藥質(zhì)量為mp=(m0-mf)(1-ε)=15 394.7kg.通常導(dǎo)彈設(shè)計(jì)時(shí)需對(duì)主動(dòng)段重推比(ν0=m0g/(mpIp/t1))加以選擇.參考一般固體火箭取值,重推比可取為ν0=0.5,據(jù)此可求出助推器工作時(shí)間t1=s,燃料質(zhì)量流量平均推力Fp=Ip=353.08kN.

    2 助推-滑翔導(dǎo)彈彈道優(yōu)化模型

    彈道優(yōu)化設(shè)計(jì)的目的是獲取一條或一類滿足約束和性能指標(biāo)條件的最優(yōu)彈道.彈道優(yōu)化問題可描述為一般的最優(yōu)控制問題,即在時(shí)間區(qū)間[t0,tf]中,尋找最優(yōu)控制變量u(t),最小化性能指標(biāo)J,并使?fàn)顟B(tài)變量x(t)、初始時(shí)間t0、終端時(shí)間tf滿足運(yùn)動(dòng)微分方程約束以及端點(diǎn)約束和過程約束.下面結(jié)合助推-滑翔導(dǎo)彈的特點(diǎn),建立運(yùn)動(dòng)微分方程,確定約束條件及性能指標(biāo).

    1)性能指標(biāo).

    助推-滑翔導(dǎo)彈彈道復(fù)雜,宜分段進(jìn)行彈道優(yōu)化.滑翔彈頭的主要優(yōu)勢(shì)即在于增大射程,總體設(shè)計(jì)時(shí)比較關(guān)心其射程覆蓋范圍,以Lf表示彈道終點(diǎn)處射程,可將滑翔彈頭飛行段的性能指標(biāo)取為J=min(-Lf),即使射程最大.

    主動(dòng)段的主要作用是加速滑翔彈頭,其射程在全彈道射程中所占比例很小,可將性能指標(biāo)取為J=-min[(gH+v2/2)|tf],即主動(dòng)段終點(diǎn)機(jī)械能最大,以利于滑翔彈頭飛行段增大射程.

    2)狀態(tài)方程.

    總體性能分析時(shí)主要關(guān)注導(dǎo)彈縱平面的運(yùn)動(dòng),忽略地球自轉(zhuǎn)和非球形影響,建立助推-滑翔導(dǎo)彈縱平面運(yùn)動(dòng)方程:

    式中:H為飛行器距離地面的高度;L為射程;Θ為當(dāng)?shù)厮俣葍A角;v為飛行器飛行速率;m為飛行器質(zhì)量;e為發(fā)動(dòng)機(jī)燃料質(zhì)量流量;FX,F(xiàn)Y分別為氣動(dòng)阻力、升力;r為飛行器到地心的距離;μ為地球引力系數(shù);φp為推力與速度方向的夾角.

    在優(yōu)化問題中取狀態(tài)變量x=(HLvΘ m)T,主動(dòng)段通過改變推力方向控制彈道,取控制量為φp;滑翔彈頭段通過氣動(dòng)舵來改變彈道,而在僅考慮縱平面運(yùn)動(dòng)的情況下,飛行器動(dòng)力取決于飛行器的速度、高度及攻角α,則取控制量為攻角α.

    3)約束條件.

    彈道設(shè)計(jì)時(shí)主要考慮飛行過程約束、終端參數(shù)約束及控制量約束.飛行過程約束,也就是飛行器在飛行過程中彈道參數(shù)必須滿足的約束.為確保飛行器安全飛行,通常要考慮熱流、法向過載、動(dòng)壓等因素對(duì)彈道的約束.具體有,法向過載約束:n=,動(dòng)壓約束駐點(diǎn)熱流約束=ksρ0.5v3.15<max.其中,ny,max,qmax,max分別為飛行中所允許的法向過載、動(dòng)壓和熱流的最大值.

    終端參數(shù)約束是指飛行器在彈道端點(diǎn)需滿足的條件,即最優(yōu)控制問題中的邊界條件.主動(dòng)段彈道起點(diǎn)參數(shù)取為導(dǎo)彈垂直上升段終點(diǎn)狀態(tài),滑翔彈頭飛行段起點(diǎn)參數(shù)應(yīng)當(dāng)與主動(dòng)段終點(diǎn)參數(shù)一致,滑翔彈頭飛行段末端通常有落地速度與落地彈道傾角的要求,可取-Θ(tf)=-60°,v(tf)=1 000m/s.

    另外,攻角、推力與速度方向的夾角等控制量的幅值和變化率不能超過限制值,可?。粒?5°,|φp|<5°,||<1 (°)/s,|p|<1 (°)/s.

    3 Gauss偽譜法優(yōu)化基本原理

    助推-滑翔導(dǎo)彈彈道形式復(fù)雜,飛行約束多,其彈道優(yōu)化是一個(gè)難點(diǎn).Gauss偽譜法(Gauss Pseudospectral Method,GPM)是一種適合解決多階段多約束最優(yōu)控制問題的數(shù)值優(yōu)化方法,它將狀態(tài)變量和控制變量在一系列Legendre Gauss(LG)點(diǎn)上離散,構(gòu)造全局插值多項(xiàng)式來近似狀態(tài)變量和控制變量,經(jīng)過一系列處理,最終將最優(yōu)控制問題轉(zhuǎn)化為具有一系列代數(shù)約束的參數(shù)優(yōu)化問題,即非線性規(guī)劃問題(NLP),并進(jìn)行數(shù)值求解.GPM離散最優(yōu)控制問題的主要流程如下[7,8].

    1)時(shí)域變換.

    由于Gauss偽譜法的配點(diǎn)都分布在區(qū)間[-1,1],因此求解中一般對(duì)時(shí)間區(qū)間進(jìn)行變換,通過引入變換:

    將最優(yōu)控制問題的時(shí)間區(qū)間由t∈[t0,tf]轉(zhuǎn)換至τ∈[-1,1],并將原問題按時(shí)間歸一化.

    2)全局插值多項(xiàng)式近似狀態(tài)與控制變量.

    Gauss偽譜法的離散點(diǎn)κ={τ1,…,τK}為K階Legendre-Gauss點(diǎn),即K階Legendre多項(xiàng)式的根PK(τ),其中,

    Legendre-Gauss點(diǎn)分布在區(qū)間(-1,1),再增加一個(gè)點(diǎn),τ0=-1,作為區(qū)間[-1,1)上的K+1個(gè)插值點(diǎn),構(gòu)成K+1個(gè)Lagrange插值多項(xiàng)式li(τ)(i=0,…,K),并以此為基函數(shù)近似狀態(tài)變量和控制變量:

    3)動(dòng)力學(xué)方程約束轉(zhuǎn)換為代數(shù)約束.

    偽譜方法中,狀態(tài)變量由全局插值多項(xiàng)式近似,狀態(tài)變量的導(dǎo)數(shù)可通過求導(dǎo)來近似,從而將動(dòng)力學(xué)微分方程約束轉(zhuǎn)換為代數(shù)約束,即

    式中,τk為L(zhǎng)G點(diǎn).微分矩陣D∈RK×(K+1)可離線確定,將式(8)代入動(dòng)力學(xué)方程:

    于是,將最優(yōu)控制問題的動(dòng)力學(xué)微分方程約束轉(zhuǎn)換為代數(shù)約束:

    4)離散條件下的性能指標(biāo).

    將Bolza型性能指標(biāo)函數(shù)J=Φ(X0,t0,Xf,tf)+中的積分項(xiàng)用Gauss積分來近似,得到Gauss偽譜方法中的性能指標(biāo)函數(shù):

    5)離散條件下的終端狀態(tài)約束.

    最優(yōu)控制問題往往包含終端狀態(tài)約束,根據(jù)動(dòng)力學(xué)方程有:

    將終端狀態(tài)約束條件離散并用Gauss積分來近似,可得:

    6)多階段不連續(xù)問題的處理.

    對(duì)于助推-滑翔導(dǎo)彈彈道優(yōu)化這樣的分段不連續(xù)最優(yōu)控制問題,提出了多范圍偽譜技術(shù),即將問題分成若干子區(qū)間,在每個(gè)小區(qū)間上分別進(jìn)行狀態(tài)變量和控制變量的離散,在分段點(diǎn)處設(shè)置連接點(diǎn),同時(shí)加上連接點(diǎn)約束條件.例如在某一點(diǎn)τ1∈[τ0,τ2]不連續(xù),記

    則連接點(diǎn)約束條件為

    基于上述的數(shù)值近似方法,原連續(xù)最優(yōu)控制問題被離散,并轉(zhuǎn)換為非線性規(guī)劃(NLP)問題.轉(zhuǎn)換所得的NLP問題可采用數(shù)值方法求解,其中序列二次規(guī)劃(SQP)是一種較好的選擇,SQP算法較為成熟,在此不做介紹.BENSON和 HUNTINGTON等人的研究[9,10]證明了經(jīng)GPM參數(shù)化得到的NLP問題的 KKT(Karush-Kuhn-Tucker)條件與離散的Hamiltonian邊值問題(HBVP)最優(yōu)性條件之間的等價(jià)性.除此之外,GPM的主要優(yōu)點(diǎn)是算法以指數(shù)速率收斂,收斂速率快,且能在離散點(diǎn)較少時(shí)得到較好的精度.

    4 優(yōu)化算例

    1)最大射程彈道優(yōu)化.

    助推-滑翔彈道的全彈道可劃分為主動(dòng)段和滑翔彈頭飛行段,導(dǎo)彈經(jīng)主動(dòng)段彈道飛行,于關(guān)機(jī)點(diǎn)處滑翔彈頭與助推級(jí)分離,分離后的彈頭經(jīng)自由飛行段飛行,再入大氣層后依靠彈頭的高升阻比氣動(dòng)性能進(jìn)行滑翔機(jī)動(dòng)飛行,最終以某落速和落角攻擊目標(biāo).利用GPM對(duì)多階段不連續(xù)問題的處理辦法,主動(dòng)段以關(guān)機(jī)點(diǎn)機(jī)械能最大為性能指標(biāo),以主動(dòng)段關(guān)機(jī)點(diǎn)參數(shù)作為滑翔彈頭飛行段起點(diǎn)參數(shù),滑翔彈頭段以射程最大為性能指標(biāo),分段進(jìn)行彈道優(yōu)化.

    由于滑翔彈頭飛行段射程占總射程的絕大部分,全彈道射程主要取決于主動(dòng)段終點(diǎn)的運(yùn)動(dòng)參數(shù)及滑翔彈頭的氣動(dòng)特性.在助推器與滑翔彈頭參數(shù)確定的情況下,主動(dòng)段終端參數(shù)由主動(dòng)段控制規(guī)律決定.假設(shè)Hf,vf,Θf分別為主動(dòng)段結(jié)束時(shí)刻的高度、速度和彈道傾角.經(jīng)初步仿真發(fā)現(xiàn),vf主要由助推器的總沖與推力程序決定,受主動(dòng)段控制規(guī)律影響較小;Hf對(duì)總射程影響有限;Θf既對(duì)總射程影響大,又對(duì)主動(dòng)段的控制較敏感,可取Θf為主動(dòng)段彈道與滑翔彈頭飛行彈道之間的銜接量.給定參數(shù)Θf的取值,性能指標(biāo)與彈道約束不變,可得到不同Θf下的彈道曲線以及Θf與全彈道射程Lf的關(guān)系,如圖1、圖2所示.

    圖1 不同Θf對(duì)應(yīng)的彈道曲線

    圖2 全彈道射程與Θf的關(guān)系

    可見,Θf對(duì)彈道影響較大,且射程隨Θf的變化趨勢(shì)并不是簡(jiǎn)單的遞增或遞減.滑翔彈頭飛行段彈道呈波狀跳躍形式,當(dāng)Θf小于18°時(shí),全彈道曲線有4個(gè)波峰;當(dāng)Θf在18°~25°之間時(shí),彈道曲線有3個(gè)波峰;當(dāng)Θf大于25°時(shí),彈道曲線有2個(gè)波峰.波峰越多,射程越大.但在以上3種彈道形式中,分別存在3個(gè)射程的局部極大值,對(duì)應(yīng)的主動(dòng)段終端傾角為17°,23°,40°,其中Θf=17°時(shí)全彈道射程最大,可將此最優(yōu)傾角記為Θf,opt,其物理意義是保證在同樣的主動(dòng)段終點(diǎn)機(jī)械能條件下,最利于飛行器的能量利用,使射程達(dá)到最大值的最優(yōu)傾角.

    最大射程優(yōu)化結(jié)果如圖3~圖6所示,可見主動(dòng)段終點(diǎn)速度略大于4 500m/s,總射程3 300km,滑翔彈頭飛行段射程占總射程的95.5%.落角、落速均嚴(yán)格滿足端點(diǎn)約束,法向過載滿足最大過載約束.最大射程彈道呈波狀跳躍形式,大部分時(shí)間的飛行高度在30km以上,彈道較平滑.其中主動(dòng)段結(jié)束至再入大氣層(認(rèn)為大氣層邊界高度為90km)之前的彈道,認(rèn)為無控制作用,取控制量為零.

    圖3 彈道曲線

    圖4 速度曲線

    圖5 彈道傾角曲線

    圖6 控制量變化曲線

    2)起飛質(zhì)量與射程的關(guān)系.

    無論是彈道導(dǎo)彈還是助推-滑翔導(dǎo)彈,其全彈道都可分為主動(dòng)段和被動(dòng)段,其中被動(dòng)段又可分為自由飛行段、再入段兩部分.對(duì)于彈道導(dǎo)彈而言,再入段射程在整個(gè)被動(dòng)段射程所占比例很小,可近似將該段彈道看成是自由飛行段橢圓彈道的延續(xù),從而整個(gè)被動(dòng)段射程都用橢圓彈道來計(jì)算,其被動(dòng)段最大射程以及對(duì)應(yīng)的最優(yōu)傾角可估算為

    式中,vs為第一宇宙速度,R0為平均地球半徑.結(jié)合式(2)、式(15)、式(16)即可得到彈道導(dǎo)彈起飛質(zhì)量與最大射程的關(guān)系.

    對(duì)于助推-滑翔導(dǎo)彈而言,再入段射程在整個(gè)被動(dòng)段射程中占大部分,而再入段彈道形式復(fù)雜,射程很難用解析式表達(dá).當(dāng)起飛質(zhì)量已知時(shí),取助推火箭燃料比沖Ip=2 800m/s,火箭結(jié)構(gòu)比ε=0.1,重推比ν0=0.5,則可得到助推器基本參數(shù).通過前面的方法,利用彈道優(yōu)化技術(shù),即得到此起飛質(zhì)量下的最大射程與對(duì)應(yīng)的主動(dòng)段終端傾角.多次計(jì)算,即可得出起飛質(zhì)量與最大射程的關(guān)系.計(jì)算結(jié)果如圖7~圖10所示.

    圖7 不同起飛質(zhì)量下的高度-射程曲線

    圖8 不同起飛質(zhì)量下的速度曲線

    由圖7、圖8可見,起飛質(zhì)量越大則主動(dòng)段終點(diǎn)速度越大,全彈道射程也就越大.由圖9可見,在達(dá)到同等射程的情況下,助推-滑翔導(dǎo)彈的起飛質(zhì)量要遠(yuǎn)比彈道導(dǎo)彈小.由于采用單級(jí)火箭助推,在射程大于4 000km以上時(shí),助推-滑翔導(dǎo)彈的起飛質(zhì)量過大,此時(shí)宜采用多級(jí)火箭助推,以減小消極質(zhì)量.由圖10可見,與彈道導(dǎo)彈相比,助推-滑翔導(dǎo)彈最大射程對(duì)應(yīng)的主動(dòng)段終端傾角要小得多,這就要求導(dǎo)彈主動(dòng)段具有較強(qiáng)的機(jī)動(dòng)能力,以便完成轉(zhuǎn)彎.

    圖9 質(zhì)量比m0/mf隨全彈道射程變化曲線

    圖10 最優(yōu)傾角Θf,opt隨m0/mf變化曲線

    5 結(jié)論

    本文推導(dǎo)了單級(jí)固體火箭助推飛行器起飛質(zhì)量與主動(dòng)段終點(diǎn)速度解析關(guān)系式,給出了助推器參數(shù)的估算方法.使用Gauss偽譜法對(duì)助推-滑翔導(dǎo)彈的全彈道進(jìn)行了優(yōu)化計(jì)算.采用分段優(yōu)化的策略,以主動(dòng)段終端傾角Θf為彈道銜接量,分析了Θf對(duì)射程的影響,結(jié)果表明存在使射程達(dá)到最大值的最優(yōu)傾角Θf,opt,且此最優(yōu)傾角與彈道導(dǎo)彈相應(yīng)值相比要小得多.最大射程彈道的優(yōu)化結(jié)果表明,Gauss偽譜法處理此類多階段多約束的彈道優(yōu)化問題效果較好,最優(yōu)彈道起伏較小,攻角變化平滑,各項(xiàng)約束都得到滿足.利用彈道優(yōu)化計(jì)算,分析了助推-滑翔導(dǎo)彈起飛質(zhì)量與射程的關(guān)系,并與彈道導(dǎo)彈進(jìn)行了比較,結(jié)果表明助推-滑翔導(dǎo)彈在提高射程方面具有明顯優(yōu)勢(shì).

    [1]劉欣,楊濤.滑翔導(dǎo)彈再入拉起段彈道優(yōu)化設(shè)計(jì)與制導(dǎo)[J].國(guó)防科技大學(xué)學(xué)報(bào),2012,34(1):67-71.LIU Xin,YANG Tao.Trajectory optimization and guidance in reentry phase for glide missile[J].Journal of National University of Denfence Technology,2012,34(1):67-71.(in Chinese)

    [2]PERSTON H C,DARRY J P.Approximate performance of periodic hypersonic cruise trajectories for global reach[C].Proceedings of AIAA 8th International Space Planes and Hypersonics and Technologies Conference.Norfolk,VA:AIAA,1998:1-6.

    [3]李瑜,楊志紅,崔乃剛.助推-滑翔導(dǎo)彈彈道優(yōu)化研究[J].宇航學(xué)報(bào),2008,29(1):66-71.LI Yu,YANG Zhi-h(huán)ong,CUI Nai-gang.A study of optimal trajectory for boost-glide missle[J].Journal of Astronautics,2008,29(1):66-71.(in Chinese)

    [4]趙吉松,谷良賢,龔春林.高超聲速飛行器起飛質(zhì)量的解析估算[J].固體火箭技術(shù),2008,31(6):548-551.ZHAO Ji-song,GU Liang-xian,GONG Chun-lin.Analytical evaluation of takeoff mass for hypersonic vehicles[J].Journal of Solid Rocket Technology,2008,31(6):548-551.(in Chinese)

    [5]賈沛然,陳克俊,何力.遠(yuǎn)程火箭彈道學(xué)[M].長(zhǎng)沙:國(guó)防科技大學(xué)出版社,1993:89-96,195-205.JIA Pei-ran,CHEN Ke-jun,HE Li.Ballistic for long-range rocket[M].Changsha:National University of Denfence Technology Press,1993:89-96,195-205.(in chinese)

    [6]TIMOTHY R J.Common aero vehicle autonomous reentry trajectory optimization satisfying waypoint and no-fly zone constraints[D].Ohio,USA:Air Force Insitute of Technology Graduate School of Engineering and Management,2007:88-89.

    [7]雍恩米.高超聲速滑翔式再入飛行器軌跡優(yōu)化與制導(dǎo)方法研究[D].長(zhǎng)沙:國(guó)防科學(xué)技術(shù)大學(xué),2008:49-53.YONG En-mi.Study on trajectory optimization and guidance approach for hypersonic glide-reentry[D].Changsha:National University of Defense Technology,2008:49-53.(in Chinese)

    [8]楊希祥,張為華.基于Gauss偽譜法的固體運(yùn)載火箭上升段軌跡快速優(yōu)化研究[J].宇航學(xué)報(bào),2011,32(1):15-21.YANG Xi-xiang,ZHANG Wei-h(huán)ua.Rapid optimization of ascent trajectory for solid launch vehicles based on Gauss pseudospectral method[J].Journal of Astronautics.2011,32(1):15-21.(in Chinese)

    [9]HUNTINGTON G T.Advancement and analysis of a Gauss pseudospectral transcription for optimal control problems[D].Massachusetts,USA:Massachusetts Institute of Technology,2007:51-57,115-143.

    [10]BENSON D A,HUNTINGTON G T,THORVALDSEN T P.Direct trajectory optimization and costate estimation via an orthogonal collocation method[J].Journal of Guidance,Control,and Dynamics,2006,29(6):1 435-1 440.

    猜你喜歡
    彈頭滑翔射程
    斜拋射程問題的多種求解方法
    蘭花螳螂會(huì)滑翔
    熱身管下小口徑槍彈彈頭殼材料對(duì)其膛內(nèi)運(yùn)動(dòng)的影響規(guī)律
    攻天掠地的先鋒武器——滑翔導(dǎo)彈
    求解斜上拋運(yùn)動(dòng)“射高”和“射程”的兩個(gè)小妙招
    導(dǎo)彈突防概率解析與仿真計(jì)算結(jié)果差異分析
    讓你心跳加速的“滑翔飛板”
    空中滑翔大比拼(下)——滑翔傘
    射程的故事
    “彈頭”太空玩跳傘
    黑人巨大精品欧美一区二区mp4| 午夜免费成人在线视频| 久久久精品免费免费高清| 国产精品一区二区在线观看99| 国产黄色免费在线视频| 精品久久蜜臀av无| 91九色精品人成在线观看| 日韩中文字幕欧美一区二区| 国产在线观看jvid| 亚洲欧美一区二区三区黑人| 亚洲人成电影免费在线| 国产精品一区二区免费欧美 | 亚洲一码二码三码区别大吗| 久久久国产一区二区| 精品国产国语对白av| 在线观看www视频免费| 欧美黑人精品巨大| 老司机在亚洲福利影院| 日韩 欧美 亚洲 中文字幕| 男女无遮挡免费网站观看| 一本一本久久a久久精品综合妖精| 永久免费av网站大全| 99久久国产精品久久久| 正在播放国产对白刺激| 97人妻天天添夜夜摸| 自线自在国产av| 国产一区有黄有色的免费视频| 国产亚洲欧美在线一区二区| 精品国产乱码久久久久久小说| 亚洲av电影在线观看一区二区三区| 男男h啪啪无遮挡| 一本综合久久免费| 国产亚洲精品第一综合不卡| 亚洲av电影在线进入| 国产av精品麻豆| 黄片大片在线免费观看| svipshipincom国产片| 久久久久精品国产欧美久久久 | 中文字幕av电影在线播放| 国产精品免费视频内射| 97精品久久久久久久久久精品| 男女午夜视频在线观看| 无遮挡黄片免费观看| 黑丝袜美女国产一区| 亚洲熟女精品中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av片天天在线观看| 老司机靠b影院| 在线永久观看黄色视频| 丝袜人妻中文字幕| 香蕉国产在线看| 精品国产乱码久久久久久小说| 婷婷丁香在线五月| 男男h啪啪无遮挡| 亚洲专区中文字幕在线| 黄色视频,在线免费观看| 两性夫妻黄色片| 性少妇av在线| 国产精品一区二区免费欧美 | 免费在线观看完整版高清| 成人免费观看视频高清| 午夜精品久久久久久毛片777| 国产精品一区二区免费欧美 | 亚洲国产av影院在线观看| 淫妇啪啪啪对白视频 | 亚洲一码二码三码区别大吗| 久久久久精品国产欧美久久久 | 日韩三级视频一区二区三区| 丝袜美腿诱惑在线| 麻豆国产av国片精品| 中国国产av一级| 中文字幕av电影在线播放| 国产一级毛片在线| 久久久国产精品麻豆| 中文字幕人妻熟女乱码| 亚洲色图综合在线观看| 韩国高清视频一区二区三区| 考比视频在线观看| 18在线观看网站| 一本—道久久a久久精品蜜桃钙片| av又黄又爽大尺度在线免费看| 一本大道久久a久久精品| 18禁观看日本| 欧美日韩中文字幕国产精品一区二区三区 | 久久国产精品大桥未久av| 99香蕉大伊视频| 亚洲色图综合在线观看| 国产亚洲av片在线观看秒播厂| 精品国产乱码久久久久久男人| 国产区一区二久久| 亚洲午夜精品一区,二区,三区| 午夜福利视频在线观看免费| 欧美精品一区二区大全| 极品人妻少妇av视频| 国产精品一区二区在线不卡| 亚洲国产欧美日韩在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 黄片大片在线免费观看| 在线观看免费高清a一片| 悠悠久久av| 亚洲久久久国产精品| 久久九九热精品免费| 两个人看的免费小视频| 亚洲精品在线美女| 精品久久久久久久毛片微露脸 | 纯流量卡能插随身wifi吗| 男女床上黄色一级片免费看| 少妇人妻久久综合中文| 欧美国产精品va在线观看不卡| 久久国产精品男人的天堂亚洲| 亚洲七黄色美女视频| 18禁裸乳无遮挡动漫免费视频| 久久精品国产综合久久久| 新久久久久国产一级毛片| 久久人妻福利社区极品人妻图片| 亚洲精品国产一区二区精华液| 久热爱精品视频在线9| 日韩欧美一区二区三区在线观看 | 80岁老熟妇乱子伦牲交| 中文欧美无线码| 黑人欧美特级aaaaaa片| 欧美亚洲日本最大视频资源| 啦啦啦 在线观看视频| 精品熟女少妇八av免费久了| 91精品国产国语对白视频| 亚洲专区中文字幕在线| 免费观看人在逋| 久久久国产成人免费| 亚洲国产精品一区三区| 日韩制服丝袜自拍偷拍| 欧美激情高清一区二区三区| 2018国产大陆天天弄谢| av不卡在线播放| 中文字幕制服av| 90打野战视频偷拍视频| 夫妻午夜视频| 亚洲精品第二区| 80岁老熟妇乱子伦牲交| 成人黄色视频免费在线看| 日日夜夜操网爽| 国产在线一区二区三区精| 黄色视频,在线免费观看| 精品欧美一区二区三区在线| 亚洲精品国产区一区二| 午夜福利在线免费观看网站| 一级片'在线观看视频| 国产精品一区二区精品视频观看| 国产一区二区激情短视频 | 亚洲第一av免费看| 69av精品久久久久久 | 久久久久国产精品人妻一区二区| 欧美黄色片欧美黄色片| 日韩一卡2卡3卡4卡2021年| www日本在线高清视频| 精品人妻在线不人妻| 亚洲一区二区三区欧美精品| 久久国产亚洲av麻豆专区| 黑人巨大精品欧美一区二区蜜桃| 交换朋友夫妻互换小说| 9色porny在线观看| 亚洲伊人色综图| 国产亚洲精品久久久久5区| 熟女少妇亚洲综合色aaa.| 日韩中文字幕欧美一区二区| 中文字幕另类日韩欧美亚洲嫩草| 欧美中文综合在线视频| 一边摸一边做爽爽视频免费| 国产三级黄色录像| √禁漫天堂资源中文www| 99热全是精品| 欧美激情极品国产一区二区三区| 国产区一区二久久| 超碰97精品在线观看| 国产一区二区三区在线臀色熟女 | 99精品欧美一区二区三区四区| 久久久精品国产亚洲av高清涩受| 亚洲精品成人av观看孕妇| 欧美激情极品国产一区二区三区| √禁漫天堂资源中文www| 老司机靠b影院| 国内毛片毛片毛片毛片毛片| 丰满迷人的少妇在线观看| 黑人巨大精品欧美一区二区蜜桃| 18在线观看网站| 热re99久久国产66热| 一个人免费看片子| 免费人妻精品一区二区三区视频| 少妇猛男粗大的猛烈进出视频| 久久国产精品男人的天堂亚洲| 在线观看免费日韩欧美大片| 午夜日韩欧美国产| 热99久久久久精品小说推荐| 极品人妻少妇av视频| 国产精品1区2区在线观看. | 国产高清国产精品国产三级| 两性午夜刺激爽爽歪歪视频在线观看 | 99re6热这里在线精品视频| 人妻一区二区av| 国产有黄有色有爽视频| 国产成人影院久久av| 精品欧美一区二区三区在线| 亚洲欧美清纯卡通| 老熟女久久久| 老汉色av国产亚洲站长工具| 欧美亚洲 丝袜 人妻 在线| 午夜久久久在线观看| 91九色精品人成在线观看| 黑人巨大精品欧美一区二区蜜桃| 精品人妻在线不人妻| 亚洲一码二码三码区别大吗| 国产激情久久老熟女| 欧美精品啪啪一区二区三区 | 日本精品一区二区三区蜜桃| 国产黄色免费在线视频| 在线av久久热| 老司机深夜福利视频在线观看 | 老鸭窝网址在线观看| 一级毛片电影观看| 丝袜脚勾引网站| 久久精品熟女亚洲av麻豆精品| 亚洲激情五月婷婷啪啪| 9191精品国产免费久久| 午夜福利,免费看| 老司机午夜十八禁免费视频| 丝袜脚勾引网站| 国产亚洲欧美在线一区二区| 9色porny在线观看| 精品一品国产午夜福利视频| 国产成人av教育| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品久久午夜乱码| 97人妻天天添夜夜摸| 亚洲精品国产精品久久久不卡| 日韩一卡2卡3卡4卡2021年| 在线观看免费高清a一片| 精品亚洲成a人片在线观看| 久久国产精品大桥未久av| 久久国产精品人妻蜜桃| videos熟女内射| 黄色视频不卡| 国产熟女午夜一区二区三区| 国产av国产精品国产| 欧美日韩精品网址| 国产一区二区激情短视频 | 久久精品国产亚洲av香蕉五月 | 成人免费观看视频高清| 丰满人妻熟妇乱又伦精品不卡| 亚洲全国av大片| 高清黄色对白视频在线免费看| 一级毛片女人18水好多| 亚洲第一青青草原| √禁漫天堂资源中文www| 美女大奶头黄色视频| 国产成人av激情在线播放| 十八禁高潮呻吟视频| 亚洲黑人精品在线| 咕卡用的链子| 亚洲男人天堂网一区| 操美女的视频在线观看| 亚洲欧美激情在线| 欧美日本中文国产一区发布| 可以免费在线观看a视频的电影网站| 超碰97精品在线观看| 成年美女黄网站色视频大全免费| 少妇猛男粗大的猛烈进出视频| 欧美黑人精品巨大| 亚洲免费av在线视频| 日韩 亚洲 欧美在线| 成年人免费黄色播放视频| 在线av久久热| 激情视频va一区二区三区| 免费观看av网站的网址| 精品国产国语对白av| 国产欧美日韩一区二区三 | 日本av手机在线免费观看| 女人被躁到高潮嗷嗷叫费观| 欧美亚洲 丝袜 人妻 在线| 色婷婷久久久亚洲欧美| 叶爱在线成人免费视频播放| 一区二区三区四区激情视频| 国产免费av片在线观看野外av| 久久性视频一级片| 精品福利观看| 肉色欧美久久久久久久蜜桃| 色播在线永久视频| 日本黄色日本黄色录像| 青青草视频在线视频观看| h视频一区二区三区| 国产成人免费观看mmmm| 18禁观看日本| 一区在线观看完整版| 欧美日韩av久久| 香蕉丝袜av| 午夜福利视频精品| 国产在线视频一区二区| 青青草视频在线视频观看| 欧美黑人精品巨大| 免费久久久久久久精品成人欧美视频| 精品人妻一区二区三区麻豆| 三级毛片av免费| 欧美亚洲日本最大视频资源| 成年人黄色毛片网站| 久久国产精品影院| 久久久久久久久久久久大奶| 欧美日韩av久久| 黄频高清免费视频| 亚洲av成人一区二区三| 亚洲中文av在线| 久久人人爽av亚洲精品天堂| 男人添女人高潮全过程视频| 首页视频小说图片口味搜索| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 欧美成人午夜精品| 亚洲人成电影观看| 亚洲熟女毛片儿| 最近最新免费中文字幕在线| 一二三四在线观看免费中文在| 免费观看av网站的网址| 两性午夜刺激爽爽歪歪视频在线观看 | 免费av中文字幕在线| 天天影视国产精品| 亚洲人成电影免费在线| 亚洲精品国产区一区二| 国产黄色免费在线视频| 国产老妇伦熟女老妇高清| 色94色欧美一区二区| 国产成人一区二区三区免费视频网站| 免费在线观看影片大全网站| 视频在线观看一区二区三区| 高清欧美精品videossex| 国产三级黄色录像| 丰满人妻熟妇乱又伦精品不卡| 美女中出高潮动态图| 久久精品久久久久久噜噜老黄| 韩国高清视频一区二区三区| 少妇精品久久久久久久| 操出白浆在线播放| 午夜视频精品福利| 国产成人影院久久av| 精品国产国语对白av| 久久精品人人爽人人爽视色| 高清在线国产一区| 新久久久久国产一级毛片| 69精品国产乱码久久久| 欧美精品亚洲一区二区| 午夜福利在线观看吧| 亚洲第一青青草原| 男女国产视频网站| 欧美老熟妇乱子伦牲交| 日韩制服骚丝袜av| 一区二区三区激情视频| 69av精品久久久久久 | 99热网站在线观看| av在线老鸭窝| 久久精品亚洲av国产电影网| 蜜桃国产av成人99| tocl精华| 一边摸一边抽搐一进一出视频| 肉色欧美久久久久久久蜜桃| 一级毛片精品| 欧美日韩中文字幕国产精品一区二区三区 | 丁香六月欧美| 国产av一区二区精品久久| 欧美中文综合在线视频| 亚洲精品在线美女| 69av精品久久久久久 | 国产成人欧美| 久久免费观看电影| 欧美日韩精品网址| 亚洲国产精品一区二区三区在线| 久久综合国产亚洲精品| 亚洲全国av大片| 狠狠狠狠99中文字幕| 日日摸夜夜添夜夜添小说| 亚洲熟女毛片儿| 欧美日韩亚洲综合一区二区三区_| 精品第一国产精品| 男女无遮挡免费网站观看| 亚洲精品国产一区二区精华液| av免费在线观看网站| 一级片免费观看大全| 伦理电影免费视频| 亚洲国产精品999| 国产一区二区在线观看av| 成人国产一区最新在线观看| 久9热在线精品视频| 女人被躁到高潮嗷嗷叫费观| 香蕉国产在线看| 每晚都被弄得嗷嗷叫到高潮| 国产精品自产拍在线观看55亚洲 | av视频免费观看在线观看| 日本一区二区免费在线视频| 亚洲 国产 在线| 久久性视频一级片| 亚洲国产成人一精品久久久| 国产av又大| 亚洲 国产 在线| 欧美日韩亚洲综合一区二区三区_| 制服人妻中文乱码| 免费在线观看完整版高清| 91麻豆av在线| 少妇被粗大的猛进出69影院| 老司机午夜福利在线观看视频 | 九色亚洲精品在线播放| 日韩人妻精品一区2区三区| 久久 成人 亚洲| 人人妻,人人澡人人爽秒播| 久久99一区二区三区| 欧美日韩国产mv在线观看视频| 国产精品二区激情视频| 国产1区2区3区精品| 91成年电影在线观看| 亚洲精品乱久久久久久| 亚洲 国产 在线| 久久久久精品人妻al黑| 中亚洲国语对白在线视频| 热99国产精品久久久久久7| 亚洲av日韩在线播放| 亚洲精品国产av成人精品| 亚洲精品乱久久久久久| 成年女人毛片免费观看观看9 | 秋霞在线观看毛片| 欧美日韩中文字幕国产精品一区二区三区 | av欧美777| 精品一区二区三卡| 精品人妻熟女毛片av久久网站| 一个人免费在线观看的高清视频 | 精品少妇内射三级| 在线观看免费视频网站a站| 亚洲国产欧美一区二区综合| 一区二区三区乱码不卡18| 午夜福利一区二区在线看| 满18在线观看网站| 黑人欧美特级aaaaaa片| 国产精品影院久久| 国产在线视频一区二区| 最近最新中文字幕大全免费视频| 日本wwww免费看| 十分钟在线观看高清视频www| 女性被躁到高潮视频| 在线天堂中文资源库| 成年女人毛片免费观看观看9 | 国产在视频线精品| 精品亚洲成a人片在线观看| 国产成人精品无人区| 日韩中文字幕视频在线看片| 99久久99久久久精品蜜桃| 国产一区有黄有色的免费视频| 亚洲一码二码三码区别大吗| 多毛熟女@视频| 一区二区三区四区激情视频| 国产精品影院久久| 男女下面插进去视频免费观看| 久久久久久久久免费视频了| 婷婷成人精品国产| 国产男女超爽视频在线观看| 欧美国产精品一级二级三级| 午夜激情av网站| 两个人免费观看高清视频| 免费看十八禁软件| 久久ye,这里只有精品| 999久久久精品免费观看国产| 欧美激情高清一区二区三区| xxxhd国产人妻xxx| 日本五十路高清| av网站免费在线观看视频| 精品卡一卡二卡四卡免费| 亚洲视频免费观看视频| 免费av中文字幕在线| 欧美成人午夜精品| 99国产精品一区二区蜜桃av | 男女边摸边吃奶| 欧美乱码精品一区二区三区| 精品久久久久久久毛片微露脸 | 青春草亚洲视频在线观看| 人人妻,人人澡人人爽秒播| 99久久精品国产亚洲精品| 成人av一区二区三区在线看 | 国产欧美日韩一区二区三区在线| 久久久国产成人免费| 日本91视频免费播放| 99久久精品国产亚洲精品| 十八禁高潮呻吟视频| av网站在线播放免费| 亚洲国产看品久久| 性少妇av在线| 亚洲精品一区蜜桃| 免费在线观看黄色视频的| 亚洲成人免费av在线播放| 一本色道久久久久久精品综合| 亚洲全国av大片| 亚洲免费av在线视频| 老熟女久久久| 色综合欧美亚洲国产小说| 一本一本久久a久久精品综合妖精| 精品高清国产在线一区| 中文字幕人妻丝袜制服| 少妇 在线观看| 99国产精品一区二区蜜桃av | 日韩视频一区二区在线观看| 国产精品熟女久久久久浪| 精品久久久久久久毛片微露脸 | 精品国产超薄肉色丝袜足j| 成在线人永久免费视频| 大片免费播放器 马上看| 91av网站免费观看| 亚洲五月色婷婷综合| 亚洲人成电影观看| 男女午夜视频在线观看| 色综合欧美亚洲国产小说| 免费在线观看完整版高清| 黑人操中国人逼视频| 美女脱内裤让男人舔精品视频| 曰老女人黄片| 色精品久久人妻99蜜桃| 伊人久久大香线蕉亚洲五| 最黄视频免费看| 女人被躁到高潮嗷嗷叫费观| 色婷婷av一区二区三区视频| 国产欧美日韩一区二区三 | 久久久久久亚洲精品国产蜜桃av| 久久中文看片网| 国产高清videossex| 国产亚洲欧美在线一区二区| 亚洲av片天天在线观看| 国产亚洲av片在线观看秒播厂| 成人三级做爰电影| 少妇精品久久久久久久| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一二三| 欧美乱码精品一区二区三区| 黄色视频,在线免费观看| 亚洲专区字幕在线| 视频在线观看一区二区三区| 亚洲精品国产区一区二| 亚洲国产欧美在线一区| avwww免费| 久久国产亚洲av麻豆专区| 国产成人精品在线电影| 国产视频一区二区在线看| 99久久国产精品久久久| 黑人巨大精品欧美一区二区蜜桃| 伦理电影免费视频| 亚洲精品美女久久av网站| 色老头精品视频在线观看| 在线观看免费午夜福利视频| 久久久久国产精品人妻一区二区| 后天国语完整版免费观看| 欧美在线黄色| 窝窝影院91人妻| 亚洲人成77777在线视频| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| 美女午夜性视频免费| 黄色怎么调成土黄色| 最近中文字幕2019免费版| 久久人人爽人人片av| 这个男人来自地球电影免费观看| 亚洲少妇的诱惑av| 欧美中文综合在线视频| 久久久精品免费免费高清| 日日爽夜夜爽网站| 男女边摸边吃奶| 亚洲欧美一区二区三区久久| 中国美女看黄片| 丁香六月天网| 人人妻人人澡人人看| 国产xxxxx性猛交| 黄色a级毛片大全视频| 黄片小视频在线播放| 91精品伊人久久大香线蕉| 五月开心婷婷网| 亚洲成人手机| 女人久久www免费人成看片| 中文字幕av电影在线播放| 成人亚洲精品一区在线观看| 啦啦啦视频在线资源免费观看| 久久久欧美国产精品| 精品一区在线观看国产| 蜜桃在线观看..| 正在播放国产对白刺激| 久久精品aⅴ一区二区三区四区| 久久久久视频综合| 婷婷成人精品国产| av福利片在线| 久久人人爽人人片av| 亚洲av日韩精品久久久久久密| 色综合欧美亚洲国产小说| 美女午夜性视频免费| av天堂久久9| 99精品欧美一区二区三区四区| 黄色毛片三级朝国网站| 在线精品无人区一区二区三| 国产男女内射视频| 在线永久观看黄色视频| 欧美精品av麻豆av| 国产欧美日韩精品亚洲av| 午夜久久久在线观看| 黄片小视频在线播放| 国产亚洲午夜精品一区二区久久| 亚洲成人免费av在线播放| 国产成人欧美在线观看 | 嫩草影视91久久| 女警被强在线播放| 欧美激情 高清一区二区三区| 黑人操中国人逼视频| 久久久国产成人免费| 免费黄频网站在线观看国产| 窝窝影院91人妻| 超色免费av| 嫩草影视91久久| 色视频在线一区二区三区|