Télécharger deveq4.eso

Retour à la liste

Numérotation des lignes :

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

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