交河故城崖體拉裂-傾倒破壞模式

袁炳祥，諶文武,，滕 軍，陳 銳

(1. 哈爾濱工業大學 深圳研究生院，廣東 深圳 518055；2. 蘭州大學 西部災害與環境力學教育部重點實驗室，蘭州 730000；3. 蘭州大學 土木工程與力學學院，蘭州 730000)

1 前 言

位于新疆吐魯番的交河故城遺址是第一批全國重點文物保護單位，在自然和人為環境的影響下，已經遭受了嚴重的破壞[1－2]，尤其是近些年來，崖體的破壞已經危害到本體的保護，引起了國家文物局的高度重視。為了避免整體破壞，國家文物局組織了敦煌研究院、蘭州大學、西北大學等單位進行了大量的科研工作[1－3]，并實施了相應的搶險加固工程。

通過大量的現場勘察發現，砂層嚴重風蝕和次生卸荷裂隙發育部位的崖體破壞較嚴重，是崖體加固的重要部分，要從根本上解決崖體的病害，必須深入研究其變形破壞機制。有限元軟件具有能夠考慮土體的非線性彈塑性本構關系、動態模擬土體的失穩過程及滑移面形狀等優點，被廣泛應用到邊坡巖土體工程[4－6]。本文運用有限元軟件ANSYS模擬崖體在差異風蝕和次生卸荷裂隙作用下的變形破壞模式。為了對比分析數值模擬結果，通過底面摩擦試驗模擬崖體拉裂-傾倒破壞演化過程。底面摩擦法是模擬重力的最常用方法之一，可以模擬整個破壞過程，并可以根據情況實時觀測試驗過程的各個細節，因此被應用于邊坡穩定性的研究[7－9]，但這種方法極少應用于土遺址保護研究。為了檢驗數值模擬的可靠性，依托交河故城典型崖體模型開展崖體變形破壞演化過程的物理模擬試驗，剖析砂層風蝕和次生卸荷裂隙對崖體變形破壞的影響，提出有效的崖體加固方案。

2 崖體拉裂-傾倒破壞模式

交河崖體主要由黃色粉質黏土層、淡黃色粉土層和黑色砂層組成，具有明顯的分層特征。對不同土層取原狀土樣進行物理性質試驗，得到3種土的粒度成分和力學參數，見表1。研究表明[10－11]，顆粒的粒徑在 0.5～0.075 mm 區間是最容易被風蝕的，圖1表明3種土的粒度成分在抗風蝕能力存在著較大的差異性。

表1 崖體土的力學參數Table1 Mechanic parameters of cliff soil

圖1 崖體粉質黏土、粉土和砂的累積曲線圖Fig.1 Cumulative curves of silty clay,silt and sand of cliff mass

吐魯番盆地春季多風，以西北風為主，最大風速可達 40 m/s以上。根據吐魯番 40年（1951－1991年）的氣象資料統計表明，8～10級大風年平均為36.2 d，沙暴年平均為11.2 d，大風持續時間最長為8 d，平均每天刮8.35 h[2]。

風是交河故城破壞的主要外動力，自始至終都參與各種病害的產生、發展，風蝕已成為交河故城最主要病害之一。為了定量地研究風蝕對交河故城崖體不同土層的影響，利用西部災害與環境力學教育部重點實驗室的多功能風洞實驗室對交河崖體粉質黏土、粉土和砂3種相關土進行風蝕試驗。依據吐魯番常見大風風速，依次選取18、28、30、40 m/s等4個風速的攜砂風，對3組5 cm ×5 cm ×5 cm立方體樣品試塊進行300、600、900 s連續吹蝕。試驗用砂采自甘肅民勤縣境內騰格里沙漠邊緣，平均粒徑為 0.228 mm。在某一風速下，其攜沙能力為一定值，風速越大攜沙能力越強，實驗結果表明，粉質黏土風蝕量最小接近于0。由圖2可以看出，在相同風速和時間下砂的風蝕量為粉土的數倍，為粉質黏土的數10倍。試驗所測得3組樣品在不同風速、不同吹蝕時間后的風蝕量，為交河故城崖體的差異性風蝕效應研究提供了有利的依據。圖中，1為粉土18 m/s風速的風蝕量；2為砂18 m/s風速的風蝕量；3為粉土28 m/s風速的風蝕量；4為砂 28 m/s風速的風蝕量；5為粉土30 m/s風速的風蝕量；6為砂30 m/s風速的風蝕量；7為粉土40 m/s風速的風蝕量；8為砂 40 m/s風速的風蝕量。

