FLAC3D 中Cable 數(shù)量對(duì)錨桿拉拔模擬試驗(yàn)結(jié)果的影響

杜學(xué)領(lǐng)

（貴州理工學(xué)院 礦業(yè)工程學(xué)院，貴州 貴陽(yáng)550003）

錨桿是支護(hù)工程中的主要材料之一，近年來(lái)，D錨桿、恒阻大變形錨桿、中空注漿錨桿等新型錨桿的研發(fā)、錨桿監(jiān)測(cè)技術(shù)的進(jìn)步及支護(hù)理論的深入研究[1-5]，進(jìn)一步提升了錨桿種類的豐富性和支護(hù)設(shè)計(jì)的科學(xué)性。FLAC3D是巖土、礦業(yè)領(lǐng)域被廣泛應(yīng)用的數(shù)值模擬軟件之一，斷層影響、沖擊地壓評(píng)價(jià)及支護(hù)系統(tǒng)特性等問題均可在FLAC3D中進(jìn)行研究[6-10]。但在FLAC3D中，依托Cable 建立錨桿的模擬效果與其數(shù)量有關(guān)。不同的Cable 數(shù)量，可能影響到模擬的準(zhǔn)確性。盡管Itasca 公司提出了2 條建議[11]：“1. Try to provide approximately one node in each FLAC3Dzone.The reasoning here is that since the zones are constant-stress regions, it is not necessary to have more than one interaction point within a zone; 2. Try to provide two to three cable element within the development length of the cable. The “development length of the cable is determined by dividing the specified yield strength, Ft, by the grout cohesive strength, cg. By following this procedure, failure by“pull-out”can occur if such conditions arise. If the cable elements are too long, then only the yield failure mode of each element is possible”，上述建議同樣適用于使用Pile 結(jié)構(gòu)單元建立巖石錨桿的情況。但在已有文獻(xiàn)中，一些研究并沒有應(yīng)用上述原則或者未提及數(shù)值模擬中Cable 的數(shù)量問題，例如李為騰等在拉拔模擬和支護(hù)模擬中分別設(shè)定單個(gè)Cable 長(zhǎng)度為0.05、0.1 m，但其并未指明Cable 密度與圍巖網(wǎng)格尺寸的關(guān)系[12]；劉少虹對(duì)單個(gè)Cable 長(zhǎng)度的設(shè)定為0.11 m，其Cable與圍巖為非均勻接觸，非均勻條件下的錨桿受力是否存在尺寸效應(yīng)仍需要進(jìn)一步評(píng)價(jià)[13]；黃慶享等對(duì)軟巖大變形巷道支護(hù)方案的模擬中，并沒有給出數(shù)值模擬中錨桿所采用的物理參數(shù)和設(shè)置方案[14]；類似的未給出數(shù)值模擬中錨桿設(shè)置詳情的研究還大量存在[15-17]，以至于降低了數(shù)值模擬的可信度和可重復(fù)性。尺寸效應(yīng)廣泛存在于PFC、UDEC 等數(shù)值模擬及物理實(shí)驗(yàn)中，并已有大量尺寸效應(yīng)報(bào)道[18-21]。但FLAC3D中Cable 數(shù)量與圍巖網(wǎng)格的匹配性問題，尚未引起足夠重視。對(duì)于煤礦巷道支護(hù)工程而言，頂板、幫部的錨桿長(zhǎng)度一般為1.6～2.4 m 左右，與邊長(zhǎng)為30～50 m 的數(shù)值模擬研究模型相比，Cable 數(shù)量是否影響錨桿性能還需要進(jìn)一步評(píng)價(jià)。在固定網(wǎng)格尺寸、錨桿長(zhǎng)度的情況下，通過在拉拔試驗(yàn)中改變Cable 的數(shù)量，評(píng)價(jià)FLAC3D中不同Cable數(shù)量對(duì)模擬結(jié)果的影響。

1 研究方案

借鑒已有研究成果，在FLAC3D中利用Cable 結(jié)構(gòu)單元建立錨桿模型，模型尺寸及測(cè)點(diǎn)布置如圖1。

圖1 模型尺寸及測(cè)點(diǎn)布置示意圖Fig.1 Schematic diagram of model size and position of the measurement points

