亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于空氣流量預測的發(fā)動機空燃比三步非線性控制

        2017-12-04 09:27:54林佳眉
        農(nóng)業(yè)機械學報 2017年11期
        關(guān)鍵詞:發(fā)動機模型

        王 萍 林佳眉 陳 虹

        (1.吉林大學汽車仿真與控制國家重點實驗室, 長春 130025; 2.吉林大學控制科學與工程系, 長春 130025)

        基于空氣流量預測的發(fā)動機空燃比三步非線性控制

        王 萍1,2林佳眉2陳 虹1,2

        (1.吉林大學汽車仿真與控制國家重點實驗室, 長春 130025; 2.吉林大學控制科學與工程系, 長春 130025)

        為降低排放并保證發(fā)動機的動力性,空燃比(Air-fuel ratio, AFR)控制成為人們?nèi)找骊P(guān)注的焦點??杖急鹊乃矐B(tài)精確控制可同時兼顧扭矩輸出和尾氣排放。為達到瞬態(tài)空燃比的精確控制,提出了空燃比三步非線性控制方法,包含穩(wěn)態(tài)控制、參考前饋控制和誤差反饋控制。當空燃比期望值為恒值時,穩(wěn)態(tài)控制起到主導作用;當車輛扭矩需求發(fā)生速變時,空燃比期望值也發(fā)生跳變,參考前饋控制此時將及時做出反應;誤差反饋控制以廢氣氧傳感器(EGO)測得的空燃比作為反饋量進行修正??紤]噴油器執(zhí)行機構(gòu)的延遲,基于模型實時預測未來進氣歧管的壓力,通過前饋控制進行延遲補償。利用仿真平臺en-DYNA中的四缸發(fā)動機模型進行仿真驗證,證明算法瞬態(tài)工況和參考輸入快變時的有效性。

        發(fā)動機; 空燃比; 三步非線性控制; 均值模型; 空氣流量預測

        引言

        為了使汽油機在不損失動力性的前提下,盡量降低燃油的消耗,減少有害物的排放,空燃比對扭矩控制和排放控制具有重要的意義[1]。汽油發(fā)動機的空燃比控制是根據(jù)進入氣缸的空氣量來調(diào)節(jié)噴油器的噴油量,從而滿足空氣和汽油的化學當量比,使氣缸中的混合氣體能完全燃燒。空燃比大于理論值的混合氣稱為稀混合氣,氣多油少,燃燒完全,油耗低,污染小,但功率較小。空燃比小于理論值的混合氣稱為濃混合氣,氣少油多,功率較大,但燃燒不完全,油耗高,污染大。在工程上,根據(jù)車輛運行工況的變化,空燃比期望值不單純設置成理論空燃比,而是在理論空燃比附近一個區(qū)域內(nèi)變化[2-4]。當汽油機工作在瞬態(tài)工況時,由于發(fā)動機轉(zhuǎn)速的瞬態(tài)變化對進氣量影響很大,而氧傳感器反饋的缸內(nèi)空燃比具有較大傳輸時滯,油路的噴油器也存在執(zhí)行延遲,這些因素都為瞬態(tài)空燃比精確控制帶來很大困難。

        目前工程上常采用查MAP表的方式來控制空燃比[5]。穩(wěn)態(tài)工況下,發(fā)動機空燃比的控制可以通過查取MAP表作為前饋并結(jié)合反饋PID控制來實現(xiàn)。瞬態(tài)工況下,采用基于查MAP表的開環(huán)控制策略。然而由于發(fā)動機較強的非線性,MAP表存在誤差,導致控制精度低[6]。文獻[7-8]提出了基于神經(jīng)網(wǎng)絡的模型預測控制,實現(xiàn)了多步預測,提高了前饋的控制精度,但在瞬態(tài)工況下,仍存在較大超調(diào)。為了提高空燃比控制精度,需要從氣路和油路2個方面來進行研究。文獻[9]采用非線性卡爾曼濾波器對瞬態(tài)工況的進氣量進行估計,文獻[10]利用輸入估計算法估計氣缸的進氣量,降低了空燃比的跟蹤誤差。燃油濕壁效應和油路部分的執(zhí)行延遲也是需要重點考慮的因素[11-13]??紤]油路傳輸?shù)臐癖谛?,文獻[14-15]提出了基于擴張卡爾曼濾波的自適應估計器和自適應模糊滑膜控制器來觀測油膜動態(tài)特性,建立濕壁效應和油膜蒸發(fā)的模型,進一步提高了空燃比的控制精度。文獻[16]中使用最小二乘的方法預測進氣量,補償了油路部分的時滯。文獻[17]中采用魯棒控制策略控制空燃比,考慮了不確定因素的影響,實現(xiàn)了瞬態(tài)工況下空燃比控制,但設計過程相對復雜,工程上不易實現(xiàn)。

        文獻[18]提出了一種易于工程實現(xiàn)的三步非線性控制器設計方法,并應用于汽油缸內(nèi)直噴發(fā)動機的軌壓控制,它采用工程上常見的前饋加反饋的結(jié)構(gòu),設計過程簡單,控制器參數(shù)物理意義明確。考慮空燃比控制的前饋加反饋經(jīng)典結(jié)構(gòu),本文提出空燃比三步非線性控制方法,兼顧穩(wěn)態(tài)工況、瞬態(tài)工況和空燃比期望值的變化。另外,針對油路執(zhí)行延遲問題,通過預測進氣歧管壓力來預測氣缸進氣量,在穩(wěn)態(tài)控制器中對延遲進行補償。

        1 系統(tǒng)建模

        發(fā)動機空燃比系統(tǒng)分為3部分:進氣歧管子模型、氣缸進氣量子模型和燃油傳輸子模型。 系統(tǒng)模型如圖1所示,均采用均值模型[19-20]。

        圖1 發(fā)動機系統(tǒng)模型Fig.1 System model of engine1.節(jié)氣門 2.進氣歧管 3.進氣閥 4.火花塞 5.排氣閥 6.排氣歧管 7.活塞 8.曲軸

        1.1 進氣歧管子模型

        (1)

        根據(jù)理想氣體方程pmVm=maRTm(R為氣體常數(shù),Tm(K)和Vm(m3)為進氣歧管的溫度和體積), 對進氣歧管壓力pm(Pa)求導, 可得

        (2)

        (3)

        1.2 氣缸進氣量子模型

        (4)

        式中Vc——發(fā)動機排量,m3

        η——發(fā)動機容積效率

        ω——根據(jù)負載扭矩變化的發(fā)動機轉(zhuǎn)速,rad/s

        (5)

        1.3 燃油傳輸子模型

        在汽油機中,噴油器噴出的燃油由2部分組成,一部分直接形成燃油蒸汽進入氣缸,另一部分則附著在進氣管管壁上,形成油膜。考慮到濕壁效應的影響,基于微分方程,建立了燃油傳輸子模型

        (6)

        式中uf——噴油器噴出的燃油質(zhì)量,kg

        mfc——進入氣缸的燃油質(zhì)量,kg

        mfw——油膜部分的燃油質(zhì)量,kg

        τw——油膜揮發(fā)時間常數(shù)

        εw——燃油沉積系數(shù)

        根據(jù)空燃比λ的定義,λ表示為

        (7)

        (8)

        2 三步非線性控制器設計

        空燃比三步法非線性控制系統(tǒng)如圖2所示。圖中θ為曲軸轉(zhuǎn)角,Te為發(fā)動機指示扭矩。

        圖2 空燃比三步法非線性控制系統(tǒng)結(jié)構(gòu)示意圖Fig.2 Block diagram of AFR control

        “三步法”控制策略包含3個階段:穩(wěn)態(tài)控制、參考前饋控制和誤差反饋控制?;凇叭椒ā钡目刂瓶驁D如圖3所示。根據(jù)前面建立的系統(tǒng)模型,控制目標為空燃比能跟蹤期望值λs。

        圖3 基于“三步法”的控制框圖Fig.3 Block diagram of triple-step control scheme

        2.1 穩(wěn)態(tài)控制

        首先設計穩(wěn)態(tài)控制器,根據(jù)式(8),對控制器的控制輸出λ求導

        (9)

        (10)

        (11)

        將式(11)整理為

        (12)

        其中

        (13)

        (14)

        2.2 參考前饋控制

        考慮發(fā)動機的扭矩需求,實際發(fā)動機空燃比期望值是在理想空燃比附近區(qū)域變化的,為了能夠及時對期望值的變化做出反應,在穩(wěn)態(tài)控制的基礎上,加入?yún)⒖记梆伩刂苪ff。假設此時的控制律為

        u=us+uff

        (15)

        (16)

        2.3 誤差反饋控制

        由于外部干擾和建模誤差的存在,穩(wěn)態(tài)控制和參考前饋控制也并不能保證輸出達到目標期望值,需要加入誤差反饋控制ue。此時的控制律為

        u=us+uff+ue

        (17)

        將式(17)代入到式(12)中,得到

        (18)

        定義跟蹤誤差為e=λs-λ,對誤差e求導

        (19)

        將式(18)整理為

        (20)

        定義李雅普諾夫函數(shù)V為

        (21)

        其中

        式中k1——待調(diào)整參數(shù),k1gt;0

        對李雅普諾夫函數(shù)V求導

        (22)

        (23)

        式中k0——待調(diào)整參數(shù),k0gt;0

        綜合式(14)、(16)、(23),得到最終的虛擬控制律為

        (24)

        3 進氣歧管壓力預測

        在實際的發(fā)動機控制系統(tǒng)中,從發(fā)動機噴油指令的更新到實際噴出燃油的執(zhí)行存在一定的時間延時。根據(jù)式(8),期望進入氣缸的燃油量可以近似表示為

        (25)

        從式(25)中可以看到,油路部分中進入氣缸的燃油量和進入氣缸的空氣量有著直接關(guān)系。假定噴油指令的延時是2個采樣周期,為了補償油路部分的執(zhí)行延遲,則需要對進入氣缸的空氣量進行2步預測。根據(jù)式(5),對進入氣缸的進氣量的預測可以轉(zhuǎn)換為對進氣歧管壓力的預測。此時的空燃比控制系統(tǒng)結(jié)構(gòu)示意圖如圖4所示。

        圖4 帶有進氣歧管壓力預測的空燃比控制系統(tǒng)示意圖Fig.4 Block diagram of AFR control with intake manifold prediction

        本文研究的是四缸發(fā)動機,將式(2)在曲軸轉(zhuǎn)角域內(nèi)每180°進行離散化,采樣周期T為

        (26)

        (27)

        其中

        因為要預測2個采樣周期的進氣歧管壓力, 將式(27)中的pm(k)用pm(k+1)代替,則pm(k+2)可以表示為

        (28)

        在實際工況中, 參數(shù)a1、a2會根據(jù)不同工況而變化,需要通過最小二乘法來辨識這2個參數(shù),辨識值可用1、2來表示。因此, 兩步的進氣歧管壓力預測可整理為

        (29)

        (30)

        (31)

        (32)

        為了得到最終的控制律,將式(24)離散化為

        (33)

        u=us+uff+ue=

        (34)

        4 仿真驗證

        在TESIS公司開發(fā)的enDYNA虛擬發(fā)動機軟件中對所設計的控制系統(tǒng)進行仿真驗證。為驗證在發(fā)動機瞬態(tài)工況下,基于“三步法”設計的空燃比控制器的控制效果,選取四缸汽油機模型。在仿真結(jié)果中,期望空燃比的變化曲線如圖5所示,對期望空燃比做歸一化處理,當發(fā)動機轉(zhuǎn)速或節(jié)氣門開度較大時,期望空燃比為1.05;當發(fā)動機轉(zhuǎn)速或節(jié)氣門開度較小時,期望空燃比為0.95;其他工況下期望空燃比為1。

        圖5 期望空燃比變化曲線Fig.5 Changing curve of reference AFR

        4.1 進氣歧管壓力預測

        進氣歧管壓力預測結(jié)果如圖6所示。由圖6可以看到,經(jīng)過預測后的進氣歧管壓力比原有的進氣歧管壓力超前一些。

        圖6 進氣歧管壓力預測結(jié)果Fig.6 Prediction of intake manifold pressure

        4.2 變發(fā)動機轉(zhuǎn)速

        為了更直觀地表現(xiàn)出“三步法”控制策略的控制效果,采取相同的工況,與文獻[21]中的空燃比反饋調(diào)節(jié)控制方法做了對比,并模擬了2種發(fā)動機工況。

        在變發(fā)動機轉(zhuǎn)速工況下,設定節(jié)氣門開度維持在30°不變,發(fā)動機轉(zhuǎn)速在1 200~2 000 r/min之間變化,模擬加速和減速的過程,如圖7所示??杖急雀櫱€、進入氣缸進氣量曲線以及燃油噴射量曲線如圖8~10所示。

        圖7 轉(zhuǎn)速變化曲線Fig.7 Changing curve of engine speed

        圖8 空燃比跟蹤曲線Fig.8 Tracking curves of AFR

        圖9 氣缸進氣量估計Fig.9 Estimation of air mass flow into cylinder

        圖10 噴油嘴的噴油量Fig.10 Fuel injection

        從圖8中可以看到,“三步法”控制器的穩(wěn)態(tài)誤差比反饋調(diào)節(jié)器的穩(wěn)態(tài)誤差要小,可以將空燃比控制在期望空燃比附近一個很窄的范圍內(nèi),同時在轉(zhuǎn)速產(chǎn)生變化時,空燃比的超調(diào)量也比較小。

        4.3 變節(jié)氣門開度

        在變節(jié)氣門開度工況下, 設定發(fā)動機轉(zhuǎn)速維持1 200 r/min不變,節(jié)氣門開度在15°~30°之間變化,如圖11所示??杖急雀櫱€、進入氣缸進氣量曲線以及燃油噴射量曲線圖12~14所示。

        圖11 節(jié)氣門開度變化曲線Fig.11 Changing curve of throttle angle

        由圖12可以看出,2個控制器都能將空燃比控制在期望空燃比附近,但當節(jié)氣門開度突變時,普通反饋調(diào)節(jié)器的超調(diào)更加明顯,“三步法”控制器能夠使空燃比經(jīng)過一個很小的超調(diào)后,快速回到穩(wěn)態(tài)。由此可見,基于“三步法”設計的控制器不僅設計步驟簡單,而且控制效果良好。

        圖12 變節(jié)氣門開度空燃比跟蹤曲線Fig.12 Tracking curves of AFR

        圖13 變節(jié)氣門開度氣缸進氣量估計Fig.13 Estimation of air mass flow into cylinder

        圖14 變節(jié)氣門開度噴油嘴的噴油量Fig.14 Fuel injection

        5 結(jié)束語

        本文首先針對油路部分的延遲,對進氣歧管壓力作了2步預測,從而使進入氣缸進氣量的估計得到改善,補償了油路部分的時滯。其次,為了更好地跟蹤瞬態(tài)工況下的空燃比,基于“三步法”控制方法設計了空燃比控制器。在仿真平臺enDYNA四缸發(fā)動機模型下,驗證了控制器的跟蹤性能,并與普通反饋調(diào)節(jié)控制器作了對比。結(jié)果表明“三步法”控制器的穩(wěn)態(tài)誤差和瞬態(tài)誤差都比較小,響應迅速,跟蹤效果好。

        1 劉麗萍,王祝煒,張振東,等. 電控汽油機滿足未來排放和經(jīng)濟性要求的研究[J]. 農(nóng)業(yè)機械學報,2000,31(1):119-122.

        LIU Liping, WANG Zhuwei, ZHANG Zhendong, et al. Study on required economy and mession of futrue electronic-controlled gasoline engines[J]. Transactions of the Chinese Society for Agricultural Machinery,2000,31(1):119-122.(in Chinese)

        2 PAK. Model predictive engine air-ratio control using online sequential extreme learning machine[J]. Neural Computing and Applications, 2016, 27(1): 79-92.

        3 GILLES T. Automotive service: inspection, maintenance, repair [M]. Cengage Learning, 2012.

        4 ZMUDKA Z, POSTRZEDINK S. Inverse aspects of the three-way catalytic converter operation in the spark ignition engine[J]. Journal of KONES, 2011, 18: 509-516.

        5 萬里平,蔣炎坤,郭巖,等. 基于油耗和空燃比的汽油機充氣效率測定方法[J/OL]. 農(nóng)業(yè)機械學報,2014,45(3):32-36,25. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1amp;file_no=20140306amp;journal_id=jcsam. DOI:10.6041/j.issn.1000-1298.2014.03.006.

        WAN Liping, JIANG Yankun, GUO Yan, et al. Gasolin engine charge-efficiency measuring based on fuel consumption and AFR test[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014,45(3):32-36,25. (in Chinese)

        6 劉曉亮,王生昌,李茂月. 汽油機空燃比控制和時間延遲的研究[J]. 公路與汽運, 2008(2):11-14.

        LIU Xiaoliang, WANG Shengchang, LI Maoyue. Study on air fuel ratio control and time delay of gasoline engine[J]. Highways and Automotive Applications, 2008(2):11-14. (in Chinese)

        7 WANG S W, YU D L, GOMM J B, et al. Adaptive neural network model based predictive control for air-fuel ratio of SI engines[J]. Engineering Applications of Artificial Intelligence, 2006, 19(2): 189-200.

        8 王莉,劉德新,李萬眾. 稀燃汽油機瞬態(tài)空燃比的滑模-神經(jīng)網(wǎng)絡控制[J]. 農(nóng)業(yè)機械學報,2007,38(10):1-4.

        WANG Li, LIU Dexin, LI Wanzhong. Sliding mode-neural network control for transient air-fuel ratio of lean burn gasoline engine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2007,38(10):1-4. (in Chinese)

        9 馮煜,焦曉紅. 汽油機非線性卡爾曼濾波暫態(tài)進氣量估計及空燃比控制[J]. 控制理論與應用, 2015,32(4):546-553.

        FENG Yu, JIAO Xiaohong. Transient cylinder air charge estimation with nonlinear Kalman filter for air-fuel ratio control of gasoline engines[J]. Control Theory and Applications, 2015,32(4):546-553. (in Chinese)

        10 ALEXANDER S. Application of input estimation techniques to charge estimation and control in automotive engines [J]. Control Engineering Practice, 2002, 10(12): 1371-1383.

        11 洪木南,李建秋,歐陽明高. 汽油機油膜模型參數(shù)辨識[J]. 農(nóng)業(yè)機械學報,2010,41(11):26-30.

        HONG Munan, LI Jianqiu, OUYANG Minggao. Parameters identification of the fuel film model for spark ignition engines[J]. Transactions of the Chinese Society for Agricultural Machinery,2010,41(11):26-30.(in Chinese)

        12 KHAJOMTRAIDET C, ITO K. Simple adaptive air-fuel ratio control of a port injection SI engine with a cylinder pressure sensor [J]. Journal of Control Theory and Applications, 2015, 13(2): 141-150.

        13 ZHANG J,SHEN T,XU G,et al. Wall-wetting model based method for air-fuel ratio transient control in gasoline engines with dual injection system [J]. International Journal of Automotive Technology, 2013, 14(6): 867-873.

        14 ARSIE I, PIANESE C, RIZZO G, et al. An adaptive estimator of fuel film dynamics in the intake port of a spark ignition engine[J]. Control Engineering Practice, 2003, 11(3): 303-309.

        15 王立標,李岳林. 瞬態(tài)空燃比的自適應模糊滑模補償控制[J]. 臺州學院學報,2016,38(3):23-29.

        WANG Libiao,LI Yuelin. Adaptive fuzzy sliding mode compensation control for engine transient air fuel ratio[J]. Journal of Taizhou University, 2016,38(3):23-29. (in Chinese)

        16 趙亞楠. 發(fā)動機暫態(tài)進氣量估計研究及實驗驗證[D].秦皇島:燕山大學, 2013.

        ZHAO Yanan. Air charge estimation under transient mode for SI engines and its experimental verification[D]. Qinhuangdao:Yanshan University, 2013. (in Chinese)

        17 石向南,鄧元望,朱浩. 電噴發(fā)動機過渡工況空燃比魯棒控制研究[J]. 湖南大學學報:自然科學版,2010,37(9):19-23.

        SHI Xiangnan,DENG Yuanwang, ZHU Hao. Research on the transient air fuel ratio control of SI engine based on robust control theory[J]. Journal of Hunan University:Natural Sciences, 2010,37(9):19-23. (in Chinese)

        18 CHEN Hong, GONG Xun, LIU Qifang, et al. Triple-step method to design non-linear controller for rail pressure of gasoline direct injection engines[J]. Control Theory and Applications Iet, 2014, 8(11): 948-959.

        19 HENDRICKS E, SORENSON S C. Mean value modelling of spark ignition engines[C]. SAE Technical Paper,1990.

        20 ELBERT H,VESTERHOLM H.The analysis of mean value SI engine models[C]. SAE Technical Paper,1992.

        21 JIAO Xiaohong, SHEN Tielong. Lyapunov-design of adaptive air-fuel ratio control for gasoline engines based on mean-value model[C]∥Control Conference (CCC), 2011 30th Chinese. IEEE, 2011: 6146-6150.

        AFRTriple-stepControlofGasolineEngineBasedonAirMassFlowPrediction

        WANG Ping1,2LIN Jiamei2CHEN Hong1,2

        (1.StateKeyLaboratoryofAutomotiveSimulationandControl,JilinUniversity,Changchun130025,China2.DepartmentofControlScienceandEngineering,JilinUniversity,Changchun130025,China)

        Control of the air-fuel ratio (AFR) in gasoline engines is of imminent importance when aiming at minimizing calibration effort and meeting performance requirements. People have higher demands on the gasoline engine, which has less exhaust emission, better economic efficiency and favorable engine power performance. In order to keep the air-fuel ratio close to the stoichiometric value under transient conditions, AFR precise control was achieved by employing triple-step method which was easily to be implemented in engineering. The structure of the designed controller consisted of three parts: steady-state control, feed-forward control concerning the reference variations and error feedback control. When the desired AFR was a constant, the steady-state control would play a dominant role. And the feed-forward control would react immediately when the desired AFR was changed on account of the torque requirements changed. The feedback control would amendment AFR measured by the exhaust gas oxygen (EGO) sensor which can enhance the close-loop performance and rearranged into a state-dependent PI. A straightforward design process was provided, and the structure of the designed non-linear controller was easily achieved, which was comparable to those widely used in current automotive control. Furthermore, taking the implement delay of the injection into consideration, it can be compensated by feed-forward control based on predicting the intake manifold pressure. Finally, the simulation results in the environment of en-DYNA with a reasonable common four-cylinder engine model showed the efficiency of the proposed method. And the predicted intake manifold pressure was visibly advanced to that without prediction one in the simulation result.

        engine; air-fuel ratio; triple-step method; mean-value model; air mass flow prediction

        10.6041/j.issn.1000-1298.2017.11.049

        TP273

        A

        1000-1298(2017)11-0398-07

        2017-06-21

        2017-08-16

        國家自然科學基金項目(61403159)和吉林省教育廳“十三五”科學研究規(guī)劃項目(2016431)

        王萍(1982—),女,副教授,主要從事預測控制和發(fā)動機控制研究,E-mail: wangping12@jlu.edu.cn

        陳虹(1963—),女,教授,博士生導師,主要從事先進控制和汽車電子控制研究,E-mail: chenh@jlu.edu.cn

        猜你喜歡
        發(fā)動機模型
        一半模型
        2015款寶馬525Li行駛中發(fā)動機熄火
        重要模型『一線三等角』
        2012年奔馳S600發(fā)動機故障燈偶爾點亮
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        3D打印中的模型分割與打包
        FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
        新一代MTU2000發(fā)動機系列
        發(fā)動機的怠速停止技術(shù)i-stop
        新型1.5L-Eco-Boost發(fā)動機
        日韩啪啪精品一区二区亚洲av| 中国老熟女重囗味hdxx| 99久久久无码国产精品试看| 一本大道久久a久久综合| 最新日本免费一区二区三区| 国产一区二区视频免费在线观看| 日本黄网站三级三级三级| 狠狠色狠狠色综合| 手机av男人天堂免费网址| 男人的天堂手机版av| 最近2019年好看中文字幕视频| chinese国产乱在线观看| 女同成片av免费观看| 亚洲色图专区在线视频| 蜜臀av午夜一区二区三区| 久久久精品免费观看国产 | 国产精品99久久精品女同| 亚洲av日韩一区二区| 日本精品αv中文字幕| 国产成人乱色伦区小说| 国产免费一区二区三区三| 日韩经典午夜福利发布| 亚洲av永久无码天堂网手机版| 白白色免费视频一区二区| 大香蕉青青草视频在线| 台湾佬中文娱乐网22| 亚洲中文字幕第一页在线| 亚洲福利第一页在线观看| 日韩人妻精品中文字幕专区| 一本一道波多野结衣av中文| 国产精品丝袜在线不卡| 午夜视频手机在线免费观看| 日韩欧美在线综合网另类 | 精品亚洲麻豆1区2区3区| 欧美孕妇xxxx做受欧美88| 精品久久久久久国产潘金莲| 日本一区二区在线免费看| 欧美成人片在线观看| 精品无码久久久九九九AV| 国产精品后入内射日本在线观看 | 熟女少妇内射日韩亚洲|