圖2 不同風速下粉土和砂的風蝕量與吹蝕時間的關系Fig.2 Relationships of wind erosion amount and time of silts and sands in different wind velocities

交河故城崖體主要有兩類破壞裂隙：構造裂隙和次生卸荷裂隙。次生卸荷裂隙具有傾角近直立、張開度大、水平延伸長、垂直延伸深等特點。裂縫切割崖體深度以4～8 m最多，切穿深達10 m以上者占13%，少數裂隙立面出露深度達18 m，已超過崖體高度的 1/2。裂隙多切穿粉土層，止于砂層，砂層為裂隙及崩塌破壞的控制層位。

由于交河崖體差異風蝕和次生卸荷裂隙兩者的共同作用，砂層在強烈的風蝕作用下被掏空，被卸荷裂隙切割的粉質黏土層在重力作用下拉張卸荷裂隙的根部，被切割土體將發生傾斜變形，從而可能拉斷卸荷裂隙的根部而使土體傾倒失穩。交河故城西北端崖體粉質黏土層向河谷方向的傾倒變形及坡腳大塊傾倒的土體均為崖體的拉裂-傾倒破壞提供了有利依據。

3 崖體破壞模式的有限元模擬

3.1 模型建立及參數選取

為研究砂層風蝕和次生卸荷裂隙對崖體破壞的影響和崖體破壞過程，運用有限元軟件ANSYS進行模擬。將崖體模型簡化為粉質黏土層和砂層互層，粉質黏土層和砂層均采用Solid82 8節點四邊形單元，每個結點有2個自由度，即X、Y兩個方向的位移。定義材料本構模型時采用Drucker- Prager理想彈塑性模型，該準則實際是對 Mohr-Coulomb準則的修正，在π平面上的屈服軌跡為一圓，其屈服函數[12]表示為

式中：I1為應力第一不變量；J2為第二偏應力張量；α、k與巖土體的黏聚力c和內摩擦角φ的關系為

式中：粉質黏土和砂的物理力學參數均用表 1數據。

為模擬卸荷裂隙和砂層掏蝕二者對崖體破壞影響，裂隙深度和砂層掏蝕深度取統計值的較大值，在模型頂部的粉質黏土層取8 m長卸荷裂隙、砂層被掏蝕深度取3 m。對模型進行自由劃分，其優點是在卸荷裂隙附近單元劃分面積比較小，劃分比較詳細，共有541個單元，1763個節點。在模型左邊界設置X方向位移約束，底邊界設置X、Y方向位移約束，模型整體施加重力荷載，如圖3所示。

圖3 ANSYS模擬崖體模型Fig.3 The ANSYS model of the cliff

3.2 結果分析

通過模擬結果得到 X、Y方向位移、總位移和塑性應變云圖，如圖4所示。被卸荷裂隙切割崖體從底部到頂部X方向位移逐漸增大，其最大位移為323.972 mm，見圖 4(a)綠色區域；其垂直方向位移，距離卸荷裂隙越遠越大，見圖 4(b)黃色區域；其總位移云圖，以卸荷裂隙根部為中心，呈圓環狀向崖體頂部邊界點發展，表明越向頂部邊界點其總位移越大，最大值為487.326 mm，見圖4(c)綠色區域。從以上崖體的位移可以看出，崖體以卸荷裂隙根部為中心發生傾倒。通過塑性應變云圖可以看出，沿卸荷裂隙塑性應變區已經貫通粉質黏土層，其最大塑性應變為 25.62%，見圖 4(d)紅色區域，說明崖體卸荷裂隙根部被拉張。以上結果都說明被卸荷裂隙切割的粉質黏土層塊體以卸荷裂隙根部為軸向外側傾倒破壞。實際調查中，張開度大于50 mm的裂隙占總數80%以上，某些裂縫局部張開達 1200 mm，若其底部有砂層掏蝕崖體就極易破壞。

