Télécharger la plaquette Cast3M


fr en

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.  
  47. TEMPS ;
  48. *
  49. *------------------------- GEOMETRIE -----------------------------------
  50. *
  51. P0 = 0 0 0 ;
  52. P1 = 0 1 0 ;
  53. P2 = 2 1 0 ;
  54. P3 = 2 0 0 ;
  55. L1 = D 8 P0 P1 ;
  56. L2 = D 16 P1 P2 ;
  57. L3 = D 8 P2 P3 ;
  58. L4 = D 16 P3 P0 ;
  59. S = L1 L2 L3 L4 DALL PLAN ;
  60. *
  61. OEIL1 = 10 10 10 ;
  62. *
  63. SI (NEG GRAPH 'N');
  64.   TRAC 'QUAL' OEIL1 S ;
  65. FINSI;
  66. *
  67. *------------------------- MATERIAU ------------------------------------
  68. *
  69. MM = MODE S MECANIQUE COQ4 ;
  70.  
  71. MA = MATE MM YOUN 30E6 NU 0.3 RHO 7.32E-4 ;
  72. CA = CARA MM EPAI 0.01 ;
  73. MA =MA ET CA;
  74. *
  75. *-------------- CONDITIONS AUX LIMITES - CALCUL DES MATRICES -----------
  76. *
  77. INC = BLOQ DEPL ROTA L1 ;
  78. K   = RIGI MM MA ;
  79. KT  = K ET INC ;
  80. M   = MASS MM MA;
  81. *
  82. *------------------------ CALCUL ---------------------------------------
  83. PROGREF = PROG 84.6 363.8 526.6 1187  ;
  84. PROG1   = PROG 85.  360.  530.  1200. ;
  85. AUTO    = VIBR PROCHE PROG1 KT M      ;
  86.  
  87. *
  88. *-------------------------- TRACE ET SORTIES ---------------------------
  89. PROGERR = VIDE 'LISTREEL';
  90. REPE SURMOD (DIME PROG1);
  91.   INDICE = &SURMOD ;
  92.   FREQEN = AUTO . MODES . INDICE . FREQUENCE ;
  93.   FREQREF= EXTR PROGREF INDICE;
  94.  
  95.   MESS ' MODE ' INDICE 
  96.         'FREQUENCE REF :' (EXTR PROGREF INDICE)
  97.         'FREQUENCE CAL :' FREQEN ;
  98.   MESS '  ERREUR RELATIVE : ' ERREUR ;
  99.   DEPLA = AUTO . MODES . INDICE . DEFORMEE_MODALE ;
  100.   DEPLZ = EXCO DEPLA 'UZ';
  101.   SI (NEG GRAPH 'N');
  102.     TRAC OEIL1 DEPLZ (DEFO 0.5 DEPLA S);
  103.   FINSI;
  104.   SAUT 1 LIGNE ;
  105.  
  106.   ERREUR = 100 * ((ABS(FREQREF -FREQEN) / FREQREF ));
  107.   PROGERR = PROGERR ET (PROG ERREUR);
  108.  
  109. FIN SURMOD;
  110.  
  111. *---------------- CODE DE FONCTIONNEMENT -------------------------------
  112. ERRMAX = MAXI PROGERR;
  113.  
  114. SI (ERRMAX < 1);
  115.   ERRE 0;
  116. SINON;
  117.   ERRE 5;
  118. FINSI;
  119.  
  120. TEMPS ;
  121. FIN;
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  

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