Télécharger dyna9.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna9.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna9.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test dyna9.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. * DYNA9
  30. *
  31. * CALCUL DES FREQUENCES PROPRES D'UNE PLAQUE-CONSOLE
  32. *
  33. * Dans cet exemple on calcule les frequences propres d'une plaque
  34. * encastree sur un cote et libre sur les trois autres cotes.
  35. *
  36. * Les solutions calcules sont comparees a celles publiees par
  37. * M.V. BARTON dans l'article :
  38. *
  39. * VIBRATIONS OF RECTANGULAR AND SHEAR PLATES
  40. * JOU. OF APPL. MECH. , VOL.18 , 1951
  41. *
  42. ************************************************************************
  43. TITR 'FREQUENCES PROPRES D UNE PLAQUE CONSOLE' ;
  44.  
  45. OPTI DIME 3 ELEM QUA4 ;
  46. OPTI EPSI LINEAIRE;
  47.  
  48. TEMPS ;
  49. *
  50. *------------------------- GEOMETRIE -----------------------------------
  51. *
  52. P0 = 0 0 0 ;
  53. P1 = 0 1 0 ;
  54. P2 = 2 1 0 ;
  55. P3 = 2 0 0 ;
  56. L1 = D 8 P0 P1 ;
  57. L2 = D 16 P1 P2 ;
  58. L3 = D 8 P2 P3 ;
  59. L4 = D 16 P3 P0 ;
  60. S = L1 L2 L3 L4 DALL PLAN ;
  61. *
  62. OEIL1 = 10 10 10 ;
  63. *
  64. SI (NEG GRAPH 'N');
  65. TRAC 'QUAL' OEIL1 S ;
  66. FINSI;
  67. *
  68. *------------------------- MATERIAU ------------------------------------
  69. *
  70. MM = MODE S MECANIQUE COQ4 ;
  71.  
  72. MA = MATE MM YOUN 30E6 NU 0.3 RHO 7.32E-4 ;
  73. CA = CARA MM EPAI 0.01 ;
  74. MA =MA ET CA;
  75. *
  76. *-------------- CONDITIONS AUX LIMITES - CALCUL DES MATRICES -----------
  77. *
  78. INC = BLOQ DEPL ROTA L1 ;
  79. K = RIGI MM MA ;
  80. KT = K ET INC ;
  81. M = MASS MM MA;
  82. *
  83. *------------------------ CALCUL ---------------------------------------
  84. PROGREF = PROG 84.6 363.8 526.6 1187 ;
  85. PROG1 = PROG 85. 360. 530. 1200. ;
  86. AUTO = VIBR PROCHE PROG1 KT M ;
  87.  
  88. *
  89. *-------------------------- TRACE ET SORTIES ---------------------------
  90. PROGERR = VIDE 'LISTREEL';
  91. REPE SURMOD (DIME AUTO);
  92. INDICE = &SURMOD ;
  93. FREQEN = (TIRE AUTO FREQ NUME INDICE) ;
  94. FREQREF= EXTR PROGREF INDICE;
  95.  
  96. MESS ' MODE ' INDICE
  97. 'FREQUENCE REF :' (EXTR PROGREF INDICE)
  98. 'FREQUENCE CAL :' FREQEN ;
  99. MESS ' ERREUR RELATIVE : ' ERREUR ;
  100. DEPLA = (TIRE AUTO DEPL NUME INDICE) ;
  101. DEPLZ = EXCO DEPLA 'UZ';
  102. SI (NEG GRAPH 'N');
  103. TRAC OEIL1 DEPLZ (DEFO 0.5 DEPLA S);
  104. FINSI;
  105. SAUT 1 LIGNE ;
  106.  
  107. ERREUR = 100 * ((ABS(FREQREF -FREQEN) / FREQREF ));
  108. PROGERR = PROGERR ET (PROG ERREUR);
  109.  
  110. FIN SURMOD;
  111.  
  112. *---------------- CODE DE FONCTIONNEMENT -------------------------------
  113. ERRMAX = MAXI PROGERR;
  114.  
  115. SI (ERRMAX < 1);
  116. ERRE 0;
  117. SINON;
  118. ERRE 5;
  119. FINSI;
  120.  
  121. TEMPS ;
  122. FIN;
  123.  
  124.  
  125.  
  126.  
  127.  

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