Sources Esope | Cast3M

tconv5
C TCONV5    SOURCE    CHAT      05/01/13    03:34:32     5004
 
C=======================================================================
C=                            T C O N V 5                              =
C=                            -----------                              =
C=                                                                     =
C=  Fonction :                                                         =
C=  ----------                                                         =
C=  Calcul du jacobien pour un element TRI3, TRI6, QUA4 et QUA8        =
C=                                                                     =
C=  Parametres :  (E)=Entree  (S)=Sortie                               =
C=  ------------                                                       =
C=   XEL    (E)   Coordonnees GLOBALES des noeuds de l'element         =
C=   SHP    (E)   Derivees des fonctions de FORME de l'element         =
C=   NBNN   (E)   Nombre de NOEUDS de l'element                        =
C=   DJAC   (S)   Jacobien                                             =
C=                                                                     =
C=  Denis ROBERT, le 6 mai 1988.                                       =
C=======================================================================
 
      SUBROUTINE TCONV5(XEL,SHP,NBNN,DJAC)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
 
      DIMENSION XEL(3,*),SHP(6,*)
 
      A11=0.
      A12=0.
      A13=0.
      A21=0.
      A22=0.
      A23=0.
      DO i=1,NBNN
        A11=A11+SHP(2,i)*XEL(1,i)
        A12=A12+SHP(2,i)*XEL(2,i)
        A13=A13+SHP(2,i)*XEL(3,i)
        A21=A21+SHP(3,i)*XEL(1,i)
        A22=A22+SHP(3,i)*XEL(2,i)
        A23=A23+SHP(3,i)*XEL(3,i)
      ENDDO
      A1=A12*A23-A13*A22
      A2=A13*A21-A11*A23
      A3=A11*A22-A12*A21
      DJAC=SQRT(A1*A1+A2*A2+A3*A3)
 
      RETURN
      END