Télécharger cmp1de.eso

Retour à la liste

Numérotation des lignes :

  1. C CMP1DE SOURCE CHAT 05/01/12 22:11:00 5004
  2. CCC
  3. C **********************************************************************
  4. CCC
  5. SUBROUTINE CARAC_MATE_POWDER1_DENSI(XMAT,XDENSI)
  6. IMPLICIT INTEGER(I-N)
  7. real*8 xmat(*),xdensi,eta
  8. real*8 you,xnu,sigy0,siginf,kiso,velo,cpar,mpar
  9. common /elasdata/ you,xnu
  10. common /miehdata/ sigy0,kiso,siginf,velo,cpar,mpar
  11. real*8 kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
  12. logical ogden
  13. common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
  14. . ogden
  15. real*8 nn1,nn2,eta0,aa1,aa2
  16. common /poder1_data/ nn1,nn2,eta0,aa1,aa2
  17.  
  18. ogden =.false.
  19. you = xmat( 1)
  20. xnu = xmat( 2)
  21. sigy0 = xmat( 5)
  22. nn1 = xmat( 6)
  23. nn2 = xmat( 7)
  24. eta0 = xmat( 8)
  25. eta = xdensi*eta0
  26. aa1 = (1.D0 - eta**2)/(2.D0 + eta**2)
  27. if (aa1.lt.1.D-30) aa1=1.D-30
  28. aa1 = aa1**nn1
  29. if (aa1.lt.1.D-5) aa1=1.D-5
  30. aa2 = abs(eta - 0.98*eta0)/(1.D0 - 0.98*eta0)
  31. if (eta.lt.eta0) aa2=0.02D0*eta0/(1.D0 - 0.98D0*eta0)
  32. aa2 = aa2**nn2
  33.  
  34. siginf = 0.D0
  35. kiso = 0.D0
  36. velo = 0.D0
  37. cpar = 0.D0
  38. mpar = 1.D0
  39. return
  40. end
  41.  
  42.  
  43.  
  44.  

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