999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于非線性超聲調(diào)制方法的損傷識別與定位

2016-11-23 09:08:43屈文忠
振動、測試與診斷 2016年5期
關(guān)鍵詞:裂紋信號方法

屈文忠, 李 拯, 王 芝, 肖 黎

(武漢大學(xué)工程力學(xué)系 武漢,430072)

?

基于非線性超聲調(diào)制方法的損傷識別與定位

屈文忠, 李 拯, 王 芝, 肖 黎

(武漢大學(xué)工程力學(xué)系 武漢,430072)

在結(jié)構(gòu)健康監(jiān)測技術(shù)中,非線性彈性波譜方法具有對結(jié)構(gòu)微小變化敏感的特性,能夠有效地對裂紋等非線性損傷進(jìn)行識別。筆者針對采用兩個持續(xù)激勵的普通非線性彈性波譜方法不能定位損傷的問題,提出了一種能夠識別并且定位鋁板中疲勞裂紋的非線性超聲調(diào)制方法。該方法通過識別脈沖與高頻超聲波之間的調(diào)制現(xiàn)象來進(jìn)行損傷檢測。實(shí)驗(yàn)中,壓電陣列粘貼于疲勞裂紋鋁板表面,漢寧窗調(diào)制的正弦脈沖激勵和正弦持續(xù)激勵同時施加在壓電陣列上。通過采集不同的作動傳感路徑的響應(yīng),利用短時傅里葉變換對響應(yīng)進(jìn)行頻域分析,構(gòu)造損傷指數(shù)量化損傷程度,對疲勞裂紋進(jìn)行識別和定位。實(shí)驗(yàn)結(jié)果表明,所提出的方法可以成功地檢測并定位疲勞裂紋損傷。

結(jié)構(gòu)健康監(jiān)測; 疲勞裂紋; 非線性超聲調(diào)制; 短時傅里葉變換; 壓電陣列

引 言

金屬材料具有高強(qiáng)度和低成本的優(yōu)點(diǎn),在土木、機(jī)械和航空航天結(jié)構(gòu)等領(lǐng)域得到廣泛應(yīng)用。在長期使用過程中,由于撞擊、化學(xué)腐蝕以及疲勞等因素影響,結(jié)構(gòu)中有可能產(chǎn)生裂紋等損傷。缺陷結(jié)構(gòu)檢測維修不及時,最終會導(dǎo)致結(jié)構(gòu)失效。對于金屬結(jié)構(gòu)微小損傷尤其是疲勞裂紋的檢測技術(shù)的研究日益受到人們的重視[1]。傳統(tǒng)的線性檢測技術(shù)基于損傷附近超聲波的反射、散射、衰減以及模式轉(zhuǎn)換等變化進(jìn)行損傷識別,但對于疲勞裂紋等微小損傷不敏感[2]。與線性檢測方法相比,非線性檢測方法基于諧波、亞諧波和調(diào)制等非線性特征進(jìn)行損傷識別,對于疲勞裂紋等微小損傷更加敏感[3]。

非線性超聲調(diào)制方法是一種典型的基于非線性特征的檢測方法[4]。疲勞裂紋開合、復(fù)合材料層脫和螺栓松動等都會產(chǎn)生明顯的非線性特征,采用非線性調(diào)制方法可以識別這些特征[5-6]。Meo等[7]利用非線性彈性波譜方法探測復(fù)合材料中的裂紋,在各向異性材料中成功識別裂紋損傷。Kim等[8]在風(fēng)力渦輪機(jī)工作時,將風(fēng)機(jī)勻速轉(zhuǎn)動的振動作為低頻振動信號,利用振動聲調(diào)制方法探測風(fēng)機(jī)葉片的疲勞裂紋。Sohn等[9]利用非線性超聲調(diào)制方法檢測鋁板和飛機(jī)機(jī)翼與機(jī)身連接接頭試件中的疲勞裂紋。Kazakov等[10]提出非線性調(diào)幅方法,該方法需發(fā)射周期性的超聲波脈沖,通過脈沖的調(diào)幅來識別裂紋。胡海峰等[11]利用希爾伯特變換對鋁板的調(diào)制信號進(jìn)行解調(diào),研究了裂紋長度對振動聲調(diào)制信號調(diào)制強(qiáng)度的影響。Solodov等[12]研究非線性自調(diào)制技術(shù),并通過掃描定位板中的裂紋。

