李 波,方俊杰,馮志國
(國家林業(yè)和草原局昆明勘察設計院 昆明市 650216)
我國是一個疆域遼闊、地質(zhì)環(huán)境十分復雜的國家,特別是西南部地區(qū),該區(qū)域多為山地、高原,地質(zhì)作用強烈,山勢陡峻、溝壑叢生,使得該區(qū)域的基礎設施建設尤為困難。我國地質(zhì)災害種類繁多,對公路工程的建設危害也巨大。涉及山區(qū)公路工程安全的主要有邊坡坍塌、崩塌、滑坡、泥石流等,其中公路路塹邊坡的穩(wěn)定性導致的災害和次生災害已經(jīng)成為主要地質(zhì)災害之一。由于我國受太平洋板塊、印度板塊和菲律賓海板塊的擠壓,地震斷裂帶十分發(fā)育,特別是西部西南部地區(qū),地震活動頻發(fā)。隨著我國近年來對西部地區(qū)投資的加大,大量公路工程將建設在西南部地震區(qū),地震作用下公路工程的穩(wěn)定性要求顯得越來越重要。由于西南部特殊的地形條件,公路建設不可避免地會造就大量的路塹邊坡,因此地震作用下的公路路塹邊坡的穩(wěn)定問題已經(jīng)成為公路工程穩(wěn)定性研究中的重點和難點問題[1-2]。
目前地震荷載作用下公路路塹邊坡的動力響應分析方法主要:大型振動臺模型試驗分析法、滑塊分析法、擬靜力法、有限單元法、有限差分法、動力時程分析法[3]等,這些方法都有各自運用領域和優(yōu)缺點。公路路塹邊坡地震作用下的動力響應包括加速度、位移、速度、應力響應等。探討各參數(shù)響應規(guī)律對地震荷載下邊坡的穩(wěn)定性分析和抗震治理有著極為重要的意義,這些不同參數(shù)的變化也直接或間接地反應了邊坡在地震作用下的變化規(guī)律。
以西南地區(qū)某二級公路均質(zhì)路塹邊坡為例,對該邊坡分不同臺階并建立相應的邊坡模型,輸入水平和豎直耦合的實測地震波,并對該路塹邊坡在實測地震波作用的動力響應譜規(guī)律進行研究??偨Y加速度、速度、位移、應力等時程曲線的響應規(guī)律,為公路路塹邊坡的抗震設計提供理論依據(jù)。
該公路路塹邊坡位于瀾滄江沿江二級公路K90+080~K90+300段右側,右側原為耕地,路線從耕地下方挖方經(jīng)過,設計坡比為1∶1。據(jù)詳勘報告,0~2m的耕植土、素填土,2~30m均為粉質(zhì)黏土,土體重度為19.6kN/m3,彈性模量為8.6MPa,泊松比為0.365,黏聚力為31.5kPa,內(nèi)摩擦角為17.2°。段落長度約220m,目前總體穩(wěn)定,段內(nèi)邊坡最高處為20m左右,最低處為16m左右,全段平均坡高約18m左右。由于該段邊坡沿路線方向較長、土質(zhì)均勻,為簡化模型,提高軟件計算效率,可以將這樣縱向很長且土質(zhì)材料均勻的路塹邊坡的三維變形簡化為二位平面變形問題。簡化后的計算模型坡比為1∶1,坡高取平均坡高約18m。根據(jù)《中國地震動峰值加速度區(qū)劃圖》(GB 18306—2015),該段路線所經(jīng)區(qū)域的設計基本地震加速度值為 0.20g,根據(jù)《中國地震動反應譜特征周期區(qū)劃圖》(GB 18306—2015),該段路線所經(jīng)區(qū)域的反應譜特征周期為 0.4s,抗震設防烈度為8度[4]。
為探尋邊坡臺階對該路塹邊坡的動力響應規(guī)律,在邊坡高度、設計坡比不變的情況下,分別模擬設置了一個臺階的二級邊坡和兩個臺階的三級邊坡,并與不設置臺階的單級邊坡進行對比,希望由此探索邊坡設立臺階后的動力響應規(guī)律和動力破壞模式,為多臺階邊坡地震作用下錨固提供理論依據(jù)[5]。
采用有限元軟件ANSYS對以上三種情況建立二維分析計算模型,對模型劃分有限單元網(wǎng)格,如圖1~圖3所示,并在邊坡坡面和坡體內(nèi)部設置相應的監(jiān)測點。
圖1 一級邊坡網(wǎng)格模型及監(jiān)測點
圖2 二級邊坡網(wǎng)格模型及監(jiān)測點
圖3 三級邊坡網(wǎng)格模型及監(jiān)測點
路塹邊坡底部施加豎直和水平兩個方向的約束,坡前和坡后施加豎直的單向約束,邊坡坡頂和邊坡臨空面豎直和水平兩個方向均不施加約束,用于模擬坡面和坡頂土體的變形和破壞。
選擇世界上第一條成功記錄全過程的地震波—1940年5月18日位于加州南部的EI-Centro地震波,該地震波完整記錄時間為54s,計算截取前段最具代表性的20s,如圖4。輸入地震波時,采用豎直和水平的雙向耦合輸入,分別探討輸入地震波后不同臺階數(shù)均質(zhì)公路路塹邊坡的坡面方向動力響應參數(shù)的變化規(guī)律。
圖4 EI波水平加速度時程曲線
為研究路塹邊坡在耦合輸入地震波時坡面方向的動力響應規(guī)律,分別在不同臺階數(shù)邊坡的坡面設置不同監(jiān)測點,如圖1~圖3,研究其加速度、速度、位移的變化,總結其動力響應規(guī)律,由此為地震作用下多臺階邊坡的治理提供理論依據(jù)。限于文章篇幅未列出監(jiān)測點的參數(shù)時程曲線,僅將時程曲線上各分量最大值和最小值統(tǒng)計于表1~表3。
表1 無臺階邊坡不同監(jiān)測點的最值統(tǒng)計表
表2 一臺階邊坡不同監(jiān)測點的最值統(tǒng)計表
根據(jù)該公路路塹邊坡坡面各監(jiān)測點的速度、加速度、位移時程曲線,統(tǒng)計各曲線的極值,如表1~表3,分析對比數(shù)據(jù),可得出以下結論:
表3 二臺階邊坡不同監(jiān)測點的最值統(tǒng)計表
(1)對于無臺階的一級邊坡
邊坡水平位移、速度、加速度分布大致呈現(xiàn)沿坡面向上逐漸增大的趨勢,坡腳位置最小、坡中次之、坡頂點處最大,說明邊坡對所輸入的地震波有沿坡面向上放大的趨勢,稱之為邊坡地震作用下的臨空面放大效應。
(2)對于一臺階的二級邊坡
①坡面上各突出點的位移、速度、加速度相對于坡腳點都有所增大,區(qū)別在于,各參數(shù)極值點不是出現(xiàn)在坡頂點,而是出現(xiàn)在第一級邊坡頂2點處,這說明一臺階邊坡最易破壞點為第一級邊坡的頂點,在治理時應該重點治理第一級,所對應的錨桿錨索長度也應較長。
②第二級邊坡坡腳3點有相反的趨勢,此點的速度、加速度最小,說明振動反應最輕微,此點為邊坡的安全點。
③計算邊坡三量最大值和坡腳點的比值,也就是邊坡各參數(shù)的放大系數(shù),并與無臺階一級邊坡進行對比。位移放大系數(shù)由1.816降到1.427,速度放大系數(shù)由1.295降到1.210,加速度放大系數(shù)由1.251降到1.231,由此證明分臺階處理后對控制邊坡位移、速度和加速度是有利的。
(3)對于二臺階的三級邊坡
①二臺階三級邊坡每一級的頂點,即2、4、6點相對于坡腳點1,速度、加速度、位移三量都有放大趨勢,但和一臺階邊坡有所不同,最大的位移、速度、加速度均出現(xiàn)在坡頂點6,其次是2點,增量最小的是4點,說明第三級最危險,其次是第一級、第二級。錨固設計時可由此確定主次關系,并相對應的優(yōu)化錨桿錨索設計長度,提高治理的效率和經(jīng)濟性。
②后緣點3和點5相對于坡腳點1,位移有增大趨勢,但增大幅度較小,速度和加速度減弱,這是由此處應力集中導致的,滑動體可能由此出露。
③將無臺階邊坡、一臺階邊坡和二臺階邊坡的各參數(shù)放大系數(shù)[6]進行比較可以發(fā)現(xiàn),位移放大系數(shù)分別為1.816、1.427和1.285,速度放大系數(shù)分別為1.295、1.210和1.201,加速度放大系數(shù)分別為1.251、1.231和1.212。隨著臺階數(shù)的增加,放大系數(shù)呈現(xiàn)逐步減小的趨勢,由此可見增加臺階可以使邊坡受力更加均勻,有利于控制邊坡的位移。因此分臺階處理邊坡可以降低邊坡的地震動力反應,提高邊坡動力穩(wěn)定性。
為了更形象地觀察耦合輸入EI地震波后邊坡內(nèi)部應力的變化,給出了邊坡最后一個計算子步的水平向應力云圖和平均應力云圖,如圖5、圖6,通過應力云圖可以很好地驗證上述結論。
圖5 水平應力云圖
通過分析邊坡的各計算子步的水平向應力云圖和平均應力云圖的變化規(guī)律,可得出如下結論:
(1)對于應力的分布,各臺階邊坡均發(fā)生了應力集中的現(xiàn)象,差別在于集中的位置和程度不一。無臺階的一級邊坡應力首先集中在坡頂和坡腳以下,主要沿著坡體內(nèi)部一倍坡高處集中,大體呈圓弧形趨勢發(fā)展,這表明隨著地震荷載的輸入,邊坡有可能沿著此應力帶破壞,最終從負應力較為集中的坡腳出露。
圖6 平均應力云圖
(2)對于一臺階的二級邊坡:應力分布動力響應規(guī)律在坡體內(nèi)部大致與無臺階邊坡相同,應力均首先集中在坡頂和坡腳以下一倍坡高處,并由此沿坡體內(nèi)部呈圓弧形發(fā)展。應力集中區(qū)別在于坡面處,其中1、2、4點應力集中明顯,為坡面危險點。3點應力較小,為安全點。一級邊坡坡面應力整體大于二級邊坡坡面應力,邊坡錨固設計時,應重點治理一級邊坡,二級次之。
(3)對于二臺階的三級邊坡:由于設置了兩個臺階,在地震荷載作用下受力最為均勻,在輸入地震波以后沒有出現(xiàn)明顯的塑性帶,坡腳以下一倍坡高處應力最為集中,坡體內(nèi)部塑性集中不明顯,呈現(xiàn)多點分布的趨勢。坡面應力分布也有不同,隨著地震動荷載的增加,三級邊坡坡頂處6點首先出現(xiàn)應力集中,其次是一級邊坡坡頂處2點,最后是二級邊坡坡頂處4點。由此進一步證明二臺階三級邊坡在地震作用下的最危險位置依次為第三級、第一級和第二級坡,與2.1節(jié)結論一致,邊坡錨固設計可由此確定主次關系并由此調(diào)整錨桿錨索設計長度。
通過對豎直和水平雙向耦合輸入地震波時三種不同邊坡模型各位置監(jiān)測點的各參數(shù)時程曲線和應力分布的分析,可得出以下結論:
(1)對于一臺階的二級邊坡:最易破壞點為第一級邊坡的頂點,治理時應該重點治理第一級,第一級的錨桿錨索長度也應較長。應力分布規(guī)律在坡體內(nèi)部大致與無臺階邊坡相同,首先集中在坡頂和坡腳以下一倍坡高處,并由此沿坡體內(nèi)部呈圓弧形發(fā)展。
(2)對于二臺階的三級邊坡:坡面應力分布隨著地震動荷載的增加,頂處6點首先出現(xiàn)應力集中,其次是2點,最后是4點。因此第三級最危險,其次是第一級、第二級。輸入地震波以后沒有出現(xiàn)明顯的塑性帶,坡體內(nèi)部塑性集中不明顯,呈現(xiàn)多點分布的趨勢。二臺階邊坡的錨固設計時可由此確定主次關系,并由此優(yōu)化錨桿錨索的設計長度,提高治理的效率和經(jīng)濟性。
(3)將無臺階邊坡、一臺階邊坡和二臺階邊坡的各參數(shù)放大系數(shù)進行比較可以發(fā)現(xiàn),位移放大系數(shù)分別為1.816、1.427和1.285,速度放大系數(shù)分別為1.295、1.210和1.201,加速度放大系數(shù)分別為1.251、1.231和1.212。隨著臺階數(shù)的增加,放大系數(shù)呈現(xiàn)逐步減小的趨勢,由此可見增加臺階可以使邊坡受力更加均勻,有利于控制邊坡的位移。因此分臺階處理邊坡可以降低邊坡的地震動力反應,提高邊坡動力穩(wěn)定性。