Télécharger reeref.eso

Retour à la liste

Numérotation des lignes :

  1. C REEREF SOURCE CHAT 05/01/13 02:47:56 5004
  2. SUBROUTINE REEREF(NDIM,ITY1,NOEL1,IZSH,XREEL,XREF)
  3. *******************************************************************
  4. *** SP 'REEREF' : permet passage entre coord reelles d'un pt de
  5. *** l'elemt reel et coord reference de ce pt dans l'elemt reference
  6. ***
  7. *** APPELES 1 = 'TRJ302', 'TRJ402', 'TRJ803', 'TRJ603', 'TRJ403'
  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 les fctions de forme, base ainsi
  14. *** que coord reelles des noeuds de l'element considere
  15. *** 'XREEL' coordonnees reelles du pt considere
  16. ***
  17. *** S = 'XREF' coordonnees de reference du pt considere
  18. ***
  19. *** RQ : les sp appeles se basent sur des methodes iteratives avec
  20. *** un critere de convergence d'ou une erreur de precision possible.
  21. *********************************************************************
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8 (A-H,O-Z)
  24. SEGMENT IZSH
  25. REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  26. ENDSEGMENT
  27. DIMENSION XREEL(3),XREF(3)
  28. DIMENSION XNOE23(2,3),XNOE24(2,4)
  29. DIMENSION XNOE34(3,4),XNOE36(3,6),XNOE38(3,8)
  30. *** cas TRI3 (triangles)
  31. IF (ITY1.EQ.4) THEN
  32. DO 1 I=1,NDIM
  33. DO 2 K=1,NOEL1
  34. *** 'XNOE23' = format imposé par le sp 'TRJ302'
  35. XNOE23(I,K)=XYZL(I,K)
  36. 2 CONTINUE
  37. 1 CONTINUE
  38. *** sp 'TRJ302' calcule le 'XREF' associé au 'XREEL'
  39. CALL TRJ302(XNOE23,XREEL,XREF)
  40. ENDIF
  41. *** cas QUA4 (quadrangles)
  42. IF (ITY1.EQ.8) THEN
  43. DO 3 I=1,NDIM
  44. DO 4 K=1,NOEL1
  45. XNOE24(I,K)=XYZL(I,K)
  46. 4 CONTINUE
  47. 3 CONTINUE
  48. CALL TRJ402(XNOE24,XREEL,XREF,IER)
  49. ENDIF
  50. *** cas CUB8 (cubes)
  51. IF (ITY1.EQ.14) THEN
  52. DO 5 I=1,NDIM
  53. DO 6 K=1,NOEL1
  54. XNOE38(I,K)=XYZL(I,K)
  55. 6 CONTINUE
  56. 5 CONTINUE
  57. CALL TRJ803(XNOE38,XREEL,XREF,IER)
  58. ENDIF
  59. *** cas PRI6 (prismes)
  60. IF (ITY1.EQ.16) THEN
  61. DO 7 I=1,NDIM
  62. DO 8 K=1,NOEL1
  63. XNOE36(I,K)=XYZL(I,K)
  64. 8 CONTINUE
  65. 7 CONTINUE
  66. CALL TRJ603(XNOE36,XREEL,XREF,IER)
  67. ENDIF
  68. *** cas TET4 (tetraedres)
  69. IF (ITY1.EQ.23) THEN
  70. DO 9 I=1,NDIM
  71. DO 10 K=1,NOEL1
  72. XNOE34(I,K)=XYZL(I,K)
  73. 10 CONTINUE
  74. 9 CONTINUE
  75. CALL TRJ403(XNOE34,XREEL,XREF)
  76. ENDIF
  77. RETURN
  78. END
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  

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