蘇永華，廖君橙，黃騰

摘? ?要：為了研究城區(qū)塌陷土洞的發(fā)育致塌機(jī)制，通過(guò)建立三維模型試驗(yàn)研究土洞的發(fā)育破壞過(guò)程，分析了試驗(yàn)過(guò)程中地表變形數(shù)據(jù)以及土洞的最終破壞形狀，同時(shí)通過(guò)土壓力量測(cè)試驗(yàn)分析土洞發(fā)展過(guò)程中覆土應(yīng)力的變化. 以此試驗(yàn)為基礎(chǔ)建立了計(jì)算土洞極限半徑的有限元模型，通過(guò)正交試驗(yàn)分析了彈性模量、泊松比、黏聚力、內(nèi)摩擦角對(duì)極限半徑的影響. 結(jié)果表明：對(duì)于發(fā)育的土層空洞，破壞區(qū)內(nèi)土壓力顯著降低;土洞發(fā)生破壞時(shí)，洞頂沉降速率發(fā)生突變;土洞內(nèi)部破壞的高度與土洞直徑間的規(guī)律與普氏理論提出的幾何規(guī)律吻合良好;對(duì)于埋深不同的土洞，其破壞形式主要分為壇形和直筒形塌陷.

關(guān)鍵詞：土洞破壞;室內(nèi)模型;有限元計(jì)算;土壓力;沉降

中圖分類(lèi)號(hào)：U416.1? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

Study on Development and Failure Characteristics

of Soil Caves in Urban Surface Subsidence

SU Yonghua？覮，LIAO Juncheng，HUANG Teng

（College of Civil Engineering，Hunan University，Changsha 410082，China）

Abstract：In order to study the development and collapse mechanism of the collapsed soil cave in the urban area， the development and failure process of the soil cave were studied through the establishment of a three-dimensional model test. The surface deformation data during the test and the final failure shape of the soil cave were analyzed. Changes in the stress of the covering soil during the development of the cave was also analyzed through the measurement tests of earth pressure. Based on this experiment， a finite element model for calculating the limit radius of the soil cave was established， and the influence of elastic modulus， Poisson's ratio， cohesion and internal friction angle on the limit radius was analyzed through orthogonal experiments. The results show that： for the developed soil cavity， the earth pressure in the failure area is significantly reduced; when the soil cave is damaged， the settlement velocity of the cave top changes suddenly; the law between the height of the damage inside the soil cave and the diameter of the soil cave is consistent with that proposed by Platts theory. The geometric rules are in good agreement; for soil caves with different burial depths， the damage modes are mainly divided into altar-shaped and straight-tube-shaped collapses.

Key words：soil cave failure;indoor model;finite element calculation;earth pressure;settlement

城區(qū)地表塌陷是指在城市開(kāi)發(fā)過(guò)程中，建設(shè)項(xiàng)目對(duì)地層不斷造成擾動(dòng)，致使土層中形成內(nèi)部空洞并逐步發(fā)展，直至地層土體達(dá)到承載能力極限狀態(tài)并最終發(fā)生淺層地表失穩(wěn)的現(xiàn)象[1]. 該類(lèi)土洞發(fā)育隱蔽，事故后維護(hù)速度快，相關(guān)事故現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)提取難度大，因此對(duì)于地表塌陷事故的可靠測(cè)量數(shù)據(jù)較少. 根據(jù)相關(guān)統(tǒng)計(jì)，北京市2012年發(fā)生道路塌陷事故300余次;深圳市2017年報(bào)道的該類(lèi)事故高達(dá)238次，我國(guó)其他各大城市，城區(qū)地表塌陷事故多有報(bào)道. 其塌陷形式、塌陷規(guī)模與產(chǎn)生的損失各不相同.

相關(guān)理論研究表明，當(dāng)?shù)貙油炼窗l(fā)育時(shí)，土洞上方發(fā)生應(yīng)力重新分布[2-3]，最終土洞發(fā)生破壞的形式為上覆土剪切破壞[4]. 同時(shí)，文獻(xiàn)[4]指出，當(dāng)土洞產(chǎn)生于黏性土中并發(fā)生塌陷時(shí)，塌陷形式為沿著空洞外圍的豎向筒狀塌陷，并利用極限平衡法構(gòu)建了圓柱形土洞覆土安全系數(shù)計(jì)算方法. 對(duì)于黃土路基下伏土洞，研究將黃土內(nèi)土洞發(fā)育簡(jiǎn)化為土洞頂面的抬升過(guò)程[5]，基于普氏理論探討了覆土臨界厚度.

