Télécharger la plaquette Cast3M


fr en

Télécharger hccd3d.eso

Retour à la liste

Numérotation des lignes : 

hccd3d

  1. C HCCD3D    SOURCE    PV090527  23/01/27    21:15:36     11574          
  2.       SUBROUTINE HCCD3D(pfcc,ppcc,poro0,epccv,Hcc)
  3. c     Module d ecrouissage plastique tangent pour CamClay 
  4. c     hcc=-dpc/depccv
  5. c     (A.Sellier 2021/04/23)
  6.       implicit real*8 (a-h,o-z)
  7.       implicit integer (i-n)
  8.  
  9.       real*8 pfcc,ppcc,poro0,epccv,Hcc
  10.  
  11. c     module de depart
  12. C       t2 = pfcc / ppcc
  13. C       t3 = 0.1D1 / poro0
  14. C       t5 = -0.1D1 + poro0 + epccv
  15. C       t8 = t2 ** (0.1D1 / t5 * t3 * epccv)
  16. C       t12 = t5 ** 2
  17. C       t14 = log(t2)
  18. C       Hcc = -ppcc * t14 / t12 * t3 * (-0.1D1 + poro0) * t8
  19.  
  20.       t2 = pfcc / ppcc
  21.       t3 = -0.1D1 + poro0 + epccv
  22.       t6 = 0.1D1 / poro0
  23.       t8 = t2 ** (t6 / t3 * epccv)
  24.       t11 = t3 ** 2
  25.       t14 = log(t2)
  26.       Hcc = -ppcc * t14 * t6 / t11 * (-0.1D1 + poro0) * t8
  27.  
  28. c     on test l ecoulement avec le module le plus elevé pour ne jamais passer sous la courbe
  29. c      t6 = log(pfcc / ppcc)
  30. c      Hcc = -t6 / poro0 * pfcc * (-0.1D1 + poro0)
  31.  
  32.  
  33.  
  34. c      print*,'Hcc=',hcc      
  35.  
  36.       return
  37.       end
  38.  
  39.  
  40.  
  41.  
  42.  

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