
注冊
/
登錄使用Abaqus的VUMAT子程序做內聚力單元的本構,abaqus會傳入一個strainInc(nblock,ndir + nshr )應變增量,我輸出了ndir=1,nshr=2,剛好第一個是法向應變增量,后面兩個是切向應變增量,然后切向的是真實剪應變,需要乘2倍后換成工程剪應變再使用。
在VUMAT中更新應力是用stressNew=stressOld+D*strainInc,其中D是剛度矩陣。我不理解的是,這個應變增量strainInc是程序通過什么方法計算出來的,特別是在分析步剛開始的時候?







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






自己編寫材料vumat,利用平面應力單元模擬,但通過輸出應變增量發現strainInc(k,3)在過程中始終等于零,所以最后3方向的應變也始終為零。但是利用abaqus自帶材料模型模擬平面應力狀態會發現3的應變不為零。自己的子程序用三維單元驗證是沒問題。請問大佬這是哪里出問題?
跳至頁
TOP