胡啟軍+蔣晶+徐亞輝+石仁丹+張宇浩

摘要：樁巖（土）接觸面力學(xué)特性的研究是樁基承載機(jī)理研究的基礎(chǔ)。通過紅層泥巖樁巖接觸面大型直剪試驗(yàn)，研究了紅層泥巖樁巖接觸面的力學(xué)特性，結(jié)果表明：接觸面剪應(yīng)力先隨剪切位移增大而增大，在達(dá)到峰值后，剪應(yīng)力隨著剪切位移增大而降低，并最終趨于穩(wěn)定值，應(yīng)力應(yīng)變曲線呈現(xiàn)出應(yīng)變軟化的特征。根據(jù)剪切試驗(yàn)結(jié)果，推導(dǎo)出樁巖接觸面應(yīng)變軟化本構(gòu)方程。利用fish語言對(duì)FLAC3D中自帶的理想彈塑性接觸單元進(jìn)行二次開發(fā)，并應(yīng)用開發(fā)的模型對(duì)樁巖接觸面直剪試驗(yàn)進(jìn)行了數(shù)值模擬，分析剪應(yīng)力與剪切位移之間的關(guān)系，證明了該本構(gòu)能夠較好地模擬接觸面間的應(yīng)變軟化特性。

關(guān)鍵詞：紅層泥巖；樁巖接觸面；應(yīng)變軟化；本構(gòu)模型；數(shù)值模擬

中圖分類號(hào)：TU473.1文獻(xiàn)標(biāo)志碼：A文章編號(hào)：16744764（2017）03012207

Abstract：The pilerock（soil） interface mechanical behavior was the basis of pile foundation bearing mechanics. Firstly， the mechanical property of the interface was investigated through the laboratory large direct shear test of red mudstone pilerock interface. The results indicated that the shear stressstrain curve presented strain softening behavior. The shear stress was increasing with the shear strain initially， and when the shear stress came to a peak， the shear stress was decreasing with the shear strain and finally reached a steady. Subsequently， the constitutive equation was deduced based on the interface mechanical property， the ideal elasticplastic contact element in FLAC3D was further developed through fish language. Additionally， the direct shear test of red mudstone was simulated according to the interface constitutive model， and the relationship of shear stress and shear displacement between the pile and soil were analyzed which proved that the constitutive model showed good performance on simulating strain softening features of the interface.

Keywords：red mudstone； pilerock interface； constitutive models； strain softening； numerical simulation

目前，樁基礎(chǔ)在巖土工程中的各個(gè)領(lǐng)域得到了廣泛的應(yīng)用，而樁巖（土）接觸面的力學(xué)性質(zhì)對(duì)于樁基承載性能具有重要的影響。紅層泥巖是中國西南地區(qū)廣泛分布的一種軟弱巖類，它的強(qiáng)度介于巖與土之間，在建筑物荷載作用下，表現(xiàn)出變形大、承載力低、復(fù)雜易變等工程特性[12]。因此，研究紅層泥巖與樁基礎(chǔ)的相互作用機(jī)理具有重要的經(jīng)濟(jì)價(jià)值和社會(huì)意義。

