Télécharger la plaquette Cast3M


fr en

Télécharger kmkf3d.eso

Retour à la liste

Numérotation des lignes : 

kmkf3d

  1. C KMKF3D    SOURCE    FD218221  24/02/07    21:15:17     11834          
  2.         subroutine kmkf3d(xmat,nmat,vnmat,nstype,nbelas,sim3,kmecap)
  3.  
  4. c       amplification du fluage par exces de cisaillement entre rc3 et rc
  5. c       dans cette version kmecap passe de 1 pour Rc a (1+khif) pour Rc
  6. c       on a donc un fluage non lineaire mais jamais asymptotique 
  7. c       (contraitement a fluendo3d) cela pour eviter les divergence par
  8. c       fluage tertiaire qui sont gerees par la plasticite
  9.  
  10.         IMPLICIT REAL*8 (A-H,O-Z)
  11.         IMPLICIT INTEGER (I-N)
  12.  
  13.         integer nmat,nstype,nbelas
  14.         real*8 xmat(nmat),kmecap
  15.         integer vnmat(nstype)
  16.         real*8 sim3(3)
  17.  
  18.         real*8 sim3d(3),devi,trac,press
  19.         real*8 xflu,Rc,delta,coeff1
  20.  
  21.  
  22. c       coeff de fluage non lineaire
  23.         call xmat3d(xflu,xmat,nmat,vnmat,nstype,nbelas,6,1)
  24. c       coeff d amplification par fluage
  25.         if (xflu.gt.0.) then
  26. c           invariants du tenseur des contraintes 
  27.             call tp3d(sim3,simd3,devi,trac,press)
  28. c           resitance  a la compression 
  29.             call xmat3d(Rc,xmat,nmat,vnmat,nstype,nbelas,6,1)
  30. c           coeff de confinement
  31.             call xmat3d(delta,xmat,nmat,vnmat,nstype,nbelas,6,4)
  32. c           cisaillement reduit de DP
  33.             taudp=max(devi-press*delta,0.d0)
  34. c           contrainte DP au pic
  35.             coeff1=(1.d0/sqrt(3.d0)-delta/3.d0)
  36.             taulim=rc*coeff1
  37. c           coeff de fluage (variation lineaire de 1 a 1+xflu entre rc/3 et rc)
  38.             if(taudp.gt.(taulim/3.d0)) then
  39.                 if(taudp.lt.taulim) then
  40.                     kmecap=1.d0+xflu*(3.d0*taudp-taulim)/(2.d0*taulim)
  41.                 else
  42.                     kmecap=1.d0+xflu
  43.                 end if
  44.             else 
  45.                 kmecap=1.d0
  46.             end if
  47.         else
  48.             kmecap=1.d0
  49.         end if
  50.  
  51.         return
  52.         end
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  

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