此子程序用于用户自定义初始应力,由*INITIAL CONDITIONS,TYPE=STRESS卡片上的USER参数表征。头部和变量说明如下:
subroutine sigini(sigma,coords,ntens,ncrds,noel,npt,layer,
& kspt,lrebar,rebarn)
!
! 用户子程序 sigini
!
! 输入:
!
! coords 积分点的坐标
! ntens 要定义的应力分量数目
! ncrds 坐标数目
! noel 单元编号
! npt 积分点编号
! layer 当前未使用
! kspt 当前未使用
! lrebar 当前未使用 (值: 0)
! rebarn 当前未使用
!
! 输出:
!
! sigma(1..ntens) 积分点处的初始应力值。如果ntens=6,
! 分量顺序为11,22,33,12,13,23
!