Télécharger la plaquette Cast3M


fr en

Télécharger deveq1.eso

Retour à la liste

Numérotation des lignes : 

deveq1

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

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