• <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)用
    頻響阻抗法診斷變壓器繞組變形
    多麗絲·萊辛《瑪莎·奎斯特》中的邊緣人
    久久久精品免费免费高清| 18禁裸乳无遮挡动漫免费视频| 国产单亲对白刺激| 可以免费在线观看a视频的电影网站| 国产精品美女特级片免费视频播放器 | 丝袜人妻中文字幕| 人人妻人人澡人人爽人人夜夜| 亚洲av美国av| 国产欧美日韩一区二区三| 99国产极品粉嫩在线观看| 12—13女人毛片做爰片一| 一级a爱视频在线免费观看| 中文字幕精品免费在线观看视频| 在线观看免费高清a一片| 国产aⅴ精品一区二区三区波| 国产精品免费大片| 国产亚洲午夜精品一区二区久久| 成人国语在线视频| 欧美日韩福利视频一区二区| 久久精品亚洲av国产电影网| 桃红色精品国产亚洲av| 1024视频免费在线观看| av网站免费在线观看视频| 亚洲精品中文字幕一二三四区 | 91九色精品人成在线观看| 老司机深夜福利视频在线观看| 777久久人妻少妇嫩草av网站| 下体分泌物呈黄色| 悠悠久久av| 丝袜喷水一区| 国产真人三级小视频在线观看| 80岁老熟妇乱子伦牲交| av视频免费观看在线观看| tube8黄色片| 日韩欧美一区视频在线观看| 久久久久精品人妻al黑| 午夜福利免费观看在线| 一级片免费观看大全| 亚洲精品自拍成人| 免费人妻精品一区二区三区视频| 亚洲色图 男人天堂 中文字幕| 久久久精品国产亚洲av高清涩受| 黄色丝袜av网址大全| 色综合欧美亚洲国产小说| 777米奇影视久久| 久久久久精品国产欧美久久久| 精品一区二区三区av网在线观看 | 国产无遮挡羞羞视频在线观看| 韩国精品一区二区三区| 久久中文看片网| 久久国产精品大桥未久av| 亚洲av第一区精品v没综合| 最近最新中文字幕大全免费视频| 夜夜爽天天搞| 热99re8久久精品国产| 性少妇av在线| 成年女人毛片免费观看观看9 | 国产精品亚洲av一区麻豆| 天堂8中文在线网| 色婷婷久久久亚洲欧美| 麻豆乱淫一区二区| 国产日韩欧美在线精品| 日本精品一区二区三区蜜桃| 久久久久久久久久久久大奶| 欧美日韩福利视频一区二区| 91麻豆av在线| 久久热在线av| 热re99久久国产66热| 女性被躁到高潮视频| 女性被躁到高潮视频| 欧美国产精品va在线观看不卡| 又紧又爽又黄一区二区| 新久久久久国产一级毛片| 一区福利在线观看| a级片在线免费高清观看视频| 美女高潮喷水抽搐中文字幕| 91精品三级在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 在线永久观看黄色视频| 日本精品一区二区三区蜜桃| 日韩三级视频一区二区三区| 在线天堂中文资源库| 欧美激情极品国产一区二区三区| 精品国产一区二区三区久久久樱花| 黑丝袜美女国产一区| 狠狠狠狠99中文字幕| a级毛片在线看网站| 日韩精品免费视频一区二区三区| 国产在线一区二区三区精| 搡老岳熟女国产| 午夜福利,免费看| 91成年电影在线观看| 久久中文字幕人妻熟女| 啪啪无遮挡十八禁网站| 丝袜人妻中文字幕| 欧美日韩一级在线毛片| 久久久久国产一级毛片高清牌| 亚洲色图综合在线观看| 色综合婷婷激情| 免费在线观看视频国产中文字幕亚洲| 老熟妇乱子伦视频在线观看| 亚洲三区欧美一区| 男人舔女人的私密视频| 新久久久久国产一级毛片| 午夜福利,免费看| 国产免费福利视频在线观看| 日韩有码中文字幕| 欧美黄色淫秽网站| 国产亚洲欧美精品永久| 国产激情久久老熟女| 伦理电影免费视频| 成人18禁高潮啪啪吃奶动态图| 啦啦啦视频在线资源免费观看| 黑人欧美特级aaaaaa片| 欧美人与性动交α欧美软件| 两性午夜刺激爽爽歪歪视频在线观看 | 国产黄频视频在线观看| tocl精华| 婷婷成人精品国产| 在线观看人妻少妇| 国产在线视频一区二区| 男女下面插进去视频免费观看| 日本撒尿小便嘘嘘汇集6| 脱女人内裤的视频| 亚洲专区国产一区二区| 12—13女人毛片做爰片一| 日韩成人在线观看一区二区三区| 99re在线观看精品视频| 宅男免费午夜| 久久中文字幕一级| 久久午夜亚洲精品久久| 大码成人一级视频| 黄色a级毛片大全视频| 亚洲午夜理论影院| 老司机亚洲免费影院| 777久久人妻少妇嫩草av网站| 建设人人有责人人尽责人人享有的| 丁香六月天网| 一本一本久久a久久精品综合妖精| kizo精华| av有码第一页| 最近最新中文字幕大全免费视频| 国产精品成人在线| tocl精华| 欧美一级毛片孕妇| 精品少妇久久久久久888优播| 午夜日韩欧美国产| 水蜜桃什么品种好| 老司机靠b影院| 免费看a级黄色片| 国产av一区二区精品久久| 中文字幕人妻丝袜一区二区| 性高湖久久久久久久久免费观看| 性高湖久久久久久久久免费观看| 国产伦理片在线播放av一区| av线在线观看网站| 久久精品国产a三级三级三级| av在线播放免费不卡| 在线观看一区二区三区激情| 日韩欧美一区二区三区在线观看 | 久久热在线av| 欧美乱码精品一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 悠悠久久av| 精品视频人人做人人爽| 国产亚洲欧美在线一区二区| 女性生殖器流出的白浆| av网站免费在线观看视频| 中文欧美无线码| 成人黄色视频免费在线看| 久9热在线精品视频| 日本撒尿小便嘘嘘汇集6| 婷婷丁香在线五月| 久久青草综合色| 变态另类成人亚洲欧美熟女 | 变态另类成人亚洲欧美熟女 | 久热这里只有精品99| 一夜夜www| 久久精品亚洲av国产电影网| 久久香蕉激情| 亚洲精品乱久久久久久| 99国产精品99久久久久| 欧美日韩亚洲高清精品| 国产精品香港三级国产av潘金莲| 十八禁网站网址无遮挡| 久久精品国产亚洲av高清一级| 免费人妻精品一区二区三区视频| av又黄又爽大尺度在线免费看| 久久 成人 亚洲| 精品国产超薄肉色丝袜足j| 欧美日韩黄片免| 80岁老熟妇乱子伦牲交| 人人妻人人澡人人爽人人夜夜| 夜夜爽天天搞| 波多野结衣一区麻豆| 无人区码免费观看不卡 | 麻豆乱淫一区二区| 亚洲少妇的诱惑av| 久久亚洲精品不卡| 亚洲av电影在线进入| 欧美精品一区二区免费开放| 亚洲国产中文字幕在线视频| 91九色精品人成在线观看| 国产高清视频在线播放一区| 亚洲成人国产一区在线观看| 亚洲黑人精品在线| 国产成人免费观看mmmm| 国精品久久久久久国模美| 超碰97精品在线观看| 久久国产精品男人的天堂亚洲| 美女高潮到喷水免费观看| 久久天躁狠狠躁夜夜2o2o| 男人操女人黄网站| 人人澡人人妻人| 亚洲成人免费电影在线观看| 99国产极品粉嫩在线观看| 黄片小视频在线播放| 最近最新中文字幕大全免费视频| 国产精品秋霞免费鲁丝片| 99精品在免费线老司机午夜| 日韩一区二区三区影片| 丁香六月天网| 国产午夜精品久久久久久| 中亚洲国语对白在线视频| 精品一区二区三区四区五区乱码| 欧美精品av麻豆av| 亚洲精品一卡2卡三卡4卡5卡| 十八禁高潮呻吟视频| 国产亚洲欧美在线一区二区| 夜夜爽天天搞| av一本久久久久| 国产精品国产高清国产av | 18禁观看日本| 亚洲精品国产精品久久久不卡| 丝袜人妻中文字幕| 午夜激情久久久久久久| 在线观看免费高清a一片| 免费日韩欧美在线观看| 丝瓜视频免费看黄片| 国产一区二区三区综合在线观看| 最新的欧美精品一区二区| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久精品电影小说| 新久久久久国产一级毛片| 免费高清在线观看日韩| 熟女少妇亚洲综合色aaa.| 亚洲av成人一区二区三| 动漫黄色视频在线观看| 日本欧美视频一区| 日本精品一区二区三区蜜桃| 1024视频免费在线观看| 国产人伦9x9x在线观看| 成人特级黄色片久久久久久久 | 极品少妇高潮喷水抽搐| 亚洲,欧美精品.| 大陆偷拍与自拍| 黄色片一级片一级黄色片| 国产成人av激情在线播放| 少妇 在线观看| 国产在线观看jvid| 丝袜美腿诱惑在线| 99热网站在线观看| 午夜福利影视在线免费观看| 老熟女久久久| 亚洲av电影在线进入| 久久香蕉激情| 免费久久久久久久精品成人欧美视频| 啦啦啦中文免费视频观看日本| 久久久久久久大尺度免费视频| 国产91精品成人一区二区三区 | 国产成人精品久久二区二区免费| 亚洲五月色婷婷综合| 美女午夜性视频免费| 成在线人永久免费视频| 精品久久久久久电影网| 美女福利国产在线| 欧美 日韩 精品 国产| 99国产精品免费福利视频| 一级,二级,三级黄色视频| 成年女人毛片免费观看观看9 | 国产精品美女特级片免费视频播放器 | 999精品在线视频| 80岁老熟妇乱子伦牲交| cao死你这个sao货| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产综合久久久| 伊人久久大香线蕉亚洲五| 女性生殖器流出的白浆| 中文字幕人妻熟女乱码| 久久午夜综合久久蜜桃| 又紧又爽又黄一区二区| 久久毛片免费看一区二区三区| 欧美成狂野欧美在线观看| 人人妻人人添人人爽欧美一区卜| 18禁观看日本| 久久久久久久久久久久大奶| 中文字幕av电影在线播放| 精品人妻在线不人妻| 日韩有码中文字幕| 中文字幕精品免费在线观看视频| 午夜福利乱码中文字幕| 青草久久国产| 一本综合久久免费| 精品人妻在线不人妻| 99香蕉大伊视频| 久久99热这里只频精品6学生| 超碰97精品在线观看| 亚洲av日韩精品久久久久久密| 国产黄色免费在线视频| 欧美精品av麻豆av| 亚洲五月色婷婷综合| 午夜福利一区二区在线看| 人成视频在线观看免费观看| 国产精品秋霞免费鲁丝片| 国产免费视频播放在线视频| 男人舔女人的私密视频| 国产伦理片在线播放av一区| 最近最新中文字幕大全免费视频| 久久久久视频综合| 亚洲av成人一区二区三| 国产高清视频在线播放一区| 亚洲成人国产一区在线观看| 天天操日日干夜夜撸| 国产片内射在线| 麻豆乱淫一区二区| 国产欧美日韩一区二区三| 最近最新中文字幕大全免费视频| 久久久国产一区二区| 9热在线视频观看99| av视频免费观看在线观看| 久久性视频一级片| 窝窝影院91人妻| 日本五十路高清| 满18在线观看网站| 99久久精品国产亚洲精品| 亚洲成人手机| 丁香六月天网| 国产福利在线免费观看视频| 亚洲欧洲日产国产| 日韩一区二区三区影片| 国产亚洲欧美在线一区二区| 热99re8久久精品国产| 91精品国产国语对白视频| av天堂在线播放| 久久人妻熟女aⅴ| 亚洲精品乱久久久久久| 一区二区三区乱码不卡18| 丁香欧美五月| av网站免费在线观看视频| 中文字幕最新亚洲高清| 中文亚洲av片在线观看爽 | 欧美精品亚洲一区二区| 国产男靠女视频免费网站| 免费观看av网站的网址| 国产单亲对白刺激| 美国免费a级毛片| 天天影视国产精品| 51午夜福利影视在线观看| av又黄又爽大尺度在线免费看| 欧美午夜高清在线| 久久亚洲真实| 国产精品一区二区在线不卡| 亚洲欧美激情在线| av有码第一页| 99热网站在线观看| 黄色成人免费大全| 亚洲成av片中文字幕在线观看| 久久午夜亚洲精品久久| 国产精品 欧美亚洲| 女性被躁到高潮视频| 欧美黑人精品巨大| 黄片大片在线免费观看| 国产99久久九九免费精品| 不卡一级毛片| 一本一本久久a久久精品综合妖精| 国产精品成人在线| 黄片播放在线免费| av天堂在线播放| 亚洲七黄色美女视频| 色综合欧美亚洲国产小说| 色婷婷久久久亚洲欧美| 精品一区二区三卡| 欧美精品人与动牲交sv欧美| 搡老乐熟女国产| 伊人久久大香线蕉亚洲五| 女人久久www免费人成看片| 成人国产av品久久久| 免费黄频网站在线观看国产| 69av精品久久久久久 | 成年动漫av网址| 高清毛片免费观看视频网站 | 国产在线一区二区三区精| 人成视频在线观看免费观看| 淫妇啪啪啪对白视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品 国内视频| 国产精品电影一区二区三区 | 国产欧美日韩一区二区三区在线| 99国产综合亚洲精品| 成人18禁在线播放| 另类精品久久| 国内毛片毛片毛片毛片毛片| 伦理电影免费视频| 免费在线观看影片大全网站| 18禁国产床啪视频网站| 中亚洲国语对白在线视频| 国产亚洲欧美精品永久| 麻豆成人av在线观看| 国产在线免费精品| 国产一区二区在线观看av| 久久久久久免费高清国产稀缺| 成在线人永久免费视频| 欧美中文综合在线视频| 精品国产超薄肉色丝袜足j| 黄色视频,在线免费观看| 在线播放国产精品三级| 999精品在线视频| 成人亚洲精品一区在线观看| 国产片内射在线| 老司机靠b影院| 亚洲av第一区精品v没综合| 欧美黄色淫秽网站| 国产亚洲一区二区精品| 国产av精品麻豆| 青青草视频在线视频观看| 好男人电影高清在线观看| 性色av乱码一区二区三区2| 亚洲情色 制服丝袜| 天天影视国产精品| 美女主播在线视频| 一本综合久久免费| 国产91精品成人一区二区三区 | 少妇粗大呻吟视频| 中文欧美无线码| 国产老妇伦熟女老妇高清| 国产一区二区三区综合在线观看| 久久久久国产一级毛片高清牌| 婷婷成人精品国产| 人成视频在线观看免费观看| 一级毛片精品| 国产精品香港三级国产av潘金莲| 在线观看免费午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 叶爱在线成人免费视频播放| 9色porny在线观看| 巨乳人妻的诱惑在线观看| 女人久久www免费人成看片| 精品视频人人做人人爽| bbb黄色大片| 在线观看66精品国产| 成在线人永久免费视频| 人人澡人人妻人| 日韩视频在线欧美| 无人区码免费观看不卡 | 99国产综合亚洲精品| 国产在线精品亚洲第一网站| tube8黄色片| 精品久久久久久电影网| 欧美成人午夜精品| 日韩三级视频一区二区三区| 国产在线精品亚洲第一网站| 美女午夜性视频免费| 一区二区av电影网| 亚洲熟女毛片儿| 纵有疾风起免费观看全集完整版| 午夜福利影视在线免费观看| 午夜福利免费观看在线| 国产日韩一区二区三区精品不卡| 天天影视国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 一边摸一边抽搐一进一出视频| 国产欧美亚洲国产| svipshipincom国产片| 日本五十路高清| 青草久久国产| 岛国在线观看网站| 亚洲中文av在线| 欧美中文综合在线视频| 乱人伦中国视频| 香蕉国产在线看| 曰老女人黄片| 国产精品自产拍在线观看55亚洲 | 男女高潮啪啪啪动态图| 亚洲国产欧美一区二区综合| 亚洲精品av麻豆狂野| 国产不卡一卡二| 99久久精品国产亚洲精品| 久久婷婷成人综合色麻豆| 欧美亚洲 丝袜 人妻 在线| 欧美激情久久久久久爽电影 | cao死你这个sao货| 黑人操中国人逼视频| 天堂动漫精品| 国产aⅴ精品一区二区三区波| 免费少妇av软件| 99热网站在线观看| 午夜精品国产一区二区电影| 91国产中文字幕| 黄色成人免费大全| 亚洲精品中文字幕一二三四区 | 19禁男女啪啪无遮挡网站| 亚洲精品国产色婷婷电影| 99re在线观看精品视频| 午夜精品久久久久久毛片777| 国产精品99久久99久久久不卡| 国产精品久久久久久人妻精品电影 | 999久久久精品免费观看国产| 中文欧美无线码| 亚洲欧美精品综合一区二区三区| 亚洲国产精品一区二区三区在线| 我要看黄色一级片免费的| 亚洲精品粉嫩美女一区| 丰满人妻熟妇乱又伦精品不卡| 午夜福利视频在线观看免费| 久久香蕉激情| 在线观看免费视频日本深夜| 久久人妻福利社区极品人妻图片| xxxhd国产人妻xxx| 人妻一区二区av| svipshipincom国产片| 天天躁狠狠躁夜夜躁狠狠躁| 操出白浆在线播放| 国产日韩欧美视频二区| 人成视频在线观看免费观看| 久久精品国产综合久久久| 日本黄色日本黄色录像| 亚洲精品国产色婷婷电影| 亚洲第一av免费看| 精品国产一区二区久久| 亚洲国产欧美日韩在线播放| 久久免费观看电影| 99国产精品一区二区三区| 国产免费现黄频在线看| 亚洲欧美日韩高清在线视频 | 国产麻豆69| 丝袜在线中文字幕| 欧美亚洲 丝袜 人妻 在线| 国产精品国产高清国产av | 老汉色av国产亚洲站长工具| 欧美精品av麻豆av| 精品视频人人做人人爽| 两人在一起打扑克的视频| 三上悠亚av全集在线观看| 国产精品成人在线| 九色亚洲精品在线播放| 亚洲精品中文字幕一二三四区 | 欧美黑人欧美精品刺激| 亚洲欧洲精品一区二区精品久久久| 巨乳人妻的诱惑在线观看| 极品人妻少妇av视频| 三上悠亚av全集在线观看| 欧美日韩视频精品一区| 国产精品亚洲一级av第二区| 亚洲av欧美aⅴ国产| 欧美 亚洲 国产 日韩一| 极品教师在线免费播放| 国产av国产精品国产| 亚洲av美国av| 夫妻午夜视频| h视频一区二区三区| 777米奇影视久久| 无遮挡黄片免费观看| 三级毛片av免费| 男人舔女人的私密视频| 国产极品粉嫩免费观看在线| 成人av一区二区三区在线看| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 18在线观看网站| 麻豆成人av在线观看| 亚洲av国产av综合av卡| 一级a爱视频在线免费观看| 精品福利永久在线观看| a级毛片黄视频| 狠狠狠狠99中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 国内毛片毛片毛片毛片毛片| 9191精品国产免费久久| 午夜福利免费观看在线| 男男h啪啪无遮挡| 久久久国产精品麻豆| 亚洲中文日韩欧美视频| 成人免费观看视频高清| 啦啦啦在线免费观看视频4| 亚洲熟妇熟女久久| 看免费av毛片| 黄色视频不卡| 美女福利国产在线| 99精国产麻豆久久婷婷| 一本久久精品| 男女下面插进去视频免费观看| 少妇 在线观看| 国产黄频视频在线观看| 国产一区二区三区视频了| 国产国语露脸激情在线看| 色在线成人网| 丰满人妻熟妇乱又伦精品不卡| 美女高潮喷水抽搐中文字幕| 午夜91福利影院| av网站免费在线观看视频| 精品一区二区三区视频在线观看免费 | 91国产中文字幕| 蜜桃国产av成人99| 欧美日韩一级在线毛片| 亚洲少妇的诱惑av| 免费久久久久久久精品成人欧美视频| 欧美老熟妇乱子伦牲交| 丁香六月欧美| 叶爱在线成人免费视频播放| 岛国在线观看网站| 久久精品91无色码中文字幕| 又紧又爽又黄一区二区| 成人手机av| 天天躁夜夜躁狠狠躁躁| 亚洲成a人片在线一区二区|