Télécharger ecrp3d.eso

Retour à la liste

Numérotation des lignes :

ecrp3d
  1. C ECRP3D SOURCE PV090527 23/01/27 21:15:28 11574
  2. c***********************************************************************
  3. subroutine ecrp3d(epldpd,edppic,Cdp_min,Cdp_max,Cdp,Hdp)
  4.  
  5. c sous programme de traitement de l ecrouissage parabolique pre pic
  6. c Alain SELLIER 2021/05/05
  7.  
  8. implicit real*8 (a-h,o-z)
  9. implicit integer (i-n)
  10.  
  11. real*8 epldpd,edppic,Cdp_min,Cdp_max,Cdp,Hdp
  12. real*8 CdpM,cg,cg1,ldppic,ldp
  13.  
  14. c variable dp equivalente = integrale du multiplicateur
  15. ldp=epldpd
  16. ldppic = edppic
  17.  
  18. if(ldppic.gt.0.d0) then
  19. if(ldp.lt.ldppic) then
  20. c cohesion
  21. cg=Cdp_min
  22. CdpM=Cdp_max
  23. t1 = CdpM - cg
  24. t2 = ldp ** 2
  25. t7 = ldppic ** 2
  26. Cdp =0.1D1/t7 *(0.2D1 *ldppic * ldp * t1+ t7 * cg - t2 * t1)
  27. c module d ecrouissage
  28. t4 = ldppic ** 2
  29. Hdp = -0.2D1 / t4 * (ldp - ldppic) * (CdpM - cg)
  30. else
  31. Cdp=Cdp_max
  32. Hdp=0.d0
  33. end if
  34. else
  35. Cdp=Cdp_max
  36. Hdp=0.d0
  37. end if
  38.  
  39. return
  40. end
  41.  
  42.  
  43.  

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