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

    基于追逃博弈的非合作目標(biāo)接近控制

    2020-03-05 05:06:18羅建軍王明明
    宇航總體技術(shù) 2020年1期
    關(guān)鍵詞:航天器控制策略坐標(biāo)系

    柴 源,羅建軍,王明明,韓 楠

    (1. 西北工業(yè)大學(xué)航天飛行動(dòng)力學(xué)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 西安 710072;2. 西北工業(yè)大學(xué)青島研究院,青島 266200)

    0 引言

    空間自主交會(huì)是航天領(lǐng)域的一大研究熱點(diǎn)。隨著航天技術(shù)的不斷進(jìn)步,合作目標(biāo)的自主交會(huì)技術(shù)已經(jīng)比較成熟,并成功應(yīng)用于空間站、在軌服務(wù)等多種空間任務(wù)中。目前,空間非合作目標(biāo)如故障衛(wèi)星、空間垃圾等的自主接近引起了學(xué)者的廣泛關(guān)注[1]。

    非合作目標(biāo)由于其非合作性,具有如下特征:信息層面不溝通、機(jī)動(dòng)行為不配合、先驗(yàn)知識(shí)不完備等。因此在設(shè)計(jì)追蹤航天器的軌道控制方法時(shí)需要綜合考慮接近過(guò)程中的不確定性。針對(duì)非合作目標(biāo)接近問(wèn)題,目前已經(jīng)開(kāi)展了部分研究工作。根據(jù)建模方法的不同可將已有研究分為3種:基于視線坐標(biāo)系的建模與控制、基于目標(biāo)軌道坐標(biāo)系的建模與控制、基于追蹤航天器軌道坐標(biāo)系的建模與控制。在目標(biāo)軌道坐標(biāo)系下,王洪宇等[2]提出了一種全局魯棒最優(yōu)滑??刂破鱽?lái)克服非合作目標(biāo)所帶來(lái)的不確定性。但是由于目標(biāo)航天器的軌道半徑和速度無(wú)法直接得到,基于目標(biāo)坐標(biāo)系的建模具有局限性。在基于視線坐標(biāo)系的建模與控制方面,陳統(tǒng)等[3]建立了姿軌聯(lián)合運(yùn)動(dòng)模型,并結(jié)合具有魯棒性的模糊控制理論實(shí)現(xiàn)對(duì)非合作目標(biāo)的接近;殷澤陽(yáng)等[4]提出了低復(fù)雜度預(yù)設(shè)性能控制方法,實(shí)現(xiàn)在未知系統(tǒng)參數(shù)情況下的快速高精度目標(biāo)接近,但是基于視線坐標(biāo)系的建模為非線性模型,增加了控制求解難度。在追蹤航天器軌道坐標(biāo)系下,盧山等[5]設(shè)計(jì)了針對(duì)自主交會(huì)和攔截兩種接近模式的基于李雅普諾夫的控制律;郭永等[6]基于人工勢(shì)場(chǎng)法與蔓葉線理論的障礙物模型,提出了可以避障的滑模控制器。該坐標(biāo)系下,追蹤航天器可以基于自身的軌道信息及星載傳感器測(cè)量得到相對(duì)位置和速度信息等進(jìn)行控制器設(shè)計(jì),更加方便簡(jiǎn)潔。因此,本文采用基于追蹤航天器軌道坐標(biāo)系的相對(duì)運(yùn)動(dòng)模型,以便于控制律的設(shè)計(jì)。

    根據(jù)上述分析,多數(shù)設(shè)計(jì)方法都是通過(guò)提高控制器的魯棒性來(lái)克服非合作目標(biāo)的機(jī)動(dòng)以及外界干擾等。但是由于非合作目標(biāo)機(jī)動(dòng)上界的不確定,控制器的設(shè)計(jì)存在保守性,不利于燃料的優(yōu)化和接近精度的提高。

    博弈論研究的是多個(gè)參與者的最優(yōu)控制與決策問(wèn)題,其中每一位參與者通過(guò)各自目標(biāo)函數(shù)的優(yōu)化獲得控制策略[7]。近年來(lái),博弈控制方法在各種工程問(wèn)題的研究中也得到了應(yīng)用。Abouheaf等[8]、Lin[9]和Mylvaganam等[10]將博弈控制方法應(yīng)用到多智能體一致性、編隊(duì)和避障等問(wèn)題中。韓楠等[11]利用微分博弈實(shí)現(xiàn)了多顆微小衛(wèi)星對(duì)失效航天器的姿態(tài)接管控制。Innocenti等[12]利用基于狀態(tài)相關(guān)里卡提方程SDRE的非合作微分博弈控制實(shí)現(xiàn)交會(huì)任務(wù)。追逃博弈研究追捕者與逃逸者以不同的策略完成追捕任務(wù)的協(xié)調(diào)過(guò)程[13]。Bardhan等[14]基于追逃博弈設(shè)計(jì)了導(dǎo)彈攔截導(dǎo)引律,Li等[15]將近圓軌道上的兩個(gè)航天器追逃問(wèn)題轉(zhuǎn)化為兩點(diǎn)邊值優(yōu)化問(wèn)題進(jìn)行求解。因此,本文將非合作目標(biāo)視為理性的博弈參與者,設(shè)計(jì)追蹤航天器的追逃博弈控制方法,從而實(shí)現(xiàn)非合作目標(biāo)的精確接近。為了簡(jiǎn)化納什均衡的求解,追逃博弈模型選擇線性二次型微分博弈模型[16],以得到控制策略的顯式表達(dá)式,便于在線應(yīng)用。

    本文介紹了追蹤航天器的追逃博弈控制器的設(shè)計(jì)思路,基于追蹤航天器和非合作目標(biāo)的軌道相對(duì)運(yùn)動(dòng)模型,設(shè)計(jì)了與相對(duì)距離和燃耗有關(guān)的目標(biāo)函數(shù),并建立了二者的追逃博弈模型,推導(dǎo)了追逃博弈的均衡策略,并給出了策略求解算法,通過(guò)數(shù)值仿真驗(yàn)證了非合作目標(biāo)接近的追逃博弈控制方法的有效性。

    1 設(shè)計(jì)思路

    本文面向非合作目標(biāo)接近的軌道控制問(wèn)題,試圖提出一種基于線性二次型追逃博弈的控制方法。追逃博弈將參與雙方定義為追蹤者和逃脫者,在博弈過(guò)程中各博弈方均以己方最大利益為目標(biāo),一方的得益必然導(dǎo)致另一方的損失,二者的得失總和為0[14]。將非合作目標(biāo)接近問(wèn)題描述為追逃博弈問(wèn)題,其中追蹤航天器扮演追蹤者,非合作目標(biāo)扮演逃避者。要接近非合作目標(biāo),一方面,追蹤航天器要選擇其控制策略以調(diào)節(jié)到某種狀態(tài),并盡量減少博弈過(guò)程中的燃料消耗。另一方面,理性的非合作目標(biāo)試圖調(diào)節(jié)到使追蹤航天器難以追上的狀態(tài),選擇其控制策略的同時(shí)將自身燃料消耗降至最低。將二者互相沖突的目標(biāo)歸納為追逃博弈的目標(biāo)函數(shù)

    (1)

    式中,X為相對(duì)狀態(tài)量,具體含義在下文給出。u為追蹤者的控制量,v為逃逸者的控制量。Q>0,Ru>0,Rv>0均為對(duì)稱矩陣。追蹤航天器的目標(biāo)是最小化J,而非合作目標(biāo)則期望最大化J。

    在考慮二者動(dòng)力學(xué)約束的情況下,通過(guò)優(yōu)化二者的目標(biāo)函數(shù),建立追逃博弈模型:

    (2)

    通過(guò)建立在追蹤航天器上的軌道相對(duì)運(yùn)動(dòng)方程,將非合作目標(biāo)軌道接近的任務(wù)要求轉(zhuǎn)化為追逃博弈控制優(yōu)化問(wèn)題中的動(dòng)力學(xué)約束。通過(guò)最優(yōu)化問(wèn)題的求解得到納什均衡控制策略,追蹤航天器盡可能在燃耗最小的情況下實(shí)現(xiàn)非合作目標(biāo)的接近。

    2 非合作目標(biāo)接近的追逃博弈建模

    空間非合作目標(biāo)接近問(wèn)題涉及兩個(gè)近距離航天器間的軌道運(yùn)動(dòng),本節(jié)先給出追蹤航天器軌道坐標(biāo)系下追蹤航天器和非合作目標(biāo)的相對(duì)運(yùn)動(dòng)模型,之后建立二者的追逃博弈模型。

    2.1 相對(duì)運(yùn)動(dòng)建模

    本文中下標(biāo)e和p分別代指非合作目標(biāo)和追蹤航天器。在慣性坐標(biāo)系下,非合作目標(biāo)追蹤航天器的軌道運(yùn)動(dòng)方程分別為

    (3)

    式中,rp和re分別為追蹤航天器和非合作目標(biāo)在慣性坐標(biāo)系下的位置矢量;up和ue分別為追蹤航天器和非合作目標(biāo)的控制加速度;μ為地球引力常數(shù),μ=3.986×1014m3/s2。

    定義追蹤航天器和非合作目標(biāo)的相對(duì)位置為

    r=re-rp

    (4)

    則慣性坐標(biāo)系下的相對(duì)運(yùn)動(dòng)方程為

    (5)

    將式(5)投影在追蹤航天器本體坐標(biāo)系中可得

    (6)

    式中,ωe和ωp分別表示二者的軌道角速度,r表示慣性系下的位置矢量。

    在二者相對(duì)距離和非合作目標(biāo)地心距之比足夠小,即r?re的條件下,re=r+rp的2階及高階泰勒展開(kāi)項(xiàng)可忽略不計(jì),則相對(duì)軌道運(yùn)動(dòng)方程寫成狀態(tài)空間形式

    (7)

    其中

    其中

    式(7)為非合作目標(biāo)和追蹤航天器之間的相對(duì)運(yùn)動(dòng)方程。由于目標(biāo)航天器為非合作目標(biāo),無(wú)法得到其軌道信息,因此坐標(biāo)系建立在追蹤航天器上。追蹤航天器可以基于自身的軌道信息及星載傳感器測(cè)量得到相對(duì)位置和速度信息,以便進(jìn)行博弈問(wèn)題的建模和求解。

    2.2 追逃博弈建模

    追逃博弈由以下3個(gè)要素構(gòu)成:博弈參與者N={p,e}、各參與者容許策略集Ui、參與者目標(biāo)函數(shù)J[7]。為滿足非合作目標(biāo)接近的任務(wù)要求,設(shè)計(jì)如下目標(biāo)函數(shù)

    (8)

    其中,Q∈R6×6、Rp∈R3×3、Re∈R3×3為對(duì)稱正定的加權(quán)矩陣。

    非合作目標(biāo)和追蹤航天器進(jìn)行追逃博弈時(shí),二者通過(guò)獨(dú)立優(yōu)化各自目標(biāo)函數(shù)(8)來(lái)獲得控制策略。該策略稱為納什均衡,其定義如下:

    (9)

    其中Ui為參與者i的容許控制策略集。

    對(duì)應(yīng)于可行控制策略集ui∈Ui的值函數(shù)為

    (10)

    則非合作目標(biāo)接近的追逃博弈問(wèn)題可由下式描述

    (11)

    追蹤航天器通過(guò)求解上述優(yōu)化問(wèn)題,得到追逃博弈的納什均衡策略,以實(shí)現(xiàn)對(duì)非合作目標(biāo)的接近。

    3 非合作目標(biāo)接近的博弈控制策略

    本節(jié)先給出追逃博弈控制策略的狀態(tài)反饋表達(dá)式,之后給出李亞普諾夫迭代法進(jìn)行狀態(tài)反饋矩陣的求解。

    3.1 追逃博弈控制策略

    值函數(shù)的微分等價(jià)為

    (12)

    其中,V(0)=0。

    定義哈密爾頓函數(shù)為

    (13)

    對(duì)應(yīng)最優(yōu)值函數(shù)的反饋控制策略為

    (14)

    將其代哈密爾頓函數(shù)中可得HJ方程為

    (15)

    其中,V*(0)=0。

    假設(shè)最優(yōu)值函數(shù)在狀態(tài)x(t)下有線性二次型形式的解

    (16)

    則追蹤航天器和非合作目標(biāo)對(duì)應(yīng)的納什均衡反饋控制策略為

    (17)

    則HJ方程可以整理為

    (18)

    整理得

    (19)

    通過(guò)對(duì)上述代數(shù)黎卡提方程(19)進(jìn)行求解,可以得到對(duì)稱正定矩陣P,從而根據(jù)式(17)得到狀態(tài)反饋控制策略。

    本文控制策略與傳統(tǒng)的線性二次型調(diào)節(jié)器(LQR)方法有相似之處。LQR方法是現(xiàn)代控制理論中較成熟的一種狀態(tài)空間設(shè)計(jì)法,針對(duì)線性系統(tǒng),設(shè)計(jì)與系統(tǒng)狀態(tài)和控制輸入相關(guān)的二次型目標(biāo)函數(shù)

    (20)

    利用動(dòng)態(tài)規(guī)劃推導(dǎo)得到代數(shù)黎卡提方程

    (21)

    從而得到狀態(tài)反饋的最優(yōu)控制律

    (22)

    但是本文的控制策略是基于追逃博弈得到的,考慮最優(yōu)性的同時(shí),比傳統(tǒng)的LQR控制有更好的魯棒性。

    3.2 控制策略求解

    代數(shù)黎卡提方程(19)的求解已有豐富的研究成果[17],本文采用李雅普諾夫迭代法進(jìn)行計(jì)算。該方法將代數(shù)黎卡提方程解耦為李雅普諾夫方程來(lái)獨(dú)立運(yùn)算,算法速度快,準(zhǔn)確性高。

    迭代算法

    (A-SPP(k))TP(k+1)+P(k+1)(A-SPP(k))=
    -(Q+P(k+1)SpP(k+1)+P(k+1)SeP(k+1)),
    k=0,1,2,…

    (23)

    初值選擇

    0=ATP(0)+P(0)A+Q-P(0)SpP(0)

    (24)

    通過(guò)迭代求解李亞普諾夫方程式(20)和式(21)可以得到矩陣P。

    4 仿真算例及結(jié)果分析

    為了突出本控制器的優(yōu)勢(shì),本節(jié)將基于追逃博弈的控制方法與傳統(tǒng)LQR控制進(jìn)行對(duì)比,通過(guò)3組數(shù)值仿真算例驗(yàn)證基于追逃博弈的控制方法應(yīng)用于非合作目標(biāo)接近問(wèn)題的有效性。假設(shè)追蹤航天器初始時(shí)刻相對(duì)于非合作目標(biāo)的位置為r=[300,150,-100]Tm,追蹤航天器進(jìn)行非合作目標(biāo)逼近,最終二者的相對(duì)運(yùn)動(dòng)狀態(tài)為0。追蹤航天器的控制加速度幅值約束為umax=5m/s2。仿真軌道初始值如表1所示。

    表1 初始軌道參數(shù)

    算例1假設(shè)空間非合作目標(biāo)不存在機(jī)動(dòng)。該工況相對(duì)理想,用于檢驗(yàn)所提出控制方法的可行性。本文方法選擇權(quán)重矩陣為:Q=10-5I6,Rp=0.01I3,Re=0.02I3。LQR方法選擇權(quán)重矩陣為:Q=10-5I6,R=0.01I3。仿真時(shí)間為200s,仿真步長(zhǎng)為0.1s。

    圖1和圖3分別為本文提出的方法在接近過(guò)程中,非合作目標(biāo)和追蹤航天器相對(duì)距離和相對(duì)速度隨時(shí)間的變化曲線。圖2和圖4分別為L(zhǎng)QR提出的方法在接近過(guò)程中,非合作目標(biāo)和追蹤航天器相對(duì)距離和相對(duì)速度隨時(shí)間的變化曲線。經(jīng)過(guò)約40s,兩種控制器均使追蹤航天器與非合作目標(biāo)的相對(duì)距離穩(wěn)定在0.5m左右。

    圖5和圖6分別為兩種控制器下追蹤航天器的控制加速度隨時(shí)間變化曲線??梢钥闯觯谡麄€(gè)非合作目標(biāo)接近過(guò)程中,初始相對(duì)距離較遠(yuǎn),接近非合作目標(biāo)所需控制力較大,隨著相對(duì)距離的減小,控制力逐漸減少并趨于0。通過(guò)上述分析,在非合作目標(biāo)無(wú)機(jī)動(dòng)的理想情況下,兩種方法均可實(shí)現(xiàn)對(duì)非合作目標(biāo)的接近。

    算例2假設(shè)非合作目標(biāo)的未知機(jī)動(dòng)為納什均衡策略。該工況下,非合作目標(biāo)為理性的博弈參與者,有意識(shí)地與追蹤航天器對(duì)抗。假設(shè)非合作目標(biāo)的最大控制加速度umax=2m/s2。本文方法選擇加權(quán)矩陣為:Q=10-5I6,Rp=0.01I3,Re=0.02I3。LQR方法選擇目標(biāo)函數(shù)中的矩陣為:Q=10-5I6,Rp=0.01I3。仿真時(shí)間為200s,仿真步長(zhǎng)為0.1s。

    圖1 相對(duì)距離隨時(shí)間變化曲線(本文)Fig.1 Relative distance by game

    圖2 相對(duì)距離隨時(shí)間變化曲線(LQR)Fig.2 Relative distance by LQR

    圖3 相對(duì)速度隨時(shí)間變化曲線(本文)Fig.3 Relative velocity by game

    圖4 相對(duì)速度隨時(shí)間變化曲線(LQR)Fig.4 Relative velocity by LQR

    圖5 控制加速度隨時(shí)間變化曲線(本文)Fig.5 Control acceleration by game

    圖6 控制加速度隨時(shí)間變化曲線(LQR)Fig.6 Control acceleration by LQR

    圖7和圖8為追逃博弈方法與LQR方法分別作用下的相對(duì)位置變化圖。圖9和圖10為兩種控制分別作用下的相對(duì)速度變化圖??梢钥闯觯诜呛献髂繕?biāo)采取納什均衡策略時(shí),本文所提出的方法能夠快速平滑地使相對(duì)距離收斂到0.5m左右的穩(wěn)定值。而LQR方法則是震蕩收斂狀態(tài),所需時(shí)間較長(zhǎng)。

    圖11和圖12分別為兩種控制器下追蹤航天器的控制加速度隨時(shí)間變化曲線??梢钥闯?,基于追逃博弈的控制方法可以在燃耗較少的情況下快速收斂到0。對(duì)比二者的目標(biāo)函數(shù),在二者都采取納什均衡策略的情況下,即基于追逃博弈的控制下,J*=104;而在LQR控制下,J*=141,由此也可以驗(yàn)證式(9)的右不等式成立。

    算例3假設(shè)非合作目標(biāo)存在未知機(jī)動(dòng)[4]:

    本文方法選擇權(quán)重矩陣為:Q=10-5I6,Rp=0.01I3,Re=0.008I3。LQR方法選擇權(quán)重矩陣為:Q=10-5I6,Rp=0.01I3。仿真時(shí)間為200s,仿真步長(zhǎng)為0.1s。

    在本工況下,圖13和圖14為追逃博弈方法與LQR方法分別作用下的相對(duì)位置變化圖。圖15和圖16為兩種控制方法下的相對(duì)速度變化圖。可以看出,盡管非合作目標(biāo)存在未知機(jī)動(dòng),追逃博弈的控制方法仍可以實(shí)現(xiàn)狀態(tài)的收斂,精度在1m左右。而LQR方法魯棒性不足,無(wú)法實(shí)現(xiàn)非合作目標(biāo)的接近。

    圖17和圖18分別為兩種控制器下追蹤航天器的控制加速度隨時(shí)間變化曲線??梢钥闯觯刂萍铀俣瘸掷m(xù)并不為0,而是隨著非合作目標(biāo)的運(yùn)動(dòng)震蕩。

    圖7 相對(duì)距離隨時(shí)間變化曲線(本文)Fig.7 Relative distance by game

    圖8 相對(duì)距離隨時(shí)間變化曲線(LQR)Fig.8 Relative distance by LQR

    圖9 相對(duì)速度隨時(shí)間變化曲線(本文)Fig.9 Relative velocity by game

    圖10 相對(duì)速度隨時(shí)間變化曲線(LQR)Fig.10 Relative velocity by LQR

    圖11 控制加速度隨時(shí)間變化曲線(本文)Fig.11 Control acceleration by game

    圖12 控制加速度隨時(shí)間變化曲線(LQR)Fig.12 Control acceleration by LQR

    圖13 相對(duì)距離隨時(shí)間變化曲線(本文)Fig.13 Relative distance by game

    圖14 相對(duì)距離隨時(shí)間變化曲線(LQR)Fig.14 Relative distance by LQR

    圖15 相對(duì)速度隨時(shí)間變化曲線(本文)Fig.15 Relative velocity by game

    圖16 相對(duì)速度隨時(shí)間變化曲線(LQR)Fig.16 Relative velocity by LQR

    圖17 控制加速度隨時(shí)間變化曲線(本文)Fig.17 Control acceleration by game

    圖18 控制加速度隨時(shí)間變化曲線(LQR)Fig.18 Control acceleration by LQR

    5 結(jié)論

    本文針對(duì)空間非合作目標(biāo)的接近控制問(wèn)題,基于追逃博弈方法設(shè)計(jì)了追蹤航天器的軌道控制器。面向非合作目標(biāo)接近的任務(wù)要求,合理設(shè)計(jì)了博弈的目標(biāo)函數(shù),并結(jié)合二者的動(dòng)力學(xué)約束,實(shí)現(xiàn)了對(duì)非合作目標(biāo)和追蹤航天器之間追逃博弈的數(shù)學(xué)描述。結(jié)合線性化動(dòng)力學(xué),通過(guò)優(yōu)化二次型目標(biāo)函數(shù),得到線性二次型追逃博弈的納什均衡解策略?;谧诽硬┺牡目刂撇呗跃哂酗@式表達(dá)式,方便工程應(yīng)用。數(shù)值仿真驗(yàn)證了本文設(shè)計(jì)的追逃博弈控制方法對(duì)于存在未知機(jī)動(dòng)的非合作目標(biāo)的有效性。本文未考慮接近過(guò)程中的姿態(tài)運(yùn)動(dòng),后續(xù)研究將進(jìn)一步考慮能夠?qū)崿F(xiàn)非合作目標(biāo)接近的姿軌聯(lián)合博弈控制。

    猜你喜歡
    航天器控制策略坐標(biāo)系
    2022 年第二季度航天器發(fā)射統(tǒng)計(jì)
    考慮虛擬慣性的VSC-MTDC改進(jìn)下垂控制策略
    能源工程(2020年6期)2021-01-26 00:55:22
    2019 年第二季度航天器發(fā)射統(tǒng)計(jì)
    工程造價(jià)控制策略
    山東冶金(2019年3期)2019-07-10 00:54:04
    2018 年第三季度航天器發(fā)射統(tǒng)計(jì)
    2018年第二季度航天器發(fā)射統(tǒng)計(jì)
    現(xiàn)代企業(yè)會(huì)計(jì)的內(nèi)部控制策略探討
    解密坐標(biāo)系中的平移變換
    坐標(biāo)系背后的故事
    基于重心坐標(biāo)系的平面幾何證明的探討
    深爱激情五月婷婷| 国产免费一级a男人的天堂| 国产真实乱freesex| 久久精品久久久久久噜噜老黄 | 免费看日本二区| 亚洲精品影视一区二区三区av| 国产主播在线观看一区二区| 久久久久久久久久成人| 久久久久精品国产欧美久久久| 在线观看美女被高潮喷水网站| 国产老妇女一区| av黄色大香蕉| 能在线免费观看的黄片| 我要看日韩黄色一级片| 无遮挡黄片免费观看| 成年版毛片免费区| 麻豆av噜噜一区二区三区| av女优亚洲男人天堂| 国产精品,欧美在线| 99国产极品粉嫩在线观看| 国产免费男女视频| 天堂影院成人在线观看| 99精品久久久久人妻精品| 国产高清有码在线观看视频| 亚洲 国产 在线| 精品久久久久久久久久久久久| 久久中文看片网| 成人特级黄色片久久久久久久| 日日夜夜操网爽| 极品教师在线视频| 久久久久九九精品影院| 亚洲av免费在线观看| 日韩欧美国产一区二区入口| 一个人看视频在线观看www免费| 午夜精品在线福利| 久久人人精品亚洲av| 少妇熟女aⅴ在线视频| 久久精品国产亚洲av天美| 毛片一级片免费看久久久久 | 亚洲av免费在线观看| 99视频精品全部免费 在线| 国产成人a区在线观看| 深爱激情五月婷婷| 欧美+亚洲+日韩+国产| 韩国av一区二区三区四区| 亚洲美女搞黄在线观看 | 一级毛片久久久久久久久女| 久久婷婷人人爽人人干人人爱| 精品一区二区三区视频在线观看免费| 亚洲美女搞黄在线观看 | 九九久久精品国产亚洲av麻豆| 全区人妻精品视频| 天堂av国产一区二区熟女人妻| 日韩 亚洲 欧美在线| 69av精品久久久久久| 综合色av麻豆| 国产av麻豆久久久久久久| www.www免费av| 亚洲熟妇中文字幕五十中出| 亚洲第一区二区三区不卡| 午夜爱爱视频在线播放| 免费看美女性在线毛片视频| 在线免费观看的www视频| 色在线成人网| 亚洲中文字幕日韩| 91狼人影院| 精品一区二区三区av网在线观看| 一进一出抽搐动态| 日韩av在线大香蕉| 国产精品一区二区性色av| 国产av一区在线观看免费| 精品欧美国产一区二区三| 听说在线观看完整版免费高清| 亚洲av.av天堂| 国产三级在线视频| av.在线天堂| 亚洲av中文字字幕乱码综合| 午夜爱爱视频在线播放| 97超视频在线观看视频| 免费看a级黄色片| 一区二区三区激情视频| 伊人久久精品亚洲午夜| 亚洲三级黄色毛片| 一卡2卡三卡四卡精品乱码亚洲| 性色avwww在线观看| 精品乱码久久久久久99久播| 内射极品少妇av片p| 亚洲男人的天堂狠狠| 久久久国产成人免费| 亚洲精华国产精华精| 日本精品一区二区三区蜜桃| 搡老岳熟女国产| 97碰自拍视频| 两个人的视频大全免费| 日韩欧美在线乱码| 两个人视频免费观看高清| 国产男靠女视频免费网站| 在线观看66精品国产| 可以在线观看毛片的网站| 蜜桃亚洲精品一区二区三区| 国产精品嫩草影院av在线观看 | 国产精品久久视频播放| 亚洲av日韩精品久久久久久密| 欧美日本亚洲视频在线播放| 男女啪啪激烈高潮av片| 婷婷亚洲欧美| 人人妻人人看人人澡| 精品人妻一区二区三区麻豆 | 97人妻精品一区二区三区麻豆| 中国美白少妇内射xxxbb| 三级毛片av免费| 黄色视频,在线免费观看| 国产男人的电影天堂91| 国产不卡一卡二| 久久人人精品亚洲av| 悠悠久久av| 午夜福利在线观看免费完整高清在 | 国产免费av片在线观看野外av| 熟女电影av网| 午夜日韩欧美国产| 亚洲图色成人| 如何舔出高潮| 最新在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 一本久久中文字幕| 91久久精品国产一区二区成人| 久久久成人免费电影| 人妻制服诱惑在线中文字幕| 免费观看精品视频网站| 久久人妻av系列| 日本在线视频免费播放| 欧美丝袜亚洲另类 | 国内精品一区二区在线观看| 不卡视频在线观看欧美| 国产在线男女| 精品99又大又爽又粗少妇毛片 | 九九在线视频观看精品| 日韩欧美国产在线观看| 亚洲国产精品久久男人天堂| 成年女人毛片免费观看观看9| 成人鲁丝片一二三区免费| 成人鲁丝片一二三区免费| 亚洲人与动物交配视频| 男人的好看免费观看在线视频| 我要搜黄色片| 久久久久久久久大av| 久久精品夜夜夜夜夜久久蜜豆| 国产男靠女视频免费网站| 老熟妇仑乱视频hdxx| 99久久久亚洲精品蜜臀av| 亚洲第一区二区三区不卡| 99视频精品全部免费 在线| 动漫黄色视频在线观看| 日本熟妇午夜| 国产精品野战在线观看| 无人区码免费观看不卡| 1024手机看黄色片| a在线观看视频网站| 日韩av在线大香蕉| 久久亚洲真实| 99国产精品一区二区蜜桃av| 久久6这里有精品| a在线观看视频网站| 免费观看人在逋| 国产精品久久久久久精品电影| 成年女人看的毛片在线观看| 三级男女做爰猛烈吃奶摸视频| 精品福利观看| 精品福利观看| 亚洲欧美日韩高清在线视频| 国产精品伦人一区二区| 亚洲欧美日韩卡通动漫| 免费人成视频x8x8入口观看| 亚洲精品国产成人久久av| 久久亚洲精品不卡| 久久热精品热| 亚洲第一区二区三区不卡| 国产视频一区二区在线看| 国产乱人伦免费视频| 久久久久久久久久成人| 看黄色毛片网站| 成人欧美大片| 在线播放无遮挡| 国产av麻豆久久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 午夜爱爱视频在线播放| 日本一二三区视频观看| 午夜免费激情av| 国产精华一区二区三区| 一本一本综合久久| 成人性生交大片免费视频hd| 熟女人妻精品中文字幕| 一级黄片播放器| 在线看三级毛片| 两人在一起打扑克的视频| 99在线视频只有这里精品首页| 久久6这里有精品| 夜夜爽天天搞| 国产精品久久视频播放| 久久久久九九精品影院| 在现免费观看毛片| 精品欧美国产一区二区三| 国产淫片久久久久久久久| 韩国av在线不卡| 亚洲欧美激情综合另类| 久久午夜亚洲精品久久| 欧美激情在线99| 欧美一区二区亚洲| 亚洲欧美日韩卡通动漫| a在线观看视频网站| 国产一区二区三区在线臀色熟女| 亚洲无线在线观看| 极品教师在线免费播放| 女人十人毛片免费观看3o分钟| 欧美zozozo另类| 日韩强制内射视频| 婷婷六月久久综合丁香| 亚洲欧美日韩无卡精品| 99久久精品热视频| 欧美精品国产亚洲| or卡值多少钱| 99国产精品一区二区蜜桃av| 精品久久久久久久末码| 亚洲国产精品成人综合色| a在线观看视频网站| 国模一区二区三区四区视频| av在线天堂中文字幕| 99热只有精品国产| 天堂动漫精品| 国内精品宾馆在线| 日本爱情动作片www.在线观看 | 我的老师免费观看完整版| 国产精品久久久久久亚洲av鲁大| 嫩草影视91久久| 久久精品国产自在天天线| 国产精品人妻久久久影院| 精品久久国产蜜桃| 噜噜噜噜噜久久久久久91| 国产成人福利小说| 欧美精品国产亚洲| 听说在线观看完整版免费高清| 99国产精品一区二区蜜桃av| 亚洲最大成人手机在线| 久久精品国产亚洲网站| 亚洲男人的天堂狠狠| 久久久久久久亚洲中文字幕| 亚洲av中文字字幕乱码综合| 两性午夜刺激爽爽歪歪视频在线观看| 老司机深夜福利视频在线观看| 91狼人影院| 亚洲最大成人中文| 国产毛片a区久久久久| 国产老妇女一区| 大型黄色视频在线免费观看| 男人舔女人下体高潮全视频| 国产欧美日韩精品亚洲av| av在线天堂中文字幕| 一边摸一边抽搐一进一小说| 在线免费观看不下载黄p国产 | 精品福利观看| 亚洲国产精品合色在线| 91狼人影院| h日本视频在线播放| 一区二区三区高清视频在线| 五月伊人婷婷丁香| 日韩欧美一区二区三区在线观看| 久久久久久久午夜电影| 亚洲aⅴ乱码一区二区在线播放| 一本久久中文字幕| 麻豆av噜噜一区二区三区| 久久精品国产亚洲av天美| 可以在线观看毛片的网站| 窝窝影院91人妻| 真人做人爱边吃奶动态| 一个人看视频在线观看www免费| 在线观看免费视频日本深夜| 日本五十路高清| 久久精品国产清高在天天线| 露出奶头的视频| 国产精品久久久久久av不卡| 成年人黄色毛片网站| 91麻豆av在线| 十八禁国产超污无遮挡网站| 日本黄大片高清| 91久久精品国产一区二区三区| 热99re8久久精品国产| 亚洲无线观看免费| 精华霜和精华液先用哪个| 亚洲性夜色夜夜综合| 日日干狠狠操夜夜爽| 国模一区二区三区四区视频| 久99久视频精品免费| 亚洲久久久久久中文字幕| 日韩一区二区视频免费看| 变态另类成人亚洲欧美熟女| 欧美日本亚洲视频在线播放| 长腿黑丝高跟| 最近中文字幕高清免费大全6 | 亚洲av不卡在线观看| 国产精品不卡视频一区二区| 日韩国内少妇激情av| 三级男女做爰猛烈吃奶摸视频| 国产高清不卡午夜福利| 成人特级黄色片久久久久久久| 成人二区视频| 观看免费一级毛片| av在线亚洲专区| 国产高潮美女av| 国产精品国产三级国产av玫瑰| 床上黄色一级片| 免费大片18禁| 亚洲av第一区精品v没综合| 中亚洲国语对白在线视频| aaaaa片日本免费| 成人综合一区亚洲| av天堂中文字幕网| 美女大奶头视频| 国产亚洲91精品色在线| 制服丝袜大香蕉在线| 精品一区二区三区视频在线| 在线播放国产精品三级| 免费av不卡在线播放| 波多野结衣高清无吗| 婷婷六月久久综合丁香| 免费不卡的大黄色大毛片视频在线观看 | 国产av在哪里看| 久久久久久伊人网av| 校园人妻丝袜中文字幕| 中文字幕精品亚洲无线码一区| 人妻丰满熟妇av一区二区三区| 国产精品亚洲一级av第二区| 亚洲中文字幕一区二区三区有码在线看| 国产一区二区亚洲精品在线观看| 少妇猛男粗大的猛烈进出视频 | 哪里可以看免费的av片| 99久久无色码亚洲精品果冻| 一级av片app| 国产伦在线观看视频一区| 色哟哟哟哟哟哟| 在线免费观看的www视频| 性插视频无遮挡在线免费观看| 免费看av在线观看网站| 天堂av国产一区二区熟女人妻| 亚洲最大成人中文| 亚洲中文字幕一区二区三区有码在线看| 99国产极品粉嫩在线观看| 日韩中字成人| 国产成人福利小说| 午夜精品在线福利| 国产大屁股一区二区在线视频| 欧美绝顶高潮抽搐喷水| 亚洲成人精品中文字幕电影| 听说在线观看完整版免费高清| 久9热在线精品视频| 国产在线精品亚洲第一网站| 久久精品91蜜桃| 欧美bdsm另类| 免费大片18禁| 亚洲精品影视一区二区三区av| 狂野欧美白嫩少妇大欣赏| 亚洲男人的天堂狠狠| 久久精品国产亚洲av涩爱 | 制服丝袜大香蕉在线| 别揉我奶头~嗯~啊~动态视频| 久久国产精品人妻蜜桃| 麻豆成人av在线观看| 久久精品国产清高在天天线| 性插视频无遮挡在线免费观看| 好男人在线观看高清免费视频| 老师上课跳d突然被开到最大视频| 久久精品国产亚洲av涩爱 | 最近视频中文字幕2019在线8| av女优亚洲男人天堂| 日韩大尺度精品在线看网址| 中文亚洲av片在线观看爽| 国产私拍福利视频在线观看| 97热精品久久久久久| 内地一区二区视频在线| 少妇高潮的动态图| 一级a爱片免费观看的视频| 99九九线精品视频在线观看视频| 99riav亚洲国产免费| 日本三级黄在线观看| 国产探花极品一区二区| 亚洲精品日韩av片在线观看| 搡老妇女老女人老熟妇| 中出人妻视频一区二区| 欧美+日韩+精品| 少妇人妻一区二区三区视频| 国产伦人伦偷精品视频| 婷婷丁香在线五月| 又黄又爽又免费观看的视频| 少妇被粗大猛烈的视频| 69av精品久久久久久| 婷婷丁香在线五月| 在线播放无遮挡| 国产爱豆传媒在线观看| 亚洲自偷自拍三级| 麻豆成人午夜福利视频| 啦啦啦韩国在线观看视频| 悠悠久久av| 国产一区二区亚洲精品在线观看| 国产成人aa在线观看| 欧美高清成人免费视频www| 最好的美女福利视频网| 嫩草影院新地址| 久久久国产成人精品二区| 免费无遮挡裸体视频| 搡女人真爽免费视频火全软件 | 亚洲av二区三区四区| 色噜噜av男人的天堂激情| 欧美成人a在线观看| 日韩欧美国产在线观看| 国产午夜福利久久久久久| 欧美精品国产亚洲| 在现免费观看毛片| 国产伦精品一区二区三区四那| 日韩欧美国产一区二区入口| 噜噜噜噜噜久久久久久91| 干丝袜人妻中文字幕| 亚洲色图av天堂| 他把我摸到了高潮在线观看| 国模一区二区三区四区视频| 国产v大片淫在线免费观看| 欧美潮喷喷水| 22中文网久久字幕| 听说在线观看完整版免费高清| 久久草成人影院| 精品福利观看| 九九久久精品国产亚洲av麻豆| 成年人黄色毛片网站| 免费人成在线观看视频色| 日韩欧美国产一区二区入口| 日本爱情动作片www.在线观看 | 亚洲国产色片| 国产精品一区二区免费欧美| 三级男女做爰猛烈吃奶摸视频| 香蕉av资源在线| 给我免费播放毛片高清在线观看| 日日夜夜操网爽| 国产乱人视频| 日本爱情动作片www.在线观看 | 成年免费大片在线观看| 色综合色国产| 国产成年人精品一区二区| 国产成人aa在线观看| 欧美日韩国产亚洲二区| eeuss影院久久| 一边摸一边抽搐一进一小说| 一本精品99久久精品77| 亚洲精品在线观看二区| 色在线成人网| 免费看美女性在线毛片视频| 欧美高清性xxxxhd video| 欧美三级亚洲精品| 国产精品三级大全| 中出人妻视频一区二区| 变态另类成人亚洲欧美熟女| 直男gayav资源| 身体一侧抽搐| 久久久久久久午夜电影| 成熟少妇高潮喷水视频| 欧美日韩黄片免| 最近视频中文字幕2019在线8| 亚洲国产精品久久男人天堂| 人妻夜夜爽99麻豆av| 国产黄色小视频在线观看| av黄色大香蕉| 久久久午夜欧美精品| 亚洲黑人精品在线| 大又大粗又爽又黄少妇毛片口| 国产一区二区三区av在线 | 在线国产一区二区在线| 日本 av在线| 能在线免费观看的黄片| 欧美不卡视频在线免费观看| 国产精品日韩av在线免费观看| 女的被弄到高潮叫床怎么办 | 欧美激情在线99| 日本一本二区三区精品| 欧美+亚洲+日韩+国产| 亚洲av熟女| 免费观看精品视频网站| 一级a爱片免费观看的视频| 国内精品宾馆在线| 亚洲狠狠婷婷综合久久图片| 精品人妻一区二区三区麻豆 | 国产精品一区二区免费欧美| 黄片wwwwww| 亚洲av不卡在线观看| 人人妻人人澡欧美一区二区| 狂野欧美白嫩少妇大欣赏| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| av在线亚洲专区| 免费av毛片视频| 日本黄色片子视频| 午夜福利在线观看免费完整高清在 | 日韩欧美免费精品| 国产av一区在线观看免费| 日日啪夜夜撸| 成年人黄色毛片网站| 人妻制服诱惑在线中文字幕| 国产精品人妻久久久影院| 亚洲av电影不卡..在线观看| 大型黄色视频在线免费观看| 97热精品久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 国产成年人精品一区二区| 免费看日本二区| 日本免费一区二区三区高清不卡| 最近中文字幕高清免费大全6 | 黄色一级大片看看| 久久久久久国产a免费观看| 直男gayav资源| 免费看a级黄色片| 在现免费观看毛片| 成年版毛片免费区| 国产精品自产拍在线观看55亚洲| 久久久久久久亚洲中文字幕| 婷婷色综合大香蕉| 少妇人妻精品综合一区二区 | 国产不卡一卡二| 国产中年淑女户外野战色| 夜夜看夜夜爽夜夜摸| 免费看美女性在线毛片视频| 亚洲在线自拍视频| 男女下面进入的视频免费午夜| 国内精品久久久久久久电影| 国产精品久久久久久久电影| 欧美激情在线99| 搡老熟女国产l中国老女人| 国产国拍精品亚洲av在线观看| 久久久久久久久久久丰满 | 哪里可以看免费的av片| 久久精品国产鲁丝片午夜精品 | 国产av不卡久久| 欧美极品一区二区三区四区| 日韩一本色道免费dvd| 俺也久久电影网| 不卡视频在线观看欧美| 狂野欧美激情性xxxx在线观看| 真实男女啪啪啪动态图| 国内揄拍国产精品人妻在线| 久久久国产成人精品二区| 欧美+亚洲+日韩+国产| 午夜精品在线福利| 女同久久另类99精品国产91| 一区二区三区四区激情视频 | 亚洲自偷自拍三级| 亚洲美女搞黄在线观看 | 欧美日韩黄片免| 亚洲精华国产精华精| 老司机深夜福利视频在线观看| 成人特级av手机在线观看| 又黄又爽又刺激的免费视频.| 国产精品久久久久久亚洲av鲁大| 两个人视频免费观看高清| 91在线观看av| 久久久久久国产a免费观看| 欧美一级a爱片免费观看看| 日日撸夜夜添| 国产高清三级在线| 久久久精品大字幕| 日韩欧美一区二区三区在线观看| 国产久久久一区二区三区| 九九久久精品国产亚洲av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 久久人人精品亚洲av| 一级黄片播放器| 日韩人妻高清精品专区| 久久九九热精品免费| 麻豆成人av在线观看| 亚洲欧美日韩卡通动漫| 国产精品亚洲美女久久久| 狂野欧美激情性xxxx在线观看| 成人一区二区视频在线观看| 亚洲自偷自拍三级| 国产在线精品亚洲第一网站| 国产精品一及| 天堂影院成人在线观看| 久久久成人免费电影| 欧美高清性xxxxhd video| 少妇被粗大猛烈的视频| 人妻久久中文字幕网| 亚洲精品成人久久久久久| 欧美成人免费av一区二区三区| 日本色播在线视频| 女人被狂操c到高潮| 国产 一区精品| 啦啦啦啦在线视频资源| 色在线成人网| 麻豆av噜噜一区二区三区| 嫩草影视91久久| 国产成人av教育| 日本免费一区二区三区高清不卡| 麻豆国产av国片精品| 身体一侧抽搐| 国产欧美日韩一区二区精品| 国产午夜福利久久久久久| 男人和女人高潮做爰伦理| 国产精品一区二区性色av| 一本一本综合久久| 身体一侧抽搐| 99精品久久久久人妻精品| 欧美一区二区精品小视频在线| 能在线免费观看的黄片| 国内少妇人妻偷人精品xxx网站| 99精品久久久久人妻精品| 在线观看av片永久免费下载| 久久精品国产亚洲网站| 神马国产精品三级电影在线观看| 欧美xxxx性猛交bbbb| 国产女主播在线喷水免费视频网站 | 搡老妇女老女人老熟妇|