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    CB215821  26/08/24    21:17:23     12622          
  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 101 I=1,6
  19.        DO 10 J=1,6
  20.         ROME(I,J)=0.D0
  21.   10  CONTINUE
  22.   101 CONTINUE
  23. C
  24. C Signe choisi pour etre coherent avec Genta....
  25. C
  26.       ROME(1,5) =VROT(2)
  27.       ROME(5,1) =VROT(2)
  28.       ROME(2,6) = (-1.)*VROT(2)
  29.       ROME(6,2) = (-1.)*VROT(2)
  30. C
  31. C
  32.       DO   104 IN=1,NNOEU
  33. C
  34.         DO   103 IM=1,NNOEU
  35. C
  36.           DO   102 IDD1=1,NDDL
  37. C
  38.             DO   100 IDD2=1,NDDL
  39. C
  40.             RELB((6*IN)-7+2*IDD1,(6*IM)-7+2*IDD2)=
  41.      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2-1)
  42.             RELB((6*IN)-6+2*IDD1,(6*IM)-6+2*IDD2)=
  43.      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2)
  44.             RELB((6*IN)-6+2*IDD1,(6*IM)-7+2*IDD2)=
  45.      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2-1)
  46.             RELB((6*IN)-7+2*IDD1,(6*IM)-6+2*IDD2)=
  47.      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2)
  48. C
  49.   100  CONTINUE
  50.   102     CONTINUE
  51.   103   CONTINUE
  52.   104 CONTINUE
  53. C
  54.       RETURN
  55.       END
  56.  
  57.  
  58.  
  59.  
  60.  

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