Télécharger la plaquette Cast3M


fr en

Télécharger devlc0.eso

Retour à la liste

Numérotation des lignes : 

devlc0

  1. C DEVLC0    SOURCE    BP208322  15/07/22    21:15:22     8586
  2. C DEVLC0
  3.       SUBROUTINE DEVLC0(Q2,XASM,FAMOR,NA1,NB1C,IND)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. *                                                                    *
  8. *     Operateur DYNE                                                 *
  9. *     ________________________________________________               *
  10. *                                                                    *
  11. *     Calcul des forces d'amortissement sur la base A                *
  12. *     FAMOR = FAMOR - K*q                                            *
  13. *     Calcule le produit matrice vecteur                             *
  14. *     Distingue le cas matrice diagonal et pleine                    *
  15. *                                                                    *
  16. *     Parametres:                                                    *
  17. *                                                                    *
  18. * e   Q2(,)     Vitesse generalise                                   *
  19. * e   XASM      Matrice d'amortissement generalisee                  *
  20. * s   FAMOR     Forces  d'amortissement generalisee                  *
  21. * e   NA1,NB1C  Dimension de la matrice                              *
  22. *               (si NB1C=1, alors matrice diagonale)                 *
  23. * e   IND       Indice du pas                                        *
  24. *                                                                    *
  25. *    -Auteur, date de creation:                                      *
  26. *     Benoit PRABEL, 11/02/2015                                      *
  27. *                                                                    *
  28. *--------------------------------------------------------------------*
  29.       REAL*8 Q2(NA1,*),XASM(NA1,*),FAMOR(NA1,*)
  30. *
  31. *
  32.       IF (NB1C.EQ.1) THEN
  33.          DO 100 I=1,NA1
  34.            FAMOR(I,IND) = XASM(I,1)*Q2(I,IND)
  35.  100     CONTINUE
  36. c       ELSEIF(NB1C.NE.NA1) THEN
  37. c         CALL ERREUR(832)
  38. c on ne va pas tester a chaque fois que cette subroutine est appelee
  39.       ELSE
  40.          DO 200 I=1,NA1
  41.            FAMOR(I,IND) = 0.D0
  42.          DO 200 J=1,NB1C
  43.            FAMOR(I,IND) = FAMOR(I,IND) + XASM(I,J)*Q2(J,IND)
  44.  200     CONTINUE
  45.       ENDIF
  46.  
  47.       END
  48.  
  49.  
  50.  

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