Télécharger deveq1.eso

Retour à la liste

Numérotation des lignes :

  1. C DEVEQ1 SOURCE BP208322 15/07/22 21:15:15 8586
  2. SUBROUTINE DEVEQ1(Q1,Q2,NA1,XASM,XM,XDT,NPAS,FTOTA,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. * Calcul des deplacements generalises pour le premier demi-pas *
  11. * de temps. *
  12. * *
  13. * Param}tres: *
  14. * *
  15. * es Q1(.,.) Vecteur des deplacements generalises *
  16. * es Q2(.,.) Vecteur des vitesses generalisees *
  17. * e NA1 Nombre total d'inconnues dans la base A *
  18. * e XASM Vecteur des amortissements generalises *
  19. * e XM Vecteur des masses generalisees *
  20. * e XDT Valeurs des pas de temps *
  21. * e NPAS Numero du pas de temps considere *
  22. * e FTOTA Forces exterieures totalisees, sur la base A *
  23. * *
  24. * Auteur, date de creation: *
  25. * *
  26. * Denis ROBERT-MOUGIN, le 19 mai 1989. *
  27. * *
  28. *--------------------------------------------------------------------*
  29. REAL*8 Q1(NA1,*),Q2(NA1,*),FINERT(NA1,*)
  30. REAL*8 XASM(NA1,*),XM(NA1,*),XDT(*),FTOTA(NA1,*)
  31. *
  32. PDT = XDT(NPAS)
  33. PDTS2 = PDT * 0.5D0
  34. AUX4 = PDT * PDT / 24.D0
  35. *
  36. * . 2
  37. * Q = Q + h/2 Q + h /24 ( 4 G - G
  38. * i,m+1/2 i,m i,m i,m i,m-1/2
  39. * . .
  40. * - A ( 4 Q - Q ) )
  41. * i i,m i,m-1/2
  42. *
  43. *
  44. * Boucle sur les inconnues de la base A
  45. *
  46. DO 10 IA = 1,NA1
  47. UNSM3 = 1.D0 / ( XM(IA,1) - FINERT(IA,3) )
  48. UNSM4 = 1.D0 / ( XM(IA,1) - FINERT(IA,4) )
  49. AUX1 = XASM(IA,1) * ( ( 4.D0 * Q2(IA,3) * UNSM3 )
  50. & - ( Q2(IA,4) * UNSM4 ) )
  51. AUX2 = (4.D0*FTOTA(IA,3)*UNSM3) - (FTOTA(IA,4)*UNSM4)
  52. AUX3 = PDTS2 * Q2(IA,3)
  53. Q1(IA,2) = Q1(IA,3) + AUX3 + ( AUX4 * ( AUX2 - AUX1) )
  54. 10 CONTINUE
  55. * end do
  56. *
  57. END
  58.  
  59.  
  60.  
  61.  

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