Télécharger la plaquette Cast3M


fr en

Télécharger prandt.eso

Retour à la liste

Numérotation des lignes : 

prandt

  1. C PRANDT    SOURCE    CHAT      05/01/13    02:21:23     5004
  2. c*****************************************************************************
  3. c       sous rouine utilisee dans l'ecoulement plastique de la loi
  4. c       de GURSON
  5. c*********************************************************************
  6.         subroutine prandt (wrkgur)
  7. c
  8. c       this subroutine is used whenever phin < phi0 and Pn+1(phin)>0
  9. c       it performs a classical radial return algorithm for the Prandtl Reuss model
  10. c       rho and sqrtj2 must have been updated for the trial state before the call
  11. c
  12. c---    variables
  13.       IMPLICIT INTEGER(I-N)
  14.         IMPLICIT REAL*8 (A-H,O-Z)
  15.         segment wrkgur
  16.          real*8 sigbar, sy0,phi0,rho0,g,b,h
  17.          real*8 epn,phin,sqrtj2,rho,sig(6)
  18.          real*8 e(7),dt
  19.          real*8 conv,tol1,tol2
  20.         endsegment
  21.  
  22. *        common/prop/sigbar,sy0,phi0,rho0,g,b,h
  23. *        common/state/epn,phin,sqrtj2,rho,sig(6)
  24. c               e is the strain rate:
  25. c                       e(1-6) is the deviatoric part
  26. c                       e(7) is the trace/3
  27. c               dt ist the time increment
  28. *        common /delta/ e(7),dt
  29.  
  30. c
  31. c---    equivalent strain increment
  32.         dep=(sqrtj2-sy0-h*epn)/(h+3.d0*g)
  33.         epn=epn+dep
  34.         sqrtj2=sqrtj2-3.d0*g*dep
  35. c
  36.         fac=1.d0+3.d0*g*dep/sqrtj2
  37.         do 100 i=1,6
  38.         sig(i)=sig(i)/fac
  39.  100    continue
  40. c
  41.         return
  42.         end
  43. c*****************************************************************************
  44.  
  45.  
  46.  
  47.  
  48.  

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