巖土工程大變形問題（如滑坡、泥石流、靜力觸探試驗、打樁等）是巖土工程理論和實踐中一直以來的難點問題。近日，河海大學(xué)力學(xué)與材料學(xué)院袁維海副研究員以“Particle finite element method implementation for large deformation analysis using Abaqus”為題在巖土工程領(lǐng)域知名期刊Acta Geotechnica發(fā)表最新研究成果。該研究通過對Abaqus進行二次開發(fā)實現(xiàn)了粒子有限元法（PFEM）這一大變形數(shù)值分析手段，彌補了傳統(tǒng)商業(yè)軟件中大變形分析功能的不足，為巖土工程大變形分析提供了簡易的計算工具，拓展了巖土工程大變形分析的新途徑，有利于PFEM方法的推廣應(yīng)用。

圖1 粒子有限元法（PFEM）計算流程

粒子有限元法將有限元網(wǎng)格中的節(jié)點視為自由粒子，根據(jù)分析步的位移增量更新粒子的位置，并通過周期性的計算域重新識別和網(wǎng)格重新剖分，解決傳統(tǒng)有限元法存在的網(wǎng)格扭曲問題，能有效的實現(xiàn)巖土大變形數(shù)值模擬。研究采用Abaqus的內(nèi)置函數(shù)進行常規(guī)增量FEM分析，并使用mesh-to-mesh mapping方法實現(xiàn)場變量的映射，使用單個Python腳本就能實現(xiàn)PFEM數(shù)值模擬，顯著減少了編碼工作量，為PFEM在巖土工程中的推廣應(yīng)用奠定了基礎(chǔ)。

圖2 粒子有限元法（PFEM）在Abaqus中實現(xiàn)

研究通過分析彈性懸臂梁彎曲，剛性基礎(chǔ)貫入試驗、T-bar貫入試驗以及海床-管線相互作用等問題，驗證了該方法的準(zhǔn)確性與穩(wěn)定性。結(jié)果表明，研究所提出的Abaqus-PFEM方法能夠很好解決巖土大變形及土-結(jié)構(gòu)相互作用問題。

Abaqus-PFEM方法結(jié)合了PFEM處理大變形問題的能力及Abaqus軟件強大的有限元求解和前后處理功能。該方法具有計算精度和計算效率高、應(yīng)力結(jié)果光滑、處理接觸非線性和材料非線性簡潔等優(yōu)點，為解決巖土工程大變形問題提供了新的思路和方法。

圖3 彈性懸臂梁彎曲模擬結(jié)果

圖4 剛性基礎(chǔ)貫入模擬結(jié)果

圖5 T-bar貫入模擬結(jié)果

圖6 海床-管線相互作用模擬結(jié)果

論文研究工作由國家重點研發(fā)計劃項目和國家自然科學(xué)基金項目等提供經(jīng)費支持。 該論文于2021年7月入選ESI高被引論文。

論文鏈接：

Yuan W.-H., Wang H.-C., Zhang W., Dai B.-B., Liu K., Wang Y., Particle finite element method implementation for large deformation analysis using Abaqus. Acta Geotechnica, 2021, 16(8): 2449–2462.

https://link.springer.com/article/10.1007/s11440-020-01124-2

END

素材來源：河海大學(xué)官網(wǎng)、力材院

編輯：小海螺記者 郭一帆；張然、柏蓉

推薦閱讀

第八次祭奠 | 河海大學(xué)舉行南京大屠殺死難者國家公祭日悼念活動

銘記歷史 | 河海大學(xué)的校園里，有一座侵華日軍南京大屠殺清涼山遇難同胞紀(jì)念碑

江蘇省委宣講團在河海大學(xué)宣講黨的十九屆六中全會精神

河海大學(xué)思想政治工作側(cè)記 | 堅持立德樹人根本 為黨育人為國育 才

重要提醒！河海大學(xué)2022年全國碩士研究生招生考試自命題科目答題紙條碼粘貼要求