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

行星探測(cè)器GNC系統(tǒng)自主安全模式設(shè)計(jì)

2010-12-11 01:31:54黃江川王曉磊

楊 巍,黃江川,王曉磊

(1.北京控制工程研究所,北京 100190; 2.空間智能控制技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,北京 100190)

行星探測(cè)器GNC系統(tǒng)自主安全模式設(shè)計(jì)

楊 巍1,2,黃江川1,王曉磊1

(1.北京控制工程研究所,北京 100190; 2.空間智能控制技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,北京 100190)

針對(duì)行星探測(cè)器GNC系統(tǒng)自主生存能力要求高的特點(diǎn),設(shè)計(jì)具有黃極捕獲功能的慣性安全模式和無(wú)慣性系定向支持的應(yīng)急安全模式.描述安全模式下姿態(tài)確定與控制任務(wù),提出相應(yīng)的方案.通過(guò)對(duì)安全模式的任務(wù)分析可知,該模式下探測(cè)器可以自主進(jìn)行姿態(tài)機(jī)動(dòng),保證能源供給和星地通信鏈路建立.

安全模式; 自主生存; 黃極捕獲; 非主軸姿態(tài)機(jī)動(dòng)

安全模式作為探測(cè)器或衛(wèi)星的最后保障十分重要,國(guó)內(nèi)外大部分飛行器都設(shè)計(jì)了安全模式.文獻(xiàn)[1]中對(duì)安全模式進(jìn)行了描述,將安全模式作為獨(dú)立的備份控制模式,自主將衛(wèi)星定向在能源和熱控安全的姿態(tài)上.文獻(xiàn)[2]中將衛(wèi)星的初始捕獲與安全模式作為一個(gè)模式,其下分為3個(gè)子模式:速率阻尼、地球捕獲、偏航捕獲以及穩(wěn)定狀態(tài).文獻(xiàn)[3]中的安全模式包括:太陽(yáng)定向和長(zhǎng)期保持兩個(gè)子模式.由于行星探測(cè)器距離地球遠(yuǎn)、空間衰減大、通信延遲時(shí)間長(zhǎng),同時(shí)還面臨著十分復(fù)雜的空間環(huán)境,因此其安全模式不但要保證探測(cè)器的能源安全,還要在星地通信中斷時(shí)重新建立星地通信鏈路,等待地面干預(yù)處理.結(jié)合參考文獻(xiàn)中的安全模式設(shè)計(jì)方法,針對(duì)行星探測(cè)器的特點(diǎn)提出帶有黃極捕獲和地球捕獲功能的慣性安全模式以及無(wú)慣性系定向支持的應(yīng)急安全模式,用以保證行星探測(cè)器能源供給和通信的基本需要.

1 安全模式任務(wù)要求

行星探測(cè)任務(wù)由于距離地球遠(yuǎn)、信道衰減大,測(cè)控天線無(wú)法設(shè)計(jì)為全天球覆蓋,其指向受到限制,因此安全模式下不僅要保證能源供給,還要自主完成星地通信鏈路的建立.考慮到以上兩個(gè)約束條件以及整星故障的不確定性,行星探測(cè)器安全模式應(yīng)滿足以下幾個(gè)條件:

1)帆板對(duì)日定向,保證整星能源安全;

2)天線對(duì)地,保證測(cè)控鏈路安全;

3)整個(gè)過(guò)程不依賴(lài)星歷計(jì)算;

4)安全模式可長(zhǎng)期堅(jiān)持,從安全模式建立到退出,整個(gè)過(guò)程的燃料消耗要盡可能少.

從以上的任務(wù)要求可以看出,行星探測(cè)器GNC系統(tǒng)的安全模式任務(wù)復(fù)雜,需要進(jìn)行一系列的姿態(tài)機(jī)動(dòng)任務(wù)才能最終建立安全模式.此外,由于行星探測(cè)器攜帶的燃料有限,對(duì)燃料消耗也有嚴(yán)格要求.當(dāng)出現(xiàn)星地通信中斷時(shí),為了恢復(fù)星地通信,需要設(shè)計(jì)慣性安全模式和應(yīng)急安全模式以建立星地通信鏈路.

2 安全模式設(shè)計(jì)

