999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于DSC和MLP的欠驅動船舶自適應滑模軌跡跟蹤控制

2018-03-14 02:28:33沈智鵬
系統(tǒng)工程與電子技術 2018年3期
關鍵詞:船舶設計

沈智鵬, 王 茹

(大連海事大學信息科學技術學院, 遼寧 大連 116026)

0 引 言

船舶運動控制是控制理論應用于航海領域的一個重要的研究課題,船舶運動的目的在于提高船舶自動化水平,進而能夠保證船舶在航行過程中的安全性、經濟性以及舒適性[1]。目前,海上航行的大多數(shù)船舶僅裝備螺旋槳主推進器和舵裝置用以船舶推進和操縱,是一類典型的欠驅動系統(tǒng)。欠驅動船舶的鎮(zhèn)定控制、路徑和軌跡跟蹤控制等問題一直是不少船舶運動控制理論研究者關注的熱點。

近年來,對船舶航跡跟蹤控制問題的研究,已有諸多成果。如文獻[2-3]基于Lyapunov直接法和坐標轉換,在保證航向角速度滿足標準持續(xù)激勵條件下,設計船舶軌跡跟蹤控制律。這些控制律的設計都必須保證船舶模型參數(shù)是已知的。此外,文獻[4-5]基于backstepping思想設計船舶航跡跟蹤控制律,實現(xiàn)了對船舶軌跡跟蹤的有效控制。文獻[6]在已知外界干擾有界的情況下,利用逆推法將船舶運動的控制律分解為兩個子回路進行設計,提出了一種基于干擾界已知的滑模控制策略。但是采用傳統(tǒng)的backstepping方法設計控制器存在“微分爆炸”問題,不利于工程實踐,文獻[7-9]采用動態(tài)面控制(dynamic surface control,DSC)思想,在以往backstepping技術的基礎上加入了一階低通濾波器,避免狀態(tài)量在進行微分計算時存在“微分爆炸”問題,簡化控制律設計,易于工程實踐。

然而船舶在實際航行中,除了遭受風、浪、流等外界環(huán)境干擾外,其自身的質量、速度、慣性矩等性能也會隨之變化,導致船舶模型參數(shù)也跟著發(fā)生一定的變化。為了克服船舶模型參數(shù)的不確定性和外界風、浪、流等干擾的影響,文獻[10]在控制器的設計中利用自適應神經網絡的萬能逼近特性估計船舶動態(tài)不確定,并根據(jù)Lyapunov穩(wěn)定性原理證明船舶閉環(huán)系統(tǒng)中所有信號一致最終有界。文獻[11]針對欠驅動船舶的路徑跟蹤問題,設計了神經網絡穩(wěn)定自適應控制器,并對風、浪、流等外界環(huán)境干擾進行自適應補償。文獻[12]為實現(xiàn)在參數(shù)攝動和外界環(huán)境干擾影響下對船舶路徑跟蹤的控制,運用滑模和神經網絡這兩種技術設計了一種基于神經滑模魯棒控制器。文獻[13]在欠驅動船舶的路徑跟蹤問題上,利用神經網絡算法對強化學習的參數(shù)進行優(yōu)化,增強了所設計控制器的自適應性。文獻[14]針對無人水面艇的航跡跟蹤問題,在無人艇存在建模誤差以及外界干擾的情況下,利用神經網絡技術,設計出一種自適應徑向基函數(shù)(radial basis function,RBF)神經網絡控制器,該控制器具有較強的魯棒性。不過,采用神經網絡技術需要在線估計神經網絡的權值向量,不可避免的增加了控制算法的計算負載,也就是所謂的“維數(shù)災難”問題。這一問題可以采用最小學習參數(shù)(minimal learning parameter,MLP)方法[15]予以解決。MLP法是通過在線估計神經網絡權值向量范數(shù),從而使得所提出的自適應律不依賴于神經網絡節(jié)點數(shù)。文獻[16]在對大氣變化非常敏感的超音速飛行器中,采用DSC和MLP結合起來設計控制器,利用簡單的自適應算法就能夠實現(xiàn)對動態(tài)不確定的估計,在線更新計算負擔大大減少。文獻[17]將采用MLP代替神經網絡控制算法,應用于三自由度欠驅動船舶上,從而降低了該控制算法的計算負擔,最后,通過Lyapunov原理證明了所設計控制律的有效性。

