Télécharger dyna8.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna8.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna8.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test dyna8.dgibi: jeux de données *
  10. * --------------------------------- *
  11. *******************************************************
  12.  
  13. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  14. * SINON SI GRAPH DIFFERENT DE N TOUS
  15. * LES GRAPHIQUES SONT AFFICHES
  16.  
  17. GRAPH = 'N' ;
  18.  
  19. SAUT PAGE;
  20. SI (NEG GRAPH 'N') ;
  21. OPTI ECHO 1 ;
  22. OPTI TRAC PSC ;
  23. SINO ;
  24. OPTI ECHO 0 ;
  25. FINSI ;
  26.  
  27. SAUT PAGE;
  28. *************************************************************************
  29. * DYNA8
  30. *
  31. * CALCUL DES FREQUENCES PROPRES D'UNE POUTRE-CONSOLE
  32. * (COMPARAISON AVEC TIMOSHENKO)
  33. *
  34. * On recherche les frequences propres d'une poutre encastree a une
  35. * extremite (surface ST) et libre a l'autre extremite.
  36. *
  37. * Ce test est similaire au test DYNA6, mais au lieu d'utiliser des
  38. * elements POUT on utilise ici des elements massifs CU20.
  39. *
  40. * Les resultats sont compares aux resultats obtenus par methode analyti-
  41. * que et cites dans l'ouvrage de Timoshenko.
  42. *
  43. *************************************************************************
  44. TITR 'FREQUENCES PROPRES D UNE POUTRE-CONSOLE' ;
  45. OPTI DIME 3 ELEM CU20 ;
  46. TEMPS;
  47. *
  48. *----------------------- GEOMETRIE --------------------------------------
  49. *
  50. P0 = 0 0 0 ;
  51. P10 = 10 0 0 ;
  52. PY = 0 1 0 ;
  53. PZ = 0 0 2 ;
  54. LYZ = P0 DROI 1 PY ;
  55. ST = LYZ TRAN 1 PZ ;
  56. TRA = ST VOLU 10 TRAN P10 ;
  57. *
  58. OEIL1 = 15 15 15 ;
  59. *
  60. SI (NEG GRAPH 'N');
  61. TRAC 'QUAL' OEIL1 TRA ;
  62. FINSI;
  63. *
  64. *----------------------- MATERIAU ---------------------------------------
  65. *
  66. MM = MODE TRA MECANIQUE CU20 ;
  67.  
  68. MA = MATE MM YOUN 2.1E11 NU 0.3 RHO 7800 ;
  69. *
  70. *---------------- CONDITIONS AUX LIMITES, RIGIDITE ET MASSE -------------
  71. *
  72. INC = BLOQ DEPL ROTA ST ;
  73. K = RIGI MM MA ;
  74. KT = K ET INC ;
  75. MT = MASS MM MA ;
  76. *
  77. *---------------- CALCUL DES FREQUENCES PAR INTERVALLE ------------------
  78. *
  79. * Cette methode de recherche en trois intervalles est plus rapide que
  80. * la recherche en une seule fois sur l'intervalle 0-150.
  81. *
  82. AUTO0 = VIBR INTERVALLE 0. 50. BASSE 2 KT MT ;
  83. AUTO1 = VIBR INTERVALLE 50. 100. BASSE 2 KT MT ;
  84. AUTO2 = VIBR INTERVALLE 100. 150. BASSE 2 KT MT ;
  85. *
  86. *----------------- TRACE ET SORTIES -------------------------------------
  87.  
  88. *
  89. F1 = AUTO0 . 'MODES' . 1 . 'FREQUENCE';
  90. F2 = AUTO0 . 'MODES' . 2 . 'FREQUENCE';
  91. F3 = AUTO1 . 'MODES' . 1 . 'FREQUENCE';
  92. F4 = AUTO1 . 'MODES' . 2 . 'FREQUENCE';
  93. F5 = AUTO2 . 'MODES' . 1 . 'FREQUENCE';
  94. F6 = AUTO2 . 'MODES' . 2 . 'FREQUENCE';
  95. *
  96. SAUT PAGE ; SAUT 2 LIGNE ;
  97. *
  98. MESS ' MODE 1 FREQUENCE THEORIQUE 8.391 CALCULEE ' F1 ;
  99. MESS ' MODE 2 FREQUENCE THEORIQUE 16.78 CALCULEE ' F2 ;
  100. MESS ' MODE 3 FREQUENCE THEORIQUE 52.52 CALCULEE ' F3 ;
  101. MESS ' MODE 4 FREQUENCE THEORIQUE 59.75 CALCULEE ' F4 ;
  102. MESS ' MODE 5 FREQUENCE THEORIQUE 105.0 CALCULEE ' F5 ;
  103. MESS ' MODE 6 FREQUENCE THEORIQUE 129.6 CALCULEE ' F6 ;
  104. *
  105. *------------- CODE DE FONCTIONNEMENT ---------------------------------*
  106. ERR=TABLE;
  107. ERR.1=100*(ABS(8.391 -F1)/8.391);
  108. ERR.2=100*(ABS(16.78 -F2)/16.78);
  109. ERR.3=100*(ABS(52.52 -F3)/52.52);
  110. ERR.4=100*(ABS(59.75 -F4)/59.75);
  111. ERR.5=100*(ABS(105.0 -F5)/105.0);
  112. ERR.6=100*(ABS(129.6 -F6)/129.6);
  113. I=0;
  114. REPETER BOUC 6;
  115. I=I+1;
  116. SI (ERR.I < 1);
  117. ERRE 0;
  118. SINON;
  119. ERRE 5;
  120. FINSI;
  121. FIN;
  122.  
  123. TEMPS;
  124.  
  125. FIN;
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  

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