基于四維離散數(shù)值方法的巖石圓環(huán)試樣動態(tài)破壞及耗能規(guī)律

張奔 趙高峰

摘 要：利用最新提出的四維離散彈簧元法（4D-LSM）對巖石圓環(huán)的動態(tài)抗變形及耗能特性進(jìn)行研究。首先對4D-LSM用于描述巖石圓環(huán)試樣破壞的適用性進(jìn)行了驗(yàn)證，發(fā)現(xiàn)4D-LSM可以較好地再現(xiàn)巖石試樣的孔徑比與破壞形態(tài)和強(qiáng)度間的影響規(guī)律。在此基礎(chǔ)上，對孔隙率、非均質(zhì)性、厚徑比等因素與抗變形及耗能間的規(guī)律進(jìn)行研究，得到了對應(yīng)的數(shù)學(xué)公式。利用4D-LSM模擬大變形方面的優(yōu)勢，發(fā)現(xiàn)當(dāng)巖石材料的變形抵抗能力較大時，巖石圓環(huán)將表現(xiàn)出與傳統(tǒng)僅考慮小變形情況下實(shí)驗(yàn)和數(shù)值計(jì)算不同的破壞形態(tài)，巖石圓環(huán)的材料抗拉強(qiáng)度與抗變形能力也分別呈現(xiàn)非線性關(guān)系。通過建立圓環(huán)陣列模型，研究了由多個巖環(huán)構(gòu)成的組合結(jié)構(gòu)的抗變形和能量抵抗能力以及對應(yīng)的破壞形態(tài)。結(jié)果表明，組合結(jié)構(gòu)中的圓環(huán)單元與單圓環(huán)受力的破裂形態(tài)有所不同，但單圓環(huán)的抗變形及耗能規(guī)律仍適用于巖環(huán)陣列結(jié)構(gòu)的分析。

關(guān)鍵詞：四維離散彈簧元法;巖石圓環(huán);耗能;變形

中圖分類號：TU457

文獻(xiàn)標(biāo)志碼：A 文章編號：2096-6717（2019）02-0020-09

Abstract：In this work， the deformation resistance and energy dissipation of rock rings were investigated by using the 4D lattice spring model （4D-LSM）. Firstly， the 4D-LSM was verified against existing experimental and numerical simulation results in terms of the influence of rock ring size on the failure pattern and strength. On this basis， empirical relationships of the influence of various factors， e.g. the porosity， heterogeneity， and thickness-diameter ratio， on the deformation resistance and energy dissipation of the rock ring were derived from numerical tests with the 4D-LSM. When considering the large deformation， it was found that the failure pattern of the rock ring would be different from the classical experimental and numerical observations. Moreover， different empirical relationships between the deformation resistance and energy dissipation as well as the material ultimate deformation were derived compared to these numerical tests considering small deformation only. Finally， the failure pattern， deformation resistance and energy dissipation of a composed structure made up from an array of rock rings were studied. Numerical results reveal that these relationships obtained for the single rock ring are also applicable to the composed structure.

Keywords：4D lattice spring model; rock ring; energy dissipation; deformation

近年來，地下空間不斷得到新的拓展和開發(fā)，許多重大基礎(chǔ)設(shè)施，如地鐵、人防工程、核廢料處置庫等在地下相繼修建。長期以來，人們對地下工程的抗震性能都持較為樂觀的態(tài)度，認(rèn)為地下工程抗震能力要遠(yuǎn)優(yōu)于地上結(jié)構(gòu)，然而，最近一些地震活動導(dǎo)致的地下工程嚴(yán)重破壞乃至坍塌的案例[1]使人們開始對地下工程在地震作用下的穩(wěn)定性進(jìn)行更深入的研究。與地面建筑以抗震設(shè)計(jì)為主不同，地下重大基礎(chǔ)設(shè)施多采用更為有效的隔震手段[2]。地震屬于振動形式的一種，隔振技術(shù)主要通過有效的防護(hù)措施將建筑物和振動環(huán)境進(jìn)行適當(dāng)隔離，通過降低能量的輸入保證建筑結(jié)構(gòu)在載荷沖擊過程中盡量保持正常狀態(tài)[3]。實(shí)際工程中，常采用在重要地下硐室旁開挖輔助硐室的方式來隔離動力災(zāi)害，通過鉆孔的方式來改變巖體的動力學(xué)阻抗也是一種可行的隔震方式。

