Télécharger cmatem.eso

Retour à la liste

Numérotation des lignes :

  1. C CMATEM SOURCE CB215821 16/04/21 21:15:47 8920
  2. c MODELO MRSLade-Macari
  3. CCC
  4. C **********************************************************************
  5. CCC
  6. SUBROUTINE CARACMATEMAC(XMAT)
  7. IMPLICIT INTEGER(I-N)
  8. real*8 xmat(*)
  9. real*8 you,xnu,exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  10. . exl,xnn,aaa,bbb,pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,
  11. . exm2,eww2,cpr,cpm,cfr
  12. common /elasdata/ you,xnu
  13. common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  14. . exl,xnn,aaa,bbb
  15. common /MRScapp/ pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,exm2,eww2,
  16. . cpr,cpm,cfr
  17. you = xmat(1)
  18. xnu = xmat(2)
  19. pcc = xmat(5)
  20. paa = xmat(6)
  21. qaa = xmat(7)
  22. exm = xmat(8)
  23. eww = xmat(9)
  24. xk1 = xmat(10)
  25. xk2 = xmat(11)
  26. etb = xmat(12)
  27. exv = xmat(13)
  28. xep = xmat(14)
  29. xnn = xmat(15)
  30. cco = xmat(16)
  31. exl = xmat(17)
  32. pc0 = xmat(18)
  33. exr = xmat(19)
  34. cca = xmat(20)
  35. phi = xmat(21)
  36. alp = xmat(22)
  37. c exv < 0 cone fijo
  38. if (exv.le.0.D0) then
  39. bbb=1.D0/(abs(exv)*(1.D0+xk1))+xk2*xep/(1.D0+xep)/(1.D0+xep-xk2)
  40. aaa=(1.D0+xk1)**(-1.D0/abs(exv))*
  41. . (1.D0-xk2/(1.D0+xep))*etb*EXP(bbb)
  42. else
  43. bbb=1.D0/(exv*(1.D0+xk1))+xk2*xep/(1.D0+xep)/(1.D0+xep-xk2)
  44. aaa=(1.D0+xk1)**(-1.D0/exv)*(1.D0-xk2/(1.D0+xep))*etb*EXP(bbb)
  45. endif
  46. cpr=(1.D0-alp)*(phi*(1.D0-alp)+alp)/(2.D0*phi*(1.D0-alp)+alp)
  47. cpm=(alp*alp+phi*(1.D0-alp*alp))/(2.D0*phi*(1.D0-alp)+alp)
  48. cfr=(phi*(1.D0-alp)+alp)*SQRT(alp/(2.D0*phi*(1.D0-alp)+alp))
  49. pcc2=pcc
  50. paa2=paa
  51. qaa2=qaa
  52. exm2=exm
  53. eww2=eww
  54. return
  55. end
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  

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