Télécharger dpfunc.eso

Retour à la liste

Numérotation des lignes :

  1. C DPFUNC SOURCE CB215821 16/04/21 21:16:30 8920
  2. FUNCTION DPfunc(sigma,alpha)
  3.  
  4. c Drucker Prager function to describe the compressive part of a concrete material.
  5.  
  6. c sigmaeq = SQRT(3 J2) + alpha I1 with J2 : 2nd invariant of the deviatoric stress
  7. c and I1 : 1st invariant of the stress
  8.  
  9. IMPLICIT REAL*8 (A-B,D-H,O-Z)
  10. implicit integer (I-K,M,N)
  11. implicit logical (L)
  12. implicit character*10 (C)
  13.  
  14. dimension sigma(6)
  15.  
  16. sigmaxx = sigma(1)
  17. sigmayy = sigma(2)
  18. sigmazz = sigma(3)
  19. sigmaxy = sigma(4)
  20. sigmaxz = sigma(5)
  21. sigmayz = sigma(6)
  22.  
  23. rloc = (sigmaxx-sigmayy)**2 + (sigmayy-sigmazz)**2 +
  24. . (sigmazz-sigmaxx)**2
  25. rloc2 = ((1d0)/(6d0))*rloc
  26. rj2 = rloc2 + sigmaxy**2 + sigmayz**2 + sigmaxz**2
  27. ri1 = sigmaxx + sigmayy + sigmazz
  28. DPfunc = SQRT(3d0*rj2) + alpha*ri1
  29.  
  30. return
  31. end
  32.  
  33.  
  34.  
  35.  
  36.  
  37.  
  38.  

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