例子如下：

*INCLUDE,INPUT=initial.ipe500

*INITIAL CONDITIONS,TYPE=STRESS,USER

*USER SUBROUTINES

SUBROUTINE SIGINI(SIGMA,COORDS,NTENS,NCRDS,NOEL,NP,LA,M,LR,RE)

INCLUDE ‘ABA_PARAM.INC’

DIMENSION SIGMA(NTENS),COORDS(NCRDS)

CHARACTER*8 RE

IF(NOEL.GE.1.AND.NOEL.LE.100)THEN

IF(M.EQ.1.OR.M.EQ.17.OR.M.EQ.33.OR.M.EQ.49)THEN

SIGMA(1)=-7.050

ELSE IF(M.EQ.2.OR.M.EQ.16.OR.M.EQ.34.OR.M.EQ.48)THEN

SIGMA(1)=-5.2875

ELSE IF(M.EQ.3.OR.M.EQ.15.OR.M.EQ.35.OR.M.EQ.47)THEN

SIGMA(1)=-3.525

ELSE IF(M.EQ.4.OR.M.EQ.14.OR.M.EQ.36.OR.M.EQ.46)THEN

SIGMA(1)=-1.7625

ELSE IF(M.EQ.5.OR.M.EQ.13.OR.M.EQ.37.OR.M.EQ.45)THEN

SIGMA(1)=0.

ELSE IF(M.EQ.6.OR.M.EQ.12.OR.M.EQ.38.OR.M.EQ.44)THEN

SIGMA(1)=1.7625

ELSE IF(M.EQ.7.OR.M.EQ.11.OR.M.EQ.39.OR.M.EQ.43)THEN

SIGMA(1)=3.525

ELSE IF(M.EQ.8.OR.M.EQ.10.OR.M.EQ.40.OR.M.EQ.42)THEN

SIGMA(1)=5.2875

ELSE IF(M.EQ.9.OR.M.EQ.41)THEN

SIGMA(1)=7.050

ELSE IF(M.EQ.18.OR.M.EQ.32)THEN

SIGMA(1)=5.3439

ELSE IF(M.EQ.19.OR.M.EQ.31)THEN

SIGMA(1)=3.6378

ELSE IF(M.EQ.20.OR.M.EQ.30)THEN

SIGMA(1)=1.9317

ELSE IF(M.EQ.21.OR.M.EQ.29)THEN

SIGMA(1)=0.2256

ELSE IF(M.EQ.22.OR.M.EQ.28)THEN

SIGMA(1)=-1.4805

ELSE IF(M.EQ.23.OR.M.EQ.27)THEN

SIGMA(1)=-3.1866

ELSE IF(M.EQ.24.OR.M.EQ.26)THEN

SIGMA(1)=-4.8927

ELSE IF(M.EQ.25)THEN

SIGMA(1)=-6.5988

ENDIF

ENDIF

RETURN

END

在分析步之前添加*INITIAL CONDITIONS,TYPE=STRESS,USER

楼主的问题解决了嘛