非線性調(diào)制方法通常采用兩個不同頻率的持續(xù)正弦信號作為激勵,一個是低頻波動或振動信號,另一個是高頻波動信號。將激勵施加在結(jié)構(gòu)的作動器上,利用傳感器采集響應(yīng)并進(jìn)行頻域分析。若頻域中存在調(diào)制成分則說明結(jié)構(gòu)中存在非線性損傷。然而,由于采用持續(xù)激勵,上述的非線性調(diào)制方法只能識別損傷而不能定位。筆者提出了一種能夠定位鋁板中疲勞裂紋損傷的非線性超聲調(diào)制方法,對存在疲勞裂紋的鋁板進(jìn)行了實(shí)驗(yàn)研究。實(shí)驗(yàn)驗(yàn)證了筆者所提出的非線性超聲調(diào)制方法能夠有效、準(zhǔn)確地識別并且定位結(jié)構(gòu)中的損傷,并給出直觀的損傷定位圖。

1 非線性超聲調(diào)制原理

在大幅值激勵狀態(tài)下,結(jié)構(gòu)中疲勞裂紋開合引起剛度變化是一種典型的非線性特征,非線性超聲調(diào)制方法通過測量非線性特征有效地探測損傷。由兩個正弦信號組成的激勵同時施加在結(jié)構(gòu)上,其中一個為低頻激勵,頻率為f1;另一個為高頻激勵,頻率為f2。采集響應(yīng)并進(jìn)行頻域分析。在頻域內(nèi),非線性系統(tǒng)的響應(yīng)與線性系統(tǒng)的響應(yīng)不同。線性系統(tǒng)內(nèi),應(yīng)力應(yīng)變關(guān)系滿足胡克定律

(1)

其中:σ為應(yīng)力;ε為應(yīng)變;E為彈性模量。

當(dāng)系統(tǒng)受到高、低頻信號同時激勵,應(yīng)變表示為

ε=A1sin2πf1t+A2sin2πf2t

(2)

其中:A1為低頻激勵的應(yīng)變幅值;A2為高頻激勵的應(yīng)變幅值。

將式(2)代入式(1),得

σ=A1Esin2πf1t+A2Esin2πf2t

(3)

線性系統(tǒng)的響應(yīng)為兩個激勵的線性疊加。若試件出現(xiàn)疲勞裂紋或其他非線性損傷,系統(tǒng)將變?yōu)榉蔷€性。非線性系統(tǒng)的應(yīng)力應(yīng)變關(guān)系可以表示[8]為

(4)

其中:α1和α2為非線性參數(shù)。

聯(lián)立式(4)和式(2),可得應(yīng)力表達(dá)式為

(5)

由式(5)可以看出,非線性系統(tǒng)響應(yīng)的頻率成分包含諧波成分(2f1和2f2)和非線性調(diào)制成分(f2±f1)。

非線性超聲調(diào)制方法采用的激勵通常是兩個持續(xù)的正弦信號,該方法可以準(zhǔn)確識別鋁板結(jié)構(gòu)中疲勞裂紋等損傷,但是不能定位損傷。因?yàn)楫?dāng)采用兩個持續(xù)激勵時,超聲波會在鋁板中朝各個方向不停傳播。只要鋁板中的任何位置出現(xiàn)疲勞裂紋,超聲波都會經(jīng)過損傷位置,傳感網(wǎng)絡(luò)中的任一條作動傳感路徑采集的響應(yīng)都含有非線性特征。即使在沒有損傷的路徑上也會有調(diào)制現(xiàn)象,因此無法找出疲勞裂紋的位置。針對上述問題,筆者提出一種改進(jìn)的非線性超聲調(diào)制方法,用來定位鋁板中的疲勞裂紋。該方法的關(guān)鍵是采用脈沖激勵,其原理如圖1所示。若采用短周期的脈沖信號,由于信號周期太少,導(dǎo)致頻域成分不明顯。綜合考慮傳感網(wǎng)絡(luò)分布和頻域分析效果,同時為了減少頻譜泄露,筆者采用漢寧窗調(diào)制的30周脈沖信號和持續(xù)正弦信號作為激勵[13]。在一條作動傳感路徑中,脈沖激勵和持續(xù)正弦激勵同時施加在作動器上,只采集傳感器特定時間段的響應(yīng),以保證采集的響應(yīng)中只包含超聲波激勵從最短路徑傳播過來的信號而沒有邊界或遠(yuǎn)離本條路徑的損傷反射的信號。若該條路徑存在疲勞裂紋損傷,則響應(yīng)的頻域中會有調(diào)制現(xiàn)象。組成作動傳感網(wǎng)絡(luò),通過該方法比較各條路徑的非線性特征,可以準(zhǔn)確地識別損傷的位置。

