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

一種新的航天器姿態(tài)快速測量方法

2017-06-27 08:14:18范生宏魯利剛
無線電通信技術(shù) 2017年4期
關(guān)鍵詞:測量產(chǎn)品

張 杰,張 巍,范生宏,王 順,魯利剛

(1. 北京衛(wèi)星制造廠,北京 100080;2.北京控制工程研究所,北京 100190;3. 北京普達(dá)迪泰科技有限公司,北京 100083)

一種新的航天器姿態(tài)快速測量方法

張 杰1,張 巍2,范生宏3,王 順3,魯利剛1

(1. 北京衛(wèi)星制造廠,北京 100080;2.北京控制工程研究所,北京 100190;3. 北京普達(dá)迪泰科技有限公司,北京 100083)

針對航天器組件裝調(diào)及實(shí)驗(yàn)階段,姿態(tài)準(zhǔn)直測量過程需要頻繁建站的問題,提出了一種航天器組件姿態(tài)快速測量方法。該方法采用調(diào)整裝置,將衛(wèi)星本體坐標(biāo)系XOY面調(diào)整到與大地水平,根據(jù)整星立方鏡和組件立方鏡的標(biāo)定數(shù)據(jù),快速反算出用于組件準(zhǔn)直測量的經(jīng)緯儀的理論指向,指導(dǎo)產(chǎn)品組件姿態(tài)的快速調(diào)整。經(jīng)實(shí)驗(yàn)驗(yàn)證,該快速調(diào)整方法精度滿足設(shè)計(jì)要求,測量效率大幅提升。

航天器姿態(tài);準(zhǔn)直測量;快速調(diào)整

0 引言

為了保證航天器的正常飛行和工作,地面總裝階段必須將航天器組件(星敏、天線、推進(jìn)及控制裝置等)調(diào)整到必要的幾何精度[1]。由于產(chǎn)品結(jié)構(gòu)復(fù)雜、尺寸較小或被測要素不易復(fù)現(xiàn),較多情況下采用準(zhǔn)直測量法予以解決[2]。即在該組件上安裝立方鏡(簡稱產(chǎn)品立方鏡),通過經(jīng)緯儀測量被測立方鏡與基準(zhǔn)立方鏡之間的關(guān)系,來確定被測組件的位置和姿態(tài)[3]。

立方鏡在航天器組件的位姿測量應(yīng)用比較廣泛[4],國內(nèi)眾多科研院所和高校都進(jìn)行過相關(guān)的研究和實(shí)踐,大都采用經(jīng)緯儀工業(yè)測量系統(tǒng)建站和準(zhǔn)直測量結(jié)合的方法[5-9],但該方法實(shí)際測量過程中需要多次互瞄建站,效率較低。提出了一種新的姿態(tài)快速測量方法,并驗(yàn)證了其精度和可行性。

1 傳統(tǒng)方法

傳統(tǒng)的測量方法一般使用經(jīng)緯儀工業(yè)系統(tǒng),通過多臺經(jīng)緯儀互瞄定向和尺度標(biāo)定構(gòu)建空間測量系統(tǒng)[7];對2組立方鏡準(zhǔn)直測量,分別計(jì)算立方鏡坐標(biāo)系與經(jīng)緯儀坐標(biāo)系間的關(guān)系;通過關(guān)系傳遞和數(shù)據(jù)計(jì)算,得到2組立方鏡坐標(biāo)系之間的關(guān)系,進(jìn)而得到產(chǎn)品組件姿態(tài)參數(shù)。

應(yīng)用該方法,至少需要使用4臺經(jīng)緯儀,保證每組立方鏡的2個(gè)準(zhǔn)直面可被觀測。當(dāng)組件實(shí)際姿態(tài)與理想狀態(tài)偏離較大或組件經(jīng)調(diào)整后對立方鏡再次準(zhǔn)直時(shí),可能超出經(jīng)緯儀準(zhǔn)直測量范圍,需要重新調(diào)整儀器,再次構(gòu)建空間測量系統(tǒng),然后重新測量、分析;至組件姿態(tài)偏離量在誤差允許范圍內(nèi)時(shí),終止建站操作。整個(gè)流程如圖1所示。