考慮欠驅動船舶軌跡跟蹤存在模型參數(shù)不確定以及在實際航行中不可避免的遭受外部環(huán)境干擾的問題,設計了一種自適應神經網絡滑模控制方法。通過引入DSC技術避免傳統(tǒng)backstepping方法帶來“微分爆炸”問題,采用MLP方法用于避免“維數(shù)災難”問題,并以單參數(shù)在線學習代替所有權值在線學習,從而使其計算負載降低,更符合工程要求。利用Lyapunov穩(wěn)定性原理證明所設計控制律可使欠驅動船舶沿期望軌跡航行,并保證欠驅動船舶軌跡跟蹤閉環(huán)系統(tǒng)中軌跡跟蹤誤差信號一致最終有界,最后利用一艘“BAY CLASS”遠程巡邏船的仿真實驗對所設計控制器的有效性進行分析和驗證。

1 問題描述

考慮船舶前向、橫蕩和艏搖3個自由度的水平面運動,則存在模型參數(shù)不確定性和外界環(huán)境干擾影響的欠驅動水面船舶的運動學和動力學模型可表示為

(1)

式中

為實現(xiàn)欠驅動船舶軌跡控制律的設計,引入以下假設。

控制目標:針對三自由欠驅動船舶非線性數(shù)學模型(1),在滿足假設1和假設2的條件下,考慮船舶存在模型不確定和外界環(huán)境干擾的情況,設計前向推進力τu和轉向力矩τr使得所設計船舶能夠沿期望軌跡航行,并保證船舶系統(tǒng)中軌跡跟蹤誤差信號一致最終有界,實現(xiàn)對欠驅動船舶軌跡跟蹤的有效控制。

2 自適應滑模軌跡跟蹤控制律設計

為后續(xù)設計船舶系統(tǒng)控制律的需要,定義誤差變量

(2)

式中,(xd,yd)為船舶的參考軌跡;xe、ye分別為縱向位置誤差和橫向位置誤差;ue、ve分別為前向速度誤差和橫漂速度誤差,根據(jù)船舶的參考軌跡和實際位置信息可設計虛擬控制律。

為了實現(xiàn)對船舶軌跡跟蹤位置誤差趨于零,對式(2)中的xe和ye求導可得

(3)

并將ue和ve代入式(3)得

(4)

設計虛擬控制律,即

(5)

式中,k1、k2>0為設計常數(shù)。

將式(5)代入到式(4)可得

(6)

(7)

由式(7)可知,當ue、ve收斂于零時

(8)

由式(8)可知,位置誤差xe、ye將收斂于零,即實現(xiàn)對船舶軌跡跟蹤的控制。

2.1 縱向推力的最小參數(shù)自適應滑模控制律設計

(9)

式中,λ1>0為設計常數(shù)。

對一階滑模面(9)求導可得

(10)

從式(10)可以看出,如果采用傳統(tǒng)的backstepping對虛擬控制量αu直接求導,將會增加計算復雜程度。為了解決這個問題,依據(jù)動態(tài)面控制思想[18],對αu引入一階低通濾波器,并設新的狀態(tài)量Xu∈R為一階濾波器的輸出,其數(shù)學表達式為

(11)

式中,T1>0為濾波時間常數(shù)。

(12)

所產生的濾波誤差為y1=Xu-αu。

對y1關于時間求導,結合式(5)和式(12)得

(13)

式中,β1(·)為非線性連續(xù)函數(shù)。

(14)

式中,ε1>0,η1>0為設計常數(shù)。

由于式(14)中的fu存在有船舶動態(tài)不確定項,不能直接用于控制律的設計,故采用RBF神經網絡控制算法對未知項進行逼近。RBF神經網絡是一種3層前饋網絡,包括輸入層、隱含層和輸出層。RBF神經網絡自學習能力較強,可以對任意函數(shù)進行逼近,避免了對未知函數(shù)進行復雜的數(shù)學理論分析,簡化控制結構。圖1為多輸入單輸出的RBF神經網絡結構圖。

圖1 RBF神經網絡拓撲結構Fig.1 RBF neural networks topology architecture

RBF神經網絡輸入層、隱含層節(jié)點數(shù)量分別為m、n,輸出為1。x=[x1x2…xm]T∈Rm為網絡的輸入向量,則神經網絡高斯基函數(shù)的輸出表達式為

(15)

式中,bj>0為高斯基函數(shù)的寬度,j為神經網絡隱含層的第j個節(jié)點;cj=[c1,c2,…,cm]T∈Rm為第j個隱層神經元的中心點向量值,與輸入向量x的維數(shù)相同。