建立探測(cè)器本體坐標(biāo)系oxbybzb,其中本體軸xbybzb為星體慣量主軸.由于太陽(yáng)系內(nèi)的行星與地球繞日運(yùn)行的軌道面基本在一個(gè)平面內(nèi),即黃道面.黃極為黃道面的法線,在慣性坐標(biāo)系中是固定的,如圖1所示建立地心赤道慣性系和地心黃道慣性系,其中OXYZ為地心赤道慣性坐標(biāo)系,OXeYeZe為地心黃道慣性坐標(biāo)系.北黃極在天龍座(赤經(jīng)180°,赤緯66.5°),南黃極則在玉夫座(赤經(jīng)0°,赤緯-66.5°).

圖1 黃極指向示意圖

2.1慣性安全模式任務(wù)描述

(1)太陽(yáng)捕獲

為了保證探測(cè)器的能源供應(yīng),在安全模式下要進(jìn)行太陽(yáng)捕獲任務(wù),使太陽(yáng)帆板指向太陽(yáng).帆板驅(qū)動(dòng)機(jī)構(gòu)旋轉(zhuǎn)軸與±yb平行,機(jī)械零位指向星體±xb軸.天線驅(qū)動(dòng)機(jī)構(gòu)分A、B兩軸,A軸以xb為軸線轉(zhuǎn)動(dòng),當(dāng)A軸處于零位時(shí),B軸以星體yb為軸線轉(zhuǎn)動(dòng),即天線零位指向在本體oxbzb平面內(nèi).具體過(guò)程是首先控制天線和帆板歸零并鎖定;進(jìn)行速率阻尼;然后利用太陽(yáng)敏感器測(cè)量太陽(yáng)敏感器方位,根據(jù)陀螺和太陽(yáng)敏感器數(shù)據(jù)進(jìn)行姿態(tài)確定;控制探測(cè)器本體+xb軸指向太陽(yáng).

(2)黃極捕獲

從行星探測(cè)器的安全模式任務(wù)要求中可知,該模式要自主完成通信鏈路建立.為了保證實(shí)現(xiàn)天線指向地球設(shè)計(jì)了黃極捕獲功能,黃極指向如圖1所示.

根據(jù)星敏感器的導(dǎo)航星表可以測(cè)量確定黃極在赤道慣性坐標(biāo)系的方位.當(dāng)完成太陽(yáng)捕獲后,控制探測(cè)器繞本體xb軸轉(zhuǎn)動(dòng),使yb軸與黃極指向平行完成黃極捕獲.

(3)地球捕獲

黃極捕獲完成后,控制探測(cè)器繞本體yb軸慢旋搜索上行信號(hào)峰值,同時(shí)帆板解除鎖定并反向跟蹤太陽(yáng).當(dāng)搜索到地面測(cè)控站峰值信號(hào)后停旋,此時(shí)天線指向地球,實(shí)現(xiàn)地球捕獲建立星地通信鏈路.

慣性安全模式的進(jìn)入條件是當(dāng)星地通信中斷超過(guò)一定時(shí)間時(shí),視為探測(cè)器嚴(yán)重故障狀態(tài),自主轉(zhuǎn)入慣性安全模式;模式的退出條件是當(dāng)?shù)厍虿东@成功時(shí),由地面遙控退出;當(dāng)連續(xù)多次地球捕獲不成功時(shí),則自主轉(zhuǎn)入應(yīng)急安全模式.

2.2應(yīng)急安全模式任務(wù)描述

應(yīng)急安全模式是在無(wú)慣性系定向支持,僅依靠陀螺和太陽(yáng)敏感器進(jìn)行姿態(tài)測(cè)量定姿的工作模式.具體任務(wù)如下:

1)應(yīng)急安全模式下,自主切換至中增益天線.中增益天線具有±16°的天線波束寬度,有利于捕獲地球;

2)重新進(jìn)行天線、帆板歸零并鎖定,太陽(yáng)捕獲;

3)太陽(yáng)捕獲后,根據(jù)中增益天線安裝位置,控制探測(cè)器俯仰軸轉(zhuǎn)動(dòng)一定角度,使中增益天線電軸方向與太陽(yáng)矢量方向夾角為20°;

4)控制本體繞太陽(yáng)矢量方向旋轉(zhuǎn),若地球、探測(cè)器和太陽(yáng)角小于20°則轉(zhuǎn)動(dòng)一圈可以捕獲地球;若24h內(nèi)未實(shí)現(xiàn)捕獲,則繼續(xù)繞探測(cè)器的俯仰軸旋轉(zhuǎn),使中增益天線電軸方向與太陽(yáng)矢量方向夾角增至40°,重復(fù)上述轉(zhuǎn)動(dòng).兩種地球捕獲方式如圖2所示,α為地球、探測(cè)器、太陽(yáng)矢量夾角,β為中增益天線波束寬度,γ為中增益天線電軸太陽(yáng)矢量偏置角(假設(shè)電軸方向與探測(cè)器本體軸方向一致).

