• <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一片| 一本综合久久免费| 母亲3免费完整高清在线观看| 女人被躁到高潮嗷嗷叫费观| 高清欧美精品videossex| 在线av久久热| 国产精品一区二区在线观看99| 亚洲专区中文字幕在线| 99九九在线精品视频| 高清黄色对白视频在线免费看| 亚洲精品久久成人aⅴ小说| 成人国产一区最新在线观看| 欧美日韩av久久| 香蕉久久夜色| 国产免费视频播放在线视频| 欧美成狂野欧美在线观看| 久久久久久久精品吃奶| 日韩三级视频一区二区三区| 欧美黑人欧美精品刺激| 天天躁狠狠躁夜夜躁狠狠躁| 国产视频一区二区在线看| 99国产精品一区二区三区| 欧美国产精品一级二级三级| 亚洲avbb在线观看| 国产单亲对白刺激| 久久99一区二区三区| 十八禁网站免费在线| 一个人免费看片子| 成人精品一区二区免费| 搡老乐熟女国产| 午夜精品久久久久久毛片777| 在线看a的网站| 国产三级黄色录像| 9191精品国产免费久久| 中文字幕色久视频| a在线观看视频网站| 欧美精品一区二区大全| 岛国毛片在线播放| 日本欧美视频一区| 久久久久网色| 一本综合久久免费| 成人精品一区二区免费| 中文字幕最新亚洲高清| 亚洲专区中文字幕在线| 免费av中文字幕在线| 十八禁人妻一区二区| 99国产综合亚洲精品| 亚洲欧美激情在线| 99精品久久久久人妻精品| 男女边摸边吃奶| 欧美日本中文国产一区发布| 国产又爽黄色视频| 日本五十路高清| 久久久久久久国产电影| 久久精品亚洲熟妇少妇任你| 国产成人免费无遮挡视频| 久久久国产一区二区| 天堂8中文在线网| 免费女性裸体啪啪无遮挡网站| 美女高潮喷水抽搐中文字幕| 国产xxxxx性猛交| 12—13女人毛片做爰片一| 午夜福利影视在线免费观看| 国产高清国产精品国产三级| 极品人妻少妇av视频| 天天操日日干夜夜撸| 国产男女超爽视频在线观看| 日韩欧美一区二区三区在线观看 | 久久精品国产亚洲av高清一级| 宅男免费午夜| 色在线成人网| 国产亚洲一区二区精品| 亚洲自偷自拍图片 自拍| 中文字幕最新亚洲高清| 久久青草综合色| 色婷婷av一区二区三区视频| 中文字幕人妻熟女乱码| 嫁个100分男人电影在线观看| 91av网站免费观看| 国产一卡二卡三卡精品| 欧美人与性动交α欧美软件| 亚洲精品久久午夜乱码| 在线播放国产精品三级| 2018国产大陆天天弄谢| 午夜福利在线观看吧| 精品国产一区二区三区久久久樱花| 久久精品国产亚洲av香蕉五月 | 精品一区二区三区视频在线观看免费 | 18禁美女被吸乳视频| 亚洲午夜理论影院| 久久久久久亚洲精品国产蜜桃av| 如日韩欧美国产精品一区二区三区| 肉色欧美久久久久久久蜜桃| 美女高潮到喷水免费观看| 捣出白浆h1v1| 国产xxxxx性猛交| 午夜视频精品福利| 黄色视频在线播放观看不卡| 中文字幕最新亚洲高清| 国产男女内射视频| 麻豆国产av国片精品| 丝袜人妻中文字幕| 亚洲成a人片在线一区二区| 久久精品国产a三级三级三级| 欧美性长视频在线观看| 欧美人与性动交α欧美软件| 亚洲中文日韩欧美视频| 日韩视频在线欧美| 免费黄频网站在线观看国产| 狠狠精品人妻久久久久久综合| 如日韩欧美国产精品一区二区三区| av有码第一页| 99在线人妻在线中文字幕 | 成年人午夜在线观看视频| 啦啦啦在线免费观看视频4| 丰满迷人的少妇在线观看| 99久久精品国产亚洲精品| 人成视频在线观看免费观看| cao死你这个sao货| 亚洲七黄色美女视频| 国产成人精品在线电影| 视频区欧美日本亚洲| 天堂动漫精品| 久久久水蜜桃国产精品网| 日韩欧美一区二区三区在线观看 | 最新在线观看一区二区三区| 成人国语在线视频| 久久久久视频综合| 久久久精品国产亚洲av高清涩受| 国产精品偷伦视频观看了| 丝袜美腿诱惑在线| 亚洲人成77777在线视频| 丝袜美腿诱惑在线| 91av网站免费观看| 午夜福利免费观看在线| 久久久欧美国产精品| 免费在线观看完整版高清| 亚洲国产看品久久| 精品国产一区二区三区四区第35| 午夜福利视频在线观看免费| 涩涩av久久男人的天堂| 免费看a级黄色片| 国产男靠女视频免费网站| 999久久久国产精品视频| 在线观看一区二区三区激情| 日韩 欧美 亚洲 中文字幕| 韩国精品一区二区三区| 天堂俺去俺来也www色官网| 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美激情在线| 精品人妻1区二区| 国产不卡一卡二| 国产免费av片在线观看野外av| 1024香蕉在线观看| 精品乱码久久久久久99久播| 国产精品九九99| 亚洲熟妇熟女久久| 精品免费久久久久久久清纯 | 久久亚洲精品不卡| 99精品欧美一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| kizo精华| 久久毛片免费看一区二区三区| 亚洲精品美女久久久久99蜜臀| 免费在线观看完整版高清| 90打野战视频偷拍视频| 亚洲九九香蕉| 69av精品久久久久久 | 高清在线国产一区| 中文字幕最新亚洲高清| 成人18禁高潮啪啪吃奶动态图| 色尼玛亚洲综合影院| 最新美女视频免费是黄的| 亚洲欧洲精品一区二区精品久久久| 国产精品免费大片| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕色久视频| 国产亚洲午夜精品一区二区久久| 欧美成人免费av一区二区三区 | 国产亚洲欧美在线一区二区| 国产三级黄色录像| 黑丝袜美女国产一区| 国产伦人伦偷精品视频| 在线看a的网站| 99riav亚洲国产免费| 欧美乱妇无乱码| 国产精品亚洲av一区麻豆| 亚洲精品国产一区二区精华液| 久久久久久久久免费视频了| 欧美人与性动交α欧美软件| 黑人巨大精品欧美一区二区mp4| 下体分泌物呈黄色| 日韩一区二区三区影片| 精品人妻1区二区| 国产又爽黄色视频| 精品人妻在线不人妻| 美女高潮喷水抽搐中文字幕| 午夜福利免费观看在线| 纵有疾风起免费观看全集完整版| 久久亚洲真实| 久久 成人 亚洲| 久久中文看片网| 日韩成人在线观看一区二区三区| 国产深夜福利视频在线观看| 男人舔女人的私密视频| 丁香六月天网| 91麻豆精品激情在线观看国产 | 在线观看免费日韩欧美大片| 18禁黄网站禁片午夜丰满| 91大片在线观看| 欧美性长视频在线观看| 欧美精品亚洲一区二区| 激情在线观看视频在线高清 | 下体分泌物呈黄色| 久久久精品区二区三区| 黄片小视频在线播放| 久热爱精品视频在线9| 一级a爱视频在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 日韩免费av在线播放| 免费少妇av软件| 他把我摸到了高潮在线观看 | 99久久精品国产亚洲精品| 国产激情久久老熟女| 亚洲av片天天在线观看| 亚洲午夜理论影院| av电影中文网址| 成人特级黄色片久久久久久久 | 19禁男女啪啪无遮挡网站| 啪啪无遮挡十八禁网站| 一区二区av电影网| 免费久久久久久久精品成人欧美视频| 亚洲精品av麻豆狂野| 一边摸一边做爽爽视频免费| 波多野结衣av一区二区av| 人人妻,人人澡人人爽秒播| 美女午夜性视频免费| 热99久久久久精品小说推荐| 日韩欧美免费精品| 欧美日韩成人在线一区二区| av又黄又爽大尺度在线免费看| 中文字幕最新亚洲高清| 精品国产国语对白av| 老司机靠b影院| 另类亚洲欧美激情| 高清视频免费观看一区二区| 亚洲人成77777在线视频| 一区二区三区精品91| 日本vs欧美在线观看视频| 中文字幕最新亚洲高清| 成人18禁高潮啪啪吃奶动态图| 我要看黄色一级片免费的| 欧美国产精品va在线观看不卡| 动漫黄色视频在线观看| 日韩视频在线欧美| 欧美日本中文国产一区发布| 视频在线观看一区二区三区| 99热网站在线观看| 黄色怎么调成土黄色| 国产精品98久久久久久宅男小说| 丰满饥渴人妻一区二区三| 久久精品aⅴ一区二区三区四区| 国产成人精品在线电影| 精品国产乱码久久久久久小说| 久久久久国产一级毛片高清牌| 久久国产精品人妻蜜桃| 欧美日韩av久久| 国产日韩欧美亚洲二区| 成年动漫av网址| 中文字幕另类日韩欧美亚洲嫩草| 肉色欧美久久久久久久蜜桃| 麻豆国产av国片精品| 久久久久久久久免费视频了| 亚洲精品中文字幕在线视频| 精品人妻在线不人妻| 人成视频在线观看免费观看| 欧美日韩亚洲国产一区二区在线观看 | 黄色成人免费大全| 男人舔女人的私密视频| 久久毛片免费看一区二区三区| 国产男女超爽视频在线观看| 高清毛片免费观看视频网站 | 精品久久久精品久久久| 九色亚洲精品在线播放| 男女无遮挡免费网站观看| 亚洲国产看品久久| 丝袜美腿诱惑在线| 性少妇av在线| 淫妇啪啪啪对白视频| 成年版毛片免费区| 乱人伦中国视频| 香蕉丝袜av| 在线观看免费午夜福利视频| 黄色丝袜av网址大全| 欧美精品亚洲一区二区| 国产男靠女视频免费网站| 女性生殖器流出的白浆| 99精国产麻豆久久婷婷| av又黄又爽大尺度在线免费看| 最新在线观看一区二区三区| 精品国内亚洲2022精品成人 | av网站在线播放免费| 免费一级毛片在线播放高清视频 | 成人特级黄色片久久久久久久 | 一个人免费在线观看的高清视频| 色精品久久人妻99蜜桃| 亚洲成a人片在线一区二区| 国精品久久久久久国模美| 欧美国产精品va在线观看不卡| 老熟妇仑乱视频hdxx| 蜜桃在线观看..| 黄色视频,在线免费观看| 亚洲熟女精品中文字幕| 亚洲avbb在线观看| 一本大道久久a久久精品| 成人国产一区最新在线观看| 日本精品一区二区三区蜜桃| 麻豆av在线久日| 亚洲国产精品一区二区三区在线| 亚洲男人天堂网一区| 在线看a的网站| 最黄视频免费看| 黄色a级毛片大全视频| 在线观看66精品国产| 久久中文看片网| 欧美亚洲 丝袜 人妻 在线| 成人三级做爰电影| 不卡av一区二区三区| 一个人免费看片子| 午夜福利,免费看| 高清黄色对白视频在线免费看| 色综合欧美亚洲国产小说| 国产aⅴ精品一区二区三区波| 日日摸夜夜添夜夜添小说| 好男人电影高清在线观看| 90打野战视频偷拍视频| 国产国语露脸激情在线看| 国产精品电影一区二区三区 | 高清在线国产一区| 水蜜桃什么品种好| 好男人电影高清在线观看| 午夜久久久在线观看| 一本大道久久a久久精品| 国产精品麻豆人妻色哟哟久久| 一夜夜www| 欧美成人免费av一区二区三区 | 国产成人免费观看mmmm| 久久精品亚洲精品国产色婷小说| 国产一区有黄有色的免费视频| 亚洲精品中文字幕一二三四区 | 亚洲三区欧美一区| 热99久久久久精品小说推荐| 一边摸一边抽搐一进一小说 | 桃花免费在线播放| 亚洲一区中文字幕在线| 18禁国产床啪视频网站| 黄色片一级片一级黄色片| 欧美人与性动交α欧美软件| 午夜福利一区二区在线看| 嫁个100分男人电影在线观看| 免费av中文字幕在线| 亚洲色图av天堂| 1024香蕉在线观看| 免费在线观看黄色视频的| 日韩免费高清中文字幕av| 色在线成人网| 母亲3免费完整高清在线观看| 中文字幕最新亚洲高清| 日韩制服丝袜自拍偷拍| 国产成人精品无人区| 久久影院123| 美女视频免费永久观看网站| 亚洲精品美女久久久久99蜜臀| av福利片在线| 我要看黄色一级片免费的| 国产日韩欧美亚洲二区| 久久人妻福利社区极品人妻图片| 成人国语在线视频| 老熟妇乱子伦视频在线观看| 国产精品二区激情视频| 99香蕉大伊视频| 美国免费a级毛片| a级片在线免费高清观看视频| 久久天躁狠狠躁夜夜2o2o| 人人妻人人添人人爽欧美一区卜| 男人舔女人的私密视频| 最近最新中文字幕大全免费视频| 黑人猛操日本美女一级片| 麻豆乱淫一区二区| 老司机影院毛片| 久久久久久亚洲精品国产蜜桃av| 国产精品久久久久成人av| avwww免费| av福利片在线| 丰满少妇做爰视频| 日本欧美视频一区| 亚洲人成电影观看| 搡老熟女国产l中国老女人| 久久精品aⅴ一区二区三区四区| 久久午夜亚洲精品久久| 亚洲av片天天在线观看| 日韩精品免费视频一区二区三区| 桃红色精品国产亚洲av| 午夜福利影视在线免费观看| 国产xxxxx性猛交| 一区二区三区精品91| 午夜福利视频精品| 日韩成人在线观看一区二区三区| 纵有疾风起免费观看全集完整版| 制服人妻中文乱码| 国产深夜福利视频在线观看| 狠狠婷婷综合久久久久久88av| www.精华液| 曰老女人黄片| 人人妻人人爽人人添夜夜欢视频| 大香蕉久久成人网| 精品国产国语对白av| 国产精品久久久av美女十八| 丝袜人妻中文字幕| 午夜91福利影院| 久久天堂一区二区三区四区| bbb黄色大片| 久热爱精品视频在线9| 90打野战视频偷拍视频| 国产成人啪精品午夜网站| 伊人久久大香线蕉亚洲五| 交换朋友夫妻互换小说| 国产亚洲欧美精品永久| 亚洲欧美精品综合一区二区三区| 久久久国产精品麻豆| 日本wwww免费看| 一本综合久久免费| 久久久欧美国产精品| 视频区图区小说| 免费日韩欧美在线观看| 久久久水蜜桃国产精品网| 亚洲国产欧美在线一区| 久久久久久久国产电影| 精品国产乱码久久久久久男人| 在线播放国产精品三级| 自线自在国产av| 免费人妻精品一区二区三区视频| 天天躁日日躁夜夜躁夜夜| 国产91精品成人一区二区三区 | av免费在线观看网站| 成人免费观看视频高清| 亚洲国产中文字幕在线视频| 另类精品久久| 日本a在线网址| 热re99久久国产66热| 热99久久久久精品小说推荐| 久久这里只有精品19| 欧美日韩精品网址| 性高湖久久久久久久久免费观看| 色综合欧美亚洲国产小说| 日韩制服丝袜自拍偷拍| 国产黄频视频在线观看| netflix在线观看网站| 久久av网站| 首页视频小说图片口味搜索| 亚洲精品成人av观看孕妇| 国产成人精品在线电影| 久久久久精品国产欧美久久久| 老熟妇仑乱视频hdxx| 交换朋友夫妻互换小说| 男男h啪啪无遮挡| 一区福利在线观看| 国产成人免费观看mmmm| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩精品亚洲av| 在线av久久热| 日韩免费av在线播放| 汤姆久久久久久久影院中文字幕| 精品国产国语对白av| 美女高潮喷水抽搐中文字幕| 狂野欧美激情性xxxx| 伊人久久大香线蕉亚洲五| 亚洲av国产av综合av卡| 成人永久免费在线观看视频 | 亚洲国产av新网站| 国产老妇伦熟女老妇高清| 精品国产一区二区三区久久久樱花| 久久久久精品国产欧美久久久| av线在线观看网站| 中亚洲国语对白在线视频| 欧美激情久久久久久爽电影 | 成年人免费黄色播放视频| 国产精品一区二区在线不卡| 亚洲美女黄片视频| 亚洲天堂av无毛| 黑人操中国人逼视频| 天堂动漫精品| 久久久久精品国产欧美久久久| 国产精品影院久久| av福利片在线| 在线观看www视频免费| 午夜福利在线观看吧| 老司机亚洲免费影院| 免费高清在线观看日韩| 黄色怎么调成土黄色| 国产一区有黄有色的免费视频| 亚洲人成电影观看| 手机成人av网站| 亚洲精品中文字幕一二三四区 | 91成年电影在线观看| 亚洲欧美一区二区三区久久| 黄色视频不卡| 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 欧美精品av麻豆av| 国产真人三级小视频在线观看| av福利片在线| 久久久久久免费高清国产稀缺| 99国产精品一区二区蜜桃av | 午夜免费鲁丝| 成年动漫av网址| 成年女人毛片免费观看观看9 | 丝袜美足系列| 99在线人妻在线中文字幕 | 久久久精品94久久精品| avwww免费| 久久久国产一区二区| 嫩草影视91久久| 久久久国产成人免费| 女同久久另类99精品国产91| 狠狠婷婷综合久久久久久88av| 精品少妇久久久久久888优播| 久热这里只有精品99| 无人区码免费观看不卡 | 久久人妻福利社区极品人妻图片| 精品国产乱子伦一区二区三区| 国产主播在线观看一区二区| 在线av久久热| 韩国精品一区二区三区| 久9热在线精品视频| 国产视频一区二区在线看| 一级,二级,三级黄色视频| 手机成人av网站| 久久99一区二区三区| netflix在线观看网站| 亚洲av电影在线进入| 丁香六月欧美| 日本黄色日本黄色录像| av又黄又爽大尺度在线免费看| 亚洲专区中文字幕在线| 国产欧美日韩一区二区精品| 亚洲熟女毛片儿| 黑人巨大精品欧美一区二区蜜桃| 老鸭窝网址在线观看| 大片免费播放器 马上看| 日韩欧美三级三区| svipshipincom国产片| a级片在线免费高清观看视频| 首页视频小说图片口味搜索| 97在线人人人人妻| 欧美人与性动交α欧美精品济南到| 少妇 在线观看| 久久99一区二区三区| 免费日韩欧美在线观看| 国产伦人伦偷精品视频| 王馨瑶露胸无遮挡在线观看| 国产一区二区三区综合在线观看| 伦理电影免费视频| av在线播放免费不卡| 亚洲天堂av无毛| 日韩免费高清中文字幕av| 人人妻人人澡人人爽人人夜夜| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇一区二区三区视频日本电影| 男女免费视频国产| 黑人操中国人逼视频| 俄罗斯特黄特色一大片| 在线观看舔阴道视频| 国产男女超爽视频在线观看| 日日爽夜夜爽网站| 国产欧美亚洲国产| 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 国产亚洲午夜精品一区二区久久| 99国产极品粉嫩在线观看| 日韩一卡2卡3卡4卡2021年| 男女午夜视频在线观看| 亚洲第一欧美日韩一区二区三区 | 丁香欧美五月| 777米奇影视久久| 国产精品久久久av美女十八| 一本色道久久久久久精品综合| 自线自在国产av| 久久这里只有精品19| 午夜福利视频在线观看免费| 成人三级做爰电影| 久久精品aⅴ一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 欧美av亚洲av综合av国产av| 大片电影免费在线观看免费| 涩涩av久久男人的天堂| 日本一区二区免费在线视频| 亚洲 国产 在线| 色婷婷av一区二区三区视频| 好男人电影高清在线观看| 热99re8久久精品国产| 日本五十路高清| 国产单亲对白刺激| 亚洲国产看品久久| 久久久精品94久久精品| 久久av网站| 黄片播放在线免费| 亚洲精品美女久久av网站| 久久 成人 亚洲| 五月开心婷婷网| 99国产精品99久久久久| 免费女性裸体啪啪无遮挡网站| 超碰成人久久| 亚洲黑人精品在线| 女同久久另类99精品国产91|