Télécharger la plaquette Cast3M


fr en

Télécharger rotor2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  rotor2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ************************************************************************      
  6. *                                                                      *
  7. *  Mots-clés : Vibrations, calcul modal, machines tournantes,          *
  8. *              poutre, modes complexes, reponse frequentielle          *
  9. *                                                                      *
  10. *                                                                      *
  11. *  Test de GYROSCOPIQUE, CAMPBELL et BALOURD pour les elements de poutre*
  12. *  Etude d'une machine tournante dans le repere inertiel (ou fixe)     *
  13. *  Elements: TIMO                                                      *
  14. *                                                                      *
  15. *  Volant d'inertie sur 2 appuis simples                               *
  16. *  et 1 appui élastique isotropes ou anisotropes                       *
  17. *  Exemple Lalanne P.13 +  cas test 14 de ROTOR                        *
  18. *                                                                      *
  19. *  Auteur: Didier Combescure Avril 2003                                *
  20. *                                                                      *                                                                         
  21. ************************************************************************      
  22. *GRAPH = 'Y';
  23. GRAPH = 'N';
  24. *
  25. *
  26. * Caractéristiques de l'arbre et du volant d'inertie
  27. *
  28. Ltot = 0.40 ;  
  29. L1 = Ltot/3.;
  30. L2 = 2.*Ltot/3.;
  31. L1 = 0.133;
  32. L2 = 0.266;
  33. *                                                        
  34. Re1 = 0.01 ;
  35. Ri2 = 0.01 ;
  36. Ri2 = 0.0 ;
  37. Re2 = 0.15 ;
  38. hvol = 0.03;
  39. *                                                             
  40. E1 = 2.e+11 ;                                                           
  41. Nu1 = 0.3 ;                                                             
  42. Rho1 = 7800. ;
  43. *
  44. Sarbre = pi*(Re1**2);
  45. Inxarb = pi*(Re1**4)/2.;
  46. Inyarb = pi*(Re1**4)/4.;
  47. *
  48. Svolant = pi*((Re2**2) - (Ri2**2));
  49. Inxvol = pi*((Re2**4) - (Ri2**4))/2.;
  50. Inyvol = pi*((Re2**4) - (Ri2**4))/4.;
  51. *
  52. * Valeur Lalanne
  53. *
  54. *Mvol = Svolant*hvol*rho1;
  55. *Inyvol = (Mvol/12.)*((3*(Re2**2)) - (3*(Ri2**2)) + (hvol**2));
  56. *Inyvol = 0.09427/(hvol*rho1);
  57. *
  58. *
  59. * Raideur et amortissement des supports
  60. *
  61. * Cas 1
  62. *
  63. *Ky = 0.;
  64. *Kz = Ky;
  65. *Cy = 0.;
  66. *Cz = Cy;
  67. *
  68. * Cas 2
  69. *
  70. Ky = 500000.;
  71. Kz = 0.;
  72. Cy = 0.;
  73. Cz = 0.;
  74. *
  75. * Cas 3
  76. *
  77. *bet= 0.0002;
  78. bet= 0.015;
  79. *bet= 0.026;
  80. *Ky = 200000.;
  81. *Kz = 500000.;
  82. *Cy = 200000.*bet;
  83. *Cz = 500000.*bet;
  84. *                                                                       
  85. OPTI 'DIME' 3 'ELEM' SEG2;                                                     
  86. *
  87. P0 = 0. 0. 0.;
  88. P1volm = (L1 - (0.5*hvol))  0.  0. ;
  89. P1 = (L1)  0.  0. ;
  90. P1volp = (L1 + (0.5*hvol))  0.  0. ;
  91. P1 = L1  0.  0. ;
  92. P2 = L2  0.  0. ;
  93. P3 = Ltot  0.  0. ;
  94. GeoArbre = (D 20 P0 P1volm) et (D 20 P1volp P2) et (D 20 P2 P3); 
  95. GeoVolan = (D 1 P1volm P1) et (D 1 P1 P1volp);
  96. *
  97. Geotot = GeoArbre et GeoVolan;
  98. *
  99. FacUnite = (2.*pi/60.);    
  100. *
  101. * Elément POUT
  102. *    
  103. *        
  104. Mod1a = MODELI GeoArbre 'MECANIQUE' TIMO;                                          
  105. Mat1a = MATE Mod1a 'YOUNG' E1 'NU' Nu1 'RHO' Rho1
  106.         'SECT' Sarbre 'INRY' Inyarb 'INRZ' Inyarb 'TORS' Inxarb
  107.         'OMEG' 1.;    
  108. Mod1b = MODELI GeoVolan 'MECANIQUE' TIMO;                                          
  109. Mat1b = MATE Mod1b 'YOUNG' E1 'NU' Nu1 'RHO' Rho1
  110.         'SECT' Svolant 'INRY' Inyvol 'INRZ' Inyvol 'TORS' Inxvol
  111.         'OMEG' 1.;    
  112. *
  113. mod1 = mod1a et mod1b;
  114. mat1 = mat1a et mat1b;
  115. *
  116. SI (EGA GRAPH 'Y');
  117. *
  118.  TABR = TABLE;
  119.  TABR. 'POUTRE' = TABLE;
  120.  (TABR. 'POUTRE'). 'CIRCULAIRE' = VRAI;
  121.  (TABR. 'POUTRE'). 'NCIRC' = 20;
  122.  TABR. 'TUYAU' = TABLE;
  123.  (TABR. 'TUYAU'). 'NCIRC' = 20;
  124.  (TABR. 'TUYAU'). 'NEPAI' = 1;
  125. *
  126.  MESHVOL = POUT2MAS mod1 mat1 GAUSS TABR;
  127.  trac meshvol cach;
  128.  MESHVOL = POUT2MAS mod1 mat1 MASSIF  TABR;;
  129.  elim 0.0001 MESHVOL;
  130.  trac meshvol cach;                                                                
  131. FINSI;
  132. *                                                                      
  133. * Cond. limites : 2 appuis 
  134. *
  135. Blo2 = (BLOQ DEPL P0);   
  136. Blo3 = (BLOQ DEPL P3);                             
  137. Blo1 = BLOQ 'RX' Geotot;  
  138. *
  139. Kxp= (APPUI UY KY p2);
  140. Kyp= (APPUI UZ KZ p2);
  141. Cxp= (APPUI UY CY p2);
  142. Cyp= (APPUI UZ CZ p2);
  143. *
  144. Kappui = Kxp et Kyp;     
  145. Cappui = Cxp et Cyp;                       
  146. *                                                                      
  147. MMass1 = MASS Mod1 (Mat1) ;                                   
  148. MRigi1 = RIGI Mod1 (Mat1) ;   
  149. *
  150. *
  151. * Matrice de couplage gyroscopique
  152. *
  153.  CG1 = GYROSCOPIQUE Mod1 Mat1;
  154. *                                   
  155. *                                                                       
  156. TBasR1 = VIBR 'PROCHE' (prog 100.) (lect 2)
  157.            (MRigi1 ET Kappui et Blo1 ET Blo2 et Blo3) 
  158.              MMass1  ;
  159. *
  160. TMod1 = TBasR1 . 'MODES' ;
  161. NbMod1 = (DIME TMod1) - 2 ;
  162. TInd1 = INDEX (TMod1 . 1) ;
  163. *
  164. SI (EGA GRAPH 'Y');
  165.  oeil = 1. 10. 0. ;
  166.  REPETER bbb1 NbMod1 ;
  167.   i = &bbb1 ;
  168.   Frq1 = TMod1 . i . 'FREQUENCE' ;
  169.   Def1 = TMod1 . i . 'DEFORMEE_MODALE' ;
  170.   MESS 'Mode ' i ':' Frq1 'Hz' ;
  171.   TITRE Frq1 ;
  172. *  trac oeil ((defo geotot def1) et (defo geotot def1 0.));
  173.  FIN bbb1 ;
  174. *
  175.  TABR = TABLE;
  176.  TABR. 'POUTRE' = TABLE;
  177.  (TABR. 'POUTRE'). 'CIRCULAIRE' = VRAI;
  178.  (TABR. 'POUTRE'). 'NCIRC' = 20;
  179.  TABR. 'TUYAU' = TABLE;
  180.  (TABR. 'TUYAU'). 'NCIRC' = 20;
  181.  (TABR. 'TUYAU'). 'NEPAI' = 1;
  182.  TABR.'DEPLACEMENTS' = TABLE;
  183.  (TABR.'DEPLACEMENTS'). 1 = (TMod1 . 1) . 'DEFORMEE_MODALE' ;
  184. *
  185.  MESHVOL = POUT2MAS mod1 mat1 GAUSS TABR;
  186.  trac (defo MESHVOL ((TABR.'DEPLACEMENTS'). 1)) cach;
  187. *
  188. FINSI;
  189. *
  190. *-Projection des matrices assemblées sur la base réelle
  191. *
  192.   MMasP1 = PJBA TBasR1 MMass1 ;
  193.   Cg1P = PJBA TBasR1 Cg1 ;
  194.   MRigP1 = PJBA TBasR1  (MRigi1) ;
  195.   KApp1p =  PJBA TBasR1 (Kappui et Blo1 ET Blo2 et Blo3);
  196.   CApp1p =  PJBA TBasR1 Cappui;
  197. *
  198. * Calcul du diagramme de Campbell - Modèle 1 (POUTRE)
  199. *
  200. TAB1 = TABLE;
  201. TAB1.'BASE_MODALE' = TBasR1 ;
  202. TAB1.'MASS_PROJ' = MMasP1;
  203. TAB1.'RIGI_PROJ' = MRigP1 et KApp1p ;
  204. TAB1.'GYRO_PROJ' = FacUnite*Cg1P;
  205. TAB1.'AMOR_PROJ' = CApp1p;
  206. TAB1.'AFFICHAGE' = VRAI;
  207. PRFREQ = (prog 0. pas 1. 120.)*100.;
  208. CAMPBELL TAB1 PRFREQ;
  209. *
  210. TMod1 = (TAB1.'BASE_MODALE') . 'MODES' ;
  211. NbModR = (DIME ((TAB1.'BASE_MODALE').'MODES')) - 2 ;
  212. NbModC = 2*NbModR;
  213. repeter lab1 NbModC;
  214.   i1= &lab1;
  215.   si (i1 ega 1);
  216.      EVTOT = (TAB1. i1). 'FREQUENCE_REELLE' ;
  217.       EVAMOT = (TAB1. i1). 'AMORTISSEMENT';
  218.    sinon;
  219.       EVTOT = EVTOT et 
  220.                ((TAB1. i1). 'FREQUENCE_REELLE' );
  221.       EVAMOT = EVAMOT et
  222.                ((TAB1. i1). 'AMORTISSEMENT');
  223.    finsi;
  224. FIN lab1;
  225. *
  226. SI (EGA GRAPH 'Y');
  227.  Titre 'Diagramme de Campbell';
  228.  dess (evtot et
  229.    (evol rouge manu (prog 0. 12000.) ((prog 0. 12000.)/60.))
  230.  et (evol rouge manu (prog 0. 12000.) ((prog 0. 12000.)/120.)) )
  231.      xbord 0. 12000.  ybord 0. 80. ;
  232.  dess (EVAMOT)  xbord 0. 12000.;
  233. FINSI;
  234. *-----------------------------
  235. *     Réponse à un balourd   *
  236. *-----------------------------
  237. md = 0.0001*(0.15);
  238. FbalReel = FORCE ( 0. (md) 0. ) P1;
  239. FbalImag = FORCE ( 0. 0. ((-1.)*md) ) P1;
  240. prfreq = FacUnite*(prog 100. pas 100. 2500. pas 2. 3500 
  241.     pas 100. 12000.);
  242. *
  243. TAB1 = TABLE;
  244. TAB1.'BASE_MODALE' = TBasR1 ;
  245. TAB1.'MASS_PROJ' = MMasP1;
  246. TAB1.'RIGI_PROJ' = MRigP1 et KApp1p;
  247. TAB1.'GYRO_PROJ' = Cg1P;
  248. TAB1.'AMOR_PROJ' = CApp1p;
  249. TAB1.'AFFICHAGE' = VRAI;
  250. TAB1.'FBALOURD' = FbalReel;
  251. *TAB1.'FBALREEL' = FbalReel;
  252. *TAB1.'FBALIMAG' = FbalImag;
  253. TAB1.'SAUVDEFO' = FAUX;
  254. TAB1.'REPONSE' = TABLE;
  255. (TAB1.'REPONSE'). 1  = TABLE;
  256. ((TAB1.'REPONSE'). 1 ). 'POINT' = P1;
  257. BALOURD TAB1 PRFREQ;
  258. *
  259. EvUYbal = ((TAB1.'REPONSE'). 1) . 'UY' ;
  260. EvUZbal = ((TAB1.'REPONSE'). 1) . 'UZ' ;
  261. EvGabal = ((TAB1.'REPONSE'). 1) . 'GAXE';
  262. EvPabal = ((TAB1.'REPONSE'). 1) . 'PAXE' ;
  263. *
  264. SI (EGA GRAPH 'Y');
  265.  DESS EvUYbal LOGY ybord 1.E-7 1E-3; 
  266.  DESS EvUZbal LOGY ybord 1.E-7 1E-3; 
  267.  dess (EvPabal et EvUYbal et EvUZbal et EvGabal) 
  268.                         LOGY ybord 1.E-7 1.E-3; 
  269.  dess  EvGabal LOGY ybord 1.E-7 1E-3; 
  270. FINSI;
  271. *
  272. AmpBal=extr (ipol (prog 377.) (extr EvGabal absc) 
  273.               (extr EvGabal ordo)) 1;
  274. AmpBaTh = 7.07D-5;
  275. SI (ABS((AmpBal - AmpBaTh)/AmpBaTh) < 5.D-2);
  276.   ERRE 0;
  277. SINON;
  278.   ERRE 5;
  279. FINSI;
  280. *
  281. FIN;
  282. * 
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  

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