Esope source files

   1 : C BJO4C     SOURCE    CHAT      05/01/12    21:40:33     5004           
   2 :       SUBROUTINE BJO4C(IGAU,XEL,BPSS,SHPTOT,SHP,BGENE,DJAC,IRRT)
   3 : C=======================================================================
   4 : C
   5 : C PPj POUR LE JOINT CISAILLEMENT JCI4 UNIQUEMENT 2D   
   6 : C     CALCUL DE LA MATRICE B = TETA * ( N , -N )
   7 : C     ET DU JACOBIEN EN IGAU
   8 : C          ROUTINE FORTRAN PUR
   9 : C          P. PEGON 2001 d'apres CODE S. FELIX MAI 92
  10 : C=======================================================================
  11 : C  INPUT
  12 : C     IGAU  = NUMERO DU POINT DE GAUSS
  13 : C     XEL   = COORDONNEES DES NOEUDS DE L'ELEMENT
  14 : C     BPSS  = MATRICE DE PASSAGE
  15 : C             BPSS(,) = AXE S1
  16 : C             BPSS(,) = AXE S2
  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 : ATTENTION : DANS LES CAS CONTRAINTES PLANES ET
  36 : C             AXISYMETRIQUE, LA MATRICE BPSS ( RESPECTIVEMENT BGENE )
  37 : C             N'A PLUS LA DIMENSION (3X3) ( RESPECTIVEMENT (3X24) ).
  38 : C             TROUVER LEURS DIMENSIONS CORRECTES ET MODIFIER LES
  39 : C             PARAMETRES LRE, NST, NBNO. CEUX CI CORRESPONDRONT A
  40 : C             CEUX DU CAS BIDIMENSIONNELS. AINSI PAR EXEMPLE, NBNO
  41 : C             SERA EGAL A 4, NST SERA EGAL A 2 ETC ...
  42 : C
  43 : C=======================================================================
  44 :       IMPLICIT INTEGER(I-N)
  45 :       IMPLICIT REAL*8(A-H,O-Z)
  46 :       PARAMETER (LRE=2*8,NST=2,NBNO=8)
  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=4 !!!!!!
  61 : C        IL NE FAUT SURTOUT PAS METTRE NBNO=8 CAR CA FAUSSE ALORS
  62 : C        LES RESULTATS ... NBNO DOIT ETRE EGAL A 8 CAR IL YA 8 NOEUDS
  63 : C
  64 :       NBNONN=NBNO/2
  64 :       CALL JACOBI(XEL,SHP,2,NBNONN,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,4
  75 :         DO 4 K=1,2
  76 :          L=2*(J-1)+K
  77 :          M=L+8
  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 :