Télécharger refree.eso

Retour à la liste

Numérotation des lignes :

  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