
注冊
/
登錄
在編寫子程序的時候遇到了點問題…想用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) 剛開始接觸子程序,
), coordMp(nblock,*), 1 direct(nblock,3,3), T(nblock,3,3), 2 charLength(nblock), props(nprops), 3 stateOld(nblock,nstatev), 4 stateNew(nblock




But when using single precision, these values are correct nblock: 144, props: 5769 3846.



extern "C" void vumat( long *nblock, long *ndir, long *nshr, long *nstatev, long *nfieldv, long *nprops, long *lanneal, double *stepTime, double *totalTime, double *dt, char* cmname, double
使用Abaqus的VUMAT子程序做內聚力單元的本構,abaqus會傳入一個strainInc(nblock,ndir + nshr )應變增量,我輸出了ndir=1,nshr=2,剛好第一個是法向應變增量,后面兩個是切向應變增量,然后切向的是真實剪應變,需要乘2倍后換成工程剪應變再使用。
跳至頁
TOP