該類(lèi)事故的試驗(yàn)研究可歸納為以下兩大類(lèi)：假設(shè)土洞為理想球體，基于三維模型的試驗(yàn)[6-7];設(shè)置土洞為二維平面模型，探究其發(fā)展破壞機(jī)理[8-9];將土洞假設(shè)為三維豎直筒形，研究覆土在該條件下的極限覆土厚度[10].

相關(guān)理論與試驗(yàn)研究結(jié)果表明，工程中發(fā)生的地表塌陷事故為三維土洞通過(guò)逐步塌落發(fā)展引發(fā)的事故，二維試驗(yàn)對(duì)于其發(fā)展過(guò)程的研究?jī)r(jià)值有限;三維模型試驗(yàn)使用黏土為填料[11]，假設(shè)土洞為標(biāo)準(zhǔn)球形且均勻發(fā)育. 對(duì)于城區(qū)地表塌陷，事故地層為較為松散的填土地層，對(duì)于類(lèi)似土洞發(fā)育過(guò)程，在工程實(shí)際中，影響因素較多，實(shí)際工程中難以預(yù)測(cè).

鑒于以上的土洞發(fā)育致塌特點(diǎn)，基于目前的研究現(xiàn)狀，模型試驗(yàn)研究計(jì)劃采用以下路線(xiàn)：選取砂土混合物為試驗(yàn)材料，通過(guò)改變填土高度研究不同埋深下土洞的破壞特征，綜合考慮土洞在外力影響下的發(fā)展，試驗(yàn)中逐漸擴(kuò)大土洞的半徑，并通過(guò)監(jiān)測(cè)地面沉降量分析其發(fā)育過(guò)程中的覆土變形特征.

1? ?試驗(yàn)步驟

1.1? ?試驗(yàn)?zāi)Ｐ拖湓O(shè)計(jì)

試驗(yàn)?zāi)Ｐ拖錇殚L(zhǎng)寬高均為1.2 m的土箱，使用螺栓連接角鋼為其外框架，為方便沉降觀(guān)測(cè)，模型箱體側(cè)面均采用鋼化玻璃材料. 在模型箱底部設(shè)置直徑為18 cm的圓形開(kāi)孔，作為模擬土洞的初始直徑，在底部布置兩道型鋼梁控制模型箱底部變形.

1.2? ?測(cè)量點(diǎn)位布置

針對(duì)土洞內(nèi)部坍塌發(fā)育過(guò)程中地面沉降情況以及土洞擴(kuò)展引發(fā)的地面塌陷形式. 試驗(yàn)沉降監(jiān)測(cè)點(diǎn)與土壓力測(cè)點(diǎn)布置如圖1和圖2所示. 試驗(yàn)過(guò)程中，隨著土洞半徑擴(kuò)展，洞頂覆土發(fā)生沉降，最終產(chǎn)生突塌，以覆土沉降為觀(guān)測(cè)對(duì)象，采用攝影測(cè)量技術(shù)，在洞頂對(duì)稱(chēng)設(shè)置間隔為10 cm的標(biāo)記物進(jìn)行地表沉降觀(guān)測(cè)，測(cè)點(diǎn)標(biāo)記為C1～C5，該測(cè)試系統(tǒng)通過(guò)不同發(fā)育狀態(tài)下測(cè)點(diǎn)的位置確定該狀態(tài)下模型箱內(nèi)表土的相對(duì)沉降量.

試驗(yàn)土壓力測(cè)量系統(tǒng)包括5個(gè)量程為30 kPa的CS-TYL050土壓力盒，以土洞垂直中心線(xiàn)為軸心，在同一水平面上對(duì)稱(chēng)布置于土洞上方，編號(hào)為P1～P5，直徑28.0 mm，厚度6.5 mm. 為保證土壓力盒不對(duì)土體應(yīng)力分布產(chǎn)生影響，參考已有研究中同樣大小土壓力傳感器的設(shè)置方法[12]，將土壓力傳感器間距設(shè)置為150 mm.

1.3? ?填土選擇與基本物理性質(zhì)測(cè)定

