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

    基于改進(jìn)牛頓法的航空發(fā)動機(jī)起動建模技術(shù)

    2018-06-20 10:19:34王繼強(qiáng)何小龍胡忠志
    航空發(fā)動機(jī) 2018年1期
    關(guān)鍵詞:慢車起動機(jī)部件

    姜 波,王繼強(qiáng),何小龍,胡忠志

    (南京航空航天大學(xué)江蘇省航空動力系統(tǒng)重點實驗室,南京210016)

    0 引言

    航空發(fā)動機(jī)由靜止逐漸加速至慢車狀態(tài)這一過程稱為起動過程。這個過程發(fā)動機(jī)氣動熱力參數(shù)變換劇烈,如何保證這個過程成功、快速的完成對延長發(fā)動機(jī)壽命、保證飛機(jī)安全有著重要的意義。對該過程進(jìn)行數(shù)學(xué)抽象,建立起動過程的數(shù)值計算模型不僅能夠為設(shè)計人員提供許多有價值的信息,而且能夠讓設(shè)計人員在發(fā)動機(jī)設(shè)計的初期就展開合理的起動過程控制規(guī)律優(yōu)化研究。因而建立能夠準(zhǔn)確反映發(fā)動機(jī)起動過程的數(shù)值仿真模型對發(fā)動機(jī)的理論研究及工程應(yīng)用都有著極大地意義。

    在20世紀(jì)80年代,國外已經(jīng)展開了發(fā)動機(jī)慢車以下狀態(tài)的氣動熱力學(xué)建模研究,M.Yunis和R.K.A-grawal在1980年基于流量平衡的假設(shè),通過外推方法獲得低轉(zhuǎn)速部件特性從而建立起了單轉(zhuǎn)子發(fā)動機(jī)的起動模型[1]。在90年代NASA的Chappell采用慢車以上建模的方法,建立了包含高空起動、風(fēng)車起動和地面起動的AVEST-V3模型[2]。此外在2001年和2007年,英國弗吉尼亞理工大學(xué)[3]英國克蘭菲爾德大學(xué)[4]的學(xué)者在上述研究的基礎(chǔ)上針對部件特性外推公式及零轉(zhuǎn)速附近特性修正展開了深入的研究。國內(nèi)學(xué)者對此研究展開的時間較晚,但也取得了大量的研究成果,如:基于大量發(fā)動機(jī)試車數(shù)據(jù),使用系統(tǒng)辨識的方法建立的實時起動模型[5-8];基于容積動力學(xué)和部件特性外推的方法建立雙軸混合排氣的渦扇發(fā)動機(jī)的地面及空中起動模型[9-10],利用壓氣機(jī)及渦輪特性數(shù)據(jù)的脊柱特性建立能夠準(zhǔn)確模擬單軸發(fā)動機(jī)風(fēng)車起動的部件級模型[11]等。但這些研究大都基于大量的試驗數(shù)據(jù)并且在起動初始點需要根據(jù)經(jīng)驗反復(fù)的試湊出十分準(zhǔn)確的初始值,工作量大且過程較為復(fù)雜。

    本文利用NASA公開的PW公司JT9D大涵道比雙轉(zhuǎn)子渦扇發(fā)動機(jī)慢車以上部件級模型,通過部件特性外推、起動附件建模、修正部件總壓恢復(fù)系數(shù)和改進(jìn)模型求解算法建立了JT9D發(fā)動機(jī)起動至最大轉(zhuǎn)速整個過程的仿真模型,并通過仿真結(jié)果驗證了模型的準(zhǔn)確性及改進(jìn)求解算法的有效性。

    1 JT9D發(fā)動機(jī)部件級模型

    建立的起動過程模型是基于美國PW公司于20世紀(jì)60年代年設(shè)計制造的大涵道比雙轉(zhuǎn)子渦扇發(fā)動機(jī)JT9D的慢車以上狀態(tài)的部件級模型,具體結(jié)構(gòu)如圖1所示,具體的建模過程這里不再贅述,參照文獻(xiàn)[12]。

    圖1 JT9D發(fā)動機(jī)部件結(jié)構(gòu)

    2 雙軸渦扇發(fā)動機(jī)起動過程

    雙軸渦扇發(fā)動機(jī)在起動時沒有空氣流經(jīng)整臺發(fā)動機(jī),因此燃燒室無法點火。只有當(dāng)轉(zhuǎn)子達(dá)到一定的轉(zhuǎn)速后,燃燒室內(nèi)的氣流才能建立穩(wěn)定燃燒所需的氣流壓力和溫度。因而發(fā)動機(jī)地面起動過程必須依靠起動機(jī)帶轉(zhuǎn),一般情況下,雙軸渦扇發(fā)動機(jī)的起動過程(如圖2所示)分為以下3個階段。

    圖2 雙軸渦扇發(fā)動機(jī)起動過程扭矩特性

    第1階段:燃燒室未供油,起動機(jī)帶動高壓轉(zhuǎn)子,起動機(jī)扭矩MST大于起動過程阻力矩MC/ηm,高壓轉(zhuǎn)子轉(zhuǎn)速逐漸增加到n1。

    第2階段:燃燒室開始供油并點火,渦輪開始產(chǎn)生力矩,高壓轉(zhuǎn)子轉(zhuǎn)速逐漸增加到n2。

    第3階段:起動機(jī)脫開,渦輪產(chǎn)生足夠的力矩并將高壓轉(zhuǎn)子轉(zhuǎn)速獨立帶轉(zhuǎn)至n3。

    3 起動建模關(guān)鍵問題及技術(shù)

    3.1 低轉(zhuǎn)速部件特性獲取

    發(fā)動機(jī)的部件特性在高轉(zhuǎn)速時差別較大,但在低轉(zhuǎn)速范圍內(nèi)卻有著很大的相似性[1]。因而,采用部件特性指數(shù)外推的方法,能夠在缺少試驗特性的情況下得到較為合理的低轉(zhuǎn)速部件特性。采用的指數(shù)外推公式如下

    式中:p、q、r為外推指數(shù);cor表示換算狀態(tài);ref表示參考狀態(tài);Wc為換算流量;N為轉(zhuǎn)速;PR為增壓比(渦輪特性為膨脹比ER);η為效率。以風(fēng)扇為例,選取換算轉(zhuǎn)速為1的特性線為參考線。利用原有的特性數(shù)據(jù)通過不斷變換外推指數(shù)的數(shù)值,使得不同換算轉(zhuǎn)速線上最高效率狀態(tài)點連成的線近似為線性(效率僅在低轉(zhuǎn)速Ncor<0.8時成線性)進(jìn)而確定外推指數(shù)。不同外推指數(shù)下風(fēng)扇流量、壓比、效率、隨換算轉(zhuǎn)速變化分別如圖3~5所示,不同外推指數(shù)下的外推結(jié)果對比如圖6所示。

    圖3 不同外推指數(shù)下風(fēng)扇流量隨換算轉(zhuǎn)速變化

    圖4 不同外推指數(shù)下風(fēng)扇壓比隨換算轉(zhuǎn)速變化

    圖5 不同外推指數(shù)下風(fēng)扇效率隨換算轉(zhuǎn)速變化

    在確定了外推指數(shù)后,以換算轉(zhuǎn)速為0.6的試驗特性線為參考,外推換算轉(zhuǎn)速為0.5的特性線并與原有試驗特性對比,結(jié)果顯示:按照不同換算轉(zhuǎn)速線上的最高效率點連成的線近似為線性,所選取的外推指數(shù)能更好的反映低轉(zhuǎn)速的部件特性。據(jù)此確定外推指數(shù)見表1。

    圖6 不同外推指數(shù)下的外推結(jié)果對比

    表1 部件特性外推指數(shù)

    在確定部件特性的外推指數(shù)后,以各部件試驗的最低狀態(tài)換算轉(zhuǎn)速線為外推參考線,外推各部件的風(fēng)扇換算流量-壓比特性、風(fēng)扇換算流量-效率特性、高壓渦輪換算流量-膨脹比、高壓渦輪換算流量效率特性特性分別如圖7~10所示。

    圖7 風(fēng)扇換算流量-壓比特性

    3.2 起動機(jī)建模

    JT9D發(fā)動機(jī)使用美國漢密爾頓公司提供的PS-700空氣渦輪起動機(jī)[13],主要參數(shù)見表2。

    空氣渦輪起動機(jī)在發(fā)動機(jī)起動過程的輸出扭矩與高壓轉(zhuǎn)子轉(zhuǎn)速近似為線性關(guān)系[14]。起動機(jī)扭矩特性如圖11所示。

    圖8 風(fēng)扇換算流量-效率特性

    圖9 高壓渦輪換算流量-膨脹比特性

    表2 PS-700空氣渦輪起動機(jī)主要參數(shù)

    根據(jù)扭矩特性按照如下公式建立起動機(jī)模型

    式中:P為起動機(jī)輸出功率;n2為起動機(jī)脫開時高壓轉(zhuǎn)子轉(zhuǎn)速;N為高壓轉(zhuǎn)子轉(zhuǎn)速;T為起動機(jī)輸出扭矩。

    圖11 起動機(jī)扭矩特性

    3.3 改進(jìn)牛頓法及其應(yīng)用

    在求解JT9D發(fā)動機(jī)共同工作非線性方程組時NASA提供了1種傳統(tǒng)的Newton-Raphson迭代求解方法。此方法在求解發(fā)動機(jī)慢車以上狀態(tài)共同工作非線性方程組時能夠在選定較為準(zhǔn)確的迭代初始值情況下達(dá)到快速收斂的效果[15-17]。但發(fā)動機(jī)起動過程部件特性的相對值變化劇烈,即使在給出相對準(zhǔn)確的迭代初始值情況下也需要多次求解才能達(dá)到收斂精度的要求。并且每次迭代求解共同工作非線性方程組時計算機(jī)都需要對模型進(jìn)行氣動熱力計算,然后再對平衡方程偏導(dǎo)數(shù)組成的Jacobian矩陣進(jìn)行中心差分,龐大的計算量降低了模型的執(zhí)行效率,使得計算時間變長,模型的實時性降低。

    本文所述的改進(jìn)牛頓法是參考廣泛用于國外先進(jìn)數(shù)值仿真程序,如美國航空航天局NASA的NPSS(Numerical Propulsion Simulation System)[18]、荷蘭國家航空航天實驗室(NLR)的GSP(Gas turbine Simulation Program)[19]等在求解非線性方程組時使用的Broyden法[20]:即在傳統(tǒng)Newton-Raphson方法基礎(chǔ)上改變Jacobian矩陣的迭代更新方式,進(jìn)而減少模型的氣動熱力計算次數(shù),提高模型的仿真求解速度。計算方法如下:

    對于發(fā)動機(jī)的共同工作非線性方程組

    式中:x=(x1,x2,…,xn)T,n 為平衡方程數(shù),xi為發(fā)動機(jī)的某參數(shù);εi為平衡方程的誤差量;fi(x)代表不同的平衡方程,且F(x)=(f1(x1),f2(x2),…,fn(xn))T,Jacobian矩陣為

    其中:偏導(dǎo)數(shù)項在模型求解計算時需根據(jù)中心差分求得。根據(jù)Newton-Raphson求解方法

    式中:k為迭代次數(shù),此時構(gòu)造如下矩陣

    根據(jù)Broyden法求解原理[20],Jacobian矩陣迭代更新的修正項為

    使用Sherman-Morrison公式解析的求Jacobian矩陣的逆,從而可得Jacobian矩陣在更新時

    當(dāng) |F(xk)|<δ(設(shè)定誤差允許量)時,認(rèn)為迭代收斂。

    這種改進(jìn)的牛頓法在求解發(fā)動機(jī)共同工作非線性方程組時第1步仍采用中心差分的方法得到Jacobian矩陣,而后續(xù)的迭代更新只需通過計算ΔJk就可完成。這樣模型的動態(tài)仿真過程就避免了每次通過中心差分方法重新計算Jacobian矩陣而帶來的龐大計算量,從而大大提高了需要多次迭代求解發(fā)動機(jī)共同工作方程的起動過程模型的執(zhí)行效率。

    3.4 冷轉(zhuǎn)動過程模型

    發(fā)動機(jī)在點火前的工作狀態(tài)即起動過程的第一階段稱為冷轉(zhuǎn)動過程。這一過程渦輪和壓氣機(jī)都是負(fù)載,建模方法如下:

    首先計算起動機(jī)和壓氣機(jī)的扭矩,然后根據(jù)轉(zhuǎn)子動力學(xué)計算高壓轉(zhuǎn)子轉(zhuǎn)速

    式中:ΔT為起動機(jī)扭矩與負(fù)載扭矩差;J為高壓轉(zhuǎn)子轉(zhuǎn)動慣量;N為高壓轉(zhuǎn)子轉(zhuǎn)速;t為時間。

    之后,確定1個合適的點火狀態(tài),根據(jù)此時的轉(zhuǎn)速 nR,WR和壓比 PR利用式(13)、(14)計算對應(yīng)時刻下的參數(shù)

    最后,低壓轉(zhuǎn)子轉(zhuǎn)速及其他截面狀態(tài)參數(shù)按高壓轉(zhuǎn)子轉(zhuǎn)速進(jìn)行線性插值即可得到。

    3.5 部件總壓恢復(fù)系數(shù)修正

    NASA公開的JT9D發(fā)動機(jī)模型各部件的總壓恢復(fù)系數(shù)都取了定值,這在慢車以上狀態(tài)建模時是適用的。但在發(fā)動機(jī)的起動過程中由于各部件的流量、壓比、轉(zhuǎn)速都很小,真實的總壓恢復(fù)系數(shù)更接近于1,并且隨著轉(zhuǎn)速的增大而逐漸減小[9]。因而在建立起動模型時,各部件的總壓恢復(fù)系數(shù)在點火前取近似略小于1的常數(shù),而在點火狀態(tài)和慢車狀態(tài)之間時按照轉(zhuǎn)速進(jìn)行線性插值

    式中:a為點火狀態(tài);b為慢車狀態(tài);σ為總壓恢復(fù)系數(shù)。

    4 仿真結(jié)果分析

    NASA公開的JT9D發(fā)動機(jī)模型的仿真范圍為:70%NHmax~100%NHmax。按照上述建模方法,補充建立包含起動過程的全狀態(tài)性能模型,并在地面標(biāo)準(zhǔn)狀況下給定供油規(guī)律(如圖12所示),仿真結(jié)果分別如圖13~15所示。

    圖12 供油曲線

    圖13 轉(zhuǎn)速變化曲線

    圖14 推力變化曲線

    圖15 燃燒室出口溫度變化曲線

    由仿真結(jié)果可以看出,起動過程的第1階段燃燒室沒有供油。高壓轉(zhuǎn)子由起動機(jī)帶動,轉(zhuǎn)速緩慢上升并帶動低壓轉(zhuǎn)子轉(zhuǎn)速上升,此時推力也緩慢上升,燃燒室出口溫度基本不變。在第16 s時開始供油,轉(zhuǎn)子轉(zhuǎn)速和燃燒室出口溫度急劇上升至慢車狀態(tài)。在第40 s以后增加供油,發(fā)動機(jī)由慢車狀態(tài)逐漸運行到最大轉(zhuǎn)速狀態(tài),轉(zhuǎn)子轉(zhuǎn)速、推力及燃燒室出口溫度都增加至100%。整個模型能很好的反映起動過程的各階段并且與原有慢車以上狀態(tài)模型有效地連接,驗證了建模方法的有效性。

    為了進(jìn)一步驗證建模方法的準(zhǔn)確性,將模型的仿真結(jié)果與僅有的ICAO(國際民航組織)公布的JT9D-7J發(fā)動機(jī)數(shù)據(jù)進(jìn)行對比,結(jié)果見表3。

    表3 7%推力狀態(tài)下(52%NHmax)模型仿真結(jié)果驗證

    各部件總壓恢復(fù)系數(shù)修正前(修正前采用固定值σb)與修正后高壓轉(zhuǎn)子轉(zhuǎn)速的仿真結(jié)果對比如圖16所示,燃燒室點火至慢車狀態(tài)過程的高壓轉(zhuǎn)子轉(zhuǎn)速最大相對誤差為1%。

    圖16 總壓恢復(fù)系數(shù)修正前后仿真結(jié)果對比

    同時,模型在點火后由各部件特性迭代求解的初始點計算時間最長,最難收斂。因而在此處用改進(jìn)牛頓法進(jìn)行驗證效果最為明顯。

    單個步長Jacobian矩陣計算時間見表4。表4計算的硬件環(huán)境為:AMD A6-3400M(CPU),主頻1.4 G Hz,內(nèi)存4 GB。由結(jié)果可以看出改進(jìn)牛頓法的求解速度要明顯高于傳統(tǒng)的Newton-Raphson方法。并且,基于改進(jìn)牛頓法的原理,當(dāng)發(fā)動機(jī)共同工作方程和獨立變量個數(shù)增加以及迭代初始值偏離真值程度越大時,這種效果將會變得更加明顯。

    表4 單個步長的Jacobian矩陣計算時間

    5 結(jié)論

    通過對JT9D發(fā)動機(jī)部件級起動過程建模、使用改進(jìn)牛頓法替換傳統(tǒng)的Newton-Raphson方法仿真求解,得到以下結(jié)論:

    (1)采用部件特性指數(shù)外推法并使用不同換算轉(zhuǎn)速線上最高效率狀態(tài)點連成的線,近似為線性來確定外推指數(shù),可以得到合理的各部件低轉(zhuǎn)速部件特性;

    (2)起動附件建模及轉(zhuǎn)子冷轉(zhuǎn)動過程建模能夠有效地幫助建立完整的起動過程模型。

    (3)發(fā)動機(jī)起動過程各部件總壓恢復(fù)系數(shù)按照轉(zhuǎn)速變化修正能夠有效的提高模型的精度。

    (4)采用改進(jìn)牛頓法替代傳統(tǒng)的Newton-Raphson方法可以改善模型的執(zhí)行效率,大幅提高模型的仿真求解速度。

    [1]Agrawal R K,Yunis M.A generalized mathematical model to estimate gas turbine starting characteristics[J].Physical Review B,1981(1):5668-5687.

    [2]Chappell M A,Mclaughlin P W.Approach of modeling continuous turbine engine operation from startup to shut down[J].Journal of Propulsion&Power,2015,9(3):466-471.

    [3]Sexton W R.A new method to control turbofan engine starting by varying compressor surge valve bleed[D].Virginia:Virginia Polytechnic Institute and State University,2001.

    [4]Howard,Jason.Sub-idle modeling of gas turbines:altitude relight and wind milling.[D].Cranfield:Cranfield University,2007.

    [5]于達(dá)仁,郭鈺鋒,牛軍,等.渦噴發(fā)動機(jī)風(fēng)車啟動工況的神經(jīng)網(wǎng)絡(luò)建模[J].推進(jìn)技術(shù),2001,22(3):183-186.YU Daren,GUO Yufeng,NIU Jun,et al.Turbojet modeling in wind milling based on radial basis function networks[J].Journal of Propulsion Technology,2001,22(3):183-186.(in Chinese)

    [6]卓剛.航空發(fā)動機(jī)智能建模與故障診斷研究[D].南京:南京航空航天大學(xué),2004.ZHUO Gang.The research on aeroengine intelligent modeling and fault diagnosis[D].Nanjing:Nanjing University of Aeronautics and Astronautics.( in Chinese)

    [7]劉建勛,王劍影,李應(yīng)紅,等.某型發(fā)動機(jī)起動模型的支持向量機(jī)辨識及應(yīng)用[J].推進(jìn)技術(shù),2004,25(5):401-404.LIU Jianxun,WANG Jianying,LI Yinghong,et al.An identification model of aero engine starting based on support vector machine and its application[J].Journal of Propulsion Technology,2004,25(5):401-404.(in Chinese)

    [8]李永進(jìn),張海波,張?zhí)旌?一種考慮非線性余項的機(jī)載發(fā)動機(jī)自適應(yīng)模型建立及其在尋優(yōu)控制中的應(yīng)用 [J].推進(jìn)技術(shù),2016,37(1):172-180.LI Yongjin,ZHANG Haibo,ZHANG Tianhong.Establishment and application in performance seeking control of an on-board adaptive aero-engine model considering nonlinear remainders[J].Journal of Propulsion Technology,2016,37(1):172-180.(in Chinese)

    [9]周文祥,黃金泉,竇建平.渦扇發(fā)動機(jī)部件級起動模型[J].航空動力學(xué)報,2006,21(2):248-253.ZHOU Wenxiang,HUANG Jinquan,DOU Jianping.Development of component-level startup model for a turbofan engine[J].Journal of Aerospace Power.2006,21(2):248-253.(in Chinese)

    [10]周文祥,黃金泉.渦扇發(fā)動機(jī)高空起動模型研究[J].航空動力學(xué)報,2007,22(8):1384-1390.ZHOU Wenxiang,HUANG Jinquan.Research on the startup model of turbofan engine at high altitude[J].Journal of Aerospace Power,2007,22(8):1384-1390.(in Chinese)

    [11]施洋,屠秋野,嚴(yán)紅明,等.一種航空發(fā)動機(jī)全狀態(tài)性能模型[J].航空動力學(xué)報,2017,32(2):373-381.SHI Yang,TU Qiuye,YAN Hongming,et al.A full states performance model for aero engine[J].Journal of Aerospace Power,2007,22(8):1384-1390.(in Chinese)

    [12]AIAA.A process for the creation of T-MATS propulsion system models from NPSS data[C]//.AIAA Joint Propulsion Conference.Cleveland,OH.2014.

    [13]胡曉煜.世界中小型航空發(fā)動機(jī)手冊[M].北京.航空工業(yè)出版社,2006:443-445.HU Xiaoyu.World minitype aeroengine manual[M]Beijing:Aviation Industry Press,2006:443-445.( in Chinese)

    [14]楊帆,樊丁,彭凱,等.基于試車數(shù)據(jù)的航空發(fā)動機(jī)起動過程建模[J].空軍工程大學(xué)學(xué)報:自然科學(xué)版,2013(6):1-4.YANG Fan,FAN Ding,PENG Kai,et al.Modeling of an aero-engine in starting process based on the test data[J].Journal of Air Force Engineering University:Natural Science Edition,2013(6):1-4.(in Chinese)

    [15]Sellers J F,Daniele C J.Dyngen:a program for calculating steady-state and transient performance of turbojet and turbofan engines[R].NASA-TN-D-7901,1975.

    [16]Daniele C J,Krosel S M,Szuch J R,et al.Digital computer program for generating dynamic turbofan engine models [R].NASA-TM-83446,1983.

    [17]李慶揚,莫孜中,祁力群.非線性方程組的數(shù)值解法[M].北京:科學(xué)出版社,1999:108-113.LI Qingyang,MO Zizhong,QI Liqun.Numerical solutions of nonlinear equations[M].Beijing:Science Press,1999:108-113.(in Chinese)

    [18]Claus R W,Evans A L,Lylte J K,et al.Numerical propulsion system simulation [J].Computing Systems in Engineering,1991,2(4):357-364.

    [19]Visser W P J,Broomhead M J.GSP,a generic object-oriented gas turbine simulation environment[C]//ASME Turbo Expo 2000:Power for Land,Sea,and Air.2000:V001T01A002.

    [20]Broyden C G.Quasi-Newton methods and their application to function minimization[J].Mathematics of Computation,1967,21(99):368-381.

    猜你喜歡
    慢車起動機(jī)部件
    小慢車“尋鮮記”
    云南畫報(2021年9期)2021-12-02 05:07:00
    2009款配備智能進(jìn)入和起動系統(tǒng)的漢蘭達(dá)車起動機(jī)控制系統(tǒng)簡介
    基于Siemens NX和Sinumerik的銑頭部件再制造
    高鐵時代的“小慢車”
    部件拆分與對外漢字部件教學(xué)
    起動機(jī)進(jìn)口壓力對某型發(fā)動機(jī)地面起動影響分析
    航空浮子式汽化器慢車系統(tǒng)工作原理與典型故障分析
    廣州地鐵14號線“快、慢車”運行方案分析
    水輪機(jī)過流部件改造與節(jié)能增效
    簡述永磁減速起動機(jī)
    汽車電器(2014年8期)2014-02-28 12:14:29
    18禁美女被吸乳视频| 看片在线看免费视频| 最新的欧美精品一区二区| 国产精品国产av在线观看| 一进一出抽搐动态| 男女之事视频高清在线观看| 国产精品 国内视频| 如日韩欧美国产精品一区二区三区| 亚洲av熟女| 久久久久久久精品吃奶| 中文欧美无线码| 超色免费av| 又大又爽又粗| 黄片大片在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | 99热网站在线观看| 热99re8久久精品国产| 国产亚洲欧美精品永久| 丰满人妻熟妇乱又伦精品不卡| 国产在线精品亚洲第一网站| 美女福利国产在线| 国产精品乱码一区二三区的特点 | 久久草成人影院| 亚洲欧洲精品一区二区精品久久久| 国产1区2区3区精品| 成人精品一区二区免费| 新久久久久国产一级毛片| 亚洲精品国产精品久久久不卡| 视频在线观看一区二区三区| 亚洲精品国产色婷婷电影| 九色亚洲精品在线播放| 精品一品国产午夜福利视频| 免费看十八禁软件| 9色porny在线观看| 欧美 日韩 精品 国产| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩综合在线一区二区| 亚洲国产欧美网| 精品电影一区二区在线| 久久精品人人爽人人爽视色| e午夜精品久久久久久久| 免费一级毛片在线播放高清视频 | 国产精品免费视频内射| 1024香蕉在线观看| 国产熟女午夜一区二区三区| 啦啦啦视频在线资源免费观看| 一进一出抽搐gif免费好疼 | 国产亚洲欧美98| 精品国产国语对白av| 亚洲 国产 在线| 久久久水蜜桃国产精品网| 欧美成狂野欧美在线观看| 久久热在线av| 久久久久久久国产电影| 国产av精品麻豆| 日本撒尿小便嘘嘘汇集6| 热re99久久国产66热| 91字幕亚洲| 亚洲精品粉嫩美女一区| 色婷婷av一区二区三区视频| 老司机亚洲免费影院| 99国产极品粉嫩在线观看| 91成年电影在线观看| 18禁黄网站禁片午夜丰满| 啦啦啦在线免费观看视频4| 身体一侧抽搐| 国产免费现黄频在线看| 18禁裸乳无遮挡动漫免费视频| netflix在线观看网站| av线在线观看网站| 91成年电影在线观看| 精品久久久久久电影网| 女性被躁到高潮视频| 国产xxxxx性猛交| 18禁美女被吸乳视频| 成人三级做爰电影| 国产精品亚洲一级av第二区| 五月开心婷婷网| 精品熟女少妇八av免费久了| 黄色毛片三级朝国网站| 日韩成人在线观看一区二区三区| 中国美女看黄片| 高清在线国产一区| 自线自在国产av| 91成年电影在线观看| 亚洲成a人片在线一区二区| 美国免费a级毛片| 村上凉子中文字幕在线| 欧美精品人与动牲交sv欧美| 久久青草综合色| 久久久久久人人人人人| 狠狠婷婷综合久久久久久88av| 亚洲精品乱久久久久久| 国产精品一区二区精品视频观看| 在线免费观看的www视频| 国产色视频综合| 免费日韩欧美在线观看| 国产又色又爽无遮挡免费看| 国产成人影院久久av| av天堂久久9| 丰满饥渴人妻一区二区三| 两个人免费观看高清视频| 视频在线观看一区二区三区| 两个人看的免费小视频| 久久久久久人人人人人| 免费观看人在逋| 人妻久久中文字幕网| 在线国产一区二区在线| 9191精品国产免费久久| 一区二区日韩欧美中文字幕| 亚洲九九香蕉| 国产成人精品无人区| 免费看十八禁软件| 色综合欧美亚洲国产小说| 在线观看舔阴道视频| 99国产精品一区二区三区| 国产男靠女视频免费网站| 欧美激情 高清一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 99精品久久久久人妻精品| 亚洲欧洲精品一区二区精品久久久| 国产成+人综合+亚洲专区| 精品国产亚洲在线| 亚洲美女黄片视频| 美女高潮喷水抽搐中文字幕| 国产亚洲欧美在线一区二区| 亚洲七黄色美女视频| 成人黄色视频免费在线看| 国产成人免费观看mmmm| 久久人人爽av亚洲精品天堂| 精品卡一卡二卡四卡免费| 777久久人妻少妇嫩草av网站| 身体一侧抽搐| 女人高潮潮喷娇喘18禁视频| 精品一品国产午夜福利视频| 一进一出好大好爽视频| 国产亚洲欧美在线一区二区| 最近最新中文字幕大全免费视频| 最新的欧美精品一区二区| 女警被强在线播放| 国产国语露脸激情在线看| 自拍欧美九色日韩亚洲蝌蚪91| 国产片内射在线| 丰满迷人的少妇在线观看| 丰满迷人的少妇在线观看| 美女福利国产在线| 超碰成人久久| 一边摸一边做爽爽视频免费| 美女视频免费永久观看网站| 久久久久久人人人人人| 老熟妇仑乱视频hdxx| 女人爽到高潮嗷嗷叫在线视频| 精品无人区乱码1区二区| 男人的好看免费观看在线视频 | 久久久国产一区二区| 69精品国产乱码久久久| 国产av又大| 久久久国产成人免费| 久久久国产成人免费| 69av精品久久久久久| 国产精品 国内视频| 两人在一起打扑克的视频| 日日爽夜夜爽网站| 国产亚洲av高清不卡| 欧美在线黄色| 91精品三级在线观看| 国产欧美日韩精品亚洲av| 精品国产美女av久久久久小说| 人成视频在线观看免费观看| 亚洲精品粉嫩美女一区| 日韩精品免费视频一区二区三区| 免费黄频网站在线观看国产| 午夜影院日韩av| 亚洲成人国产一区在线观看| 久久婷婷成人综合色麻豆| 超碰成人久久| 下体分泌物呈黄色| 中文字幕精品免费在线观看视频| 免费在线观看视频国产中文字幕亚洲| 欧美日韩瑟瑟在线播放| 麻豆成人av在线观看| 曰老女人黄片| 变态另类成人亚洲欧美熟女 | 欧美在线黄色| 丁香欧美五月| 国产蜜桃级精品一区二区三区 | 精品卡一卡二卡四卡免费| 亚洲 欧美一区二区三区| 99热国产这里只有精品6| 热99re8久久精品国产| 精品无人区乱码1区二区| 三级毛片av免费| 成人国语在线视频| 后天国语完整版免费观看| 夜夜爽天天搞| 最近最新免费中文字幕在线| 电影成人av| 欧美精品一区二区免费开放| 欧美在线黄色| avwww免费| 啦啦啦免费观看视频1| 免费看a级黄色片| 精品福利观看| 一级毛片女人18水好多| 黄片播放在线免费| 久久午夜综合久久蜜桃| 久久国产精品人妻蜜桃| 亚洲在线自拍视频| 9色porny在线观看| 国产99久久九九免费精品| 夜夜躁狠狠躁天天躁| 免费在线观看日本一区| 国产精品电影一区二区三区 | 国产主播在线观看一区二区| 日韩欧美三级三区| 人妻一区二区av| 亚洲精品粉嫩美女一区| 中文字幕人妻丝袜制服| 欧美乱妇无乱码| 精品熟女少妇八av免费久了| 亚洲精品乱久久久久久| 黄色丝袜av网址大全| 精品久久久久久久久久免费视频 | 俄罗斯特黄特色一大片| 1024视频免费在线观看| 一区在线观看完整版| 两人在一起打扑克的视频| 欧美丝袜亚洲另类 | 欧美 日韩 精品 国产| 黄色 视频免费看| 欧美日韩黄片免| 伊人久久大香线蕉亚洲五| 国产精品久久久av美女十八| 少妇被粗大的猛进出69影院| 黄色视频不卡| 日韩视频一区二区在线观看| 国内毛片毛片毛片毛片毛片| 免费人成视频x8x8入口观看| 久久精品亚洲av国产电影网| 亚洲成国产人片在线观看| 丰满迷人的少妇在线观看| 亚洲中文日韩欧美视频| 女性被躁到高潮视频| 激情视频va一区二区三区| 老熟妇仑乱视频hdxx| 精品视频人人做人人爽| 在线天堂中文资源库| 黄频高清免费视频| 香蕉国产在线看| 中文字幕色久视频| av天堂久久9| 中文字幕精品免费在线观看视频| 国产视频一区二区在线看| 亚洲av成人不卡在线观看播放网| av在线播放免费不卡| 最近最新免费中文字幕在线| 高清在线国产一区| 亚洲午夜理论影院| 免费日韩欧美在线观看| 超碰97精品在线观看| 99re在线观看精品视频| 欧美激情极品国产一区二区三区| 亚洲欧美一区二区三区黑人| 中亚洲国语对白在线视频| 中文字幕人妻熟女乱码| svipshipincom国产片| 我的亚洲天堂| a级毛片黄视频| 久久精品国产a三级三级三级| 午夜福利一区二区在线看| 精品乱码久久久久久99久播| 免费高清在线观看日韩| a在线观看视频网站| 18禁裸乳无遮挡动漫免费视频| 成人18禁在线播放| 亚洲午夜理论影院| 美女视频免费永久观看网站| 在线免费观看的www视频| 欧美激情高清一区二区三区| 国产精品美女特级片免费视频播放器 | 亚洲专区字幕在线| 黑人操中国人逼视频| 亚洲国产欧美日韩在线播放| 亚洲少妇的诱惑av| 丰满的人妻完整版| 叶爱在线成人免费视频播放| 免费看十八禁软件| 色综合欧美亚洲国产小说| 91在线观看av| 纯流量卡能插随身wifi吗| 午夜免费观看网址| cao死你这个sao货| 在线观看免费午夜福利视频| 精品第一国产精品| 黄色成人免费大全| 午夜福利欧美成人| 一区在线观看完整版| 黄色片一级片一级黄色片| 国内久久婷婷六月综合欲色啪| 亚洲欧美一区二区三区久久| 亚洲精品在线美女| 日本欧美视频一区| 亚洲av电影在线进入| 国产欧美日韩综合在线一区二区| 最近最新中文字幕大全电影3 | 又大又爽又粗| 少妇的丰满在线观看| 男男h啪啪无遮挡| 12—13女人毛片做爰片一| 国产高清激情床上av| 亚洲九九香蕉| 免费在线观看亚洲国产| 国产激情欧美一区二区| 亚洲欧美精品综合一区二区三区| 自线自在国产av| 男女床上黄色一级片免费看| 久久久久久久久久久久大奶| 欧美黄色淫秽网站| 人人妻人人添人人爽欧美一区卜| 一二三四在线观看免费中文在| 国产在线观看jvid| 国产av精品麻豆| 女警被强在线播放| 91大片在线观看| 老司机亚洲免费影院| 午夜福利视频在线观看免费| 亚洲成av片中文字幕在线观看| 精品福利永久在线观看| x7x7x7水蜜桃| 在线观看免费高清a一片| 欧美 日韩 精品 国产| 精品人妻熟女毛片av久久网站| 亚洲精华国产精华精| 日本一区二区免费在线视频| 美女福利国产在线| 国产在线一区二区三区精| 欧美日本中文国产一区发布| 一边摸一边抽搐一进一小说 | 国产深夜福利视频在线观看| 在线天堂中文资源库| 99国产精品一区二区蜜桃av | 国产野战对白在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 1024香蕉在线观看| 五月开心婷婷网| 美女 人体艺术 gogo| 黄片播放在线免费| 精品国产超薄肉色丝袜足j| 他把我摸到了高潮在线观看| 母亲3免费完整高清在线观看| 国产亚洲欧美98| 在线观看舔阴道视频| 欧美+亚洲+日韩+国产| 五月开心婷婷网| 国产精品美女特级片免费视频播放器 | 又大又爽又粗| 国产男女内射视频| 亚洲国产精品合色在线| 国产一区有黄有色的免费视频| 精品免费久久久久久久清纯 | 久久精品国产a三级三级三级| 1024视频免费在线观看| 在线国产一区二区在线| 超碰97精品在线观看| 岛国毛片在线播放| 三级毛片av免费| 中文字幕高清在线视频| 男人舔女人的私密视频| 国产免费现黄频在线看| 人人妻人人澡人人看| 在线国产一区二区在线| 久久久国产欧美日韩av| 在线看a的网站| 狂野欧美激情性xxxx| 国产精品国产高清国产av | x7x7x7水蜜桃| 国产在视频线精品| 久久精品国产亚洲av香蕉五月 | 9热在线视频观看99| 欧美黑人精品巨大| 两性夫妻黄色片| 热re99久久国产66热| 男女午夜视频在线观看| 满18在线观看网站| www.熟女人妻精品国产| 日本黄色视频三级网站网址 | 法律面前人人平等表现在哪些方面| 777久久人妻少妇嫩草av网站| 一进一出好大好爽视频| 免费观看a级毛片全部| 亚洲成人手机| 十八禁高潮呻吟视频| 男女下面插进去视频免费观看| 可以免费在线观看a视频的电影网站| 91九色精品人成在线观看| 麻豆成人av在线观看| 国产主播在线观看一区二区| 国产精品一区二区在线观看99| a级毛片在线看网站| 久久青草综合色| 国产欧美日韩一区二区三区在线| 国产黄色免费在线视频| 王馨瑶露胸无遮挡在线观看| 亚洲av成人av| 正在播放国产对白刺激| av在线播放免费不卡| 18禁裸乳无遮挡免费网站照片 | 亚洲七黄色美女视频| 美女扒开内裤让男人捅视频| 变态另类成人亚洲欧美熟女 | av超薄肉色丝袜交足视频| 久久热在线av| 纯流量卡能插随身wifi吗| 91av网站免费观看| 老司机福利观看| 动漫黄色视频在线观看| 91成年电影在线观看| 中文字幕色久视频| 国产亚洲精品久久久久久毛片 | 亚洲国产中文字幕在线视频| 国产亚洲精品久久久久久毛片 | 人人澡人人妻人| 嫩草影视91久久| 亚洲精品中文字幕一二三四区| 国产高清videossex| 亚洲精品中文字幕在线视频| 国产成人精品在线电影| 国产亚洲欧美精品永久| 一边摸一边做爽爽视频免费| 亚洲精品国产精品久久久不卡| 夜夜夜夜夜久久久久| 水蜜桃什么品种好| 久久精品人人爽人人爽视色| 老司机福利观看| 999精品在线视频| 国产有黄有色有爽视频| 国产男女内射视频| 久久久久久免费高清国产稀缺| 国产成人精品在线电影| 法律面前人人平等表现在哪些方面| 久久久久久免费高清国产稀缺| 精品人妻熟女毛片av久久网站| 性少妇av在线| 国产精品免费视频内射| 久久中文看片网| 丝袜人妻中文字幕| e午夜精品久久久久久久| bbb黄色大片| 国产xxxxx性猛交| 久久中文字幕一级| 国产高清国产精品国产三级| 飞空精品影院首页| 国产精品98久久久久久宅男小说| 一进一出抽搐gif免费好疼 | 99国产综合亚洲精品| 老司机福利观看| 中文字幕制服av| 69av精品久久久久久| 国产1区2区3区精品| 久久狼人影院| 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看| 久久久久精品人妻al黑| 精品午夜福利视频在线观看一区| 热re99久久精品国产66热6| 男女之事视频高清在线观看| 精品国产乱子伦一区二区三区| 无遮挡黄片免费观看| 亚洲色图 男人天堂 中文字幕| 国产在视频线精品| 亚洲久久久国产精品| 成熟少妇高潮喷水视频| 国产欧美日韩一区二区三区在线| 国产精华一区二区三区| 国产精品一区二区精品视频观看| 亚洲av日韩在线播放| 午夜91福利影院| 另类亚洲欧美激情| 亚洲免费av在线视频| 亚洲精品国产区一区二| 精品免费久久久久久久清纯 | 久久精品aⅴ一区二区三区四区| 日本一区二区免费在线视频| 9热在线视频观看99| 亚洲人成77777在线视频| 欧美国产精品一级二级三级| 男人舔女人的私密视频| 亚洲一区二区三区欧美精品| 亚洲欧美一区二区三区久久| 欧美黄色片欧美黄色片| 国产97色在线日韩免费| 精品一区二区三区视频在线观看免费 | 国产aⅴ精品一区二区三区波| 女人被躁到高潮嗷嗷叫费观| 亚洲五月天丁香| 亚洲欧美一区二区三区久久| 欧美黄色片欧美黄色片| 成熟少妇高潮喷水视频| 久久久久久久国产电影| 亚洲国产欧美一区二区综合| 午夜91福利影院| 久久精品国产清高在天天线| 大码成人一级视频| 嫩草影视91久久| 国产精华一区二区三区| 久久精品亚洲熟妇少妇任你| 啦啦啦 在线观看视频| 伦理电影免费视频| 国产精品久久久av美女十八| 在线观看日韩欧美| 国产精品98久久久久久宅男小说| 午夜福利免费观看在线| 搡老乐熟女国产| 国产真人三级小视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲成人免费电影在线观看| 最近最新中文字幕大全免费视频| 国产成人精品久久二区二区免费| 天堂动漫精品| 成人18禁在线播放| 无遮挡黄片免费观看| 操出白浆在线播放| 一区二区三区激情视频| 亚洲aⅴ乱码一区二区在线播放 | 自线自在国产av| 成人手机av| 久久 成人 亚洲| 欧美av亚洲av综合av国产av| 一进一出抽搐动态| 欧美乱码精品一区二区三区| 亚洲精品国产色婷婷电影| 国产野战对白在线观看| 无人区码免费观看不卡| 国产一区二区三区视频了| 亚洲一区二区三区欧美精品| 女性被躁到高潮视频| 99国产综合亚洲精品| 黄色毛片三级朝国网站| 国产精品二区激情视频| 免费久久久久久久精品成人欧美视频| 国产野战对白在线观看| 国产精品秋霞免费鲁丝片| 少妇粗大呻吟视频| 十八禁人妻一区二区| 久热爱精品视频在线9| 欧美精品人与动牲交sv欧美| 欧美日韩瑟瑟在线播放| 国产精品av久久久久免费| 亚洲成人手机| 亚洲熟妇中文字幕五十中出 | 国产单亲对白刺激| 国产精品影院久久| 美女 人体艺术 gogo| 午夜视频精品福利| 别揉我奶头~嗯~啊~动态视频| 在线观看www视频免费| 中国美女看黄片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品一区二区三区在线| 香蕉久久夜色| 精品国产一区二区三区久久久樱花| 啦啦啦视频在线资源免费观看| 捣出白浆h1v1| 久久青草综合色| 亚洲av成人av| 久久青草综合色| 日本精品一区二区三区蜜桃| 久久热在线av| 日韩欧美一区二区三区在线观看 | 国产免费男女视频| 欧美日韩乱码在线| 亚洲一区二区三区不卡视频| 成在线人永久免费视频| 色播在线永久视频| 国产午夜精品久久久久久| 精品少妇久久久久久888优播| 男女床上黄色一级片免费看| 黑人欧美特级aaaaaa片| 亚洲男人天堂网一区| 黄片播放在线免费| 中文字幕人妻丝袜一区二区| 免费在线观看日本一区| 国产欧美亚洲国产| 真人做人爱边吃奶动态| 成人免费观看视频高清| 亚洲av片天天在线观看| 国产精品一区二区在线观看99| 精品卡一卡二卡四卡免费| 久久国产精品大桥未久av| 咕卡用的链子| 性色av乱码一区二区三区2| 在线十欧美十亚洲十日本专区| 精品免费久久久久久久清纯 | 欧美日韩中文字幕国产精品一区二区三区 | 啦啦啦在线免费观看视频4| 日韩三级视频一区二区三区| 免费观看精品视频网站| 国产欧美亚洲国产| 午夜日韩欧美国产| 极品教师在线免费播放| 中文字幕精品免费在线观看视频| 男女免费视频国产| 777久久人妻少妇嫩草av网站| 久久精品成人免费网站| 亚洲五月色婷婷综合| 一区二区三区精品91| 黄色 视频免费看| 99香蕉大伊视频| 天天躁夜夜躁狠狠躁躁| e午夜精品久久久久久久| 国产成人系列免费观看| 国产精品免费大片| 91麻豆精品激情在线观看国产 | 国产高清国产精品国产三级| 日本精品一区二区三区蜜桃|