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

    隨機(jī)激勵(lì)下非線性振動(dòng)系統(tǒng)特性的定性分析

    2010-06-08 05:03:50孔憲仁楊正賢徐大富
    航天器環(huán)境工程 2010年4期
    關(guān)鍵詞:奎斯特量級(jí)頻響

    孔憲仁,廖 俊,楊正賢,徐大富

    (哈爾濱工業(yè)大學(xué) 衛(wèi)星技術(shù)研究所,哈爾濱 150086)

    0 引言

    結(jié)構(gòu)在隨機(jī)振動(dòng)試驗(yàn)時(shí),有時(shí)會(huì)超出線彈性范圍,此時(shí)線性模型已經(jīng)不能滿足需要,為正確地研究結(jié)構(gòu)特性,需要建立非線性模型。非線性動(dòng)力學(xué)系統(tǒng)的辨識(shí)通常有3個(gè)步驟:一是識(shí)別非線性的存在;二是判別非線性的性質(zhì);三是建立系統(tǒng)的非線性模型進(jìn)行參數(shù)識(shí)別。目前判別非線性的存在已經(jīng)比較成熟,主要有幅值域、頻域、時(shí)域方法[1]。

    定性研究包括系統(tǒng)非線性的定位、非線性類型、非線性方程幾種,本文集中研究非線性類型這部分。Feldmen 研究了時(shí)域中應(yīng)用Hilbert變換判別隨機(jī)激勵(lì)下系統(tǒng)非線性類型的方法[2];Stasze 應(yīng)用小波變換、Franco 應(yīng)用Gabor 變換研究了系統(tǒng)非線性定性方法[3];K.Worden[4]等研究了正弦激勵(lì)下非線性系統(tǒng)頻響函數(shù)圖及奈奎斯特圖隨外激勵(lì)量級(jí)變化而發(fā)生的變形,根據(jù)系統(tǒng)分別存在非線性阻尼和非線性剛度時(shí)變形的規(guī)律不同而由此判斷系統(tǒng)中非線性的類別,并給出了非線性定性分析的準(zhǔn)則。本文參考文獻(xiàn)[4]的準(zhǔn)則,旨在給出隨機(jī)激勵(lì)下,通過(guò)系統(tǒng)頻響函數(shù)圖及奈奎斯特圖隨外激勵(lì)量級(jí)的變化進(jìn)行定性分析的方法,并將系統(tǒng)擴(kuò)展到多自由度上。

    本文對(duì)幾種帶有非線性因素的系統(tǒng)模型進(jìn)行了近似求解,并分別畫出了它們的等效頻響函數(shù)圖及奈奎斯特圖,得到隨機(jī)激勵(lì)下系統(tǒng)頻響函數(shù)隨外激勵(lì)量級(jí)的變化規(guī)律,從而給出了適用于隨機(jī)激勵(lì)的,通過(guò)等效頻響函數(shù)隨激勵(lì)量級(jí)變化的趨勢(shì)來(lái)判別系統(tǒng)非線性特性的方法。

    1 解李雅普諾夫方程求隨機(jī)響應(yīng)方差

    對(duì)于線性時(shí)不變系統(tǒng)可以寫成方程[2]

    式中:sx為狀態(tài)向量;w為輸入向量;z為輸出向量;A、B、C為響應(yīng)維數(shù)的常量矩陣;下標(biāo)s、w、z表示矩陣維數(shù)。

    假設(shè)輸入為白噪聲,則輸入的功率譜矩陣為對(duì)稱正定的常量矩陣Wuu。狀態(tài)變量的穩(wěn)態(tài)協(xié)方差矩陣表示為

    其中Xss為協(xié)方差矩陣。

    容易證明式(2)滿足李雅普諾夫方程

    其中W為激勵(lì)w的功率譜矩陣乘以2π得到的矩陣。

    對(duì)ssA進(jìn)行特征值分解,則得

    式中:ssE為特征向量構(gòu)成的矩陣;ssΛ為特征值構(gòu)成的對(duì)角線矩陣。

    經(jīng)過(guò)矩陣推導(dǎo)求得李雅普諾夫方程的解的矩陣形式為

    式中:1矩陣為合適維數(shù)的矩陣元素全為1的矩陣;符號(hào)“./”表示矩陣之間對(duì)應(yīng)元素的相乘運(yùn)算。

    得到輸出z的均方根(RMS)值為

    這樣僅通過(guò)式(5)和式(6)矩陣運(yùn)算即可以求得所有響應(yīng)點(diǎn)均方根值(或方差)。

    對(duì)于一般線性結(jié)構(gòu)的動(dòng)力方程

    可令

    式中:Mnn、Cnn、Knn分別為系統(tǒng)的質(zhì)量、阻尼和剛度矩陣;I為n乘以n維的單位矩陣;Onn為元素為零的n乘以n維矩陣。

    將式(8)代入式(5)、(6)即可很方便地得到隨機(jī)激勵(lì)作用下響應(yīng)的方差。采用適當(dāng)Czs矩陣可以獲得響應(yīng)的位移、應(yīng)力的方差值以及von mises屈服應(yīng)力等。令Czs=[InnOnn],可以得到響應(yīng)位移的方差值。

    上述過(guò)程也可以使用模態(tài)疊加方法進(jìn)行縮減,再采用適當(dāng)?shù)腃zs將方差值從模態(tài)坐標(biāo)返回到原坐標(biāo)下。

    將式(7)通過(guò)模態(tài)變換到模態(tài)坐標(biāo)下[3],則得

    其中Mmm、Zmm、Ωmm、Φnm分別為模態(tài)質(zhì)量、模態(tài)阻尼比、固有頻率和模態(tài)的矩陣。經(jīng)過(guò)推導(dǎo)可以得到模態(tài)坐標(biāo)下的李雅普諾夫方程及其解法。

    可令

    當(dāng)Czs=[ΦnmOnm],通過(guò)式(4)、(5)、(6)即可解出各自由度位移響應(yīng)的方差值。

    2 等效線性化方法[1,5]

    從非線性隨機(jī)振動(dòng)工程分析的角度來(lái)看,如果非線性程度不是太強(qiáng),不考慮分叉、跳躍等本質(zhì)非線性現(xiàn)象,則等效線性化法已經(jīng)取得了較大的進(jìn)展,是相對(duì)有效且簡(jiǎn)單可行的方法,是工程實(shí)際中最具有應(yīng)用潛質(zhì)的預(yù)測(cè)非線性系統(tǒng)隨機(jī)響應(yīng)的近似方法[6]。其基本思想為用有精確解的線性系統(tǒng)代替給定的非線性系統(tǒng),并使兩方程之差在某種統(tǒng)計(jì)意義上最小。文獻(xiàn)[7]比較了以能量差為基礎(chǔ)以及以方程差為基礎(chǔ)的等效線性化方法計(jì)算的隨機(jī)振動(dòng)響應(yīng),認(rèn)為以能量差為等效準(zhǔn)則的方法效果要好。

    考慮具有隨機(jī)激勵(lì)的單自由度非線性系統(tǒng),其動(dòng)力學(xué)方程為

    其中: (,)gxx˙為非線性項(xiàng),m、c、k分別為系統(tǒng)質(zhì)量、阻尼、剛度。建立等效的線性化方程

    其中ce、ke為等效的阻尼系數(shù)和剛度系數(shù)。以能量差最小為準(zhǔn)則可求得ce、ke分別為

    以方程差為基礎(chǔ)的等效阻尼系數(shù)和剛度系數(shù)可通過(guò)下式獲得[7]:

    先設(shè)定ce、ke的一個(gè)初始值,求解式(11)得到x的某統(tǒng)計(jì)值,代入式(13)或(14)求得新的等效系數(shù)。經(jīng)過(guò)反復(fù)迭代后得到所求的值。一般將式(11)的等效系數(shù)的結(jié)果化成E(x?xT)的形式或其中的元素來(lái)表示,這正是式(2)得到的解,用它來(lái)參與迭代,效率和精度將得到很大的提高。只需通過(guò)式(8)或者式(10)轉(zhuǎn)換成狀態(tài)空間表示的運(yùn)動(dòng)方程,其協(xié)方差矩陣即可求得,并參與到等效線性化的迭代計(jì)算中,從而獲得穩(wěn)定的等效線性化的系數(shù)。本文采用能量差為準(zhǔn)則的等效線性化方法對(duì)系統(tǒng)進(jìn)行處理。

    3 幾種情況下的非線性特性研究

    根據(jù)線性等效原理,等效頻響函數(shù)為

    將幾種情況下的近似解分別代入式(15)即可得到存在于不同非線性特性系統(tǒng)的等效頻響函數(shù)。

    3.1 單自由度情況

    令m=1、c=5、k=100,功率譜S0分別為1、4、8、12,ε=1時(shí)為硬彈簧,ε=-1時(shí)為軟彈簧;則不同量級(jí)隨機(jī)激勵(lì)下的硬彈簧和軟彈簧的等效頻響函數(shù)圖與奈奎斯特圖如圖1~圖4所示。

    由圖1~圖4可以看出,當(dāng)存在非線性剛度時(shí),系統(tǒng)頻響函數(shù)隨著激勵(lì)量級(jí)增大而發(fā)生改變。當(dāng)K3>0時(shí),圖1中系統(tǒng)固有頻率隨激勵(lì)量級(jí)增加而增加,發(fā)生向上移頻現(xiàn)象;圖2中奈奎斯特圖向內(nèi)收縮,系統(tǒng)呈現(xiàn)典型的單自由度硬彈簧性質(zhì)[8]。由圖3、圖4可知,當(dāng)K3<0時(shí),圖3中系統(tǒng)固有頻率隨外加激勵(lì)量級(jí)增加而降低,發(fā)生向下移頻現(xiàn)象;圖4中奈奎斯特圖向外膨脹,呈現(xiàn)軟彈簧性質(zhì)。在隨機(jī)激勵(lì)試驗(yàn)中,可應(yīng)用不同量級(jí)的頻響函數(shù)曲線與外激勵(lì)的變化規(guī)律定性地分析系統(tǒng)中存在的非線性類型,從而指導(dǎo)系統(tǒng)建模。

    圖1 不同量級(jí)隨機(jī)激勵(lì)下硬彈簧等效線性化頻響函數(shù)圖(k3>0)Fig.1 Linearized FRF of a cubic stiffness model for different levels of random excitation(k3>0)

    圖2 不同量級(jí)隨機(jī)激勵(lì)下硬彈簧等效奈奎斯特圖(k3>0)Fig.2 Nyquist plot of a cubic stiffness model for different levels of random excitation(k3>0)

    圖3 不同量級(jí)隨機(jī)激勵(lì)下軟彈簧等效線性化頻響函數(shù)圖(k3<0)Fig.3 Linearized FRF of a cubic stiffness model for different levels of random excitation(k3<0)

    圖4 不同量級(jí)隨機(jī)激勵(lì)下軟彈簧等效奈奎斯特圖(k3<0)Fig.4 Nyquist plot of a cubic stiffness model for different levels of random excitation(k3<0)

    3.2 二自由度情況

    圖5為一個(gè)二自由度的系統(tǒng),兩個(gè)質(zhì)量塊之間存在非線性彈簧。

    圖5 二自由度系統(tǒng)模型Fig.5 A 2 DOF model

    3.2.1 非線性硬彈簧在兩質(zhì)量塊之間非線性硬彈簧的動(dòng)力學(xué)方程為

    式中:

    其中m=1,k=3 000,k3=2×108,c=20。激勵(lì)功率譜S1=S2,分別為1、1.5、2、2.5。該系統(tǒng)經(jīng)過(guò)變量替換可變換成一個(gè)線性單自由度系統(tǒng)和一個(gè)帶非線性硬剛度單自由度系統(tǒng)的疊加[4],這在系統(tǒng)頻響函數(shù)中得到體現(xiàn),見圖6、圖7。

    圖6 不同量級(jí)二自由度H11等效線性化頻響函數(shù)圖Fig.6 Linearized FRF of H11 for a 2 DOF system under different levels of random excitation

    圖7 不同量級(jí)二自由度頻響函數(shù)H11等效奈奎斯特圖(k3>0)Fig.7 Nyquist plot of H11for the 2 DOF system under different levels of random excitation (k3>0)

    圖6顯示,隨著激勵(lì)量級(jí)的增加,頻響函數(shù)在第一等效共振頻率處未發(fā)生頻率漂移,在第二等效共振頻率處頻率向上漂移;圖7顯示,隨著隨機(jī)激勵(lì)量級(jí)增大,奈奎斯特圖的小圓出現(xiàn)明顯向內(nèi)收縮。這說(shuō)明在第二階模態(tài)上存在硬彈簧的非線性特征。這也證明了該方法對(duì)于非線性特性分析的正確性。對(duì)于其他3個(gè)頻響函數(shù)也有類似特征,限于篇幅不予詳述。

    3.2.2 非線性硬彈簧與地面相連

    將圖 5中非線性的彈簧與連接地面和質(zhì)量塊的彈簧進(jìn)行互換,其運(yùn)動(dòng)方程形式見式(16)。此時(shí)的系統(tǒng)兩個(gè)自由度的運(yùn)動(dòng)方程都包含了非線性。限于篇幅只取了4個(gè)頻響函數(shù)中的H11和H12進(jìn)行分析計(jì)算,其他情況也可采用類似方法進(jìn)行分析。式(16)中各常量的定義分別為

    其中m=1,k=3 000,k3=2×108,c=20。激勵(lì)功率譜S1=S2,分別為 1、2、4、8。

    分析計(jì)算結(jié)果見圖8~圖11,其中:圖8、圖9分別為頻響函數(shù)H11的幅值圖和奈奎斯特圖;圖10、圖11分別為頻響函數(shù)H12的幅值圖和奈奎斯特圖。

    圖8 不同量級(jí)二自由度H11等效線性化頻響函數(shù)圖Fig.8 Linearized FRF of H11 for the 2 DOF system under different levels of random excitation

    圖 9 不同量級(jí)二自由度頻響函數(shù)H11等效奈奎斯特圖Fig.9 Nyquist plot of H11 for 2 DOFs system under different levels of random excitation

    圖10 不同量級(jí)二自由度H12等效線性化頻響函數(shù)圖Fig.10 Linearized FRF of H12 for the 2 DOF system under different levels of random excitation

    圖11 不同量級(jí)二自由度頻響函數(shù)H12等效奈奎斯特圖Fig.11 Nyquist plot of H12 for the 2 DOF system under different levels of random excitation

    圖8和圖10顯示,在兩個(gè)自由度上,隨著激勵(lì)量級(jí)的增大,都有固有頻率向上移動(dòng),這與硬彈簧的特性相同;隨著激勵(lì)量級(jí)的增大,圖9中的小圓是向外擴(kuò)張的,圖11中的兩個(gè)圓都有向外擴(kuò)張的現(xiàn)象,這是單自由度情況下軟彈簧的特點(diǎn)??梢?,多自由度因?yàn)榇嬖谀B(tài)之間的能量傳遞,使得其非線性特性有了自身的特點(diǎn),不能簡(jiǎn)單用單自由度的方法來(lái)分析。

    4 結(jié)論

    本文采用隨機(jī)激勵(lì)源,應(yīng)用求解李雅普諾夫方程與等效線性化法對(duì)非線性特性進(jìn)行了定性的研究。分別對(duì)單自由度與二自由度振動(dòng)系統(tǒng)在隨機(jī)激勵(lì)下的非線性特性進(jìn)行了分析,并給出了等效頻響函數(shù)及奈奎斯特圖隨激勵(lì)量級(jí)變化的規(guī)律,根據(jù)變化規(guī)律的不同給出隨機(jī)激勵(lì)下系統(tǒng)非線性定性研究的方法,對(duì)非線性動(dòng)力學(xué)系統(tǒng)的定性研究及建模有一定的工程應(yīng)用價(jià)值。

    (References)

    [1]朱位秋.隨機(jī)振動(dòng)[M].北京: 科學(xué)出版社, 1998

    [2]Feldoman M.Nonlinear system vibration analysis using Hilbert transform-II forced vibration analysis method‘Forcevib’[J].Mechanical Systems and Signal Processing,1994, 8(3): 309-318

    [3]Ga?tan Kerschen.Past, present and future of nonlinear system identification in structural dynamics[J].Mechanical Systems and Signal Processing, 2006, 20: 505-592

    [4]Worden K, Tomlinson G R.Nonlinearity in structural dynamics: detection, identification and modelling[M].Insitute of Physics Publishing, 2001

    [5]Elishakoff I, Zhang X.An appraisal of different stochastic linearization techniques[J].Journal of Sound and Vibration, 1992, 153(2): 370-375

    [6]朱位秋.非線性隨機(jī)振動(dòng)理論的近期進(jìn)展[J].力學(xué)進(jìn)展, 1994, 24(2): 163-173

    [7]羅務(wù)揆, 邢譽(yù)峰, 楊蓉.辛算法在隨機(jī)振動(dòng)響應(yīng)求解中的應(yīng)用[J].航空動(dòng)力學(xué)報(bào), 2008, 23(1): 37-43

    [8]何蕊, 羅文波, 王本利, 等.基于Lyapunov指數(shù)的鋁蜂窩板頻率漂移機(jī)理分析[J].宇航學(xué)報(bào), 2009, 30(2):463-467

    猜你喜歡
    奎斯特量級(jí)頻響
    含例鄰域邏輯的薩奎斯特對(duì)應(yīng)理論
    基于分塊化頻響函數(shù)曲率比的砌體房屋模型損傷識(shí)別研究
    美團(tuán)外賣哥
    論《瑪莎·奎斯特》中瑪莎的女性主義意識(shí)
    21連勝
    頻響函數(shù)殘差法在有限元模型修正中的應(yīng)用
    頻響阻抗法診斷變壓器繞組變形
    多麗絲·萊辛《瑪莎·奎斯特》中的邊緣人
    婷婷精品国产亚洲av| 国内精品一区二区在线观看| 黑人巨大精品欧美一区二区mp4| 国产精品亚洲一级av第二区| 欧美绝顶高潮抽搐喷水| 国产精品精品国产色婷婷| 亚洲精品粉嫩美女一区| 午夜福利欧美成人| 成人国产一区最新在线观看| 少妇人妻一区二区三区视频| 国产一区在线观看成人免费| 99国产极品粉嫩在线观看| 三级男女做爰猛烈吃奶摸视频| 成年版毛片免费区| 国产精品99久久99久久久不卡| 中文亚洲av片在线观看爽| 精品久久久久久,| 搡老岳熟女国产| 亚洲中文av在线| 一级片免费观看大全| 99久久国产精品久久久| 亚洲中文日韩欧美视频| 日本五十路高清| 亚洲精品美女久久av网站| 中国美女看黄片| 精品久久久久久成人av| www日本在线高清视频| 日韩精品青青久久久久久| www.999成人在线观看| 一级片免费观看大全| 欧美一级毛片孕妇| 欧美日韩精品网址| 18禁国产床啪视频网站| 禁无遮挡网站| 丁香六月欧美| 日本五十路高清| 亚洲无线在线观看| 热99re8久久精品国产| 成人三级做爰电影| 女警被强在线播放| 色老头精品视频在线观看| 色综合婷婷激情| 又黄又爽又免费观看的视频| ponron亚洲| 精品少妇一区二区三区视频日本电影| 18禁观看日本| 丰满人妻一区二区三区视频av | 男人的好看免费观看在线视频 | 国产精品98久久久久久宅男小说| 999精品在线视频| 在线观看一区二区三区| 露出奶头的视频| 日韩av在线大香蕉| 一本一本综合久久| 欧美中文日本在线观看视频| 久久午夜综合久久蜜桃| 国产熟女xx| 一区二区三区激情视频| 村上凉子中文字幕在线| 亚洲激情在线av| 成在线人永久免费视频| 十八禁人妻一区二区| 少妇的丰满在线观看| 久久人妻av系列| 成人永久免费在线观看视频| 在线观看www视频免费| 欧美av亚洲av综合av国产av| 亚洲精品色激情综合| 51午夜福利影视在线观看| 制服诱惑二区| 99久久综合精品五月天人人| 窝窝影院91人妻| 日本成人三级电影网站| 亚洲美女黄片视频| 男女床上黄色一级片免费看| 久久精品国产亚洲av高清一级| 国产真实乱freesex| 欧美中文综合在线视频| 久热爱精品视频在线9| 悠悠久久av| 国产真人三级小视频在线观看| 亚洲中文av在线| 色播亚洲综合网| 欧美色欧美亚洲另类二区| 日本一本二区三区精品| 欧美中文日本在线观看视频| 亚洲成av人片在线播放无| 亚洲精品在线美女| 色综合亚洲欧美另类图片| 久久精品国产亚洲av高清一级| 国产成人aa在线观看| 精品欧美一区二区三区在线| 国产v大片淫在线免费观看| 欧美性猛交黑人性爽| 老汉色av国产亚洲站长工具| 九九热线精品视视频播放| 91成年电影在线观看| 日本 欧美在线| 亚洲人成伊人成综合网2020| 欧美中文综合在线视频| 亚洲av中文字字幕乱码综合| www日本在线高清视频| 999久久久国产精品视频| 久久亚洲精品不卡| 舔av片在线| 日本一二三区视频观看| 精品午夜福利视频在线观看一区| 中文字幕高清在线视频| 精品不卡国产一区二区三区| 一a级毛片在线观看| 88av欧美| 久久久国产成人免费| 免费一级毛片在线播放高清视频| 日本免费一区二区三区高清不卡| 午夜亚洲福利在线播放| 99热这里只有精品一区 | 精品福利观看| 五月玫瑰六月丁香| 99热6这里只有精品| 香蕉久久夜色| 久久精品国产清高在天天线| 免费在线观看黄色视频的| 亚洲中文日韩欧美视频| 婷婷六月久久综合丁香| 成人手机av| 在线观看免费日韩欧美大片| 麻豆国产av国片精品| 看免费av毛片| 国产成人系列免费观看| www.www免费av| 精品久久久久久久久久免费视频| 国产三级中文精品| 免费搜索国产男女视频| 日日夜夜操网爽| 亚洲 国产 在线| 久久久精品欧美日韩精品| 99久久综合精品五月天人人| 亚洲电影在线观看av| 九色成人免费人妻av| 精品一区二区三区四区五区乱码| 小说图片视频综合网站| 琪琪午夜伦伦电影理论片6080| 最新在线观看一区二区三区| 国产亚洲精品久久久久久毛片| 亚洲国产欧洲综合997久久,| 欧美中文日本在线观看视频| 在线a可以看的网站| 免费av毛片视频| 最近最新中文字幕大全免费视频| 午夜福利免费观看在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成av人片在线播放无| 此物有八面人人有两片| 国产黄a三级三级三级人| 久9热在线精品视频| 国产成人av激情在线播放| 精品福利观看| 亚洲国产精品999在线| 国产一区二区三区在线臀色熟女| 热99re8久久精品国产| 国产成人精品久久二区二区免费| 国产三级中文精品| 欧美日韩中文字幕国产精品一区二区三区| 禁无遮挡网站| 国产99久久九九免费精品| 麻豆久久精品国产亚洲av| 免费高清视频大片| 毛片女人毛片| 淫妇啪啪啪对白视频| 成年版毛片免费区| 精品人妻1区二区| 亚洲av日韩精品久久久久久密| 制服诱惑二区| 观看免费一级毛片| 在线观看66精品国产| 欧美日本亚洲视频在线播放| 久久这里只有精品中国| 亚洲狠狠婷婷综合久久图片| 国产欧美日韩一区二区精品| 级片在线观看| 欧美日韩精品网址| 成人三级黄色视频| 丝袜美腿诱惑在线| 精品无人区乱码1区二区| 国产激情欧美一区二区| 动漫黄色视频在线观看| 少妇被粗大的猛进出69影院| 日韩av在线大香蕉| 国产精品亚洲av一区麻豆| 国产亚洲精品av在线| 亚洲av成人不卡在线观看播放网| 国产一区二区在线av高清观看| 欧美色视频一区免费| 久久精品国产99精品国产亚洲性色| 国产精品一区二区精品视频观看| 久久久久久大精品| 国产av一区在线观看免费| 伊人久久大香线蕉亚洲五| 久久精品国产亚洲av高清一级| 99热这里只有是精品50| 久久人人精品亚洲av| 男女视频在线观看网站免费 | 男女下面进入的视频免费午夜| 好男人在线观看高清免费视频| 欧美绝顶高潮抽搐喷水| 国产成年人精品一区二区| 成人国语在线视频| ponron亚洲| 99久久无色码亚洲精品果冻| 中文字幕av在线有码专区| 两个人免费观看高清视频| 欧美日韩福利视频一区二区| 亚洲中文字幕日韩| 69av精品久久久久久| 国产激情欧美一区二区| 亚洲激情在线av| 国产亚洲欧美98| 在线a可以看的网站| 黄色视频,在线免费观看| 国产又黄又爽又无遮挡在线| av欧美777| 欧美日韩一级在线毛片| 国产区一区二久久| 又黄又粗又硬又大视频| 麻豆国产97在线/欧美 | av中文乱码字幕在线| 99国产精品一区二区蜜桃av| 在线国产一区二区在线| 国产精品久久久久久精品电影| 欧美色欧美亚洲另类二区| 亚洲中文字幕日韩| 午夜亚洲福利在线播放| 一区二区三区国产精品乱码| 岛国视频午夜一区免费看| 久久草成人影院| 99久久久亚洲精品蜜臀av| 亚洲中文av在线| 99国产综合亚洲精品| 一进一出抽搐gif免费好疼| 亚洲乱码一区二区免费版| 亚洲狠狠婷婷综合久久图片| 午夜福利成人在线免费观看| 国产精品免费视频内射| 欧美黑人欧美精品刺激| 亚洲欧美日韩高清专用| 亚洲av片天天在线观看| 亚洲国产精品sss在线观看| 人成视频在线观看免费观看| 亚洲第一电影网av| 欧美三级亚洲精品| 19禁男女啪啪无遮挡网站| 免费在线观看亚洲国产| 国产伦在线观看视频一区| 国产一区在线观看成人免费| 日本 欧美在线| 亚洲av成人精品一区久久| 18禁观看日本| 成在线人永久免费视频| 三级毛片av免费| 无限看片的www在线观看| 在线永久观看黄色视频| 国产精华一区二区三区| 日韩欧美在线乱码| 日韩精品青青久久久久久| 国产三级黄色录像| 国产av又大| 久久国产精品人妻蜜桃| 白带黄色成豆腐渣| 色噜噜av男人的天堂激情| 久久香蕉激情| 日本黄大片高清| 日本 欧美在线| 日日夜夜操网爽| 国产精品一区二区精品视频观看| 悠悠久久av| 麻豆国产97在线/欧美 | 欧美色欧美亚洲另类二区| 欧美大码av| 国产精品乱码一区二三区的特点| 18禁观看日本| 亚洲av成人av| 欧美成狂野欧美在线观看| 最近最新中文字幕大全电影3| 国语自产精品视频在线第100页| АⅤ资源中文在线天堂| e午夜精品久久久久久久| 久久 成人 亚洲| 久久国产乱子伦精品免费另类| 午夜亚洲福利在线播放| 免费在线观看亚洲国产| 中文字幕久久专区| 少妇粗大呻吟视频| 亚洲av日韩精品久久久久久密| 久久久精品国产亚洲av高清涩受| 亚洲av成人一区二区三| 午夜免费成人在线视频| 亚洲国产精品成人综合色| 欧美中文综合在线视频| 999久久久精品免费观看国产| 床上黄色一级片| 欧美黄色淫秽网站| 成人高潮视频无遮挡免费网站| 久久精品国产亚洲av香蕉五月| tocl精华| 免费观看精品视频网站| 亚洲乱码一区二区免费版| 91国产中文字幕| 精品福利观看| 亚洲精华国产精华精| 久久久久免费精品人妻一区二区| 国产一区二区三区在线臀色熟女| 波多野结衣巨乳人妻| 深夜精品福利| 在线播放国产精品三级| 黄色女人牲交| 久久伊人香网站| 精华霜和精华液先用哪个| bbb黄色大片| 不卡一级毛片| 久久性视频一级片| 国产99白浆流出| a在线观看视频网站| 久久久久国产精品人妻aⅴ院| 精品国产乱码久久久久久男人| 可以在线观看毛片的网站| 深夜精品福利| 亚洲在线自拍视频| 99精品欧美一区二区三区四区| 日本三级黄在线观看| 美女大奶头视频| aaaaa片日本免费| 伊人久久大香线蕉亚洲五| 亚洲成av人片免费观看| 欧美三级亚洲精品| 真人一进一出gif抽搐免费| 久久国产精品影院| 99riav亚洲国产免费| 成人国产综合亚洲| 69av精品久久久久久| 无限看片的www在线观看| av在线天堂中文字幕| 亚洲av中文字字幕乱码综合| 露出奶头的视频| 国产激情欧美一区二区| 欧美日韩中文字幕国产精品一区二区三区| 国产激情偷乱视频一区二区| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 亚洲av成人av| 99精品在免费线老司机午夜| 欧美一区二区国产精品久久精品 | 国产精品久久久av美女十八| 性欧美人与动物交配| 女人高潮潮喷娇喘18禁视频| 夜夜夜夜夜久久久久| 久久久久久九九精品二区国产 | 色在线成人网| 美女 人体艺术 gogo| 三级国产精品欧美在线观看 | 久久久精品大字幕| 亚洲精品色激情综合| 亚洲人成网站高清观看| 18禁国产床啪视频网站| 夜夜躁狠狠躁天天躁| 国产精品久久久久久精品电影| 国产精品永久免费网站| 香蕉av资源在线| 日本在线视频免费播放| 国产精品av视频在线免费观看| 美女扒开内裤让男人捅视频| 日韩精品免费视频一区二区三区| 国产激情久久老熟女| avwww免费| 校园春色视频在线观看| 亚洲美女视频黄频| 中文字幕最新亚洲高清| 五月玫瑰六月丁香| 亚洲 欧美一区二区三区| 国产欧美日韩精品亚洲av| 免费搜索国产男女视频| 在线观看舔阴道视频| 国产精品自产拍在线观看55亚洲| 搞女人的毛片| 久久久久久久精品吃奶| 在线观看66精品国产| 日韩高清综合在线| 亚洲国产高清在线一区二区三| 国产亚洲av嫩草精品影院| 中文亚洲av片在线观看爽| 日本一二三区视频观看| 我的老师免费观看完整版| 99在线视频只有这里精品首页| 欧美三级亚洲精品| 午夜精品久久久久久毛片777| 黄色 视频免费看| 一本一本综合久久| 国产三级中文精品| 久久久久久亚洲精品国产蜜桃av| 国产又色又爽无遮挡免费看| 国产成人啪精品午夜网站| 色综合亚洲欧美另类图片| 欧美中文综合在线视频| 欧美黄色片欧美黄色片| 看片在线看免费视频| 床上黄色一级片| 看黄色毛片网站| 亚洲性夜色夜夜综合| 欧美乱码精品一区二区三区| 国产亚洲精品久久久久5区| 亚洲自偷自拍图片 自拍| 大型av网站在线播放| 成人国产综合亚洲| 午夜免费激情av| 亚洲 欧美 日韩 在线 免费| 免费在线观看黄色视频的| 久久精品人妻少妇| www国产在线视频色| 午夜福利欧美成人| 天堂影院成人在线观看| 长腿黑丝高跟| 麻豆一二三区av精品| a级毛片a级免费在线| 成人欧美大片| 亚洲人成77777在线视频| 久久久久免费精品人妻一区二区| 人成视频在线观看免费观看| 亚洲国产高清在线一区二区三| 青草久久国产| 国产精品久久久久久亚洲av鲁大| 成年女人毛片免费观看观看9| 国产一区二区在线观看日韩 | 最近最新免费中文字幕在线| 精品国产亚洲在线| av福利片在线| 2021天堂中文幕一二区在线观| 久久久久国产一级毛片高清牌| 国产又色又爽无遮挡免费看| 男女那种视频在线观看| 久久久国产精品麻豆| 国产亚洲精品综合一区在线观看 | 日韩中文字幕欧美一区二区| 久久午夜亚洲精品久久| 日韩精品青青久久久久久| 少妇的丰满在线观看| 男女做爰动态图高潮gif福利片| 午夜久久久久精精品| 99在线人妻在线中文字幕| 中文字幕高清在线视频| 亚洲一码二码三码区别大吗| 在线观看www视频免费| 午夜精品久久久久久毛片777| 亚洲精品在线观看二区| 特大巨黑吊av在线直播| 国产黄a三级三级三级人| 国产激情欧美一区二区| 午夜亚洲福利在线播放| 国产精品98久久久久久宅男小说| 国产熟女午夜一区二区三区| 久久久国产欧美日韩av| svipshipincom国产片| 亚洲精品色激情综合| 日日爽夜夜爽网站| 午夜亚洲福利在线播放| 久久香蕉精品热| 久久香蕉激情| 丁香欧美五月| 制服人妻中文乱码| 一个人免费在线观看的高清视频| 欧美不卡视频在线免费观看 | 天天躁夜夜躁狠狠躁躁| 久久久久国产一级毛片高清牌| 久久精品国产亚洲av香蕉五月| 精品国产乱子伦一区二区三区| 久久久国产精品麻豆| www日本在线高清视频| 俺也久久电影网| 欧美一区二区精品小视频在线| 欧美黄色淫秽网站| 色综合亚洲欧美另类图片| 国产99久久九九免费精品| 叶爱在线成人免费视频播放| 国产精品国产高清国产av| 99riav亚洲国产免费| 久久久久久人人人人人| 性欧美人与动物交配| 欧美午夜高清在线| 成人欧美大片| 欧美zozozo另类| 国产99白浆流出| 亚洲人成伊人成综合网2020| 久久这里只有精品中国| 国产探花在线观看一区二区| 51午夜福利影视在线观看| 亚洲成人免费电影在线观看| 午夜激情福利司机影院| 久久亚洲精品不卡| 国产1区2区3区精品| 搡老妇女老女人老熟妇| 亚洲 欧美 日韩 在线 免费| 非洲黑人性xxxx精品又粗又长| avwww免费| 啦啦啦韩国在线观看视频| av在线天堂中文字幕| 一级黄色大片毛片| 亚洲人成伊人成综合网2020| 久久精品夜夜夜夜夜久久蜜豆 | 久久婷婷人人爽人人干人人爱| 高潮久久久久久久久久久不卡| 床上黄色一级片| 少妇熟女aⅴ在线视频| 嫁个100分男人电影在线观看| 伊人久久大香线蕉亚洲五| 欧美日韩亚洲国产一区二区在线观看| 亚洲18禁久久av| 成人永久免费在线观看视频| 久久性视频一级片| 欧美日本亚洲视频在线播放| 久久国产精品人妻蜜桃| 90打野战视频偷拍视频| 久热爱精品视频在线9| 热99re8久久精品国产| 搡老岳熟女国产| 欧美一级a爱片免费观看看 | 精品国产超薄肉色丝袜足j| 亚洲人成伊人成综合网2020| 国产欧美日韩一区二区三| 国产av麻豆久久久久久久| 欧美zozozo另类| 99久久久亚洲精品蜜臀av| 又爽又黄无遮挡网站| 免费在线观看完整版高清| 亚洲成av人片在线播放无| 国模一区二区三区四区视频 | videosex国产| 国产黄片美女视频| 久久久精品国产亚洲av高清涩受| 国产免费av片在线观看野外av| 免费一级毛片在线播放高清视频| 久久亚洲真实| 18禁观看日本| 国产成人啪精品午夜网站| 欧美日韩中文字幕国产精品一区二区三区| 国产精品1区2区在线观看.| 操出白浆在线播放| 国产精品av视频在线免费观看| 啦啦啦观看免费观看视频高清| 精品免费久久久久久久清纯| 每晚都被弄得嗷嗷叫到高潮| 老司机福利观看| 国产一区二区在线av高清观看| 成年版毛片免费区| 国产精品美女特级片免费视频播放器 | 午夜日韩欧美国产| 麻豆av在线久日| 亚洲av成人不卡在线观看播放网| 国产精品永久免费网站| 黄色视频,在线免费观看| 亚洲国产高清在线一区二区三| 夜夜夜夜夜久久久久| 亚洲精品色激情综合| 亚洲欧美日韩东京热| 好男人在线观看高清免费视频| 国产一区二区在线观看日韩 | 操出白浆在线播放| 一区福利在线观看| 国产亚洲精品综合一区在线观看 | 人人妻,人人澡人人爽秒播| 伦理电影免费视频| 亚洲av电影在线进入| 国产视频一区二区在线看| xxxwww97欧美| 男女床上黄色一级片免费看| 欧美性长视频在线观看| 日韩有码中文字幕| 丰满的人妻完整版| 搡老熟女国产l中国老女人| 两性夫妻黄色片| 国产不卡一卡二| 免费一级毛片在线播放高清视频| 国模一区二区三区四区视频 | 男女之事视频高清在线观看| 身体一侧抽搐| 熟女少妇亚洲综合色aaa.| 国产三级黄色录像| 一边摸一边做爽爽视频免费| aaaaa片日本免费| 中文字幕人妻丝袜一区二区| 亚洲午夜精品一区,二区,三区| 国产精品久久久人人做人人爽| 一区二区三区国产精品乱码| a在线观看视频网站| 搡老岳熟女国产| 99国产精品99久久久久| 男人的好看免费观看在线视频 | 淫妇啪啪啪对白视频| 欧美在线一区亚洲| 此物有八面人人有两片| 亚洲av成人精品一区久久| 亚洲国产精品成人综合色| 男人舔奶头视频| 99热这里只有是精品50| 免费人成视频x8x8入口观看| 日本五十路高清| 91老司机精品| 亚洲黑人精品在线| 国产精品影院久久| 国产在线精品亚洲第一网站| 色播亚洲综合网| 最好的美女福利视频网| 国产精品一区二区免费欧美| 亚洲黑人精品在线| 国产亚洲精品第一综合不卡| 国内精品久久久久精免费| 免费人成视频x8x8入口观看| 亚洲精品一区av在线观看|