Télécharger camarh.eso

Retour à la liste

Numérotation des lignes :

  1. C CAMARH SOURCE CB215821 16/04/21 21:15:31 8920
  2. c MODELO RHMC
  3. CCC
  4. C **********************************************************************
  5. CCC
  6. SUBROUTINE CARAC_MATE_RHMC (XMAT)
  7. IMPLICIT INTEGER(I-N)
  8. real*8 xmat(*)
  9. real*8 you,xnu,coh,phi,ar,tmax,aap,aan,bbp,bbn
  10. common /elasdata/ you,xnu
  11. common /rhmcdata/ coh,phi,ar,tmax,aap,aan,bbp,bbn
  12. real*8 kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
  13. logical ogden
  14. common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
  15. . ogden
  16. ogden=.false.
  17. you = xmat( 1)
  18. xnu = xmat( 2)
  19. c en anteriores versiones el dato era realmente coh*sqrt(2/3)
  20. c pq se multiplicaba luego por ese valor para usarlo en todos lados !!
  21. coh = xmat( 5)
  22. c *SQRT(2.D0/3.D0)
  23. phi = xmat( 6)
  24. ar = 0.05D0*cos(phi)
  25. tmax = 25.D0*atan(1.D0)/45.D0
  26. aap = cos(tmax)/3.D0*(3.D0+tan(tmax)*tan(3.D0*tmax)+
  27. . sin(phi)*(tan(3.D0*tmax)-3.D0*tan(tmax))/sqrt(3.D0))
  28. aan = cos(tmax)/3.D0*(3.D0+tan(tmax)*tan(3.D0*tmax)-
  29. . sin(phi)*(tan(3.D0*tmax)-3.D0*tan(tmax))/sqrt(3.D0))
  30. bbp = (sin(phi)*cos(tmax)/sqrt(3.D0)+sin(tmax))
  31. . /3.D0/cos(3.D0*tmax)
  32. bbn = (sin(phi)*cos(tmax)/sqrt(3.D0)-sin(tmax))
  33. . /3.D0/cos(3.D0*tmax)
  34. return
  35. end
  36.  
  37.  
  38.  
  39.  
  40.  
  41.  

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