Sources Esope | Cast3M

vfsgp2
C VFSGP2    SOURCE    CB215821  16/04/21    21:18:39     8920
 
      SUBROUTINE VFLSIGPOWDER2(SIG,VECM)
      IMPLICIT INTEGER(I-N)
      integer i
      real*8  sig(3),vecm(3),aux1,aux2,aux3,xi1,xj2,aux4
      real*8               nn1,nn2,eta0,aa1,aa2,bb1,bb2,
     .                     pia,pib,sigy0,pcc,rrr
      common /poder2_data/ nn1,nn2,eta0,aa1,aa2,bb1,bb2,
     .                     pia,pib,sigy0,pcc,rrr
      call Invari_I1(sig,3,xi1)
      call desviador(sig,vecm,3)
      if (xi1.lt.pia) then
         aux2 = 2.D0*(xi1-pia)/rrr**2
         do i=1,3
           vecm(i)= vecm(i)+aux2
         enddo
      else
         aux1 = (pia-xi1)/(pia-pcc)
         if (aux1.gt.1.D0) aux1=1.D0
         aux2 = bb1/3.D0
         call Invari_J2(sig,3,xj2)
         aux3 = 1.D0 / SQRT(2.D0*xj2)
         aux4 = xj2/(aa2*sigy0**2/3.D0)
         if (aux4.gt.1.D0) aux4=1.D0
****** D-P Meschke examples
*         aux1 = 1.D0
*         aux4 = 1.D0
****** D-P Meschke examples
         do i=1,3
           vecm(i)= aux4*aux3*vecm(i)+aux2*aux1
         enddo
      endif
      return
      end