Télécharger pccd3d.eso

Retour à la liste

Numérotation des lignes :

pccd3d
  1. C PCCD3D SOURCE PV090527 23/02/13 21:15:12 11592
  2. SUBROUTINE PCCD3D(precision3d,poro0,epccv,ppcc,pfcc,pc,rt,pt)
  3. c pressions limites actuelle de CamClay
  4. c (A.Sellier 2021/04/23)
  5. implicit real*8 (a-h,o-z)
  6. implicit integer (i-n)
  7.  
  8. real*8 precision3d,poro0,EPLCCV,ppcc,pfcc,pc,rt,pt
  9. real*8 pc1,pc2
  10.  
  11. c pression de consolidation actuelle
  12. if(epccv.le.0.d0) then
  13. c consolidation
  14. if(epccv+poro0*(1.d0-2.d0*precision3d).gt.0.d0) then
  15. c la pression de consolidation peut encore évoluer
  16. t8 = (pfcc / ppcc) ** (epccv / (-0.1D1 + poro0 + epccv) / poro0)
  17. pc1 = ppcc * t8
  18. pc2=max(ppcc,pc,pc1)
  19. pc=min(pfcc,pc2)
  20. if(pc.ne.pc) then
  21. print*,'pb dans pccd3d appele par endo3d'
  22. end if
  23. else
  24. c on a atteint la pression limite
  25. pc=pfcc
  26. end if
  27. else
  28. c deconsolidation
  29. c interdite en dessous valeur initiale
  30. pc=max(ppcc,pc)
  31. pc=min(pfcc,pc)
  32. end if
  33.  
  34. c pression limite de traction triaxiale constante
  35. pt = -rt
  36.  
  37.  
  38. return
  39. end
  40.  
  41.  
  42.  
  43.  
  44.  

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