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

    基于貝葉斯優(yōu)化BiLSTM模型的輸電塔損傷識(shí)別

    2023-02-01 07:08:32魏佳恒郭惠勇
    振動(dòng)與沖擊 2023年1期
    關(guān)鍵詞:貝葉斯準(zhǔn)確率有限元

    魏佳恒,郭惠勇

    (重慶大學(xué) 土木工程學(xué)院,重慶 400045)

    特高壓輸電塔穿越高海拔地區(qū),所處環(huán)境復(fù)雜、惡劣,由于其結(jié)構(gòu)高聳、剛度弱的特點(diǎn),在長(zhǎng)期服役期間易產(chǎn)生結(jié)構(gòu)損傷[1]。因此對(duì)輸電塔結(jié)構(gòu)進(jìn)行運(yùn)行狀態(tài)監(jiān)測(cè)和損傷識(shí)別至關(guān)重要。結(jié)構(gòu)健康監(jiān)測(cè)系統(tǒng)在運(yùn)營(yíng)期間會(huì)積累海量監(jiān)測(cè)數(shù)據(jù),基于大數(shù)據(jù)及深度學(xué)習(xí)思想的結(jié)構(gòu)健康監(jiān)測(cè)研究是當(dāng)前該領(lǐng)域的研究熱點(diǎn)[2-3]。

    深度學(xué)習(xí)具備精度高、魯棒性強(qiáng)和抗噪能力強(qiáng)等優(yōu)點(diǎn),利用大量數(shù)據(jù)對(duì)多層網(wǎng)絡(luò)模型進(jìn)行訓(xùn)練以得到精確、足夠的特征信息,而不受傳統(tǒng)神經(jīng)網(wǎng)絡(luò)的限制[4]。目前深度學(xué)習(xí)已廣泛運(yùn)用于各領(lǐng)域的研究,其中在結(jié)構(gòu)損傷識(shí)別領(lǐng)域主要是基于結(jié)構(gòu)振動(dòng)響應(yīng)的方法,利用損傷結(jié)構(gòu)的振動(dòng)響應(yīng)信號(hào)提取損傷特征,并以此建立信號(hào)與結(jié)構(gòu)損傷狀態(tài)的映射關(guān)系[5]。由于循環(huán)神經(jīng)網(wǎng)絡(luò)(recurren neural network,RNN)和長(zhǎng)短時(shí)記憶神經(jīng)網(wǎng)絡(luò)(long and short-term memory,LSTM)對(duì)于處理時(shí)間序列數(shù)據(jù)具有明顯優(yōu)勢(shì),目前已有不少運(yùn)用于機(jī)械故障診斷領(lǐng)域的研究。Talebi等[6]運(yùn)用角速度作為的輸入向量,建立了某風(fēng)能轉(zhuǎn)換結(jié)構(gòu)的故障監(jiān)測(cè)系統(tǒng),研究發(fā)現(xiàn)該方法在短期監(jiān)測(cè)中具有較好的檢測(cè)能力;趙志宏等[7]等提出一種基于雙向長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)模型BiLSTM的故障診斷的方法,將機(jī)器原始振動(dòng)信號(hào)作為模型輸入向量,訓(xùn)練提取不同故障類型的故障特征,結(jié)果驗(yàn)證了該方法的有效性;杜小磊等[8]將深層小波卷積自編碼分析與長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)相結(jié)合運(yùn)用于軸承故障診斷之中,結(jié)果表明該方法提取特征和識(shí)別故障能力優(yōu)于人工神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)等傳統(tǒng)方法;Cabrera等[9]將貝葉斯方法與LSTM模型相結(jié)合,提出一種壓縮機(jī)的故障檢測(cè)方法,結(jié)果故障識(shí)別的準(zhǔn)確率達(dá)到93%以上。

    超參數(shù)選取也是目前深度學(xué)習(xí)面臨的較為棘手的問題,模型超參數(shù)組合與網(wǎng)絡(luò)模型性能密切相關(guān)。貝葉斯優(yōu)化是一種十分有效的全局優(yōu)化算法,其目標(biāo)是找到優(yōu)化問題的全局最優(yōu)解[10]。目前,已有貝葉斯優(yōu)化運(yùn)用于各類模型參數(shù)優(yōu)選的研究,如Zhong等[11]提出通過應(yīng)用貝葉斯優(yōu)化搜索支持向量機(jī)模型的C和γ超參數(shù),用于工業(yè)測(cè)量數(shù)據(jù)分類。Wang等[12]提出了貝葉斯優(yōu)化結(jié)合粒子濾波的軸承故障特征提取方法等。

    損傷結(jié)構(gòu)的加速度響應(yīng)數(shù)據(jù)是時(shí)間序列數(shù)據(jù),可以反映結(jié)構(gòu)的狀態(tài)信息,蘊(yùn)含結(jié)構(gòu)的損傷特征,通過BiLSTM模型可以有效提取其中特征,解決損傷識(shí)別的問題。在輸電塔的健康監(jiān)測(cè)中,可以較為容易地測(cè)得結(jié)構(gòu)的加速度響應(yīng)序列,而雙向長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)可以有效地處理加速度響應(yīng)數(shù)據(jù),故本文采用了BiLSTM模型進(jìn)行輸電塔的損傷識(shí)別研究,并采用貝葉斯優(yōu)化選取模型的超參數(shù)組合,通過輸電塔仿真模擬和試驗(yàn)研究驗(yàn)證了該方法的有效性,并且相對(duì)其他模型具有更好的損傷識(shí)別效果。

    1 BiLSTM模型

    LSTM由RNN神經(jīng)網(wǎng)絡(luò)演變而來,在傳統(tǒng)RNN神經(jīng)網(wǎng)絡(luò)輸入層、隱藏層和輸出層的結(jié)構(gòu)基礎(chǔ)上,改進(jìn)了隱藏層的結(jié)構(gòu)。通過引入門機(jī)制來控制信息傳遞的路徑,并通過門對(duì)經(jīng)過網(wǎng)絡(luò)的信息進(jìn)行有選擇的記憶或刪除[13]。

    it=σ(Wixt+Uiht-1+bi)

    (1)

    ft=σ(Wfxt+Uiht-1+bf)

    (2)

    ot=σ(Woxt+Uoht-1+bo)

    (3)

    (4)

    (5)

    ht=ot?tanh(Ct)

    (6)

    式中:Wi、Wf和Wo表示從輸入門、遺忘門和輸出門到下一個(gè)輸入門的權(quán)重矩陣;Ui、Uf和Uo表示隱含層的權(quán)重矩陣;bi、bf和bo分別表示各個(gè)門結(jié)構(gòu)的偏置矩陣。σ為Sigmoid激活函數(shù)。

    圖1 LSTM結(jié)構(gòu)示意圖Fig.1 The structure of LSTM

    (7)

    (8)

    (9)

    圖2 BiLSTM結(jié)構(gòu)示意圖Fig.2 The structure of BiLSTM

    2 貝葉斯優(yōu)化的BiLSTM模型

    2.1 貝葉斯優(yōu)化

    應(yīng)用BiLSTM模型進(jìn)行損傷識(shí)別時(shí),需要確定模型最優(yōu)超參數(shù)組合。構(gòu)造一種優(yōu)秀的算法需要大量參數(shù)調(diào)節(jié)試驗(yàn),耗時(shí)耗力也不一定取得好效果。本文引入一種全局優(yōu)化算法——貝葉斯優(yōu)化(Bayesian optimization,BO)選取模型超參數(shù)組合。貝葉斯優(yōu)化基于貝葉斯定理,運(yùn)用概率代理模型擬合真實(shí)目標(biāo)函數(shù),并根據(jù)擬合的結(jié)果選擇接下來最有可能的點(diǎn)進(jìn)行評(píng)估,利用歷史信息合理減少評(píng)估次數(shù),提高搜索效率。

    貝葉斯優(yōu)化的框架主要包含兩個(gè)核心部分:概率代理模型和采集函數(shù)。概率代理模型包含先驗(yàn)概率模型和觀測(cè)模型;采集函數(shù)是根據(jù)后驗(yàn)概率分布構(gòu)造的函數(shù)[16],選用合適的概率代理模型和采集函數(shù)才能獲得更好的優(yōu)化效果。

    2.1.1 概率代理模型

    概率代理模型通常分為參數(shù)模型與非參數(shù)模型,非參數(shù)模型更為靈活,且不易發(fā)生“過擬合”,其中以高斯過程(Gaussian process,GP)應(yīng)用最為廣泛[17]。

    高斯過程是多元高斯概率分布的泛化[18],由一個(gè)均值函數(shù)m和一個(gè)協(xié)方差函數(shù)k構(gòu)成。在高斯過程中,有限個(gè)隨機(jī)變量都滿足一個(gè)高斯聯(lián)合分布。首先假設(shè)一個(gè)均值為0的先驗(yàn)分布

    p(f|X,θ)=N(0,Σ)

    (10)

    式中:X為訓(xùn)練集;f為未知函數(shù)f函數(shù)值的集合;Σ為k(x,x′)構(gòu)成的協(xié)方差矩陣;θ為超參數(shù)。

    假設(shè)存在的噪聲ε滿足獨(dú)立同分布的高斯分布:p(ε)=N(0,σ2),則似然分布為

    p(y|f)=N(f,σ2I)

    (11)

    式中,y=f(x)+ε,y表示觀測(cè)值y的集合。

    可得邊際似然分布為

    p(y|X,θ)=N(0,Σ+σ2I)

    (12)

    因此,存在聯(lián)合分布為

    (13)

    p(f*|X,y,X*)=N(〈f*〉,cov(f*))

    (14)

    (15)

    (16)

    式中:〈f*〉表示預(yù)測(cè)均值;cov(f*)是預(yù)測(cè)協(xié)方差。

    2.1.2 采集函數(shù)

    在貝葉斯優(yōu)化中,采集函數(shù)是用來選擇下一個(gè)最有可能的模型最佳性能評(píng)估點(diǎn)的方法。該函數(shù)通過已經(jīng)觀測(cè)到的數(shù)據(jù)集D1:t,得到后驗(yàn)分布構(gòu)造,并通過對(duì)其最大化來指導(dǎo)下一個(gè)評(píng)估點(diǎn)xt+1。期望提升(expected improvement,EI)的參數(shù)少且具有整合提升的概率,平衡深度和寬度的關(guān)系的優(yōu)點(diǎn),故本文采用EI作為采集函數(shù),其具體公式如下

    αt(x,D1:t)=

    (17)

    式中:v*為當(dāng)前最優(yōu)函數(shù)值;φ為標(biāo)準(zhǔn)正態(tài)分布概率密度函數(shù)。

    2.2 貝葉斯優(yōu)化的BiLSTM模型

    利用貝葉斯優(yōu)化對(duì)BiLSTM模型進(jìn)行超參數(shù)選取時(shí),可以考慮成以下公式進(jìn)行求解

    x*=arg minf(x),x∈X

    (18)

    式中:f(x)為目標(biāo)函數(shù);x表示一組超參數(shù)組合;X表示超參數(shù)組合空間;x*則表示X中的一組x,使目標(biāo)函數(shù)f(x)取得最優(yōu)解。本文以BiLSTM模型識(shí)別損傷的誤差率valError作為貝葉斯優(yōu)化的評(píng)價(jià)函數(shù)f(x),公式如下所示

    (19)

    式中:N為所有待檢測(cè)樣本組數(shù);R為N組樣本數(shù)據(jù)中所預(yù)測(cè)損傷類型與真實(shí)損傷類型相一致的樣本組數(shù)。

    本文采用高斯過程代理超參數(shù)組合x到BiLSTM模型的函數(shù)關(guān)系,根據(jù)已觀測(cè)的數(shù)據(jù)集得到后驗(yàn)分布構(gòu)造,然后通過EI函數(shù)選取下一個(gè)評(píng)估點(diǎn),迭代修正先驗(yàn)信息,逐步提高代理模型準(zhǔn)確性,尋找使目標(biāo)函數(shù)取得最優(yōu)解的超參數(shù)組合。

    BiLSTM模型超參數(shù)眾多,本文重點(diǎn)選取其中兩個(gè)超參數(shù)(初始學(xué)習(xí)率、小批量尺寸)進(jìn)行優(yōu)化選取。根據(jù)BiLSTM反向傳播傳播算法可知,其權(quán)重更新公式為

    (20)

    式中:α為學(xué)習(xí)率;n為小批量尺寸。由式(20)可知,除了梯度本身,學(xué)習(xí)率和小批量尺寸直接決定模型的權(quán)重更新,從優(yōu)化本身來看都是影響性能收斂最重要的參數(shù)。小批量尺寸決定模型每次參數(shù)更新?lián)p失函數(shù)的數(shù)據(jù)量,而學(xué)習(xí)率決定了權(quán)重迭代的步長(zhǎng)。對(duì)于不同的數(shù)據(jù)與模型,其最適的學(xué)習(xí)率和小批量尺寸也不同,沒有合適的選取法則,故本文采用貝葉斯優(yōu)化的方法選取。

    2.3 基于貝葉斯優(yōu)化的BiLSTM模型損傷識(shí)別方法

    2.3.1 模型執(zhí)行步驟

    運(yùn)用BO-BiLSTM模型進(jìn)行損傷識(shí)別的具體步驟如圖3所示。

    圖3 基于貝葉斯優(yōu)化BiLSTM模型的結(jié)構(gòu)損傷檢測(cè)流程Fig.3 Structural damage detection process based on BO- BiLSTM model

    步驟1首先獲取不同損傷工況下結(jié)構(gòu)每一層的加速度數(shù)據(jù),可以通過模型試驗(yàn)或者數(shù)值模擬獲得,使用z-score函數(shù)對(duì)數(shù)據(jù)標(biāo)準(zhǔn)化處理,建立樣本數(shù)據(jù)集,劃分訓(xùn)練集、測(cè)試集和驗(yàn)證集。

    步驟2建立BiLSTM模型,確定網(wǎng)絡(luò)模型的超參數(shù),并設(shè)置需要優(yōu)化的超參數(shù)(初始學(xué)習(xí)率、小批量尺寸)區(qū)間。

    步驟3進(jìn)行貝葉斯優(yōu)化,以BiLSTM網(wǎng)絡(luò)模型作為優(yōu)化的目標(biāo)函數(shù),以valError值作為評(píng)價(jià)函數(shù),使其到達(dá)最小值。

    步驟4在當(dāng)前超參數(shù)組合下,以訓(xùn)練集和驗(yàn)證集訓(xùn)練BiLSTM模型,計(jì)算當(dāng)前超參數(shù)下模型評(píng)價(jià)函數(shù)的值,完成后返回貝葉斯優(yōu)化,根據(jù)概率模型采集函數(shù)選取下一組超參數(shù)進(jìn)行新一輪訓(xùn)練,直至到達(dá)迭代次數(shù)。

    步驟5輸出貝葉斯優(yōu)化后模型性能最優(yōu)一組超參數(shù)及其BO-BiLSTM網(wǎng)絡(luò)模型。

    步驟6此時(shí)BO-BiLSTM模型已經(jīng)具備損傷檢測(cè)能力,將測(cè)試集數(shù)據(jù)輸入模型,輸出損傷狀態(tài),同時(shí)輸出網(wǎng)絡(luò)Softmax-score,驗(yàn)證模型性能。

    2.3.2 損傷評(píng)價(jià)指標(biāo)

    BO-BiLSTM模型進(jìn)行機(jī)械故障診斷或結(jié)構(gòu)損傷識(shí)別,從本質(zhì)上來說是一個(gè)多分類的問題。對(duì)于每一個(gè)待檢測(cè)的損傷數(shù)據(jù)樣本,通過網(wǎng)絡(luò)模型后可以得到一個(gè)N維數(shù)組作為輸出結(jié)果。數(shù)組中的每一個(gè)維度對(duì)應(yīng)一個(gè)結(jié)構(gòu)損傷類別,通過前向傳播算法得到輸出層的輸出緯度值代表屬于這個(gè)結(jié)構(gòu)損傷類別的可能性大小。Softmax回歸可以將神經(jīng)網(wǎng)絡(luò)的前向傳播輸出的數(shù)值分布轉(zhuǎn)化為概率分布,使每個(gè)樣本分類到不同損傷類別的概率和為1。

    (21)

    如式(21)所示,BO-BiLSTM模型原始輸出層的輸出值被作為置信度通過softmax層來生成新的輸出,該輸出滿足概率分布的所有要求。因此,使用softmax層輸出的各類損傷工況概率值Softmax-score作為本文方法中驗(yàn)證網(wǎng)絡(luò)損傷識(shí)別效果的評(píng)判指標(biāo)。

    3 輸電塔有限元模型損傷識(shí)別仿真

    3.1 輸電塔有限元模型

    為驗(yàn)證本文提出的BO-BiLSTM模型的有效性,通過ABAQUS軟件對(duì)一個(gè)小型的干字型輸電鐵塔進(jìn)行有限元模擬,以高斯白噪聲作為底部激勵(lì)獲取結(jié)構(gòu)在不同損傷工況下的加速度響應(yīng)數(shù)據(jù),并以此訓(xùn)練網(wǎng)絡(luò)模型,計(jì)算其損傷識(shí)別的準(zhǔn)確率。

    該模型采用梁?jiǎn)卧#⒌妮旊娝P图捌涑叽缭敿?xì)信息如圖4和圖5所示,輸電塔的主材、橫隔件以及橫擔(dān)采用鋼材,其斜材和交叉斜撐采用鋁材,具體材料與截面信息如表1所示。將輸電塔按分為8層子結(jié)構(gòu),損傷桿件的剛度折減系數(shù)設(shè)置為0.95,采用損傷不同桿件數(shù)量的方法模擬不同程度損傷工況,損傷工況如表2所示。本文基于Matlab2019b建立BO-BiLSTM模型,采用BiLSTM模型,BO-LSTM模型作為對(duì)比方法,驗(yàn)證其優(yōu)越性。

    圖4 輸電塔有限元模型Fig.4 The finite element model of transmission tower

    圖5 輸電塔模型尺寸平面圖Fig.5 Dimension plane of transmission tower model

    表1 構(gòu)件材料信息表Tab.1 The information of component material

    表2 損傷工況信息表Tab.2 The information of damage conditions

    3.2 數(shù)據(jù)集構(gòu)建

    本文研究其橫擔(dān)構(gòu)造之下的5層子結(jié)構(gòu)損傷識(shí)別情況。每種工況取如圖5所示6個(gè)測(cè)點(diǎn)的加速度響應(yīng)數(shù)據(jù)作為構(gòu)造特征的原始數(shù)據(jù),在不同白噪聲激勵(lì)下進(jìn)行10次模擬,每次白噪聲激勵(lì)時(shí)長(zhǎng)為25 s,取樣頻率為1 000 Hz,單一工況數(shù)據(jù)集長(zhǎng)度為250 000×6。共得到5種損傷位置(1-5層),2種不同損傷程度(損傷2桿和損傷4桿)的11種工況的數(shù)據(jù),將每種工況數(shù)據(jù)劃分為6×500的矩陣樣本,即每個(gè)樣本采樣時(shí)間為0.5 s,每種工況500個(gè)樣本。采用z-score方法對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化,11個(gè)工況共5 500個(gè)樣本。取整個(gè)數(shù)據(jù)集的前60%并打亂順序作為網(wǎng)絡(luò)的訓(xùn)練集;取后20%作為貝葉斯優(yōu)化的驗(yàn)證集,驗(yàn)證貝葉斯優(yōu)化效果;最后取數(shù)據(jù)集剩余20%作為模型的測(cè)試集,測(cè)試模型損傷識(shí)別效果,且防止模型對(duì)部分?jǐn)?shù)據(jù)過擬合,測(cè)試模型的泛化能力。

    3.3 BiLSTM模型架構(gòu)

    BiLSTM損傷識(shí)別模型架構(gòu)如圖6所示,網(wǎng)絡(luò)模型共分為5層,包含1個(gè)輸入層,1個(gè)BiLSTM隱含層,1個(gè)全連接層,1個(gè)softmax層以及1個(gè)分類輸出層。輸入數(shù)據(jù)首先進(jìn)入BiLSTM神經(jīng)單元中進(jìn)行計(jì)算,將其輸出值作為全連接層的輸入,計(jì)算后輸出值作為softmax層輸入值,輸出樣本數(shù)據(jù)分類至每種工況的概率分布,最后由分類輸出層輸出樣本對(duì)應(yīng)的損傷工況。

    圖6 BiLSTM損傷識(shí)別模型架構(gòu)圖Fig.6 The structure of BiLSTM damage identification model

    根據(jù)輸電塔有限元模型結(jié)構(gòu)實(shí)際情況以及所采集的加速度響應(yīng)數(shù)據(jù),對(duì)三個(gè)重要超參數(shù)(初始學(xué)習(xí)率、小批量尺寸、L2正則化系數(shù))設(shè)置優(yōu)化區(qū)間,其他超參數(shù)依據(jù)調(diào)參經(jīng)驗(yàn)進(jìn)行設(shè)置,具體參數(shù)設(shè)置如表3所示。

    表3 BiLSTM參數(shù)設(shè)置Tab.3 BiLSTM parameter settings

    對(duì)BiLSTM模型進(jìn)行貝葉斯優(yōu)化,迭代次數(shù)設(shè)置為30次,達(dá)到計(jì)算次數(shù)后停止優(yōu)化,輸出最優(yōu)超參數(shù)組合。尋優(yōu)過程如圖7所示,優(yōu)化至第30次迭代時(shí)模型性能達(dá)到最優(yōu),超參數(shù)組合為InitalLearnRate=0.014 511,MiniBatchsize=255,即模型在經(jīng)過貝葉斯優(yōu)化選取不同超參數(shù)組合后,BiLSTM模型整體識(shí)別準(zhǔn)確率由63.9%提升至92.9%。使用測(cè)試集數(shù)據(jù)對(duì)BO-BiLSTM模型進(jìn)行驗(yàn)證,結(jié)果發(fā)現(xiàn)該模型對(duì)測(cè)試集樣本損傷識(shí)別準(zhǔn)確率為94.20%,相對(duì)驗(yàn)證集提高1.34%,表明BO-BiLSTM模型具有很好的數(shù)據(jù)泛化能力,具備損傷識(shí)別能力。

    圖7 貝葉斯優(yōu)化尋優(yōu)過程Fig.7 Bayesian optimization process

    3.4 損傷識(shí)別結(jié)果分析

    使用測(cè)試集1 100個(gè)樣本對(duì)上節(jié)最優(yōu)BO-LSTM模型進(jìn)行損傷識(shí)別測(cè)試,分別從損傷位置以及損傷程度兩方面對(duì)模型損傷識(shí)別能力進(jìn)行進(jìn)一步分析。

    3.4.1 損傷位置分析

    輸電塔有限元試驗(yàn)設(shè)置的損傷情況為1~5層每層分別損傷4根桿件和損傷兩根桿件,故每層有100×2個(gè)測(cè)試樣本?,F(xiàn)將測(cè)試集分別經(jīng)過BiLSTM、BO-LSTM和BO-BiLSTM網(wǎng)絡(luò)模型識(shí)別,各層的識(shí)別結(jié)果及定位準(zhǔn)確率,如表4所示。

    表4 數(shù)值仿真損傷位置識(shí)別準(zhǔn)確率Tab.4 Damage location identification accuracy rate of numerical simulation

    根據(jù)表4結(jié)果顯示,BO-BiLSTM網(wǎng)絡(luò)模型對(duì)于輸電塔有限元模型層間損傷位置定位的整體準(zhǔn)確率在94.36%,整體識(shí)別效果相對(duì)其他兩個(gè)模型更加穩(wěn)定:BiLSTM模型對(duì)高層損傷識(shí)別準(zhǔn)確,而對(duì)底層損傷識(shí)別效果較差;BO-LSTM總體識(shí)別效果不錯(cuò),但是對(duì)第四層損傷識(shí)別效果不太理想。為進(jìn)一步分析模型性能,故提取出不同模型每層損傷工況樣本Softmax-score均值,其結(jié)果如圖8所示。

    圖8 數(shù)值仿真損傷位置識(shí)別結(jié)果Fig.8 Damage location identification results of numerical simulation

    Softmax-score反映模型判斷待檢測(cè)樣本到各種損傷位置的概率分布。根據(jù)圖分析可知,各層損傷的Softmax-score值隨損傷層數(shù)的增加而上升,表明各模型對(duì)高層損傷的損傷特征提取更為準(zhǔn)確,對(duì)高層損傷的定位更為精準(zhǔn)。對(duì)比不同模型Softmax-score,BO-BiLSTM模型基本處于各模型最高值,而BO-LSTM模型對(duì)第4層損傷的識(shí)別出現(xiàn)誤差,易與第3層損傷產(chǎn)生混淆。綜合以上對(duì)比結(jié)果可知,經(jīng)過訓(xùn)練后的BO-BiLSTM模型可以準(zhǔn)確定位出輸電塔模型層間損傷的位置,定位效果相比另兩個(gè)模型更為穩(wěn)定。

    3.4.2 損傷模式分析

    主要考慮了輸電塔每層損傷的兩種損傷模式,即損傷2根桿件和4根桿件。故五層塔架結(jié)構(gòu)考慮了10種損傷模式,具體如表5所示。該表包含10種損傷模式1~10,而模式0代表無損傷模式。表5還包含各長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)模型對(duì)不同損傷工況的識(shí)別結(jié)果與準(zhǔn)確率。

    由上述結(jié)果可以看出,BO-BiLSTM模型可以在定位輸電塔層間損傷的基礎(chǔ)上,區(qū)分相同損傷位置下不同模式的損傷。單層損傷4根桿件模式的識(shí)別準(zhǔn)確率基本達(dá)到100%;損傷2桿模式的準(zhǔn)確率有所下降,但也維持在84%以上。同時(shí),隨著損傷所在層數(shù)的增加,模型對(duì)損傷模式識(shí)別準(zhǔn)確率也隨之上升。對(duì)比三個(gè)模型結(jié)果可知,各模型對(duì)強(qiáng)損傷模式均保有較高的識(shí)別準(zhǔn)確率;對(duì)于弱損傷模式,BO-BiLSTM模型識(shí)別效果更好,未進(jìn)行貝葉斯優(yōu)化的BiLSTM模型對(duì)低層弱損傷模式失去判斷能力,提取特征能力較差。為進(jìn)一步分析各模型性能,將不同程度模式下Softmax-score提取出來如圖9所示,其中模式0表示結(jié)構(gòu)未損傷狀態(tài)。

    表5 數(shù)值仿真損傷模式識(shí)別準(zhǔn)確率Tab.5 Damage pattern identification accuracy rate of numerical simulation

    圖9 數(shù)值仿真損傷模式結(jié)果Fig.9 Damage pattern identification results of numerical simulation

    整體分析,三種模型對(duì)弱損傷模式識(shí)別效果相對(duì)較差,如模式6(第1層損傷2桿)情況下,模型對(duì)該模式損傷判定較為模糊,分類至模式6的概率值僅略微高于其他模式,表明模型對(duì)該模式下?lián)p傷的特征提取效果不夠好,容易出現(xiàn)誤判現(xiàn)象;而對(duì)強(qiáng)損傷模式的Softmax-score均接近于1,基本都可以準(zhǔn)確判別損傷情況。結(jié)合識(shí)別準(zhǔn)確率分析,BO-BiLSTM模型對(duì)輸電塔模型層間損傷的損傷特征提取更為準(zhǔn)確,特別是對(duì)小損傷模式識(shí)別效果優(yōu)于另兩個(gè)模型,經(jīng)過訓(xùn)練后對(duì)輸電塔損傷具有良好的診斷能力。

    4 輸電塔損傷識(shí)別試驗(yàn)研究

    4.1 試驗(yàn)介紹

    為驗(yàn)證本文所提出模型對(duì)實(shí)際結(jié)構(gòu)的損傷識(shí)別效果,本節(jié)嘗試將有限元模擬與試驗(yàn)相結(jié)合,即先通過有限元模型獲取結(jié)構(gòu)在各種損傷工況下的加速度響應(yīng)數(shù)據(jù)訓(xùn)練BO-BiLSTM模型,再通過試驗(yàn)獲取的損傷數(shù)據(jù)組成測(cè)試集,以此測(cè)試BO-BiLSTM對(duì)實(shí)際輸電塔結(jié)構(gòu)損傷提取的魯棒性。本文將以第3章的輸電塔有限元模型為例,設(shè)計(jì)制作輸電塔結(jié)構(gòu)試驗(yàn)對(duì)本文所提方法進(jìn)行研究。

    輸電塔試驗(yàn)?zāi)P统叽缗c材料均與有限元模型一致,輸電塔共分為8層,塔桿總高度為2.6 m;第1~5層高度均為400 mm;第6~8層每層高度均為200 mm,塔身四肢主材為空心鋼管,直徑為10 mm,所有橫隔件為實(shí)心圓鋼,直徑為6 mm;塔身四面斜材以及每層橫截面方向上交叉支撐由正方形截面實(shí)心鋁桿組成。如圖10所示,塔身主材每層節(jié)點(diǎn)上焊接有空的薄鋼片,通過螺栓將斜材、交叉支撐相連接。橫擔(dān)件通過焊接的方法與主材相連接;輸電塔的每個(gè)柱腳焊接在一個(gè)10 mm厚的鋼板上,每個(gè)鋼板通過4顆螺栓與振動(dòng)臺(tái)相連接,以確保固定方式為固接。

    (a) 試驗(yàn)?zāi)P?/p>

    本次試驗(yàn)利用小型振動(dòng)臺(tái)對(duì)結(jié)構(gòu)進(jìn)行加載,在y方向上施加10段隨機(jī)生成帶寬為20~50 Hz的25 s高斯白噪聲作為結(jié)構(gòu)的激勵(lì)。通過東華動(dòng)態(tài)信號(hào)測(cè)試分析系統(tǒng)以及加速度傳感器采集結(jié)構(gòu)第1~5層和第7層的加速度響應(yīng)數(shù)據(jù)。加速度傳感器的布置位置與有限元模型中所選取的節(jié)點(diǎn)保持一致,易于模型的識(shí)別匹配。

    通過拆除單根或多根塔身斜材來模擬結(jié)構(gòu)層間損傷,而由于鋁材輕質(zhì)的特點(diǎn),可以減少由于拆除構(gòu)件后對(duì)結(jié)構(gòu)整體質(zhì)量的影響。本次試驗(yàn)一共設(shè)置了六種損傷工況和對(duì)應(yīng)模式,分別考慮第3層、第4層和第5層發(fā)生不同程度損傷的情況,損傷工況與有限元模擬對(duì)應(yīng)工況如表6所示。

    表6 試驗(yàn)的損傷工況設(shè)置Tab.6 States of transmission tower model experiment

    為避免振動(dòng)臺(tái)剛啟動(dòng)時(shí)和即將結(jié)束時(shí)產(chǎn)生的誤差對(duì)結(jié)果的影響,本文對(duì)每個(gè)通道獲得的加速度數(shù)據(jù)做截?cái)嗵幚?,將試?yàn)采集的加速度數(shù)據(jù)去掉前后,截取中間100 000個(gè)數(shù)據(jù)進(jìn)行研究。圖11是試驗(yàn)工況6下通道1、3、5采集的加速度響應(yīng)數(shù)據(jù)節(jié)選。將試驗(yàn)每種損傷工況采集的數(shù)據(jù)劃分為6×500的矩陣樣本,采用z-score對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化構(gòu)成驗(yàn)證集。

    圖11 試驗(yàn)工況6下通道1、3、5加速度響應(yīng)數(shù)據(jù)節(jié)選Fig.11 Part of the accelerations data of the channel 1,3 and 5 at test case 6

    4.2 試驗(yàn)結(jié)果分析

    為驗(yàn)證3.3節(jié)訓(xùn)練所得模型對(duì)實(shí)際結(jié)構(gòu)的損傷識(shí)別效果,將試驗(yàn)所得樣本直接作為測(cè)試集樣本輸入上節(jié)訓(xùn)練完成BO-BiLSTM模型進(jìn)行損傷識(shí)別,整體識(shí)別準(zhǔn)確率達(dá)到了79.6%。下面將具體從損傷位置和損傷程度對(duì)結(jié)果進(jìn)行分析。

    4.2.1 損傷位置識(shí)別分析

    現(xiàn)將測(cè)試集經(jīng)過BO-BiLSTM、BiLSTM和BO-LSTM網(wǎng)絡(luò)模型識(shí)別,各層的識(shí)別結(jié)果及定位準(zhǔn)確率,如表7所示,各層Softmax-score如圖12所示。

    表7 試驗(yàn)損傷位置識(shí)別準(zhǔn)確率Tab.7 Damage location identification accuracy rate of experiment

    根據(jù)表7結(jié)果顯示,使用輸電塔試驗(yàn)數(shù)據(jù)作為測(cè)試集,測(cè)試通過有限元模擬數(shù)據(jù)訓(xùn)練出的BO-BiLSTM模型,整體的損傷定位準(zhǔn)確率為82.8%,相對(duì)有限元結(jié)果有所下降,但仍保持有較高準(zhǔn)確率。層間損傷的定位準(zhǔn)確率從第3層的80.75%逐漸上升至第5層84.5%,準(zhǔn)確率隨層數(shù)的增加而上升。對(duì)比三個(gè)模型識(shí)別結(jié)果可以發(fā)現(xiàn),BO-BiLSTM模型整體識(shí)別效果更穩(wěn)定,而BiLSTM和BO-LSTM模型在某些層均出現(xiàn)了較大波動(dòng),準(zhǔn)確率較低。

    分析各模型的Softmax-score可知,各模型對(duì)試驗(yàn)層間損傷的定位表現(xiàn)結(jié)果與有限元模型結(jié)果基本相同,經(jīng)過貝葉斯優(yōu)化的BO-BiLSTM模型整體保持較高性能,可以較為準(zhǔn)確提取不同位置損傷工況的定位特征;BiLSTM模型則對(duì)較低層的損傷特征提取不夠明顯,識(shí)別效果較差;BO-LSTM模型則對(duì)第4層損傷特征提取出現(xiàn)較大誤差,易與第3層損傷混淆。

    圖12 試驗(yàn)的損傷位置識(shí)別結(jié)果Fig.12 Damage location identification results of experiment

    4.2.2 損傷模式識(shí)別分析

    在試驗(yàn)中每層損傷也設(shè)置了兩種損傷模式,即2桿損傷和4桿損傷。在本文第3章中,輸電塔有限元模型的桿件損傷通過剛度折減0.95的方法進(jìn)行模擬。由于試驗(yàn)條件限制,試驗(yàn)中采取拆除桿件的方法模擬結(jié)構(gòu)的剛度變化,故存在一定誤差。試驗(yàn)中設(shè)置的2種損傷模式,與有限元模型模擬工況一致。表8所示為三種長(zhǎng)短記憶網(wǎng)絡(luò)模型對(duì)不同損傷工況的識(shí)別結(jié)果準(zhǔn)確率,圖13為各試驗(yàn)工況Softmax-score。

    表8 試驗(yàn)的損傷模式識(shí)別準(zhǔn)確率Tab.8 Damage pattern identification accuracy rate of experiment

    從識(shí)別準(zhǔn)確率不難看出,各模型損傷四根桿件桿工況準(zhǔn)確率普遍高于損傷兩個(gè)桿件準(zhǔn)確率,且其Softmax- score也更為突出,一定程度上表明模型對(duì)強(qiáng)損傷工況更為敏感,但整體而言BO-BiLSTM模型識(shí)別精度依舊保持良好。具體從各模型結(jié)果分析,雖然BO-LSTM與BiLSTM模型對(duì)損傷4根桿件識(shí)別效果略優(yōu)于BO-BiLSTM模型,但實(shí)際準(zhǔn)確率上差別不大;而對(duì)于損傷兩根桿件的工況而言,BiLSTM模型和BO-LSTM模型對(duì)部分工況識(shí)別效果很差,甚至完全失去判斷能力,而BO-BiLSTM模型仍保有較為不錯(cuò)的識(shí)別精度,識(shí)別更加穩(wěn)定。

    對(duì)比試驗(yàn)與有限元模型準(zhǔn)確率,由于訓(xùn)練BO-BiLSTM模型的有限元數(shù)據(jù)并未包含噪聲,因此有限元模擬數(shù)據(jù)集的識(shí)別精度都高于試驗(yàn)數(shù)據(jù)集,但兩者精度并未存在過大差距。在實(shí)際試驗(yàn)過程中,試驗(yàn)場(chǎng)地具有較大且復(fù)雜的環(huán)境噪聲(如人員走動(dòng)、液壓機(jī)的振動(dòng)、開關(guān)門、空調(diào)等等),BO-BiLSTM模型對(duì)試驗(yàn)測(cè)試集識(shí)別精度仍保持在一個(gè)較高水平,由此可知BO-BiLSTM模型在結(jié)構(gòu)的損傷識(shí)別的實(shí)際應(yīng)用過程中具有較好的抗噪能力和魯棒性。

    圖13 試驗(yàn)的損傷模式識(shí)別結(jié)果Fig.13 Damage pattern identification results of experiment

    5 結(jié) 論

    本文提出了一種基于貝葉斯優(yōu)化的BiLSTM模型的損傷識(shí)別方法,該方法運(yùn)用貝葉斯優(yōu)化選取BiLSTM模型的超參數(shù)組合,通過結(jié)構(gòu)加速度響應(yīng)數(shù)據(jù)訓(xùn)練模型,使其具備損傷識(shí)別能力。采用輸電塔有限元模型損傷試驗(yàn)驗(yàn)證了提出方法的有效性,并對(duì)一個(gè)實(shí)際輸電塔結(jié)構(gòu)進(jìn)行了試驗(yàn)研究,結(jié)果表明:

    (1) 本文提出的損傷識(shí)別方法能夠有效識(shí)別輸電塔結(jié)構(gòu)的層間損傷。以結(jié)構(gòu)不同損傷狀態(tài)下加速度響應(yīng)作為輸入數(shù)據(jù),經(jīng)過訓(xùn)練后的BO-BiLSTM模型可以判斷其損傷的位置與損傷模式。

    (2) BiLSTM模型的超參數(shù)與模型性能密切相關(guān)。從輸電塔有限元試驗(yàn)中得知,貝葉斯優(yōu)化可以有效的運(yùn)用于BiLSTM模型的超參數(shù)選取,充分利用歷史信息選取最優(yōu)超參數(shù)組合,提高搜索效率。

    (3) BO-BiLSTM模型可以采用異源數(shù)據(jù)(有限元模擬和振動(dòng)臺(tái)試驗(yàn))進(jìn)行損傷識(shí)別。有限元模擬的方法可以獲得大量不同損傷工況下的加速度樣本作為訓(xùn)練集,訓(xùn)練用于輸電塔損傷識(shí)別的BO-BiLSTM模型,提取結(jié)構(gòu)損傷特征信息,再將之應(yīng)用于實(shí)際結(jié)構(gòu)的損傷診斷之中。

    本研究通過有限元模型模擬部分損傷工況,以此作為訓(xùn)練集訓(xùn)練長(zhǎng)短時(shí)記憶網(wǎng)絡(luò),并以試驗(yàn)數(shù)據(jù)驗(yàn)證該方法的有效性。由于實(shí)際可能出現(xiàn)的工況種類眾多,模型模擬的數(shù)據(jù)量龐大,對(duì)網(wǎng)絡(luò)模型的要求更高,受試驗(yàn)條件限制本文模型試驗(yàn)中無法涵蓋所有可能出現(xiàn)的損傷工況,故應(yīng)用于實(shí)際中仍有待進(jìn)一步研究。

    猜你喜歡
    貝葉斯準(zhǔn)確率有限元
    乳腺超聲檢查診斷乳腺腫瘤的特異度及準(zhǔn)確率分析
    健康之家(2021年19期)2021-05-23 11:17:39
    不同序列磁共振成像診斷脊柱損傷的臨床準(zhǔn)確率比較探討
    2015—2017 年寧夏各天氣預(yù)報(bào)參考產(chǎn)品質(zhì)量檢驗(yàn)分析
    高速公路車牌識(shí)別標(biāo)識(shí)站準(zhǔn)確率驗(yàn)證法
    貝葉斯公式及其應(yīng)用
    基于貝葉斯估計(jì)的軌道占用識(shí)別方法
    一種基于貝葉斯壓縮感知的說話人識(shí)別方法
    電子器件(2015年5期)2015-12-29 08:43:15
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    IIRCT下負(fù)二項(xiàng)分布參數(shù)多變點(diǎn)的貝葉斯估計(jì)
    精品一区二区三区视频在线| 国产永久视频网站| 亚洲国产精品专区欧美| 日韩欧美精品免费久久| 又爽又黄无遮挡网站| 日韩一区二区视频免费看| 插逼视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 欧美日韩国产mv在线观看视频 | 综合色丁香网| 啦啦啦啦在线视频资源| 欧美激情国产日韩精品一区| 蜜臀久久99精品久久宅男| 国产精品人妻久久久影院| 一区二区三区免费毛片| 亚洲美女搞黄在线观看| 国产男女超爽视频在线观看| 亚洲av男天堂| 亚洲婷婷狠狠爱综合网| 国产91av在线免费观看| 亚洲最大成人中文| 国产极品天堂在线| 亚洲精品乱码久久久久久按摩| 嘟嘟电影网在线观看| 又爽又黄无遮挡网站| 亚洲av电影不卡..在线观看| 亚洲成人一二三区av| 又爽又黄a免费视频| 免费看不卡的av| 日本-黄色视频高清免费观看| 精华霜和精华液先用哪个| 国产精品久久久久久精品电影| 欧美日韩视频高清一区二区三区二| 日韩亚洲欧美综合| 蜜桃亚洲精品一区二区三区| 亚洲精品456在线播放app| 亚洲天堂国产精品一区在线| 亚洲性久久影院| 你懂的网址亚洲精品在线观看| 国产精品一区二区在线观看99 | 亚洲精品国产成人久久av| 日韩av免费高清视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产大屁股一区二区在线视频| 久久精品国产亚洲av天美| 男人狂女人下面高潮的视频| 熟妇人妻不卡中文字幕| 高清在线视频一区二区三区| 中文在线观看免费www的网站| 99久久人妻综合| 免费黄频网站在线观看国产| 国产一区有黄有色的免费视频 | 亚洲欧美精品专区久久| 国产成年人精品一区二区| 一区二区三区乱码不卡18| 极品少妇高潮喷水抽搐| av网站免费在线观看视频 | 丝袜喷水一区| 纵有疾风起免费观看全集完整版 | 欧美极品一区二区三区四区| 九九在线视频观看精品| 国产亚洲91精品色在线| 亚洲欧美中文字幕日韩二区| 国产精品不卡视频一区二区| 少妇被粗大猛烈的视频| 国产在视频线精品| 免费少妇av软件| 天堂中文最新版在线下载 | 免费大片18禁| 国产v大片淫在线免费观看| 十八禁网站网址无遮挡 | 日本一本二区三区精品| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 建设人人有责人人尽责人人享有的 | 有码 亚洲区| 简卡轻食公司| 成人午夜精彩视频在线观看| freevideosex欧美| 欧美zozozo另类| 国产91av在线免费观看| 国产午夜福利久久久久久| 五月玫瑰六月丁香| 国产精品一区www在线观看| 波多野结衣巨乳人妻| 午夜日本视频在线| 少妇裸体淫交视频免费看高清| 男人狂女人下面高潮的视频| 熟女电影av网| 国产亚洲最大av| 天堂俺去俺来也www色官网 | 一级片'在线观看视频| 简卡轻食公司| 人妻制服诱惑在线中文字幕| eeuss影院久久| av一本久久久久| 欧美成人一区二区免费高清观看| 99久久精品国产国产毛片| 国产乱来视频区| 秋霞伦理黄片| 网址你懂的国产日韩在线| 午夜福利在线观看免费完整高清在| 国产一区亚洲一区在线观看| 少妇熟女欧美另类| 成年版毛片免费区| 亚洲熟妇中文字幕五十中出| 亚洲精品第二区| 国产午夜精品论理片| 精品久久久久久久久亚洲| 一个人观看的视频www高清免费观看| av线在线观看网站| 男的添女的下面高潮视频| 综合色av麻豆| 超碰97精品在线观看| 国产精品一区二区三区四区免费观看| xxx大片免费视频| 国产精品一区www在线观看| 亚洲欧美成人精品一区二区| 国产一区二区三区av在线| 一夜夜www| 欧美成人午夜免费资源| 在线观看av片永久免费下载| 69av精品久久久久久| av线在线观看网站| 神马国产精品三级电影在线观看| 国产成人a∨麻豆精品| 国产亚洲av片在线观看秒播厂 | 国产成人a∨麻豆精品| 国产色爽女视频免费观看| 午夜福利网站1000一区二区三区| 日韩成人av中文字幕在线观看| 蜜桃久久精品国产亚洲av| 最近中文字幕2019免费版| 又粗又硬又长又爽又黄的视频| 99热6这里只有精品| 国产成人freesex在线| 色综合亚洲欧美另类图片| 在现免费观看毛片| 亚洲欧美日韩东京热| 99久久精品热视频| 国产中年淑女户外野战色| 国产探花极品一区二区| 中文精品一卡2卡3卡4更新| 免费黄色在线免费观看| 亚洲最大成人中文| 五月玫瑰六月丁香| 晚上一个人看的免费电影| 国产高清不卡午夜福利| 一级a做视频免费观看| 免费在线观看成人毛片| 三级经典国产精品| 淫秽高清视频在线观看| 99久久人妻综合| 寂寞人妻少妇视频99o| 国产高清国产精品国产三级 | 日本免费在线观看一区| 一级二级三级毛片免费看| 久久久久久久国产电影| 亚洲激情五月婷婷啪啪| 欧美xxxx黑人xx丫x性爽| 99热这里只有是精品50| 亚洲精品视频女| 一级毛片久久久久久久久女| 丝袜喷水一区| 亚洲精品日韩在线中文字幕| 精品久久久久久久末码| 91精品一卡2卡3卡4卡| 国内揄拍国产精品人妻在线| 成年人午夜在线观看视频 | 一级a做视频免费观看| 91精品国产九色| 高清午夜精品一区二区三区| 校园人妻丝袜中文字幕| 成年女人在线观看亚洲视频 | 国精品久久久久久国模美| 国产综合懂色| 日本色播在线视频| 国产美女午夜福利| 精品国产一区二区三区久久久樱花 | 亚洲精品一二三| 日韩不卡一区二区三区视频在线| 人妻少妇偷人精品九色| 国产老妇女一区| 国产精品爽爽va在线观看网站| 国产91av在线免费观看| 久久99精品国语久久久| 在线 av 中文字幕| 高清视频免费观看一区二区 | 中文字幕av在线有码专区| 26uuu在线亚洲综合色| 精品欧美国产一区二区三| 麻豆成人av视频| 免费av毛片视频| 三级毛片av免费| 岛国毛片在线播放| 搡老乐熟女国产| 午夜福利成人在线免费观看| 国产极品天堂在线| 午夜老司机福利剧场| 亚洲av中文字字幕乱码综合| 丰满人妻一区二区三区视频av| 97精品久久久久久久久久精品| 国产成人精品福利久久| 国产精品日韩av在线免费观看| 女的被弄到高潮叫床怎么办| 日韩精品青青久久久久久| 亚洲精品影视一区二区三区av| 日韩欧美一区视频在线观看 | 联通29元200g的流量卡| 美女国产视频在线观看| 欧美成人精品欧美一级黄| 国产午夜福利久久久久久| 亚洲av一区综合| 国产色爽女视频免费观看| 中文欧美无线码| 国内精品宾馆在线| 亚洲在久久综合| 亚洲av福利一区| 女的被弄到高潮叫床怎么办| 国产av在哪里看| 男插女下体视频免费在线播放| 久久亚洲国产成人精品v| 亚洲国产精品专区欧美| 18禁动态无遮挡网站| 成人毛片a级毛片在线播放| 69av精品久久久久久| 午夜免费观看性视频| 国产精品av视频在线免费观看| 亚洲人与动物交配视频| 身体一侧抽搐| 丰满少妇做爰视频| 亚洲性久久影院| 国产一区亚洲一区在线观看| 国产在线一区二区三区精| 日产精品乱码卡一卡2卡三| 一区二区三区乱码不卡18| 欧美人与善性xxx| 国产 亚洲一区二区三区 | 一本久久精品| 男女那种视频在线观看| 99视频精品全部免费 在线| av免费观看日本| 少妇裸体淫交视频免费看高清| 久久99热6这里只有精品| 久久久精品免费免费高清| 永久免费av网站大全| 日日撸夜夜添| 欧美日韩国产mv在线观看视频 | 日韩av免费高清视频| 精品人妻熟女av久视频| av专区在线播放| 亚洲精品aⅴ在线观看| 亚洲av成人av| 一本久久精品| 久久久久久久亚洲中文字幕| 边亲边吃奶的免费视频| 嫩草影院精品99| 亚洲一区高清亚洲精品| 欧美激情久久久久久爽电影| 亚洲熟妇中文字幕五十中出| 亚洲一区高清亚洲精品| 欧美不卡视频在线免费观看| 美女大奶头视频| 美女xxoo啪啪120秒动态图| 国产乱来视频区| 特级一级黄色大片| 婷婷六月久久综合丁香| 18禁裸乳无遮挡免费网站照片| 久久久a久久爽久久v久久| 欧美最新免费一区二区三区| 亚洲av成人精品一二三区| 精品久久久精品久久久| 午夜久久久久精精品| 你懂的网址亚洲精品在线观看| 亚洲欧美成人综合另类久久久| 大片免费播放器 马上看| 亚洲图色成人| av.在线天堂| 97超视频在线观看视频| 亚洲国产精品成人久久小说| 亚洲丝袜综合中文字幕| 人人妻人人澡欧美一区二区| 欧美zozozo另类| 丝袜喷水一区| 国产极品天堂在线| 亚洲精品日本国产第一区| 国产精品1区2区在线观看.| 一个人免费在线观看电影| 欧美潮喷喷水| 26uuu在线亚洲综合色| 一级毛片我不卡| 看免费成人av毛片| 国产黄片视频在线免费观看| av在线亚洲专区| 夫妻午夜视频| 欧美潮喷喷水| 日本猛色少妇xxxxx猛交久久| 日韩欧美国产在线观看| 身体一侧抽搐| 日韩亚洲欧美综合| 精品一区二区三区人妻视频| 国产亚洲精品久久久com| 亚洲成色77777| 啦啦啦中文免费视频观看日本| 赤兔流量卡办理| 亚洲av免费在线观看| 欧美成人午夜免费资源| 热99在线观看视频| 狠狠精品人妻久久久久久综合| 自拍偷自拍亚洲精品老妇| 国产黄色视频一区二区在线观看| 欧美精品一区二区大全| 午夜爱爱视频在线播放| 精品人妻偷拍中文字幕| 日韩 亚洲 欧美在线| 午夜福利在线观看免费完整高清在| 18禁在线无遮挡免费观看视频| 97超碰精品成人国产| 久久久久久国产a免费观看| 天堂网av新在线| 99热网站在线观看| 国产在视频线在精品| 国产不卡一卡二| 成人鲁丝片一二三区免费| 成人毛片60女人毛片免费| 久久久久九九精品影院| 人体艺术视频欧美日本| 插阴视频在线观看视频| 高清午夜精品一区二区三区| 国产成人a∨麻豆精品| 特级一级黄色大片| 免费观看性生交大片5| 亚洲性久久影院| 久久久精品欧美日韩精品| 久久这里有精品视频免费| 久久久国产一区二区| 亚洲在线自拍视频| 亚洲18禁久久av| 十八禁国产超污无遮挡网站| xxx大片免费视频| 国产色爽女视频免费观看| 亚洲av男天堂| 国产精品麻豆人妻色哟哟久久 | 国产乱人偷精品视频| 亚洲欧洲国产日韩| 日本三级黄在线观看| 在线观看一区二区三区| 欧美高清性xxxxhd video| av国产免费在线观看| 又爽又黄无遮挡网站| av播播在线观看一区| 麻豆精品久久久久久蜜桃| 国产黄色免费在线视频| 成年免费大片在线观看| 国产午夜精品久久久久久一区二区三区| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 99久国产av精品| 黄片无遮挡物在线观看| 极品教师在线视频| 成人国产麻豆网| 亚洲国产成人一精品久久久| 亚洲美女视频黄频| 国产大屁股一区二区在线视频| 亚洲av.av天堂| 老师上课跳d突然被开到最大视频| 国产成人91sexporn| 久久久久久久久久久丰满| 九九久久精品国产亚洲av麻豆| 亚洲图色成人| 97超视频在线观看视频| 成年女人在线观看亚洲视频 | 最近的中文字幕免费完整| 一区二区三区高清视频在线| 51国产日韩欧美| 亚洲av中文av极速乱| 少妇熟女欧美另类| 亚洲国产精品sss在线观看| 久热久热在线精品观看| 亚洲av成人精品一二三区| 亚洲精品一区蜜桃| 亚洲精品色激情综合| 亚洲在线观看片| 七月丁香在线播放| 亚洲av成人精品一区久久| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 久久久久久国产a免费观看| 欧美变态另类bdsm刘玥| 美女xxoo啪啪120秒动态图| 久久久久久久久久成人| 成人午夜高清在线视频| 人人妻人人澡人人爽人人夜夜 | 国产成人精品福利久久| 欧美日韩国产mv在线观看视频 | 看十八女毛片水多多多| 国产在视频线精品| 男的添女的下面高潮视频| 欧美xxxx性猛交bbbb| 亚洲精品视频女| 99热这里只有是精品在线观看| 日韩欧美精品v在线| 91精品国产九色| 免费播放大片免费观看视频在线观看| 午夜久久久久精精品| 国产一区亚洲一区在线观看| 国产激情偷乱视频一区二区| 欧美性感艳星| 99久久九九国产精品国产免费| 亚洲精品影视一区二区三区av| 国产有黄有色有爽视频| 国产精品久久久久久精品电影小说 | 国产成人精品一,二区| av又黄又爽大尺度在线免费看| 久久精品熟女亚洲av麻豆精品 | 国产成人a区在线观看| 中文精品一卡2卡3卡4更新| av专区在线播放| 最近视频中文字幕2019在线8| 熟妇人妻久久中文字幕3abv| 久久6这里有精品| 久久久久久久久久久免费av| 99久久精品国产国产毛片| 边亲边吃奶的免费视频| 国产亚洲一区二区精品| 亚洲欧美精品专区久久| 国产成人aa在线观看| 久久久久久久亚洲中文字幕| 免费av观看视频| 久久久久久伊人网av| 免费看av在线观看网站| 日韩人妻高清精品专区| 在线免费十八禁| 午夜福利在线观看免费完整高清在| 精品久久久久久久末码| 欧美激情久久久久久爽电影| 欧美日韩视频高清一区二区三区二| 久久久欧美国产精品| 国产av不卡久久| 中文字幕亚洲精品专区| 日日啪夜夜撸| 三级国产精品片| 国产av在哪里看| 国产人妻一区二区三区在| 岛国毛片在线播放| 国产综合精华液| av专区在线播放| 免费不卡的大黄色大毛片视频在线观看 | 国产乱人偷精品视频| 日韩中字成人| 国产精品国产三级国产专区5o| av在线蜜桃| 最近手机中文字幕大全| 免费av不卡在线播放| 日本爱情动作片www.在线观看| 亚洲精品一二三| 97在线视频观看| 一级毛片aaaaaa免费看小| 日韩三级伦理在线观看| 最近最新中文字幕大全电影3| 免费电影在线观看免费观看| 亚洲av中文字字幕乱码综合| 国产av国产精品国产| 国产精品久久久久久久电影| 三级毛片av免费| 美女主播在线视频| 国产精品麻豆人妻色哟哟久久 | 国产精品.久久久| 久久国内精品自在自线图片| 日日啪夜夜撸| 色尼玛亚洲综合影院| 一级片'在线观看视频| 91狼人影院| 嫩草影院入口| 久久久久久九九精品二区国产| av女优亚洲男人天堂| 神马国产精品三级电影在线观看| 99热这里只有精品一区| 国产av不卡久久| 欧美人与善性xxx| 一区二区三区免费毛片| 欧美精品一区二区大全| 极品教师在线视频| 亚洲国产欧美在线一区| 日日摸夜夜添夜夜添av毛片| 欧美成人一区二区免费高清观看| 日韩av不卡免费在线播放| 久久久久久久大尺度免费视频| 我的女老师完整版在线观看| 中文字幕av在线有码专区| 免费观看无遮挡的男女| 国产免费一级a男人的天堂| 青春草国产在线视频| 久久这里只有精品中国| 欧美潮喷喷水| 日日干狠狠操夜夜爽| 日本-黄色视频高清免费观看| 人妻系列 视频| 99久久精品热视频| 小蜜桃在线观看免费完整版高清| 国产成人a∨麻豆精品| 边亲边吃奶的免费视频| 亚洲av不卡在线观看| 亚洲国产欧美在线一区| 亚洲在久久综合| 禁无遮挡网站| 日本与韩国留学比较| 久久精品国产鲁丝片午夜精品| 国产 亚洲一区二区三区 | 十八禁国产超污无遮挡网站| 欧美性感艳星| 91午夜精品亚洲一区二区三区| 国产淫片久久久久久久久| 国产精品国产三级专区第一集| 99久久人妻综合| 欧美另类一区| 伦理电影大哥的女人| 男人爽女人下面视频在线观看| 97人妻精品一区二区三区麻豆| 久久久久久久亚洲中文字幕| 午夜亚洲福利在线播放| 肉色欧美久久久久久久蜜桃 | 国产精品人妻久久久久久| 欧美成人a在线观看| 国产永久视频网站| 亚洲av电影不卡..在线观看| 精品久久久久久久久av| 欧美不卡视频在线免费观看| 水蜜桃什么品种好| 成人国产麻豆网| 久久久精品免费免费高清| 日日摸夜夜添夜夜爱| 亚洲av日韩在线播放| 国产精品一二三区在线看| 国产 亚洲一区二区三区 | 成年女人在线观看亚洲视频 | 内地一区二区视频在线| 成人特级av手机在线观看| 波野结衣二区三区在线| 麻豆乱淫一区二区| 最近最新中文字幕免费大全7| 午夜免费男女啪啪视频观看| 亚洲国产最新在线播放| 国产av在哪里看| 精品久久久久久成人av| 精品久久国产蜜桃| 99re6热这里在线精品视频| 久久韩国三级中文字幕| 国产精品一区二区性色av| 久久精品久久久久久久性| 色综合亚洲欧美另类图片| 美女被艹到高潮喷水动态| 久久99热这里只频精品6学生| 国产熟女欧美一区二区| 一级黄片播放器| 亚洲欧美成人综合另类久久久| 国产精品麻豆人妻色哟哟久久 | 97超碰精品成人国产| 免费观看无遮挡的男女| 身体一侧抽搐| 一级av片app| 免费不卡的大黄色大毛片视频在线观看 | 22中文网久久字幕| 99久国产av精品| 美女国产视频在线观看| 一级毛片aaaaaa免费看小| 亚洲av一区综合| 夜夜看夜夜爽夜夜摸| 国产精品美女特级片免费视频播放器| 免费观看a级毛片全部| 又黄又爽又刺激的免费视频.| 国产精品爽爽va在线观看网站| 色综合色国产| 简卡轻食公司| 黄片wwwwww| 国产成人精品一,二区| 国产综合精华液| 精品国产一区二区三区久久久樱花 | 99久久人妻综合| 日韩av免费高清视频| 国产伦理片在线播放av一区| 国模一区二区三区四区视频| 亚洲av中文av极速乱| 别揉我奶头 嗯啊视频| 久久久精品94久久精品| 国产精品.久久久| 亚洲国产精品成人久久小说| 婷婷色综合大香蕉| 国产淫语在线视频| 国产高清有码在线观看视频| 亚洲欧美精品自产自拍| 午夜激情欧美在线| 蜜臀久久99精品久久宅男| 久久久色成人| 国产精品1区2区在线观看.| 国产精品女同一区二区软件| 中文在线观看免费www的网站| 婷婷色综合大香蕉| 欧美高清性xxxxhd video| 搡女人真爽免费视频火全软件| 欧美高清性xxxxhd video| 国产成人福利小说| 少妇熟女aⅴ在线视频| 婷婷色综合www| 国产极品天堂在线| 欧美成人午夜免费资源| 国产成人免费观看mmmm| 国产色爽女视频免费观看| 成人特级av手机在线观看| 人妻一区二区av| 亚洲欧美日韩卡通动漫| 大陆偷拍与自拍| 亚洲精品日本国产第一区| 婷婷色麻豆天堂久久| 中文字幕免费在线视频6|