近年來(lái)，巖土工程大變形數值仿真分析方法迅猛發(fā)展。該方法既能模擬評估巖土工程的穩定性，為實(shí)際工程提供設計支撐的基礎，同時(shí)能進(jìn)一步模擬工程結構破壞形式、失效過(guò)程等，評價(jià)工程失穩所造成的危害。為推動(dòng)大變形數值仿真分析方法在實(shí)際工程的應用，中國科學(xué)院武漢巖土力學(xué)研究所王頔助理研究員、陳昊博士后及其團隊成員開(kāi)發(fā)了一套自主化國產(chǎn)軟件，該軟件基于物質(zhì)點(diǎn)法（MPM）、離散單元法（DEM）和傳統有限元法（FEM），結合圖形加速處理技術(shù)（GPU），構建了高效的巖土工程大變形數值仿真分析平臺，實(shí)現了諸如新磨村滑坡等大變形失效過(guò)程的高效模擬（如圖1）。

　　目前，巖土工程大變形數值模擬方法仍存在一些缺陷，特別是對巖土體失效機制的模擬和分析方面。已有數值模擬方法多注重于宏觀(guān)評估巖土工程問(wèn)題的穩定性、破壞形式、失效過(guò)程等，通常采用唯象本構模型和連續性介質(zhì)假設對巖土體進(jìn)行模擬。這些方法雖然能把握巖土體的物理力學(xué)響應，但缺乏對巖土體內部作用機制的捕捉，難以充分刻畫(huà)內在失效機制；另一方面，基于離散介質(zhì)的數值模擬方法能較好地描述土體顆粒層面的微觀(guān)機制，但在處理實(shí)際工程問(wèn)題時(shí)受到計算成本的限制。因此，工程尺度邊坡失效破壞過(guò)程仍缺乏有效的模擬方法和深入的機制研究。

　　針對這一現狀，該套自主化國產(chǎn)軟件新增了基于MPM-DEM耦合的多尺度數值分析模塊，開(kāi)發(fā)了一種基于MPM-DEM耦合的層級多尺度數值模擬方法（如圖2），并將其應用于三軸實(shí)驗模擬、邊坡破壞分析等方面，旨在全面了解邊坡失效和破壞后的巖土體宏微觀(guān)響應機制。一方面，該方法在傳統物質(zhì)點(diǎn)法（MPM）的基礎上引入Affine-particle-in-cell (APIC)速度更新格式和B樣條基函數，提高了宏觀(guān)模擬的穩定性；另一方面，引入了代表性體積單元（RVE）替代本構模型，捕捉巖土體微觀(guān)物理力學(xué)演化機制。通過(guò)在顆粒間使用Coulomb摩擦定律和Johnson-Kendall-Robert（JKR）粘聚力模型再現微觀(guān)顆粒的摩擦和粘結效應。

　　基于該層級多尺度數值模擬方法，模擬再現了不同摩擦、粘結、級配特性下邊坡的大變形失效過(guò)程，探究了邊坡大變形失效過(guò)程中滑帶內部土體的微觀(guān)失效機制。多尺度分析顯示：滑帶破壞中，坡趾處的土體受到剪切作用強烈，而頂部破壞作用則受到拉伸作用（如圖3）；基于該方法對漸進(jìn)式邊坡破壞過(guò)程進(jìn)行了更深入的分析，發(fā)現這類(lèi)失效往往由坡趾處微小的破壞引起，但可以持續發(fā)展至數公里的范圍。層級多尺度方法模擬分析結果顯示：土壤的流動(dòng)破壞與顆粒材料的剪切破壞相關(guān)，而宏觀(guān)剪切帶的發(fā)展通常伴隨著(zhù)拉伸破壞（如圖4）。由于剪切和拉伸破壞是摩擦和粘聚材料的典型破壞機制，可以推斷摩擦和粘聚效應在不同的漸進(jìn)式破壞模式中起到關(guān)鍵作用。因此，通過(guò)進(jìn)一步探究微觀(guān)顆粒粘結和摩擦作用力的分布規律，提出了一個(gè)微觀(guān)描述參數 來(lái)量化摩擦和粘聚之間的相互作用。進(jìn)一步的分析證明： 可以捕捉摩擦和粘聚的效應，并區分不同的退步性破壞模式（如圖5）。最后，建立了一個(gè)顆粒邊坡的一系列漸進(jìn)式破壞模式云圖，其中破壞機制通過(guò)不同的主導效應來(lái)解釋?zhuān)赐寥乐械哪Σ粱蛘尘郏ㄈ鐖D6）。

　　相關(guān)研究成果發(fā)表在Journal of Geophysical Research: Solid Earth和Computers and Geotechnics等國際期刊上，研究工作獲得了中國科學(xué)院青年人才項目、國家自然科學(xué)基金面上項目（51979270、51709258）資助。

　　論文鏈接：

　　https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2022JB026008

　　https://www.sciencedirect.com/science/article/abs/pii/S0266352X22002300?via%3Dihub　

圖1 新磨村滑坡堆積模擬結果與勘探及其他模擬方法對比　　  

圖2 基于物質(zhì)點(diǎn)法-離散單元法耦合的多尺度數值模擬方法　　 

圖3 邊坡失效滑帶多尺度分析　　 

圖4 漸進(jìn)式邊坡破壞機制多尺度分析　　 

圖5 破壞模式量化指標 與分布模式

圖6 邊坡漸進(jìn)式破壞模式分布云圖