基于Abaqus的修正劍橋模型的vumat子程序開發(fā)

巖土工程中廣泛采用的莫爾–庫侖準(zhǔn)則,能較好地描述巖土材料的強(qiáng)度特性和破壞行為,而涉及變形和固結(jié)分析問題時(shí)該模型從理論和實(shí)踐中均有一定的不足,基于等向加載試驗(yàn)建立起來的劍橋模型則能較好地描述黏土的彈塑性變形特性。

傳統(tǒng)相關(guān)聯(lián)修正劍橋本構(gòu)模型的屈服面方程、硬化準(zhǔn)則和流動(dòng)法則分別為

基于Abaqus的修正劍橋模型的vumat子程序開發(fā)的圖1

由于式(1)所確定的屈服軌跡在p平面是一個(gè) 圓,不能反映巖土介質(zhì)拉壓不等(S-D)效應(yīng),而且劍 橋模型是基于正常固結(jié)狀態(tài)試驗(yàn)推導(dǎo)而來,子午面上臨界狀態(tài)線通過應(yīng)力坐標(biāo)原點(diǎn),表現(xiàn)為不考慮土 體黏聚力的純摩擦型本構(gòu),而大多數(shù)的巖土介質(zhì)具 有一定的黏聚力,屬摩擦–黏聚型材料。考慮到莫 爾–庫侖準(zhǔn)則有此特征,故將莫爾–庫侖準(zhǔn)則與劍橋模型相結(jié)合,得

基于Abaqus的修正劍橋模型的vumat子程序開發(fā)的圖2

聯(lián)合式(1),(4),(5),有改進(jìn)的屈服函數(shù):

基于Abaqus的修正劍橋模型的vumat子程序開發(fā)的圖3

根據(jù)上述理論,可以通過以下流程圖完成VUMAT子程序編寫

基于Abaqus的修正劍橋模型的vumat子程序開發(fā)的圖4

通過單胞模型計(jì)算得到得結(jié)果如下圖所示。

基于Abaqus的修正劍橋模型的vumat子程序開發(fā)的圖5

[1]袁克闊,陳衛(wèi)忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應(yīng)的修正劍橋模型及數(shù)值實(shí)現(xiàn)[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(08):1574-1579.

最后,有需要?dú)g迎通過微信公眾號(hào)聯(lián)系我們。

微信公眾號(hào):320科技工作室。

登錄后免費(fèi)查看全文
立即登錄
App下載
技術(shù)鄰APP
工程師必備
  • 項(xiàng)目客服
  • 培訓(xùn)客服
  • 平臺(tái)客服

TOP

22
8
16