圖2 應(yīng)急安全模式地球捕獲示意圖

以火星探測(cè)器為例,根據(jù)火星和地球繞日運(yùn)行的軌道關(guān)系可知,火星探測(cè)器的α角不超過(guò)50°,此時(shí)兩種搜索必有一種能捕獲地球.由于該模式下探測(cè)器轉(zhuǎn)軸不是星體主軸,且需每天調(diào)整一次轉(zhuǎn)軸,燃料消耗較大,因此該模式只能作為最后的安全措施.

應(yīng)急安全模式的進(jìn)入條件是慣性安全模式多次未能捕獲地球,自主轉(zhuǎn)入應(yīng)急安全模式;退出條件是當(dāng)完成地球捕獲后,等待地面干預(yù),由地面遙控退出應(yīng)急安全模式.

2.3姿態(tài)測(cè)量方案

(1)太陽(yáng)捕獲

為了實(shí)現(xiàn)太陽(yáng)捕獲采用0-1太陽(yáng)敏感器、數(shù)字太陽(yáng)敏感器搜索太陽(yáng);由數(shù)字太陽(yáng)敏感器和陀螺精確定姿.

(2)陀螺-星敏感器濾波定姿

在黃極捕獲階段,由于黃極的指向相對(duì)慣性空間是不變的.因此可以利用星敏感器測(cè)量黃極方位,從而得到探測(cè)器本體Y軸相對(duì)黃極方位的姿態(tài)偏差,進(jìn)行慣性空間姿態(tài)機(jī)動(dòng)使Y軸指向黃極方向,完成黃極捕獲.

在太陽(yáng)捕獲后、黃極慣性姿態(tài)機(jī)動(dòng)前、后,都需要利用陀螺和星敏感器進(jìn)行精確定姿.采用陀螺數(shù)據(jù)進(jìn)行姿態(tài)角預(yù)估,同時(shí)利用星敏感器修正陀螺漂移[4].

(3)陀螺預(yù)估計(jì)

在故障情況下需要快速的建立慣性定向姿態(tài),因此在黃極捕獲過(guò)程中繞X軸旋轉(zhuǎn)機(jī)動(dòng)采用推力器進(jìn)行姿態(tài)機(jī)動(dòng)使yb軸指向黃極.由于姿態(tài)機(jī)動(dòng)的角度較大、角速率較高,星敏感器無(wú)法連續(xù)地進(jìn)行觀測(cè),因此只能采用高精度的速率積分陀螺測(cè)量估計(jì)探測(cè)器的慣性姿態(tài).利用陀螺對(duì)姿態(tài)角和姿態(tài)角速度在短時(shí)間內(nèi)進(jìn)行預(yù)估計(jì),為姿態(tài)機(jī)動(dòng)控制提供姿態(tài)測(cè)量數(shù)據(jù).

2.4姿態(tài)控制方案

(1)慣性安全模式姿控方案

從節(jié)省燃料的角度,在對(duì)機(jī)動(dòng)角度和角速度要求不高的情況下,應(yīng)采用動(dòng)量輪控制;當(dāng)需要大角度姿態(tài)機(jī)動(dòng)時(shí)需采用相平面噴氣姿態(tài)控制.該模式下的太陽(yáng)捕獲過(guò)程中,可以采用動(dòng)量輪和噴氣聯(lián)合控制;黃極捕獲過(guò)程,由于姿態(tài)機(jī)動(dòng)角度和機(jī)動(dòng)角速度比較大,因此采用噴氣控制,同時(shí)動(dòng)量輪保持額定轉(zhuǎn)速.

(2)應(yīng)急安全模式姿態(tài)機(jī)動(dòng)控制方案

由應(yīng)急安全模式的任務(wù)分析可知,在地球捕獲時(shí)探測(cè)器姿態(tài)機(jī)動(dòng)轉(zhuǎn)軸不是本體主軸.在機(jī)動(dòng)過(guò)程由于主軸慣量不相等而出現(xiàn)角速度耦合的情況,給姿態(tài)機(jī)動(dòng)控制帶來(lái)了困難.目前星上常用的姿態(tài)機(jī)動(dòng)控制算法是相平面噴氣控制,但計(jì)算量相對(duì)較大.本文提出的另一種姿態(tài)機(jī)動(dòng)方案,在保證足夠精度要求的基礎(chǔ)上,具有計(jì)算量小的特點(diǎn).

