
注冊
/
登錄

使用Abaqus的VUMAT子程序做內(nèi)聚力單元的本構(gòu),abaqus會傳入一個strainInc(nblock,ndir + nshr )應變增量,我輸出了ndir=1,nshr=2,剛好第一個是法向應變增量,后面兩個是切向應變增量,然后切向的是真實剪應變,需要乘2倍后換成工程剪應變再使用。


在VUMAT中更新應力是用stressNew=stressOld+D*strainInc,其中D是剛度矩陣。我不理解的是,這個應變增量strainInc是程序通過什么方法計算出來的,特別是在分析步剛開始的時候?

請問一下,如果需要在VUMAT子程序構(gòu)建滯回曲線,如下圖所示。那么該如何判斷加卸載的轉(zhuǎn)折點1以及轉(zhuǎn)折點2?我看別人實現(xiàn)的時候是用應變增量小于0時為卸載,應變增量大于0時為加載。但是在計算過程中發(fā)現(xiàn),顯式分析的應變會存在震蕩現(xiàn)象,就算是純加載階段,應變增量也會頻繁出現(xiàn)負值,好像應變的增減是完全隨機的。即使添加了瑞利阻尼也只是減小了這種震蕩現(xiàn)象,不能消除。




本人水利碩士,使用Abaqus對土石壩進行模擬,經(jīng)常遇到不收斂的問題,想詳細學習MSG文件,來知道模型不收斂的原因,但苦于找不到詳細的MSG文件說明教程,還請看到的大佬賜教,感激不盡。





有無做過的大佬分享一下子程序嗎???
跳至頁
TOP