Télécharger pourea.eso

Retour à la liste

Numérotation des lignes :

  1. C POUREA SOURCE CHAT 05/01/13 02:19:34 5004
  2. SUBROUTINE POUREA(XX,YY,ZZ,XDDL,WORK,P,VECT,KERRE)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE REACTUALISATION DU VECTEUR VECT
  5. C-----------------------------------------------------------------------
  6. C ENTREE
  7. C XX,YY,ZZ LES COORDONNEES DES 2 NOEUDS
  8. C XDDL LES DEPLACEMENTS EN AXES GLOBAUX
  9. C WORK TABLEAU DE TRAVAIL
  10. C P MATRICE DE PASSAGE
  11. C
  12. C SORTIE
  13. C VECT NOUVEAU VECTEUR
  14. C KERRE INDICE D'ERREUR ( 0 SI TOUT EST OK )
  15. C
  16. C-----------------------------------------------------------------------
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19. DIMENSION WORK(*),P(3,*),VECT(*),VN(3),VNP(3),XDDL(*)
  20. DIMENSION OMEG(3),PP(3,3),XX(*),YY(*),ZZ(*)
  21. C
  22. C CALCUL DE VNP
  23. C
  24. VNP(1)=XX(2)+XDDL(7)-XX(1)-XDDL(1)
  25. VNP(2)=YY(2)+XDDL(8)-YY(1)-XDDL(2)
  26. VNP(3)=ZZ(2)+XDDL(9)-ZZ(1)-XDDL(3)
  27. RLVNP=SQRT(VNP(1)*VNP(1)+VNP(2)*VNP(2)+VNP(3)*VNP(3))
  28. IF(RLVNP.EQ.0.D0) THEN
  29. KERRE=1
  30. RETURN
  31. ENDIF
  32. DO 1 I=1,3
  33. VNP(I)=VNP(I)/RLVNP
  34. 1 OMEG(I)=VNP(I)
  35. C
  36. C CALCUL DE LA MATRICE DE PASSAGE
  37. C
  38. CALL POUPAS(XX,YY,ZZ,OMEG,PP,KERRE)
  39. IOMEG=KERRE
  40. KERRE=0
  41. C
  42. C PASSAGE DES DEPLACEMENTS EN AXES LOCAUX - ROTATION PSI
  43. C
  44. CALL POUVEC(XDDL,WORK,P,1)
  45. PSI=(WORK(4)+WORK(10))*0.5D0
  46. C2=PP(1,1)*VNP(1)+PP(1,2)*VNP(2)+PP(1,3)*VNP(3)
  47. S2=PP(2,1)*VNP(1)+PP(2,2)*VNP(2)+PP(2,3)*VNP(3)
  48. C3=COS(PSI)
  49. S3=SIN(PSI)
  50. C
  51. C REACTUALISATION DU VECTEUR
  52. C
  53. IF(IOMEG.EQ.0) THEN
  54. DO 2 I=1,3
  55. VN(I)=0.D0
  56. DO 2 J=1,3
  57. 2 VN(I)=VN(I)+PP(I,J)*VECT(J)
  58. OMEG(1)=C2*VN(1)-S2*C3*VN(2)+S2*S3*VN(3)
  59. OMEG(2)=S2*VN(1)+C2*C3*VN(2)-C2*S3*VN(3)
  60. OMEG(3)=S3*VN(2)+C3*VN(3)
  61. DO 3 I=1,3
  62. VECT(I)=0.D0
  63. DO 3 J=1,3
  64. 3 VECT(I)=VECT(I)+PP(J,I)*OMEG(J)
  65. ELSE
  66. VN(1)=0.D0
  67. VN(2)=C3
  68. VN(3)=S3
  69. DO 6 I=1,3
  70. VECT(I)=0.D0
  71. DO 6 J=1,3
  72. 6 VECT(I)=VECT(I)+P(J,I)*VN(J)
  73. ENDIF
  74. RETURN
  75. END
  76.  
  77.  

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