作為一種新型多功能材料，孔隙材料在能量吸收方面展現(xiàn)了良好的性能[4-8]，近年來備受關(guān)注。其中，類巖石孔隙材料由于質(zhì)量輕、耗能好、耐久性高等特點(diǎn)而被用于土建工程中，其吸能減振的機(jī)理受到了科研工作者的廣泛重視。王智等[9]對多孔混凝土的制備方法進(jìn)行探究，在保證孔隙率的條件下有效提高了混凝土的強(qiáng)度。黃勝等[10]、趙武勝等[11]將自制的泡沫混凝土作為西藏扎墨公路嘎隆拉隧道的隔震層進(jìn)行模擬計(jì)算，破壞區(qū)的范圍得到有效縮小，證明了孔隙材料保護(hù)隧道免遭地震波沖擊的良好性能。王代華等[12]對含有泡沫混凝土層的復(fù)合結(jié)構(gòu)在爆破中的能量狀態(tài)進(jìn)行模擬，結(jié)果顯示，孔隙結(jié)構(gòu)對于改變復(fù)合結(jié)構(gòu)中的能量分布效果明顯，能量穿過孔隙結(jié)構(gòu)后明顯衰減。劉海燕等[13]對泡沫混凝土高應(yīng)變率狀態(tài)下的吸能特征進(jìn)行了實(shí)驗(yàn)探究，總結(jié)了泡沫混凝土在抗沖擊過程中的變形規(guī)律。葉燕華等[14]在空心砌塊墻體中分別注入、不注入泡沫混凝土，并反復(fù)加、卸壓，對比發(fā)現(xiàn)，前者在載荷作用下抗剪強(qiáng)度和極限承載力得到提高，墻體破壞時裂紋擴(kuò)展分散、速度慢，顯示了優(yōu)良的抗倒塌性能。

巖石類材料以脆性為主，具有較高的強(qiáng)度，但抗變形能力很小?？紫兜囊雽r石類材料特性改變非常大，能提高其隔振能力，但動力學(xué)數(shù)值模擬中對剛性轉(zhuǎn)角大變形的分析一直比較困難，因?yàn)閿?shù)學(xué)上需要進(jìn)行近似處理。最近提出的四維離散彈簧元法（4D-LSM）通過四維空間的相互作用，一定程度上解決了這個問題[15]，該方法可以較好地處理動態(tài)大變形問題。筆者先對4D-LSM用于單個巖石圓環(huán)模型的適用性進(jìn)行了驗(yàn)證;研究了單個巖環(huán)模型在不同條件下的破裂形態(tài)以及動態(tài)變形和能量抵抗能力;對含有隨機(jī)孔隙結(jié)構(gòu)的巖環(huán)陣列模型的動態(tài)能量吸收特性和尺寸效應(yīng)進(jìn)行了分析。

1 四維離散數(shù)值方法4D-LSM

1.1 模型簡介

4D-LSM是在三維離散彈簧元法（DLSM）的基礎(chǔ)上發(fā)展的。在DLSM[16]中，模型由彈簧和離散球形顆粒連接構(gòu)成，顆粒之間的連接彈簧包括法向和切向兩個方向。DLSM使用相同的顆粒模型可以產(chǎn)生不同的晶格結(jié)構(gòu)，從而表征不同的宏觀力學(xué)行為。與傳統(tǒng)離散方法相比，DLSM無需宏觀參數(shù)校準(zhǔn)，自由度是離散單元法（DEM）的一半，且容易進(jìn)行并行處理。DLSM中切向彈簧根據(jù)局部應(yīng)變進(jìn)行計(jì)算，具有旋轉(zhuǎn)不變性。4D-LSM[15]基于平行世界的概念，在DLSM的基礎(chǔ)上引入第四維度構(gòu)造“四維超薄膜”，通過三維模型與平行世界中的模型建立相互作用關(guān)系，有效解決了傳統(tǒng)LSM方法中的泊松比限制問題，4D-LSM的原理圖見文獻(xiàn)[15]。相對于經(jīng)典連續(xù)介質(zhì)力學(xué)，4D-LSM的離散性使其更適合模擬巖石或類巖石等非均質(zhì)材料的破壞問題。

