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

    基于試飛數(shù)據(jù)的民用發(fā)動機(jī)部件特性修正方法

    2024-02-05 21:44:06周超李密
    科技創(chuàng)新與應(yīng)用 2024年5期
    關(guān)鍵詞:航空發(fā)動機(jī)

    周超 李密

    摘? 要:由于民用航空發(fā)動機(jī)外涵尾噴管容易不臨界,不同馬赫數(shù)下發(fā)動機(jī)工作線不同,僅利用地面臺架試驗(yàn)數(shù)據(jù)修正發(fā)動機(jī)部件特性,整機(jī)模型精度較差。因此,該文利用發(fā)動機(jī)試飛數(shù)據(jù),提出一種改進(jìn)的部件特性刪除法,將發(fā)動機(jī)非設(shè)計(jì)點(diǎn)計(jì)算轉(zhuǎn)化為設(shè)計(jì)點(diǎn)計(jì)算,減少目標(biāo)方程數(shù)量,簡化算法,同時(shí)考慮對不同轉(zhuǎn)速和不同馬赫數(shù)下部件特性修正,以此建立整個(gè)飛行包線內(nèi)發(fā)動機(jī)整機(jī)數(shù)學(xué)模型,計(jì)算結(jié)果與試飛數(shù)據(jù)進(jìn)行對比,發(fā)動機(jī)整機(jī)模型在整個(gè)飛行包線內(nèi)各參數(shù)誤差均小于3%,發(fā)動機(jī)推力誤差更是小于1%。

    關(guān)鍵詞:航空發(fā)動機(jī);發(fā)動機(jī)模型;部件特性;修正因子;部件特性刪除法

    中圖分類號:V235.13? ? ? 文獻(xiàn)標(biāo)志碼:A? ? ? ? ? 文章編號:2095-2945(2024)05-0052-05

    Abstract: Due to the fact that the nozzle is easy to be non-critical, and the working line of the engine is different under different Mach numbers, the accuracy of the whole engine model is poor only by using the ground test data to modify the characteristics of the engine components. So this paper based on the data of engine flight test, an improved method of deleting the characteristics of components is proposed, which transforms the engine non-design point calculation into the design point calculation. The number of target equations is reduced and the algorithm is simplified. As well the modification of component characteristics at different speeds and Mach numbers is considered. Based on this the mathematical model of the whole engine in the whole flight envelope is established. Compared with the flight test data ,the error of all parameters in the whole flight envelope is less than 3%, and thrust error is less than 1%.

    Keywords: Aero Engine; engine model; component characteristic; correction factor; deletion of part characteristics

    民用航空發(fā)動機(jī)在適航審定過程中,《航空發(fā)動機(jī)適航規(guī)定》第33.5(b)(2)條中明確規(guī)定必須提供發(fā)動機(jī)推力在非標(biāo)準(zhǔn)大氣條件下的修正程序,基于部件特性的航空發(fā)動機(jī)的總體性能計(jì)算中,發(fā)動機(jī)部件特性圖的準(zhǔn)確程度直接決定了發(fā)動機(jī)特性的計(jì)算精度。一方面在發(fā)動機(jī)試飛過程中,往往得不到真實(shí)的部件特性,使用通用部件特性發(fā)動機(jī)的特性計(jì)算精度得不到保證;另一方面已經(jīng)定型量產(chǎn)的發(fā)動機(jī),制造和安裝存在一定的偏差,長時(shí)間工作性能的衰退,真實(shí)的部件特性會產(chǎn)生偏離。因此,為了提高發(fā)動機(jī)模型的精度,必須對發(fā)動機(jī)最初的部件特性或者通用部件特性進(jìn)行修正。

    文獻(xiàn)[1-2]提出了常用的部件特性修正方法,通常使用設(shè)計(jì)點(diǎn)的修正因子,修正發(fā)動機(jī)部件完整特性,當(dāng)偏離設(shè)計(jì)點(diǎn)較遠(yuǎn)時(shí),模型的精度不能得到有效保證。文獻(xiàn)[3-4]以發(fā)動機(jī)測量參數(shù)和性能參數(shù)為目標(biāo)函數(shù),選取修正因子為變量,利用單純形法獲得發(fā)動機(jī)真實(shí)部件特性。文獻(xiàn)[5]利用“部件特性刪除法”計(jì)算發(fā)動機(jī)不同轉(zhuǎn)速下的部件特性的修正因子。文獻(xiàn)[6-10]基于試驗(yàn)數(shù)據(jù)的優(yōu)化算法:遺傳算法、粒子群法和最小二乘法,獲得發(fā)動機(jī)部件特性。

    本文基于文獻(xiàn)[5]采用了改進(jìn)的部件特性刪除法對通用部件特性進(jìn)行修正,其中并未介紹當(dāng)尾噴管面積固定時(shí)特性的修正方法,且僅獲得不同轉(zhuǎn)速下的修正因子,然后使用同一修正因子對該轉(zhuǎn)速下的特性進(jìn)行修正,而對于民用大涵道比分排渦扇發(fā)動機(jī)而言,內(nèi)外涵尾噴管常處于不臨界狀態(tài),不同馬赫數(shù)下,同樣轉(zhuǎn)速發(fā)動機(jī)工作點(diǎn)不同,使用單點(diǎn)因子修正同一轉(zhuǎn)速特性勢必會影響發(fā)動機(jī)性能模型的精度。因此,本文在文獻(xiàn)[5]方法的基礎(chǔ)上,加以改進(jìn)利用設(shè)計(jì)點(diǎn)和非設(shè)計(jì)點(diǎn),綜合計(jì)算方法對部件特性進(jìn)行修正,獲得整個(gè)包線范圍內(nèi)發(fā)動機(jī)部件特性的修正因子,以此獲得發(fā)動機(jī)真實(shí)部件特性,并建立發(fā)動機(jī)整機(jī)計(jì)算模型并與試飛數(shù)據(jù)對比,結(jié)果表明,該方法修正的部件特性使得發(fā)動機(jī)整機(jī)模型在整個(gè)包線范圍內(nèi)具有較高的精度。

    1? 整機(jī)數(shù)學(xué)模型建立

    采用變比熱法建立基于部件特性的發(fā)動機(jī)數(shù)學(xué)模型,具有較高的精度,部件級的氣動熱力學(xué)模型參考文獻(xiàn)[11-12]。發(fā)動機(jī)設(shè)計(jì)點(diǎn)計(jì)算的氣動熱力學(xué)模型從進(jìn)氣道至尾噴管各個(gè)部件依次計(jì)算。發(fā)動機(jī)非設(shè)計(jì)點(diǎn)計(jì)算與設(shè)計(jì)點(diǎn)計(jì)算不同,部件級氣動熱力學(xué)模型需要根據(jù)部件特性和修正因子插值求出部件工作點(diǎn),然后依照部件順序依次計(jì)算。發(fā)動機(jī)非設(shè)計(jì)點(diǎn)計(jì)算已知量較少,計(jì)算過程中遇到未知量給定一個(gè)初始值,稱為獨(dú)立變量,依據(jù)共同工作條件建立相應(yīng)的平衡方程,將工作點(diǎn)計(jì)算過程轉(zhuǎn)化為非線性方程的求解過程。

    在發(fā)動機(jī)非設(shè)計(jì)點(diǎn)求解過程中,利用平衡方程殘差的計(jì)算結(jié)果,反復(fù)修正獨(dú)立變量的數(shù)值,使得方程的殘差越來越小,直到滿足一定的精度要求,則迭代成功,完成民用大涵道比分排渦扇發(fā)動機(jī)非設(shè)計(jì)點(diǎn)計(jì)算,獲得發(fā)動機(jī)各個(gè)部件的參數(shù)。

    2? 部件特性修正方法

    2.1? 通用部件特性

    大涵道比分排渦扇發(fā)動機(jī)非設(shè)計(jì)點(diǎn)計(jì)算過程中,是以渦輪發(fā)動機(jī)部件特性曲線為基礎(chǔ)進(jìn)行插值計(jì)算。通用部件特性風(fēng)扇、壓氣機(jī)和渦輪的部件特性表達(dá)方式如式(1)和式(2)所示,選取的通用部件特性圖如圖1和圖2所示。

    換算流量:Wcor=F(ncor,π)? ,? (1)

    效率:η=F(ncor,π)? ,? ?(2)

    式中:ncor表示相對換算轉(zhuǎn)速,π表示壓比。

    2.2? 部件特性修正因子

    建立高精度的發(fā)動機(jī)模型的前提是必須獲得實(shí)際發(fā)動機(jī)的部件特性,但實(shí)際部件特性受到裝機(jī)條件、使用時(shí)長等因素的影響。因此,本文提出了利用試飛數(shù)據(jù),通過通用部件特性和計(jì)算的修正因子獲得發(fā)動機(jī)裝機(jī)條件下的實(shí)際部件特性

    Xact=Fm×Xref , (3)

    式中:Xact為發(fā)動機(jī)實(shí)際部件特性,Xref為發(fā)動機(jī)通用部件特性,F(xiàn)m為修正因子。

    如前文所述,基于部件特性建立的發(fā)動機(jī)仿真模型,通常使用設(shè)計(jì)點(diǎn)的修正因子,修正整個(gè)部件特性,該方法僅能保證設(shè)計(jì)點(diǎn)附近的精度。對于民用大涵道比分排渦扇發(fā)動機(jī),發(fā)動機(jī)工作線隨著進(jìn)口馬赫數(shù)的變化而發(fā)生偏離,文獻(xiàn)[5]雖然獲得了不同轉(zhuǎn)速下的部件特性修正因子,依然不能保證發(fā)動機(jī)模型的精度。因此,需要基于試飛數(shù)據(jù)在發(fā)動機(jī)整個(gè)飛行包線內(nèi),不同馬赫數(shù)、不同轉(zhuǎn)速下所有的工作點(diǎn)計(jì)算真實(shí)的部件特性修正因子,才能獲得具有高精度的發(fā)動機(jī)仿真數(shù)學(xué)模型。

    2.3? 改進(jìn)的部件特性刪除法

    傳統(tǒng)的部件特性修正方法是一個(gè)將測試參數(shù)方程和平衡方程組成優(yōu)化目標(biāo)函數(shù),并將其最小化的過程

    ΔY=Ym-Y′m=Ym-F(Xn) ,(4)

    式中:m代表方程個(gè)數(shù)(目標(biāo)參數(shù)和發(fā)動機(jī)模型的平衡方程),n為變量(修正因子和發(fā)動機(jī)模型獨(dú)立變量),修正因子的求解過程即是求解非線性方程組的過程,同時(shí)方程具有解的條件為m≥n。

    對于大涵道比分排渦扇發(fā)動機(jī)而言,部件特性修正風(fēng)扇、增壓級、壓氣機(jī)、低壓渦輪的換算流量和效率共8個(gè)變量。加上發(fā)動機(jī)自身的7個(gè)平衡方程,共有15個(gè)目標(biāo)方程。因此在使用傳統(tǒng)的部件特性修正方法時(shí),如果m=n,目標(biāo)方程過多,方程收斂性很差。當(dāng)m>n時(shí),部件特性的修正精度無法得到滿足。

    因此本文基于部件特性刪除法,即將各部件進(jìn)出口試飛測量參數(shù)代替發(fā)動機(jī)部件特性,直接參與到發(fā)動機(jī)氣動熱力學(xué)計(jì)算,以發(fā)動機(jī)共同工作條件為約束,建立共同工作方程,減少了方程數(shù)量,增加模型的收斂性。同時(shí)以試飛直接測量參數(shù)為輸入量,模型計(jì)算結(jié)果代表了發(fā)動機(jī)實(shí)際工作點(diǎn),計(jì)算得到的部件的流量、壓比、效率代表了發(fā)動機(jī)部件的實(shí)際工作特性,繼而針對各部件特性圖獲得工作點(diǎn)的修正因子。該方法建立的修正模型滿足發(fā)動機(jī)各部件共同工作條件,利用修正因子結(jié)合原始特性還原發(fā)動機(jī)真實(shí)部件特性,并替代原始特性進(jìn)行發(fā)動機(jī)整機(jī)性能計(jì)算,計(jì)算結(jié)果與試飛直接測量參數(shù)一一對應(yīng),具有非常高的精度。

    考慮到大涵道比分排渦扇發(fā)動機(jī)內(nèi)外涵尾噴管在中等轉(zhuǎn)速和低轉(zhuǎn)速下不臨界,使得低壓轉(zhuǎn)子的工作線發(fā)生偏移,而高壓轉(zhuǎn)子工作線幾乎不變,因此,需要針對不同馬赫數(shù)建立不同修正模型。具體工作思路為,先修正設(shè)計(jì)點(diǎn)馬赫數(shù)下高低壓轉(zhuǎn)子不同轉(zhuǎn)速下的部件特性,然后在修正非設(shè)計(jì)點(diǎn)馬赫數(shù)部件特性時(shí),將獲得的設(shè)計(jì)點(diǎn)馬赫數(shù)高壓轉(zhuǎn)子特性作為已知量帶入模型中,僅修正低壓轉(zhuǎn)子特性,減少了修正特性的工作量,并能保證模型的精度。

    3? 設(shè)計(jì)點(diǎn)馬赫數(shù)下特性圖修正

    基于上述理論分析,在不同轉(zhuǎn)速下部件特性修正過程中壓氣機(jī)、燃燒室、渦輪轉(zhuǎn)化為設(shè)計(jì)點(diǎn)計(jì)算,但是內(nèi)外涵尾噴管面積固定,需要建立內(nèi)外涵尾噴管流量平衡方程;高壓渦輪處于臨界狀態(tài),渦輪的換算流量和效率幾乎不變,可以將設(shè)計(jì)點(diǎn)修正的高壓渦輪特性作為已知條件帶入模型中,因此需要建立高壓渦輪的流量平衡方程;利用測量得到的發(fā)動機(jī)燃油流量,建立相應(yīng)的目標(biāo)方程。因此部件特性修正過程中平衡方程如下。

    利用牛頓法,獲得不同轉(zhuǎn)速下部件的效率和換算流量,通過該狀態(tài)點(diǎn)在原有特性圖上的位置,計(jì)算不同轉(zhuǎn)速下部件特性修正因子,繼而獲得設(shè)計(jì)點(diǎn)馬赫數(shù)下部件的真實(shí)特性。如圖3—圖6為修正前后風(fēng)扇、增壓級、壓氣機(jī)和低壓渦輪特性圖。

    4? 不同馬赫數(shù)下特性圖修正

    與軍用小涵道比渦扇發(fā)動機(jī)不同的是,大涵道比分排渦扇發(fā)動機(jī)內(nèi)外涵尾噴管很容易不臨界,這就造成了發(fā)動機(jī)的低壓轉(zhuǎn)子工作線隨著馬赫數(shù)的變化而不再是同一條直線,如圖7—圖9所示。

    圖7—圖9不同馬赫數(shù)下風(fēng)扇和增壓級工作線的偏移量較大,由于渦輪臨界馬赫數(shù)對高壓轉(zhuǎn)子的工作線幾乎無影響,因此,需要利用試飛數(shù)據(jù)對低壓轉(zhuǎn)子部件特性在不同馬赫數(shù)下進(jìn)行修正。平衡方程與第2章節(jié)相似,壓氣機(jī)特性需要作為已知量帶入模型中,則需要增加壓氣機(jī)和增壓級流量平衡方程,同時(shí)增加壓氣機(jī)壓比作為獨(dú)立變量,進(jìn)行求解,獲得低壓部件特性修正因子,計(jì)算得到不同馬赫數(shù)下低壓轉(zhuǎn)子的真實(shí)部件特性。

    如圖10和圖11所示,同一轉(zhuǎn)速下部件特性分為了3段,中間淺灰色圍成的區(qū)域是利用不同馬赫數(shù)下的試飛數(shù)據(jù)修正的部件特性,其他兩段由于發(fā)動機(jī)工作點(diǎn)未在此處,未進(jìn)行修正,可以看出修正后的特性與未修正的特性走勢完全不同,表征了本文完成了發(fā)動機(jī)整個(gè)包線內(nèi)不同轉(zhuǎn)速下各部件特性的修正。

    5? 方法驗(yàn)證

    用修正后的特性圖將原來的通用特性圖替換掉,在不同高度和不同馬赫數(shù)下進(jìn)行發(fā)動機(jī)整機(jī)模型計(jì)算,并與試飛數(shù)據(jù)對比,如圖12和表1所示。

    由圖12和表1表明,在工作包線范圍內(nèi)模型計(jì)算結(jié)果試飛數(shù)據(jù)的誤差很小,低壓渦輪出口總溫和耗油量誤差相對較大,但都在3%的誤差范圍之內(nèi),完全滿足工程需求,其他各截面和性能參數(shù)誤差均在2%以內(nèi),尤其是計(jì)算的推力誤差更是在1%以內(nèi)。由此表明本文使用的改進(jìn)的“部件特性刪除法”完全有效,并且具有較高的精度。

    6? 結(jié)論

    為充分利用發(fā)動機(jī)試飛數(shù)據(jù),提高部件特性修正的精度,本文提出了一種改進(jìn)的“部件特性刪除法”,完成部件特性的修正,研究結(jié)果表明如下。

    1)本文使用的部件特性刪除法,將已知的試飛數(shù)據(jù)的截面參數(shù)直接參與到發(fā)動機(jī)氣動熱力學(xué)計(jì)算中,大大減少了目標(biāo)方程的數(shù)量,解決了特性計(jì)算的不收斂的問題,并且保證了模型的精度。

    2)民用大涵道比分排渦扇發(fā)動機(jī)內(nèi)外涵尾噴管不臨界,低壓轉(zhuǎn)子工作線偏移,說明僅在不同轉(zhuǎn)速下修正部件特性是不夠的,本文利用試飛數(shù)據(jù),在不同馬赫數(shù)不同發(fā)動機(jī)轉(zhuǎn)速下對部件特性進(jìn)行修正,結(jié)果表明,能夠完全修正出發(fā)動機(jī)整個(gè)包線范圍內(nèi)工作點(diǎn)的實(shí)際部件特性。

    3)本文使用的部件特性刪除法修正的部件特性,帶入到整機(jī)模型進(jìn)行計(jì)算,計(jì)算結(jié)果與實(shí)際試飛數(shù)據(jù)具有較高的貼合度,誤差較小,滿足工程需求,解決了基于試飛數(shù)據(jù)未知實(shí)際部件特性的情況下建立高精度發(fā)動機(jī)數(shù)學(xué)模型最棘手的問題。

    參考文獻(xiàn):

    [1] STAMATIS A, MATHIOUDAKIS K. Adaptive simulation of gas turbine performance[J].ASME J. Eng. Gas Turbines Power,1990,112(2):168-175.

    [2] KONG C, KHO S, KI J. Component map generation of a gas turbine using genetic algorithms[R].ASME Paper GT, 2004:2004-53736.

    [3] 吳虎,肖洪,蔣建軍.渦扇發(fā)動機(jī)部件特性自適應(yīng)模擬[J].推進(jìn)技術(shù),2005(5):430-433.

    [4] 蔣建軍,魯偉,尹洪舉,等.航空發(fā)動機(jī)部件特性自適應(yīng)計(jì)算方法[J].航空計(jì)算技術(shù),2008(4):11-13,17.

    [5] 陳玉春,黃興,徐思遠(yuǎn),等.渦輪發(fā)動機(jī)部件特性自適應(yīng)模型的確定方法[J].推進(jìn)技術(shù),2008(2):214-218.

    [6] LI Y G, MARINAI L, LO G E, et al. Multiple point adaptive performance simulation tuned to aerospace test-bed data[J]. Journal of Propulsion and Power,2009,25(3):635-641.

    [7] 劉永葆,賀星,黃樹紅.基于改進(jìn)遺傳算法的燃?xì)廨啓C(jī)自適應(yīng)建模[J].航空動力學(xué)報(bào),2012,27(3):695-700.

    [8] 劉盾,江和甫,都慶福.部件特性圖的非線性最小二乘擬合修正[J].燃?xì)鉁u輪試驗(yàn)與研究,2000(4):18-21.

    [9] 楊欣毅,沈偉,劉海峰,等.一種應(yīng)用滑動最小二乘求取壓氣機(jī)特性的方法[J].航空動力學(xué)報(bào),2009,24(8):1741-1746.

    [10] 楊欣毅,王永華,賀孝濤,等.基于改進(jìn)量子粒子群的發(fā)動機(jī)部件特性修正[J].航空計(jì)算技術(shù),2012,42(2):4-8.

    [11] 廉小純,吳虎.航空燃?xì)廨啓C(jī)原理[M].北京:國防工業(yè)出版社,2001:45-74.

    [12] KONG C,KI J. A new scaling method for component maps of gasturbine using system identification[J]. Journal of Engineering for GasTurbines and Power,2003,125(1):958-979.

    猜你喜歡
    航空發(fā)動機(jī)
    關(guān)于《航空發(fā)動機(jī)結(jié)構(gòu)與系統(tǒng)》課程改革的探索
    基于小波分析的航空發(fā)動機(jī)氣動失穩(wěn)信號分析方法
    整體葉盤葉片加工技術(shù)研究
    某型發(fā)動機(jī)喘振故障的研究
    “航空發(fā)動機(jī)強(qiáng)度與振動”教學(xué)改革的探索與思考
    基于PDM的航空發(fā)動機(jī)技術(shù)狀態(tài)管理
    航空發(fā)動機(jī)振動監(jiān)測硬件系統(tǒng)振動通道的設(shè)計(jì)
    航空發(fā)動機(jī)小零件的標(biāo)準(zhǔn)化實(shí)踐與數(shù)字化管理
    科技視界(2016年6期)2016-07-12 09:59:21
    航空發(fā)動機(jī)管理信息系統(tǒng)研究
    科技視界(2016年14期)2016-06-08 12:05:58
    淺談航空發(fā)動機(jī)起動與點(diǎn)火系統(tǒng)
    科技視界(2016年14期)2016-06-08 00:41:57
    热99re8久久精品国产| 精品福利观看| 一区二区三区激情视频| 国产一区在线观看成人免费| 亚洲av美国av| 夜夜看夜夜爽夜夜摸| 亚洲精品一区av在线观看| 国产又色又爽无遮挡免费看| 亚洲人成77777在线视频| 欧美黑人巨大hd| 亚洲成av人片免费观看| 一本久久中文字幕| 正在播放国产对白刺激| 欧美绝顶高潮抽搐喷水| 国产一区二区三区视频了| 亚洲美女黄片视频| 国产精品二区激情视频| 久久精品91无色码中文字幕| 十分钟在线观看高清视频www| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美激情综合另类| 人人澡人人妻人| 国产伦一二天堂av在线观看| 午夜福利高清视频| 日韩欧美三级三区| 夜夜爽天天搞| 亚洲人成伊人成综合网2020| 中文字幕久久专区| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 日韩欧美三级三区| 人成视频在线观看免费观看| 久久久国产成人精品二区| 久久久久久久久免费视频了| 国产亚洲精品一区二区www| 女同久久另类99精品国产91| 搡老岳熟女国产| 国产欧美日韩精品亚洲av| 亚洲欧美一区二区三区黑人| 亚洲欧美激情综合另类| 中文字幕高清在线视频| 亚洲精品国产区一区二| 久久精品夜夜夜夜夜久久蜜豆 | 999精品在线视频| 亚洲美女黄片视频| 女性生殖器流出的白浆| 一区二区三区精品91| 在线av久久热| 精品国产国语对白av| 成人特级黄色片久久久久久久| 亚洲国产毛片av蜜桃av| 99国产综合亚洲精品| 午夜精品在线福利| 法律面前人人平等表现在哪些方面| 亚洲av中文字字幕乱码综合 | 变态另类成人亚洲欧美熟女| 无遮挡黄片免费观看| 波多野结衣高清作品| 99国产精品一区二区蜜桃av| 很黄的视频免费| 一本精品99久久精品77| 777久久人妻少妇嫩草av网站| 国产一级毛片七仙女欲春2 | 国产区一区二久久| 中文字幕av电影在线播放| 亚洲色图av天堂| 国产亚洲精品第一综合不卡| 十八禁网站免费在线| 日本a在线网址| 99国产精品一区二区蜜桃av| 亚洲五月天丁香| 免费av毛片视频| 成人欧美大片| 亚洲电影在线观看av| 无遮挡黄片免费观看| 午夜影院日韩av| 日韩欧美国产一区二区入口| 天堂动漫精品| 一边摸一边抽搐一进一小说| 观看免费一级毛片| 欧美日韩精品网址| 亚洲国产看品久久| 国产精品av久久久久免费| 久久久久国内视频| 黄色片一级片一级黄色片| 久久午夜综合久久蜜桃| 精品不卡国产一区二区三区| 国产精品影院久久| 免费电影在线观看免费观看| 亚洲精华国产精华精| 久久精品国产99精品国产亚洲性色| 国产高清videossex| 日韩国内少妇激情av| 老熟妇仑乱视频hdxx| 国产av一区二区精品久久| 亚洲狠狠婷婷综合久久图片| 精品一区二区三区四区五区乱码| 男女那种视频在线观看| 免费一级毛片在线播放高清视频| 身体一侧抽搐| 免费高清视频大片| 午夜福利视频1000在线观看| 国内久久婷婷六月综合欲色啪| 国产免费男女视频| 久久久国产欧美日韩av| 两个人视频免费观看高清| 热99re8久久精品国产| 一卡2卡三卡四卡精品乱码亚洲| 丁香六月欧美| 午夜视频精品福利| 国产精品日韩av在线免费观看| 好男人在线观看高清免费视频 | 成人亚洲精品一区在线观看| 99久久精品国产亚洲精品| 国产成人啪精品午夜网站| 国产精品 欧美亚洲| 女警被强在线播放| 男人操女人黄网站| 久久香蕉激情| 99在线视频只有这里精品首页| 国产熟女午夜一区二区三区| aaaaa片日本免费| 午夜激情福利司机影院| www国产在线视频色| 91九色精品人成在线观看| 欧美av亚洲av综合av国产av| 免费看十八禁软件| 怎么达到女性高潮| 久久国产精品男人的天堂亚洲| 母亲3免费完整高清在线观看| 伦理电影免费视频| 久久久久久人人人人人| 欧美激情极品国产一区二区三区| 在线观看日韩欧美| 我的亚洲天堂| 久久香蕉精品热| 亚洲国产毛片av蜜桃av| 精品国产乱码久久久久久男人| 别揉我奶头~嗯~啊~动态视频| 九色国产91popny在线| 国产精品1区2区在线观看.| 国产又爽黄色视频| 亚洲国产精品成人综合色| 在线看三级毛片| 欧美最黄视频在线播放免费| 国产午夜福利久久久久久| 黄色女人牲交| 久久久久九九精品影院| 欧美久久黑人一区二区| 精品国产超薄肉色丝袜足j| 精品乱码久久久久久99久播| 精品欧美国产一区二区三| 亚洲国产精品久久男人天堂| 在线十欧美十亚洲十日本专区| 无遮挡黄片免费观看| 国产一区二区三区在线臀色熟女| 中文字幕精品亚洲无线码一区 | 亚洲国产精品合色在线| 身体一侧抽搐| 亚洲熟妇熟女久久| 国产真实乱freesex| 欧美一级a爱片免费观看看 | 亚洲精品在线观看二区| 久久天躁狠狠躁夜夜2o2o| 黑人欧美特级aaaaaa片| 午夜福利18| 夜夜爽天天搞| 波多野结衣巨乳人妻| 97人妻精品一区二区三区麻豆 | 9191精品国产免费久久| 欧美午夜高清在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产三级在线视频| 精品不卡国产一区二区三区| 国内少妇人妻偷人精品xxx网站 | 亚洲国产精品合色在线| 又紧又爽又黄一区二区| 亚洲av电影不卡..在线观看| 伊人久久大香线蕉亚洲五| 美女午夜性视频免费| 2021天堂中文幕一二区在线观 | 亚洲精品国产精品久久久不卡| 国产精品久久久av美女十八| 午夜免费观看网址| 欧美激情高清一区二区三区| 久久精品91无色码中文字幕| 久久香蕉激情| 欧美国产日韩亚洲一区| 国产又黄又爽又无遮挡在线| www日本在线高清视频| 午夜影院日韩av| 一区二区日韩欧美中文字幕| 精品国产亚洲在线| 久久天堂一区二区三区四区| 性色av乱码一区二区三区2| 老鸭窝网址在线观看| 久久久国产成人精品二区| 91在线观看av| 久久久久亚洲av毛片大全| 视频区欧美日本亚洲| 亚洲欧美激情综合另类| 欧美乱码精品一区二区三区| av天堂在线播放| 老汉色av国产亚洲站长工具| 男人操女人黄网站| 一区二区三区国产精品乱码| 欧美不卡视频在线免费观看 | 999久久久国产精品视频| 黄色视频不卡| 欧美日韩福利视频一区二区| 国产亚洲精品综合一区在线观看 | 亚洲欧洲精品一区二区精品久久久| 波多野结衣巨乳人妻| 亚洲精品一区av在线观看| 麻豆av在线久日| 亚洲av片天天在线观看| 在线永久观看黄色视频| 欧美黄色淫秽网站| 免费在线观看影片大全网站| 久久亚洲真实| 91九色精品人成在线观看| 国产欧美日韩一区二区三| 18禁黄网站禁片午夜丰满| 白带黄色成豆腐渣| 国产精品av久久久久免费| 伦理电影免费视频| 国产男靠女视频免费网站| 丝袜人妻中文字幕| 在线观看舔阴道视频| 91麻豆精品激情在线观看国产| 久久这里只有精品19| 国产精品99久久99久久久不卡| 啦啦啦 在线观看视频| 欧美一级毛片孕妇| 老司机在亚洲福利影院| 日韩精品免费视频一区二区三区| 日韩av在线大香蕉| 午夜成年电影在线免费观看| 一边摸一边做爽爽视频免费| 免费观看人在逋| 人人澡人人妻人| 亚洲美女黄片视频| 国产91精品成人一区二区三区| 黄色成人免费大全| 免费在线观看影片大全网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲精品久久久久5区| 日本撒尿小便嘘嘘汇集6| netflix在线观看网站| 日韩精品青青久久久久久| 亚洲熟妇中文字幕五十中出| 最近最新免费中文字幕在线| 亚洲国产精品久久男人天堂| 欧美日本亚洲视频在线播放| 禁无遮挡网站| 99国产精品一区二区三区| 19禁男女啪啪无遮挡网站| 亚洲国产中文字幕在线视频| 嫩草影视91久久| 啪啪无遮挡十八禁网站| 精品久久久久久久久久久久久 | 一本大道久久a久久精品| ponron亚洲| 亚洲美女黄片视频| 十分钟在线观看高清视频www| 日日爽夜夜爽网站| 国产又爽黄色视频| 亚洲欧美激情综合另类| 老司机靠b影院| 男女视频在线观看网站免费 | 老司机深夜福利视频在线观看| 99热只有精品国产| 美女扒开内裤让男人捅视频| 欧美又色又爽又黄视频| 久久久国产成人免费| 曰老女人黄片| 一夜夜www| 亚洲av中文字字幕乱码综合 | 国产伦一二天堂av在线观看| 久久久久久大精品| 欧美最黄视频在线播放免费| 波多野结衣高清作品| 午夜久久久在线观看| 好男人电影高清在线观看| 嫁个100分男人电影在线观看| 欧美性长视频在线观看| 久久久国产精品麻豆| 欧美久久黑人一区二区| 草草在线视频免费看| 桃色一区二区三区在线观看| 亚洲av成人不卡在线观看播放网| 女警被强在线播放| 精品国产亚洲在线| 黄片大片在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 少妇裸体淫交视频免费看高清 | 亚洲av五月六月丁香网| 欧美日韩亚洲综合一区二区三区_| 啦啦啦观看免费观看视频高清| 午夜福利一区二区在线看| 国产成人欧美在线观看| 午夜影院日韩av| 国产国语露脸激情在线看| 欧美三级亚洲精品| 亚洲五月婷婷丁香| 亚洲中文字幕日韩| 这个男人来自地球电影免费观看| 又紧又爽又黄一区二区| 又大又爽又粗| 久久国产精品影院| 中国美女看黄片| av电影中文网址| 国产精品九九99| 精品第一国产精品| 欧美激情 高清一区二区三区| 欧美另类亚洲清纯唯美| 一a级毛片在线观看| 很黄的视频免费| 久久久久久九九精品二区国产 | 脱女人内裤的视频| 国产欧美日韩精品亚洲av| 国产aⅴ精品一区二区三区波| 久久亚洲精品不卡| 99热这里只有精品一区 | 欧美精品啪啪一区二区三区| 久久人妻福利社区极品人妻图片| 9191精品国产免费久久| 精品久久久久久久毛片微露脸| 啦啦啦韩国在线观看视频| 18禁国产床啪视频网站| 91大片在线观看| av欧美777| 精品欧美一区二区三区在线| 99re在线观看精品视频| 人妻久久中文字幕网| 女人高潮潮喷娇喘18禁视频| 夜夜夜夜夜久久久久| 好男人电影高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看影片大全网站| 两性夫妻黄色片| 亚洲熟女毛片儿| 在线观看午夜福利视频| 久久香蕉精品热| 麻豆av在线久日| 一二三四在线观看免费中文在| 欧美zozozo另类| 不卡av一区二区三区| 露出奶头的视频| 一区二区三区国产精品乱码| 国产精品爽爽va在线观看网站 | 欧美黑人欧美精品刺激| 黄色丝袜av网址大全| 久久婷婷人人爽人人干人人爱| 人妻久久中文字幕网| 巨乳人妻的诱惑在线观看| 久久草成人影院| 国产成人系列免费观看| 老司机午夜福利在线观看视频| 村上凉子中文字幕在线| 在线av久久热| 午夜福利欧美成人| 午夜成年电影在线免费观看| 国产免费av片在线观看野外av| 亚洲国产中文字幕在线视频| 两个人免费观看高清视频| 午夜福利成人在线免费观看| 日韩欧美一区视频在线观看| 日韩欧美在线二视频| avwww免费| 欧美中文综合在线视频| 国产精品久久久人人做人人爽| 精品久久久久久成人av| 免费av毛片视频| 一本精品99久久精品77| 十八禁网站免费在线| 中文字幕久久专区| 亚洲无线在线观看| 好看av亚洲va欧美ⅴa在| 黄色视频,在线免费观看| 好看av亚洲va欧美ⅴa在| 亚洲九九香蕉| 免费在线观看视频国产中文字幕亚洲| 欧美zozozo另类| 老熟妇乱子伦视频在线观看| 在线观看午夜福利视频| 欧美亚洲日本最大视频资源| 黄色片一级片一级黄色片| 波多野结衣高清无吗| 免费在线观看成人毛片| 久久久国产欧美日韩av| 亚洲精品国产一区二区精华液| 美女大奶头视频| 中出人妻视频一区二区| 男女做爰动态图高潮gif福利片| 亚洲国产中文字幕在线视频| 久久精品国产清高在天天线| 中文字幕久久专区| 精品久久久久久久毛片微露脸| 搡老熟女国产l中国老女人| 亚洲九九香蕉| 搞女人的毛片| 69av精品久久久久久| 成年版毛片免费区| 俺也久久电影网| 欧美三级亚洲精品| 国产麻豆成人av免费视频| 午夜福利在线在线| 可以免费在线观看a视频的电影网站| 国产在线观看jvid| 欧美另类亚洲清纯唯美| 国产一卡二卡三卡精品| 啦啦啦韩国在线观看视频| 国产亚洲欧美在线一区二区| 欧美乱码精品一区二区三区| 国产高清有码在线观看视频 | 一区二区三区高清视频在线| 亚洲人成伊人成综合网2020| 精华霜和精华液先用哪个| 99久久综合精品五月天人人| 美女大奶头视频| 美女扒开内裤让男人捅视频| 久久精品国产清高在天天线| 国产99久久九九免费精品| 亚洲中文av在线| 麻豆一二三区av精品| 欧美性猛交黑人性爽| 久久久久国产一级毛片高清牌| 一二三四在线观看免费中文在| 亚洲成人国产一区在线观看| 久久精品成人免费网站| 一本综合久久免费| 亚洲aⅴ乱码一区二区在线播放 | 亚洲,欧美精品.| 免费观看精品视频网站| 免费高清在线观看日韩| 国产精品乱码一区二三区的特点| 亚洲黑人精品在线| 国产亚洲精品久久久久久毛片| 国产亚洲av嫩草精品影院| 一进一出抽搐gif免费好疼| 精品无人区乱码1区二区| 波多野结衣av一区二区av| 最新美女视频免费是黄的| 宅男免费午夜| 国产精品野战在线观看| 欧美日韩乱码在线| 变态另类成人亚洲欧美熟女| 麻豆成人午夜福利视频| 男女午夜视频在线观看| 熟女少妇亚洲综合色aaa.| 国产av一区二区精品久久| 人人澡人人妻人| 欧美性猛交黑人性爽| 亚洲自偷自拍图片 自拍| 69av精品久久久久久| 99国产极品粉嫩在线观看| 亚洲免费av在线视频| 国产不卡一卡二| 日韩欧美国产在线观看| 国产一区二区激情短视频| 搡老熟女国产l中国老女人| 99国产精品一区二区蜜桃av| 亚洲精品在线美女| 99久久国产精品久久久| 丁香欧美五月| 国产高清videossex| 久久精品影院6| 中文资源天堂在线| 日本 欧美在线| 老汉色∧v一级毛片| 老司机午夜十八禁免费视频| 99re在线观看精品视频| 亚洲在线自拍视频| 久久中文字幕人妻熟女| 国产精品一区二区精品视频观看| 久久香蕉激情| 丁香欧美五月| 亚洲人成网站高清观看| 久久久精品欧美日韩精品| 欧美精品啪啪一区二区三区| 91国产中文字幕| 伊人久久大香线蕉亚洲五| 91麻豆av在线| 久久中文看片网| 久久 成人 亚洲| 日本a在线网址| 亚洲三区欧美一区| 人人妻人人澡人人看| 国产一区在线观看成人免费| 91老司机精品| 免费观看人在逋| 精品不卡国产一区二区三区| 三级毛片av免费| 国产又色又爽无遮挡免费看| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩亚洲综合一区二区三区_| av免费在线观看网站| 国产真实乱freesex| 在线视频色国产色| 黑人欧美特级aaaaaa片| 久久中文字幕一级| 亚洲三区欧美一区| e午夜精品久久久久久久| 少妇被粗大的猛进出69影院| 中文字幕高清在线视频| 国产精品 国内视频| 久久婷婷成人综合色麻豆| 最新在线观看一区二区三区| 国产精品久久久人人做人人爽| 脱女人内裤的视频| а√天堂www在线а√下载| 熟女电影av网| 国产视频内射| 亚洲成人国产一区在线观看| 韩国av一区二区三区四区| 久久精品91无色码中文字幕| 好男人电影高清在线观看| 亚洲avbb在线观看| 亚洲久久久国产精品| 女人爽到高潮嗷嗷叫在线视频| a级毛片在线看网站| 日本免费一区二区三区高清不卡| 中出人妻视频一区二区| www.999成人在线观看| 国产熟女xx| 久久精品成人免费网站| 国产精品久久久久久亚洲av鲁大| 国产精品永久免费网站| 久久久久久国产a免费观看| 女性被躁到高潮视频| 91国产中文字幕| 午夜福利18| 男女那种视频在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲精品一卡2卡三卡4卡5卡| 久99久视频精品免费| 99re在线观看精品视频| 脱女人内裤的视频| 亚洲黑人精品在线| 久久精品人妻少妇| 日本三级黄在线观看| 美女高潮到喷水免费观看| 少妇熟女aⅴ在线视频| 午夜福利视频1000在线观看| 久久人人精品亚洲av| 中文字幕最新亚洲高清| 成在线人永久免费视频| 国产精品亚洲一级av第二区| 精品一区二区三区四区五区乱码| 久热这里只有精品99| 日日摸夜夜添夜夜添小说| 亚洲aⅴ乱码一区二区在线播放 | 国产精品日韩av在线免费观看| 欧美日韩精品网址| 中文字幕精品亚洲无线码一区 | 成人亚洲精品av一区二区| 变态另类丝袜制服| 国产乱人伦免费视频| 午夜激情福利司机影院| 91麻豆精品激情在线观看国产| 91九色精品人成在线观看| 在线观看免费日韩欧美大片| 美女大奶头视频| 成人一区二区视频在线观看| 久久久久久久精品吃奶| 欧美+亚洲+日韩+国产| 国产高清videossex| 久久中文字幕人妻熟女| 亚洲一区二区三区色噜噜| 91麻豆av在线| 国产亚洲欧美在线一区二区| 亚洲五月色婷婷综合| 日韩一卡2卡3卡4卡2021年| 日韩欧美三级三区| 欧美精品啪啪一区二区三区| 日韩有码中文字幕| 婷婷六月久久综合丁香| 欧美乱码精品一区二区三区| 久久午夜亚洲精品久久| 亚洲国产看品久久| 国产成人精品久久二区二区免费| 麻豆av在线久日| 久久人人精品亚洲av| 人妻丰满熟妇av一区二区三区| cao死你这个sao货| 久久久久亚洲av毛片大全| 两个人看的免费小视频| 男女之事视频高清在线观看| 欧美一区二区精品小视频在线| 午夜成年电影在线免费观看| 成在线人永久免费视频| 亚洲一区二区三区色噜噜| 亚洲男人的天堂狠狠| 亚洲无线在线观看| 日本黄色视频三级网站网址| 精品欧美一区二区三区在线| 欧美日韩中文字幕国产精品一区二区三区| 黄色成人免费大全| 国产97色在线日韩免费| 99在线人妻在线中文字幕| 一级a爱视频在线免费观看| 欧美丝袜亚洲另类 | 免费电影在线观看免费观看| 日韩国内少妇激情av| 可以在线观看的亚洲视频| 欧美午夜高清在线| 欧美国产精品va在线观看不卡| av福利片在线| 亚洲成人国产一区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产麻豆成人av免费视频| 久久久久久久久久黄片| 老熟妇乱子伦视频在线观看| e午夜精品久久久久久久|