Télécharger la plaquette Cast3M


fr en

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