對于系統(tǒng)中的含有參數(shù)不確定函數(shù)fu的神經網絡輸出表達式為

(16)

式中,Z=[u,v,r]T為神經網絡的輸入,h(Z)=[h1(Z),h2(Z),…,hj(Z)]T為神經網絡的高斯基函數(shù)輸出;Wu∈Rj×1為神經網絡的理想權值,δ1為神經網絡的逼近誤差,且滿足|δ1|≤δU,δU為誤差的界值。

假設用于逼近fu的神經網絡權值Wu有界,即‖Wu‖≤WU,令

(17)

設計參數(shù)自適應律為

(18)

式中,Γu>0,ξ1>0為設計常數(shù)。根據(jù)式(17)和式(18)設計縱向推進力控制律,即

(19)

避免參數(shù)漂移的的方法有死區(qū)法和σ-修正法等,由于采用死區(qū)法要設置死區(qū)寬度,寬度過大時,會導致系統(tǒng)誤差趨于一個較大的死區(qū),影響船舶閉環(huán)系統(tǒng)的精度,故為了克服死區(qū)問題,設計帶有“σ-修正”的單參數(shù)自適應律

(20)

式中,γ1>0,σ1>0為設計常數(shù)。故縱向推進力控制律式(16)變換為

(21)

證明定義Lyapunov函數(shù)為

(22)

對于給定的正數(shù)B0、?0,考慮如下緊集:

注意到Ωd×Ω1也是緊集,且非線性連續(xù)函數(shù)β1(·)在Ωd×Ω1緊集內有最大值Nu。所以

(23)

式中,α1為正常數(shù)。

對式(22)求導,并將式(10)、式(21)、式(23)代入可得

(24)

根據(jù)文獻[19]的Young’s不等式,可得式(25)和式(26)成立

(25)

(26)

可得

(27)

(28)

將式(28)代入到式(27)可得

(29)

式中

μ1=min{(2η1-1),2(1/T1-α1,σ1}

2η1-1>0

1/T1-α1>0

(30)

證畢

2.2 轉向力矩的最小參數(shù)自適應滑模控制律設計

(31)

對二階滑模面式(31)求導可得

(32)

由式(5)得

(33)

為方便后續(xù)表達式簡潔,令

(34)

為防止對P直接求導出現(xiàn)項數(shù)滋升,對P引入一階低通濾波器,并設新的狀態(tài)量XP∈R為一階濾波器的輸出,其數(shù)學表達式為

(35)

式中,T2>0為濾波時間常數(shù)。

(36)

定義濾波誤差為y2=XP-P。

(37)

(38)

對y2關于時間求導,結合式(5)、式(36)得

(39)

式中,β2(·)為非線性連續(xù)函數(shù)。

與前面縱向推力的自適應滑模控制律設計類似,式(38)中存在含有動態(tài)不確定的fr,采用RBF神經網絡進行逼近,表達式為

(40)

式中,Wr為神經網絡的理想權值;δ2為神經網絡的逼近誤差且|δ2|≤δR。

設計帶有“σ-修正”的參數(shù)自適應律為

(41)

式中,γ2>0、σ2>0為設計常數(shù)。

(42)

式中,ε2>0、η2>0為設計常數(shù)。

證明定義Lyapunov函數(shù):

(43)

對于給定的正數(shù)B0、?2,考慮如下緊集

注意到Ωd×Ω2也是緊集,且非線性連續(xù)函數(shù)β2(·)在Ωd×Ω2緊集內有最大值NP。所以結合式(39)可得

(44)

式中,α2為正常數(shù)。

對式(43)求導并將式(38)、式(41)、式(42)、式(44)代入得

(45)

根據(jù)文獻[19]的Young’s不等式,可得式(46)和式(47)成立

(46)

(47)

代入式(45)可得

(48)

與式(28)相同,考慮不等式:

(49)

代入式(48)得

(50)

式中

μ2=min{(2η2-1),2(1/T2-α2),σ2}

2η2-1>0

1/T2-α2>0

(51)

證畢

3 穩(wěn)定性分析

綜上所述,通過設計縱向推力控制律τu和轉向力矩控制律τr,分別滿足了式(30)和式(50)可達條件,保證系統(tǒng)運動狀態(tài)可在有限時間內趨向并保持在相應的滑模面s1和s2,同時所設計的兩個滑模面是漸進收斂于原點的。所以系統(tǒng)運動點也會以指數(shù)速率收斂于原點,即滿足

(52)

當速度誤差ue、ve趨于0時,由式(8)可知,位置誤差xe、ye也將趨于0,故船舶運動的實際軌跡能夠很好沿參考軌跡航行。

為證明船舶的艏搖角速度r是有界輸入有界輸出,定義如下Lyapunov函數(shù):

(53)

由船舶模型式(1)可知

(54)

式中,Δfr為船舶動態(tài)未知不確定項,其他項為已知項,本文對fr進行了逼近,故Δfr可由式(54)得到,后續(xù)證明可直接應用Δfr。

對式(53)求導得

m33Δfr+br+τr]

