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

    基于帶置信區(qū)間穩(wěn)定圖的鐵路連續(xù)梁橋模態(tài)參數(shù)識別

    2021-07-06 10:40:12王秋萍秦世強
    鐵道建筑 2021年6期
    關(guān)鍵詞:階次置信區(qū)間阻尼比

    王秋萍 秦世強

    武漢理工大學土木工程與建筑學院,武漢430070

    高速鐵路橋梁運營過程中,列車和橋梁形成一個車橋耦合系統(tǒng),系統(tǒng)的動力響應(yīng)十分復雜,會影響行車安全。為評價車橋耦合動力學行為,準確識別結(jié)構(gòu)的模態(tài)參數(shù)十分重要[1-2]。模態(tài)參數(shù)包括頻率、阻尼比和模態(tài)振型,是結(jié)構(gòu)基本的動力參數(shù),也是進行結(jié)構(gòu)動力分析、損傷識別的基本參數(shù)。與傳統(tǒng)測力法相比,環(huán)境激勵的模態(tài)參數(shù)識別方法通過測試結(jié)構(gòu)在風、水流、地脈動等環(huán)境荷載下的微幅振動響應(yīng),結(jié)合時域、頻域算法識別結(jié)構(gòu)模態(tài)參數(shù),具有測試簡便、無須額外激勵設(shè)備,不會對結(jié)構(gòu)造成損傷等優(yōu)點。因此,環(huán)境激勵的模態(tài)參數(shù)識別方法在橋梁工程領(lǐng)域得到廣泛關(guān)注和應(yīng)用[3-4]。

    隨機子空間識別(Stochastic Subspace Identification,SSI)是一種時域識別方法,是環(huán)境激勵模態(tài)參數(shù)識別中較為可靠的一種算法。西寧北川河橋[5]、岸橋[6]、Millau高架橋[7]等已經(jīng)驗證了SSI算法的實用性。基于環(huán)境激勵下的結(jié)構(gòu)動力響應(yīng),通過SSI算法識別結(jié)構(gòu)模態(tài)參數(shù)分為四個步驟:①數(shù)據(jù)采集;②數(shù)據(jù)預處理;③系統(tǒng)矩陣識別;④模態(tài)參數(shù)識別。第4步也稱作模態(tài)分析后處理,主要是利用穩(wěn)定圖方法來確定各個測試組數(shù)據(jù)識別的頻率和阻尼比,并整合各個測試組的信息識別模態(tài)振型。對于分多個測試組完成的橋梁環(huán)境振動測試,這種后處理十分繁瑣,不僅要處理多個測試組的數(shù)據(jù),而且須要區(qū)分穩(wěn)定圖中虛假模態(tài),使得整個模態(tài)分析過程效率大幅降低。

    許多學者針對這個問題開展了研究。樊江玲等[8]在穩(wěn)定圖中引入了一種可表征各階模態(tài)貢獻量的分量能量指標作為穩(wěn)定判據(jù),從而剔除虛假模態(tài)。陳永高等[9]引入余弦相似定理,提出頻率相似度和振型相似度,從而剔除穩(wěn)定圖中虛假模態(tài)。Reynders等[10]提出模態(tài)轉(zhuǎn)移范數(shù)ni的概念來區(qū)分系統(tǒng)真實模態(tài)和虛假模態(tài),ni越小,該階模態(tài)是虛假模態(tài)的可能性就越高。常軍等[11]提出一種兩階段穩(wěn)定圖的方法,首先對輸出信號進行分段構(gòu)造穩(wěn)定圖,然后對各段識別的模態(tài)參數(shù)再一次用穩(wěn)定圖方法分析,從而去除虛假模態(tài)。Cara等[12]通過研究各階模態(tài)對實測響應(yīng)的貢獻,確定了系統(tǒng)階次的選定方法,避免系統(tǒng)階次過高帶來虛假模態(tài)或過低引起模態(tài)遺漏。

    在上述研究基礎(chǔ)上,針對穩(wěn)定圖中存在的系統(tǒng)階次確定較為隨意和虛假模態(tài)問題,本文提出基于奇異熵增量的系統(tǒng)階次計算方法和帶置信區(qū)間的穩(wěn)定圖,并通過一座高速鐵路連續(xù)梁橋?qū)е眯艆^(qū)間穩(wěn)定圖的效果進行驗證。

    1 隨機子空間識別

    結(jié)構(gòu)的隨機狀態(tài)空間模型為

    式中:x k、x k+1分別為系統(tǒng)在第k、k+1個時間點處的狀態(tài)向量;y k為測試到的輸出響應(yīng)向量;w k、v k分別為測試過程誤差和系統(tǒng)建模誤差,均假定為互不相關(guān)的白噪聲;A、C分別為系統(tǒng)的狀態(tài)矩陣和輸出矩陣。

    隨機子空間識別的主要思路是通過實測輸出響應(yīng)向量y k識別系統(tǒng)矩陣(AC),結(jié)構(gòu)的模態(tài)參數(shù)可通過系統(tǒng)矩陣特征值分解獲取。有兩種方式識別系統(tǒng)矩陣:①直接從響應(yīng)數(shù)據(jù)出發(fā)構(gòu)建Hankel矩陣,引入卡爾曼濾波、正交投影和QR分解,稱為數(shù)據(jù)驅(qū)動的SSI(SSI-data);②從響應(yīng)數(shù)據(jù)的協(xié)方差出發(fā)構(gòu)建輸出協(xié)方差矩陣(Toeplitz矩陣),引入奇異值分解,稱為協(xié)方差驅(qū)動的SSI(SSI-cov)。兩種方式得到的結(jié)果近似,但計算效率不同;另外,SSI-cov識別過程更容易量化模態(tài)參數(shù)的方差[13]。因此,本文選擇SSI-cov作為試驗模態(tài)分析算法。SSI算法較成熟,不再贅述其理論過程。詳細的SSI算法實現(xiàn)過程可以參考文獻[6]。

    2 系統(tǒng)階次的確定

    系統(tǒng)階次n是SSI算法中須要確定的一個參數(shù)。由于模態(tài)以共軛對形式出現(xiàn),因此系統(tǒng)階次取偶數(shù)。系統(tǒng)階次取值過小會導致結(jié)構(gòu)真實模態(tài)遺漏,而取值過大會出現(xiàn)虛假模態(tài)。因此,合理地估計系統(tǒng)階次十分重要。

    確定系統(tǒng)階次的方法主要有兩種:

    1)奇異值分解法。通過觀察Toeplitz矩陣奇異值分解(Singular Value Decomposition,SVD)的結(jié)果來確定系統(tǒng)階次。當系統(tǒng)階次逐漸增大時,奇異值會逐漸減小至0。這種方法在數(shù)值模擬中比較精確,但實際工程中Toeplitz矩陣的奇異值并不是突變到0,而是緩緩地減小,通過觀察奇異值的變化很難確定系統(tǒng)階次。

    2)穩(wěn)定圖法?;舅悸肥墙o系統(tǒng)階次賦予一個較大的值nmax,循環(huán)求解每個系統(tǒng)階次對應(yīng)的模態(tài)參數(shù),然后以頻率為橫軸、系統(tǒng)階次為縱軸繪制穩(wěn)定圖,穩(wěn)定圖中的散點稱為穩(wěn)定點,須滿足相鄰兩階頻率的相對誤差dfi不超過1%、相鄰兩階阻尼比的相對誤差dξi不超過5%,相鄰兩階模態(tài)振型的模態(tài)置信準則(Modal Assurance Criterion,MAC)值的相對誤差dMACi不超過1%。結(jié)構(gòu)真實模態(tài)的穩(wěn)定點會逐漸形成一條穩(wěn)定軸,而虛假模態(tài)一般只會形成一些離散點。

    目前,穩(wěn)定圖方法已經(jīng)成為SSI算法的通用工具,其有效性已經(jīng)在許多工程實例中得到證明。合理地確定穩(wěn)定圖中nmax具有一定的價值。當nmax取值較大時,循環(huán)求解各階模態(tài)參數(shù)的計算時間會顯著增加,且虛假模態(tài)會逐漸形成穩(wěn)定軸,影響系統(tǒng)模態(tài)的拾取。本文在已有研究的基礎(chǔ)上,提出一種奇異熵增量加權(quán)法來估計系統(tǒng)階次。

    熵是反映系統(tǒng)不確定度的一個指標。對奇異值而言,奇異熵Ek定義為

    式中:k為奇異熵的階次;ΔEi為奇異熵階次i時的奇異熵增量,可以通過下式得到。

    式中:λi為Toeplitz矩陣的奇異值。

    奇異熵增量在低階時λi較大,在高階時λi較小,且λi隨系統(tǒng)階次的增加逐漸減小,這種變化規(guī)律不隨信號的噪聲情況而改變。已有研究利用奇異熵增量的一階導數(shù)來判定系統(tǒng)階次,當奇異熵增量的一階導數(shù)趨于穩(wěn)定時,對應(yīng)的階次即為系統(tǒng)階次。然而,這種方法仍須要作圖觀察。本文利用奇異熵增量的特性構(gòu)造奇異熵階次i處的權(quán)重系數(shù)εi。εi定義為

    計算得到系統(tǒng)階次n,即

    系統(tǒng)階次n應(yīng)取其鄰近的較大偶數(shù)。該方法簡單易行,避免了作圖觀察,且便于模態(tài)參數(shù)自動化識別。

    3 帶置信區(qū)間的穩(wěn)定圖

    穩(wěn)定圖方法中另一個問題是虛假模態(tài)。引起虛假模態(tài)的原因包括實測過程產(chǎn)生的信號噪聲以及計算假定帶來的模擬誤差。信號噪聲引起的虛假模態(tài)可以通過信號預處理、合理選定系統(tǒng)階次避免,而模擬誤差則須要量化。量化的誤差稱為模態(tài)參數(shù)的不確定度σ,能夠衡量模態(tài)參數(shù)的可信程度。在SSI算法中,模態(tài)參數(shù)最終從狀態(tài)矩陣A和輸出矩陣C獲得。因此,模態(tài)參數(shù)的不確定度由A和C的誤差導致,即識別的A和C僅僅是其真實值的估計

    為了分析系統(tǒng)狀態(tài)的變化,首先分析可觀矩陣O i的一階擾動ΔO i, 即

    式中:U、S分別為輸出協(xié)方差矩陣奇異值分解后得到的正交矩陣和對角矩陣;ΔU、ΔS分別為正交矩陣和對角矩陣的一階擾動。

    由于狀態(tài)矩陣A可以寫成可觀矩陣O i的表達式,因此,狀態(tài)矩陣A的一階擾動ΔA可以表達為

    式中:(·)?表示(·)的廣義逆。

    此時,問題轉(zhuǎn)化為如何求解正交矩陣和對角矩陣的一階擾動ΔU和ΔS。具體步驟:基于矩陣敏感性分析推導狀態(tài)矩陣的一階擾動,并量化單個測試組識別的模態(tài)參數(shù)的不確定度;構(gòu)建頻率置信區(qū)間(f±σf)和阻尼比置信區(qū)間(ξ±σξ),σf、σξ分別為頻率f和阻尼比ξ的不確定度,并將其以誤差圖的形式繪制于穩(wěn)定圖中,用于輔助判斷虛假模態(tài)。置信區(qū)間越大,表明模態(tài)參數(shù)的不確定度越高,是虛假模態(tài)的可能性就越大;反之,表明是結(jié)構(gòu)真實模態(tài)的可能性越大。

    4 連續(xù)梁橋模態(tài)參數(shù)識別

    4.1 橋梁概況

    Jalon高架橋位于西班牙境內(nèi),是一座高速鐵路連續(xù)梁橋,主橋跨徑布置為(35+45×4+35)m。主橋全長250 m,采用單箱單室等高截面雙線梁,箱梁頂寬12.94 m,底寬5.90 m,截面高3.22 m,見圖1。

    圖1 橋梁立面及橫截面

    4.2 環(huán)境振動試驗

    為研究該橋的動力特性,須對其進行環(huán)境振動測試。設(shè)置12個GeoSIG無線加速度傳感器,其中4個作為參考點傳感器,其余8個作為移動測點傳感器。每個測試組采樣時間約為15 min,采樣頻率為200 Hz;測點均布置于箱梁內(nèi)部,每個截面共布置3個測點(參見圖1),沿橋縱向每隔2.5 m設(shè)置1個測試截面,共101個測試截面,303個測點,見圖2。

    圖2 測點布置及參考點位置(單位:m)

    測試組4中參考點27處橫向加速度響應(yīng)見圖3??芍?,在環(huán)境激勵作用下橋梁呈微幅振動狀態(tài),最大橫向加速度未超過1×10-3m/s2。

    圖3 參考點27處橫向加速度響應(yīng)

    4.3 識別結(jié)果分析

    利用測試組4中加速度數(shù)據(jù)構(gòu)建輸出協(xié)方差矩陣得到奇異值分解結(jié)果和對應(yīng)的奇異熵增量,見圖4??芍浩娈愔惦S系統(tǒng)階次變化無明顯突變,很難通過觀察奇異值變化情況來確定系統(tǒng)階次;隨著系統(tǒng)階次的增加,奇異熵增量逐漸減小。按照前文方法計算得到系統(tǒng)階次為49.3,將系統(tǒng)階次設(shè)置為50時即可得到結(jié)構(gòu)的低階模態(tài)。

    圖4 奇異值與奇異熵增量隨系統(tǒng)階次變化曲線

    通過分析測試組4的數(shù)據(jù)得到帶置信區(qū)間的穩(wěn)定圖,見圖5。圖中頻率誤差線顯示了各個穩(wěn)定點的準確程度,尤其頻率為0~1 Hz時,存在較多不確定度較高的穩(wěn)定點,影響了系統(tǒng)模態(tài)的拾取。當頻率不確定度σf≤0.25、阻尼比不確定度σξ≤0.25時,圖5(a)中的虛假模態(tài)和離散點會被剔除,得到圖5(b)所示的穩(wěn)定圖,其更利于人工拾取真實模態(tài)。此外,圖5(a)中頻率為4~6 Hz時穩(wěn)定軸發(fā)生了分裂,另一個穩(wěn)定軸頻率不確定度相對較高,須進一步判定是否為結(jié)構(gòu)真實模態(tài)。通過引入置信區(qū)間,能夠較好地區(qū)分部分虛假模態(tài)和真實模態(tài),有效提高拾取穩(wěn)定軸的效率。

    圖5 帶置信區(qū)間的穩(wěn)定圖

    按照上述方法對其余測試組的數(shù)據(jù)進行處理,并獲取與圖5(b)類似的帶置信區(qū)間穩(wěn)定圖進行各測試組模態(tài)參數(shù)的融合,以獲得模態(tài)參數(shù)的試驗值。為了校核試驗值,采用ANSYS建立橋梁有限元模型,通過理論模態(tài)分析獲取模態(tài)頻率的有限元計算值。采用實體單元模擬混凝土箱梁,采用梁單元模擬橋墩;邊界條件按圖紙設(shè)置,在E-1處設(shè)置移動鉸支座,在E-2處設(shè)置固定支座;橋墩底部采用固結(jié)處理。

    橋梁模態(tài)參數(shù)識別結(jié)果見表1,其中頻率和阻尼比的試驗值為各測試組識別結(jié)果的均值。

    表1 橋梁模態(tài)參數(shù)識別結(jié)果

    由表1可知:①頻率試驗值與有限元計算值吻合良好,橫向振動模態(tài)中,第1、2階頻率有限元計算值略低于試驗值,第3、7、8階頻率有限元計算值略高于試驗值;②除第8階外,其余階次頻率有限元計算值和試驗值的相對誤差未超過5.00%;③第8階頻率有限元計算值與試驗值相對誤差為8.36%,相對低階模態(tài)較高,但未超過10.00%;④由于環(huán)境激勵對橋梁的激勵程度有限,高階模態(tài)的識別誤差相對低階模態(tài)較高。豎向模態(tài)頻率的試驗值與有限元計算值基本一致。⑤最大阻尼比不超過1.37%,阻尼較低。

    5 結(jié)論

    1)基于奇異熵增量的系統(tǒng)階次計算方法簡單易行,避免通過觀察奇異值的變化情況來確定系統(tǒng)階次,且能夠在程序編制中直接實現(xiàn),利于模態(tài)識別過程自動化。

    2)帶置信區(qū)間的穩(wěn)定圖能夠有效區(qū)分系統(tǒng)真實模態(tài)和虛假模態(tài)。由于信號受實測過程誤差等因素影響,對有些虛假模態(tài)須進一步判定,但帶置信區(qū)間的穩(wěn)定圖仍能去除部分虛假模態(tài)。改進的穩(wěn)定圖方法能夠有效協(xié)助用戶進行模態(tài)參數(shù)識別后處理,提高SSI算法后處理的效率。

    3)帶置信區(qū)間的穩(wěn)定圖識別出的橋梁模態(tài)參數(shù)試驗值與有限元計算值吻合良好,表明帶置信區(qū)間的穩(wěn)定圖能有效應(yīng)用于工程結(jié)構(gòu)模態(tài)參數(shù)識別。

    猜你喜歡
    階次置信區(qū)間阻尼比
    定數(shù)截尾場合三參數(shù)pareto分布參數(shù)的最優(yōu)置信區(qū)間
    p-范分布中參數(shù)的置信區(qū)間
    基于細觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    多個偏正態(tài)總體共同位置參數(shù)的Bootstrap置信區(qū)間
    階次分析在驅(qū)動橋異響中的應(yīng)用
    基于Vold-Kalman濾波的階次分析系統(tǒng)設(shè)計與實現(xiàn)*
    黏滯阻尼器在時程分析下的附加有效阻尼比研究
    振動與沖擊(2019年4期)2019-02-22 02:33:34
    列車定位中置信區(qū)間的確定方法
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    基于齒輪階次密度優(yōu)化的變速器降噪研究
    價值工程(2017年28期)2018-01-23 20:48:29
    最近手机中文字幕大全| 国产1区2区3区精品| 在线观看免费高清a一片| 亚洲精品第二区| 亚洲精品一区蜜桃| 国产在线视频一区二区| 国产精品无大码| √禁漫天堂资源中文www| 热re99久久国产66热| 97精品久久久久久久久久精品| 男人舔女人的私密视频| 黄色一级大片看看| 亚洲熟女精品中文字幕| 久久99热这里只频精品6学生| 成年人午夜在线观看视频| 精品亚洲乱码少妇综合久久| 少妇高潮的动态图| 免费在线观看黄色视频的| 成人18禁高潮啪啪吃奶动态图| 亚洲国产色片| 日韩视频在线欧美| 熟女电影av网| 少妇的丰满在线观看| 午夜免费鲁丝| 国产一级毛片在线| 成年人免费黄色播放视频| 2022亚洲国产成人精品| 国产精品女同一区二区软件| 另类精品久久| 亚洲在久久综合| 丝袜人妻中文字幕| 久久久久人妻精品一区果冻| 国产日韩欧美视频二区| 亚洲丝袜综合中文字幕| 国产永久视频网站| 国产视频首页在线观看| 欧美xxxx性猛交bbbb| 午夜影院在线不卡| 亚洲一码二码三码区别大吗| 日韩av免费高清视频| 欧美 亚洲 国产 日韩一| 高清在线视频一区二区三区| 亚洲三级黄色毛片| 免费人成在线观看视频色| 亚洲av国产av综合av卡| 熟女av电影| 久久狼人影院| 久久久久久久亚洲中文字幕| 丁香六月天网| 黑人巨大精品欧美一区二区蜜桃 | 久热这里只有精品99| 成年女人在线观看亚洲视频| kizo精华| 宅男免费午夜| 18禁国产床啪视频网站| 在线天堂中文资源库| 人成视频在线观看免费观看| 丰满乱子伦码专区| 99热6这里只有精品| 欧美变态另类bdsm刘玥| 亚洲一区二区三区欧美精品| 免费播放大片免费观看视频在线观看| 国产男人的电影天堂91| 欧美国产精品va在线观看不卡| 青春草国产在线视频| 亚洲国产精品国产精品| 精品视频人人做人人爽| 一级爰片在线观看| 午夜福利乱码中文字幕| 国产不卡av网站在线观看| 欧美日本中文国产一区发布| 国产欧美亚洲国产| 五月开心婷婷网| 国产白丝娇喘喷水9色精品| 国产av精品麻豆| 在线天堂最新版资源| 亚洲av日韩在线播放| 亚洲av国产av综合av卡| 久久久久视频综合| 免费观看无遮挡的男女| 日韩 亚洲 欧美在线| 国产69精品久久久久777片| 午夜影院在线不卡| 免费av不卡在线播放| 免费高清在线观看视频在线观看| 少妇的逼好多水| 日韩免费高清中文字幕av| 美女国产高潮福利片在线看| www日本在线高清视频| 免费黄色在线免费观看| 免费黄色在线免费观看| 超碰97精品在线观看| 三级国产精品片| 内地一区二区视频在线| 亚洲图色成人| 欧美97在线视频| 国产精品不卡视频一区二区| 日韩av在线免费看完整版不卡| 欧美精品国产亚洲| av片东京热男人的天堂| 色婷婷久久久亚洲欧美| 两个人免费观看高清视频| 在线天堂最新版资源| 精品久久久久久电影网| 最近手机中文字幕大全| 久久精品熟女亚洲av麻豆精品| 国产精品.久久久| 欧美人与性动交α欧美精品济南到 | 国产精品一二三区在线看| 国产成人精品一,二区| 在线看a的网站| av免费在线看不卡| 欧美国产精品va在线观看不卡| 久久这里有精品视频免费| 欧美最新免费一区二区三区| 高清av免费在线| 久久人人爽人人片av| 国产淫语在线视频| 2022亚洲国产成人精品| av在线播放精品| 亚洲少妇的诱惑av| 在线天堂最新版资源| 久久久久久伊人网av| 美女福利国产在线| 国产在线一区二区三区精| 欧美激情极品国产一区二区三区 | 午夜精品国产一区二区电影| 宅男免费午夜| 久久精品久久久久久噜噜老黄| 久久精品国产综合久久久 | 亚洲精品色激情综合| 亚洲综合色网址| 成人无遮挡网站| 2021少妇久久久久久久久久久| 国产爽快片一区二区三区| 久久久久精品人妻al黑| 亚洲国产av影院在线观看| 久久青草综合色| 韩国高清视频一区二区三区| 人人妻人人爽人人添夜夜欢视频| 69精品国产乱码久久久| 亚洲婷婷狠狠爱综合网| 免费少妇av软件| 大话2 男鬼变身卡| 国产亚洲午夜精品一区二区久久| 99热6这里只有精品| 久久久久国产网址| 亚洲国产精品专区欧美| a级毛色黄片| 熟女av电影| 视频区图区小说| 久久 成人 亚洲| 午夜精品国产一区二区电影| 香蕉丝袜av| 午夜久久久在线观看| av女优亚洲男人天堂| 精品视频人人做人人爽| 欧美日韩综合久久久久久| 18在线观看网站| 亚洲三级黄色毛片| 免费黄网站久久成人精品| 一级爰片在线观看| 久久毛片免费看一区二区三区| 国产成人精品福利久久| 岛国毛片在线播放| 欧美日韩一区二区视频在线观看视频在线| 日本黄大片高清| 汤姆久久久久久久影院中文字幕| 波多野结衣一区麻豆| 亚洲美女视频黄频| 99热国产这里只有精品6| 亚洲内射少妇av| 国精品久久久久久国模美| 亚洲精品中文字幕在线视频| 国产免费现黄频在线看| xxx大片免费视频| 嫩草影院入口| 国产精品麻豆人妻色哟哟久久| 亚洲第一区二区三区不卡| 18禁裸乳无遮挡动漫免费视频| 大香蕉97超碰在线| 国产精品久久久av美女十八| 婷婷色麻豆天堂久久| 制服人妻中文乱码| 性色av一级| 亚洲国产日韩一区二区| 最新中文字幕久久久久| 中文字幕精品免费在线观看视频 | 国产精品偷伦视频观看了| 中国国产av一级| 妹子高潮喷水视频| 大陆偷拍与自拍| 少妇的丰满在线观看| 老司机亚洲免费影院| h视频一区二区三区| 男女边摸边吃奶| 男女边吃奶边做爰视频| 美女视频免费永久观看网站| 母亲3免费完整高清在线观看 | 国产国语露脸激情在线看| 春色校园在线视频观看| 亚洲,欧美,日韩| 色吧在线观看| 一二三四在线观看免费中文在 | 五月伊人婷婷丁香| 午夜福利影视在线免费观看| 国产成人精品一,二区| 国产精品一国产av| 日日爽夜夜爽网站| 久久久久久久亚洲中文字幕| 亚洲精品久久久久久婷婷小说| 久久人人爽av亚洲精品天堂| 伊人久久国产一区二区| 满18在线观看网站| 欧美国产精品一级二级三级| 18禁动态无遮挡网站| 亚洲经典国产精华液单| 国产国拍精品亚洲av在线观看| 一本久久精品| 亚洲国产日韩一区二区| 亚洲av男天堂| 51国产日韩欧美| 赤兔流量卡办理| 午夜av观看不卡| 观看美女的网站| 欧美97在线视频| 日产精品乱码卡一卡2卡三| 国产免费福利视频在线观看| 国产精品欧美亚洲77777| 国产精品一二三区在线看| 欧美激情极品国产一区二区三区 | 80岁老熟妇乱子伦牲交| 免费高清在线观看视频在线观看| 咕卡用的链子| 国产精品国产av在线观看| 成人黄色视频免费在线看| 欧美日韩视频高清一区二区三区二| 亚洲美女搞黄在线观看| 中文天堂在线官网| 国精品久久久久久国模美| 丁香六月天网| av网站免费在线观看视频| 国产精品一区二区在线不卡| 亚洲精品一区蜜桃| 久久精品久久久久久噜噜老黄| 欧美 亚洲 国产 日韩一| 成人午夜精彩视频在线观看| 看非洲黑人一级黄片| videossex国产| 久久精品熟女亚洲av麻豆精品| 国产精品偷伦视频观看了| 亚洲av电影在线观看一区二区三区| 国产亚洲av片在线观看秒播厂| 午夜福利网站1000一区二区三区| 免费观看无遮挡的男女| 男女无遮挡免费网站观看| 欧美亚洲 丝袜 人妻 在线| 欧美性感艳星| 女人被躁到高潮嗷嗷叫费观| av播播在线观看一区| 精品少妇内射三级| 亚洲精华国产精华液的使用体验| 国产精品秋霞免费鲁丝片| 国产精品嫩草影院av在线观看| 两个人免费观看高清视频| 免费在线观看黄色视频的| 高清欧美精品videossex| 不卡视频在线观看欧美| 亚洲av.av天堂| a 毛片基地| 亚洲久久久国产精品| 亚洲伊人久久精品综合| 成年人午夜在线观看视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久精品人人爽人人爽视色| www.熟女人妻精品国产 | av片东京热男人的天堂| 久久99热6这里只有精品| 最近最新中文字幕大全免费视频 | 日本91视频免费播放| 免费观看a级毛片全部| 大香蕉97超碰在线| 欧美日本中文国产一区发布| 天堂8中文在线网| 乱码一卡2卡4卡精品| 男人操女人黄网站| 在线 av 中文字幕| 成年人免费黄色播放视频| 黑人高潮一二区| 亚洲成人一二三区av| 亚洲成人av在线免费| 国产一区亚洲一区在线观看| 日韩成人av中文字幕在线观看| 午夜老司机福利剧场| 热re99久久精品国产66热6| 亚洲人成网站在线观看播放| 精品酒店卫生间| 超色免费av| 亚洲欧洲日产国产| 一区二区三区乱码不卡18| 熟女人妻精品中文字幕| 久久国产精品男人的天堂亚洲 | 制服人妻中文乱码| 水蜜桃什么品种好| 亚洲国产毛片av蜜桃av| 人人妻人人澡人人爽人人夜夜| 在线观看免费日韩欧美大片| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品一二三| 国产国语露脸激情在线看| 亚洲人与动物交配视频| www.av在线官网国产| 国产亚洲最大av| 极品少妇高潮喷水抽搐| 欧美日韩成人在线一区二区| 色婷婷久久久亚洲欧美| 亚洲人成77777在线视频| 精品少妇久久久久久888优播| 成人综合一区亚洲| 精品久久久久久电影网| 韩国精品一区二区三区 | 飞空精品影院首页| 婷婷色麻豆天堂久久| 99九九在线精品视频| 亚洲中文av在线| 男女国产视频网站| 国产av一区二区精品久久| 国产成人91sexporn| 亚洲国产成人一精品久久久| av在线观看视频网站免费| 亚洲四区av| 国产av码专区亚洲av| 亚洲情色 制服丝袜| 午夜福利视频在线观看免费| 日韩三级伦理在线观看| 男人舔女人的私密视频| 国产欧美日韩综合在线一区二区| 国产精品一二三区在线看| 亚洲av中文av极速乱| 性色avwww在线观看| 18禁国产床啪视频网站| 国产精品女同一区二区软件| 在线观看免费视频网站a站| 自拍欧美九色日韩亚洲蝌蚪91| 狠狠精品人妻久久久久久综合| 免费观看a级毛片全部| 欧美日韩av久久| 边亲边吃奶的免费视频| 国产在视频线精品| av.在线天堂| 国产免费视频播放在线视频| 国产成人精品福利久久| 天堂8中文在线网| 日韩视频在线欧美| 亚洲欧洲日产国产| 日韩视频在线欧美| 五月玫瑰六月丁香| 国产乱人偷精品视频| 亚洲美女视频黄频| 国产精品蜜桃在线观看| 男女边摸边吃奶| 久久精品夜色国产| 男女免费视频国产| 精品一区二区三区视频在线| 99视频精品全部免费 在线| 黑人欧美特级aaaaaa片| 久久久精品免费免费高清| 日本vs欧美在线观看视频| 亚洲伊人色综图| av播播在线观看一区| 妹子高潮喷水视频| 亚洲av国产av综合av卡| 99热全是精品| 国产永久视频网站| 国语对白做爰xxxⅹ性视频网站| 国产白丝娇喘喷水9色精品| 免费大片18禁| 亚洲在久久综合| 男女高潮啪啪啪动态图| 在线亚洲精品国产二区图片欧美| 日本色播在线视频| 久久精品夜色国产| 国产成人精品无人区| 好男人视频免费观看在线| 巨乳人妻的诱惑在线观看| 一级片'在线观看视频| 新久久久久国产一级毛片| 亚洲天堂av无毛| 97人妻天天添夜夜摸| 成人亚洲欧美一区二区av| 成人无遮挡网站| 亚洲精品美女久久久久99蜜臀 | 中文乱码字字幕精品一区二区三区| 天堂中文最新版在线下载| 精品少妇内射三级| 看免费成人av毛片| 在线观看美女被高潮喷水网站| 国产一区二区三区综合在线观看 | 亚洲国产色片| 欧美日韩一区二区视频在线观看视频在线| 午夜91福利影院| 久久久欧美国产精品| 成人18禁高潮啪啪吃奶动态图| 夜夜爽夜夜爽视频| 日韩三级伦理在线观看| 乱码一卡2卡4卡精品| 99热国产这里只有精品6| av线在线观看网站| 欧美日韩精品成人综合77777| 亚洲图色成人| 免费高清在线观看视频在线观看| 如何舔出高潮| 亚洲天堂av无毛| a级毛色黄片| 欧美日韩视频精品一区| 国产精品久久久久久精品古装| 日韩熟女老妇一区二区性免费视频| 亚洲内射少妇av| 免费少妇av软件| 国产一区二区三区av在线| 在现免费观看毛片| 免费高清在线观看日韩| 免费观看无遮挡的男女| 九草在线视频观看| 久久综合国产亚洲精品| 久久狼人影院| 亚洲精品国产色婷婷电影| 精品视频人人做人人爽| 国产国拍精品亚洲av在线观看| 国产色爽女视频免费观看| 精品午夜福利在线看| 伊人久久国产一区二区| 国产精品秋霞免费鲁丝片| 亚洲国产欧美日韩在线播放| 国产精品偷伦视频观看了| 国产精品蜜桃在线观看| 人人妻人人澡人人看| 日本91视频免费播放| 久久人人爽人人爽人人片va| 天天影视国产精品| 丰满乱子伦码专区| 宅男免费午夜| 免费观看无遮挡的男女| 国产成人精品福利久久| 欧美 亚洲 国产 日韩一| a级毛片在线看网站| 日韩 亚洲 欧美在线| 久久99蜜桃精品久久| 久久99热6这里只有精品| 视频在线观看一区二区三区| 中文字幕制服av| 欧美精品一区二区免费开放| 女人被躁到高潮嗷嗷叫费观| 国产精品偷伦视频观看了| 久久久久久久国产电影| 亚洲天堂av无毛| 97在线人人人人妻| 少妇的丰满在线观看| 欧美精品高潮呻吟av久久| 视频区图区小说| 午夜影院在线不卡| 午夜激情久久久久久久| 亚洲精品久久成人aⅴ小说| 青青草视频在线视频观看| 十八禁网站网址无遮挡| 亚洲欧美日韩卡通动漫| 久久av网站| 免费观看性生交大片5| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜久久久在线观看| 少妇 在线观看| 日本午夜av视频| 欧美成人精品欧美一级黄| 菩萨蛮人人尽说江南好唐韦庄| 日日啪夜夜爽| 国产精品三级大全| 国产精品国产av在线观看| 久久97久久精品| 亚洲国产日韩一区二区| 热99国产精品久久久久久7| 国产av码专区亚洲av| 51国产日韩欧美| 2022亚洲国产成人精品| 国产精品国产三级国产专区5o| 日本-黄色视频高清免费观看| 2021少妇久久久久久久久久久| 在线天堂最新版资源| 免费黄色在线免费观看| 婷婷色麻豆天堂久久| 啦啦啦啦在线视频资源| 免费高清在线观看视频在线观看| 建设人人有责人人尽责人人享有的| 久久99精品国语久久久| 欧美成人午夜精品| 亚洲国产精品专区欧美| 午夜福利视频在线观看免费| 亚洲av电影在线观看一区二区三区| 宅男免费午夜| 精品久久久精品久久久| 最近最新中文字幕大全免费视频 | 欧美亚洲 丝袜 人妻 在线| 免费黄网站久久成人精品| 又粗又硬又长又爽又黄的视频| 伦精品一区二区三区| 美国免费a级毛片| 熟女电影av网| 一边摸一边做爽爽视频免费| 波多野结衣一区麻豆| 欧美亚洲日本最大视频资源| 亚洲美女视频黄频| 亚洲图色成人| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| 丰满迷人的少妇在线观看| 男女无遮挡免费网站观看| 久久99热这里只频精品6学生| 少妇高潮的动态图| av天堂久久9| 97超碰精品成人国产| 日韩成人av中文字幕在线观看| www.av在线官网国产| 一二三四中文在线观看免费高清| 寂寞人妻少妇视频99o| 曰老女人黄片| 亚洲精品乱码久久久久久按摩| 99热这里只有是精品在线观看| 亚洲av中文av极速乱| 另类精品久久| 草草在线视频免费看| 97在线人人人人妻| 日韩大片免费观看网站| 免费在线观看完整版高清| 啦啦啦中文免费视频观看日本| av有码第一页| 嫩草影院入口| 亚洲国产精品一区二区三区在线| 黄色视频在线播放观看不卡| 看非洲黑人一级黄片| 乱码一卡2卡4卡精品| 精品久久久久久电影网| 国产精品人妻久久久久久| 综合色丁香网| 亚洲av日韩在线播放| 亚洲国产精品成人久久小说| 观看av在线不卡| 亚洲少妇的诱惑av| 久久影院123| 亚洲在久久综合| 亚洲高清免费不卡视频| 久久国产精品大桥未久av| 欧美日韩视频精品一区| 国产欧美日韩一区二区三区在线| 一边亲一边摸免费视频| 欧美激情 高清一区二区三区| 最近手机中文字幕大全| 天天躁夜夜躁狠狠久久av| 99久久人妻综合| 色婷婷av一区二区三区视频| 国产精品麻豆人妻色哟哟久久| 久久久久久久精品精品| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 国产黄频视频在线观看| 内地一区二区视频在线| 国产精品欧美亚洲77777| 国产在线视频一区二区| 男人爽女人下面视频在线观看| 久久精品国产自在天天线| 97在线人人人人妻| 国产综合精华液| 国产免费现黄频在线看| 麻豆乱淫一区二区| 我的女老师完整版在线观看| 少妇的逼水好多| 一本久久精品| 欧美精品人与动牲交sv欧美| 在线观看www视频免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人午夜福利电影在线观看| 人人妻人人添人人爽欧美一区卜| 永久免费av网站大全| 欧美丝袜亚洲另类| 亚洲欧洲国产日韩| 国产精品不卡视频一区二区| 日日爽夜夜爽网站| 午夜福利在线观看免费完整高清在| 香蕉精品网在线| 日本vs欧美在线观看视频| 日韩欧美一区视频在线观看| 在线观看免费高清a一片| av.在线天堂| 免费日韩欧美在线观看| 男的添女的下面高潮视频| 国产在线一区二区三区精| 午夜久久久在线观看| 欧美激情 高清一区二区三区| 老司机亚洲免费影院| 黄片无遮挡物在线观看| 日韩制服丝袜自拍偷拍| 久久影院123| 亚洲av日韩在线播放| 久热久热在线精品观看| 97在线视频观看| 在线天堂中文资源库| 好男人视频免费观看在线| 如日韩欧美国产精品一区二区三区| 涩涩av久久男人的天堂| 亚洲欧洲日产国产| 亚洲成av片中文字幕在线观看 | 在线观看国产h片| 中文字幕另类日韩欧美亚洲嫩草| 狂野欧美激情性bbbbbb| 中国三级夫妇交换| 久久这里有精品视频免费| 卡戴珊不雅视频在线播放|