Esope source files

   1 : C MTC_J2    SOURCE    CHAT      05/01/13    01:53:41     5004           
   2 : CCC
   3 : C **********************************************************************
   4 : CCC
   5 :       SUBROUTINE MTC_J2 (KMAT,NDIMK,SIG,NDIMS,VINT,LAM,XMAT,
   6 :      .                   nescri,ues,kerre)
   7 :       IMPLICIT INTEGER(I-N)
   8 :       integer ndims,kerre,nescri,ues,ndimk,ifl,ndimx,ndimv
   9 :       real*8    sig(ndims),vint,Kmat(*),xmat(*),x(7),void,lam
  10 :       real*8            sigy0,kiso,siginf,velo
  11 :       common /j2j2data/ sigy0,kiso,siginf,velo
  12 :       kerre=0
  13 :       do i=1,7
  14 :         x(i)=0.D0
  15 :       enddo
  15 :       call carac_mate_j2(xmat)
  16 :       if (lam.le.0.D0) then
  17 :          if (nescri.eq.1) write(ues,*)' Paso elastico'
  18 :          ifl=1
  18 :          call MatHok(Kmat,ndimk,ndims,ifl)
  19 :          return      
  20 :       endif
  21 :       if (nescri.eq.1) write(ues,*)' Paso plastico'
  22 :       if ((kiso.eq.0.D0).and.(velo.eq.0.D0)) then
  23 :          if (nescri.eq.1) write(ues,*)' Plasticidad sin endurecimiento'
  24 :          ndimv=0
  24 :          call MTC0J2(Kmat,ndimk,sig,ndims,ndimv,lam,
  25 :      .               nescri,ues,kerre)
  26 :          return
  27 :       else
  28 :          if (nescri.eq.1) write(ues,*)' Plasticidad con endurecimiento'
  29 :          ndimv=1
  30 :          ndimx=ndims+1
  31 :          do i=1,ndims
  32 :               x(i)=sig(i)
  33 :          enddo
  34 :          x(ndimx)=vint
  35 :          void=0.D0
  35 :          call MTC1J2(Kmat,ndimk,x,ndimx,ndimv,lam,ndims,
  36 :      .               nescri,ues,kerre)
  37 :       endif
  38 :       return
  39 :       end
  40 :  
  41 :