Esope source files

   1 : C ECO_J2    SOURCE    CHAT      05/01/12    23:25:03     5004           
   2 : CCC
   3 : C **********************************************************************
   4 : CCC
   5 :       SUBROUTINE ECO_J2 (SIGINI,VARINI,DSIGTR,SIGFIN,VARFIN,
   6 :      .                   ddefpl,iplast,ndims,xmat,kerre,
   7 :      .                   precis,nitmax,nescri,ues,iterlocal)
   8 :       IMPLICIT INTEGER(I-N)
   9 :       integer ndims,kerre,nitmax,nescri,ues,iterlocal,ione
  10 :       real*8  sigini(ndims),varini(2),ddefpl(ndims),dsigtr(ndims),
  11 :      .        sigfin(ndims),varfin(2),xmat(*),precis
  12 :       integer iplast,ndimx
  13 :       real*8  sigtri(6),xtri(7),xini(7),lini,void,resu
  14 :       real*8            sigy0,kiso,siginf,velo
  15 :       common /j2j2data/ sigy0,kiso,siginf,velo
  16 :       void  =0.D0
  17 :       ione  =1
  18 :       kerre =0
  19 :       iplast=1
  20 :       do i=1,6
  21 :         sigtri(i)=0.D0
  22 :       enddo
  23 :       do i=1,7
  24 :         xtri(i)=0.D0
  25 :         xini(i)=0.D0
  26 :       enddo
  26 :       call carac_mate_j2(xmat)
  27 :       do i=1,ndims
  28 :         sigtri(i)=sigini(i)+dsigtr(i)
  29 :       enddo
  29 :       call yielddJ2(sigtri,ndims,varini,ione,resu)
  30 :       if (resu.le.0.D0) then
  31 :          if (nescri.eq.1) write(ues,*)' Paso elastico'
  32 :          iplast=0
  33 :          do i=1,ndims
  34 :            sigfin(i)=sigtri(i)
  35 :          enddo
  36 :          varfin(1)=varini(1)
  37 :          varfin(2)=-1.D0
  38 :          iterlocal=0
  39 :          return      
  40 :       endif
  41 :       if (nescri.eq.1) write(ues,*)' Paso plastico'
  42 :       do i=1,ndims
  43 :          xtri(i)=sigtri(i)
  44 :          xini(i)=sigtri(i)
  45 :       enddo
  46 :       lini=0.D0
  47 :       if ((kiso.eq.0.D0).and.(velo.eq.0.D0)) then
  48 :          if (nescri.eq.1) write(ues,*)' Plasticidad sin endurecimiento'
  48 :          call Integra0J2(xtri,xini,ndims,lini,ddefpl,ndims,
  49 :      .                   precis,nitmax,nescri,ues,kerre,iterlocal)
  50 :          varfin(1)=varini(1)+lini
  51 :          varfin(2)=lini
  52 :       else
  53 :          if (nescri.eq.1) write(ues,*)' Plasticidad con endurecimiento'
  54 :          ndimx=ndims+1
  55 :          xtri(ndimx)=varini(1)
  56 :          xini(ndimx)=varini(1)
  56 :          call Integra1J2(xtri,xini,ndimx,ione,lini,ddefpl,ndims,
  57 :      .                   precis,nitmax,nescri,ues,kerre,iterlocal)
  58 :          varfin(1)=xini(ndimx)
  59 :          varfin(2)=lini
  60 :       endif
  61 :         do i=1,ndims
  62 :           sigfin(i)=xini(i)
  63 :         enddo
  64 :       return
  65 :       end
  66 :  
  67 :