1.2 理論公式

與DLSM相比，4D-LSM的力和法向矢量具有4個分量。在4D-LSM中，代表三維相互作用的彈簧剛度（k3D）全部相同，而代表四維相互作用的彈簧剛度出現(xiàn)差異，差異特征由四維剛度比λ4D表征。為了獲得彈性的各向同性，四維相互作用中彈簧剛度通過式（1）進(jìn)行調(diào)整。

2 巖石圓環(huán)模型計(jì)算驗(yàn)證

2.1 計(jì)算模型

建立半徑R為50 mm、厚度T為50 mm的圓盤，組成圓盤的球形顆粒直徑為1 mm。探究4D-LSM對于研究圓環(huán)劈裂的適用性，以10 mm為一個梯度，在圓盤模型中分別挖去半徑為10、20、30 mm的同心圓孔，如圖1所示。用參數(shù)δ來表征內(nèi)外徑之比，即δ=r/R（δ=0、0.2、0.4、0.6）。模型彈性模量、泊松比、密度的設(shè)置參考卡拉拉大理石，系統(tǒng)具體輸入?yún)?shù)如表1所示。

因巴西圓盤劈裂實(shí)驗(yàn)中加載端存在應(yīng)力集中問題，有學(xué)者對試樣進(jìn)行改進(jìn)[17]，將加載端打磨成兩個相互平行的平面，使線載荷轉(zhuǎn)化為面載荷，有效地保證了圓盤試樣的中間起裂。對于4D-LSM模擬劈裂問題，模型采用速度控制的加載方式，速度作用在模型最上層顆粒面和最底層顆粒面，加載角度φ為16.1°。對所有模型最上層的顆粒面設(shè)置大小為4 mm/s的垂直向下速度;最底層顆粒面沿垂直方向速度設(shè)為0，表示對模型的底端的固定。

2.2 驗(yàn)證

將δ=0時的圓盤模型考慮成圓環(huán)模型的特殊情況。由圖1可知，當(dāng)δ值為0、0.2時，模型表現(xiàn)出單一的縱向拉伸破壞，試樣破壞成兩部分;δ值等于0.4、0.6時，模型同時表現(xiàn)出縱向和橫向拉伸破壞，模型破壞成4部分。這表明，隨著δ值的增大，模型破壞裂紋數(shù)量表現(xiàn)出增多的趨勢。由圖1可知，通過4D-LSM模擬巖石圓環(huán)模型的加載，得到的破壞形態(tài)與前人[18-19]得出的物理實(shí)驗(yàn)結(jié)果相吻合。

將模型加載峰值載荷除去πRL進(jìn)行標(biāo)準(zhǔn)化處理，即抗拉強(qiáng)度σt=Pmax/πRL。為對4D-LSM求解的圓環(huán)模型的抗拉強(qiáng)度與現(xiàn)有文獻(xiàn)中的結(jié)果進(jìn)行對比，將各試樣的σt除去δ=0時的σt，進(jìn)行歸一化處理。圖2為歸一化處理后的實(shí)驗(yàn)數(shù)據(jù)[18]、FEM/DEM模擬數(shù)據(jù)[20]以及4D-LSM模擬數(shù)據(jù)和擬合曲線，由圖2可知數(shù)據(jù)相差不大。4D-LSM數(shù)據(jù)擬合曲線的公式為

巖環(huán)模型隨著孔徑比δ的增加抗拉強(qiáng)度呈指數(shù)下降的趨勢，δ對模型抗拉強(qiáng)度的影響隨著δ的增加逐漸減弱。擬合曲線有效驗(yàn)證了尤明慶等[21]指出的巖環(huán)劈裂載荷隨內(nèi)徑增加而呈現(xiàn)出指數(shù)下降的結(jié)論。通過破壞形態(tài)、模擬數(shù)據(jù)兩方面證實(shí)了4D-LSM對于研究巖石圓環(huán)力學(xué)特性的可行性。

