
注冊
/
登錄



在編寫子程序的時候遇到了點問題…想用vsprinc來提取最大主應力的…但是得到的結果很怪 dimension s(nblock, ndir+nshr), eigVal(nblock, 3) call vsprinc(nblock, s, eigVal, ndir, nshr) do k =1, nblock stateNew(k,2)=eigVal(k,1) 剛開始接觸子程序,





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






跳至頁
TOP