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

    一種改進(jìn)的特征系統(tǒng)實(shí)現(xiàn)算法及其在海洋平臺上的應(yīng)用

    2016-01-28 03:44:50辛峻峰張永波盛進(jìn)路
    振動與沖擊 2015年24期

    辛峻峰, 張永波, 盛進(jìn)路

    (1.青島科技大學(xué) 機(jī)電學(xué)院,山東 青島 266061;2.青島國家海洋科學(xué)研究中心,山東 青島 266071;3.重慶交通大學(xué) 交通運(yùn)輸學(xué)院,重慶 400074;4.大連海事大學(xué) 船機(jī)修造工程交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,大連 116026)

    ?

    一種改進(jìn)的特征系統(tǒng)實(shí)現(xiàn)算法及其在海洋平臺上的應(yīng)用

    辛峻峰1, 張永波2, 盛進(jìn)路3,4

    (1.青島科技大學(xué) 機(jī)電學(xué)院,山東 青島266061;2.青島國家海洋科學(xué)研究中心,山東 青島266071;3.重慶交通大學(xué) 交通運(yùn)輸學(xué)院,重慶400074;4.大連海事大學(xué) 船機(jī)修造工程交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,大連116026)

    有效的模態(tài)識別方法能精確地得到結(jié)構(gòu)的相關(guān)參數(shù),從而可以準(zhǔn)確地掌握大型結(jié)構(gòu)的健康狀況。國內(nèi)外許多學(xué)者提出了基于時域響應(yīng)的模態(tài)參數(shù)識別方法,如復(fù)指數(shù)法(Complex Exponential Method, CEM)、時間序列法(Time Series Method,TSM)、隨機(jī)減量法(Random Decrement Technique, RDT)、隨機(jī)子空間法(Stochastic Subspace Identification Method, SSI)等。其中,ERA是目前為止沖擊和張拉激勵下較為先進(jìn)高效的模態(tài)參數(shù)識別方法之一。它能夠更精確地提取海洋平臺等大型結(jié)構(gòu)的模態(tài)參數(shù)[1-3]。

    ERA是一種多輸入多輸出的時域模態(tài)參數(shù)識別技術(shù)。1984年Juang[1]首先將其應(yīng)用到結(jié)動力學(xué)領(lǐng)域,目前在國外航天領(lǐng)域廣泛使用。但是當(dāng)數(shù)據(jù)的信噪比較低時,ERA方法的運(yùn)用受到限制。針對于此,本文首先基于矩陣映射理論提出了改進(jìn)的ERA算法-C/ERA,繼而通過五自由度阻尼彈簧質(zhì)量系統(tǒng)仿真和導(dǎo)管架平臺縮尺模型的物理模型試驗(yàn)將ERA和C/ERA兩種方法進(jìn)行了對比。

    1C/ERA基本理論

    結(jié)構(gòu)的狀態(tài)可用下面的微分方程表示(m表示輸入數(shù)量、n表示輸出數(shù)量)[1,4]:

    y(k)=Cx(k)

    (1)

    式中:x(k)x(k)為kΔt時刻系統(tǒng)的狀態(tài)矢量(2n維); y(k)為kΔt時刻系統(tǒng)的響應(yīng)矢量(n維);p(k)為kΔt時刻系統(tǒng)的激勵矢量(m維)。A是系統(tǒng)矩陣, B是控制矩陣,C是輸出矩陣;其中Δt為采樣時間間隔。

    基于脈沖響應(yīng)函數(shù)矩陣序列(Markov參數(shù))構(gòu)建廣義Hankel矩陣

    H(k-1)=

    (2)

    (3)

    A=S-1/2UTH(1)VS-1/2

    B=S-1/2VTE

    C=ETmUS1/2

    (4)

    2數(shù)值仿真

    基于Matlab建立了一個五自由度質(zhì)量-彈簧-阻尼系統(tǒng)數(shù)值模型[7],如圖1所示。

    圖1 自由度彈簧-質(zhì)量-阻尼系統(tǒng)Fig.1 A 5-DOF mass-spring-dashpot system

    每個單元的質(zhì)量、剛度和阻尼系數(shù)分別取為:mn= 50 kg,kn= 2.9 × 107N/m,cn= 1 000 N·s/m。位移記為xn,其中n=1,…,5。通過特征值分析,得到5階模態(tài)頻率的理論值為:34.499、100.700、158.730、203.880、232.520 Hz;5階模態(tài)阻尼比的理論值為:0.003 737、0.010 909、0.017 197、0.022 092、0.025 198[3]。

    沖擊荷載是用Matlab impulse函數(shù)生成的單位脈沖(見圖2),繼而將其加載到靠固定端的質(zhì)量塊上 (見圖1)。

    圖2 單位脈沖Fig.2 Unit impulse

    本文采用靠固定端質(zhì)量塊輸出的脈沖響應(yīng)作為待用數(shù)據(jù),并加入20%方差為1.0103的隨機(jī)噪聲(見圖3)。

    圖3 靠固定端質(zhì)量塊的響應(yīng)Fig.3 Response signal of the first mass

    在方法的使用方式上,本文將ERA和C/ERA的Hankel矩陣維數(shù)都設(shè)定為i=512。

    在展示工具的選擇上,本文采用了穩(wěn)定圖,其中采用了信號的傅里葉變換作為背景,顯示模態(tài)階數(shù)在0到20或25之間模態(tài)識別結(jié)果的情況。對于同一模態(tài),如果前后兩次識別結(jié)果同一模態(tài)的模態(tài)頻率誤差在1%,同時響應(yīng)阻尼的識別結(jié)果的誤差在5%,即

    則此次識別的結(jié)果標(biāo)識為穩(wěn)定,否則為不穩(wěn)定。

    圖4 ERA(a)和 C/ERA (b)沖擊激勵下穩(wěn)定圖(*穩(wěn)定,o不穩(wěn)定)Fig.4 Stability Diagram for ERA (a) and C/ERA(b) under impulse loading(*stable,o unstable)

    從圖4中能夠很容易發(fā)現(xiàn),C/ERA(圖4b)的識別結(jié)果明顯優(yōu)于ERA(圖4a)的識別結(jié)果:

    (1) C/ERA(圖4b)的穩(wěn)定圖非常干凈,而ERA(圖4a)的識別結(jié)果非常雜亂(橫軸50~100 Hz之間出現(xiàn)了虛假模態(tài),橫軸200 Hz之后縱軸12~20的區(qū)域內(nèi)非?;靵y),這說明C/ERA比ERA有更強(qiáng)的消噪能力。

    (2) ERA(圖4a)只能識別出前三階模態(tài)(橫軸在34.499 Hz、100.700 Hz、158.730 Hz處對應(yīng)定的紅色直線,203.880 Hz處對應(yīng)的直線在縱軸17開始發(fā)生混亂,232.520 Hz處沒有對應(yīng)的穩(wěn)定的直線),而C/ERA(圖b)清晰地識別出了0~250 Hz之間的5個模態(tài)(橫軸在34.499 Hz、100.700 Hz、158.730 Hz、203.880 Hz、232.520 Hz處對應(yīng)穩(wěn)定的紅色直線)。這說明C/ERA能夠比ERA識別出更多的模態(tài),尤其對高階模態(tài)有更強(qiáng)的識別能力。

    穩(wěn)定圖的對比結(jié)果清晰地展示了:較ERA而言,C/ERA消噪能力更強(qiáng)從而能夠識別出更多的高階模態(tài)。

    表1 頻率估計值(定階為10)

    表1和2分別羅列了兩種方法模態(tài)頻率和模態(tài)阻尼的估計結(jié)果,可以看到C/ERA不但能識別出第4和第5階高階模態(tài)(203.880和232.520 Hz),而且在低階模態(tài)阻尼的估計上有較大改善(相對誤差超過5%):第一階模態(tài)阻尼相對誤差由-42.52%降到-2.19%,第二階模態(tài)阻尼由7.48%降到-1.6%。

    表2 阻尼估計值的均值(定階為10)

    總之,五自由度質(zhì)量-彈簧-阻尼系統(tǒng)沖擊激勵仿真證明:較/ERA 而言,C/ERA消噪能力更強(qiáng),從而導(dǎo)致該算法對高階模態(tài)有更強(qiáng)的識別能力,并且對低階模態(tài)的阻尼有更精確的識別結(jié)果。

    3導(dǎo)管架平臺試驗(yàn)

    3.1目標(biāo)平臺

    圖5 目標(biāo)平臺Fig.5 Target platform

    某井組平臺位于東經(jīng)118°50′57″、北緯38°14′10″,是一座6井式采油、修井一體化平臺,工藝平臺為四腿式樁基鋼制導(dǎo)管架結(jié)構(gòu)。井口平臺與計量平臺之間以鋼制衍架結(jié)構(gòu)棧橋連接。海洋平臺如圖5所示。

    3.2物理模型

    圖6 海洋平臺物理模型Fig.6 Physical model of jacket platform

    導(dǎo)管架物理模型采用鋼管制作,在Y向設(shè)3根主梁,X向設(shè)2根,有3層水平橫撐,分別位于-88 cm、142 cm和188 cm處設(shè)水平橫撐,在88 cm至-142 cm處設(shè)豎向斜撐,主要構(gòu)件的尺寸為:主腿25×2.5 mm,水平撐桿16×1.5 mm,斜撐桿16×1.5 mm。模型的構(gòu)件尺寸見表3,平臺物理模型見圖6~圖7所示,對應(yīng)數(shù)值模型的模態(tài)信息見表4所示。

    圖7 海洋平臺物理模型實(shí)際比例Fig.7 Actual proportion of physical model of jacket platform

    主要構(gòu)件原型尺寸/mm模型尺寸/mm主腿1350×2425×2.5水平撐桿700×2216×1.5斜撐桿500×1816×1.5

    表4 有限元結(jié)果

    3.3沖擊激勵加載

    用橡皮錘在選定目標(biāo)平臺甲板上左側(cè)節(jié)點(diǎn)(見圖8),在第0 s施加6 kN的沖擊力,記錄平臺各個節(jié)點(diǎn)的自由振動加速度響應(yīng)。

    3.4響應(yīng)數(shù)據(jù)獲取

    以500 Hz采樣頻率用傳感器A采集了1 024個數(shù)據(jù)點(diǎn),約2 s(見圖9)。同樣將ERA和C/ERA的Hankel矩陣維數(shù)都設(shè)定為i=512。

    從穩(wěn)定圖的分析結(jié)果中(圖10)能夠發(fā)現(xiàn):

    (1) C/ERA(圖b)的穩(wěn)定圖更清晰,而ERA(圖a)的識別結(jié)果在橫軸130Hz之后縱軸10-25的區(qū)域內(nèi)非常雜亂,這說明C/ERA比ERA有更強(qiáng)的消噪能力。

    (2) ERA(圖a)和C/ERA(圖b)都能識別出前五階模態(tài)(橫軸大約在22 Hz、27 Hz、46 Hz、59 Hz、137 Hz處對應(yīng)穩(wěn)定的紅色直線),但是在高階模態(tài)的識別上(橫軸150 Hz之后),C/ERA能非常穩(wěn)定地識別出三個模態(tài)(橫軸大約在160 Hz、180 Hz、190 Hz處),而ERA在同樣的區(qū)域相對混亂,這證明C/ERA對高階模態(tài)有更強(qiáng)的識別能力。

    圖9 靠固定端質(zhì)量塊響應(yīng)Fig.9Responseofmassclosetofixedpoint圖10 物理模型試驗(yàn)ERA(a)和C/ERA(b)沖擊激勵下穩(wěn)定圖(*穩(wěn)定,o不穩(wěn)定)Fig.10StabilityDiagramforERA(a)andC/ERA(b)fromphysicalmodelofjacketplatformunderimpulseloading(*stable,ounstable)

    穩(wěn)定圖的對比結(jié)果清晰地展示了:較ERA而言,C/ERA消噪能力能強(qiáng),由此該方法能夠識別出更多的模態(tài),尤其對高階模態(tài)有更強(qiáng)的識別能力。

    在模態(tài)頻率(表5)的估計上,跟前6階模態(tài)的識別精度相比較,C/ERA對高階模態(tài) (y-4和θ-4兩個模態(tài))的識別誤差(-7.02%和3.91%)并沒有出現(xiàn)較大變化,保持了相對穩(wěn)定,由此可見該方法有能力穩(wěn)定地估計高階模態(tài)頻率。

    表5 頻率估計值(定階為20).

    總之,導(dǎo)管架平臺物理模型試驗(yàn)沖擊激勵試驗(yàn)與數(shù)值仿真的研究結(jié)論基本一致:C/ERA消噪能力更強(qiáng)從而導(dǎo)致該方法對高階模態(tài)更加敏感,能夠識別出更多的高階模態(tài)。

    4結(jié)論

    本文提出了傳統(tǒng)ERA的改進(jìn)方法-C/ERA,該算法在ERA的奇異值截斷后,重建Hankel矩陣,繼而將矩陣中各元素由其所在的反對角線上的元素數(shù)學(xué)平均值代替并引入Frobenius范數(shù)準(zhǔn)則控制迭代次數(shù)。最終通過五自由度阻尼彈簧質(zhì)量系統(tǒng)仿真和導(dǎo)管架平臺縮尺模型的物理模型錘擊試驗(yàn)驗(yàn)證了該方法的優(yōu)勢:①更強(qiáng)的消噪能力;②對高階模態(tài)更敏感;③對低階模態(tài)阻尼有更精確的識別能力。由此,本文提出的C/ERA方法比ERA更適合作為今后海洋平臺等大型結(jié)構(gòu)沖擊激勵下的模態(tài)參數(shù)識別工具。

    參 考 文 獻(xiàn)

    [1] Juang J N,Pappa R S. An eigensystem realizationm algorithm (ERA) for modal parameter identification and model reduction[J]. Journal of Guidance, Control, and Dynamics. 1985, 8(5): 620-627.

    [2] Hu S L J, Li P, Vincent H, et al. Modal parameter estimation for jacket-type platforms using free-vibration data[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 2011:1943-1460.

    [3] 辛峻峰,盛進(jìn)路,張永波.數(shù)據(jù)驅(qū)動隨機(jī)子空間法矩陣維數(shù)選擇與噪聲問題研究[J]. 振動與沖擊,2013,32(16):152-157.

    XIN Jun-feng,SHENG Jin-lu,ZHANG Yong-bo.Relation between noise and matrix dimension of data-driven stochostic subspace identification method[J]. Journal of Vibration and Shock, 2013,32(16):152-157.

    [4] Cunha A, Caetano E. Experimental modal analysis of civil engineering structures[J]. Sound and Vibration,2006:12-20.

    [5] Juang J N,Pappa R S. Effects of noise on modal parameters identified by the eigensystem realization algorithm[J]. Journal of Guidance, Control, and Dynamics,1986, 9(3):294-303.

    [6] Cadzow J A. Signal enhancement-A composite property mapping algorithm[J]. IEEE Trans Acoust, Speech, Signal Proeess,1988,36(1):49-62.

    [7] Chu M T, Funderlic R E, Plemmons R J. Structured low rank approximation[J].Linear Algebra Appl.June2003,366:157-172.

    [8] ANSYS.Inc,ed.,2004. ANSYS Academic Research.Release 9.0. Help System. Coupled Field Analysis Guide10.

    第一作者 辛峻峰 男,博士,講師,1982年12月生

    摘要:基于矩陣映射理論,提出一種ERA的改進(jìn)算法--C/ERA。該算法在傳統(tǒng)的ERA方法奇異值截斷過程之后重建Hankel矩陣,繼而將矩陣各元素由其所在反對角線上的元素數(shù)學(xué)平均值代替并引入Frobenius第二范數(shù)準(zhǔn)則控制迭代次數(shù)。五自由度阻尼彈簧質(zhì)量系統(tǒng)仿真以及海洋平臺縮尺模型物理模型的試驗(yàn)結(jié)果表明: C/ERA方法能更強(qiáng)的消噪能力,對高階模態(tài)更敏感而且在低階模態(tài)阻尼的識別上有更高的識別精度。

    關(guān)鍵詞:模態(tài)參數(shù)識別; C/ERA;消噪

    Improved eigen system realization algorithm and its application in the analysis of ocean platform

    XINJun-feng1,ZHANGYong-bo2,SHENGJin-lu3,4(1. Qingdao University of Science and Technology, Qingdao 266061, China; 2. Qingdao National Oceanographic center, Qingdao 266071, China;3. Traffic and Transportation Colleage, Chongqing Jiaotong University, Chongqing 400074, China;4. Key Laboratory of Ship Mechanics Made Traffic Engineering Industry, Dalian Maritime University, Dalian 116026, China)

    Abstract:Based on the matrix mapping theory, an improved eigensystem realization algorithm(ERA), called C/ERA, was proposed. In the algorithm, a Hankel matrix was rebuilt by replacing all the elements on each anti-subdiagonal by the arithmetic average of the elements on the anti-subdiagonal, and the concept of using the Frobenius norm (L2-norm) to control iterations number was introduced for enhancing the denoising ability of ERA. In accordance with the simulation results of a 5-DOF mass-spring-dashpot system and the experimental data of the model of a jacket-type platform under impact loading, it is proved that C/ERA has a better capacity of de-noising, higher accuracy for low order modes and can identify more high order modes than ERA does.

    Key words:modal parameter identification; C/ERA; noising

    中圖分類號:O32

    文獻(xiàn)標(biāo)志碼:A DOI:10.13465/j.cnki.jvs.2015.24.033

    通信作者盛進(jìn)路 男,博士,副教授,1978年生

    收稿日期:2014-08-07修改稿收到日期:2015-03-31

    基金項(xiàng)目:青島市建設(shè)科技計劃項(xiàng)目(JK2015-15);青島市青年專項(xiàng)(14-2-4-116-jch);大連海事大學(xué)船機(jī)修造工程交通行業(yè)重點(diǎn)實(shí)驗(yàn)室開放課題

    怎么达到女性高潮| 制服人妻中文乱码| а√天堂www在线а√下载| 999久久久精品免费观看国产| 人妻丰满熟妇av一区二区三区| 国产精品 国内视频| 日韩欧美一区二区三区在线观看| 免费在线观看视频国产中文字幕亚洲| xxxwww97欧美| 变态另类丝袜制服| 精品国产超薄肉色丝袜足j| 亚洲一区高清亚洲精品| av视频在线观看入口| 免费看a级黄色片| 国产精品一区二区三区四区免费观看 | 国产男靠女视频免费网站| 黑人操中国人逼视频| 少妇的丰满在线观看| 久久久久久大精品| 伊人久久大香线蕉亚洲五| 麻豆一二三区av精品| 日韩有码中文字幕| 免费看美女性在线毛片视频| 热99在线观看视频| 999久久久精品免费观看国产| 99久久99久久久精品蜜桃| 99久久成人亚洲精品观看| 后天国语完整版免费观看| 成在线人永久免费视频| 91麻豆精品激情在线观看国产| 日韩免费av在线播放| 亚洲自偷自拍图片 自拍| 真实男女啪啪啪动态图| 精品无人区乱码1区二区| 日韩高清综合在线| 亚洲狠狠婷婷综合久久图片| 免费电影在线观看免费观看| 97碰自拍视频| 亚洲人成网站在线播放欧美日韩| 国产高清三级在线| 久久久久国产一级毛片高清牌| 免费看十八禁软件| 国产伦人伦偷精品视频| 亚洲av免费在线观看| 啦啦啦免费观看视频1| 亚洲专区中文字幕在线| 在线观看免费视频日本深夜| 人人妻人人看人人澡| 国内精品久久久久久久电影| 欧美成人性av电影在线观看| 久久久水蜜桃国产精品网| 国产三级在线视频| 麻豆久久精品国产亚洲av| 国产精品亚洲美女久久久| 一本精品99久久精品77| 国内精品久久久久久久电影| 欧美乱妇无乱码| 国产黄色小视频在线观看| 日韩欧美三级三区| 日韩成人在线观看一区二区三区| 成人特级黄色片久久久久久久| 久久久久免费精品人妻一区二区| 午夜精品一区二区三区免费看| 校园春色视频在线观看| 国内精品久久久久久久电影| 国内精品久久久久久久电影| 国产一区二区在线av高清观看| 美女午夜性视频免费| 亚洲人成网站在线播放欧美日韩| 亚洲成av人片在线播放无| 欧美日韩亚洲国产一区二区在线观看| 在线观看日韩欧美| 日本免费a在线| 国产精品日韩av在线免费观看| 丰满人妻一区二区三区视频av | 亚洲成人精品中文字幕电影| 亚洲精品中文字幕一二三四区| 日韩欧美在线乱码| 麻豆成人午夜福利视频| 亚洲精品久久国产高清桃花| 香蕉丝袜av| 成年女人永久免费观看视频| 黄片小视频在线播放| avwww免费| 伊人久久大香线蕉亚洲五| 身体一侧抽搐| 国产亚洲精品久久久久久毛片| 又黄又爽又免费观看的视频| 精品久久久久久久久久免费视频| 日本三级黄在线观看| 老司机午夜福利在线观看视频| x7x7x7水蜜桃| 亚洲黑人精品在线| 我要搜黄色片| 亚洲狠狠婷婷综合久久图片| 99在线人妻在线中文字幕| 精品国产美女av久久久久小说| 午夜福利18| 日本五十路高清| av天堂在线播放| 99久久久亚洲精品蜜臀av| 99riav亚洲国产免费| 国产精品久久久久久亚洲av鲁大| 手机成人av网站| 国产激情偷乱视频一区二区| 成人18禁在线播放| 欧美日韩亚洲国产一区二区在线观看| 成人av在线播放网站| 在线观看午夜福利视频| 搡老妇女老女人老熟妇| 欧美一区二区国产精品久久精品| 草草在线视频免费看| 99国产极品粉嫩在线观看| 国产午夜福利久久久久久| 亚洲欧美精品综合一区二区三区| 国产熟女xx| xxx96com| 国产又黄又爽又无遮挡在线| 国产激情久久老熟女| 国产成人欧美在线观看| 国产成人一区二区三区免费视频网站| 欧美成狂野欧美在线观看| 久久久色成人| 欧美成狂野欧美在线观看| 给我免费播放毛片高清在线观看| 麻豆av在线久日| 欧美zozozo另类| 一个人免费在线观看电影 | 成人18禁在线播放| 欧美中文综合在线视频| 欧美中文综合在线视频| 亚洲av电影不卡..在线观看| 成人亚洲精品av一区二区| 亚洲av电影不卡..在线观看| 国产单亲对白刺激| 成年免费大片在线观看| e午夜精品久久久久久久| 色视频www国产| 狂野欧美白嫩少妇大欣赏| 国内精品久久久久久久电影| x7x7x7水蜜桃| 免费在线观看成人毛片| 无人区码免费观看不卡| 国产麻豆成人av免费视频| 99在线视频只有这里精品首页| 精品欧美国产一区二区三| 午夜福利视频1000在线观看| 一本久久中文字幕| 亚洲色图 男人天堂 中文字幕| 亚洲色图 男人天堂 中文字幕| 午夜激情福利司机影院| 女人被狂操c到高潮| 看免费av毛片| 国产在线精品亚洲第一网站| 亚洲av成人精品一区久久| 夜夜躁狠狠躁天天躁| 国产av麻豆久久久久久久| 91九色精品人成在线观看| 亚洲成av人片在线播放无| 韩国av一区二区三区四区| 日韩欧美在线乱码| 欧美乱色亚洲激情| 熟妇人妻久久中文字幕3abv| 免费在线观看视频国产中文字幕亚洲| 在线十欧美十亚洲十日本专区| 在线观看66精品国产| 麻豆成人av在线观看| av在线蜜桃| 亚洲美女黄片视频| 少妇丰满av| 国产伦精品一区二区三区视频9 | 国产午夜精品论理片| 国内精品一区二区在线观看| 搞女人的毛片| 亚洲国产日韩欧美精品在线观看 | 黄色片一级片一级黄色片| 日本免费a在线| 免费大片18禁| 亚洲av第一区精品v没综合| 免费搜索国产男女视频| 狠狠狠狠99中文字幕| 我要搜黄色片| 久久久久久久精品吃奶| 亚洲黑人精品在线| 亚洲av熟女| 岛国在线免费视频观看| 琪琪午夜伦伦电影理论片6080| 久久香蕉国产精品| or卡值多少钱| 日韩三级视频一区二区三区| 宅男免费午夜| 国产91精品成人一区二区三区| 国产高清激情床上av| 一a级毛片在线观看| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲av香蕉五月| 一进一出好大好爽视频| 日韩欧美国产在线观看| 国产亚洲精品一区二区www| 国产免费av片在线观看野外av| 亚洲国产色片| 国产精品亚洲一级av第二区| 18禁美女被吸乳视频| 观看免费一级毛片| 欧美av亚洲av综合av国产av| 亚洲国产色片| 99视频精品全部免费 在线 | 国语自产精品视频在线第100页| 亚洲精品美女久久av网站| 久久草成人影院| 日本黄色视频三级网站网址| 免费看美女性在线毛片视频| 国产av麻豆久久久久久久| 午夜a级毛片| av女优亚洲男人天堂 | 男人和女人高潮做爰伦理| 韩国av一区二区三区四区| 18禁国产床啪视频网站| 人妻久久中文字幕网| 国产高清视频在线播放一区| www.www免费av| 在线国产一区二区在线| 久久精品影院6| 欧美一级毛片孕妇| 亚洲真实伦在线观看| 丝袜人妻中文字幕| 激情在线观看视频在线高清| 国产精品精品国产色婷婷| 天堂影院成人在线观看| 99国产极品粉嫩在线观看| a在线观看视频网站| 亚洲一区二区三区不卡视频| 变态另类丝袜制服| 真实男女啪啪啪动态图| 欧美在线黄色| 两人在一起打扑克的视频| 亚洲中文字幕日韩| 色播亚洲综合网| 中文字幕高清在线视频| 男女床上黄色一级片免费看| 国产精品久久电影中文字幕| 午夜影院日韩av| 美女cb高潮喷水在线观看 | 丝袜人妻中文字幕| 国产亚洲精品一区二区www| 亚洲av美国av| 久久伊人香网站| 日本五十路高清| 欧美乱妇无乱码| 激情在线观看视频在线高清| 午夜亚洲福利在线播放| 少妇的逼水好多| 18禁国产床啪视频网站| 亚洲精品乱码久久久v下载方式 | 国产成人精品无人区| 国产欧美日韩精品一区二区| 国产成年人精品一区二区| 成人三级做爰电影| 欧美av亚洲av综合av国产av| 熟女人妻精品中文字幕| 亚洲精品在线美女| 俺也久久电影网| www日本在线高清视频| 日韩三级视频一区二区三区| 精品一区二区三区av网在线观看| 久久这里只有精品中国| 两人在一起打扑克的视频| 亚洲一区二区三区色噜噜| 免费搜索国产男女视频| 成人高潮视频无遮挡免费网站| 久久精品综合一区二区三区| 啦啦啦韩国在线观看视频| 久久久成人免费电影| 狂野欧美白嫩少妇大欣赏| 夜夜看夜夜爽夜夜摸| 日日干狠狠操夜夜爽| 午夜日韩欧美国产| 久久久水蜜桃国产精品网| 99久久成人亚洲精品观看| 国产成人系列免费观看| 美女高潮的动态| 国产高潮美女av| 丝袜人妻中文字幕| 三级男女做爰猛烈吃奶摸视频| 国产伦人伦偷精品视频| 中文在线观看免费www的网站| 精华霜和精华液先用哪个| 久久这里只有精品19| 国产极品精品免费视频能看的| 母亲3免费完整高清在线观看| 亚洲一区二区三区不卡视频| av在线天堂中文字幕| 国产欧美日韩精品一区二区| 老汉色av国产亚洲站长工具| 99热这里只有是精品50| 精品乱码久久久久久99久播| 国产又黄又爽又无遮挡在线| 国产精品一及| 天堂影院成人在线观看| 国产精品久久久久久精品电影| 99久久综合精品五月天人人| 亚洲国产精品sss在线观看| 国产成人一区二区三区免费视频网站| 亚洲专区中文字幕在线| 成年人黄色毛片网站| 黄色成人免费大全| 久久天躁狠狠躁夜夜2o2o| 亚洲最大成人中文| 欧美成人性av电影在线观看| 免费看光身美女| 亚洲性夜色夜夜综合| 一级毛片精品| 哪里可以看免费的av片| 精品国产乱子伦一区二区三区| 久久天堂一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 99re在线观看精品视频| 色av中文字幕| 一区福利在线观看| 校园春色视频在线观看| 男女下面进入的视频免费午夜| 美女 人体艺术 gogo| 午夜两性在线视频| 午夜福利视频1000在线观看| 在线观看66精品国产| 亚洲aⅴ乱码一区二区在线播放| 国产成+人综合+亚洲专区| 90打野战视频偷拍视频| 一区二区三区国产精品乱码| 露出奶头的视频| 免费在线观看日本一区| 久久久水蜜桃国产精品网| 99热这里只有是精品50| 波多野结衣高清作品| xxx96com| 亚洲国产欧洲综合997久久,| 看免费av毛片| 女生性感内裤真人,穿戴方法视频| 天堂网av新在线| 精品乱码久久久久久99久播| 亚洲av成人精品一区久久| 色吧在线观看| 天堂网av新在线| 丰满人妻一区二区三区视频av | 脱女人内裤的视频| 国产午夜精品久久久久久| 亚洲欧美日韩高清在线视频| 九色成人免费人妻av| 精品国产乱子伦一区二区三区| 97超视频在线观看视频| 国内精品久久久久精免费| 美女高潮喷水抽搐中文字幕| 国产精品永久免费网站| 亚洲激情在线av| 精品久久久久久久久久久久久| 在线观看一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 亚洲美女视频黄频| 国产精品永久免费网站| 国产精品,欧美在线| 精品久久久久久久久久久久久| 国产蜜桃级精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 香蕉av资源在线| 高清在线国产一区| 在线免费观看的www视频| 级片在线观看| 亚洲,欧美精品.| 九色国产91popny在线| 丰满人妻一区二区三区视频av | 色综合婷婷激情| 久久久久久大精品| 亚洲一区二区三区不卡视频| 日韩欧美免费精品| 老司机深夜福利视频在线观看| 日韩欧美免费精品| 搞女人的毛片| 99久久无色码亚洲精品果冻| 国产精品亚洲一级av第二区| 久久午夜亚洲精品久久| 午夜福利高清视频| 中文字幕精品亚洲无线码一区| 日韩欧美国产一区二区入口| 亚洲一区高清亚洲精品| 少妇丰满av| 欧美中文综合在线视频| 国产91精品成人一区二区三区| 中文字幕最新亚洲高清| 欧美黄色片欧美黄色片| 可以在线观看毛片的网站| 久久中文看片网| 制服人妻中文乱码| 国产av不卡久久| 亚洲 欧美一区二区三区| 天堂√8在线中文| 免费一级毛片在线播放高清视频| 成在线人永久免费视频| 757午夜福利合集在线观看| 亚洲国产精品成人综合色| 免费看光身美女| 欧美在线一区亚洲| 亚洲成人免费电影在线观看| 国产欧美日韩一区二区精品| 国产乱人视频| 狂野欧美白嫩少妇大欣赏| 久久久久亚洲av毛片大全| 日本撒尿小便嘘嘘汇集6| 免费看a级黄色片| 少妇的丰满在线观看| 中文字幕人成人乱码亚洲影| 国内久久婷婷六月综合欲色啪| 在线永久观看黄色视频| 亚洲av成人精品一区久久| 99在线人妻在线中文字幕| 一区二区三区国产精品乱码| 午夜亚洲福利在线播放| 亚洲精品一卡2卡三卡4卡5卡| 十八禁人妻一区二区| 亚洲一区二区三区色噜噜| 欧美日韩福利视频一区二区| 婷婷精品国产亚洲av在线| 国产单亲对白刺激| 国产蜜桃级精品一区二区三区| 亚洲天堂国产精品一区在线| 亚洲精品粉嫩美女一区| 久久精品国产综合久久久| 精品国产乱码久久久久久男人| www日本在线高清视频| 日本 欧美在线| 九九在线视频观看精品| 一级作爱视频免费观看| 波多野结衣高清无吗| 又爽又黄无遮挡网站| 最新中文字幕久久久久 | 国产野战对白在线观看| 欧美极品一区二区三区四区| 久久精品国产清高在天天线| 脱女人内裤的视频| 极品教师在线免费播放| 日韩精品中文字幕看吧| 高清毛片免费观看视频网站| 亚洲 国产 在线| 国产亚洲欧美98| 97人妻精品一区二区三区麻豆| 欧美另类亚洲清纯唯美| 男女下面进入的视频免费午夜| 高清在线国产一区| 19禁男女啪啪无遮挡网站| 国产精品 欧美亚洲| avwww免费| 天堂影院成人在线观看| 九色成人免费人妻av| 在线观看舔阴道视频| 51午夜福利影视在线观看| 精品福利观看| 精品一区二区三区视频在线 | 叶爱在线成人免费视频播放| 国产精品一区二区精品视频观看| 久久九九热精品免费| 国产三级中文精品| 精品99又大又爽又粗少妇毛片 | 91九色精品人成在线观看| 久久久久久久精品吃奶| 成人特级黄色片久久久久久久| 亚洲成人久久爱视频| 亚洲国产精品sss在线观看| 国产精品久久电影中文字幕| 日本精品一区二区三区蜜桃| 日韩欧美在线乱码| 母亲3免费完整高清在线观看| 国产99白浆流出| 久久久久性生活片| 国内揄拍国产精品人妻在线| 看免费av毛片| 搡老妇女老女人老熟妇| 亚洲 欧美 日韩 在线 免费| aaaaa片日本免费| 国产一区二区在线观看日韩 | 黑人操中国人逼视频| 国产成人精品久久二区二区91| 99视频精品全部免费 在线 | 久久这里只有精品中国| 国产激情欧美一区二区| 日本免费a在线| 在线观看免费视频日本深夜| 午夜激情欧美在线| 久久亚洲精品不卡| 中文字幕久久专区| 18禁美女被吸乳视频| 午夜两性在线视频| 亚洲欧美激情综合另类| 国产蜜桃级精品一区二区三区| 成人一区二区视频在线观看| 日韩欧美免费精品| 亚洲国产高清在线一区二区三| 欧美日韩福利视频一区二区| 美女午夜性视频免费| 亚洲av成人精品一区久久| 成在线人永久免费视频| 又紧又爽又黄一区二区| 一进一出抽搐动态| 亚洲成人中文字幕在线播放| 国产成人av激情在线播放| 亚洲九九香蕉| 日韩精品中文字幕看吧| av女优亚洲男人天堂 | 久久午夜综合久久蜜桃| 中文字幕熟女人妻在线| 国语自产精品视频在线第100页| av女优亚洲男人天堂 | 男女那种视频在线观看| 久久中文字幕一级| 日韩中文字幕欧美一区二区| 两人在一起打扑克的视频| avwww免费| 日日夜夜操网爽| 两个人视频免费观看高清| 成人特级黄色片久久久久久久| 欧美色欧美亚洲另类二区| 他把我摸到了高潮在线观看| 欧美高清成人免费视频www| a级毛片在线看网站| 999久久久国产精品视频| 熟女少妇亚洲综合色aaa.| 久久久水蜜桃国产精品网| 国产免费男女视频| 一本一本综合久久| 五月伊人婷婷丁香| 国产精品电影一区二区三区| 亚洲国产精品成人综合色| 我的老师免费观看完整版| 99久久综合精品五月天人人| АⅤ资源中文在线天堂| 国产毛片a区久久久久| 国产真实乱freesex| 97人妻精品一区二区三区麻豆| 亚洲在线观看片| 亚洲欧美一区二区三区黑人| 1000部很黄的大片| 亚洲欧美精品综合久久99| 哪里可以看免费的av片| 亚洲av五月六月丁香网| 精品熟女少妇八av免费久了| 一a级毛片在线观看| 国产伦精品一区二区三区四那| 一个人观看的视频www高清免费观看 | 波多野结衣高清无吗| www.精华液| 视频区欧美日本亚洲| av福利片在线观看| 久久这里只有精品19| 免费av毛片视频| 欧美日本视频| 黄片大片在线免费观看| 日本与韩国留学比较| 久久久久久久久免费视频了| 亚洲七黄色美女视频| 一区二区三区高清视频在线| 午夜福利成人在线免费观看| 国产精品日韩av在线免费观看| 成年版毛片免费区| 很黄的视频免费| 国产精品免费一区二区三区在线| 特级一级黄色大片| 别揉我奶头~嗯~啊~动态视频| 久久精品夜夜夜夜夜久久蜜豆| 成人精品一区二区免费| 日韩高清综合在线| 少妇人妻一区二区三区视频| 亚洲 欧美一区二区三区| 久久亚洲真实| 欧美又色又爽又黄视频| xxx96com| 最近最新中文字幕大全免费视频| 18禁观看日本| 亚洲第一电影网av| 我要搜黄色片| 一区二区三区高清视频在线| 丰满人妻熟妇乱又伦精品不卡| 听说在线观看完整版免费高清| 亚洲精品在线美女| 欧美极品一区二区三区四区| 综合色av麻豆| 久久香蕉精品热| 久久久久国产一级毛片高清牌| 国产精品九九99| 在线观看66精品国产| 欧美黄色淫秽网站| 床上黄色一级片| 一级黄色大片毛片| 又爽又黄无遮挡网站| 丰满的人妻完整版| 一本综合久久免费| 一级作爱视频免费观看| 在线永久观看黄色视频| 国产免费男女视频| 可以在线观看的亚洲视频| 久久精品国产清高在天天线| 日韩国内少妇激情av| 女警被强在线播放| 欧美午夜高清在线| 国产精品99久久久久久久久| 欧美色视频一区免费| 亚洲片人在线观看| 好看av亚洲va欧美ⅴa在| tocl精华| 久久99热这里只有精品18| 熟妇人妻久久中文字幕3abv| 麻豆成人午夜福利视频| 国产高清激情床上av| 91在线观看av| 久久久国产成人免费| 日韩 欧美 亚洲 中文字幕| 小蜜桃在线观看免费完整版高清| 国产亚洲av嫩草精品影院|