3 巖石圓環(huán)的耗能規(guī)律

3.1 厚度與外直徑之比

將模型厚度T與外直徑D的比值定義為w，即w=T/D，選取δ=0.6的模型，通過設(shè)置不同厚度值，探究w（w=0.25、0.50、0.75、1.00）對模擬結(jié)果的影響，如圖3所示。對孔隙率n、均質(zhì)度m、顆粒間彈簧抗拉伸長量u 3個因素以δ=0.6、w=0.50的圓環(huán)模型為基礎(chǔ)進(jìn)行研究。

圓環(huán)破壞的載荷加載點(diǎn)位移曲線中，峰值載荷對應(yīng)的位移值為模型的動態(tài)變形抵抗能力;在動態(tài)載荷的作用下，破壞后吸收的能量為模型的動態(tài)能量抵抗能力，數(shù)值等于載荷加載點(diǎn)位移曲線中開始加載到峰值載荷的曲線段與橫坐標(biāo)軸圍成圖形的面積。因該曲線段接近直線，圖形近似看成三角形，則模型的動態(tài)能量抵抗能力計(jì)算式為

式中：Q為動態(tài)能量抵抗能力;z為動態(tài)變形抵抗能力。

由圖4（a）可知，不同w值模型的峰值載荷差異較大，但整理數(shù)據(jù)后發(fā)現(xiàn)，圓環(huán)模型抗拉強(qiáng)度并未因?yàn)槟Ｐ秃穸炔煌霈F(xiàn)太大變化，均值為0.230 0 MPa。這表明厚度的增加對于抗拉強(qiáng)度不會帶來較大的影響。模型的峰值載荷與模型厚度在一定范圍內(nèi)成正比例關(guān)系，可以通過調(diào)節(jié)厚度的大小對模型峰值載荷進(jìn)行調(diào)整。隨著w值的改變，動態(tài)變形抵抗能力在均值13.360 0 μm上下波動，而能量抵抗能力與w值呈正比例關(guān)系，函數(shù)表達(dá)式分別為

3.2 孔隙率

巖石中各類孔洞、裂隙的總體積與巖石總體積的比值稱為巖石的孔隙率[22]。分別生成孔隙率為0.05、0.10、0.15、0.20的模型，研究孔隙率對圓環(huán)劈裂的影響，如圖5所示。從圖6（b）中可知，隨著孔隙率的增加，抗拉強(qiáng)度呈二次函數(shù)下降，擬合曲線為

對圓環(huán)模型的動態(tài)變形抵抗能力和動態(tài)能量抵抗能力與孔隙率的關(guān)系進(jìn)行擬合，如圖6（c）、（d）所示，圓環(huán)模型抵抗變形的能力與動態(tài)能量抵抗能力隨著孔隙率的增加分別呈冪函數(shù)和對數(shù)關(guān)系下降，表達(dá)式分別為

3.3 非均質(zhì)性

由于巖石礦物組成成分、粒徑以及分布方式的不同，再加上裂隙、節(jié)理、孔隙等結(jié)構(gòu)的存在，導(dǎo)致自然界中巖石具有非均質(zhì)的特性。在4D-LSM系統(tǒng)中，模型材料的均質(zhì)性符合韋布爾分布函數(shù)，用參數(shù)m表征材料的均質(zhì)度，m值越大均質(zhì)度越高。對圓環(huán)模型m值分別設(shè)置為5、10、20、30進(jìn)行圓環(huán)劈裂實(shí)驗(yàn)，如圖7所示。載荷加載點(diǎn)位移曲線如圖8（a）所示，峰值載荷隨均質(zhì)度的增加逐漸增大。抗拉強(qiáng)度與均質(zhì)度的擬合式為

模型的動態(tài)變形抵抗能力和動態(tài)能量抵抗能力隨著m值的增加分別呈對數(shù)關(guān)系增加，表達(dá)式分別為

