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. OPTI EPSI LINEAIRE;
  47. TEMPS;
  48. *
  49. *----------------------- GEOMETRIE --------------------------------------
  50. *
  51. P0 = 0 0 0 ;
  52. P10 = 10 0 0 ;
  53. PY = 0 1 0 ;
  54. PZ = 0 0 2 ;
  55. LYZ = P0 DROI 1 PY ;
  56. ST = LYZ TRAN 1 PZ ;
  57. TRA = ST VOLU 10 TRAN P10 ;
  58. *
  59. OEIL1 = 15 15 15 ;
  60. *
  61. SI (NEG GRAPH 'N');
  62. TRAC 'QUAL' OEIL1 TRA ;
  63. FINSI;
  64. *
  65. *----------------------- MATERIAU ---------------------------------------
  66. *
  67. MM = MODE TRA MECANIQUE CU20 ;
  68.  
  69. MA = MATE MM YOUN 2.1E11 NU 0.3 RHO 7800 ;
  70. *
  71. *---------------- CONDITIONS AUX LIMITES, RIGIDITE ET MASSE -------------
  72. *
  73. INC = BLOQ DEPL ROTA ST ;
  74. K = RIGI MM MA ;
  75. KT = K ET INC ;
  76. MT = MASS MM MA ;
  77. *
  78. *---------------- CALCUL DES FREQUENCES PAR INTERVALLE ------------------
  79. *
  80. * Cette methode de recherche en trois intervalles est plus rapide que
  81. * la recherche en une seule fois sur l'intervalle 0-150.
  82. *
  83. AUTO0 = VIBR INTERVALLE 0. 50. BASSE 2 KT MT ;
  84. AUTO1 = VIBR INTERVALLE 50. 100. BASSE 2 KT MT ;
  85. AUTO2 = VIBR INTERVALLE 100. 150. BASSE 2 KT MT ;
  86. *
  87. *----------------- TRACE ET SORTIES -------------------------------------
  88.  
  89. *
  90. F1 = AUTO0 . 'MODES' . 1 . 'FREQUENCE';
  91. F2 = AUTO0 . 'MODES' . 2 . 'FREQUENCE';
  92. F3 = AUTO1 . 'MODES' . 1 . 'FREQUENCE';
  93. F4 = AUTO1 . 'MODES' . 2 . 'FREQUENCE';
  94. F5 = AUTO2 . 'MODES' . 1 . 'FREQUENCE';
  95. F6 = AUTO2 . 'MODES' . 2 . 'FREQUENCE';
  96. *
  97. SAUT PAGE ; SAUT 2 LIGNE ;
  98. *
  99. MESS ' MODE 1 FREQUENCE THEORIQUE 8.391 CALCULEE ' F1 ;
  100. MESS ' MODE 2 FREQUENCE THEORIQUE 16.78 CALCULEE ' F2 ;
  101. MESS ' MODE 3 FREQUENCE THEORIQUE 52.52 CALCULEE ' F3 ;
  102. MESS ' MODE 4 FREQUENCE THEORIQUE 59.75 CALCULEE ' F4 ;
  103. MESS ' MODE 5 FREQUENCE THEORIQUE 105.0 CALCULEE ' F5 ;
  104. MESS ' MODE 6 FREQUENCE THEORIQUE 129.6 CALCULEE ' F6 ;
  105. *
  106. *------------- CODE DE FONCTIONNEMENT ---------------------------------*
  107. ERR=TABLE;
  108. ERR.1=100*(ABS(8.391 -F1)/8.391);
  109. ERR.2=100*(ABS(16.78 -F2)/16.78);
  110. ERR.3=100*(ABS(52.52 -F3)/52.52);
  111. ERR.4=100*(ABS(59.75 -F4)/59.75);
  112. ERR.5=100*(ABS(105.0 -F5)/105.0);
  113. ERR.6=100*(ABS(129.6 -F6)/129.6);
  114. I=0;
  115. REPETER BOUC 6;
  116. I=I+1;
  117. SI (ERR.I < 1);
  118. ERRE 0;
  119. SINON;
  120. ERRE 5;
  121. FINSI;
  122. FIN;
  123.  
  124. TEMPS;
  125.  
  126. FIN;
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  

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