Esope source files

   1 : C BJT3C     SOURCE    CHAT      05/01/12    21:40:43     5004           
   2 :       SUBROUTINE BJT3C(IGAU,MFR,IFOU,NIFOU,XEL,BPSS,SHPTOT,SHP,
   3 :      .                                  BGENE,DJAC,IRRT)
   4 : C=======================================================================
   5 : C
   6 : C PPj POUR LE JOINT CISAILLEMENT JCT3 UNIQUEMENT 2D 
   7 : C     CALCUL DE LA MATRICE B = TETA * ( N , -N )
   8 : C     ET DU JACOBIEN EN IGAU
   9 : C          ROUTINE FORTRAN PUR
  10 : C          P. PEGON 2001 d'apres DERIVEE DE LA ROUTINE BJO4 PAR S. FELIX
  11 : C=======================================================================
  12 : C  INPUT
  13 : C     IGAU  = NUMERO DU POINT DE GAUSS
  14 : C     XEL   = COORDONNEES DES NOEUDS DE L'ELEMENT
  15 : C     BPSS  = MATRICE DE PASSAGE
  16 : C             BPSS(,) = AXE S1
  17 : C             BPSS(,) = AXE S2
  18 : C             BPSS(,) = AXE SN
  19 : C     SHPTOT= FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE DE REFERENCE
  20 : C             SHPTOT(1, ) = FONCTION DE FORME
  21 : C             SHPTOT(2, ) = DERIVEES PAR RAPPORT A QSI
  22 : C             SHPTOT(3, ) = DERIVEES PAR RAPPORT A ETA
  23 : C  OUTPUT
  24 : C     SHP   = FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE GEOMETRIQUE
  25 : C             SHP(1, ) = FONCTION DE FORME
  26 : C             SHP(2, ) = DERIVEES PAR RAPPORT A X LOCAL
  27 : C             SHP(3, ) = DERIVEES PAR RAPPORT A Y LOCAL
  28 : C     DJAC  = JACOBIEN AU POINT D'INTEGRATION IGAU
  29 : C     BGENE = MATRICE B AU POINT D'INTEGRATION IGAU
  30 : C     IRRT  = DIFFERENT DE ZERO SI ERREUR
  31 : C  CONSTANTES
  32 : C     NBNO  = NOMBRE DE NOEUDS
  33 : C     LRE   = NOMBRE DE COLONNES DE LA MATRICE B
  34 : C     NST   = NOMBRE DE COMPOSANTES DE CONTRAINTES
  35 : C
  36 : C  REMARQUE : ATTENTION : DANS LES CAS CONTRAINTES PLANES ET
  37 : C             AXISYMETRIQUE, LA MATRICE BPSS ( RESPECTIVEMENT BGENE )
  38 : C             N'A PLUS LA DIMENSION (3X3) ( RESPECTIVEMENT (3X18) ).
  39 : C             TROUVER LEURS DIMENSIONS CORRECTES ET MODIFIER LES
  40 : C             PARAMETRES LRE, NST, NBNO. CEUX CI CORRESPONDRONT A
  41 : C             CEUX DU CAS BIDIMENSIONNELS.
  42 : C
  43 : C=======================================================================
  44 :       IMPLICIT INTEGER(I-N)
  45 :       IMPLICIT REAL*8(A-H,O-Z)
  46 :       PARAMETER (LRE=2*6,NST=2,NBNO=6)
  47 :       DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)
  48 :       DIMENSION BPSS(3,3)
  49 : C
  50 :       IRRT = 0
  51 : C
  52 : C     MATRICE JACOBIENNE
  53 : C
  54 :       DO 1 I=1,NBNO
  55 :        SHP(1,I) = SHPTOT(1,I,IGAU)
  56 :        SHP(2,I) = SHPTOT(2,I,IGAU)
  57 :        SHP(3,I) = SHPTOT(3,I,IGAU)
  58 :     1 CONTINUE
  59 : C
  60 : C !!!!!! ATTENTION : IL FAUT CALCULER LE JACOBIEN AVEC NBNO=3 !!!!!!
  61 : C        IL NE FAUT SURTOUT PAS METTRE NBNO=6 CAR CA FAUSSE ALORS
  62 : C        LES RESULTATS ... NBNO DOIT ETRE EGAL A 6 CAR IL YA 6 NOEUDS
  63 : C
  64 :       NBNONN=NBNO/2
  64 :       CALL DEVOLU(XEL,SHP,MFR,NBNONN,IFOU,NIFOU,2,1.D0,RR,DJAC)
  65 :       IF (DJAC.LT.0.0D0) THEN
  66 :          IRRT = 1
  67 :       ELSE IF (DJAC.EQ.0.0D0) THEN
  68 :          IRRT = 2
  69 :       ENDIF
  70 : C
  71 : C     MATRICE B
  72 : C
  72 :       CALL ZERO(BGENE,NST,LRE)
  73 :       DO 2 I=1,2
  74 :        DO 3 J=1,3
  75 :         DO 4 K=1,2
  76 :          L=2*(J-1)+K
  77 :          M=L+6
  78 :          BGENE(I,L)=BPSS(I,K)*SHP(1,J)
  79 :          BGENE(I,M)=-BGENE(I,L)
  80 :     4   CONTINUE
  81 :     3  CONTINUE
  82 :     2 CONTINUE
  83 : C
  84 :       RETURN
  85 :       END
  86 : 
  87 : 
  88 : 
  89 :  
  90 :