call vgetvrm('PE',strain,jData,cData,jStatus)?
e11=strain(k)
e22=strain(nblock+k)
e33=strain(nblock*2+k)
e12=strain(nblock*3+k)
e23=strain(nblock*4+k)
e13=strain(nblock*5+k)
stateNew(k,2)=e11
stateNew(k,3)=e22
stateNew(k,4)=e33
stateNew(k,5)=e12
stateNew(k,6)=e23
stateNew(k,7)=e13
do r=1,6
STR(k,r)=strain(k+(r-1)*nblock)
end do
li913 发表于 2022-9-12 17:42
do r=1,6
STR(k,r)=strain(k+(r-1)*nblock)
end do
欢迎光临 Fortran Coder (http://bbs.fcode.cn/) | Powered by Discuz! X3.2 |