圖1 傳統(tǒng)方法測量流程圖

2 改進(jìn)方法

由于經(jīng)緯儀準(zhǔn)直量程有限,組件初始狀態(tài)未知,產(chǎn)品姿態(tài)調(diào)整控制過程中,反復(fù)的系統(tǒng)建站操作,繁復(fù)的數(shù)據(jù)處理環(huán)節(jié),導(dǎo)致方法實(shí)際應(yīng)用效果不太理想。此處對傳統(tǒng)方法提出改進(jìn)。

如圖2所示,衛(wèi)星裝配過程中,將衛(wèi)星本體置放于支撐裝置上,支撐裝置的位姿可通過調(diào)整裝置進(jìn)行調(diào)節(jié)。

圖2 產(chǎn)品立方鏡準(zhǔn)直示意圖

經(jīng)緯儀的儀器坐標(biāo)系(也稱為測站坐標(biāo)系),建立方式定義如下:儀器調(diào)置呈水平狀態(tài),儀器三軸中心為原點(diǎn),過原點(diǎn)沿鉛垂向上為+Z軸,零度方位角視準(zhǔn)方向?yàn)?X軸,由右手規(guī)則確定Y軸方向。儀器坐標(biāo)系的XOY平面為水平面。假設(shè)衛(wèi)星坐標(biāo)系Z軸向上,XOY面與大地近似水平。則調(diào)整過程如下:

① 根據(jù)衛(wèi)星坐標(biāo)系和衛(wèi)星立方鏡坐標(biāo)系的標(biāo)定數(shù)據(jù),計(jì)算出衛(wèi)星坐標(biāo)系XOY面與大地水平時(shí),準(zhǔn)直衛(wèi)星立方鏡的經(jīng)緯儀T3、T4的水平角和天頂距;

② 根據(jù)步驟①計(jì)算出的理論水平角和天頂距,經(jīng)緯儀對衛(wèi)星立方鏡姿態(tài)進(jìn)行實(shí)時(shí)觀測,指導(dǎo)調(diào)整裝置使整星XOY面與大地嚴(yán)格水平,調(diào)整水平后,將準(zhǔn)直整星X軸的水平角傳遞給T0,并將經(jīng)緯儀T0水平度盤零方向置為與X軸平行;

③ 根據(jù)設(shè)計(jì)參數(shù)、標(biāo)定數(shù)據(jù)等,計(jì)算得到產(chǎn)品立方鏡與衛(wèi)星坐標(biāo)系的關(guān)系,即在衛(wèi)星坐標(biāo)系下,準(zhǔn)直產(chǎn)品立方鏡的經(jīng)緯儀T1、T2的水平角和天頂距;

④ 將經(jīng)緯儀搬到準(zhǔn)直產(chǎn)品組件立方鏡的位置T1和T2,根據(jù)步驟③計(jì)算出的理論值調(diào)整產(chǎn)品姿態(tài),對產(chǎn)品立方鏡準(zhǔn)直圖像予以實(shí)時(shí)觀測、判讀,對調(diào)整操作予以即時(shí)指導(dǎo);

⑤ 產(chǎn)品姿態(tài)逼近理論值后,根據(jù)觀測儀器參數(shù)值,計(jì)算得到產(chǎn)品立方鏡和衛(wèi)星坐標(biāo)系的關(guān)系,進(jìn)而得到在衛(wèi)星坐標(biāo)系下產(chǎn)品的最終姿態(tài)。

2.1 準(zhǔn)直衛(wèi)星立方鏡的經(jīng)緯儀角度計(jì)算