錨桿長(zhǎng)度為0.4 m，起點(diǎn)及拉拔端均位于坐標(biāo)原點(diǎn)O 處。圍巖尺寸為x×y×z=寬×長(zhǎng)×高=0.3 m×0.6 m×0.3 m，圍巖網(wǎng)格密度為0.01 m 1 個(gè)網(wǎng)格。圍巖采用應(yīng)變軟化模型，其參數(shù)與文獻(xiàn)[22]中表1、表2 中所采用的煤的參數(shù)相同；錨桿采用全長(zhǎng)錨固，Cable 的參數(shù)與文獻(xiàn)[22]中表3、表5 中的參數(shù)相同。文獻(xiàn)[22]未提及的grout-perimeter、grout-stiffness，本文分別設(shè)定為0.094 2 m、28 GPa。與文獻(xiàn)[22]有所不同的是，本文僅研究Cable 數(shù)量對(duì)模擬結(jié)果的影響，因此不對(duì)錨桿施加預(yù)緊力，且由于錨桿長(zhǎng)度、模型尺寸及網(wǎng)格劃分與文獻(xiàn)[22]有所不同，本文采用的拉拔速度為10-6m/step，所有模型均運(yùn)行10 000步以便于對(duì)比模擬結(jié)果。拉拔結(jié)束時(shí)，拉拔端的結(jié)點(diǎn)（node）位移為10 mm，根據(jù)物理拉拔試驗(yàn)中，常規(guī)的錨桿在此位移下已達(dá)到峰值載荷并已進(jìn)入峰后階段[23-24]。模型的邊界條件與FLAC3D手冊(cè)中的案例“Simulation of Pull-Tests for Fully Bonded Rock Reinforcement”使用Cable 建立錨桿的邊界條件相同，且不考慮密度及重力的影響[11]。沿錨桿軸向每隔0.05 m 布置1 個(gè)測(cè)點(diǎn)監(jiān)測(cè)拉拔過程中的錨桿軸力演化，并分別編號(hào)為C0、C1、…、C8。在z=0 平面，在x=y=0.01 m 處開始沿y 軸間隔0.2 m 布置3 個(gè)測(cè)點(diǎn)監(jiān)測(cè)拉拔過程中圍巖中的應(yīng)力響應(yīng)，此3 個(gè)測(cè)點(diǎn)可反映錨桿兩端及中部所引用的圍巖應(yīng)力變化。

根據(jù)文獻(xiàn)[22]的參數(shù)及文獻(xiàn)[11]的建議，本文的研究條件下錨桿的development length 為0.133 m，該長(zhǎng)度上應(yīng)有2～3 個(gè)Cable，即：約0.04 m 1 個(gè)Cable。但文獻(xiàn)[11]同時(shí)提出，應(yīng)對(duì)1 個(gè)zone 提供接近1個(gè)而無(wú)需過多的結(jié)點(diǎn)，圍巖網(wǎng)格密度為0.01 m 1 個(gè)網(wǎng)格，按此計(jì)算則約0.01 m 1 個(gè)Cable。二者存在一定矛盾性，這種矛盾也會(huì)因研究問題的差異而出現(xiàn)在不同的數(shù)值模擬中。為了研究Cable 數(shù)量對(duì)模擬結(jié)果，本文以0.04 m 1 個(gè)Cable 為基準(zhǔn)，并標(biāo)記為“標(biāo)準(zhǔn)組”；分別以0.02、0.01、0.005 m 1 個(gè)Cable 為基準(zhǔn)建立加密的Cable 研究方案；分別以0.05、0.08、0.1、0.2 m 1 個(gè)Cable 為基準(zhǔn)建立相對(duì)疏松的Cable研究方案，此時(shí)部分軸力測(cè)點(diǎn)的數(shù)據(jù)完全相同，在后續(xù)分析中僅選取部分測(cè)點(diǎn)用于說(shuō)明問題。

2 數(shù)值模擬結(jié)果及分析

2.1 拉拔端的軸力-時(shí)步演化特征

