Télécharger la plaquette Cast3M


fr en

Télécharger trjvpo.eso

Retour à la liste

Numérotation des lignes : 

trjvpo

  1. C TRJVPO    SOURCE    CHAT      05/01/13    03:52:04     5004
  2.       SUBROUTINE TRJVPO(UNOEU,UE,XY,ITYP,IDIM,NOEL1,SHP)
  3. C*******************************************************************
  4. C      CALCUL DE LA VITESSE AU POINT XY
  5. C
  6. C********************************************************************
  7. C   ENTREES
  8. C     UNOEU VITESSES AUX NOEUDS DE L ELEMENT
  9. C     XY COORDONNEES DE REFERENCES DU POINT DE CALCUL
  10. C     ITYP TYPE DE L ELEMENT
  11. C     NOEL1  NOMBRE DE NOEUDS DE L ELEMENT POSITION
  12. C   SORTIES
  13. C     UE VITESSE DANS L ELEMENT DE REFERENCE POUR LE POINT
  14. C        DE COORDONNEES XY
  15. C
  16. C     SHP(1,I) FONCTIONS DE FORMES AU POINT DE CALCUL
  17. C     SHP(2,I) DERIVEE DES FONCTIONS DE FORME / XSI AU POINT DE CALCUL
  18. C     SHP(3,I)    "     "      "     "    "   / ETA "    "   "    "
  19. C     SHP(4,I)    "     "      "      "   "   /DZETA "   "    "   "
  20. C ********************************************************************
  21. C
  22. C
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8 (A-H,O-Z)
  25. C
  26. C
  27. C
  28.       DIMENSION UNOEU(*),XY(*),UE(*)
  29.       DIMENSION SHP(6,*)
  30. C
  31. C
  32. C***  CALCUL DES XI ET DES DXI
  33. C
  34.       CALL SHAPE(XY(1),XY(2),XY(3),ITYP,SHP,IRET)
  35. C
  36. C***  VITESSE
  37. C
  38.       DO 51 I=1,IDIM
  39.       UE(I)=0.D0
  40.       DO 52 K=1,NOEL1
  41.       J=(K-1)*IDIM+I
  42.       UE(I)=UE(I)+SHP(1,K)*UNOEU(J)
  43.   52  CONTINUE
  44.   51  CONTINUE
  45. C     write(6,*) 'vitesse',(ue(i),i=1,idim)
  46. C
  47.   50  CONTINUE
  48. C
  49.       RETURN
  50.       END
  51.  
  52.  
  53.  

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