C REFREE SOURCE CHAT 05/01/13 02:48:04 5004 SUBROUTINE REFREE(NDIM,ITY1,NOEL1,IZSH,XREF,XREEL) *********************************************************************** *** SP 'REFREE' : permet le passage entre coord reference d'un pt de *** l'element de reference et coord reelles de ce pt dans l'element reel *** *** APPELES 1 = 'SHAPE' *** APPELES 2 = aucun *** *** E = 'NDIM' dimension de l'espace *** 'ITY1' entier caracterisant le type de l'element *** 'NOEL1' nombre de noeuds de l'element considere *** 'IZSH' segment contenant fonctions de forme et de base ainsi *** que coord reelles des noeuds de l'element considere *** 'XREF' coordonnees reference du pt considere *** *** S = 'XREEL' coordonnees reelles du pt considere *********************************************************************** IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) SEGMENT IZSH REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9) ENDSEGMENT DIMENSION XREF(3),XREEL(3) *** calcul des fonctions de forme au pt de calcul de l'element reference CALL SHAPE(XREF(1),XREF(2),XREF(3),ITY1,SHP,IRET) *** passage en coord reelles via les fonctions de forme et valeurs aux noeuds DO 10 I=1,NDIM XREEL(I)=0.D0 DO 20 K=1,NOEL1 XREEL(I)=XREEL(I)+SHP(1,K)*XYZL(I,K) 20 CONTINUE 10 CONTINUE RETURN END