圖1 非線性調(diào)制方法原理圖Fig.1 Schematic of the nonlinear modulation method

2 疲勞實(shí)驗(yàn)

利用MB公司的MODAL50激振器進(jìn)行疲勞實(shí)驗(yàn),在鋁板上生成疲勞裂紋。如圖2所示,將鋁板四周固定在框架上,在板中間位置開個小槽,將激振器頂桿置于小槽附近,利用應(yīng)力集中生成初始疲勞裂紋。

圖2 鋁板疲勞實(shí)驗(yàn)Fig.2 Fatigue testing arrangements

激振器由頻率為5 Hz的正弦信號驅(qū)動,信號峰-峰值為2 V,功率放大器增益調(diào)到最大。在激勵大約535 000周后,小槽底端出現(xiàn)疲勞裂紋,然后將正弦信號峰-峰值降到1.5 V,保持增益不變。再經(jīng)歷43 200周激勵后,疲勞裂紋長度達(dá)到15 mm,停止疲勞實(shí)驗(yàn)。產(chǎn)生疲勞裂紋后,可進(jìn)行非線性調(diào)制實(shí)驗(yàn)定位裂紋損傷。

3 非線性超聲調(diào)制實(shí)驗(yàn)

3.1 實(shí)驗(yàn)裝置與過程

實(shí)驗(yàn)裝置如圖3所示,其中鋁板尺寸為900 mm×600 mm×1 mm。鋁板的幾何形狀和壓電片(PZT)的位置如圖4所示。采用STEMINC公司的圓形壓電片,直徑為12 mm,厚度為0.6 mm。在鋁板表面粘貼兩列壓電片,每列6個。作動傳感網(wǎng)絡(luò)由這兩列壓電片組成,其中一列壓電片作為作動器,另一列作為傳感器。壓電片的坐標(biāo)如表1所示。疲勞裂紋長度為15 mm,其中心坐標(biāo)為(283.5 mm,450 mm)。利用Agilent公司的33522A任意函數(shù)發(fā)生器產(chǎn)生激勵信號;由TEGAM公司的2350電壓放大器將信號放大;利用Agilent公司的DSO-X 3014A示波器采集信號。

圖3 實(shí)驗(yàn)裝置Fig.3 Experimental setup

圖4 壓電陣列和損傷位置示意圖(單位:mm)Fig.4 Geometry of plate showing the PZT array and damage position (unit: mm)

傳感器Lx,y/mm作動器Lx,y/mmACT1(125,250)SEN1(125,650)ACT2(195,250)SEN2(195,650)ACT3(265,250)SEN3(265,650)ACT4(335,250)SEN4(335,650)ACT5(405,250)SEN5(405,650)ACT6(475,250)SEN6(475,650)

漢寧窗調(diào)制的30周50 kHz正弦脈沖和170 kHz的正弦波作為激勵信號,電壓峰-峰值均為50 V。利用BNC三通連接頭,將兩個激勵信號同時施加在1個壓電片(作動器)上。170 kHz的正弦波信號持續(xù)施加在作動器上,以脈沖信號開始施加在作動器上作為傳感器測量響應(yīng)的時間零點(diǎn)。同樣的激勵信號輪流激勵6個作動器(ACT1~ACT6),每當(dāng)激勵1個作動器時,6個傳感器(SEN1~SEN6)同時記錄響應(yīng),稱這些響應(yīng)為測量信號。由于每條作動傳感路徑長度不同,這必然會影響各條路徑非線性特征的比較。在頻域內(nèi),距離較短的完整路徑的響應(yīng)在調(diào)制頻率220 kHz處的值有可能比距離較長的損傷路徑的值大。因此,直接比較各條作動傳感路徑的調(diào)制成分會導(dǎo)致錯誤的判斷。為了避免路徑長度不同的影響,必須引入?yún)⒖夹盘枴,F(xiàn)將只有脈沖信號施加在作動器上時傳感器的響應(yīng)作為參考信號。由于各條路徑的參考信號在頻率220 kHz值的大小同樣受到路徑長度不同的影響,故將各條路徑測量信號在頻率220 kHz處的值減去參考信號在頻率220 kHz處的值便可去除路徑長度不同的影響。相減后的值可以準(zhǔn)確表征各條路徑的非線性特征。

