錨桿錨固系統(tǒng)應(yīng)力波分形機制及其應(yīng)用研究

馮劍鋒，趙金昌，梁騰飛，賈少華，尹志強，閆 東，卞德存

(太原理工大學(xué) 礦業(yè)工程學(xué)院，太原 030024)

?

錨桿錨固系統(tǒng)應(yīng)力波分形機制及其應(yīng)用研究

馮劍鋒，趙金昌，梁騰飛，賈少華，尹志強，閆 東，卞德存

(太原理工大學(xué) 礦業(yè)工程學(xué)院，太原 030024)

為研究錨桿振動系統(tǒng)應(yīng)力波分形特性，采用分形理論的盒維數(shù)方法，建立了適用于錨桿振動信號的盒維數(shù)模型。結(jié)合錨桿工作荷載實驗，計算得到了錨桿振動信號分形盒維數(shù)值D，從理論上探討了錨桿振動信號波形的物理力學(xué)機制，解釋了D與錨桿工作荷載之間的內(nèi)在聯(lián)系。研究表明，D與錨桿振動基頻及錨固介質(zhì)阻尼關(guān)系密切，振動基頻反映了振動信號的不規(guī)則度和復(fù)雜度，而錨固介質(zhì)阻尼直接影響振動信號的衰減，從而影響D的數(shù)值大小。將錨桿無損檢測與分形理論結(jié)合，為錨桿無損檢測分析領(lǐng)域提供了全新的研究思路。

錨桿錨固；應(yīng)力波；盒維數(shù)；工作荷載；分形機制；基頻；無損檢測

目前,在地下工程施工中,采用錨桿加固圍巖的方法已經(jīng)被普遍應(yīng)用。然而由于多方面原因,在工程現(xiàn)場對錨桿錨固質(zhì)量進行檢測時,大多仍采取破壞性試驗的手段,這種檢測方法不僅效率低下,而且會對周邊圍巖形成比較強烈的擾動,破壞圍巖結(jié)構(gòu),大大弱化了錨桿的支護效果,這種弊端在巖層比較軟或者遇到破碎巖層時尤為突出。因此，亟待探求一種方法,在不破壞圍巖的前提下,對錨桿錨固質(zhì)量進行快速、準確的判斷。

隨著對錨桿錨固系統(tǒng)振動特性的深入研究,發(fā)現(xiàn)錨桿錨固系統(tǒng)應(yīng)力波作為各種頻率成分簡諧波的組合,其波形表現(xiàn)相當(dāng)復(fù)雜；但在相同的錨固參數(shù)條件下,進行重復(fù)、多次實驗所獲取的振動信號波形參數(shù),如頻率、反射時間等參數(shù)基本相近,且滿足統(tǒng)計自相似性隨機分型特性[1]。分形理論作為一種新理論,近年來已經(jīng)被引入到爆破振動信號、地震波信號的處理等研究領(lǐng)域。為了拓展分形理論在錨桿無損檢測領(lǐng)域的應(yīng)用,筆者提出將錨桿錨固系統(tǒng)振動信號的分形盒維數(shù)與錨桿工作荷載相結(jié)合,通過盒維數(shù)模型的建立、實驗研究與理論分析,給出判斷錨桿工作荷載的一個新思路。

1 盒維數(shù)模型定義

設(shè)錨桿錨固系統(tǒng)振動信號W(t)?F,F是n維歐式空間Rn上的閉集。Nδ是用來覆蓋F所需邊長為r的n維立方體最小數(shù)目,則W(t)的盒維數(shù)定義為[2-3]:

本次實驗的數(shù)據(jù)采集儀在采集信號時是每間隔一段時間Δt采集一次信號,最后將采集到的所有數(shù)據(jù)點連成光滑曲線。因此，嚴格意義上來說,在一個采樣時間間隔內(nèi),兩個數(shù)據(jù)點之間是一段直線,因此在用長度為δ的盒子對集合F進行覆蓋時,δ>Δt。在實際計算過程中,一般不根據(jù)盒維數(shù)的定義直接得到D,而是根據(jù)分形對象在其無標(biāo)度區(qū),將一系列尺寸為kδ(k為網(wǎng)格放大倍數(shù))的方形網(wǎng)格對分形對象進行覆蓋,從而計算出不同尺度下的有效覆蓋網(wǎng)格數(shù)量Nkδ,采用最小二乘法得出lg(1/kδ)-lgNkδ的擬合直線,其斜率就是該波形的盒維數(shù)D[3]。

2 錨固系統(tǒng)振動信號盒維數(shù)計算模型

如圖1所示的波形就是典型的錨桿振動信號波形。錨桿錨固系統(tǒng)振動信號包括兩個尺度,縱向表示其振幅,橫向表示時差。要對振動信號曲線進行較為準確的覆蓋,必須建立與其相匹配的矩形盒模型,從而表征其特性[4]。本文采用的矩形盒模型尺度為δ1×δ2, 即橫向×縱向。

