Télécharger la plaquette Cast3M


fr en

Télécharger deveq7.eso

Retour à la liste

Numérotation des lignes : 

deveq7

  1. C DEVEQ7    SOURCE    BP208322  18/12/20    21:15:28     10048          
  2. C
  3. c       SUBROUTINE DEVEQ7(Q1,Q2,NA1,XK,XASM,XM,PDT,NPAS,FTOTA,NB1)
  4.       SUBROUTINE DEVEQ7(Q1,Q2,NA1,XM,PDT,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   PDT     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,*),FTOTA(NA1,*)
  35.       REAL*8 XOPM1(NB1,NB1,*),FAMOR(NA1,*)
  36. *
  37.       AUX4 = PDT * PDT / 6.D0
  38. *
  39. *                       .       2
  40. *     Q      = Q    + h Q    + h /6 M^-1 ( F    + 2 F
  41. *      i,m+1    i,m      i,m                i,m      i,m+1/2
  42. *                      .        .
  43. *              - C   ( Q    + 2 Q        ) )
  44. *                 ij    j,m      j,m+1/2
  45. *
  46. *
  47. C     CALCUL DU DEPLACEMENT
  48.  
  49. C     MATRICE M PLEINE
  50.       IF(NB1M.GT.1) THEN
  51.  
  52.         DO 21 IA = 1,NA1
  53.           AUX2 = 0.D0
  54.           DO 22 IB = 1,NB1
  55.             AUX2 = AUX2 + XOPM1(IA,IB,3) *
  56.      &      ((FTOTA(IB,3)-FAMOR(IB,3))+2.D0*(FTOTA(IB,2)-FAMOR(IB,2)))
  57.  22       CONTINUE
  58.           Q1(IA,1) = Q1(IA,3) + (PDT * Q2(IA,3)) + (AUX4 * AUX2)
  59.  21     CONTINUE
  60.  
  61. C     MATRICE M DIAGO
  62.       ELSE
  63.  
  64.         DO 20 I = 1,NA1
  65.            AUX2 = ( FTOTA(I,3)-FAMOR(I,3)
  66.      &      + 2.D0*(FTOTA(I,2)-FAMOR(I,2)) ) / XM(I,1)
  67.            Q1(I,1) = Q1(I,3) + (PDT * Q2(I,3)) + ( AUX4 * AUX2)
  68.  20     CONTINUE
  69.  
  70.       ENDIF
  71. *
  72.       END
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  

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