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

    基于應(yīng)變模態(tài)響應(yīng)重構(gòu)的損傷識(shí)別方法

    2022-09-03 03:56:56鄒云峰盧玄東陽勁松何旭輝
    工程力學(xué) 2022年9期
    關(guān)鍵詞:幅值重構(gòu)模態(tài)

    鄒云峰,盧玄東,陽勁松,何旭輝

    (1. 中南大學(xué)土木工程學(xué)院,湖南,長沙 410075;2. 中南大學(xué)交通運(yùn)輸工程學(xué)院,湖南,長沙 410075;3. 軌道交通工程結(jié)構(gòu)防災(zāi)減災(zāi)湖南省重點(diǎn)實(shí)驗(yàn)室,湖南,長沙 410075)

    土木結(jié)構(gòu)的性能在其運(yùn)營階段會(huì)隨著時(shí)間推移而下降,從而影響結(jié)構(gòu)的安全性和耐久性,土木結(jié)構(gòu)健康監(jiān)測技術(shù)一直是國內(nèi)外的研究熱點(diǎn)[1-2]。土木結(jié)構(gòu)損傷形式多樣、損傷位置隱蔽,由于土木結(jié)構(gòu)形式復(fù)雜、尺度大等特點(diǎn),很難對(duì)土木結(jié)構(gòu)損傷進(jìn)行準(zhǔn)確識(shí)別,若放任損傷不斷發(fā)展,很有可能導(dǎo)致重大安全事故,因此開展土木結(jié)構(gòu)損傷識(shí)別研究具有重要意義。

    現(xiàn)有結(jié)構(gòu)損傷識(shí)別方法的基本原理是,損傷會(huì)導(dǎo)致結(jié)構(gòu)物理特性的變化,該變化會(huì)反映在實(shí)際測量的結(jié)構(gòu)動(dòng)態(tài)響應(yīng)中,通過分析結(jié)構(gòu)動(dòng)態(tài)響應(yīng)的變化實(shí)現(xiàn)損傷識(shí)別[3-5]。損傷識(shí)別方法可大致分為3 種:一是頻域方法[6-8],主要識(shí)別參數(shù)包括頻率、振型和曲率模態(tài)等,該方法需要在頻域上識(shí)別模態(tài)參數(shù),具有一定的復(fù)雜性且難以保證計(jì)算精度;二是時(shí)域方法,如基于響應(yīng)靈敏度的識(shí)別方法[9-10],此種方法需要通過外部激勵(lì)和結(jié)構(gòu)動(dòng)態(tài)響應(yīng)的關(guān)系來進(jìn)行損傷識(shí)別,具有數(shù)據(jù)量充足、識(shí)別精度高的特點(diǎn),但需要準(zhǔn)確的外荷載信息,而實(shí)際結(jié)構(gòu)的外荷載通常難以獲?。蝗菚r(shí)頻域方法,如基于小波包靈敏度的識(shí)別方法[11-12],此種方法同樣需要獲取外荷載信息,而且在時(shí)頻域上進(jìn)行的正向運(yùn)算和反向運(yùn)算,其計(jì)算效率較低,過程繁復(fù)。

    針對(duì)以上問題,本文將動(dòng)態(tài)響應(yīng)重構(gòu)方法引入損傷識(shí)別中,通過已有的測量信息去獲取全局的結(jié)構(gòu)動(dòng)態(tài)響應(yīng),從而解決需要已知外部載荷的局限。Ribeiro 等[13]針對(duì)多自由度系統(tǒng)提出一種基于廣義傳遞矩陣的響應(yīng)重構(gòu)法,在頻域內(nèi)建立兩組響應(yīng)的傳遞矩陣,進(jìn)而實(shí)現(xiàn)響應(yīng)待測點(diǎn)的重構(gòu)工作;Law 等[14-16]基于頻響函數(shù)將上述位移傳遞率推廣至加速度傳遞率,完成了全模型或模型子結(jié)構(gòu)上的響應(yīng)重構(gòu)。He 等[17-19]提出一種基于經(jīng)驗(yàn)?zāi)B(tài)分解的時(shí)域響應(yīng)重構(gòu)方法,該方法通過結(jié)構(gòu)模態(tài)特性構(gòu)建傳遞矩陣,結(jié)合經(jīng)驗(yàn)?zāi)B(tài)分解法和模態(tài)疊加法得到未安裝傳感器位置的動(dòng)力響應(yīng),相對(duì)于頻域重構(gòu)方法,該方法具有更高的計(jì)算效率。目前,時(shí)域重構(gòu)方法僅用于無損傷結(jié)構(gòu)的響應(yīng)重構(gòu),若結(jié)構(gòu)出現(xiàn)損傷而引起結(jié)構(gòu)參數(shù)改變,使用原本的傳遞函數(shù)重構(gòu)得到的響應(yīng)必然與真實(shí)響應(yīng)存在差別,通過利用這種差異性結(jié)合有限元更新方法可以實(shí)現(xiàn)對(duì)土木結(jié)構(gòu)的損傷識(shí)別[20-21]。

    本文結(jié)合基于經(jīng)驗(yàn)?zāi)B(tài)分解的響應(yīng)重構(gòu)方法,提出一種基于應(yīng)變模態(tài)響應(yīng)重構(gòu)的損傷識(shí)別方法。該方法根據(jù)傳感器采集的響應(yīng)信息直接分解得到的模態(tài)響應(yīng)和使用無損傷模型的應(yīng)變模態(tài)重構(gòu)得到的模態(tài)響應(yīng)的差異作為有限元模型修正的依據(jù),通過應(yīng)變模態(tài)比值構(gòu)建的傳遞率的靈敏度矩陣進(jìn)行迭代運(yùn)算,求得損傷位置及損傷程度。由于該方法不需要獲取外部激勵(lì)信息,且在時(shí)域中進(jìn)行,無需進(jìn)行頻域的相關(guān)計(jì)算和時(shí)頻域轉(zhuǎn)換,關(guān)鍵信息的獲取方便快捷。最后,通過連續(xù)梁的單一損傷識(shí)別和多損傷識(shí)別驗(yàn)證了該方法的正確性和有效性,并討論了不同等級(jí)的測量噪聲和模態(tài)階次的選取對(duì)識(shí)別結(jié)果的影響。

    1 基于應(yīng)變模態(tài)響應(yīng)重構(gòu)的損傷識(shí)別理論推導(dǎo)

    1.1 基于經(jīng)驗(yàn)?zāi)B(tài)分解的應(yīng)變重構(gòu)方法

    使用經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)法處理時(shí)域信號(hào)的方法分解傳感器的測量響應(yīng),可以將測量響應(yīng)近似表示為:

    1.2 基于應(yīng)變模態(tài)響應(yīng)重構(gòu)的有限元模型修正

    將損傷識(shí)別問題表述為使實(shí)際結(jié)構(gòu)和有限元模型之間有關(guān)結(jié)構(gòu)特性的信息差異最小值問題:

    基于應(yīng)變模態(tài)響應(yīng)重構(gòu)的損傷識(shí)別方法,其基本流程如圖1 所示。

    圖1 損傷識(shí)別流程圖Fig. 1 Flowchart of damage identification method

    2 數(shù)值模擬分析

    2.1 連續(xù)梁損傷識(shí)別

    通過連續(xù)梁仿真案例以驗(yàn)證本文提出的基于應(yīng)變模態(tài)響應(yīng)重構(gòu)的損傷識(shí)別方法。用于數(shù)值模擬的連續(xù)梁有限元模型如圖2 所示,該連續(xù)梁彈性模量為210 GPa,密度為7850 kg/m3,長、寬、高分別為3 m、0.06 m、0.03 m。采用ANSYS 建立有限元模型,單元類型為beam3 單元,共有31 個(gè)節(jié)點(diǎn),30 個(gè)單元,每個(gè)單元長為0.1 m,從左往右分別給節(jié)點(diǎn)編號(hào)1~31,單元編號(hào)1~30。無損傷狀態(tài)下連續(xù)梁的前5 階自振頻率分別為20.85 Hz、32.48 Hz、83.38 Hz、105.32 Hz、187.61 Hz。對(duì)該連續(xù)梁模型設(shè)置多組損傷工況,損傷工況主要考慮連續(xù)梁容易發(fā)生損傷的部位,如跨中和支座附近部位,以此來驗(yàn)證所提出的損傷識(shí)別方法的可行性,損傷工況的設(shè)置如表1 所示。

    表1 損傷工況Table 1 Damage condition

    圖2 連續(xù)梁模型Fig. 2 Model of continuous beam

    由理論推導(dǎo)可知,本文損傷識(shí)別方法與外界激勵(lì)無關(guān),將幅值在一定范圍內(nèi)的豎向沖擊荷載施加于連續(xù)梁的8 號(hào)節(jié)點(diǎn)上。應(yīng)變片布置在每個(gè)單元上部中心位置,采樣頻率為1 kHz,如圖2 所示。選取結(jié)構(gòu)自由振動(dòng)期間1 s 的響應(yīng)信號(hào),提取一階模態(tài)響應(yīng)并根據(jù)式(9)~式(11)得到差異向量。在這里需要注意的是,對(duì)于多組應(yīng)變片采集的數(shù)據(jù),可以有不同種重構(gòu)方式組合,在該數(shù)值模擬中,經(jīng)過多次對(duì)比發(fā)現(xiàn),重構(gòu)方式組合只會(huì)在存在噪聲的情況下對(duì)識(shí)別結(jié)果產(chǎn)生影響,且該影響較小,在此不予討論。

    為了提高識(shí)別精度,迭代求解加入了自適應(yīng)Tikhonov 正則化方法[22],識(shí)別結(jié)果如圖3 所示。

    圖3 損傷識(shí)別結(jié)果Fig. 3 Damage identification results

    從圖3 可以看出,預(yù)設(shè)損傷的單元的識(shí)別結(jié)果與實(shí)際預(yù)設(shè)損傷量較吻合度很高,并且未預(yù)設(shè)損傷的單元上的錯(cuò)誤識(shí)別很小,可以忽略不計(jì)。無論是單個(gè)損失還是多個(gè)損失的情況,識(shí)別結(jié)果都令人滿意。

    2.2 噪聲對(duì)損傷識(shí)別的影響

    結(jié)構(gòu)響應(yīng)信號(hào)的測量不可避免地存在噪聲,在本研究中,將白噪聲添加到每一個(gè)單元的應(yīng)變中,以研究噪聲對(duì)損傷識(shí)別的影響。測量噪聲是高斯脈沖過程,其均方根設(shè)置為應(yīng)變最大均方根的百分比。這里使用噪聲等級(jí)來表征噪聲程度,即添加到應(yīng)變中的噪聲是應(yīng)變最大均方根值等同于噪聲等級(jí)的高斯過程。結(jié)合連續(xù)梁數(shù)值模擬案例,分析不同噪聲等級(jí)對(duì)損傷識(shí)別的影響(噪聲等級(jí)分別為5%、10%、和15%)。

    為了更好地對(duì)比說明噪聲對(duì)該案例中連續(xù)梁的損傷識(shí)別的影響,對(duì)無損傷的連續(xù)梁作不同噪聲等級(jí)下的損傷識(shí)別??紤]到高斯白噪聲的隨機(jī)性,在研究噪聲對(duì)損傷識(shí)別的影響時(shí),對(duì)損傷工況1、2、5 和6 各設(shè)置1000 個(gè)對(duì)照組,每個(gè)對(duì)照組隨機(jī)生成不同的應(yīng)變測量噪聲。由于迭代求解中加入了自適應(yīng)Tikhonov 正則化方法,各個(gè)單元損傷因子的識(shí)別值均小于0。用伽馬分布擬合全部1000 個(gè)對(duì)照組識(shí)別結(jié)果中的各個(gè)非預(yù)設(shè)損傷單元損傷因子識(shí)別值,用高斯分布擬合預(yù)設(shè)損傷單元的損傷因子識(shí)別值,分別計(jì)算其95%置信區(qū)間。損傷工況1、2、5 和6 在不同噪聲等級(jí)下的各單元識(shí)別損傷因子的95%置信區(qū)間如圖4~圖6所示。

    圖4 5%噪聲下的損傷識(shí)別結(jié)果Fig. 4 Damage identification results under 5% noise

    圖5 10%噪聲下的損傷識(shí)別結(jié)果Fig. 5 Damage identification results under 10% noise

    圖6 15%噪聲下的損傷識(shí)別結(jié)果Fig. 6 Damage identification results under 15% noise

    從圖4~圖6 可以看出,在不同的損傷工況中,噪聲引起識(shí)別結(jié)果的波動(dòng)也不相同,但對(duì)于預(yù)設(shè)損傷值為-0.05 和-0.2 而言,噪聲的影響程度差別不大,并不會(huì)隨著預(yù)設(shè)損傷的增大而增大;由損傷工況1、5、6 的識(shí)別結(jié)果可以看出,隨著噪聲等級(jí)的增大,識(shí)別結(jié)果的波動(dòng)明顯增大,可見噪聲對(duì)識(shí)別存在一定影響。但即使噪聲等級(jí)提升到15%,識(shí)別結(jié)果的95%置信區(qū)間范圍也不超過0.015,對(duì)于損傷工況1、5、6 而言是完全可以接受的,識(shí)別結(jié)果依舊能很好地反映實(shí)際損傷情況,而對(duì)于損傷工況2 而言,15%噪聲的影響可以忽略不計(jì),可見所提出的損傷識(shí)別方法對(duì)噪聲有較好的魯棒性。

    2.3 不同模態(tài)對(duì)損傷識(shí)別的影響

    前文對(duì)連續(xù)梁的損傷識(shí)別都依據(jù)結(jié)構(gòu)自由振動(dòng)期間的響應(yīng)進(jìn)行,對(duì)于該連續(xù)梁而言,自由振動(dòng)表現(xiàn)為單頻振動(dòng),而有時(shí)候土木結(jié)構(gòu)響應(yīng)可能受到連續(xù)加持的外荷載作用(如風(fēng)荷載、流水荷載等),此時(shí)結(jié)構(gòu)響應(yīng)通常包含不同的模態(tài)成分,且本文所提出的方法是基于結(jié)構(gòu)的模態(tài)信息進(jìn)行損傷識(shí)別的,所以要研究不同模態(tài)對(duì)識(shí)別結(jié)果的影響。為了使結(jié)構(gòu)響應(yīng)能包含多個(gè)模態(tài)成分,將之前施加的沖擊荷載改為隨機(jī)荷載,該隨機(jī)荷載由6 階低通巴特沃斯濾波器過濾的均值大于0 的高斯過程模擬,結(jié)構(gòu)響應(yīng)的各階模態(tài)響應(yīng)可以反映在結(jié)構(gòu)響應(yīng)的頻譜圖上,以16 號(hào)單元應(yīng)變?yōu)槔涓道锶~譜如圖7 所示。根據(jù)圖7,這里主要選取有著較明顯幅值的前4 階模態(tài)響應(yīng)來討論模態(tài)階次的選取對(duì)識(shí)別的影響。

    圖7 16 號(hào)單元應(yīng)變數(shù)據(jù)的傅里葉譜Fig. 7 Fourier spectrum of strain data of Element 16

    需要注意的是,對(duì)于一些應(yīng)變模態(tài)接近0 的響應(yīng),使用EMD 方法提取其模態(tài)響應(yīng)將產(chǎn)生很大的誤差,從而無法按照式(5)和式(11)進(jìn)行近似替代,因此,對(duì)于結(jié)構(gòu)中某一階應(yīng)變模態(tài)為0 或接近0 的位置的損傷情況,無法選取這一階應(yīng)變模態(tài)用于識(shí)別該位置的損傷。在避免出現(xiàn)上述情況的前提下,對(duì)表1 中的6 個(gè)損傷工況選用不同的模態(tài)階次進(jìn)行識(shí)別以研究選取的模態(tài)對(duì)損傷識(shí)別的影響。對(duì)比選用各階模態(tài)的識(shí)別結(jié)果,多個(gè)模態(tài)組合如表2 所示。在表2 中,組合1 和組合2 選取的單元模態(tài)響應(yīng)幅值都相對(duì)較大,而組合3 中包含了多個(gè)幅值較小的模態(tài)響應(yīng),這里模態(tài)響應(yīng)的幅值可以通過各單元同階應(yīng)變模態(tài)值的相對(duì)大小體現(xiàn)出來。例如,圖7 中16 號(hào)單元應(yīng)變的第一、三、四階模態(tài)響應(yīng)相對(duì)于第二階模態(tài)響應(yīng)而言較小。損傷識(shí)別結(jié)果如圖8 所示。

    表2 模態(tài)組合Table 2 Damage condition

    在圖8 中,模態(tài)組合1、2 的識(shí)別結(jié)果區(qū)別不大。而模態(tài)組合3 包含了多個(gè)幅值較小的模態(tài)響應(yīng),其識(shí)別結(jié)果較差,因此應(yīng)盡量選擇幅值較大的模態(tài)進(jìn)行損傷識(shí)別。此外,相較于圖3 中的識(shí)別結(jié)果而言,在隨機(jī)荷載作用下的識(shí)別誤差有所增大,但在本案例中,即使對(duì)于預(yù)設(shè)損傷因子為-0.05 的單元而言,這樣的誤差依然是可以接受的。

    圖8 不同模態(tài)組合下的損傷識(shí)別結(jié)果Fig. 8 Damage identification results under different modal combinations

    3 結(jié)論

    本文提出了一種基于應(yīng)變模態(tài)響應(yīng)重構(gòu)的損傷識(shí)別方法,開展了連續(xù)梁單一損傷和多損傷識(shí)別仿真計(jì)算研究,并探討了測量噪聲和模態(tài)階次選取對(duì)識(shí)別結(jié)果的影響,得主要結(jié)論如下:

    (1) 該識(shí)別方法在時(shí)域上實(shí)現(xiàn),不需要進(jìn)行模態(tài)參數(shù)識(shí)別和時(shí)頻分析,也不需要獲取結(jié)構(gòu)外部激勵(lì)信息,通過將傳感器采集的應(yīng)變進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解分解,得到的模態(tài)時(shí)域響應(yīng)信息可直接用于損傷識(shí)別,通過連續(xù)梁仿真案例驗(yàn)證了該方法的高效性和可行性。

    (2) 噪聲影響的研究表明,噪聲對(duì)損傷識(shí)別影響很小,對(duì)于該案例中的連續(xù)梁結(jié)構(gòu)而言,噪聲等級(jí)達(dá)到15%,依然能得到較好的識(shí)別結(jié)果。在實(shí)際工程中,應(yīng)變測量的噪聲等級(jí)通常不超過10%,該損傷識(shí)別方法在實(shí)際運(yùn)用中對(duì)測量噪聲具有較強(qiáng)的魯棒性。

    (3) 應(yīng)變模態(tài)接近0 的模態(tài)響應(yīng)無法進(jìn)行響應(yīng)重構(gòu),從而不能用于損傷識(shí)別。除了上述情況外,模態(tài)階次選取影響的研究表明,幅值較大的模態(tài)更有利于識(shí)別損傷;并且對(duì)比使用隨機(jī)荷載下結(jié)構(gòu)振動(dòng)響應(yīng)進(jìn)行損傷識(shí)別,使用自由振動(dòng)下的結(jié)構(gòu)響應(yīng)識(shí)別能得到更好的結(jié)果。

    (4) 本文提出的方法,通過對(duì)比重構(gòu)信號(hào)與實(shí)際信號(hào)的差異,通過更新有限元模型獲取結(jié)構(gòu)的損傷。在損傷識(shí)別過程中,需要在損傷單元布貼有傳感器,才能實(shí)現(xiàn)對(duì)損傷的識(shí)別,因此,本方法更適用于對(duì)大型結(jié)構(gòu)敏感區(qū)域的檢測。

    猜你喜歡
    幅值重構(gòu)模態(tài)
    長城敘事的重構(gòu)
    攝影世界(2022年1期)2022-01-21 10:50:14
    北方大陸 重構(gòu)未來
    北京的重構(gòu)與再造
    商周刊(2017年6期)2017-08-22 03:42:36
    基于S變換的交流電網(wǎng)幅值檢測系統(tǒng)計(jì)算機(jī)仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    論中止行為及其對(duì)中止犯的重構(gòu)
    正序電壓幅值檢測及諧波抑制的改進(jìn)
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    低壓電力線信道脈沖噪聲的幅值與寬度特征
    基于零序電壓幅值增量的消弧線圈調(diào)諧新方法
    電測與儀表(2015年7期)2015-04-09 11:40:10
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
    欧美日韩亚洲国产一区二区在线观看 | 男女午夜视频在线观看| 在线观看www视频免费| 亚洲欧美日韩另类电影网站| 黄色视频在线播放观看不卡| 成年人午夜在线观看视频| 在线 av 中文字幕| a 毛片基地| 国产黄色视频一区二区在线观看| 老司机在亚洲福利影院| 欧美成人午夜精品| 飞空精品影院首页| 亚洲av成人不卡在线观看播放网 | 亚洲精品国产色婷婷电影| 搡老乐熟女国产| 看免费av毛片| 大片免费播放器 马上看| 婷婷色综合www| 欧美日韩av久久| 国产精品一二三区在线看| 国产欧美日韩一区二区三区在线| 成年av动漫网址| av在线老鸭窝| 80岁老熟妇乱子伦牲交| 高清视频免费观看一区二区| 女人久久www免费人成看片| 国产精品av久久久久免费| 精品福利永久在线观看| 久久久亚洲精品成人影院| 国产精品成人在线| 久久久久精品久久久久真实原创| 国产97色在线日韩免费| 国产精品三级大全| 综合色丁香网| 亚洲国产欧美网| 老司机亚洲免费影院| 亚洲欧美激情在线| 免费黄网站久久成人精品| 国产精品麻豆人妻色哟哟久久| 少妇被粗大的猛进出69影院| 亚洲精品美女久久av网站| 亚洲成人av在线免费| 国产福利在线免费观看视频| 亚洲人成电影观看| 日本一区二区免费在线视频| 国产一区亚洲一区在线观看| 成人18禁高潮啪啪吃奶动态图| 国产99久久九九免费精品| 大片免费播放器 马上看| 亚洲欧美中文字幕日韩二区| 免费观看性生交大片5| avwww免费| 国产av一区二区精品久久| 国产野战对白在线观看| 韩国高清视频一区二区三区| 久久精品亚洲av国产电影网| 一本一本久久a久久精品综合妖精| 精品一品国产午夜福利视频| 一二三四在线观看免费中文在| 99国产综合亚洲精品| 国产精品久久久久久精品古装| 国产深夜福利视频在线观看| 美女中出高潮动态图| 亚洲,欧美,日韩| 天天添夜夜摸| 色精品久久人妻99蜜桃| 欧美日韩视频精品一区| 久久久久视频综合| 亚洲情色 制服丝袜| 美国免费a级毛片| 高清欧美精品videossex| 欧美成人午夜精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品国产av蜜桃| 韩国精品一区二区三区| 久久精品熟女亚洲av麻豆精品| a级片在线免费高清观看视频| 男人操女人黄网站| xxx大片免费视频| av在线播放精品| 老司机亚洲免费影院| 亚洲成人一二三区av| 成年动漫av网址| 欧美日韩综合久久久久久| 国产精品国产三级专区第一集| 99久久精品国产亚洲精品| 久久久精品94久久精品| 午夜福利视频精品| 免费日韩欧美在线观看| 青春草国产在线视频| 999久久久国产精品视频| 久久久国产欧美日韩av| 久久久欧美国产精品| 国产高清不卡午夜福利| 成年美女黄网站色视频大全免费| 这个男人来自地球电影免费观看 | 国精品久久久久久国模美| 国产片特级美女逼逼视频| 久热爱精品视频在线9| 狠狠精品人妻久久久久久综合| 99久久人妻综合| 午夜免费男女啪啪视频观看| 亚洲av电影在线进入| 制服丝袜香蕉在线| 99精品久久久久人妻精品| 一本色道久久久久久精品综合| 国产精品一区二区在线不卡| 啦啦啦在线免费观看视频4| 中文乱码字字幕精品一区二区三区| 久久久久久久精品精品| 欧美中文综合在线视频| 视频区图区小说| 超色免费av| 男女之事视频高清在线观看 | 婷婷色综合www| 日本av手机在线免费观看| 免费久久久久久久精品成人欧美视频| 男女免费视频国产| 精品少妇黑人巨大在线播放| 免费少妇av软件| 婷婷成人精品国产| 下体分泌物呈黄色| 亚洲婷婷狠狠爱综合网| 国产亚洲最大av| 久久这里只有精品19| 欧美另类一区| 亚洲综合精品二区| 国产亚洲av高清不卡| 精品亚洲成国产av| 宅男免费午夜| av.在线天堂| 国产熟女午夜一区二区三区| 亚洲av成人精品一二三区| 亚洲精品国产一区二区精华液| 人人妻人人添人人爽欧美一区卜| 午夜免费鲁丝| 综合色丁香网| 又黄又粗又硬又大视频| 韩国av在线不卡| 久久精品国产a三级三级三级| 欧美97在线视频| 久久久国产精品麻豆| 成人手机av| 777米奇影视久久| 国产亚洲最大av| 一级毛片我不卡| 久久精品久久精品一区二区三区| 国产精品免费大片| 一级毛片 在线播放| 亚洲欧洲国产日韩| 91aial.com中文字幕在线观看| 男的添女的下面高潮视频| 国产av精品麻豆| 黄色一级大片看看| 在线精品无人区一区二区三| 欧美变态另类bdsm刘玥| 新久久久久国产一级毛片| 国产爽快片一区二区三区| 蜜桃国产av成人99| 秋霞在线观看毛片| 曰老女人黄片| 一级毛片我不卡| 9191精品国产免费久久| 免费观看性生交大片5| 久久久久网色| 97精品久久久久久久久久精品| 日韩制服丝袜自拍偷拍| 秋霞伦理黄片| 黄色一级大片看看| 最新的欧美精品一区二区| 久久久精品国产亚洲av高清涩受| 中文字幕最新亚洲高清| 久久久久精品性色| 久久久久久人妻| 免费黄频网站在线观看国产| 一区福利在线观看| 亚洲av日韩精品久久久久久密 | 捣出白浆h1v1| 天堂中文最新版在线下载| 欧美人与性动交α欧美精品济南到| 亚洲精品国产色婷婷电影| 好男人视频免费观看在线| 欧美精品人与动牲交sv欧美| 男女午夜视频在线观看| 国产av码专区亚洲av| 爱豆传媒免费全集在线观看| 国产免费一区二区三区四区乱码| 久久久久久人妻| 捣出白浆h1v1| 香蕉国产在线看| 国产精品一区二区在线不卡| 亚洲三区欧美一区| 久久毛片免费看一区二区三区| 一级,二级,三级黄色视频| 亚洲成国产人片在线观看| 少妇猛男粗大的猛烈进出视频| 一边摸一边做爽爽视频免费| 亚洲美女搞黄在线观看| 免费看不卡的av| 美女大奶头黄色视频| 国产午夜精品一二区理论片| 久久久久精品人妻al黑| 汤姆久久久久久久影院中文字幕| 国产精品无大码| xxx大片免费视频| 老司机亚洲免费影院| 国产在线一区二区三区精| 女性生殖器流出的白浆| 超色免费av| 亚洲精品久久午夜乱码| 欧美 亚洲 国产 日韩一| 精品亚洲乱码少妇综合久久| 国产成人精品福利久久| 成人黄色视频免费在线看| 激情五月婷婷亚洲| 又黄又粗又硬又大视频| 丝袜脚勾引网站| 日本91视频免费播放| 一区福利在线观看| 一本久久精品| 亚洲成人手机| 超碰成人久久| av不卡在线播放| 在线免费观看不下载黄p国产| 日本vs欧美在线观看视频| 少妇人妻 视频| 狂野欧美激情性xxxx| 两个人看的免费小视频| 天堂8中文在线网| 久久免费观看电影| 国产成人精品久久二区二区91 | www.av在线官网国产| 一区二区三区精品91| 男女边吃奶边做爰视频| 久久精品亚洲熟妇少妇任你| 国产在线一区二区三区精| 丝袜美腿诱惑在线| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲av高清不卡| 夫妻性生交免费视频一级片| 亚洲av电影在线观看一区二区三区| 精品亚洲乱码少妇综合久久| 亚洲av国产av综合av卡| 热99国产精品久久久久久7| 国产伦理片在线播放av一区| 久久久久国产一级毛片高清牌| 少妇人妻精品综合一区二区| 久久毛片免费看一区二区三区| 丰满迷人的少妇在线观看| 亚洲精品av麻豆狂野| 一边亲一边摸免费视频| 精品国产乱码久久久久久小说| 性色av一级| 亚洲少妇的诱惑av| 99热全是精品| 一二三四中文在线观看免费高清| 日韩欧美精品免费久久| 热99国产精品久久久久久7| 亚洲专区中文字幕在线 | 一本一本久久a久久精品综合妖精| 久久99精品国语久久久| 久久久久人妻精品一区果冻| 国产成人欧美| 亚洲,欧美,日韩| 成人黄色视频免费在线看| 精品免费久久久久久久清纯 | 一个人免费看片子| 性色av一级| 国产xxxxx性猛交| 老汉色av国产亚洲站长工具| 成人手机av| 毛片一级片免费看久久久久| 中文字幕亚洲精品专区| 亚洲,欧美,日韩| 最近最新中文字幕免费大全7| 丝袜脚勾引网站| 欧美日本中文国产一区发布| 中文字幕色久视频| 亚洲成人一二三区av| 精品少妇黑人巨大在线播放| 亚洲熟女精品中文字幕| 亚洲综合精品二区| 亚洲人成网站在线观看播放| 免费在线观看黄色视频的| 女性被躁到高潮视频| 搡老乐熟女国产| 久久99热这里只频精品6学生| 亚洲图色成人| 久久免费观看电影| av卡一久久| 另类亚洲欧美激情| 亚洲精品,欧美精品| 好男人视频免费观看在线| 久久毛片免费看一区二区三区| 精品国产一区二区久久| 久久久久国产精品人妻一区二区| 欧美日本中文国产一区发布| 国产亚洲欧美精品永久| 国产一级毛片在线| 美女脱内裤让男人舔精品视频| 国产亚洲最大av| 伊人久久大香线蕉亚洲五| 欧美最新免费一区二区三区| 啦啦啦在线观看免费高清www| 美女国产高潮福利片在线看| 人人妻人人爽人人添夜夜欢视频| 成人国产麻豆网| 人妻人人澡人人爽人人| 另类亚洲欧美激情| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产av成人精品| 国产精品久久久人人做人人爽| 观看av在线不卡| 亚洲欧美一区二区三区久久| av视频免费观看在线观看| av免费观看日本| 婷婷色综合www| 国产一卡二卡三卡精品 | 青春草亚洲视频在线观看| 老司机在亚洲福利影院| 精品少妇黑人巨大在线播放| 日本wwww免费看| 亚洲欧美一区二区三区久久| 女人高潮潮喷娇喘18禁视频| 国产一区二区激情短视频 | 亚洲av国产av综合av卡| 只有这里有精品99| 免费久久久久久久精品成人欧美视频| 亚洲自偷自拍图片 自拍| 十八禁高潮呻吟视频| 久久av网站| 久久性视频一级片| 搡老乐熟女国产| 国产亚洲一区二区精品| 精品一品国产午夜福利视频| 亚洲成av片中文字幕在线观看| 丰满乱子伦码专区| 亚洲欧美激情在线| 午夜影院在线不卡| 18在线观看网站| 人人妻人人添人人爽欧美一区卜| 亚洲精品日本国产第一区| 欧美变态另类bdsm刘玥| 日本午夜av视频| 国产免费视频播放在线视频| 在线亚洲精品国产二区图片欧美| 精品酒店卫生间| 日本wwww免费看| 欧美 亚洲 国产 日韩一| 欧美在线黄色| 男女免费视频国产| 美女扒开内裤让男人捅视频| 日韩 亚洲 欧美在线| 国产一区二区三区综合在线观看| 丰满乱子伦码专区| 精品视频人人做人人爽| 国产97色在线日韩免费| 老熟女久久久| 操出白浆在线播放| 天堂俺去俺来也www色官网| 最近最新中文字幕大全免费视频 | 精品久久久久久电影网| 国产精品蜜桃在线观看| 街头女战士在线观看网站| 亚洲精品日本国产第一区| videosex国产| 国产成人av激情在线播放| 国产日韩欧美在线精品| 桃花免费在线播放| 亚洲伊人色综图| 成人黄色视频免费在线看| 9热在线视频观看99| 国产精品一区二区在线观看99| 搡老乐熟女国产| avwww免费| 天天躁夜夜躁狠狠久久av| 日日撸夜夜添| 天天躁夜夜躁狠狠久久av| 一级毛片 在线播放| 好男人视频免费观看在线| 精品人妻一区二区三区麻豆| 亚洲欧美精品自产自拍| 亚洲伊人久久精品综合| 中国三级夫妇交换| 一级a爱视频在线免费观看| 宅男免费午夜| av在线观看视频网站免费| 丰满少妇做爰视频| 日韩欧美一区视频在线观看| 成人黄色视频免费在线看| 久久久久精品国产欧美久久久 | 亚洲一卡2卡3卡4卡5卡精品中文| 综合色丁香网| 国产精品99久久99久久久不卡 | 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| 欧美日韩亚洲综合一区二区三区_| 丰满迷人的少妇在线观看| 十八禁高潮呻吟视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av电影在线进入| 黑人猛操日本美女一级片| 欧美成人午夜精品| 亚洲成人av在线免费| 成人影院久久| 飞空精品影院首页| 欧美成人精品欧美一级黄| 一本大道久久a久久精品| 国产精品一区二区在线观看99| 亚洲av电影在线观看一区二区三区| 青春草视频在线免费观看| 亚洲美女视频黄频| 亚洲欧美中文字幕日韩二区| www.av在线官网国产| 国产成人免费观看mmmm| 欧美亚洲日本最大视频资源| 最新的欧美精品一区二区| 国产精品久久久久久精品电影小说| 岛国毛片在线播放| 亚洲欧美成人综合另类久久久| 99热全是精品| 大话2 男鬼变身卡| 午夜91福利影院| 日本黄色日本黄色录像| 国产成人啪精品午夜网站| 成年av动漫网址| 香蕉国产在线看| 少妇 在线观看| 肉色欧美久久久久久久蜜桃| 女人精品久久久久毛片| 女人被躁到高潮嗷嗷叫费观| 国产女主播在线喷水免费视频网站| 国产精品女同一区二区软件| 国产亚洲av高清不卡| a级毛片黄视频| 一本一本久久a久久精品综合妖精| 久久久久精品久久久久真实原创| 亚洲欧美中文字幕日韩二区| 国产成人精品在线电影| 国产一区二区 视频在线| 国产色婷婷99| 国产国语露脸激情在线看| 久热爱精品视频在线9| av.在线天堂| 精品少妇内射三级| a级片在线免费高清观看视频| 色综合欧美亚洲国产小说| 亚洲国产精品国产精品| 在线亚洲精品国产二区图片欧美| 大片电影免费在线观看免费| 国产一级毛片在线| 亚洲精品自拍成人| xxxhd国产人妻xxx| 伊人久久国产一区二区| 91国产中文字幕| 一二三四中文在线观看免费高清| 亚洲精品av麻豆狂野| 欧美精品高潮呻吟av久久| 69精品国产乱码久久久| 99久久综合免费| 国产精品久久久人人做人人爽| 国产探花极品一区二区| 成人免费观看视频高清| 国产97色在线日韩免费| 热99国产精品久久久久久7| 少妇人妻精品综合一区二区| 国产在线一区二区三区精| 国产精品久久久久久人妻精品电影 | 免费人妻精品一区二区三区视频| 国产激情久久老熟女| 国产亚洲一区二区精品| 欧美激情高清一区二区三区 | 欧美av亚洲av综合av国产av | 妹子高潮喷水视频| 男女下面插进去视频免费观看| 99九九在线精品视频| 老汉色av国产亚洲站长工具| 免费av中文字幕在线| 国产亚洲av片在线观看秒播厂| 夫妻性生交免费视频一级片| 中文欧美无线码| 欧美日韩国产mv在线观看视频| 日韩欧美一区视频在线观看| 国产一区二区在线观看av| 蜜桃国产av成人99| 亚洲七黄色美女视频| av网站在线播放免费| 大片电影免费在线观看免费| 久久久久久久大尺度免费视频| 在线看a的网站| 啦啦啦在线免费观看视频4| 考比视频在线观看| 成人三级做爰电影| 国产成人系列免费观看| 精品一区二区三区四区五区乱码 | 国产激情久久老熟女| 国产淫语在线视频| 国产精品人妻久久久影院| 韩国精品一区二区三区| 日韩伦理黄色片| 国产成人欧美在线观看 | 超色免费av| 国产片特级美女逼逼视频| 亚洲国产精品一区二区三区在线| 女性被躁到高潮视频| 99国产综合亚洲精品| av卡一久久| 美女高潮到喷水免费观看| 久久 成人 亚洲| 亚洲,欧美精品.| 久久综合国产亚洲精品| 久久亚洲国产成人精品v| 中国国产av一级| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美成人综合另类久久久| 自线自在国产av| 视频区图区小说| 乱人伦中国视频| 亚洲免费av在线视频| 中文字幕制服av| 亚洲人成电影观看| 97精品久久久久久久久久精品| 成年美女黄网站色视频大全免费| 女人爽到高潮嗷嗷叫在线视频| 色94色欧美一区二区| 日日爽夜夜爽网站| 国产伦理片在线播放av一区| 精品一区在线观看国产| 赤兔流量卡办理| 久久久精品94久久精品| av视频免费观看在线观看| 男女国产视频网站| 国产深夜福利视频在线观看| 精品视频人人做人人爽| 桃花免费在线播放| www.自偷自拍.com| 极品人妻少妇av视频| 亚洲欧美精品综合一区二区三区| 成人免费观看视频高清| 黄色毛片三级朝国网站| 国产在线视频一区二区| 一区二区三区四区激情视频| 欧美激情极品国产一区二区三区| 亚洲色图综合在线观看| 欧美日本中文国产一区发布| 国精品久久久久久国模美| 大香蕉久久网| 亚洲成人手机| 国产av精品麻豆| 日本欧美国产在线视频| 亚洲熟女精品中文字幕| 搡老岳熟女国产| 国产精品女同一区二区软件| av.在线天堂| 久久久久久免费高清国产稀缺| 最近最新中文字幕免费大全7| 熟女av电影| 美国免费a级毛片| 只有这里有精品99| 日本91视频免费播放| 久久人人爽av亚洲精品天堂| 男女之事视频高清在线观看 | 亚洲欧美清纯卡通| 亚洲七黄色美女视频| 色婷婷久久久亚洲欧美| av在线老鸭窝| 一级,二级,三级黄色视频| 少妇人妻久久综合中文| 涩涩av久久男人的天堂| 一级片'在线观看视频| 亚洲成人国产一区在线观看 | 交换朋友夫妻互换小说| 一本大道久久a久久精品| 久久女婷五月综合色啪小说| 青青草视频在线视频观看| 日韩,欧美,国产一区二区三区| 男男h啪啪无遮挡| 国产一区二区三区综合在线观看| 日本av手机在线免费观看| 成人亚洲精品一区在线观看| 欧美乱码精品一区二区三区| 国产有黄有色有爽视频| 一本久久精品| 免费女性裸体啪啪无遮挡网站| 99九九在线精品视频| 纵有疾风起免费观看全集完整版| 午夜日韩欧美国产| 亚洲人成77777在线视频| 精品卡一卡二卡四卡免费| 天天操日日干夜夜撸| 热re99久久精品国产66热6| 国产麻豆69| 丝袜喷水一区| 男女国产视频网站| 黄色视频不卡| 人妻人人澡人人爽人人| 中文天堂在线官网| 亚洲三区欧美一区| 久久久精品国产亚洲av高清涩受| 精品一区二区三卡| 91老司机精品| 18禁国产床啪视频网站| 黄色毛片三级朝国网站| 亚洲精品美女久久av网站| 精品一区二区三区四区五区乱码 | 建设人人有责人人尽责人人享有的| 婷婷色综合大香蕉| 欧美日韩国产mv在线观看视频| 中文天堂在线官网| 另类精品久久| 你懂的网址亚洲精品在线观看| 色婷婷久久久亚洲欧美| 国产精品国产三级专区第一集| 精品少妇久久久久久888优播| 老司机影院毛片|