目前對(duì)土洞坍塌的試驗(yàn)研究主要采用的材料有以下兩種：其一，以干砂作為填料，該方案試驗(yàn)效果顯著，體現(xiàn)了松散介質(zhì)中土體的沉降破壞過(guò)程，但是其忽略了實(shí)際工程中的黏聚力參數(shù);其二，以黏土為試驗(yàn)材料，考慮了實(shí)際工程材料的黏聚力，但是其試驗(yàn)現(xiàn)象不夠明顯. 為了兼顧工程實(shí)際中土體的材料參數(shù)特性，試驗(yàn)填料采用與路面路基材料成分相近的工程砂土混合材料作為模型填料，使用Geotest室內(nèi)大型剪切儀測(cè)試材料，結(jié)果如圖3所示.

試驗(yàn)分別在20 kPa、30 kPa和50 kPa的側(cè)壓力下進(jìn)行剪切，直至剪應(yīng)力無(wú)法增加. 取曲線(xiàn)水平段平均值為該壓力下材料抗剪強(qiáng)度. 直線(xiàn)擬合其黏聚力為10 kPa，摩擦角為25°.

土樣級(jí)配曲線(xiàn)如圖4所示，其不均勻系數(shù)Cu = 6.05，曲率系數(shù)Cc = 1.05，試驗(yàn)材料為良好級(jí)配砂質(zhì)填土.

1.4? ?不同埋深土洞發(fā)育模擬

針對(duì)試驗(yàn)砂土進(jìn)行3組對(duì)比試驗(yàn)，每組試驗(yàn)中均以試驗(yàn)箱底部9 cm半徑圓形開(kāi)口作為初始土洞直徑，通過(guò)逐步擴(kuò)大土洞半徑模擬土洞發(fā)育過(guò)程，設(shè)定每次將半徑擴(kuò)大1 cm，測(cè)試其穩(wěn)定后的土壓力與表面沉降. 試驗(yàn)安排見(jiàn)表1.

試驗(yàn)在紅磚砌筑的下部懸空的操作臺(tái)上進(jìn)行，具體操作為：1）將模型箱在操作臺(tái)上水平放置;2）封閉試驗(yàn)箱下部開(kāi)孔;3）分層填土，試驗(yàn)材料模型采用分層壓實(shí)的方法，每10 cm壓實(shí)一次;4）擊實(shí)后靜壓24 h進(jìn)行模型土填筑;5）布置試驗(yàn)測(cè)點(diǎn)標(biāo)記物;6）打開(kāi)試驗(yàn)箱下部開(kāi)口，形成初始塌落面;7）使用每一步對(duì)應(yīng)長(zhǎng)度刀具分步擴(kuò)大土洞半徑，待土洞發(fā)育穩(wěn)定后測(cè)量塌落尺寸并使用攝影測(cè)量技術(shù)監(jiān)測(cè)沉降點(diǎn)位的位移直至土洞發(fā)生貫穿塌孔.

模型土填筑過(guò)程中要保證試驗(yàn)土的均勻性，具體操作如圖5所示. 為了保證試驗(yàn)土體的均勻，在試驗(yàn)箱中埋設(shè)取土盒，對(duì)試驗(yàn)土體進(jìn)行取樣，試驗(yàn)結(jié)束后取出埋設(shè)的取樣盒并且測(cè)量試驗(yàn)取樣參數(shù)，由此測(cè)量試驗(yàn)用土的密度參數(shù)γ與含水率w，試驗(yàn)土體物理力學(xué)參數(shù)見(jiàn)表2.

2? ?不同埋深土洞發(fā)育破壞

2.1? ?土洞發(fā)育致塌過(guò)程

對(duì)于土洞發(fā)育的研究表明，當(dāng)土體內(nèi)部產(chǎn)生懸空時(shí)，上覆土層逐步發(fā)生坍塌，其坍塌產(chǎn)生類(lèi)似拱形的界面[13-14]，試驗(yàn)過(guò)程顯示，在該類(lèi)坍塌過(guò)程中，最終坍塌面為橢球面. 因此，將空洞坍塌形狀參數(shù)簡(jiǎn)化為開(kāi)孔半徑R與其對(duì)應(yīng)的坍塌高度h來(lái)描述.

試驗(yàn)過(guò)程中每一步發(fā)育土洞底部的半徑R不斷擴(kuò)大，同時(shí)對(duì)應(yīng)的土洞高度h也逐漸增大，最終地面塌陷，產(chǎn)生半徑為r的塌陷窗口. 針對(duì)以上三組不同填土高度試驗(yàn)，土洞發(fā)生驟發(fā)塌陷時(shí)的土洞尺寸各不相同：埋深為400 mm時(shí)覆土塌穿時(shí)土洞的開(kāi)口半徑為140 mm;埋深為500 mm時(shí)，覆土塌穿時(shí)土洞開(kāi)口半徑為190 mm;埋深為600 mm時(shí)，土洞的塌孔半徑為160 mm，其塌陷高度隨土洞開(kāi)口半徑的變化見(jiàn)表3.

對(duì)于不同填土高度的模型，試驗(yàn)結(jié)果表明在相同填土材料的條件下，土洞處于極限狀態(tài)時(shí)，對(duì)于埋深較淺土洞，其直徑隨埋深增加而增加，而埋深較深的土洞由于上覆荷載的增加，其直徑有所減小. 為了分析土洞發(fā)育過(guò)程中半徑R與塌落高度h的關(guān)系，塌落高度與開(kāi)孔直徑之比即為土洞發(fā)育過(guò)程中的無(wú)量綱化高跨比n，高跨比隨半徑R變化如圖6所示.

對(duì)于覆土下存在土洞的情況，普氏理論針對(duì)巖溶土洞的發(fā)育過(guò)程將土洞直徑與覆土塌陷后自然平衡拱高間建立的關(guān)系為：

hmax = ■? ? ? ? ?（1）

式中：fk為堅(jiān)固系數(shù)，《巖土工程手冊(cè)》建議在松散土中取該值為0.5. 實(shí)際工程應(yīng)用中將式（1）乘以經(jīng)驗(yàn)系數(shù)0.828進(jìn)行計(jì)算，即：

hmax = 0.828■? ? ? ? ?（2）

將試驗(yàn)結(jié)果與理論計(jì)算值相對(duì)照，已有理論研究主要針對(duì)埋藏較深巖溶或隧道上方，其覆土強(qiáng)度指標(biāo)相比于淺層的松散土體更高，引發(fā)的塌陷范圍相比試驗(yàn)結(jié)果更小. 對(duì)于淺層空洞而言，采用普氏理論的計(jì)算結(jié)果更符合實(shí)際工程情況.

2.2? ?土洞塌陷類(lèi)型分析

對(duì)以上三組試驗(yàn)進(jìn)一步分析其塌陷模式，不同埋深情況下，土洞的極限直徑D、地面塌陷直徑d及其比值k如表4所示.

試驗(yàn)結(jié)果表明，相同材料條件下，當(dāng)不同埋深的土洞發(fā)生破壞時(shí)，其覆土塌陷模式可以分為壇形塌陷與直筒形塌陷.

2.3? ?覆土沉降模式

采用攝影測(cè)量的方法記錄N1和N2每一次土洞發(fā)育完成后洞頂?shù)某两盗縎與其半徑方向坐標(biāo)X，繪制為如圖7、圖8所示沉降分布圖. 同時(shí)，將C3號(hào)沉降監(jiān)測(cè)點(diǎn)的沉降與相應(yīng)直徑R進(jìn)行對(duì)比分析，如圖9所示.

對(duì)比以上兩組試驗(yàn)的沉降測(cè)量結(jié)果，該結(jié)果表明位于塌陷區(qū)域的土體沉降隨著土洞破壞的加劇，上覆土的沉降槽由近似圓弧形逐漸發(fā)生變化，土洞到達(dá)極限平衡狀態(tài)時(shí)其沉降可近似看作拋物線(xiàn)型. 分析其沉降過(guò)程，發(fā)現(xiàn)在土洞發(fā)生坍塌事故前，地面由于土洞發(fā)育引發(fā)的沉降速度逐漸減小，在土洞覆土發(fā)生完全坍塌時(shí)，其沉降量發(fā)生突變從而引發(fā)地面塌陷，該現(xiàn)象符合實(shí)際工程[7]中城市路面下土洞發(fā)育時(shí)地面沉降發(fā)展規(guī)律.

2.4? ?土壓力分布

