Télécharger mtogyf.eso

Retour à la liste

Numérotation des lignes :

  1. C MTOGYF SOURCE DC 05/02/15 21:15:51 5030
  2. SUBROUTINE MTOGYF(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 3D
  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. *
  23. ROME(1,6) = (-1.)*VROT(2)
  24. ROME(2,5) = (-1.)*VROT(2)
  25. ROME(6,1) = VROT(2)
  26. ROME(5,2) = VROT(2)
  27. C
  28. C
  29. DO 100 IN=1,NNOEU
  30. C
  31. DO 100 IM=1,NNOEU
  32. C
  33. DO 100 IDD1=1,NDDL
  34. C
  35. DO 100 IDD2=1,NDDL
  36. C
  37. RELB((6*IN)-7+2*IDD1,(6*IM)-7+2*IDD2)=
  38. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2-1)
  39. RELB((6*IN)-6+2*IDD1,(6*IM)-6+2*IDD2)=
  40. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2)
  41. RELB((6*IN)-6+2*IDD1,(6*IM)-7+2*IDD2)=
  42. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2-1)
  43. RELB((6*IN)-7+2*IDD1,(6*IM)-6+2*IDD2)=
  44. . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2)
  45. C
  46. 100 CONTINUE
  47. C
  48. RETURN
  49. END
  50.  
  51.  
  52.  
  53.  

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