Télécharger la plaquette Cast3M


fr en

Télécharger mtogy2.eso

Retour à la liste

Numérotation des lignes : 

mtogy2

  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