當(dāng)土洞在發(fā)育擴(kuò)展的過(guò)程中，使用土壓力傳感器記錄土壓力隨土洞開(kāi)口直徑變化過(guò)程. 試驗(yàn)設(shè)計(jì)綜合考慮土壓力盒的靈敏度以及設(shè)置土壓力盒對(duì)試驗(yàn)結(jié)果的擾動(dòng)情況，僅對(duì)于600 mm埋深的土洞模型，在埋深100 mm的位置布置一排土壓力盒，土壓力盒平面位置如圖2所示，將土壓力盒從左至右編號(hào)為1～5號(hào). 土洞中心測(cè)點(diǎn)為3號(hào)測(cè)點(diǎn)，發(fā)育過(guò)程中記錄各土壓力盒的測(cè)試值P，在豎直方向上土洞各步發(fā)育后土壓力測(cè)量值如圖10所示. 根據(jù)不同點(diǎn)位的測(cè)試結(jié)果，土洞發(fā)育過(guò)程中，洞頂土壓力明顯小于洞周的土壓力. 為排除土壓力監(jiān)測(cè)過(guò)程中的無(wú)關(guān)因素，對(duì)稱(chēng)位置上各測(cè)點(diǎn)的數(shù)據(jù)取平均，測(cè)量值與發(fā)育步驟間關(guān)系如圖11所示.

土壓力測(cè)試試驗(yàn)結(jié)果表明，當(dāng)填土下方產(chǎn)生土洞時(shí)，土洞軸心處土壓力變化最為明顯，隨著半徑逐漸擴(kuò)大，土壓力值顯著減小. 其次當(dāng)土洞發(fā)育至直徑為32 cm時(shí)，2、4號(hào)土壓力盒位置在水平方向位于該范圍內(nèi)，其土壓力降低速度加快. 試驗(yàn)1、5號(hào)土壓力盒在整個(gè)發(fā)育過(guò)程中位于土洞水平投影范圍之外，第1步發(fā)育完成后其豎向土壓力上升26%，并且在其后的試驗(yàn)過(guò)程中變化較小.

3? ?數(shù)值模型計(jì)算分析

3.1? ?數(shù)值模型建立方法

對(duì)于試驗(yàn)所取得的實(shí)測(cè)數(shù)據(jù)，選取N1組試驗(yàn)結(jié)果，建立相對(duì)應(yīng)的數(shù)值模型，其模型相關(guān)參數(shù)如表2所示. 考慮三維模型的對(duì)稱(chēng)性，取1/4模型進(jìn)行有限元計(jì)算，取摩爾庫(kù)侖模型為計(jì)算本構(gòu)模型. 建立模型如圖12所示. 該模型共計(jì)17 846個(gè)三維計(jì)算單元，其內(nèi)部側(cè)面設(shè)置為對(duì)稱(chēng)邊界條件，外部側(cè)面則限制水平方向位移，模型下部采用固定邊界. 將表3中N1組別的坍塌數(shù)據(jù)代入有限元模型中，對(duì)照試驗(yàn)過(guò)程改變模型，對(duì)比有限元計(jì)算模型沉降計(jì)算結(jié)果與試驗(yàn)實(shí)測(cè)沉降數(shù)據(jù).

以土洞為中心，取室內(nèi)模型試驗(yàn)C3、C4點(diǎn)沉降測(cè)量值與數(shù)值模型試驗(yàn)對(duì)應(yīng)點(diǎn)沉降值進(jìn)行對(duì)比，數(shù)值模擬結(jié)果與實(shí)際測(cè)量結(jié)果對(duì)比如圖13所示. 結(jié)果表明，對(duì)于試驗(yàn)用砂土材料，數(shù)值模擬方法對(duì)于橢球形土洞正中心處其模擬效果較好，C3點(diǎn)位試驗(yàn)與模擬結(jié)果誤差不超過(guò)3%，C4點(diǎn)的誤差為7%. 根據(jù)模型試驗(yàn)與試驗(yàn)結(jié)果對(duì)比，在土洞發(fā)生破壞前各步計(jì)算中，數(shù)值模擬較好地反映了土洞發(fā)育過(guò)程中塌陷土體的變形特征，在工程中具有一定參考價(jià)值.

根據(jù)計(jì)算模型進(jìn)一步選取長(zhǎng)沙一處路面塌陷事故為例，其現(xiàn)場(chǎng)塌坑深度3 m，地表產(chǎn)生近似圓形塌孔，塌孔直徑為2.8 m. 由臨近地質(zhì)資料推斷塌陷事故土洞在地表填土層中，物理力學(xué)參數(shù)見(jiàn)表5.

