丁瑋 柴華友 劉思浩 陳燦 聶田 胡哲
摘 要: 敲擊 回波法是基樁完整性檢測的有效方法,基于一維應(yīng)力波理論,可以根據(jù)樁頂反射波質(zhì)點振動速度響應(yīng)曲線(包括走時、幅值及相位)來確定樁體波阻抗變化范圍及程度。然而,在實際中,樁是柱體而非桿件,截面變形不均勻。通過點源激發(fā)的波分量及在樁頂側(cè)面的反射波,分析了樁頂近場三維波動現(xiàn)象,并基于圓柱中縱向波對遠場波進行分析。結(jié)果表明,在測點位置及激振頻率滿足一定的要求時,三維波動特性及縱向波頻散特性可以削弱,樁的動測信號可近似用一維應(yīng)力波理論來分析。在樁土相互作用下,樁體中應(yīng)力波隨傳播衰減,分析了樁土剪切波速差異對應(yīng)力波衰減影響規(guī)律,結(jié)果表明,應(yīng)力波衰減不僅與樁土剪切波速比值有關(guān),而且還與樁土剪切波速絕對值有關(guān),樁體剪切波速越低、土體剪切波速越高,應(yīng)力波衰減越快。
關(guān)鍵詞: 強度準則;粗粒土;強度參數(shù);卸載;加載
中圖分類號:TU473.1? ? 文獻標志碼:A? ?文章編號:2096-6717(2020)02-0073-07
Analysis of influence of source and pile-soil interaction in low strain pile integrity testing
Ding Wei, Chai Huayou, Liu Sihao, Chen Can, Nie Tian, Hu Zhe
(School of Civil Engineering and Architecture, Wuhan Institute of Technology, Wuhan 430073, P.R.China)
Abstract:? The impact-echo method is an effective method for pile integrity testing.Based on the one dimensional (1-D) stress wave theory, the range and the extent of the changes in pile impedance can be determined from the vibration curve of reflected wave at the top of pile (including travel time, amplitude and phase).However, in practice, a pile is a cylinder rather than a rod for the induced waves, and the section deformation is not uniform.In this study, the three-dimensional wave phenomena exhibited in the wave-field near the pile top are analyzed from the components of the wave excited by the point sources and the multiple reflections from the side of the pile top.The propagation behavior of the waves in the far-field, which is far from the top, is analyzed based on the longitudinal waves in cylinders.The results show that the three-dimensional wave in the near field and the dispersive characteristics of the waves in the far field can be mitigated, and the signals measured at the top can be approximately analyzed using 1-D stress wave theory when the receiver position and source frequency content meet some conditions.Under the interaction between the pile and the soil, the stress waves in the pile decay with propagation.The effects of shear wave velocities of pile and soil on the attenuation of stress wave are analyzed.It is shown that the stress wave attenuation is not only related to the shear wave velocity ratio of the pile and soil, but also related to the absolute values of the shear wave velocity of the pile and soil.The lower the shear wave velocity of the pile is, the higher the shear wave velocity of the soil is, and the faster the stress waves decay.
Keywords: strength criterion; coarse-grained soils; strength parameter; unloading; loading
低應(yīng)變反射波法又稱錘擊 回波法,其基本原理為一維彈性桿應(yīng)力波理論。通過在樁頂施加一敲擊荷載,激發(fā)的彈性波沿樁身傳播,當遇到樁身波阻抗變化截面時(如擴頸、縮頸、斷樁等),應(yīng)力波會發(fā)生反射及透射現(xiàn)象。通過安裝在樁頂位置處的加速度傳感器可以接收到反射波質(zhì)點速度響應(yīng),基于一維波動理論可以得到應(yīng)力波沿樁身傳播的平均速度、缺陷類型及位置,從而對樁身完整性進行判定。
實際上樁體是柱體,在應(yīng)力波作用下,截面變形不均勻[1],同時,應(yīng)力波在傳播過程中受樁土相互作用影響,這樣,基于一維波動理論分析就有局限性,一些學(xué)者探討了樁身應(yīng)力波傳播三維效應(yīng)及樁周土對其傳播的影響。Liao等[2]通過對比一維波動理論和三維數(shù)值模擬結(jié)果,研究了樁頂?shù)娜S效應(yīng),結(jié)果表明三維效應(yīng)主要受頻率影響,高頻時表現(xiàn)更為強烈。陳凡等[3]分析了疊加在響應(yīng)中的振蕩,并將這種振蕩解釋為由于S波和R波在樁側(cè)面的多次反射所致。Chow等[4]用軸對稱有限元分析研究了源半徑與樁半徑比率對結(jié)果的影響,認為源和接收器之間的距離應(yīng)大于0.5R,可以最大限度地減少第1個峰值后的反向(負向)過沖分量。Zheng等[5]通過理論分析認為,在0.6R處側(cè)面反射波干擾影響最小。王雪峰等[6]基于應(yīng)力波理論,系統(tǒng)地研究了樁頂不同部位質(zhì)點速度響應(yīng)的振蕩程度以及應(yīng)力波視速度等隨錘樁比、波長樁徑比、長徑比的變化規(guī)律,給出了三維效應(yīng)的動測規(guī)避方式。陳安國等[7]建立了以樁的連續(xù)桿件模型和土的線彈性模型為基礎(chǔ)的低應(yīng)變數(shù)值計算模型,并利用行波理論求解出波動方程。岳向紅[8]基于理論、數(shù)值分析和工程實踐,對基樁中應(yīng)力波傳播三維效應(yīng)和衰減做出了探討。陳輝等[9]對不同缺陷基樁的低應(yīng)變反射波曲線進行數(shù)值模擬,得出反射波質(zhì)點速度響應(yīng)曲線,為實際基樁檢測提供了理論參考。Chai等[10-12]分析了測點布置、振源特征波長與樁徑比值對測試信號的影響。此外,樁土相互作用會影響波傳播衰減及樁土系統(tǒng)動力響應(yīng)[13],研究樁土相互作用對應(yīng)力波衰減的影響對基樁動測也非常重要。
由于樁是柱體,波在樁中傳播不滿足一維波動理論,因此,研究樁中波場對提高基樁結(jié)構(gòu)完整性分析的精度非常必要。筆者將樁體波場分為距樁頂較近的近場及距樁頂較遠的遠場,基于波場成份及傳播路徑分析樁側(cè)多次反射波影響較小的位置,基于柱體中縱向波傳播特性,分析頻率成份對遠場波傳播特性的影響,探討了樁頂測試信號的一維近似條件。在樁土相互作用下,樁體應(yīng)力波部分能量會向周圍土體輻射,導(dǎo)致沿樁體向下能量不斷衰減,筆者分析了樁土剪切波速等參數(shù)對樁土相互作用的影響。
1 樁中波場分析
1.1 樁頂近場
在均勻半空間表面上施加豎直向沖擊源時,向下的P波和S波將會以半球形波前的形式傳播,在介質(zhì)體內(nèi)遵循r-1的幾何衰減定律,其中,r是距離源中心的徑向距離。沿表面?zhèn)鞑サ闹边_P波和S波以r-2形式衰減;R波以圓柱形波陣面?zhèn)鞑ィ⑶易裱璻-1/2的幾何衰減定律。在P波、S波和R波中,R波的波速最慢,幾何衰減最慢,所以,R波將主導(dǎo)遠離源的表面場。S波和R波速度之間的回歸關(guān)系可以近似表示為[14-15]
cR≈ 0.87+1.12v 1+v cS (1)
式中:cS、cR分別為S波速度和R波速度;v為泊松比。
樁頂中心振源作用下近場及遠場波傳播如圖1所示,在到達樁的圓周側(cè)面之前,樁中波與均勻無限體中波有相同的傳播特性。在完整樁樁頂面中心區(qū)域用半正弦脈沖激振,數(shù)值模擬得到的不同時刻應(yīng)力波豎直向質(zhì)點速度幅值云圖如圖2所示,圖中中間黑色矩形框表示脈沖作用區(qū),白色虛線表示波前??梢钥闯?,激發(fā)的球面P、S波與柱面R波相互疊加,無法識別出各類型波。
當波遇到樁側(cè)面沿頂面?zhèn)鞑サ牟ㄔ趫A周側(cè)來回反射,直到它們的能量消散。由于P波的垂直運動在樁頂?shù)牟▓鲋邢鄬^小,因此,垂直表面運動由S波和R波主導(dǎo),多次反射波的基階共振頻率在式(2)所示范圍[3]。
cR 2R 式中:R為樁半徑,由于cR接近剪切波速cS,多次反射的共振頻率約等于cS/2R。 1.2 樁頂遠場 樁頂激振產(chǎn)生的下行波一部分在樁側(cè)邊界來回反射,反射波和向下傳播波的相互疊加會在距樁頂較遠處形成一種新的波,即導(dǎo)波。圖3給出了P波和S波在邊界反射,相互疊加在圓柱體形成導(dǎo)波的過程。取樁徑為0.2 m、樁長為5 m的完整單樁,進行低應(yīng)變數(shù)值模擬分析,圖4為樁頂以下不同位置豎直向質(zhì)點速度幅值云圖。由質(zhì)點速度沿截面分布可以看出,距樁頂較近波陣面是曲面(圖中白色虛線),樁軸線附近質(zhì)點速度相對兩側(cè)較大,隨著與樁頂距離增加,波陣面趨于平面,波陣面上質(zhì)點速度分布也較均勻,表明波在遠場近似一維桿波,在遠場,下行波傳播特性可以基于圓柱中縱向波傳播分析。 圓柱體中的導(dǎo)波有3種可能的傳播模式:縱向、扭轉(zhuǎn)和彎曲??v向模式的位移縱向軸對稱;扭轉(zhuǎn)模式的位移主要在圓周方向上;彎曲模式是非軸對稱的。當振源沿軸向施加于樁頂時,扭轉(zhuǎn)和彎曲模式可以被忽略。圓柱體中,縱向?qū)Рǔ尸F(xiàn)出多階模態(tài)和頻散特征??v向模式的頻散由Pochhammer頻率方程表示[1]。 2α R (β2+ξ2)J1(βR)-(β2-ξ2)J0(αR)×J1(βR)-? ?4ξ2αβJ1(αR)J0(βR)=0?(3) 式中:J0為第一類零階貝塞爾函數(shù);J1為第1類一階貝塞爾函數(shù)。 ξ、α和β定義為 ξ=2π/λL (4) α2=ω2/c2P-ξ2 (5) β2=ω2/c2S-ξ2 (6) 式中:ω和λL分別為縱向?qū)Рǖ膱A頻率和波長。模態(tài)截面位移分布與頻率有關(guān)。在給定頻率下,較高模態(tài)以較高的速度行進并具有較復(fù)雜的位移分布。 引入無量綱波數(shù)k=Rf/cP,其中,R是圓柱半徑;f是頻率;cP是縱波波速。對v=0.29,歸一化基階模態(tài)的頻散和截面變形形狀如圖5所示[16]。 圖中顏色較深的部分表示縱向位移??梢钥闯?,當k足夠小時(譬如k<0.1),導(dǎo)波大致表現(xiàn)為平面波。也就是說,不同頻率的波以幾乎相同的速度行進,并且每個橫截面變形近似均勻(參見圖5中的陰影區(qū)域)。 2 一維波動近似條件 在實際低應(yīng)變檢測中,敲擊脈沖可以近似用半正弦脈沖模擬,假設(shè)Td是半正弦脈沖的持續(xù)時間,脈沖的最大頻率fmax≈1/Td。樁頂敲擊激發(fā)的近場是三維的,波場含P、S、R波成份,可以用這些波的傳播特性及在樁頂側(cè)面反射來分析樁頂激發(fā)波響應(yīng)。在遠場,波傳播特性可以用基階模態(tài)縱向?qū)Рǚ治?,由于圖5頻散僅與無量綱波數(shù)有關(guān),為了建立敲擊脈沖寬度Td與無量綱波數(shù)的聯(lián)系,引入脈沖特征波長參數(shù),脈沖的特征波長定義為 Wp=Td×c0 (7) 式中:c0= E/ρ 是桿波速,E、ρ分別為楊氏模量及密度。半正弦脈沖的無量綱波數(shù)kmax的最大值可近似重寫為 kmax=Rfmax/cRc0/(WpcP) (8) 通過Ls-Dyna數(shù)值模擬方法分析敲擊脈沖寬度及樁側(cè)面反射波對樁頂不同位置質(zhì)點速度響應(yīng)及遠場波傳播的影響,并驗證以上理論分析??s頸樁模型如圖6所示,樁直徑(D=2R,R為半徑)為0.5 m,樁長L為10 m,樁材料的楊氏模量E、密度ρ、泊松比v、分別取為40 GPa、2 500 kg/m3、0.25,側(cè)面及樁底取為自由邊界。a是樁頂圓形加載區(qū)域的半徑,取a/R=1/8。在3個不同位置r/R=0.4、0.6、0.8處,無量綱脈沖特征波長Wp/R取4和8時,質(zhì)點響應(yīng)分別如圖7(a)、(b)所示,r為測點與樁中心軸的徑向距離。 圖7中,第1個波峰是對應(yīng)于沖擊脈沖的響應(yīng)(即直達S波和R波),波峰后的振蕩是來自樁圓周側(cè)的多次反射,多次反射波共振頻率滿足式(2)。樁底反射后強烈的振蕩是由基階縱向波不同頻率的波以不同的速度傳播引起,即頻散現(xiàn)象,圓柱中基階縱向波頻散曲線如圖5所示。通過對比圖7(a)、(b)可以發(fā)現(xiàn),增大Wp/R的比例,可以抑制多次反射強度,特別當傳感器放到離中心約0.6R的位置時,多次反射波在此位置范圍出現(xiàn)相消干涉現(xiàn)象,多次反射波影響較小。 對于Wp/R=4的情況,kmax≈0.25,由脈沖產(chǎn)生的波分布在較寬的波數(shù)k范圍內(nèi),頻散現(xiàn)象較為明顯,基階模態(tài)在該范圍內(nèi)的頻散導(dǎo)致尖端反射之后的強烈振蕩,當比率增加到8時,kmax≈0.125,在此范圍內(nèi),頻散曲線較為平坦,相速度趨于一維桿中波速,截面變形后仍近似平面,在此情況下,圓柱中遠場波近似為桿中波。圖8為Wp/R=4、8情況下,位置r/R=0.6處樁頂質(zhì)點速度響應(yīng)與基于一維桿理論得到的響應(yīng)比較??梢钥闯?,當樁頂測點布置在r/R=0.6附近且敲擊脈沖特征波長與樁徑比值Wp/R較大時,遠場反射信號可近似基于一維波動理論分析。 3 樁土相互作用影響分析 上述分析表明,當敲擊脈沖特征波長相對半徑較大時,遠場波傳播可近似用一維應(yīng)力波來分析,樁體向下傳播應(yīng)力波在樁土相互作用下,部分能量向周圍土體散射,部分能量沿樁體向上傳播,這導(dǎo)致向下傳播的波不斷衰減,衰減與樁幾何參數(shù)及樁土剪切波速有關(guān)。通過數(shù)值模擬側(cè)重研究樁土剪切波速對波傳播影響,數(shù)值模擬樁長為5 m,樁身直徑為0.2 m,樁周土外徑為1 m,為模擬水平無限土體,避免人工土邊界反射波影響,在土的外側(cè)施加無反射邊界。對低應(yīng)變測試,樁土之間采用彈性接觸。根據(jù)以上分析,為避免近場三維波動現(xiàn)象,對距樁軸線2R/3、樁頂1.5 m以下不同位置質(zhì)點縱向振動速度幅值進行分析,為了比較樁土相互作用影響,以1.5 m位置的質(zhì)點速度幅值為基礎(chǔ),對其他位置點質(zhì)點速度幅值進行歸一化處理。 3.1 樁身波速與土體波速 3.1.1 樁身波速 取土的剪切波速為50 m/s,密度為1 580 kg/m3,泊松比為0.44。樁身材料波速分別取為2 500、2 000、1 500 m/s,通過對數(shù)值模擬計算得到的響應(yīng)分析,得到樁身質(zhì)點縱向振動速度衰減曲線,如圖9所示。 由圖9可以看出,樁身剪切波速不同,其樁身質(zhì)點縱向振動速度衰減速率也隨之發(fā)生變化,且隨著樁體材料剪切波速的減小,衰減越快。即認為樁體材料剪切波速越小,應(yīng)力波在樁身傳播時的衰減速率越快。 3.1.2 土體波速 取土的密度為1 580 kg/m3,泊松比為0.44,樁體材料波速為3 000 m/s,改變樁周土體剪切波速,分別取50、60、100 m/s。衰減規(guī)律如圖10所示。 樁周土體剪切波速不同,其樁身質(zhì)點縱向振動速度衰減速率也不同,且隨著土體剪切波速的增大,衰減越快。即認為樁周土體剪切波速越大,應(yīng)力波在樁身傳播時的衰減越快。 3.1.3 保持剪切波速比值不變 取土的密度為1 580 kg/m3,泊松比為0.44,樁身剪切波速分別取2 500、2 000、1 500 m/s,對應(yīng)的樁周土體剪切波速分別取為100、80、60 m/s,從而保證3種情況下剪切波速比值為25,進行數(shù)值模擬分析,如圖11。 從圖11中可以看出,雖然樁土剪切波速比值不變,但不同樁土剪切波速導(dǎo)致衰減不同,樁體剪切波速越低,則衰減越快,這表明樁土剪切波速比值并不是應(yīng)力波衰減的決定性參數(shù)。 3.2 土體泊松比 取土的剪切波速為50 m/s,密度為1 580 kg/m3; 樁身材料波速為3 000 m/s。改變樁周土體泊松比,分別為0.43、0.40、0.35,進行低應(yīng)變數(shù)值模擬分析,得到曲線如圖12所示。 從圖12可以看出,樁周土體泊松比不同,其樁身質(zhì)點縱向振動速度衰減速率幾乎一致。即認為樁周土體泊松比對應(yīng)力波在樁中的傳播衰減速率幾乎無影響。 以上分析表明,樁體剪切波速與應(yīng)力波的衰減速率負相關(guān),樁周土體剪切波速與應(yīng)力波的衰減速率正相關(guān),土體泊松比對應(yīng)力波衰減速率幾乎無影響。 4 模型實驗 縮徑模型樁長3.3 m,幾何參數(shù)見圖13,縮徑處完整系數(shù)為0.36。按設(shè)計尺寸將模型樁制作好,然后埋入預(yù)先挖好的坑中,經(jīng)過一段時間固結(jié)后,進行低應(yīng)變測試,測試信號見圖14,樁體平均波速為3 654 m/s,為了考慮樁土相互作用的影響,取土的密度1 800 kg/m3,剪切波速150 m/s?;谝痪S波動理論,當完整系數(shù)取0.45時,測試與計算曲線達到最佳匹配,如圖14所示。由于土體擠入模型樁縮徑處導(dǎo)致完整系數(shù)比預(yù)設(shè)值大。這表明,可近似采用一維波動理論對遠場反射波進行分析。 5 結(jié)論 1)在樁頂處施加脈沖作用,距源較近,樁頂附近波場呈三維波動效應(yīng),波場中含P、S、R成份,通過對這些波的波陣面、能量以及波在樁側(cè)面反射分析可以研究近場波動響應(yīng);距振源較遠的遠場,可以用圓桿中縱向波傳播特性進行分析,當敲擊脈沖特征波長相對樁徑較大時,遠場縱向波相速度趨于桿波速,截面變形均勻,在此情況下,遠場波近似于一維桿中波。 2)樁頂距振源約0.6R處,側(cè)面反射波相消相干,側(cè)面多次反射波影響程度最小,當測點布置在此位置,且敲擊脈沖特征波長相對樁徑較大,可以近似基于一維波動理論對遠場反射波信號進行分析,但不適合對近場樁體反射波進行分析。 3)應(yīng)力波在樁土相互作用下衰減不只與樁土剪切波速比值有關(guān),更重要的是與樁體、土體剪切波速具體值有關(guān),樁體剪切波速越低,土體剪切波速越高,應(yīng)力波衰減越快。 參考文獻: [1]? ? GRAFF K F.Wave motion in elastic solids [M].Columbus: Ohio State University Press, 1976: 245. [2]? ?LIAO S T, ROESSET J M.Dynamic response of intact piles to impulse loads [J].International Journal for Numerical and Analytical Methods in Geomechanics, 1997, 21(4): 255-275. [3]? ?陳凡, 王仁軍.尺寸效應(yīng)對基樁低應(yīng)變完整性檢測的影響[J].巖土工程學(xué)報, 1998,20(5): 95-99. CHEN F, WANG R J.Dimension effect on low strain integrity testing of piles [J].Chinese Journal of Geotechnical Engineering, 1998, 20(5): 95-99.(in Chinese) [4]? ? CHOW Y K, PHOON K K, CHOW W F, et al.Low strain integrity testing of piles: Three-dimensional effects? ?[J].Journal of Geotechnical and Geoenvironmental Engineering, 2003, 129(11): 1057-1062. [5]? ?ZHENG C J, KOURETZIS G P, DING X M, et al.Three-dimensional effects in low-strain integrity testing of piles: analytical solution [J].Canadian Geotechnical Journal, 2016, 53(2): 225-235. [6]? ?王雪峰, 陳培, 甘軼.樁頂動測三維效應(yīng)分析[C]//第二屆全國巖土與工程學(xué)術(shù)大會論文集,武漢, 2006: 506-513. WANG X F, CHEN P, GAN Y.Three-dimensional effect analysis of pile top dynamic measurement[C] // Proceedings of the 2nd National Conference on Geotechnical and Engineering, Wuhan, 2006: 506-513. [7]? ?陳安國, 劉東甲.基樁低應(yīng)變檢測的波動方程法數(shù)值計算[J].合肥工業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2007, 30(10): 1336-1340. CHENA G, LIU D J.Numerical computation of the wave equation for low strain test piles [J].Journal of Hefei University of Technology(Natural Science), 2007, 30(10): 1336-1340.(in Chinese) [8]? ?岳向紅.基于三維導(dǎo)波理論的基樁和錨桿無損檢測技術(shù)研究[D].武漢: 中國科學(xué)院研究生院(武漢巖土力學(xué)研究所), 2008. YUE X H.Research of pile and bolts non-destructive test technology based on 3D guided-wave theory [D].Wuhan: Graduate School of Chinese Academy of Sciences (Wuhan Institute of Rock and Soil Mechanics), 2008.(in Chinese) [9]? ?陳輝, 董承全, 張佰戰(zhàn), 等.低應(yīng)變法檢測中幾種典型缺陷樁的理論模型分析[J].工程地球物理學(xué)報, 2012, 9(3): 342-345. CHENH, DONG C Q, ZHANG B Z, et al.Theoretical model analysis of defective piles in low strain method [J].Chinese Journal of Engineering Geophysics, 2012, 9(3): 342-345.(in Chinese) [10]? ?CHAI H Y, PHOON K K.Detection of shallow anomalies in pile integrity testing [J].International Journal of Geomechanics, 2013, 13(5): 672-677. [11]? ?CHAI H Y, WEI C F, PHOON K K, et al.Some observations on the performance of the signal matching technique in assessment of pile integrity [J].Journal of Nondestructive Evaluation, 2011, 30(4): 246-258. [12]? ?CHAI H Y, PHOON K K, ZHANG D J.Effects of the source on wave propagation in pile integrity testing [J].Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2010, 136(9): 1200-1208. [13]? ?吳文兵, 蔣國盛, 王奎華, 等.土塞效應(yīng)對管樁縱向動力特性的影響研究[J].巖土工程學(xué)報, 2014, 36(6): 1129-1141. WU W B, JIANG G S, WANG K H, et al.Influence of soil plug effect on vertical dynamic response of pipe piles [J].Chinese Journal of Geotechnical Engineering, 2014, 36(6): 1129-1141.(in Chinese) [14]? ?CHAMMAS R, ABRAHAM O, COTE P, et al.Characterization of heterogeneous soils using surface waves: Homogenization and numerical modeling [J].International Journal of Geomechanics, 2003, 3(1): 55-63. [15]? ?UENISHI K.On a possible role of Rayleigh surface waves in dynamic slope failures [J].International Journal of Geomechanics, 2010, 10(4): 153-160. [16]? ?PAVLAKOVIC N B.Leaky guided ultrasonic waves in NDT [D].London: Imperial College London, UK, 1998. (編輯 王秀玲) 收稿日期:2019-07-01 基金項目:? 國家自然科學(xué)基金(41474113) 作者簡介:? 丁瑋(1995- ),女,主要從事表面波及應(yīng)力波的測試研究,E-mail:yueyinding@outlook.com。 柴華友(通信作者),博士,副研究員,E-mail:chy_rsm@hotmail.com。 Received: 2019-07-01 Foundation items:? National Natural Science Foundation of China (No. 41474113) Author brief:? Ding Wei (1995- ), main research interest: test research on surface wave and stress wave, E-mail: yueyinding@outlook.com. Chai Huayou (corresponding author), PhD, associate professor, E-mail: chy_rsm@hotmail.com.