Télécharger detjac.eso

Retour à la liste

Numérotation des lignes :

  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