湯 真
(重慶交通大學(xué),重慶 400074)
20世紀(jì)30年代,人們開始研究混凝土無損檢測方法,80年代,美國學(xué)者Mary Sansalone利用機械波反射法進行無損檢測,近年來,測試內(nèi)容也由強度推定、內(nèi)部缺陷探測等擴展到混凝土結(jié)構(gòu)應(yīng)力范疇,2012年,宋麗莉等[1]對水泥和巖石樣品開展散射波試驗研究,證明多次散射波形成的尾波比直達波對介質(zhì)性質(zhì)變化更為敏感,2018年,黎志謀[2]開展了超聲尾波對混凝土梁三點彎曲應(yīng)力場超聲測試,利用泰勒展開式,建立了幅值一階變化系數(shù)與應(yīng)力變化近似線性曲線。
本文制定了混凝土工字梁偏心受壓超聲尾波應(yīng)力測試方案,以奇異值分解為基礎(chǔ)理論,建立了超聲尾波波形的幅值特征向量與混凝土結(jié)構(gòu)荷載等級關(guān)系曲線。
本方案對混凝土工字梁空載及20個荷載的偏心壓力作用下進行超聲測試,設(shè)計好荷載取值范圍及施加荷載的具體位置,利用Midas Civil建立試驗加載模型提取出相應(yīng)的應(yīng)力計算值,由RSM-SY5(T)型非金屬超聲波檢測儀發(fā)射超聲脈沖信號,通過接收多次散射波形成的尾波記錄及保存數(shù)據(jù)。
本次試驗采用的梁體構(gòu)件為鋼筋混凝土工字型梁,跨徑2 m,整體梁高45 cm,其中翼緣板寬30 cm,翼緣板高10 cm,腹板寬10 cm,腹板高25 cm??v向鋼筋和箍筋均采用直徑為10 mm的HPB300級鋼筋。
1.2.1 超聲測試裝置
超聲測試裝置由RSM-SY5(T)型非金屬超聲波檢測儀和換能器組成,平面換能器中心頻率50 kHz,直徑為4 cm,布置方式為一發(fā)雙收,換能器布設(shè)位置如圖1所示。
1.2.2 試驗加載系統(tǒng)
利用Midas Civil對不同荷載等級下工字型梁體的應(yīng)力大小進行計算,提取出相應(yīng)的應(yīng)力計算值,表1所示。
圖1 換能器布設(shè)位置示意圖(單位:cm)
表1 鋼筋混凝土工字型梁部分荷載等級下的應(yīng)力值
加載方式采用單級加載,采用千斤頂對混凝土工字梁端面施加偏心壓力,千斤頂力作用點距梁體側(cè)表面形心軸下方7.5 cm。為保持千斤頂對梁體所施加的壓力不回縮并實時監(jiān)控施加偏心壓力,在梁體兩側(cè)鋼筋貼應(yīng)變片并配套應(yīng)變檢測儀,加載系統(tǒng)如圖2所示。
圖2 加載系統(tǒng)示意圖
1.2.3 試驗流程
首先對超聲測試儀器進行30 min左右的預(yù)熱,然后利用千斤頂對梁體進行偏心加載。每天進行一次完整混凝土工字梁偏心加載試驗,每次完整試驗共7組超聲測試,每組測試共21個荷載狀態(tài)(含自重),荷載范圍(1~20 kN),每一級荷載測試將記錄超聲波形11條,共進行5次完整試驗,如圖3所示。
圖3 現(xiàn)場加載測試示意圖
由于本試驗有2處接收端,故偏心加載試驗數(shù)據(jù)可構(gòu)成7×21×11×1 024×2的五維矩陣,為建立了超聲尾波波形幅值特征向量與荷載等級間關(guān)系曲線,在數(shù)據(jù)分析中以偏心荷載等級為自變量,把數(shù)據(jù)進行平均處理,原始波形進行奇異值分解,利用SVD矩陣性質(zhì)建立幅值與荷載等級關(guān)系,實現(xiàn)采用幅值特征向量來表征荷載等級的數(shù)據(jù)處理目標(biāo)。
選用任一端一次完整試驗測試數(shù)據(jù)得到7×21×11×1 024的4維矩陣W′,將一次采集數(shù)據(jù)進行預(yù)處理,物理意義容易理解每一級荷載有一個波形,將W矩陣按11次重復(fù)移點測試平均為一次,每一天重復(fù)測試8組平均為1組,得到21×1 024的二維矩陣。
由于每次超聲測試采集到的波形能量不盡相同,波形數(shù)據(jù)不具備可比性,為了使波形數(shù)據(jù)各類指標(biāo)位于同一數(shù)量級上,對其進行了能量歸一化處理,具體過程如下:一條完整的波形Wj由1 024個點構(gòu)成,Xj為幅值。
Wj=Wj(x1,x2……x1024)
(1)
波性能量:
(2)
波形能量歸一化處理:
(3)
數(shù)據(jù)經(jīng)預(yù)處理后得到21×1 024二維矩陣W,矩陣W為實數(shù)矩陣,W∈Rm×n,則一定存在正交矩陣U∈Rm×m,正交矩陣V∈Rn×n和對角矩陣Σ∈Rm×n,本文中m=1024,n=21,進行奇異值分解矩陣W可表示為[3]:
W=UΣVT
(4)
Σ可表示如下:
(5)
式中,0為零矩陣;p為min(m,n);σp為矩陣W的奇異值,且σ1≥σ2≥…≥σp≥0。
正交矩陣U與VT滿足:
UHU=Em,m
(6)
(VT)HVT=En,n
(7)
式中,UH、(VT)H分別為U、VT的共軛轉(zhuǎn)置為Em,m、En,n分別為m階、n階單位矩陣。
為建立不同荷載等級下超聲尾波波形數(shù)據(jù)特征值與荷載等級間的函數(shù)關(guān)系,設(shè)存在特征向量wx,使其與奇異值分解后矩陣相乘等于W與線性向量[ax+b]相乘可表示為:
wxUΣVT=[ax+b]W
(8)
式(8)左右兩側(cè)同時右乘:
wxEm,mEm,nEn,n=[ax+b]W(VT)HΣ-1UH
(9)
由于矩陣W′經(jīng)過多次平均處理,為向工程實際過渡,需要確保數(shù)據(jù)在組內(nèi)、組間的穩(wěn)定性、重復(fù)性,下面考察各組內(nèi)及組間數(shù)據(jù)結(jié)果,逐一進行分析。
1)組內(nèi)數(shù)據(jù)結(jié)果分析。從組內(nèi)角度考察各天的測試數(shù)據(jù):①單天測試數(shù)據(jù)構(gòu)成的11×1 024×7×21 四維數(shù)據(jù)矩陣;②7 個 11×1 024×21 的三維數(shù)據(jù)矩陣;③進行能量歸一化處理使得波形數(shù)據(jù)各類指標(biāo)位于同一數(shù)量級上;④采用基于奇異值分解的數(shù)據(jù)處理方式求取單組組內(nèi) 11 條波形各自的幅值特征向量Wx;⑤建立幅值特征向量與荷載等級間的近似線性關(guān)系曲線。
基于奇異值分解原理求取不同組數(shù)下組內(nèi) 11 條波形的幅值特征向量,建立了幅值特征向量與荷載等級間線性關(guān)系,如圖4所示:
圖4 組內(nèi)數(shù)據(jù)特征向量與不同荷載等級間的關(guān)系曲線
接收端1、接收端2單組組內(nèi)11條波形的特征向量Wx與荷載等級間呈現(xiàn)近似線性關(guān)系,并且特征向量Wx的離散程度極小。
2)組間數(shù)據(jù)結(jié)果分析。單天測試數(shù)據(jù)構(gòu)成的11×1 024×7×21四維數(shù)據(jù)矩陣按測試組數(shù)的不同拆分為7個11×1 024×21的三維數(shù)據(jù)矩陣,每個三維數(shù)據(jù)矩陣代表單組測試過程中采集到的波形數(shù)據(jù),首先求取11條波形的平均,可將7個三維數(shù)據(jù)矩陣整合成7個21×1 024的二維矩陣。
進行能量歸一化處理使得波形數(shù)據(jù)各類指標(biāo)位于同一數(shù)量級上,以單天整體測試數(shù)據(jù)所構(gòu)成的二維矩陣W作為開展奇異值分解的基準(zhǔn),分別帶入單組測試數(shù)據(jù),求解出各組測試數(shù)據(jù)在不同荷載等級下的特征向量,建立不同組數(shù)據(jù)的幅值特征向量與荷載等級間的近似線性關(guān)系曲線,如圖5所示。
圖5 各組特征向量與不同荷載等級間的關(guān)系曲線
接收端各組測試數(shù)據(jù)任然存在特征向量Wx與荷載等級間近似線性關(guān)系曲線,但同一荷載等級下各組測試數(shù)據(jù)的特征向量Wx上下位置隨機分布,未出現(xiàn)按加載測試先后順序分布。剔除第一組(黑色線),其余六組測試數(shù)據(jù)的特征向量Wx間離散程度較小,能達到利用特征向量Wx來識別荷載等級的數(shù)據(jù)處理目標(biāo)。
第一組測試數(shù)據(jù)特征向量Wx明顯偏離于其余六組,主要考慮兩方面,①儀器設(shè)備考慮預(yù)熱及梁體本身對荷載反應(yīng)過程記憶了之前的受力歷史擁有了“凱撒效應(yīng)”;②耦合劑與梁體耦合狀態(tài)穩(wěn)定時間。
對偏心受壓應(yīng)力場下構(gòu)件及混凝土梁開展了超聲尾波試驗,基于奇異值分解后的超聲尾波波形幅值特征向量來表征偏心荷載等級的結(jié)構(gòu)測試模型,建立了超聲尾波波形幅值特征向量與荷載等級間的近似線性關(guān)系曲線。
從組內(nèi)的角度對單天測試數(shù)據(jù)展開處理及分析,單組組內(nèi) 11條波形數(shù)據(jù)的幅值特征向量與荷載等級間呈現(xiàn)出近似線性關(guān)系,并且不同荷載等級下單組組內(nèi) 11 條波形數(shù)據(jù)的幅值特征向量離散程度極??;從組間的角度對單天測試數(shù)據(jù)展開處理及分析,線性關(guān)系明顯,雖然不同組的波形幅值特征向量間出現(xiàn)了隨機的偏移,但識別荷載等級仍然有效。
[ID:009085]