對(duì)于拉拔試驗(yàn)而言，常規(guī)錨桿的軸力沿拉拔反方向快速衰減，最大軸力出現(xiàn)在拉拔端附近[1]，因此，拉拔端的軸力在一定程度上能反映錨桿在拉拔過程中的特性。不同Cable 數(shù)量時(shí)錨桿拉拔端軸力-時(shí)步演化過程如圖2，由圖2 可知，盡管圍巖屬性、網(wǎng)格密度、邊界條件、錨桿長(zhǎng)度等條件均相同，但對(duì)相同長(zhǎng)度錨桿劃分為不同Cable 數(shù)量時(shí)，模擬結(jié)果有明顯不同。Cable 數(shù)量與單元網(wǎng)格密度具有重要相關(guān)性，當(dāng)單一Cable 長(zhǎng)度與網(wǎng)格長(zhǎng)度劃分一致或比網(wǎng)格更短時(shí)，錨桿能獲得較大的拉拔軸力，如0.01 m 1 個(gè)Cable 時(shí)（單一Cable 長(zhǎng)度與網(wǎng)格長(zhǎng)度劃分一致）錨桿最大軸力可達(dá)到114 kN，0.005 m 1 個(gè)Cable 時(shí)（單一Cable 長(zhǎng)度是網(wǎng)格長(zhǎng)度的2 倍）錨桿最大軸力可達(dá)到錨桿的抗拉強(qiáng)度266 kN；而當(dāng)單一Cable 長(zhǎng)度低于網(wǎng)格長(zhǎng)度時(shí)，錨桿拉拔端最大軸力要低于50 kN，該數(shù)值遠(yuǎn)小于錨桿的抗拉強(qiáng)度，說(shuō)明此時(shí)錨桿并沒有發(fā)揮出高承載的性能。由于不同的錨固條件下所獲得的錨固效果是不同的，可以根據(jù)現(xiàn)場(chǎng)的拉拔試驗(yàn)確定錨桿的極限拉拔力，然后在數(shù)值模擬中設(shè)定單一Cable 長(zhǎng)度為網(wǎng)格長(zhǎng)度的0.5～1 倍之間，尋找與現(xiàn)場(chǎng)試驗(yàn)相符合的模擬方案。

圖2 錨桿拉拔端軸力-時(shí)步演化過程圖Fig.2 Evolution of axial force at the pull-out end and step

2.2 拉拔引起的圍巖應(yīng)力演化

不同模擬方案中圍巖不同位置沿拉拔方向的應(yīng)力演化過程如圖3。

圖3 圍巖不同位置沿拉拔方向的應(yīng)力演化過程Fig.3 Stress evolution of different positions of surrounding rocks along the pulling direction

從圖3 可知，因?qū)ο嗤L(zhǎng)度的錨桿設(shè)置不同的Cable 數(shù)量，導(dǎo)致圍巖中的應(yīng)力演化也會(huì)出現(xiàn)較大差異。當(dāng)1 個(gè)Cable 的長(zhǎng)度顯著大于網(wǎng)格長(zhǎng)度時(shí)（1 個(gè)Cable 長(zhǎng)度在0.05～0.5 m 時(shí)），造成拉拔端附近的圍巖應(yīng)力較低（圖3（a）），而錨桿尾部的應(yīng)力則相對(duì)較大（圖3（c）），這種情況與物理實(shí)驗(yàn)是不符合的。在物理實(shí)驗(yàn)中，錨桿尾部的軸力較小，因此其引起的附加應(yīng)力場(chǎng)也相對(duì)較小，從這個(gè)意義而言，對(duì)于相同長(zhǎng)度的錨桿，當(dāng)Cable 數(shù)量較少時(shí)，無(wú)法正確反映出錨桿與圍巖的相互作用。在單純拉拔條件下，錨桿中部附近的圍巖應(yīng)力變化應(yīng)體現(xiàn)出從拉拔端到錨桿尾部的漸變特性，其數(shù)量級(jí)一般應(yīng)低于拉拔端附近的應(yīng)力量級(jí)，由此觀之，圖3（b）中1 個(gè)Cable 長(zhǎng)度為0.04 時(shí)的情況也與物理實(shí)驗(yàn)不符。