3.4 脆性小變形破壞

將圓環(huán)模型顆粒間彈簧抗拉伸長量u分別設(shè)為0.10、0.15、0.20、0.30 μm，進(jìn)一步探究脆性階段圓環(huán)破壞與u的關(guān)系。圓環(huán)模型載荷加載點(diǎn)位移曲線如圖9（a）所示，4條曲線表現(xiàn)出了極高的相似性，說明顆粒間彈簧抗拉伸長量的改變對載荷加載點(diǎn)位移曲線的走勢影響不大。抗拉強(qiáng)度與顆粒間彈簧抗拉伸長量的擬合線如圖9（b）所示，直線方程為

模型的動態(tài)變形抵抗能力和動態(tài)能量抵抗能力隨著u值的增加分別呈正比例函數(shù)和二次函數(shù)關(guān)系增加，表達(dá)式分別為

3.5 大變形破壞

將u分別設(shè)為0.04、0.08、0.12、0.16 mm，對圓環(huán)模型大變形階段的破壞模式和耗能特征進(jìn)行研究，如圖10所示。結(jié)果表明，大變形階段圓環(huán)模型抗拉強(qiáng)度擬合曲線區(qū)別于小變形階段的正比例函數(shù)，更接近二次函數(shù)，函數(shù)關(guān)系式為

大變形階段圓環(huán)模型動態(tài)變形抵抗能力與小變形階段正比例函數(shù)相比也轉(zhuǎn)化為冪函數(shù)關(guān)系，擬合公式為

與小變形階段相比大變形階段圓環(huán)模型的動態(tài)能量抵抗能力表現(xiàn)出了相同的二次函數(shù)規(guī)律，函數(shù)表達(dá)式為Q=481.670 0u2+20.356 0u-0.648 5，

圖11為大變形階段（u=0.16 mm）圓環(huán)模型破壞過程圖。通過對u的設(shè)置提高模型變形量，以對圓環(huán)的大變形破壞過程進(jìn)行模擬，研究發(fā)現(xiàn)，與脆性小變形階段相比，不僅圓環(huán)模型動態(tài)變形抵抗能力和動態(tài)能量抵抗能力增加，圓環(huán)模型的破壞模式也發(fā)生根本轉(zhuǎn)變：從小變形階段破壞成4部分，演變成僅縱向發(fā)生拉伸破壞。

3.6 巖環(huán)組合體

以δ=0.6、w=0.5、n=0.15的圓環(huán)模型為單元，分別建立N×N（N=1、2、3、4）的圓環(huán)陣列模型，如圖12所示，通過N對組合體模型尺度進(jìn)行表征。對模型抗變形能力進(jìn)行研究，對比不同組合體模型的載荷加載點(diǎn)位移曲線發(fā)現(xiàn)，隨著N的增加，圓環(huán)陣列模型的峰值載荷也逐漸變大，但第2峰值載荷與第1峰值載荷相比逐漸弱化。圓環(huán)陣列模型的第1峰值載荷Pmax與N呈正比例函數(shù)關(guān)系，表達(dá)式為

對巖環(huán)組合體模型的動態(tài)變形抵抗能力和動態(tài)能量抵抗能力進(jìn)行分析，如圖13（c）、（d）所示，發(fā)現(xiàn)巖環(huán)組合體的變形抵抗能力隨著N的增加呈正比例增加，而動態(tài)變形抵抗能力與N呈冪函數(shù)關(guān)系。函數(shù)表達(dá)式分別為

對組成陣列材料圓環(huán)單元的動態(tài)變形抵抗能力和動態(tài)能量抵抗能力進(jìn)行分析，如圖13（e）、（f）所示。圓環(huán)單元的動態(tài)變形抵抗能力和動態(tài)能量抵抗能力隨著N的變化基本沒有改變。單位圓環(huán)的變形抵抗能力均值為10.799 0 μm，能量抵抗能力均值為6.129 3×10-3 J，可知圓環(huán)單元在能量吸收方面幾乎不受陣列模型尺度的影響。

4 討論

