Télécharger devlc0.eso

Retour à la liste

Numérotation des lignes :

  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