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

    火箭助推器翼傘回收動力學仿真與試驗分析

    2017-08-01 11:10:34蔣萬松榮偉滕海山劉濤李春
    航天返回與遙感 2017年3期
    關(guān)鍵詞:有效載荷助推器氣動力

    蔣萬松 榮偉 滕海山 劉濤 李春

    (北京空間機電研究所,北京 100094)

    火箭助推器翼傘回收動力學仿真與試驗分析

    蔣萬松 榮偉 滕海山 劉濤 李春

    (北京空間機電研究所,北京 100094)

    翼傘具有良好的滑翔性、操縱性和穩(wěn)定性,能夠解決火箭助推器落點散布大導(dǎo)致的安全性問題。為對翼傘回收系統(tǒng)和控制系統(tǒng)進行優(yōu)化設(shè)計,以基于可控翼傘回收技術(shù)的火箭助推器-控制平臺—翼傘多體飛行系統(tǒng)為研究對象,采用拉格朗日乘子法建立了三體組合10自由度多體動力學仿真模型,考慮了翼傘的表觀質(zhì)量特性和火箭助推器的氣動力影響,對某次空投飛行試驗進行了動力學過程仿真重建,通過仿真與試驗的對比分析飛行機理和系統(tǒng)性能。分析結(jié)果表明,翼傘系統(tǒng)間存在多體相對運動,必須采用多體動力學模型進行研究;機動飛行時,火箭助推器與翼傘間的相對運動角度并不大,而且火箭助推器的大尺寸對相對偏航運動的影響也并不顯著;航向跟蹤誤差主要來自操縱效率低,在小角度航向跟蹤誤差時,需提高操縱控制增益。研究成果可以為翼傘系統(tǒng)的工程設(shè)計與應(yīng)用提供參考。

    翼傘 多體動力學 拉格朗日乘子法 火箭助推器 航天回收

    0 引言

    隨著我國火箭發(fā)射活動日益頻繁,火箭助推器殘骸及廢液等的回落對一些村莊與小鎮(zhèn)居民的生命和財產(chǎn)造成不同程度的威脅。目前,火箭及其分離體的回收技術(shù)已成為國內(nèi)研究熱點。

    翼傘不僅像普通降落傘一樣折疊包裝方便、體積小和質(zhì)量輕,而且滑翔性、操縱性和穩(wěn)定性好,能夠解決火箭助推器落點散布大的問題,在航空航天飛行器和分離體的精確定點回收方面有著很高的應(yīng)用價值,但目前國內(nèi)外尚無采用翼傘系統(tǒng)回收火箭助推器的先例。

    針對火箭助推器的定點回收問題,開發(fā)了縮比質(zhì)量的翼傘回收系統(tǒng),并進行了空投飛行試驗,如圖1(a)所示。系統(tǒng)由翼傘(包括傘衣、傘繩和吊帶等)、操縱控制系統(tǒng)和控制平臺(系統(tǒng)集成結(jié)構(gòu))組成,操縱控制系統(tǒng)設(shè)備布置在控制平臺上,如圖1(b)所示,火箭助推器作為其有效載荷。在火箭發(fā)射時,翼傘系統(tǒng)折疊包裝在火箭助推器的頭錐部位。

    為優(yōu)化系統(tǒng)總體設(shè)計和歸航控制律,對某次空投試驗進行了飛行過程動力學仿真重建和機理分析。然而,翼傘系統(tǒng)與傳統(tǒng)航空飛行器不同,翼傘與有效載荷或控制平臺之間存在明顯相對運動,尤其是開傘、轉(zhuǎn)彎機動或雀降機動過程中更加明顯,要從多體系統(tǒng)動力學的角度對翼傘-載荷系統(tǒng)進行動力學建模和分析;翼傘是一種超輕結(jié)構(gòu),需要考慮翼傘的表觀質(zhì)量的影響;火箭助推器的外形對系統(tǒng)的影響尚不明確;試驗表明翼傘的操縱響應(yīng)有較大的遲滯。本文在飛行過程動力學建模時考慮了上述因素,并研究了響應(yīng)遲滯和載荷體外形氣動力對系統(tǒng)性能的影響。

    1 翼傘系統(tǒng)動力學建模

    本文主要研究火箭助推器翼傘回收系統(tǒng)從開傘充氣張滿到著陸過程的動力學行為。目前用于一般翼傘系統(tǒng)多體動力學特性研究的仿真模型均將傘體和有效載荷分別看作兩個獨立的剛體,不同之處在于連接約束模型,歸納起來主要分為三類:

    第一類是“兩體+彈簧”。兩體間通過彈簧相連,連接點的空間位置顯式求解,彈簧參數(shù)設(shè)置依賴于設(shè)計和經(jīng)驗,增加了求解自由度且難于收斂;以文獻[1-3]等為代表。

    第二類是“兩體+約束方程”。兩體間通過約束方程建立關(guān)系,保證了連接點空間位置的一致性,建模過程程式化,便于模型的擴展,但增加了求解自由度,以文獻[4-11]等為代表。

    第三類采用“共鉸點兩體+扭簧”。兩體間共用一個坐標系原點,相對姿態(tài)關(guān)系通過扭簧建立關(guān)系,其參數(shù)的設(shè)置依賴于實際系統(tǒng)和仿真經(jīng)驗,為最小解集系統(tǒng),以文獻[12-15]等為代表。

    1.1 基本假設(shè)

    本文采用擴展性好的拉格朗日乘子法對一般可控翼傘—控制平臺—載荷系統(tǒng)建立了三體10DOF(自由度)動力學仿真模型,翼傘表觀質(zhì)量計算采用Barrows方法[16-17],類似的兩體8DOF模型已得到了驗證[18]。在建立系統(tǒng)動力學仿真模型時,做如下基本假設(shè):

    1)大地是水平的,忽略地球自轉(zhuǎn),且重力加速度恒定;

    2)翼傘完成充氣張滿后幾何形狀不變,左右對稱,后緣操縱只影響氣動力;

    3)翼傘、控制平臺和火箭助推器均看作6DOF剛體,兩體間通過剛性吊帶連接約束;

    4)控制平臺和火箭助推器的質(zhì)量特性恒定不變,但翼傘含有表觀質(zhì)量,且隨大氣密度變化;

    5)火箭助推器的氣動力模型簡化為圓柱體擾流模型,忽略控制平臺的氣動力影響。

    1.2 坐標系選用

    1)坐標系

    慣性坐標系:固定于地球表面的東-北-天坐標系,坐標原點在著陸目標點地面位置,x軸指向當?shù)貣|向,y軸指向當?shù)乇毕?,z軸由右手定則確定;

    翼傘坐標系:固定翼傘上的本體坐標系,坐標原點位于翼傘等效c/4(c為平均氣動弦長),x軸在翼傘縱向?qū)ΨQ面內(nèi)與平均氣動弦線平行并指向開口方向,z軸在翼傘縱向?qū)ΨQ面上并與 x軸正交并指向上表面,y軸由右手定則確定,如圖1所示;

    控制平臺坐標系:固定在有效載荷平臺上的本體坐標系,坐標原點位于有效載荷平臺的質(zhì)心位置,x軸與穩(wěn)定飛行前方一致,z軸與x軸垂直并指向上方,y軸由右手定則確定;

    有效載荷坐標系:固定在火箭助推器上的本體坐標系,坐標原點位于火箭助推器的質(zhì)心位置,x軸指向頭錐,z軸與x軸垂直并處于吊點一側(cè),y軸由右手定則確定,如圖1所示;

    翼傘氣流坐標系:坐標原點位于翼傘等效c/4處,x軸指向氣流的來流方向,z軸在翼傘的縱向?qū)ΨQ面內(nèi)且與x軸垂直,y軸由右手定則確定。

    2)坐標系間關(guān)系

    從慣性坐標系到本體坐標系采用z-y-x轉(zhuǎn)序的歐拉角(偏航、俯仰和滾轉(zhuǎn))來描述姿態(tài),記為(ψ, θ, φ),通過下角標區(qū)分(翼傘為 p;控制平臺為 C;有效載荷為 B),慣性系到本體系的轉(zhuǎn)換矩陣定義[19]為S =Sx(φ) Sy(θ)Sz(ψ),且有 S-1=ST其( )T為轉(zhuǎn)置矩陣。其中,

    從氣流坐標系到翼傘本體坐標系采用z-y-x轉(zhuǎn)序的攻角α和側(cè)滑角β來描述,氣流系到本體系的轉(zhuǎn)換矩陣定義[19]為 Sa=Sy(- α) Sz(β)。

    1.3 系統(tǒng)動力學方程

    翼傘系統(tǒng)中,翼傘、控制平臺和火箭助推器各看作6DOF剛體,多體系統(tǒng)動力學方程為

    式中 q為廣義坐標;˙q為廣義速度;λ為拉格朗日乘子;Q為廣義力;f為約束方程向量;T為系統(tǒng)動能,表示為

    式中 M為質(zhì)量矩陣;J為慣量矩陣;H為慣量交叉矩陣;S本體系到慣性系的坐標轉(zhuǎn)換矩陣;ω為本體角速度;ix為慣性系下位置;i˙x為慣性系下速度。

    根據(jù)方程(1),三體系統(tǒng)的動力學方程可以表示為

    即式(3)簡寫為

    式中 轉(zhuǎn)換矩陣 STT=diag(Sp,E, SC,E, SB,E),為氣動力;L?為氣動力矩;E為單位矩陣;V為本體系下平動速度,

    下面給出方程(4)的約束方程及拉格朗日乘子的求解方法。翼傘系統(tǒng)多體間的約束模型可以看作一個復(fù)合約束副,如圖2所示。其中,A1和A2代表翼傘連接點;A0為A1A2的中點,在翼傘縱向?qū)ΨQ面內(nèi);B1和B2代表控制平臺的2個上連接點;B0為B1B2的中點,在控制平臺縱向?qū)ΨQ面內(nèi);C1,C2,C3和C4為控制平臺4個下連接點;D1和D2為助推器上兩個連接點。

    1)翼傘—控制平臺的連接特征線A0B0與B1B2(平行于控制平臺y軸)正交,即約束f1

    式中 ε=(- s in θ,0,-cosθ )T; θ為翼傘安裝角; e ( j=1,2,3)為單位矢量基,滿足 E =(e, e, e )。

    2)翼傘—控制平臺連接特征線A0B0總是在翼傘縱向?qū)ΨQ面內(nèi),但沿A0B0線可移動,即f2和f3分別為

    式中 rB0為B0點在控制平臺本體坐標系下的位置向量。

    3)翼傘—控制平臺間沿A0B0線移動受單個吊帶長度約束,假設(shè)其不發(fā)生變化,即f4表示為

    式中 rA1和 rB1分別為連接吊索兩端點在翼傘和控制平臺坐標系下的向量;l1為連接吊索長度。

    4)控制平臺與有效載荷間約束與式(7)和式(8)類似,即f5~ f8表示為

    式中1Dr ,2Dr ,2Cr 和3Cr 為連接吊索兩端點在有效載荷和控制平臺本體坐標系下向量;l2和 l3為控制平臺與有效載荷(助推器)間連接吊索長度,如圖2所示。

    為得到方程(1)第三項,并直接求解出拉格朗日乘子,對約束方程f求導(dǎo)

    式中 拉格朗日方程廣義坐標 q =(xip,θp,xiC,θC,xiB,θB)T;=(Vp, ωp, VC,ωC, VB,ωB)T為本體系下表達的廣義速度; f=(f1, f2,f3, f4,f5, f6, f7,f8)T為約束方程;N為約束方程對本體系廣義速度的導(dǎo)數(shù)矩陣;K為余項列陣。

    式中 S =STS; S =STS ; S˙ =S ω×-ω×S ; S˙ =S ω×-ω×S ;( )×為反對稱矩陣。1Cp2 BC 11pC122CB2ρ x S r x S r ;Br和Cr分別為ρ x S r x S r;121CCpp= + --i BiA11ρ x S r x S r ;232 BBCC= + --i DiC3BBCC= + --i DiC單根吊帶兩端點在載荷和控制平臺本體坐標系下的位置向量;

    應(yīng)用式(4)和式(13),可求解出拉格朗日乘子

    從而最終得到的系統(tǒng)動力學方程為

    以上方程可通過通常數(shù)值積分方法進行求解。

    1.4 表觀質(zhì)量估算

    本文所述翼傘沿展向展開后為矩形的圓弧形翼傘,采用Barrows估算方法[16]給出的翼傘的表觀質(zhì)量矩陣如式(18)所示,計算時加到方程(3)中與翼傘相關(guān)質(zhì)量項中。

    =0.055;*h為翼傘比高度(翼傘頂點到翼傘兩側(cè)邊緣的高度與展長之比);b為翼傘展長;c為翼傘弦長;w為傘繩匯聚點;0e為翼型剖面百分比厚度;下標k為俯仰中心;下標r為滾動中心;γ為展弦比;kA和kB為系數(shù);aρ為大氣密度;R0為系統(tǒng)長度;SA=bc為翼傘參考面積。

    1.5 氣動模型

    翼傘氣動力計算采用經(jīng)典飛行力學方法,計算式及相關(guān)氣動系數(shù)如表1所示。升阻力氣動系數(shù)參考了 lingard[20]文獻,靜導(dǎo)數(shù)系數(shù)和橫向氣動參數(shù)參考了 Prakash[13]和 Iacomini[21]文獻,瞬時飛行狀態(tài)大氣密度(GJB 365.1-1987)和氣動參數(shù)均通過樣條插值計算。

    表1 翼傘氣動力模型Tab.1 Aerodynamic model for parafoil

    在考慮有效載荷氣動力時,將火箭助推器簡化為圓柱體擾流模型,不考慮兩端圓形端面的氣動力,圓柱面受到氣動阻力和氣動升力作用,如圖2所示,采用經(jīng)典飛行力學計算方法,主要與來流總攻角有關(guān),阻力系數(shù)和升力系數(shù)如圖3所示[22]。

    1.6 歸航控制律說明

    翼傘系統(tǒng)完整的飛行過程包括定向飛行段、盤旋管理段、逆風調(diào)整段和雀降機動段。定向飛行段根據(jù)翼傘系統(tǒng)的當前位置及目標落點位置,自動調(diào)整航向,使系統(tǒng)向著目標點飛行;盤旋管理段在落點的目標上空進行盤旋機動,對飛行高度進行管理,使之在落點上空附近飛行;逆風調(diào)整段根據(jù)當前飛行高度、位置、風向等信息,計算欲落點位置,并在該段結(jié)束時使系統(tǒng)飛至預(yù)落點位置,且飛行方向處于逆風飛行狀態(tài);雀降機動段通過雙邊操縱制動實施精確軟著陸。

    本文采用具有基準下偏量的簡單線性下偏操縱歸航控制律,即翼傘后緣下偏量輸出量指令與航向偏差量 Δψ 成正比關(guān)系,如表2所示。其中h(t)和R(t)分別為系統(tǒng)在任意時刻t,離目標點高度和水平距離,h1、h2、R1和R2為模式轉(zhuǎn)換設(shè)置的定值,下標m和n分別表示翼傘左、右后緣。

    根據(jù)事先預(yù)置或地面通過遙控上傳的風向數(shù)據(jù)(優(yōu)先級高于計算值),計算實際航向與風向的夾角Δψ 。飛行試驗中,逆風調(diào)整段和雀降著陸段因地著陸點高度偏差沒有實施。

    表2 各飛行模式歸航控制算法Tab.2 Homing control algorithm in all types of flight modes

    1.7 操縱響應(yīng)延遲模型

    翼傘操縱機構(gòu)對操縱控制指令執(zhí)行需要一定的響應(yīng)時間,大型翼傘更加明顯,長達數(shù)秒。動力學建模時有必要考慮操縱響應(yīng)的延遲,建立與操縱率和操縱指令有關(guān)的簡單模型

    式中max˙ 為最大操縱速率;Δt為步長;cδ 根據(jù)表2計算得到。

    2 動力學仿真

    2.1 仿真條件說明

    翼型Clark-Y,弦長5.55m,翼展14.4m,厚度0.15m,系統(tǒng)長度11.5m,安裝角6°(不含剖面4°),翼傘凈質(zhì)量35kg;助推器模型長5.6m,直徑0.9m,質(zhì)量800kg,重心距噴口1.8m;航向定義為東偏北為正,角度范圍為(-180°,180°)。

    仿真初始條件:東向477m,北向1 002m,高度3 980m,航向-160°,飛行時間258s,盤旋半徑550m,無風。著陸目標點(0,0,1 575)。

    2.2 仿真與試驗結(jié)果對比

    仿真與試驗結(jié)果的對比如圖4所示。

    2.3 結(jié)果分析

    通過以上仿真與試驗結(jié)果的對比,經(jīng)分析得到以下基本結(jié)論:

    從圖 4(a)和(b)可以看出,翼傘系統(tǒng)整體上至少經(jīng)歷了定向飛行和盤旋管理段,是經(jīng)過飛行模式轉(zhuǎn)變的復(fù)合運動。從圖4(c)看出,翼傘系統(tǒng)的下降速度基本穩(wěn)定在一個恒值,約9.5m/s;仿真得到的水平速度與試驗結(jié)果有所差別,仿真結(jié)果是逐漸減小的,約20~23m/s,而試驗測得水平速度有波動,系統(tǒng)滑翔比為2.2左右;結(jié)合圖4(d),水平速度與操縱量存在關(guān)聯(lián)性。圖4(d)中仿真得到的操縱量比試驗記錄的操縱量略小,試驗曲線沒有體現(xiàn)定向飛行到盤旋管理段的轉(zhuǎn)換過程,而仿真結(jié)果在60~70s左右出現(xiàn)了零操縱狀態(tài)到50%操縱量峰值突變,明顯體現(xiàn)了該轉(zhuǎn)換過程。

    圖5給出了仿真分析結(jié)果,包括氣流角、姿態(tài)角和相對姿態(tài)角。圖示均表明了在60~70s出現(xiàn)了飛行模式轉(zhuǎn)變,攻角、側(cè)滑角、俯仰角、滾轉(zhuǎn)角及相對俯仰角和滾轉(zhuǎn)角均有所變化,在進入盤旋管理段后趨于穩(wěn)定。圖 5(c)和(d)中,翼傘與平臺間存在相對運動,在機動時表現(xiàn)得比較明顯,說明系統(tǒng)間存在相對運動,但是角度比較小,進一步說明了采用多體動力學仿真模型的必要性,也說明至少應(yīng)采用 7或8自由度以上的模型才能較真實反映物理過程。從圖5(c)初步得到,火箭助推器的大尺寸對相對偏航運動的影響并不大。

    結(jié)合歸航控制律和飛行過程,定向飛行階段存在大角度的跟蹤誤差,始終沒有朝向目標點(原點),以致于定向飛行到盤旋管理段的過度過程比較平緩,但按照控制律邏輯應(yīng)有明顯的界限,原因可能是轉(zhuǎn)彎操縱控制效率低、較大的相對偏航運動或航向測量誤差。從仿真結(jié)果來看,翼傘與控制平臺間的相對偏航運動較小(<4°);航向測量誤差不超過±10°,都不是航向跟蹤誤差的主因,因此航向大角度跟蹤誤差更可能是轉(zhuǎn)彎操縱控制效率低導(dǎo)致,這就要求在小角度航向誤差時,增大操縱控制增益,但同時需改進航向飛行到盤旋管理段的過度過程控制律。

    3 結(jié)束語

    本文針對火箭助推器翼傘回收多體飛行動力學系統(tǒng)為研究對象,采用拉格朗日乘子法建立了三體10DOF動力學仿真模型,考慮了翼傘的表觀質(zhì)量特性。對某次空投飛行試驗過程進行了動力學仿真重建,驗證了仿真模型的有效性,分析了飛行過程機理和性能,提出小角度操縱時應(yīng)增大操縱控制增益。

    References)

    [1]VISHNYAK A. Simulation of the Payload-parachute-wing System Flight Dynamics[C]//Aerospace Design Conference, Irvine, CA, U.S.A. AIAA 93-1250.

    [2]熊菁. 翼傘系統(tǒng)動力學與歸航方案研究[D]. 長沙: 國防科技大學博士學位論文, 2005. XIONG Jing. Research on Dynamics and Homming Control of Parafoil System[D]. Changsha: National University of DefenseTechnology, 2005. (in Chinese)

    [3]熊菁, 宋旭民, 秦子增. 翼傘系統(tǒng)兩體相對運動分析[J]. 航天返回與遙感, 2004, 25(2): 10-16. XIONG Jing, SONG Xuming, QIN Zizeng. Study on Relative Motion of Two-body Parafoil System[J]. Spacecraft Recovery & Remote Sensing, 2004, 25(2): 10-16. (in Chinese)

    [4]YAKIMENKO O A. On the Development of a Scalable 8-DOF Model for a Generic Parafoil-payload Delivery System[C]// 18th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Munich, Germany. AIAA 2005-1665.

    [5]PILLASCH D W, SHEN Y C, VALERO N. Parachute/Submunition System Coupled Dynamics[C]//8th Aerodynamic Decelerator and Balloon Technology Conference, Hyannis, MA, U.S.A. AIAA 84-0784.

    [6]WISE K A. Dynamics of a UAV with Parafoil under Powered Flight[C]//AIAA Guidance, Navigation, and Control Conference and Exhibit, Keystone, Colorado. AIAA 2006-6795.

    [7]STRICKER G, WITTE L. Analysis of the Relative Motion in a Parafoil-load System[C]//16th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Boston, MA, U.S.A. AIAA 2001-2013.

    [8]CHRISTIAAN R. A Flight Simulation Algorithm for a Parafoil Suspending an Air Vehicle[J]. Journal of Guidance, Control and Dynamics, 2007, 30(3): 791-803.

    [9]MOOIJ E, WIJNANDS Q G, et al. 9 DOF Parafoil/ Payload Simulator Development and Validation[C]//AIAA Modeling and Simulation Technologies Conference and Exhibit, Austin, Texas. AIAA 2003-5459.

    [10]陳建平, 寧雷鳴, 張紅英, 等. 基于多體動力學的大型翼傘系統(tǒng)飛行仿真分析[J]. 飛行力學, 2015, 33(6): 486-490. CHEN Jianping, NING Leiming, ZHANG Hongying, et al. Flight Simulation and Analysis of Large Parafoil System Based on Multibody Dynamics[J]. Flight Dynamics, 2015, 33(6): 486-490. (in Chinese)

    [11]陳建平, 張紅英, 童明波, 等. 翼傘系統(tǒng)縱向飛行性能仿真[J]. 中國空間科學技術(shù), 2015, 35(2): 25-32. CHEN Jianping, ZHANG Hongying, TONG Mingbo, et al. Longitudinal Flight Performances Simulation for Parafoil System[J]. Chinese Space Science and Technology, 2015, 35(2): 25-32. (in Chinese)

    [12]BARROWS T M. Multibody Parafoil Model[C]//20th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Seattle, Washington. AIAA 2009-2945.

    [13]PRAKASH O, ANANTHKRISHNAN N. Modeling and Simulation of 9-DOF Parafoil-payload System Flight Dynamics[C]// AIAA Atmospheric Flight Mechanics Conference and Exhibit, Keystone, Colorado. AIAA 2006-6130.

    [14]MULLER S, WAGNER O, SACHS G. A High-fidelity Nonlinear Multibody Simulation Model for Parafoil System[C]//17th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Monterey, California. AIAA 2003-2120.

    [15]GORMAN C M, SLEGERS N J. Comparison and Analysis of Multi-body Parafoil Models with Varying Degrees of Freedom[C]//21st AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Dublin, Ireland. AIAA 2011-2615.

    [16]BARROWS T M. Apparent Mass of Parafoils with Spanwise Camber[C]//16th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Boston, MA, U.S.A. AIAA 2001-2006.

    [17]LISSAMAN B S, BROWN G J. Apparent Mass Effects on Parafoil Dynamics[C]//Aerospace Design Conference, Irvine, CA, U.S.A. AIAA 93-1236.

    [18]蔣萬松, 榮偉, 滕海山, 等. 翼傘載荷系統(tǒng)多體動力學仿真分析[J]. 南京航空航天大學學報, 2016, 48(4): 474-481. JIANG Wansong, RONG Wei, TENG Haishan, et al. Multibody Dynamical Simulation Analysis for Parafoil-payload System[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2016, 48(4): 474-481. (in Chinese)

    [19]趙育善, 師鵬. 航天器飛行動力學建模理論與方法[M]. 北京: 北京航空航天大學出版社, 2012. ZHAO Yushan, SHI Peng. Spacecraft Flight Dynamics Modeling Theory and Methods[M]. Beijing: Beijing University of Aeronautics & Astronautics Press, 2012. (in Chinese)

    [20]LINGARD J S. Ram Airparachute Design[C]//13th AIAA Aerodynamic Decelerator Systems Technology Conference, Clearwater Beach. AIAA 1995.

    [21]IACOMINI C S, CERIMELE C J. Lateral-directional Aerodynamics from a Large Scale Parafoil Test Program[C]//15th Aerodynamic Decelerator Systems Technology Conference, Toulouse, France. AIAA 99-1731.

    [22]夏雪湔, 鄧學鎣. 工程分離流動力學[M]. 北京: 北京航空航天大學出版社, 1991. XIA Xuejian, DENG Xueying. Engineering Separation Flow Dynamics[M]. Beijing: Beijing University of Aeronautics & Astronautics Press, 1991. (in Chinese)

    Dynamical Simulation and Test Analysis for Booster Recovery with Parafoil System

    JIANG Wansong RONG Wei TENG Haishan LIU Tao LI Chun
    (Beijing Institute of Space Mechanics and Electricity, Beijing 100094, China)

    Controllable parafoil system could solve the security problem resulted from large dispersion of deserted boosters with its favorable gliding, manueverability and stability. To optimize the parafoil recovery system and its control system, taking a multibody dynamic system of booster-controller-parafoil based on parafoil recovery technology as a researching object, a three-body simulation model with 10-degree-of-freedom is established with Lagrange multiplier method, which considers the apparent mass characteristics of parafoil and aerodynamic impact of the booster. The process of an airdrop test is reconstructed with this model, and flight mechanism and performances are analyzed. The results show that multibody relative motion does exist among parafoil-payload systems, thus it is essential to establish a mutltibody model for simulation. While in maneuvering, the relative attitudes are not so large, and the aerodynamic of large size of rocket on parafoil system is not so profound. The tracking errors mainly come from low steering efficiency of parafoil, thus it is necessary to increase the gain of steering. The results provide a reference for the parafoil system engineering and application.

    parafoil; multibody dynamics; Lagrange multiplier method; booster; spacecraft recovery

    V275

    A

    1009-8518(2017)03-0013-11

    10.3969/j.issn.1009-8518.2017.03.002

    蔣萬松,男,1981年生,2010年獲中國空間技術(shù)研究院飛行器設(shè)計專業(yè)碩士學位,高級工程師,中國空間技術(shù)研究院博士在讀。研究方向為航天返回與著陸技術(shù)。E-mail: skylingy24@sina.com。

    (編輯:劉穎)

    2017-03-20

    猜你喜歡
    有效載荷助推器氣動力
    理念牽引 機制創(chuàng)新 人才驅(qū)動 做有效載荷創(chuàng)新發(fā)展領(lǐng)跑者
    助推器殘骸被成功回收
    美國SLS重型運載火箭助推器測試
    軍事文摘(2021年16期)2021-11-05 08:49:20
    面向有效載荷數(shù)字化研制的標準化工作轉(zhuǎn)型初探
    飛行載荷外部氣動力的二次規(guī)劃等效映射方法
    衛(wèi)星有效載荷研制流程的策劃與推進
    檢驗——提高分數(shù)的助推器
    導(dǎo)讀案,英語啟發(fā)式閱讀教學的“助推器”
    側(cè)風對拍動翅氣動力的影響
    高速鐵路接觸線覆冰后氣動力特性的風洞試驗研究
    最近在线观看免费完整版| 亚洲一区二区三区不卡视频| 亚洲av中文字字幕乱码综合| 国产精品亚洲美女久久久| 亚洲欧美日韩东京热| 精品不卡国产一区二区三区| 99久久久亚洲精品蜜臀av| 天天躁夜夜躁狠狠躁躁| 男人舔女人下体高潮全视频| 母亲3免费完整高清在线观看| 亚洲免费av在线视频| 老司机午夜福利在线观看视频| 在线国产一区二区在线| 久久人人精品亚洲av| 国产亚洲欧美在线一区二区| 黄色片一级片一级黄色片| 啦啦啦观看免费观看视频高清| 欧美最黄视频在线播放免费| 97碰自拍视频| 在线观看美女被高潮喷水网站 | 午夜久久久久精精品| 日本熟妇午夜| 亚洲美女视频黄频| 一本久久中文字幕| 日韩大尺度精品在线看网址| 男人舔女人的私密视频| 精品久久久久久久毛片微露脸| 亚洲人成电影免费在线| 天堂影院成人在线观看| 久久精品91蜜桃| 国产欧美日韩一区二区三| 大型av网站在线播放| 可以免费在线观看a视频的电影网站| 国产亚洲精品综合一区在线观看 | 成年女人毛片免费观看观看9| 久久精品影院6| 精品国产乱码久久久久久男人| 狂野欧美白嫩少妇大欣赏| 脱女人内裤的视频| 国产精品久久久久久亚洲av鲁大| 日韩欧美免费精品| 美女扒开内裤让男人捅视频| 欧美日韩精品网址| 18禁观看日本| 伊人久久大香线蕉亚洲五| 身体一侧抽搐| 伊人久久大香线蕉亚洲五| 国产成人系列免费观看| 精品人妻1区二区| 亚洲国产精品999在线| 老司机在亚洲福利影院| 夜夜看夜夜爽夜夜摸| 身体一侧抽搐| 欧美人与性动交α欧美精品济南到| av在线天堂中文字幕| 久久久久久久午夜电影| 日韩欧美免费精品| 桃红色精品国产亚洲av| 久久久久久亚洲精品国产蜜桃av| 国产片内射在线| 色综合婷婷激情| 国产精品一区二区精品视频观看| 在线观看www视频免费| 国产野战对白在线观看| 日韩 欧美 亚洲 中文字幕| 男女床上黄色一级片免费看| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色| 国产1区2区3区精品| 国产一区二区三区视频了| 天堂动漫精品| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲精品综合一区在线观看 | 成年免费大片在线观看| 不卡av一区二区三区| 啦啦啦免费观看视频1| 大型av网站在线播放| www.精华液| 麻豆成人av在线观看| 国产精品久久久久久精品电影| 1024视频免费在线观看| 窝窝影院91人妻| 久久精品人妻少妇| 黄频高清免费视频| 国产高清videossex| 精品不卡国产一区二区三区| 岛国视频午夜一区免费看| 禁无遮挡网站| 一区二区三区国产精品乱码| www日本黄色视频网| 日本a在线网址| 精品久久久久久成人av| 在线观看www视频免费| 亚洲精品美女久久av网站| 在线观看日韩欧美| 99国产精品99久久久久| 99国产极品粉嫩在线观看| 国内揄拍国产精品人妻在线| 精品久久蜜臀av无| 黄色丝袜av网址大全| 精品国内亚洲2022精品成人| 男女那种视频在线观看| 精品第一国产精品| 精品国产乱子伦一区二区三区| 午夜日韩欧美国产| 国产精品亚洲一级av第二区| 久久人人精品亚洲av| 老汉色av国产亚洲站长工具| 岛国在线免费视频观看| 黄色片一级片一级黄色片| 中文字幕最新亚洲高清| netflix在线观看网站| 国产99白浆流出| 国产精品香港三级国产av潘金莲| 长腿黑丝高跟| 国产精品自产拍在线观看55亚洲| 老司机午夜福利在线观看视频| 国产精品一区二区三区四区免费观看 | 桃红色精品国产亚洲av| 日日爽夜夜爽网站| 久久香蕉精品热| svipshipincom国产片| 国产成人精品无人区| 久9热在线精品视频| 亚洲国产中文字幕在线视频| 国产1区2区3区精品| 变态另类成人亚洲欧美熟女| 久久久国产欧美日韩av| 久久精品影院6| 搡老妇女老女人老熟妇| 免费观看人在逋| 日韩大码丰满熟妇| 精品欧美一区二区三区在线| 看黄色毛片网站| 三级国产精品欧美在线观看 | 亚洲美女黄片视频| 国产成年人精品一区二区| 母亲3免费完整高清在线观看| 欧美黄色淫秽网站| 窝窝影院91人妻| 亚洲熟女毛片儿| 午夜福利欧美成人| 国产伦人伦偷精品视频| 欧美久久黑人一区二区| 日韩欧美 国产精品| 久久久精品欧美日韩精品| 精品久久蜜臀av无| ponron亚洲| 日韩欧美免费精品| 欧美一级a爱片免费观看看 | 精品久久久久久久久久久久久| 一区二区三区激情视频| 精品国产乱码久久久久久男人| 亚洲国产精品久久男人天堂| 久久精品人妻少妇| 18禁裸乳无遮挡免费网站照片| 99精品在免费线老司机午夜| 国产主播在线观看一区二区| 成人手机av| 成人精品一区二区免费| 国产精品久久久久久精品电影| 久9热在线精品视频| 国产成+人综合+亚洲专区| 国产蜜桃级精品一区二区三区| 久久香蕉激情| 免费无遮挡裸体视频| 深夜精品福利| 欧美日韩亚洲国产一区二区在线观看| 成人av一区二区三区在线看| 人成视频在线观看免费观看| 中文字幕av在线有码专区| 天天躁狠狠躁夜夜躁狠狠躁| 一级作爱视频免费观看| 91九色精品人成在线观看| av欧美777| 久久久久久免费高清国产稀缺| 亚洲精品国产精品久久久不卡| 性欧美人与动物交配| 三级男女做爰猛烈吃奶摸视频| 亚洲熟妇熟女久久| 少妇的丰满在线观看| 欧美三级亚洲精品| 一夜夜www| 男女那种视频在线观看| 美女免费视频网站| 操出白浆在线播放| 国产视频一区二区在线看| 视频区欧美日本亚洲| 亚洲av电影不卡..在线观看| www.自偷自拍.com| 国产精品 国内视频| 丁香六月欧美| 国产私拍福利视频在线观看| 99久久综合精品五月天人人| 欧美+亚洲+日韩+国产| 国产成人精品久久二区二区91| 高潮久久久久久久久久久不卡| 日韩欧美精品v在线| www.www免费av| 欧美高清成人免费视频www| 国产一区在线观看成人免费| 午夜精品久久久久久毛片777| 搞女人的毛片| 久久久久久久午夜电影| 男女那种视频在线观看| 国产主播在线观看一区二区| 欧美黄色淫秽网站| 精品国产超薄肉色丝袜足j| 香蕉丝袜av| 日本一本二区三区精品| 精品国产超薄肉色丝袜足j| 他把我摸到了高潮在线观看| 成在线人永久免费视频| 舔av片在线| 国产69精品久久久久777片 | 欧美不卡视频在线免费观看 | aaaaa片日本免费| 九色国产91popny在线| 香蕉av资源在线| 亚洲aⅴ乱码一区二区在线播放 | 夜夜看夜夜爽夜夜摸| 美女午夜性视频免费| 日韩av在线大香蕉| 欧美日韩亚洲国产一区二区在线观看| 岛国在线观看网站| 五月伊人婷婷丁香| 日韩欧美免费精品| 亚洲成人久久爱视频| 精品一区二区三区av网在线观看| 一级黄色大片毛片| 日韩精品免费视频一区二区三区| 欧美成人性av电影在线观看| 欧美日韩亚洲国产一区二区在线观看| 久久久久国产精品人妻aⅴ院| 亚洲国产看品久久| 欧美在线一区亚洲| 国产精品 欧美亚洲| 日本成人三级电影网站| 色综合亚洲欧美另类图片| 国产成人影院久久av| 日韩欧美国产在线观看| 亚洲精品在线美女| 波多野结衣高清作品| 久久草成人影院| 两个人的视频大全免费| 一本久久中文字幕| 国产精品免费视频内射| 两个人视频免费观看高清| 看片在线看免费视频| 两个人免费观看高清视频| x7x7x7水蜜桃| 999久久久国产精品视频| 亚洲精华国产精华精| 欧美乱码精品一区二区三区| 天堂动漫精品| 观看免费一级毛片| 久久精品国产亚洲av香蕉五月| 国产成人一区二区三区免费视频网站| 老司机靠b影院| a级毛片在线看网站| 国产熟女xx| 长腿黑丝高跟| 国产在线精品亚洲第一网站| 少妇被粗大的猛进出69影院| 成人亚洲精品av一区二区| 国产精品久久久久久亚洲av鲁大| 在线观看免费午夜福利视频| 久久久久久免费高清国产稀缺| 大型av网站在线播放| 9191精品国产免费久久| 久久香蕉国产精品| 久久久久精品国产欧美久久久| 色在线成人网| 国内精品久久久久久久电影| 天天躁夜夜躁狠狠躁躁| av福利片在线| 人妻夜夜爽99麻豆av| 老鸭窝网址在线观看| 国产亚洲精品一区二区www| 老汉色av国产亚洲站长工具| 黑人欧美特级aaaaaa片| 久久天堂一区二区三区四区| 国产精品 欧美亚洲| 成人三级黄色视频| 亚洲免费av在线视频| ponron亚洲| 国产亚洲av嫩草精品影院| 国产精品一区二区三区四区免费观看 | 国产精品1区2区在线观看.| 亚洲精品色激情综合| www国产在线视频色| 女同久久另类99精品国产91| 国产黄a三级三级三级人| 香蕉国产在线看| 亚洲成av人片在线播放无| 天天躁狠狠躁夜夜躁狠狠躁| 岛国在线免费视频观看| 人妻丰满熟妇av一区二区三区| 国产三级中文精品| 欧美成人性av电影在线观看| xxx96com| 国产精品久久久久久亚洲av鲁大| 国产高清有码在线观看视频 | 国产精品国产高清国产av| 色噜噜av男人的天堂激情| 免费看十八禁软件| 50天的宝宝边吃奶边哭怎么回事| 久久久久久九九精品二区国产 | 狂野欧美激情性xxxx| 久久精品亚洲精品国产色婷小说| 人妻夜夜爽99麻豆av| 午夜免费观看网址| 真人一进一出gif抽搐免费| 男女视频在线观看网站免费 | 国内少妇人妻偷人精品xxx网站 | 嫁个100分男人电影在线观看| 亚洲人成网站高清观看| 婷婷精品国产亚洲av在线| 欧美zozozo另类| 亚洲 欧美一区二区三区| 亚洲片人在线观看| 国产麻豆成人av免费视频| 成人高潮视频无遮挡免费网站| 亚洲av美国av| 中出人妻视频一区二区| 一区二区三区高清视频在线| 嫩草影院精品99| 日本在线视频免费播放| 一区二区三区高清视频在线| 黄频高清免费视频| 黄色a级毛片大全视频| 97人妻精品一区二区三区麻豆| 欧美人与性动交α欧美精品济南到| www日本在线高清视频| 床上黄色一级片| 高清在线国产一区| 亚洲精品在线观看二区| 日韩大码丰满熟妇| 好看av亚洲va欧美ⅴa在| 99久久精品热视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美日韩东京热| 久久天堂一区二区三区四区| 一进一出抽搐动态| 国产精品98久久久久久宅男小说| 久久久久久久午夜电影| svipshipincom国产片| 丁香欧美五月| 嫁个100分男人电影在线观看| 欧美乱色亚洲激情| 黄片小视频在线播放| 欧美另类亚洲清纯唯美| 国产精品99久久99久久久不卡| 国产精品一及| 特大巨黑吊av在线直播| 日韩中文字幕欧美一区二区| 极品教师在线免费播放| 日韩欧美在线二视频| 午夜视频精品福利| 最近最新免费中文字幕在线| 国产精品亚洲美女久久久| 99精品欧美一区二区三区四区| 成人三级做爰电影| 动漫黄色视频在线观看| 搡老妇女老女人老熟妇| 日日夜夜操网爽| 99re在线观看精品视频| 可以在线观看的亚洲视频| 久久香蕉精品热| 亚洲av熟女| 国产精品亚洲美女久久久| 亚洲精品美女久久av网站| 精品久久久久久成人av| 人人妻人人澡欧美一区二区| 国产欧美日韩一区二区精品| 欧美性长视频在线观看| 精品久久久久久久毛片微露脸| 久久伊人香网站| 日韩免费av在线播放| 黄片大片在线免费观看| 一二三四在线观看免费中文在| 亚洲中文av在线| 老司机午夜十八禁免费视频| 国产1区2区3区精品| 国产av一区二区精品久久| 一级作爱视频免费观看| 精品第一国产精品| 两个人免费观看高清视频| 99久久国产精品久久久| 婷婷亚洲欧美| 国产精品日韩av在线免费观看| 久久人妻av系列| 每晚都被弄得嗷嗷叫到高潮| 日韩大码丰满熟妇| 日本熟妇午夜| 久久久国产精品麻豆| 国产精品永久免费网站| 好男人电影高清在线观看| 欧美乱妇无乱码| 亚洲国产精品合色在线| 琪琪午夜伦伦电影理论片6080| 中出人妻视频一区二区| 久久久久久人人人人人| 母亲3免费完整高清在线观看| 搞女人的毛片| 成在线人永久免费视频| 久久久久久久久免费视频了| av欧美777| 亚洲免费av在线视频| 妹子高潮喷水视频| 国产黄色小视频在线观看| 天堂av国产一区二区熟女人妻 | 色播亚洲综合网| 成人永久免费在线观看视频| 淫秽高清视频在线观看| 国产aⅴ精品一区二区三区波| 国产精品久久久av美女十八| 99精品久久久久人妻精品| 色哟哟哟哟哟哟| 中国美女看黄片| 免费无遮挡裸体视频| 俄罗斯特黄特色一大片| 亚洲av电影不卡..在线观看| av中文乱码字幕在线| 亚洲色图 男人天堂 中文字幕| 曰老女人黄片| 成在线人永久免费视频| www.熟女人妻精品国产| 午夜福利18| 三级毛片av免费| 欧美日韩黄片免| 色精品久久人妻99蜜桃| 熟妇人妻久久中文字幕3abv| 国产一区二区三区在线臀色熟女| 欧美一级a爱片免费观看看 | 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美国产一区二区入口| 啪啪无遮挡十八禁网站| 又黄又爽又免费观看的视频| 最近视频中文字幕2019在线8| 真人做人爱边吃奶动态| 中文字幕人妻丝袜一区二区| 法律面前人人平等表现在哪些方面| 免费看a级黄色片| 亚洲欧美日韩高清专用| 国产精品1区2区在线观看.| 两个人视频免费观看高清| 日韩国内少妇激情av| 51午夜福利影视在线观看| www.熟女人妻精品国产| 欧美日本亚洲视频在线播放| 三级毛片av免费| 又紧又爽又黄一区二区| 亚洲av第一区精品v没综合| 搡老熟女国产l中国老女人| 黄片小视频在线播放| 久久草成人影院| av中文乱码字幕在线| 亚洲成人中文字幕在线播放| 国产精华一区二区三区| 精品少妇一区二区三区视频日本电影| 久久人妻福利社区极品人妻图片| 久久 成人 亚洲| a级毛片a级免费在线| 国产1区2区3区精品| 亚洲精品国产精品久久久不卡| 久久久久久久久免费视频了| 国产蜜桃级精品一区二区三区| 午夜福利成人在线免费观看| 无遮挡黄片免费观看| 国产69精品久久久久777片 | 亚洲美女黄片视频| 亚洲第一欧美日韩一区二区三区| 久久国产乱子伦精品免费另类| 免费在线观看亚洲国产| 97碰自拍视频| 美女 人体艺术 gogo| 亚洲国产高清在线一区二区三| 色精品久久人妻99蜜桃| 国产精品乱码一区二三区的特点| 一个人免费在线观看的高清视频| 久久久久久久久免费视频了| 久久99热这里只有精品18| 成人三级做爰电影| 一级a爱片免费观看的视频| 久久久久久亚洲精品国产蜜桃av| 欧美日本亚洲视频在线播放| 国产69精品久久久久777片 | 久久精品夜夜夜夜夜久久蜜豆 | 亚洲va日本ⅴa欧美va伊人久久| 欧美性猛交╳xxx乱大交人| 亚洲成av人片在线播放无| 淫秽高清视频在线观看| 国产激情偷乱视频一区二区| 国产一区二区在线av高清观看| 亚洲在线自拍视频| 日韩欧美三级三区| 亚洲人成网站高清观看| 精品乱码久久久久久99久播| 亚洲欧洲精品一区二区精品久久久| 国产精品av视频在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 欧美av亚洲av综合av国产av| 久久精品国产99精品国产亚洲性色| 12—13女人毛片做爰片一| 午夜精品在线福利| 亚洲中文字幕日韩| 18禁裸乳无遮挡免费网站照片| 又紧又爽又黄一区二区| 国模一区二区三区四区视频 | 久久香蕉精品热| tocl精华| 一二三四在线观看免费中文在| 最近最新免费中文字幕在线| 黄色成人免费大全| 亚洲国产欧美一区二区综合| 久久国产精品影院| 高清毛片免费观看视频网站| 久久久久久免费高清国产稀缺| 久久精品成人免费网站| 精品少妇一区二区三区视频日本电影| 1024视频免费在线观看| 小说图片视频综合网站| 日韩高清综合在线| 亚洲人成网站在线播放欧美日韩| 欧美成人一区二区免费高清观看 | 老司机在亚洲福利影院| 欧美日韩亚洲综合一区二区三区_| 国产成+人综合+亚洲专区| 欧美性猛交╳xxx乱大交人| 天堂av国产一区二区熟女人妻 | 很黄的视频免费| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 丰满人妻熟妇乱又伦精品不卡| 男女下面进入的视频免费午夜| 亚洲avbb在线观看| 国产主播在线观看一区二区| 丰满人妻熟妇乱又伦精品不卡| 精品国产美女av久久久久小说| 精品第一国产精品| 国产av不卡久久| 国产av麻豆久久久久久久| 国产伦人伦偷精品视频| 老汉色∧v一级毛片| 亚洲欧美精品综合久久99| 久久精品影院6| 色老头精品视频在线观看| 在线观看日韩欧美| 久久精品aⅴ一区二区三区四区| 天堂av国产一区二区熟女人妻 | 又黄又粗又硬又大视频| 毛片女人毛片| 欧美av亚洲av综合av国产av| 中文字幕久久专区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲男人的天堂狠狠| 亚洲国产欧美人成| 国产精品免费视频内射| 国产日本99.免费观看| 丰满人妻一区二区三区视频av | 色综合欧美亚洲国产小说| 五月玫瑰六月丁香| 欧美日韩瑟瑟在线播放| 黄色片一级片一级黄色片| 他把我摸到了高潮在线观看| 好看av亚洲va欧美ⅴa在| 两性夫妻黄色片| 99re在线观看精品视频| 99国产精品99久久久久| 日韩大码丰满熟妇| 国产在线观看jvid| 日日摸夜夜添夜夜添小说| 久久天躁狠狠躁夜夜2o2o| ponron亚洲| 精品久久久久久久久久久久久| 叶爱在线成人免费视频播放| 亚洲精品中文字幕在线视频| 精品久久久久久成人av| 久久久精品欧美日韩精品| 亚洲 欧美 日韩 在线 免费| 国产单亲对白刺激| 久久久水蜜桃国产精品网| 国产精品一区二区精品视频观看| 老司机午夜十八禁免费视频| 亚洲精品色激情综合| 亚洲一区中文字幕在线| 精品一区二区三区四区五区乱码| 亚洲,欧美精品.| 岛国在线观看网站| 中文字幕精品亚洲无线码一区| 国产精品免费一区二区三区在线| 免费在线观看视频国产中文字幕亚洲| 好男人在线观看高清免费视频| 日本 av在线| 搡老熟女国产l中国老女人| 精品欧美国产一区二区三| 久久精品aⅴ一区二区三区四区| 亚洲国产日韩欧美精品在线观看 | 亚洲五月天丁香| 午夜影院日韩av| 啦啦啦韩国在线观看视频| 亚洲五月天丁香| 无遮挡黄片免费观看| 久久 成人 亚洲| 午夜免费激情av| 色噜噜av男人的天堂激情| netflix在线观看网站| 1024手机看黄色片| 少妇被粗大的猛进出69影院| 欧美绝顶高潮抽搐喷水| 99久久国产精品久久久| 国产成年人精品一区二区| 国产乱人伦免费视频| 首页视频小说图片口味搜索|