初始应力场 (sigini.f)

此子程序用于用户自定义初始应力,由*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
!