1)動(dòng)力學(xué)模型

將探測(cè)器作為剛體建立其動(dòng)力學(xué)模型:

將式(2)展開(kāi)可得:

假設(shè)I1>I2>I3,且當(dāng)無(wú)外力矩作用時(shí),剛體能量和角動(dòng)量守恒.

式(4)為能量橢球和角動(dòng)量橢球,式中h為角動(dòng)量大小,ε為能量大小.非對(duì)稱(chēng)自由剛體的角速度矢量為角動(dòng)量橢球和能量橢球的交線如圖3所示.由圖3可以看出無(wú)外力矩作用時(shí),探測(cè)器無(wú)法繞固定方向即太陽(yáng)矢量方向轉(zhuǎn)動(dòng).因此,需要采用非慣性主軸的姿態(tài)機(jī)動(dòng)算法,實(shí)現(xiàn)應(yīng)急安全模式下的姿態(tài)機(jī)動(dòng)要求,同時(shí)滿足燃料消耗的要求.

圖3 非對(duì)稱(chēng)剛體自由轉(zhuǎn)動(dòng)的角速度矢量

2)燃料消耗分析

由圖3可以看出要維持機(jī)動(dòng)方案中的角速度則需要相應(yīng)的力矩維持標(biāo)稱(chēng)轉(zhuǎn)動(dòng).由于應(yīng)急安全模式的運(yùn)行時(shí)間長(zhǎng),需要對(duì)燃料消耗進(jìn)行分析.

假設(shè)三軸慣量為:

脈沖噴氣力矩均為T(mén)im=8N·m,

推力器推力為F=10N,

推力器比沖為Isp=2746N·s/kg,

自旋軸與xb軸夾角為θ=20°,

則三軸角速度為[0.282 0 0.1026]T(°)/s.

由式(3)得:

在一個(gè)自旋周期內(nèi)對(duì)T積分:

假設(shè)噴氣控制為脈沖力矩,令Δt為等效的脈沖噴氣作用時(shí)間則有:

令P為燃料消耗,則:

代入具體參數(shù)得探測(cè)器自旋一周20min的燃料消耗量:P=1.444×10-3kg,則12h消耗燃料量為0.052kg.

3)非主軸自旋姿態(tài)機(jī)動(dòng)算法

星體控制力矩與噴氣點(diǎn)火指令有以下關(guān)系:

將探測(cè)器相對(duì)于慣性空間的本體角速度表示為歐拉軸-角式:

將式(10)化為矩陣形式并代入式(2)得

為了便于分析,將式(11)改寫(xiě)為:

式中,|δi|為[0,1]之間的連續(xù)量,p=Ir,u=-r×Ir.

進(jìn)一步得

式中,g=A-1p,h=A-1u.

圖4所示為非主軸旋轉(zhuǎn)姿態(tài)機(jī)動(dòng)算法框圖.

圖中,E為姿態(tài)誤差,ΔθC為每個(gè)采樣周期的指令角度增量,Δθ為每個(gè)采樣周期陀螺測(cè)量的角度增量,D為姿態(tài)角誤差和角速度誤差的合并誤差,K為姿態(tài)誤差增益,S為噴氣指令向量,ωC為指令三軸角速度,ω為估計(jì)三軸角速度.

根據(jù)前面的動(dòng)力學(xué)分析有:

4)角速度估計(jì)器設(shè)計(jì)

由圖4可知陀螺的高頻測(cè)量噪聲可能會(huì)導(dǎo)致錯(cuò)誤噴氣,因此需要對(duì)陀螺測(cè)得的角度增量經(jīng)角速度估計(jì)器進(jìn)行處理,這樣得到的估計(jì)角速度與指令角速度進(jìn)行比較得到角速度誤差.

考慮ω的離散z變換:

由式(17)可以看出,該式具有低通特性,因此可得:

由于

可得

將式(19)、(21)代入式(18)可得:

將式(9)化為矩陣形式后代入式(2)得:

當(dāng)實(shí)際角速度與指令角速度相等時(shí)并結(jié)合式(11),則式(20)可化為:

式(18)、(22)和(25)組成了角速度估計(jì)器的計(jì)算公式.