3.2 損傷指數(shù)

構(gòu)造損傷指數(shù)量化各條路徑的非線性特征。以ACT4~SEN6(完整路徑)和ACT5~SEN1(損傷路徑)這兩條路徑為例說明。利用短時傅里葉變換(short-time Fourier transform,簡稱STFT)計算這兩條路徑測量信號和參考信號的頻譜[14],計算結(jié)果如圖5所示。根據(jù)式(5)中的理論分析,損傷路徑響應(yīng)頻譜圖中應(yīng)包含脈沖激勵倍頻信號頻譜100 kHz和非線性調(diào)制信號差頻頻譜120 kHz以及和頻頻譜220 kHz。由于實(shí)驗(yàn)中鋁板上的疲勞裂紋 非線性程度較弱以及實(shí)驗(yàn)過程中噪聲干擾,響應(yīng)信號頻譜中,脈沖激勵倍頻信號頻譜100 kHz和非線性調(diào)制信號差頻頻譜120 kHz較弱,無法與背景噪聲區(qū)分,非線性調(diào)制信號和頻頻譜220 kHz較強(qiáng)。因此,選擇220 kHz的調(diào)制成分構(gòu)造損傷指數(shù)。需要說明的是,筆者主要研究脈沖激勵50 kHz與持續(xù)正弦激勵170 kHz之間的非線性調(diào)制現(xiàn)象,選取的分析頻率范圍為0~270 kHz,故頻譜圖中沒有顯示出持續(xù)正弦激勵倍頻信號頻譜340 kHz。

圖5 響應(yīng)的頻譜Fig.5 The spectrograms of responses

(6)

由于每次激勵一個作動器,6個作動器輪流激勵,因此將各條路徑的損傷指數(shù)分為6組,每組包含一個作動器到6個傳感器的6條路徑的損傷指數(shù),各組之間相互沒有影響。為保證各組損傷指數(shù)在定位損傷時的權(quán)重相同,將每組的6個損傷指數(shù)進(jìn)行歸一化,最終得到每條路徑的損傷指數(shù)。所有路徑的損傷指數(shù)如表2所示,損傷指數(shù)值越大,表明該路徑上存在損傷的可能性越大,或者損傷距離路徑越近。

表2 各條路徑的損傷指數(shù)

3.3 損傷定位方法

傳感網(wǎng)絡(luò)中各檢測點(diǎn)存在損傷的概率由各條作動傳感路徑的損傷指數(shù)決定,檢測點(diǎn)與某條路徑的距離決定了該條路徑對該點(diǎn)計算損傷存在概率時所起作用的大小。路徑與檢測點(diǎn)的距離越大,所起的作用越小,于是令作動傳感路徑對損傷存在概率計算的影響符合正態(tài)分布[15]

(7)

(8)

利用36條作動傳感路徑定位損傷結(jié)果如圖6所示,圖中紅色圓圈表示作動器和傳感器的位置,短線表示疲勞裂紋的位置。像素點(diǎn)紅色越深表明疲勞裂紋位于該點(diǎn)的概率越大。從圖中可以看出,疲勞裂紋位于紅色最深的區(qū)域。實(shí)驗(yàn)結(jié)果表明,非線性超聲調(diào)制方法能準(zhǔn)確定位損傷。

圖6 利用36條路徑定位損傷Fig.6 Damage imaging constructed with 36 actuator-sensor paths

4 結(jié)束語