樁巖（土）接觸面上兩種材料的變形不一致常常會(huì)引起接觸面上產(chǎn)生較大的剪應(yīng)力，而接觸面上的力學(xué)特性與巖（土）體的性質(zhì)有著密切的關(guān)系[34]。現(xiàn)階段，主要通過樁巖（土）接觸面的剪切試驗(yàn)來研究樁巖相互作用的力學(xué)特性[515]。文獻(xiàn)[5]用盒式直剪儀研究了土與混凝土接觸面的力學(xué)特性，建立了雙曲線本構(gòu)模型，因參數(shù)較少，得到了廣泛的應(yīng)用。文獻(xiàn)[610]通過細(xì)粒土與混凝土的剪切試驗(yàn)，對(duì)細(xì)粒土與混凝土接觸面的本構(gòu)模型進(jìn)行了研究。文獻(xiàn)[1114]根據(jù)粗粒土與混凝土的剪切試驗(yàn)結(jié)果，提出了粗粒土與混凝土接觸面的本構(gòu)模型。文獻(xiàn)[15]通過對(duì)泥巖與混凝土進(jìn)行剪切試驗(yàn)，修正了泥巖與混凝土的剪應(yīng)力退化模型。以上對(duì)樁巖（土）接觸面研究的試驗(yàn)對(duì)象主要還是針對(duì)砂土、粘土以及粗粒土，對(duì)于泥巖等軟巖的研究較少。本文通過紅層泥巖有無泥皮樁巖接觸面大型直剪試驗(yàn)，研究紅層泥巖樁巖接觸面力學(xué)特性。

有限差分軟件FLAC3D自帶的Interface接觸單元，可用于巖層接觸、樁土接觸等物體間的接觸、滑移、分離等問題[1618]。但它是一種基于理想彈塑性Coulomb滑動(dòng)本構(gòu)關(guān)系的無厚度接觸單元，不能反映接觸面復(fù)雜的力學(xué)特性[19]。因此，本文對(duì)FLAC3D中自帶接觸單元進(jìn)行二次開發(fā)，提高其對(duì)接觸面復(fù)雜力學(xué)特性的模擬能力。

本文進(jìn)行了樁巖接觸面大型直剪試驗(yàn)，并對(duì)接觸面的力學(xué)特性進(jìn)行研究，在此基礎(chǔ)上建立本構(gòu)方程，利用fish語言對(duì)FLAC3D中自帶的理想彈塑性接觸單元進(jìn)行二次開發(fā)，并運(yùn)用開發(fā)的模型對(duì)樁巖接觸面直剪試驗(yàn)進(jìn)行數(shù)值模擬，分析剪應(yīng)力與剪切位移之間的關(guān)系，對(duì)本構(gòu)模型的準(zhǔn)確性進(jìn)行驗(yàn)證。

1室內(nèi)樁巖接觸面剪切試驗(yàn)

1.1試驗(yàn)設(shè)備及試驗(yàn)材料

建立合理的接觸面本構(gòu)模型，需要深入地揭示接觸面在受荷過程中的主要力學(xué)規(guī)律和受力變形機(jī)理[12]，文獻(xiàn)[1115]利用樁巖（土）接觸面大型直剪試驗(yàn)，對(duì)樁巖（土）之間的力學(xué)規(guī)律進(jìn)行了研究。因此，本文進(jìn)行紅層泥巖樁巖接觸面直剪試驗(yàn)對(duì)接觸面的力學(xué)規(guī)律進(jìn)行研究。

試驗(yàn)設(shè)備采用西南石油大學(xué)土建院實(shí)驗(yàn)中心的QXZ5000微機(jī)控制電液伺服壓剪試驗(yàn)機(jī)，如圖1所示。該試驗(yàn)機(jī)豎向液壓加壓系統(tǒng)最大垂直荷載5 000 kN，水平剪力加載裝置最大水平荷載1 000 kN，采用電液伺服控制。計(jì)算機(jī)及分析軟件可繪制力位移曲線、位移時(shí)間曲線、力時(shí)間曲線等，自行設(shè)置控制參數(shù)，程序自動(dòng)運(yùn)行。

為了模擬符合現(xiàn)場(chǎng)實(shí)際的樁巖相互作用，在試驗(yàn)機(jī)基礎(chǔ)上自行設(shè)計(jì)并制作了上下剪切盒，在混凝土試塊和下壓板之間設(shè)置基槽，基槽由一個(gè)整塊鋼板鑲嵌在下壓板上，然后在鋼板上面處理得到5條凹槽，上面放置直徑12 mm的鋼珠。每次試驗(yàn)之前，在鋼板凹槽及滾珠上面涂上一層凡士林以減小剪切時(shí)的摩擦阻力。整個(gè)裝置示意圖如圖2所示。