按照上述方法建立有限元計(jì)算模型，地層為雜填土層，經(jīng)驗(yàn)取值認(rèn)為對(duì)于砂質(zhì)土完整系數(shù)取0.5，砂質(zhì)黏土為0.6，粉質(zhì)黏土以及礫石層則取為0.8. 因填土層粒徑分布復(fù)雜，考慮其類(lèi)似礫石土，取其完整系數(shù)為0.8，以確定其發(fā)育模式. 考慮塌陷土體對(duì)土洞的填埋，取土洞的近似埋深為3.5 m.

當(dāng)土洞在水平方向半徑逐漸發(fā)育時(shí)，取塑性區(qū)貫通至地表面時(shí)土洞發(fā)生塌陷，最終塌陷時(shí)塑性區(qū)分布如圖14所示. 由其地面塑性區(qū)分布，可以預(yù)測(cè)其地面塌陷形式為土洞頂部發(fā)生塌穿，周邊地面產(chǎn)生環(huán)狀塑性區(qū)，該環(huán)狀塑性區(qū)產(chǎn)生拉裂破壞. 數(shù)值模擬結(jié)果與實(shí)際工程均顯示地面發(fā)生塌穿的范圍為圓形，數(shù)值計(jì)算塌陷半徑為1.5 m，實(shí)際事故發(fā)生的塌穿范圍半徑約為1.4 m，模擬結(jié)果較好地反映了工程中的事故現(xiàn)場(chǎng).

3.2? ?數(shù)值模擬參數(shù)分析

以上針對(duì)三維試驗(yàn)?zāi)Ｐ蜏y(cè)試結(jié)果與實(shí)際工程實(shí)例建立了相關(guān)的數(shù)值計(jì)算模型，結(jié)果表明了本文數(shù)值模擬計(jì)算模型的可靠性. 在實(shí)際工程中，探討土體材料參數(shù)對(duì)塌陷事故的影響具有重要的工程意義. 因此，以地表填土的彈性模量E、泊松比μ、黏聚力c、內(nèi)摩擦角φ作為變量進(jìn)行分析，填土物理參數(shù)的經(jīng)驗(yàn)取值范圍見(jiàn)表6.

正交試驗(yàn)設(shè)計(jì)是探討材料參數(shù)對(duì)土洞坍塌半徑影響的一種有效分析方法[15]，該方法針對(duì)填土參數(shù)的經(jīng)驗(yàn)取值，將填土參數(shù)在取值范圍內(nèi)進(jìn)行變量水平劃分，以此作為變量的編碼依據(jù)，將因素水平劃分為上星號(hào)臂γ、上水平、零水平、下水平、下星號(hào)臂-γ五個(gè)等級(jí)，其具體劃分如表7所示.

為了分析各因素對(duì)土洞最終塌陷時(shí)半徑的影響，采用二次正交組合試驗(yàn)進(jìn)行分析. 選用正交表L8（27） 為二因素設(shè)計(jì)依據(jù)，其星號(hào)試驗(yàn)次數(shù)為8次，零水平試驗(yàn)次數(shù)為1，總試驗(yàn)次數(shù)為17次. 試驗(yàn)所得的回歸方程包含七項(xiàng)，其中Z1～Z4為因素X1～X4的一次項(xiàng)，交互項(xiàng)為Z1Z2、Z1Z3、Z2Z3，平方項(xiàng)為Z ′1、Z ′2、Z ′3、Z ′4，則擬合公式為：

y=a+■bjZj+b12Z1Z2+b13Z1Z3+b23Z2Z3+■bjjZ ′j （3）

將上述變量水平依據(jù)二次正交計(jì)算模型輸入有限元計(jì)算軟件，利用湘江路路面塌陷事故建立的模型進(jìn)行有限元模擬，分別得到不同條件下土洞發(fā)育的極限半徑R的值，計(jì)算結(jié)果見(jiàn)表8.

將試驗(yàn)結(jié)果進(jìn)行回歸計(jì)算，忽略參數(shù)中影響極小的項(xiàng)，得到土洞極限半徑關(guān)于以上表6中的4個(gè)參數(shù)的回歸方程如式（4），將上述回歸變量代回可得到極限半徑相對(duì)于各參數(shù)值的回歸方程如式（5）.

