Esope source files

   1 : C BJO3      SOURCE    CHAT      05/01/12    21:40:29     5004           
   2 :       SUBROUTINE BJO3(IGAU,MFR,IFOU,NIFOU,XEL,BPSS,SHPTOT,SHP,
   3 :      .                                  BGENE,DJAC,IRRT)
   4 : C=======================================================================
   5 : C
   6 : C     CALCUL DE LA MATRICE B = TETA * ( N , -N )
   7 : C     ET DU JACOBIEN EN IGAU POUR UN JOINT 2D
   8 : C          ROUTINE FORTRAN PUR
   9 : C     Pierre Pegon (ISPRA) 28/10/93
  10 : C=======================================================================
  11 : C  INPUT
  12 : C     IGAU  = NUMERO DU POINT DE GAUSS
  13 : C     MFR   = NUMERO DE LA FORMULATION 
  14 : C     XEL   = COORDONNEES DES NOEUDS DE L'ELEMENT
  15 : C     BPSS  = MATRICE DE PASSAGE
  16 : C             BPSS(,) = AXE S1
  17 : C             BPSS(,) = AXE SN
  18 : C     SHPTOT= FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE DE REFERENCE
  19 : C             SHPTOT(1, ) = FONCTION DE FORME
  20 : C             SHPTOT(2, ) = DERIVEES PAR RAPPORT A QSI
  21 : C             SHPTOT(3, ) = DERIVEES PAR RAPPORT A ETA
  22 : C  OUTPUT
  23 : C     SHP   = FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE GEOMETRIQUE
  24 : C             SHP(1, ) = FONCTION DE FORME
  25 : C             SHP(2, ) = DERIVEES PAR RAPPORT A X LOCAL
  26 : C             SHP(3, ) = DERIVEES PAR RAPPORT A Y LOCAL
  27 : C     DJAC  = JACOBIEN AU POINT D'INTEGRATION IGAU
  28 : C     BGENE = MATRICE B AU POINT D'INTEGRATION IGAU
  29 : C     IRRT  = DIFFERENT DE ZERO SI ERREUR
  30 : C  CONSTANTES
  31 : C     NBNO  = NOMBRE DE NOEUDS
  32 : C     LRE   = NOMBRE DE COLONNES DE LA MATRICE B
  33 : C     NST   = NOMBRE DE COMPOSANTES DE CONTRAINTES
  34 : C
  35 : C  REMARQUE : CE SOUS-PROGRAMME EST VALABLE EN CONTRAINTES PLANES,
  36 : C             DEFORMATIONS PLANES ET AXISYMETRIQUES. LA MATRICE
  37 : C                     BGENE = TETA ( Ntop , - Nbot )
  38 : C             EST UNE MATRICE (2x12). EN EFFET, TETA EST UNE MATRICE
  39 : C             (2x2), ET ( Ntop , - Nbot ) EST UNE MATRICE (2x12).
  40 : C             LA MATRICE TETA EST ICI APPELEE BPSS.
  41 : C
  42 : C=======================================================================
  43 :       IMPLICIT INTEGER(I-N)
  44 :       IMPLICIT REAL*8(A-H,O-Z)
  45 :       PARAMETER (LRE=12,NST=2,NBNO=6)
  46 :       DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)
  47 :       DIMENSION BPSS(3,3)
  48 : C
  49 :       IRRT = 0
  50 : C
  51 : C     JACOBIEN
  52 : C
  53 :       DO 1 I=1,NBNO
  54 :        SHP(1,I) = SHPTOT(1,I,IGAU)
  55 :        SHP(2,I) = SHPTOT(2,I,IGAU)
  56 :     1 CONTINUE
  57 : C
  58 :       NBNONN=NBNO/2
  58 :       CALL DEVOLU(XEL,SHP,MFR,NBNONN,IFOU,NIFOU,86,1.D0,RR,DJAC)
  59 : C
  60 :       IF (DJAC.EQ.0.0D0) THEN
  61 :          IRRT = 2
  62 :       ENDIF
  63 : C
  64 : C     MATRICE B ATTENTION AU RESPECT DE LA NUMEROTATION LOCALE 
  65 : C     DES NOEUDS !!!
  66 : C
  66 :       CALL ZERO(BGENE,NST,LRE)
  67 :       DO 2 I=1,2
  68 :        DO 3 J=1,3
  69 :         DO 4 K=1,2
  70 :          L=2*(J-1)+K
  71 : C CE N'EST PAS M=L+6 ....
  72 :          M=L+2*(2*(3-J)+1)
  73 :          BGENE(I,L)=BPSS(I,K)*SHP(1,J)
  74 :          BGENE(I,M)=-BGENE(I,L)
  75 :     4   CONTINUE
  76 :     3  CONTINUE
  77 :     2 CONTINUE
  78 : C
  79 :       RETURN
  80 :       END
  81 :  
  82 :