紅層泥巖為樁側(cè)巖層，因此，在室內(nèi)樁巖接觸面剪切試驗(yàn)中，采用現(xiàn)場(chǎng)采集紅層泥巖試樣與混凝土試塊相互作用來進(jìn)行試驗(yàn)的研究以確保物質(zhì)材料的一致性。泥巖試塊尺寸設(shè)計(jì)為50 cm×50 cm×20 cm，試驗(yàn)采用室內(nèi)配置C35混凝土，混凝土試塊制作尺寸為60 cm×60 cm×20 cm，大于紅層泥巖試樣尺寸，以使在剪切過程中保持剪切面積不變。

1.2混凝土表面粗糙度設(shè)計(jì)

在實(shí)際樁基工程中，受樁側(cè)孔壁土體強(qiáng)度以及旋挖鉆機(jī)的影響，混凝土表面粗糙度有很大的隨機(jī)性。試驗(yàn)根據(jù)現(xiàn)場(chǎng)樁身表面粗糙度分布的隨機(jī)性，采用人工鑿毛法，對(duì)混凝土試塊表面進(jìn)行初步粗糙度處理，保證其表面有一定的粗糙度以符合工程實(shí)際。

參考已有學(xué)者[2021]對(duì)混凝土表面粗糙度進(jìn)行的劃分，定義混凝土表面粗糙度，y為灌砂平均深度，D為鑿毛最大深度限值，本實(shí)驗(yàn)設(shè)定為10 mm。灌砂平均深度y可由灌砂法求得，。

依據(jù)所得粗糙度n，將混凝土表面粗糙度劃分為3個(gè)等級(jí)：粗糙度Ⅰ為混凝土表面粗糙度n<01；粗糙度Ⅱ?yàn)榛炷帘砻娲植诙?.1≤n；粗糙度Ⅲ為混凝土表面粗糙度n>0.2。

試驗(yàn)人工鑿毛法對(duì)混凝土表面粗糙度處理結(jié)果如圖3所示。

e11.3試驗(yàn)工況及模型示意圖

試驗(yàn)工況如表1所示。

1.4試驗(yàn)結(jié)果與分析

1.4.1接觸面直剪試驗(yàn)曲線特征接觸面粗糙程度以及接觸面上作用的法向應(yīng)力對(duì)接觸面剪切力學(xué)特性都具有重要影響。對(duì)室內(nèi)樁巖接觸面剪切試驗(yàn)結(jié)果進(jìn)行整理，作出在不同接觸面粗糙程度、不同法向應(yīng)力作用下樁巖接觸面剪應(yīng)力剪切位移曲線如圖4～6所示。

1從樁巖接觸面剪應(yīng)力剪切位移曲線圖4～6可以看出，3種粗糙度條件下，在加載初期，接觸面剪應(yīng)力先隨剪切位移增加而增大，在達(dá)到峰值應(yīng)力后，隨著剪切位移增大，剪應(yīng)力反而降低，并最終趨于穩(wěn)定值。說明3種粗糙度條件下樁巖接觸面剪應(yīng)力剪切位移曲線均表現(xiàn)為應(yīng)變軟化特征。同時(shí)，隨著接觸面粗糙程度的增大，剪應(yīng)力峰值也隨之增大。

當(dāng)法向應(yīng)力較大時(shí)，剪應(yīng)力剪切位移曲線出現(xiàn)明顯拐點(diǎn)；當(dāng)法向應(yīng)力較小時(shí)，剪應(yīng)力剪切位移曲線達(dá)到峰值后出現(xiàn)下降段，最后趨于穩(wěn)定，但拐點(diǎn)不明顯。法向應(yīng)力越大，接觸面初始切向剛度越高，且在接觸面粗糙度相同條件下，法向應(yīng)力愈大，峰值應(yīng)力及殘余應(yīng)力也越大。

