Sources Esope | Cast3M

delskm
C DELSKM    SOURCE    CB215821  16/04/21    21:16:22     8920
 
      SUBROUTINE DETERMINA_LS_KMA(SIG,NESCRI,KMA,IZONE)
      IMPLICIT INTEGER(I-N)
      integer kma,nescri,izone
      real*8  sig(*),xi1,xj2,valor
      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_J2(sig,3,xj2)
      call Invari_I1(sig,3,xi1)
      if (xi1.lt.pia) then
         kma=2
         valor=rrr**2*xj2+(xi1-pia)**2-(pia-pib)**2
         valor=xj2+((xi1-pia)/rrr)**2-aa2*sigy0**2/3.D0
         if (nescri.eq.1)
     .   write(199,*)' CAP ',xi1,xj2,
     .   pia,pib,rrr,((xi1-pia)/rrr)**2,aa2,valor
         izone = 0
      else
         kma=3
         valor=SQRT(2.D0*xj2)+bb1*xi1/3.D0-bb2
         if (nescri.eq.1)
     .   write(199,*)' CONE',xi1,pia,valor,kma
         izone = 1
      endif
      return
      end