Télécharger la plaquette Cast3M


fr en

Télécharger arete2.eso

Retour à la liste

Numérotation des lignes : 

arete2

  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