Télécharger deveq7.eso

Retour à la liste

Numérotation des lignes :

  1. C DEVEQ7 SOURCE BP208322 15/07/22 21:15:19 8586
  2. C
  3. c SUBROUTINE DEVEQ7(Q1,Q2,NA1,XK,XASM,XM,XDT,NPAS,FTOTA,NB1)
  4. SUBROUTINE DEVEQ7(Q1,Q2,NA1,XM,XDT,NPAS,FTOTA,FAMOR,
  5. & XOPM1,NB1,NB1M)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. *--------------------------------------------------------------------*
  9. * *
  10. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  11. * ________________________________________________ *
  12. * *
  13. * Calcul des deplacements generalises pour le second demi-pas *
  14. * de temps. *
  15. * *
  16. * Param}tres: *
  17. * *
  18. * es Q1(.,.) Vecteur des deplacements generalises *
  19. * e Q2(.,.) Vecteur des vitesses generalisees *
  20. * e NA1 Nombre total d'inconnues dans la base A *
  21. * e XK Vecteur des raideurs generalisees *
  22. * e XASM Vecteur des amortissements generalises *
  23. * e XM Vecteur des masses generalisees *
  24. * e XDT Valeurs des pas de temps *
  25. * e NPAS Numero du pas de temps *
  26. * e FTOTA Forces exterieures totalisees, sur la base A *
  27. * *
  28. * Auteur, date de creation: *
  29. * *
  30. * Denis ROBERT-MOUGIN, le 19 mai 1989. *
  31. * *
  32. *--------------------------------------------------------------------*
  33. REAL*8 Q1(NA1,*),Q2(NA1,*)
  34. REAL*8 XM(NA1,*),XDT(*),FTOTA(NA1,*)
  35. REAL*8 XOPM1(NB1,NB1,*),FAMOR(NA1,*)
  36. *
  37. PDT = XDT(NPAS)
  38. AUX4 = PDT * PDT / 6.D0
  39. *
  40. * . 2
  41. * Q = Q + h Q + h /6 M^-1 ( F + 2 F
  42. * i,m+1 i,m i,m i,m i,m+1/2
  43. * . .
  44. * - C ( Q + 2 Q ) )
  45. * ij j,m j,m+1/2
  46. *
  47. *
  48. C CALCUL DU DEPLACEMENT
  49.  
  50. C MATRICE M PLEINE
  51. IF(NB1M.GT.1) THEN
  52.  
  53. DO 21 IA = 1,NA1
  54. AUX2 = 0.D0
  55. DO 22 IB = 1,NB1
  56. AUX2 = AUX2 + XOPM1(IA,IB,3) *
  57. & ((FTOTA(IB,3)-FAMOR(IB,3))+2.D0*(FTOTA(IB,2)-FAMOR(IB,2)))
  58. 22 CONTINUE
  59. Q1(IA,1) = Q1(IA,3) + (PDT * Q2(IA,3)) + (AUX4 * AUX2)
  60. 21 CONTINUE
  61.  
  62. C MATRICE M DIAGO
  63. ELSE
  64.  
  65. DO 20 I = 1,NA1
  66. AUX2 = ( FTOTA(I,3)-FAMOR(I,3)
  67. & + 2.D0*(FTOTA(I,2)-FAMOR(I,2)) ) / XM(I,1)
  68. Q1(I,1) = Q1(I,3) + (PDT * Q2(I,3)) + ( AUX4 * AUX2)
  69. 20 CONTINUE
  70.  
  71. ENDIF
  72. *
  73. END
  74.  
  75.  
  76.  
  77.  

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