y=-0.122Z1-0.131Z2+0.228Z3+0.209Z1Z2+

0.169Z1Z3+0.174Z2Z3-0.23Z ′1-0.211Z ′2-

0.211Z ′3-0.222Z ′4+2.252? ? ? ? （4）

y=-5.726-0.241X1+6.109X2+0.483X3+0.542X4+

0.957X1X2+0.031X1X3+1.274X2X3-0.027X21-

38.632X22-0.062X23-0.011X24? ? ? ? （5）

針對(duì)計(jì)算結(jié)果，當(dāng)E=6 MPa，μ=0.3時(shí)，黏聚力與內(nèi)摩擦角與土洞極限半徑Rmax繪制為圖15. 圖15分別顯示了土洞的極限半徑受內(nèi)摩擦角和黏聚力的影響，當(dāng)內(nèi)摩擦角增大時(shí)，在表6的取值范圍內(nèi)，土洞的半徑呈現(xiàn)先增大后減小的關(guān)系，在φ = 24°時(shí)其半徑取得極大值;當(dāng)填土黏聚力增大時(shí)，土洞極限半徑在c = 8.5 kPa處取得極大值. 該擬合結(jié)果顯示在一定范圍內(nèi)增加填土材料參數(shù)有助于增強(qiáng)土洞的穩(wěn)定，從而使塌陷時(shí)土洞的半徑更大，但若繼續(xù)提高填土的強(qiáng)度，其土體的塌陷半徑反而有所降低. 由此，在道路建設(shè)中應(yīng)當(dāng)針對(duì)材料進(jìn)行多次試驗(yàn)以選取合適強(qiáng)度的填土材料以防止該類(lèi)事故的發(fā)生.

4? ?結(jié)? ?論

本文通過(guò)室內(nèi)模型試驗(yàn)，重點(diǎn)研究了土洞發(fā)育致塌過(guò)程中土洞形狀、地表沉降量與其對(duì)應(yīng)土壓力變化數(shù)據(jù)，分析試驗(yàn)結(jié)果得出以下結(jié)論：

1）對(duì)于同一填土材料，當(dāng)土洞的水平向直徑發(fā)展時(shí)，其豎向高度與水平向直徑的比值符合普氏理論提出的幾何規(guī)律.

2）對(duì)于松散填土中埋深不同的土洞，其塌陷形式主要分為壇形與筒式塌陷兩類(lèi).

3）土洞產(chǎn)生引發(fā)地面沉降的監(jiān)測(cè)數(shù)據(jù)結(jié)果表明，在土洞坍塌之前沉降發(fā)育較緩慢，塌陷為突發(fā)式失穩(wěn).

4）土洞發(fā)育過(guò)程中，上覆土發(fā)生明顯的應(yīng)力轉(zhuǎn)移現(xiàn)象，土洞外側(cè)豎向土壓力發(fā)生明顯升高，其內(nèi)側(cè)豎直土壓力明顯降低.

5）利用有限單元法結(jié)合試驗(yàn)結(jié)果建立了模擬土洞發(fā)育的模型，以此為基礎(chǔ)通過(guò)正交回歸試驗(yàn)分析了各參數(shù)項(xiàng)對(duì)塌陷半徑的影響.

參考文獻(xiàn)

[1]? ? 萬(wàn)志清，秦四清，李志剛，等. 土洞形成的機(jī)理及起始條件[J]. 巖石力學(xué)與工程學(xué)報(bào)，2003，22（8）：1377—1382.

