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

    基于FTESO和漂角補(bǔ)償?shù)拇昂较蚧?刂?/h1>
    2022-03-19 08:40:14儲(chǔ)瑞婷劉志全
    中國艦船研究 2022年1期
    關(guān)鍵詞:船舶設(shè)計(jì)

    儲(chǔ)瑞婷,劉志全

    上海海事大學(xué) 航運(yùn)技術(shù)與控制工程交通運(yùn)輸行業(yè)重點(diǎn)實(shí)驗(yàn)室,上海 201306

    0 引 言

    隨著全球海洋事業(yè)的拓展,船舶運(yùn)動(dòng)控制領(lǐng)域愈發(fā)受到關(guān)注,其中,航向跟蹤從始至終都是不可忽略的一項(xiàng)重要性能。除海面環(huán)境擾動(dòng)影響外,船舶操縱本身具有的復(fù)雜性和不確定性給船舶航向控制帶來了巨大挑戰(zhàn)。因此,探討魯棒性強(qiáng)的控制算法對船舶航向控制具有實(shí)際意義。目前,很多控制方法已在船舶控制領(lǐng)域得到廣泛應(yīng)用,例如,狀態(tài)反饋線性化[1]、反步法[2]、滑??刂疲⊿MC)[3]及自適應(yīng)[4]、模糊算法[5]和模型預(yù)測算法[6]等。為了解決反步法中輸出能量高和非線性函數(shù)難以消除的實(shí)際應(yīng)用問題,Zhang等[2]通過使用正弦函數(shù)對控制誤差進(jìn)行調(diào)制,在不改變控制器結(jié)構(gòu)的情況下提高了船舶自動(dòng)駕駛儀的控制性能。但在一般情況下,非線性船舶轉(zhuǎn)向系統(tǒng)的控制設(shè)計(jì)過程還會(huì)受到舵角、舵速率和航向角的非線性影響。Perera等[7]為了簡化非線性船舶的轉(zhuǎn)向系統(tǒng),對基于李雅普諾夫(Lyapunov)、赫爾維茨(Hurwitz)和PID的輸入輸出線性化控制方法進(jìn)行了研究,但尚未解決欠驅(qū)動(dòng)水面船舶路徑跟蹤控制中有關(guān)復(fù)雜漂角的問題[8]。

    為了補(bǔ)償漂角,最直接的方法是使用GPS、加速度計(jì)和其他傳感器等進(jìn)行測量[9],但傳感器的噪聲和高成本又使得這種方法不可行。因此,Wang等[10]提出了一種基于濾波擴(kuò)張狀態(tài)觀測器(ESO)的視線法制導(dǎo)律,其利用有限時(shí)間擴(kuò)張狀態(tài)觀測器(FTESO)來估計(jì)由海流、風(fēng)浪擾動(dòng)引起的時(shí)變漂角,從而實(shí)現(xiàn)路徑跟蹤,但未考慮舵飽和的問題。李蕓等[11]設(shè)計(jì)了一種新型的擴(kuò)張狀態(tài)觀測器,通過結(jié)合滑模控制實(shí)現(xiàn)了船舶航向有效跟蹤。與傳統(tǒng)的線性擴(kuò)張狀態(tài)觀測器(LESO)相比,F(xiàn)TESO非線性項(xiàng)的引入可以保障擴(kuò)張狀態(tài)的有限時(shí)間估計(jì)。為了提高觀測性能,Xiong等[12]提出了一種通過非線性項(xiàng)和切換項(xiàng)來反饋輸出估計(jì)誤差的新型擴(kuò)張狀態(tài)觀測器,基于該擴(kuò)張狀態(tài)觀測器的結(jié)構(gòu),將干擾觀測器設(shè)計(jì)的問題轉(zhuǎn)化為存在外部干擾時(shí)狀態(tài)觀測器的設(shè)計(jì)問題。基于機(jī)械性特性,船舶控制中存在舵飽和的問題,故不少研究人員針對帶有舵飽和特性的控制系統(tǒng)分析及設(shè)計(jì)方法進(jìn)行了廣泛研究。Liang等[13]利用FTESO估計(jì)動(dòng)態(tài)定位船在有限時(shí)間內(nèi)的不可測狀態(tài),構(gòu)造了一個(gè)輔助系統(tǒng)用來處理控制飽和狀態(tài)。此外,也可以通過飽和函數(shù)限制進(jìn)入觀測器之前的控制信號,減少執(zhí)行器約束飽和的影響[14]。

    本文將針對船舶航向控制的非線性特性和外部干擾問題,在考慮時(shí)變漂角和輸入飽和問題的基礎(chǔ)上,利用FTESO對不可測量狀態(tài)進(jìn)行估計(jì)并得到漂角估計(jì),用以修正航向誤差,然后通過設(shè)計(jì)帶有積分項(xiàng)的滑模面并結(jié)合趨近律,研究一種考慮漂角修正的欠驅(qū)動(dòng)船舶航向魯棒自適應(yīng)控制方法,并通過Lyapunov理論證明所提控制器系統(tǒng)的穩(wěn)定性。

    1 定義與引理

    定義1:齊次性。對于連續(xù)標(biāo)量函數(shù)V(x):Rn→R,若滿足以下條件:對于任意實(shí)數(shù) λ>0,存在數(shù)σ>0 和 向 量 (r1,···,rn)∈Rn>0, 其 中ri>0(i=1,2,···,n), 使 得V(λr1x1,···,λrn xn)=λdV(x1,···,xn),則 稱V關(guān)于 (r1,···rn) 具 有 齊次度 σ[15]。若向量 函 數(shù)f(x):Rn→Rn, 對于任意的 λ>0, 存在(r1,···,rn)∈Rfi(nλ>r1x01,,·其··,中λrnx rni)>=0λ(rii+=d f1i(,x21,,· ·· ·· ·,,nx)n,),i使=得1,2,f·(·x·),滿n,足其中d>?min{ri,i=1,2,···,n}, 則稱向量函數(shù)f關(guān)于(r1,···,rn)∈Rn>0 具 有齊次度d[16]。

    引理1:假設(shè)存在連續(xù)可微函數(shù)V(x)>0:U1→R,U1?U∈Rn,滿足

    式中:U表示整個(gè)狀態(tài)空間;集合U1為自變量x的定義域;Vα為V的α次冪;t為時(shí)間;常數(shù)c>0,0<α<1,則表示該系統(tǒng)局部有限時(shí)間穩(wěn)定。對于任意給定的初始條件x(t0)∈U1,收斂時(shí)間滿足T≤V1?α(x(t0),t0)/c(1?α)[17],其中t0為初始時(shí)刻。

    引理2:假設(shè)存在連續(xù)可微函數(shù)V(x)>0:U1→R,U1?U∈Rn,滿足

    式中,常數(shù)c1>0,c2>0 , 0<α<1,則表示該系統(tǒng)局部有限時(shí)間穩(wěn)定。對于任意給定的初始條件x(t0)∈{U1∩U2}, 收斂時(shí)間滿足T≤ln(1?(c2/c1)V1?α(x0,t0))/(c2α?c2), 其中U2={x|V1?α(x,t)≤c1/c2}[18]。

    引理3:對于任意的xi∈R(i=1,2,···,n),以及實(shí)數(shù) 0

    引理4:對于任意的xi∈R(i=1,2,···,n),以及實(shí)數(shù)p>1,有[20]

    2 模型建立

    基于欠驅(qū)動(dòng)船舶的六自由度模型,忽略縱蕩、橫蕩、艏搖運(yùn)動(dòng)與垂蕩、橫搖、縱搖運(yùn)動(dòng)之間的耦合,建立水面船舶縱蕩、橫蕩和艏搖方向的三自由度欠驅(qū)動(dòng)船舶模型。假定模型中慣性矩陣、附加質(zhì)量矩陣和阻尼矩陣為對角陣,且當(dāng)船舶具有3個(gè)對稱面時(shí)上述設(shè)定成立。大多數(shù)船舶的左右舷具有對稱性,但是,若船舶作水平面運(yùn)動(dòng)時(shí),則可無需考慮上下對稱性。船舶的前后非對稱性意味著慣性矩陣和阻尼矩陣的非對角項(xiàng)非零。然而,這些非對角項(xiàng)與主對角項(xiàng)相比非常小?;谏鲜鲈O(shè)定,水面船舶的數(shù)學(xué)模型可以簡化表示為[21]:

    式中:ψ為定義于固定坐標(biāo)系下的實(shí)際航向,(°);u,v分別為船舶的縱蕩和橫蕩速度, m/s;r為艏搖角速度, rad/s;mjj(1≤j≤3)為正常數(shù),表示包括附加質(zhì)量的船舶慣性系數(shù);du,dv,dr,dui,dvi,dri分別為縱蕩、橫蕩和艏搖方向上的水動(dòng)力阻尼系數(shù);未 知 時(shí) 變 項(xiàng) τwu(t), τwv(t), τwr(t)分 別 為 由 風(fēng)、浪、流引起的環(huán)境擾動(dòng);τu為 水面船舶推進(jìn)力,N,其由螺旋槳或噴水器提供;τr為艏搖力矩, N·m,其通過改變每個(gè)螺旋槳或噴水器的速度而產(chǎn)生,這表明僅能在縱蕩速度為非零時(shí)控制航向。因此,本文設(shè)計(jì)了獨(dú)立的控制系統(tǒng)來控制縱蕩速度。

    在舵位于中間位置時(shí),若沿縱剖面方向直行,因流體的對稱性,船舶將不會(huì)受到側(cè)向力的作用。在舵偏轉(zhuǎn)一個(gè)角度時(shí),因改變了水流的對稱性,舵將受到一個(gè)側(cè)向力的作用,且側(cè)向力作用點(diǎn)會(huì)對船舶重心產(chǎn)生一個(gè)繞船舶重心的力矩。在力矩的作用下,船體相對于水流會(huì)發(fā)生偏轉(zhuǎn),其縱剖面與水流速度方向?qū)⑿纬梢粋€(gè)漂角β,取值定義如下:

    為了簡化控制器的設(shè)計(jì),本文將利用擴(kuò)張狀態(tài)觀測器對模型中艏搖方向的總擾動(dòng)進(jìn)行估計(jì)。因此,針對航向控制器設(shè)計(jì),將式(5)的三自由度欠驅(qū)動(dòng)船舶模型簡化為

    假設(shè)1:外界干擾有界,即 |τwu(t)|≤τwumax<∞,|τwv(t)|≤τwvmax<∞, |τwr(t)|≤τwrmax<∞。

    假設(shè)2:因外部干擾的能量有限,故合理假設(shè)w的導(dǎo)數(shù)有界,即 ‖w˙‖≤Δ,Δ為w˙的上界。

    3 控制系統(tǒng)

    3.1 航向預(yù)濾波與漂角補(bǔ)償

    為了避免較大的航向變化率,本文在控制系統(tǒng)中加入了預(yù)濾波器以平滑期望航向,可保證滑??刂破鞯聂敯粜浴nA(yù)濾波器的主要作用是對原期望航向ψr進(jìn)行濾波,用于實(shí)現(xiàn)濾波后參考航向ψd和航向變化率ψ˙d的平滑過渡,從而避免更高的控制增益要求,提高控制器性能。本文采用的二階預(yù)濾波器形式如下:

    式中, λi(i=1,2,3)為待確定的預(yù)濾波器參數(shù)。

    因外部干擾引起的時(shí)變漂角較小,故在船舶航向控制中通常忽略漂角的影響而假定其為0,但實(shí)際上漂角對跟蹤誤差具有一定的影響。因此,本文考慮了帶漂角修正的航向控制,如圖1所示。圖中,ψh為實(shí)際航向,ψd為預(yù)濾波后的參考航向,初始航向誤差為 ψh?ψd??紤]漂角存在時(shí),船舶實(shí)際運(yùn)動(dòng)方向與期望切線方向相差了一個(gè)漂角[8],則基于漂角修正后的期望航向ψda與 航向誤差ea分別為:

    圖1 漂角補(bǔ)償后的期望航向Fig.1 The desired heading with sideslip angle compensation

    3.2 控制器

    3.2.1 滑??刂坡稍O(shè)計(jì)

    為了保證船舶航向的自動(dòng)保持與跟蹤特點(diǎn),使用閉環(huán)反饋控制系統(tǒng)研究了航向的跟蹤問題。根據(jù)滑??刂圃?,利用誤差反饋思想進(jìn)行了滑模控制器設(shè)計(jì)。本文設(shè)計(jì)的帶有積分項(xiàng)的滑模面s如下:

    式中: 0

    在滑??刂浦?,會(huì)不可避免地出現(xiàn)抖振現(xiàn)象,之所以如此,主要是因?yàn)椴贿B續(xù)的符號函數(shù)項(xiàng)的增益需要足夠大來保證魯棒性。為了減小滑模抖振的影響,采用趨近律方法改善系統(tǒng)的抖振。本文使用基于指數(shù)趨近律的滑??刂疲员WC運(yùn)動(dòng)點(diǎn)在有限的時(shí)間內(nèi)到達(dá)滑動(dòng)模態(tài)。

    3.2.2 擴(kuò)張狀態(tài)觀測器設(shè)計(jì)

    為保障跟蹤精度,采用FTESO對由不確定項(xiàng)和外界干擾構(gòu)成的總擾動(dòng)項(xiàng)進(jìn)行估計(jì)和補(bǔ)償。結(jié)合滑??刂破?,得到系統(tǒng)控制結(jié)構(gòu)如圖2所示。

    圖2 船舶航向控制系統(tǒng)示意圖Fig.2 Schematic diagram of the heading control system for ships

    定義 η=[u,v,r]T,w? 為系統(tǒng)擴(kuò)張的狀態(tài)向量,包含有未知外界干擾和內(nèi)部未建模動(dòng)態(tài)項(xiàng);觀測器估計(jì)誤差z1=η?η?, 擴(kuò)張狀態(tài)誤差z2=w?w?,構(gòu)造如下形式的FTESO:

    令σ =α1α2。由定義1,可推斷出系統(tǒng)(式(17))關(guān)于權(quán)值 (1,α1) 具有齊次度 α1?1。定義可微函數(shù)[[sigVσ1α(z=1)Z]TT,P[sZig,σ1α1且(z2)V]Tα]T>,0,P為其L中yapuZno=v[ZT1等,ZT2式]T=ATP+PA=?I6的 解,其 中Ii(i=2,3,4,···)為i階 單位矩陣,定義系統(tǒng)矩陣

    可知A為Hurwitz矩陣。根據(jù)文獻(xiàn)[15],可知Vα即 為系統(tǒng)(式(17))的Lyapunov函數(shù),并令fα為系統(tǒng)(式(17))的向量場,則LfαVα為Vα沿 著向量場fα的李導(dǎo)數(shù)。由此可以推斷出,Vα,LfαVα關(guān) 于權(quán)值(1,α1)分別具有齊次度 2 /σ, 2/σ+(α1?1)。由文獻(xiàn)[22],可以得到以下不等式:

    設(shè)計(jì)誤差系統(tǒng)(式(16))的Lyapunov函數(shù)如下:

    對式(19)求導(dǎo),得

    將式(18)代入式(20),得

    將式(22)代入式(21),得

    其中,

    根據(jù)引理3,有

    最終,速度估計(jì)誤差和擴(kuò)張狀態(tài)觀測誤差可以收斂到緊致集Ω內(nèi),Ω的表達(dá)式如下:

    3.3 系統(tǒng)穩(wěn)定性分析

    基于 FTESO 的滑??刂破髟O(shè)計(jì)表述為如下定理。

    定理1: 對于船舶航向控制系統(tǒng)(式(7)),設(shè)計(jì)擴(kuò)張狀態(tài)觀測器(式(15)),利用估計(jì)漂角β?=arctan(v?/u?)選取滑模面為式(11) ,設(shè)計(jì)滑??刂破魅缦拢?/p>

    由于FTESO保證了各項(xiàng)估計(jì)誤差都足夠小,因此,只要選取適當(dāng)大的k2就 能保證k2≥(wr?

    為了得到更好的控制效果,選用一個(gè)陡峭的飽和函數(shù),定義為

    來逼近控制器中的符號函數(shù),其中ε1為極小值。在極限情況下,當(dāng)ε1趨于0時(shí),飽和函數(shù)近似為符號函數(shù)。即控制律取為

    4 仿真結(jié)果分析

    為了驗(yàn)證所提控制器的有效性,考慮了整個(gè)航向系統(tǒng)。本文使用Matlab軟件中的Simulink建模工具搭建了整個(gè)系統(tǒng)模型,并進(jìn)行了仿真應(yīng)用。模型(式(5))中的相關(guān)參數(shù)見文獻(xiàn)[22]:m21151=×1 10220,×d1v0=3,1m4272×=1 10737,.9d×r=1 0830,m23×31=064;3 6d×u21=050;.d2ud=u,du3=0.1du;dv2=0.2dv,dv3=0.1dv;dr2=0.2dr,dr3=0.1dr。本文中所需設(shè)計(jì)和選擇的參數(shù)見表1??v蕩速度由獨(dú)立的控制系統(tǒng)控制,本文設(shè)計(jì)的基于擴(kuò)張狀態(tài)觀測器的PID控制器控制船舶保持7 m/s的縱蕩速度,另外,選擇艏搖控制輸入約束為τm=7×108N·m。 設(shè)定初始航速為 6 m/s,期望航向角分別為20?, ?20?和 0?,仿真步長為0.01 s,仿真時(shí)間為 1 50 s。 假設(shè)仿真中的環(huán)境擾動(dòng)為τwu=150×103dr, τwr=900×105dr,τwv=2×103×(1+0.3 cos(0.4t)), 其中dr為三級海浪模型,描述如下:

    表1 不同控制方法下的參數(shù)設(shè)計(jì)Table 1 The parameters of different controllers

    式中: 為三級海浪模型的表達(dá)式; 指功率譜密度為0.1的零均值高斯白噪聲過程; 為一個(gè)二階波的傳遞函數(shù),y(s) ω(s)h(s)

    為了保證仿真結(jié)果的可靠性,本文選取基于線性擴(kuò)張狀態(tài)觀測器的反步法與本文設(shè)計(jì)的基于FTESO的滑??刂品椒ㄟM(jìn)行了仿真對比,仿真結(jié)果如圖3~圖7所示。在不同的仿真模型中,相同的初始條件和相關(guān)參數(shù)值保持了一致。圖中,“SMC+FTESO?漂角修正”為考慮漂角修正的滑??刂品椒ǎ癝MC+FTESO?無漂角修正”為不考慮漂角修正的滑??刂品椒?,“反步+LESO”為結(jié)合了線性擴(kuò)張狀態(tài)觀測器的反步法 。

    由圖3(a)可以看出,考慮了漂角修正后的船舶航向能更快地到達(dá)期望航向(在航向發(fā)生變化后的4 s內(nèi)即達(dá)到期望值),而未考慮漂角修正的方法和傳統(tǒng)的反步法則需要8 s甚至是更長的時(shí)間才能到達(dá)穩(wěn)定狀態(tài)。圖3(b)所示為未考慮輸入飽和的控制器艏搖力矩。從圖中可見,當(dāng)航向發(fā)生變化時(shí),力矩遠(yuǎn)遠(yuǎn)超出了舵能正常使用的極限值。為解決此問題,本文在設(shè)計(jì)帶漂角修正的航向控制器時(shí)增加了輸入約束來進(jìn)行相應(yīng)的處理。

    圖3 無約束下控制器的航向和控制力矩對比Fig.3 Comparison of heading angles and yaw torque of the controllers without input constraint

    圖4所示為考慮輸入飽和后(有約束)3種不同控制器的航向及其誤差對比。從圖中可以看出,對控制輸入進(jìn)行約束后,航向及其誤差較約束前未發(fā)生變化。

    圖4 有輸入約束下控制器的航向及其誤差對比Fig.4 Comparison of heading angles and their errors of the controllers with input constraint

    圖5所示為通過觀測器估計(jì)得到的漂角值。從圖中可以看出,漂角一直存在于船舶航行過程中,雖然其值較小,但當(dāng)航向發(fā)生變化時(shí)漂角也會(huì)發(fā)生突變。因此,考慮漂角修正的船舶航向控制能夠改善航向控制性能,減小航向誤差。約束后的艏搖力矩如圖6所示,從圖中可以看出,艏搖力矩得到了相應(yīng)的減小。

    圖5 漂角估計(jì)值Fig.5 The estimation of sideslip angle

    圖6 有輸入約束的艏搖控制力矩對比Fig.6 Comparison of yaw torque of the controllers with input constraint

    從最終的航向跟蹤效果來看,本文設(shè)計(jì)的滑??刂品椒ㄅc基于線性擴(kuò)張狀態(tài)觀測器的傳統(tǒng)反步法相比優(yōu)勢并不明顯,但從觀測器的角度分析,本文設(shè)計(jì)的FTESO觀測性能更好。圖7示出了2種觀測器對縱蕩、橫蕩、艏搖這3個(gè)方向上的速度估計(jì)值及其誤差。FTESO對艏搖角速度的估計(jì)誤差穩(wěn)定在 ±4×10?5rad/s內(nèi),而傳統(tǒng)LESO在初始時(shí)間其誤差最大達(dá)1 .4×10?3rad/s ,經(jīng)過0.98 s的調(diào)節(jié)時(shí)間后才穩(wěn)定至 ±6×10?5rad/s內(nèi)波動(dòng)。因考慮了漂角修正,觀測器對縱蕩速度u和橫蕩速度v的估計(jì)性能也極為重要。FTESO對縱蕩速度的估計(jì)誤差始終保持在 ±4×10?5m/s內(nèi),傳統(tǒng)LESO在初始時(shí)間其誤差最大達(dá)1 .12×10?3m/s,收斂時(shí)間為 0.99 s, 誤差穩(wěn)定在 ±6×10?5m/s內(nèi)波動(dòng)。盡管在對橫蕩速度v的估計(jì)上LESO的觀測誤差最終能穩(wěn)定在 ±2×10?7m/s內(nèi) ,但1 .15 s的收斂時(shí)間相對FTESO的 0.06 s仍較長,且初始時(shí)刻的誤差達(dá)到4.3×10?3m/s , 而FTESO的最大估計(jì)誤差僅6.68×10?5m/s。

    圖7 有輸入約束時(shí)縱蕩、橫蕩及艏搖方向上的速度估計(jì)值及其誤差Fig.7 The estimations and errors of velocity in surge, sway and yaw directions with input constraint

    5 結(jié) 語

    本文以水面三自由度欠驅(qū)動(dòng)船舶運(yùn)動(dòng)模型為對象,在航向跟蹤控制器設(shè)計(jì)中考慮了漂角修正,不同于利用相對速度對漂角的直接計(jì)算,其采用FTESO對漂角進(jìn)行實(shí)時(shí)估計(jì),從而極大地改善了船舶的航向跟蹤性能。仿真結(jié)果表明,相對于傳統(tǒng)的線性觀測器,F(xiàn)TESO的觀測性能良好;而對于總擾動(dòng)項(xiàng)的估計(jì),則在基于趨近律的滑模控制器設(shè)計(jì)中予以補(bǔ)償,削弱了抖振影響,使得航向跟蹤誤差在收斂至0的同時(shí),還可保證航向控制系統(tǒng)有較強(qiáng)的抗干擾能力。此外,本文利用飽和函數(shù)對艏搖力矩進(jìn)行了約束,具有較好的實(shí)際意義。雖然本文對輸入力矩進(jìn)行了一定的約束,但不難看出仍有改進(jìn)的空間,下一步的工作將繼續(xù)研究輸入飽和問題,并通過改進(jìn)算法來提高控制精度。

    猜你喜歡
    船舶設(shè)計(jì)
    計(jì)算流體力學(xué)在船舶操縱運(yùn)動(dòng)仿真中的應(yīng)用
    基于改進(jìn)譜分析法的船舶疲勞強(qiáng)度直接計(jì)算
    何為設(shè)計(jì)的守護(hù)之道?
    《豐收的喜悅展示設(shè)計(jì)》
    流行色(2020年1期)2020-04-28 11:16:38
    船舶!請加速
    BOG壓縮機(jī)在小型LNG船舶上的應(yīng)用
    瞞天過?!律O(shè)計(jì)萌到家
    設(shè)計(jì)秀
    海峽姐妹(2017年7期)2017-07-31 19:08:17
    有種設(shè)計(jì)叫而專
    Coco薇(2017年5期)2017-06-05 08:53:16
    船舶壓載水管理系統(tǒng)
    中國船檢(2017年3期)2017-05-18 11:33:09

    久久午夜亚洲精品久久| 夫妻午夜视频| 超色免费av| 成年人免费黄色播放视频| 国产精品自产拍在线观看55亚洲 | 在线av久久热| 国产精品二区激情视频| 一级,二级,三级黄色视频| 又紧又爽又黄一区二区| 色婷婷av一区二区三区视频| 欧美日韩精品网址| 大香蕉久久成人网| 亚洲国产精品sss在线观看 | 午夜福利免费观看在线| 国产精品电影一区二区三区 | 18在线观看网站| 日韩中文字幕欧美一区二区| 免费在线观看亚洲国产| 国产精品国产av在线观看| 亚洲中文日韩欧美视频| 嫁个100分男人电影在线观看| 日本a在线网址| 久久精品亚洲熟妇少妇任你| 中文字幕最新亚洲高清| av网站免费在线观看视频| 精品视频人人做人人爽| 黄色成人免费大全| 男男h啪啪无遮挡| 女人被狂操c到高潮| 1024视频免费在线观看| 国产欧美日韩精品亚洲av| 在线视频色国产色| 亚洲成人免费电影在线观看| av不卡在线播放| 国精品久久久久久国模美| 麻豆成人av在线观看| 亚洲欧美色中文字幕在线| 人人妻人人澡人人爽人人夜夜| 中文欧美无线码| 国产精品欧美亚洲77777| 国产成+人综合+亚洲专区| 三上悠亚av全集在线观看| 亚洲第一青青草原| 这个男人来自地球电影免费观看| 久久ye,这里只有精品| 欧美成狂野欧美在线观看| 91麻豆av在线| 下体分泌物呈黄色| 日韩 欧美 亚洲 中文字幕| 成人免费观看视频高清| 国产精品秋霞免费鲁丝片| 午夜免费鲁丝| 成人av一区二区三区在线看| 捣出白浆h1v1| bbb黄色大片| 黄色a级毛片大全视频| 亚洲欧美精品综合一区二区三区| 国产免费av片在线观看野外av| 亚洲一码二码三码区别大吗| 免费久久久久久久精品成人欧美视频| av免费在线观看网站| 欧美精品高潮呻吟av久久| 深夜精品福利| 久久精品亚洲av国产电影网| 91老司机精品| 悠悠久久av| 亚洲第一青青草原| 欧美精品人与动牲交sv欧美| 一级a爱片免费观看的视频| 9191精品国产免费久久| aaaaa片日本免费| 亚洲,欧美精品.| 中文字幕制服av| 色老头精品视频在线观看| 亚洲 欧美一区二区三区| www.精华液| 国产成人av教育| 国产精品成人在线| 国产成人精品无人区| 大香蕉久久网| 精品人妻熟女毛片av久久网站| 国产欧美日韩一区二区三| 国产男靠女视频免费网站| 美女高潮到喷水免费观看| 午夜精品在线福利| 啦啦啦视频在线资源免费观看| 丝袜在线中文字幕| 女性生殖器流出的白浆| 两人在一起打扑克的视频| 亚洲av成人一区二区三| www.自偷自拍.com| 男女之事视频高清在线观看| 亚洲精品国产一区二区精华液| 欧美乱色亚洲激情| 在线观看免费视频网站a站| 亚洲欧美一区二区三区黑人| 国产精品欧美亚洲77777| av免费在线观看网站| 国产一卡二卡三卡精品| 日韩大码丰满熟妇| 免费一级毛片在线播放高清视频 | 成人18禁在线播放| 老司机靠b影院| 在线免费观看的www视频| 超色免费av| 久久久久久久久免费视频了| 老熟女久久久| а√天堂www在线а√下载 | 国产亚洲欧美在线一区二区| 久久天躁狠狠躁夜夜2o2o| 久久这里只有精品19| 久久久国产成人免费| 日本撒尿小便嘘嘘汇集6| 欧美成人免费av一区二区三区 | 久久精品亚洲熟妇少妇任你| 身体一侧抽搐| 一级a爱片免费观看的视频| 久久久久精品人妻al黑| cao死你这个sao货| 国产欧美日韩综合在线一区二区| 亚洲综合色网址| 国产精品电影一区二区三区 | 国产精品 国内视频| 国产精品美女特级片免费视频播放器 | 三级毛片av免费| 性少妇av在线| 久久草成人影院| 两个人看的免费小视频| 亚洲精品av麻豆狂野| 91精品国产国语对白视频| 又黄又爽又免费观看的视频| 中文字幕制服av| 两个人免费观看高清视频| 黄色视频,在线免费观看| 欧美日韩成人在线一区二区| 国产精品永久免费网站| 美女福利国产在线| 午夜两性在线视频| 日韩免费av在线播放| 99re在线观看精品视频| 久久久国产精品麻豆| 欧美日韩成人在线一区二区| 久久精品国产99精品国产亚洲性色 | 99精品久久久久人妻精品| 制服诱惑二区| 亚洲avbb在线观看| 亚洲 欧美一区二区三区| 国产精品免费视频内射| 久久久国产欧美日韩av| 国产精品免费一区二区三区在线 | 如日韩欧美国产精品一区二区三区| 欧美黄色淫秽网站| 丁香六月欧美| 麻豆乱淫一区二区| 亚洲欧美色中文字幕在线| 免费在线观看完整版高清| 最近最新免费中文字幕在线| 看片在线看免费视频| 极品少妇高潮喷水抽搐| 一边摸一边抽搐一进一小说 | а√天堂www在线а√下载 | 久久99一区二区三区| 黄色丝袜av网址大全| 久久婷婷成人综合色麻豆| 久久国产精品大桥未久av| 国产亚洲欧美精品永久| 纯流量卡能插随身wifi吗| 久久人妻福利社区极品人妻图片| 久久热在线av| 日韩成人在线观看一区二区三区| 欧美色视频一区免费| 在线观看舔阴道视频| 久久狼人影院| 波多野结衣一区麻豆| 久久精品国产综合久久久| 精品一区二区三卡| 中文字幕人妻熟女乱码| a级毛片黄视频| 99re6热这里在线精品视频| 亚洲欧美激情综合另类| 亚洲精品在线美女| 69av精品久久久久久| 黄片播放在线免费| 美女午夜性视频免费| 又大又爽又粗| 国产免费男女视频| 大香蕉久久网| 伦理电影免费视频| 国产精品久久久av美女十八| 午夜精品久久久久久毛片777| 搡老熟女国产l中国老女人| 亚洲片人在线观看| 99热网站在线观看| 亚洲欧美激情在线| 满18在线观看网站| 少妇猛男粗大的猛烈进出视频| 一级a爱片免费观看的视频| 成人手机av| 一本综合久久免费| 国产一区二区三区综合在线观看| 国产精品免费大片| 婷婷精品国产亚洲av在线 | 悠悠久久av| а√天堂www在线а√下载 | av超薄肉色丝袜交足视频| 欧美最黄视频在线播放免费 | 一本一本久久a久久精品综合妖精| 日本撒尿小便嘘嘘汇集6| 九色亚洲精品在线播放| 少妇裸体淫交视频免费看高清 | 咕卡用的链子| 亚洲成av片中文字幕在线观看| 亚洲精品乱久久久久久| 老熟妇仑乱视频hdxx| 亚洲一码二码三码区别大吗| 9191精品国产免费久久| 亚洲精品美女久久av网站| 国产成人欧美在线观看 | 日本vs欧美在线观看视频| 亚洲五月色婷婷综合| 91成年电影在线观看| 在线播放国产精品三级| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 色婷婷久久久亚洲欧美| 91成年电影在线观看| 伦理电影免费视频| 99re在线观看精品视频| 国产成人免费观看mmmm| 老熟妇仑乱视频hdxx| av网站免费在线观看视频| 成人精品一区二区免费| 高清视频免费观看一区二区| 免费av中文字幕在线| 丝瓜视频免费看黄片| av超薄肉色丝袜交足视频| 午夜日韩欧美国产| 久久久久精品人妻al黑| 高清在线国产一区| 啦啦啦免费观看视频1| 欧美日韩成人在线一区二区| 看黄色毛片网站| 中文亚洲av片在线观看爽 | 欧美日韩亚洲国产一区二区在线观看 | 黄片播放在线免费| 亚洲欧美激情在线| 亚洲av美国av| 后天国语完整版免费观看| 人人澡人人妻人| 中文字幕av电影在线播放| 国产视频一区二区在线看| 一级a爱片免费观看的视频| 久久中文字幕人妻熟女| 最新在线观看一区二区三区| 日日摸夜夜添夜夜添小说| 飞空精品影院首页| 捣出白浆h1v1| 如日韩欧美国产精品一区二区三区| 久久热在线av| 黄网站色视频无遮挡免费观看| 国产真人三级小视频在线观看| 国产亚洲精品久久久久5区| 精品午夜福利视频在线观看一区| 不卡av一区二区三区| 久久亚洲精品不卡| 午夜精品在线福利| 变态另类成人亚洲欧美熟女 | 亚洲熟妇熟女久久| 91大片在线观看| 波多野结衣av一区二区av| 在线观看日韩欧美| 免费在线观看影片大全网站| 久久性视频一级片| 亚洲五月色婷婷综合| 亚洲av第一区精品v没综合| 亚洲情色 制服丝袜| 国产淫语在线视频| 成年女人毛片免费观看观看9 | 午夜亚洲福利在线播放| 丝袜美足系列| 亚洲精品粉嫩美女一区| 我的亚洲天堂| 免费在线观看完整版高清| 亚洲avbb在线观看| 国产又色又爽无遮挡免费看| 少妇的丰满在线观看| 嫩草影视91久久| 香蕉国产在线看| 九色亚洲精品在线播放| 亚洲一码二码三码区别大吗| 免费日韩欧美在线观看| 久久人妻福利社区极品人妻图片| 亚洲欧美一区二区三区久久| 男女床上黄色一级片免费看| 搡老乐熟女国产| 在线视频色国产色| aaaaa片日本免费| 欧美日韩亚洲高清精品| 91国产中文字幕| 国产精品久久久久成人av| 欧美激情极品国产一区二区三区| xxxhd国产人妻xxx| 人人妻人人爽人人添夜夜欢视频| 欧美大码av| 精品国产一区二区三区四区第35| bbb黄色大片| 日韩免费高清中文字幕av| 欧美黑人欧美精品刺激| 免费久久久久久久精品成人欧美视频| 免费在线观看影片大全网站| 亚洲精华国产精华精| 在线十欧美十亚洲十日本专区| 日韩视频一区二区在线观看| 啦啦啦视频在线资源免费观看| 人人妻人人爽人人添夜夜欢视频| 亚洲情色 制服丝袜| 另类亚洲欧美激情| 国产精品 欧美亚洲| 国产真人三级小视频在线观看| 亚洲免费av在线视频| 一进一出抽搐gif免费好疼 | 国产精品电影一区二区三区 | 国产亚洲精品一区二区www | 欧美成人午夜精品| 中国美女看黄片| 色婷婷久久久亚洲欧美| 99国产精品一区二区蜜桃av | 两性夫妻黄色片| 国产亚洲av高清不卡| 欧美精品人与动牲交sv欧美| 欧美精品亚洲一区二区| 国产精品 国内视频| 美女国产高潮福利片在线看| av网站免费在线观看视频| 国产精品乱码一区二三区的特点 | 国产一区二区激情短视频| 精品国产乱子伦一区二区三区| 欧美激情 高清一区二区三区| 色尼玛亚洲综合影院| 日本欧美视频一区| 欧美另类亚洲清纯唯美| 久久精品国产亚洲av香蕉五月 | 在线天堂中文资源库| 91老司机精品| 日本黄色日本黄色录像| 日日爽夜夜爽网站| 美女高潮到喷水免费观看| av国产精品久久久久影院| 亚洲熟妇中文字幕五十中出 | 亚洲成人免费电影在线观看| 人成视频在线观看免费观看| 少妇猛男粗大的猛烈进出视频| 欧美午夜高清在线| 国产区一区二久久| 亚洲熟女毛片儿| 国产三级黄色录像| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 黑丝袜美女国产一区| 国产片内射在线| 亚洲av成人不卡在线观看播放网| 夜夜爽天天搞| 在线国产一区二区在线| 亚洲视频免费观看视频| 日本wwww免费看| 国产精品1区2区在线观看. | 性少妇av在线| 黄网站色视频无遮挡免费观看| 叶爱在线成人免费视频播放| 两个人看的免费小视频| 国产又爽黄色视频| 欧美日韩福利视频一区二区| 久热爱精品视频在线9| 老司机影院毛片| 真人做人爱边吃奶动态| 成在线人永久免费视频| 国产成+人综合+亚洲专区| 少妇裸体淫交视频免费看高清 | 黄网站色视频无遮挡免费观看| a在线观看视频网站| 中文亚洲av片在线观看爽 | 国产日韩欧美亚洲二区| 18禁裸乳无遮挡动漫免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 婷婷丁香在线五月| 国产高清激情床上av| 法律面前人人平等表现在哪些方面| 一a级毛片在线观看| 人妻久久中文字幕网| 午夜福利欧美成人| 777米奇影视久久| www.自偷自拍.com| tocl精华| 黄色怎么调成土黄色| av电影中文网址| 色婷婷久久久亚洲欧美| 成人免费观看视频高清| 欧洲精品卡2卡3卡4卡5卡区| 女人被躁到高潮嗷嗷叫费观| 在线看a的网站| 热re99久久国产66热| 亚洲男人天堂网一区| 国产单亲对白刺激| 18禁观看日本| 大香蕉久久成人网| 两性午夜刺激爽爽歪歪视频在线观看 | 王馨瑶露胸无遮挡在线观看| 老汉色av国产亚洲站长工具| 亚洲精品乱久久久久久| 热re99久久精品国产66热6| 亚洲九九香蕉| 黑人欧美特级aaaaaa片| 精品久久久久久,| 日本一区二区免费在线视频| 一区二区三区激情视频| 亚洲国产欧美一区二区综合| 成人精品一区二区免费| 精品乱码久久久久久99久播| 夫妻午夜视频| 99久久人妻综合| 欧美亚洲 丝袜 人妻 在线| 国产成+人综合+亚洲专区| 免费高清在线观看日韩| 亚洲av成人一区二区三| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 久久久久久免费高清国产稀缺| 亚洲七黄色美女视频| 窝窝影院91人妻| 一进一出好大好爽视频| 操出白浆在线播放| 国产乱人伦免费视频| 国产欧美日韩精品亚洲av| 夫妻午夜视频| 嫩草影视91久久| 国产黄色免费在线视频| 午夜91福利影院| 欧美日本中文国产一区发布| 国产高清国产精品国产三级| 久久国产精品男人的天堂亚洲| 一区在线观看完整版| 一级黄色大片毛片| 久久人妻福利社区极品人妻图片| 免费一级毛片在线播放高清视频 | 黄色怎么调成土黄色| 日韩 欧美 亚洲 中文字幕| 日本vs欧美在线观看视频| 欧美亚洲日本最大视频资源| 国产极品粉嫩免费观看在线| 在线视频色国产色| 最近最新中文字幕大全电影3 | www.精华液| 免费在线观看完整版高清| 免费一级毛片在线播放高清视频 | 亚洲欧洲精品一区二区精品久久久| 亚洲片人在线观看| 黄色a级毛片大全视频| 午夜福利免费观看在线| 亚洲一区中文字幕在线| 国产精华一区二区三区| 老熟女久久久| 久久久久久亚洲精品国产蜜桃av| ponron亚洲| 国产欧美日韩一区二区精品| 成人特级黄色片久久久久久久| 一边摸一边抽搐一进一小说 | 亚洲精华国产精华精| 亚洲精品一卡2卡三卡4卡5卡| 一级,二级,三级黄色视频| 亚洲欧美精品综合一区二区三区| 久久久久久久精品吃奶| 欧美激情久久久久久爽电影 | 一级a爱视频在线免费观看| 久久中文字幕一级| 在线观看免费高清a一片| 亚洲成av片中文字幕在线观看| 国产高清videossex| 国产免费av片在线观看野外av| www.精华液| 人妻 亚洲 视频| 熟女少妇亚洲综合色aaa.| 亚洲精品久久成人aⅴ小说| 亚洲综合色网址| 亚洲第一av免费看| 欧美成人免费av一区二区三区 | 一个人免费在线观看的高清视频| 黄色a级毛片大全视频| 亚洲精品一二三| √禁漫天堂资源中文www| 啦啦啦免费观看视频1| 久久久国产成人精品二区 | 国产又爽黄色视频| 十分钟在线观看高清视频www| 日韩欧美国产一区二区入口| 99热国产这里只有精品6| 在线观看www视频免费| 久久久久国内视频| 国产精品1区2区在线观看. | 国产高清videossex| 日韩欧美一区视频在线观看| 亚洲avbb在线观看| 久久狼人影院| 多毛熟女@视频| 久久精品aⅴ一区二区三区四区| 成人三级做爰电影| 久久狼人影院| 午夜成年电影在线免费观看| 一级a爱片免费观看的视频| 欧美丝袜亚洲另类 | av在线播放免费不卡| 亚洲精品乱久久久久久| 亚洲精品在线观看二区| 精品久久久久久电影网| 亚洲专区中文字幕在线| 在线观看免费午夜福利视频| 每晚都被弄得嗷嗷叫到高潮| 国产色视频综合| 亚洲精品国产区一区二| 黄片大片在线免费观看| 午夜两性在线视频| 在线观看免费日韩欧美大片| 午夜两性在线视频| 久久性视频一级片| 国产亚洲精品第一综合不卡| xxx96com| 欧美日韩乱码在线| 热re99久久精品国产66热6| 首页视频小说图片口味搜索| 国产91精品成人一区二区三区| netflix在线观看网站| 老鸭窝网址在线观看| 大型黄色视频在线免费观看| √禁漫天堂资源中文www| 女人久久www免费人成看片| 老司机在亚洲福利影院| 美女 人体艺术 gogo| av视频免费观看在线观看| 色综合婷婷激情| 最近最新免费中文字幕在线| 亚洲精品在线观看二区| 制服人妻中文乱码| 成在线人永久免费视频| 久久久精品区二区三区| 老司机午夜福利在线观看视频| 亚洲精品美女久久久久99蜜臀| 久久精品亚洲熟妇少妇任你| 99精品久久久久人妻精品| avwww免费| 丁香六月欧美| 1024视频免费在线观看| 国产亚洲精品第一综合不卡| 淫妇啪啪啪对白视频| 狠狠狠狠99中文字幕| 欧美日韩亚洲高清精品| 欧美国产精品一级二级三级| 国产一区二区三区在线臀色熟女 | 人人妻人人爽人人添夜夜欢视频| 超碰成人久久| 精品熟女少妇八av免费久了| 国产蜜桃级精品一区二区三区 | 亚洲一区二区三区欧美精品| av线在线观看网站| 精品第一国产精品| 亚洲七黄色美女视频| ponron亚洲| 成人18禁在线播放| 国产精品av久久久久免费| 中文亚洲av片在线观看爽 | 精品乱码久久久久久99久播| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩一卡2卡3卡4卡2021年| a级片在线免费高清观看视频| 欧美日韩国产mv在线观看视频| 91国产中文字幕| 啦啦啦免费观看视频1| 国产成人免费无遮挡视频| 欧美精品人与动牲交sv欧美| 99国产极品粉嫩在线观看| 久久国产精品人妻蜜桃| 黄色视频,在线免费观看| 国产无遮挡羞羞视频在线观看| 大片电影免费在线观看免费| 90打野战视频偷拍视频| 亚洲熟女毛片儿| 欧美一级毛片孕妇| videosex国产| 亚洲欧美激情在线| 777米奇影视久久| 久久精品人人爽人人爽视色| 久久香蕉精品热| 精品午夜福利视频在线观看一区| 久久久久久久午夜电影 | 人人妻人人爽人人添夜夜欢视频| 久久久久国产精品人妻aⅴ院 | 欧美黑人精品巨大| 搡老熟女国产l中国老女人| 波多野结衣一区麻豆| 久久精品国产亚洲av香蕉五月 | 超色免费av| 99re在线观看精品视频| 在线国产一区二区在线| 成年人黄色毛片网站| 久久精品熟女亚洲av麻豆精品| 亚洲av片天天在线观看| 国产色视频综合| 国产主播在线观看一区二区| 国产亚洲精品一区二区www | 国产xxxxx性猛交| 看黄色毛片网站| 午夜精品在线福利| 80岁老熟妇乱子伦牲交| 欧美人与性动交α欧美软件| 黄色片一级片一级黄色片| 欧美黑人欧美精品刺激| 国产精品一区二区精品视频观看|