對(duì)上述不合理方案進(jìn)一步檢視發(fā)現(xiàn)，當(dāng)Cable的數(shù)量較少時(shí)，圍巖中的應(yīng)力呈現(xiàn)出分散性分布的特點(diǎn)，標(biāo)準(zhǔn)組方案中圍巖應(yīng)力分布圖如圖4，Cable數(shù)量少于標(biāo)準(zhǔn)組時(shí)，其圍巖中的應(yīng)力分布均與此類似。由圖4 可知，當(dāng)Cable 數(shù)量較少時(shí)，圍巖中的高應(yīng)力以靠近結(jié)點(diǎn)所在位置為中心，而與拉拔位置關(guān)系不大。這顯然并不符合常規(guī)全長(zhǎng)錨固錨桿的拉拔特性。這不僅說(shuō)明Cable 數(shù)量較少時(shí)所表征的工程問題極有可能是錯(cuò)誤的，還說(shuō)明按照FLAC3D幫助手冊(cè)中根據(jù)development length 設(shè)置Cable 數(shù)量的建議也是不可靠的。

圖4 標(biāo)準(zhǔn)組方案圍巖中分散性應(yīng)力分布Fig.4 Dispersion stress distribution in surrounding rock of the standard group

3 數(shù)值模擬討論

3.1 單位長(zhǎng)度錨固劑剛度問題

前述研究設(shè)定單位長(zhǎng)度錨固劑剛度grout-stiffness（以下簡(jiǎn)稱“錨固劑剛度”）為28 GPa，由于文獻(xiàn)[22]并未明確給出錨固劑剛度的數(shù)值，本文分別設(shè)定錨固劑剛度為5、10、15、20、35、40 GPa，以單一Cable 長(zhǎng)度與網(wǎng)格尺寸相一致的模型為基礎(chǔ)，進(jìn)一步觀察模擬結(jié)果的差異。不同錨固劑剛度時(shí)的軸力-時(shí)步演化圖如圖5。

圖5 不同錨固劑剛度時(shí)的軸力-時(shí)步演化圖Fig.5 Axial force-step evolution with different grout-stiffness

由圖5 可知，在本文的研究條件下，錨固劑剛度在5～40 GPa 之間變化時(shí)，雖然錨桿軸力隨著拉拔時(shí)步的演化不盡相同，但從軸力的峰值而言錨固劑剛度變化的影響并不顯著，峰值軸力約在110～130 kN 之間變化，與Cable 數(shù)量的改變相比，該變化幅值對(duì)錨桿峰值軸力的影響是相對(duì)較小的。這也說(shuō)明，盡管本文指定了錨固劑剛度為28 GPa，但其對(duì)峰值軸力的分析是相對(duì)可靠的。但錨固劑剛度改變后，會(huì)改變軸力的峰后變化特征。在本文的研究條件下，錨固劑剛度低于35 GPa 時(shí)，錨桿峰后軸力波動(dòng)趨穩(wěn)，而錨固劑剛度不小于35 GPa 時(shí)，錨桿峰后軸力波動(dòng)下降。

3.2 圍巖軟化特性對(duì)模擬結(jié)果的影響

本文的圍巖軟化特性參照文獻(xiàn)[22]設(shè)定，以前述參數(shù)為基礎(chǔ)，對(duì)標(biāo)準(zhǔn)試樣進(jìn)行單軸壓縮的數(shù)值模擬測(cè)試，單軸壓縮模擬結(jié)果如圖6。

圖6 單軸壓縮數(shù)值模擬結(jié)果Fig.6 Numerical simulation results of uniaxial compression test

