Télécharger capa4.eso

Retour à la liste

Numérotation des lignes :

  1. C CAPA4 SOURCE FANDEUR 11/07/19 21:15:20 7042
  2.  
  3. C=======================================================================
  4. C= C A P A 4 =
  5. C= --------- =
  6. C= =
  7. C= Fonction : =
  8. C= ---------- =
  9. C= Calcul, au point de Gauss (iGau) considere, du jacobien (DJAC), =
  10. C= des fonctions de forme (FORME) et de leurs derivees (SHP) =
  11. C= =
  12. C= Parametres : (E)=Entree (S)=Sortie =
  13. C= ------------ =
  14. C= iGau (E) Numero du point de Gauss considere =
  15. C= NBNN (E) Nombre de NOEUDS de l'element =
  16. C= XE (E) COORDONNEES des noeuds de l'element (repere GLOBAL) =
  17. C= SHPTOT (E) Fonctions de FORME et leurs DERIVEES par rapport aux =
  18. C= coordonnees de REFERENCE (Qsi,Eta,Dzeta) =
  19. C= SHP (S) Fonctions de FORME et leurs DERIVEES par rapport aux =
  20. C= coordonnees REELLES (X,Y,Z) =
  21. C= FORME (S) Fonctions de FORME dans la geometrie REELLE =
  22. C= DJAC (S) JACOBIEN au point de Gauss considere =
  23. C= =
  24. C= Denis ROBERT, le 17 fevrier 1988. =
  25. C=======================================================================
  26.  
  27. SUBROUTINE CAPA4 (NEF,iGau,NBNN,XE,SHPTOT,SHP,FORME,DJAC)
  28.  
  29. IMPLICIT INTEGER(I-N)
  30. IMPLICIT REAL*8 (A-H,O-Z)
  31.  
  32.  
  33. -INC PPARAM
  34. -INC CCOPTIO
  35. -INC CCREEL
  36.  
  37. C= Quelques constantes (2.Pi et 4.Pi)
  38. PARAMETER (X2Pi= 6.283185307179586476925286766559D0)
  39. PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  40.  
  41. DIMENSION XE(3,*),FORME(*),SHPTOT(6,NBNN,*),SHP(6,*)
  42.  
  43. C ====================
  44. C 1 - Initialisation
  45. C ====================
  46. C*OF CALL ZERO(FORME,1,NBNN)
  47.  
  48. C ============================================
  49. C 2 - Cas des elements MASSIFS 2D PLAN ou 3D
  50. C ============================================
  51. IF (IFOMOD.EQ.-1.OR.IFOMOD.EQ.2) THEN
  52. iFin=IDIM+1
  53. DO j=1,NBNN
  54. FORME(j)=SHPTOT(1,j,iGau)
  55. SHP(1,j)=FORME(j)
  56. DO i=2,iFin
  57. SHP(i,j)=SHPTOT(i,j,iGau)
  58. ENDDO
  59. ENDDO
  60. CALL JACOBI(XE,SHP,IDIM,NBNN,DJAC)
  61.  
  62. C ================================================
  63. C 3 - Cas des elements MASSIFS 2D AXISYMETRIQUES
  64. C ================================================
  65. ELSE IF (IFOMOD.EQ.0) THEN
  66. DO i=1,NBNN
  67. FORME(i)=SHPTOT(1,i,iGau)
  68. SHP(1,i)=FORME(i)
  69. SHP(2,i)=SHPTOT(2,i,iGau)
  70. SHP(3,i)=SHPTOT(3,i,iGau)
  71. ENDDO
  72. CALL JACOBI(XE,SHP,IDIM,NBNN,DJAC)
  73. CALL DISTRR(XE,SHP,NBNN,RR)
  74. DJAC=X2Pi*DJAC*RR
  75.  
  76. C =================================
  77. C 4 - Cas des elements MASSIFS 1D
  78. C =================================
  79. ELSE IF (IFOMOD.EQ.3.OR.IFOMOD.EQ.4.OR.IFOMOD.EQ.5) THEN
  80. C =====
  81. C 4.1 - Cas generique
  82. C =====
  83. DO i=1,NBNN
  84. FORME(i)=SHPTOT(1,i,iGau)
  85. SHP(1,i)=FORME(i)
  86. SHP(2,i)=SHPTOT(2,i,iGau)
  87. ENDDO
  88. CALL JACOBI(XE,SHP,IDIM,NBNN,DJAC)
  89. C =====
  90. C 4.2 - Cas particulier des modes axisymetriques et spherique 1D
  91. C =====
  92. IF (IFOMOD.EQ.4 .OR. IFOMOD.EQ.5) THEN
  93. CALL DISTRR(XE,SHP,NBNN,RR)
  94. IF (IFOMOD.EQ.5) THEN
  95. DJAC=X4Pi*DJAC*RR*RR
  96. ELSE
  97. DJAC=X2Pi*DJAC*RR
  98. ENDIF
  99. ENDIF
  100.  
  101. C ==============================================
  102. C 5 - Cas des elements MASSIFS en mode FOURIER -> ERREUR
  103. C ==============================================
  104. ELSE IF (IFOMOD.EQ.1) THEN
  105. CALL ERREUR(251)
  106. ENDIF
  107.  
  108. RETURN
  109. END
  110.  
  111.  
  112.  

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