設(shè)衛(wèi)星坐標(biāo)系系O-XYZ先平移(xo,yo,zo),再旋轉(zhuǎn) (εx,εy,εz),最后縮放k倍后,轉(zhuǎn)換到衛(wèi)星立方鏡坐標(biāo)O′-X′Y′Z′。點(diǎn)P在O-XYZ中的坐標(biāo)為(x,y,z),在O′-X′Y′Z′中的坐標(biāo)為(x′,y′,z′),則衛(wèi)星坐標(biāo)系與衛(wèi)星立方鏡坐標(biāo)系間的關(guān)系如下所示[10]:

(1)

經(jīng)緯儀T3準(zhǔn)直立方鏡坐標(biāo)系的X軸,T4準(zhǔn)直立方鏡坐標(biāo)系的Y軸,則衛(wèi)星坐標(biāo)系XOY面與大地水平時(shí),經(jīng)緯儀T3、T4準(zhǔn)直衛(wèi)星立方鏡工作面時(shí)的角度H3、V3和H4、V4分別為:

(2)

2.2 指導(dǎo)調(diào)整衛(wèi)星坐標(biāo)系XOY面與大地水平

如圖3所示,觀測經(jīng)緯儀目鏡視場中準(zhǔn)直圖像,通過儀器瞄準(zhǔn)十字線與準(zhǔn)直圖像水平線的位置關(guān)系,即可以準(zhǔn)確判斷出衛(wèi)星立方鏡工作面法線相對理論狀態(tài)的偏離方向,并計(jì)算出偏離量;調(diào)整經(jīng)緯儀角度調(diào)節(jié)鈕,使兩線重合,即可由儀器表盤讀數(shù)準(zhǔn)確得到立方鏡工作面的法向偏差。如儀器目鏡視場中不出現(xiàn)準(zhǔn)直圖像,說明在該方向上,衛(wèi)星坐標(biāo)系XOY平面與水平面夾角過大,超出經(jīng)緯儀準(zhǔn)直量程??梢韵葘πl(wèi)星立方鏡準(zhǔn)直,通過準(zhǔn)直狀態(tài)下儀器俯仰角與90°的差值,判斷出衛(wèi)星立方鏡工作面法向偏離方向及大小,對支撐調(diào)節(jié)裝置進(jìn)行初步調(diào)整,直到偏離量逼近允許誤差區(qū)間(-k,k)。其中可能需要對經(jīng)緯儀進(jìn)行數(shù)次位姿調(diào)整,確??蓪πl(wèi)星立方鏡實(shí)現(xiàn)準(zhǔn)直。

圖3 經(jīng)緯儀準(zhǔn)直示意圖

經(jīng)緯儀T3、T4俯仰角調(diào)整至理論值,實(shí)時(shí)監(jiān)測準(zhǔn)直圖像,根據(jù)支撐調(diào)節(jié)裝置調(diào)節(jié)點(diǎn)和衛(wèi)星立方鏡工作面法向的空間位置關(guān)系,逐步調(diào)整各調(diào)節(jié)點(diǎn)高度。經(jīng)過精確調(diào)整,可將衛(wèi)星立方鏡單個(gè)方向的法向偏差值控制在0.001°以內(nèi)。

2.3 衛(wèi)星坐標(biāo)系X軸方位傳遞

經(jīng)緯儀T3、T4空間互瞄,尺度標(biāo)定,構(gòu)建空間測站,建立測量坐標(biāo)系T-XYZ;對衛(wèi)星立方鏡準(zhǔn)直測量,得到衛(wèi)星立方鏡坐標(biāo)系C-XYZ與T-XYZ的關(guān)系;根據(jù)衛(wèi)星立方鏡C-XYZ與衛(wèi)星坐標(biāo)系S-XYZ的關(guān)系,可得到衛(wèi)星坐標(biāo)系與測量坐標(biāo)系T-XYZ的關(guān)系。由此,可計(jì)算出衛(wèi)星坐標(biāo)系X軸與測量坐標(biāo)系的空間位置關(guān)系,得到其與經(jīng)緯儀T3(或T4)儀器坐標(biāo)系X軸的空間夾角β。

