Esope source files

   1 : C BMAT28    SOURCE    FANDEUR   10/08/31    21:15:16     6735           
   2 : 
   3 : C=======================================================================
   4 : C=                          B M A T 2 8                                =
   5 : C=                          -----------                                =
   6 : C=  Fonction :                                                         =
   7 : C=  ----------                                                         =
   8 : C=   Calcul de la matrice B reliant les deformations en un point d'un  =
   9 : C=   element fini aux ddls de deplacement aux noeuds de cet element    =
  10 : C=   Le jacobien est egalement evalue au point de Gauss pour verifier  =
  11 : C=   ulterieurement si l'element fini n'est pas trop distordu.         =
  12 : C=                                                                     =
  13 : C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  14 : C=  ------------                                                       =
  15 : C=   iGau     (E)   Numero du point de Gauss considere                 =
  16 : C=   NBPGAU   (E)   Nombre de points de Gauss de l'element fini        =
  17 : C=   POIGAU,QSIGAU   (E)   |  Poids et coordonnees de reference des    =
  18 : C=   ETAGAU,DZEGAU   (E)   |  differents points de Gauss de l'element  =
  19 : C=   ITEL     (E)   Type de l'element fini (cf. NOMTP dans bdata.eso)  =
  20 : C=   MFR      (E)   Formulation associee a l'element fini              =
  21 : C=   NBNO     (E)   Nombre de noeuds de l'element fini                 =
  22 : C=   LRE      (E)   Nombre de DDL de l'element fini                    =
  23 : C=   IFOU     (E)   Mode de calcul utilise (cf. IFOUR dans CCOPTIO)    =
  24 : C=   NST      (E)   Nombre de composantes de deformations              =
  25 : C=   NN       (E)   Numero du mode de Fourier (si IFOU=1)              =
  26 : C=   DIM3     (E)   Epaisseur dans l'hypothese des contraintes planes  =
  27 : C=   XEL      (E)   Coordonnees des noeuds de l'element fini etudie    =
  28 : C=   SHPTOT   (E)   Fonctions de forme et leurs derivees               =
  29 : C=   SHP      (S)   Fonctions de forme et leurs derivees actuelles     =
  30 : C=   BGENE    (S)   Matrice de gradients (B) calcule au point de Gauss =
  31 : C=   DJAC     (S)   Jacobien au point de Gauss etudie                  =
  32 : C=   XDPGE,YDPGE   (E)   Coordonnees du noeud support (mode PLAN GENE) =
  33 : C=                                                                     =
  34 : C=  Remarque :                                                         =
  35 : C=  ----------                                                         =
  36 : C=   Lors de l'entree dans le sous-programme, SHPTOT(2 a 4,*) contient =
  37 : C=   les DERIVEES des fonctions de forme par rapport aux coordonnees   =
  38 : C=   de REFERENCE Qsi,Eta,Dzeta.                                       =
  39 : C=   En sortie du sous-programme, SHP(2 a 4,*) contient les DERIVEES   =
  40 : C=   des fonctions de FORME par rapport aux coordonnees REELLES x,y,z. =
  41 : C=======================================================================
  42 : 
  43 :       SUBROUTINE BMAT28 (iGau,NBPGAU,POIGAU,QSIGAU,ETAGAU,DZEGAU,
  44 :      &                   ITEL,MFR,NBNO,LRE,IFOU,NST,NN,DIM3,XEL,
  45 :      &                   SHPTOT,SHP,BGENE,DJAC,XDPGE,YDPGE)
  46 : 
  47 :       IMPLICIT INTEGER(I-N)
  48 :       IMPLICIT REAL*8 (A-H,O-Z)
  49 : 
  50 : -INC CCREEL
  51 : 
  52 :       DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*),
  53 :      &          POIGAU(*),QSIGAU(*),ETAGAU(*),DZEGAU(*)
  54 : 
  55 :       DIMENSION BB1(3,9),GEOM(20),XX(3),YY(3)
  56 : 
  57 :       DATA XX / .5D0,.0D0,.5D0 /
  58 :       DATA YY / .0D0,.5D0,.5D0 /
  59 : 
  60 :       IF (ITEL.NE.28.AND.ITEL.NE.45) RETURN
  61 : 
  62 :       RR=XZero
  62 :       CALL ZERO(BGENE,NST,LRE)
  63 : 
  64 : C ===
  65 : C  Element COQUE DKT
  66 : C ===
  67 :       DO i=1,NBNO
  68 :         SHP(1,i)=SHPTOT(1,i,iGau)
  69 :         SHP(2,i)=SHPTOT(2,i,iGau)
  70 :         SHP(3,i)=SHPTOT(3,i,iGau)
  71 :       ENDDO
  71 :       CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC)
  72 :       k = 1
  73 :       DO i = 1, NBNO
  74 :         BGENE(1,k) = SHP(2,i)
  75 :         BGENE(2,k) = XZero
  76 :         BGENE(3,k) = SHP(3,i)
  77 :         BGENE(1,k+1) = XZero
  78 :         BGENE(2,k+1) = SHP(3,i)
  79 :         BGENE(3,k+1) = SHP(2,i)
  80 :         K=K+6
  81 :       ENDDO
  81 :       CALL GEOCST(XEL,GEOM)
  81 :       CALL BFDKT(XX(iGau),YY(iGau),GEOM,BB1)
  82 :       DJAC=GEOM(17)
  83 :       K=2
  84 :       KK=0
  85 :       DO i=1,3
  86 :         DO IX=1,3
  87 :           DO IY=1,3
  88 :             BGENE(3+IX,K+IY)=BB1(IX,IY+KK)
  89 :           ENDDO
  90 :         ENDDO
  91 :         K=K+6
  92 :         KK=KK+3
  93 :       ENDDO
  94 : 
  95 :       RETURN
  96 :       END
  97 : 
  98 :