1.4.2接觸面直剪試驗(yàn)結(jié)果分析根據(jù)室內(nèi)樁巖接觸面剪切試驗(yàn)曲線的特征，可以看出剪應(yīng)力剪切位移曲線呈現(xiàn)一定的應(yīng)變軟化特征，大致表現(xiàn)為彈性階段、塑性階段、應(yīng)變軟化階段和殘余摩擦階段4個(gè)階段。

在加載初期，樁巖接觸面的剪應(yīng)力隨剪切位移的增大近似呈線性增加，抗剪強(qiáng)度主要由紅層泥巖與混凝土結(jié)構(gòu)粗糙接觸面上的咬合力承擔(dān)，此階段為彈性階段。

隨著水平荷載的持續(xù)增加，超過極限彈性剪應(yīng)力之后進(jìn)入塑性階段。在此階段，接觸面抗剪強(qiáng)度主要由兩方面構(gòu)成：一方面是由于紅層泥巖部分土體嵌入到混凝土粗糙表面的凸凹縫中，而發(fā)生的混凝土在剪切方向?qū)ν馏w產(chǎn)生的阻力；另一方面為剪切面上破碎的土顆粒之間的滑動(dòng)摩擦和紅層泥巖與混凝土結(jié)構(gòu)接觸面的滑動(dòng)摩擦。

隨著接觸面滑移增加，接觸面處紅層泥巖的原結(jié)構(gòu)被破壞，越來越多的破碎的土體顆粒嵌入到混凝土表面的凸凹縫中，從而使剪切過程中紅層泥巖與混凝土接觸面逐漸向平面剪切過度，表現(xiàn)為應(yīng)變軟化特征，此階段為軟化階段。

最后，混凝土的凸凹縫逐漸被破碎的土體填滿，而使得剪切過程中紅層泥巖與混凝土接觸面接近平面剪切，此時(shí)剪切強(qiáng)度主要由破壞后的紅層泥巖與混凝土接觸面存在的滑動(dòng)摩擦提供，因此，剪應(yīng)力隨著剪切位移的增加而最終趨于穩(wěn)定，此階段為殘余階段。

2接觸面本構(gòu)模型

2.1接觸面本構(gòu)模型

根據(jù)室內(nèi)樁巖接觸面剪切試驗(yàn)曲線呈現(xiàn)出：彈性階段、塑性階段、應(yīng)變軟化階段和殘余摩擦階段4個(gè)階段的特征，基于經(jīng)典雙曲線接觸面本構(gòu)模型[5]，建立以相對(duì)位移為變量的樁巖接觸面應(yīng)變軟化本構(gòu)模型，其關(guān)系曲線如圖7所示。

2.2參數(shù)的確定

令殘余強(qiáng)度τr=β τsu，則：當(dāng)β=1時(shí)，接觸面剪應(yīng)力剪切位移（τωs）關(guān)系為雙曲線接觸面模型；當(dāng)β<1時(shí)，接觸面剪應(yīng)力剪切位移（τωs）關(guān)系為應(yīng)變軟化接觸面本構(gòu)模型。由式（5）可求得參數(shù)a、b、c，見式（6）。b1c2=βτsu

綜上可知，確定a、b、c試驗(yàn)參數(shù)，需要用到極限剪切位移ωsu，極限剪應(yīng)力τsu，殘余強(qiáng)度τr，通過直剪試驗(yàn)結(jié)果，可以得到此3個(gè)參數(shù)，然后通過公式求出參數(shù)a、b、c。

3數(shù)值分析