經(jīng)過上步調(diào)整,衛(wèi)星坐標(biāo)系XOY平面與經(jīng)緯儀T3(或T4)儀器坐標(biāo)系XOY平面呈平行狀態(tài),將經(jīng)緯儀T3(或T4)方位角調(diào)置呈β(或360°-β),即可使儀器視準(zhǔn)方向與衛(wèi)星坐標(biāo)系X軸方向重合;將儀器該視準(zhǔn)方向設(shè)為零度方位角,傳遞于經(jīng)緯儀T0,則儀器坐標(biāo)系T0-XYZ與衛(wèi)星坐標(biāo)系S-XYZ方向完全一致。

因衛(wèi)星上需要調(diào)整的組件數(shù)量多,裝調(diào)過程中始終保持經(jīng)緯儀T0位置,將其作為基準(zhǔn)儀器,傳遞基準(zhǔn)角度。

2.4 準(zhǔn)直產(chǎn)品立方鏡的經(jīng)緯儀水平角和天頂距計(jì)算

根據(jù)產(chǎn)品坐標(biāo)系與衛(wèi)星坐標(biāo)系之間理論關(guān)系、產(chǎn)品立方鏡與產(chǎn)品坐標(biāo)系標(biāo)定關(guān)系,可計(jì)算得到產(chǎn)品立方鏡相對衛(wèi)星坐標(biāo)系的關(guān)系。計(jì)算過程如下:

產(chǎn)品坐標(biāo)系下的坐標(biāo)表示為(X,Y,Z),產(chǎn)品立方鏡坐標(biāo)系下的坐標(biāo)為(X′,Y′,Z′),衛(wèi)星坐標(biāo)系下的坐標(biāo)為(X″,Y″,Z″)。以產(chǎn)品坐標(biāo)系為當(dāng)前坐標(biāo)系,其參數(shù)表示為(0,0,0,0,0,0,1),產(chǎn)品立方鏡坐標(biāo)系的坐標(biāo)參數(shù)為(X01,Y01,Z01,εx01,εy01,εz01,k01),衛(wèi)星坐標(biāo)系的坐標(biāo)系參數(shù)表示為(X02,Y02,Z02,εx02,εy02,εz02,k02),則產(chǎn)品立方鏡坐標(biāo)系在衛(wèi)星坐標(biāo)系下的坐標(biāo)系的平移量和旋轉(zhuǎn)矩陣推導(dǎo)如下[11-12]:

(3)

令M=R2-1″,則可由式(3)推導(dǎo)得出用于準(zhǔn)直的經(jīng)緯儀水平角和天頂距值。

2.5 產(chǎn)品姿態(tài)測量與調(diào)整

如圖2對現(xiàn)場進(jìn)行布局,保持經(jīng)緯儀T0初始位置。經(jīng)緯儀T1、T2對產(chǎn)品立方鏡兩相鄰工作面進(jìn)行準(zhǔn)直,其中經(jīng)緯儀T1瞄準(zhǔn)立方鏡+X方向,經(jīng)緯儀T2瞄準(zhǔn)立方鏡+Y方向;經(jīng)緯儀T0將零度方位角傳遞于經(jīng)緯儀T1、T2。則經(jīng)緯儀T1、T2儀器坐標(biāo)系與衛(wèi)星坐標(biāo)系S-XYZ方向完全一致。

對產(chǎn)品立方鏡精確準(zhǔn)直測量,由2臺儀器得到4個(gè)參數(shù):經(jīng)緯儀T1俯仰角V1、方位角H1及經(jīng)緯儀T2俯仰角V2、方位角H2。

因而,在產(chǎn)品調(diào)整階段,可根據(jù)經(jīng)緯儀讀數(shù)對產(chǎn)品姿態(tài)予以快速判斷:① 當(dāng)其中一項(xiàng)參數(shù)偏差超過最大允許誤差k時(shí),可判定產(chǎn)品姿態(tài)超差;② 當(dāng)所有參數(shù)偏差均小于最大允許誤差k的0.5倍時(shí),可判定產(chǎn)品姿態(tài)合格。