圖1 典型的錨桿振動信號Fig.1 Typical bolt vibration signal

如圖2所示,試件全長L=1.6 m,錨固段長度L2=0.5 m,則固端反射時間

底端反射時間

取應(yīng)力波在桿體中的傳播速度v1=5 070 m/s,在錨固系統(tǒng)中的傳播速度近似為v1=4 000 m/s,則振動信號包含固端反射和底端反射的時間可以確定。實驗中的采樣頻率為500 kHz,所以時間間隔為

圖2 試件示意圖Fig.2 Schematic diagram of pieces

關(guān)于縱向尺度,由于其與振幅相關(guān),而振幅又受到激發(fā)裝置的激發(fā)沖量大小影響,所以沒有一個固定的值。在實驗過程中,我們盡可能地讓每次激發(fā)的沖量相等。通過觀察對比,得到最大振幅(振幅用電荷Q表示)差值為ΔQ<11 nC,故取ΔQ=11 nC.同時,為保證選取的數(shù)據(jù)點數(shù)能夠完整包含一次固端反射和底端反射及顯示首波,我們選取n=500個數(shù)據(jù)點,得到

在用矩形盒覆蓋振動信號時,-lg(kδi)與lgNkδi滿足線性回歸方程:

盒維數(shù)計算式為[4-5]:

圖3為本次實驗中的一條擬合直線。

3 實驗研究

實驗選用3組試件。為研究錨固介質(zhì)的阻尼對盒維數(shù)的影響,用介質(zhì)的不同密實程度來模擬介質(zhì)的不同阻尼,其他錨固參數(shù)相同。外徑90 mm的鋼管中的端錨錨桿,硅酸鹽水泥型號425,用水泥與細沙質(zhì)量比為1∶2的水泥砂漿作為錨固劑,水灰質(zhì)量比1∶0.5,養(yǎng)護28 d后進行測試,試件參數(shù)見表1所示。

表1 試件參數(shù)

在MATLAB中,根據(jù)上述計算模型,對得到的加速度響應(yīng)信號進行處理,如圖3所示。

圖3 盒維數(shù)擬合結(jié)果Fig.3 Box dimension fitting results

利用DAP軟件提取信號基頻[5]。DAP7.1瞬態(tài)測試分析軟件是由成都泰斯特電子信息有限責(zé)任公司開發(fā)的一款用于瞬態(tài)信號測試、采集、分析和存儲的軟件,在數(shù)據(jù)處理方面有數(shù)字濾波器、時域處理、頻域處理、特征值顯示等功能。

針對此次實驗數(shù)據(jù),選擇該軟件的頻域處理功能,對信號進行頻域處理,處理類型選擇幅值線性譜,其計算式為:

fft為快速傅立葉變換,其公式為:

加載函數(shù)選擇Hanning,其公式為:

頻域處理函數(shù)如圖4所示。盒維數(shù)計算結(jié)果及基頻見表2和圖5所示。根據(jù)實驗結(jié)果可得出如下結(jié)論。

圖4 頻域處理函數(shù)Fig.4 Frequency domain processing window

荷載/kN試件1試件2試件3盒維數(shù)D基頻f/kHz盒維數(shù)D基頻f/kHz盒維數(shù)D基頻f/kHz01.13063.3201.09953.0271.10222.9351.13183.5161.10583.2231.11534.004101.14313.7111.11143.8091.11654.199151.14944.1021.11294.0041.12114.297201.15994.2971.11514.6881.12724.590251.16395.0781.11725.1761.13724.883301.16585.2731.11935.3711.14174.980351.16885.4691.12665.5661.14815.078401.17155.6641.13285.7621.15355.176451.15475.8591.14155.9571.16635.237501.16585.0781.11694.9801.13404.883551.13854.6881.11494.8831.12414.492601.13924.4921.11244.5901.11424.199

圖5 盒維數(shù)、基頻與荷載對應(yīng)關(guān)系Fig.5 Relationship between box-dimension,the fundamental frequency and load

1) 從表2測試數(shù)據(jù)可以發(fā)現(xiàn),3組試件基頻值及盒維數(shù)值均以45 kN錨固失效點為拐點，如圖5所示。荷載在0～45 kN階段時,隨著荷載的增加,錨桿振動基頻值及振動信號盒維數(shù)值增大；荷載增加到45 kN以后,錨固體松動,水泥砂漿開始脫落,相應(yīng)地錨桿振動基頻值及振動信號盒維數(shù)值減小。

2) 由于3組試件錨固密實程度有差異,而錨桿的固有頻率又受周圍錨固介質(zhì)抗剪強度的影響較大,故而其基頻值也有相應(yīng)差異。