有限差分軟件FLAC3D自帶的Interface接觸面單元，它是基于理想彈塑性Coulomb滑動(dòng)本構(gòu)關(guān)系的無厚度接觸單元，無法模擬樁巖接觸面在彈性階段表現(xiàn)出的非線性，以及在塑性階段表現(xiàn)出的應(yīng)變軟化等特征[19]。因此，利用fish語言對(duì)FLAC3D中自帶的接觸單元進(jìn)行二次開發(fā)，并應(yīng)用開發(fā)的模型對(duì)樁巖接觸面直剪試驗(yàn)進(jìn)行了數(shù)值模擬，分析剪應(yīng)力與剪切位移之間的關(guān)系，對(duì)樁巖接觸面應(yīng)變軟化本構(gòu)模型的準(zhǔn)確性進(jìn)行驗(yàn)證。

3.1FLAC3D接觸單元的二次開發(fā)

利用fish語言將樁巖接觸面應(yīng)變軟化本構(gòu)模型植入到FLAC3D中，由此建立新接觸面單元計(jì)算流程圖如圖8所示。

樁巖接觸面本構(gòu)模型計(jì)算循環(huán)步驟圖如下：

1）先在指定表面節(jié)點(diǎn)施加集中荷載，產(chǎn)生的表面位移轉(zhuǎn)化為節(jié)點(diǎn)速度，再通過節(jié)點(diǎn)速度計(jì)算新的應(yīng)變率，從而得到接觸面單元節(jié)點(diǎn)剪應(yīng)力，剪切位移；

2）利用接觸面剪應(yīng)力的取值判定接觸面是否發(fā)生破壞，即若，通過接觸面應(yīng)變軟化本構(gòu)模型計(jì)算新的接觸面單元節(jié)點(diǎn)剪應(yīng)力和剪切位移；若，則樁巖接觸面發(fā)生破壞，自動(dòng)賦予法向剛度較小值；

3）由上一步得到新的應(yīng)力和荷載，通過動(dòng)量方程計(jì)算新的節(jié)點(diǎn)速度和應(yīng)變率。

不斷重復(fù)上述3個(gè)步驟，在計(jì)算過程中，如果記錄的最大不平衡力接近零或者為常數(shù)，則停止計(jì)算。

3.2數(shù)值驗(yàn)證

參照室內(nèi)剪切試驗(yàn)條件，利用FLAC3D建立模型，選取粗糙度Ⅰ條件，法向應(yīng)力分別為50、150、250、350 kPa進(jìn)行數(shù)值模擬。所建模型如圖9所示。

對(duì)數(shù)值模擬結(jié)果進(jìn)行整理，作出不同法向應(yīng)力作用下，樁巖接觸面剪應(yīng)力剪切位移曲線如圖10所示。

從圖10中可以看出，在數(shù)值模擬中，加載初期，剪應(yīng)力增加到一定程度后產(chǎn)生剪切位移，之后，剪應(yīng)力隨著剪切位移呈近似線性增加，在達(dá)到應(yīng)力峰值后，隨著剪切位移增大，剪應(yīng)力反而降低，并最終趨于穩(wěn)定值。

通過分析比較室內(nèi)試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果，可以看出樁巖接觸面剪應(yīng)力剪切位移曲線趨勢(shì)大致相同，都經(jīng)歷了彈性階段、塑性階段、應(yīng)變軟化階段和殘余摩擦階段4個(gè)階段。表明剪應(yīng)力剪切位移曲線呈現(xiàn)出一定的應(yīng)變軟化特征，證明了該本構(gòu)能夠較好地模擬接觸面間的應(yīng)變軟化特性。

4結(jié)論

通過室內(nèi)樁巖接觸面剪切試驗(yàn)，結(jié)合理論分析和有限差分?jǐn)?shù)值模擬等方法，針對(duì)我國西南地區(qū)廣泛分布的紅層泥巖地層，從樁巖相互作用出發(fā)，對(duì)樁巖接觸面力學(xué)性質(zhì)進(jìn)行研究，得出以下結(jié)論：

