Esope source files

   1 : C JO3LOC    SOURCE    CHAT      05/01/13    00:49:44     5004           
   2 :       SUBROUTINE JO3LOC(XE,SHPTOT,IGAU,NBNO,BPSS)
   3 : C=======================================================================
   4 : C
   5 : C    -CALCUL DE LA MATRICE DE PASSAGE BPSS
   6 : C          ROUTINE FORTRAN PUR
   7 : C    Pierre Pegon (ISPRA) OCTOBER 93
   8 : C=======================================================================
   9 : C  INPUT
  10 : C     XE     = COORDONNEES  DE L ELEMENT
  11 : C     SHPTOT = FONCTIONS DE FORME
  12 : C            = SHPTOT(1,...) = FONCTIONS DE FORME
  13 : C            = SHPTOT(2,...) = DERIVEE PAR RAPPORT A QSI
  14 : C            = SHPTOT(3,...) = DERIVEE PAR RAPPORT A ETA
  15 : C     IGAU   = NUMERO DU POINT DE GAUSS
  16 : C     NBNO   = NOMBRE DE NOEUDS DE L'ELEMENT
  17 : C  OUTPUT
  18 : C     BPSS   = MATRICE DE PASAGE REPERE GLOBAL/REPERE LOCAL
  19 : C
  20 : C=======================================================================
  21 :       IMPLICIT INTEGER(I-N)
  22 :       IMPLICIT REAL*8(A-H,O-Z)
  23 :       DIMENSION XE(3,*),BPSS(3,3),SHPTOT(6,NBNO,*)
  24 :       DIMENSION XD(2,4)
  25 :       DIMENSION S1(2),SN(2)
  26 : C
  27 :       NOQUAL = 0
  28 : C
  29 : C----------     CALCUL DE LA MATRICE DE  PASSAGE
  30 : C
  31 :       DO 6 I=1,2
  32 :          S1(I)=0.0D0
  33 :          SN(I)=0.0D0
  34 :     6 CONTINUE
  35 : C
  36 :       DO 7 I=1,NBNO
  37 : C
  38 : C-------------------TANGENTE AU POINT DE GAUSS IGAU SELON QSI
  39 : C
  40 :          S1(1) = S1(1) + ( SHPTOT(2,I,IGAU)*XE(1,I) )
  41 :          S1(2) = S1(2) + ( SHPTOT(2,I,IGAU)*XE(2,I) )
  42 :     7  CONTINUE
  43 : C
  44 : C-------------------NORMALISATION DE S1
  45 : C
  46 :        XNORME = SQRT((S1(1)**2) + (S1(2)**2))
  47 :        S1(1) = S1(1) / XNORME
  48 :        S1(2) = S1(2) / XNORME
  49 : C
  50 : C-------------------NORMALE AU JOINT (PAR ROTATION DE 90 DEGRES)
  51 : C
  52 :        SN(1) =-S1(2)
  53 :        SN(2) = S1(1)
  54 : C
  55 : C-------------------STOCKAGE DE LA MATRICE DE PASSAGE
  56 : C
  57 :        DO 10 I=1,2
  58 :           BPSS(1,I) = S1(I)
  59 :           BPSS(2,I) = SN(I)
  60 :    10  CONTINUE
  61 : C
  62 :       RETURN
  63 :       END
  64 :  
  65 :