Télécharger la plaquette Cast3M


fr en

Télécharger eco_j2.eso

Retour à la liste

Numérotation des lignes : 

eco_j2

  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
  27.       call carac_mate_j2(xmat)
  28.       do i=1,ndims
  29.         sigtri(i)=sigini(i)+dsigtr(i)
  30.       enddo
  31.       call yielddJ2(sigtri,ndims,varini,ione,resu)
  32.       if (resu.le.0.D0) then
  33.          if (nescri.eq.1) write(ues,*)' Paso elastico'
  34.          iplast=0
  35.          do i=1,ndims
  36.            sigfin(i)=sigtri(i)
  37.          enddo
  38.          varfin(1)=varini(1)
  39.          varfin(2)=-1.D0
  40.          iterlocal=0
  41.          return
  42.       endif
  43.       if (nescri.eq.1) write(ues,*)' Paso plastico'
  44.       do i=1,ndims
  45.          xtri(i)=sigtri(i)
  46.          xini(i)=sigtri(i)
  47.       enddo
  48.       lini=0.D0
  49.       if ((kiso.eq.0.D0).and.(velo.eq.0.D0)) then
  50.          if (nescri.eq.1) write(ues,*)' Plasticidad sin endurecimiento'
  51.          call Integra0J2(xtri,xini,ndims,lini,ddefpl,ndims,
  52.      .                   precis,nitmax,nescri,ues,kerre,iterlocal)
  53.          varfin(1)=varini(1)+lini
  54.          varfin(2)=lini
  55.       else
  56.          if (nescri.eq.1) write(ues,*)' Plasticidad con endurecimiento'
  57.          ndimx=ndims+1
  58.          xtri(ndimx)=varini(1)
  59.          xini(ndimx)=varini(1)
  60.          call Integra1J2(xtri,xini,ndimx,ione,lini,ddefpl,ndims,
  61.      .                   precis,nitmax,nescri,ues,kerre,iterlocal)
  62.          varfin(1)=xini(ndimx)
  63.          varfin(2)=lini
  64.       endif
  65.         do i=1,ndims
  66.           sigfin(i)=xini(i)
  67.         enddo
  68.       return
  69.       end
  70.  
  71.  
  72.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales