strain(nblock*6)这个数组的大小是这个,我想讲其变为矩阵STR(nblock,6),不知道这样的操作是否正确
[Fortran] 纯文本查看 复制代码 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
|