1）通過室內(nèi)樁巖接觸面剪切試驗(yàn)，無論有無泥皮、何種粗糙度，接觸面剪應(yīng)力剪切位移曲線都出現(xiàn)拐點(diǎn)，在剪應(yīng)力達(dá)到峰值強(qiáng)度后開始降低并最終趨于穩(wěn)定值。3種混凝土表面粗糙度情況下接觸面的剪應(yīng)力剪切位移曲線都表現(xiàn)出應(yīng)變軟化的特征。

2）根據(jù)接觸面剪切試驗(yàn)成果，建立樁巖接觸面應(yīng)變軟化本構(gòu)模型，并推導(dǎo)其本構(gòu)方程。

3）利用fish語言對(duì)FLAC3D中自帶的理想彈塑性接觸單元進(jìn)行二次開發(fā)，并應(yīng)用開發(fā)的模型對(duì)樁巖接觸面剪切試驗(yàn)進(jìn)行了數(shù)值模擬，分析剪應(yīng)力與剪切位移之間的關(guān)系，證明了該本構(gòu)能夠較好地模擬接觸面間的應(yīng)變軟化特性。

參考文獻(xiàn)：

[1] 楊宗才，張俊云，周德培.紅層泥巖邊坡快速風(fēng)化特性研究[J].巖石力學(xué)與工程學(xué)報(bào)，2006，25（2）：275283.

YANGZ C，ZHANG J Y，ZHOU D P.Study on fast weathering characteristics of red bed mudstone slope [J].Chinese Journal of Rock Mechanics and Engineering， 2006，25（2）：275283.（in Chinese）

[2] 馮強(qiáng).四川紅層泥巖的分布及其路用性能研究[D].成都：西南交通大學(xué)，2011.

FENG Q. The distribution and road behavior of Sichuan red mudstone [D]. Chengdu： Southwest Jiaotong University，2011.（in Chinese）

[3] POTYONDY J G. Skin friction between various soils and construction materials [J]. Géotechnique， 1961， 11（4）：339353.

[4] MAYORAL J M， PESTANA J M， SEED R B. Modeling claypile interface during multidirectional loading [J].Computers and Geotechnics， 2016， 74： 163173.

[5] CLOUGH G W， DUNCAN J M. Finite element analyses of retaining wall behavior [J]. Asce Soil Mechanics & Foundation Division Journal， 1973， 97（12）：16571673.

[6] BRANDT J R T. Behavior of soilconcrete interfaces [D].Edmonton， Canada： University of Alberta， 1985.

[7] BOULON M. Basic features of soil structure interface behavior [J]. Computers and Geotechnics， 1989， 7（1）： 115131.

[8] 張冬霽，盧廷浩.一種土與結(jié)構(gòu)接觸面模型的建立及其應(yīng)用[J].巖土工程學(xué)報(bào)，1998（6）：6569.

ZHANG D J， LU T H. Establishment and application of an interface modelbetween soil and structure [J].Chinese Journal of Geotechnical Engineering，1998（6）：6569. （in Chinese）

[9] DAGUIAR S C， MODARESSIFARAHMANDRAZAVI A， SANTOS J A D， et al. Elastoplastic constitutive modelling of soilstructure interfaces under monotonic and cyclic loading [J]. Computers & Geotechnics， 2011， 38（4）：430447.

[10] 王偉，盧廷浩，周愛兆.土與混凝土接觸面反向單剪本構(gòu)模型[J].土木建筑與環(huán)境工程，2011，33（1）：5964.

WANG W，LU T H，ZHOU A Z.Constitutive model for negative simplesheared soilconcrete interface [J].Journal of Civil， Architectural & Environmental Engineering，2011，33（1）：5964. （in Chinese）

[11] 殷宗澤，朱泓，許國華. 土與結(jié)構(gòu)材料接觸面的變形及其數(shù)學(xué)模擬[J].巖土工程學(xué)報(bào)，1994，16（3）：1422.

YIN Z Z， ZHU H， XU G H.Numerical simulation of the deformation in the interface between soil and structural material [J].Chinese Journal of Geotechnical Engineering， 1994， 16（3）： 1422. （in Chinese）

