Télécharger matlas.eso

Retour à la liste

Numérotation des lignes :

matlas
  1. C MATLAS SOURCE AM 09/12/08 21:15:22 6582
  2. SUBROUTINE MATLAS(CMAT,NMAT,NSTRS,IFOUR,DEP,KERRE)
  3. C
  4. C ROUTINE APPELE DANS CAMECA
  5. C_______________________________________________________________________
  6. C
  7. C FORMATION DE LA MATRICE D ELASTICITE DE HOOKE
  8. C_______________________________________________________________________
  9. C
  10. IMPLICIT INTEGER(I-N)
  11. IMPLICIT REAL*8(A-H,O-Z)
  12. DIMENSION DEP(NSTRS,NSTRS),CMAT(NMAT)
  13. C
  14.  
  15. YOUN = CMAT(1)
  16. XNU = CMAT(2)
  17. C
  18. CALL ZERO(DEP,NSTRS,NSTRS)
  19. C
  20. IF(IFOUR.EQ.-2) THEN
  21. C
  22. C CALCUL DES MATRICES D ELASTICITE SUIVANT LE TYPE DE CALCUL : DEP
  23. C - CONTRAINTES PLANES
  24. C
  25. E1=YOUN/(1.D0-XNU*XNU)
  26. DEP(1,1)=E1
  27. DEP(1,2)=XNU*E1
  28. DEP(2,1)=XNU*E1
  29. DEP(2,2)=E1
  30. DEP(3,3)=E1*(1.D0-XNU)/2.D0
  31. GOTO 100
  32. C
  33. ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  34. C
  35. C - DEFORMATION PLANE OU AXISYMETRIQUE
  36. C
  37. E1=YOUN/(1.D0+XNU)/(1.D0-2.D0*XNU)
  38. DEP(1,1)=E1*(1.D0-XNU)
  39. DEP(1,2)=E1*XNU
  40. DEP(1,3)=E1*XNU
  41. DEP(2,1)=E1*XNU
  42. DEP(2,2)=E1*(1.D0-XNU)
  43. DEP(2,3)=E1*XNU
  44. DEP(3,1)=E1*XNU
  45. DEP(3,2)=E1*XNU
  46. DEP(3,3)=E1*(1.D0-XNU)
  47. DEP(4,4)=E1*(1.D0-2.D0*XNU)/2.D0
  48. GOTO 100
  49. ELSE
  50. C
  51. C - CALCUL TRIDIMENSIONEL
  52. C
  53. E1=YOUN/(1.D0+XNU)/(1.D0-2.D0*XNU)
  54. DEP(1,1)=E1*(1.D0-XNU)
  55. DEP(1,2)=E1*XNU
  56. DEP(1,3)=E1*XNU
  57. DEP(2,1)=E1*XNU
  58. DEP(2,2)=E1*(1.D0-XNU)
  59. DEP(2,3)=E1*XNU
  60. DEP(3,1)=E1*XNU
  61. DEP(3,2)=E1*XNU
  62. DEP(3,3)=E1*(1.D0-XNU)
  63. DEP(4,4)=E1*(1.D0-2.D0*XNU)/2.D0
  64. DEP(5,5)=E1*(1.D0-2.D0*XNU)/2.D0
  65. DEP(6,6)=E1*(1.D0-2.D0*XNU)/2.D0
  66. GOTO 100
  67. ENDIF
  68. C
  69. WRITE(*,*) 'ERREUR DANS MATLAS : PB DE DIMENSION'
  70. KERRE=99
  71. 100 CONTINUE
  72. C
  73. END
  74.  
  75.  
  76.  

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