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

    基于干擾觀測器的撓性衛(wèi)星姿態(tài)滑模變結(jié)構(gòu)控制*

    2017-05-03 01:04:34王雪冰
    關(guān)鍵詞:帆板撓性姿態(tài)控制

    王雪冰,吳 忠

    (北京航空航天大學(xué)儀器科學(xué)與光電工程學(xué)院,北京 100191)

    ?

    基于干擾觀測器的撓性衛(wèi)星姿態(tài)滑模變結(jié)構(gòu)控制*

    王雪冰,吳 忠

    (北京航空航天大學(xué)儀器科學(xué)與光電工程學(xué)院,北京 100191)

    為抑制各種內(nèi)外干擾因素對撓性衛(wèi)星姿態(tài)控制性能的影響,設(shè)計基于干擾觀測器的滑模變結(jié)構(gòu)控制器.該控制器采用干擾觀測器對系統(tǒng)中存在的內(nèi)外干擾進(jìn)行估計,并對估計值加以前饋補(bǔ)償.在此基礎(chǔ)上,采用滑模變結(jié)構(gòu)控制器對未補(bǔ)償?shù)母蓴_進(jìn)一步抑制,實(shí)現(xiàn)衛(wèi)星姿態(tài)與姿態(tài)角速度的漸近收斂.與單純的滑模變結(jié)構(gòu)控制器相比,本文的控制器已根據(jù)干擾估計值對干擾進(jìn)行了前饋補(bǔ)償,采用較小的切換增益即可抑制剩余干擾,顫振現(xiàn)象減弱.仿真結(jié)果表明,本文設(shè)計的控制器是可行的.

    撓性航天器;姿態(tài)控制;干擾觀測器;滑模變結(jié)構(gòu)控制

    0 引 言

    為滿足未來航天任務(wù)的復(fù)雜需求,撓性衛(wèi)星的姿態(tài)必須具有足夠高的精度和穩(wěn)定度[1].然而,在撓性衛(wèi)星姿態(tài)控制系統(tǒng)中,存在多種內(nèi)外干擾,如太陽帆板撓性振動、太陽帆板驅(qū)動機(jī)構(gòu)(SADA)干擾和環(huán)境干擾等,嚴(yán)重影響了其控制性能.因此,抗干擾姿態(tài)控制器設(shè)計一直是高精度高穩(wěn)定度衛(wèi)星所要解決的核心關(guān)鍵問題之一.

    針對撓性衛(wèi)星姿態(tài)控制問題,文獻(xiàn)[2]設(shè)計了一種H∞魯棒控制器,并以ETS-VI衛(wèi)星為對象,進(jìn)行了驗證,但其控制器階數(shù)偏高,且性能保守.文獻(xiàn)[3]設(shè)計了基于姿態(tài)誤差四元數(shù)的自適應(yīng)姿態(tài)跟蹤控制律,保證了衛(wèi)星的姿態(tài)跟蹤性能,但該方法處理快變參數(shù)和未建模動態(tài)的能力較弱.文獻(xiàn)[4]針對撓性附件不確定干擾,設(shè)計了內(nèi)??刂破?,提高了控制精度,但是只考慮了正弦周期性干擾.文獻(xiàn)[5]提出了一種雙閉環(huán)自抗擾姿態(tài)控制器,具有較強(qiáng)的魯棒性和干擾抑制能力,但其調(diào)試參數(shù)較多.滑模變結(jié)構(gòu)控制具有設(shè)計方法簡單、響應(yīng)速度快、魯棒性強(qiáng)的特點(diǎn),如文獻(xiàn)[6-9]均研究了撓性航天器的滑模變結(jié)構(gòu)控制方案,但顫振現(xiàn)象是其固有缺陷.作為有效的干擾抑制策略,干擾觀測器可以估計干擾并通過控制器進(jìn)行反饋補(bǔ)償,文獻(xiàn)[10]設(shè)計了一種基于干擾觀測器的前饋補(bǔ)償控制系統(tǒng),姿態(tài)精度達(dá)到了要求,但存在一定靜差.

    本文綜合滑模變結(jié)構(gòu)控制和干擾觀測器的優(yōu)點(diǎn),提出了一種基于干擾觀測器的滑模變結(jié)構(gòu)控制器,通過干擾觀測器有效估計干擾并進(jìn)行前饋補(bǔ)償,滑模變結(jié)構(gòu)控制器對未補(bǔ)償?shù)母蓴_進(jìn)一步抑制,既提高了姿態(tài)控制性能,又能有效減弱顫振現(xiàn)象.

    1 動力學(xué)模型及問題描述

    1.1 衛(wèi)星姿態(tài)動力學(xué)模型

    參照文獻(xiàn)[11-12],考慮帶有單翼旋轉(zhuǎn)撓性太陽帆板的受擾衛(wèi)星,其動力學(xué)方程可表示為

    (1)

    式中,J為撓性衛(wèi)星系統(tǒng)轉(zhuǎn)動慣量,ω為衛(wèi)星本體相對慣性系的角速度在本體系中的投影,Rs,sa為帆板驅(qū)動與星體驅(qū)動的耦合系數(shù)陣,l=[0 1 0]T為太陽帆板的驅(qū)動軸方向,ωsa為帆板轉(zhuǎn)動角速度,F(xiàn)s,sa為帆板振動與星體轉(zhuǎn)動的耦合系數(shù)陣,u為衛(wèi)星的控制力矩,η(t)=[η1(t)η2(t) …ηN(t)]T為太陽帆板的模態(tài)坐標(biāo)矢量,N為模態(tài)階數(shù),ud為衛(wèi)星受到的環(huán)境干擾力矩矢量,Isa為帆板關(guān)于鉸鏈點(diǎn)的轉(zhuǎn)動慣量,F(xiàn)sa為帆板轉(zhuǎn)動與帆板振動的耦合系數(shù)陣,Tsa為帆板的驅(qū)動力矩,Td為太陽帆板驅(qū)動機(jī)構(gòu)的擾動力矩,x×表示x的3個分量構(gòu)成的斜對稱矩陣,矩陣C和K分別為阻尼和剛度矩陣,具體形式為

    (2)

    其中,ξ1,ξ2,…,ξN為阻尼比,Λ1,Λ2,…,ΛN為模態(tài)頻率.

    如設(shè)ωb為衛(wèi)星本體系相對于軌道系的角速度,那么ω與ωb之間的關(guān)系可以表示為

    ω=ωb+ωo

    (3)

    (4)

    1.2 干擾力矩分析

    在式(1)表示的動力學(xué)模型中,衛(wèi)星不僅受到環(huán)境干擾力矩ud的影響,還受到SADA運(yùn)動帶來干擾的影響,下面將逐一分析.

    (1)環(huán)境干擾力矩

    衛(wèi)星受到的空間環(huán)境干擾力矩主要包括:氣動力矩up、重力梯度力矩ug、地磁力矩um和太陽光壓力矩us[13].那么衛(wèi)星受到總的環(huán)境干擾力矩可表示為:

    ud=up+ug+um+us

    (5)

    需要說明的是,環(huán)境干擾力矩對衛(wèi)星姿態(tài)的影響與其運(yùn)動的軌道高度、結(jié)構(gòu)分布情況和空間環(huán)境條件等有關(guān).

    (2)SADA干擾力矩

    設(shè)SADA采用步進(jìn)電機(jī)驅(qū)動,其諧波力矩及導(dǎo)電環(huán)的摩擦力矩是其干擾力矩的主要因素[14-17].

    電機(jī)諧波力矩模型可以表示為[17]:

    Tcog=Tksin(4Ntθsa)

    (6)

    其中,Tk為四階齒槽力矩幅值,Nt為步進(jìn)電機(jī)轉(zhuǎn)子齒數(shù),θsa為太陽帆板的轉(zhuǎn)角.

    導(dǎo)電環(huán)的摩擦力矩模型可以表示為[18]:

    Tf=σ2ωsa+TCsign(ωsa)+

    (TS-TC)e-(ωsa/ωS)2sign(ωsa)

    (7)

    其中,σ2是粘性摩擦阻尼系數(shù),TC為庫倫摩擦力矩,TS為最大靜摩擦力矩,ωS為Stribeck轉(zhuǎn)速.這樣,SADA干擾力矩可表示為:

    Td=Tcog+Tf

    (8)

    (3)轉(zhuǎn)動慣量不確定性

    帆板轉(zhuǎn)動將引起整星轉(zhuǎn)動慣量變化,設(shè)J0為衛(wèi)星轉(zhuǎn)動慣量中的確定部分,ΔJ為變化部分,那么整星轉(zhuǎn)動慣量可表示為:

    J=J0+ΔJ

    (9)

    將其代入式(1)中的第一式,可得:

    (10)

    1.3 問題描述

    將式(3)代入式(10),可得

    (11)

    定義等價干擾d為:

    (12)

    則式(11)可重寫為

    (13)

    由此可見,如果把環(huán)境干擾、帆板驅(qū)動機(jī)構(gòu)干擾、撓性附件振動等干擾因素統(tǒng)一視為等價干擾d,式(1)的動力學(xué)模型可以降階為式(13)的形式.同時,對于由式(4)、式(13)描述的衛(wèi)星系統(tǒng),作如下假設(shè):

    (14)

    假設(shè)2.角速度ω和四元數(shù)[q0q]T均可直接測量或間接計算得到.

    那么,撓性衛(wèi)星姿態(tài)控制問題可以描述為,針對滿足假設(shè)1、假設(shè)2的衛(wèi)星系統(tǒng)式(4)、式(13),考慮系統(tǒng)中存在的干擾,設(shè)計控制力矩u,使得

    (15)

    2 基于干擾觀測器的滑模變結(jié)構(gòu)控制器設(shè)計

    2.1 干擾觀測器設(shè)計

    為抑制干擾d對姿態(tài)控制性能的影響,可對其進(jìn)行前饋補(bǔ)償.為此,需要獲得干擾d的估計信息.根據(jù)文獻(xiàn)[19-20],設(shè)計干擾觀測器如下:

    (16)

    (17)

    將式(17)代入式(16)的第二式,并整理可得

    (18)

    根據(jù)式(13),式(18)可重寫為

    (19)

    (20)

    對式(19)求解,可得:

    (21)

    (22)

    根據(jù)假設(shè)1及式(22),可得:

    (23)

    可見,在式(16)的干擾觀測器下,干擾估計誤差是有界的.

    2.2 控制器設(shè)計

    在干擾觀測器(16)的基礎(chǔ)上,可以對滑模變結(jié)構(gòu)控制器進(jìn)行設(shè)計.定義滑動模態(tài)為

    s=ωb+cq

    (24)

    其中,c>0為常數(shù).

    選取Lyapunov函數(shù)為

    (25)

    對式(25)兩端求導(dǎo),并將式(24)、式(13)和式(4)依次代入,整理可得

    (26)

    設(shè)計控制器u為

    (27)

    其中,k1、k2為大于0的常數(shù),且

    在式(27)的控制器下,式(26)可變?yōu)?/p>

    (28)

    (29)

    此時,Lyapunov函數(shù)的一階導(dǎo)數(shù)式(26)將變?yōu)?/p>

    (30)

    3 仿真分析

    為驗證本文所設(shè)計控制器的有效性,本部分以某帶有單翼旋轉(zhuǎn)撓性太陽帆板的衛(wèi)星為例,進(jìn)行了仿真,并對滑模變結(jié)構(gòu)控制器和基于干擾觀測器的滑模變結(jié)構(gòu)控制器的性能進(jìn)行了對比.在仿真中,采用的衛(wèi)星參數(shù)如下[21]:

    根據(jù)以上參數(shù),可計算得到:

    考慮帆板的前五階模態(tài),其模態(tài)頻率為:[Λ1,Λ2,Λ3,Λ4,Λ5]=[1.084 6,4.351 4,6.739 6,11.514 2,15.756 7]rad/s,阻尼比為:[ξ1,ξ2,ξ3,ξ4,ξ5]= [0.004,0.005,0.006 4,0.008,0.008 5],初始模態(tài)坐標(biāo)設(shè)為0.

    設(shè)與SADA相關(guān)的干擾力矩參數(shù)為:Tk=0.06N·m,Nt=300,σ2=46N·m·s/rad,TC=0.01N·m,TS=0.02N·m,ωS=0.03(°)/s.設(shè)環(huán)境干擾力矩為:ud=10-3×[7.8cos(nt),1.6sin(nt)-2.6cos(nt),-1.2sin(nt)]TN·m,n=0.060(°)/s.

    設(shè)衛(wèi)星三軸姿態(tài)的初值為[1°,1°,-1°],角速度初值為0(°)/s,并采用表1中的控制器參數(shù),對衛(wèi)星姿態(tài)控制系統(tǒng)進(jìn)行了仿真,仿真結(jié)果如圖1~圖11所示.其中,SMC表示滑模變結(jié)構(gòu)控制器,DO+SMC表示基于干擾觀測器的滑模變結(jié)構(gòu)控制器.

    圖1~3為衛(wèi)星三軸姿態(tài)角變化曲線,可看出,進(jìn)入穩(wěn)態(tài)后,與SMC相比,基于干擾觀測器的滑模變結(jié)構(gòu)控制器使衛(wèi)星姿態(tài)控制精度和穩(wěn)定度明顯提高,三軸姿態(tài)誤差均小于2.1×10-6(°),而SMC最大為3.4×10-5(°).圖4~6為衛(wèi)星三軸角速度曲線,可以看出,DO+SMC的角速度誤差小于4.8×10-6(°)/s,而SMC角速度誤差最大為6.5×10-5(°)/s.

    表1 控制器參數(shù)Tab.1 Parameters of the controller

    圖1 滾轉(zhuǎn)角變化曲線Fig.1 Time responses of roll angle

    圖3 偏航角變化曲線Fig.3 Time responses of yaw angle

    圖4 滾轉(zhuǎn)角速度變化曲線Fig.4 Time responses of roll angular velocity

    圖5 俯仰角速度變化曲線Fig.5 Time responses of pitch angular velocity

    圖6 偏航角速度變化曲線Fig.6 Time responses of yaw angular velocity

    圖7 DO+SMC滑模面變化曲線Fig.7 Time responses of sliding mode surface by DO+SMC

    圖7~8表示滑模面變化,對比看出,DO+SMC滑模面收斂速度更快、精度更高,顫振明顯減弱.

    圖9~11分別為衛(wèi)星三軸干擾力矩、干擾力矩估計和估計誤差曲線.可以看出,所設(shè)計的干擾觀測器能對環(huán)境干擾、帆板驅(qū)動機(jī)構(gòu)干擾、撓性附件振動等干擾因素進(jìn)行有效估計.其中滾轉(zhuǎn)軸干擾估計誤差小于實(shí)際值的2.8%,俯仰軸干擾估計誤差小于實(shí)際值的2.6%,偏航軸干擾估計誤差小于實(shí)際值的1.5%.

    圖8 SMC滑模面變化曲線Fig.8 Time responses of sliding mode surface by SMC

    圖9 干擾力矩變化曲線Fig.9 Time responses of disturbance

    圖10 干擾力矩估計曲線Fig.10 Time responses of disturbance estimation

    圖11 干擾力矩估計誤差曲線Fig.11 Time responses of disturbance estimation error

    4 結(jié) 論

    本文針對帶撓性太陽帆板的衛(wèi)星姿態(tài)控制系統(tǒng),設(shè)計了基于干擾觀測器的滑模變結(jié)構(gòu)控制器.干擾觀測器用于估計系統(tǒng)中存在的內(nèi)外干擾,并進(jìn)行前饋補(bǔ)償,滑模變結(jié)構(gòu)控制器進(jìn)一步抑制干擾,既減弱了顫振現(xiàn)象,又能有效實(shí)現(xiàn)衛(wèi)星姿態(tài)與角速度的漸進(jìn)收斂.仿真結(jié)果表明,所設(shè)計的控制系統(tǒng)有效提高了撓性衛(wèi)星姿態(tài)控制性能.

    [1] WANG Z, WU Z. Nonlinear attitude control scheme with disturbance observer for flexible spacecrafts[J]. Nonlinear Dynamics, 2015, 81(1-2): 257-264.

    [2] KIDA T, YAMAGUCHI I, CHIDA Y, et al. On-orbit robust control experiment of flexible spacecraft ETS-VI[J]. Journal of Guidance, Control and Dynamics. 1997, 20(5): 865-872.

    [3] 宗紅, 張洪華, 徐福祥. 撓性衛(wèi)星自適應(yīng)姿態(tài)跟蹤控制[J]. 中國空間科學(xué)技術(shù), 2007, 27(6): 25-30. ZONG H, ZHANG H H, XU F X. Adaptive attitude tracking control of flexible satellite[J]. Chinese Space Science and Technology, 2007, 27(6): 25-30.

    [4] 張洪華, 張國峰. 帶有撓性附件衛(wèi)星的自適應(yīng)內(nèi)??刂芠J]. 空間控制技術(shù)與應(yīng)用, 2001,27(3): 10-16. ZHANG H H, ZHANG G F. Adaptive internal model control for satellite with flexible appendages[J]. Aerospace Control and Application, 2001(3): 10-16.

    [5] 朱承元, 楊滌, 楊旭. 大撓性多體衛(wèi)星的自抗擾姿態(tài)控制系統(tǒng)設(shè)計[J]. 航天控制, 2004, 22(6): 25-31. ZHU C Y, YANG D, YANG X. Design of active disturbance rejection attitude control system for large flexible multi-body satellite[J]. Aerospace Control, 2004, 22(6): 25-31.

    [6] 李新國, 陳士櫓. 撓性衛(wèi)星太陽帆板的變結(jié)構(gòu)主動控制[J]. 西北工業(yè)大學(xué)學(xué)報, 1997, 15(1): 84-87. LI X G, CHEN S L. Variable structure active control of flexible satellites with solar paddles[J]. Journal of Northwestern Polytechnical University, 1997, 15(1): 84-87.

    [7] HU Q, MA G. Variable structure control and active vibration suppression of flexible spacecraft during attitude maneuver[J]. Aerospace Science and Technology, 2005, 9(4): 307-317.

    [8] 秦浩. 帶柔性太陽帆板的航天器姿態(tài)控制實(shí)驗研究[D]. 北京: 北京理工大學(xué), 2015. QIN H. Experimental study on the attitude control of spacecraft with flexible solar arrays[D]. Beijing: Beijing Institute of Technology, 2015.

    [9] 呂建婷, 李傳江, 馬廣富. 衛(wèi)星姿態(tài)調(diào)節(jié)的滑模PID控制器設(shè)計[J]. 哈爾濱工業(yè)大學(xué)學(xué)報, 2008, 40(7):1009-1012. LV J T, LI C J, MA G F. A sliding mode PID controller design for satellite attitude regulation[J]. Journal of Harbin Institute of Technology, 2008, 40(7):1009-1012.

    [10] 畢開波, 孫劍, 隋先輝. 基于干擾觀測器的天基觀測航天器姿態(tài)控制設(shè)計[J]. 海軍航空工程學(xué)院學(xué)報, 2016, 31(1): 17-21. BI K B, SUN J, SUI X H. Design of attitude control system of space-based observation spacecraft based on disturbance observer[J]. Journal of Naval Aeronautical and Astronautical University, 2016, 31(1): 17-21.

    [11] 陸棟寧, 劉一武. 基于內(nèi)模原理的復(fù)雜撓性衛(wèi)星姿態(tài)控制研究[J]. 宇航學(xué)報, 2014, 35(3): 306-314. LU D N, LIU Y W. Research on internal model principle based attitude controlfor complex flexible satellites[J]. Journal of Astronautics, 2014, 35(3): 306-314.

    [12] 章仁為.衛(wèi)星軌道姿態(tài)動力學(xué)與控制[M]. 北京: 北京航空航天大學(xué)出版社, 1998.

    [13] 靳瀟瀟. 撓性航天器姿態(tài)抗干擾控制[D]. 北京: 北京航空航天大學(xué), 2015. JIN X X. Anti-disturbance attitude control of flexible spacecrafts[D]. Beijing: Beihang University, 2015.

    [14] 張猛, 祝曉麗, 陸嬌娣,等. 一種高穩(wěn)定度太陽帆板驅(qū)動機(jī)構(gòu)控制方法[J]. 空間控制技術(shù)與應(yīng)用, 2010, 36(4): 46-49. ZHANG M, ZHU X L, LU J D, et al. A high stability control method for solar array drive mechanism[J]. Aerospace Control and Application, 2010, 36(4): 46-49.

    [15] 斯祝華, 劉一武, 黎康. 太陽帆板驅(qū)動裝置建模及其驅(qū)動控制研究[J]. 空間控制技術(shù)與應(yīng)用, 2010, 36(2): 13-19. SI Z H, LIU Y W, LI K. Research on modeling and driver design of solar array drive assembly[J]. Aerospace Control and Application, 2010, 36(2): 13-19.

    [16] 朱仕堯, 雷勇軍, 謝燕,等. 高精度航天器太陽能電池陣減擾驅(qū)動方法[J]. 振動與沖擊, 2014, 33(21): 26-30. ZHU S Y,LEI Y J, X Y, et al. Disturbance-mitigation drive methods for solar array of a high-precision spacecraft[J]. Journal of Vibration and Shock, 2014, 33(21): 26-30.

    [17] ZHU S, LEI Y. Disturbance analysis and feed-forward compensation for the flexible solar array sun-tracking drive[J]. Institution of Mechanical Engineers Part G Journal of Aerospace Engineering, 2015, 229(14):71-76.

    [18] CANUDASD W C, OLSSON H, ASTROM K J, et al. A new model for control of systems with friction[J]. IEEE Transactions on automatic control, 1995, 40(3): 419-425.

    [19] LIU H, GUO L, ZHANG Y M. An anti-disturbance PD control scheme for attitude control and stabilization of flexible spacecrafts[J]. Nonlinear Dynamics, 2012(67): 2081-2088.

    [20] 王振. 空天飛行器復(fù)合控制技術(shù)研究[D]. 北京: 北京航空航天大學(xué), 2015. WANG Z. On composite control of aerospace vehicles[D]. Beijing: Beihang University, 2015.

    [21] TADIKONDA S S K. Articulated, flexible multibody dynamics modeling: geostationary operational environmental satellite-8 case study[J]. Journal of Guidance, Control, and Dynamics, 1997, 20(2): 276-283.

    Attitude Control of Flexible Satellites via Disturbance-Observer-BasedSliding Mode Controller

    WANG Xuebing, WU Zhong

    (SchoolofInstrumentationScienceandOpto-ElectronicsEngineering,BeihangUniversity,Beijing100191,China)

    To mitigate the effect of diverse internal and external disturbances of flexible spacecraft on attitude control performance, a composite control scheme is proposed by combining disturbance-observer (DO) and sliding mode control (SMC). The disturbance observer can observe the disturbances and output the estimation, which can be compensated through feedforward. Then, SMC is used to suppress the uncompensated disturbances, so that the attitude and rate of the spacecraft are asymptotically converged. Compared with conventional SMC, the proposed DO+SMC approach just need a smaller switching gain to restrain the remaining disturbances that the DO does not compensate, and the chattering is attenuated effectively. Simulation results demonstrate the feasibility of the designed controller.

    flexible spacecraft; attitude control; disturbance observer; sliding mode control

    資助項目(10772011). 收稿日期:2017-01-04

    V448.2

    A

    1674-1579(2017)02-0036-07

    10.3969/j.issn.1674-1579.2017.02.006

    王雪冰(1990—),女,碩士研究生,研究方向為航天器姿態(tài)控制技術(shù);吳 忠(1970—),男,教授,研究方向為航天器姿態(tài)控制、伺服控制技術(shù)等.

    猜你喜歡
    帆板撓性姿態(tài)控制
    提高帆板競技能力的教學(xué)與訓(xùn)練思考
    風(fēng)擾動下空投型AUV的飛行姿態(tài)控制研究
    多星發(fā)射上面級主動抗擾姿態(tài)控制技術(shù)研究
    什么是撓性轉(zhuǎn)子軸承
    剛性轉(zhuǎn)子和撓性轉(zhuǎn)子的區(qū)別
    基于UC/OS-II四旋翼姿態(tài)控制系統(tǒng)設(shè)計
    帆板比賽
    一種帆板驅(qū)動機(jī)構(gòu)用永磁同步電機(jī)微步控制方法
    撓性航天器退步自適應(yīng)姿態(tài)機(jī)動及主動振動控制
    國內(nèi)首臺擺動式太陽帆板驅(qū)動裝置交付衛(wèi)星使用
    久久精品国产亚洲av涩爱| 精品一品国产午夜福利视频| 捣出白浆h1v1| 午夜av观看不卡| 777米奇影视久久| 91精品伊人久久大香线蕉| 精品国产超薄肉色丝袜足j| 一个人免费看片子| 日本a在线网址| 精品国产乱码久久久久久小说| 在线观看国产h片| 97人妻天天添夜夜摸| 国产亚洲av片在线观看秒播厂| 亚洲成av片中文字幕在线观看| 久久精品国产亚洲av高清一级| 国产成人91sexporn| av国产精品久久久久影院| 9热在线视频观看99| 精品亚洲成a人片在线观看| 女性生殖器流出的白浆| 亚洲欧洲国产日韩| 精品国产国语对白av| 国产高清视频在线播放一区 | 久久久久精品人妻al黑| 精品国产乱码久久久久久男人| 成在线人永久免费视频| 免费观看av网站的网址| 国产成人精品久久二区二区91| 大陆偷拍与自拍| 免费人妻精品一区二区三区视频| 首页视频小说图片口味搜索 | 巨乳人妻的诱惑在线观看| 欧美日本中文国产一区发布| 亚洲黑人精品在线| 日韩制服骚丝袜av| 91老司机精品| a级毛片黄视频| 亚洲图色成人| 精品视频人人做人人爽| 男女无遮挡免费网站观看| 欧美激情 高清一区二区三区| 国产在线视频一区二区| 免费在线观看完整版高清| 九草在线视频观看| 伊人亚洲综合成人网| 亚洲精品自拍成人| 嫁个100分男人电影在线观看 | 黑人巨大精品欧美一区二区蜜桃| 97在线人人人人妻| 国产在线观看jvid| 久久亚洲精品不卡| 欧美人与善性xxx| 少妇裸体淫交视频免费看高清 | 久久精品久久久久久噜噜老黄| 国产成人一区二区三区免费视频网站 | 久热这里只有精品99| 国产淫语在线视频| 一区二区日韩欧美中文字幕| 十八禁网站网址无遮挡| 99精国产麻豆久久婷婷| 欧美精品人与动牲交sv欧美| 久9热在线精品视频| 亚洲精品乱久久久久久| 亚洲精品国产色婷婷电影| 中文字幕高清在线视频| 亚洲国产精品成人久久小说| 超碰97精品在线观看| 亚洲成人免费av在线播放| 蜜桃在线观看..| 成在线人永久免费视频| 亚洲天堂av无毛| 中文字幕色久视频| 国产精品免费大片| 亚洲精品日本国产第一区| 国产免费又黄又爽又色| 久久免费观看电影| av国产精品久久久久影院| 999久久久国产精品视频| 黄色视频不卡| 一本—道久久a久久精品蜜桃钙片| 日韩免费高清中文字幕av| 欧美+亚洲+日韩+国产| 亚洲人成网站在线观看播放| 国产不卡av网站在线观看| 亚洲,欧美精品.| 成年人免费黄色播放视频| 免费观看人在逋| 男人操女人黄网站| 欧美中文综合在线视频| 9色porny在线观看| 啦啦啦视频在线资源免费观看| av不卡在线播放| 国产精品一二三区在线看| 午夜福利视频在线观看免费| 男女床上黄色一级片免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产免费一区二区三区四区乱码| 女人高潮潮喷娇喘18禁视频| 久久精品成人免费网站| 高清不卡的av网站| 免费少妇av软件| 国产精品久久久久久人妻精品电影 | videosex国产| xxx大片免费视频| 国产老妇伦熟女老妇高清| 久久亚洲精品不卡| 亚洲人成电影免费在线| 色婷婷久久久亚洲欧美| 99国产精品免费福利视频| av福利片在线| 久久精品国产综合久久久| 国产野战对白在线观看| 欧美日韩成人在线一区二区| 日日夜夜操网爽| 嫩草影视91久久| 国产又色又爽无遮挡免| 欧美黑人欧美精品刺激| 黄色视频在线播放观看不卡| 久久久久久久大尺度免费视频| 久久av网站| 免费观看av网站的网址| 男人操女人黄网站| 久久热在线av| 丰满饥渴人妻一区二区三| 欧美激情 高清一区二区三区| 青草久久国产| 久久99热这里只频精品6学生| 国产一级毛片在线| 成年人免费黄色播放视频| 国产精品偷伦视频观看了| 国产片内射在线| 日本欧美国产在线视频| 巨乳人妻的诱惑在线观看| videosex国产| 亚洲人成电影免费在线| 日日摸夜夜添夜夜爱| 午夜精品国产一区二区电影| 中文字幕人妻丝袜一区二区| 在线观看免费日韩欧美大片| 亚洲精品乱久久久久久| 啦啦啦视频在线资源免费观看| 黄色 视频免费看| 在线精品无人区一区二区三| 久9热在线精品视频| 精品少妇久久久久久888优播| 国产亚洲一区二区精品| 欧美老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 日本vs欧美在线观看视频| 男男h啪啪无遮挡| 中文字幕亚洲精品专区| 蜜桃在线观看..| 日韩av在线免费看完整版不卡| 大码成人一级视频| 一级毛片我不卡| 欧美成人精品欧美一级黄| 久久天堂一区二区三区四区| 亚洲国产欧美一区二区综合| 菩萨蛮人人尽说江南好唐韦庄| 丰满迷人的少妇在线观看| 久久久久网色| 国产日韩欧美亚洲二区| 亚洲精品美女久久av网站| 亚洲欧美一区二区三区黑人| 日本午夜av视频| 欧美中文综合在线视频| 国产精品久久久久成人av| 亚洲av成人不卡在线观看播放网 | 精品亚洲成a人片在线观看| 看十八女毛片水多多多| 大片电影免费在线观看免费| 激情视频va一区二区三区| 成人免费观看视频高清| 色综合欧美亚洲国产小说| 激情视频va一区二区三区| 波野结衣二区三区在线| www.精华液| 老司机靠b影院| 国产有黄有色有爽视频| 男女国产视频网站| 老熟女久久久| 国产成人影院久久av| 免费在线观看日本一区| 久久午夜综合久久蜜桃| 久久精品国产综合久久久| 两个人免费观看高清视频| 欧美日韩一级在线毛片| 久热爱精品视频在线9| 在线天堂中文资源库| 各种免费的搞黄视频| 国产视频一区二区在线看| av有码第一页| 欧美久久黑人一区二区| 巨乳人妻的诱惑在线观看| 中文欧美无线码| 一区二区三区乱码不卡18| 1024香蕉在线观看| 色综合欧美亚洲国产小说| 老汉色∧v一级毛片| 欧美激情高清一区二区三区| 免费在线观看日本一区| 国产又色又爽无遮挡免| 大码成人一级视频| 999精品在线视频| 黄频高清免费视频| 国产av国产精品国产| 欧美精品高潮呻吟av久久| 嫁个100分男人电影在线观看 | 国产精品 国内视频| 亚洲 欧美一区二区三区| 天天影视国产精品| 日韩制服骚丝袜av| 蜜桃国产av成人99| 久久免费观看电影| 欧美人与性动交α欧美软件| 亚洲av日韩精品久久久久久密 | 丝袜喷水一区| 日韩人妻精品一区2区三区| 亚洲国产av影院在线观看| 熟女av电影| 99re6热这里在线精品视频| 久久久亚洲精品成人影院| 久久影院123| 欧美+亚洲+日韩+国产| 色网站视频免费| 亚洲伊人久久精品综合| 国产精品一区二区免费欧美 | 午夜日韩欧美国产| 人人妻,人人澡人人爽秒播 | xxx大片免费视频| 亚洲专区国产一区二区| 大话2 男鬼变身卡| 亚洲国产欧美网| 国产91精品成人一区二区三区 | 亚洲av国产av综合av卡| 亚洲午夜精品一区,二区,三区| 女人久久www免费人成看片| 久久久精品94久久精品| 搡老岳熟女国产| av福利片在线| 久热爱精品视频在线9| 欧美日韩综合久久久久久| 久久久久久久久免费视频了| 热99久久久久精品小说推荐| 又粗又硬又长又爽又黄的视频| 欧美 日韩 精品 国产| 电影成人av| 色播在线永久视频| 啦啦啦 在线观看视频| 日本猛色少妇xxxxx猛交久久| www.自偷自拍.com| 亚洲色图 男人天堂 中文字幕| 日本欧美国产在线视频| 精品久久久久久久毛片微露脸 | 精品国产国语对白av| 亚洲av美国av| 国产片内射在线| 亚洲国产毛片av蜜桃av| 亚洲国产av影院在线观看| 建设人人有责人人尽责人人享有的| 天堂中文最新版在线下载| av在线播放精品| 成人午夜精彩视频在线观看| 大片电影免费在线观看免费| 亚洲欧美色中文字幕在线| 99香蕉大伊视频| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区黑人| 日韩免费高清中文字幕av| 国产一区二区三区综合在线观看| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区| 国产一区二区在线观看av| 久久久精品区二区三区| 欧美亚洲 丝袜 人妻 在线| 在现免费观看毛片| 精品国产乱码久久久久久男人| 人妻一区二区av| 欧美日韩亚洲国产一区二区在线观看 | av视频免费观看在线观看| 日本vs欧美在线观看视频| 国产一级毛片在线| 亚洲人成电影观看| 亚洲精品一卡2卡三卡4卡5卡 | 在线观看免费日韩欧美大片| 亚洲五月色婷婷综合| 亚洲欧美精品自产自拍| 国产在视频线精品| 看免费成人av毛片| 啦啦啦中文免费视频观看日本| 蜜桃在线观看..| 欧美日韩亚洲国产一区二区在线观看 | 免费高清在线观看日韩| 亚洲中文日韩欧美视频| av有码第一页| 久久 成人 亚洲| 尾随美女入室| 亚洲av电影在线观看一区二区三区| 男女无遮挡免费网站观看| 一级毛片 在线播放| 女人被躁到高潮嗷嗷叫费观| 亚洲天堂av无毛| 午夜两性在线视频| 欧美日韩成人在线一区二区| 精品少妇一区二区三区视频日本电影| 成人亚洲精品一区在线观看| √禁漫天堂资源中文www| 秋霞在线观看毛片| 成年人免费黄色播放视频| 91麻豆av在线| 亚洲精品中文字幕在线视频| 欧美老熟妇乱子伦牲交| 精品国产乱码久久久久久小说| 水蜜桃什么品种好| 欧美黑人欧美精品刺激| 黄色片一级片一级黄色片| 2021少妇久久久久久久久久久| 国产又爽黄色视频| 51午夜福利影视在线观看| 久久久精品94久久精品| 欧美亚洲日本最大视频资源| 女性生殖器流出的白浆| 大片电影免费在线观看免费| 亚洲成人免费电影在线观看 | 免费在线观看日本一区| 亚洲成人免费av在线播放| 1024视频免费在线观看| 亚洲,一卡二卡三卡| 免费少妇av软件| 国产精品偷伦视频观看了| 99久久综合免费| 赤兔流量卡办理| 久久久久久亚洲精品国产蜜桃av| 亚洲国产精品国产精品| 少妇人妻久久综合中文| 亚洲国产精品国产精品| 国产97色在线日韩免费| 国产亚洲av片在线观看秒播厂| 亚洲国产av新网站| 精品人妻一区二区三区麻豆| 国产精品亚洲av一区麻豆| 亚洲少妇的诱惑av| 一级毛片 在线播放| 精品人妻1区二区| 精品国产一区二区久久| av天堂在线播放| 午夜老司机福利片| 老熟女久久久| 久久久久久亚洲精品国产蜜桃av| 国产成人系列免费观看| 高清av免费在线| 免费少妇av软件| 91麻豆av在线| 一区二区日韩欧美中文字幕| 久久ye,这里只有精品| 国产主播在线观看一区二区 | a级片在线免费高清观看视频| 亚洲欧洲日产国产| 久久毛片免费看一区二区三区| 国产激情久久老熟女| 精品一区二区三卡| 黑人猛操日本美女一级片| 久久亚洲精品不卡| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| 久久99热这里只频精品6学生| www.999成人在线观看| 国产精品一区二区免费欧美 | 69精品国产乱码久久久| 在线亚洲精品国产二区图片欧美| 国产精品秋霞免费鲁丝片| 欧美+亚洲+日韩+国产| 国产成人av教育| 亚洲av成人不卡在线观看播放网 | 一本综合久久免费| 亚洲欧美精品自产自拍| 日本欧美国产在线视频| 久久精品国产亚洲av高清一级| 欧美精品啪啪一区二区三区 | 国产精品一国产av| 悠悠久久av| 国产熟女午夜一区二区三区| 伊人久久大香线蕉亚洲五| 午夜老司机福利片| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲综合色网址| 成年动漫av网址| 国产精品三级大全| 亚洲成人国产一区在线观看 | 丝袜美腿诱惑在线| 久久精品亚洲熟妇少妇任你| a级毛片在线看网站| 久久国产精品影院| 欧美日韩视频精品一区| 午夜精品国产一区二区电影| av片东京热男人的天堂| 精品国产一区二区三区四区第35| 丰满人妻熟妇乱又伦精品不卡| 久久久欧美国产精品| 亚洲国产精品一区三区| 一二三四社区在线视频社区8| 亚洲,一卡二卡三卡| videosex国产| 欧美 亚洲 国产 日韩一| 男女免费视频国产| 精品少妇一区二区三区视频日本电影| 制服诱惑二区| 国产在线一区二区三区精| 亚洲免费av在线视频| 免费少妇av软件| 男女边吃奶边做爰视频| 人人妻人人爽人人添夜夜欢视频| 久久久久精品国产欧美久久久 | 19禁男女啪啪无遮挡网站| 国产一区二区在线观看av| 亚洲国产精品成人久久小说| 亚洲人成77777在线视频| 夫妻午夜视频| 伊人久久大香线蕉亚洲五| 母亲3免费完整高清在线观看| 久久99热这里只频精品6学生| 少妇精品久久久久久久| 日韩伦理黄色片| 侵犯人妻中文字幕一二三四区| 国产精品九九99| 欧美另类一区| 国产日韩欧美亚洲二区| av福利片在线| 国产精品一二三区在线看| 亚洲一码二码三码区别大吗| netflix在线观看网站| 99re6热这里在线精品视频| av片东京热男人的天堂| 免费看av在线观看网站| 国产免费福利视频在线观看| 日日夜夜操网爽| 国产又色又爽无遮挡免| 可以免费在线观看a视频的电影网站| 亚洲国产成人一精品久久久| 最近中文字幕2019免费版| 2021少妇久久久久久久久久久| 国产精品一国产av| 老司机影院成人| 黑丝袜美女国产一区| 黄色a级毛片大全视频| av欧美777| 18禁裸乳无遮挡动漫免费视频| 18禁观看日本| 最黄视频免费看| 丰满饥渴人妻一区二区三| a级毛片在线看网站| 在现免费观看毛片| 精品国产一区二区三区久久久樱花| 国产av国产精品国产| 亚洲,欧美精品.| 夫妻午夜视频| 精品国产一区二区三区四区第35| 亚洲免费av在线视频| 日韩制服丝袜自拍偷拍| 美女主播在线视频| 天堂俺去俺来也www色官网| 母亲3免费完整高清在线观看| 一本综合久久免费| 日韩,欧美,国产一区二区三区| 国产亚洲午夜精品一区二区久久| 亚洲伊人色综图| 只有这里有精品99| 国产在视频线精品| 午夜福利,免费看| 男女国产视频网站| 免费黄频网站在线观看国产| 久9热在线精品视频| 久久久久国产精品人妻一区二区| 亚洲九九香蕉| www日本在线高清视频| 免费高清在线观看视频在线观看| 欧美成人午夜精品| 亚洲激情五月婷婷啪啪| 亚洲欧美中文字幕日韩二区| 视频区欧美日本亚洲| 久久午夜综合久久蜜桃| 日本91视频免费播放| 操出白浆在线播放| 丰满饥渴人妻一区二区三| 久久国产精品大桥未久av| 一级毛片女人18水好多 | 亚洲人成电影观看| 久久精品亚洲熟妇少妇任你| 黑丝袜美女国产一区| 国产成人免费无遮挡视频| 看免费成人av毛片| 黄色a级毛片大全视频| 国产一卡二卡三卡精品| 欧美精品av麻豆av| 国产精品久久久人人做人人爽| 日本猛色少妇xxxxx猛交久久| 久久国产精品男人的天堂亚洲| 日本vs欧美在线观看视频| 老司机靠b影院| a级毛片在线看网站| 欧美 亚洲 国产 日韩一| 看十八女毛片水多多多| 亚洲欧美日韩高清在线视频 | 欧美精品av麻豆av| 久久国产精品男人的天堂亚洲| 首页视频小说图片口味搜索 | 精品一区在线观看国产| 大香蕉久久网| 国产精品久久久久久精品电影小说| 在线观看免费日韩欧美大片| 97在线人人人人妻| 在线观看国产h片| 久久国产精品大桥未久av| 亚洲欧美精品综合一区二区三区| 黄色片一级片一级黄色片| 成年动漫av网址| 国精品久久久久久国模美| 满18在线观看网站| 晚上一个人看的免费电影| 日韩一卡2卡3卡4卡2021年| 亚洲少妇的诱惑av| 两人在一起打扑克的视频| 午夜免费成人在线视频| 日本午夜av视频| 国产无遮挡羞羞视频在线观看| 性色av乱码一区二区三区2| 91精品伊人久久大香线蕉| www.av在线官网国产| 青青草视频在线视频观看| 国产又色又爽无遮挡免| 熟女少妇亚洲综合色aaa.| 巨乳人妻的诱惑在线观看| 可以免费在线观看a视频的电影网站| 欧美变态另类bdsm刘玥| 精品高清国产在线一区| 亚洲精品自拍成人| 精品人妻1区二区| 亚洲国产中文字幕在线视频| 80岁老熟妇乱子伦牲交| 制服人妻中文乱码| 欧美人与善性xxx| 国产精品一区二区免费欧美 | 精品亚洲成a人片在线观看| 女警被强在线播放| 中文字幕制服av| 最新的欧美精品一区二区| 日韩 亚洲 欧美在线| 久久久久久久久久久久大奶| 九色亚洲精品在线播放| av网站免费在线观看视频| 丝袜脚勾引网站| 亚洲精品自拍成人| 激情视频va一区二区三区| 日本午夜av视频| 欧美+亚洲+日韩+国产| 色婷婷av一区二区三区视频| 午夜福利视频在线观看免费| 中文精品一卡2卡3卡4更新| 国产一卡二卡三卡精品| 一本大道久久a久久精品| 麻豆av在线久日| 国产成人av教育| 婷婷色av中文字幕| 日本午夜av视频| 中文字幕亚洲精品专区| 别揉我奶头~嗯~啊~动态视频 | 无限看片的www在线观看| 精品少妇黑人巨大在线播放| 天堂8中文在线网| www.999成人在线观看| 啦啦啦在线免费观看视频4| 大香蕉久久网| 午夜91福利影院| 国产亚洲精品久久久久5区| 热99国产精品久久久久久7| 天天操日日干夜夜撸| 国产日韩欧美视频二区| 欧美人与性动交α欧美精品济南到| 成人国产一区最新在线观看 | 三上悠亚av全集在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品国产精品| 少妇人妻久久综合中文| 新久久久久国产一级毛片| 成年女人毛片免费观看观看9 | 欧美 亚洲 国产 日韩一| 亚洲欧美精品自产自拍| 久久久久网色| 嫁个100分男人电影在线观看 | 别揉我奶头~嗯~啊~动态视频 | 一区福利在线观看| 欧美久久黑人一区二区| 国产精品.久久久| 最近手机中文字幕大全| av网站免费在线观看视频| 婷婷色av中文字幕| 日韩一卡2卡3卡4卡2021年| 成年动漫av网址| 91老司机精品| 亚洲成色77777| 欧美 亚洲 国产 日韩一| 亚洲国产欧美一区二区综合| 午夜两性在线视频| 国产欧美日韩综合在线一区二区| 亚洲五月色婷婷综合| 国产精品偷伦视频观看了| 18禁观看日本| 在线av久久热| avwww免费| 丰满迷人的少妇在线观看| 午夜福利,免费看| 人人妻人人添人人爽欧美一区卜| 19禁男女啪啪无遮挡网站| 咕卡用的链子| 高清黄色对白视频在线免费看|