Télécharger la plaquette Cast3M


fr en

Télécharger trjnoe.eso

Retour à la liste

Numérotation des lignes : 

trjnoe

  1. C TRJNOE    SOURCE    CHAT      05/01/13    03:50:16     5004
  2.       SUBROUTINE TRJNOE(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL1,IEL,
  3.      * INOELO,IZPART,IZUN,IZCOU,ITP,IFORML,XY)
  4. C
  5. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  6. C
  7. C     LORSQU UNE PARTICULE EST LACHEE SUR UN NOEUD DU MAILLAGE
  8. C     GEOMETRIQUE ON DETERMINE SI ELLE EST OU NON DANS L ELEMENT IEL1
  9. C
  10. C     IZUN VITESSE OU FLUX
  11. C     XY COORDONNEES DE REFERENCE DU POINT CONSIDERE
  12. C
  13. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  14. C
  15. C
  16.       IMPLICIT INTEGER(I-N)
  17.       IMPLICIT REAL*8 (A-H,O-Z)
  18. C
  19.  
  20. -INC PPARAM
  21. -INC CCOPTIO
  22. C
  23. C
  24.       SEGMENT IZPART
  25.            INTEGER NLEPA(NPART),NUMPA(NPART)
  26.            REAL*8  COORPA(NDIM,NPART)
  27.       ENDSEGMENT
  28.       POINTEUR IZREF.IZPART
  29.       SEGMENT  IZCOU
  30.            REAL*8 DTCO(NEL),COU
  31.       ENDSEGMENT
  32. C
  33.       SEGMENT IZTRAV
  34.            REAL*8  COOR(NDIM,NPART)
  35.       ENDSEGMENT
  36.       SEGMENT  IZNOEU
  37.             REAL*8 XELE(IDIM,NOEL)
  38.             INTEGER NOEGLO(NOEL)
  39.       ENDSEGMENT
  40.       SEGMENT  IZAPAR
  41.             INTEGER IAPAR(4,NPART)
  42.       ENDSEGMENT
  43.       SEGMENT IZUN
  44.            REAL*8 UN(I1,I2,I3)
  45.       ENDSEGMENT
  46.       SEGMENT HYTRAV
  47.             REAL*8 XEL(3,NOEL),SHP(6,NOEL),SHY(12,NOEL)
  48.             REAL*8 ZJAC(NDIM,NDIM)
  49.       ENDSEGMENT
  50. C
  51. C
  52.       DIMENSION XY(3),UELEM(3),VECT(3)
  53.       DIMENSION TLI(4)
  54. C
  55. C
  56.       NDIM=COORPA(/1)
  57.       COUR=0.001D0
  58.       INOEL1= INOELO
  59.       IF(IFORML.EQ.1)THEN
  60. C                        FORMULATION ELEMENTS FINIS
  61.       DO 10 ID=1,NDIM
  62.       UELEM(ID)=UN(ID,INOELO,IEL)
  63.    10 CONTINUE
  64.       ELSEIF(IFORML.EQ.2)THEN
  65. C                     FORMULATION HYBRIDE ON VA CALCULER LA VITESSE
  66. C                   REELLE AU NOEUD DU MAILLAGE
  67.  
  68.          CALL VHYNOE(UELEM,IZNOEU,IZUN,NDIM,XY,ITP,IEL)
  69.       ENDIF
  70.       DO 70 ID=1,NDIM
  71.       COOR(ID,IPART)=COORPA(ID,IPART)+
  72.      *              UELEM(ID)*DTCO(IEL1)*COUR
  73. C     write(6,*)' coor trjnoe ',COOR(ID,IPART),COORPA(ID,IPART),
  74. C    *              UELEM(ID),DTCO(IEL1),COUR
  75.   70  CONTINUE
  76.       CALL INITI(IAPAR(1,IPART),4,0)
  77. C***  TRIANGLES
  78. C
  79.             IF(ITP.EQ.4.OR.ITP.EQ.6.OR.ITP.EQ.7)THEN
  80.                 CALL TRJTRI(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL1,
  81.      *                    INOELO,TLI)
  82. C
  83. C***  QUADRANGLES
  84. C
  85.             ELSEIF(ITP.EQ.8.OR.ITP.EQ.11)THEN
  86.                 CALL TRJQUA(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL1,
  87.      *                    INOELO,TLI)
  88. C***  PRISMES
  89. C
  90.             ELSEIF(ITP.EQ.16)THEN
  91.                 CALL TRJPRI(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL,
  92.      *                                  INOELO,TLI,ITRI)
  93. C
  94. C***  CUBES
  95. C
  96.             ELSEIF(ITP.EQ.14)THEN
  97.                 CALL TRJCUB(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL,
  98.      *                                               INOELO,TLI)
  99. C
  100. C***  TETRAEDRE
  101. C
  102.             ELSEIF(ITP.EQ.23)THEN
  103.                 CALL TRJTET(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL,
  104.      *                                               INOELO,TLI)
  105.           ELSE
  106.              CALL ERREUR(16)
  107.           ENDIF
  108.       DO 80 ID=1,NDIM
  109.             COOR(ID,IPART)=COORPA(ID,IPART)
  110.   80  CONTINUE
  111.       INOELO=INOEL1
  112. C
  113.       END
  114.  
  115.  
  116.  

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