根據(jù)不同因素下巖石圓環(huán)模型的破裂形態(tài)可知，孔隙率、非均質(zhì)度均能引起圓環(huán)模型裂紋擴(kuò)展形態(tài)的改變，但總體來說豎向裂紋差異較小，橫向裂紋大體在圓環(huán)腰部水平線附近上下浮動;厚徑比對裂紋擴(kuò)展有一定影響，主要表現(xiàn)為橫向次生裂紋的數(shù)量差異。通過對孔隙度或均質(zhì)度設(shè)置合理的參數(shù)，圓環(huán)模型破壞形態(tài)表現(xiàn)出了較規(guī)則的“四扇葉形”，如圖5、圖7所示，與前人實(shí)驗(yàn)中所得破壞形式吻合[18-19，23]。由此可見，通過對圓環(huán)模型的均質(zhì)度或孔隙率設(shè)置合理的參數(shù)，4D-LSM可以對巖石圓環(huán)模型的破壞模式進(jìn)行良好的再現(xiàn)。

圓環(huán)模型厚度的增加并不會帶來抗拉強(qiáng)度和變形抵抗能力的增加，但能量抵抗能力隨著w值呈正比例增長。孔隙率的增加會導(dǎo)致圓環(huán)模型抗拉強(qiáng)度、變形抵抗能力與能量抵抗能力的不同規(guī)律下降，但總體而言，對變形抵抗能力影響幅值較小，抗拉強(qiáng)度、能量抵抗能力降低幅值較大?？估瓘?qiáng)度、變形抵抗能力、能量抵抗能力隨著模型均質(zhì)度的增加而出現(xiàn)不同規(guī)律增大，由此可知，非均質(zhì)性是影響巖石強(qiáng)度的一個重要因素，通過均質(zhì)度的提高可以有效改善圓環(huán)模型的隔振性能。與小變形階段抗拉強(qiáng)度與變形抵抗能力的線性增長相比，大變形階段得到了非線性的規(guī)律和不同于現(xiàn)有實(shí)驗(yàn)和數(shù)值的破壞結(jié)果。由此可知，當(dāng)對類巖石材料（如混凝土等）通過加強(qiáng)材料韌性如添加納米材料進(jìn)行改性時，其對應(yīng)圓環(huán)模型的耗能與抗變形能力將會有大的提升，破壞模式與脆性階段相比也將表現(xiàn)出較大差異。

巖環(huán)組合體模型破壞形式主要表現(xiàn)為45°斜向破壞帶貫穿模型，導(dǎo)致模型喪失承載能力。組合體模型雖然總體上承受壓力，但本質(zhì)上仍表現(xiàn)出圓環(huán)單元的受拉破壞，圓環(huán)單元破壞成4部分。不同尺度組合體模型中，各圓環(huán)單元表現(xiàn)出不同的破壞形態(tài)，但圓環(huán)單元的變形特征和耗能規(guī)律卻與單個圓環(huán)受載荷作用時表現(xiàn)一致。由此可見，對巖環(huán)組合體的耗能規(guī)律可以通過單個圓環(huán)進(jìn)行較好的表征。基于以上特點(diǎn)，當(dāng)通過增加孔隙對隔振超材料進(jìn)行設(shè)計(jì)時，可先對單個圓環(huán)著手研究。

5 結(jié)論

1）4D-LSM用于分析巖石圓環(huán)破壞的適用性得到驗(yàn)證。模擬結(jié)果表明，圓環(huán)試樣的的破壞形態(tài)由多個因素相互耦合共同作用而成。非均質(zhì)性、孔隙率對圓環(huán)模型的破壞形式有著重要影響。

2）分別擬合得出了抗拉強(qiáng)度、動態(tài)變形抵抗能力、動態(tài)能量抵抗能力與厚徑比、孔隙率、非均質(zhì)度、顆粒間彈簧抗拉伸長量的函數(shù)關(guān)系式?？偨Y(jié)了圓環(huán)模型不同因素的抗變形特征和耗能規(guī)律。

3）N×N巖石圓環(huán)陣列模型的承載強(qiáng)度與N呈正比例關(guān)系，組合體模型中圓環(huán)單元的動態(tài)變形抵抗能力與動態(tài)能量抵抗能力不受陣列模型尺度的影響。