3) 錨固介質(zhì)的阻尼直接影響到錨桿振動信號的衰減,介質(zhì)阻尼越大,信號越趨于平緩,表現(xiàn)在盒維數(shù)值上,即介質(zhì)阻尼越大,盒維數(shù)值越小。故而3組試件的盒維數(shù)值大小關(guān)系為:

試件1>試件3>試件2 .

4 機理研究

由于錨桿錨固系統(tǒng)受力及振動情況相當(dāng)復(fù)雜,為了便于計算,做出下列假設(shè)。

1)錨桿側(cè)介質(zhì)是均勻的。

2)用一個線性彈簧和一個與速度有關(guān)的阻尼器以平行方式耦合,以模仿錨桿周圍介質(zhì)與錨桿之間的相互作用。

3)圍巖深度對錨桿周圍介質(zhì)的剪切應(yīng)力不構(gòu)成影響。

根據(jù)上述假設(shè),結(jié)合理論分析,得出錨桿的工作荷載F與錨桿錨固系統(tǒng)振動基頻f呈以下關(guān)系,即[6-10]:

式中:L為錨桿的長度;f為錨桿的振動基頻;vc為激發(fā)應(yīng)力波在錨桿桿體中的傳播速度,取5 070m/s;E為錨桿桿體彈性模量；A為錨桿的橫截面面積。錨桿為周圍介質(zhì)系統(tǒng)共同作用時得錨桿縱向波動方程[11-12]:

錨桿底端固定,頂端自由,則其振動邊界條件為:

解得,加速度響應(yīng):

在實際振動測試時,加速度傳感器固定在錨桿頂端,把x=L代入上式中可得:

綜上,對于錨桿振動信號,其曲線復(fù)雜程度可以根據(jù)它的頻率成分反映出來[10]。由上述的理論分析可知,應(yīng)力波在激發(fā)狀態(tài)下,隨著錨桿工作荷載的增加,錨桿頂端反射加劇,錨桿振動基頻隨之增大；表現(xiàn)在振動信號上即信號趨于“復(fù)雜”,對應(yīng)的信號變化加快；表現(xiàn)在波形圖像上,即信號波形所占據(jù)的平面比例進一步增大；而盒維數(shù)恰恰是反映圖形的復(fù)雜程度或不規(guī)則程度[13-15],相應(yīng)地盒維數(shù)值增大。隨著荷載的不斷增加,當(dāng)錨桿接近“拉脫”狀態(tài)時,此時錨桿相當(dāng)于錨固劣質(zhì)的錨桿。由于桿體與周圍介質(zhì)粘結(jié)力減小及周圍介質(zhì)松動,周圍介質(zhì)的抗剪強度降低,錨桿振動基頻減小[16-17],從而振動信號趨于“簡單”,即對應(yīng)的信號變化減慢,信號波形所占據(jù)的平面比例進一步減小,相應(yīng)的盒維數(shù)值減小。

對于錨固介質(zhì)阻尼不同的試件,由于介質(zhì)阻尼直接影響到信號的衰減,故介質(zhì)阻尼越大,信號衰減越快,表現(xiàn)在波形圖像上,即信號波形所占據(jù)的平面比例進一步縮小,相應(yīng)的盒維數(shù)值越小。

5 結(jié)論

1) 研究及實驗表明,根據(jù)錨桿振動信號的特點與規(guī)律,結(jié)合分形理論對錨桿振動波形進行雙尺度矩形盒維數(shù)建模可行、有效。

2) 對于同一根錨桿,隨著錨桿工作荷載的增大,錨桿振動信號盒維數(shù)值增大,當(dāng)錨固體出現(xiàn)松動,錨桿接近“拉脫”狀態(tài),即將要失效時,錨桿振動信號盒維數(shù)值隨著工作荷載的增大而減小。

3) 以上盒維數(shù)值的變化是因為隨著荷載的增加,錨桿振動基頻增大,振動信號趨于復(fù)雜,盒維數(shù)值相應(yīng)增大。直至錨固體出現(xiàn)松動,導(dǎo)致錨桿振動基頻減小,振動信號趨于簡單,盒維數(shù)值相應(yīng)減小。

4) 對于不同的錨桿,振動信號盒維數(shù)值大小受錨固介質(zhì)阻尼的影響較大,介質(zhì)阻尼越大,振動信號盒維數(shù)值越小。

[1] 張濟忠.分形[M].北京:清華大學(xué)出版社,1995.

[2] 朱華，姬翠翠.分形理論及其應(yīng)用[M].北京：科學(xué)出版社，2011.

[3] 婁建武，龍源，徐全軍，等.爆破地震信號的分形盒維數(shù)值分析[J].爆炸與沖擊，2004,24(4):364-369.

[4] 晏俊偉，龍源，婁建武，等.爆破地震波的分形機制[J].解放軍理工大學(xué)學(xué)報，2008,9(4)：367-372.

