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

    考慮環(huán)境因素影響的海洋平臺(tái)結(jié)構(gòu)損傷檢測(cè)研究

    2021-09-27 08:18:02吳文開(kāi)徐明強(qiáng)王樹(shù)青蔣玉峰王國(guó)興
    振動(dòng)與沖擊 2021年16期
    關(guān)鍵詞:協(xié)整殘差閾值

    吳文開(kāi),徐明強(qiáng),王樹(shù)青,蔣玉峰,王國(guó)興

    (1.中國(guó)海洋大學(xué) 工程學(xué)院,山東 青島 266100;2.山東交通學(xué)院 船舶與港口工程學(xué)院,山東 威海 264200)

    海洋平臺(tái)體積龐大、結(jié)構(gòu)復(fù)雜,且長(zhǎng)期服役于惡劣的海洋環(huán)境下,結(jié)構(gòu)損傷逐漸累積。為保證結(jié)構(gòu)作業(yè)的安全性與耐久性,對(duì)海洋平臺(tái)進(jìn)行健康監(jiān)測(cè)是必不可少的[1]。由于具有全局、自動(dòng)化檢測(cè)的能力,基于結(jié)構(gòu)振動(dòng)特性的損傷檢測(cè)技術(shù)在結(jié)構(gòu)健康監(jiān)測(cè)領(lǐng)域得到了廣泛重視,學(xué)者們發(fā)展出了一系列的損傷檢測(cè)方法[2]。其基本思想是:損傷會(huì)改變結(jié)構(gòu)的物理屬性(例如剛度),進(jìn)而影響結(jié)構(gòu)的動(dòng)力特性,反過(guò)來(lái),利用結(jié)構(gòu)的動(dòng)力特性構(gòu)建損傷靈敏度指標(biāo),即可實(shí)現(xiàn)損傷判定[3]。然而,常用的損傷指標(biāo),如結(jié)構(gòu)頻率,極易受到海洋環(huán)境因素(如溫度、潮汐、海生物附著等)變化的影響,以致在實(shí)際應(yīng)用中難以準(zhǔn)確反映出結(jié)構(gòu)損傷[4]。為了將損傷檢測(cè)方法更好地推廣到工程實(shí)踐,發(fā)展能夠消除環(huán)境因素影響的損傷檢測(cè)技術(shù)至關(guān)重要。

    溫度是影響結(jié)構(gòu)剛度變化的一種典型環(huán)境要素[5]。Askegaard等[6]在對(duì)一座三跨鋼筋混凝土人行橋?yàn)槠?年的監(jiān)測(cè)中發(fā)現(xiàn),季節(jié)性的溫度變化對(duì)橋梁頻率的改變達(dá)到了10%。Doebling等[7]對(duì)美國(guó)阿拉莫斯峽谷大橋進(jìn)行的模態(tài)測(cè)試則表明,氣溫日變化引起大橋1階頻率的改變接近5%。然而,F(xiàn)arrar等[8]對(duì)I-40橋進(jìn)行破壞性試驗(yàn)時(shí)發(fā)現(xiàn),將橋梁一側(cè)的工字梁沿橫截面切割1/2,其1階頻率的改變僅為8%。不難看出,環(huán)境因素的變化對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)的影響較為明顯,在一定程度上會(huì)掩蓋結(jié)構(gòu)真實(shí)損傷,從而導(dǎo)致?lián)p傷誤判和損傷漏判的發(fā)生。

    已有的考慮環(huán)境因素影響的損傷檢測(cè)方法主要分為兩大類(lèi):一類(lèi)是建立環(huán)境因素與結(jié)構(gòu)動(dòng)力響應(yīng)之間的相關(guān)性模型,Peeters等[9]建立了ARX模型,通過(guò)測(cè)量溫度的變化預(yù)測(cè)結(jié)構(gòu)頻率,當(dāng)預(yù)測(cè)頻率與識(shí)別頻率不吻合時(shí)認(rèn)為結(jié)構(gòu)發(fā)生了損傷,然而此類(lèi)方法要求環(huán)境因素是可測(cè)的,但對(duì)海洋結(jié)構(gòu)而言,由于其服役環(huán)境相當(dāng)惡劣,通常不具備可測(cè)量條件;第二類(lèi)方法考慮了環(huán)境因素不可測(cè)量或難以測(cè)量的情況,此類(lèi)方法將結(jié)構(gòu)動(dòng)力響應(yīng)分解成結(jié)構(gòu)損傷和環(huán)境因素影響兩個(gè)部分,僅需要響應(yīng)信息,是近年來(lái)的研究熱點(diǎn)。應(yīng)用較為廣泛的有主成分分析(principal component analysis,PCA)和協(xié)整分析(cointegration analysis,CA)等。

    Yan等[10]考慮了線(xiàn)性或弱非線(xiàn)性環(huán)境因素變化的影響,首次將通過(guò)PCA降維處理后得到的殘差作為損傷指標(biāo),以一個(gè)三跨橋的有限元模型和一個(gè)木橋的物理模型試驗(yàn)對(duì)方法的有效性進(jìn)行了驗(yàn)證。隨后,Yan等[11]融合了兩種新型聚類(lèi)策略,將PCA方法進(jìn)一步推廣到處理環(huán)境因素的非線(xiàn)性影響。吳森等[12]先利用PCA消除一鋼結(jié)構(gòu)平臺(tái)動(dòng)態(tài)響應(yīng)中的溫度影響,繼而以小波包系數(shù)節(jié)點(diǎn)能量譜計(jì)算結(jié)構(gòu)損傷敏感特征來(lái)識(shí)別結(jié)構(gòu)損傷。常鵬等[13]采用結(jié)構(gòu)響應(yīng)的小波包能量譜作為特征參數(shù)的輸入,以主成分殘差作為損傷指標(biāo),通過(guò)一個(gè)藏式古建筑兩年的實(shí)測(cè)數(shù)據(jù)驗(yàn)證了該方法可以剔除溫度變化的影響。Wang等[14]采用PCA方法消除振型數(shù)據(jù)中的環(huán)境影響,進(jìn)而構(gòu)造反映結(jié)構(gòu)真實(shí)狀態(tài)的殘差應(yīng)變能,通過(guò)多變量假設(shè)檢驗(yàn)進(jìn)行損傷判定,對(duì)某一海上風(fēng)機(jī)現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)的分析表明,該方法可有效避免誤判問(wèn)題。

    Cross等[15]首先提出采用響應(yīng)數(shù)據(jù)的協(xié)整殘差作為損傷判定指標(biāo),并以一個(gè)溫度變化條件下的層合板損傷試驗(yàn)驗(yàn)證了該方法的有效性。由于實(shí)際信號(hào)可能在不同的時(shí)間尺度下具有不同的共同趨勢(shì),Worden等[16]先對(duì)信號(hào)進(jìn)行多分辨分析,進(jìn)而求解各時(shí)間尺度下分解信號(hào)的協(xié)整殘差,在一定程度上提高了協(xié)整方法的損傷檢測(cè)靈敏度。梁亞斌等[17]采用EG(engle-granger)兩步協(xié)整求解結(jié)構(gòu)前2階頻率的協(xié)整殘差,通過(guò)鋼筋混凝土梁和鋼桁架橋的數(shù)值模擬表明該方法可有效消除溫度影響且具有一定的噪聲魯棒性。刁延松等[18]以測(cè)點(diǎn)響應(yīng)的AR(autoregressive model)模型系數(shù)作為協(xié)整變量,通過(guò)一個(gè)海洋平臺(tái)模型的實(shí)驗(yàn)證明了所提方法可以有效消除溫度變化的影響。Huang等[19]發(fā)展了一種基于卡爾曼濾波和協(xié)整的損傷識(shí)別方法,先通過(guò)協(xié)整系數(shù)構(gòu)建卡爾曼濾波的狀態(tài)向量,而后利用遞歸過(guò)程在線(xiàn)估計(jì)結(jié)構(gòu)狀態(tài)的變化,將該方法應(yīng)用于天津永和大橋,成功地識(shí)別出了兩處結(jié)構(gòu)損傷。

    PCA和CA均為多元統(tǒng)計(jì)分析的重要內(nèi)容,兩者頗有類(lèi)似之處。然而,在健康監(jiān)測(cè)領(lǐng)域尚未對(duì)兩種方法進(jìn)行過(guò)比較研究。本文以一個(gè)海洋平臺(tái)模型為例,考慮空氣、海水和海底泥土溫度變化的聯(lián)合影響,以結(jié)構(gòu)的模態(tài)頻率作為樣本數(shù)據(jù),以X-bar控制圖作為損傷判別標(biāo)準(zhǔn),分析比較兩種方法在損傷識(shí)別效果以及噪聲魯棒性方面的性能差異,為海洋平臺(tái)損傷檢測(cè)方法的選取提供參考借鑒。

    1 主成分分析(PCA)

    1.1 基本原理

    PCA又稱(chēng)“Karhunen-Loeve變換”、“本征正交分解”,是一種多元統(tǒng)計(jì)分析方法。PCA通過(guò)求解樣本數(shù)據(jù)協(xié)方差矩陣的特征值問(wèn)題,以前幾階特征向量構(gòu)造的主成分代替樣本數(shù)據(jù),達(dá)到了簡(jiǎn)化問(wèn)題和減少計(jì)算量的目的。

    取結(jié)構(gòu)n階模態(tài)頻率的N個(gè)觀測(cè)值,對(duì)頻率作中心化處理,組成N×n階測(cè)量頻率矩陣X,其協(xié)方差矩陣可表示為

    (1)

    則主成分分析的實(shí)質(zhì)是求解以下特征值問(wèn)題

    CΦ=ΛΦ

    (2)

    式中:Λ=diag(λ1,λ2,…,λn)為由n個(gè)特征值組成的對(duì)角矩陣,且λ1≥λ2≥…≥λn;Φ=[φ1,φ2,…,φn]為與特征值對(duì)應(yīng)的特征向量所組成的特征向量矩陣。

    設(shè)定能量閾值ε,使得對(duì)于前k(k≤n)個(gè)特征值,滿(mǎn)足

    (3)

    則前k個(gè)特征向量為主成分向量,可構(gòu)成投影矩陣

    T=[φ1,φ2,…,φk]

    (4)

    進(jìn)而得到原頻率矩陣基于低維空間的估計(jì)

    (5)

    定義殘差矩陣

    (6)

    則第t個(gè)觀測(cè)值對(duì)應(yīng)的殘差向量為et=[et1,et2,…,etn],對(duì)其求Euclid范數(shù),可得到描述觀測(cè)值與估計(jì)值之間偏離程度的主成分殘差

    (7)

    1.2 X-bar控制圖

    (8)

    由于主成分分析已消除了環(huán)境因素的影響,因此正態(tài)分布假設(shè)下基準(zhǔn)數(shù)據(jù)對(duì)應(yīng)的主成分殘差位于區(qū)間[NI-3σ,NI+3σ]內(nèi)的概率為99.74%。

    (9)

    此外,測(cè)試數(shù)據(jù)和基準(zhǔn)數(shù)據(jù)的殘差均值比

    (10)

    也可以作為損傷判別標(biāo)準(zhǔn),該比值趨于1,則結(jié)構(gòu)正常,當(dāng)殘差均值比較大時(shí)可以認(rèn)為結(jié)構(gòu)發(fā)生了損傷。

    考慮到樣本數(shù)目的有限性,用于定量描述結(jié)構(gòu)損傷的異常值比例不宜設(shè)置過(guò)小,本文以異常值比例超過(guò)10%,殘差均值比超過(guò)1.5作為損傷判別標(biāo)準(zhǔn)。

    2 協(xié)整分析

    2.1 協(xié)整理論概述

    從時(shí)間序列分析的角度看,協(xié)整就是將一組具有長(zhǎng)期共同趨勢(shì)的時(shí)間序列,通過(guò)線(xiàn)性組合的方式規(guī)整為一個(gè)時(shí)間序列,這個(gè)新的序列在反映了原始序列特征的同時(shí),已經(jīng)消除了共同趨勢(shì)。對(duì)于結(jié)構(gòu)頻率時(shí)間序列(即一系列采集的結(jié)構(gòu)頻率數(shù)據(jù))而言,這種長(zhǎng)期的共同趨勢(shì)通常是由環(huán)境因素變化引起的。

    為了更好地理解協(xié)整的概念,首先需要定義單整過(guò)程。如果一個(gè)非平穩(wěn)隨機(jī)時(shí)間序列{yt,t=1,2,…,N}經(jīng)過(guò)d次差分后為平穩(wěn)過(guò)程,則{yt}為d階單整過(guò)程,記為yt~I(xiàn)(d)?,F(xiàn)假設(shè)一組d階單整的非平穩(wěn)序列yt=[y1t,y2t,…,ynt]T,如果其線(xiàn)性組合

    ζt=β1y1t+β2y2t+…+βnynt

    (11)

    為d-1階單整序列,即為I(d-1)過(guò)程,則yt中的各序列存在協(xié)整關(guān)系。其中β=[β1,β2,…,βn]T稱(chēng)為協(xié)整向量。

    2.2 平穩(wěn)性檢驗(yàn)

    由2.1節(jié)可知,協(xié)整與時(shí)間序列的單整階數(shù),即平穩(wěn)性密切相關(guān),因此在協(xié)整分析前,需要對(duì)時(shí)間序列進(jìn)行平穩(wěn)性判斷,通常用增廣單位根(augmented Dickey-Fuller,ADF)進(jìn)行檢驗(yàn)[20]。構(gòu)造時(shí)間序列{yt}的p階自回歸模型AR(p)的誤差修正形式

    (12)

    式中:Δ為差分算子;ρ和γj為模型系數(shù);εt為高斯白噪聲序列。

    若{yt}是非平穩(wěn)過(guò)程,則模型至少有一單位根,對(duì)應(yīng)ρ=0。故可作以下假設(shè)檢驗(yàn)

    H0∶ρ=0,H1∶ρ<0

    (13)

    檢驗(yàn)統(tǒng)計(jì)量為

    (14)

    若拒絕H0,則{yt}是平穩(wěn)序列,即yt~I(xiàn)(0)。若接受H0,則{yt}是非平穩(wěn)序列,需構(gòu)造原序列差分{Δyt}的自回歸模型,并根據(jù)式(13)進(jìn)行假設(shè)檢驗(yàn)。以此類(lèi)推,如果對(duì)于{Δkyi(t)}的自回歸模型拒絕H0,則yt~I(xiàn)(k)。由于實(shí)際觀測(cè)序列通常不是發(fā)散的,因此多為I(1)或I(2)過(guò)程。

    2.3 Johansen方法

    對(duì)于一個(gè)n維多元時(shí)間序列yt=[y1t,y2t,…,ynt]T,若其各序列同為1階單整,即?yit~I(xiàn)(1),i=1,2,…,n,t=1,2,…,N,則可通過(guò)Johansen方法估計(jì)協(xié)整向量β,其實(shí)質(zhì)是對(duì)yt的向量自回歸模型(vector auto regression,VAR)的參數(shù)作極大似然估計(jì)[21]。

    首先構(gòu)造多元時(shí)間序列yt的p階向量自回歸模型VAR(p)

    (15)

    式中:Πi和Φ為系數(shù)矩陣;εt為多元高斯白噪聲過(guò)程,有εt~N(0,Ω);dt為確定性趨勢(shì)。通過(guò)適當(dāng)?shù)淖儞Q,可以得到模型的誤差修正形式

    (16)

    Johansen方法認(rèn)為,如果n個(gè)時(shí)間序列之間存在協(xié)整關(guān)系,則系數(shù)矩陣Π必不滿(mǎn)秩,設(shè)rank(Π)=m,則原序列yt有m個(gè)協(xié)整向量,且存在n×m的滿(mǎn)秩矩陣α和β,滿(mǎn)足

    Π=αβT

    (17)

    z0t=αβTz1t+Ψz2t+εt

    (18)

    忽略常數(shù)系數(shù)(2π)-nN/2,則似然函數(shù)為

    (19)

    從而得到Ψ的極大似然估計(jì)為

    (20)

    r0t=αβTr1t+εt

    (21)

    (22)

    記m個(gè)特征值λ1≥λ2≥…≥λm及對(duì)應(yīng)的特征向量為φ1,φ2,…,φm,通常取最大的特征值對(duì)應(yīng)的特征向量作為協(xié)整向量。

    2.4 模型定階

    在進(jìn)行Johansen協(xié)整檢驗(yàn)之前,需要先確定模型的滯后階數(shù)p,在正態(tài)分布假設(shè)下,通常可采用AIC,BIC和HQ信息準(zhǔn)則。文獻(xiàn)[22]認(rèn)為,由BIC準(zhǔn)則確定的模型得到的極大似然估計(jì)的效果更好,因此本文采用BIC準(zhǔn)則進(jìn)行定階,其表達(dá)式為

    (23)

    2.5 協(xié)整秩檢驗(yàn)

    通過(guò)求解式(22),我們可以找到協(xié)整向量,但是我們并不能保證各變量確實(shí)存在協(xié)整關(guān)系,為此,Johansen提出了似然比統(tǒng)計(jì)量。

    (24)

    (25)

    且觀察式(19)發(fā)現(xiàn)似然函數(shù)有上界

    (26)

    考慮假設(shè)檢驗(yàn)

    H0∶r=m,H1∶r>m

    (27)

    式中,r為系數(shù)矩陣Π的秩,即為協(xié)整的秩。

    引入似然函數(shù)比Q并結(jié)合式(22)、式(24)~式(26),有

    (28)

    構(gòu)造負(fù)對(duì)數(shù)似然比統(tǒng)計(jì)量LR

    (29)

    其漸近分布為

    LR? tr{f(Wr)[f(Wr)]-1f(Wr)}

    (30)

    式中,tr{·}為矩陣的跡,且有

    (31)

    式中,Wr為n-r維維納過(guò)程。LR統(tǒng)計(jì)量的臨界值與協(xié)整秩r有關(guān),可以通過(guò)數(shù)值模擬得到。檢驗(yàn)的流程按照假設(shè)r=0,1,…,n-1依次進(jìn)行。當(dāng)LRm大于臨界值時(shí),拒絕H0,表明協(xié)整關(guān)系的個(gè)數(shù)大于m,檢驗(yàn)繼續(xù),直至不能拒絕H0為止。若此時(shí)統(tǒng)計(jì)量為L(zhǎng)Rm*,則說(shuō)明原時(shí)間序列存在m*個(gè)協(xié)整關(guān)系。

    2.6 損傷判定

    根據(jù)式(22)求得的結(jié)果,取最大的特征值對(duì)應(yīng)的特征向量作為協(xié)整向量,即β=φ1,則第t個(gè)樣本對(duì)應(yīng)的協(xié)整殘差為

    (32)

    為方便與PCA方法作對(duì)比,對(duì)協(xié)整殘差作中心化處理

    ξt=|ζt-μξ|

    (33)

    式中,μξ為基準(zhǔn)數(shù)據(jù)協(xié)整殘差的均值。

    對(duì)比式(2)和式(22)可以看出:主成分分析的核心算法是求解樣本數(shù)據(jù)協(xié)方差矩陣的特征值問(wèn)題,再取原始數(shù)據(jù)與估計(jì)數(shù)據(jù)的殘差作損傷判定指標(biāo),而協(xié)整分析的核心算法是求解殘差乘積矩的特征值問(wèn)題,直接取協(xié)整殘差作為損傷判定指標(biāo)。因此,兩種方法具有一定的相似性。

    3 穩(wěn)定性評(píng)價(jià)

    基于主成分分析和協(xié)整分析的損傷判定方法本質(zhì)上是一種二元分類(lèi)器,即將結(jié)構(gòu)響應(yīng)數(shù)據(jù)分類(lèi)為“陰性”和“陽(yáng)性”,其中“陰性”代表正常狀態(tài),“陽(yáng)性”代表?yè)p傷狀態(tài)。分類(lèi)器的魯棒性通??捎檬茉囌吖ぷ魈卣髑€(xiàn)(receiver operating characteristic curve,ROC)來(lái)衡量。給定某一閾值水平η,分別記錄損傷狀態(tài)下殘差超出閾值范圍的樣本數(shù)目,即真陽(yáng)性樣本數(shù)目TP;損傷狀態(tài)下位于閾值范圍之內(nèi)的樣本數(shù)目,即假陰性樣本數(shù)目FN;正常狀態(tài)下殘差超出閾值范圍的樣本數(shù)目,即假陽(yáng)性樣本數(shù)目FP;正常狀態(tài)下位于閾值范圍之內(nèi)的樣本數(shù)目,即真陰性樣本數(shù)目TN,最終組成如表1所示的混淆矩陣?;诖?,我們可以計(jì)算該閾值水平下陰性樣本被判定為陽(yáng)性的比例,即假陽(yáng)性率FPR(η);以及陽(yáng)性樣本被判定為陽(yáng)性的比例,即真陽(yáng)性率TPR(η)

    表1 混淆矩陣Tab.1 Confusion matrix

    (34)

    (35)

    使閾值水平η在從0變化到+∞,我們可以得到體現(xiàn)該分類(lèi)器性能的ROC。為了便于描述,通常通過(guò)計(jì)算ROC的面積進(jìn)行說(shuō)明。當(dāng)ROC位于曲線(xiàn)y=x附近,即當(dāng)其所圍面積約為0.5時(shí),表明該分類(lèi)器等同于隨機(jī)試驗(yàn),已不具備判別能力;ROC所圍面積越接近1,表明該分類(lèi)器魯棒性越好。

    然而,直接使用ROC對(duì)分類(lèi)器的性能優(yōu)劣進(jìn)行評(píng)價(jià)是不嚴(yán)謹(jǐn)?shù)?,不同的測(cè)試數(shù)據(jù)將產(chǎn)生不同的ROC,因此通過(guò)方差進(jìn)行度量才能更準(zhǔn)確地反映出分類(lèi)器性能的優(yōu)劣[23]。方差的獲取可通過(guò)對(duì)多個(gè)測(cè)試數(shù)據(jù)集取平均實(shí)現(xiàn),常用的有垂直平均法和閾值平均法。

    3.1 垂直平均法

    垂直平均法通過(guò)給定假陽(yáng)性率FPR的采樣間隔,對(duì)多條ROC在對(duì)應(yīng)的假陽(yáng)性率FPR采樣點(diǎn)處的真陽(yáng)性率TPR取平均,并記錄相應(yīng)的標(biāo)準(zhǔn)差。對(duì)于某一特定的FPR采樣點(diǎn),當(dāng)ROC上不存在對(duì)應(yīng)的真陽(yáng)性率TPR時(shí),可通過(guò)相鄰的兩個(gè)假陽(yáng)性率FPR對(duì)應(yīng)的TPR插值獲取。

    3.2 閾值平均法

    垂直平均法的優(yōu)點(diǎn)在于其平均值由因變量TPR組成,這簡(jiǎn)化了置信區(qū)間的計(jì)算。然而,ROC的自變量FPR往往不是人為控制的。而閾值平均法使獨(dú)立變量人工可控。該方法根據(jù)產(chǎn)生ROC上的點(diǎn)的閾值進(jìn)行采樣,然后為每個(gè)閾值找到每條ROC的對(duì)應(yīng)點(diǎn),再分別關(guān)于對(duì)應(yīng)ROC點(diǎn)的FPR和TPR取平均。因此,閾值平均法將產(chǎn)生垂直和水平兩個(gè)方向的方差。

    4 數(shù)值模擬

    4.1 海洋平臺(tái)模型

    本文以一個(gè)導(dǎo)管架平臺(tái)模型為例,比較主成分分析與協(xié)整分析的性能差異。主要考慮了環(huán)境溫度對(duì)結(jié)構(gòu)頻率的影響。結(jié)構(gòu)頻率通過(guò)MATLAB軟件求取結(jié)構(gòu)的特征方程獲得。

    導(dǎo)管架平臺(tái)有限元模型簡(jiǎn)圖,如圖1所示。此模型劃分為80個(gè)單元,48個(gè)節(jié)點(diǎn)。結(jié)構(gòu)的材料密度ρ=7 850 kg/m3,泊松比μ=0.3。平臺(tái)總質(zhì)量為550 t,其中上部荷載為250 t,為簡(jiǎn)化模擬,將上部載荷簡(jiǎn)化為質(zhì)量單元平均分配到平臺(tái)頂層的4個(gè)節(jié)點(diǎn)上。平臺(tái)頂層和底層尺寸分別為6.097 m×5.284 m和10.112 m×8.764 m。平臺(tái)各桿件尺寸,如表2所示。該平臺(tái)作業(yè)海域水深22 m,取泥面以下9 m為簡(jiǎn)化固定端。

    圖1 導(dǎo)管架平臺(tái)模型(m)Fig.1 Sketch of the offshore platform structure(m)

    表2 平臺(tái)各桿件尺寸Tab.2 Member dimensions of the offshore platform

    4.2 溫度對(duì)結(jié)構(gòu)頻率的影響

    溫度對(duì)結(jié)構(gòu)頻率的影響是通過(guò)改變結(jié)構(gòu)材料的彈性模量實(shí)現(xiàn)的。根據(jù)Woon等[24]的研究結(jié)果,鋼材彈性模量與環(huán)境溫度之間存在以下線(xiàn)性關(guān)系

    En(κν)=(1+cv)En(κ0)

    (36)

    式中,En(κ0)和En(κν)分別為參考環(huán)境溫度κ0和當(dāng)前環(huán)境溫度κν下單元n的彈性模量。本文取參考環(huán)境溫度κ0為20 ℃,對(duì)應(yīng)的材料彈性模量為206 GPa。變化系數(shù)取為cv=c0(κν-κ0)/En(κ0),c0=-1.0×108N/m2·°C。環(huán)境溫度與鋼材彈性模量的關(guān)系,如圖2所示。

    圖2 環(huán)境溫度與鋼材彈性模量的關(guān)系Fig.2 Relationship between steel elastic modulus and ambient temperature

    4.3 工況分析

    海洋平臺(tái)結(jié)構(gòu)位于空氣、海水和海底泥土3種介質(zhì)中,3種介質(zhì)的傳熱性質(zhì)差異較大,其溫度變化不盡相同。本文采用了某海洋觀測(cè)站實(shí)測(cè)的溫度變化數(shù)據(jù),數(shù)據(jù)記錄了該海域每小時(shí)的氣溫、水溫和海底以下5 m的泥溫,截取了4 000組溫度變化的數(shù)據(jù)以模擬結(jié)構(gòu)頻率的變化,3種環(huán)境溫度的變化曲線(xiàn),如圖3所示。其中,1~3 500組溫度用于模擬正常狀態(tài)下結(jié)構(gòu)的頻率變化,作為基準(zhǔn)數(shù)據(jù);3 501~4 000組溫度用于模擬各損傷工況下結(jié)構(gòu)的頻率變化,作為測(cè)試數(shù)據(jù)。損傷工況設(shè)置,如表3所示,其中工況A為健康工況,用于探究方法是否會(huì)發(fā)生損傷誤判,工況B~工況E為損傷工況,分別模擬平臺(tái)內(nèi)水平撐單元、外水平撐單元、立面斜撐單元以及立柱單元的損傷,用于探究方法是否會(huì)發(fā)生損傷漏判,結(jié)構(gòu)損傷通過(guò)對(duì)相應(yīng)單元進(jìn)行剛度削減來(lái)模擬。圖4給出了結(jié)構(gòu)前3階頻率的模擬結(jié)果,其中1~3 500組作為基準(zhǔn)數(shù)據(jù),3 501~6 000組以每500組為一個(gè)工況,分別對(duì)應(yīng)工況A~工況E。可以看出,溫度變化在一定程度上掩蓋了結(jié)構(gòu)損傷導(dǎo)致的結(jié)構(gòu)頻率變化,為了避免損傷誤判和損傷漏判的發(fā)生,有必要消除溫度變化的影響。

    圖3 實(shí)測(cè)溫度變化曲線(xiàn)Fig.3 Measured temperature history series

    表3 模擬的海洋平臺(tái)損傷工況Tab.3 Simulated damage cases of the offshore platform

    圖4 結(jié)構(gòu)前3階頻率Fig.4 First three natural frequencies under temperature changing

    提取前3階頻率的3 500組基準(zhǔn)數(shù)據(jù),分別進(jìn)行主成分分析和協(xié)整分析,所得主成分轉(zhuǎn)換矩陣T和協(xié)整向量β將用于各損傷工況的判定。其中,主成分分析的能量閾值取為0.95,由表4可知選取的主成分階數(shù)為1。協(xié)整秩檢驗(yàn)的顯著性水平取為0.01,參見(jiàn)表5。首先假設(shè)r=0,根據(jù)式(29)計(jì)算得到LR統(tǒng)計(jì)量為75.641 0>41.072 2,拒絕原假設(shè)。進(jìn)一步假設(shè)r=1,此時(shí)統(tǒng)計(jì)量為16.075 5<23.157 4,接受原假設(shè),因此頻率數(shù)據(jù)存在一個(gè)協(xié)整關(guān)系。可見(jiàn),溫度影響這一“共同趨勢(shì)”起了主要作用。

    表4 主成分階數(shù)的選取Tab.4 Selection number of principle components

    表5 協(xié)整秩檢驗(yàn)Tab.5 Hypothesis test of cointegration rank

    為方便比較,將主成分殘差dt與協(xié)整殘差ξt均作最大值歸一化處理。由圖5及表6、表7可以看出:對(duì)于工況A,兩種算法的異常值比例均遠(yuǎn)小于10%,殘差均值比均小于1.5,其判定結(jié)果為結(jié)構(gòu)沒(méi)有發(fā)生損傷,與實(shí)際情況一致。對(duì)于工況B,兩種算法的異常值比例均小于10%,殘差均值比均小于1.5,未能準(zhǔn)確識(shí)別出損傷,主要原因是單元71為內(nèi)水平撐,其剛度損失對(duì)結(jié)構(gòu)整體剛度的變化幾乎沒(méi)有影響,以至于損傷工況B對(duì)結(jié)構(gòu)頻率的變化影響甚微,因此導(dǎo)致了兩種算法均出現(xiàn)了漏判。對(duì)于工況C~工況E,兩種算法的異常值比例均為100%,殘差均值比遠(yuǎn)大于1.5,完全準(zhǔn)確地識(shí)別出結(jié)構(gòu)損傷。對(duì)比來(lái)看,工況A下兩種算法的殘差均值比均在1.1附近,差別并不明顯。對(duì)于工況C~工況E,協(xié)整分析對(duì)應(yīng)的殘差均值比明顯較高,說(shuō)明該算法在察覺(jué)結(jié)構(gòu)頻率的異常數(shù)據(jù)方面具有更高的靈敏度。

    圖5 各損傷工況判定結(jié)果Fig.5 The damage detection results of damage case A to damage case E

    表6 各工況異常值比例Tab.6 Outliners ratio in each case %

    表7 各工況殘差均值比Tab.7 Mean value ratio of residuals

    由于工況B~工況E對(duì)應(yīng)不同的桿件類(lèi)型發(fā)生相同程度的損傷,因此其殘差均值比的大小代表了兩種算法對(duì)不同桿件類(lèi)型損傷的靈敏度。由表7可以看出,對(duì)于主成分分析,不同桿件類(lèi)型的損傷檢測(cè)靈敏度從高到低依次為:立面斜撐>立柱>外水平撐>內(nèi)水平撐;對(duì)于協(xié)整分析,不同桿件類(lèi)型的損傷檢測(cè)靈敏度從高到低依次為:立面斜撐>外水平撐>立柱>內(nèi)水平撐??梢?jiàn),兩種算法均對(duì)立面斜撐的損傷最為靈敏。對(duì)比立柱和外水平撐兩種桿件的損傷,若采用主成分分析,則立柱的損傷更容易識(shí)別;若采用協(xié)整分析,則外水平撐的損傷更容易識(shí)別。因此,實(shí)際檢測(cè)過(guò)程中可根據(jù)目標(biāo)檢測(cè)桿件的不同采用不同的方法。

    4.4 算法的噪聲魯棒性

    實(shí)際頻率識(shí)別過(guò)程中還不可避免的受到系統(tǒng)和量測(cè)噪聲等不確定性因素的干擾,參考梁亞斌等的做法,本文采用以下噪聲模型

    (37)

    仍以結(jié)構(gòu)前3階頻率作為樣本數(shù)據(jù),設(shè)定某一噪聲水平,通過(guò)一次模擬,得到噪聲影響下的一組頻率樣本,通過(guò)改變X-bar圖的閾值控制線(xiàn),得到該頻率樣本下的ROC(PCA和CA各對(duì)應(yīng)一條ROC)。將以上模擬重復(fù)1 000次,得到1 000條ROC,對(duì)曲線(xiàn)取垂直平均,則可得到某一噪聲水平下的平均ROC。以下取工況A的頻率樣本作為陰性樣本,工況C的頻率樣本作為陽(yáng)性樣本,測(cè)試兩種算法的噪聲魯棒性。其中假陽(yáng)性率FPR的采樣間隔取為0.002。由圖6和表8可以看出,當(dāng)噪聲水平為5%時(shí),協(xié)整分析對(duì)應(yīng)的平均ROC面積更接近于1,且其1 000次模擬得到的ROC面積的標(biāo)準(zhǔn)差相對(duì)更小,說(shuō)明其魯棒性更好。但隨著噪聲水平的增加,主成分分析對(duì)應(yīng)的平均ROC面積更接近于1,且其ROC面積的標(biāo)準(zhǔn)差相對(duì)更小,因此其魯棒性更好。

    圖6 工況C平均ROCFig.6 The averaging ROC curves of case C

    表8 工況D平均ROC所圍面積及標(biāo)準(zhǔn)差Tab.8 Standard deviation of the area under ROC

    4.5 基準(zhǔn)樣本數(shù)目的影響

    由于主成分分析的投影矩陣T和協(xié)整分析的協(xié)整向量β是通過(guò)基準(zhǔn)數(shù)據(jù)得到的,因此基準(zhǔn)數(shù)據(jù)的樣本數(shù)目對(duì)算法的性能會(huì)造成一定的影響。取結(jié)構(gòu)前3階頻率作為樣本數(shù)據(jù),設(shè)定某一噪聲水平,改變基準(zhǔn)數(shù)據(jù)的樣本數(shù)目N,求解該基準(zhǔn)樣本數(shù)目對(duì)應(yīng)的投影矩陣T和協(xié)整向量β以用于損傷工況的判定。仍取工況A的頻率樣本作為陰性樣本,工況C的頻率樣本作為陽(yáng)性樣本,以1.2節(jié)定義的上控制線(xiàn)UCL作為閾值控制線(xiàn),記錄相應(yīng)的假陽(yáng)性率FPR。由于假陽(yáng)性率代表了損傷誤判的概率,因此應(yīng)盡可能的小??紤]到噪聲的隨機(jī)性,將以上模擬重復(fù)1 000次,并對(duì)相應(yīng)的FPR取平均,則可得到FPR隨基準(zhǔn)樣本數(shù)目N變化的曲線(xiàn)。由圖8可以看出,對(duì)于噪聲水平為10%的情況,當(dāng)基準(zhǔn)樣本數(shù)目為500~1 200時(shí),主成分分析對(duì)應(yīng)的假陽(yáng)性率大于0.96,誤判較為嚴(yán)重,而協(xié)整分析對(duì)應(yīng)的假陽(yáng)性率小于0.21,誤判率較低。對(duì)于噪聲水平為15%的情況,當(dāng)基準(zhǔn)樣本數(shù)目為500~1 500時(shí)也有類(lèi)似的現(xiàn)象。以上說(shuō)明了在一定的噪聲水平下,當(dāng)基準(zhǔn)樣本數(shù)目較小時(shí),協(xié)整分析可以更加有效地抑制損傷誤判的問(wèn)題。

    圖7 基準(zhǔn)頻率的選取(以第1階頻率為例)Fig.7 Baseline natural frequencies selection(the 1st frequency as an example)

    圖8 假陽(yáng)性率隨基準(zhǔn)樣本數(shù)目的變化曲線(xiàn)Fig.8 Curves of FPR changing with sample number of baseline natural frequencies

    5 結(jié) 論

    本文以一個(gè)海洋平臺(tái)模型為例,考慮了空氣、海水和海底泥土溫度變化的影響。以結(jié)構(gòu)的模態(tài)頻率作為樣本數(shù)據(jù),分別利用主成分分析和協(xié)整分析消除溫度變化影響,提取能夠表達(dá)結(jié)構(gòu)真實(shí)特性的主成分殘差和協(xié)整殘差,繼而根據(jù)X-bar控制圖并結(jié)合結(jié)構(gòu)損傷前后的殘差均值比進(jìn)行損傷判定。通過(guò)改變頻率數(shù)據(jù)的噪聲水平以及基準(zhǔn)樣本的數(shù)目,比較分析了兩種方法的性能差異。結(jié)果表明:

    (1)兩種方法均能有效消除環(huán)境溫度的影響,避免損傷誤判和損傷漏判的發(fā)生。在不考慮噪聲的條件下,各損傷工況對(duì)應(yīng)的協(xié)整殘差均值比更大,說(shuō)明協(xié)整分析具有更高的損傷檢測(cè)靈敏度。

    (2)小噪聲水平下(噪聲水平小于5%),協(xié)整分析對(duì)應(yīng)的平均ROC所圍面積更接近1,且其ROC面積的標(biāo)準(zhǔn)差相對(duì)更小,因此其魯棒性更好。隨著噪聲水平的增加,主成分分析的噪聲魯棒性更好。

    (3)在一定的噪聲水平下,當(dāng)基準(zhǔn)頻率樣本數(shù)目較小時(shí),主成分分析對(duì)應(yīng)的假陽(yáng)性率較高,誤判嚴(yán)重;協(xié)整分析對(duì)應(yīng)的假陽(yáng)性率較低,可更加有效地抑制損傷誤判的問(wèn)題。

    猜你喜歡
    協(xié)整殘差閾值
    基于雙向GRU與殘差擬合的車(chē)輛跟馳建模
    基于殘差學(xué)習(xí)的自適應(yīng)無(wú)人機(jī)目標(biāo)跟蹤算法
    小波閾值去噪在深小孔鉆削聲發(fā)射信號(hào)處理中的應(yīng)用
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    外商直接投資對(duì)我國(guó)進(jìn)出口貿(mào)易影響的協(xié)整分析
    基于自適應(yīng)閾值和連通域的隧道裂縫提取
    河南金融發(fā)展和城鄉(xiāng)居民收入差距的協(xié)整分析
    比值遙感蝕變信息提取及閾值確定(插圖)
    河北遙感(2017年2期)2017-08-07 14:49:00
    室內(nèi)表面平均氡析出率閾值探討
    平穩(wěn)自相關(guān)過(guò)程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    中文字幕久久专区| 亚洲人与动物交配视频| 琪琪午夜伦伦电影理论片6080| 国产午夜精品久久久久久一区二区三区 | av黄色大香蕉| 欧美丝袜亚洲另类 | 在线a可以看的网站| 亚洲国产欧美网| 亚洲人成网站在线播| 高清毛片免费观看视频网站| 国产精品三级大全| 久久99热这里只有精品18| 麻豆久久精品国产亚洲av| 三级男女做爰猛烈吃奶摸视频| 日日夜夜操网爽| 亚洲在线自拍视频| 成年女人毛片免费观看观看9| 一级作爱视频免费观看| 69人妻影院| 搡女人真爽免费视频火全软件 | 99精品久久久久人妻精品| 国产高清有码在线观看视频| 国产黄a三级三级三级人| 黄色女人牲交| 精品人妻1区二区| 老汉色∧v一级毛片| 99视频精品全部免费 在线| 欧美乱色亚洲激情| 午夜免费男女啪啪视频观看 | 校园春色视频在线观看| 久久久久性生活片| 午夜福利在线观看吧| 欧美成人性av电影在线观看| 国产一区二区三区在线臀色熟女| 哪里可以看免费的av片| 97超视频在线观看视频| 国产男靠女视频免费网站| 又粗又爽又猛毛片免费看| 激情在线观看视频在线高清| 一个人看视频在线观看www免费 | 91麻豆精品激情在线观看国产| 一a级毛片在线观看| 久久国产乱子伦精品免费另类| 伊人久久大香线蕉亚洲五| 一区二区三区高清视频在线| 国产精品久久久久久久电影 | 国产欧美日韩精品亚洲av| 在线免费观看的www视频| 男人舔女人下体高潮全视频| 99国产极品粉嫩在线观看| 欧美三级亚洲精品| 一个人看的www免费观看视频| 色老头精品视频在线观看| 伊人久久精品亚洲午夜| 日日摸夜夜添夜夜添小说| 法律面前人人平等表现在哪些方面| 蜜桃亚洲精品一区二区三区| 美女 人体艺术 gogo| 一级作爱视频免费观看| 欧美色视频一区免费| 一级作爱视频免费观看| 熟女人妻精品中文字幕| 欧美乱妇无乱码| 毛片女人毛片| 麻豆国产av国片精品| 久久久国产成人免费| 亚洲最大成人手机在线| 午夜老司机福利剧场| 一个人观看的视频www高清免费观看| 动漫黄色视频在线观看| 欧美最新免费一区二区三区 | 国产成人欧美在线观看| 桃红色精品国产亚洲av| 18禁裸乳无遮挡免费网站照片| 亚洲五月天丁香| 亚洲五月天丁香| 久久亚洲真实| 国产伦人伦偷精品视频| 午夜福利免费观看在线| 久久久久久久久久黄片| 中文字幕人妻丝袜一区二区| 波多野结衣高清作品| 亚洲精品美女久久久久99蜜臀| 九九久久精品国产亚洲av麻豆| 亚洲美女视频黄频| 国产主播在线观看一区二区| 黄色女人牲交| 99在线人妻在线中文字幕| a级毛片a级免费在线| 久久精品影院6| 美女黄网站色视频| 高清日韩中文字幕在线| 国产精品香港三级国产av潘金莲| 午夜影院日韩av| 国产爱豆传媒在线观看| 91字幕亚洲| 亚洲专区中文字幕在线| 老熟妇仑乱视频hdxx| 国产免费一级a男人的天堂| 国产精品综合久久久久久久免费| 男人舔奶头视频| 精品久久久久久,| avwww免费| 久久婷婷人人爽人人干人人爱| 婷婷六月久久综合丁香| 美女黄网站色视频| 天堂√8在线中文| 国产精品日韩av在线免费观看| 19禁男女啪啪无遮挡网站| 精品久久久久久久毛片微露脸| 日韩欧美免费精品| 69av精品久久久久久| 18美女黄网站色大片免费观看| 久久久久国内视频| 久久午夜亚洲精品久久| 国产真实乱freesex| 51国产日韩欧美| 日本黄色片子视频| 色av中文字幕| 级片在线观看| 看免费av毛片| 观看免费一级毛片| 国内精品久久久久精免费| 国产激情偷乱视频一区二区| 人人妻人人看人人澡| 天堂av国产一区二区熟女人妻| 亚洲国产日韩欧美精品在线观看 | 精品国产美女av久久久久小说| 色吧在线观看| 有码 亚洲区| 嫩草影院精品99| 日本a在线网址| av视频在线观看入口| 午夜免费激情av| 国产欧美日韩一区二区精品| 欧美色视频一区免费| 精品无人区乱码1区二区| 99久国产av精品| 亚洲天堂国产精品一区在线| 啦啦啦韩国在线观看视频| 免费在线观看日本一区| 国产精品野战在线观看| 欧美精品啪啪一区二区三区| 在线观看免费视频日本深夜| 此物有八面人人有两片| 五月伊人婷婷丁香| 男女做爰动态图高潮gif福利片| 岛国视频午夜一区免费看| 中文字幕熟女人妻在线| 欧美成人性av电影在线观看| 日韩高清综合在线| 成人性生交大片免费视频hd| 可以在线观看的亚洲视频| 日本黄色片子视频| 国内久久婷婷六月综合欲色啪| 国产精品野战在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 内射极品少妇av片p| 黄色成人免费大全| 精品国产超薄肉色丝袜足j| 在线观看av片永久免费下载| 欧美日本亚洲视频在线播放| 久久婷婷人人爽人人干人人爱| 亚洲国产欧美人成| 午夜亚洲福利在线播放| 久久久久九九精品影院| 国产精品美女特级片免费视频播放器| 蜜桃亚洲精品一区二区三区| 国产又黄又爽又无遮挡在线| 看片在线看免费视频| 国产毛片a区久久久久| 久久精品91蜜桃| 天堂av国产一区二区熟女人妻| 中文资源天堂在线| 国产亚洲精品综合一区在线观看| 久久久久国产精品人妻aⅴ院| av福利片在线观看| 久久久久亚洲av毛片大全| 男女下面进入的视频免费午夜| 深夜精品福利| 两个人视频免费观看高清| 国产亚洲av嫩草精品影院| 国产成人欧美在线观看| 欧美大码av| 欧美乱色亚洲激情| 日本黄大片高清| 最好的美女福利视频网| 久久中文看片网| 18禁在线播放成人免费| 久久天躁狠狠躁夜夜2o2o| 国产成人av教育| 亚洲天堂国产精品一区在线| 午夜免费男女啪啪视频观看 | 久久精品亚洲精品国产色婷小说| 熟女电影av网| 亚洲精品在线观看二区| 欧美bdsm另类| 女生性感内裤真人,穿戴方法视频| 男女之事视频高清在线观看| 最近最新免费中文字幕在线| 51国产日韩欧美| 高清在线国产一区| 久久久久免费精品人妻一区二区| 成人特级黄色片久久久久久久| 免费看a级黄色片| 丰满人妻一区二区三区视频av | 亚洲最大成人手机在线| 在线看三级毛片| 亚洲在线自拍视频| 婷婷亚洲欧美| 老汉色∧v一级毛片| 一区福利在线观看| 久久久久久大精品| 中文字幕熟女人妻在线| 国产亚洲av嫩草精品影院| 久久精品国产99精品国产亚洲性色| 亚洲第一欧美日韩一区二区三区| 日本一本二区三区精品| 在线播放无遮挡| 国产老妇女一区| 国产真人三级小视频在线观看| 九色成人免费人妻av| 最近最新免费中文字幕在线| 亚洲成a人片在线一区二区| 99久久成人亚洲精品观看| 久久久久久久久中文| 熟女电影av网| av国产免费在线观看| 亚洲熟妇熟女久久| 久久亚洲真实| 国产精品亚洲av一区麻豆| 欧美黄色淫秽网站| 美女被艹到高潮喷水动态| 成年版毛片免费区| av中文乱码字幕在线| 国产精品免费一区二区三区在线| 我要搜黄色片| 国产麻豆成人av免费视频| 中文字幕人妻熟人妻熟丝袜美 | 听说在线观看完整版免费高清| a在线观看视频网站| 久久草成人影院| 国产亚洲精品综合一区在线观看| 日韩中文字幕欧美一区二区| 日韩成人在线观看一区二区三区| 久久香蕉国产精品| www.www免费av| 亚洲天堂国产精品一区在线| 久久精品国产亚洲av涩爱 | 好男人在线观看高清免费视频| 成人亚洲精品av一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 成人特级黄色片久久久久久久| 精品国产亚洲在线| 午夜a级毛片| 淫秽高清视频在线观看| 国产精品久久视频播放| 久久6这里有精品| 日本黄色片子视频| 亚洲人成伊人成综合网2020| 亚洲天堂国产精品一区在线| 亚洲国产精品合色在线| 欧美日本亚洲视频在线播放| 亚洲精品国产精品久久久不卡| 十八禁人妻一区二区| 欧美日韩亚洲国产一区二区在线观看| 波多野结衣巨乳人妻| 国产高潮美女av| 婷婷六月久久综合丁香| 啪啪无遮挡十八禁网站| 夜夜夜夜夜久久久久| 国产伦精品一区二区三区四那| 一区二区三区激情视频| 国产精品av视频在线免费观看| 最近最新中文字幕大全免费视频| 99热6这里只有精品| 夜夜躁狠狠躁天天躁| 欧美不卡视频在线免费观看| 精品人妻偷拍中文字幕| 亚洲av美国av| 国模一区二区三区四区视频| 久久天躁狠狠躁夜夜2o2o| 51国产日韩欧美| 一级毛片女人18水好多| 欧美日韩综合久久久久久 | 亚洲性夜色夜夜综合| 成人无遮挡网站| 性欧美人与动物交配| 国产探花极品一区二区| 国产一区二区三区视频了| 精品乱码久久久久久99久播| 日日干狠狠操夜夜爽| 女人被狂操c到高潮| 18+在线观看网站| 精品国产超薄肉色丝袜足j| 色综合亚洲欧美另类图片| 51午夜福利影视在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品电影一区二区三区| 999久久久精品免费观看国产| 婷婷六月久久综合丁香| 母亲3免费完整高清在线观看| www国产在线视频色| 18+在线观看网站| 成人国产一区最新在线观看| 麻豆国产97在线/欧美| 好男人在线观看高清免费视频| 国产 一区 欧美 日韩| 在线观看免费视频日本深夜| 日本 av在线| 欧美成人a在线观看| 亚洲久久久久久中文字幕| 白带黄色成豆腐渣| 午夜免费成人在线视频| 久久久久久久午夜电影| 国产成人av教育| 国产激情偷乱视频一区二区| 久久久久久九九精品二区国产| 99热只有精品国产| 亚洲成av人片免费观看| 精品久久久久久久久久免费视频| 国产一级毛片七仙女欲春2| 美女高潮喷水抽搐中文字幕| 国产亚洲av嫩草精品影院| 观看免费一级毛片| 久久国产精品影院| 久久久久久久久大av| 十八禁人妻一区二区| 最近最新中文字幕大全电影3| 亚洲成人免费电影在线观看| 宅男免费午夜| 男人的好看免费观看在线视频| 国产在线精品亚洲第一网站| 女同久久另类99精品国产91| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久久精品吃奶| 日韩欧美在线二视频| 国语自产精品视频在线第100页| 日韩欧美一区二区三区在线观看| 亚洲黑人精品在线| 国产又黄又爽又无遮挡在线| 动漫黄色视频在线观看| 亚洲18禁久久av| 成人性生交大片免费视频hd| 天堂av国产一区二区熟女人妻| 久久欧美精品欧美久久欧美| 亚洲性夜色夜夜综合| 欧美3d第一页| 岛国在线观看网站| 日本 av在线| www.www免费av| 国产色爽女视频免费观看| 熟女电影av网| 中文字幕人妻丝袜一区二区| 99精品久久久久人妻精品| 一级毛片高清免费大全| 女人高潮潮喷娇喘18禁视频| 天美传媒精品一区二区| 亚洲av电影不卡..在线观看| 国语自产精品视频在线第100页| 久久久成人免费电影| 日韩高清综合在线| 啦啦啦观看免费观看视频高清| 欧美+亚洲+日韩+国产| 欧美性猛交╳xxx乱大交人| 久久草成人影院| 免费无遮挡裸体视频| 黄色女人牲交| 久久午夜亚洲精品久久| 午夜视频国产福利| 亚洲专区中文字幕在线| 国产精品自产拍在线观看55亚洲| 中国美女看黄片| 在线播放国产精品三级| 在线天堂最新版资源| 美女黄网站色视频| 成人性生交大片免费视频hd| 一区二区三区国产精品乱码| 女人高潮潮喷娇喘18禁视频| 日日干狠狠操夜夜爽| 国产综合懂色| 人妻夜夜爽99麻豆av| 亚洲av电影在线进入| 色综合欧美亚洲国产小说| 欧美日韩综合久久久久久 | www.www免费av| 看免费av毛片| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 色尼玛亚洲综合影院| 一区二区三区激情视频| a级毛片a级免费在线| 精品午夜福利视频在线观看一区| 一级黄色大片毛片| 国产一区二区亚洲精品在线观看| 国产精品亚洲av一区麻豆| 免费av不卡在线播放| 哪里可以看免费的av片| x7x7x7水蜜桃| 国产精品乱码一区二三区的特点| 亚洲美女视频黄频| 给我免费播放毛片高清在线观看| 99久久精品一区二区三区| 97碰自拍视频| 久久亚洲真实| 免费人成在线观看视频色| av天堂在线播放| 母亲3免费完整高清在线观看| 亚洲av成人精品一区久久| 国产高清视频在线播放一区| 少妇裸体淫交视频免费看高清| 在线播放国产精品三级| 国产爱豆传媒在线观看| 3wmmmm亚洲av在线观看| 国产亚洲精品久久久com| 欧美绝顶高潮抽搐喷水| 久久久久久久久大av| 中文字幕精品亚洲无线码一区| 成人高潮视频无遮挡免费网站| 91av网一区二区| 午夜福利视频1000在线观看| 波多野结衣巨乳人妻| 男女下面进入的视频免费午夜| 好男人在线观看高清免费视频| 亚洲国产日韩欧美精品在线观看 | 两人在一起打扑克的视频| 白带黄色成豆腐渣| 老汉色av国产亚洲站长工具| 色综合欧美亚洲国产小说| 国产一级毛片七仙女欲春2| 日本免费a在线| 国产aⅴ精品一区二区三区波| 亚洲av成人不卡在线观看播放网| 亚洲国产精品成人综合色| 国产欧美日韩精品一区二区| a级一级毛片免费在线观看| 999久久久精品免费观看国产| 亚洲av电影在线进入| 成人18禁在线播放| 午夜福利高清视频| 亚洲av第一区精品v没综合| 99热这里只有是精品50| 国产av在哪里看| 性色avwww在线观看| 国产亚洲精品一区二区www| 宅男免费午夜| 国产日本99.免费观看| 国产精品98久久久久久宅男小说| 亚洲精品国产精品久久久不卡| 国产单亲对白刺激| 国产99白浆流出| а√天堂www在线а√下载| 国产高清三级在线| 亚洲乱码一区二区免费版| 高清毛片免费观看视频网站| 国产伦在线观看视频一区| 99热只有精品国产| 免费看日本二区| 人妻夜夜爽99麻豆av| 亚洲国产精品成人综合色| 美女被艹到高潮喷水动态| 国产精品久久久人人做人人爽| 熟女人妻精品中文字幕| 国产精品一区二区三区四区免费观看 | 少妇的逼水好多| 欧美午夜高清在线| 亚洲一区高清亚洲精品| 免费av不卡在线播放| 香蕉丝袜av| 淫妇啪啪啪对白视频| 国产一区二区在线观看日韩 | 尤物成人国产欧美一区二区三区| 国内少妇人妻偷人精品xxx网站| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品综合一区在线观看| 亚洲激情在线av| 狂野欧美激情性xxxx| 老汉色∧v一级毛片| 男女视频在线观看网站免费| 一级毛片高清免费大全| 久久中文看片网| 久久草成人影院| 中文字幕av成人在线电影| 欧美在线一区亚洲| 91九色精品人成在线观看| 丰满的人妻完整版| 国产单亲对白刺激| 99久久久亚洲精品蜜臀av| 成人无遮挡网站| 国产成+人综合+亚洲专区| 亚洲内射少妇av| 欧美大码av| 午夜福利欧美成人| 少妇人妻精品综合一区二区 | 久久中文看片网| 婷婷丁香在线五月| 在线观看av片永久免费下载| 国产精华一区二区三区| 国产乱人伦免费视频| 国产亚洲欧美在线一区二区| 午夜久久久久精精品| 在线视频色国产色| 一级作爱视频免费观看| 一个人看的www免费观看视频| 天美传媒精品一区二区| 天堂网av新在线| 黄色女人牲交| 首页视频小说图片口味搜索| 伊人久久大香线蕉亚洲五| 性色avwww在线观看| 国产毛片a区久久久久| 亚洲无线在线观看| 高潮久久久久久久久久久不卡| 国产主播在线观看一区二区| 夜夜夜夜夜久久久久| 欧美一区二区精品小视频在线| 性色avwww在线观看| 成人国产综合亚洲| 日韩欧美免费精品| 人人妻人人看人人澡| 国产精品一区二区三区四区久久| 亚洲人成网站在线播| 欧美日本亚洲视频在线播放| 午夜精品久久久久久毛片777| 国语自产精品视频在线第100页| 我要搜黄色片| 在线国产一区二区在线| 18禁国产床啪视频网站| 色综合欧美亚洲国产小说| 亚洲精品在线观看二区| 久99久视频精品免费| 精品国产超薄肉色丝袜足j| 99久久无色码亚洲精品果冻| 天天一区二区日本电影三级| 90打野战视频偷拍视频| 中文在线观看免费www的网站| 一级a爱片免费观看的视频| 老司机午夜十八禁免费视频| 老汉色∧v一级毛片| 最近最新免费中文字幕在线| 少妇高潮的动态图| 在线观看美女被高潮喷水网站 | 精品欧美国产一区二区三| 99视频精品全部免费 在线| 欧美zozozo另类| 国产精品亚洲美女久久久| 乱人视频在线观看| 麻豆一二三区av精品| 久久中文看片网| 天堂影院成人在线观看| 欧美bdsm另类| 少妇的丰满在线观看| 精品电影一区二区在线| 给我免费播放毛片高清在线观看| 免费av不卡在线播放| 国产精品久久久久久人妻精品电影| 九色国产91popny在线| 最新中文字幕久久久久| 欧美中文日本在线观看视频| 此物有八面人人有两片| 18禁黄网站禁片午夜丰满| 亚洲专区中文字幕在线| 国产欧美日韩一区二区精品| 91字幕亚洲| x7x7x7水蜜桃| 国产三级中文精品| 亚洲人与动物交配视频| 动漫黄色视频在线观看| 18禁裸乳无遮挡免费网站照片| av天堂中文字幕网| 久久久久国内视频| 亚洲人成网站在线播放欧美日韩| 69av精品久久久久久| 国产免费男女视频| 少妇裸体淫交视频免费看高清| 国产伦在线观看视频一区| 欧美性感艳星| 国产亚洲精品综合一区在线观看| 欧美日韩中文字幕国产精品一区二区三区| 黑人欧美特级aaaaaa片| 日韩欧美一区二区三区在线观看| 精品99又大又爽又粗少妇毛片 | 国产视频内射| 日本五十路高清| 亚洲人成伊人成综合网2020| 波多野结衣高清无吗| 欧美日韩一级在线毛片| 日韩免费av在线播放| 热99re8久久精品国产| 日韩免费av在线播放| 热99re8久久精品国产| 此物有八面人人有两片| 非洲黑人性xxxx精品又粗又长| 精品电影一区二区在线| 别揉我奶头~嗯~啊~动态视频| 国产高清视频在线播放一区| 成人av在线播放网站| 亚洲专区中文字幕在线| 美女 人体艺术 gogo| 日本五十路高清| 性欧美人与动物交配| 啦啦啦免费观看视频1| 内地一区二区视频在线| 国产精品av视频在线免费观看| 亚洲精品乱码久久久v下载方式 | 首页视频小说图片口味搜索| 婷婷精品国产亚洲av在线| 欧美丝袜亚洲另类 | 在线观看日韩欧美| 久久精品人妻少妇| 99久久精品国产亚洲精品| 在线视频色国产色| 成人亚洲精品av一区二区| 亚洲美女视频黄频| 亚洲天堂国产精品一区在线|