參考文獻(xiàn)：

[1] 陳國興，陳蘇，杜修力，等. 城市地下結(jié)構(gòu)抗震研究進(jìn)展[J]. 防災(zāi)減災(zāi)工程學(xué)報，2016，36（1）：1-23.

CHEN G X，CHEN S，DU X L，et al. Review of seismic damage， model test， available design and analysis methods of urban underground structures： Retrospect and prospect [J]. Journal of Disaster Prevention and Mitigation Engineering， 2016，36（1）： 1-23. （in Chinese）

[2] 吳畏，蔣美蓉，石磊，等. 地下工程隔震技術(shù)應(yīng)用探討[J]. 地下空間與工程學(xué)報，2007，3（Sup2）： 1583-1586.

WU W， JIANG M R， SHI L， et al. Discussion on application of shock-isolation technology in underground engineering [J]. Chinese Journal of Underground Space Engineering， 2007， 3（Sup2）： 1583-1586. （in Chinese）

[3] 孫鐵成，高波，葉朝良. 地下結(jié)構(gòu)抗震減震措施與研究方法探討[J]. 現(xiàn)代隧道技術(shù)，2007， 44（3）： 1-5.

SUN T C， GAO B， YE C L. Discussion on anti-seismic and seismic-relieving measures and corresponding research methods for underground structures [J]. Modern Tunnelling Technology， 2007， 44（3）： 1-5. （in Chinese）

[4] 喻寅， 賀紅亮， 王文強(qiáng)，等. 多孔脆性材料對高能量密度脈沖的吸收和抵抗能力[J]. 物理學(xué)報， 2015， 64（12）： 263-271.

YU Y， HE H L，WANG W Q， et al. The ability of porous brittle materials to absorb and withstand high energy density pulse [J]. Acta Physica Sinica， 2015， 64（12）： 263-271. （in Chinese）

[5] HANSSEN A G， ENSTOCK L， LANGSETH M. Close-range blast loading of aluminium foam panels [J]. International Journal of Impact Engineering， 2002， 27（6）： 593-618.

[6] RATHBUN H J， RADFORD D D， XUE Z， et al. Performance of metallic honeycomb-core sandwich beams under shock loading [J]. International Journal of Solids and Structures， 2006， 43（6）： 1746-1763.

[7] RADFORD D D， MCSHANE G J， DESHPANDE V S， et al. The response of clamped sandwich plates with metallic foam cores to simulated blast loading [J]. International Journal of Solids and Structures， 2006， 43（7/8）： 2243-2259.

[8] 王建軍，郭偉國，趙融，等. 輕質(zhì)泡沫混凝土吸能特性及其壓潰流動應(yīng)力方程[J]. 土木建筑與環(huán)境工程，2013， 35（6）： 96-102.

WANG J J， GUO W G， ZHAO R， et al. Energy-absorbing properties and crushing flow stress equation of lightweight foamed concrete [J]. Journal of Civil， Architectural & Environmental Engineering， 2013， 35（6）： 96-102. （in Chinese）

[9] 王智， 錢覺時， 張朝輝，等. 多孔混凝土配合比設(shè)計(jì)方法初探[J]. 重慶建筑大學(xué)學(xué)報， 2008， 30（3）： 121-124.

WANG Z， QIAN J S， ZHANG Z H， et al. Preliminary discussion on mixing design of porous concrete [J]. Journal of Chongqing Jianzhu University， 2008， 30（3）： 121-124. （in Chinese）

[10] 黃勝，陳衛(wèi)忠，楊建平，等. 地下工程地震動力響應(yīng)及抗震研究[J]. 巖石力學(xué)與工程學(xué)報，2009，28（3）： 483-490.

HUANG S， CHEN W Z， YANG J P， et al. Research on earthquake-induced dynamic responses and aseismic measures for underground engineering [J]. Chinese Journal of Rock Mechanics and Engineering， 2009， 28（3）： 483-490. （in Chinese）

