Télécharger la plaquette Cast3M


fr en

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