當(dāng)產(chǎn)品嚴(yán)重偏離理想狀態(tài)時(shí),與衛(wèi)星坐標(biāo)系調(diào)置水平狀態(tài)相同,根據(jù)經(jīng)緯儀T1、T2儀器讀數(shù)指導(dǎo)初調(diào)過程,直至2個(gè)方向的3個(gè)姿態(tài)參數(shù)偏差均位于經(jīng)緯儀準(zhǔn)直量程內(nèi),其中涉及對經(jīng)緯儀T1、T2的多次位姿調(diào)整。

2.6 姿態(tài)參數(shù)計(jì)算

在判定產(chǎn)品姿態(tài)合格后,需根據(jù)儀器讀數(shù),精確計(jì)算產(chǎn)品的實(shí)際狀態(tài),即確定產(chǎn)品立方鏡坐標(biāo)系和衛(wèi)星坐標(biāo)系的關(guān)系。

經(jīng)對產(chǎn)品立方鏡嚴(yán)格準(zhǔn)直后,經(jīng)緯儀T1俯仰角V1、方位角H1,經(jīng)緯儀T2俯仰角V2、方位角H2。由于T1測站準(zhǔn)直向量的反方向?yàn)榕c衛(wèi)星坐標(biāo)系+X軸平行的向量,T2測站準(zhǔn)直向量的反方向?yàn)榕c衛(wèi)星坐標(biāo)系+Y軸平行的向量。設(shè)2個(gè)坐標(biāo)系的轉(zhuǎn)換矩陣如下:

(4)

(5)

則觀測角度與坐標(biāo)系旋轉(zhuǎn)矩陣間的關(guān)系如下:

(6)

a11、a21、a31即為產(chǎn)品立方鏡X軸在衛(wèi)星坐標(biāo)系下的向量,a12、a22、a32為產(chǎn)品立方鏡Y軸在衛(wèi)星坐標(biāo)系下的向量,X軸與Y軸叉乘可得到產(chǎn)品坐標(biāo)系Z軸在衛(wèi)星坐標(biāo)系下的向量,即旋轉(zhuǎn)矩陣的a13、a23、a33。根據(jù)準(zhǔn)直解算出產(chǎn)品立方鏡和衛(wèi)星坐標(biāo)系的實(shí)際關(guān)系,對比根據(jù)設(shè)計(jì)和標(biāo)定數(shù)據(jù)計(jì)算出的產(chǎn)品立方鏡參數(shù),即可解獲得產(chǎn)品的實(shí)際狀態(tài)。

3 實(shí)驗(yàn)驗(yàn)證

為了驗(yàn)證該方法的可行性,采用傳統(tǒng)方法先進(jìn)行準(zhǔn)直測量,再采用快速測量方法進(jìn)行3組準(zhǔn)直測試,實(shí)驗(yàn)結(jié)果如表1所示。

表1 標(biāo)定實(shí)驗(yàn)數(shù)據(jù)/(°)

測量方法RxRyRz傳統(tǒng)方法359.3040112.2878837.99078快速方法第1次359.3032512.2879437.99032快速方法第2次359.3039812.2886237.99116快速方法第3次359.3038112.2869237.99103

從表1測試結(jié)果可以看出,3組快速測量實(shí)驗(yàn)自身重復(fù)性較好,其測試結(jié)果與傳統(tǒng)方法標(biāo)定結(jié)果各軸的差均小于3 s,可以滿足大部分航天器組件的總裝精測要求。

4 結(jié)束語

