Télécharger la plaquette Cast3M


fr en

Télécharger mtogyf.eso

Retour à la liste

Numérotation des lignes : 

mtogyf

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

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