Télécharger la plaquette Cast3M


fr en

Télécharger vibr10.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  vibr10.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. ************************************************************************
  6. *
  7. * Mots-clés : Vibrations, calcul modal, machines tournantes, 
  8. *             reponse frequentielle
  9. * 
  10. * Calcul de la reponse au balourd   
  11. *
  12. * Creation : J Charpentier ?
  13. *            Test de l'option 'BALOURD' de l'operateur VIBC
  14. *
  15. * Modif : B Prabel, 2015-09-09 : plutot que d'utiliser VIBC 'BALOU',
  16. *         on utilisera desormais IMPE, RESOU et RECO
  17. *                                                                       
  18. ************************************************************************
  19.  
  20. OPTI 'DIME' 3 'ELEM' SEG2 ;
  21. *
  22. *donnees
  23. P1 =  1. 0. 0. ;
  24. m1 = 40. ;
  25. k1 = 4.e5 ;
  26. f1 = ((k1 / m1) ** 0.5) / (2. * pi) ;
  27. MESS 'f1 =' f1 ;
  28. *
  29. TR1 = MANU 'POI1' P1 ;
  30. *
  31. *Matrices de masse et de rigidite :
  32. MMass = MANU 'RIGI' TR1 (MOTS 'UY' 'UZ') 
  33.         (PROG m1 0. 0. m1) ;
  34. MRigi = MANU 'RIGI' TR1 (MOTS 'UY' 'UZ') 
  35.         (PROG k1 0. 0. k1) ;
  36. *
  37. *deformees modales :
  38. Phi1 = MANU 'CHPO' TR1 2 'UY' 1. 'UZ' 0. ;
  39. Phi2 = MANU 'CHPO' TR1 2 'UY' 0. 'UZ' 1. ;
  40. *
  41. *points reperes des modes
  42. PR1 = 0. 0. 0. ;
  43. PR2 = 0. 0. 0. ;
  44. *
  45. *base des mods reels :
  46. TBasR = TABLE 'BASE_MODALE' ;
  47. TBasR.MODES = TABLE 'BASE_DE_MODES' ;
  48. TBasR.MODES . 'MAILLAGE' = TR1 ;
  49. TBasR.MODES . 1 = TABLE 'MODE' ;
  50. TBasR.MODES . 1 . 'NUMERO_MODE' = 1 ;
  51. TBasR.MODES . 1 . 'POINT_REPERE' = PR1 ;
  52. TBasR.MODES . 1 . 'FREQUENCE' = f1 ;
  53. TBasR.MODES . 1 . 'MASSE_GENERALISEE' = m1 ;
  54. TBasR.MODES . 1 . 'DEPLACEMENTS_GENERALISES' =
  55.   k1 - (2 * pi * f1) ;
  56. TBasR.MODES . 1 . 'DEFORMEE_MODALE' = Phi1 ;
  57. TBasR.MODES . 2 = TABLE 'MODE' ;
  58. TBasR.MODES . 2 . 'NUMERO_MODE' = 2 ;
  59. TBasR.MODES . 2 . 'POINT_REPERE' = PR2 ;
  60. TBasR.MODES . 2 . 'FREQUENCE' = f1 ;
  61. TBasR.MODES . 2 . 'MASSE_GENERALISEE' = m1 ;
  62. TBasR.MODES . 2 . 'DEPLACEMENTS_GENERALISES' =
  63.   k1 - (2 * pi * f1) ;
  64. TBasR.MODES . 2 . 'DEFORMEE_MODALE' = Phi2 ;
  65. *
  66. *vitesse (rad/s)
  67. omeg1 = 19.0985 ;
  68. *
  69. *-projection des matrices sur la base réelle :
  70.   MMasP = PJBA TBasR MMass ;
  71.   MRigP = PJBA TBasR MRigi ;
  72. mess 'MMasP='; list MMasP;
  73. mess 'MRigP='; list MRigP;
  74.  
  75. *-ecriture impedance
  76.   MMasP2 = IMPE MMasP 'MASSE';
  77.   MRigP2 = IMPE MRigP 'RAIDEUR';
  78. mess 'MMasP2='; list MMasP2;
  79. mess 'MRigP2='; list MRigP2;
  80.   MdynP2 = MRigP2 et ((omeg1 ** 2) * MMasP2);
  81. *
  82. *-calcul de la force d'excitation du balourd
  83.   MxD = 1e-4 ;
  84.   Fy1s = MxD * (omeg1 ** 2) ;
  85.   Fz1s = (-1) * MxD * (omeg1 ** 2) ;
  86.   Fy1c = MxD * (omeg1 ** 2) ;
  87.   Fz1c = MxD * (omeg1 ** 2) ;
  88. *
  89.   F1s = MANU 'CHPO' P1 2 'FY' Fy1s 'FZ' Fz1s ;
  90.   F1c = MANU 'CHPO' P1 2 'FY' Fy1c 'FZ' Fz1c ;
  91.   PjF1s = PJBA F1s TBasR ;
  92.   PjF1c = (PJBA F1c TBasR) 'EXCO' 'FALF' 'IFAL' ;
  93. *
  94.   PjFBal1 = PjF1s et PjF1c;
  95. mess 'PjFBal1='; list PjFBal1;
  96. *
  97. *-calcul de la reponse :
  98.   PjQBal1 = RESOU MdynP2 PjFBal1;
  99. mess 'PjQBal1='; list PjQBal1;
  100.  
  101. *-recombinaison modale
  102.   QRe1 = EXCO PjQBal1 'ALFA' 'ALFA';
  103.   QIm1 = EXCO PjQBal1 'IALF' 'ALFA';
  104.   RepRe1 = RECO QRe1 TBasR;
  105.   RepIm1 = RECO QIm1 TBasR;
  106. mess 'RepRe1='; list RepRe1;
  107. mess 'RepIm1='; list RepIm1;
  108. *
  109. *-calcul des pt et gd axes de l'ellipse :
  110.    YR1 = EXTR RepRe1 'UY' P1 ;
  111.    ZR1 = EXTR RepRe1 'UZ' P1 ;
  112.    MESS YR1 ZR1 ;
  113. * systeme isotrope --> ellipse = cercle
  114.  
  115. *-resultats attendus :
  116. An1= (MxD * (omeg1 ** 2)) / (k1 - ((omeg1 ** 2) * m1));
  117. MESS 'resultats attendus : ' An1 ;
  118. *
  119. Si (((Abs(ZR1+An1)/An1) < 0.001)
  120.  et ((Abs(YR1-An1)/An1) < 0.001));
  121.   Erre 0;
  122. Sinon;
  123.   Erre 5;
  124. Finsi;
  125.  
  126. FIN ; 
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  

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