Télécharger la plaquette Cast3M


fr en

Télécharger deveq2.eso

Retour à la liste

Numérotation des lignes : 

deveq2

  1. C DEVEQ2    SOURCE    BP208322  18/12/20    21:15:24     10048          
  2.       SUBROUTINE DEVEQ2(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 premier 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 XASM(NA1,*),XM(NA1,*),FTOTA(NA1,*),Q2(NA1,*)
  29.       REAL*8 FINERT(NA1,*)
  30. *
  31.       PDTS4 = PDT / 4.D0
  32. *
  33. *     .             4       .                                 .
  34. *     Q        = ________ ( Q    + h/4 ( G    + G        - A  Q   ) )
  35. *      i,m+1/2   4 + A  h    i,m          i,m    i,m+1/2    i  i,m
  36. *                     I
  37. *
  38. *     Boucle sur les inconnues de la base A
  39. *
  40.       DO 10 IA = 1,NA1
  41.          UNSM2 = 1.D0 / ( XM(IA,1) - FINERT(IA,2) )
  42.          UNSM3 = 1.D0 / ( XM(IA,1) - FINERT(IA,3) )
  43.          AUX1 = 4.D0 / ( 4.D0 + ( XASM(IA,1) * UNSM2 * PDT ) )
  44.          AUX2 = ( FTOTA(IA,3) * UNSM3 ) + ( FTOTA(IA,2) * UNSM2 )
  45.      &                        - ( XASM(IA,1) * UNSM3 * Q2(IA,3) )
  46.          AUX3 = Q2(IA,3) + ( PDTS4 * AUX2 )
  47.          Q2(IA,2) = AUX1 * AUX3
  48.  10      CONTINUE
  49. *     end do
  50. *
  51.       END
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  

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