Esope source files

   1 : C ECORHM    SOURCE    CHAT      05/01/12    23:25:17     5004           
   2 : CCC
   3 : C **********************************************************************
   4 : CCC
   5 :       SUBROUTINE ECO_RHMC (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,ione,
  10 :      .          iplast,i,iterlocal
  11 :       real*8    sigini(    6),varini(2),ddefpl(    6),sigtr(    6),
  12 :      .          sigfin(    6),varfin(2),xmat(*),precis,sigtri(6),
  13 :      .          xtri(7),xini(7),lini,void,resu,dsigtr(    6)
  14 :       void  =0.D0
  15 :       ione  =1
  16 :       kerre =0
  17 :       iplast=1
  18 :       do i=1,7
  19 :         xtri(i)=0.D0
  20 :         xini(i)=0.D0
  21 :       enddo
  21 :       call carac_mate_rhmc(xmat)
  22 :       do i=1,ndims
  23 :         sigtri(i)=sigini(i)+dsigtr(i)
  24 :       enddo
  24 :       call yielddRHMC(sigtri,ndims,void,ione,resu)
  25 :       if (resu.le.0.D0) then
  26 :          if (nescri.eq.1) write(ues,*)' Paso elastico'
  27 :          iplast=0
  28 :          do i=1,ndims
  29 :            sigfin(i)=sigtri(i)
  30 :          enddo
  31 :          varfin(1)=varini(1)
  32 :          varfin(2)=-1.D0
  33 :          iterlocal=0
  34 :          return
  35 :       endif
  36 :       if (nescri.eq.1) write(ues,*)' Paso plastico'
  37 :       do i=1,ndims
  38 :          xtri(i)=sigtri(i)
  39 :          xini(i)=sigtri(i)
  40 :       enddo
  41 :       lini=0.D0
  41 :       call Integra0RHMC(xtri,xini,ndims,lini,ddefpl,ndims,
  42 :      .                  precis,nitmax,nescri,ues,kerre,iterlocal)
  43 :       varfin(1)=varini(1)+lini
  44 :       varfin(2)=lini
  45 :         do i=1,ndims
  46 :           sigfin(i)=xini(i)
  47 :         enddo
  48 :       return
  49 :       end
  50 :  
  51 :  
  52 :