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

    基于AE-BP模型的楊木膠合板應(yīng)力損傷識(shí)別

    2024-07-02 00:00:00劉佳于孟言高珊陳昱龍馮蔓萱杜鑫宇

    摘 要:【目的】利用聲發(fā)射(AE)技術(shù)對(duì)應(yīng)力損傷全過程中的楊木膠合板進(jìn)行無損檢測,并利用BP神經(jīng)網(wǎng)絡(luò)對(duì)AE檢測結(jié)果進(jìn)行識(shí)別,以提高膠合板損傷檢測精度?!痉椒ā恳允袌稣加辛枯^高的托盤用楊木膠合板作為研究對(duì)象,在聯(lián)合AE和應(yīng)力損傷試驗(yàn)過程中,提取6個(gè)AE特征參數(shù),利用聲發(fā)射RA-AF聯(lián)合分析法區(qū)分楊木膠合板產(chǎn)生裂紋的類型,采用K-均值聚類分析方法確定損傷演化程度和AE特征參數(shù)之間的對(duì)應(yīng)關(guān)系,利用BP神經(jīng)網(wǎng)絡(luò)建立損傷識(shí)別模型,并對(duì)識(shí)別網(wǎng)絡(luò)進(jìn)行測試訓(xùn)練。【結(jié)果】AE信號(hào)幅度和上升時(shí)間可有效地表征楊木膠合板應(yīng)力損傷從微裂紋萌生、產(chǎn)生宏觀裂紋至完全斷裂的損傷演化過程;通過RA-AF聯(lián)合分析發(fā)現(xiàn):楊木膠合板在應(yīng)力損傷試驗(yàn)第一階段主要為剪切破壞損傷,第二、三階段主要為拉伸破壞損傷;通過K-均值聚類分析發(fā)現(xiàn)損傷類型與AE峰值頻率之間的存在較強(qiáng)對(duì)應(yīng)關(guān)系,可有效的表征不同的損傷類型:在31 kHz內(nèi)為基體開裂,在31~100 kHz內(nèi)為脫膠分層,大于100 kHz為纖維斷裂;構(gòu)建AE-BP神經(jīng)網(wǎng)絡(luò)模型對(duì)應(yīng)力損傷類型訓(xùn)練樣本的擬合優(yōu)度是95.94%,測試集的擬合優(yōu)度是98.89%,模型總擬合優(yōu)度是96.51%,網(wǎng)絡(luò)訓(xùn)練效果較好。【結(jié)論】在應(yīng)力承載AE監(jiān)測過程中,通過構(gòu)建AE-BP模型,可對(duì)楊木膠合板產(chǎn)生的未知損傷進(jìn)行有效檢測并準(zhǔn)確識(shí)別。

    關(guān)鍵詞:楊木膠合板;聲發(fā)射;BP神經(jīng)網(wǎng)絡(luò);損傷識(shí)別

    中圖分類號(hào):S781.3 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1673-923X(2024)04-0169-11

    基金項(xiàng)目:國家自然科學(xué)基金面上項(xiàng)目(32071685);國家自然科學(xué)基金青年項(xiàng)目(31600453)。

    Identification of stress damage in poplar plywood based on AE-BP model

    LIU Jia, YU Mengyan, GAO Shan, CHEN Yulong, FENG Manxuan, DU Xinyu

    (College of Civil Engineering and Transportation, Northeast Forestry University, Harbin 150040, Heilongjiang, China)

    Abstract:【Objective】Acoustic emission (AE) is used to detect the damage of poplar plywood in the whole process of stress damage, and BP neural network is applied to identify the results of AE, so as to improve the damage detection accuracy of plywood.【Method】Poplar plywood used in pallets with high market share was taken as the research object. During the joint AE and stress damage test, six AE characteristic parameters were extracted, the crack types of plywood were distinguished by acoustic emission RAAF joint analysis method, and the corresponding relationship between damage evolution degree and AE characteristic parameters was determined by K-means clustering analysis method. The damage identification model was established by BP neural network, and the identification network was trained by test.【Result】AE signal amplitude and rise time effectively characterized the evolution process of stress damage from microcrack initiation, macroscopic crack to complete fracture; Through RA-AF analysis, it was found that in the first stage of bending test, the main damage of poplar plywood is shear failure. In the second and third stages, the main damage was tensile failure; Based on cluster analysis, it was found that there was a strong corresponding relationship between damage types and AE Peak frequencies,the different damage patterns could be effectively characterized: matrix cracking within 31 kHz, debond and delamination within 31-100 kHz, and fiber fracture above 100 kHz; The AE-BP neural network model could identify the fuzzy damage types with the goodness of fit of the training samples was 95.94%, the goodness of fit of the test set was 98.89%, and the total goodness of fit of the model was 96.51%, and the network training was more effective.【Conclusion】The damage types of poplar plywood during AE monitoring can be effectively detected and accurately identified by constructing AE-BP model.

    Keywords: poplar plywood; acoustic emission; BP neural network; damage identification

    膠合板具有強(qiáng)度高、剛度大、穩(wěn)定性好、加工性能好、耐磨性好、耐腐蝕性好、使用壽命長等優(yōu)點(diǎn),廣泛用于家具、建筑、車船制造業(yè)及包裝等領(lǐng)域[1,2]。在包裝領(lǐng)域,膠合板主要用于制作木箱、木托盤等。其中膠合板托盤,具有成本低、生產(chǎn)工藝簡單、出口無熏蒸等優(yōu)點(diǎn),是替代實(shí)木托盤的最佳選擇。托盤在使用過程中因外載荷的變化導(dǎo)致內(nèi)部出現(xiàn)應(yīng)力集中而產(chǎn)生裂痕,造成的材料局部開裂會(huì)縮短木質(zhì)托盤的壽命,嚴(yán)重時(shí)會(huì)造成經(jīng)濟(jì)損失[3,4]。依照木質(zhì)平托盤通用技術(shù)要求國家標(biāo)準(zhǔn)(2022)[5],木質(zhì)托盤性能檢測主要包括整體抗壓實(shí)驗(yàn)、鋪板抗壓實(shí)驗(yàn)等破壞性實(shí)驗(yàn),無法實(shí)時(shí)監(jiān)測托盤鋪板內(nèi)部的形變和開裂,不適用于生產(chǎn)線上的連續(xù)快速檢測。為了科學(xué)合理地使用膠合板托盤,提高其生產(chǎn)質(zhì)量,擴(kuò)大其適用范圍,需掌握膠合板的力學(xué)性能及實(shí)時(shí)在線檢測技術(shù)。了解托盤用膠合板的性能和性質(zhì)對(duì)于提高整個(gè)托盤結(jié)構(gòu)的耐久性和穩(wěn)定性有著重要的作用。因此,發(fā)展一種有效的膠合板托盤在線無損檢測方法具有十分重要的意義。

    聲發(fā)射(Acoustic Emission,AE)可提供應(yīng)力作用下的木質(zhì)材料缺陷動(dòng)態(tài)信息,實(shí)現(xiàn)在線實(shí)時(shí)檢測,且對(duì)檢測材料形狀、大小等無特定要求等優(yōu)點(diǎn)[6],應(yīng)用領(lǐng)域比其他檢測方法具有更廣泛的良好前景,已廣泛應(yīng)用于人造板領(lǐng)域。Brunner等[7]使用參數(shù)分析法對(duì)兩種不同類型的層壓木材樣品在準(zhǔn)靜態(tài)拉伸測試過程中產(chǎn)生的AE信號(hào)數(shù)據(jù)進(jìn)行分析;郭曉磊等[8]在三點(diǎn)彎曲載荷過程中利用AE特征參數(shù)累計(jì)能量和累計(jì)振鈴計(jì)數(shù)分析了木基復(fù)合材料損傷過程及特點(diǎn);Niemz等[9]應(yīng)用AE信號(hào)能量、幅度等參數(shù)來研究膠合板的斷裂行為,發(fā)現(xiàn)非均質(zhì)材料(顆粒膠合板)比均質(zhì)的中密度纖維板更容易預(yù)測宏觀斷裂位置;Ritschel等[10]通過AE和數(shù)字圖像相關(guān)(DIC)方法研究了云杉制成的膠合板在拉伸載荷過程中的亞宏觀損傷演變,發(fā)現(xiàn)AE參數(shù)能量與膠合板損傷內(nèi)部有一致對(duì)應(yīng)關(guān)系;徐峰等[11]基于EMD-SVD的特征提取方法對(duì)膠合板損傷進(jìn)行了識(shí)別;Wang等[12]通過單參數(shù)分析和K均值聚類分析,在四點(diǎn)加載試驗(yàn)中比較了竹木、竹膠合板、SPF木材的力學(xué)性能及AE活動(dòng);Clerc等[13]將無監(jiān)督模式識(shí)別法應(yīng)用于山毛櫸木膠合木的AE信號(hào)識(shí)別,發(fā)現(xiàn)AE信號(hào)可區(qū)分木材斷裂和膠合界面斷裂。

    綜合上述研究可見,AE參數(shù)分析法在木基材料損傷研究中具有很強(qiáng)的應(yīng)用性,但僅能對(duì)加載損傷過程進(jìn)行階段分析,近年來學(xué)者們開始將聚類分析和人工神經(jīng)網(wǎng)絡(luò)等模式識(shí)別法與AE技術(shù)結(jié)合起來[14-18],研究多用于研究混凝土、巖石等損傷識(shí)別,用于人造板研究較為鮮見。

    本研究以市場上占有量較高的楊木膠合板為研究對(duì)象,以三點(diǎn)彎曲加載實(shí)驗(yàn)?zāi)M實(shí)際應(yīng)力損傷產(chǎn)生過程,探索聲發(fā)射特征參數(shù)對(duì)膠合板內(nèi)部基體開裂、脫膠分層、纖維斷裂不同損傷階段和類型的響應(yīng)規(guī)律,運(yùn)用參數(shù)分析法、RA-AF聯(lián)合法和K-means聚類法,對(duì)AE檢測結(jié)果進(jìn)行分析,并結(jié)合BP神經(jīng)網(wǎng)絡(luò)建立應(yīng)力損傷識(shí)別模型,實(shí)現(xiàn)對(duì)損傷演化過程檢測和損傷類型識(shí)別,為實(shí)際應(yīng)用中的膠合板托盤鋪板損傷在線無損檢測提供基礎(chǔ)理論參考。

    1 材料與方法

    1.1 試驗(yàn)材料及設(shè)備

    1)試驗(yàn)材料。楊木膠合板樣品采用3層膠合板,其內(nèi)部為楊木夾芯,厚度為3 mm,外覆桃花芯面板,厚度為1 mm,由環(huán)保脲醛膠粘合而成。尺寸選擇膠合板托盤鋪板國標(biāo)尺寸1/4比例[19],即300 mm×250 mm×5 mm,試件數(shù)為35,表面無視覺缺陷。測試溫度為常溫, 樣品含水率為氣干狀態(tài)(圖1)。

    2)試驗(yàn)設(shè)備。阿爾泰科技有限公司USB-3133A高速采集卡、LabVIEW軟件、PXR15I單端諧振聲發(fā)射傳感器及PXDBlock信號(hào)分離器(長沙鵬翔電子科技有限公司),兩通道聲發(fā)射信號(hào)采集系統(tǒng);PXR15I單端諧振聲發(fā)射傳感器(長沙鵬翔電子科技有限公司),頻率帶寬為100~400 kHz,諧振頻率為150 kHz,內(nèi)置前置放大器增益為40 dB;PXD Block信號(hào)分離器(長沙鵬翔電子科技有限公司),用來給傳感器內(nèi)置放大器模塊電路供電。根據(jù)香農(nóng)采樣定理,設(shè)置每個(gè)通道的采樣頻率為500 kHz,輸出電壓范圍為(-5 V,5 V);深圳市新三思計(jì)量技術(shù)有限公司生產(chǎn)的SANS-CMT5504力學(xué)試驗(yàn)機(jī)。試驗(yàn)裝置見圖2。

    1.2 AE聯(lián)合三點(diǎn)彎曲荷載試驗(yàn)

    將兩個(gè)AE傳感器安裝于托盤用楊木膠合板距離中心位置30 mm處左右兩側(cè),然后在試件與傳感器接觸位置涂抹凡士林耦合劑使二者充分耦合接觸,用絕緣膠帶將傳感器固定在鋪板表面。支點(diǎn)跨距為200 mm,進(jìn)給速度為5 mm/min。首先,將試件加載至預(yù)估極限荷載的10%,持續(xù)2 min,然后卸載2 min,清零所有數(shù)值,正式加載;設(shè)備啟動(dòng)后持續(xù)加壓,直至托盤用楊木膠合板試件斷裂為止。當(dāng)膠合板破裂時(shí),停止試驗(yàn)并記錄該點(diǎn)的極限載荷[20]。施加荷載過程中AE同步動(dòng)態(tài)跟蹤監(jiān)測托盤用楊木膠合板的損傷演化全過程,載荷快速下降時(shí)停止加載,同時(shí)停止AE信號(hào)采集。

    2 結(jié)果與分析

    2.1 楊木膠合板應(yīng)力損傷演化過程AE特征參數(shù)分析

    托盤用楊木膠合板應(yīng)力損傷破壞形態(tài)圖見圖3。將能量、振鈴計(jì)數(shù)、幅度、持續(xù)時(shí)間、上升時(shí)間、峰值頻率6個(gè)AE特征參數(shù)及累計(jì)參數(shù)歷程圖與荷載曲線關(guān)聯(lián)起來對(duì)于整個(gè)實(shí)驗(yàn)過程進(jìn)行分析,見圖4。根據(jù)試件應(yīng)力損傷破壞過程的AE特征參數(shù)損傷演化可以分為三個(gè)階段:微裂紋萌生、產(chǎn)生宏觀裂紋和完全斷裂。研究發(fā)現(xiàn)所有測定試樣的AE特征參數(shù)隨荷載變化的規(guī)律基本一致(圖4)。

    1)微裂紋萌生階段:0~30 s,荷載隨著加載時(shí)間以一定的斜率呈上升狀態(tài),AE信號(hào)在試件開始受力后立即產(chǎn)生,該階段產(chǎn)生的AE信號(hào)微弱且不連續(xù)(幅度<45 dB,振鈴計(jì)數(shù)<10,能量<170 mV·mS,上升時(shí)間<1 200 μs,持續(xù)時(shí)間<1 800 μs),峰值頻率主要分布在0~20 kHz和125 kHz范圍內(nèi),6個(gè)累計(jì)參數(shù)曲線持續(xù)緩慢上升。AE信號(hào)的產(chǎn)生是由加載過程中的機(jī)械噪聲以及試件內(nèi)部產(chǎn)生原始微損傷并持續(xù)產(chǎn)生新的損傷引起的,在此階段托盤用膠合板表面無肉眼可見裂紋。

    2)產(chǎn)生宏觀裂紋階段:30~59 s,隨著荷載的增加,此階段振鈴計(jì)數(shù)、能量、持續(xù)時(shí)間多次出現(xiàn)較高峰值,峰值頻率中低頻信號(hào)明顯增加,累計(jì)振鈴計(jì)數(shù)、累計(jì)能量、累計(jì)上升時(shí)間、累計(jì)持續(xù)時(shí)間曲線總體變化趨勢一致,在該階段后期斜率增大,同時(shí)幅度、上升時(shí)間達(dá)到最大峰值(幅度達(dá)到60 dB左右,上升時(shí)間達(dá)到3 100 μs左右),聲發(fā)射累計(jì)幅度曲線的斜率逐漸增大,膠合板表面產(chǎn)生肉眼可見的宏觀裂紋,并伴有輕微開裂聲,說明托盤用楊木膠合板內(nèi)部結(jié)構(gòu)正在發(fā)生破壞,內(nèi)部活動(dòng)劇烈,釋放大量的AE信號(hào)。

    3)完全斷裂階段:59~80 s,第三階段中期達(dá)到最大載荷,伴隨較大的響聲,托盤用楊木膠合板纖維持續(xù)拉脫斷開,載荷呈階躍式下降,直至完全斷裂。在完全斷裂時(shí)刻,能量、振鈴計(jì)數(shù)和持續(xù)時(shí)間達(dá)到最大峰值,分別為24 000 mV·mS、550次和66 000 μs左右;峰值頻率高頻信號(hào)占比增加,累計(jì)振鈴計(jì)數(shù)、累計(jì)能量、累計(jì)持續(xù)時(shí)間相關(guān)曲線總體變化趨勢一致,斜率在第三階段前期出現(xiàn)增大現(xiàn)象,該階段試件舊裂紋失穩(wěn)擴(kuò)展并伴有新裂紋的不斷產(chǎn)生,材料發(fā)生破壞后,釋放出的AE信號(hào)逐漸減少。

    選用AE參數(shù)分析法研究了楊木膠合板彎曲損傷實(shí)驗(yàn)全過程,并將損傷演化分為三個(gè)階段。在第二階段,即產(chǎn)生宏觀裂紋階段,幅度、上升時(shí)間達(dá)到峰值,而能量、振鈴計(jì)數(shù)及持續(xù)時(shí)間在第三階段即裂紋失穩(wěn)直至斷裂階段達(dá)到峰值。累計(jì)幅度、累計(jì)上升時(shí)間損傷一、二階段的分界點(diǎn)較清晰,但依舊無法作為損傷階段劃分的精確判定方式,能量、振鈴計(jì)數(shù)、持續(xù)時(shí)間、累計(jì)幅度、累計(jì)振鈴計(jì)數(shù)、累計(jì)能量、累計(jì)持續(xù)時(shí)間、累計(jì)上升時(shí)間難以直觀確定裂紋失穩(wěn)斷裂臨界點(diǎn)。幅度、上升時(shí)間隨時(shí)間變化曲線中,各階段分界點(diǎn)較明顯,可以作為楊木膠合板應(yīng)力損傷各階段的劃分標(biāo)準(zhǔn)。

    2.2 RA-AF關(guān)聯(lián)分析

    聲發(fā)射參數(shù)RA為上升時(shí)間與幅度的比值(μs/V),AF為振鈴計(jì)數(shù)與持續(xù)時(shí)間的比值(kHz),RA-AF關(guān)聯(lián)分析法可有效區(qū)分張拉型裂紋和剪切型裂紋[21]。已有研究表明張拉型裂紋主要以縱波的形式傳播,其引起的裂紋側(cè)向擴(kuò)展具有高頻短時(shí)(高AF)的特征,而剪切型裂紋主要以橫波形式傳播且具有低頻長時(shí)(低AF)的特點(diǎn)[22],因此認(rèn)為高RA和低AF的信號(hào)對(duì)應(yīng)剪切型裂紋,而低RA和高AF的信號(hào)對(duì)應(yīng)張拉型裂紋。

    托盤用楊木膠合板應(yīng)力損傷試驗(yàn)過程中RA、AF以及荷載隨時(shí)間變化的趨勢如圖5所示。第一階段(0~30 s)開始時(shí),AF保持較小的波動(dòng),走勢較為平緩,在30 s左右出現(xiàn)較小峰值,RA在20 s左右出現(xiàn)較大峰值;第二階段(30~59 s)與第一階段相比,AF波動(dòng)范圍明顯增大,RA波動(dòng)范圍明顯減小,AF在第二階段中期出現(xiàn)峰值,RA在第二階段后期出現(xiàn)峰值;第三階段(59~80 s),AF值連續(xù)幾次出現(xiàn)峰值,RA值波動(dòng)范圍較前兩個(gè)階段減小。第一階段,托盤用楊木膠合板主要為剪切破壞,第二、三階段主要為拉伸破壞。

    2.3 AE特征參數(shù)聚類分析

    針對(duì)膠合板的四種類型損傷類型:基體開裂、脫膠、分層、纖維斷裂[23],使用DBi指標(biāo)和SC系數(shù)對(duì)6個(gè)AE特征參數(shù)的檢測結(jié)果進(jìn)行聚類有效性評(píng)價(jià),聚類數(shù)目K的范圍取2~10。其中,SC輪廓系數(shù)越大,DBi指標(biāo)越小,聚類數(shù)越好。圖6a展示了試件在DBi指標(biāo)和SC輪廓系數(shù)下的曲線圖,當(dāng)聚類數(shù)K=3時(shí),DBi指數(shù)達(dá)到最小值,SC系數(shù)達(dá)到最大值,選擇K=3進(jìn)行初步聚類,K=3時(shí)聚類結(jié)果理想,試件的AE信號(hào)最終被聚類為3類。

    托盤用膠合板應(yīng)力損傷實(shí)驗(yàn)過程中聲發(fā)射參數(shù)信號(hào)的K-均值聚類分析(圖6b—6f),可以分為3個(gè)聲發(fā)射信號(hào)類別。A類信號(hào)(Cluster A)的振鈴計(jì)數(shù)分布在0~20次,持續(xù)時(shí)間分布在0~9 000 μs,上升時(shí)間分布在0~300 μs,能量分布在0~2 000 mV·mS ,峰值頻率分布在0~40 kHz。B類信號(hào)(Cluster B)的振鈴計(jì)數(shù)分布在5~70次,持續(xù)時(shí)間分布在9 000~22 000 μs,上升時(shí)間分布在300~1 100 μs,能量分布在2 000~6 200 mV·mS,峰值頻率分布在40~100 kHz,C類信號(hào)(Cluster C)的振鈴計(jì)數(shù)、持續(xù)時(shí)間、上升時(shí)間、能量、峰值頻率分布分別大于70次、22 000 μs、1 100 μs、6 200 mV·mS、100 kHz。

    A類信號(hào)具有低能量、低振鈴計(jì)數(shù)、低持續(xù)時(shí)間、低上升時(shí)間、低峰值頻率的特征,幅度分布為35~62 dB,托盤用楊木膠合板發(fā)生損傷時(shí)基體開裂最先出現(xiàn),并且覆蓋整個(gè)應(yīng)力損傷過程,結(jié)合圖4聲發(fā)射參數(shù)圖分析,發(fā)現(xiàn)A類信號(hào)在應(yīng)力損傷試驗(yàn)過程中一直出現(xiàn),因此A類信號(hào)認(rèn)為是基體開裂信號(hào);C類信號(hào)具有高能量、高振鈴計(jì)數(shù)、高持續(xù)時(shí)間、高上升時(shí)間、高峰值頻率的特征,結(jié)合圖4聲發(fā)射參數(shù)圖分析,C類信號(hào)認(rèn)為是托盤用楊木膠合板纖維斷裂信號(hào);B類信號(hào)都為中頻信號(hào),B類信號(hào)認(rèn)為是托盤用膠合板脫膠或分層信號(hào)。觀察圖6發(fā)現(xiàn)峰值頻率很好地把3類聲發(fā)射信號(hào)區(qū)分開,認(rèn)為是表征試件不同損傷模式效果最佳的聲發(fā)射參數(shù)。應(yīng)力損傷試驗(yàn)全過程中托盤用楊木膠合板各損傷階段對(duì)應(yīng)AE特征參數(shù)值見表1。

    3 AE-BP模型的損傷識(shí)別

    3.1 試驗(yàn)準(zhǔn)備及激活函數(shù)的選擇

    BP神經(jīng)網(wǎng)絡(luò)是經(jīng)過反復(fù)迭代訓(xùn)練形成的具有識(shí)別與預(yù)測能力的神經(jīng)網(wǎng)絡(luò)模型[24-25],拓?fù)浣Y(jié)構(gòu)如圖7所示。

    圖7中,輸入層參數(shù)(m)為振鈴計(jì)數(shù)、能量、絕對(duì)能量、幅度、上升計(jì)數(shù)、上升時(shí)間、平均頻率、峰值頻率、持續(xù)時(shí)間9個(gè)AE參數(shù)值,設(shè)定輸出層參數(shù)(n)為膠合板托盤鋪板的損傷類型:基體開裂、脫膠分層和纖維斷裂,即m=9,n=3。本研究中隱含層節(jié)點(diǎn)數(shù)量L采用下式(1)進(jìn)行初步確定,a為0~10的調(diào)節(jié)常數(shù),以網(wǎng)絡(luò)訓(xùn)練的均方根誤差最小為指標(biāo),確定L=10。

    隱含層、輸出層激勵(lì)函數(shù)分別選取S型正切函數(shù)Tansig、Purelin函數(shù),并采用Trainrp算法訓(xùn)練網(wǎng)絡(luò),性能函數(shù)選用均方誤差函數(shù),學(xué)習(xí)速率為0.01,最大迭代次數(shù)設(shè)為1 000次,目標(biāo)誤差為0.001。

    3.2 識(shí)別結(jié)果

    利用BP神經(jīng)網(wǎng)絡(luò)建立AE信號(hào)與損傷類型的對(duì)應(yīng)關(guān)系,將AE特征參數(shù)代入模型,得到相應(yīng)的損傷類型,并與實(shí)際值相對(duì)比。選取編號(hào)為3的托盤用楊木膠合板AE參數(shù)數(shù)據(jù),前199組數(shù)據(jù)進(jìn)行網(wǎng)絡(luò)訓(xùn)練,最后70組數(shù)據(jù)進(jìn)行測試,將測試結(jié)果和真實(shí)值相對(duì)比。表2展示了其中10組有代表性的BP神經(jīng)網(wǎng)絡(luò)輸出數(shù)據(jù)。

    圖8a—c為輸出層3個(gè)節(jié)點(diǎn)訓(xùn)練集與預(yù)測集結(jié)果對(duì)比圖,圖中均方根誤差RMSE都小于0.2,該指標(biāo)越小,表示模型的預(yù)測能力越強(qiáng)。圖8d為 3層BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練過程中的收斂曲線圖,隨著訓(xùn)練次數(shù)的增加,均方差逐漸降低且趨于極小值,第11次網(wǎng)絡(luò)誤差為1.32×10-5,達(dá)到了網(wǎng)絡(luò)誤差設(shè)置的要求,發(fā)現(xiàn)訓(xùn)練后BP神經(jīng)網(wǎng)絡(luò)具有較好的識(shí)別效果。圖8e為BP神經(jīng)網(wǎng)絡(luò)預(yù)測值和目標(biāo)值的擬合優(yōu)度圖,橫坐標(biāo)為樣本原目標(biāo)值,縱坐標(biāo)為網(wǎng)絡(luò)輸出預(yù)測值,二者相關(guān)度用系數(shù)R表示,若R越接近1,則表示線性化程度越高,結(jié)果越好。訓(xùn)練樣本的擬合優(yōu)度是95.94%,測試集的擬合優(yōu)度是98.89%,模型總擬合優(yōu)度是96.51%,網(wǎng)絡(luò)訓(xùn)練效果較好,表明AE聯(lián)合BP神經(jīng)網(wǎng)絡(luò)模型可以精準(zhǔn)有效的進(jìn)行膠合板損傷檢測與識(shí)別。

    4 結(jié)論與討論

    4.1 結(jié) 論

    本研究基于AE技術(shù)研究了楊木膠合板不同損傷程度下AE參數(shù)特征值變化規(guī)律,并利用AE參數(shù)和BP神經(jīng)網(wǎng)絡(luò)建立了損傷類型識(shí)別模型,得出以下結(jié)論:

    1)通過AE特征參數(shù)變化可有效的判斷楊木膠合板微裂紋萌生、產(chǎn)生宏觀裂紋和完全斷裂三個(gè)階段。通過RA-AF聯(lián)合分析發(fā)現(xiàn)拉伸破壞損傷是楊木膠合板的損傷的主要形式。

    2)通過K-均值聚類分析發(fā)現(xiàn)楊木膠合板AE峰值頻率與損傷類型之間的對(duì)應(yīng)關(guān)系較好,對(duì)損傷模式區(qū)分表征效果最佳。

    3)聯(lián)合AE-BP神經(jīng)網(wǎng)絡(luò)模型可準(zhǔn)確有效地對(duì)楊木膠合板未知損傷類型進(jìn)行有效檢測與識(shí)別。

    4.2 討 論

    楊木膠合板試件在抗彎試驗(yàn)加載過程中,會(huì)逐漸出現(xiàn)多種損傷以及對(duì)應(yīng)的聲發(fā)射現(xiàn)象。根據(jù)AE基本參數(shù)隨時(shí)間的變化圖和荷載位移曲線,可將試驗(yàn)分為三個(gè)階段,即微裂紋產(chǎn)生階段、產(chǎn)生宏觀裂紋階段和完全斷裂階段,這與涂郡成等[26]根據(jù)木梁尖端損傷演變試驗(yàn)劃分階段類似。其中材料損傷程度越嚴(yán)重,聲發(fā)射活動(dòng)越劇烈。AE信號(hào)幅度、上升時(shí)間隨時(shí)間變化曲線中,各階段分界點(diǎn)較明顯,可以作為托盤用楊木膠合板應(yīng)力損傷各階段的劃分標(biāo)準(zhǔn)(圖4)。在之前的研究中,郭曉磊等[8]研究了膠合板(11層,相鄰單板層木紋方向相互垂直)和三聚氰胺貼面刨花板兩種材料,發(fā)現(xiàn)聲發(fā)射信號(hào)的幅值、累計(jì)能量、振鈴計(jì)數(shù)等參數(shù)可有效表征兩種材料內(nèi)部損傷演化過程。楊洲凌等[27]研究發(fā)現(xiàn)櫸木三點(diǎn)彎曲試驗(yàn)中所產(chǎn)生的累計(jì)AE振鈴計(jì)數(shù)是反映櫸木損傷的重要特征參數(shù)。猜測選取不同AE特征參數(shù)的原因是研究對(duì)象不同,AE特征參數(shù)對(duì)于不同材料的敏感性不同。

    本研究RA-AF關(guān)聯(lián)分析法可有效區(qū)分楊木膠合板張拉型裂紋和剪切型裂紋,而RA-AF關(guān)聯(lián)分析方法在該方面的應(yīng)用較少。王婷[22]曾利用RA-AF方法來區(qū)分SPF、楊木和杉木壓縮剪切和彎曲剪切試件的張拉型裂紋與剪切型裂紋,結(jié)果發(fā)現(xiàn),低幅度的信號(hào)簇以張拉型裂紋為主,中高幅度的信號(hào)簇則為張剪復(fù)合破壞。本研究在此基礎(chǔ)上將該方法用于托盤用楊木膠合板的張拉型裂紋與剪切型裂紋的區(qū)分中,拉伸破壞占楊木膠合板托盤鋪板的60%以上,在材料破壞中起主導(dǎo)作用,因此在使用托盤用楊木膠合板時(shí)更應(yīng)關(guān)注其抗拉伸性能。

    經(jīng)K-均值處理后的楊木膠合板AE各簇信號(hào)與損傷類型之間的對(duì)應(yīng)關(guān)系,這個(gè)結(jié)果與目前研究[16,28]得出的結(jié)果較為吻合,即低峰值頻率信號(hào)對(duì)應(yīng)基體開裂,中峰值頻率信號(hào)對(duì)應(yīng)脫膠分層,高峰值頻率信號(hào)對(duì)應(yīng)纖維斷裂。其中峰值頻率表征托盤用膠合板不同損傷模式效果最佳,與曹枝軍等[29]對(duì)熱障涂層三點(diǎn)彎曲實(shí)驗(yàn)表征失效模式的AE參數(shù)的選擇結(jié)果具有一致性。

    本研究利用AE設(shè)備全程監(jiān)測楊木膠合板應(yīng)力損傷試驗(yàn)過程,基于AE技術(shù)和BP神經(jīng)網(wǎng)絡(luò)建立的損傷識(shí)別模型對(duì)于膠合板的應(yīng)力損傷檢測效果優(yōu)異,可為楊木膠合板托盤鋪板在實(shí)際使用過程中的損傷狀況提供建議與參考,也可用于托盤用楊木膠合板的實(shí)時(shí)監(jiān)測和養(yǎng)護(hù),防止膠合板失效導(dǎo)致經(jīng)濟(jì)損傷。

    但本研究仍存在以下局限性:

    1)本研究環(huán)境溫度和膠合板的水分含量單一,未考慮溫度和水分含量是否會(huì)影響膠合板損傷演化過程及過程中AE參數(shù)的變化。

    2)本研究是以膠合板三點(diǎn)彎曲加載實(shí)驗(yàn)?zāi)M實(shí)際應(yīng)力損傷產(chǎn)生過程,而膠合板在實(shí)際應(yīng)用中損傷失效模式情形復(fù)雜,之后應(yīng)針對(duì)實(shí)際應(yīng)用中應(yīng)力損傷、跌落損傷等不同方式損傷展開研究,并逐步嘗試將AE技術(shù)直接用于膠合板托盤或其他膠合板結(jié)構(gòu)的現(xiàn)場監(jiān)測。

    3)本研究后期可結(jié)合先進(jìn)的數(shù)字圖像處理技術(shù)(DIC)、計(jì)算機(jī)斷層掃描(CT)和掃描電子顯微鏡技術(shù)(SEM)開展膠合板多尺度損傷檢測,將AE技術(shù)與微觀檢測技術(shù)結(jié)合起來,實(shí)現(xiàn)從表面到內(nèi)部監(jiān)測膠合板損傷演變過程。

    參考文獻(xiàn):

    [1] 易錦,祝彬,賀國京,等.CFRP增強(qiáng)膠合木梁剛度退化試驗(yàn)研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2021,41(4):165-172. YI J, ZHU B, HE G J, et al. Experimental study on stiffness degradation of CFRP reinforced glulam beams[J]. Journal of Central South University of Forestry Technology,2021,41(4): 165-172.

    [2] 易錦,楊帆,梁承昊,等.基于剛度退化的落葉松膠合木梁疲勞壽命預(yù)測[J].中南林業(yè)科技大學(xué)報(bào),2022,42(6):135-143,151. YI J, YANG F, LIANG C H, et al. Fatigue life prediction of larch glulam beams based on stiffness degradation[J]. Journal of Central South University of Forestry Technology,2022,42(6): 135-143,151.

    [3] 鞏桂芬,劉隨強(qiáng),田洋洋.鋼釘數(shù)量對(duì)托盤彎曲力學(xué)性能的影響[J].包裝工程,2022,43(1):141-147. GONG G F, LIU S Q, TIAN Y Y. Effects of the number of steel nails on bending mechanical properties of pallets[J]. Packaging Engineering,2022,43(1):141-147.

    [4] 高珊,周欣宇,高歌,等.基于壓電效應(yīng)的木質(zhì)托盤跌落損傷無損檢測[J].中南林業(yè)科技大學(xué)報(bào),2022,42(12):164-174. GAO S, ZHOU X Y, GAO G, et al. Nondestructive detection for falling damage of wooden pallet based on piezoelectric effect[J]. Journal of Central South University of Forestry Technology, 2022,42(12):164-174.

    [5] 國家市場監(jiān)督管理總局.GB/T 31148—2022,木質(zhì)平托盤-通用技術(shù)要求[S].北京:中國標(biāo)準(zhǔn)出版社,2022. State Administration for Market Regulation. GB/T 31148-2022, Wooden flat pallet: General technical requirement[S]. Beijing: China Standards Press,2022.

    [6] NASIR V, AYANLEYE S, KAZEMIRAD S, et al. Acoustic emission monitoring of wood materials and timber structures: a critical review[J]. Construction and Building Materials,2022,350: 128877.

    [7] BRUNNER A J, HOWALD M T, NIEMZ P. Acoustic emission rate behavior of laminated wood specimens under tensile loading[J]. Progress in Acoustic Emission,2006,24:104-110.

    [8] 郭曉磊,郭勇,胡偉,等.木基復(fù)合材料破壞過程中聲發(fā)射特性的研究[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,35(3): 97-100. GUO X L, GUO Y, HU W, et al. Acoustic emission characterization of wood base composite materials fracture[J]. Journal of Nanjing Forestry University (Natural Sciences Edition), 2011,35(3):97-100.

    [9] NIEMZ P, BRUNNER A J, WALTER O. Investigation of the mechanism of failure behaviour of wood based materials using acoustic emission analysis and image processing[J]. Wood Research,2009,54(2):49-62.

    [10] RITSCHEL F, ZHOU Y, BRUNNER A J, et al. Acoustic emission analysis of industrial plywood materials exposed to destructive tensile load[J]. Wood Science Technology,2014,48:611-631.

    [11] 徐鋒,劉云飛.基于EMD-SVD的聲發(fā)射信號(hào)特征提取及分類方法[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2014,22(6):1238-1247. XU F, LIU Y F. Feature extraction and classification method of acoustic emission signal based on EMD-SVD[J]. Journal of Basic Science and Engineering,2014,22(6):1238-1247.

    [12] WANG T, WANG Z Q, YANG Y, et al. Acoustic emission characteristics of different bamboo and wood materials in bending failure process[J]. Journal of Renewable Materials, 2022,10(2):527-540.

    [13] CLERC G, SAUSE M G R, BRUNNER A J, et al. Fractography combined with unsupervised pattern recognition of acoustic emission signals for a better understanding of crack propagation in adhesively bonded wood[J]. Wood Science and Technology, 2019,53(6):1235-1253.

    [14] LAXMAN K C, ROSS A, AI L, et al. Determination of vehicle loads on bridges by acoustic emission and an improved ensemble artificial neural network[J]. Construction and Building Materials,2023,364:129844.

    [15] LIMA R A A, DROBIAZKO M, BERNASCONI A, et al. On crack tip localisation in quasi-statically loaded, adhesively bonded double cantilever beam specimens by acoustic emission[J]. Theoretical and Applied Fracture Mechanics,2022,118:103286.

    [16] HOSEINLAGHAB S, FARAHANI M, SAFARABADI M, et al. Tension-after-impact analysis and damage mechanism evaluation in laminated composites using AE monitoring[J]. Mechanical Systems and Signal Processing,2023,186:109844.

    [17] ZHOU L, WANG P Z, ZHANG C J, et al. Multi-mode fusion BP neural network model with vibration and acoustic emission signals for process pipeline crack location[J]. Ocean Engineering,2022,264:112384.

    [18] KARIMIAN S F, MODARRES M. Acoustic emission signal clustering in cfrp laminates using a new feature set based on waveform analysis and information entropy analysis[J]. Composite Structures,2021,268:113987.

    [19] 國家市場監(jiān)督管理總局.GB/T 2934—2007,聯(lián)運(yùn)通用平托盤主要尺寸及公差[S].北京:中國標(biāo)準(zhǔn)出版社,2007. State Administration for Market Regulation. GB/T 2934-2007, General-purpose flat pallets for through transit of goods: Principal dimensions and tolerances[S]. Beijing: China Standards Press,2007.

    [20] 邢碧瀅,龐燕.小徑木通用平托盤力學(xué)性能檢測實(shí)驗(yàn)研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2019,39(8):131-138. XING B Y, PANG Y. Experimental research on mechanical properties detectionof small diameter wood universal flat pallet[J]. Journal of Central South University of Forestry Technology, 2019,39(8):131-138.

    [21] WANG S Q, FU L X, CHEN X H, et al. Characteristic analysis of cement grouted asphalt mixture cracking based on acoustic emission[J]. Construction and Building Materials,2023,375: 130927.

    [22] 王婷.平面內(nèi)外加載下木質(zhì)材料滾動(dòng)剪切聲發(fā)射特性研究[D].南京:南京林業(yè)大學(xué),2022. WANG T. Acoustic emission characteristics of rolling shear failure of wood materials under in-plane and out-plane shear loads[D]. Nanjing: Nanjing Forestry University,2022.

    [23] 徐鋒,劉云飛.基于EMD的膠合板損傷聲發(fā)射信號(hào)特征提取及神經(jīng)網(wǎng)絡(luò)模式識(shí)別[J].振動(dòng)與沖擊,2012,31(15):30-35. XU F, LIU Y F. Feature extraction and pattern recognition of acoustic emission signals generated from plywood damage based on EMD and neural network[J]. Journal of Vibration and Shock,2012,31(15): 30-35.

    [24] 趙奎,宋林,曾鵬,等.單軸壓縮條件下充填體—砂巖組合體破壞過程及前兆特征聲發(fā)射試驗(yàn)研究[J].金屬礦山, 2022(5):70-76. ZHAO K, SONG L, ZENG P, et al. Experimental study on failure process and precursor characteristic acoustic emission of filling body-sandstone combination under uniaxial compression[J]. Metal Mine,2022(5):70-76.

    [25] 張美林,李俊萩,張晴暉,等.基于熵和波形特征的木材損傷斷裂過程聲發(fā)射信號(hào)處理[J].林業(yè)工程學(xué)報(bào),2022,7(2): 159-166.ZHANG M L, LI J Q, ZHANG Q H, et al. Acoustic emission signal processing and analysis of wood damage and fracture process based on entropy and waveform characteristics[J]. Journal of Forestry Engineering,2022,7(2):159-166.

    [26] 涂郡成,趙東,趙健.應(yīng)用AE和DIC原位監(jiān)測含橫紋裂紋木構(gòu)件的裂紋演化規(guī)律試驗(yàn)研究[J].北京林業(yè)大學(xué)學(xué)報(bào), 2020,42(1):142-148. TU J C, ZHAO D, ZHAO J. Experimental study on in situ monitoring of the evolution law of cracks in wood components with transverse cracks based on acoustic emission and image correlation[J]. Journal of Beijing Forestry University,2020,42(1): 142-148.

    [27] 楊洲凌,黃長林,王明華,等.櫸木橫紋三點(diǎn)彎曲受壓聲發(fā)射損傷模型研究[J].西北林學(xué)院學(xué)報(bào),2022,37(4):238-242,256. YANG Z L, HUANG C L, WANG M H, et al. Acoustic emission damage model of Zelkova schneideriana transverse grain compression three-point bending[J]. Journal of Northwest Forestry University,2022,37(4):238-242,256.

    [28] LI L, SWOLFS Y, STRAUMIT I, et al. Cluster analysis of acoustic emission signals for 2D and 3D woven carbon fiber/ epoxy composites[J]. Journal of Composite Materials,2016,50(14): 1921-1935.

    [29] 曹枝軍,袁建輝,蘇懷宇,等.基于原位聲發(fā)射信號(hào)的熱障涂層三點(diǎn)彎曲失效模式的機(jī)器學(xué)習(xí)[J].硅酸鹽學(xué)報(bào),2023,51(2): 373-388. CAO Z J, YUAN J H, SU H Y, et al. Failure mode of thermal barrier coatings based on acoustic emission under three-point bending via machine learning based on in-situ acoustic emission signals[J]. Journal of the Chinese Ceramic Society,2023,51(2): 373-388.

    [本文編校:羅 列]

    午夜激情久久久久久久| 下体分泌物呈黄色| 成人二区视频| 久久毛片免费看一区二区三区| 亚洲成人av在线免费| 天天躁夜夜躁狠狠久久av| 亚洲精品aⅴ在线观看| 人体艺术视频欧美日本| 99视频精品全部免费 在线| 夫妻性生交免费视频一级片| av福利片在线观看| 精品国产露脸久久av麻豆| 日本-黄色视频高清免费观看| 老司机影院毛片| 国产有黄有色有爽视频| 精品国产露脸久久av麻豆| 我的老师免费观看完整版| 中文字幕人妻丝袜制服| 国产精品一区二区在线观看99| 欧美区成人在线视频| 三级国产精品欧美在线观看| 国产成人精品久久久久久| 国产av一区二区精品久久| 少妇精品久久久久久久| 精品国产露脸久久av麻豆| 日韩制服骚丝袜av| 99久久综合免费| 亚洲精品中文字幕在线视频 | 国产成人精品久久久久久| 日本wwww免费看| 久久97久久精品| 另类精品久久| 18+在线观看网站| 精品人妻偷拍中文字幕| 精品午夜福利在线看| 亚洲va在线va天堂va国产| 久久久久久久大尺度免费视频| 国产精品99久久久久久久久| 久久国产乱子免费精品| 精品少妇久久久久久888优播| 三级国产精品片| 国产午夜精品一二区理论片| 一级,二级,三级黄色视频| 亚洲自偷自拍三级| 精品久久久精品久久久| 精品亚洲成a人片在线观看| 国产白丝娇喘喷水9色精品| 69精品国产乱码久久久| 久久午夜综合久久蜜桃| 丰满饥渴人妻一区二区三| 97超碰精品成人国产| 亚洲激情五月婷婷啪啪| 精品一区在线观看国产| 国产日韩欧美视频二区| 亚洲精华国产精华液的使用体验| 免费高清在线观看视频在线观看| 成人漫画全彩无遮挡| 99热全是精品| 国产精品99久久久久久久久| 欧美bdsm另类| 精品国产国语对白av| 乱人伦中国视频| 亚洲av二区三区四区| 免费播放大片免费观看视频在线观看| 毛片一级片免费看久久久久| 蜜桃久久精品国产亚洲av| 日韩熟女老妇一区二区性免费视频| 亚洲欧洲国产日韩| 在线观看免费高清a一片| 九九爱精品视频在线观看| 女性生殖器流出的白浆| 亚洲怡红院男人天堂| 久久午夜福利片| 涩涩av久久男人的天堂| 九九在线视频观看精品| av又黄又爽大尺度在线免费看| 亚洲精品国产av成人精品| 在线观看人妻少妇| 在线看a的网站| 大片电影免费在线观看免费| 少妇 在线观看| 熟女av电影| 免费人妻精品一区二区三区视频| 在线观看美女被高潮喷水网站| 国产无遮挡羞羞视频在线观看| 国产伦在线观看视频一区| 91aial.com中文字幕在线观看| 一级毛片久久久久久久久女| 欧美成人精品欧美一级黄| www.av在线官网国产| 人人妻人人澡人人看| 精品少妇内射三级| 午夜影院在线不卡| 欧美精品人与动牲交sv欧美| 久久ye,这里只有精品| 午夜激情福利司机影院| 午夜福利,免费看| 国产一区二区三区av在线| 亚洲精品久久午夜乱码| 高清毛片免费看| 精品国产露脸久久av麻豆| 亚洲精品乱久久久久久| a 毛片基地| 少妇被粗大的猛进出69影院 | 午夜免费男女啪啪视频观看| 久久99精品国语久久久| 黄色怎么调成土黄色| 国产成人精品一,二区| 色94色欧美一区二区| 色吧在线观看| a 毛片基地| 欧美老熟妇乱子伦牲交| 国产熟女午夜一区二区三区 | av线在线观看网站| 免费观看在线日韩| 国产精品三级大全| 青春草国产在线视频| 国产伦理片在线播放av一区| 欧美bdsm另类| 美女内射精品一级片tv| 乱码一卡2卡4卡精品| 两个人的视频大全免费| 久久综合国产亚洲精品| 亚洲不卡免费看| 熟女电影av网| 黑人猛操日本美女一级片| 国产精品国产三级国产专区5o| 国产一区二区三区综合在线观看 | 在线观看免费日韩欧美大片 | 国产深夜福利视频在线观看| 熟女电影av网| 国产视频首页在线观看| 女的被弄到高潮叫床怎么办| 肉色欧美久久久久久久蜜桃| av一本久久久久| 欧美丝袜亚洲另类| 99精国产麻豆久久婷婷| 国产精品免费大片| 超碰97精品在线观看| 黄色一级大片看看| 国产一区亚洲一区在线观看| 免费在线观看成人毛片| av有码第一页| 国产亚洲91精品色在线| 欧美丝袜亚洲另类| 国产精品久久久久久久电影| av一本久久久久| 美女视频免费永久观看网站| 最近中文字幕2019免费版| 777米奇影视久久| 少妇熟女欧美另类| 久久久精品免费免费高清| 成人亚洲欧美一区二区av| 国产有黄有色有爽视频| 男女啪啪激烈高潮av片| 极品人妻少妇av视频| 青春草亚洲视频在线观看| 在线观看免费高清a一片| a 毛片基地| 国产深夜福利视频在线观看| 免费看不卡的av| 精品国产露脸久久av麻豆| av福利片在线| 国国产精品蜜臀av免费| 性色avwww在线观看| 亚洲精品成人av观看孕妇| 黄色欧美视频在线观看| 在线天堂最新版资源| 欧美少妇被猛烈插入视频| 欧美 日韩 精品 国产| 精品熟女少妇av免费看| 成人特级av手机在线观看| 国产精品蜜桃在线观看| 国产av码专区亚洲av| 欧美另类一区| 久久久久久人妻| 欧美日韩一区二区视频在线观看视频在线| 久久午夜综合久久蜜桃| videos熟女内射| 男女啪啪激烈高潮av片| 精品少妇内射三级| 午夜视频国产福利| 少妇人妻 视频| 曰老女人黄片| 欧美激情极品国产一区二区三区 | av免费观看日本| 亚洲精品国产av成人精品| 欧美高清成人免费视频www| 日韩伦理黄色片| 欧美 亚洲 国产 日韩一| 国产精品不卡视频一区二区| 菩萨蛮人人尽说江南好唐韦庄| 精品一区二区三区视频在线| 精品酒店卫生间| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 日本黄大片高清| 国产精品一区二区在线不卡| 精品亚洲成国产av| 免费人妻精品一区二区三区视频| 国产在视频线精品| 97超碰精品成人国产| 日韩制服骚丝袜av| 韩国av在线不卡| 99精国产麻豆久久婷婷| 成人漫画全彩无遮挡| av在线app专区| 久久99一区二区三区| 成人漫画全彩无遮挡| 精品99又大又爽又粗少妇毛片| av.在线天堂| 免费大片18禁| 免费人成在线观看视频色| 久久国产乱子免费精品| 亚洲国产精品成人久久小说| 黄色日韩在线| 久久久久人妻精品一区果冻| 日本欧美视频一区| 久久久久精品久久久久真实原创| 最近手机中文字幕大全| 国产淫语在线视频| 欧美精品高潮呻吟av久久| 久久青草综合色| 少妇 在线观看| 久久亚洲国产成人精品v| 蜜臀久久99精品久久宅男| 国产综合精华液| a 毛片基地| 亚洲欧美成人精品一区二区| 亚洲熟女精品中文字幕| 国产在视频线精品| av视频免费观看在线观看| 欧美精品一区二区免费开放| 午夜免费鲁丝| 热re99久久精品国产66热6| 女性生殖器流出的白浆| 99久久精品热视频| 精品一区二区三区视频在线| 欧美老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 久久国内精品自在自线图片| 青春草视频在线免费观看| 自拍偷自拍亚洲精品老妇| 亚洲av电影在线观看一区二区三区| 日本黄色片子视频| 亚洲国产日韩一区二区| 我的老师免费观看完整版| 久久人人爽人人爽人人片va| 国产乱来视频区| 国产黄片美女视频| 久久久久视频综合| 高清欧美精品videossex| 国产日韩欧美在线精品| 在线看a的网站| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品成人久久小说| 久热这里只有精品99| 日韩人妻高清精品专区| 伦理电影免费视频| 天堂俺去俺来也www色官网| 人妻一区二区av| 一区二区三区精品91| 久久综合国产亚洲精品| 国产精品99久久99久久久不卡 | 男人和女人高潮做爰伦理| 黄色配什么色好看| 青春草亚洲视频在线观看| 国产日韩欧美在线精品| 久久久久久久久大av| 妹子高潮喷水视频| 91午夜精品亚洲一区二区三区| 能在线免费看毛片的网站| 亚洲,欧美,日韩| 午夜久久久在线观看| 涩涩av久久男人的天堂| 日韩不卡一区二区三区视频在线| 99热这里只有是精品在线观看| 久久99一区二区三区| 好男人视频免费观看在线| 国产探花极品一区二区| 久久久精品94久久精品| 国产免费又黄又爽又色| 97超视频在线观看视频| 十八禁高潮呻吟视频 | 久久久a久久爽久久v久久| 久久精品国产自在天天线| 少妇人妻精品综合一区二区| 午夜福利网站1000一区二区三区| 高清欧美精品videossex| 女性生殖器流出的白浆| 啦啦啦中文免费视频观看日本| 80岁老熟妇乱子伦牲交| 九色成人免费人妻av| 在线观看人妻少妇| 日本vs欧美在线观看视频 | 日韩中文字幕视频在线看片| 国产日韩欧美亚洲二区| 26uuu在线亚洲综合色| 99热全是精品| 亚洲一区二区三区欧美精品| a 毛片基地| 精品熟女少妇av免费看| 夜夜骑夜夜射夜夜干| 久久久国产精品麻豆| 偷拍熟女少妇极品色| 久热这里只有精品99| 大码成人一级视频| av在线老鸭窝| 中国三级夫妇交换| 国产亚洲av片在线观看秒播厂| 国产极品粉嫩免费观看在线 | 老女人水多毛片| 又大又黄又爽视频免费| 香蕉精品网在线| 久久99蜜桃精品久久| 亚洲人与动物交配视频| 久久久亚洲精品成人影院| av不卡在线播放| 国产精品一区二区在线观看99| 最近中文字幕高清免费大全6| 日日摸夜夜添夜夜添av毛片| 日韩强制内射视频| 久热这里只有精品99| 欧美精品国产亚洲| 男人添女人高潮全过程视频| 精品午夜福利在线看| 国产亚洲av片在线观看秒播厂| 亚洲av成人精品一二三区| 亚洲精品456在线播放app| 亚洲精品亚洲一区二区| 秋霞在线观看毛片| 蜜臀久久99精品久久宅男| 狂野欧美白嫩少妇大欣赏| 天天操日日干夜夜撸| 热99国产精品久久久久久7| 成人亚洲欧美一区二区av| 午夜福利视频精品| 成人免费观看视频高清| 卡戴珊不雅视频在线播放| 十八禁网站网址无遮挡 | 免费av不卡在线播放| 亚洲欧美日韩卡通动漫| 国产白丝娇喘喷水9色精品| 亚洲,一卡二卡三卡| 高清视频免费观看一区二区| 天美传媒精品一区二区| 国产欧美亚洲国产| 免费观看无遮挡的男女| 人人妻人人澡人人爽人人夜夜| 嫩草影院新地址| 亚洲国产精品999| 久久影院123| 日韩在线高清观看一区二区三区| 搡老乐熟女国产| 另类亚洲欧美激情| 亚洲国产精品一区三区| 亚洲性久久影院| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 久久久久人妻精品一区果冻| 久久亚洲国产成人精品v| 我要看日韩黄色一级片| 伦精品一区二区三区| 偷拍熟女少妇极品色| 中国美白少妇内射xxxbb| 久久久精品94久久精品| 免费看不卡的av| 成人影院久久| 色5月婷婷丁香| 久久久久人妻精品一区果冻| 欧美激情国产日韩精品一区| 久久6这里有精品| 免费看光身美女| 国产综合精华液| 欧美激情国产日韩精品一区| 高清不卡的av网站| 国产白丝娇喘喷水9色精品| 热re99久久国产66热| 色94色欧美一区二区| 亚洲高清免费不卡视频| 人体艺术视频欧美日本| 老司机影院成人| 边亲边吃奶的免费视频| 午夜激情福利司机影院| 国产精品福利在线免费观看| 女的被弄到高潮叫床怎么办| 精品一区二区免费观看| 久久久久久久亚洲中文字幕| 有码 亚洲区| 少妇的逼水好多| 久久久久人妻精品一区果冻| 精品人妻熟女毛片av久久网站| 久久精品久久精品一区二区三区| 女人久久www免费人成看片| 久久国产亚洲av麻豆专区| 成人免费观看视频高清| 婷婷色综合www| 亚洲国产精品一区二区三区在线| 狂野欧美白嫩少妇大欣赏| 日日啪夜夜撸| 国产精品久久久久久精品古装| 日本av免费视频播放| 国产一级毛片在线| 日韩中字成人| 亚洲图色成人| 免费不卡的大黄色大毛片视频在线观看| videos熟女内射| 在线观看三级黄色| 另类亚洲欧美激情| 丰满迷人的少妇在线观看| 人妻一区二区av| 中国国产av一级| 久久久久久久久久成人| 国产成人免费无遮挡视频| .国产精品久久| 一本大道久久a久久精品| 亚洲欧美一区二区三区黑人 | 精华霜和精华液先用哪个| 午夜精品国产一区二区电影| 桃花免费在线播放| 精品视频人人做人人爽| 9色porny在线观看| 中国国产av一级| 一级爰片在线观看| 五月伊人婷婷丁香| 最近2019中文字幕mv第一页| 肉色欧美久久久久久久蜜桃| 卡戴珊不雅视频在线播放| 久久av网站| 久久人人爽人人爽人人片va| 国产在线一区二区三区精| 午夜av观看不卡| 日韩 亚洲 欧美在线| 久久人妻熟女aⅴ| 久久韩国三级中文字幕| 蜜桃在线观看..| 午夜激情久久久久久久| 精品亚洲乱码少妇综合久久| 日本av免费视频播放| 少妇熟女欧美另类| av网站免费在线观看视频| 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 卡戴珊不雅视频在线播放| 欧美日本中文国产一区发布| 久久午夜综合久久蜜桃| 国产欧美日韩精品一区二区| 美女国产视频在线观看| 一级av片app| 欧美精品一区二区大全| 中文字幕精品免费在线观看视频 | 亚洲成人一二三区av| 午夜日本视频在线| 麻豆成人午夜福利视频| 欧美日韩一区二区视频在线观看视频在线| 日本爱情动作片www.在线观看| 在线播放无遮挡| 久久久久久久精品精品| av天堂中文字幕网| 丝袜脚勾引网站| 欧美亚洲 丝袜 人妻 在线| 欧美3d第一页| 日日撸夜夜添| 少妇的逼水好多| 久久青草综合色| 99视频精品全部免费 在线| 简卡轻食公司| 午夜av观看不卡| 高清在线视频一区二区三区| 天堂中文最新版在线下载| 久久久久网色| 日本与韩国留学比较| 蜜桃久久精品国产亚洲av| www.色视频.com| 亚洲丝袜综合中文字幕| 最近中文字幕高清免费大全6| 最近中文字幕2019免费版| 亚洲欧美中文字幕日韩二区| 哪个播放器可以免费观看大片| 涩涩av久久男人的天堂| 亚洲精品国产色婷婷电影| 黄色配什么色好看| 五月天丁香电影| 国产黄频视频在线观看| 极品少妇高潮喷水抽搐| 精品一区在线观看国产| 日韩熟女老妇一区二区性免费视频| 99久国产av精品国产电影| 免费高清在线观看视频在线观看| 乱人伦中国视频| 日本猛色少妇xxxxx猛交久久| 在线精品无人区一区二区三| 性色avwww在线观看| 日韩av在线免费看完整版不卡| 欧美 日韩 精品 国产| 七月丁香在线播放| 在线播放无遮挡| 大片免费播放器 马上看| 国产精品一区二区性色av| 欧美3d第一页| 搡女人真爽免费视频火全软件| 美女中出高潮动态图| 9色porny在线观看| 欧美激情极品国产一区二区三区 | 在线观看免费高清a一片| 亚洲av日韩在线播放| 另类精品久久| 各种免费的搞黄视频| 中文资源天堂在线| videossex国产| 亚洲精品乱码久久久久久按摩| 中文字幕亚洲精品专区| 高清av免费在线| 国产一区二区三区综合在线观看 | 六月丁香七月| 男女边摸边吃奶| 亚洲美女搞黄在线观看| 欧美另类一区| 亚洲经典国产精华液单| 狂野欧美激情性bbbbbb| 国产精品秋霞免费鲁丝片| 日韩欧美一区视频在线观看 | 寂寞人妻少妇视频99o| 国产视频内射| 97超视频在线观看视频| 欧美人与善性xxx| 中文字幕精品免费在线观看视频 | 欧美精品一区二区大全| 亚洲精品国产av蜜桃| 一级片'在线观看视频| 99re6热这里在线精品视频| 美女cb高潮喷水在线观看| 国产黄频视频在线观看| 国产一区有黄有色的免费视频| 国产综合精华液| 免费观看av网站的网址| 免费看日本二区| xxx大片免费视频| 亚洲在久久综合| 老司机亚洲免费影院| 免费观看a级毛片全部| 国产色婷婷99| 久久97久久精品| 中文天堂在线官网| 九九爱精品视频在线观看| 欧美少妇被猛烈插入视频| 国产精品一区二区在线观看99| 在线免费观看不下载黄p国产| 少妇人妻 视频| 亚洲人成网站在线观看播放| 成人无遮挡网站| 蜜臀久久99精品久久宅男| 91成人精品电影| 天堂8中文在线网| 中文字幕人妻丝袜制服| 国产一级毛片在线| 国产精品免费大片| 亚洲经典国产精华液单| 日韩一区二区三区影片| 各种免费的搞黄视频| 不卡视频在线观看欧美| 在线 av 中文字幕| 最黄视频免费看| 日韩精品免费视频一区二区三区 | 久久久久精品久久久久真实原创| 午夜视频国产福利| 精品久久久精品久久久| 欧美日韩在线观看h| 欧美精品高潮呻吟av久久| 久久精品久久精品一区二区三区| 男的添女的下面高潮视频| 国产深夜福利视频在线观看| 国产精品三级大全| 国产黄片美女视频| 国产乱来视频区| 国产欧美另类精品又又久久亚洲欧美| 久久鲁丝午夜福利片| 黄色视频在线播放观看不卡| 精品一区二区免费观看| 亚洲性久久影院| 成年人午夜在线观看视频| 美女主播在线视频| 久久这里有精品视频免费| av国产久精品久网站免费入址| 亚洲欧美精品专区久久| 精品亚洲成a人片在线观看| 丰满迷人的少妇在线观看| 日本黄大片高清| 蜜桃久久精品国产亚洲av| 久久久久人妻精品一区果冻| 亚洲电影在线观看av| 久热久热在线精品观看| 欧美国产精品一级二级三级 | 一本久久精品| 美女主播在线视频| 午夜免费男女啪啪视频观看| 我要看黄色一级片免费的| 看非洲黑人一级黄片| 两个人免费观看高清视频 | 国内少妇人妻偷人精品xxx网站| 精品亚洲乱码少妇综合久久| 欧美日韩视频高清一区二区三区二| 两个人的视频大全免费| 亚洲av在线观看美女高潮| 亚洲av电影在线观看一区二区三区| 亚洲欧美一区二区三区黑人 | 有码 亚洲区| 国产av精品麻豆| 超碰97精品在线观看| 夫妻午夜视频| 插逼视频在线观看| 中文字幕亚洲精品专区| 国产精品一区二区在线不卡| 高清午夜精品一区二区三区| 亚洲国产欧美日韩在线播放 |