Esope source files

   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
  17 :       CALL ZERO(DEP,NSTRS,NSTRS)
  18 : C
  19 :       IF(IFOUR.EQ.-2) THEN
  20 : C
  21 : C CALCUL DES MATRICES D ELASTICITE SUIVANT LE TYPE DE CALCUL : DEP
  22 : C  - CONTRAINTES PLANES 
  23 : C
  24 :          E1=YOUN/(1.D0-XNU*XNU)
  25 :          DEP(1,1)=E1
  26 :          DEP(1,2)=XNU*E1
  27 :          DEP(2,1)=XNU*E1
  28 :          DEP(2,2)=E1
  29 :          DEP(3,3)=E1*(1.D0-XNU)/2.D0
  30 :          GOTO 100
  31 : C
  32 :       ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  33 : C
  34 : C  - DEFORMATION PLANE OU AXISYMETRIQUE 
  35 : C
  36 :          E1=YOUN/(1.D0+XNU)/(1.D0-2.D0*XNU)
  37 :          DEP(1,1)=E1*(1.D0-XNU)
  38 :          DEP(1,2)=E1*XNU
  39 :          DEP(1,3)=E1*XNU
  40 :          DEP(2,1)=E1*XNU
  41 :          DEP(2,2)=E1*(1.D0-XNU)
  42 :          DEP(2,3)=E1*XNU
  43 :          DEP(3,1)=E1*XNU
  44 :          DEP(3,2)=E1*XNU
  45 :          DEP(3,3)=E1*(1.D0-XNU)
  46 :          DEP(4,4)=E1*(1.D0-2.D0*XNU)/2.D0
  47 :          GOTO 100
  48 :       ELSE
  49 : C
  50 : C  - CALCUL TRIDIMENSIONEL 
  51 : C
  52 :          E1=YOUN/(1.D0+XNU)/(1.D0-2.D0*XNU)
  53 :          DEP(1,1)=E1*(1.D0-XNU)
  54 :          DEP(1,2)=E1*XNU
  55 :          DEP(1,3)=E1*XNU
  56 :          DEP(2,1)=E1*XNU
  57 :          DEP(2,2)=E1*(1.D0-XNU)
  58 :          DEP(2,3)=E1*XNU
  59 :          DEP(3,1)=E1*XNU
  60 :          DEP(3,2)=E1*XNU
  61 :          DEP(3,3)=E1*(1.D0-XNU)
  62 :          DEP(4,4)=E1*(1.D0-2.D0*XNU)/2.D0
  63 :          DEP(5,5)=E1*(1.D0-2.D0*XNU)/2.D0
  64 :          DEP(6,6)=E1*(1.D0-2.D0*XNU)/2.D0
  65 :          GOTO 100
  66 :       ENDIF
  67 : C
  68 :       WRITE(*,*) 'ERREUR DANS MATLAS : PB DE DIMENSION'
  69 :       KERRE=99
  70 :   100 CONTINUE
  71 : C
  72 :       END
  73 :  
  74 :