本文對航天器準(zhǔn)直測量的傳統(tǒng)方法進(jìn)行了介紹,針對傳統(tǒng)方法存在的問題,提出了航天器姿態(tài)測量的快速測量方法,對方法的原理、使用范圍和實(shí)際效果進(jìn)行了詳細(xì)介紹。相對于傳統(tǒng)的測量流程,該方法實(shí)現(xiàn)了對組件姿態(tài)的快速測量和調(diào)整,簡化了經(jīng)緯儀搬站和定向過程,極大地提高準(zhǔn)直測量的效率和簡易型,應(yīng)用前景廣闊。

[1] 徐福祥.衛(wèi)星工程概論[M].北京:中國宇航出版社,2004.

[2] 張杰,王京海,石磊,等. 一種經(jīng)緯儀自準(zhǔn)直姿態(tài)實(shí)時(shí)測量方法[J]. 光電工程,2015,42(5):39-44.

[3] 楊再華. 航天器總裝精度測量方法分析[J]. 航天器環(huán)境工程,2007,24(6):390-392.

[4] 柏宏武,冀有志. 立方鏡在航天器天線總裝測量中的應(yīng)用[J]. 空間電子技術(shù),2013,10(2):62-66,82.

[5] 楊振,李廣云,黃桂平,等. 基于自準(zhǔn)直的立方鏡姿態(tài)測量方法研究[J]. 宇航計(jì)測技術(shù),2010,30(2):14-17.

[6] 黃桂平,葉聲華,李廣云,等. 經(jīng)緯儀非接觸大尺寸三坐標(biāo)測量系統(tǒng)的開發(fā)及其在航天器檢測中的應(yīng)用[J]. 上海計(jì)量測試,2002,29(3):9-12.

[7] 鄧勇,李宗春,范百興,等. MetroIn在飛船精度測量中的應(yīng)用[J]. 導(dǎo)彈與航天運(yùn)載技術(shù),2009(4):47-51.

[8] 楊振,李廣云,范百興,等. 多種互瞄條件下立方鏡準(zhǔn)直傳遞方法研究[J]. 宇航計(jì)測技術(shù),2013(6):10-12.

[9] 沈兆欣,陳曉暉. 電子經(jīng)緯儀測量系統(tǒng)中立方鏡坐標(biāo)系建立技術(shù)探討[J]. 宇航計(jì)測技術(shù),2006,26(4):73-75.

[10]黃桂平. 數(shù)字近景工業(yè)攝影測量關(guān)鍵技術(shù)研究與應(yīng)用[D].天津:天津大學(xué),2005.

[11]魏新國,張廣軍,樊巧云,等. 太陽敏感器測量坐標(biāo)系與立方鏡坐標(biāo)系轉(zhuǎn)換矩陣的標(biāo)定[J]. 光學(xué)精密工程,2009,17(4):849-853.

[12]王智,吳國棟. 測繪相機(jī)立方鏡與星敏立方鏡轉(zhuǎn)換矩陣的標(biāo)定[J]. 光學(xué)精密工程,2012,20(1):96-101.

A New Method for Rapid Measurement of Spacecraft Attitude

ZHANG Jie1,ZHANG Wei2,FAN Sheng-hong3,WANG Shun3,LU Li-gang1

(1. Beijing Spacecrafts,Beijing 100194,China;2. Beijing Institute of Control Engineering,Beijing 100190,China;3. Beijing Prodetec Technology Co.,Ltd,Beijing 100083,China)

Aiming at the problem of frequent alignment in the process of spacecraft alignment and experiment,a new method is proposed to measure the attitude of spacecraft components. This method uses the adjusting device to adjust the XOY surface of satellite coordinate system parallel with the Geodetic plane,then according to the calibration data of entire satellite cubic mirror and component cubic mirror,calculates the theodolite angle quickly,and guide the rapid adjustment of product component attitude. The experimental results show that the precision of this method can meet the design requirements,and the measurement efficiency is significantly improved.

spacecraft attitude;alignment measurement;rapid adjustment

10. 3969/j.issn. 1003-3114. 2017.04.16

