GitLab

シェル要素の層別結果出力において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