4 底面摩擦模擬試驗

4.1 物理模擬基本原理

底面摩擦模擬試驗采用二維物理模型，具有試驗簡單直觀，效果顯著的特點，試驗主要考慮研究對象受重力作用的影響，利用物理模型底面所受到的摩擦力模擬原型自然邊坡所受重力[7]。其原理如圖5所示，即將研究對象的剖面制成模型平放在環形活動橡皮帶的平直段上，并使原剖面的深度方向與X方向一致，在橡皮帶轉動方向有一固定框架，當模型受其阻擋時，則模型與橡皮帶接觸面上每一點就形成摩擦阻力。根據圣維南原理[7－9]，當模型足夠薄時，認為摩擦力均勻作用在整個厚度上，相當于原型物體在天然狀態下受到重力作用。

圖4 ANSYS模擬結果云圖Fig.4 The contour images of ANSYS simulation results

圖5 底面摩擦模型試驗原理示意圖Fig.5 Sketch of basal contact friction

4.2 試驗過程及現象

按模型比例，在模型頂部的粉質黏土層做 80 mm卸荷裂隙，并在砂層掏蝕30 mm，見圖6(a)。試驗開始后，由于模型層與層之間有空隙，可以觀察到整個模型向皮帶轉動方向稍有壓密現象；經過1 h沿卸荷裂隙產生約15 mm的裂紋，1.5 h裂紋長約20 mm，2 h裂紋約30 mm，2.4 h裂隙貫通粉質黏土層，并且崖體頂部的裂隙張開度約6 mm，2.7 h崖體頂部的裂隙張開度約20 mm，土體已經傾倒破壞，試驗結束，如圖6所示。

底面摩擦試驗表明，在卸荷裂隙和砂層風蝕同時作用下，崖體沿卸荷裂隙根部產生拉裂，逐漸貫通粉質黏土層，并最終傾倒破壞。通過底面摩擦試驗較好地驗證了數值模擬結果，說明崖體在次生卸荷裂隙和差異風蝕作用下的破壞模式是拉裂-傾倒破壞。

圖6 底面摩擦試驗照片Fig.6 Basal contact friction test pictures

5 加固措施

數值計算和物理模擬結果均表明，在砂層風蝕和次生卸荷裂隙發育部位的崖體極易發生拉裂-傾倒破壞，崖體的穩定性主要由次生卸荷裂隙根部來維持，卸荷裂隙根部的應力和塑性應變隨著拉裂的發展而逐漸的積累。因此，加固的重點應是通過工程手段阻止崖體拉裂的繼續發展，加強臨空土體底部的支頂，減小土體的下滑力；給被卸荷裂隙切割土體足夠的連接，加強土體的抗拉張能力。在風蝕砂層部位采用土坯砌補，給上部粉質黏土層較好的支撐，并在砂層噴灑膠凝劑加強其抗風蝕能力；在次生卸荷裂隙及其根部采用楠竹加筋復合錨桿加固，并采用膠凝劑+粉煤灰+黏土漿液對裂隙灌漿，加強裂隙兩側土體的連接。交河故城搶險加固工程實踐證明這種加固措施科學有效。

6 結 論

（1）交河故城崖體沉積作用具有成層性，通過對不同土層的物理試驗和風洞試驗發現，粉質黏土、粉土和砂的粒度成分、力學性質和風蝕量都有較大差別，為崖體差異性風蝕提供了有利的依據。

（2）在卸荷裂隙和差異性風蝕二者同時作用情況下，被卸荷裂隙切割粉質黏土層塊體底部臨空。在重力作用下，拉張卸荷裂隙的根部被切割土體，以卸荷裂隙根部為軸將發生傾斜變形，可能拉斷卸荷裂隙的根部而使土體傾倒失穩。

（3）在砂層風蝕和次生卸荷裂隙作用下，崖體極易發生拉裂-傾倒破壞。按照砂層風蝕和次生卸荷裂隙對崖體破壞的影響，提出了崖體加固方法，對崖體的搶險加固具有重要意義。

[1]李最雄,王旭東,孫滿利. 交河故城保護加固技術研究[M]. 北京: 科學出版社,2008.