(55)

4 仿真研究

為了驗證本文所設計的帶“σ-修正”的DSC和MLP自適應滑模控制方法的有效性,應用文獻[20]中一艘單槳單舵的遠程巡邏船進行仿真研究,該船長為38 m,質量為m=118×103kg,其他參數(shù)為m11=120×103kg,m22=177.9×103kg,m33=636×105kg,du=215×102kg/s,dv=147×103kg/s,dr=802×104kg/s,du2=0.2,du,du3=0.1,du,dv2=0.2dv,dv3=0.1dv,dr2=0.2dr,dr3=0.1dr,系統(tǒng)的不確定性部分取為

外界風浪流產生的干擾力和力矩取為

應用本文設計的基于DSC和MLP的自適應滑模控制器進行船舶軌跡跟蹤控制,仿真結果如圖2~圖8所示。圖2表示選用兩種不同的算法欠驅動船舶參考軌跡曲線xd、yd和實際運動軌跡曲線x、y。從圖2中曲線可以看出,本文算法可以使船舶快速的跟蹤上參考軌跡,并且保持穩(wěn)定跟蹤狀態(tài),具有較好的跟蹤能力,而采用文獻[6]中的算法基本上能跟蹤上參考軌跡,并且所用的時間較長。圖3和圖4分別為位置誤差收斂曲線和速度誤差收斂曲線,從圖3、圖4中可以看出,位置誤差和速度誤差都基本收斂于零,但采用滑模控制算法的位置誤差xe和速度誤差ue、ve在零附近波動較大,從位置誤差ye曲線圖可以看出采用滑模控制算法需200 s能夠跟蹤預期的軌跡,采用本文算法40 s就可以使船舶運動在預期軌跡上航行。

圖2 初始位置在圓內的軌跡跟蹤控制Fig.2 Trajectory tracking control of the vessel with initial position inside the circle

圖3 位置誤差收斂曲線Fig.3 Convergence of position errors

圖4 速度誤差收斂曲線Fig.4 Convergence of speed errors

圖5 軌跡跟蹤控制輸入τu變化曲線Fig.5 Trajectory tracking control input τu curves

圖6 軌跡跟蹤控制輸入τr變化曲線Fig.6 Trajectory tracking control input τr curves

圖7 縱向動態(tài)阻尼的逼近曲線 Fig.7 Approximation curve of surge dynamic damping

圖8 轉向動態(tài)阻尼的逼近曲線Fig.8 Approximation curve of yaw dynamic damping

圖5和圖6為船舶的控制輸入,表明采用本文算法轉向力矩τr在前10 s有抖振現(xiàn)象是由于船舶起始點偏離參考軌跡初始位置較大造成的,經過一段時間的調節(jié)后變?yōu)楣饣那€,符合實際的控制要求。采用滑模控制算法使船舶系統(tǒng)的控制輸入出現(xiàn)抖振現(xiàn)象,而且不可消除。圖7和圖8表示不確定函數(shù)fu和fr的真實值和估計值,從圖中可以看出帶有“σ-修正”的MLP能很好地估計船舶模型參數(shù)中的動態(tài)變量。

圖9 初始位置在圓外的軌跡跟蹤控制Fig.9 Trajectory tracking control of the vessel with initial position inside the circle

5 結 論

針對欠驅動船舶存在模型參數(shù)不確定和外界環(huán)境干擾的軌跡跟蹤控制問題,提出了一種基于DSC和MLP的自適應滑模控制方法。該方法克服了文獻[6]模型參數(shù)先驗知識的局限性,并消除了符號函數(shù)帶來的抖振問題。引入DSC的目的是為了避免傳統(tǒng)的backstepping方法對虛擬控制直接求導存在的“微分爆炸”問題,采用帶有“σ-修正”的MLP以單參數(shù)在線學習代替神經網絡控制算法中所有權值在線學習,既減少控制律的計算量,避免出現(xiàn)“維數(shù)災難”問題,又能防止參數(shù)漂移,易于工程實踐,仿真結果證明了該算法的有效性。

