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

    直升機(jī)-艦動(dòng)態(tài)配合的艦面流場(chǎng)環(huán)境研究綜述

    2014-03-11 06:23:49馬鴻儒孫文勝穆志韜李旭東
    裝備環(huán)境工程 2014年5期
    關(guān)鍵詞:尾流甲板艦船

    馬鴻儒,孫文勝,穆志韜,李旭東

    (海軍航空工程學(xué)院青島校區(qū),山東青島266041)

    直升機(jī)-艦動(dòng)態(tài)配合的艦面流場(chǎng)環(huán)境研究綜述

    馬鴻儒,孫文勝,穆志韜,李旭東

    (海軍航空工程學(xué)院青島校區(qū),山東青島266041)

    研究了制約直升機(jī)-艦動(dòng)態(tài)配合性能的首要原因——直升機(jī)的載艦在海上風(fēng)場(chǎng)、艦面特殊結(jié)構(gòu)、艦艇運(yùn)動(dòng)等因素作用下形成的非定常艦面紊流場(chǎng)。在總結(jié)早期直升機(jī)-艦動(dòng)態(tài)配合試驗(yàn)研究特點(diǎn)的基礎(chǔ)上,分析了艦面流場(chǎng)海上實(shí)測(cè)試驗(yàn)的方法及其發(fā)展;針對(duì)艦面流場(chǎng)模擬普遍采用的CFD技術(shù),分析了艦面流場(chǎng)CFD仿真的發(fā)展過(guò)程及現(xiàn)狀。評(píng)價(jià)了新型CFD方法——DES在處理大分離渦外流方面的優(yōu)勢(shì),指出了DES在艦面流場(chǎng)仿真領(lǐng)域的光明前景。

    直升機(jī);動(dòng)態(tài)配合;艦面流場(chǎng);DES

    艦載直升機(jī)與載艦的動(dòng)態(tài)配合性能是直升機(jī)-艦適配性的重要組成部分,意指直升機(jī)在近船或艦面操縱時(shí)直升機(jī)與艦的運(yùn)動(dòng)和其復(fù)雜流場(chǎng)相適應(yīng)完成任務(wù)的能力。機(jī)-艦動(dòng)態(tài)配合是一個(gè)復(fù)雜的空氣動(dòng)力學(xué)過(guò)程,包含大氣紊流和突風(fēng)的海面風(fēng)場(chǎng)、流經(jīng)上層建筑等結(jié)構(gòu)產(chǎn)生畸變的艦面流場(chǎng)、直升機(jī)旋翼尾流場(chǎng)以及它們之間的耦合作用構(gòu)成了直升機(jī)起飛、降落所處的非定常的流場(chǎng)環(huán)境。在直升機(jī)-艦動(dòng)態(tài)配合諸多的影響因素中,艦面紊流場(chǎng)對(duì)于直升機(jī)發(fā)揮其良好的操縱性和穩(wěn)定性,完成艦上起降任務(wù)的能力具有重要作用,因此艦面流場(chǎng)環(huán)境測(cè)試在動(dòng)態(tài)配合測(cè)試程序中居于首要地位。文中結(jié)合直升機(jī)-艦動(dòng)態(tài)配合研究的發(fā)展,對(duì)艦面流場(chǎng)測(cè)試及仿真研究的現(xiàn)狀進(jìn)行綜述,并對(duì)其未來(lái)發(fā)展方向作出展望。

    1 直升機(jī)-艦動(dòng)態(tài)配合試驗(yàn)

    1.1 國(guó)外研究

    在機(jī)-艦配合試驗(yàn)研究的早期階段,國(guó)外對(duì)特定的直升機(jī)/艦組合進(jìn)行了針對(duì)性極強(qiáng)的試驗(yàn),最主要的目的就是確定安全起降包絡(luò)線(xiàn)。直升機(jī)在艦上的安全起降包絡(luò)線(xiàn)是直升機(jī)與載艦配合執(zhí)行起降任務(wù)能力的體現(xiàn),一般包括以下幾個(gè)方面[1]:特定的直升機(jī)和艦組合;起飛和進(jìn)場(chǎng)的過(guò)程或航向;不同風(fēng)向允許的最大相對(duì)風(fēng)速;艦的最大橫搖角、縱搖角。如圖1所示。

    圖1 典型的直升機(jī)甲板操縱包絡(luò)線(xiàn)[3]Fig.1 Typical envelope of helicopter operating on the deck

    經(jīng)過(guò)這一階段的研究,一方面形成了艦載直升機(jī)艦上操縱的規(guī)范和對(duì)載艦的合格性要求;另一方面,經(jīng)過(guò)這一階段的發(fā)展,動(dòng)態(tài)配合試驗(yàn)法已經(jīng)成為了一門(mén)專(zhuān)門(mén)的學(xué)科,一些實(shí)驗(yàn)室比如美國(guó)海軍航空兵實(shí)驗(yàn)中心NATC和荷蘭海軍研究實(shí)驗(yàn)室NRL已經(jīng)將其標(biāo)準(zhǔn)化[2]。總結(jié)這一階段的試驗(yàn),可以得出一些研究的共性。

    在這些試驗(yàn)中,飛行器上通常都裝有電子儀表來(lái)記錄動(dòng)態(tài)配合時(shí)飛行器的運(yùn)動(dòng)和姿態(tài)參數(shù)。雖然試驗(yàn)方案具有多重目標(biāo),其當(dāng)時(shí)試驗(yàn)的重點(diǎn)都放在起飛/降落包線(xiàn)的研究和對(duì)已有包線(xiàn)的擴(kuò)展上。試驗(yàn)過(guò)程中要求直升機(jī)在某一特定的外界環(huán)境條件以及一定的直升機(jī)狀態(tài)下起飛和降落,之后,由飛行員利用飛行員等級(jí)評(píng)定標(biāo)準(zhǔn)對(duì)整個(gè)過(guò)程進(jìn)行評(píng)級(jí)。評(píng)定等級(jí)的依據(jù)是艦載環(huán)境下直升機(jī)的控制裕度、飛行品質(zhì)和性能,以及由此引起的飛行員工作負(fù)荷的變化。最后,以上述依據(jù)繪制各種環(huán)境下的晝夜間操縱包線(xiàn)。

    雖然上述海上試驗(yàn)在一定程度上確保了直升機(jī)艦上操縱安全的問(wèn)題,但是,要對(duì)此過(guò)程中影響因素進(jìn)行量化分析,還要進(jìn)行精確測(cè)定才行。總結(jié)以往飛行試驗(yàn)的結(jié)果,人們逐漸發(fā)現(xiàn),受海上紊流、艦面設(shè)備及上層建筑引起的渦流、艦面溫度場(chǎng)分布不均等干擾形成的復(fù)雜艦尾流對(duì)直升機(jī)操縱影響巨大,于是進(jìn)行了一系列艦尾流評(píng)估和測(cè)試。

    美國(guó)海軍研究生學(xué)校(NPS)對(duì)驅(qū)逐艦的艦面流場(chǎng)進(jìn)行了一系列的模型試驗(yàn)研究。其中,Gustav A. Anderson[4]在DD963模型艦上進(jìn)行了直升機(jī)特定飛行航線(xiàn)上尾流場(chǎng)的繪圖,并與之前該校進(jìn)行的低速風(fēng)洞可視化流場(chǎng)研究結(jié)果進(jìn)行對(duì)比,指出了該航線(xiàn)上非定常的艦尾流效應(yīng)。M.M.Rhoades[5]在另一型艦——AOR級(jí)上進(jìn)行了類(lèi)似試驗(yàn),給出了測(cè)得的飛行甲板上方流場(chǎng)示意圖,之后他又利用改進(jìn)的速度測(cè)量工具進(jìn)行試驗(yàn),在甲板上方發(fā)現(xiàn)了非常規(guī)的速度分量和強(qiáng)烈紊流的存在[6]。這些試驗(yàn)結(jié)果一方面證明了前期各國(guó)海上試驗(yàn)總結(jié)得出的規(guī)律——艦面上方存在紊亂的艦尾流場(chǎng)會(huì)給直升機(jī)操縱帶來(lái)顯著影響,同時(shí)也為之后的艦尾流仿真提供了可參照的依據(jù)。

    澳大利亞國(guó)防部也在20世紀(jì)90年代初對(duì)直升機(jī)載艦艦面流場(chǎng)特性進(jìn)行了研究。首先進(jìn)行的是FFG7驅(qū)逐艦甲板流場(chǎng)海上測(cè)量試驗(yàn)[7],記錄了艦船運(yùn)動(dòng)、參照風(fēng)速、以及甲板風(fēng)速等信息。A.M.Arney[8]對(duì)該試驗(yàn)實(shí)測(cè)數(shù)據(jù)的采集和處理的細(xì)節(jié)進(jìn)行了描述,針對(duì)數(shù)據(jù)采集時(shí)存在的艦上高頻發(fā)射裝置干擾問(wèn)題,提出了還原數(shù)據(jù)剔除干擾的方法。Lin-coln P.Erm[9]則對(duì)SH-60B/FFG-7機(jī)艦仿真程序所用的艦尾流模型進(jìn)行了詳細(xì)研究。該模型基于另一型艦的風(fēng)洞試驗(yàn)數(shù)據(jù),但針對(duì)FFG-7艦外形特點(diǎn)對(duì)數(shù)據(jù)進(jìn)行了修正,結(jié)果顯示風(fēng)洞中速度測(cè)量值與實(shí)測(cè)數(shù)據(jù)有很多不同,證明風(fēng)洞試驗(yàn)數(shù)據(jù)在一定程度上的不可靠性。

    1.2 國(guó)內(nèi)研究

    我國(guó)自20世紀(jì)80年代以來(lái),隨著海軍艦載直升機(jī)的逐步發(fā)展,開(kāi)始直升機(jī)-艦配合方面的研究工作。陳碧云結(jié)合國(guó)外艦載直升機(jī)發(fā)展的特點(diǎn),討論了直升機(jī)裝艦若干問(wèn)題[10],詳細(xì)指出了影響直升機(jī)/艦船動(dòng)態(tài)配合的有關(guān)因素,重點(diǎn)介紹了艦載直升機(jī)起降限制和飛行包線(xiàn)的確定方法。林河泉、傅百先等關(guān)于動(dòng)態(tài)配合的討論也多是在參考國(guó)外動(dòng)態(tài)配合試驗(yàn)的基礎(chǔ)上對(duì)飛行包線(xiàn)的確定進(jìn)行論述,據(jù)目前掌握資料尚無(wú)公開(kāi)進(jìn)行的機(jī)艦動(dòng)態(tài)配合試驗(yàn)。模型研究方面,劉航[11]在南航NH-2閉口回流式低速風(fēng)洞中進(jìn)行了艦?zāi)N擦鲌?chǎng)風(fēng)洞試驗(yàn)、流譜試驗(yàn)和地面模擬試驗(yàn)。試驗(yàn)中觀察到了美國(guó)海軍研究院試驗(yàn)中出現(xiàn)的流譜圖像和流動(dòng)特性,結(jié)果表明,飛行甲板上方的流動(dòng)極為復(fù)雜,表現(xiàn)出十分強(qiáng)烈的三維性和隨機(jī)性,特別是經(jīng)過(guò)機(jī)庫(kù)后發(fā)生了大規(guī)模分離,形成了一定的渦流區(qū)。由于該試驗(yàn)在國(guó)內(nèi)尚屬首次,船模中未考慮復(fù)雜的上層建筑,也未模擬艦船各方向運(yùn)動(dòng),因此流場(chǎng)模擬結(jié)果的精度還較低。

    2 艦尾流仿真計(jì)算

    隨著計(jì)算流體力學(xué)(CFD)技術(shù)的發(fā)展,動(dòng)態(tài)配合的計(jì)算機(jī)仿真為目前的直升機(jī)/艦系統(tǒng)開(kāi)發(fā)、測(cè)試和性能預(yù)測(cè)包括預(yù)測(cè)飛行包線(xiàn),提供了新的解決思路。

    在已經(jīng)進(jìn)行的直升機(jī)/艦動(dòng)態(tài)配合仿真系統(tǒng)中,最有代表性的當(dāng)屬美國(guó)的 JSHIP(Joint Shipboard Helicopter Integration Process)項(xiàng)目。該項(xiàng)目于1998年7月開(kāi)始實(shí)施,由美國(guó)軍方資助,對(duì)一系列機(jī)-艦系統(tǒng)進(jìn)行測(cè)試和評(píng)估,用以解決在聯(lián)合軍事行動(dòng)中陸軍、空軍直升機(jī)在艦上使用與載艦的配合能力,項(xiàng)目為期4~5年。要達(dá)到JSHIP中提出的增強(qiáng)聯(lián)合行動(dòng)中諸軍種直升機(jī)在艦艇上操縱能力的目標(biāo),項(xiàng)目主要是借助 DIMSS(Dynamic Interface Modeling and Simulation System)實(shí)現(xiàn)的。有了海上試驗(yàn)的經(jīng)驗(yàn)和初步結(jié)論,美國(guó)軍方一開(kāi)始就對(duì)DIMSS提出了明確要求,文獻(xiàn)[12]對(duì)DIMSS的目標(biāo)和結(jié)構(gòu)進(jìn)行了說(shuō)明,著重研究了整個(gè)模型與動(dòng)態(tài)配合試驗(yàn)相比所要達(dá)到的逼真度,給出了DI配合仿真逼真度的算法,如果某仿真系統(tǒng)達(dá)不到DIMSS逼真度的要求,文獻(xiàn)[13]介紹了一種基于逼真度反饋的方法改進(jìn)原仿真系統(tǒng)。

    2.1 艦面流場(chǎng)仿真的發(fā)展

    鑒于流場(chǎng)控制方程的直接求解在大多數(shù)情況下仍然難以實(shí)現(xiàn),通常需針對(duì)不同流體類(lèi)型、特征尺度、流場(chǎng)幾何特性等采用合適的求解方程及格式。

    一般地,描述艦面流場(chǎng)的流體力學(xué)方程可以表示為三維定常/非定常粘性不可壓N-S方程[14]。有關(guān)艦面流場(chǎng)的定常N-S方程的計(jì)算,國(guó)內(nèi)外進(jìn)行了大量的研究。Tsze C.Tai[15—16]對(duì)層流和紊流情況都進(jìn)行了考慮,描述了利用多區(qū)域、薄層N-S方程對(duì)DD-963艦尾流進(jìn)行仿真。所得的流場(chǎng)包括含有大面積自由渦的氣流分離區(qū),風(fēng)洞試驗(yàn)中表現(xiàn)的包括渦粘干擾的流場(chǎng)特征也有很好的體現(xiàn)。在某個(gè)特定的飛行路徑上計(jì)算結(jié)果與風(fēng)洞試驗(yàn)數(shù)據(jù)吻合良好,艦船尺寸越接近全尺寸,在風(fēng)速給定的情況下其對(duì)應(yīng)的雷諾數(shù)越大,而使艦船尾流的影響性相對(duì)減小,同時(shí)流場(chǎng)隨風(fēng)速的變化而產(chǎn)生的變化差異性減小。利物浦大學(xué)前期的研究[17—18]表明,使用定常CFD尾流模型與FLGHTLAB直升機(jī)模型相結(jié)合在模擬器上進(jìn)行甲板著艦?zāi)M,比海上實(shí)測(cè)執(zhí)行任務(wù)時(shí)飛行員的工作負(fù)荷有所增加。Sharma等[19]利用改進(jìn)的并行流體求解器PUMA計(jì)算了LPD-17艦周?chē)鸁o(wú)粘、穩(wěn)態(tài)、時(shí)準(zhǔn)的流場(chǎng)計(jì)算解。他們的穩(wěn)態(tài)解和風(fēng)洞試驗(yàn)的頻譜之間取得了很好的一致性。Polsky與Bruner[20—21]利用并行非結(jié)構(gòu)流動(dòng)求解器COBALT以不同的湍流模型Laminar N-S,MILES,k-ε和SST研究了LHA艦船尾流。結(jié)果表明,穩(wěn)態(tài)CFD解不能準(zhǔn)確預(yù)測(cè)紊流場(chǎng)的時(shí)均值,并發(fā)現(xiàn)湍流模型在計(jì)算中計(jì)入了太多耗散,因?yàn)橥糠蛛x格式有很大的耗散性。此外,Tattersall[22],Liu J[23],Guillot[24]以及Reddy[25]等人也分別用不同的湍流模型對(duì)N-S方程進(jìn)行了求解。

    2.2 國(guó)內(nèi)研究現(xiàn)狀

    我國(guó)研究人員也進(jìn)行了相關(guān)研究。非航空艦船方面:孫文勝利用有限元法求解二維定常粘性不可壓N-S方程,對(duì)甲板上方艦船縱向?qū)ΨQ(chēng)面內(nèi)氣流場(chǎng)進(jìn)行了仿真[26],表明了采用有限元法求解艦面流場(chǎng)的可行性。吳裕平[27]求解三維N-S方程對(duì)艦面流場(chǎng)了仿真,指出甲板后方機(jī)庫(kù)產(chǎn)生下洗氣流,將對(duì)直升機(jī)起降產(chǎn)生不利影響。為了滿(mǎn)足在實(shí)時(shí)仿真中直升機(jī)飛行仿真計(jì)算對(duì)于流場(chǎng)數(shù)據(jù)的實(shí)時(shí)性需要,孫文勝[28]在求解N-S方程的基礎(chǔ)上,利用神經(jīng)網(wǎng)絡(luò)的誤差反向傳播算法,可仿真計(jì)算甲板流場(chǎng)任一點(diǎn)的流場(chǎng)數(shù)據(jù)。結(jié)果顯示,神經(jīng)網(wǎng)絡(luò)對(duì)流場(chǎng)的模擬結(jié)果與實(shí)驗(yàn)結(jié)果及采用高精度網(wǎng)格的N-S方程模擬結(jié)果相比有良好的一致性,同時(shí)計(jì)算時(shí)間大為縮短。機(jī)艦動(dòng)態(tài)配合實(shí)時(shí)仿真及直升機(jī)該課目訓(xùn)練模擬器的開(kāi)發(fā)對(duì)于流場(chǎng)數(shù)據(jù)的實(shí)時(shí)性都有較高要求,在提高仿真計(jì)算硬件之外,該方法提供了一條切實(shí)可行的思路。

    航空艦船方面對(duì)LHA艦的研究較多,中國(guó)艦船研究設(shè)計(jì)中心進(jìn)行了LHA型艦縮比模型的艦面空氣流場(chǎng)多工況條件的CFD數(shù)值模擬計(jì)算[29],該數(shù)值模擬試驗(yàn)在定常條件下進(jìn)行。同時(shí)進(jìn)行了2種飛行甲板形式的艦船空氣流場(chǎng)特性比較研究[30],認(rèn)為具備采用全通平甲板的各類(lèi)技術(shù)條件時(shí),不推薦使用滑躍式飛行甲板的設(shè)計(jì)形式。曲飛等[31]通過(guò)CFD計(jì)算得到了LHA飛行甲板不同著艦點(diǎn)區(qū)域的定常氣流場(chǎng)數(shù)據(jù),并將CFX定常流數(shù)值模擬所得的速度時(shí)均值與風(fēng)洞試驗(yàn)結(jié)果及實(shí)船測(cè)量數(shù)據(jù)均進(jìn)行了對(duì)比。

    上述研究反映了利用CFD求解艦面流場(chǎng)的發(fā)展趨勢(shì)——從二維流到三維流,從有限元法到有限體積法等多種方法,求解網(wǎng)格尺寸的精細(xì),在求解中關(guān)注了越來(lái)越多的外部因素等。要達(dá)到為機(jī)艦動(dòng)態(tài)配合提供速度場(chǎng)的目的,還更高的求解速度。

    3 DES方法在艦面流場(chǎng)仿真中的應(yīng)用

    考慮湍流模型的精度,Hodge S.J.等在文獻(xiàn)[18]的基礎(chǔ)上引入了非定常模型[32],增加了整個(gè)仿真的逼真度,將計(jì)算結(jié)果應(yīng)用于直升機(jī)艦載起降飛行仿真時(shí)飛行員工作載荷有所增加。該研究首先采用的是簡(jiǎn)單驅(qū)逐艦?zāi)P?如圖2a所示),來(lái)生成CFD尾流數(shù)據(jù)。近期的研究使用了更切合實(shí)際的艦船模型——皇家海軍23型護(hù)衛(wèi)艦[33](如圖2b所示),使得模型精度進(jìn)一步增加。

    圖2 計(jì)算采用的艦幾何外形[33]Fig.2 Geometries used for computations[33]

    包括上述研究在內(nèi),很多利用CFD計(jì)算艦船尾流的文章都采用了無(wú)粘流方法或雷諾平均的N-S方法(RANS)和一些其他形式的湍流模型[34],然而,多數(shù)情況下RANS方法和標(biāo)準(zhǔn)湍流模型在應(yīng)用于陡壁體產(chǎn)生大分離尾流時(shí)給出的結(jié)果較差[35]。在高分離區(qū)域湍流渦各向異性,并且經(jīng)常表現(xiàn)出非定常特性,因此湍流模型作用下的結(jié)果與他們?cè)谠囼?yàn)條件下得到的參數(shù)相去甚遠(yuǎn)。由于RANS在處理大范圍分離流時(shí)的局限性,人們嘗試用大渦模擬(LES)方法解決問(wèn)題。大渦模擬的本質(zhì)是用非穩(wěn)態(tài)的N-S方程直接模擬大尺度渦運(yùn)動(dòng),而通過(guò)模型來(lái)模擬小尺度渦運(yùn)動(dòng),即大渦直接求解,小渦用模型求解。Camelli等[36]利用Smagorinsky湍流模型對(duì)LPD-17艦尾流進(jìn)行了LES模擬且加入了溫度場(chǎng)。這種情況下的近壁區(qū)網(wǎng)格處理是一個(gè)難題,另外LES多年來(lái)已經(jīng)被證明對(duì)于高雷諾數(shù)的外流是不適用的[37]。

    由此,綜合RANS法和LES方法優(yōu)點(diǎn)的分離渦仿真(DES)建模技術(shù)[37]應(yīng)運(yùn)而生。DES將分離渦區(qū)域湍流結(jié)構(gòu)網(wǎng)格進(jìn)行單獨(dú)的精確劃分,而不是對(duì)整個(gè)流場(chǎng)整體建模,從而克服了傳統(tǒng)湍流模型的限制,保證了尾流的不穩(wěn)定特性并因此得到了更切合實(shí)際的湍流譜。一些驗(yàn)證性研究表明DES生成的艦船尾流頻譜比非定常RANS方法更接近實(shí)驗(yàn)數(shù)據(jù)。DES也被應(yīng)用于多種高雷諾數(shù)外部氣動(dòng)流場(chǎng)中,經(jīng)證明在氣動(dòng)預(yù)測(cè)中,結(jié)果等于或在某些情況下好于RANS模型[38]。例如,James S Forrest[39]采用DES方法,對(duì)于兩種艦?zāi)P汀?jiǎn)單護(hù)衛(wèi)艦和23型護(hù)衛(wèi)艦,以一系列甲板風(fēng)向角產(chǎn)生了不同的流場(chǎng), CFD數(shù)據(jù)與風(fēng)洞數(shù)據(jù)和全尺寸海上試驗(yàn)數(shù)據(jù)有較好的一致性,提高了CFD模擬流場(chǎng)的精度,從而改善了直升機(jī)-艦動(dòng)態(tài)配合仿真精度。

    采用DES等更精確CFD方法的重要前提是需要?jiǎng)澐指?xì)的計(jì)算網(wǎng)格。隨著網(wǎng)格更加精細(xì),計(jì)算單元數(shù)量和計(jì)算尾流所需的時(shí)間將會(huì)增加。當(dāng)模型復(fù)雜時(shí),問(wèn)題更加嚴(yán)重,因?yàn)閺?fù)雜的外形需要更多的網(wǎng)格來(lái)精確描述。對(duì)CAD模型劃分網(wǎng)格時(shí),采用23型驅(qū)逐艦比SFS2所需的計(jì)算時(shí)間要高出一個(gè)數(shù)量級(jí)。

    國(guó)內(nèi)有關(guān)該方法的研究相對(duì)國(guó)外較晚,楊春蕾等[40]對(duì)RANS和DES方法的船體繞流模擬及其不確定度進(jìn)行了分析??諝鈩?dòng)力學(xué)領(lǐng)域,目前國(guó)內(nèi)對(duì)DES的工程應(yīng)用案例較少,少數(shù)的研究如應(yīng)用于汽車(chē)外流場(chǎng)模擬研究[41]和基本起落架的非定常分離流[42],已經(jīng)顯示出了在處理大范圍分離流方面良好的計(jì)算性能。隨著計(jì)算機(jī)計(jì)算能力的提高,該優(yōu)勢(shì)還將使其在空氣動(dòng)力學(xué)領(lǐng)域內(nèi)具有廣闊的應(yīng)用前景。

    4 結(jié)論

    無(wú)論是直升機(jī)-艦動(dòng)態(tài)配合海上試驗(yàn)還是動(dòng)態(tài)配合特性仿真研究抑或是訓(xùn)練模擬器的開(kāi)發(fā),艦面流場(chǎng)的測(cè)量和計(jì)算都是直升機(jī)-艦動(dòng)態(tài)配合的前提條件,為直升機(jī)飛行動(dòng)力學(xué)仿真提供了輸入。艦面流場(chǎng)數(shù)據(jù)的獲取方式也隨動(dòng)態(tài)配合研究的需要不斷變化,分析其發(fā)展歷程和趨勢(shì),可得以下結(jié)論。

    1)艦面流場(chǎng)的海上實(shí)測(cè)由于其耗資大、周期長(zhǎng),應(yīng)用正逐漸減少,但其可為仿真計(jì)算數(shù)據(jù)的對(duì)比提供重要參考依據(jù),具有不可替代的作用,因此隨數(shù)據(jù)采集和處理技術(shù)的進(jìn)步,該方法還將進(jìn)一步得到應(yīng)用。

    2)艦面流場(chǎng)的CFD仿真隨更高精度湍流模型和求解方法的出現(xiàn)得到了更廣泛的應(yīng)用,DES方法由于其自身的優(yōu)勢(shì)提高了大范圍分離渦仿真的精度,將在艦面外流場(chǎng)計(jì)算方面發(fā)揮更多的作用。

    3)機(jī)艦動(dòng)態(tài)配合實(shí)時(shí)仿真對(duì)艦面流場(chǎng)計(jì)算耗時(shí)提出了很高的要求,使得流場(chǎng)求解在滿(mǎn)足精度要求的前提下更加關(guān)注計(jì)算效率,同時(shí)差值方法的應(yīng)用也是一種可行的思路。

    [1] 林河泉,高正.直升機(jī)/艦動(dòng)態(tài)配合問(wèn)題的探討[C]//第五屆全國(guó)直升機(jī)年會(huì)論文集.1989:168—171.(余不詳) LIN Quan-he,GAO Zheng.Investigation on Helicopter/ ship Dynamic Interface[C]//Proceedings of the 5th National Helicopter Society.1989:168—171.(余不詳)

    [2] 中國(guó)直升機(jī)設(shè)計(jì)研究所.艦載直升機(jī)文集(上卷) [M].2006:94—95.(余不詳) China Helicopter Design Research Institute.Corpus of Shipboard Helicopter[M].2006:94—95.(余不詳)

    [3] WILKINSON C H,VANDER G M,ROSCOE M F.Modeling and Simulation of The Ship-helicopter Environment [C]//AIAA.2000:4583.(余不詳)

    [4] GUSTAV A A.Mapping the Airwake of a Model DD-963 Along Specific Heicopter Flight Paths[C]//Monterey, California,1989.(余不詳)

    [5] MARK M R.A Study of the Airwake Aerodynamics over the Flight Deck of an AOR Model Ship[C]//Monterey, California,1989.(余不詳)

    [6] RHOADES M M,HEALEY J Val.Flight Deck Aerodynamics of a Nonaviation Ship[J].Journal of Aircraft,1992,29 (4):619—626.

    [7] HOURIGAN D T,SUTTON C W,BIRD F J.Instrumentation For Airwake Measurements on the Flight Deck of a FFG-7[R].Australia,Department of Defence Defence Science and Technology Organisation Aeronautical Research Laboratory.1991:AR-006-623.

    [8] ARNEY A M.FFG-7 Ship Motion and Airwake Trial:Part 1-Data Processing Procedures[R].Australia,Department of Defence,Defence Science And Technology Organisation Aeronautical Research Laboratory,1994:AR-008-349.

    [9] LINCOLN P E.A Preliminary Study of the Airwake Model Used in an Existing SH-60B/FFG-7 Helicopter/Ship Simulation Program[R].Australia,Department of Defence, Defence Science And Technology Organisation Aeronautical Research Laboratory.1994:AR-008-644.

    [10]陳碧云.直升機(jī):艦載與海上作戰(zhàn)[M].北京:海洋出版社,2003. CHEN Bi-yun.Helicopter:Shipboard&Maritime Conduct [M].Beijing:Ocean Press,2003.

    [11]劉航,趙維義.直升機(jī)載艦空氣尾流特性試驗(yàn)方法[J].流體力學(xué)實(shí)驗(yàn)與測(cè)量,2002,16(4):57—63. LIU Hang,ZHAO Wei-yi.Experimental Methods on Air Wake Characteristics of Shipborne Helicopter[J].Experiments and Measurements in Fluid Mechanics,2002,16 (4):57—63.

    [12]WILKINSON C H,ROSCOE M F,VANDERVLIET G M. Determining Fidelity Standards for the Shipboard Launch and Recovery Task[C]//AIAA Modeling and Simulation Technologies Conference and Exhibit.Montreal,Canada, 2001:4583.

    [13]ADVANI S K,WILKINSON C H.Dynamic Interface Mod-eling and Simulation-A Unique Challenge[C]//The Royal Aeronautical Society Conference on Helicopter Flight Simulation.London,UK,2001.

    [14]呂紅.艦船周?chē)鷼饬鲌?chǎng)數(shù)值模擬方法及分布規(guī)律的研究[D].哈爾濱:哈爾濱工程大學(xué),2008. LYU Hong.The Study on the Numeric Simulation Method and the Distributed Law of Airflow Field Around Naval Ship[D].Harbin:Harbin Engineering University,2008.

    [15]TAI T C.Simulation of DD-963 Ship Airwake by Navier-Stokes Method[J].Journal of Aircraft.1995,32(6): 1399—1401.

    [16]TAI T C.Simulation and Analysis of LHD Ship Airwake by Navier-Stokes Method[C]//Presented at the NATO RTO Symposium on Fluid Dynamics Problems of Vehicles Operating near or in the Air-Sea Interface.Amsterdam, The Netherlands,1998.

    [17]ZAN S J.On Aerodynamic Modeling and Simulation of the Dynamic Interface,Proceedings of the Institution of Mechanical Engineers[J].Part G:Aerospace Engineering, 2005,219:393—410.

    [18]ROPER D M,HODGE S J,OWEN I,et al.Integrating CFD and Piloted Simulation to Quantify Ship-Helicopter Operating Limits[J].The Aeronautical Journal,2006,110 (1109):419—428.

    [19]SHARMA A,LONG L N.Airwake Simulations on an LPD 17 Ship[C]//AIAA Computational Fluid Dynamics Conference.Anaheim,CA,2001.

    [20]POLSKY S A,BRUNER C W S.Time-Accurate Computational Simulations of an LHA Ship Airwake[C]//18th AIAA Applied Aerodynamics Conference.Denver,CO,2000.

    [21]POLSKY S A.A Computational Study of Unsteady Ship Airwake[C]//40th AIAA Aerospace Sciences Meeting& Exhibit.Reno,Nevada,2002.

    [22]TATTERSALL P,ALBONE C M,SOLIMAN M M,et al. Prediction of Ship Air Wakes over Flight Decks Using CFD[C]//In AGARD Fluid Dynamics Symposium on Fluid dynamics of Vehicles Operating Near or In the Air-Sea Interface.Amsterdam,1998.

    [23]LIU J,LONG L N.Higher Order Accurate Ship Airwake Predictions for the Helicopter-Ship Interface Problem [C]//Presented at AHS 54th Annual Forum.Washington D C,1998.

    [24]GUILLOT M J,WALKER M A.Unsteady Analysis of the Air Wake over the LPD-17[C]//AIAA Applied Aerodynamics Conference and Exhibit.Denver,CO,2000.

    [25]REDDY K R,TOFFOLETTO R,JONES K R W.Numerical Simulation of Ship Airwake[J].Computers and Fluids,2000,129:451—465.

    [26]孫文勝.直升機(jī)艦載起降環(huán)境的氣動(dòng)特性仿真[J].系統(tǒng)仿真學(xué)報(bào),2003,15(2):254—256. SUN Wen-sheng.Aerodynamic Simulation of Helicopter Landing and Take-off Environment on Shipboard[J]. Journal of System Simulation,2003,15(2):254—256.

    [27]孫文勝,林明.基于神經(jīng)網(wǎng)絡(luò)模型的艦面流場(chǎng)仿真算法[J].中國(guó)工程科學(xué),2003,5(5):76—79. SUN Wen-sheng,LIN Ming.Simulation Algorithm of Flightdeck Airflow Based on Neural Network[J].Engineering Science,2003,5(5):76—79.

    [28]吳裕平.CFD在艦船甲板流場(chǎng)計(jì)算與特性研究中的應(yīng)用[J].直升機(jī)技術(shù),2006,147(3):1—4. WU Yu-ping.The Research and Applications of CFD to Calculation the Flow Field on Warship Dec[J].Helicopter Technique,2006,147(3):1—4.

    [29]陸超,姜治芳,王濤.不同工況條件對(duì)艦船艦面空氣流場(chǎng)的影響[J].艦船科學(xué)技術(shù),2009,31(9):38—42. LU Chao,JIANG Zhi-fang,WANG Tao.Influence of Different Airflow Situation for Ship Airwake[J].Ship Science and Technology,2009,31(9):38—42.

    [30]陸超,姜治芳,王濤.兩種飛行甲板形式的艦船空氣流場(chǎng)特性比較[J].艦船科學(xué)技術(shù),2009,31(7):29—31. LU Chao,JIANG Zhi-fang,WANG Tao.A Comparison of Ship's Airwakes With Different Flight Decks[J].Ship Science and Technology,2009,31(7):29—31.

    [31]曲飛,陸超,姜治芳,等.艦船艦面空氣流場(chǎng)的CFD數(shù)值模擬探討[J].中國(guó)艦船研究,2009,4(5):23—27. QU Fei,LU Chao,JIANG Zhi-fang,et al.CFD Numerical Simulation of Ship Air-wake[J].Chinese Journal of Ship Research,2009,4(5):23—27.

    [32]HODGE S J,PADFIELD G D,SOUTHWORTH M R.Helicopter-Ship Dynamic Interface Simulation:Fidelity at Low-Cost,Cutting Costs in Flight Simulation-Balancing Quality and Capability[C]//RaeS Conference.London, UK,2006.

    [33]FORREST J S,OWEN I,PADFIELD G D,et al.Detached-Eddy Simulation of Ship Airwakes for Piloted Helicopter Flight Simulation[C]//1st International Aerospace CFD Conference.Paris,France,2007.

    [34]REDDY K R,TOFFOLETTO R,JONES K R.Numerical Simulation of Ship Airwake[J].Computers&Fluids, 2000,29:451—465.

    [35]SPALART P R.Strategies for Turbulence Modelling and Simulations[J].International Journal of Heat and Fluid Flow,2000,21:252—263.

    [36]CAMELLI F E,SOTO O,LOHNER R,et al.Topside LPD17 Flow and Temperature Study with an Implicit Monolithic Scheme[C]//41st AIAA Aerospace Sciences Meeting& Exhibit.Reno,Nevada,2003.

    [37]SPALART P R,JOU W-H,STRELETS M,et al.Comments on the Feasibility of LES For Wings and on a Hybrid RANS/LES Approach[C]//1st AFOSR International Conference on DNS/LES.Ruston,LA,1997.

    [38]STRELETS M.Detached Eddy Simulation of Massively Separated Flows[C]//39th AIAA Aerospace Sciences Meeting and Exhibit.Reno,NV,2001.

    [39]JAMES S F,STEVEN J H,IEUAN O,et al.Towards Fully Simulated Ship-Helicopter Operating Limits[C]//The Importance of Ship Airwake Fidelity the American Helicopter Society 64th Annual Forum.Montréal,Canada,2008.

    [40]楊春蕾,朱仁傳,繆國(guó)平,等.基于RANS和DES法船體繞流模擬及不確定度分析[J].上海交通大學(xué)學(xué)報(bào), 2012,46(3):430—435. YANG Chun-lei,ZHU Ren-chuan,MOU Guo-ping,et al. Uncertainty Analysis in CFD for Flow Simulation Around Ship Using RANS and DES[J].Journal of Shanghai Jiaotong University,2012,46(3):430—435.

    [41]林鐵平.汽車(chē)外流場(chǎng)DES模擬研究[D].長(zhǎng)沙:湖南大學(xué),2010. LIN Tie-ping.Detached Eddy Simulation of Vehicle External Flow Field[D].Changsha:Hunan University,2010.

    [42]劉建.采用DES類(lèi)方法研究基本起落架非定常大范圍分離流動(dòng)[D].北京:清華大學(xué),2011. LIU Jian.Study of Unsteady and Massively Separated Flows Past Rudimentary Landing Gear Using DES-type Methods[D].Beijing:Tsinghua University,2011.

    Review on Shipboard Flow Study for Dynamic Interface of Helicopter and Ship

    MA Hong-ru,SUN Wen-sheng,MU Zhi-tao,LI Xu-dong
    (Naval Aeronautical Engineering Academy Qingdao Campus,Qingdao 266041,China)

    To study the main factor restraining the dynamic interface performance of helicopter and ship-an unsteady flow field formed on the helicopter carrier by sea breeze,special structure of the ship and ship movement.Based on the feature of dynamic interface test of helicopter and ship in early stage,the methods of shipboard flow measurement as well as its development were analyzed.As for the CFD technique-usually used in flow simulation,its development and status were discussed.The advantage of DES,a new method of CFD,in handling Large-separation external flow was introduced and its bright perspective in simulation of shipboard flow was pointed out.

    helicopter;dynamic interface;shipboard flow;DES

    MU Zhi-tao(1963—),Male,Ph.D.,Professor,Research focus:corrosion fatigue and fatigue life evaluation of materials.

    10.7643/issn.1672-9242.2014.05.015

    V212.4

    :A

    1672-9242(2014)05-0075-06

    2014-06-05;

    2014-07-17

    Received:2014-06-05;Revised:2014-07-17

    馬鴻儒(1988—),男,河北人,博士研究生,主要研究方向?yàn)楹Q蟓h(huán)境與飛行器性能研究。

    Biography:MA Hong-ru(1988—),Male,from Hebei,Doctorial student,Research focus:ocean environment&aircraft performance.

    穆志韜(1963—),男,博士,教授,主要研究方向?yàn)椴牧细g疲勞及壽命評(píng)估。

    猜你喜歡
    尾流甲板艦船
    艦船通信中的噪聲消除研究
    艦船測(cè)風(fēng)傳感器安裝位置數(shù)值仿真
    科考船木甲板安裝工藝
    殲-15艦載戰(zhàn)斗機(jī)在遼寧艦飛行甲板準(zhǔn)備起飛
    國(guó)防(2017年12期)2017-04-09 07:52:09
    飛機(jī)尾流的散射特性與探測(cè)技術(shù)綜述
    替代甲板結(jié)合負(fù)壓吸引治療修復(fù)甲床缺損
    艦船腐蝕預(yù)防與控制系統(tǒng)工程
    錐形流量計(jì)尾流流場(chǎng)分析
    甲板前置治療甲床部分缺損
    水面艦船風(fēng)尾流效應(yīng)減弱的模擬研究
    在线观看一区二区三区| 麻豆久久精品国产亚洲av| 中文在线观看免费www的网站| 亚洲va在线va天堂va国产| 丝袜美腿在线中文| 91久久精品国产一区二区三区| 亚洲熟妇中文字幕五十中出| 少妇被粗大猛烈的视频| 国产精品一区二区性色av| 人人妻人人澡人人爽人人夜夜 | 国产亚洲精品av在线| 日韩欧美国产在线观看| 色视频www国产| 18禁裸乳无遮挡免费网站照片| 91久久精品国产一区二区三区| 97超视频在线观看视频| 亚洲精品乱码久久久久久按摩| 熟妇人妻久久中文字幕3abv| 蜜桃久久精品国产亚洲av| 亚洲国产精品成人久久小说| 亚洲国产精品成人综合色| 日本猛色少妇xxxxx猛交久久| 日韩av不卡免费在线播放| 日本-黄色视频高清免费观看| 精品酒店卫生间| 国产精品人妻久久久久久| 国产色爽女视频免费观看| 干丝袜人妻中文字幕| 久久久久网色| 久久久久九九精品影院| 国产黄色小视频在线观看| 国国产精品蜜臀av免费| 爱豆传媒免费全集在线观看| 亚洲av中文av极速乱| 看十八女毛片水多多多| 国产一区二区三区综合在线观看 | 老女人水多毛片| 免费看日本二区| 三级毛片av免费| 性插视频无遮挡在线免费观看| 视频中文字幕在线观看| 精品午夜福利在线看| 男人舔女人下体高潮全视频| 边亲边吃奶的免费视频| 一个人看的www免费观看视频| 久久精品国产亚洲av涩爱| 99久久精品一区二区三区| 亚洲美女搞黄在线观看| 亚洲精品456在线播放app| 男人和女人高潮做爰伦理| 在线观看免费高清a一片| 搡老妇女老女人老熟妇| 久久久精品欧美日韩精品| 爱豆传媒免费全集在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品久久久久久婷婷小说| 两个人的视频大全免费| 在线观看av片永久免费下载| 久久精品人妻少妇| 国产一级毛片七仙女欲春2| 插逼视频在线观看| 久久久久网色| 精品久久久久久久久亚洲| 男女啪啪激烈高潮av片| 午夜日本视频在线| 身体一侧抽搐| 夫妻性生交免费视频一级片| 熟妇人妻久久中文字幕3abv| 高清视频免费观看一区二区 | 汤姆久久久久久久影院中文字幕 | 国产精品国产三级国产av玫瑰| 亚洲国产精品专区欧美| 99久久中文字幕三级久久日本| 免费观看a级毛片全部| 日韩av不卡免费在线播放| 99视频精品全部免费 在线| 在线天堂最新版资源| 伊人久久精品亚洲午夜| 日韩欧美国产在线观看| 熟妇人妻不卡中文字幕| 国产精品蜜桃在线观看| a级毛色黄片| 人妻少妇偷人精品九色| 女人十人毛片免费观看3o分钟| av在线天堂中文字幕| 色哟哟·www| 亚洲综合色惰| 国产亚洲5aaaaa淫片| 日日啪夜夜撸| 国产av码专区亚洲av| 日本黄大片高清| 亚洲色图av天堂| 欧美人与善性xxx| 80岁老熟妇乱子伦牲交| 国产单亲对白刺激| a级毛色黄片| 国产精品久久久久久精品电影小说 | 免费大片18禁| 噜噜噜噜噜久久久久久91| 国产欧美另类精品又又久久亚洲欧美| 国产精品一区二区三区四区久久| xxx大片免费视频| 一边亲一边摸免费视频| 亚洲欧洲日产国产| 亚洲内射少妇av| 国产精品爽爽va在线观看网站| 亚洲最大成人av| 中国国产av一级| 久久久久久国产a免费观看| 国产亚洲午夜精品一区二区久久 | 国产精品久久久久久精品电影| 免费播放大片免费观看视频在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲精品亚洲一区二区| 欧美日韩在线观看h| 美女脱内裤让男人舔精品视频| 永久免费av网站大全| 国产极品天堂在线| 日韩一本色道免费dvd| 亚洲欧美日韩无卡精品| 亚洲精品日韩在线中文字幕| 一区二区三区高清视频在线| 神马国产精品三级电影在线观看| 久久久亚洲精品成人影院| 亚洲av中文av极速乱| ponron亚洲| 高清午夜精品一区二区三区| 黄色一级大片看看| 精品酒店卫生间| 身体一侧抽搐| 国产乱人视频| 国产极品天堂在线| 美女大奶头视频| 美女大奶头视频| 亚洲美女视频黄频| 七月丁香在线播放| 国产精品爽爽va在线观看网站| 亚洲高清免费不卡视频| 日日干狠狠操夜夜爽| 91久久精品国产一区二区成人| 亚洲四区av| 亚洲精品影视一区二区三区av| 午夜爱爱视频在线播放| 国产亚洲精品久久久com| 男插女下体视频免费在线播放| 国产精品爽爽va在线观看网站| 亚洲激情五月婷婷啪啪| 免费大片黄手机在线观看| 天天一区二区日本电影三级| 高清午夜精品一区二区三区| 老司机影院毛片| 成年版毛片免费区| av.在线天堂| 99视频精品全部免费 在线| 一区二区三区高清视频在线| 午夜福利网站1000一区二区三区| 久久草成人影院| 特大巨黑吊av在线直播| 国产乱人偷精品视频| 精品久久久久久久人妻蜜臀av| a级一级毛片免费在线观看| 久久久欧美国产精品| 国产精品人妻久久久久久| 精品不卡国产一区二区三区| 寂寞人妻少妇视频99o| 国产一区二区三区av在线| 久久久久精品久久久久真实原创| 色综合亚洲欧美另类图片| 亚洲综合色惰| 一级二级三级毛片免费看| 一夜夜www| 大香蕉97超碰在线| 午夜精品在线福利| 欧美丝袜亚洲另类| 欧美成人精品欧美一级黄| 老女人水多毛片| 国产探花极品一区二区| 99视频精品全部免费 在线| 男人舔奶头视频| 久久久成人免费电影| 乱人视频在线观看| 成人午夜高清在线视频| 性插视频无遮挡在线免费观看| 禁无遮挡网站| 成年人午夜在线观看视频 | 国产综合懂色| 亚洲,欧美,日韩| 国产av码专区亚洲av| 欧美最新免费一区二区三区| 我的老师免费观看完整版| 亚洲精品自拍成人| 26uuu在线亚洲综合色| 亚洲综合色惰| 亚洲av福利一区| 久久综合国产亚洲精品| 久久99热这里只频精品6学生| videossex国产| 国产亚洲精品久久久com| a级毛色黄片| 国产中年淑女户外野战色| 免费看av在线观看网站| 亚洲一区高清亚洲精品| 欧美高清性xxxxhd video| 国产精品人妻久久久久久| 少妇高潮的动态图| 人妻制服诱惑在线中文字幕| 亚洲经典国产精华液单| 亚洲成人av在线免费| 亚洲精品亚洲一区二区| 深夜a级毛片| 国产精品国产三级国产av玫瑰| 欧美性猛交╳xxx乱大交人| 女人十人毛片免费观看3o分钟| 五月天丁香电影| 国产片特级美女逼逼视频| 亚洲精品456在线播放app| 国产亚洲av片在线观看秒播厂 | 国产精品不卡视频一区二区| 色综合色国产| 午夜福利在线在线| 久久久久久伊人网av| 激情 狠狠 欧美| 国产精品精品国产色婷婷| 好男人在线观看高清免费视频| 国产乱来视频区| a级毛色黄片| 精品久久久噜噜| 亚洲av二区三区四区| 街头女战士在线观看网站| 日韩三级伦理在线观看| 乱人视频在线观看| 久久午夜福利片| 一区二区三区免费毛片| 免费在线观看成人毛片| 成年av动漫网址| 丝袜喷水一区| 亚洲精品乱码久久久v下载方式| 午夜免费男女啪啪视频观看| 亚洲经典国产精华液单| 哪个播放器可以免费观看大片| 日韩av免费高清视频| videossex国产| 最近的中文字幕免费完整| 亚洲人成网站高清观看| 国产一区二区三区综合在线观看 | 天天躁夜夜躁狠狠久久av| 久久久久精品久久久久真实原创| 免费电影在线观看免费观看| 亚洲精品久久午夜乱码| 大香蕉97超碰在线| 亚洲av电影不卡..在线观看| 免费少妇av软件| 亚洲经典国产精华液单| 91午夜精品亚洲一区二区三区| 91av网一区二区| 少妇裸体淫交视频免费看高清| 最近2019中文字幕mv第一页| videos熟女内射| 人人妻人人澡人人爽人人夜夜 | 搡老妇女老女人老熟妇| 亚洲国产精品sss在线观看| 男女啪啪激烈高潮av片| 国产一区亚洲一区在线观看| 亚洲国产最新在线播放| 成人亚洲精品av一区二区| 小蜜桃在线观看免费完整版高清| 高清午夜精品一区二区三区| 国产成人免费观看mmmm| av卡一久久| 18禁动态无遮挡网站| 99re6热这里在线精品视频| 精品人妻偷拍中文字幕| 我的女老师完整版在线观看| 亚洲高清免费不卡视频| 国产在视频线在精品| 特级一级黄色大片| 亚洲熟女精品中文字幕| 亚洲国产欧美在线一区| 一个人免费在线观看电影| 一级黄片播放器| 精品久久久久久久人妻蜜臀av| 午夜福利网站1000一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 插阴视频在线观看视频| 狠狠精品人妻久久久久久综合| 亚洲无线观看免费| 精品一区二区三区人妻视频| 高清日韩中文字幕在线| 在线观看美女被高潮喷水网站| 韩国高清视频一区二区三区| 中国美白少妇内射xxxbb| 2022亚洲国产成人精品| 青春草亚洲视频在线观看| 国产视频内射| 99热这里只有是精品50| 欧美日韩一区二区视频在线观看视频在线 | 国产人妻一区二区三区在| 国产大屁股一区二区在线视频| 女人久久www免费人成看片| 草草在线视频免费看| 精品久久久久久久末码| 丝袜喷水一区| 三级国产精品欧美在线观看| 亚洲一级一片aⅴ在线观看| 老司机影院成人| 国产淫语在线视频| 国产乱人视频| 国产综合精华液| 久热久热在线精品观看| 国产精品一区二区在线观看99 | 嫩草影院入口| 久久精品国产鲁丝片午夜精品| 看十八女毛片水多多多| eeuss影院久久| 成人午夜高清在线视频| av线在线观看网站| 91狼人影院| 欧美日韩精品成人综合77777| 男女视频在线观看网站免费| 久久精品国产亚洲av天美| 欧美日韩综合久久久久久| xxx大片免费视频| 亚洲精品成人av观看孕妇| 国产伦在线观看视频一区| 美女主播在线视频| 国精品久久久久久国模美| 国产v大片淫在线免费观看| 成人亚洲精品一区在线观看 | 一本一本综合久久| 亚洲在线观看片| 国产爱豆传媒在线观看| 国产高清国产精品国产三级 | 黄片无遮挡物在线观看| 偷拍熟女少妇极品色| 免费黄网站久久成人精品| 夜夜看夜夜爽夜夜摸| 国产黄a三级三级三级人| 能在线免费看毛片的网站| 久久精品国产亚洲av天美| 九九在线视频观看精品| 国产精品国产三级国产av玫瑰| 亚洲国产精品成人久久小说| 国产精品人妻久久久久久| 美女主播在线视频| 国产精品久久视频播放| 国产一区有黄有色的免费视频 | 我的女老师完整版在线观看| 午夜福利视频1000在线观看| 日韩一区二区视频免费看| 汤姆久久久久久久影院中文字幕 | 亚洲精品成人av观看孕妇| 麻豆成人av视频| 国产一区二区在线观看日韩| 18禁裸乳无遮挡免费网站照片| 成人欧美大片| 2022亚洲国产成人精品| 久久精品综合一区二区三区| 欧美 日韩 精品 国产| 久久久久久久久久久免费av| 日韩欧美精品v在线| 女人十人毛片免费观看3o分钟| 最近最新中文字幕大全电影3| 亚洲精华国产精华液的使用体验| 国产精品一区二区在线观看99 | 久久国产乱子免费精品| 美女内射精品一级片tv| 国产免费又黄又爽又色| 亚洲精品亚洲一区二区| 久久精品久久精品一区二区三区| 亚洲精品国产av蜜桃| 熟女电影av网| 精品久久久久久久久av| 国产69精品久久久久777片| 国产成人freesex在线| 成人国产麻豆网| 国产人妻一区二区三区在| 亚洲性久久影院| 亚洲精品日韩av片在线观看| 嫩草影院入口| 免费少妇av软件| 久久99精品国语久久久| 18禁在线无遮挡免费观看视频| 直男gayav资源| 国产高潮美女av| 青青草视频在线视频观看| 春色校园在线视频观看| 美女被艹到高潮喷水动态| 亚洲精品亚洲一区二区| 亚洲精品久久午夜乱码| 免费黄频网站在线观看国产| 午夜日本视频在线| 性插视频无遮挡在线免费观看| 成人综合一区亚洲| 亚洲乱码一区二区免费版| 久久午夜福利片| 免费av观看视频| 久久久成人免费电影| videossex国产| 天堂av国产一区二区熟女人妻| 亚洲精品,欧美精品| 好男人视频免费观看在线| 国产色婷婷99| 国产亚洲5aaaaa淫片| 一区二区三区乱码不卡18| 赤兔流量卡办理| 久久精品国产自在天天线| 亚洲精品日本国产第一区| 欧美变态另类bdsm刘玥| 日韩制服骚丝袜av| 99热6这里只有精品| 岛国毛片在线播放| 亚洲18禁久久av| 狂野欧美白嫩少妇大欣赏| 三级经典国产精品| 夫妻午夜视频| 特大巨黑吊av在线直播| 纵有疾风起免费观看全集完整版 | 国产精品一二三区在线看| 久久久久久久久久人人人人人人| 淫秽高清视频在线观看| 直男gayav资源| 国产黄频视频在线观看| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| 亚洲伊人久久精品综合| 99热全是精品| 高清午夜精品一区二区三区| 乱人视频在线观看| av黄色大香蕉| 亚洲乱码一区二区免费版| 99久国产av精品国产电影| 国产亚洲精品av在线| 亚洲美女视频黄频| 天天躁日日操中文字幕| 亚洲国产av新网站| 日产精品乱码卡一卡2卡三| av卡一久久| 网址你懂的国产日韩在线| 天堂√8在线中文| 午夜免费观看性视频| 又粗又硬又长又爽又黄的视频| 免费黄网站久久成人精品| 午夜福利网站1000一区二区三区| 中文乱码字字幕精品一区二区三区 | 国产成人freesex在线| 欧美日韩视频高清一区二区三区二| 国产伦在线观看视频一区| 最近中文字幕高清免费大全6| 不卡视频在线观看欧美| 精品久久国产蜜桃| 尤物成人国产欧美一区二区三区| 一级毛片 在线播放| 老女人水多毛片| 国产成人精品婷婷| 日韩在线高清观看一区二区三区| 一级毛片久久久久久久久女| 不卡视频在线观看欧美| 一本久久精品| 男的添女的下面高潮视频| 亚洲欧美一区二区三区国产| 亚洲成人精品中文字幕电影| 成年av动漫网址| 亚洲一级一片aⅴ在线观看| 高清毛片免费看| 色吧在线观看| 午夜免费观看性视频| 激情五月婷婷亚洲| 亚洲,欧美,日韩| 日本免费a在线| 欧美不卡视频在线免费观看| 卡戴珊不雅视频在线播放| 三级男女做爰猛烈吃奶摸视频| 久久国内精品自在自线图片| 亚洲精品亚洲一区二区| 成人亚洲精品av一区二区| 一级毛片 在线播放| 欧美97在线视频| 日韩强制内射视频| 黄色配什么色好看| 最近手机中文字幕大全| 日本免费在线观看一区| 最近的中文字幕免费完整| 91aial.com中文字幕在线观看| 亚洲久久久久久中文字幕| 美女黄网站色视频| 国产成人午夜福利电影在线观看| 国产黄片视频在线免费观看| 精品国产一区二区三区久久久樱花 | 国内精品美女久久久久久| 最近的中文字幕免费完整| 婷婷色综合大香蕉| 久久久a久久爽久久v久久| 综合色av麻豆| 精品酒店卫生间| 午夜福利视频1000在线观看| 亚洲av国产av综合av卡| 中文字幕免费在线视频6| 99久久中文字幕三级久久日本| 成年女人在线观看亚洲视频 | 日本黄色片子视频| 草草在线视频免费看| 精品久久久噜噜| 中国国产av一级| 成年版毛片免费区| 国产老妇伦熟女老妇高清| 男人舔奶头视频| 午夜福利成人在线免费观看| 国产一区亚洲一区在线观看| 欧美日本视频| 久久6这里有精品| 不卡视频在线观看欧美| 国精品久久久久久国模美| 精品不卡国产一区二区三区| 麻豆成人午夜福利视频| 国产免费一级a男人的天堂| 亚洲四区av| 亚洲av中文av极速乱| 午夜精品国产一区二区电影 | 国产精品女同一区二区软件| 白带黄色成豆腐渣| av女优亚洲男人天堂| 国产一区二区三区av在线| 高清在线视频一区二区三区| 欧美另类一区| 欧美不卡视频在线免费观看| 能在线免费观看的黄片| 亚洲精品成人久久久久久| 最近最新中文字幕大全电影3| 男女国产视频网站| 女人久久www免费人成看片| 老司机影院成人| 亚洲成人精品中文字幕电影| 国产成人a区在线观看| 观看美女的网站| 全区人妻精品视频| 青春草国产在线视频| www.色视频.com| 欧美变态另类bdsm刘玥| 一个人观看的视频www高清免费观看| 精品一区二区三区人妻视频| 91精品国产九色| 又粗又硬又长又爽又黄的视频| 国产一区二区三区综合在线观看 | 午夜福利在线观看吧| 久久热精品热| 久久精品国产鲁丝片午夜精品| 国产亚洲午夜精品一区二区久久 | 九九久久精品国产亚洲av麻豆| 免费人成在线观看视频色| 久久精品久久精品一区二区三区| 国产一级毛片七仙女欲春2| 干丝袜人妻中文字幕| 在线免费观看的www视频| 中文资源天堂在线| 午夜老司机福利剧场| freevideosex欧美| 国产综合懂色| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看| 久久久久性生活片| 亚洲av男天堂| 欧美日本视频| 亚洲四区av| 少妇被粗大猛烈的视频| 成人欧美大片| 91精品伊人久久大香线蕉| 亚洲精品自拍成人| 久久精品国产鲁丝片午夜精品| 青春草亚洲视频在线观看| 99视频精品全部免费 在线| 麻豆精品久久久久久蜜桃| 国产精品一及| 不卡视频在线观看欧美| 国产精品一区www在线观看| 热99在线观看视频| 少妇的逼水好多| 一区二区三区免费毛片| 街头女战士在线观看网站| 国产精品麻豆人妻色哟哟久久 | 插逼视频在线观看| 天堂俺去俺来也www色官网 | 99久久人妻综合| 嫩草影院精品99| 国产精品久久久久久久电影| 晚上一个人看的免费电影| 九色成人免费人妻av| 成年版毛片免费区| 国产精品国产三级国产专区5o| 亚洲av男天堂| 伊人久久精品亚洲午夜| 麻豆成人午夜福利视频| 精品一区二区免费观看| 日日摸夜夜添夜夜添av毛片| 国产亚洲av嫩草精品影院| 亚洲精品日韩av片在线观看| 日韩电影二区| 天天躁日日操中文字幕| 久久99蜜桃精品久久| 亚洲在线自拍视频| 久久97久久精品| 在线观看av片永久免费下载| 亚洲精品一区蜜桃| 99九九线精品视频在线观看视频| 亚洲电影在线观看av| 精品99又大又爽又粗少妇毛片| a级毛片免费高清观看在线播放| 你懂的网址亚洲精品在线观看| 久久久久久久大尺度免费视频| 国产亚洲午夜精品一区二区久久 | 亚洲av福利一区| 精品酒店卫生间| 久久99热6这里只有精品| xxx大片免费视频| 欧美不卡视频在线免费观看| 久久99热6这里只有精品|