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

    幾何形狀對十字傘充氣和滑翔性能的影響

    2022-11-15 01:57:32黃云堯武士輕張揚(yáng)
    航天返回與遙感 2022年5期
    關(guān)鍵詞:傘衣風(fēng)洞試驗(yàn)滑翔

    黃云堯 武士輕 張揚(yáng),*

    幾何形狀對十字傘充氣和滑翔性能的影響

    黃云堯1,3武士輕2,3張揚(yáng)1,3,*

    (1 南京航空航天大學(xué)航天學(xué)院,南京 210016)(2北京空間機(jī)電研究所,北京 100094)(3中國航天科技集團(tuán)有限公司航天進(jìn)入、減速與著陸技術(shù)實(shí)驗(yàn)室,北京 100094)

    為了研究幾何形狀對十字傘充氣和滑翔性能的影響,文章采用LS-DYNA軟件中的結(jié)構(gòu)化任意拉格朗日-歐拉(S-ALE)流固耦合求解器,利用數(shù)值仿真對5種不同幾何形狀的十字傘進(jìn)行了動(dòng)態(tài)開傘過程模擬,得到了十字傘展開過程中不同時(shí)刻的傘衣形狀,并對不同傘型的充氣性能進(jìn)行了對比分析。此外,對充氣穩(wěn)定后的十字傘,通過在其傘繩端施加恒定載荷使其具備滑翔能力,并對比分析了5種傘型的滑翔性能。數(shù)值結(jié)果表明,異形結(jié)構(gòu)對十字傘的充氣性能影響較大,但對滑翔性能的影響因傘型而異,十字傘的長寬比對充氣和滑翔性均具有較大影響。

    充氣性能 滑翔性能 十字傘 流固耦合 航天返回

    0 引言

    降落傘因具有質(zhì)量小、減速效果好等優(yōu)點(diǎn),已被廣泛應(yīng)用于載人航天、登月返回、空降兵訓(xùn)練以及定點(diǎn)空投等重大航空航天任務(wù)。傳統(tǒng)的空投系統(tǒng)容易受到風(fēng)場等環(huán)境因素的影響,一經(jīng)投放便無法對降落傘進(jìn)行后續(xù)控制。在惡劣的天氣條件下,空投行動(dòng)要么無法實(shí)施,要么嚴(yán)重偏離預(yù)定著陸區(qū),不僅降低了空投效果而且容易對地面人員及設(shè)施造成損壞[1]。

    沖壓式翼傘由于具備優(yōu)良的滑翔性能和可操縱性,不僅可以通過自主歸航控制實(shí)現(xiàn)載荷的精確定點(diǎn)回收,還可通過“雀降”操縱實(shí)現(xiàn)高效的安全無損著陸[2]。以翼傘為代表的精確空投系統(tǒng)(JPADS)[1]落點(diǎn)非常準(zhǔn)確,但成本較高,原因是為翼傘提供高機(jī)動(dòng)性的高性能減速器的制造需高昂的人工成本,同時(shí)為了能夠承受傘衣打開過程中的巨大沖擊力并提供足夠的滑翔和轉(zhuǎn)彎特性,所使用的材料也十分昂貴,加工難度大,這些因素導(dǎo)致基于翼傘的精確空投系統(tǒng)的研制成本相對過高[1]。傳統(tǒng)的不受控系統(tǒng)雖然研制成本相對較低,但落點(diǎn)誤差較大。因此,為了兼顧空投的準(zhǔn)確性和經(jīng)濟(jì)性,一種利用現(xiàn)有技術(shù)來提供較為精確落點(diǎn)位置的低成本可控空投系統(tǒng)(Affordable Guided Airdrop System,AGAS)應(yīng)運(yùn)而生[1]。

    最初的AGAS研究基于圓形傘,當(dāng)前則主要針對十字傘[1]。與翼傘相比,十字傘的構(gòu)造既簡單又經(jīng)濟(jì),且氣動(dòng)特性和機(jī)理的研究較為成熟;與標(biāo)準(zhǔn)的圓形傘、環(huán)帆傘以及其他類型的降落傘相比,十字傘的高幾何孔隙率使其在自由來流和尾流區(qū)中能夠更好地維持穩(wěn)定性。高幾何孔隙率雖然可能導(dǎo)致較大的下降速率,但由于減少了總的下降時(shí)間,可以在下降后迅速清理空域,并快速穿越未知或變化劇烈的風(fēng)況[2]。雖然十字傘的控制能力與翼傘相比仍較為有限,但最近一些基于十字傘的研究發(fā)展出了各種各樣的控制策略和手段,在應(yīng)用于AGAS方面展現(xiàn)出了巨大的潛力[2-7]。

    目前針對AGAS的研究主要有風(fēng)洞試驗(yàn)、飛行試驗(yàn)以及數(shù)值仿真等手段。文獻(xiàn)[4-5]利用NASA的6.096m垂直風(fēng)洞對可控十字傘進(jìn)行試驗(yàn),發(fā)展出了一種可對十字傘航向穩(wěn)定性進(jìn)行控制和調(diào)整的控制器,并對5種不同幾何形狀傘衣的十字傘進(jìn)行了風(fēng)洞試驗(yàn),比較了不同傘型的滑翔特性和落點(diǎn)精度,確定出能提供最佳落點(diǎn)精度的十字傘型;文獻(xiàn)[6]通過飛行測試研究了標(biāo)準(zhǔn)型十字傘和混合型十字傘的滑翔性能,通過對十字傘的形狀以及傘繩長度進(jìn)行調(diào)整,達(dá)到控制十字傘的滑翔比和轉(zhuǎn)彎性能的目的;文 獻(xiàn)[7]利用LS-DYNA軟件模擬了操縱十字傘的開傘過程,并通過對傘繩的長度進(jìn)行調(diào)整實(shí)現(xiàn)控制的效果;文獻(xiàn)[8]通過ALE(Arbitrary Lagrange-Euler)算法對十字型阻力傘的充氣過程進(jìn)行數(shù)值模擬,討論了不同開傘速度、計(jì)算域網(wǎng)格密度等對阻力傘開傘動(dòng)載荷的影響;文獻(xiàn)[9]采用S-ALE(Structured ALE)算法對半折疊狀態(tài)的十字型阻力傘的開傘過程進(jìn)行仿真計(jì)算,分析了傘衣充氣張滿過程,獲得了十字傘的阻力特性。

    相對于風(fēng)洞測試和飛行試驗(yàn),數(shù)值仿真受不可控因素(如風(fēng)速、氣候)的影響較小,更改參數(shù)更為靈活,研究成本較低。因此,通過數(shù)值仿真研究十字傘的工作過程已逐漸成為一個(gè)熱點(diǎn),但目前大多是基于標(biāo)準(zhǔn)十字傘外形進(jìn)行充氣過程模擬[7-9],針對不同形狀十字傘的充氣和滑翔特性的相關(guān)數(shù)值仿真研究仍較為匱乏。本文將參照文獻(xiàn)[5]中的垂直風(fēng)洞試驗(yàn)研究,針對5種不同幾何形狀的十字傘,采用S-ALE方法仿真分析十字傘傘衣幾何形狀對其充氣和滑翔性能的影響。

    1 數(shù)值仿真方法

    1.1 方法介紹

    ALE方法最早出現(xiàn)于數(shù)值模擬流體動(dòng)力學(xué)問題的有限差分方法中[10],這種方法兼具Lagrange方法和Euler方法的優(yōu)勢:首先,在物體邊界運(yùn)動(dòng)的處理上引進(jìn)了Lagrange方法的特點(diǎn),能夠使內(nèi)部網(wǎng)格單元有效地跟蹤物體邊界的運(yùn)動(dòng);其次,在內(nèi)部網(wǎng)格的劃分上,吸收了Euler方法的長處,使內(nèi)部網(wǎng)格與物質(zhì)材料相互獨(dú)立存在,網(wǎng)格可以根據(jù)定義的參數(shù)在求解過程中適當(dāng)調(diào)整位置,使得網(wǎng)格不至于出現(xiàn)嚴(yán)重的畸變。

    S-ALE方法與ALE方法的理論完全一致,二者的區(qū)別在于S-ALE方法進(jìn)行流場建模的方式與ALE方法直接進(jìn)行網(wǎng)格劃分不同,S-ALE流場域只需要設(shè)置三個(gè)維度的位置和網(wǎng)格數(shù)量,然后通過關(guān)鍵字對網(wǎng)格進(jìn)行物質(zhì)填充即可。S-ALE求解器是在ALE求解器的基礎(chǔ)上發(fā)展而來,適用于結(jié)構(gòu)化網(wǎng)格流固耦合問題的求解,在理論和算法相同的情況下,利用結(jié)構(gòu)網(wǎng)格節(jié)點(diǎn)排序的規(guī)律性對ALE算法進(jìn)行簡化,通過關(guān)鍵字控制自動(dòng)生成結(jié)構(gòu)化網(wǎng)格,大大減少了計(jì)算時(shí)間和內(nèi)存開銷。

    1.2 控制方程

    在ALE方法的描述中,引入拉格朗日和歐拉坐標(biāo)外的第三個(gè)任意參照坐標(biāo),則拉格朗日坐標(biāo)下物體的運(yùn)動(dòng)可以描述為[11]:

    因此ALE算法的控制方程可以由質(zhì)量、動(dòng)量以及能量守恒方程給定,具體為:

    十字傘結(jié)構(gòu)的控制方程[12]為:

    2 模型設(shè)置

    2.1 背景介紹

    文獻(xiàn)[5]的垂直風(fēng)洞試驗(yàn)在不嚴(yán)重影響十字傘穩(wěn)定性的基礎(chǔ)上,為了提高十字傘的滑翔能力,研究了5種傘衣面積相同而外形不同的十字傘,如圖1所示。風(fēng)洞試驗(yàn)中將A型傘設(shè)置為長寬比為3.0的基礎(chǔ)傘型;B、C兩種傘型針對十字傘的長寬比進(jìn)行了調(diào)整,B型傘有著較大的長寬比,C型傘的長寬比較??;D、E型傘為異形結(jié)構(gòu),在長寬比為3.0的基礎(chǔ)上提出對十字傘的形狀進(jìn)行修改,以期望獲得較好的滑翔能力;其中D型傘采用的是文獻(xiàn)[6]中首次提及到的雜交形狀,E型傘采用非對稱設(shè)置。對5種不同形狀的十字傘在垂直風(fēng)洞中相同的環(huán)境下進(jìn)行了測試,試驗(yàn)比較分析了5種傘型在風(fēng)洞試驗(yàn)中的穩(wěn)定性,并通過控制器對十字傘傘繩進(jìn)行操控使其具有一定的滑翔能力,然后通過傳感器測得的數(shù)據(jù)來計(jì)算比較各個(gè)傘型的滑翔比。試驗(yàn)結(jié)果表明:與基礎(chǔ)的傘型A相比,B型傘具有較好的穩(wěn)定性,但滑翔能力較差;C型傘滑翔能力較強(qiáng)但穩(wěn)定性較差;D型傘滑翔能力較強(qiáng),但穩(wěn)定性較差難以控制;E型傘的滑翔能力和穩(wěn)定性能都較差。

    圖1 不同傘型十字傘形狀示意

    2.2 計(jì)算模型和參數(shù)設(shè)置

    本文將參照文獻(xiàn)[5]中提及到的5種十字傘傘衣模型,通過仿真研究十字傘的充氣性能和滑翔特性,其中充氣性能主要通過充氣后十字傘是否穩(wěn)定和十字傘俯仰特性表現(xiàn),而滑翔特性主要通過十字傘的滑翔比G來體現(xiàn),滑翔比越高則滑翔性能越優(yōu)秀。5種傘型的傘衣相關(guān)參數(shù)如表1所示。

    表1 5種傘型的傘衣尺寸參數(shù)

    Tab.1 Dimensions of the five canopy types

    為了簡化計(jì)算,對十字傘模型作如下假設(shè):1)傘衣為如圖2所示的簡單交叉結(jié)構(gòu),傘衣的每條傘輻上各有4條傘繩,4條傘繩通過連接繩與重物相連;2)模型只考慮傘繩的連接,對于傘衣徑向帶、頂孔繩等不作具體處理;3)傘衣采用殼單元,傘繩采用梁單元;4)重物為剛體,下表面通過連接繩與固定端相連,傘衣平面到重物的垂直距離為4m,重物通過一根連接繩與固定端相連,重物下表面與固定端的距離為1.6m。A型傘系統(tǒng)的幾何模型,如圖2所示。試驗(yàn)將通過5種傘型在展開過程中的表現(xiàn)及俯仰角來判斷各種傘型的穩(wěn)定性,試驗(yàn)中俯仰角即為十字傘中心點(diǎn)和重物中心點(diǎn)連線與豎直方向的夾角[3]。

    在實(shí)際風(fēng)洞試驗(yàn)中,重物端有兩個(gè)執(zhí)行器用于動(dòng)態(tài)控制傘繩的長短。當(dāng)試驗(yàn)中兩個(gè)執(zhí)行器控制的靜態(tài)傘繩長度和動(dòng)態(tài)傘繩長度均小于標(biāo)稱長度時(shí),十字傘開始滑翔[4]。

    為了比較各種十字傘型的滑翔特性,本文在充氣穩(wěn)定后的十字傘傘繩端施加一段時(shí)間的恒定載荷,用于模擬實(shí)際風(fēng)洞實(shí)驗(yàn)中執(zhí)行器對十字傘的控制作用,使十字傘具備一定的滑翔能力,具體操作為:通過關(guān)鍵字*DEFINE_CURVE來定義時(shí)間載荷曲線,按照定義的曲線,通過*LOAD_NODE將大小恒定的1kN載荷施加在其中的兩根傘繩上,如圖2中紅色箭頭所示。

    傘衣材料密度為74.00kg/m3,彈性模量為0.064GPa,泊松比為0.14,傘衣厚度為0.002 5m。傘繩材料密度為623.00kg/m3,彈性模量為10.90GPa,傘繩截面為圓形,圓截面的面積為7.068 6×10–6m2。

    圖3 側(cè)向力和軸向力示意

    3 關(guān)鍵字設(shè)置

    LS-DYNA軟件中S-ALE求解器的流場解算設(shè)置步驟為:1)通過ALE_STRUCTURED_CONTROL_POINTS關(guān)鍵字控制流場在三維方向上的位置和網(wǎng)格數(shù)量;2)通過ALE_ STRUCTURED_MESH創(chuàng)建網(wǎng)格;3)最后通過INITIAL_ VOLUME_FRACTION_GEOMETRY關(guān)鍵字對流場區(qū)域進(jìn)行物質(zhì)填充,對傘衣附近的流體域網(wǎng)格進(jìn)行加密,加密區(qū)域網(wǎng)格尺度與傘衣表面網(wǎng)格尺度相同,加密區(qū)域以外至遠(yuǎn)場區(qū)域的網(wǎng)格逐漸稀疏,流場計(jì)算域網(wǎng)格劃分如圖4所示,流場入口處施加20m/s的恒定來流速度。

    圖4 流場計(jì)算域網(wǎng)格劃分

    除上述S-ALE特有的關(guān)鍵字外,其他的關(guān)鍵字設(shè)置[14-16]與ALE求解器相同,具體設(shè)置如表2所示。

    表2 流固耦合過程關(guān)鍵字設(shè)置

    Tab.2 Major keywords for fluid-structure interaction (FSI) setup

    4 仿真結(jié)果及討論

    4.1 傘的充氣性能分析

    十字傘的充氣性能主要通過充氣穩(wěn)定后的傘衣形狀和俯仰特性來表現(xiàn),這里選取傘衣變化幅度較大的幾個(gè)時(shí)間節(jié)點(diǎn)(0.7s、1s、1.5s)記錄圖像,進(jìn)行比較,并將整個(gè)試驗(yàn)過程劃分為預(yù)充氣階段(0~0.7s)、充分充氣階段(0.7~1.5s)、充滿后工作階段(1.5~2.0s)第3個(gè)階段來分析十字傘的充氣性能,同時(shí)記錄仿真結(jié)束時(shí)(2s時(shí))十字傘傘衣和重物的相對位置,便于直觀地比較十字傘的偏轉(zhuǎn)情況,各傘型傘衣的具體狀態(tài)如圖6所示,對應(yīng)的流固耦合現(xiàn)象如圖7所示。俯仰特性主要通過最大俯仰角表現(xiàn),通過仿真得到傘型A~E對應(yīng)的最大俯仰角數(shù)據(jù)分別為0.86°、0.96°、10.7°、10.7°、16.0°。

    結(jié)合圖5~6可以清楚地看到十字傘充分展開的完整過程,這一過程大致可分為3個(gè)階段:預(yù)充氣狀態(tài)、充分充氣狀態(tài)、充滿氣后工作階段[17]。

    1)預(yù)充氣階段(0~0.7s)。自流場入口處氣流開始垂直上升,在充氣后0.4s左右開始與傘衣接觸,隨著氣流在傘衣部分的聚集,導(dǎo)致傘衣頂部逐漸形成高壓區(qū),并從傘頂逐漸擴(kuò)散到傘衣底部,傘頂變化比較明顯;當(dāng)空氣向十字傘頂部聚集時(shí),傘衣和傘繩連接處因有束縛較少出現(xiàn)“卷邊”現(xiàn)象;充氣達(dá)0.7s時(shí),十字傘初具雛形,此時(shí)5種傘型都呈較為穩(wěn)定狀態(tài)。

    圖5 不同傘型傘衣變化

    2)充分充氣階段(0.7~1.5s)。從開始充氣后的0.7s到1.0s,傘衣的變化幅度最為劇烈,由于氣流對十字傘的持續(xù)作用,傘衣繼續(xù)膨脹,到達(dá)1.0s左右時(shí),傘衣變形達(dá)到最大,隨后,聚集在傘頂端的氣流向四周流動(dòng),內(nèi)部氣流速度逐漸變小并繞傘衣底邊處流出,傘臂逐漸向外張開,且降落傘開始旋轉(zhuǎn)和擺動(dòng);到達(dá)1.5s時(shí),傘衣形狀經(jīng)過多次調(diào)整后趨于穩(wěn)定,基本不再發(fā)生大幅度的變化。在這一階段可以明顯觀察到五種傘型的變化趨勢略有不同,根據(jù)圖6中5種傘型的速度分布云圖可以看出:A、B型傘的穩(wěn)定性較強(qiáng),C型傘在氣流通過后傘衣變化較為劇烈,D、E型傘由于非對稱性的影響,氣流通過時(shí),傘衣表現(xiàn)出了更為嚴(yán)重的失穩(wěn)現(xiàn)象。

    圖6 不同形狀十字傘的流固耦合現(xiàn)象

    3)充滿后工作階段(1.5~2.0s)。這個(gè)階段的十字傘傘衣形狀變化幅度較小,表現(xiàn)出有規(guī)律的“呼吸”現(xiàn)象。原先在傘頂上方的對稱渦旋變得極為不穩(wěn)定,在不對稱的渦旋和壓差的作用下,傘衣開始向其他方向擺動(dòng),擺動(dòng)的位置也會(huì)在相應(yīng)部位產(chǎn)生非對稱渦旋,但大小和發(fā)展趨勢不確定,傘因此產(chǎn)生繞固定端打轉(zhuǎn)、前后左右隨機(jī)擺動(dòng)等現(xiàn)象[16],十字傘的偏轉(zhuǎn)主要發(fā)生在這個(gè)階段。因此,此階段5種傘型的區(qū)別最為明顯。觀察圖5~6中充氣2s時(shí)各傘型的位置和速度分布云圖,可以明顯看出,A、B型傘偏轉(zhuǎn)較小,C、D、E型傘偏轉(zhuǎn)較大。

    綜合上述分析以及最大俯仰角數(shù)據(jù),可以得出:A、B型傘最大俯仰角較小,有較強(qiáng)的穩(wěn)定性,充氣性能較好;C、D、E型傘的最大俯仰角較大,穩(wěn)定性較差,充氣性能較差。數(shù)值模擬結(jié)果與文獻(xiàn)[5]中垂直風(fēng)洞試驗(yàn)得到的結(jié)論一致。

    4.2 傘的滑翔特性分析

    通過LS-DYNA對有無施加載荷得到的受力情況進(jìn)行比較分析,5種傘型施加載荷前后的軸向力與側(cè)向力如圖7所示。記錄仿真結(jié)束時(shí)各個(gè)傘型的形狀如圖8所示。

    圖7 不同形狀十字傘受力曲線(1.2s施加載荷)

    圖8 2s時(shí)不同傘型受力后形狀變化(從左到右依次為A、B、C、D、E型)

    結(jié)合圖7~8進(jìn)行分析發(fā)現(xiàn),A型傘在施加控制之后,側(cè)向力提升較為明顯,且相較未施加載荷(圖5)有明顯的滑翔能力;B型傘和E型傘即使在施加控制之后側(cè)向力也無明顯提升,滑翔能力較差。而C、D型傘在施加控制之前就已經(jīng)有較大的側(cè)向力,滑翔能力相較于其他傘型較強(qiáng)。

    從圖7曲線中提取3組數(shù)據(jù)計(jì)算施加控制后5種十字傘的滑翔比,結(jié)果如表3所示:

    表3 5種十字傘不同時(shí)刻的滑翔比

    Tab.3 Glide ratios of cruciform parachutes at different time steps

    綜上,滑翔能力的數(shù)值模擬結(jié)果表明:與基準(zhǔn)的A型傘相比,C、D型傘的滑翔能力有顯著提高,B、E型傘的滑翔能力較差,即使在施加載荷后也無明顯提升,該項(xiàng)模擬結(jié)果與文獻(xiàn)[5]中的垂直風(fēng)洞試驗(yàn)結(jié)論一致。

    5 結(jié)束語

    本文利用LS-DYNA中的S-ALE求解器對5種不同傘衣形狀的十字傘進(jìn)行了數(shù)值模擬,分析比較幾何外形對十字傘充氣和滑翔性能的影響,數(shù)值模擬結(jié)果顯示,不同幾何形狀的十字傘充氣展開后的充氣性能和滑翔特性有著較大的差異,具體表現(xiàn)為:

    1)十字傘的長寬比A會(huì)影響到俯仰特性和穩(wěn)定性,對滑翔能力的影響較大:與A型基準(zhǔn)傘相比,A較大的傘(B型傘)展開后的俯仰特性和穩(wěn)定性較好,而A較小的傘(C型傘)俯仰特性和穩(wěn)定性較差。A較小的傘的滑翔能力較強(qiáng),而A較大的傘滑翔能力表現(xiàn)較差。

    2)異形結(jié)構(gòu)會(huì)嚴(yán)重影響十字傘的俯仰特性和穩(wěn)定性,但對滑翔能力的影響因傘型不同而有所區(qū)別,其中D型傘相較于A型基準(zhǔn)傘而言,俯仰特性和穩(wěn)定性較差,但滑翔能力有一定的提升;而非對稱結(jié)構(gòu)的E型傘在數(shù)值模擬中既沒有很好的穩(wěn)定性和俯仰特性,也無法在施加載荷控制后使滑翔能力得到顯著提升。

    3)通過在十字傘的傘繩端施加恒定載荷的方式來模擬垂直風(fēng)洞試驗(yàn)中對傘繩的控制效果,由于風(fēng)洞試驗(yàn)中未明確施加載荷的大小,且載荷大小可能隨著時(shí)間而變化,因此無法定量比較風(fēng)洞試驗(yàn)和數(shù)值模擬得到的滑翔比結(jié)果,后續(xù)將分析降落傘在重物上的吊點(diǎn)分布、重物的幾何尺寸、質(zhì)量、慣量,固定端與重物的連接等因素對結(jié)果產(chǎn)生的影響,更加準(zhǔn)確地模擬風(fēng)洞試驗(yàn)中的控制效果,得到更為精確的結(jié)果。

    [1] CIVELEK B, KIVRAK S. A Review on the Precision Guided Airdrop Systems[J]. International Journal of Latest Technology in Engineering, Management & Applied Science, 2019, 8(1): 15-17.

    [2] FAGLEV C P, SEIDEL J, MCLAUGHLIN T E, et al. Computational Study of Air Drop Control Mechanisms for Cruciform Parachutes[C]//24th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar 2017, June 5-9, 2017, Denver, Colorado, USA. Reston, VA: AIAA, 2017: 5-9.

    [3] NOETSCHER G, ROSE T M, BERGERON K. Modeling and Experimental Efforts towards Robust Low-Cost Cruciform Canopy Control[C]//Aerodynamic Decelerator Systems Technology: AIAA Aviation Forum 2019, June 17-21, 2019, Dallas, Texas, USA. Reston, VA: AIAA, 2019.

    [4] HERRINGTON S M, RENZELMAN J T, FIELDS T D, et al.Vertical Wind-Tunnel Testing of Steerable Cruciform Parachute System[J].Journal of Aircraft, 2018, 56(2):747-757.

    [5] HERRINGTON S M, RENZELMAN J T, FIELDS T D, et al. Experimental Investigation into the Effects of Geometry on the Glide Performance of Cruciform Parachutes[C]//Aerodynamic Decelerator Systems Technology: AIAA Aviation Forum 2019, June 17-21, 2019, Dallas, Texas, USA. Reston, VA: AIAA, 2019.

    [6] JEAN P, JASON P, ERIC B, et al. Glide Performance Study of Standard and Hybrid Cruciform Parachutes[C]//17th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, May 19-22, 2003, Monterey, California. [S.l.]: AIAA,2003.

    [7] GAO Xinglong, ZHANG Qingbin, CHEN Qin, et al. Fluid-structure Interactions on Steerable Cruciform Parachute Inflation Dynamics[C]//5th International Conference on Mechanical and Aeronautical Engineering: ICMAE 2019, Sanya, China, December 12-15, 2019. Bristol, UK: Institute of Physics, 2020.

    [8] 連紅博, 喻東明, 周昊, 等. 十字型阻力傘充氣過程的數(shù)值模擬[J]. 航空兵器, 2021, 28(6): 83-87.

    LIAN Hongbo, YU Dongming, ZHOU Hao, et al. Numerical Simulation on Inflation Process of Cruciform Drag Parachute[J]. Aero Weaponry, 2021, 28(6): 83-87. (in Chinese)

    [9] 楊品, 胡雪鵬, 金樓. 基于LS-DYNA的某型阻力傘半折疊數(shù)值仿真研究[C]//第八屆中國航空學(xué)會(huì)青年科技論壇, 2018年11月4-6日, 江門. 2018: 774-781.

    YANG Pin, HU Xuepeng, JIN Lou. Numerical Simulation Study of a Certain Drag Parachute Based on LS-DYNA[C]//The 8th Youth Forum of Chinese Society of Aeronautics, November 4-6, 2018, Jiangmen. 2018: 774-781. (in Chinese)

    [10] 房明. 火箭橇-阻力傘系統(tǒng)力學(xué)特性與流場分析[D].南京: 南京航空航天大學(xué), 2018.

    FANG Ming. Dynamic Characteristics and Flow Field Analyses on a Rocket Sled-Drag Parachute System[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2018. (in Chinese)

    [11] 李裕春, 時(shí)黨勇, 趙遠(yuǎn). ANSYS11.0/LS-DYNA基礎(chǔ)理論與工程實(shí)踐[M]. 北京: 中國水利水電出版社, 2008.

    LI Yuchun, SHI Dangyong, ZHAO Yuan. ANSYS11.0/LS-DYNA Foundation Theory and Engineering Practice[M]. Beijing: China Water & Power Press, 2008. (in Chinese)

    [12] 徐琳, 梁建, 段麗華, 等. 基于CSD/CFD舵面氣動(dòng)力流固耦合仿真分析[J]. 航空兵器, 2020, 27(2): 47-52.

    XU Lin, LIANG Jian, DUAN Lihua, et al. Fluid Solid Coupling Simulation Analysis of Rudder Surface Aerodynamic Based on CSD/CFD[J]. Aero Weaponry, 2020, 27(2): 47-52. (in Chinese)

    [13] JORGENSEN D S. Cruciform Parachute Aerodynamics[D]. Leicester: University of Leicester, 1982.

    [14] 王中陽. 降落傘充氣過程動(dòng)力學(xué)數(shù)值模擬[D]. 南京: 南京航空航天大學(xué), 2013.

    WANG Zhongyang. The Dynamic Numerical Simulation of the Parachute Inflation[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2013. (in Chinese)

    [15] 賈賀, 榮偉, 陳國良. 基于LS-DYNA的降落傘傘衣織物透氣性參數(shù)仿真驗(yàn)證[J]. 航天返回與遙感, 2009, 30(1): 15-20.

    JIA He, RONG Wei, CHEN Guoliang. The Use of LS-DYNA to Simulate the Permeability Parameters of the Parachute Canopy[J]. Spacecraft Recovery & Remote Sensing, 2009, 30(1): 15-20. (in Chinese)

    [16] 賈賀, 榮偉, 陳國良. 基于LS-DYNA軟件的降落傘充氣過程仿真研究[J]. 航天器環(huán)境工程, 2010, 27(3): 367-373.

    JIA He, RONG Wei, CHEN Guoliang. The Simulation of Parachute Inflation Process Based on LS-DYNA Software[J]. Spacecraft Environment Engineering, 2010, 27(3): 367-373. (in Chinese)

    [17] 程涵, 余莉, 李勝全. 基于ALE的降落傘充氣過程數(shù)值仿真[J]. 南京航空航天大學(xué)學(xué)報(bào), 2012, 44(3): 290-293.

    CHENG Han, YU Li, LI Shengquan. Numerical Simulation of Parachute Inflation Process Based on ALE[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2012, 44(3): 290-293. (in Chinese)

    Effects of Geometry on the Inflation and Gliding Performance of Cruciform Parachutes

    HUANG Yunyao1,3WU Shiqing2,3ZHANG Yang1,3,*

    (1 College of Astronautics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)(2 Beijing Institute of Space Mechanics & Electricity, Beijing 100094, China)(3 Laboratory of Aerospace Entry, Descent and Landing Technology, CASC, Beijing 100094, China)

    In order to study the effect of geometry on the inflation and glide performance of cruciform parachute, the S-ALE (Structured Arbitrary Lagrange-Euler) fluid-structure interaction solver in LS-DYNA is employed. Cruciform parachutes with five different canopy shapes are taken and the corresponding dynamic opening process are simulated. The inflated shapes of the cruciform parachutes at different time instants are obtained and the inflation performances are compared. Furthermore, loads are applied to these stably inflated cruciform parachutes to test and compare the glide performance. Numerical results show that the abnormal structure has a great influence on the inflation performance but varies with the type of the parachute on the glide performance, whereas the aspect ratio impacts a lot on both of the inflation and glide performance.

    inflation performance; glide performance; cruciform parachutes; fluid-structure interaction; spacecraft recovery

    O35

    A

    1009-8518(2022)05-0048-11

    10.3969/j.issn.1009-8518.2022.05.005

    2022-07-05

    航天進(jìn)入減速與著陸技術(shù)實(shí)驗(yàn)室開發(fā)基金(EDL19092128)

    黃云堯, 武士輕, 張揚(yáng). 幾何形狀對十字傘充氣和滑翔性能的影響[J]. 航天返回與遙感, 2022, 43(5): 48-58.

    HUANG Yunyao, WU Shiqing, ZHANG Yang. Effects of Geometry on the Inflation and Gliding Performance of Cruciform Parachutes[J]. Spacecraft Recovery & Remote Sensing, 2022, 43(5): 48-58. (in Chinese)

    黃云堯,男,1997年生,南京航空航天大學(xué)機(jī)械專業(yè)在讀碩士研究生,主要研究方向?yàn)榱鞴恬詈?、航天器返回技術(shù)。E-mail:huangyy1415@qq.com。

    武士輕,男,1982年生,2008年獲長春理工大學(xué)機(jī)械電子工程專業(yè)碩士學(xué)位,高級(jí)工程師。主要研究方向?yàn)楹教炱鳠o損著陸技術(shù)。E-mail:wushi508@foxmail.con。

    張揚(yáng),男,1987年生,2016年獲南京航空航天大學(xué)流體力學(xué)專業(yè)博士學(xué)位,副研究員。主要研究方向?yàn)榱鞴恬詈?、航天器返回技術(shù)。E-mail:zhy@nuaa.edu.cn。

    (編輯:夏淑密)

    猜你喜歡
    傘衣風(fēng)洞試驗(yàn)滑翔
    蘭花螳螂會(huì)滑翔
    小型降落傘充氣和傘衣塌陷過程
    攻天掠地的先鋒武器——滑翔導(dǎo)彈
    設(shè)計(jì)參數(shù)及大氣參數(shù)對降落傘充氣性能的影響
    低風(fēng)壓架空導(dǎo)線的風(fēng)洞試驗(yàn)
    電線電纜(2017年5期)2017-10-18 00:52:03
    讓你心跳加速的“滑翔飛板”
    傘衣透氣性對翼傘氣動(dòng)特性的影響
    滾轉(zhuǎn)機(jī)動(dòng)載荷減緩風(fēng)洞試驗(yàn)
    空中滑翔大比拼(下)——滑翔傘
    遮擋條件下超高層建筑風(fēng)洞試驗(yàn)研究
    重慶建筑(2014年12期)2014-07-24 14:00:32
    午夜久久久在线观看| 老司机深夜福利视频在线观看| 国产亚洲av高清不卡| 久久久国产欧美日韩av| 免费女性裸体啪啪无遮挡网站| 99热国产这里只有精品6| 国产高清videossex| 一本大道久久a久久精品| 12—13女人毛片做爰片一| 亚洲精品一二三| 久久亚洲真实| 法律面前人人平等表现在哪些方面| av电影中文网址| 韩国精品一区二区三区| 19禁男女啪啪无遮挡网站| 中国美女看黄片| 国产亚洲午夜精品一区二区久久| 国产精品免费视频内射| 国产一卡二卡三卡精品| 欧美乱妇无乱码| 国产高清videossex| 国产精品98久久久久久宅男小说| 99re在线观看精品视频| 亚洲一码二码三码区别大吗| 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲| 免费少妇av软件| 久久中文字幕一级| 首页视频小说图片口味搜索| 国产精品熟女久久久久浪| 999久久久国产精品视频| 高清视频免费观看一区二区| 69精品国产乱码久久久| 男男h啪啪无遮挡| 每晚都被弄得嗷嗷叫到高潮| 视频区欧美日本亚洲| 如日韩欧美国产精品一区二区三区| 正在播放国产对白刺激| 欧美日韩黄片免| 亚洲精品久久成人aⅴ小说| 99国产综合亚洲精品| 欧美精品一区二区免费开放| 欧美黄色淫秽网站| 五月天丁香电影| av免费在线观看网站| 久久精品国产a三级三级三级| 精品人妻熟女毛片av久久网站| 一区二区三区精品91| 国产人伦9x9x在线观看| 老汉色∧v一级毛片| 亚洲男人天堂网一区| 十八禁人妻一区二区| 黑人欧美特级aaaaaa片| 久久久久久久大尺度免费视频| 91精品三级在线观看| 国产精品偷伦视频观看了| 欧美精品av麻豆av| 飞空精品影院首页| 成在线人永久免费视频| 国产精品秋霞免费鲁丝片| 亚洲国产欧美网| 亚洲全国av大片| 99精品欧美一区二区三区四区| 一区二区三区精品91| 亚洲欧洲精品一区二区精品久久久| 国产精品98久久久久久宅男小说| 美女高潮喷水抽搐中文字幕| 波多野结衣一区麻豆| 国产在线精品亚洲第一网站| 每晚都被弄得嗷嗷叫到高潮| 90打野战视频偷拍视频| 日本vs欧美在线观看视频| 咕卡用的链子| 亚洲第一av免费看| 纵有疾风起免费观看全集完整版| 成人三级做爰电影| 国产精品1区2区在线观看. | 免费黄频网站在线观看国产| av网站免费在线观看视频| 国产精品免费大片| 国产精品亚洲一级av第二区| 国产午夜精品久久久久久| 两性夫妻黄色片| 淫妇啪啪啪对白视频| 国产精品99久久99久久久不卡| av在线播放免费不卡| 满18在线观看网站| 中文字幕最新亚洲高清| 国产日韩欧美视频二区| 国产精品久久久人人做人人爽| 久久99一区二区三区| 国产成人欧美在线观看 | 丁香六月天网| 一夜夜www| 久久久久视频综合| 老熟妇乱子伦视频在线观看| 久久中文看片网| 嫁个100分男人电影在线观看| 成人18禁在线播放| 国产精品1区2区在线观看. | 老汉色av国产亚洲站长工具| 日本撒尿小便嘘嘘汇集6| 成人精品一区二区免费| 久久国产精品男人的天堂亚洲| 亚洲欧美一区二区三区黑人| 熟女少妇亚洲综合色aaa.| 欧美国产精品va在线观看不卡| 丁香六月天网| 黄频高清免费视频| 伦理电影免费视频| 大型黄色视频在线免费观看| 亚洲九九香蕉| 黄色丝袜av网址大全| 亚洲欧美激情在线| 一个人免费在线观看的高清视频| 国产成人啪精品午夜网站| 97人妻天天添夜夜摸| 欧美激情 高清一区二区三区| 欧美人与性动交α欧美精品济南到| 啦啦啦在线免费观看视频4| 嫩草影视91久久| 国产精品av久久久久免费| 老司机午夜福利在线观看视频 | 天天添夜夜摸| 法律面前人人平等表现在哪些方面| 黄网站色视频无遮挡免费观看| 99国产精品免费福利视频| 最黄视频免费看| 亚洲av片天天在线观看| 国产主播在线观看一区二区| 欧美人与性动交α欧美精品济南到| 亚洲av片天天在线观看| 另类亚洲欧美激情| 丝袜人妻中文字幕| 亚洲 欧美一区二区三区| 国产在线一区二区三区精| 少妇裸体淫交视频免费看高清 | 亚洲熟妇熟女久久| 最新在线观看一区二区三区| 天堂俺去俺来也www色官网| 亚洲精品美女久久av网站| 国产麻豆69| 欧美 亚洲 国产 日韩一| h视频一区二区三区| 嫁个100分男人电影在线观看| 国产免费视频播放在线视频| 精品视频人人做人人爽| 久9热在线精品视频| 91老司机精品| 国产一区有黄有色的免费视频| 1024视频免费在线观看| 美国免费a级毛片| av在线播放免费不卡| 国产精品偷伦视频观看了| 国产成人精品久久二区二区免费| 国产高清激情床上av| 国产黄频视频在线观看| 美国免费a级毛片| av视频免费观看在线观看| 嫁个100分男人电影在线观看| 成人精品一区二区免费| 国产不卡av网站在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 99国产精品一区二区三区| 亚洲精品一二三| 无人区码免费观看不卡 | 99riav亚洲国产免费| 国产99久久九九免费精品| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 99九九在线精品视频| 久热这里只有精品99| tocl精华| 在线 av 中文字幕| 色94色欧美一区二区| 久久久久国产一级毛片高清牌| 在线观看免费视频日本深夜| 亚洲男人天堂网一区| 老鸭窝网址在线观看| 国产一区二区 视频在线| 在线十欧美十亚洲十日本专区| 无遮挡黄片免费观看| 午夜福利在线观看吧| 国精品久久久久久国模美| 国产又爽黄色视频| 亚洲一卡2卡3卡4卡5卡精品中文| av免费在线观看网站| av有码第一页| 国产高清videossex| www.自偷自拍.com| 国产色视频综合| 国产成+人综合+亚洲专区| av视频免费观看在线观看| 人人妻人人澡人人看| 亚洲午夜精品一区,二区,三区| 国产在线视频一区二区| 国产成人欧美在线观看 | 国产精品二区激情视频| 国产区一区二久久| 欧美亚洲日本最大视频资源| 精品免费久久久久久久清纯 | 亚洲三区欧美一区| 18禁美女被吸乳视频| av又黄又爽大尺度在线免费看| 国产亚洲一区二区精品| 日本a在线网址| 国产精品一区二区在线不卡| 999久久久国产精品视频| 国产精品av久久久久免费| 国产精品偷伦视频观看了| 自线自在国产av| 妹子高潮喷水视频| 久久精品熟女亚洲av麻豆精品| 99国产综合亚洲精品| 日本wwww免费看| 亚洲成人手机| 国产精品1区2区在线观看. | 久久影院123| 高清欧美精品videossex| 蜜桃在线观看..| 亚洲免费av在线视频| 麻豆乱淫一区二区| 51午夜福利影视在线观看| 久久久久久久久免费视频了| 亚洲美女黄片视频| 国产熟女午夜一区二区三区| 97人妻天天添夜夜摸| 久久精品人人爽人人爽视色| 亚洲情色 制服丝袜| 午夜免费成人在线视频| 欧美日韩成人在线一区二区| 色播在线永久视频| 日韩欧美一区视频在线观看| av片东京热男人的天堂| 成人av一区二区三区在线看| 日韩大码丰满熟妇| 老司机深夜福利视频在线观看| 十八禁网站网址无遮挡| 高清在线国产一区| 国产淫语在线视频| 日韩一卡2卡3卡4卡2021年| 成人特级黄色片久久久久久久 | 啦啦啦视频在线资源免费观看| 国产高清激情床上av| 极品人妻少妇av视频| 又紧又爽又黄一区二区| 国产免费福利视频在线观看| 99久久99久久久精品蜜桃| 一进一出好大好爽视频| 午夜激情久久久久久久| 男女午夜视频在线观看| 黄色片一级片一级黄色片| 一边摸一边做爽爽视频免费| 久久久久久久久久久久大奶| 在线观看舔阴道视频| a级片在线免费高清观看视频| 欧美黄色淫秽网站| 纯流量卡能插随身wifi吗| 老鸭窝网址在线观看| 丝袜美足系列| www.999成人在线观看| 亚洲成国产人片在线观看| 精品国产超薄肉色丝袜足j| 国产免费福利视频在线观看| 亚洲中文日韩欧美视频| 日韩视频在线欧美| 最新在线观看一区二区三区| 99国产极品粉嫩在线观看| 制服诱惑二区| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 丁香六月欧美| 亚洲熟妇熟女久久| 9色porny在线观看| 日韩欧美一区二区三区在线观看 | 欧美国产精品一级二级三级| 在线观看一区二区三区激情| 日本vs欧美在线观看视频| 亚洲国产欧美日韩在线播放| 一级,二级,三级黄色视频| 咕卡用的链子| 午夜福利在线免费观看网站| 中文字幕高清在线视频| 美女高潮到喷水免费观看| 久久久精品94久久精品| 变态另类成人亚洲欧美熟女 | 黄色视频在线播放观看不卡| 看免费av毛片| 日韩欧美国产一区二区入口| 日本a在线网址| 欧美精品av麻豆av| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成77777在线视频| 午夜激情久久久久久久| av免费在线观看网站| 十八禁网站免费在线| 日本vs欧美在线观看视频| 成年版毛片免费区| 99re6热这里在线精品视频| 久久久精品国产亚洲av高清涩受| 欧美性长视频在线观看| 亚洲成人国产一区在线观看| 久久久国产一区二区| 精品亚洲成国产av| 性少妇av在线| 日韩欧美国产一区二区入口| 精品国产一区二区久久| 午夜精品国产一区二区电影| 久久香蕉激情| 欧美亚洲日本最大视频资源| 波多野结衣一区麻豆| 国产日韩欧美视频二区| 欧美亚洲 丝袜 人妻 在线| 757午夜福利合集在线观看| 国产男女超爽视频在线观看| 国产99久久九九免费精品| 蜜桃在线观看..| 欧美日韩亚洲国产一区二区在线观看 | 最新美女视频免费是黄的| 日本欧美视频一区| 狂野欧美激情性xxxx| 老熟女久久久| 精品国产一区二区久久| 天天躁夜夜躁狠狠躁躁| 操出白浆在线播放| 嫩草影视91久久| 精品国产乱码久久久久久小说| 亚洲熟女毛片儿| 精品福利永久在线观看| 亚洲午夜精品一区,二区,三区| 久久人妻熟女aⅴ| 精品卡一卡二卡四卡免费| 美女高潮到喷水免费观看| 日本av手机在线免费观看| 91麻豆精品激情在线观看国产 | 免费不卡黄色视频| 国产亚洲精品第一综合不卡| 色婷婷av一区二区三区视频| 色在线成人网| 18禁裸乳无遮挡动漫免费视频| 国产深夜福利视频在线观看| 国产麻豆69| 伊人久久大香线蕉亚洲五| 超色免费av| 老司机在亚洲福利影院| 大型av网站在线播放| 又黄又粗又硬又大视频| 狠狠婷婷综合久久久久久88av| 一区二区三区国产精品乱码| 亚洲精品国产精品久久久不卡| 青青草视频在线视频观看| 国产不卡av网站在线观看| 国产片内射在线| 久久亚洲精品不卡| 欧美国产精品一级二级三级| 91麻豆精品激情在线观看国产 | 精品欧美一区二区三区在线| 丝袜喷水一区| 成人国语在线视频| 巨乳人妻的诱惑在线观看| 人人妻人人澡人人看| 午夜精品国产一区二区电影| 黄色片一级片一级黄色片| 91字幕亚洲| 无人区码免费观看不卡 | 熟女少妇亚洲综合色aaa.| 亚洲国产毛片av蜜桃av| 熟女少妇亚洲综合色aaa.| 日韩成人在线观看一区二区三区| 色综合欧美亚洲国产小说| 国产又爽黄色视频| 国产精品 欧美亚洲| 国产日韩欧美在线精品| av网站在线播放免费| 国产淫语在线视频| 亚洲avbb在线观看| 国产精品香港三级国产av潘金莲| 国产片内射在线| 考比视频在线观看| 成人影院久久| 深夜精品福利| 两个人看的免费小视频| 国产伦理片在线播放av一区| 国产在线精品亚洲第一网站| 国产午夜精品久久久久久| 在线观看免费高清a一片| 成年人免费黄色播放视频| 久久影院123| 欧美老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 一区福利在线观看| 欧美另类亚洲清纯唯美| 国产精品久久久人人做人人爽| 亚洲欧洲精品一区二区精品久久久| 91麻豆精品激情在线观看国产 | 国产精品久久电影中文字幕 | 岛国毛片在线播放| 色综合婷婷激情| 亚洲av电影在线进入| 亚洲午夜理论影院| 国产av精品麻豆| 国产成人精品久久二区二区免费| 国产亚洲精品久久久久5区| 国产精品熟女久久久久浪| 男女下面插进去视频免费观看| 九色亚洲精品在线播放| 欧美精品一区二区大全| 久热爱精品视频在线9| 中国美女看黄片| 丰满饥渴人妻一区二区三| 麻豆国产av国片精品| 999精品在线视频| 国产极品粉嫩免费观看在线| 女人高潮潮喷娇喘18禁视频| 熟女少妇亚洲综合色aaa.| 一进一出好大好爽视频| 精品国内亚洲2022精品成人 | 亚洲av日韩在线播放| 色老头精品视频在线观看| 久久免费观看电影| 一区二区三区激情视频| 视频区图区小说| 国产在线一区二区三区精| e午夜精品久久久久久久| tube8黄色片| 亚洲综合色网址| 在线观看免费视频网站a站| 搡老乐熟女国产| www.999成人在线观看| 午夜福利视频精品| 亚洲欧美一区二区三区黑人| 国产精品亚洲一级av第二区| 日韩有码中文字幕| 丁香六月天网| 色尼玛亚洲综合影院| 91麻豆精品激情在线观看国产 | 高清在线国产一区| 天堂动漫精品| 搡老熟女国产l中国老女人| 满18在线观看网站| 中文亚洲av片在线观看爽 | 国产av国产精品国产| 久久国产精品大桥未久av| 夜夜骑夜夜射夜夜干| 黄色a级毛片大全视频| 男女无遮挡免费网站观看| 久久ye,这里只有精品| 免费在线观看完整版高清| 日日摸夜夜添夜夜添小说| 丝袜美腿诱惑在线| 九色亚洲精品在线播放| 国产熟女午夜一区二区三区| 午夜福利视频在线观看免费| 成人国产av品久久久| 菩萨蛮人人尽说江南好唐韦庄| 91成人精品电影| 成人三级做爰电影| 一本一本久久a久久精品综合妖精| 久久国产亚洲av麻豆专区| 黑丝袜美女国产一区| 热99国产精品久久久久久7| 电影成人av| 国产亚洲一区二区精品| 精品久久久久久电影网| 欧美黑人精品巨大| 色视频在线一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩精品网址| 成人手机av| 老熟女久久久| 欧美中文综合在线视频| 亚洲成av片中文字幕在线观看| 亚洲国产欧美网| 狠狠狠狠99中文字幕| 少妇 在线观看| 久久性视频一级片| 女人高潮潮喷娇喘18禁视频| 黑人猛操日本美女一级片| 久久久水蜜桃国产精品网| 亚洲中文日韩欧美视频| 精品免费久久久久久久清纯 | 久久中文看片网| 真人做人爱边吃奶动态| 51午夜福利影视在线观看| 午夜两性在线视频| 一进一出抽搐动态| 免费少妇av软件| 国产人伦9x9x在线观看| 午夜福利免费观看在线| 777久久人妻少妇嫩草av网站| 久久性视频一级片| 91成人精品电影| 国产片内射在线| 女性生殖器流出的白浆| 曰老女人黄片| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品第一综合不卡| 欧美大码av| 五月开心婷婷网| 不卡av一区二区三区| 免费人妻精品一区二区三区视频| 亚洲av成人不卡在线观看播放网| 2018国产大陆天天弄谢| 久久久水蜜桃国产精品网| 欧美 亚洲 国产 日韩一| 超色免费av| 十八禁高潮呻吟视频| 久久久精品94久久精品| 色94色欧美一区二区| 9热在线视频观看99| 成人黄色视频免费在线看| 极品少妇高潮喷水抽搐| av一本久久久久| 国产在线免费精品| 下体分泌物呈黄色| 久9热在线精品视频| 欧美精品高潮呻吟av久久| 中文字幕最新亚洲高清| 黄色视频,在线免费观看| 国产精品一区二区精品视频观看| bbb黄色大片| 一区二区日韩欧美中文字幕| 日日夜夜操网爽| 国产成人精品无人区| 999久久久精品免费观看国产| 精品国产一区二区久久| 日韩欧美一区二区三区在线观看 | 美女国产高潮福利片在线看| 男人操女人黄网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美日韩高清在线视频 | 亚洲三区欧美一区| 亚洲免费av在线视频| 亚洲第一青青草原| 一本大道久久a久久精品| 久久性视频一级片| 丁香欧美五月| 午夜免费鲁丝| 国产成人精品久久二区二区91| 国产一区二区三区视频了| 少妇的丰满在线观看| 亚洲,欧美精品.| 丰满少妇做爰视频| 国产精品偷伦视频观看了| 色在线成人网| 高潮久久久久久久久久久不卡| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩欧美视频二区| 日韩大片免费观看网站| 日韩欧美三级三区| 日韩免费高清中文字幕av| 成人手机av| 亚洲第一av免费看| 一级毛片精品| 久久亚洲真实| av天堂在线播放| 久久99一区二区三区| 美女视频免费永久观看网站| 亚洲国产欧美在线一区| 亚洲精品国产精品久久久不卡| 黄色视频在线播放观看不卡| 男人操女人黄网站| 国产精品香港三级国产av潘金莲| 在线观看免费高清a一片| 亚洲精品自拍成人| 精品亚洲乱码少妇综合久久| 中文字幕最新亚洲高清| 狠狠婷婷综合久久久久久88av| 久久婷婷成人综合色麻豆| 美女扒开内裤让男人捅视频| 午夜福利在线免费观看网站| 国产亚洲一区二区精品| 欧美成人免费av一区二区三区 | 99国产精品一区二区蜜桃av | 亚洲,欧美精品.| 男男h啪啪无遮挡| 日本a在线网址| 成人特级黄色片久久久久久久 | 水蜜桃什么品种好| 久久久久国内视频| 日韩一区二区三区影片| 国产麻豆69| 国产免费av片在线观看野外av| 亚洲精品av麻豆狂野| 在线看a的网站| 亚洲精品成人av观看孕妇| 精品熟女少妇八av免费久了| 精品亚洲成a人片在线观看| 久久久久精品国产欧美久久久| 捣出白浆h1v1| 午夜福利影视在线免费观看| 最近最新免费中文字幕在线| 日韩免费av在线播放| 欧美另类亚洲清纯唯美| 精品乱码久久久久久99久播| 亚洲av成人一区二区三| 夜夜爽天天搞| 欧美 日韩 精品 国产| 久久香蕉激情| 免费少妇av软件| 青青草视频在线视频观看| 日韩欧美免费精品| 一区二区三区激情视频| 国产成人免费观看mmmm| 精品一品国产午夜福利视频| av国产精品久久久久影院| 欧美日韩亚洲综合一区二区三区_| 欧美激情久久久久久爽电影 | 国产在线精品亚洲第一网站| 久久中文字幕人妻熟女| 美女国产高潮福利片在线看| 大香蕉久久成人网| 日韩免费av在线播放| 女人精品久久久久毛片| 考比视频在线观看|