Télécharger arete2.eso

Retour à la liste

Numérotation des lignes :

  1. C ARETE2 SOURCE PV 09/01/08 21:15:13 6254
  2. SUBROUTINE ARETE2 (MELEME,TS,ISEG,TN,INOR,TI,ICPR,LCPR)
  3.  
  4. * APPELE PAR ARETE1
  5. * REMPLIT LES TABLES TSEG ET TNOR POUR LES CALCULS
  6. * D'ELIMINATION D'ARETES .
  7. *
  8. *-------------------------------------------------------------------
  9. *
  10. * ARBRE D'APPEL :
  11. * --------------
  12. * ARETE2 : PREPARATION DU MAILLAGE PAR CALCUL DES NORMALES
  13. * AUX FACETTES , ET CLASSEMENT DES SEGMENTS APRES
  14. * ECLATES EN SEG2
  15. *
  16. *-------------------------------------------------------------------
  17. *
  18. * VARIABLES PRINCIPALES :
  19. * ---------------------
  20. *
  21. * TSEG : SEGMENT 2 DIMENSIONS CONTENANT
  22. *
  23. * | A TRACER (= 0 NON, =1 OUI) |
  24. * | COULEUR |
  25. * | N› NOEUD MIN |
  26. * | N› NOEUD MAX |
  27. * | ENTREES DANS TSEG |
  28. *
  29. * TNOR : SEGMENT 2 DIMENSIONS CONTENANT LES NORMALES A COMPARER
  30. * POUR CHAQUE FACE ELEMENTAIRE
  31. *
  32. * ISEG : TAILLE REELLE DU SEGMENT TSEG
  33. * INOR : TAILLE REELLE DU SEGMENT TNOR
  34. *
  35. *--------------------------------------------------------------------
  36. * AUTEUR : J.BRUN (JUIN 90)
  37.  
  38. IMPLICIT INTEGER(I-N)
  39. IMPLICIT REAL*8 (A-H,O-Y)
  40. IMPLICIT LOGICAL (Z)
  41. SEGMENT TN
  42. REAL*8 TNOR (N1,3)
  43. ENDSEGMENT
  44. SEGMENT TS
  45. INTEGER TSEG (M1,M2)
  46. ENDSEGMENT
  47. SEGMENT TI
  48. INTEGER TINDIC(MI1,MI2,2)
  49. ENDSEGMENT
  50. SEGMENT ICPR(0)
  51. -INC SMELEME
  52. REAL*8 VN(3)
  53.  
  54. *---------------------------------------------------------------------
  55. *
  56. * ELIMINATION ELEMENT POI1 ,SEG2 ,SEG3
  57. *
  58. *---------------------------------------------------------------------
  59. 100 IF (ITYPEL.LE.3) GOTO 200
  60. *---------------------------------------------------------------------
  61. *
  62. * BOUCLE SUR LES ELEMENTS
  63. *
  64. *---------------------------------------------------------------------
  65. DO 2 NEL=1,NUM(/2)
  66. *
  67. * ELEMENT PLAN TRI3 OU TRI4
  68. *
  69. IF (ITYPEL.LT.5) THEN
  70. CALL NORFAC(VN,NUM(1,NEL),NUM(2,NEL),NUM(3,NEL))
  71. DO 3 I=1,3
  72.  
  73. IF (I.EQ.3) THEN
  74. N1=NUM(3,NEL)
  75. N2=NUM(1,NEL)
  76. N3=ICOLOR(NEL)
  77. CALL TRISEG(TS,ISEG,TN,INOR,TI,ICPR,LCPR,VN,N1,N2,N3)
  78. ELSE
  79. N1=NUM(I,NEL)
  80. N2=NUM(I+1,NEL)
  81. N3=ICOLOR(NEL)
  82. CALL TRISEG(TS,ISEG,TN,INOR,TI,ICPR,LCPR,VN,N1,N2,N3)
  83. ENDIF
  84. 3 CONTINUE
  85. ELSE
  86. *
  87. * ELEMENT POUVANT ETRE NON PLAN ( A REDECOUPER EN TRI3 )
  88. *
  89. NB=NUM(/1)
  90. * TRI7 -> TRI6
  91. IF (ITYPEL.EQ.7) NB=6
  92. * QUA5 -> QUA4
  93. IF (ITYPEL.EQ.9) NB=4
  94. * QUA9 -> QUA8
  95. IF (ITYPEL.EQ.11) NB=8
  96. DO 4 I=1,NB,2
  97. IF (I.EQ.1) THEN
  98. CALL NORFAC(VN,NUM(NB,NEL),NUM(1,NEL),NUM(2,NEL))
  99. N1=NUM(NB,NEL)
  100. N2=NUM(1,NEL)
  101. N3=ICOLOR(NEL)
  102. CALL TRISEG(TS,ISEG,TN,INOR,TI,ICPR,LCPR,VN,N1,N2,N3)
  103. N1=NUM(1,NEL)
  104. N2=NUM(2,NEL)
  105. N3=ICOLOR(NEL)
  106. CALL TRISEG(TS,ISEG,TN,INOR,TI,ICPR,LCPR,VN,N1,N2,N3)
  107. ELSE
  108. CALL NORFAC(VN,NUM(I-1,NEL),NUM(I,NEL),NUM(I+1,NEL))
  109. N1=NUM(I-1,NEL)
  110. N2=NUM(I,NEL)
  111. N3=ICOLOR(NEL)
  112. CALL TRISEG(TS,ISEG,TN,INOR,TI,ICPR,LCPR,VN,N1,N2,N3)
  113. N1=NUM(I,NEL)
  114. N2=NUM(I+1,NEL)
  115. N3=ICOLOR(NEL)
  116. CALL TRISEG(TS,ISEG,TN,INOR,TI,ICPR,LCPR,VN,N1,N2,N3)
  117. ENDIF
  118. 4 CONTINUE
  119. ENDIF
  120. 2 CONTINUE
  121. 200 RETURN
  122. END
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  

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