[11] 趙武勝，陳衛(wèi)忠，譚賢君，等. 高性能泡沫混凝土隧道隔震材料研究[J]. 巖土工程學(xué)報，2013，35（8）： 1544-1552.

ZHAO W S， CHEN W Z， TAN X J， et al. High-performance foam concrete for seismic-isolation materials of tunnels [J]. Chinese Journal of Geotechnical Engineering， 2013， 35（8）： 1544-1552. （in Chinese）

[12] 王代華，劉殿書，杜玉蘭，等. 含泡沫吸能層防護(hù)結(jié)構(gòu)爆炸能量分布的數(shù)值模擬研究[J]. 爆炸與沖擊，2006，26（6）： 562-567.

WANG D H， LIU D S， DU Y L， et al. Numerical simulation of anti-blasting mechanism and energy distribution of composite protective structure with foam concrete [J]. Explosion and Shock Waves， 2006， 26（6）： 562-567. （in Chinese）

[13] 劉海燕，李然. 泡沫混凝土吸能機(jī)理試驗(yàn)研究[J]. 成都大學(xué)學(xué)報（自然科學(xué)版）， 2010， 29（2）： 166-167，180.

LIU H Y， LI R. Experimental study on endergonic mechanism of foam concrete [J]. Journal of Chengdu University （Natural Science Edition）， 2010， 29（2）： 166-167，180. （in Chinese）

[14] 葉燕華， 李利群， 孫偉民，等. 內(nèi)注泡沫混凝土空心砌塊墻體抗震性能的試驗(yàn)研究[J]. 地震工程與工程振動， 2004， 24（5）： 154-158.

YE Y H， LI L Q， SUN W M， et al. Experimental study on seismic behaviors of hollow block wall filled with foaming concrete [J]. Earthquake Engineering and Engineering Vibration， 2004， 24（5）： 154-158. （in Chinese）

[15] ZHAO G F. Developing a four-dimensional lattice spring model for mechanical responses of solids [J]. Computer Methods in Applied Mechanics and Engineering， 2017， 315： 881-895.

[16] ZHAO G F. Development of micro-macro continuum-discontinuum coupled numerical method [D]. EPFL， Switzerland， 2010.

[17] WANG Q Z， XING L. Determination of fracture toughness KIC by using the flattened Brazilian disk specimen for rocks [J]. Engineering Fracture Mechanics， 1999， 64（2）： 193-201.

[18] LI D Y， WANG T， CHENG T J， et al. Static and dynamic tensile failure characteristics of rock based on splitting test of circular ring [J]. Transactions of Nonferrous Metals Society of China， 2016， 26（7）： 1912-1918.

[19] MELLOR M， HAWKES I. Measurement of tensile strength by diametral compression of discs and annuli [J]. Engineering Geology， 1971， 5（3）： 173-225.

[20] LI X， FENG F， LI D. Numerical simulation of rock failure under static and dynamic loading by splitting test of circular ring [J]. Engineering Fracture Mechanics， 2018，188：184-201.

[21] 尤明慶， 陳向雷， 蘇承東. 干燥及飽水巖石圓盤和圓環(huán)的巴西劈裂強(qiáng)度[J]. 巖石力學(xué)與工程學(xué)報， 2011， 30（3）： 464-472.

YOU M Q， CHEN X L， SU C D. Brazilian splitting strengths of discs and rings of rocks in dry and saturated conditions [J]. Chinese Journal of Rock Mechanics and Engineering， 2011， 30（3）： 464-472. （in Chinese）

[22] 王作棠，周華強(qiáng)，謝耀社. 礦山巖體力學(xué)[M]. 江蘇 徐州： 中國礦業(yè)大學(xué)出版社， 2007.

WANG Z T， ZHOU H Q， XIE Y S. Mine rock mechanics [M]. Xuzhou， Jiangsu： China University of Mining and Technology Press， 2007. （in Chinese）

[23] LI X B， WU Q H， TAO M， et al. Dynamic Brazilian splitting test of ring-shaped specimens with different hole diameters [J]. Rock Mechanics and Rock Engineering， 2016， 49（10）： 4143-4151.

（編輯 胡英奎）