Télécharger la plaquette Cast3M


fr en

Télécharger mtokce.eso

Retour à la liste

Numérotation des lignes : 

mtokce

  1. C MTOKCE    SOURCE    BP208322  12/03/08    21:15:04     7294
  2. c       SUBROUTINE MTOKCE(LRE,NDDL,REL,VROT,RELB)
  3.       SUBROUTINE MTOKCE(LRE,NDDL,REL,ROME,RELB)
  4. C=======================================================================
  5. C
  6. C  Calcul de la matrice de raideur centrifuge à partir dela matrice de masse
  7. C  pour les éléments massifs en 3D
  8. C
  9. C     Didier COMBESCURE Mars 2003
  10. C     BP on sort le remplissage de ROME (cst) de la boucle (cf KCENT3.eso)
  11. C
  12. C=======================================================================
  13.       IMPLICIT INTEGER(I-N)
  14.       IMPLICIT REAL*8 (A-H,O-Z)
  15. c       DIMENSION REL(LRE,1),RELB(LRE,1),VROT(*),ROME(3,3)
  16.       DIMENSION REL(LRE,1),RELB(LRE,1),ROME(3,3)
  17. C
  18.       NNOEU = LRE/NDDL
  19. C
  20. c       ROME(1,1) = (-1.)*((VROT(2)**2) + (VROT(3)**2))
  21. c       ROME(2,2) = (-1.)*((VROT(1)**2) + (VROT(3)**2))
  22. c       ROME(3,3) = (-1.)*((VROT(1)**2) + (VROT(2)**2))
  23. c       ROME(1,2) = VROT(1)*VROT(2)
  24. c       ROME(1,3) = VROT(1)*VROT(3)
  25. c       ROME(2,3) = VROT(2)*VROT(3)
  26. c       ROME(2,1) = ROME(1,2)
  27. c       ROME(3,1) = ROME(1,3)
  28. c       ROME(3,2) = ROME(2,3)
  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((3*IN)-3+IDD1,(3*IM)-3+IDD2)=
  40.      . REL((3*IN)-3+IDD1,(3*IM)-3+IDD1)*ROME(IDD1,IDD2)
  41. C
  42.   100      CONTINUE
  43. C
  44.       RETURN
  45.       END
  46.  
  47.  
  48.  
  49.  
  50.  

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