Télécharger mtogy2.eso

Retour à la liste

Numérotation des lignes :

  1. C MTOGY2 SOURCE DC 05/02/15 21:15:50 5030
  2. SUBROUTINE MTOGY2(LRE,NDDL,REL,VROT,RELB)
  3. C=======================================================================
  4. C
  5. C Calcul de la matrice du couplage de Coriolis à partir de la matrice de masse
  6. C pour les éléments massifs en Fourier
  7. C
  8. C Didier COMBESCURE Janvier 2005
  9. C
  10. C=======================================================================
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8 (A-H,O-Z)
  13. DIMENSION REL(LRE,1),RELB(2*LRE,1),VROT(*),ROME(6,6)
  14. C
  15. NNOEU = LRE/NDDL
  16. C
  17. C
  18. DO 10 I=1,6
  19. DO 10 J=1,6
  20. ROME(I,J)=0.D0
  21. 10 CONTINUE
  22. C
  23. C Signe choisi pour etre coherent avec Genta....
  24. C
  25. ROME(1,5) =VROT(2)
  26. ROME(5,1) =VROT(2)
  27. ROME(2,6) = (-1.)*VROT(2)
  28. ROME(6,2) = (-1.)*VROT(2)
  29. C
  30. C
  31. DO 100 IN=1,NNOEU
  32. C
  33. DO 100 IM=1,NNOEU
  34. C
  35. DO 100 IDD1=1,NDDL
  36. C
  37. DO 100 IDD2=1,NDDL
  38. C
  39. RELB((6*IN)-7+2*IDD1,(6*IM)-7+2*IDD2)=
  40. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2-1)
  41. RELB((6*IN)-6+2*IDD1,(6*IM)-6+2*IDD2)=
  42. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2)
  43. RELB((6*IN)-6+2*IDD1,(6*IM)-7+2*IDD2)=
  44. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2-1)
  45. RELB((6*IN)-7+2*IDD1,(6*IM)-6+2*IDD2)=
  46. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2)
  47. C
  48. 100 CONTINUE
  49. C
  50. RETURN
  51. END
  52.  
  53.  
  54.  
  55.  

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