[12] 張嘎，張建民.粗粒土與結(jié)構(gòu)接觸面的靜動(dòng)本構(gòu)規(guī)律[J].巖土工程學(xué)報(bào)，2005（5）：516520.

ZHANG G， ZHANG J M. Monotonic and cyclicconstitutive law of interface between structure and coarsegrained soil [J]. Chinese Journal of Geotechnical Engineering， 2005， 27（5）： 516520. （in Chinese）

[13] 張嘎，張建民. 粗粒土與結(jié)構(gòu)接觸面三維本構(gòu)關(guān)系及數(shù)值模型[J]. 巖土力學(xué)，2007，28（2）：288292.

ZHANG G， ZHANG J M. Threedimensional model of interface between structure and coarse grained soil [J].Rock and Soil Mechanics， 2007， 28（2）： 288292. （in Chinese）

[14] 馮大闊，侯文峻，張建民. 粗粒土與結(jié)構(gòu)接觸面三維力學(xué)特性的直剪試驗(yàn)研究[J]. 土木工程學(xué)報(bào)，2012（5）：169175.

FENG D K， HOU W J， ZHANG J M. Large scale direct shear test investigation of the 3D behavior of a gravel structure interfaces [J].China Civil Engineering Journal， 2012（5）： 169175. （in Chinese）

[15] 舒丹.內(nèi)河軟巖地基框架碼頭樁基承載特性研究[D].重慶：重慶交通大學(xué)，2015.

SHU D. Study on bearing characteristics of pile in soft rockfoundation with frame wharf in inland river [D].Chongqing：Chongqing Jiaotong University，2015. （in Chinese）

[16] WU H， SHU Y， ZHU J. Implementation and verification of interface constitutive model in FLAC 3D [J]. Water Science and Engineering， 2011， 4（3）： 305316.

[17] 金長宇，張春生，馮夏庭. 非線性彈性硬化接觸本構(gòu)模型在層間錯(cuò)動(dòng)帶中的研究[J]. 水力發(fā)電學(xué)報(bào)，2013（1）：206211，217.

JIN C Y， ZHANG C S， FENG X T. Study and application of nonlinear elasticharden contact constitutive model to disturbed belt [J].Journal of Hydroelectric Engineering，2013（1）：206211，217. （in Chinese）

[18] 鄭文華，劉德軍. 基于FLAC3D二次開發(fā)的不連續(xù)變形模擬方法探討[J]. 巖土工程學(xué)報(bào)，2013（Sup1）：197202.

ZHENG W H， LIU D J. Discussion of discontinuous deformation simulation method based on secondary development of FLAC3D [J].Chinese Journal of Geotechnical Engineering，2013（Sup1）：197202. （in Chinese）

[19] 劉瑩駿，張運(yùn)良. 一種薄層接觸單元的開發(fā)及FLAC3D實(shí)現(xiàn)[J]. 水利與建筑工程學(xué)報(bào)，2014（4）：18，49.

LIU Y J， ZHANG Y L. Development of a thinlayer interface element and its implementation in FLAC3D [J]. Journal of Water Resources and Architecture Engineering，2014（4）：18，49. （in Chinese）

[20] 劉蓓.紅粘土與混凝土結(jié)構(gòu)接觸面力學(xué)特性的試驗(yàn)研究[D].長沙：中南大學(xué)，2013.

LIU P. The test study on mechanical properties of interface between red clay and concrete structure [D]. Changsha：Central South University，2013. （in Chinese）

[21] 曾衛(wèi)，何淼. 混凝土表面粗糙度對(duì)梁耐熱性能影響的實(shí)驗(yàn)研究[J]. 四川建筑，2010（5）：160161，164.

ZENG W， HE M. Experiment study on thermal resistant property of beam effected by concrete surface roughness [J]. Sichuan Architecture，2010（5）：160161，164. （in Chinese）

（編輯王秀玲）