[2]孫滿利,王旭東,李最雄,等. 交河故城瞭望臺保護加固技術[J]. 巖土力學,2007,28(1): 163－168.SUN Man-li,WANG Xu-dong,LI Zui-xiong,et al.Technology of protection and reinforcement for observation platform in Ruins of Jiaohe[J]. Rock and Soil Mechanics,2007,28(1): 163－168.

[3]崔凱,諶文武,張景科,等. 多元層狀邊坡土體風蝕速率與微結構參數關系[J]. 巖土力學,2009,30(9): 2741－2746.CUI Kai,CHEN Wen-wu,ZHANG Jing-ke,et al.Relationships between microstructure parameters and wind erosion rate of multivariate layered soil in slope[J].Rock and Soil Mechanics,2009,30(9): 2741－2746.

[4]馬乙一. 應力重塑方法在治理邊坡中的應用研究[碩士學位論文D]. 合肥: 合肥工業大學,2006.

[5]徐力群,沈振中,俞凱加,等. 思林水電站通航建筑物高邊坡穩定分析和評價[J]. 巖土力學,2008,29(增刊):42－46.XU Li-qun,SHEN Zhen-zhong,YU Kai-jia,et al.Stability analysis and estimation of slope excavation of Silin hydropower station’s navigation lock[J]. Rock and Soil Mechanics,2008,29(Supp.): 42－46.

[6]李榮建,于玉貞,李廣信. 抗滑樁加固非飽和土邊坡三維穩定性分析[J]. 巖土力學,2008,29(4): 968－972.LI Rong-jian,YU Yu-zhen,LI Guang-xin. 3D global stability analysis of unsaturated soil slope reinforced with piles[J]. Rock and Soil Mechanics,2008,29(4): 968－972.

[7]蔡國軍,黃潤秋,嚴明,等. 反傾向邊坡開挖變形破裂響應的物理模擬研究[J]. 巖石力學與工程學報,2008,27(4): 811－817.CAI Guo-jun,HUANG Run-qiu,YAN Ming,et al.Physical simulation study on deformation and failure response of an anti-inclined slope during excavation[J].Chinese Journal of Rock Mechanics and Engineering,2008,27(4): 811－817.

[8]張登項,許強. 基于底摩擦試驗的錦屏一級水電站左岸巖石高邊坡變形機制研究[J]. 地質災害與環境保護,2008,19(1): 71－75.ZHANG Deng-xiang,XU Qiang. Research on deformation mechanism of the high rock-slope at left of Jingping first stage hydropower station by experiment of basal contacting friction[J]. Journal of Geological Hazards and Environment Preservation,2008,19(1):71－75.

[9]馮文凱,石豫川,柴賀軍,等. 緩傾角層狀高邊坡變形破壞機制物理模擬研究[J]. 中國公路學報,2004,17(3):32－36.FENG Wen-kai,SHI Yu-chuan,CHAI He-jun,et al.Study on mechanism of deformation failure of a low-angle bedded high slope with physical simulation method[J]. China Journal of Highway and Transport,2004,17(3): 32－36.

[10]李曉麗,申向東. 裸露耕地土壤風蝕躍移顆粒分布特征的試驗研究[J]. 農業工程學報,2006,22(5): 74－77.LI Xiao-li,SHEN Xiang-dong. Experimental study on the distribution characteristics of the saltation particle of aeolian sediment in bare tillage[J]. Transactions of the Chinese Society of Agricultural Engineering,2006,22(5): 74－77.

[11]李萬源,沈志寶,呂世華,等. 風蝕影響因子的敏感性試驗[J]. 中國沙漠,2007,(6): 984－992.LI Wan-yuan,SHEN Zhi-bao,Lü Shi-hua,et al.Sensitivity tests of factors influencing wind erosion[J].Journal of Desert Research,2007,(6): 984－992.

[12]王桂堯,胡振南,匡希龍. 紅砂巖路基強夯處理大變形數值模擬方法與效果分析[J]. 巖土力學,2008,29(9):2451－2456.WANG Gui-yao,HU Zhen-nan,KUANG Xi-long.Large-strain numerical simulation and experimental result research about improving red-sandstone embankment by dynamic compaction[J]. Rock and Soil Mechanics,2008,29(9): 2451－2456.