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

    基于攝像測量法的在軌柔性結(jié)構(gòu)模態(tài)參數(shù)辨識

    2014-09-03 01:53:32高晶波張春芳
    哈爾濱工業(yè)大學學報 2014年11期
    關鍵詞:標定坐標系柔性

    許 暢,王 聰,高晶波,張春芳

    (哈爾濱工業(yè)大學航天學院,150001哈爾濱)

    在軌柔性結(jié)構(gòu)是指衛(wèi)星主體以外的柔性附件,包括太陽能帆板,天線等.這些柔性附件的自由振動具有低頻、密頻的特點[1-3].在軌時飛輪或推力器工作會引起柔性附件長時間振動,產(chǎn)生低頻擾動力矩,給衛(wèi)星姿態(tài)穩(wěn)定性帶來了較大的影響[4-6].柔性附件的模態(tài)參數(shù)對衛(wèi)星姿態(tài)控制和振動控制至關重要,但由于其尺寸較大,且在軌空間環(huán)境復雜,地面模態(tài)測試難以獲得所需精度,因此,在軌模態(tài)辨識方法受到越來越多的重視[7-9].

    傳統(tǒng)的位移/加速度傳感器由于設備復雜,會大大降低在軌結(jié)構(gòu)的可靠性,也增加了有效載荷的質(zhì)量和成本.同時,接觸性測量也會影響柔性結(jié)構(gòu)的動態(tài)性能,數(shù)字攝像測量能在實現(xiàn)非接觸測量的同時保證較高的精度,而且能降低傳感器帶來的在軌設備復雜度和成本,是非常有前景的一種在軌響應測量方法.

    攝影測量是基于專業(yè)攝像設備和數(shù)字圖像處理算法的測量技術,被廣泛應用于物體幾何尺寸和空間位置、姿態(tài)等測量.Chang等[10]對建筑結(jié)構(gòu)振動進行了攝像測量,采用數(shù)字圖像相關法進行測點特征匹配;Barrows[11]將攝像測量法用于風洞實驗模型形變的測量;Burner等[12]利用攝影測量技術進行了機翼形變測量;Ryall等[13]采用單相機測量法實現(xiàn)了對機翼的振動測量.在航天領域,NASA早在1996年就利用攝像測量技術獲得了和平號空間站太陽能帆板的在軌動力學特性[14],1998年NASA還利用攝像測量法獲得了不同環(huán)境下哈勃望遠鏡太陽帆版的變形信息[15],這些測量任務都是在航天飛機上完成的.對于大型星載天線,由于其結(jié)構(gòu)較太陽能帆板復雜,對攝像測量系統(tǒng)要求更高.國內(nèi)利用攝像測量進行結(jié)構(gòu)振動測試的文獻較少,一般只針對大型柔性附件進行有限元分析以獲得其在軌動態(tài)特性,其精度往往難以滿足要求.

    本文提出一種基于數(shù)字攝像測量的柔性星載天線的在軌模態(tài)辨識方法,利用兩臺或以上攝像機測量柔性天線在軌動態(tài)響應,如圖1所示.其關鍵技術包括相機標定,像素提取和三維重構(gòu)及辨識算法等.利用天線模型和兩臺高速攝像機進行了實驗驗證,測得了各個測點的振動曲線,并成功地辨識出了天線模型前5階固有頻率和阻尼比,為在軌柔性結(jié)構(gòu)模態(tài)參數(shù)辨識提供了一種有效的方法.

    圖1 星載柔性天線在軌響應的攝像測量

    1 雙目三維重構(gòu)技術

    1.1 三維動態(tài)重構(gòu)

    圖像三維重構(gòu)是指利用兩個或兩個以上相機對被測物進行成像,然后從圖像中獲得物體三維坐標的過程,過程中涉及4個坐標系,分別為:

    1)世界坐標系——根據(jù)自然環(huán)境所選定的坐標系,坐標值用(Xw,Yw,Zw)表示;

    2)圖像坐標系——坐標原點在CCD圖像平面的中心,X軸、Y軸分別為平行于圖像平面的兩條垂直邊,坐標值用(x,y)表示;

    3)像素坐標系——坐標原點在CCD圖像平面的左上角,X軸、Y軸分別平行于圖像坐標系的X軸和Y軸,坐標值用(u,v)來表示;

    4)光心坐標系——以相機的光心為坐標原點,X軸、Y軸分別平行于圖像坐標系的X軸和Y軸,相機的光軸為Z軸,坐標值用(Xc,Yc,Zc)表示.三維重構(gòu)的步驟是將像素坐標系轉(zhuǎn)換到圖像坐標系,再從圖像坐標系轉(zhuǎn)換到光心坐標系,最后由光心坐標系轉(zhuǎn)換到世界坐標系.最終像素坐標和世界坐標可由如下關系表述:

    式中:(u0,v0)是圖像坐標系原點在像素坐標系中的坐標;(dx,dy)分別是像素坐標系在x方向和y方向相鄰像素間的距離;f為相機焦距;R和t是旋轉(zhuǎn)矩陣和平移矩陣.可將式(1)簡化成

    式中:N為相機的內(nèi)部參數(shù)矩陣,包含了相機的焦距,像素單元距離等自身參數(shù);M為相機的外部參數(shù)矩陣,包括了相機安裝位置和角度等信息;Q為投影矩陣.通過兩個相機同時拍攝被測點,可以建立以下方程組:

    當相機內(nèi)外參數(shù)Q1,Q2已知,則上式消去Zc1和Zc2后得到3個未知數(shù)4個方程,可解得被測點的世界坐標值,理論上采用的攝像機數(shù)目越多,得到的結(jié)果越精確.利用不同時刻拍攝得到的圖像序列,就可以得到測點的位移振動曲線,完成對測點的三維動態(tài)重構(gòu).

    1.2 相機標定

    將式(2)改寫成

    相機的內(nèi)部參數(shù)包括f、(dx,dy)及(u0,v0),內(nèi)參數(shù)標定方法有很多,本文采用張正友標定法[16].當內(nèi)部參數(shù)已經(jīng)確定時,可以進行外部參數(shù)的標定,以下重點敘述用于三維重構(gòu)的相機外參數(shù)標定方法.

    從式(4)可見內(nèi)部參數(shù)矩陣N為滿秩矩陣,令:αx=dx/f;αy=dy/f,用N-1左乘式(4)的兩邊并整理得

    令k=αx(u-u0),l=αy(v-v0),M=Rt[ ],進一步將式(5)改寫成

    其中:(Xwi,Ywi,Zwi)是世界坐標系中的第i點坐標,(ki,li)是由對應的圖像像素坐標(ui,vi)的計算值.

    將式(6)展開可得到如下3個方程:

    將式(7)中的第3式分別代入上兩式消去Zci可得到如下兩個方程:

    若在世界坐標系中取n個點,將會產(chǎn)生2n個方程,用矩陣形式寫出這些方程如下:

    從式(9)可見,若已知世界坐標系中n個點的坐標,同時也知道各對應點在像素坐標系中的坐標,將式(9)兩邊同時除以m34(m34≠0)或取m34=1不影響方程的求解.式(9)可簡寫成

    式中h,U分別為式(9)等號左側(cè)向量和右側(cè)向量除以m34后的向量.當2n>11時,上式成為一個超定方程,其最小二乘解可用下式表示:

    則數(shù)碼相機的外部參數(shù)矩陣與其有如下關系:

    因R是正交矩陣,所以有

    2 圖像角點檢測

    在相機標定和三維重構(gòu)中,如何準確地得到測點的像素坐標是攝像測量的關鍵.當測點的灰度值與背景灰度值相差較大時,利用角點檢測算法可以準確地獲得測點的像素坐標,本文采用Harris角點檢測算法[17]對每一測點進行像素點提取.Harris角點檢測算法將每一像素點的相似度抽象為

    式中:Ⅰx,y表示像素點(x,y)的灰度值,x,y分別為像素點在u,v方向的變化量,且

    其中?為卷積算子.

    因此,對于比較小的位移,E可表示為

    式中:

    為了解決響應噪聲的問題,算法引入一個平滑的高斯窗,通過對w的設定實現(xiàn),其中

    由式(9)Harris角點檢測區(qū)域變化式,E可以通過矩陣形式表現(xiàn)出來:

    式中:

    Harris算法的核心是計算矩陣M的特征值.M的兩個特征值分別指示了該像素點在X,Y兩個方向上的強度.若兩個值都很大,則該像素點可認為是角點,只有一個較大而另一個較小時則認為是邊緣點.在實際操作中,算法將M的兩個特征值進行混合運算得出單一的強度值,并取局部強度值的最大點為特征點.

    綜上所述,利用數(shù)字攝像測量法測量在軌柔性結(jié)構(gòu)振動的步驟可以總結(jié)如下:

    1)分別對兩個相機進行內(nèi)外參數(shù)的標定;

    2)選擇合適的幀率對振動結(jié)構(gòu)進行動態(tài)拍攝;

    3)利用角點檢測算法獲得測點的像素坐標;

    4)利用三維重構(gòu)算法獲得測點的世界坐標;

    5)利用測點不同時刻的世界坐標得到其位移振動曲線.

    3 星載天線模型實驗

    為驗證上述方法的可行性,對一個星載柔性桁架天線模型進行了振動測試,并利用測試結(jié)果,采用特征系統(tǒng)實現(xiàn)法(ERA)[18]對其進行了模態(tài)辨識.

    星載天線模型示意圖如圖2所示,反射面直徑為750 mm,上下兩層之間距離為150 mm,兩層之間采用四根均勻分布直桿支撐,長度為150 mm.兩根長支撐桿之間距離為50 mm,支撐桿長度為700 mm,所用材料直徑均為3 mm,為了使其具有高柔性,低頻的特性,采用了低剛度材料.在天線模型外圈測點上粘貼標志點便于圖像處理,將天線支撐桿懸掛起來,兩臺高速攝像機一左一右放置于天線模型約1.5 m處,調(diào)節(jié)相機視場使其覆蓋模型上所有測點,實驗裝置相對位置關系如圖3所示.

    圖2 星載天線模型

    圖3 雙目攝像測量裝置

    實驗步驟如下:

    1)對相機進行標定時采用黑白相間的棋盤圖,如圖4所示.相機參數(shù)標定時采集標定板的多個位置的靜態(tài)圖片,找出每個圖像中的角點,根據(jù)每個角點已知的物理坐標關系,利用1.2節(jié)中的標定方法得到相機的內(nèi)外參數(shù).

    圖4 相機內(nèi)外參數(shù)標定

    2)保持相機位置不變,設置相機觸發(fā)模式與幀率,保證兩個相機同步觸發(fā).本實驗中,采用幀率150 fps進行圖像采集,即對應采樣頻率為150 Hz.

    3)給定天線模型一個初始位移,同時觸發(fā)相機進行記錄,直至振動衰減.

    4)保存圖像數(shù)據(jù),用于分析處理.

    利用相機標定方法求解出左側(cè)相機內(nèi)外參數(shù)矩陣分別為:

    右側(cè)相機內(nèi)外參數(shù)矩陣分別為:

    兩個相機在某一時刻得到的圖像如圖5、圖6所示.利用角點檢測算法獲得圖中18個測點的像素坐標.以4號測點為例,兩個相機圖像中的像素位移曲線分別如圖7~10所示.

    圖5 某一時刻左側(cè)相機圖像

    圖6 某一時刻右側(cè)相機圖像

    圖7 左側(cè)相機4號側(cè)點X方向像素振動

    圖8 左側(cè)相機4號側(cè)點Y方向像素振動

    圖9 右側(cè)相機4號側(cè)點X方向像素振動

    圖10 右側(cè)相機4號側(cè)點Y方向像素振動

    采用三維動態(tài)重構(gòu)算法得到每個測點三自由度位移振動曲線,其中4號點三軸位移響應分別如圖11~13所示.

    圖11 4號測點X方向振動位移

    圖12 4號測點Y方向振動位移

    圖13 4號測點Z方向振動位移

    采用測得的天線模型上8個測點的振動響應,利用特征系統(tǒng)實現(xiàn)(ERA)法對天線模型進行模態(tài)參數(shù)辨識,得到其前5階固有頻率和阻尼比如表1所示,振型如圖14(a)~14(e)所示,從圖中可以看出,前三階表現(xiàn)為反射器整體變形,是連接桿的彎曲和扭轉(zhuǎn)所致,從第四階開始出現(xiàn)了反射器本身的扭轉(zhuǎn)變形.從表1中還可以看出,該天線模型固有頻率較低,且第二、三階頻率接近.因此,采用攝像測量法能有效地辨識低頻、密頻在軌柔性結(jié)構(gòu)的模態(tài)參數(shù).

    表1 天線模型ERA辨識結(jié)果

    圖14 天線模型前5階振型

    4 誤差與精度分析

    對天線進行模態(tài)分析的精度取決于成像系統(tǒng)分辨力、圖像處理算法和辨識算法等眾多因素.這里主要考慮攝像測量法特有的成像系統(tǒng)分辨力.成像系統(tǒng)對空間光學圖像的最小細節(jié)的分辨能力稱為成像分辨力,一般用圖像物面分辨率或成像系統(tǒng)角分辨率來衡量.對于具體的測量對象,鏡頭焦距、光敏傳感器芯片面積和成像物距決定了成像視場大小.當待測物體物距確定時,圖像物面分辨率給出了圖像中一個像素所代表的感興趣的物面尺寸大小,因此這時用圖像物面分辨率來衡量成像系統(tǒng)分辨率真具有直觀性.而對于成像物距不定或變化的情況,對于鏡頭焦距一定的成像系統(tǒng),其視場角是確定的,此時則適合用角分辨率來衡量成像系統(tǒng)分辨率.

    每個像素所代表物面空間的尺寸稱為圖像物面分辨率.設拍攝視場大小為W×H,攝像機分辨率為M×N(稱為硬件分辨率),則圖像物面分辨率為

    本次實驗拍攝視場大小為100 mm×100 mm,分辨率為512×512,由上式可得圖像物面分辨率為0.20 mm.角點檢測算法提取像素點精度為1個像素,因此,本實驗系統(tǒng)位移測量精度0.20/2=0.1 mm.為了提高拍攝精度,可以減小拍攝視場,或者提高相機分辨率,在相機與被測物距離受限時也可以通過長焦鏡頭減小視場.

    測量中的誤差主要由圖像中存在的噪聲引起,這種噪聲不但影響測量過程,同時也會給標定帶來誤差.由于在軌和溫度環(huán)境的復雜性,圖像噪聲會比地面試驗中更大,需要采用合適的圖像處理算法濾除以獲得更精確地結(jié)果.

    5 結(jié)語

    針對大型柔性附件在軌模態(tài)辨識問題,提出了一種基于攝像測量的三維動態(tài)重構(gòu)方法.首先對相機進行內(nèi)外參數(shù)的標定,然后利用兩個相機同時記錄被測物的動態(tài)圖像,通過角點檢測算法獲得測點的像素坐標,再根據(jù)三維重構(gòu)算法得到測點的三維坐標,每一時刻測點三維坐標時間序列即是測點的三自由度位移振動曲線.利用天線模型和高速攝像機進行了地面驗證實驗,獲得了模型上各個測點的振動曲線并據(jù)此辨識出了模型的前5階模態(tài)參數(shù),包括固有頻率,阻尼比和振型.本文所述攝像測量法設備簡單,可以根據(jù)實際測量要求選擇合適的相機參數(shù),以達到精度要求.實驗表明,攝像測量法能有效地對大型柔性天線進行在軌模態(tài)辨識,且易于在軌實現(xiàn),保證了衛(wèi)星系統(tǒng)的可靠性,是一種具有良好應用前景的在軌響應測量方法.

    [1]司宏偉,李東旭,陳衛(wèi)東.大撓性航天桁架結(jié)構(gòu)動力學及其主動控制研究進展[J].力學進展,2008,38(2):167-176.

    [2]周舟,陸秋海,任革學,等.低密頻太陽能帆板動力學參數(shù)在軌辨識和振動控制[J].工程力學,2004,21(3):84-89.

    [3]黃文虎,曹登慶,韓增堯.航天器動力學與控制的研究進展與展望[J].力學進展,2012,42(4):367-394.

    [4]張逸群,段寶巖,李團結(jié).基于濾波的柔性可展開天線展開過程控制方法[J].機械工程學報,2012,48(3):180-188.

    [5]游斌弟,趙志剛,李文博,等.空間熱載荷作用下星載天線耦合動態(tài)影響分析[J].振動與沖擊,2012,31(17):61-66.

    [6]GENG L H,XIAO D Y,WANG Q,et al.Attitudecontrol model identification of on-orbit satellites actuated by reaction wheels[J].Acta Astronautica,2010,66(5):714-721.

    [7]ANTHONY T,ANDERSEN G.On-orbit modal identification of the hubble space telescope[C]//American Control Conference, 1995. Proceedings of the IEEE.Piscataway:IEEE,1995,1:402-406.

    [8]趙壽根,程偉,孫國江,等.航天器動力學特性參數(shù)在軌辨識技術[J].北京航空航天大學學報,2005,31(9):999-1003.

    [9]黨朝輝,項軍華.超低軌衛(wèi)星氣動參數(shù)及轉(zhuǎn)動慣量在軌實時辨識[J].上海航天,2012,29(5):20-24.

    [10]CHANG C C,JI Y F.Flexible videogrammetric technique for three-dimensional structural vibration measurement[J].Journal of Engineering Mechanics,2007,133(6):656-664.

    [11]BARROWS D A.Videogrammetric model deformation measurement technique for wind tunnel applications[C]//45th AIAA Aerospace Sciences Meeting and Exhitbit.Reston VA:AIAA,2007:AIAA 2007-1163.

    [12]BURNER A W,LIU T.Videogrammetric model deformation measurement technique[J].Journal of Aircraft,2001,38(4):745-754.

    [13]RYALL T G,F(xiàn)RASER C S.Determination of structural modes of vibration using digital photogrammetry[J].Journal of Aircraft,2002,39(1):114-119.

    [14]GILBERT M G,WELCH S S.STS-74/MIR photogrammetric appendage structural dynamics experiment[R].Washington DC:NASA,1996.

    [15]SAPP C A,DRAGG J L,SNYDER M W,et al.Photogrammetric assessment of the hubble space telescope solar arrays during the second servicing mission:final report[R].Washington DC:NASA,1998.

    [16]ZHANG Z.A flexible new technique for camera calibration[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2000,22(11):1330-1334.

    [17]HARRIS C,STEPHENS M.A combined corner and edge detector[C]//Proceedings of the Fourth Alvey Vision Conference.Sheftield:University of Sheffield,1988:147-151.

    [18]JUANG J N,PAPPA R S.An eigensystem realization algorithm for modal parameter identification and model reduction[J].Journal of Guidance, Control, and Dynamics,1985,8(5):620-627.

    猜你喜歡
    標定坐標系柔性
    一種柔性拋光打磨頭設計
    灌注式半柔性路面研究進展(1)——半柔性混合料組成設計
    石油瀝青(2021年5期)2021-12-02 03:21:18
    高校學生管理工作中柔性管理模式應用探索
    使用朗仁H6 Pro標定北汽紳寶轉(zhuǎn)向角傳感器
    解密坐標系中的平移變換
    坐標系背后的故事
    基于勻速率26位置法的iIMU-FSAS光纖陀螺儀標定
    基于重心坐標系的平面幾何證明的探討
    船載高精度星敏感器安裝角的標定
    基于Harris-張正友平面標定法的攝像機標定算法
    午夜视频精品福利| 亚洲伊人色综图| 97人妻天天添夜夜摸| 午夜精品国产一区二区电影| 高清视频免费观看一区二区| 色精品久久人妻99蜜桃| 日韩中文字幕视频在线看片| 亚洲中文日韩欧美视频| 狠狠精品人妻久久久久久综合| 亚洲视频免费观看视频| 天天躁夜夜躁狠狠久久av| 99热网站在线观看| 亚洲精品久久午夜乱码| 在线观看www视频免费| 国产成人系列免费观看| 精品国产一区二区三区久久久樱花| 精品人妻熟女毛片av久久网站| 叶爱在线成人免费视频播放| 最近最新中文字幕大全免费视频 | 侵犯人妻中文字幕一二三四区| 97在线人人人人妻| 色网站视频免费| 午夜av观看不卡| 国产精品欧美亚洲77777| 一本综合久久免费| 欧美精品高潮呻吟av久久| 亚洲av日韩精品久久久久久密 | 久久精品国产综合久久久| 亚洲精品中文字幕在线视频| 欧美黑人欧美精品刺激| 日本av免费视频播放| svipshipincom国产片| 2021少妇久久久久久久久久久| 男人舔女人的私密视频| 午夜福利免费观看在线| 国产精品秋霞免费鲁丝片| 青青草视频在线视频观看| 九色亚洲精品在线播放| 91麻豆精品激情在线观看国产 | 视频区图区小说| 欧美日韩视频高清一区二区三区二| 狂野欧美激情性bbbbbb| 精品国产一区二区三区四区第35| 国产精品av久久久久免费| 麻豆av在线久日| 十八禁人妻一区二区| 亚洲av美国av| av欧美777| 啦啦啦视频在线资源免费观看| 这个男人来自地球电影免费观看| 久久亚洲国产成人精品v| 国产无遮挡羞羞视频在线观看| 亚洲欧美日韩另类电影网站| 欧美国产精品一级二级三级| 久久久久久久久久久久大奶| 欧美成人午夜精品| 久久人人爽av亚洲精品天堂| 可以免费在线观看a视频的电影网站| 看免费成人av毛片| 97人妻天天添夜夜摸| 两个人免费观看高清视频| av片东京热男人的天堂| 最近最新中文字幕大全免费视频 | 最黄视频免费看| 美女中出高潮动态图| 飞空精品影院首页| 久久久久久久大尺度免费视频| 欧美日韩黄片免| 一边摸一边做爽爽视频免费| 亚洲精品国产av蜜桃| 国产精品免费视频内射| 飞空精品影院首页| 搡老岳熟女国产| 免费少妇av软件| 在线观看一区二区三区激情| 精品一品国产午夜福利视频| 国产精品久久久久久人妻精品电影 | 日韩中文字幕欧美一区二区 | 久久青草综合色| 国产亚洲一区二区精品| 亚洲视频免费观看视频| 免费高清在线观看视频在线观看| 久久久久久久久免费视频了| 老司机靠b影院| 欧美黄色片欧美黄色片| 亚洲av国产av综合av卡| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品国产区一区二| 日本黄色日本黄色录像| 视频区图区小说| 免费在线观看日本一区| 老司机靠b影院| 亚洲成av片中文字幕在线观看| 99国产精品一区二区蜜桃av | 午夜影院在线不卡| 99国产精品一区二区蜜桃av | 国产成人系列免费观看| 免费女性裸体啪啪无遮挡网站| 亚洲av成人精品一二三区| 看免费av毛片| 国产深夜福利视频在线观看| 精品国产乱码久久久久久男人| 日韩一卡2卡3卡4卡2021年| 国产成人a∨麻豆精品| 美女视频免费永久观看网站| 性高湖久久久久久久久免费观看| 国产国语露脸激情在线看| 国产精品欧美亚洲77777| 人人妻人人爽人人添夜夜欢视频| 久久亚洲精品不卡| 欧美性长视频在线观看| 在线观看免费高清a一片| 狂野欧美激情性xxxx| 欧美日韩福利视频一区二区| 日韩中文字幕视频在线看片| 久久人人爽av亚洲精品天堂| 人妻人人澡人人爽人人| 侵犯人妻中文字幕一二三四区| 欧美亚洲日本最大视频资源| av在线app专区| 免费看不卡的av| 久久精品人人爽人人爽视色| 欧美在线一区亚洲| 国产色视频综合| 可以免费在线观看a视频的电影网站| av天堂久久9| av国产久精品久网站免费入址| 一边摸一边做爽爽视频免费| 啦啦啦 在线观看视频| 亚洲av国产av综合av卡| av在线播放精品| 国产爽快片一区二区三区| 大香蕉久久网| 久久久久久亚洲精品国产蜜桃av| 国产精品一区二区精品视频观看| 欧美人与善性xxx| 亚洲熟女精品中文字幕| 91九色精品人成在线观看| 亚洲图色成人| 免费女性裸体啪啪无遮挡网站| 中文字幕色久视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲 国产 在线| 欧美黄色片欧美黄色片| 中文欧美无线码| av福利片在线| xxx大片免费视频| 国产精品一区二区在线不卡| 王馨瑶露胸无遮挡在线观看| 国产亚洲一区二区精品| 国产精品香港三级国产av潘金莲 | www.精华液| 中国国产av一级| 高清不卡的av网站| 狂野欧美激情性bbbbbb| 18禁观看日本| 亚洲国产精品999| 午夜福利免费观看在线| 电影成人av| 99香蕉大伊视频| 首页视频小说图片口味搜索 | 丁香六月天网| 久久鲁丝午夜福利片| 亚洲欧美精品自产自拍| 国产精品秋霞免费鲁丝片| 高潮久久久久久久久久久不卡| 狠狠精品人妻久久久久久综合| 欧美日韩综合久久久久久| xxx大片免费视频| 97精品久久久久久久久久精品| 午夜福利一区二区在线看| 成人午夜精彩视频在线观看| 99久久精品国产亚洲精品| 日本午夜av视频| 18禁裸乳无遮挡动漫免费视频| 一区二区三区四区激情视频| 又紧又爽又黄一区二区| 国产精品久久久av美女十八| av有码第一页| 国产精品成人在线| www.精华液| 超碰成人久久| 自线自在国产av| 国产在线视频一区二区| av网站免费在线观看视频| 一本—道久久a久久精品蜜桃钙片| 91精品伊人久久大香线蕉| 超碰97精品在线观看| 国产成人精品在线电影| 亚洲精品国产色婷婷电影| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品一区二区三区在线| 91老司机精品| 日本欧美视频一区| 一区在线观看完整版| 天天躁日日躁夜夜躁夜夜| 在线观看免费视频网站a站| 欧美成狂野欧美在线观看| 啦啦啦视频在线资源免费观看| tube8黄色片| 免费看十八禁软件| 啦啦啦在线免费观看视频4| 久久久国产欧美日韩av| 久久精品亚洲av国产电影网| 老司机深夜福利视频在线观看 | 女人久久www免费人成看片| 首页视频小说图片口味搜索 | 别揉我奶头~嗯~啊~动态视频 | 国产亚洲av片在线观看秒播厂| 在线看a的网站| 久久久久久久久免费视频了| 少妇精品久久久久久久| 成人黄色视频免费在线看| 女人爽到高潮嗷嗷叫在线视频| 精品国产一区二区三区四区第35| 女警被强在线播放| 麻豆乱淫一区二区| 国产亚洲精品久久久久5区| 国产高清videossex| 777久久人妻少妇嫩草av网站| 欧美日韩国产mv在线观看视频| 午夜福利在线免费观看网站| 母亲3免费完整高清在线观看| 亚洲黑人精品在线| 欧美日韩视频精品一区| 99热国产这里只有精品6| 99国产精品一区二区三区| 黄片小视频在线播放| 午夜影院在线不卡| 欧美国产精品一级二级三级| cao死你这个sao货| 久久久国产欧美日韩av| 女性被躁到高潮视频| 色婷婷av一区二区三区视频| 赤兔流量卡办理| 国产老妇伦熟女老妇高清| 性少妇av在线| 久久中文字幕一级| 涩涩av久久男人的天堂| 国产日韩一区二区三区精品不卡| 日韩人妻精品一区2区三区| videos熟女内射| 中文字幕最新亚洲高清| 香蕉丝袜av| 91字幕亚洲| 老司机亚洲免费影院| 国产一区亚洲一区在线观看| 精品少妇黑人巨大在线播放| 中文字幕精品免费在线观看视频| 欧美成狂野欧美在线观看| 久久av网站| 91九色精品人成在线观看| 性色av乱码一区二区三区2| 亚洲五月色婷婷综合| 美女主播在线视频| 亚洲国产精品999| 好男人视频免费观看在线| 日韩大码丰满熟妇| 午夜福利一区二区在线看| 飞空精品影院首页| 亚洲欧洲精品一区二区精品久久久| 欧美日本中文国产一区发布| 脱女人内裤的视频| 日本色播在线视频| 99re6热这里在线精品视频| 国产1区2区3区精品| 国产熟女欧美一区二区| 欧美性长视频在线观看| 人人妻人人澡人人看| 精品人妻熟女毛片av久久网站| 少妇猛男粗大的猛烈进出视频| 汤姆久久久久久久影院中文字幕| 亚洲国产精品999| 欧美日韩黄片免| 2018国产大陆天天弄谢| 久久鲁丝午夜福利片| 亚洲中文字幕日韩| 激情视频va一区二区三区| 2021少妇久久久久久久久久久| 嫩草影视91久久| 看十八女毛片水多多多| 国产成人av教育| 高潮久久久久久久久久久不卡| 在线亚洲精品国产二区图片欧美| 久久女婷五月综合色啪小说| 日韩大码丰满熟妇| 国产99久久九九免费精品| 国产片特级美女逼逼视频| 极品少妇高潮喷水抽搐| 99精国产麻豆久久婷婷| 国产男人的电影天堂91| av不卡在线播放| 91成人精品电影| 亚洲人成网站在线观看播放| 性色av乱码一区二区三区2| 黑丝袜美女国产一区| 制服诱惑二区| 国产一区二区 视频在线| 一本大道久久a久久精品| 免费黄频网站在线观看国产| 99国产综合亚洲精品| 亚洲av男天堂| 日本欧美国产在线视频| 黄片播放在线免费| 欧美成人午夜精品| 久久久久精品人妻al黑| 亚洲中文av在线| 男人操女人黄网站| 免费在线观看日本一区| 亚洲欧美激情在线| 婷婷成人精品国产| 在线 av 中文字幕| 又黄又粗又硬又大视频| 成人影院久久| 99国产精品99久久久久| 黄色毛片三级朝国网站| 国产视频首页在线观看| 久久久国产一区二区| 后天国语完整版免费观看| 性色av乱码一区二区三区2| 精品一品国产午夜福利视频| 久久精品亚洲熟妇少妇任你| av电影中文网址| av在线播放精品| 2021少妇久久久久久久久久久| netflix在线观看网站| 校园人妻丝袜中文字幕| e午夜精品久久久久久久| 大香蕉久久成人网| 制服诱惑二区| 精品熟女少妇八av免费久了| 国产一区二区三区av在线| 欧美成狂野欧美在线观看| 欧美日韩成人在线一区二区| 高清黄色对白视频在线免费看| 国产av精品麻豆| 爱豆传媒免费全集在线观看| 国产精品熟女久久久久浪| www日本在线高清视频| 少妇猛男粗大的猛烈进出视频| 国产在视频线精品| 日韩 欧美 亚洲 中文字幕| 汤姆久久久久久久影院中文字幕| 两个人看的免费小视频| 久久热在线av| 韩国精品一区二区三区| 亚洲精品日韩在线中文字幕| 一级毛片女人18水好多 | 嫁个100分男人电影在线观看 | av在线老鸭窝| 丝袜美腿诱惑在线| 久久影院123| 美女国产高潮福利片在线看| 免费黄频网站在线观看国产| 亚洲一区二区三区欧美精品| 日日爽夜夜爽网站| 新久久久久国产一级毛片| 十八禁人妻一区二区| 中文字幕精品免费在线观看视频| 一本一本久久a久久精品综合妖精| 老司机午夜十八禁免费视频| 老司机在亚洲福利影院| 岛国毛片在线播放| 又粗又硬又长又爽又黄的视频| 亚洲精品在线美女| 香蕉丝袜av| 久久久久久久久免费视频了| 乱人伦中国视频| 国产免费福利视频在线观看| 久久ye,这里只有精品| 丝袜在线中文字幕| 亚洲,欧美精品.| 一区在线观看完整版| 99久久99久久久精品蜜桃| 天天影视国产精品| 18在线观看网站| 久久 成人 亚洲| 国产真人三级小视频在线观看| 国产色视频综合| 国产在视频线精品| 国产精品国产av在线观看| 别揉我奶头~嗯~啊~动态视频 | 不卡av一区二区三区| a级片在线免费高清观看视频| 国产成人一区二区三区免费视频网站 | 亚洲九九香蕉| 两个人看的免费小视频| 国产真人三级小视频在线观看| 男人操女人黄网站| 日韩 亚洲 欧美在线| 美女脱内裤让男人舔精品视频| 国产高清videossex| 亚洲国产精品国产精品| 麻豆av在线久日| 亚洲精品av麻豆狂野| 国产高清视频在线播放一区 | 欧美亚洲日本最大视频资源| 人体艺术视频欧美日本| 久久影院123| 亚洲精品在线美女| 久久鲁丝午夜福利片| 精品久久久久久电影网| 一边摸一边做爽爽视频免费| 人体艺术视频欧美日本| 免费久久久久久久精品成人欧美视频| 欧美日韩亚洲高清精品| 亚洲av欧美aⅴ国产| 国产99久久九九免费精品| 午夜福利乱码中文字幕| 大香蕉久久成人网| av不卡在线播放| 日本黄色日本黄色录像| 中文精品一卡2卡3卡4更新| 亚洲一区中文字幕在线| 亚洲国产毛片av蜜桃av| 人体艺术视频欧美日本| 18禁裸乳无遮挡动漫免费视频| 欧美乱码精品一区二区三区| 精品人妻1区二区| 黑丝袜美女国产一区| 99热网站在线观看| 日韩免费高清中文字幕av| 97在线人人人人妻| 一级片'在线观看视频| 亚洲欧洲日产国产| 侵犯人妻中文字幕一二三四区| av欧美777| 纵有疾风起免费观看全集完整版| 日韩中文字幕视频在线看片| 欧美 亚洲 国产 日韩一| 国产亚洲精品第一综合不卡| 丝袜美腿诱惑在线| 国产精品亚洲av一区麻豆| 国语对白做爰xxxⅹ性视频网站| 丝瓜视频免费看黄片| 成人免费观看视频高清| 国产国语露脸激情在线看| 天天躁日日躁夜夜躁夜夜| 黄片小视频在线播放| 啦啦啦在线观看免费高清www| 免费在线观看影片大全网站 | 日本猛色少妇xxxxx猛交久久| 狠狠精品人妻久久久久久综合| 亚洲av欧美aⅴ国产| 天天添夜夜摸| 日本vs欧美在线观看视频| 亚洲专区国产一区二区| 国产男女超爽视频在线观看| 国产av国产精品国产| 国产欧美日韩一区二区三 | 亚洲七黄色美女视频| 制服人妻中文乱码| 女人高潮潮喷娇喘18禁视频| 只有这里有精品99| 最近手机中文字幕大全| 赤兔流量卡办理| 亚洲人成电影免费在线| 美女国产高潮福利片在线看| 欧美人与善性xxx| av又黄又爽大尺度在线免费看| 在线亚洲精品国产二区图片欧美| 亚洲国产欧美一区二区综合| 欧美精品亚洲一区二区| 日韩大码丰满熟妇| 亚洲成人免费电影在线观看 | 亚洲成人手机| 国产一卡二卡三卡精品| 国产97色在线日韩免费| 欧美日韩亚洲高清精品| 亚洲av国产av综合av卡| 老司机深夜福利视频在线观看 | 国产男女超爽视频在线观看| 国产精品一区二区免费欧美 | 久久精品熟女亚洲av麻豆精品| 国产麻豆69| 老汉色∧v一级毛片| 日本a在线网址| 咕卡用的链子| 成年人免费黄色播放视频| 国产亚洲午夜精品一区二区久久| av在线app专区| 后天国语完整版免费观看| 水蜜桃什么品种好| 狠狠婷婷综合久久久久久88av| 精品第一国产精品| 精品视频人人做人人爽| 日本wwww免费看| 麻豆乱淫一区二区| 纵有疾风起免费观看全集完整版| 欧美乱码精品一区二区三区| 捣出白浆h1v1| 国产视频首页在线观看| avwww免费| 精品久久蜜臀av无| 嫩草影视91久久| 美女主播在线视频| 久久人人爽人人片av| 交换朋友夫妻互换小说| 午夜91福利影院| 国产欧美日韩精品亚洲av| 精品一区二区三区av网在线观看 | 波多野结衣一区麻豆| 久久久久国产精品人妻一区二区| 国产有黄有色有爽视频| av国产久精品久网站免费入址| 久久鲁丝午夜福利片| 国产日韩一区二区三区精品不卡| 久久人人97超碰香蕉20202| 国产精品久久久久久人妻精品电影 | 老司机影院成人| 久久精品国产亚洲av高清一级| 一级毛片我不卡| 看十八女毛片水多多多| 国产成人免费无遮挡视频| 国产三级黄色录像| 一级毛片女人18水好多 | 午夜久久久在线观看| 国产亚洲精品第一综合不卡| 精品国产乱码久久久久久小说| 建设人人有责人人尽责人人享有的| 亚洲av国产av综合av卡| 亚洲午夜精品一区,二区,三区| 精品国产一区二区三区久久久樱花| 久久国产亚洲av麻豆专区| 欧美激情 高清一区二区三区| 亚洲综合色网址| 免费黄频网站在线观看国产| 国产成人影院久久av| 五月天丁香电影| 免费观看av网站的网址| 午夜久久久在线观看| 精品国产一区二区三区四区第35| 视频在线观看一区二区三区| 国产99久久九九免费精品| av欧美777| 精品久久蜜臀av无| av线在线观看网站| 午夜av观看不卡| 精品福利观看| 国产精品久久久久久人妻精品电影 | 老司机影院成人| 在线天堂中文资源库| 青春草亚洲视频在线观看| 国产高清不卡午夜福利| 精品欧美一区二区三区在线| 乱人伦中国视频| 嫁个100分男人电影在线观看 | tube8黄色片| 亚洲国产av新网站| 亚洲 国产 在线| 国产一区二区三区av在线| 黄网站色视频无遮挡免费观看| 大话2 男鬼变身卡| 咕卡用的链子| 精品久久久久久久毛片微露脸 | 一级,二级,三级黄色视频| 夜夜骑夜夜射夜夜干| 亚洲av日韩精品久久久久久密 | 天天躁日日躁夜夜躁夜夜| 国产高清不卡午夜福利| 久久青草综合色| 亚洲,欧美,日韩| 日本av手机在线免费观看| 国产精品久久久人人做人人爽| 国产精品免费大片| 精品一区二区三卡| 国产亚洲欧美精品永久| 欧美97在线视频| 国产精品人妻久久久影院| 天天影视国产精品| 亚洲国产欧美在线一区| 欧美av亚洲av综合av国产av| 大码成人一级视频| 天天影视国产精品| www.精华液| 国产精品国产av在线观看| 激情五月婷婷亚洲| 18禁国产床啪视频网站| 人人澡人人妻人| 国产欧美日韩精品亚洲av| 丰满迷人的少妇在线观看| 高清视频免费观看一区二区| 两个人看的免费小视频| 国产精品二区激情视频| 精品一区二区三区av网在线观看 | 永久免费av网站大全| 99九九在线精品视频| 日韩,欧美,国产一区二区三区| 高潮久久久久久久久久久不卡| 一区二区日韩欧美中文字幕| 国产片内射在线| 精品亚洲成国产av| 菩萨蛮人人尽说江南好唐韦庄| 欧美+亚洲+日韩+国产| 欧美精品一区二区免费开放| 国产亚洲av高清不卡| 一区二区三区精品91| 超碰97精品在线观看| 又大又黄又爽视频免费| www.自偷自拍.com| 久久久久久亚洲精品国产蜜桃av| 国产又爽黄色视频| 欧美国产精品一级二级三级| 国产精品国产av在线观看| 青春草亚洲视频在线观看| 成人亚洲欧美一区二区av| 亚洲欧美激情在线| 亚洲成国产人片在线观看| 欧美人与善性xxx| 欧美精品啪啪一区二区三区 | 女警被强在线播放| 午夜免费男女啪啪视频观看| 天天躁日日躁夜夜躁夜夜| 国产一区二区三区av在线| 男女无遮挡免费网站观看|