[1] FOSSEN T I. Handbook of marine craft hydrodynamics and motion control[M]. New York: Wiley,2011.

[2] JIANG Z P. Global tracking control of underactuated ships by Lyapunov’s direct method[J].Automatica,2002,38(2):301-309.

[3] DONG W J,GUO Y. Nonlinear tracking control of underactuated surface vessel[C]∥Proc.of the IEEE American Control Conference, 2005:4351-4356.

[4] GHOMMAM J,MNIF F,BENALI A,et al. Asymptotic back stepping stabilization of an underactuated surface vessel[J]. IEEE Trans.on Control Systems Technology,2006,14(6):1150-1157.

[5] S?RENSEN M E N,BJ?RNE E S,BREIVIK M. Performance comparison of backstepping-based adaptive controllers for marine surface vessels[C]∥Proc.of the IEEE Control Applications,2016:891-897.

[6] 朱齊丹,于瑞亭,夏桂華,等. 風浪流干擾及參數(shù)不確定欠驅動船舶航跡跟蹤的滑模魯棒控制[J]. 控制理論與應用,2012,29(7):959-964.

ZHU Q D,YU R T,XIA G H,et al. Sliding-mode robust tracking control for underactuated surface vessels with parameter uncertainties and external disturbances[J]. Control Theory & Applications,2012,29(7):959-964.

[7] AMEZQUITA S K,YAN L,BUTT W A,et al. Dynamic surface control for a class of nonlinear feedback linearizable systems with actuator failures[J]. IEEE Trans.on Neural Networks & Learning Systems,2017,28 (9):2209-2214.

[8] 杜佳璐,楊楊,胡鑫,等. 基于動態(tài)面控制的船舶動力定位控制律設計[J]. 交通運輸工程學報,2014,14(5):36-42.

DU J L,YANG Y,HU X,et al. Control law design of dynamic positioning for ship based on dynamic surface control[J]. Journal of Traffic and Transportation Engineering,2014,14(5):36-42.

[9] 王昊,王丹,彭周華,等. 多自主船協(xié)同路徑跟蹤的自適應動態(tài)面控制[J]. 控制理論與應用,2013,30(5):637-643.

WANG H,WANG D,PENG Z H,et al. Adaptive dynamic surface control for cooperative path following of multiple autonomous surface vessels[J]. Control Theory & Applications,2013,30(5):637-643.

[10] DU J L,HU X,LIU H B,et al. Adaptive robust output feedback control for a marine dynamic positioning system based on a high-gain observer[J]. IEEE Trans.on Neural Networks & Learning Systems,2015,26(11):2775-2786.

[11] 劉楊,郭晨,沈智鵬,等.欠驅動船舶路徑跟蹤的神經網絡穩(wěn)定自適應控制[J]. 控制理論與應用,2010,27(2):169-174.

LIU Y,GUO C,SHEN Z P,et al. Stable adaptive neural network control of path following for underactuated ships[J]. Control Theory & Applications,2010,27(2):169-174.

[12] 楊震,劉繁明,王巖. 欠驅動船舶路徑跟蹤的神經滑模控制[J]. 中國造船,2015,56(2):45-55.

YANG Z,LIU F M,WANG Y. Path following of underactuated surface vessels based on neural sliding mode[J]. Shipbuilding of China,2015,56(2):45-55.

[13] 沈智鵬,代昌盛. 欠驅動船舶路徑跟蹤的強化學習迭代滑模控制[J]. 哈爾濱工程大學學報,2017,38(5):697-704.

SHEN Z P,DAI C S. Iterative sliding mode control based on reinforced learning and used for path tracking of underactuated ship[J].Journal of Harbin Engineering University,2017,38(5):697-704.

[14] LI C H,ZHAO Y S, WANG G F, et al. Adaptive RBF neural network control for unmanned surface vessel course tracking[C]∥IEEE International Conference on Information Science & Technology, 2016:285-290.

[15] 劉程,李鐵山,陳納新. 帶有舵機特性的船舶航向自動舵DSC-MLP設計[J]. 哈爾濱工程大學學報,2012,33(1):9-14.

LIU C,LI T S,CHEN N X. Dynamic surface control and minimal learning parameter (DSC-MLP) design of a ship’s autopilot with rudder dynamics[J]. Journal of Harbin Engineering University,2012,33(1):9-14.

