Télécharger mtokce.eso

Retour à la liste

Numérotation des lignes :

  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