飽和土壤-結構動態相互作用中的摩擦接觸模型研究

摘要解讀
本研究針對飽和土壤與結構物之間復雜的動態接觸問題，提出了一種融合流體-固體耦合與接觸耦合的雙重控制模型。通過構建基于u-p雙場耦合的接觸控制算法，創新性地將有效應力原理與摩擦接觸狀態機相結合，有效解決了傳統接觸算法在孔隙介質中的數值不穩定問題。研究團隊通過四類典型工況的數值驗證，證實該模型在預測接觸界面"粘附-滑移-分離"狀態轉換方面具有顯著優勢，其計算結果與經典解析解、商業軟件模擬及實驗數據均呈現良好吻合。

研究背景與問題提出
在巖土工程領域，結構物與飽和土體的動態相互作用廣泛存在于邊坡滑動、樁基穿透、管道屈曲等工程場景中。由于土體固相骨架與孔隙流體的多物理場耦合特性，接觸界面同時面臨流體滲透與固體變形的雙重作用。傳統接觸模型主要針對單一介質系統開發，在處理飽和土體這類雙相耦合問題時存在顯著局限性。

具體表現為：（1）流體相（孔隙水）與固相（土顆粒）的動力學耦合尚未在接觸算法中得到充分體現；（2）現有接觸界面離散方法難以有效處理飽和土體特有的孔隙壓力傳遞問題；（3）動態加載條件下接觸狀態的瞬態轉換難以準確建模。這些問題導致現有模型在預測接觸界面能量耗散、摩擦力傳遞機制及多場耦合效應時存在較大誤差。

現有方法分析
當前接觸問題處理主要存在兩大技術路線：其一為特殊接觸單元法，通過 Goodman 單元、薄層單元等特殊構造單元直接嵌入接觸約束條件。此類方法物理意義明確，但存在單元類型單一、適用性受限的缺陷。其二為邊界約束轉化法，將接觸問題轉化為等效的邊界約束條件。這種方法通過空間離散轉換簡化了求解過程，但面臨接觸界面離散精度與算法穩定性之間的矛盾。

針對動態接觸問題，主流算法包括 Lagrange 乘子法、懲罰法和增廣拉格朗日法。Lagrange 乘子法雖能精確滿足接觸約束，但會引入額外的未知量導致計算規模激增。懲罰法則通過引入接觸罰參數實現約束，雖然計算形式簡單且無需附加變量，但其數值穩定性高度依賴罰參數的選取。增廣拉格朗日方法通過懲罰參數與約束條件的協同優化，有效緩解了傳統懲罰法的收斂問題，但算法復雜度較高。

在接觸界面離散方面，節點-線段（NTS）方法因具有簡單直觀、計算效率高等優勢被廣泛應用。然而該方法存在兩個固有缺陷：（1）接觸節點與主從段匹配存在幾何敏感性問題；（2）雙通帶算法在處理孔隙流體傳遞時存在數值振蕩。Mortar 方法通過弱形式約束有效規避了離散不匹配問題，但其需要構造雙空間基函數，導致計算成本顯著增加。

本文方法創新
研究團隊基于u-p雙場耦合理論，構建了適用于飽和土-結構系統的動態接觸模型。主要創新點體現在以下三個方面：

1. 動態接觸狀態機開發
創新性地將Mohr-Coulomb破壞準則擴展至雙場耦合系統，建立"粘附-滑移-分離"三態動態轉換模型。通過有效應力原理與摩擦強度理論的耦合分析，準確捕捉接觸界面在動態荷載下的狀態演化規律。研究特別考慮了孔隙水壓力對摩擦系數的弱化效應，當有效應力低于臨界值時自動觸發滑移狀態轉換。

2. 耦合NTS接觸元素設計
開發了具有雙場耦合特性的節點-線段接觸單元。該單元將固相位移場與液相孔隙壓力場進行統一編碼，通過共享節點實現位移連續性約束，通過獨立單元實現孔隙壓力的跳躍連續。特別設計了雙通帶算法，在位移場和孔隙壓力場分別設置獨立通帶，有效解決了傳統單通帶方法在處理雙場耦合時的數值振蕩問題。

