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

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












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



跳至頁
TOP