[5] 梁騰飛，馮劍鋒，卞德存,等.基于盒維數(shù)的錨桿工作荷載無損檢測研究[J].太原理工大學(xué)學(xué)報，2015,46(3)：327-331.

[6] 任智敏，李義.基于聲波測試的錨桿錨固質(zhì)量檢測信號分析與評價系統(tǒng)實現(xiàn)[J].煤炭學(xué)報，2011,36(1):191-196.

[7] 劉海峰，楊維武，李義.錨桿工作荷載實時檢測技術(shù)研究[J].礦業(yè)研究與開發(fā)，2008,28(2):54-57.

[8] 劉海峰.錨桿錨固系統(tǒng)及工作狀態(tài)動測技術(shù)研究[D].太原：太原理工大學(xué)，2000.

[9] 李義，張昌鎖，王成.錨桿錨固質(zhì)量無損檢測幾個關(guān)鍵問題的研究[J].巖石力學(xué)與工程學(xué)報，2008,27(1):108-116.

[10] 劉海峰，李義.應(yīng)力波反射法錨桿底端反射規(guī)律的實驗研究[J].太原理工大學(xué)學(xué)報，2000,25(2)：160-164.

[11] 李義，劉海峰，王富春.錨桿錨固狀態(tài)參數(shù)無損檢測及其應(yīng)用[J]，巖石力學(xué)與工程學(xué)報，2004,23(10):1741-1744.

[12] 李義.錨桿錨固質(zhì)量無損檢測與巷道圍巖穩(wěn)定性預(yù)測機理研究[D].太原：太原理工大學(xué)，2009.

[13] 王成.錨桿錨固質(zhì)量動測新技術(shù)的研究[D].太原：太原理工大學(xué)，1998.

[14] 謝全民，龍源，鐘明壽,等.小波與分形組合分析技術(shù)在爆破振動信號分析中的應(yīng)用[J]，振動與沖擊，2011,30(12):120-124.

[15] 劉國華，黃平捷，龔翔.基于分形維和獨立分量分析的聲發(fā)射特征提取[J].華南理工大學(xué)報(自然科學(xué)版)，2008,36(1):76-80.

[16] 謝全民，龍源，田作威.爆破振動信號時頻特征的三維分形特征研究[J].振動與沖擊，2010,29(12):122-126.

[17] 楊湖，王成.錨桿圍巖系統(tǒng)數(shù)學(xué)模型的建立及動態(tài)響應(yīng)分析[J].測試技術(shù)學(xué)報，2002,16(1):41-44.

[18] 李義，王成.應(yīng)力波反射法檢測錨桿錨固質(zhì)量的實驗研究[J].煤炭學(xué)報，2000,25(2):160-164.

(編輯：龐富祥)

Research on Fractal Mechanism of Bolt Anchor System Stress Wave and Its Application

FENG Jianfeng,ZHAO Jinchang,LIANG Tengfei,JIA Shaohua, YIN Zhiqiang,YAN Dong,BIAN Decun

(CollegeofMiningEngineering,TaiyuanUniversityofTechnology,Taiyuan030024,China)

For studying the stress wave fractal characteristics of bolt vibration system, using the fractal box-dimension method was used to, according to the propagation and characteristics of bolt vibration signal,establish the double scale rectangular box dimension model. Combined with the bolt working load experiments,calculations were carried out to obtain the bolt vibration signal fractal box dimension valuesD,the physical mechanism of anchor vibration signal waveform was theoretially discussed,and the relationship betweenDand the bolt working load was revealed.Research shows that,Dis close to the fundamental frequency and the anchoring medium damping, vibration fundamental frequency reflects the irregularity and complexity of the vibration signal, and the anchoring media damping directly affects the vibration signal attenuation,and thus affects theDnumerical size. Combining fractal theory with anchor nondestructive testing will provide a new research idea.

bolt anchor system;stress wave;box-dimension; working load;fractal mechanism;fundamental frequency;nondestructive testing

1007-9432(2016)03-0321-05

2015-10-30

國家自然科學(xué)基金資助項目：涵蓋峰后大變形過程的巷道圍巖與支護平衡規(guī)律及控制機理研究(51274145)；教育部博士點基金資助項目：群錨結(jié)構(gòu)錨桿荷載演化與圍巖穩(wěn)定性預(yù)測機理研究(20111402110003)

馮劍鋒(1989-)，男，山西澤州人，碩士生，主要從事結(jié)構(gòu)無損檢測方面的研究，(E-mail)15135151732@163.com

趙金昌(1974-)，男，副教授，主要從事結(jié)構(gòu)無損檢測及水中高壓脈沖放電方面的研究，(E-mail)zjc_8989@163.com

TD-05

A

10.16355/j.cnki.issn1007-9432tyut.2016.03.009