巖石蠕變的相場斷裂模型:理論框架與工程應用
《Computers and Geotechnics》:A phase field fracture model for rock creep: Theoretical framework and engineering applications
【字體:
大
中
小
】
時間:2026年03月02日
來源:Computers and Geotechnics 6.2
編輯推薦:
巖石蠕變裂縫擴展與損傷演化耦合機理研究。提出雙相場模型分離拉伸和壓縮-剪切裂縫,考慮蠕變與裂縫耦合作用,驗證模型能有效捕捉單雙預制缺陷巖石的裂縫演化規律及軸向應變特征,并在煤礦隧道工程中實現損傷區定量分析。
劉睿|周書偉|秦世康|張成凱|邢夢
同濟大學土木工程學院巖土工程系,中國上海200092
摘要
本研究提出了一種雙相場框架,用于模擬巖石蠕變斷裂過程,該框架考慮了拉伸和壓縮應變引起的斷裂,并探討了巖石蠕變與斷裂之間的耦合關系。通過雙相場模型,該框架能夠區分拉伸裂紋和壓縮-剪切裂紋,進而影響巖石的本構矩陣。此外,它還考慮了蠕變機制與斷裂機制之間的相互作用。實驗觀察結果證實,所提出的模型在三個方面與現有測試結果一致:蠕變行為、破壞模式以及蠕變引起的破壞。對含有單個和多個預制缺陷的試樣的研究顯示:蠕變會引發壓縮裂紋的形成;在相同的加載壓力下,蠕變速率不會改變裂紋的擴展方向;所提出的相場斷裂模型能夠有效捕捉巖石斷裂過程中的裂紋演變。進一步針對工程實踐的驗證表明,該模型在表征地下工程中圍巖損傷方面表現出了良好的性能。
引言
隨著地下空間的持續開發和利用,地質資源的利用變得越來越重要,例如壓縮空氣儲能(CAES)(Li等人,2023a;Zhou等人,2018a)、壓縮氫儲能(CHES)(Leng等人,2025;Liu等人,2025)、碳捕獲、利用與封存(CCUS)(Yu等人,2023a)以及增強型地熱系統(EGS)(Zhong等人,2023;Duwiquet等人,2021)等應用。在這些地下工程項目中,圍巖的斷裂會產生嚴重后果:在較輕的情況下,可能導致較大的洞室變形,需要維護和加固;在更嚴重的情況下,可能會導致隧道廢棄或地下空間坍塌,使其無法使用。因此,在地下設施的建設和運營過程中,準確預測并有效控制圍巖的斷裂至關重要。
與實驗測試相比,數值模擬提供了一種更為經濟的選擇。因此,已經開發了許多用于分析巖石斷裂和損傷的數值方法,如周期性動力學(PD)(Yu等人,2024;Ni等人,2021;Guo等人,2025)、邊界積分變分位移(BIVD)、有限體積法(BIVD蠕變巖石)、擴展有限元法(XFEM)(Zhang等人,2022a;Cheng等人,2024)、數值流形方法(NMM)(Yang等人,2020)、有限離散元法(FDEM)(Li等人,2025)、嵌入式離散斷裂模型(EDFM)(Kuang等人,2022)、離散元法(DEM)(Zhang等人,2022b)、粒子流代碼(PFC)(Wang等人,2023a;Wang等人,2024a)、物理信息神經網絡(PINN)(Goswami等人,2020;Samaniego等人,2020)以及真實失效過程分析(RFPA)(Wang等人,2012)。然而,這些方法在模擬裂紋擴展或圍巖損傷時需要額外的斷裂或損傷準則,并且在三維建模斷裂聚合和分支方面存在困難。相比之下,相場方法(PFM)(Navidtehrani等人,2021;Xu等人,2022a)作為一種高效的計算方法,近年來在巖石斷裂研究中得到了廣泛應用,因為它能夠準確捕捉裂紋的起始、擴展和失效過程。
巖石蠕變顯著影響地下工程的穩定性(Gu等人,2023;Gu等人,2026)。關于蠕變與斷裂建模的結合,Guo等人(2025)、Tian和Zhou(2024)以及Wang和Yin(2025)分別將PD與Bingham模型和分數階Bingham模型相結合,形成了蠕變PD模型;Wang和Cai(2021)將UDEC與Burgers模型集成;Wang等人(2023a;Wang等人,2024b)使用PFC計算巖石蠕變;Yu等人(2021)進行了NMM模擬。目前的巖石斷裂相場模型主要關注疲勞條件下的脆性斷裂和延性斷裂。然而,對于巖石的蠕變斷裂,現有的相場模型無法很好地預測這種類型的斷裂。鑒于許多隧道圍巖的破壞主要由蠕變效應引起,因此進一步發展相場斷裂理論并構建專門適用于圍巖蠕變斷裂的相場模型至關重要。
對于巖石的拉伸和剪切斷裂,裂紋的應力狀態存在顯著差異(Fei和Choo,2020;Guo等人,2025)。對于拉伸裂紋,法向應力消失,而剪切引起的裂紋則保持法向應力。Xu等人(2025)、Yu等人(2023b)、Wang等人(2024c)、Yue等人(2022)和He等人(2025)將驅動力分解為拉伸和剪切分量,以區分它們的不同效應。然而,在本構模型的模量退化過程中并未考慮這種區別。Wang等人(2023b)考慮了與不同斷裂機制相對應的本構退化,但忽略了法向應力的影響。Fei和Choo(2021)考慮了拉伸和剪切裂紋。值得注意的是,大多數當前的巖石相場斷裂模型無法同時考慮驅動力和本構建模中的不同斷裂機制。此外,采用多相場方法可以有效減輕不同斷裂驅動力之間的相互干擾(Fei和Choo,2021;Liu等人,2024),并在本構模型中區分拉伸和壓縮行為。為此,本研究提出了一種考慮本構模型中拉剪分解的蠕變斷裂相場模型。本研究采用雙相場框架,分別控制由拉伸和壓縮-剪切引起的本構斷裂行為。此外,所提出的模型還考慮了巖石蠕變過程中臨界能量釋放率的退化。
本文的結構如下:第2節回顧了相場建模和巖石蠕變的基本理論;第3節詳細闡述了巖石蠕變相場模型的構建和實現方法;第4節對所提出的模型進行了驗證;第5節研究了在不同條件下單個和雙重裂紋的蠕變行為,分析了裂紋擴展模式和軸向應變特性;第6節通過煤礦隧道驗證了該模型的適用性;最后,第7節對本文進行了總結。
節選內容
相場斷裂模型和巖石蠕變的基本理論
相場方法(PFM)基于Griffith的斷裂理論(Miehe等人,2010),其理論基礎來源于Francfort和Marigo(1998)提出的斷裂變分公式。此外,分數階本構模型(Miller和Ross,1993)能夠有效描述巖石蠕變的初始階段,并表現出良好的適應性(Kilbas,2006)。因此,本節將主要討論脆性PFM和蠕變本構模型
用于巖石蠕變斷裂的改進型PFM
第2節介紹了相場斷裂模型的基本理論方面。然而,該模型是一個彈性-拉伸模型,沒有考慮壓縮應力或壓縮的影響。在地質工程中,巖石斷裂不僅由拉伸引起,還由壓縮和剪切引起。因此,必須考慮壓縮和剪切引起的斷裂。這將通過引入壓縮應變的影響來實現。此外,目前的巖石相場模型
模型驗證
為了驗證所提出的巖石蠕變相場模型的可靠性,本節測試了三個數值案例。第一個案例是未存在預先裂紋的巖石的單軸壓縮蠕變。第二個案例是存在預先裂紋的巖石的單軸壓縮試驗,用于驗證斷裂形態。第三個案例是考慮蠕變效應的具有預先裂紋的單軸壓縮試驗。這三個案例分別進行
實驗室規模的數值示例
本節進一步展示了所提出的巖石蠕變相場模型在實驗室規模的應用。分別研究了具有單個和多個預制缺陷的巖石試樣的裂紋擴展行為。對于單缺陷試樣,考察了不同加載水平和蠕變速率對裂紋發展特性的影響。對于雙缺陷試樣,研究了不同加載水平和缺陷中心距離對裂紋擴展的影響
煤礦隧道中的工程應用
前幾節驗證了所提出的PFM在具有單個和多個缺陷的試樣中用于蠕變裂紋起始、擴展和橋接的效果,并分析了軸向應變曲線的時間演變。為了進一步證明該模型在實際工程中的實用性,本研究將所提出的模型擴展到隧道圍巖的灌漿中,該模型量化了損傷區域以優化灌漿孔深度。
結論
提出了一種考慮巖石蠕變效應的相場斷裂模型。該模型通過雙相場斷裂框架建立,考慮了拉伸和壓縮應變引起的斷裂對本構關系和巖石累積損傷特性的不同影響。主要結論如下:
(1)與傳統的相場斷裂模型相比,所提出的模型更準確地描述了裂紋的起始和擴展
CRediT作者貢獻聲明
劉睿:撰寫——原始草稿。周書偉:監督、資源支持。秦世康:軟件開發。張成凱:模型驗證。邢夢:實驗研究。
利益沖突聲明
作者聲明他們沒有已知的財務利益或個人關系可能影響本文報告的工作。
致謝
本工作得到了國家自然科學基金(52278408)和中央高;究蒲袠I務費(22120250186)的支持。
生物通微信公眾號
生物通新浪微博
今日動態 |
人才市場 |
新技術專欄 |
中國科學人 |
云展臺 |
BioHot |
云講堂直播 |
會展中心 |
特價專欄 |
技術快訊 |
免費試用
版權所有 生物通
Copyright© eBiotrade.com, All Rights Reserved
聯系信箱:
粵ICP備09063491號