Télécharger dyne16.eso

Retour à la liste

Numérotation des lignes :

dyne16
  1. C DYNE16 SOURCE BP208322 20/09/18 21:16:22 10718
  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. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  8. * ________________________________________________ *
  9. * *
  10. * Decalage ou mise a zero dans les tableaux pour le pas suivant. *
  11. * *
  12. * Parametres: *
  13. * *
  14. * es Q1( , ) Vecteur des deplacements generalises *
  15. * es Q2( , ) Vecteur des vitesses generalisees *
  16. * es Q3( , ) Vecteur des accelerations generalisees *
  17. * e NA1 Nombre total d'inconnues sur la base A *
  18. * es FTOTA Forces exterieures totalisees, sur la base A *
  19. * es XPTB Tableau des deplacements des points de liaisons *
  20. * e NPLB Nombre maxi de points de liaisons *
  21. * e IDIMB Nombre maxi de ddl par liaisons *
  22. * *
  23. *--------------------------------------------------------------------*
  24. REAL*8 Q1(NA1,*),Q2(NA1,*),Q3(NA1,*),FTOTA(NA1,*)
  25. REAL*8 XPTB(NPLB,2,*),FINERT(NA1,*)
  26. *
  27. PARAMETER ( ZERO=0.D0 )
  28. *
  29. * Boucle sur les inconnues de la base A
  30. *
  31. DO 10 I=1,NA1
  32. Q1(I,3) = Q1(I,1)
  33. Q1(I,4) = Q1(I,2)
  34. Q1(I,1) = ZERO
  35. Q1(I,2) = ZERO
  36. *
  37. Q2(I,3) = Q2(I,1)
  38. Q2(I,4) = Q2(I,2)
  39. Q2(I,1) = ZERO
  40. Q2(I,2) = ZERO
  41. *
  42. Q3(I,3) = Q3(I,1)
  43. Q3(I,4) = Q3(I,2)
  44. Q3(I,1) = ZERO
  45. Q3(I,2) = ZERO
  46. *
  47. FTOTA(I,3) = FTOTA(I,1)
  48. FTOTA(I,4) = FTOTA(I,2)
  49. FTOTA(I,1) = ZERO
  50. FTOTA(I,2) = ZERO
  51. *
  52. FINERT(I,3) = FINERT(I,1)
  53. FINERT(I,4) = FINERT(I,2)
  54. FINERT(I,1) = ZERO
  55. FINERT(I,2) = ZERO
  56. 10 CONTINUE
  57. cbp2020-09 : ci dessous devenu inutile (cf devrco)
  58. cbp2020-09 DO 20 IP = 1,NPLB
  59. cbp2020-09 DO 22 ID = 1,IDIMB
  60. cbp2020-09 XPTB(IP,3,ID) = XPTB(IP,1,ID)
  61. cbp2020-09 XPTB(IP,4,ID) = XPTB(IP,2,ID)
  62. cbp2020-09 XPTB(IP,1,ID) = ZERO
  63. cbp2020-09 XPTB(IP,2,ID) = ZERO
  64. cbp2020-09 22 CONTINUE
  65. cbp2020-09 20 CONTINUE
  66. *
  67. END
  68.  
  69.  
  70.  
  71.  

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