大哥们，sigini随机区域如何施加残余应力？ 100

Subroutine sigini(sigma, coords, ntens, ncrds, noel, npt, layer, kspt, lrebar, names)

!

  Include 'ABA_PARAM.INC'

!

  Dimension sigma(ntens), coords(ncrds)

  Character *80 names(*)

!

  Real :: x_1, y_1, z_1, wi

  Integer :: stripnum, ni, lo, ui, vi

  wi = 0

  Call random_seed()

  Call random_number(x_1)

  wi=100*x_1

  ni = int(wi)

  If (coords(2)>0 .And. coords(2)<ni) Then

    If (coords(1)>0 .And. coords(1)<200) Then

      If (layer==3 .And. kspt==1) Then

        sigma(1) = 0

      End If

      If (layer==3 .And. kspt==2) Then

        sigma(1) = -100000000

      End If

      If (layer==3 .And. kspt==3) Then

        sigma(1) = -300000000

      End If

      If (layer==3 .And. kspt==4) Then

        sigma(1) = -900000000

      End If

      If (layer==3 .And. kspt==5) Then

        sigma(1) = -1200000000

      End If

      If (layer==3 .And. kspt==6) Then

        sigma(1) = -1100000000

      End If

      If (layer==3 .And. kspt==7) Then

        sigma(1) = -1000000000

      End If

    End If

  End If

  Return

End Subroutine sigini