WAN Z Q，QIN S Q，LI Z G，et al. Formation mechanism and initial condition of soil cavity[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22（8）：1377—1382. （In Chinese）

[2]? ? KARL T. Theoretical soil mechanics[M]. Hoboken，NJ，USA：John Wiley & Sons，Inc，1943：66—76.

[3]? ? 趙明華，吳家繼，何瑋茜，等. 基于土拱效應(yīng)的樁承式路堤承載變形計(jì)算研究[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，43（3）：135—141.

ZHAO M H，WU J J，HE W X，et al. Researches on bearing mechanism and deformation of pile-supported embankments based on soil arch effect[J]. Journal of Hunan University （Natural Sciences），2016，43（3）：135—141. （In Chinese）

[4]? ? 李濤，張子真，宗露丹. 地下空洞引起土質(zhì)地層地陷的形成機(jī)制與預(yù)測(cè)研究[J]. 巖土力學(xué)，2015，36（7）：1995—2002.

LI T，ZHANG Z Z，ZONG L D. Study of formation mechanism and prediction of sinkholes in soil stratum induced by subterranean cavity[J]. Rock and Soil Mechanics，2015，36（7）：1995—2002. （In Chinese）

[5]? ? 朱才輝，李寧. 黃土高填方地基中暗穴擴(kuò)展對(duì)機(jī)場(chǎng)道面變形分析[J]. 巖石力學(xué)與工程學(xué)報(bào)，2015，34（1）：198—206.

ZHU C H，LI N. Analysis of airstrip deformation due to expansion of hidden cavities in loess filled high embankment[J]. Chinese Journal of Rock Mechanics and Engineering，2015，34（1）：198—206. （In Chinese）

[6]? ? SAGASETA C. Analysis of undraind soil deformation due to ground loss[J]. Géotechnique，1987，37（3）：301—320.

[7]? ? 周乾. 基于尖點(diǎn)突變理論的城區(qū)孕育型土洞地陷判據(jù)研究[D]. 長(zhǎng)沙：湖南大學(xué)，2019：38—48.

ZHOU Q. Study on the criterion of ground subsidence of pregnant soil cave in urban area based on cusp catastrophe theory[D]. Changsha：Hunan University，2019：38—48. （In Chinese）

[8]? ? FEHDI C，BAALI F，BOUBAYA D，et al. Detection of sinkholes using 2D electrical resistivity imaging in the Cheria Basin （north-east of Algeria）[J]. Arabian Journal of Geosciences，2011，4（1/2）：181—187.

[9]? ? 陳其志，郭生根，徐長(zhǎng)節(jié)，等. 砂土中松動(dòng)土壓力及松動(dòng)區(qū)位移破壞形式的試驗(yàn)研究[J]. 中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，50（1）：108—117.

CHEN Q Z，GUO S G，XU C J，et al. Trapdoor model tests on loosening earth pressure and failure mode of loosening zone in sand[J]. Journal of Central South University （Science and Technology），2019，50（1）：108—117. （In Chinese）

[10]? XIAO X X，XU M，DING Q Z，et al. Experimental study investigating deformation behavior in land overlying a Karst cave caused by groundwater level changes[J]. Environmental Earth Sciences，2018，77（3）：1—11.

[11]? 蘇永華，周乾，蹇宜霖. 孕育型土洞地陷判據(jù)研究[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2020，47（5）：39—47.

SU Y H，ZHOU Q，JIAN Y L. Study on the criterion of land subsidence for gestation type soil cave[J]. Journal of Hunan University （Natural Sciences），2020，47（5）：39—47. （In Chinese）

[12]? 梁波，陳弘楊，聶影，等. 不同扶壁間距的倉(cāng)扶式支擋結(jié)構(gòu)土壓力分布[J]. 中國(guó)公路學(xué)報(bào)，2020，33（4）：24—31.

LIANG B，CHEN H Y，NIE Y，et al. Distribution of earth pressure of warehouse-counterfort retaining wall with different counterfort spacing[J]. China Journal of Highway and Transport，2020，33（4）：24—31. （In Chinese）

[13]? YANG X L，HUANG F. Three-dimensional failure mechanism of a rectangular cavity in a Hoek-Brown rock medium[J]. International Journal of Rock Mechanics and Mining Sciences，2013，61：189—195.

[14]? 王濱，賀可強(qiáng). 巖溶塌陷臨界土洞的極限平衡高度公式[J]. 巖土力學(xué)，2006，27（3）：458—462.

WANG B，HE K Q. Study on limit equilibrium height expression of critical soil cave of Karst collapse[J]. Rock and Soil Mechanics，2006，27（3）：458—462. （In Chinese）

[15]? 蘇永華，孫旺，方硯兵. 基于正交試驗(yàn)的隧道開(kāi)挖面位移釋放分析[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，44（11）：156—163.

SU Y H，SUN W，F(xiàn)ANG Y B. Analysis on displacement release of tunnel face based on orthogonal tests[J]. Journal of Hunan University （Natural Sciences），2017，44（11）：156—163. （In Chinese）