Esope source files

   1 : C BDBSTX    SOURCE    CHAT      08/09/17    21:15:07     6160           
   2 : C BDBST     SOURCE    CHAT      05/01/12    21:37:58     5004           
   3 : 
   4 : C=======================================================================
   5 : C=                             B D B S T                               =
   6 : C=                             ---------                               =
   7 : C=                                                                     =
   8 : C=  Fonction :                                                         =
   9 : C=  ----------                                                         =
  10 : C=   Calcul de la contribution d'un point de Gauss a la matrice de     =
  11 : C=   RIGIDITE elastique REL d'un element fini.                         =
  12 : C=   Contribution : Jacobien.transpose(BGENE)*DHOOK*BGENE              =
  13 : C=                                                                     =
  14 : C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  15 : C=  ------------                                                       =
  16 : C=   BGENE    (E)   Matrice B                                          =
  17 : C=   DJac     (E)   "Volume" associe au point de Gauss                 =
  18 : C=   DHOOK    (E)   Matrice de Hooke au point de Gauss                 =
  19 : C=   LRE      (E)   Nombre de DDL de l'element fini                    =
  20 : C=   NSTRS    (E)   Nombre de composantes de contraintes/deformations  =
  21 : C=   REL     (E/S)  Matrice de rigidite de l'element fini              =
  22 : C=                                                                     =
  23 : C=   Remarque : Actuellement, NSTRS doit etre inferieur ou egal a 10 ! =
  24 : C=======================================================================
  25 : 
  26 :       SUBROUTINE BDBSTX(BGENE,DJAC,DHOOK,LRELOC,NSTRS,LRE,REL)
  27 : 
  28 :       IMPLICIT INTEGER(I-N)
  29 :       IMPLICIT REAL*8 (A-H,O-Z)
  30 : 
  31 : -INC CCREEL
  32 : 
  33 :       DIMENSION DHOOK(NSTRS,*),REL(LRE,*),BGENE(NSTRS,*)
  34 :       DIMENSION WORK1(10),WORK2(10)
  35 : 
  36 : C= Boucle sur les DDL
  37 :       DO IX=1,LRELOC
  38 : C== Calcul de JACobien*Poids*B pour le DDL IX
  39 :         DO i=1,NSTRS
  40 :           WORK1(i)=DJAC*BGENE(i,IX)
  41 :         ENDDO
  42 : C== Calcul de Vol.*B*D pour le DDL IX
  43 :         DO i=1,NSTRS
  44 :           CC=XZero
  45 :           DO j=1,NSTRS
  46 :             CC=CC+WORK1(j)*DHOOK(i,j)
  47 :           ENDDO
  48 :           WORK2(i)=CC
  49 :         ENDDO
  50 : C== Calcul de Vol.Bt.D.B pour les DDL IX et IY
  51 :         DO IY=1,IX
  52 :           CC=XZero
  53 :           DO i=1,NSTRS
  54 :             CC=CC+BGENE(i,IY)*WORK2(i)
  55 :           ENDDO
  56 :           REL(IX,IY)=REL(IX,IY)+CC
  57 :           REL(IY,IX)=REL(IX,IY)
  58 :         ENDDO
  59 :       ENDDO
  60 : 
  61 :       RETURN
  62 :       END
  63 : 
  64 :  
  65 :  
  66 :