Télécharger dyne16.eso

Retour à la liste

Numérotation des lignes :

  1. C DYNE16 SOURCE CHAT 05/01/12 23:14:35 5004
  2. SUBROUTINE DYNE16(Q1,Q2,Q3,NA1,FTOTA,XPTB,NPLB,IDIMB,FINERT)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. * *
  7. * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
  8. * ________________________________________________ *
  9. * *
  10. * D{calage ou mise @ z{ro dans les tableaux pour le pas suivant. *
  11. * *
  12. * Param}tres: *
  13. * *
  14. * es Q1( , ) Vecteur des d{placements g{n{ralis{s *
  15. * es Q2( , ) Vecteur des vitesses g{n{ralis{es *
  16. * es Q3( , ) Vecteur des acc{l{rations g{n{ralis{es *
  17. * e NA1 Nombre total d'inconnues sur la base A *
  18. * es FTOTA Forces ext{rieures totalis{es, sur la base A *
  19. * es XPTB Tableau des d{placements des points de liaisons *
  20. * e NPLB Nombre maxi de points de liaisons *
  21. * e IDIMB Nombre maxi de ddl par liaisons *
  22. * *
  23. * Auteur, date de cr{ation: *
  24. * *
  25. * Denis ROBERT-MOUGIN, le 16 juin 1989. *
  26. * *
  27. *--------------------------------------------------------------------*
  28. REAL*8 Q1(NA1,*),Q2(NA1,*),Q3(NA1,*),FTOTA(NA1,*)
  29. REAL*8 XPTB(NPLB,4,*),FINERT(NA1,*)
  30. *
  31. PARAMETER ( ZERO=0.D0 )
  32. *
  33. * Boucle sur les inconnues de la base A
  34. *
  35. DO 10 I=1,NA1
  36. Q1(I,3) = Q1(I,1)
  37. Q1(I,4) = Q1(I,2)
  38. Q1(I,1) = ZERO
  39. Q1(I,2) = ZERO
  40. *
  41. Q2(I,3) = Q2(I,1)
  42. Q2(I,4) = Q2(I,2)
  43. Q2(I,1) = ZERO
  44. Q2(I,2) = ZERO
  45. *
  46. Q3(I,3) = Q3(I,1)
  47. Q3(I,4) = Q3(I,2)
  48. Q3(I,1) = ZERO
  49. Q3(I,2) = ZERO
  50. *
  51. FTOTA(I,3) = FTOTA(I,1)
  52. FTOTA(I,4) = FTOTA(I,2)
  53. FTOTA(I,1) = ZERO
  54. FTOTA(I,2) = ZERO
  55. *
  56. FINERT(I,3) = FINERT(I,1)
  57. FINERT(I,4) = FINERT(I,2)
  58. FINERT(I,1) = ZERO
  59. FINERT(I,2) = ZERO
  60. 10 CONTINUE
  61. * end do
  62. DO 20 IP = 1,NPLB
  63. DO 22 ID = 1,IDIMB
  64. XPTB(IP,3,ID) = XPTB(IP,1,ID)
  65. XPTB(IP,4,ID) = XPTB(IP,2,ID)
  66. XPTB(IP,1,ID) = ZERO
  67. XPTB(IP,2,ID) = ZERO
  68. 22 CONTINUE
  69. * end do
  70. 20 CONTINUE
  71. * end do
  72. *
  73. END
  74.  
  75.  

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