Télécharger la plaquette Cast3M


fr en

Télécharger capa4.eso

Retour à la liste

Numérotation des lignes : 

capa4

  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