Télécharger la plaquette Cast3M


fr en

Télécharger bdiff1.eso

Retour à la liste

Numérotation des lignes : 

bdiff1

  1. C BDIFF1    SOURCE    PV090527  23/01/12    21:15:03     11551          
  2.  
  3. C=======================================================================
  4. C=                          B D I F F 1                                =
  5. C=                          -----------                                =
  6. C=  Fonction :                                                         =
  7. C=  ----------                                                         =
  8. C=   Calcul de la matrice Bdif reliant le gradient (au signe pres) en  =
  9. C=   un point d'un element fini aux ddls de la quantite "diffusante"   =
  10. C=                                                                     =
  11. C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  12. C=  ------------                                                       =
  13. C=   XEL     (E)  Coordonnees des noeuds de l'element fini considere   =
  14. C=   SHPREF  (E)  Fonctions de forme et leurs derivees au point etudie =
  15. C=                dans l'element fini de reference                     =
  16. C=   NBNO    (E)  Nombre de noeuds de l'element fini                   =
  17. C=   NDIM    (E)  Nombre de lignes de la matrice GRDDIF                =
  18. C=   ISGNG   (E)  Valeur negative (<0) si on veut l'oppose du gradient =
  19. C=                                                                     =
  20. C=   SHP     (S)  Fonctions de forme et leurs derivees "actuelles"     =
  21. C=   GRDDIF  (S)  Matrice de gradients Bdif calculee au point etudie   =
  22. C=                                                                     =
  23. C=  Remarque :                                                         =
  24. C=  ----------                                                         =
  25. C=   Lors de l'entree dans le sous-programme, SHPREF(2 a 4,*) contient =
  26. C=   les DERIVEES des fonctions de forme par rapport aux coordonnees   =
  27. C=   de REFERENCE Qsi,Eta,Dzeta.                                       =
  28. C=   En sortie du sous-programme, SHP(2 a 4,*) contient les DERIVEES   =
  29. C=   des fonctions de FORME par rapport aux coordonnees REELLES x,y,z. =
  30. C=======================================================================
  31.  
  32.       SUBROUTINE BDIFF1(XEL,SHPREF,NBNO,NDIM,ISGNG, SHP,GRDDIF)
  33.  
  34.       IMPLICIT INTEGER(I-N)
  35.       IMPLICIT REAL*8 (A-H,O-Z)
  36.  
  37. -INC CCREEL
  38.  
  39. -INC PPARAM
  40. -INC CCOPTIO
  41.  
  42. C= Quelques constantes (2.Pi et 4.Pi)
  43.       PARAMETER (X2Pi= 6.283185307179586476925286766559D0)
  44.       PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  45.  
  46.       DIMENSION XEL(3,*),SHPREF(6,*), SHP(6,*),GRDDIF(NDIM,*)
  47.  
  48. C ====================
  49. C  1 - Initialisation
  50. C ====================
  51.  
  52. C ===========================================
  53. C  2.4 - Cas des elements BARR, TUY2 et TUY3
  54. C ===========================================
  55.       JDIM = 1
  56.       CALL GTEMRD (XEL,SHPREF,JDIM,NBNO,DJAC)
  57.  
  58.       DO j = 1,NBNO
  59. C       SHP(1,j) = SHPREF(1,j)
  60.         SHP(2,j) = SHPREF(2,j) / DJAC
  61.       ENDDO
  62.  
  63. C=    Traitement des modes axisymetriques et spherique 1D
  64.       IF (ISGNG .GE. 0) THEN
  65.         DO j = 1,NBNO
  66.           GRDDIF(1,j) = SHP(2,j)
  67.         ENDDO
  68.       ELSE
  69.         DO j = 1, NBNO
  70.           GRDDIF(1,j) = -SHP(2,j)
  71.         ENDDO
  72.       ENDIF
  73.  
  74.       RETURN
  75.       END
  76.  
  77.  

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