筆者提出了一種基于非線性超聲調(diào)制技術(shù)的損傷檢測方法,闡釋了其基本原理并進(jìn)行相關(guān)實(shí)驗(yàn)。該方法采用脈沖信號作為激勵,克服了傳統(tǒng)非線性超聲調(diào)制方法不能定位損傷的缺點(diǎn)。實(shí)驗(yàn)中,通過在鋁板上粘貼壓電片構(gòu)成傳感網(wǎng)絡(luò),利用短時傅里葉變換對各條路徑的響應(yīng)信號進(jìn)行頻域分析,通過脈沖和持續(xù)正弦激勵的非線性調(diào)制現(xiàn)象構(gòu)造損傷指數(shù),量化各條路徑的非線性特征,通過定位算法進(jìn)行損傷定位。實(shí)驗(yàn)結(jié)果表明,筆者提出的非線性超聲調(diào)制方法可以準(zhǔn)確定位鋁板中的疲勞裂紋損傷。該方法不僅可以識別和定位疲勞裂紋,對于其他非線性損傷同樣有效。

[1] Sohn H, Farrar C R, Hemez F M. A review of structural health monitoring literature: 1996-2001[R]. [S.l.]: Los Alamos National Laboratory, LA-13976-MS, 2004.

[2] Kim J Y, Jacobs L J, Qu J, et al. Experimental characterization of fatigue damage in a nickel-base super alloy using nonlinear ultrasonic waves [J]. The Journal of the Acoustical Society of America, 2006, 120(3): 1266-1273.

[3] Zaitsev V, Nazarov V, Gusev V, et al. Novel nonlinear-modulation acoustic technique for crack detection[J]. NDT & E International, 2006, 39(3): 184-194.

[4] Van Den Abeele K E A, Johnson P A, Sutin A. Nonlinear elastic wave spectroscopy (NEWS) techniques to discern material damage, part I: nonlinear wave modulation spectroscopy (NWMS) [J]. Research in Nondestructive Evaluation, 2000, 12(1): 17-30.

[5] Duffour P, Morbidini M, Cawley P. A study of the vibro-acoustic modulation technique for the detection of cracks in metals[J]. The Journal of the Acoustical Society of America, 2006, 119(3): 1463-1475.

[6] Klepka A, Staszewski W J, Jenal R B, et al. Nonlinear acoustics for fatigue crack detection-experimental investigations of vibro-acoustic wave modulations [J]. Structural Health Monitoring, 2012, 11(2): 197-211.

[7] Meo M, Polimeno U, Zumpano G. Detecting damage in composite material using nonlinear elastic wave spectroscopy methods[J]. Applied Composite Materials, 2008, 15(3): 115-126.

[8] Kim S, Adams D E, Sohn H. Crack detection on wind turbine blades in an operating environment using vibro-acoustic modulation technique[C]∥ Aip Conference Proceedings. Denver, America: American Institute of Physics, 2013, 1511(1): 286-293.

[9] Sohn H, Lim H J, DeSimio M P, et al. Nonlinear ultrasonic wave modulation for online fatigue crack detection [J]. Journal of Sound and Vibration, 2014, 333(5): 1473-1484.

[10]Kazakov V V, Sutin A, Johnson P A. Sensitive imaging of an elastic nonlinear wave-scattering source in a solid [J]. Applied Physics Letters, 2002, 81(4): 646-648.

[11]胡海峰, 胡蔦慶, 秦國軍. 非線性振動聲調(diào)制信號耦合特征分析[J]. 機(jī)械工程學(xué)報, 2010,46(23): 68-76.

Hu Haifeng, Hu Niaoqing, Qin Guojun. Coupling Analysis of Nonlinear Vibro-acoustic Modulation Signals[J]. Journal of Mechanical Engineering, 2010,46(23): 68-76. (in Chinese)

[12]Solodov I, Wackerl J, Pfleiderer K, et al. Nonlinear self-modulation and subharmonic acoustic spectroscopy for damage detection and location [J]. Applied Physics Letters, 2004, 84(26): 5386-5388.

[13]張宇, 余龍, 閆云聚. Lamb 波結(jié)合相控陣技術(shù)的結(jié)構(gòu)損傷檢測[J]. 振動、測試與診斷, 2013, 33(6): 977-980.

