Télécharger la plaquette Cast3M


fr en

Télécharger refree.eso

Retour à la liste

Numérotation des lignes : 

refree

  1. C REFREE    SOURCE    CHAT      05/01/13    02:48:04     5004
  2.       SUBROUTINE REFREE(NDIM,ITY1,NOEL1,IZSH,XREF,XREEL)
  3. ***********************************************************************
  4. *** SP 'REFREE' : permet le passage entre coord reference d'un pt de
  5. *** l'element de reference et coord reelles de ce pt dans l'element reel
  6. ***
  7. *** APPELES 1 = 'SHAPE'
  8. *** APPELES 2 = aucun
  9. ***
  10. *** E = 'NDIM' dimension de l'espace
  11. ***     'ITY1' entier caracterisant le type de l'element
  12. ***     'NOEL1' nombre de noeuds de l'element considere
  13. ***     'IZSH' segment contenant fonctions de forme et de base ainsi
  14. ***            que coord reelles des noeuds de l'element considere
  15. ***     'XREF' coordonnees reference du pt considere
  16. ***
  17. *** S = 'XREEL' coordonnees reelles du pt considere
  18. ***********************************************************************
  19.  
  20.       IMPLICIT INTEGER(I-N)
  21.       IMPLICIT REAL*8 (A-H,O-Z)
  22.       SEGMENT IZSH
  23.           REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  24.       ENDSEGMENT
  25.       DIMENSION XREF(3),XREEL(3)
  26. ***   calcul des fonctions de forme au pt de calcul de l'element reference
  27.       CALL SHAPE(XREF(1),XREF(2),XREF(3),ITY1,SHP,IRET)
  28. ***   passage en coord reelles via les fonctions de forme et valeurs aux noeuds
  29.       DO 10 I=1,NDIM
  30.          XREEL(I)=0.D0
  31.          DO 20 K=1,NOEL1
  32.             XREEL(I)=XREEL(I)+SHP(1,K)*XYZL(I,K)
  33.  20      CONTINUE
  34.  10   CONTINUE
  35.       RETURN
  36.       END
  37.  
  38.  
  39.  

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