Esope source files

   1 : C HOOKSE    SOURCE    AM        08/02/14    21:24:13     6050           
   2 :       SUBROUTINE HOOKSE(VALMAT,IB,IGAU,MFR,CRIGI,IFOU,
   3 :      +            NBGMAT,NELMAT,SECT,LHOOK,DDHOOK,IRET)
   4 : C
   5 : C----------------------------------------------------------------------
   6 : C
   7 : C   Calcul de la matrice de HOOKE dans le cas d'un 
   8 : C   matériau homogeneise
   9 : C
  10 : C   Entrees:
  11 : C   --------
  12 : C   VALMAT   tableau de materiau
  13 : C   IB       numero de l'element
  14 : C   IGAU     numero du point de Gauss
  15 : C   MFR      numero de formulation
  16 : C   CRIGI    ELEMENT DE REDUCTION DE LA RIGIDITE
  17 : C   IFOU     numero d'harmonique de Fourier
  18 : C   NBGMAT, NELMAT tailles des tableaux
  19 : C   SECT     SECTION DE L'ELEMENT IB (<> 0 SI MFR.EQ.27)
  20 : C   LHOOK    taille de la matrice de HOOKE
  21 : C
  22 : C   Sorties:
  23 : C   --------
  24 : C   DDHOOK   matrice de HOOKE
  25 : C   IRET     = 1 si option existante, 0 sinon
  26 : C
  27 : C---------------------------------------------------------------------
  28 : C
  29 :       IMPLICIT INTEGER(I-N)
  30 :       IMPLICIT REAL*8(A-H,O-Z)
  31 : C
  32 : -INC CCHAMP                                                                     
  33 : -INC CCOPTIO                                                                    
  34 : C
  35 :       DIMENSION VALMAT(*)
  36 :       DIMENSION DDHOOK(LHOOK,*)
  37 :       DIMENSION CRIGI(12)
  38 : C                                                                               
  39 : C     INITIALISATION 
  40 : C                                                                               
  40 :         CALL ZERO(DDHOOK,LHOOK,LHOOK)
  41 : C
  42 :         IF (IGAU.LE.NBGMAT.AND.(IB.LE.NELMAT.OR.NBGMAT.GT.1)) THEN 
  43 : C
  44 : C  Cas des poutres et des tuyaux
  45 : C
  46 :          IF (MFR.EQ.7.OR.MFR.EQ.13) THEN 
  47 : C
  48 :             IF(IFOU.EQ.2) THEN
  49 : C
  50 : C  LA MATRICE EST PLEINE
  51 : C
  52 :             DDHOOK(1,1)=CRIGI(1)
  53 :             DDHOOK(1,5)=CRIGI(3)
  54 :             DDHOOK(5,1)=CRIGI(3)
  55 :             DDHOOK(1,6)=-CRIGI(2)
  56 :             DDHOOK(6,1)=-CRIGI(2)
  57 : C
  58 :             DDHOOK(5,5)=CRIGI(6)
  59 :             DDHOOK(5,6)=-CRIGI(5)
  60 :             DDHOOK(6,5)=-CRIGI(5)
  61 : C
  62 :             DDHOOK(6,6)=CRIGI(4)
  63 : C
  64 :             DDHOOK(2,2)=CRIGI(8)
  65 :             DDHOOK(2,4)=-CRIGI(10)
  66 :             DDHOOK(4,2)=-CRIGI(10)
  67 : C 
  68 :             DDHOOK(3,3)=CRIGI(7)
  69 :             DDHOOK(3,4)=CRIGI(9)
  70 :             DDHOOK(4,3)=CRIGI(9)
  71 : C
  72 :             DDHOOK(4,4)=CRIGI(11)+CRIGI(12)
  73 : C
  74 :            ELSEIF(IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.-3) THEN
  75 : C
  76 : C  LA MATRICE EST PLEINE
  77 : C
  78 :             DDHOOK(1,1)=CRIGI(1)
  79 :             DDHOOK(1,3)=-CRIGI(2)
  80 :             DDHOOK(3,1)=-CRIGI(2)
  81 : C
  82 :             DDHOOK(3,3)=CRIGI(3)
  83 : C
  84 :             DDHOOK(2,2)=CRIGI(4)
  85 : C
  86 :            ELSE
  87 :             IRET=0
  88 :            ENDIF
  89 : C
  90 : C  Cas des barres
  91 : C
  92 :          ELSEIF (MFR.EQ.27) THEN
  93 : C                                                                               
  94 :           YOU=VALMAT(1)                                                               
  95 :           DDHOOK(1,1)=YOU*SECT 
  96 : C
  97 :          ELSE
  98 : C                                                                               
  99 : C  Option non definie
 100 : C                                                                               
 101 :           IRET=0 
 102 :          ENDIF
 103 :         ENDIF
 104 : C
 105 :       RETURN             
 106 :       END
 107 : 
 108 : 
 109 : 
 110 : 
 111 :  
 112 :  
 113 :  
 114 :