Zhang Yu, Yu Long, Yan Yunju. Application of phased array technique in lamb wave based structural damage detection[J] Journal of Vibration, Measurement & Diagnosis, 2013, 33(6): 977-980. (in Chinese)

[14]Michaels J E, Michaels T E. Detection of structural damage from the local temporal coherence of diffuse ultrasonic signals[J]. Ultrasonics, Ferroelectrics and Frequency Control, IEEE Transactions, 2005, 52(10): 1769-1782.

[15]Dong Wang, Lin Ye, Ye Lu. A probabilistic diagnostic algorithm for identification of multiple notches using digital damage fingerprints (DDFs)[J]. Journal of Intelligent Material Systems and Structures, 2009, 20(12):1439-1450.

10.16450/j.cnki.issn.1004-6801.2016.05.006

*國家自然科學(xué)基金資助項目(51078293,51378402)

2014-08-09;

2014-12-15

TH113

屈文忠,男,1968年8月生,博士、教授、博士生導(dǎo)師。主要研究方向?yàn)檎駝庸こ毯徒Y(jié)構(gòu)健康監(jiān)測等。曾發(fā)表《Lamb wave damage detection using time reversal DORT method》(《Smart Materials and Structures》2013, Vol.22, No.4)等論文。

E-mail: qwz@whu.edu.cn

猜你喜歡
裂紋信號方法
裂紋長度對焊接接頭裂紋擴(kuò)展驅(qū)動力的影響
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
完形填空二則
Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
基于FPGA的多功能信號發(fā)生器的設(shè)計
電子制作(2018年11期)2018-08-04 03:25:42
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
基于LabVIEW的力加載信號采集與PID控制
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
捕魚
預(yù)裂紋混凝土拉壓疲勞荷載下裂紋擴(kuò)展速率
主站蜘蛛池模板: 浮力影院国产第一页| 国产91精品最新在线播放| 青青操国产视频| 99精品免费欧美成人小视频| 欧美日韩一区二区三区四区在线观看| 综合亚洲网| 国产精品对白刺激| 国产夜色视频| 中文字幕第1页在线播| 99热国产在线精品99| 中文字幕在线免费看| 自拍中文字幕| 成人毛片免费观看| 久久无码av一区二区三区| 国产精品青青| 久久视精品| 成人一区在线| 欧美亚洲日韩中文| 成年人视频一区二区| 制服丝袜无码每日更新| 99国产精品国产| 国产无码网站在线观看| 在线一级毛片| 国产成人免费| 亚洲欧美日韩高清综合678| 国产极品美女在线观看| 久久99精品久久久久纯品| 国产亚洲精品资源在线26u| 国产性爱网站| 国产一在线| 在线观看国产黄色| 高清乱码精品福利在线视频| 久久国产精品影院| 全部毛片免费看| 国产黄色免费看| 亚洲欧美极品| 免费高清a毛片| 国产精品亚洲一区二区三区z| 日韩一区二区在线电影| 国产主播在线一区| 这里只有精品在线| 国产精品一区二区在线播放| 国产丝袜91| 国产中文在线亚洲精品官网| 亚洲国产精品无码久久一线| 99re精彩视频| 欧美日韩国产在线播放| 女人一级毛片| 九九热精品在线视频| 一级毛片免费的| 被公侵犯人妻少妇一区二区三区| 日韩不卡高清视频| 国产无码在线调教| 久久这里只有精品23| 日韩人妻少妇一区二区| 国产一级小视频| 激情综合婷婷丁香五月尤物| 天堂成人av| 欧美激情二区三区| 久久频这里精品99香蕉久网址| 国产毛片高清一级国语 | 亚卅精品无码久久毛片乌克兰 | 99在线观看视频免费| 久久免费看片| 国产亚洲视频免费播放| 亚洲日本韩在线观看| 3p叠罗汉国产精品久久| 国产丝袜91| 又大又硬又爽免费视频| 欧美另类精品一区二区三区| 中字无码av在线电影| 国产午夜人做人免费视频| 国产成人成人一区二区| 99精品这里只有精品高清视频| 亚洲国产高清精品线久久| 国产流白浆视频| 亚洲乱码视频| 一区二区三区在线不卡免费| 亚洲国产av无码综合原创国产| 国产成人乱无码视频| 精品亚洲国产成人AV| 国产地址二永久伊甸园|