當我們乘坐地鐵或駕車穿行隧道時，很少會想到支撐在我們頭頂的巖土經歷了怎樣復雜的力學平衡。在隧道工程，無論是采用盾構法還是新奧法（NATM），開挖過程中掌子面的穩定性都是關乎工程安全與效率的核心。如果掌子面支撐壓力不足，前方土體可能發生失穩坍塌；而壓力過高，則可能造成不必要的工程浪費甚至引發地表隆起。因此，如何精確評估和設定平衡掌子面土體所需的支撐壓力，是隧道工程師們孜孜以求的目標。

長久以來，評估這一壓力的主流理論方法，大多借鑒了筒倉效應的思想，即假設破壞區由兩條垂直的滑移線圍成，形成一個類似“筒倉”的區域。然而，無論是先前的隧道模型試驗還是離心機試驗都表明，在松散地層中，實際的滑移破壞線往往是曲線而非直線。強行將其簡化為直線，無疑會引入誤差。更棘手的是，在傳統的地壓理論中，側向土壓力系數是一個關鍵的、但通常難以準確確定的參數。這兩個問題——不真實的直線滑移線假設和不確定的側向土壓力系數——共同限制了現有理論評估的精度。此外，基于上限定理的數值分析方法雖然功能強大，但需要構建復雜的靜力許可速度場并進行大規模矩陣運算，在實際工程設計中應用不便。面對這些挑戰，亟需一種既能更真實反映物理現象，又便于工程應用的新評估方法。

為了應對上述挑戰，日本埼玉大學的研究團隊在《Computers and Geotechnics》上發表了一項研究，提出了一種基于下限定理、利用曲線滑移線評估掌子面穩定性的新方法。研究團隊從基本的力學平衡原理出發，構建了一個由一條垂直滑移線（g）和一條等角對數螺旋曲線滑移線（h）圍成的二維破壞區模型。在這個破壞區內，他們提取了一個梯形非對稱微元體，對其建立豎向力、水平力以及力矩的平衡方程。通過對微元體上應力的張量方向進行分析，并結合滑移線上的應力狀態關系（基于莫爾-庫侖破壞準則），他們成功地推導出了控制該力學系統的微分方程。這個方法的巧妙之處在于，它無需預設側向土壓力系數，也無需像傳統上限分析法那樣進行復雜的空間離散化。通過數值求解這個微分方程，可以得到掌子面土壓力的下限解，從而評估其穩定性。研究還采用了對數螺旋曲線來模擬外滑移線，這與土力學中考慮三軸剪切角的常規做法一致，使得模型更能貼合實際觀測到的曲線破壞形態。

為了驗證所提方法的合理性與適用性，研究人員進行了多組參數計算，并與先前的離心模型試驗結果進行了對比。

研究首先詳細闡述了基于曲線滑移線的掌子面壓力估算理論。該方法本質上是一種極限平衡法，其核心創新是用曲線來刻畫與開挖相關的塌陷區。如圖2所示，破壞區由掌子面上方的垂直滑移線g和一條等角對數螺旋線h圍成。通過建立并求解微元體（如圖3所示）的平衡方程組，可以導出描述現象的微分方程，從而構成極限平衡法的下限解。對數螺旋線的形狀參數（m, z）由覆土厚度C、隧道直徑D以及地表破壞區寬度L等幾何條件決定，確保了曲線能真實地連接地表與掌子面下部邊緣。

在確定的滑移線框架下，研究對一個位于破壞區內部、以滑移線g為上底、曲線h為下底的梯形非對稱微元體進行了細致的力學分析。如圖3所示，分別考慮了微元體上豎向力、水平力的平衡以及繞其幾何中心的力矩平衡，并推導出相應的平衡方程（文中式6, 7, 8）。通過數學變換，這些方程最終被整理為關于角度θ的微分方程形式（文中式9, 10, 11），為后續數值求解奠定了基礎。

為了封閉方程組，研究引入了滑移線上的應力條件。基于莫爾-庫侖破壞準則，作用在滑移線g和h上的剪應力與法向應力滿足關系（文中式12）。更重要的是，通過分析滑移線g上的應力狀態（如圖4, 5, 6所示），并結合覆土壓力為γC sinθ這一條件，可以推導出該滑移線上法向應力σ_g與剪應力τ_g之間的一個重要線性關系（文中式14），其中系數ζ_g = 2 cot(π/2 - φ)。這個關系式的引入，巧妙地避免了對難以確定的側向土壓力系數的依賴，這是本方法的關鍵優勢之一。

通過數值求解上述微分方程組，研究對提出的方法進行了系統驗證。計算結果表明，所提出的解能夠正確反映摩擦材料由于覆土和摩擦阻力而產生的行為特性。尤為重要的是，計算得到的掌子面穩定性指標與先前一項離心模型試驗所獲得的結果能夠充分且準確地吻合。這強有力地證明了，采用曲線滑移線模型能夠更真實地模擬實際破壞模式，從而提高掌子面土壓力評估的精度。

本研究成功開發并驗證了一種用于評估隧道掌子面穩定性的新理論方法。該方法的核心貢獻與重要意義體現在以下幾個方面：

綜上所述，這項研究為隧道掌子面穩定性分析提供了一種新穎、嚴謹且更可靠的極限平衡方法。它不僅增強了我們對掌子面破壞機制的理論理解，而且為隧道工程，特別是盾構掘進中的泥水壓力或土壓平衡設定，提供了更具指導意義的理論工具，有助于實現更安全、更經濟的設計與施工。未來，研究者計劃將這一方法從目前的二維問題推廣到更復雜的三維情況，以拓展其應用范圍。