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

    考慮進(jìn)排氣影響的運(yùn)輸機(jī)增升構(gòu)型氣動(dòng)特性研究

    2014-09-24 13:45:37張菁張曉亮江奕廷
    航空兵器 2014年3期

    張菁+張曉亮+江奕廷

    摘 要:以某型運(yùn)輸機(jī)增升構(gòu)型為研究對(duì)象,通過(guò)數(shù)值模擬方法研究了發(fā)動(dòng)機(jī)進(jìn)排氣對(duì)全機(jī) 氣動(dòng)特性的影響。計(jì)算結(jié)果表明:在發(fā)動(dòng)機(jī)進(jìn)排氣因素影響下,全機(jī)最大升力系數(shù)明顯增加,失 速迎角有較大幅度延遲。通過(guò)對(duì)流場(chǎng)特性對(duì)比分析知:進(jìn)排氣因素不僅對(duì)短艙后方襟翼當(dāng)?shù)亓鲌?chǎng) 有較大改善,而且對(duì)主翼上表面流場(chǎng)以及平尾當(dāng)?shù)赜且灿酗@著影響?;谝陨戏治稣J(rèn)為,在翼 吊發(fā)動(dòng)機(jī)增升構(gòu)型設(shè)計(jì)過(guò)程中,進(jìn)排氣因素對(duì)各個(gè)部件當(dāng)?shù)亓鲌?chǎng)的影響需要納入設(shè)計(jì)考慮范圍。

    關(guān)鍵詞:增升系統(tǒng);氣動(dòng)特性;動(dòng)力影響;短距起降

    中圖分類(lèi)號(hào):V231.3 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1673-5048(2014)03-0028-04

    TheResearchonAerodynamicCharacteristicsofHighLift ConfigurationofTransportPlanewiththeEffectofEngineJet

    ZHANGJing1,ZHANGXiaoliang2,JIANGYiting2

    (1.ChinaAirborneMissileAcademy,Luoyang471009,China;2.NorthwesternPolytechnical University,Xian710072,China)

    Abstract:BythemethodofCFD(ComputationalFluidDynamic),theaerodynamiccharacteristics ofhighliftconfigurationwiththeeffectofenginejetisresearched.Theresultofnumericalsimulation demonstratesthatwiththeeffectofenginejet,themaximumliftcoefficientincreases,andthestallangle increases,longitudinalstabilityfactordecreasessignificantly.Conclusioncanbeobtainedbyanalyzingthe physicscharacteristicsofflowthattheenginejetnotonlychangeslocalflowfieldofflapaftertheengine, butalsochangesthelocalflowfieldontheuppersurfaceofnacelleandmainwing.Inviewoftheabove analysis,duringtheprocessofdesigninghighliftconfigurationofwinghangingengine,theeffectsofen ginejetonthelocalflowfieldofeachparthastobeconsidered.

    Keywords:highliftsystem;aerodynamiccharacteristics;dynamiceffect;shorttakingoffandland ing

    0 引 言

    由于受到起降場(chǎng)地等條件的限制,要求運(yùn)輸 機(jī)起飛、著陸滑跑距離短以及起飛、進(jìn)場(chǎng)和著陸速 度低。而能實(shí)現(xiàn)上述性能的關(guān)鍵即為全機(jī)起降時(shí) 具有很大的升力系數(shù)。以失速速度Vs=20m/s、翼載W/S=200kg/m2要求為例,在起降距離300m、 飛越15m高障礙物條件下,要求全機(jī)最大升力系 數(shù)CLmax≈7.4[1]。

    由于普通機(jī)械式增升裝置自身所具有的局限 性,在運(yùn)輸機(jī)起降時(shí)所提供升力系數(shù)有限。鑒于 此,當(dāng)代大型運(yùn)輸機(jī)有些開(kāi)始采取動(dòng)力增升技術(shù)。 動(dòng)力增升技術(shù)分為邊界層控制、噴氣襟翼、展向吹 氣、動(dòng)力轉(zhuǎn)向等。其中屬于邊界層控制技術(shù)的外吹 式襟翼直接利用翼吊發(fā)動(dòng)機(jī)排氣控制襟翼表面流動(dòng),較之于其他動(dòng)力增升技術(shù)簡(jiǎn)單而實(shí)用,并且對(duì) 高速巡航構(gòu)型影響不大。該項(xiàng)技術(shù)已非常成熟,特 別是經(jīng)過(guò)了YC-15驗(yàn)證機(jī)的飛行試驗(yàn)驗(yàn)證并在C -17大型軍用運(yùn)輸機(jī)上得到了正式應(yīng)用[2-3]。故 外吹式襟翼增升系統(tǒng)成為翼吊布局大型運(yùn)輸機(jī)的 最佳選擇。

    在采用外吹式襟翼動(dòng)力增升技術(shù)的前提下, 必須研究動(dòng)力對(duì)增升構(gòu)型氣動(dòng)特性的影響。國(guó)外 在此技術(shù)領(lǐng)域有著豐富的研究經(jīng)驗(yàn),并且取得了 很多研究成果,然而國(guó)外公開(kāi)發(fā)表的文獻(xiàn)中只公 布了研究結(jié)果以及所得出的結(jié)論[2-6],并沒(méi)有討論 動(dòng)力因素對(duì)氣動(dòng)特性產(chǎn)生影響的流動(dòng)機(jī)理。國(guó)內(nèi) 在該技術(shù)領(lǐng)域的研究主要集中于發(fā)動(dòng)機(jī)噴流對(duì)襟 翼的影響上,并沒(méi)有探究動(dòng)力因素對(duì)其他部件氣 動(dòng)特性的影響[7-9]。

    本文通過(guò)求解雷諾平均N-S方程(RANS)對(duì) 某型采用外吹式襟翼技術(shù)的大型運(yùn)輸機(jī)增升構(gòu)型 進(jìn)行數(shù)值模擬,研究了發(fā)動(dòng)機(jī)進(jìn)、排氣對(duì)整機(jī)氣動(dòng) 性能的影響并分析了相應(yīng)的流動(dòng)機(jī)理,得出了較 為有意義的結(jié)論。

    1 計(jì)算方法

    1.1 數(shù)值方法

    對(duì)于粘性起主導(dǎo)作用的增升構(gòu)型繞流問(wèn)題, 其流場(chǎng)伴隨著尾跡混合、流動(dòng)分離以及激波附面 層干擾等復(fù)雜流動(dòng)特性[10]。因此本文通過(guò)求解 RANS方程對(duì)增升構(gòu)型進(jìn)行數(shù)值模擬。采用有限體 積法求解該方程,空間離散格式為二階迎風(fēng)Roe 格式,時(shí)間推進(jìn)格式為L(zhǎng)U-SGS格式。

    綜合考慮計(jì)算效率和計(jì)算精度,流場(chǎng)模擬采 用Menter提出的k-ωSST湍流模型[11]。

    1.2 計(jì)算模型和計(jì)算網(wǎng)格

    本文研究對(duì)象為某型大型運(yùn)輸機(jī)增升構(gòu)型, 為描述方便將其命名為ModelA,具體幾何外形如 圖1所示。本文采用點(diǎn)對(duì)接結(jié)構(gòu)網(wǎng)格策略,在物面 附近應(yīng)用O型網(wǎng)格進(jìn)行加密以模擬真實(shí)流動(dòng)的附 面層,附面層第一層網(wǎng)格高度為1×10-5C(C為平 均氣動(dòng)弦長(zhǎng))量級(jí),保證較小的y+值。endprint

    1.3 數(shù)值方法可靠性驗(yàn)證

    針對(duì)本文所采用的數(shù)值模擬方法的可靠性驗(yàn) 證分為兩部分。

    首先采用一個(gè)無(wú)動(dòng)力增升構(gòu)型風(fēng)洞標(biāo)模應(yīng)用 數(shù)值模擬方法對(duì)于增升裝置流場(chǎng)模擬可靠性進(jìn)行 驗(yàn)證,該標(biāo)模為某大型客機(jī)三段增升裝置風(fēng)洞試 驗(yàn)?zāi)P?,該?gòu)型帶有大型翼吊發(fā)動(dòng)機(jī)短艙,襟翼與 主翼由滑軌艙連接,翼吊短艙內(nèi)側(cè)裝有渦流發(fā)生 器。在風(fēng)洞試驗(yàn)工況下(自由來(lái)流Ma=0.2,Re= 2.1×105),計(jì)算所得升力特性曲線(xiàn)與試驗(yàn)數(shù)據(jù)對(duì) 比如圖2所示,從圖中可以看出,計(jì)算所得升力特 性曲線(xiàn)與試驗(yàn)所得吻合良好,尤其是對(duì)最大升力 系數(shù)和失速迎角的捕捉較為準(zhǔn)確。計(jì)算所得失速 迎角比試驗(yàn)值小0.94°,最大升力系數(shù)小0.033。 該算例表明針對(duì)增升構(gòu)型流場(chǎng)的網(wǎng)格策略及數(shù)值 模擬方法是合理、可信的。

    其次對(duì)模擬發(fā)動(dòng)機(jī)噴流所采用數(shù)值方法可靠 性進(jìn)行驗(yàn)證。數(shù)值模擬中Ma=0.2,高度為海平面 標(biāo)準(zhǔn)大氣。驗(yàn)證模型為如圖3(a)所示的某型單獨(dú) 渦扇發(fā)動(dòng)機(jī)。圖3(b)分別為沿發(fā)動(dòng)機(jī)中軸線(xiàn)噴流總溫隨距離變化的數(shù)值模擬結(jié)果與試驗(yàn)值的對(duì)比。 從圖中可以看出在距核心機(jī)后緣端面最近及最遠(yuǎn) 處計(jì)算值和試驗(yàn)值符合很好,在中等距離處計(jì)算 值與試驗(yàn)值存在小的差異,但誤差范圍不超過(guò) 10%,該算例表明,本文針對(duì)發(fā)動(dòng)機(jī)帶動(dòng)力模型所 采用的數(shù)值模擬方法可以滿(mǎn)足本文研究工作的需 求。

    圖4所示為ModelA有無(wú)動(dòng)力構(gòu)型全機(jī)升力 系數(shù)對(duì)比,由圖可知,在動(dòng)力因素影響下該構(gòu)型最 大升力系數(shù)及失速迎角均有大幅提高。以該構(gòu)型 為研究對(duì)象,分別就發(fā)動(dòng)機(jī)進(jìn)、排氣對(duì)升力特性影 響機(jī)理展開(kāi)研究。

    2.1 發(fā)動(dòng)機(jī)進(jìn)氣對(duì)升力特性的影響

    圖5為20°攻角下動(dòng)力關(guān)閉和開(kāi)啟時(shí)短艙上表 面空間馬赫數(shù)云圖對(duì)比。由圖可知,在大迎角下, 進(jìn)氣時(shí)短艙上方存在較大范圍的流動(dòng)分離區(qū),而 無(wú)進(jìn)氣時(shí)該分離區(qū)消失。圖6為過(guò)短艙中軸線(xiàn)鉛垂 截面馬赫數(shù)云圖、駐點(diǎn)位置及流線(xiàn)對(duì)比圖。由圖可 得,在發(fā)動(dòng)機(jī)吸氣的作用下,短艙正前方流場(chǎng)呈現(xiàn) 收縮管道形態(tài),駐點(diǎn)位置上移,短艙上表面流場(chǎng)當(dāng) 地迎角降低,這對(duì)短艙上表面分離區(qū)域消失產(chǎn)生 有利影響。

    綜合以上分析,發(fā)動(dòng)機(jī)進(jìn)氣通過(guò)改變短艙、掛 架以及短艙后方主翼流場(chǎng)形態(tài),從而改善了短艙上表面及掛架內(nèi)側(cè)當(dāng)?shù)亓鲌?chǎng)。

    2.2 發(fā)動(dòng)機(jī)排氣對(duì)升力特性的影響

    本文所用構(gòu)型發(fā)動(dòng)機(jī)前伸量較大,發(fā)動(dòng)機(jī)射 流對(duì)短艙后方縫翼及主翼前緣流場(chǎng)有顯著影響。

    圖7為有無(wú)動(dòng)力構(gòu)型內(nèi)翼段空間流線(xiàn)及分離 區(qū)對(duì)比。由圖7(a)可知,流經(jīng)短艙及掛架內(nèi)側(cè)低 能量氣流通過(guò)掛架與內(nèi)縫翼間縫隙上翻至主翼上 表面,并且在機(jī)翼上表面沿展向內(nèi)外擴(kuò)展,形成了 范圍較大的低能量區(qū)域。沿展向向外擴(kuò)展的氣流 在流經(jīng)主翼后緣時(shí)發(fā)生物面分離。由圖7(b)可知, 在高速?lài)娏鞯囊渥饔孟拢鹘?jīng)短艙及掛架內(nèi)側(cè) 的低能量氣流并未流向主翼上表面,而是隨噴流 流向下游,降低了縫翼及主翼前緣當(dāng)?shù)赜?;短?上表面后緣流場(chǎng)有較大改善,這對(duì)短艙上表面前 緣分離的消失產(chǎn)生有利影響。流經(jīng)短艙及掛架外 側(cè)氣流通過(guò)掛架與外縫翼內(nèi)端面間隙流向主翼上 表面,并且沿翼展方向只向內(nèi)擴(kuò)展,從空間流線(xiàn)對(duì) 比可以看出,當(dāng)?shù)亓鲌?chǎng)流速較高,在流經(jīng)主翼后緣 處未發(fā)生物面分離。

    本文所研究的增升構(gòu)型采用外吹式動(dòng)力增升 系統(tǒng),噴流對(duì)增升構(gòu)型的襟翼影響較大。噴流對(duì)襟 翼的影響范圍有限,主要集中在短艙正后方區(qū)域。 鑒于此,選取兩個(gè)展向站位導(dǎo)流片及主襟翼的壓 力分布,與無(wú)動(dòng)力構(gòu)型的壓力分布結(jié)果進(jìn)行對(duì)比, 如圖8所示。從圖中可以看出,有進(jìn)排氣構(gòu)型的導(dǎo) 流片及襟翼下表面壓力普遍有所提高,上表面吸 力峰值也顯著提高,主襟翼上表面壓力普遍降低。 襟翼所提供升力大幅增加。圖9是截面馬赫數(shù)云圖對(duì)比,從圖中可以看出,有動(dòng)力構(gòu)型襟翼表面上方 的低馬赫數(shù)區(qū)明顯小于無(wú)動(dòng)力構(gòu)型。分析其原因: 一方面由于進(jìn)氣對(duì)主翼上表面流場(chǎng)的有利改善一 直延續(xù)到襟翼處,另一方面由于發(fā)動(dòng)機(jī)排氣通過(guò) 縫道對(duì)襟翼上表面的強(qiáng)加速作用,不僅使得襟翼 表面流速提高而且使得襟翼上方低馬赫數(shù)區(qū)大大 縮小,推遲襟翼表面分離的發(fā)生,增大失速迎角。

    (1)數(shù)值計(jì)算結(jié)果表明,在發(fā)動(dòng)機(jī)動(dòng)力因素的影響下,全機(jī)最大升力系數(shù)及失速迎角較無(wú)動(dòng) 力構(gòu)型均有大幅提高,升力特性曲線(xiàn)整體上移;

    (2)發(fā)動(dòng)機(jī)進(jìn)、排氣對(duì)短艙和主翼上表面、 襟翼當(dāng)?shù)亓鲌?chǎng)均有大幅改善,其綜合作用是提高 最大升力系數(shù)及失速迎角,升力特性曲線(xiàn)整體上 移;

    (3)為滿(mǎn)足運(yùn)輸機(jī)高效增升系統(tǒng)的設(shè)計(jì)要求, 要綜合考慮發(fā)動(dòng)機(jī)動(dòng)力因素對(duì)短艙、主翼、襟翼當(dāng) 地流場(chǎng)的影響。

    參考文獻(xiàn):

    [1]張錫金.飛機(jī)設(shè)計(jì)手冊(cè):6分冊(cè)[M].北京:航空工業(yè) 出版社,2003:85-86.

    [2]ThompsonJD.YC-15PowerPlantSystemDesignand Development[R].AIAA1974-973,1974.

    [3]GriffinHA,GonzalezLF,SrinivasK.Computational FluidDynamicsAnalysisofExternallyBlownFlapConfig urationforTransportAircraft[J].JournalofAircraft, 2008,45(1):172-184.

    [4]PetrovAV.AerodynamicsofAircraftwithWing-Pow eredLiftSystems[R].AIAA1993-4386,1993

    [5]SmithCC.EffectofEnginePositionandHighLiftDe vicesonAerodynamicCharacteristicsofanExternally BlownFlapJetSTOLModel[R].NASATRD-8478, 1977.

    [6]SlotnickJP,AnMY,MyskoSJ,etal.Navier-Stokes AnalysisofaHighWingTransportHigh-LiftConfigura tionwithExternallyBlownFlaps[R].AIAA2000-4219, 2000.

    [7]譚兆光,陳迎春,李杰,等.機(jī)體/動(dòng)力裝置一體化分 析中的動(dòng)力影響效應(yīng)數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào), 2009,24(8):1766-1772.

    [8]劉李濤,楊永,李喜樂(lè).外吹式動(dòng)力吹氣襟翼N-S方 程數(shù)值分析[J].航空計(jì)算技術(shù),2008,38(3):61- 64.

    [9]郭少杰,王豪杰,李杰.外吹式襟翼動(dòng)力增升數(shù)值模 擬方法研究[J].航空工程進(jìn)展,2010,1(1):49-54.

    [10]EliassonP,CatalanoP,LePapeMC,etal.Improved CFDPredictionsforHighLiftFlowsintheEuropeanPro jectEUROLIFTII[C]//25thAIAAAppliedAerody namicsConference,2007.

    [11]MenterFR.Two-EquationEddy-ViscosityTurbulence ModelsforEngineeringApplications[J].AIAAJournal, 1994,32(8):1598-1605.endprint

    1.3 數(shù)值方法可靠性驗(yàn)證

    針對(duì)本文所采用的數(shù)值模擬方法的可靠性驗(yàn) 證分為兩部分。

    首先采用一個(gè)無(wú)動(dòng)力增升構(gòu)型風(fēng)洞標(biāo)模應(yīng)用 數(shù)值模擬方法對(duì)于增升裝置流場(chǎng)模擬可靠性進(jìn)行 驗(yàn)證,該標(biāo)模為某大型客機(jī)三段增升裝置風(fēng)洞試 驗(yàn)?zāi)P停摌?gòu)型帶有大型翼吊發(fā)動(dòng)機(jī)短艙,襟翼與 主翼由滑軌艙連接,翼吊短艙內(nèi)側(cè)裝有渦流發(fā)生 器。在風(fēng)洞試驗(yàn)工況下(自由來(lái)流Ma=0.2,Re= 2.1×105),計(jì)算所得升力特性曲線(xiàn)與試驗(yàn)數(shù)據(jù)對(duì) 比如圖2所示,從圖中可以看出,計(jì)算所得升力特 性曲線(xiàn)與試驗(yàn)所得吻合良好,尤其是對(duì)最大升力 系數(shù)和失速迎角的捕捉較為準(zhǔn)確。計(jì)算所得失速 迎角比試驗(yàn)值小0.94°,最大升力系數(shù)小0.033。 該算例表明針對(duì)增升構(gòu)型流場(chǎng)的網(wǎng)格策略及數(shù)值 模擬方法是合理、可信的。

    其次對(duì)模擬發(fā)動(dòng)機(jī)噴流所采用數(shù)值方法可靠 性進(jìn)行驗(yàn)證。數(shù)值模擬中Ma=0.2,高度為海平面 標(biāo)準(zhǔn)大氣。驗(yàn)證模型為如圖3(a)所示的某型單獨(dú) 渦扇發(fā)動(dòng)機(jī)。圖3(b)分別為沿發(fā)動(dòng)機(jī)中軸線(xiàn)噴流總溫隨距離變化的數(shù)值模擬結(jié)果與試驗(yàn)值的對(duì)比。 從圖中可以看出在距核心機(jī)后緣端面最近及最遠(yuǎn) 處計(jì)算值和試驗(yàn)值符合很好,在中等距離處計(jì)算 值與試驗(yàn)值存在小的差異,但誤差范圍不超過(guò) 10%,該算例表明,本文針對(duì)發(fā)動(dòng)機(jī)帶動(dòng)力模型所 采用的數(shù)值模擬方法可以滿(mǎn)足本文研究工作的需 求。

    圖4所示為ModelA有無(wú)動(dòng)力構(gòu)型全機(jī)升力 系數(shù)對(duì)比,由圖可知,在動(dòng)力因素影響下該構(gòu)型最 大升力系數(shù)及失速迎角均有大幅提高。以該構(gòu)型 為研究對(duì)象,分別就發(fā)動(dòng)機(jī)進(jìn)、排氣對(duì)升力特性影 響機(jī)理展開(kāi)研究。

    2.1 發(fā)動(dòng)機(jī)進(jìn)氣對(duì)升力特性的影響

    圖5為20°攻角下動(dòng)力關(guān)閉和開(kāi)啟時(shí)短艙上表 面空間馬赫數(shù)云圖對(duì)比。由圖可知,在大迎角下, 進(jìn)氣時(shí)短艙上方存在較大范圍的流動(dòng)分離區(qū),而 無(wú)進(jìn)氣時(shí)該分離區(qū)消失。圖6為過(guò)短艙中軸線(xiàn)鉛垂 截面馬赫數(shù)云圖、駐點(diǎn)位置及流線(xiàn)對(duì)比圖。由圖可 得,在發(fā)動(dòng)機(jī)吸氣的作用下,短艙正前方流場(chǎng)呈現(xiàn) 收縮管道形態(tài),駐點(diǎn)位置上移,短艙上表面流場(chǎng)當(dāng) 地迎角降低,這對(duì)短艙上表面分離區(qū)域消失產(chǎn)生 有利影響。

    綜合以上分析,發(fā)動(dòng)機(jī)進(jìn)氣通過(guò)改變短艙、掛 架以及短艙后方主翼流場(chǎng)形態(tài),從而改善了短艙上表面及掛架內(nèi)側(cè)當(dāng)?shù)亓鲌?chǎng)。

    2.2 發(fā)動(dòng)機(jī)排氣對(duì)升力特性的影響

    本文所用構(gòu)型發(fā)動(dòng)機(jī)前伸量較大,發(fā)動(dòng)機(jī)射 流對(duì)短艙后方縫翼及主翼前緣流場(chǎng)有顯著影響。

    圖7為有無(wú)動(dòng)力構(gòu)型內(nèi)翼段空間流線(xiàn)及分離 區(qū)對(duì)比。由圖7(a)可知,流經(jīng)短艙及掛架內(nèi)側(cè)低 能量氣流通過(guò)掛架與內(nèi)縫翼間縫隙上翻至主翼上 表面,并且在機(jī)翼上表面沿展向內(nèi)外擴(kuò)展,形成了 范圍較大的低能量區(qū)域。沿展向向外擴(kuò)展的氣流 在流經(jīng)主翼后緣時(shí)發(fā)生物面分離。由圖7(b)可知, 在高速?lài)娏鞯囊渥饔孟?,流?jīng)短艙及掛架內(nèi)側(cè) 的低能量氣流并未流向主翼上表面,而是隨噴流 流向下游,降低了縫翼及主翼前緣當(dāng)?shù)赜?;短?上表面后緣流場(chǎng)有較大改善,這對(duì)短艙上表面前 緣分離的消失產(chǎn)生有利影響。流經(jīng)短艙及掛架外 側(cè)氣流通過(guò)掛架與外縫翼內(nèi)端面間隙流向主翼上 表面,并且沿翼展方向只向內(nèi)擴(kuò)展,從空間流線(xiàn)對(duì) 比可以看出,當(dāng)?shù)亓鲌?chǎng)流速較高,在流經(jīng)主翼后緣 處未發(fā)生物面分離。

    本文所研究的增升構(gòu)型采用外吹式動(dòng)力增升 系統(tǒng),噴流對(duì)增升構(gòu)型的襟翼影響較大。噴流對(duì)襟 翼的影響范圍有限,主要集中在短艙正后方區(qū)域。 鑒于此,選取兩個(gè)展向站位導(dǎo)流片及主襟翼的壓 力分布,與無(wú)動(dòng)力構(gòu)型的壓力分布結(jié)果進(jìn)行對(duì)比, 如圖8所示。從圖中可以看出,有進(jìn)排氣構(gòu)型的導(dǎo) 流片及襟翼下表面壓力普遍有所提高,上表面吸 力峰值也顯著提高,主襟翼上表面壓力普遍降低。 襟翼所提供升力大幅增加。圖9是截面馬赫數(shù)云圖對(duì)比,從圖中可以看出,有動(dòng)力構(gòu)型襟翼表面上方 的低馬赫數(shù)區(qū)明顯小于無(wú)動(dòng)力構(gòu)型。分析其原因: 一方面由于進(jìn)氣對(duì)主翼上表面流場(chǎng)的有利改善一 直延續(xù)到襟翼處,另一方面由于發(fā)動(dòng)機(jī)排氣通過(guò) 縫道對(duì)襟翼上表面的強(qiáng)加速作用,不僅使得襟翼 表面流速提高而且使得襟翼上方低馬赫數(shù)區(qū)大大 縮小,推遲襟翼表面分離的發(fā)生,增大失速迎角。

    (1)數(shù)值計(jì)算結(jié)果表明,在發(fā)動(dòng)機(jī)動(dòng)力因素的影響下,全機(jī)最大升力系數(shù)及失速迎角較無(wú)動(dòng) 力構(gòu)型均有大幅提高,升力特性曲線(xiàn)整體上移;

    (2)發(fā)動(dòng)機(jī)進(jìn)、排氣對(duì)短艙和主翼上表面、 襟翼當(dāng)?shù)亓鲌?chǎng)均有大幅改善,其綜合作用是提高 最大升力系數(shù)及失速迎角,升力特性曲線(xiàn)整體上 移;

    (3)為滿(mǎn)足運(yùn)輸機(jī)高效增升系統(tǒng)的設(shè)計(jì)要求, 要綜合考慮發(fā)動(dòng)機(jī)動(dòng)力因素對(duì)短艙、主翼、襟翼當(dāng) 地流場(chǎng)的影響。

    參考文獻(xiàn):

    [1]張錫金.飛機(jī)設(shè)計(jì)手冊(cè):6分冊(cè)[M].北京:航空工業(yè) 出版社,2003:85-86.

    [2]ThompsonJD.YC-15PowerPlantSystemDesignand Development[R].AIAA1974-973,1974.

    [3]GriffinHA,GonzalezLF,SrinivasK.Computational FluidDynamicsAnalysisofExternallyBlownFlapConfig urationforTransportAircraft[J].JournalofAircraft, 2008,45(1):172-184.

    [4]PetrovAV.AerodynamicsofAircraftwithWing-Pow eredLiftSystems[R].AIAA1993-4386,1993

    [5]SmithCC.EffectofEnginePositionandHighLiftDe vicesonAerodynamicCharacteristicsofanExternally BlownFlapJetSTOLModel[R].NASATRD-8478, 1977.

    [6]SlotnickJP,AnMY,MyskoSJ,etal.Navier-Stokes AnalysisofaHighWingTransportHigh-LiftConfigura tionwithExternallyBlownFlaps[R].AIAA2000-4219, 2000.

    [7]譚兆光,陳迎春,李杰,等.機(jī)體/動(dòng)力裝置一體化分 析中的動(dòng)力影響效應(yīng)數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào), 2009,24(8):1766-1772.

    [8]劉李濤,楊永,李喜樂(lè).外吹式動(dòng)力吹氣襟翼N-S方 程數(shù)值分析[J].航空計(jì)算技術(shù),2008,38(3):61- 64.

    [9]郭少杰,王豪杰,李杰.外吹式襟翼動(dòng)力增升數(shù)值模 擬方法研究[J].航空工程進(jìn)展,2010,1(1):49-54.

    [10]EliassonP,CatalanoP,LePapeMC,etal.Improved CFDPredictionsforHighLiftFlowsintheEuropeanPro jectEUROLIFTII[C]//25thAIAAAppliedAerody namicsConference,2007.

    [11]MenterFR.Two-EquationEddy-ViscosityTurbulence ModelsforEngineeringApplications[J].AIAAJournal, 1994,32(8):1598-1605.endprint

    1.3 數(shù)值方法可靠性驗(yàn)證

    針對(duì)本文所采用的數(shù)值模擬方法的可靠性驗(yàn) 證分為兩部分。

    首先采用一個(gè)無(wú)動(dòng)力增升構(gòu)型風(fēng)洞標(biāo)模應(yīng)用 數(shù)值模擬方法對(duì)于增升裝置流場(chǎng)模擬可靠性進(jìn)行 驗(yàn)證,該標(biāo)模為某大型客機(jī)三段增升裝置風(fēng)洞試 驗(yàn)?zāi)P?,該?gòu)型帶有大型翼吊發(fā)動(dòng)機(jī)短艙,襟翼與 主翼由滑軌艙連接,翼吊短艙內(nèi)側(cè)裝有渦流發(fā)生 器。在風(fēng)洞試驗(yàn)工況下(自由來(lái)流Ma=0.2,Re= 2.1×105),計(jì)算所得升力特性曲線(xiàn)與試驗(yàn)數(shù)據(jù)對(duì) 比如圖2所示,從圖中可以看出,計(jì)算所得升力特 性曲線(xiàn)與試驗(yàn)所得吻合良好,尤其是對(duì)最大升力 系數(shù)和失速迎角的捕捉較為準(zhǔn)確。計(jì)算所得失速 迎角比試驗(yàn)值小0.94°,最大升力系數(shù)小0.033。 該算例表明針對(duì)增升構(gòu)型流場(chǎng)的網(wǎng)格策略及數(shù)值 模擬方法是合理、可信的。

    其次對(duì)模擬發(fā)動(dòng)機(jī)噴流所采用數(shù)值方法可靠 性進(jìn)行驗(yàn)證。數(shù)值模擬中Ma=0.2,高度為海平面 標(biāo)準(zhǔn)大氣。驗(yàn)證模型為如圖3(a)所示的某型單獨(dú) 渦扇發(fā)動(dòng)機(jī)。圖3(b)分別為沿發(fā)動(dòng)機(jī)中軸線(xiàn)噴流總溫隨距離變化的數(shù)值模擬結(jié)果與試驗(yàn)值的對(duì)比。 從圖中可以看出在距核心機(jī)后緣端面最近及最遠(yuǎn) 處計(jì)算值和試驗(yàn)值符合很好,在中等距離處計(jì)算 值與試驗(yàn)值存在小的差異,但誤差范圍不超過(guò) 10%,該算例表明,本文針對(duì)發(fā)動(dòng)機(jī)帶動(dòng)力模型所 采用的數(shù)值模擬方法可以滿(mǎn)足本文研究工作的需 求。

    圖4所示為ModelA有無(wú)動(dòng)力構(gòu)型全機(jī)升力 系數(shù)對(duì)比,由圖可知,在動(dòng)力因素影響下該構(gòu)型最 大升力系數(shù)及失速迎角均有大幅提高。以該構(gòu)型 為研究對(duì)象,分別就發(fā)動(dòng)機(jī)進(jìn)、排氣對(duì)升力特性影 響機(jī)理展開(kāi)研究。

    2.1 發(fā)動(dòng)機(jī)進(jìn)氣對(duì)升力特性的影響

    圖5為20°攻角下動(dòng)力關(guān)閉和開(kāi)啟時(shí)短艙上表 面空間馬赫數(shù)云圖對(duì)比。由圖可知,在大迎角下, 進(jìn)氣時(shí)短艙上方存在較大范圍的流動(dòng)分離區(qū),而 無(wú)進(jìn)氣時(shí)該分離區(qū)消失。圖6為過(guò)短艙中軸線(xiàn)鉛垂 截面馬赫數(shù)云圖、駐點(diǎn)位置及流線(xiàn)對(duì)比圖。由圖可 得,在發(fā)動(dòng)機(jī)吸氣的作用下,短艙正前方流場(chǎng)呈現(xiàn) 收縮管道形態(tài),駐點(diǎn)位置上移,短艙上表面流場(chǎng)當(dāng) 地迎角降低,這對(duì)短艙上表面分離區(qū)域消失產(chǎn)生 有利影響。

    綜合以上分析,發(fā)動(dòng)機(jī)進(jìn)氣通過(guò)改變短艙、掛 架以及短艙后方主翼流場(chǎng)形態(tài),從而改善了短艙上表面及掛架內(nèi)側(cè)當(dāng)?shù)亓鲌?chǎng)。

    2.2 發(fā)動(dòng)機(jī)排氣對(duì)升力特性的影響

    本文所用構(gòu)型發(fā)動(dòng)機(jī)前伸量較大,發(fā)動(dòng)機(jī)射 流對(duì)短艙后方縫翼及主翼前緣流場(chǎng)有顯著影響。

    圖7為有無(wú)動(dòng)力構(gòu)型內(nèi)翼段空間流線(xiàn)及分離 區(qū)對(duì)比。由圖7(a)可知,流經(jīng)短艙及掛架內(nèi)側(cè)低 能量氣流通過(guò)掛架與內(nèi)縫翼間縫隙上翻至主翼上 表面,并且在機(jī)翼上表面沿展向內(nèi)外擴(kuò)展,形成了 范圍較大的低能量區(qū)域。沿展向向外擴(kuò)展的氣流 在流經(jīng)主翼后緣時(shí)發(fā)生物面分離。由圖7(b)可知, 在高速?lài)娏鞯囊渥饔孟?,流?jīng)短艙及掛架內(nèi)側(cè) 的低能量氣流并未流向主翼上表面,而是隨噴流 流向下游,降低了縫翼及主翼前緣當(dāng)?shù)赜?;短?上表面后緣流場(chǎng)有較大改善,這對(duì)短艙上表面前 緣分離的消失產(chǎn)生有利影響。流經(jīng)短艙及掛架外 側(cè)氣流通過(guò)掛架與外縫翼內(nèi)端面間隙流向主翼上 表面,并且沿翼展方向只向內(nèi)擴(kuò)展,從空間流線(xiàn)對(duì) 比可以看出,當(dāng)?shù)亓鲌?chǎng)流速較高,在流經(jīng)主翼后緣 處未發(fā)生物面分離。

    本文所研究的增升構(gòu)型采用外吹式動(dòng)力增升 系統(tǒng),噴流對(duì)增升構(gòu)型的襟翼影響較大。噴流對(duì)襟 翼的影響范圍有限,主要集中在短艙正后方區(qū)域。 鑒于此,選取兩個(gè)展向站位導(dǎo)流片及主襟翼的壓 力分布,與無(wú)動(dòng)力構(gòu)型的壓力分布結(jié)果進(jìn)行對(duì)比, 如圖8所示。從圖中可以看出,有進(jìn)排氣構(gòu)型的導(dǎo) 流片及襟翼下表面壓力普遍有所提高,上表面吸 力峰值也顯著提高,主襟翼上表面壓力普遍降低。 襟翼所提供升力大幅增加。圖9是截面馬赫數(shù)云圖對(duì)比,從圖中可以看出,有動(dòng)力構(gòu)型襟翼表面上方 的低馬赫數(shù)區(qū)明顯小于無(wú)動(dòng)力構(gòu)型。分析其原因: 一方面由于進(jìn)氣對(duì)主翼上表面流場(chǎng)的有利改善一 直延續(xù)到襟翼處,另一方面由于發(fā)動(dòng)機(jī)排氣通過(guò) 縫道對(duì)襟翼上表面的強(qiáng)加速作用,不僅使得襟翼 表面流速提高而且使得襟翼上方低馬赫數(shù)區(qū)大大 縮小,推遲襟翼表面分離的發(fā)生,增大失速迎角。

    (1)數(shù)值計(jì)算結(jié)果表明,在發(fā)動(dòng)機(jī)動(dòng)力因素的影響下,全機(jī)最大升力系數(shù)及失速迎角較無(wú)動(dòng) 力構(gòu)型均有大幅提高,升力特性曲線(xiàn)整體上移;

    (2)發(fā)動(dòng)機(jī)進(jìn)、排氣對(duì)短艙和主翼上表面、 襟翼當(dāng)?shù)亓鲌?chǎng)均有大幅改善,其綜合作用是提高 最大升力系數(shù)及失速迎角,升力特性曲線(xiàn)整體上 移;

    (3)為滿(mǎn)足運(yùn)輸機(jī)高效增升系統(tǒng)的設(shè)計(jì)要求, 要綜合考慮發(fā)動(dòng)機(jī)動(dòng)力因素對(duì)短艙、主翼、襟翼當(dāng) 地流場(chǎng)的影響。

    參考文獻(xiàn):

    [1]張錫金.飛機(jī)設(shè)計(jì)手冊(cè):6分冊(cè)[M].北京:航空工業(yè) 出版社,2003:85-86.

    [2]ThompsonJD.YC-15PowerPlantSystemDesignand Development[R].AIAA1974-973,1974.

    [3]GriffinHA,GonzalezLF,SrinivasK.Computational FluidDynamicsAnalysisofExternallyBlownFlapConfig urationforTransportAircraft[J].JournalofAircraft, 2008,45(1):172-184.

    [4]PetrovAV.AerodynamicsofAircraftwithWing-Pow eredLiftSystems[R].AIAA1993-4386,1993

    [5]SmithCC.EffectofEnginePositionandHighLiftDe vicesonAerodynamicCharacteristicsofanExternally BlownFlapJetSTOLModel[R].NASATRD-8478, 1977.

    [6]SlotnickJP,AnMY,MyskoSJ,etal.Navier-Stokes AnalysisofaHighWingTransportHigh-LiftConfigura tionwithExternallyBlownFlaps[R].AIAA2000-4219, 2000.

    [7]譚兆光,陳迎春,李杰,等.機(jī)體/動(dòng)力裝置一體化分 析中的動(dòng)力影響效應(yīng)數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào), 2009,24(8):1766-1772.

    [8]劉李濤,楊永,李喜樂(lè).外吹式動(dòng)力吹氣襟翼N-S方 程數(shù)值分析[J].航空計(jì)算技術(shù),2008,38(3):61- 64.

    [9]郭少杰,王豪杰,李杰.外吹式襟翼動(dòng)力增升數(shù)值模 擬方法研究[J].航空工程進(jìn)展,2010,1(1):49-54.

    [10]EliassonP,CatalanoP,LePapeMC,etal.Improved CFDPredictionsforHighLiftFlowsintheEuropeanPro jectEUROLIFTII[C]//25thAIAAAppliedAerody namicsConference,2007.

    [11]MenterFR.Two-EquationEddy-ViscosityTurbulence ModelsforEngineeringApplications[J].AIAAJournal, 1994,32(8):1598-1605.endprint

    在线播放国产精品三级| 色婷婷久久久亚洲欧美| av中文乱码字幕在线| 他把我摸到了高潮在线观看| 18禁美女被吸乳视频| 国产国语露脸激情在线看| 国产成人啪精品午夜网站| 国产激情欧美一区二区| 麻豆久久精品国产亚洲av| 日本欧美视频一区| 国产精品久久久av美女十八| 女生性感内裤真人,穿戴方法视频| 亚洲人成电影观看| 成人亚洲精品av一区二区| 窝窝影院91人妻| 777久久人妻少妇嫩草av网站| 好男人电影高清在线观看| 欧美成人午夜精品| 免费在线观看完整版高清| 香蕉丝袜av| 两性夫妻黄色片| 国产成人av教育| 久久精品91蜜桃| 麻豆久久精品国产亚洲av| 精品国产一区二区三区四区第35| 欧美在线一区亚洲| 18美女黄网站色大片免费观看| 免费在线观看亚洲国产| 窝窝影院91人妻| 不卡av一区二区三区| 午夜福利18| 亚洲av美国av| 极品教师在线免费播放| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产亚洲av麻豆专区| 国产精品一区二区精品视频观看| 搡老岳熟女国产| 亚洲性夜色夜夜综合| 香蕉国产在线看| 午夜福利成人在线免费观看| 国产成+人综合+亚洲专区| 日韩中文字幕欧美一区二区| 91麻豆精品激情在线观看国产| 俄罗斯特黄特色一大片| 精品熟女少妇八av免费久了| 日韩欧美国产在线观看| 999精品在线视频| 亚洲国产看品久久| 国产伦人伦偷精品视频| 亚洲男人的天堂狠狠| 亚洲av片天天在线观看| 丝袜在线中文字幕| 亚洲狠狠婷婷综合久久图片| 国产欧美日韩一区二区精品| 最新美女视频免费是黄的| 久久国产精品影院| 国产不卡一卡二| 一二三四社区在线视频社区8| 国产高清有码在线观看视频 | 一边摸一边抽搐一进一出视频| 亚洲人成网站在线播放欧美日韩| 日本精品一区二区三区蜜桃| 一区二区三区精品91| 久久草成人影院| 欧美国产精品va在线观看不卡| 亚洲全国av大片| 国产一卡二卡三卡精品| 51午夜福利影视在线观看| 日韩有码中文字幕| 淫妇啪啪啪对白视频| 成人永久免费在线观看视频| 欧美黑人精品巨大| 啪啪无遮挡十八禁网站| 可以免费在线观看a视频的电影网站| 一本综合久久免费| 午夜福利高清视频| 99在线人妻在线中文字幕| 欧美人与性动交α欧美精品济南到| 欧美中文综合在线视频| 激情在线观看视频在线高清| 村上凉子中文字幕在线| 一个人观看的视频www高清免费观看 | 国产黄a三级三级三级人| 他把我摸到了高潮在线观看| 不卡一级毛片| 亚洲欧美精品综合一区二区三区| 男女下面进入的视频免费午夜 | 亚洲中文av在线| 免费人成视频x8x8入口观看| 老司机福利观看| 一边摸一边抽搐一进一出视频| 国产亚洲精品第一综合不卡| 亚洲精品美女久久av网站| 久久精品91无色码中文字幕| 国产成人av激情在线播放| 国产精品99久久99久久久不卡| 一进一出抽搐gif免费好疼| 国产三级黄色录像| 18禁黄网站禁片午夜丰满| av天堂久久9| 在线观看www视频免费| 日韩 欧美 亚洲 中文字幕| 精品国产乱码久久久久久男人| 亚洲成av人片免费观看| 日本三级黄在线观看| 成人永久免费在线观看视频| 巨乳人妻的诱惑在线观看| 亚洲国产精品成人综合色| 少妇的丰满在线观看| 怎么达到女性高潮| 亚洲自拍偷在线| 女性生殖器流出的白浆| 无人区码免费观看不卡| 又黄又爽又免费观看的视频| 亚洲国产看品久久| 一级,二级,三级黄色视频| 国产aⅴ精品一区二区三区波| 又黄又爽又免费观看的视频| 国产av一区二区精品久久| 一本久久中文字幕| 一级毛片女人18水好多| 自线自在国产av| 亚洲熟妇熟女久久| 国产成人av教育| 黄色毛片三级朝国网站| 日本免费一区二区三区高清不卡 | x7x7x7水蜜桃| 国产精品98久久久久久宅男小说| 国产免费av片在线观看野外av| 免费女性裸体啪啪无遮挡网站| 美女 人体艺术 gogo| 中文字幕久久专区| 波多野结衣高清无吗| 人人妻人人爽人人添夜夜欢视频| 热99re8久久精品国产| 久久久国产成人精品二区| 黄色毛片三级朝国网站| 亚洲精品一区av在线观看| 国产成人精品在线电影| 亚洲男人的天堂狠狠| 一边摸一边抽搐一进一出视频| 国产真人三级小视频在线观看| 啦啦啦 在线观看视频| 国产av又大| 在线播放国产精品三级| av免费在线观看网站| 男人舔女人下体高潮全视频| 在线观看免费午夜福利视频| 欧美乱码精品一区二区三区| 国产又爽黄色视频| 制服诱惑二区| 99久久99久久久精品蜜桃| 乱人伦中国视频| 琪琪午夜伦伦电影理论片6080| 男人操女人黄网站| 国产亚洲精品av在线| 国产成人系列免费观看| 一夜夜www| 国产精品亚洲av一区麻豆| 丁香六月欧美| 97人妻精品一区二区三区麻豆 | av有码第一页| 国产精品秋霞免费鲁丝片| 久久精品亚洲熟妇少妇任你| 国产免费av片在线观看野外av| 亚洲国产欧美一区二区综合| 精品高清国产在线一区| 亚洲欧美日韩无卡精品| 国产精品 国内视频| 久久久久久久久免费视频了| 亚洲专区中文字幕在线| 免费看美女性在线毛片视频| av天堂久久9| 动漫黄色视频在线观看| 国产成年人精品一区二区| 日本撒尿小便嘘嘘汇集6| 丝袜美足系列| 巨乳人妻的诱惑在线观看| 一区二区三区精品91| 我的亚洲天堂| 亚洲伊人色综图| cao死你这个sao货| 国内久久婷婷六月综合欲色啪| 老熟妇乱子伦视频在线观看| 成人亚洲精品一区在线观看| 成年女人毛片免费观看观看9| 精品久久久精品久久久| 亚洲激情在线av| 高清毛片免费观看视频网站| 欧美一级a爱片免费观看看 | 欧美最黄视频在线播放免费| 波多野结衣av一区二区av| 久久中文字幕人妻熟女| 久久久久久久久免费视频了| 亚洲精品一区av在线观看| 国产成年人精品一区二区| 桃色一区二区三区在线观看| 一本综合久久免费| 亚洲激情在线av| 老司机午夜十八禁免费视频| 亚洲午夜精品一区,二区,三区| 国产不卡一卡二| 日韩高清综合在线| 真人做人爱边吃奶动态| 久久亚洲精品不卡| a在线观看视频网站| 亚洲欧美日韩高清在线视频| 久久国产精品男人的天堂亚洲| 免费人成视频x8x8入口观看| 脱女人内裤的视频| 一进一出好大好爽视频| 后天国语完整版免费观看| 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 美女高潮喷水抽搐中文字幕| 久久久久久久午夜电影| 午夜亚洲福利在线播放| 一边摸一边抽搐一进一小说| 日韩大码丰满熟妇| 国产欧美日韩一区二区精品| 亚洲欧美精品综合久久99| 欧美最黄视频在线播放免费| 欧美绝顶高潮抽搐喷水| 香蕉久久夜色| 欧美成人午夜精品| 午夜福利18| 久久天堂一区二区三区四区| or卡值多少钱| 制服诱惑二区| 老司机福利观看| 波多野结衣巨乳人妻| 亚洲五月色婷婷综合| 久久人人精品亚洲av| 国产黄a三级三级三级人| 别揉我奶头~嗯~啊~动态视频| 欧美日韩瑟瑟在线播放| 在线免费观看的www视频| 青草久久国产| 色老头精品视频在线观看| 国产一区在线观看成人免费| 亚洲专区中文字幕在线| 精品乱码久久久久久99久播| 亚洲黑人精品在线| 欧美成人一区二区免费高清观看 | 亚洲色图 男人天堂 中文字幕| 好看av亚洲va欧美ⅴa在| 美女大奶头视频| 国产欧美日韩一区二区精品| 欧美激情久久久久久爽电影 | 久久午夜综合久久蜜桃| www.自偷自拍.com| 18禁裸乳无遮挡免费网站照片 | 久久人人爽av亚洲精品天堂| 日韩精品青青久久久久久| 午夜免费激情av| 久久国产亚洲av麻豆专区| 午夜免费成人在线视频| 国产亚洲av嫩草精品影院| 亚洲国产欧美网| 老司机在亚洲福利影院| 免费看a级黄色片| 国产成人精品无人区| 正在播放国产对白刺激| 亚洲国产看品久久| 无人区码免费观看不卡| 女人爽到高潮嗷嗷叫在线视频| 99国产精品免费福利视频| 精品无人区乱码1区二区| 国产免费男女视频| 国产97色在线日韩免费| 999精品在线视频| 日本 欧美在线| 日日夜夜操网爽| 一区二区日韩欧美中文字幕| 国产成人一区二区三区免费视频网站| 在线观看66精品国产| 欧美日韩乱码在线| 国产成人av激情在线播放| 一边摸一边做爽爽视频免费| 亚洲av第一区精品v没综合| 久久午夜综合久久蜜桃| 激情视频va一区二区三区| 国产精品av久久久久免费| √禁漫天堂资源中文www| 午夜福利18| 日日干狠狠操夜夜爽| 国产精品一区二区在线不卡| 久久婷婷人人爽人人干人人爱 | 国产亚洲精品第一综合不卡| 精品欧美国产一区二区三| 国产1区2区3区精品| 国产欧美日韩一区二区三| 亚洲人成电影免费在线| 少妇 在线观看| 狠狠狠狠99中文字幕| 人人澡人人妻人| www.自偷自拍.com| av网站免费在线观看视频| 啦啦啦韩国在线观看视频| 美女午夜性视频免费| 日韩一卡2卡3卡4卡2021年| 国产精品98久久久久久宅男小说| 久久精品人人爽人人爽视色| 看片在线看免费视频| 亚洲片人在线观看| 成人国语在线视频| 国产一区二区三区视频了| 亚洲精品美女久久av网站| 国产成人欧美在线观看| 久久亚洲真实| 国产精品久久视频播放| 午夜精品在线福利| 成人18禁在线播放| 在线观看免费视频日本深夜| 国产一卡二卡三卡精品| 亚洲国产精品成人综合色| 精品欧美一区二区三区在线| 久久精品亚洲精品国产色婷小说| 亚洲成人国产一区在线观看| 美女扒开内裤让男人捅视频| 黑人巨大精品欧美一区二区蜜桃| 电影成人av| 欧美国产精品va在线观看不卡| 深夜精品福利| 很黄的视频免费| 女同久久另类99精品国产91| 亚洲欧美一区二区三区黑人| 国产精品亚洲美女久久久| 18禁国产床啪视频网站| 很黄的视频免费| 18禁观看日本| 亚洲人成电影观看| 国产精品1区2区在线观看.| 琪琪午夜伦伦电影理论片6080| 老司机福利观看| 欧美av亚洲av综合av国产av| 99久久国产精品久久久| 麻豆久久精品国产亚洲av| 18禁美女被吸乳视频| 免费在线观看日本一区| 欧美黑人欧美精品刺激| 午夜久久久在线观看| 亚洲国产欧美日韩在线播放| 欧美人与性动交α欧美精品济南到| 欧美国产精品va在线观看不卡| 免费观看人在逋| 日韩大尺度精品在线看网址 | 国产精品爽爽va在线观看网站 | 欧美中文日本在线观看视频| 99久久99久久久精品蜜桃| 欧美日韩福利视频一区二区| 成人手机av| 亚洲天堂国产精品一区在线| 91成人精品电影| 男人操女人黄网站| 国产一区二区三区视频了| 欧美在线一区亚洲| 精品久久久久久久人妻蜜臀av | 十八禁人妻一区二区| 亚洲专区字幕在线| 久久人妻福利社区极品人妻图片| 91精品国产国语对白视频| 亚洲国产看品久久| 亚洲性夜色夜夜综合| 亚洲欧美日韩高清在线视频| 国产蜜桃级精品一区二区三区| www.精华液| 亚洲精品一区av在线观看| 色播亚洲综合网| 欧美日韩瑟瑟在线播放| 成在线人永久免费视频| 搡老熟女国产l中国老女人| 天堂√8在线中文| 欧美另类亚洲清纯唯美| 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 久久亚洲精品不卡| 自拍欧美九色日韩亚洲蝌蚪91| 免费在线观看影片大全网站| 久久久久九九精品影院| 亚洲三区欧美一区| 久久亚洲真实| 亚洲三区欧美一区| 一级作爱视频免费观看| 成人av一区二区三区在线看| 如日韩欧美国产精品一区二区三区| 亚洲精品国产区一区二| 禁无遮挡网站| 亚洲欧美激情在线| tocl精华| 欧美一区二区精品小视频在线| 精品熟女少妇八av免费久了| 亚洲精品中文字幕在线视频| 人人妻人人澡人人看| 女人被躁到高潮嗷嗷叫费观| 真人一进一出gif抽搐免费| 少妇熟女aⅴ在线视频| 国产成人精品久久二区二区91| 欧美日本视频| 午夜亚洲福利在线播放| 亚洲国产欧美网| 国产成年人精品一区二区| 国内精品久久久久精免费| 亚洲精品久久成人aⅴ小说| 黄色丝袜av网址大全| 99re在线观看精品视频| 日本撒尿小便嘘嘘汇集6| 久久国产乱子伦精品免费另类| 国产欧美日韩综合在线一区二区| 乱人伦中国视频| 精品欧美国产一区二区三| 亚洲一区中文字幕在线| av视频免费观看在线观看| 久久性视频一级片| 一本大道久久a久久精品| 精品人妻在线不人妻| 午夜免费激情av| 一级毛片高清免费大全| 国产蜜桃级精品一区二区三区| av福利片在线| 国产亚洲av高清不卡| 欧美成人性av电影在线观看| 久热爱精品视频在线9| 男女下面插进去视频免费观看| 国语自产精品视频在线第100页| 亚洲欧美一区二区三区黑人| 国产亚洲精品一区二区www| 中文字幕另类日韩欧美亚洲嫩草| 99久久99久久久精品蜜桃| 91av网站免费观看| 亚洲一区二区三区不卡视频| 亚洲第一电影网av| 在线免费观看的www视频| 最新美女视频免费是黄的| 国产男靠女视频免费网站| 国产亚洲欧美在线一区二区| 久久久久久久久免费视频了| 男人舔女人下体高潮全视频| 午夜久久久久精精品| 国产精品香港三级国产av潘金莲| 女人被狂操c到高潮| 国产亚洲精品第一综合不卡| 精品无人区乱码1区二区| 大型黄色视频在线免费观看| 嫁个100分男人电影在线观看| 12—13女人毛片做爰片一| 青草久久国产| 91字幕亚洲| 国产av又大| 午夜激情av网站| 一本久久中文字幕| 激情视频va一区二区三区| 日本一区二区免费在线视频| 一进一出抽搐gif免费好疼| 一边摸一边抽搐一进一小说| 欧美一区二区精品小视频在线| 真人一进一出gif抽搐免费| 亚洲av成人不卡在线观看播放网| 亚洲人成伊人成综合网2020| 国产真人三级小视频在线观看| 成人国语在线视频| 丰满的人妻完整版| 亚洲欧洲精品一区二区精品久久久| 丁香欧美五月| 国产亚洲欧美精品永久| 日韩欧美三级三区| 黄片播放在线免费| 成人亚洲精品一区在线观看| 搡老熟女国产l中国老女人| cao死你这个sao货| 九色亚洲精品在线播放| 一本久久中文字幕| 色哟哟哟哟哟哟| 亚洲五月天丁香| 久久久久久免费高清国产稀缺| 久久性视频一级片| 婷婷精品国产亚洲av在线| 亚洲欧美激情综合另类| 一区福利在线观看| 日本在线视频免费播放| 变态另类丝袜制服| 亚洲精品av麻豆狂野| 搡老妇女老女人老熟妇| 国产激情欧美一区二区| 色婷婷久久久亚洲欧美| 亚洲成人免费电影在线观看| 亚洲人成77777在线视频| 最近最新中文字幕大全电影3 | 免费看美女性在线毛片视频| 国产日韩一区二区三区精品不卡| 此物有八面人人有两片| 国产av精品麻豆| 级片在线观看| 亚洲伊人色综图| 精品免费久久久久久久清纯| 在线播放国产精品三级| 91成人精品电影| 午夜精品久久久久久毛片777| 久久久国产欧美日韩av| 国产单亲对白刺激| 亚洲第一青青草原| 日日摸夜夜添夜夜添小说| 两个人视频免费观看高清| 欧美精品啪啪一区二区三区| 久久人妻熟女aⅴ| 欧美国产精品va在线观看不卡| 欧美日韩亚洲综合一区二区三区_| 国产精品免费视频内射| 99riav亚洲国产免费| 一本大道久久a久久精品| 51午夜福利影视在线观看| 激情在线观看视频在线高清| 国产精品 国内视频| 国产亚洲精品综合一区在线观看 | 午夜影院日韩av| 日本一区二区免费在线视频| 精品一品国产午夜福利视频| 首页视频小说图片口味搜索| 亚洲国产高清在线一区二区三 | 午夜福利,免费看| 无人区码免费观看不卡| 亚洲av五月六月丁香网| 欧美日本中文国产一区发布| 亚洲在线自拍视频| 亚洲国产看品久久| 欧美色欧美亚洲另类二区 | 亚洲欧美一区二区三区黑人| 国产精品99久久99久久久不卡| 国产精品美女特级片免费视频播放器 | 法律面前人人平等表现在哪些方面| 午夜福利一区二区在线看| 成人亚洲精品av一区二区| 亚洲熟妇熟女久久| 黑人操中国人逼视频| 久久精品成人免费网站| 精品一区二区三区视频在线观看免费| 制服诱惑二区| 亚洲全国av大片| 欧美黄色淫秽网站| 亚洲第一av免费看| 在线播放国产精品三级| 99精品欧美一区二区三区四区| 亚洲色图综合在线观看| 亚洲av美国av| 久久久久久久久免费视频了| 国产亚洲精品久久久久5区| 亚洲va日本ⅴa欧美va伊人久久| 真人一进一出gif抽搐免费| 久久久久久人人人人人| 午夜福利成人在线免费观看| 亚洲视频免费观看视频| 91在线观看av| 亚洲精品中文字幕一二三四区| 国产成人欧美| 亚洲欧美激情综合另类| 欧美激情 高清一区二区三区| 成人欧美大片| 电影成人av| 久久久久久免费高清国产稀缺| 可以免费在线观看a视频的电影网站| 日韩欧美一区视频在线观看| 久久久久久久久中文| 狠狠狠狠99中文字幕| 少妇熟女aⅴ在线视频| 亚洲人成电影免费在线| 黑人巨大精品欧美一区二区蜜桃| 亚洲天堂国产精品一区在线| 亚洲精品美女久久久久99蜜臀| 又黄又爽又免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人不卡在线观看播放网| 99久久综合精品五月天人人| 日韩三级视频一区二区三区| 成人永久免费在线观看视频| 99在线视频只有这里精品首页| 嫩草影视91久久| 丰满人妻熟妇乱又伦精品不卡| 长腿黑丝高跟| 国产精品98久久久久久宅男小说| 好男人在线观看高清免费视频 | 在线免费观看的www视频| 日日摸夜夜添夜夜添小说| 18美女黄网站色大片免费观看| 国产一区二区三区视频了| 欧美另类亚洲清纯唯美| 88av欧美| 欧美日韩瑟瑟在线播放| 禁无遮挡网站| 久久狼人影院| 淫妇啪啪啪对白视频| 色播亚洲综合网| 日本黄色视频三级网站网址| 成人三级做爰电影| 99热只有精品国产| 国产精品国产高清国产av| 十八禁人妻一区二区| 免费少妇av软件| 久久国产精品人妻蜜桃| 国产成人精品无人区| 一边摸一边做爽爽视频免费| 亚洲第一欧美日韩一区二区三区| 国产av一区二区精品久久| 国产成+人综合+亚洲专区| 成年人黄色毛片网站| 国产成人系列免费观看| 99国产精品99久久久久| 亚洲欧美一区二区三区黑人| 嫩草影视91久久| 国产激情久久老熟女| 亚洲人成77777在线视频| 精品福利观看| 在线观看免费午夜福利视频| 亚洲欧美日韩无卡精品| xxx96com| 久久久精品欧美日韩精品|