張杰,張巍,范生宏,等. 一種新的航天器姿態(tài)快速測量方法[J].無線電通信技術(shù),2017,43(4):67-70.

[ZHANG Jie,ZHANG Wei,FAN Shenghong,et al. A New Method for Rapid Measurement of Spacecraft Attitude[J].Radio Communications Technology,2017,43(4):67-70. ]

2017-04-20

國防科研基礎(chǔ)計(jì)劃項(xiàng)目

張 杰( 1979—) ,男,高級工程師,主要研究方向:精密測量技術(shù)。范生宏(1978—),男,博士,主要研究方向:數(shù)字近景攝影測量與精密工程測量技術(shù)。

文獻(xiàn)標(biāo)志碼:A 文章編號:1003-3114(2017) 04-67-4

猜你喜歡
測量產(chǎn)品
好產(chǎn)品,可持續(xù)
從靈感出發(fā),邂逅好產(chǎn)品
把握四個(gè)“三” 測量變簡單
滑動摩擦力的測量和計(jì)算
滑動摩擦力的測量與計(jì)算
測量的樂趣
測量
2015產(chǎn)品LOOKBOOK直擊
Coco薇(2015年1期)2015-08-13 02:23:50
新產(chǎn)品
玩具(2009年10期)2009-11-04 02:33:14
產(chǎn)品
主站蜘蛛池模板: 91外围女在线观看| 国内99精品激情视频精品| 久久人搡人人玩人妻精品一| 91福利免费视频| 国产成人1024精品下载| 无码AV日韩一二三区| 麻豆国产精品| 日韩欧美91| 国产成人精品亚洲日本对白优播| 国产亚洲精品无码专| AV熟女乱| 欧美黄色网站在线看| 九色在线观看视频| 国产永久无码观看在线| 久久久久久国产精品mv| 中文无码毛片又爽又刺激| 找国产毛片看| 亚洲成综合人影院在院播放| a级毛片免费网站| 国产乱子精品一区二区在线观看| 国产女人水多毛片18| 狠狠久久综合伊人不卡| 亚洲午夜国产精品无卡| 久久免费看片| 看你懂的巨臀中文字幕一区二区| 国产一级毛片yw| 国产精选自拍| 又爽又大又黄a级毛片在线视频| 乱系列中文字幕在线视频| 日韩免费毛片| 2019年国产精品自拍不卡| 在线毛片免费| 波多野结衣久久精品| 亚洲天天更新| 日韩亚洲高清一区二区| 人妻少妇乱子伦精品无码专区毛片| 午夜欧美理论2019理论| 伊人91视频| 亚洲愉拍一区二区精品| 色视频国产| 免费毛片网站在线观看| 国外欧美一区另类中文字幕| 免费高清毛片| 亚洲三级成人| 动漫精品啪啪一区二区三区| 毛片网站在线看| 91亚洲精品国产自在现线| 老司机久久精品视频| 亚洲欧美不卡视频| 免费黄色国产视频| 国产一区二区三区日韩精品| 国产一在线观看| 国产欧美精品一区二区| 992Tv视频国产精品| 国产女人综合久久精品视| 直接黄91麻豆网站| 亚洲精品免费网站| 亚洲人成人无码www| 精品国产成人a在线观看| 欧美全免费aaaaaa特黄在线| www.狠狠| 亚洲综合久久成人AV| 国产精品永久不卡免费视频 | 88国产经典欧美一区二区三区| 色婷婷视频在线| 精品综合久久久久久97| 国产亚洲精品yxsp| 国产精品色婷婷在线观看| 她的性爱视频| 在线色国产| 日韩国产综合精选| 国产成人a在线观看视频| 色悠久久久久久久综合网伊人| 欧美在线黄| 日韩午夜片| 亚洲色成人www在线观看| 精品偷拍一区二区| 香蕉视频在线观看www| 久久精品免费国产大片| 国产区福利小视频在线观看尤物| 欧美三级日韩三级| 99ri精品视频在线观看播放|