由圖6 可知，數(shù)值模擬中所采用的煤的單軸抗壓強(qiáng)度約為5.37 MPa，該數(shù)值相對(duì)于巖石而言屬于軟弱巖體。進(jìn)一步，在模擬中受應(yīng)變軟化參數(shù)設(shè)置影響，煤體破壞后的強(qiáng)度依然保持在2.3 MPa 左右，峰后強(qiáng)度約為峰值強(qiáng)度的42.83%。在真實(shí)的物理實(shí)驗(yàn)中，煤巖體一般表現(xiàn)為脆性或塑性破壞，但破壞后的強(qiáng)度不會(huì)維持在較高的強(qiáng)度水平，而是隨著位移的增加而逐漸衰減。從圖6 還可以看出，圍巖達(dá)到峰值載荷的位移量非常小，而在峰后則維持了較大的位移量。對(duì)于在FLAC3D中的靜力學(xué)分析而言，各種物理量的疊加計(jì)算依賴于網(wǎng)格及其節(jié)點(diǎn)，由于軟化后的網(wǎng)格依然能夠承載、拉拔過程中的破壞具有非均勻演化的特點(diǎn)（體現(xiàn)在拉伸和剪切破壞有可能在不同時(shí)步發(fā)生于同一網(wǎng)格）等原因，就導(dǎo)致了拉拔模擬的軸力-時(shí)步曲線出現(xiàn)較為明顯的波動(dòng)性，這種波動(dòng)性的波動(dòng)幅值、波動(dòng)趨勢(shì)與物理實(shí)驗(yàn)是無(wú)法做到完全一致的。因此，盡管在數(shù)值模擬中可以通過參數(shù)校核來(lái)獲得與物理實(shí)驗(yàn)較為吻合的結(jié)果，但實(shí)際上參數(shù)校核也只是在數(shù)值模擬中實(shí)現(xiàn)部分物理實(shí)驗(yàn)的重現(xiàn)，與天然煤巖體的在實(shí)驗(yàn)前的未知離散性表征還有一定差距?？梢?，在當(dāng)前的技術(shù)及理論研究條件下，數(shù)值模擬完全替代物理實(shí)驗(yàn)而獲得服務(wù)于現(xiàn)實(shí)的精確解是非常困難的。

3.3 數(shù)值模擬中的錨桿特性問題

文獻(xiàn)[22]研究了利用Fish 語(yǔ)言實(shí)現(xiàn)FLAC3D中可破斷錨桿的模擬問題，但從圖2、圖5（a）可以看出，Cable 數(shù)量、錨固劑剛度等均會(huì)對(duì)錨桿特性造成重要影響。特別是在拉拔過程中錨桿未達(dá)到極限抗拉強(qiáng)度的情況下，錨桿達(dá)到峰值軸力后軸力會(huì)出現(xiàn)衰減，此后即便發(fā)生較大的位移量錨桿也不會(huì)達(dá)到極限抗拉強(qiáng)度。按照文獻(xiàn)[22]基于錨桿自由段總伸長(zhǎng)量來(lái)進(jìn)行錨桿破斷判定的方法，如果數(shù)值模擬中錨桿與圍巖的交互作用導(dǎo)致錨桿的響應(yīng)本身就不合理，如圖5（a），此時(shí)錨桿軸力演化本身不合理，其位移也就不具備參考價(jià)值，此時(shí)應(yīng)用文獻(xiàn)[22]的判定方法，可能也是存疑的。此外，對(duì)于錨固系統(tǒng)破壞時(shí)錨桿被直接拔出的情況，如軟巖、礦井水影響等條件下就有可能因錨固不牢靠而發(fā)生錨桿的非破斷脫錨，如何在FLAC3D中更加準(zhǔn)確的反映出這些現(xiàn)象，依然需要深入研究。

4 結(jié) 論

1）對(duì)相同長(zhǎng)度錨桿劃分為不同Cable 數(shù)量時(shí)，模擬結(jié)果有明顯不同。當(dāng)單一Cable 長(zhǎng)度與網(wǎng)格長(zhǎng)度劃分一致或比網(wǎng)格更短時(shí)，錨桿能獲得較大的拉拔軸力，可以根據(jù)現(xiàn)場(chǎng)的拉拔試驗(yàn)確定錨桿的極限拉拔力，然后在數(shù)值模擬中設(shè)定單一Cable 長(zhǎng)度為網(wǎng)格長(zhǎng)度的0.5～1 倍之間，尋找與現(xiàn)場(chǎng)試驗(yàn)相符合的模擬方案。按照FLAC3D幫助手冊(cè)中根據(jù)development length 設(shè)置Cable 數(shù)量的建議并不完全可靠。

2）從錨桿軸力的峰值而言，錨固劑剛度的變化影響并不顯著，但會(huì)改變軸力的峰后變化特征。錨固劑剛度較大時(shí)，錨桿峰后軸力波動(dòng)下降。

3）當(dāng)前的技術(shù)及理論研究條件下，數(shù)值模擬完全替代物理實(shí)驗(yàn)而獲得服務(wù)于現(xiàn)實(shí)的精確解是非常困難的。對(duì)于僅有模擬結(jié)果而缺乏模擬細(xì)節(jié)的結(jié)論，應(yīng)審慎采信。