シェル要素の層別結果出力においてmises応力が出力されない
シェル要素の層別結果出力においてmises応力が出力されない
fstr_NodalStress6Dでlayer配下のmisesが計算されていないため
!C** calculate von MISES stress
do i = 1, hecMESH%n_node
fstrSOLID%MISES(i) = get_mises(fstrSOLID%STRESS(6*(i-1)+1:6*(i-1)+6))
enddo
do i = 1, hecMESH%n_elem
fstrSOLID%EMISES(i) = get_mises(fstrSOLID%ESTRESS(6*(i-1)+1:6*(i-1)+6))
enddo
deallocate( nnumber )
end subroutine fstr_NodalStress6D