5)算法仿真

a)仿真參數(shù)

采樣周期:ΔT=0.1s,

死區(qū)閾值:θD=0.8°,

噴管推力:F=10N,

脈沖噴氣力矩:Tim=8N·m,

陀螺測(cè)量噪聲方差:0.05(°)/h.

b)仿真結(jié)果

圖5 三軸估計(jì)姿態(tài)角速度

圖6 三軸姿態(tài)誤差

圖7 三軸噴氣指令

c)燃料消耗分析

噴氣控制采用10N推力器,假設(shè)單分支情況下,每軸兩個(gè)推力器分別提供該軸正向和負(fù)向的控制力矩.設(shè)推力器額定比沖為2746N·s/kg.根據(jù)噴氣次數(shù)和噴氣時(shí)間間隔,可得旋轉(zhuǎn)一周的燃料消耗:

式中:n為噴氣次數(shù);F為推力器推力;ΔT為每次噴氣時(shí)長(zhǎng);Isp為推力器額定比沖.

由仿真程序可以計(jì)算出應(yīng)急安全模式運(yùn)行12h的總噴氣次數(shù)n,再由式(26)可以得到在姿態(tài)誤差精度為1°要求下的燃料消耗為0.0579kg,基本與理論計(jì)算結(jié)果一致.

d)結(jié)果分析

通過(guò)仿真曲線可知,在姿態(tài)機(jī)動(dòng)過(guò)程中三軸姿態(tài)誤差小于1°,姿態(tài)角速度誤差小于0.01(°)/s.由于安全模式下姿態(tài)精度要求較低,因此能夠滿足姿態(tài)精度的要求.同時(shí)通過(guò)仿真還可計(jì)算出姿態(tài)機(jī)動(dòng)的燃料消耗量,與理論計(jì)算值基本一致,因此該算法可以實(shí)現(xiàn)繞非本體主軸的姿態(tài)機(jī)動(dòng)控制,滿足在應(yīng)急安全模式下進(jìn)行地球捕獲的姿態(tài)精度要求和燃料消耗要求.由算法分析過(guò)程中可知,該算法采用增量方法更新姿態(tài)誤差,不計(jì)算姿態(tài)誤差矩陣,從而使計(jì)算量減少,適合于應(yīng)急安全模式下的非本體軸快速機(jī)動(dòng)姿態(tài)控制.此外,當(dāng)選擇比較小的指令角速度時(shí),轉(zhuǎn)動(dòng)一周所消耗的燃料也相應(yīng)減少,因此可以根據(jù)對(duì)燃料消耗的限制靈活選擇自旋角速度.

3 結(jié) 論

由于一般的近地空間飛行器測(cè)控天線可以設(shè)計(jì)為全天球覆蓋,因此其GNC系統(tǒng)安全模式僅需滿足能源供給即可,而行星探測(cè)器的安全模式不僅要保證能源供給,還要實(shí)現(xiàn)故障情況下自主建立星地通信鏈路.本文針對(duì)行星探測(cè)器自主生存能力要求高的特點(diǎn),從探測(cè)器任務(wù)需求的角度出發(fā),設(shè)計(jì)了行星探測(cè)器GNC系統(tǒng)的自主安全模式,包括慣性安全模式和應(yīng)急安全模式.在慣性安全模式中,利用了太陽(yáng)系內(nèi)行星基本均運(yùn)行在黃道面內(nèi)的特點(diǎn),提出了黃極捕獲方式使探測(cè)器測(cè)控天線指向與地球處于同一平面,進(jìn)而實(shí)現(xiàn)地球捕獲;在應(yīng)急安全模式下,針對(duì)該模式姿態(tài)機(jī)動(dòng)轉(zhuǎn)軸不是本體軸,提出了非本體軸姿態(tài)機(jī)動(dòng)方案,可以滿足應(yīng)急模式下姿態(tài)機(jī)動(dòng)任務(wù)要求,同時(shí)又具有計(jì)算量小的特點(diǎn).

[1] David M,Joseph T.Guidance,navigation,and control subsystem for the EOS-AM spacecraft [C].AIAA Guidance,Navigation and Control Conference,Hilton Head Island,USA,August,1992

[2] Winkler S,Duske N,Wiedermann G,Gockel W.Earth-oriented safe mode: concept,design,and results for the GMES sentinel-2 satellite [C].AIAA Modeling and Simulation Technologies Conference,EADS Astrium,Friedrichshafen,Germany,2009