3. 增廣拉格朗日穩定性增強
將傳統增廣拉格朗日方法改進為雙場耦合特化的穩定性控制算法。通過構建位移場與孔隙壓力場的聯合罰函數，在保證接觸界面穿透約束的同時，有效抑制了孔隙壓力場的數值振蕩。算法引入動態罰參數調整機制，根據計算過程自動優化罰參數，既保證了初始階段的計算穩定性，又實現了后期高精度求解的需求。

算法實現與驗證
研究采用兩階段懲罰算法實現接觸控制，第一階段基于位移通帶進行接觸檢測，第二階段根據孔隙壓力通帶調整接觸狀態。具體實施步驟包括：

（1）接觸界面離散：采用改進的NTS元素進行接觸面離散，每個接觸單元包含3個主節點和2個從節點，通過幾何插值算法實現主從段的高精度匹配。

（2）狀態機建模：基于Mohr-Coulomb準則建立接觸強度判斷矩陣，設置有效應力閾值（σ'c/2）和孔隙壓力閾值（u/u_max）。當有效應力低于閾值時觸發滑移狀態，孔隙水壓力超過臨界值時進入分離狀態。

（3）耦合控制方程：通過雙場耦合本構方程，建立包含位移、應變率、孔隙壓力和滲透率四場耦合的動力學方程組。特別在接觸界面處，耦合了固相位移連續條件和液相孔隙壓力跳躍條件。

數值驗證部分設計了四類典型工況進行模型驗證：

（1）剛性基礎下的固結沉降：對比傳統Terzaghi固結解，驗證雙場耦合模型在考慮孔隙水壓力變化時的沉降預測精度。結果顯示模型在排水路徑超過5個孔隙半徑時仍能保持高精度。

（2）移動荷載下的接觸響應：采用軌道車輛-路基系統進行時程分析，驗證模型在接觸面滑移帶形成與消失過程中的動態響應捕捉能力。計算結果與實驗數據吻合度達92%以上。

（3）循環荷載下的疲勞接觸：模擬海底管道在波浪荷載下的接觸行為，驗證模型在接觸面反復粘附-滑移-分離循環中的預測可靠性。通過2000次循環加載的數值模擬，接觸強度衰減曲線與實驗結果誤差控制在5%以內。

（4）三維接觸問題：構建三維樁-土接觸模型，驗證算法在復雜幾何形態下的適用性。對比傳統二維模型，新增的剪切變形模式使接觸力計算誤差降低至3.8%。

工程應用價值
本成果在多個工程領域具有顯著應用價值：（1）在基礎工程中，可精確預測樁基側摩阻力分布，特別是對存在液化風險的飽和砂土層；（2）在邊坡防護工程中，能準確模擬錨桿與巖體的接觸應力傳遞，對滑移面預測精度提升達40%；（3）在地下工程中，對隧道管片與圍巖接觸壓力的動態響應預測誤差小于5%，顯著優于傳統接觸模型。

研究局限性及改進方向
當前模型仍存在兩個待優化方向：（1）接觸界面局部化處理仍需改進，當接觸面積超過總接觸面積的30%時，收斂速度下降約25%；（2）孔隙壓力場在接觸區域的跳躍連續性假設存在一定誤差。后續研究計劃引入非連續伽遼金（NAG)方法，通過構造混合有限元單元實現更精確的場變量連續性控制。

結論
本研究成功構建了雙場耦合的動態摩擦接觸模型，通過四階段驗證流程（理論推導-算法實現-數值模擬-工程驗證）確保模型可靠性。創新性地將增廣拉格朗日方法應用于孔隙介質系統，解決了傳統接觸算法在雙場耦合中的數值不穩定問題。實驗數據顯示，在接觸面存在顯著孔隙壓力變化時，模型預測的接觸力誤差小于8%，較傳統方法提升約50%。該成果為地下結構物設計、邊坡穩定性分析及地下工程監測提供了新的理論工具和計算平臺。