Télécharger tconv4.eso

Retour à la liste

Numérotation des lignes :

  1. C TCONV4 SOURCE CHAT 05/01/13 03:34:29 5004
  2.  
  3. C=======================================================================
  4. C= T C O N V 4 =
  5. C= ----------- =
  6. C= =
  7. C= Fonction : =
  8. C= ---------- =
  9. C= Calcul du jacobien et des derivees des fonctions de forme dans la =
  10. C= geometrie reelle pour les elements SEG2, SEG3, LIA3 et LIA4. =
  11. C= =
  12. C= Parametres : (E)=Entree (S)=Sortie =
  13. C= ------------ =
  14. C= XEL (E) COORDONNEES des noeuds de l'element (repere GLOBAL) =
  15. C= SHP (E/S) DERIVEES des fonctions de FORME de l'element =
  16. C= IDIM (E) DIMENSION du probleme =
  17. C= NBNN (E) Nombre de NOEUDS de l'element fini considere =
  18. C= DJAC (S) JACOBIEN =
  19. C= =
  20. C= Denis ROBERT, le 3 mai 1988. =
  21. C=======================================================================
  22.  
  23. SUBROUTINE TCONV4 (XEL,SHP,IDIM,NBNN,DJAC)
  24.  
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8 (A-H,O-Z)
  27.  
  28. -INC CCREEL
  29.  
  30. DIMENSION XEL(3,*),SHP(6,*)
  31.  
  32. C ===========================
  33. C 1 - Cas de la DIMENSION 2
  34. C ===========================
  35. IF (IDIM.EQ.2) THEN
  36. A1=XZERO
  37. A2=XZERO
  38. DO i=1,NBNN
  39. A1=A1+SHP(2,i)*XEL(1,i)
  40. A2=A2+SHP(2,i)*XEL(2,i)
  41. ENDDO
  42. DJAC=SQRT(A1*A1+A2*A2)
  43. IF (NBNN.EQ.2) THEN
  44. A1=XEL(2,2)-XEL(2,1)
  45. A2=XEL(1,2)-XEL(1,1)
  46. DLON=SQRT(A1*A1+A2*A2)
  47. COSANG=A2/DLON
  48. SINANG=A1/DLON
  49. C** Ces 4 lignes restent a verifier
  50. SHP(3,1)=SHP(2,1)*SINANG
  51. SHP(2,1)=SHP(2,1)*COSANG
  52. SHP(3,2)=SHP(2,2)*SINANG
  53. SHP(2,2)=SHP(2,2)*COSANG
  54. ENDIF
  55. C ===========================
  56. C 2 - Cas de la DIMENSION 3
  57. C ===========================
  58. ELSE IF (IDIM.EQ.3) THEN
  59. A1=XZERO
  60. A2=XZERO
  61. A3=XZERO
  62. DO i=1,NBNN
  63. A1=A1+SHP(2,i)*XEL(1,i)
  64. A2=A2+SHP(2,i)*XEL(2,i)
  65. A3=A3+SHP(2,i)*XEL(3,i)
  66. ENDDO
  67. DJAC=SQRT(A1*A1+A2*A2+A3*A3)
  68. ENDIF
  69.  
  70. RETURN
  71. END
  72.  
  73.  
  74.  
  75.  

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