Télécharger pourea.eso

Retour à la liste

Numérotation des lignes :

pourea
  1. C POUREA SOURCE CB215821 19/03/18 21:15:48 10160
  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.  
  22. C Initialisation à ZERO sinon SIGSEV parfois car utilisation de termes non initialises
  23. CALL ZERO(PP,3,3)
  24. C
  25. C CALCUL DE VNP
  26. C
  27. VNP(1)=XX(2)+XDDL(7)-XX(1)-XDDL(1)
  28. VNP(2)=YY(2)+XDDL(8)-YY(1)-XDDL(2)
  29. VNP(3)=ZZ(2)+XDDL(9)-ZZ(1)-XDDL(3)
  30. RLVNP=SQRT(VNP(1)*VNP(1)+VNP(2)*VNP(2)+VNP(3)*VNP(3))
  31. IF(RLVNP.EQ.0.D0) THEN
  32. KERRE=1
  33. RETURN
  34. ENDIF
  35. DO 1 I=1,3
  36. VNP(I)=VNP(I)/RLVNP
  37. 1 OMEG(I)=VNP(I)
  38. C
  39. C CALCUL DE LA MATRICE DE PASSAGE
  40. C
  41. CALL POUPAS(XX,YY,ZZ,OMEG,PP,KERRE)
  42. IOMEG=KERRE
  43. KERRE=0
  44. C
  45. C PASSAGE DES DEPLACEMENTS EN AXES LOCAUX - ROTATION PSI
  46. C
  47. CALL POUVEC(XDDL,WORK,P,1)
  48. PSI=(WORK(4)+WORK(10))*0.5D0
  49. C2=PP(1,1)*VNP(1)+PP(1,2)*VNP(2)+PP(1,3)*VNP(3)
  50. S2=PP(2,1)*VNP(1)+PP(2,2)*VNP(2)+PP(2,3)*VNP(3)
  51. C3=COS(PSI)
  52. S3=SIN(PSI)
  53. C
  54. C REACTUALISATION DU VECTEUR
  55. C
  56. IF(IOMEG.EQ.0) THEN
  57. DO 2 I=1,3
  58. VN(I)=0.D0
  59. DO 2 J=1,3
  60. 2 VN(I)=VN(I)+PP(I,J)*VECT(J)
  61. OMEG(1)=C2*VN(1)-S2*C3*VN(2)+S2*S3*VN(3)
  62. OMEG(2)=S2*VN(1)+C2*C3*VN(2)-C2*S3*VN(3)
  63. OMEG(3)=S3*VN(2)+C3*VN(3)
  64. DO 3 I=1,3
  65. VECT(I)=0.D0
  66. DO 3 J=1,3
  67. 3 VECT(I)=VECT(I)+PP(J,I)*OMEG(J)
  68. ELSE
  69. VN(1)=0.D0
  70. VN(2)=C3
  71. VN(3)=S3
  72. DO 6 I=1,3
  73. VECT(I)=0.D0
  74. DO 6 J=1,3
  75. 6 VECT(I)=VECT(I)+P(J,I)*VN(J)
  76. ENDIF
  77. END
  78.  
  79.  

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