[16] XU B,YU L,WANG S X,et al. Neural dynamic surface hypersonic flight control using minimal-learning-parameter technique[C]∥Proc.of the IEEE International Conference on Unmanned Aircraft Systems,2014:960-966.

[17] ZHANG G Q,ZHANG X K. Concise robust adaptive path-following control of underactuated ships using DSC and MLP[J]. IEEE Journal of Oceanic Engineering,2014,39(4):685-694.

[18] SWAROOP D,HEDRICK J K,YIP P P,et al. Dynamic surface control for a class of nonlinear systems[J]. IEEE Trans.on Automatic Control,2000,45(10):1893-1899.

[19] CAI J,WEN C,SU H,et al. Adaptive backstepping control for a class of nonlinear systems with non-triangular structural uncertainties[J].IEEE Trans.on Automatic Control,2017,62(10): 5220-5226.

[20] DO K D,JIANG Z P,PAN J. Robust adaptive path following of underactuated ships[J]. Automatic,2004,40(6):929-944.

猜你喜歡
船舶設計
計算流體力學在船舶操縱運動仿真中的應用
基于改進譜分析法的船舶疲勞強度直接計算
何為設計的守護之道?
《豐收的喜悅展示設計》
流行色(2020年1期)2020-04-28 11:16:38
船舶!請加速
BOG壓縮機在小型LNG船舶上的應用
瞞天過海——仿生設計萌到家
藝術啟蒙(2018年7期)2018-08-23 09:14:18
設計秀
海峽姐妹(2017年7期)2017-07-31 19:08:17
有種設計叫而專
Coco薇(2017年5期)2017-06-05 08:53:16
船舶壓載水管理系統(tǒng)
中國船檢(2017年3期)2017-05-18 11:33:09
主站蜘蛛池模板: 免费不卡视频| 国产极品美女在线观看| 国产男女XX00免费观看| 好久久免费视频高清| 精品国产免费观看| 欧美成人一级| 伊人国产无码高清视频| 蜜臀AVWWW国产天堂| 欧洲一区二区三区无码| 国产成人精品三级| 国产美女免费网站| 中文字幕亚洲另类天堂| 国产午夜一级淫片| 成人免费视频一区| 久久免费观看视频| 香蕉eeww99国产在线观看| 国产亚洲精品无码专| 色一情一乱一伦一区二区三区小说| 免费毛片a| 亚洲色图另类| 九色国产在线| 亚洲v日韩v欧美在线观看| 亚洲天堂成人在线观看| 九九九久久国产精品| 亚洲欧美激情另类| 女人爽到高潮免费视频大全| 伦伦影院精品一区| 毛片免费在线| 亚洲天堂视频网站| 日韩A∨精品日韩精品无码| 欧美日韩中文字幕二区三区| 人妻精品全国免费视频| swag国产精品| 国产福利微拍精品一区二区| 国产在线自在拍91精品黑人| 97视频在线精品国自产拍| 国产在线精品美女观看| 福利一区在线| 欧美一级高清片欧美国产欧美| 刘亦菲一区二区在线观看| 国产精品成人一区二区不卡| 欧美一区二区自偷自拍视频| 亚洲天堂2014| 麻豆国产原创视频在线播放 | 91在线国内在线播放老师| 网久久综合| 亚洲成人精品久久| 久久特级毛片| 四虎亚洲精品| 九九热精品视频在线| AV无码一区二区三区四区| 美女潮喷出白浆在线观看视频| 久综合日韩| 九九线精品视频在线观看| 制服丝袜在线视频香蕉| 精品国产香蕉伊思人在线| 在线观看免费黄色网址| 欧美国产成人在线| 欧美国产日韩一区二区三区精品影视 | 亚洲午夜福利精品无码| 久久精品人妻中文系列| 亚洲中文制服丝袜欧美精品| 天天躁夜夜躁狠狠躁躁88| 国产极品粉嫩小泬免费看| 亚洲美女一级毛片| 欧洲极品无码一区二区三区| 成人久久18免费网站| аv天堂最新中文在线| 久久人人97超碰人人澡爱香蕉| 亚洲国产无码有码| 国产亚洲精久久久久久久91| 97se亚洲综合在线韩国专区福利| 国内99精品激情视频精品| 91最新精品视频发布页| 日本高清成本人视频一区| 亚洲一级毛片在线播放| 国产精品无码久久久久AV| 国产主播一区二区三区| 91青青草视频| 爆乳熟妇一区二区三区| 一级高清毛片免费a级高清毛片| 免费啪啪网址|