• <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
    免费看不卡的av| 亚洲国产av新网站| 男女无遮挡免费网站观看| 久久人人97超碰香蕉20202| 男人爽女人下面视频在线观看| 亚洲精品国产av成人精品| 亚洲欧美一区二区三区久久| 人妻人人澡人人爽人人| 欧美乱码精品一区二区三区| 久久久久人妻精品一区果冻| 午夜久久久在线观看| 男女高潮啪啪啪动态图| 日韩av免费高清视频| 高清av免费在线| 青春草亚洲视频在线观看| 午夜激情久久久久久久| 亚洲国产中文字幕在线视频| 最近中文字幕高清免费大全6| 在现免费观看毛片| 午夜福利乱码中文字幕| 亚洲精品aⅴ在线观看| 各种免费的搞黄视频| 久久青草综合色| 天堂俺去俺来也www色官网| 成人午夜精彩视频在线观看| 国产成人精品久久二区二区91 | 久久综合国产亚洲精品| 日韩中文字幕欧美一区二区 | 午夜福利网站1000一区二区三区| 夫妻午夜视频| 国产精品一二三区在线看| 日韩一区二区视频免费看| 亚洲七黄色美女视频| 色婷婷av一区二区三区视频| 无限看片的www在线观看| 精品亚洲乱码少妇综合久久| 天天添夜夜摸| 十分钟在线观看高清视频www| 亚洲成人手机| 亚洲美女搞黄在线观看| 女性生殖器流出的白浆| 国产精品秋霞免费鲁丝片| 极品少妇高潮喷水抽搐| 九草在线视频观看| 亚洲欧美成人综合另类久久久| 一区二区三区四区激情视频| 国产爽快片一区二区三区| 欧美激情极品国产一区二区三区| 亚洲国产av新网站| 国产精品熟女久久久久浪| tube8黄色片| 丝袜美腿诱惑在线| 2018国产大陆天天弄谢| 老司机在亚洲福利影院| 亚洲第一青青草原| 欧美日韩福利视频一区二区| 日日爽夜夜爽网站| 国产日韩欧美亚洲二区| 欧美xxⅹ黑人| 国产片内射在线| 最近的中文字幕免费完整| 午夜福利在线免费观看网站| 少妇的丰满在线观看| 99九九在线精品视频| 熟妇人妻不卡中文字幕| 国产成人精品久久久久久| 各种免费的搞黄视频| 日本欧美国产在线视频| 亚洲欧美色中文字幕在线| 中文字幕人妻丝袜制服| 国产成人欧美在线观看 | 国产野战对白在线观看| 欧美精品人与动牲交sv欧美| 国产视频首页在线观看| 女人被躁到高潮嗷嗷叫费观| 色网站视频免费| 一本色道久久久久久精品综合| av电影中文网址| 极品少妇高潮喷水抽搐| 日韩一本色道免费dvd| 欧美精品av麻豆av| 久久人人97超碰香蕉20202| 亚洲国产最新在线播放| 日日啪夜夜爽| 日韩av在线免费看完整版不卡| svipshipincom国产片| 少妇 在线观看| 久久亚洲国产成人精品v| 亚洲精品,欧美精品| 观看美女的网站| 欧美成人精品欧美一级黄| 欧美激情极品国产一区二区三区| av片东京热男人的天堂| 久久久久国产精品人妻一区二区| 深夜精品福利| 精品一区二区免费观看| 亚洲欧美中文字幕日韩二区| 国产一卡二卡三卡精品 | 亚洲欧美中文字幕日韩二区| 欧美黄色片欧美黄色片| 欧美日韩福利视频一区二区| 人妻人人澡人人爽人人| 另类亚洲欧美激情| 欧美日韩一级在线毛片| av不卡在线播放| 男人添女人高潮全过程视频| 国产黄频视频在线观看| 精品国产国语对白av| 午夜福利视频在线观看免费| 亚洲人成电影观看| 国精品久久久久久国模美| 亚洲美女搞黄在线观看| 高清黄色对白视频在线免费看| 国产亚洲午夜精品一区二区久久| 晚上一个人看的免费电影| 亚洲av日韩精品久久久久久密 | 伦理电影免费视频| 免费看av在线观看网站| 99精国产麻豆久久婷婷| 激情五月婷婷亚洲| 18禁观看日本| 伊人亚洲综合成人网| 制服人妻中文乱码| 国产在线一区二区三区精| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲精品一区二区精品久久久 | 男女无遮挡免费网站观看| 国产黄色免费在线视频| 欧美日韩视频高清一区二区三区二| 欧美成人精品欧美一级黄| 欧美最新免费一区二区三区| 精品午夜福利在线看| 侵犯人妻中文字幕一二三四区| 亚洲国产精品999| 哪个播放器可以免费观看大片| 精品一区二区三区av网在线观看 | 如日韩欧美国产精品一区二区三区| 亚洲色图综合在线观看| 亚洲精品一区蜜桃| 丝袜人妻中文字幕| 精品少妇内射三级| 考比视频在线观看| 亚洲国产av新网站| 国产日韩欧美视频二区| 成人亚洲欧美一区二区av| xxx大片免费视频| 亚洲精品第二区| 欧美黄色片欧美黄色片| 久久天躁狠狠躁夜夜2o2o | 久久国产亚洲av麻豆专区| 妹子高潮喷水视频| 久久人妻熟女aⅴ| 久久久久精品国产欧美久久久 | 亚洲精品av麻豆狂野| av在线app专区| 男的添女的下面高潮视频| 精品第一国产精品| 国产亚洲精品第一综合不卡| 久久性视频一级片| av不卡在线播放| 亚洲成人国产一区在线观看 | 99久久99久久久精品蜜桃| 欧美97在线视频| 亚洲人成电影观看| 免费在线观看视频国产中文字幕亚洲 | 曰老女人黄片| 国产在线视频一区二区| 免费久久久久久久精品成人欧美视频| 欧美日韩精品网址| 一边摸一边做爽爽视频免费| 欧美黑人精品巨大| 久久ye,这里只有精品| 欧美国产精品一级二级三级| 亚洲一区中文字幕在线| 亚洲精品国产区一区二| 一区福利在线观看| 男女边摸边吃奶| www.熟女人妻精品国产| 毛片一级片免费看久久久久| 亚洲一区二区三区欧美精品| 国产精品一国产av| 国产精品 欧美亚洲| 亚洲综合精品二区| 久久久久久免费高清国产稀缺| 国产又色又爽无遮挡免| 国产成人精品在线电影| 亚洲美女黄色视频免费看| 精品久久蜜臀av无| 亚洲精品久久午夜乱码| 在线观看www视频免费| 男女午夜视频在线观看| 女性生殖器流出的白浆| 亚洲精品自拍成人| 亚洲国产欧美在线一区| 久久久精品区二区三区| 午夜福利网站1000一区二区三区| 9191精品国产免费久久| 十八禁网站网址无遮挡| 一级片'在线观看视频| 另类精品久久| 国产精品久久久久久人妻精品电影 | 最近2019中文字幕mv第一页| 五月开心婷婷网| 一边亲一边摸免费视频| av线在线观看网站| 黄片无遮挡物在线观看| 亚洲国产精品成人久久小说| 色网站视频免费| 在线观看国产h片| 黑人巨大精品欧美一区二区蜜桃| 久久性视频一级片| 九色亚洲精品在线播放| 午夜福利在线免费观看网站| xxxhd国产人妻xxx| 久久婷婷青草| 中国三级夫妇交换| www.自偷自拍.com| 日本色播在线视频| 9热在线视频观看99| av有码第一页| 国产精品欧美亚洲77777| 亚洲美女搞黄在线观看| 国产 一区精品| av不卡在线播放| 欧美日本中文国产一区发布| 国产片内射在线| 涩涩av久久男人的天堂| 欧美日韩福利视频一区二区| 一区二区三区激情视频| 国产成人精品无人区| 久久99热这里只频精品6学生| 国产精品久久久av美女十八| 日本av免费视频播放| 99久国产av精品国产电影| 国产色婷婷99| 免费观看人在逋| 亚洲第一av免费看| 亚洲欧美成人精品一区二区| 欧美xxⅹ黑人| 国产精品人妻久久久影院| 久久97久久精品| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 亚洲在久久综合| 精品免费久久久久久久清纯 | 午夜影院在线不卡| 精品国产乱码久久久久久男人| 日韩制服骚丝袜av| av在线播放精品| 婷婷成人精品国产| a级片在线免费高清观看视频| 日韩大码丰满熟妇| 久久影院123| 国产亚洲av高清不卡| 色网站视频免费| 观看美女的网站| 午夜免费观看性视频| 精品一区二区免费观看| 亚洲专区中文字幕在线 | 在线亚洲精品国产二区图片欧美| 91国产中文字幕| 色婷婷久久久亚洲欧美| 精品国产乱码久久久久久男人| 国产99久久九九免费精品| 一区福利在线观看| 又粗又硬又长又爽又黄的视频| 日韩制服丝袜自拍偷拍| 99香蕉大伊视频| 亚洲欧美日韩另类电影网站| 天堂中文最新版在线下载| 国产av国产精品国产| 午夜免费男女啪啪视频观看| 高清av免费在线| av线在线观看网站| 午夜老司机福利片| 免费看av在线观看网站| 91国产中文字幕| 欧美日韩精品网址| 精品福利永久在线观看| 不卡av一区二区三区| 久久女婷五月综合色啪小说| 午夜福利视频精品| www.av在线官网国产| svipshipincom国产片| 波多野结衣av一区二区av| 国产极品粉嫩免费观看在线| 久久人人爽人人片av| 亚洲欧美一区二区三区黑人| 一本久久精品| 成人亚洲精品一区在线观看| 中文字幕色久视频| 亚洲欧美一区二区三区久久| 爱豆传媒免费全集在线观看| 热99国产精品久久久久久7| 丝袜在线中文字幕| 亚洲国产精品一区二区三区在线| 欧美中文综合在线视频| 五月天丁香电影| 欧美变态另类bdsm刘玥| 中文字幕人妻熟女乱码| 一区二区三区乱码不卡18| 另类亚洲欧美激情| 久久久精品国产亚洲av高清涩受| 2021少妇久久久久久久久久久| 国产淫语在线视频| 国产精品久久久久久人妻精品电影 | 国产成人精品福利久久| 精品午夜福利在线看| 亚洲精品国产区一区二| 肉色欧美久久久久久久蜜桃| 国产日韩一区二区三区精品不卡| 欧美日韩亚洲高清精品| 一级毛片黄色毛片免费观看视频| 午夜久久久在线观看| 飞空精品影院首页| 天天躁狠狠躁夜夜躁狠狠躁| 免费观看av网站的网址| 在线观看免费日韩欧美大片| 男男h啪啪无遮挡| 一个人免费看片子| 国产精品成人在线| 亚洲精品久久午夜乱码| 午夜福利视频在线观看免费| 亚洲一卡2卡3卡4卡5卡精品中文| a级毛片黄视频| 深夜精品福利| av免费观看日本| 亚洲国产最新在线播放| 在线观看一区二区三区激情| 夫妻性生交免费视频一级片| 久久久久精品国产欧美久久久 | 久久久国产精品麻豆| 一级片'在线观看视频| 在线免费观看不下载黄p国产| 秋霞伦理黄片| 日韩伦理黄色片| 不卡av一区二区三区| 91精品三级在线观看| 狂野欧美激情性xxxx| 啦啦啦啦在线视频资源| 美女脱内裤让男人舔精品视频| 国产一区二区 视频在线| 老熟女久久久| 国产极品天堂在线| 亚洲成人手机| 看免费成人av毛片| xxx大片免费视频| 侵犯人妻中文字幕一二三四区| 一区福利在线观看| 一区在线观看完整版| 国产精品久久久久成人av| 亚洲成人国产一区在线观看 | 操美女的视频在线观看| 欧美日韩av久久| 国产在线免费精品| 各种免费的搞黄视频| 在线免费观看不下载黄p国产| 黑丝袜美女国产一区| 国产亚洲最大av| 精品人妻一区二区三区麻豆| 亚洲激情五月婷婷啪啪| 日日爽夜夜爽网站| 丝袜喷水一区| 制服诱惑二区| 久热爱精品视频在线9| 丝袜美腿诱惑在线| 免费观看性生交大片5| 久久精品国产a三级三级三级| 国产精品av久久久久免费| 黑人欧美特级aaaaaa片| 国产又爽黄色视频| 久久精品亚洲av国产电影网| 在线观看三级黄色| 国产欧美日韩一区二区三区在线| 亚洲精品aⅴ在线观看| 亚洲婷婷狠狠爱综合网| 黄色一级大片看看| 精品第一国产精品| 秋霞伦理黄片| 亚洲一级一片aⅴ在线观看| 女的被弄到高潮叫床怎么办| 1024视频免费在线观看| 欧美乱码精品一区二区三区| 欧美日本中文国产一区发布| 搡老乐熟女国产| 欧美在线黄色| 日本wwww免费看| 国产一级毛片在线| 国产女主播在线喷水免费视频网站| 色播在线永久视频| 国产一卡二卡三卡精品 | 高清欧美精品videossex| 国产在线视频一区二区| 国产亚洲最大av| 免费少妇av软件| 日韩一卡2卡3卡4卡2021年| 女人被躁到高潮嗷嗷叫费观| 欧美亚洲 丝袜 人妻 在线| 五月天丁香电影| 在线天堂中文资源库| 男女国产视频网站| 国产日韩欧美在线精品| 又大又黄又爽视频免费| 亚洲欧美激情在线| 亚洲精品aⅴ在线观看| 午夜福利影视在线免费观看| 女性生殖器流出的白浆| 一区二区三区四区激情视频| av在线老鸭窝| 十八禁人妻一区二区| 亚洲欧美一区二区三区久久| 黄片小视频在线播放| 亚洲一级一片aⅴ在线观看| 九九爱精品视频在线观看| 黄片无遮挡物在线观看| 久久久久精品性色| 香蕉丝袜av| 中文字幕av电影在线播放| h视频一区二区三区| 精品一区二区三区四区五区乱码 | 国产精品欧美亚洲77777| 中文字幕亚洲精品专区| 母亲3免费完整高清在线观看| 精品卡一卡二卡四卡免费| 在线观看一区二区三区激情| 免费观看a级毛片全部| 老司机亚洲免费影院| 免费高清在线观看日韩| 国产精品女同一区二区软件| 在线观看免费视频网站a站| 国产97色在线日韩免费| 国产一级毛片在线| 国产极品天堂在线| av不卡在线播放| 精品国产乱码久久久久久男人| 啦啦啦在线免费观看视频4| 久久婷婷青草| 超色免费av| 欧美日韩国产mv在线观看视频| 十八禁高潮呻吟视频| 99精品久久久久人妻精品| 日韩 亚洲 欧美在线| 黄色一级大片看看| av天堂久久9| 欧美av亚洲av综合av国产av | 黄色视频在线播放观看不卡| 亚洲人成77777在线视频| 久久影院123| 精品一区二区三卡| 美女主播在线视频| 国产精品国产三级国产专区5o| 大话2 男鬼变身卡| 欧美人与性动交α欧美软件| 人妻一区二区av| 搡老岳熟女国产| 国产97色在线日韩免费| 亚洲伊人色综图| 日本午夜av视频| 97精品久久久久久久久久精品| 国产福利在线免费观看视频| 中文字幕高清在线视频| 黄色视频不卡| 久久99热这里只频精品6学生| 国产免费福利视频在线观看| 国产精品一区二区在线观看99| 久久综合国产亚洲精品| 午夜福利一区二区在线看| 日本猛色少妇xxxxx猛交久久| 国产深夜福利视频在线观看| 亚洲熟女精品中文字幕| 国产老妇伦熟女老妇高清| 久久国产精品大桥未久av| 大片免费播放器 马上看| 深夜精品福利| 中文字幕最新亚洲高清| 国产精品香港三级国产av潘金莲 | 熟女少妇亚洲综合色aaa.| h视频一区二区三区| 国产片内射在线| 国产成人午夜福利电影在线观看| 免费久久久久久久精品成人欧美视频| 丁香六月天网| 精品久久久久久电影网| 免费黄色在线免费观看| 久久国产精品大桥未久av| 日韩欧美精品免费久久| 看十八女毛片水多多多| 亚洲情色 制服丝袜| 亚洲,欧美精品.| 欧美精品av麻豆av| 久热这里只有精品99| 菩萨蛮人人尽说江南好唐韦庄| 黄色 视频免费看| 成人黄色视频免费在线看| 最近2019中文字幕mv第一页| 男女免费视频国产| 最近手机中文字幕大全| 久久女婷五月综合色啪小说| 久久人人爽av亚洲精品天堂| 国产不卡av网站在线观看| 午夜久久久在线观看| 天天影视国产精品| 18禁裸乳无遮挡动漫免费视频| 悠悠久久av| 国产 一区精品| 婷婷成人精品国产| 欧美日韩视频精品一区| 韩国高清视频一区二区三区| 久久久亚洲精品成人影院| 亚洲一区中文字幕在线| 天天添夜夜摸| 久久久久久久精品精品| 在线天堂中文资源库| 成人漫画全彩无遮挡| 久久久欧美国产精品| 国产精品女同一区二区软件| 亚洲久久久国产精品| 成人18禁高潮啪啪吃奶动态图| 日韩电影二区| 搡老乐熟女国产| 亚洲精华国产精华液的使用体验| 男女国产视频网站| 亚洲国产精品999| 一区二区三区四区激情视频| 亚洲在久久综合| 成人漫画全彩无遮挡| 亚洲成人国产一区在线观看 | 又粗又硬又长又爽又黄的视频| 日韩 欧美 亚洲 中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲精品久久成人aⅴ小说| 国产男人的电影天堂91| 国产又爽黄色视频| 中文精品一卡2卡3卡4更新| 午夜影院在线不卡| 亚洲精品乱久久久久久| 国产又爽黄色视频| 大码成人一级视频| 亚洲第一av免费看| 久久久久视频综合| 亚洲成人免费av在线播放| av又黄又爽大尺度在线免费看| avwww免费| 国产精品 欧美亚洲| 亚洲,欧美精品.| 午夜免费鲁丝| 女人久久www免费人成看片| 好男人视频免费观看在线| 大陆偷拍与自拍| 精品国产国语对白av| 国产精品一区二区在线观看99| 精品卡一卡二卡四卡免费| 欧美激情高清一区二区三区 | 色网站视频免费| 老司机深夜福利视频在线观看 | 久久久久久人人人人人| 成人国产麻豆网| 免费黄频网站在线观看国产| av又黄又爽大尺度在线免费看| 欧美在线黄色| 男女边吃奶边做爰视频| 免费在线观看完整版高清| 丰满迷人的少妇在线观看| 欧美黑人欧美精品刺激| 国产成人免费无遮挡视频| 人人妻人人澡人人爽人人夜夜| 一级毛片我不卡| 男女之事视频高清在线观看 | 天美传媒精品一区二区| 看免费成人av毛片| 一级毛片我不卡| 在线看a的网站| 国产精品免费视频内射| 黄片无遮挡物在线观看| 999久久久国产精品视频| 操美女的视频在线观看| 欧美97在线视频| 午夜久久久在线观看| 午夜福利一区二区在线看| 伦理电影免费视频| 久久精品国产亚洲av涩爱| 日韩欧美精品免费久久| 日韩av不卡免费在线播放| 欧美成人精品欧美一级黄| 一本色道久久久久久精品综合| 赤兔流量卡办理| 国精品久久久久久国模美| 成人亚洲精品一区在线观看| 天天躁日日躁夜夜躁夜夜| 热re99久久精品国产66热6| 免费观看人在逋| 久久久久精品性色| 视频在线观看一区二区三区| 建设人人有责人人尽责人人享有的| 免费人妻精品一区二区三区视频| 亚洲国产欧美日韩在线播放| 一区二区av电影网| 久久国产精品男人的天堂亚洲| 国产人伦9x9x在线观看| 国产在线免费精品| 成人毛片60女人毛片免费| 亚洲欧美中文字幕日韩二区| 久久免费观看电影| 亚洲精品国产区一区二| 欧美日韩一级在线毛片| 国产野战对白在线观看| 大片电影免费在线观看免费| 国产不卡av网站在线观看| 国产1区2区3区精品| 一级毛片黄色毛片免费观看视频| 国产黄色视频一区二区在线观看| 亚洲熟女精品中文字幕| 日韩免费高清中文字幕av| 国产亚洲av片在线观看秒播厂| 爱豆传媒免费全集在线观看|