[3] Robertson B,Sabelhaus P,Mendenhall T,Fesq L.The recovery of TOMS-EP [R].AAS 00-76

[4] 李琳琳.衛(wèi)星自主軌道確定及姿態(tài)確定技術(shù)研究[D].中國(guó)科學(xué)院空間科學(xué)與應(yīng)用中心,2003

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

[6] 耿長(zhǎng)福.航天器動(dòng)力學(xué)[M].北京: 中國(guó)科學(xué)技術(shù)出版社,2001

[7] 屠善澄.衛(wèi)星姿態(tài)動(dòng)力學(xué)與控制(2)[M].北京: 宇航出版社,2002

GNCSystemAutonomousSafeModeDesignforInterplanetarySpacecraft

YANG Wei1,2,HUANG Jiangchuan1,WANG Xiaolei1

(1.BeijingInstituteofControlEngineering,Beijing100190,China;2.NationalKeyLaboratoryofScienceandTechnologyonSpaceIntelligentControl,Beijing100190,China)

Two safe modes including inertial reference frame safe mode with ecliptic pole pointing and emergency safe mode without inertial reference frame orientation are designed for the requirement of autonomous survival capability of interplanetary spacecraft.Tasks of attitude determination and control in the safe mode are described in this paper,and then the corresponding schemes are proposed.It is shown from the task analysis for safe modes that the spacecraft autonomous attitude maneuver will guarantee power supply to be provided and communication link established.

safe mode; autonomous survival; ecliptic pole pointing; non-principal axis attitude maneuver

V448.2

A

1674-1579(2010)04-0006-06

2010-04-20

楊巍(1983—),男,黑龍江人,碩士研究生,研究方向?yàn)楹教炱鲗?dǎo)航、制導(dǎo)與控制 (e-mail: wangyi1983128@163.com).

主站蜘蛛池模板: 国产午夜福利亚洲第一| 国产H片无码不卡在线视频| 日本不卡视频在线| 国产在线观看一区精品| 一级一级特黄女人精品毛片| 国产不卡在线看| 国产一级二级三级毛片| 福利一区三区| 一本大道无码日韩精品影视| 成年A级毛片| 性欧美久久| 亚洲精品中文字幕午夜 | 久久精品国产电影| 好紧好深好大乳无码中文字幕| 精品国产污污免费网站| 四虎国产成人免费观看| 国产精品自拍合集| 欧美在线伊人| 一级看片免费视频| 欧美一区二区三区不卡免费| 日本免费新一区视频| 国产原创第一页在线观看| 亚洲动漫h| AV熟女乱| 亚洲av无码牛牛影视在线二区| 日韩天堂网| 亚洲欧美另类久久久精品播放的| 亚洲欧美日韩成人高清在线一区| 91视频日本| 超清无码熟妇人妻AV在线绿巨人| 亚洲无码精品在线播放| 婷婷中文在线| 91精品小视频| 色偷偷综合网| 热这里只有精品国产热门精品| 国产99视频免费精品是看6| 亚洲有无码中文网| 国产a v无码专区亚洲av| 天天色天天综合| 国产欧美视频综合二区| 国产欧美视频在线观看| 成人亚洲视频| 71pao成人国产永久免费视频| 国产永久在线视频| 一级毛片免费观看不卡视频| 国产在线观看第二页| 日韩欧美中文| 色婷婷国产精品视频| 在线欧美a| 91精品国产自产在线观看| 欧美伦理一区| 国产成人啪视频一区二区三区| 亚洲综合二区| 婷婷激情亚洲| 精品欧美视频| 美女一区二区在线观看| 亚洲国产成熟视频在线多多| 亚洲欧洲日韩国产综合在线二区| h视频在线播放| 欧美激情视频二区| 亚洲视频影院| 丰满人妻一区二区三区视频| 日韩精品无码免费专网站| 青青草综合网| 国内精品视频| 国产精品欧美亚洲韩国日本不卡| 国产精品久久久久久久伊一| 91免费片| 国产91在线免费视频| 暴力调教一区二区三区| 亚洲欧美激情小说另类| 中文字幕亚洲电影| 怡红院美国分院一区二区| 精品无码国产自产野外拍在线| 国产美女在线观看| 久久亚洲黄色视频| 91视频首页| 国产成人精品18| 黄色福利在线| 无码精品一区二区久久久| 国产黑丝一区| 好吊妞欧美视频免费|