Télécharger la plaquette Cast3M


fr en

Télécharger detjac.eso

Retour à la liste

Numérotation des lignes : 

detjac

  1. C DETJAC    SOURCE    CHAT      05/01/12    22:44:08     5004
  2. C
  3.       SUBROUTINE DETJAC(XNOEU,XY,ITYP,IDIM,NOEL1,IFORMU,SHP,DET)
  4. C*******************************************************************
  5. C     CALCUL DU JACOBIEN POUR LE PASSAGE D'UN SCALAIRE  DE
  6. C            L'ELEMENT DE REFERENCE A L'ELEMENT REEL
  7. C********************************************************************
  8. C   ENTREES
  9. C
  10. C     XY COORDONNEES DE REFERENCES DU POINT DE CALCUL
  11. C     ITYP TYPE DE L ELEMENT
  12. C     XNOEU COORDONNEES REELLES DES NOEUDS DE L ELEMENT
  13. C     NOEL1  NOMBRE DE NOEUDS DE L ELEMENT POSITION
  14. C
  15. C           COORDONNEES DE REFERENCE AU POINT DE CALCUL
  16. C
  17. C     SHP(1,I) FONCTIONS DE FORMES AU POINT DE CALCUL
  18. C     SHP(2,I) DERIVEE DES FONCTIONS DE FORME / XSI AU POINT DE CALCUL
  19. C     SHP(3,I)    "     "      "     "    "   / ETA "    "   "    "
  20. C     SHP(4,I)    "     "      "      "   "   /DZETA "   "    "   "
  21. C
  22. C   SORTIE
  23. C     DET VALEUR DU JACOBIEN
  24. C ********************************************************************
  25. C
  26. C
  27.       IMPLICIT INTEGER(I-N)
  28.       IMPLICIT REAL*8 (A-H,O-Z)
  29. C
  30. C
  31. C
  32.       DIMENSION XNOEU(3,*),SHP(6,*),XY(*)
  33. C
  34. C
  35. C***  CALCUL DES XI ET DES DXI
  36. C
  37.       CALL SHAPE(XY(1),XY(2),XY(3),ITYP,SHP,IRET)
  38. C
  39. C***  COEFFICIENTS DE DERIVATION
  40. C
  41.       DET=0.D0
  42.       IF(IDIM.EQ.2)THEN
  43. C
  44. C***  BIDIM
  45. C
  46.             A=0.D0
  47.             B=0.D0
  48.             C=0.D0
  49.             D=0.D0
  50.             DO 63 I=1,NOEL1
  51.                  A=A+SHP(2,I)*XNOEU(1,I)
  52.                  B=B+SHP(3,I)*XNOEU(1,I)
  53.                  C=C+SHP(2,I)*XNOEU(2,I)
  54.                  D=D+SHP(3,I)*XNOEU(2,I)
  55.   63        CONTINUE
  56.             DET=A*D-B*C
  57. C
  58. C***  TRIDIM
  59. C
  60.       ELSE
  61.             A1=0.D0
  62.             A2=0.D0
  63.             A3=0.D0
  64.             B1=0.D0
  65.             B2=0.D0
  66.             B3=0.D0
  67.             C1=0.D0
  68.             C2=0.D0
  69.             C3=0.D0
  70.             DO 75 I=1,NOEL1
  71.                  A1=A1+SHP(2,I)*XNOEU(1,I)
  72.                  A2=A2+SHP(3,I)*XNOEU(1,I)
  73.                  A3=A3+SHP(4,I)*XNOEU(1,I)
  74.                  B1=B1+SHP(2,I)*XNOEU(2,I)
  75.                  B2=B2+SHP(3,I)*XNOEU(2,I)
  76.                  B3=B3+SHP(4,I)*XNOEU(2,I)
  77.                  C1=C1+SHP(2,I)*XNOEU(3,I)
  78.                  C2=C2+SHP(3,I)*XNOEU(3,I)
  79.                  C3=C3+SHP(4,I)*XNOEU(3,I)
  80.   75        CONTINUE
  81.             DET=A1*(B2*C3-C2*B3)+A2*(C1*B3-B1*C3)+A3*(B1*C2-B2*C1)
  82. C
  83.       ENDIF
  84.       RETURN
  85.       END
  86.  
  87.  
  88.  
  89.  
  90.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales