Télécharger nelec.eso

Retour à la liste

Numérotation des lignes :

  1. C NELEC SOURCE FANDEUR 10/12/17 21:18:18 6427
  2.  
  3. C=======================================================================
  4. C= N E L E C =
  5. C= --------- =
  6. C= Fonction : =
  7. C= ---------- =
  8. C= Calcul de la matrice N des fonctions de forme de l'element fini =
  9. C= en un point donne de l'element, matrice N intervenant lors du =
  10. C= calcul des forces nodales equivalentes a une densite de sources =
  11. C= volumiques. =
  12. C= =
  13. C= Parametres : (E)=Entree (S)=Sortie =
  14. C= ------------ =
  15. C= NBNO (E) Nombre de noeuds de l'element fini =
  16. C= XEL (E) Coordonnees des noeuds de l'element fini etudie =
  17. C= SHPREF (E) Fonctions de forme et leurs derivees (reference) =
  18. C= SHP (S) Fonctions de forme et leurs derivees (reelles) =
  19. C= FORME (S) Matrice des fcts de forme (N) au point considere =
  20. C= DJAC (S) Jacobien au point considere =
  21. C= =
  22. C= Remarque : =
  23. C= ---------- =
  24. C= Lors de l'entree dans le sous-programme, SHPREF(2 a 4,*) contient =
  25. C= les DERIVEES des fonctions de forme par rapport aux coordonnees =
  26. C= de REFERENCE Qsi,Eta,Dzeta. =
  27. C= En sortie du sous-programme, SHP(2 a 4,*) contient les DERIVEES =
  28. C= des fonctions de FORME par rapport aux coordonnees REELLES x,y,z. =
  29. C=======================================================================
  30.  
  31. SUBROUTINE NELEC (NBNO,XEL,SHPREF, SHP,FORME,DJAC)
  32.  
  33. IMPLICIT INTEGER(I-N)
  34. IMPLICIT REAL*8 (A-H,O-Z)
  35.  
  36. -INC CCOPTIO
  37. -INC CCREEL
  38.  
  39. C= Quelques constantes (2.Pi et 4.Pi)
  40. PARAMETER (X2Pi= 6.283185307179586476925286766559D0)
  41. PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  42.  
  43. DIMENSION XEL(3,*),SHPREF(6,*), SHP(6,*),FORME(*)
  44.  
  45. C ====================
  46. C 1 - Initialisation
  47. C ====================
  48. C* DO j = 1, NBNO
  49. C* FORME(j) = XZero
  50. C* ENDDO
  51.  
  52. C ===================================
  53. C 2.1 - Cas des elements MASSIFS 3D
  54. C ===================================
  55. IF (IFOMOD.EQ.2) THEN
  56. DO j = 1, NBNO
  57. SHP(1,j) = SHPREF(1,j)
  58. SHP(2,j) = SHPREF(2,j)
  59. SHP(3,j) = SHPREF(3,j)
  60. SHP(4,j) = SHPREF(4,j)
  61. FORME(j) = SHPREF(1,j)
  62. ENDDO
  63. CALL JACOBI(XEL,SHP,IDIM,NBNO,DJAC)
  64. C =========================================================
  65. C 2.2 - Cas des elements MASSIFS 2D PLAN ou AXISYMETRIQUE
  66. C =========================================================
  67. ELSE IF (IFOMOD.EQ.-1 .OR. IFOMOD.EQ.0) THEN
  68. DO j = 1, NBNO
  69. SHP(1,j) = SHPREF(1,j)
  70. SHP(2,j) = SHPREF(2,j)
  71. SHP(3,j) = SHPREF(3,j)
  72. FORME(j) = SHPREF(1,j)
  73. ENDDO
  74. CALL JACOBI(XEL,SHP,IDIM,NBNO,DJAC)
  75. IF (IFOMOD.EQ.0) THEN
  76. CALL DISTRR(XEL,SHP,NBNO,RR)
  77. DJAC = X2Pi*DJAC*RR
  78. ENDIF
  79. C ===================================
  80. C 2.3 - Cas des elements MASSIFS 1D
  81. C ===================================
  82. ELSE IF (IFOMOD.EQ.3 .OR. IFOMOD.EQ.4 .OR. IFOMOD.EQ.5) THEN
  83. DO j = 1, NBNO
  84. SHP(1,j) = SHPREF(1,j)
  85. SHP(2,j) = SHPREF(2,j)
  86. FORME(j) = SHPREF(1,j)
  87. ENDDO
  88. CALL JACOBI(XEL,SHP,IDIM,NBNO,DJAC)
  89. C Cas particulier des modes axisymetriques et spherique 1D
  90. IF (IFOMOD.EQ.4.OR.IFOMOD.EQ.5) THEN
  91. CALL DISTRR(XEL,SHP,NBNO,RR)
  92. IF (IFOMOD.EQ.5) THEN
  93. DJAC = X4Pi*DJAC*RR*RR
  94. ELSE
  95. DJAC = X2Pi*DJAC*RR
  96. ENDIF
  97. ENDIF
  98. C ===========================================
  99. C 2.4 - Cas des elements MASSIFS 2D FOURIER
  100. C ===========================================
  101. ELSE IF (IFOMOD.EQ.1) THEN
  102. DO j = 1,NBNO
  103. SHP(1,j) = SHPREF(1,j)
  104. SHP(2,j) = SHPREF(2,j)
  105. SHP(3,j) = SHPREF(3,j)
  106. FORME(j) = SHPREF(1,j)
  107. ENDDO
  108. CALL JACOBI(XEL,SHP,IDIM,NBNO,DJAC)
  109. CALL DISTRR(XEL,SHP,NBNO,RR)
  110. IF (NIFOUR.EQ.0) THEN
  111. DJAC = X2Pi*DJAC*RR
  112. ELSE
  113. DJAC = XPi*DJAC*RR
  114. ENDIF
  115. ENDIF
  116.  
  117. RETURN
  118. END
  119.  
  120.  
  121.  

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