Esope source files

   1 : C TCONV5    SOURCE    CHAT      05/01/13    03:34:32     5004           
   2 : 
   3 : C=======================================================================
   4 : C=                            T C O N V 5                              =
   5 : C=                            -----------                              =
   6 : C=                                                                     =
   7 : C=  Fonction :                                                         =
   8 : C=  ----------                                                         =
   9 : C=  Calcul du jacobien pour un element TRI3, TRI6, QUA4 et QUA8        =
  10 : C=                                                                     =
  11 : C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  12 : C=  ------------                                                       =
  13 : C=   XEL    (E)   Coordonnees GLOBALES des noeuds de l'element         =
  14 : C=   SHP    (E)   Derivees des fonctions de FORME de l'element         =
  15 : C=   NBNN   (E)   Nombre de NOEUDS de l'element                        =
  16 : C=   DJAC   (S)   Jacobien                                             =
  17 : C=                                                                     =
  18 : C=  Denis ROBERT, le 6 mai 1988.                                       =
  19 : C=======================================================================
  20 : 
  21 :       SUBROUTINE TCONV5(XEL,SHP,NBNN,DJAC)
  22 : 
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8 (A-H,O-Z)
  25 : 
  26 :       DIMENSION XEL(3,*),SHP(6,*)
  27 : 
  28 :       A11=0.
  29 :       A12=0.
  30 :       A13=0.
  31 :       A21=0.
  32 :       A22=0.
  33 :       A23=0.
  34 :       DO i=1,NBNN
  35 :         A11=A11+SHP(2,i)*XEL(1,i)
  36 :         A12=A12+SHP(2,i)*XEL(2,i)
  37 :         A13=A13+SHP(2,i)*XEL(3,i)
  38 :         A21=A21+SHP(3,i)*XEL(1,i)
  39 :         A22=A22+SHP(3,i)*XEL(2,i)
  40 :         A23=A23+SHP(3,i)*XEL(3,i)
  41 :       ENDDO
  42 :       A1=A12*A23-A13*A22
  43 :       A2=A13*A21-A11*A23
  44 :       A3=A11*A22-A12*A21
  45 :       DJAC=SQRT(A1*A1+A2*A2+A3*A3)
  46 : 
  47 :       RETURN
  48 :       END
  49 : 
  50 :  
  51 :