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

        ?

        基于非線性濾波環(huán)節(jié)的新型魯棒PID控制策略的研究

        2013-09-21 11:00:54萬(wàn)文軍
        動(dòng)力工程學(xué)報(bào) 2013年2期
        關(guān)鍵詞:頻率特性開(kāi)環(huán)微分

        李 軍, 萬(wàn)文軍, 張 曦

        (廣東電網(wǎng)公司電力科學(xué)研究院,廣州510080)

        魯棒特性較好的PID控制策略在工業(yè)過(guò)程控制領(lǐng)域中占主導(dǎo)地位,為提高熱工對(duì)象的控制品質(zhì),在調(diào)節(jié)回路中加入微分環(huán)節(jié)是工程實(shí)踐中常采用的典型方法之一.然而,微分環(huán)節(jié)在改善系統(tǒng)動(dòng)態(tài)性能的同時(shí),其高通濾波特性卻放大了高頻干擾,往往使微分環(huán)節(jié)的超前校正作用難以實(shí)現(xiàn),甚至在有外部高頻干擾時(shí)容易造成系統(tǒng)調(diào)節(jié)機(jī)構(gòu)大幅波動(dòng)[1].而且,隨著過(guò)程特性及操作條件的改變,PID控制系統(tǒng)難以繼續(xù)保持良好的控制品質(zhì).因此,魯棒PID控制器的設(shè)計(jì)受到了廣泛的關(guān)注.針對(duì)多模型不確定對(duì)象,Ge等[2]提出一種基于LQR-LMI的魯棒PID控制器設(shè)計(jì)方法;Goncalves等[3]提出一種魯棒兩自由度PID控制器設(shè)計(jì)方法,以滿足閉環(huán)系統(tǒng)的H2/H∞性能指標(biāo);Kim等[4]提出一種基于擴(kuò)展Lagrange粒子群算法的魯棒PID控制器設(shè)計(jì)方法,來(lái)滿足多個(gè)H∞性能指標(biāo);劉建民等[5-8]也都提出了各具特色的魯棒PID控制器設(shè)計(jì)方法.但這些方法均存在算法復(fù)雜,在DCS中較難實(shí)現(xiàn)的問(wèn)題.

        針對(duì)以上問(wèn)題,將非線性濾波環(huán)節(jié)與傳統(tǒng)的PID控制策略相結(jié)合[1],提出了一種基于非線性濾波環(huán)節(jié)的新型魯棒PID控制策略,該控制策略算法簡(jiǎn)單,易于實(shí)現(xiàn),提高了傳統(tǒng)PID控制策略的魯棒性,且控制回路能獲得優(yōu)良的控制品質(zhì).

        1 非線性濾波環(huán)節(jié)的構(gòu)造與分析

        控制系統(tǒng)的頻域特性分析方法是控制系統(tǒng)設(shè)計(jì)、性能分析、參數(shù)調(diào)整的基本方法[9-11],該方法的特點(diǎn)是易于理解.因此,筆者主要從頻率特性角度進(jìn)行計(jì)算與分析.

        1.1 非線性濾波環(huán)節(jié)構(gòu)造圖

        圖1為非線性濾波環(huán)節(jié)的構(gòu)造圖,非線性濾波環(huán)節(jié)由二階微分環(huán)節(jié)非線性壓縮模塊、預(yù)失真模塊、非線性信號(hào)重整模塊、一階慣性模塊等構(gòu)成.各模塊的作用如下:

        (1)二階微分環(huán)節(jié)非線性壓縮模塊:采用平方根特性對(duì)二階微分環(huán)節(jié)輸出信號(hào)進(jìn)行非線性壓縮.

        (2)預(yù)失真模塊[12]:對(duì)輸入信號(hào)進(jìn)行非線性預(yù)處理,按預(yù)定的(平方根特性)規(guī)律產(chǎn)生失真.

        (3)非線性信號(hào)重整模塊:對(duì)二階微分環(huán)節(jié)非線性壓縮模塊、預(yù)失真模塊輸出信號(hào)進(jìn)行重整,其作用為進(jìn)行非線性補(bǔ)償,對(duì)于正弦波或余弦波函數(shù),理論上非線性補(bǔ)償后的失真為零.另外,該模塊還可以保證輸出信號(hào)與輸入信號(hào)之間成比例變化.

        (4)一階慣性模塊:對(duì)非線性信號(hào)重整模塊的輸出進(jìn)行一階慣性濾波.

        圖1 非線性濾波環(huán)節(jié)的構(gòu)造Fig.1 Structure of the non-linear filter link

        1.2 非線性濾波環(huán)節(jié)頻域特性分析

        圖1中二階微分環(huán)節(jié)所對(duì)應(yīng)的頻域函數(shù)表達(dá)式如下:

        令s=ωj,得到相應(yīng)的頻域函數(shù)表達(dá)式為:

        對(duì)上述二階微分環(huán)節(jié)信號(hào)進(jìn)行帶符號(hào)開(kāi)根號(hào)運(yùn)算,其帶符號(hào)開(kāi)根號(hào)的定義為:

        為不失一般性,假設(shè)回路的輸入信號(hào)是幅值為A,角頻率為ω的正弦波信號(hào),即輸入信號(hào)為X(t)=Asin(ωt).結(jié)合式(2)和式(3),該正弦波信號(hào)的二階微分值帶符號(hào)開(kāi)方等式為:

        對(duì)過(guò)程輸入信號(hào)進(jìn)行絕對(duì)值開(kāi)根號(hào)運(yùn)算,獲得非線性預(yù)失真信號(hào),得到相應(yīng)的過(guò)程函數(shù)表達(dá)式為:

        將二階微分信號(hào)的非線性壓縮信號(hào)與非線性預(yù)失真信號(hào)相乘,獲得非線性重整信號(hào),得到相應(yīng)的過(guò)程函數(shù)表達(dá)式為:經(jīng)過(guò)非線性重整后的信號(hào)與輸入信號(hào)傳遞的頻域函數(shù)為:

        式(7)也可表達(dá)為:

        式(9)也可表達(dá)為:

        將式(9)的頻域函數(shù)進(jìn)行分母有理化得到式(11)

        為方便分析,令增益Kd=1,則實(shí)數(shù)項(xiàng)為虛數(shù)項(xiàng)為實(shí)數(shù)項(xiàng)和虛數(shù)項(xiàng)對(duì)應(yīng)的幅值-頻率特性數(shù)學(xué)計(jì)算結(jié)果見(jiàn)圖2(其中Kd=1、ωo=1/Td).

        從圖2可以看出,實(shí)數(shù)項(xiàng)(圖2中實(shí)線)對(duì)應(yīng)的幅頻特性具有典型的帶通濾波特性,等價(jià)為一帶通濾波器,負(fù)值表示反相位.虛數(shù)項(xiàng)(圖2中虛線)對(duì)應(yīng)的幅頻特性具有典型的高通濾波特性.實(shí)數(shù)項(xiàng)和虛數(shù)項(xiàng)的這種特性與實(shí)際微分環(huán)節(jié)的幅值特性存在本質(zhì)區(qū)別.

        圖2 非線性環(huán)節(jié)幅頻特性分析Fig.2 Frequency-domain characteristic analysis of the non-linear filter link

        2 非線性魯棒PID控制策略

        為了驗(yàn)證上述方法的有效性,進(jìn)行了大量仿真試驗(yàn).非線性濾波環(huán)節(jié)結(jié)合傳統(tǒng)PID控制策略的仿真試驗(yàn)方框圖見(jiàn)圖3,其中取非線性濾波環(huán)節(jié)的參數(shù)與微分環(huán)節(jié)相同,仿真對(duì)象為6階慣性等容對(duì)象.對(duì)應(yīng)的開(kāi)環(huán)傳遞函數(shù)為:

        閉環(huán)傳遞函數(shù)為:

        2.1 非線性濾波環(huán)節(jié)提高系統(tǒng)魯棒性的定性分析

        圖3 魯棒PID控制策略仿真框圖Fig.3 Block diagram of proposed robust PID control strategy for simulation

        經(jīng)典控制理論的Nyquist穩(wěn)定性判據(jù)是通過(guò)確定性開(kāi)環(huán)系統(tǒng)P(s)的頻率特性進(jìn)行的,控制系統(tǒng)的魯棒性能分析則是通過(guò)不確定開(kāi)環(huán)系統(tǒng)P(s)的頻率特性進(jìn)行的.理論上,系統(tǒng)穩(wěn)定性是由開(kāi)環(huán)系統(tǒng)穩(wěn)定性裕量決定的.而熱工控制系統(tǒng)穩(wěn)定性裕量又主要由系統(tǒng)在中頻段的特性決定,中頻段是指系統(tǒng)增益為0dB附近的一段頻率.在對(duì)象確定以后,系統(tǒng)穩(wěn)定性裕量便由PID控制器中頻段輸出超前虛頻幅值與實(shí)頻幅值特性決定.

        理論與實(shí)踐經(jīng)驗(yàn)表明,對(duì)于高階被控對(duì)象,相應(yīng)的PID調(diào)節(jié)器回路的響應(yīng)速度與穩(wěn)定性之間存在矛盾.而非線性濾波環(huán)節(jié)提高系統(tǒng)魯棒性的原理在于:非線性濾波環(huán)節(jié)既能夠較大幅度地降低PID控制器中頻段輸出實(shí)頻幅值,又能夠較大幅度地提高PID控制器高中頻段輸出超前虛頻幅值,在較大幅度提高PID調(diào)節(jié)器回路響應(yīng)速度的同時(shí)又具有較高的穩(wěn)定性,較好地解決了PID調(diào)節(jié)器回路的響應(yīng)速度與穩(wěn)定性之間的矛盾.

        2.2 非線性濾波環(huán)節(jié)提高穩(wěn)定性裕量的定量分析

        對(duì)于圖3給出的仿真圖,具體參數(shù)見(jiàn)3.1節(jié),其開(kāi)環(huán)調(diào)節(jié)器對(duì)應(yīng)的實(shí)頻和超前虛頻輸出特性曲線見(jiàn)圖4和圖5.選擇非線性濾波環(huán)節(jié)參數(shù)的原則是使實(shí)數(shù)項(xiàng)的峰值出現(xiàn)在系統(tǒng)頻率特性的中頻段范圍,這樣就能夠較大幅度地降低調(diào)節(jié)器在中頻段輸出實(shí)頻幅值,具體見(jiàn)圖4.非線性濾波環(huán)節(jié)中的虛數(shù)項(xiàng)能夠較大幅度地提高調(diào)節(jié)器在中高頻段輸出超前虛頻幅值,具體見(jiàn)圖5所示(正值為超前、負(fù)值為滯后).

        圖4 2種控制系統(tǒng)開(kāi)環(huán)調(diào)節(jié)器輸出實(shí)頻對(duì)比圖Fig.4 Output comparison of open loop system between traditional PID and robust PID(R)

        圖6為2種控制系統(tǒng)開(kāi)環(huán)頻率特性對(duì)比圖.由圖6可知,在加入非線性濾波環(huán)節(jié)后,相位穩(wěn)定裕量從70°略微下降到69.5°,可認(rèn)為幾乎沒(méi)有變化.然而,幅值穩(wěn)定裕量從6dB提高到10.2dB,可見(jiàn)非線性濾波環(huán)節(jié)可以明顯提高開(kāi)環(huán)系統(tǒng)穩(wěn)定性裕量.

        圖5 2種控制系統(tǒng)開(kāi)環(huán)調(diào)節(jié)器輸出虛頻對(duì)比圖Fig.5 Output comparison of open loop system between traditional PID and robust PID (I)

        圖6 6階慣性等容仿真對(duì)象的2種控制系統(tǒng)開(kāi)環(huán)頻率特性對(duì)比圖Fig.6 Comparison of open-loop frequency characteristics between two control systems by 6-order inertial isovolumetric simulation

        由圖6還可以看出,在相位穩(wěn)定裕量判定點(diǎn)頻率和幅值穩(wěn)定裕量判定點(diǎn)頻率范圍內(nèi),系統(tǒng)開(kāi)環(huán)頻率特性幅頻增益顯著下降,相頻相位滯后值顯著下降,即開(kāi)環(huán)對(duì)象滯后特性得到了改善.此外,在加入非線性濾波環(huán)節(jié)后,2個(gè)穩(wěn)定裕量判定點(diǎn)對(duì)應(yīng)的頻率寬度B2比未加入非線性濾波環(huán)節(jié)前2個(gè)穩(wěn)定裕量判定點(diǎn)對(duì)應(yīng)的頻率寬度明顯增加,這意味著魯棒性提高了.原因在于:在臨界穩(wěn)定狀態(tài),2個(gè)穩(wěn)定裕量判定點(diǎn)頻率相等,其對(duì)應(yīng)的頻率寬度為零.而在臨界穩(wěn)定邊界以內(nèi),2個(gè)穩(wěn)定裕量判定點(diǎn)之間的頻率寬度越大,其適應(yīng)不確定開(kāi)環(huán)系統(tǒng)P(s)頻率特性變化的范圍也越大.當(dāng)圖3中仿真慣性等容對(duì)象的階數(shù)從6階增加到12階時(shí),得到的2種控制系統(tǒng)開(kāi)環(huán)頻率特性的對(duì)比見(jiàn)圖7.

        圖7 12階慣性等容仿真對(duì)象的2種控制系統(tǒng)開(kāi)環(huán)頻率特性對(duì)比圖Fig.7 Comparison of open-loop frequency characteristics between two control systems by 12-order inertial isovolumetric simulation

        由圖7可知,未加入非線性濾波環(huán)節(jié)時(shí)的相位穩(wěn)定裕量和幅值穩(wěn)定裕量幾乎為零,可認(rèn)為已處于臨界穩(wěn)定狀態(tài),而加入非線性濾波環(huán)節(jié)后的相位穩(wěn)定裕量為14°、幅值穩(wěn)定裕量為1.6dB,仍處于臨界穩(wěn)定邊界內(nèi).

        非線性濾波環(huán)節(jié)必須與微分環(huán)節(jié)聯(lián)合使用,單獨(dú)使用非線性濾波環(huán)節(jié)的效果并不理想.例如對(duì)某些具體的實(shí)際過(guò)程系統(tǒng),如果通過(guò)簡(jiǎn)單提高系統(tǒng)響應(yīng)速度就能較好地解決問(wèn)題,則應(yīng)用非線性濾波環(huán)節(jié)可起到良好的效果.對(duì)于系統(tǒng)來(lái)說(shuō),通過(guò)微分作用可提高系統(tǒng)響應(yīng)速度,通過(guò)非線性濾波環(huán)節(jié)可充分保證系統(tǒng)穩(wěn)定性,將非線性濾波器與傳統(tǒng)PID聯(lián)合使用的控制方法稱為新型魯棒PID控制器.下面的仿真將充分說(shuō)明該P(yáng)ID控制回路的特性.

        3 仿真研究

        3.1 階躍響應(yīng)仿真

        選取仿真對(duì)象為6階慣性等容對(duì)象,對(duì)應(yīng)的時(shí)間常數(shù)為Ta=20s.依據(jù)傳統(tǒng)PID參數(shù)整定原則得到控制器參數(shù)為:Kp=0.85、Ti=95s、Kd=0.75、Td=33s,并取非線性濾波環(huán)節(jié)中的參數(shù)等于PID中的微分參數(shù)Kd=0.75、Td=33s.

        將提出的魯棒PID控制策略與傳統(tǒng)PID及PI控制策略進(jìn)行了對(duì)比分析,系統(tǒng)的階躍響應(yīng)曲線見(jiàn)圖8.從圖8可以看出,魯棒PID控制策略調(diào)節(jié)效果良好,與傳統(tǒng)方法相比具有更大的穩(wěn)定裕量.

        由圖8還可以看出,采用非線性濾波環(huán)節(jié)的PID控制策略后,其調(diào)節(jié)上升速率和超調(diào)量方面的品質(zhì)均優(yōu)于傳統(tǒng)PID控制策略.另外,在提高響應(yīng)速度后,采用非線性濾波環(huán)節(jié)可以充分保證調(diào)節(jié)回路的穩(wěn)定性.可見(jiàn)加入非線性濾波環(huán)節(jié)后,允許PID參數(shù)的調(diào)整范圍更寬,一定程度上增強(qiáng)了控制回路的魯棒性.

        圖8 3種控制策略階躍響應(yīng)仿真結(jié)果的對(duì)比Fig.8 Comparison of simulation results among robust PID,PI and traditional PID

        3.2 400s方波響應(yīng)仿真

        為了進(jìn)一步說(shuō)明所提出方法的有效性,對(duì)系統(tǒng)進(jìn)行了400s方波響應(yīng)測(cè)試,得到的響應(yīng)曲線見(jiàn)圖9.從圖9可以看出,所提出的控制策略抗擾動(dòng)能力較強(qiáng),具有更為優(yōu)良的控制特性.

        圖9 2種控制策略400s方波響應(yīng)仿真結(jié)果的對(duì)比Fig.9 Comparison of 400ssquare wave response between robust and traditional PID

        4 工程應(yīng)用

        將提出的非線性濾波環(huán)節(jié)應(yīng)用于某660MW機(jī)組的汽包鍋爐過(guò)熱汽溫調(diào)節(jié)系統(tǒng)中,該鍋爐過(guò)熱汽溫調(diào)節(jié)系統(tǒng)僅有一級(jí)調(diào)節(jié).優(yōu)化前,鍋爐過(guò)熱汽溫調(diào)節(jié)系統(tǒng)采用的是常規(guī)PID控制策略,系統(tǒng)的過(guò)程量在變負(fù)荷過(guò)程中經(jīng)常發(fā)生較大幅度的波動(dòng).為此,在原有的主控回路中引入本文提出的非線性濾波環(huán)節(jié),并且較大幅度地提高主控PID的比率增益(提高30%)和微分增益(提高40%),系統(tǒng)仍然很穩(wěn)定.投運(yùn)效果表明:引入非線性濾波環(huán)節(jié)后,系統(tǒng)調(diào)節(jié)品質(zhì)得到了明顯提高,變負(fù)荷過(guò)程過(guò)熱汽溫動(dòng)態(tài)波動(dòng)范圍有較大幅度降低,其中B側(cè)過(guò)熱汽溫調(diào)節(jié)系統(tǒng)優(yōu)化前后對(duì)比如圖10所示.

        圖10 B側(cè)過(guò)熱汽溫調(diào)節(jié)系統(tǒng)優(yōu)化前后過(guò)熱汽溫變化趨勢(shì)圖Fig.10 Temperature variation of B-side superheated steam before and after optimization

        從圖10可以看出,優(yōu)化前,汽溫變化速率較慢,趨勢(shì)雖然比較平滑,但汽溫上下波動(dòng)幅度較大,其中的一個(gè)原因在于PID調(diào)節(jié)回路響應(yīng)速度較慢.曾經(jīng)嘗試提高PID調(diào)節(jié)回路的響應(yīng)速度,但調(diào)節(jié)回路的穩(wěn)定性顯著變差.采用本文提出的非線性濾波環(huán)節(jié)進(jìn)行優(yōu)化后,在較大幅度提高PID調(diào)節(jié)回路響應(yīng)速度的同時(shí)又具有較高的穩(wěn)定性,對(duì)抑制變負(fù)荷過(guò)程中汽溫波動(dòng)的效果顯著.通過(guò)對(duì)比分析實(shí)際應(yīng)用效果,進(jìn)一步證實(shí)了本文提出方法的優(yōu)越性.

        5 結(jié)束語(yǔ)

        提出了一種基于非線性濾波環(huán)節(jié)的新型魯棒PID控制策略,并對(duì)非線性環(huán)節(jié)的頻域特性進(jìn)行了深入分析,給出了濾波環(huán)節(jié)提高系統(tǒng)魯棒性的原理.通過(guò)將非線性濾波環(huán)節(jié)與傳統(tǒng)PID控制相結(jié)合,較大程度提高了傳統(tǒng)PID控制策略的魯棒性,獲得了良好的控制品質(zhì),在提高系統(tǒng)響應(yīng)速度的同時(shí)也提高了系統(tǒng)的穩(wěn)定性.階躍擾動(dòng)、方波擾動(dòng)仿真試驗(yàn)以及對(duì)象特性變動(dòng)后魯棒性分析及實(shí)際應(yīng)用效果驗(yàn)證了該策略的有效性.

        [1]李軍.基于改進(jìn)微分環(huán)節(jié)的控制策略的研究及其應(yīng)用[J].熱力發(fā)電,2010,39(8):72-75.LI Jun.Study on new type control strategy based on improved differential link and application thereof[J].Thermal Power Generation,2010,39(8):72-75.

        [2]GE M,CHIU M S,WANG Q G.Robust PID controller design via LMI approach[J].Journal of Process Control,2002,12(1):3-13.

        [3]GONCALVES E N,PALHARES R M,TAKAHASHI R H C.A novel approach for H2/H∞r(nóng)obust PID synthesis for uncertain systems[J].Journal of Process Control,2008,18(1):19-26.

        [4]KIM T H,MARUTA I,SUGIE T.Robust PID controller tuning based on the constrained particle swarm optimization[J].Automatica,2008,44(4):1104-1110.

        [5]韓愷,趙均,朱豫才,等.一種擾動(dòng)自適應(yīng)的魯棒預(yù)測(cè)控制算法[J].化工學(xué)報(bào),2009,60(7):1730-1738.HAN Kai,ZHAO Jun,ZHU Yucai,et al.A robust MPC technique with adaptive disturbance model[J].CIESC Journal,2009,60(7):1730-1738.

        [6]勒其兵,葉瓊瑜,趙大力.一種改進(jìn)的魯棒分散PID控制器的設(shè)計(jì)[J].系統(tǒng)仿真學(xué)報(bào),2009,21(3):780-783.JIN Qibing,YE Qiongyu,ZHAO Dali.Modified robust decentralized PID controllers design[J].Journal of System Simulation,2009,21(3):780-783.

        [7]王傳峰,李東海,姜學(xué)智,等.基于概率魯棒性的鍋爐過(guò)熱汽溫串級(jí)PID控制器[J].清華大學(xué)學(xué)報(bào):自然科學(xué)版,2009,49(2):249-252.WANG Chuanfeng,LI Donghai,JIANG Xuezhi,et al.Cascade PID controllers for superheated steam boiler temperatures based on probabilistic robustness[J].J Tsinghua Univ:Science and Technology,2009,49(2):249-252.

        [8]劉建民,韓璞,開(kāi)平安,等.一種改進(jìn)的α-β-γ濾波器和PID控制器[J].動(dòng)力工程,2007,27(4):551-554.LIU Jianmin,HAN Pu,KAI Pingan,et al.An improvedα-β-γfilter and PID controller[J].Journal of Power Engineering,2007,27(4):551-554.

        [9]THAM M T.Why frequency response[R].Australia:University of Newcastle upon Tyne,1999:1-17.

        [10]WANG Yagang,CAI Wenjian,GE Ming.Decentralized relay-based multivariable process identification in the frequency domain[J].IEEE Trans On Automatic Control,2003,48(5):873-877.

        [11]李軍,萬(wàn)文軍,張曦.一種基于階躍響應(yīng)的理想頻率信號(hào)源及頻域分析的研究[J].動(dòng)力工程學(xué)報(bào),2012,32(4):308-314.LI Jun,WAN Wenjun,ZHANG Xi.Study on ideal frequency signal source based on step response and its utility in frequency domain analysis[J].Journal of Chinese Society of Power Engineering,2012,32(4):308-314.

        [12]張玉興,趙宏飛,向榮,等.非線性電路與系統(tǒng)[M].北京:機(jī)械工業(yè)出版社,2007:183-185.

        猜你喜歡
        頻率特性開(kāi)環(huán)微分
        擬微分算子在Hp(ω)上的有界性
        上下解反向的脈沖微分包含解的存在性
        轉(zhuǎn)速開(kāi)環(huán)恒壓頻比的交流調(diào)速系統(tǒng)的分析與仿真研究
        考慮頻率特性的變頻負(fù)荷模型研究
        借助微分探求連續(xù)函數(shù)的極值點(diǎn)
        DAM中波發(fā)射機(jī)電聲指標(biāo)頻率特性初探
        一種溫和環(huán)醚開(kāi)環(huán)成雙酯的新方法
        基于開(kāi)環(huán)補(bǔ)償?shù)娘w機(jī)偏航角控制系統(tǒng)設(shè)計(jì)及仿真
        對(duì)不定積分湊微分解法的再認(rèn)識(shí)
        基于AD9854的簡(jiǎn)易頻率特性測(cè)試儀的設(shè)計(jì)
        久久久久亚洲AV片无码乐播 | 久久精品国产av大片| 91精品亚洲熟妇少妇| 精品精品国产高清a毛片| 亚洲av无码av男人的天堂| 亚洲免费av电影一区二区三区| 人妻精品人妻一区二区三区四五| 中文字幕人妻少妇伦伦| 精品免费久久久久久久| h国产视频| 人妻在线中文字幕视频| 日本伊人精品一区二区三区| 一本一道av无码中文字幕﹣百度 | 日韩精品人妻视频一区二区三区| 无码人妻丰满熟妇啪啪网不卡 | 国产成人一区二区三中文| 成人免费毛片在线播放| 中文字幕一区二区人妻秘书| 玩弄放荡人妻少妇系列视频| 国产 中文 制服丝袜 另类| 91国产视频自拍在线观看| 国产精品成人观看视频国产奇米| 欧美猛男军警gay自慰| 国产精品白浆视频一区| 国家一级内射高清视频| 亚洲av永久无码天堂网小说区| 免费无码av片在线观看网址| 丰满人妻中文字幕乱码| 美女人妻中出日本人妻| 肥臀熟女一区二区三区| 国产精品无码无片在线观看3D| 久久精品国产亚洲av沈先生| 国产一区二区三区经典| 亚洲视频在线观看一区二区三区| 人妻体内射精一区二区三区| 又白又嫩毛又多15p| 久久精品国产屋| 可以直接在线看国产在线片网址| 香港三级日本三级a视频| 欧美性福利| 国产精品久久一区二区蜜桃|