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    MB234859  26/06/04    21:15:06     12564          
  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. -INC PPARAM
  39. -INC CCOPTIO
  40.  
  41. C= Quelques constantes (2.Pi et 4.Pi)
  42.       PARAMETER (X2Pi= 6.283185307179586476925286766559D0)
  43.       PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  44.  
  45.       DIMENSION XEL(3,*),SHPREF(6,*), SHP(6,*),GRDDIF(NDIM,*)
  46.  
  47. C ====================
  48. C  1 - Initialisation
  49. C ====================
  50.  
  51. C ===========================================
  52. C  2.4 - Cas des elements BARR, TUY2 et TUY3
  53. C ===========================================
  54.       JDIM = 1
  55.       CALL GTEMRD (XEL,SHPREF,JDIM,NBNO,DJAC)
  56.  
  57.       DO j = 1,NBNO
  58. C       SHP(1,j) = SHPREF(1,j)
  59.         SHP(2,j) = SHPREF(2,j) / ABS(DJAC)
  60.       ENDDO
  61.  
  62. C=    Traitement des modes axisymetriques et spherique 1D
  63.       XFAC=1.D0
  64.       IF (ISGNG .LT. 0) XFAC=-1.D0
  65.       DO j = 1,NBNO
  66.         GRDDIF(1,j) = XFAC * SHP(2,j)
  67.       ENDDO
  68.  
  69.       RETURN
  70.       END
  71.  
  72.  

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