Télécharger la plaquette Cast3M


fr en

Télécharger trjvlo.eso

Retour à la liste

Numérotation des lignes : 

trjvlo

  1. C TRJVLO    SOURCE    CHAT      05/01/13    03:51:59     5004
  2.       SUBROUTINE TRJVLO(IZNOEU,IZUN,IZTRAV,ITYP,UELEM,NUPOIN,IEL,IZSH)
  3. C
  4. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  5. C
  6. C     CALCUL DE LA VITESSE (REELLE) UELEM A L'INTERIEUR D'UN ELEMENT
  7. C     EN FONCTION DES VITESSES NODALES
  8. C
  9. C     IZNOEU COORDONNEES ET NUMEROS DES NOEUDS DE L'ELEMENT
  10. C     IZUN VITESSES AUX NOEUDS DE L'ELEMENT
  11. C     IZTRAV COORDONNEES DES NOEUDS (ON S'INTERESSE AU NUPOIN-IEME )
  12. C     ITYP TYPE DE L'ELEMENT
  13. C     NUPOIN NUMERO DE LA PARTICULE A LAQUELLE ON S'INTERESSE
  14. C     IEL NUMERO DE L ELEMENT
  15. C     ITYP TYPE DE LA VITESSE
  16. C     IZSH SEGMENT DE TRAVAIL POUR LES FONCTIONS DE FORME
  17. C
  18. C
  19. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  20. C
  21. C
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8 (A-H,O-Z)
  24. C
  25.  
  26.       SEGMENT  IZNOEU
  27.             REAL*8 XELE(IDIM,NOEL)
  28.             INTEGER NOEGLO(NOEL)
  29.       ENDSEGMENT
  30.       SEGMENT IZTRAV
  31.            REAL*8  COOR(NDIM,NPART)
  32.       ENDSEGMENT
  33. C
  34.       SEGMENT IZSH
  35.             REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  36.       ENDSEGMENT
  37. C
  38.       SEGMENT IZUN
  39.            REAL*8 UN(I1,I2,I3)
  40.       ENDSEGMENT
  41. C
  42.       DIMENSION UELEM(3),XYREF(3),XYZP(3)
  43. C
  44.       IDIM=XELE(/1)
  45.       NOEL=XELE(/2)
  46. C
  47.       DO 1 ID=1,IDIM
  48.         XYZP(ID)=COOR(ID,NUPOIN)
  49.    1  CONTINUE
  50. C
  51. C
  52.        IER=0
  53. C                     TRI3
  54.         IF(ITYP.EQ.4)CALL TRJ302(XELE,XYZP,XYREF)
  55. C                     QUA4
  56.         IF(ITYP.EQ.8)CALL TRJ402(XELE,XYZP,XYREF,IER)
  57. C                     TRI6
  58.         IF(ITYP.EQ.6)CALL TRJ602(XELE,XYZP,XYREF,IER)
  59. C                     TRI7
  60.         IF(ITYP.EQ.7)CALL TRJ702(XELE,XYZP,XYREF,IER)
  61. C                     QUA9
  62.         IF(ITYP.EQ.11)CALL TRJ902(XELE,XYZP,XYREF,IER)
  63. C                     PRI6
  64.         IF(ITYP.EQ.16)CALL TRJ603(XELE,XYZP,XYREF,IER)
  65. C                     CUB8
  66.         IF(ITYP.EQ.14)CALL TRJ803(XELE,XYZP,XYREF,IER)
  67. C                     TET4
  68.         IF(ITYP.EQ.23)CALL TRJ403(XELE,XYZP,XYREF)
  69.       IF(IER.NE.0)THEN
  70.                    CALL ERREUR(IER)
  71.       ENDIF
  72. C
  73. C
  74.         CALL TRJVPO(UN(1,1,IEL),UELEM,XYREF,ITYP,IDIM,NOEL,SHP)
  75.       RETURN
  76.       END
  77.  
  78.  
  79.  
  80.  

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