Télécharger dyna11.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna11.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna11.dgibi: Jeux de données *
  5. * ---------------------------------- *
  6. * *
  7. GRAPH='N';
  8. SAUT PAGE ;
  9. *******************************************************
  10. * REPONSE TRANSITOIRE A UNE PRESSION INTEREN
  11. *******************************************************
  12. *
  13. * DESCRIPTION DU PROBLEME
  14. * ------------------------
  15. * UN CYLINDRE A PAROIS EPAISSES EST SOUMISE BRUTALEMENT
  16. * A UNE PRESSION INTERNE . IL S 'AGIT DE TROUVER LES
  17. * PICS DES DEPLACEMENTS DANS CERTAINS POINTS DE LA
  18. * STRUCTURE EN QUESTION .
  19. *
  20. * REFERENCE :
  21. * ----------
  22. * S. H. CRANDALL,N.C. DAHL,AND T. J. LARDNER ,AN
  23. * INTRODUCTION TO THE MECHANICS OF SOLIDS,2ND ED.
  24. * (NEW YORK :MCGRAW-HILL BOOK COMPANY,1972) PP. 293-297
  25. *******************************************************
  26. OPTION ECHO 0 ;
  27. OPTIO DIME 2 MODE AXIS ELEM QUA4 ;
  28. *
  29. * MAILLAGE
  30. *
  31. P1 = 6. 0. ; P2 = 12. 0. ; P3 = 12. 8. ; P4 = 6. 8. ;
  32. L1 = D 6 P1 P2 ; L2 = D 8 P2 P3 ;L3 = D 6 P3 P4 ;
  33. L4 = D 8 P4 P1 ; LTOT = L1 ET L2 ET L3 ET L4 ;
  34. S1 = SURF LTOT PLAN ;
  35. *
  36. * MODELE
  37. *
  38. MOD1 = MODELE S1 MECANIQUE ELASTIQUE ;
  39. *
  40. * PROPRIETES DE MATERIAUX
  41. *
  42. RO = 0.28 * 0.002588 ;
  43. MAT1 = MATER MOD1 YOUN 30.E6 NU 0.0 RHO RO ;
  44. *
  45. * CONDITIONS AUX LIMITES
  46. *
  47. CL1 = BLOQ L1 UZ ;CL2 = BLOQ L3 UZ ;
  48. *
  49. * MATRICE DE RIGIDITE
  50. *
  51. RIG1 = RIGIFITE MOD1 MAT1 ;
  52. *
  53. * MATRICE DE MASSE
  54. *
  55. MAS1 = MASSE MOD1 MAT1 ;
  56. *
  57. * CHARGEMENT
  58. *
  59. LI1 = (PROG 0. 1.E-6 2.E-6 3.E-6 4.E-6 ) ET
  60. (PROG 5e-6 pas 1e-6 204E-6 );
  61. LI2 = (PROG 0. 0. 0. 0. 0. ) ET (PROG 200 * 1.) ;
  62. F1 = PRES MOD1 'MASS' 10. L4 ;
  63. EV = EVOL MANU T LI1 F(T) LI2 ;
  64. CHA1 = CHAR 'PRES' F1 EV ;
  65. *
  66. * PREPARATION DE LA TABLE DE DONNEES POUR DYNAMIC
  67. *
  68. LIS1 = PROG 0. PAS 1.E-6 200.E-6 ;
  69. TAB2 = TABLE ;
  70. TAB2.'DEPL' = MANU CHPO S1 2 UZ 0. UR 0. ;
  71. TAB2.'VITE' = MANU CHPO S1 2 UZ 0. UR 0. ;
  72. TAB2.'CHAR' = CHA1 ;
  73. TAB2.'RIGI' = RIG1 ;
  74. TAB2.'MASS' = MAS1 ;
  75. TAB2.'FREQ' = 250000. ;
  76. TAB2.'INST' = LIS1 ;
  77. TAB1 = DYNAMIC TAB2 ;
  78. *
  79. * EXTRAIRE LES OBJETS D'EVOLUTION DE DEPLACEMENTS
  80. *
  81. LT = PROG ; LUR1 = PROG ;
  82. LUR2 = PROG ;
  83. LUR3 = PROG ;
  84. I = 0 ;
  85. M1 = NOEUD 25 ;M2 = NOEUD 62 ; M3 = NOEUD 36 ;
  86. NBB = (DIME LIS1) - 1 ;
  87. REPETER BB1 NBB ;
  88. I = I + 1 ;
  89. LT = LT ET( PROG(TAB1 . I . 'TEMP' ));
  90. DEP1 = TAB1 . I . 'DEPL' ;
  91. LUR1 = LUR1 ET( PROG (EXTR DEP1 UR M1) );
  92. LUR2 = LUR2 ET( PROG (EXTR DEP1 UR M2) );
  93. LUR3 = LUR3 ET( PROG (EXTR DEP1 UR M3) );
  94. FIN BB1 ;
  95. EV1 = EVOL MANU 'TIME' LT 'UR ' LUR1 ;
  96. EV2 = EVOL MANU 'TIME' LT 'UR ' LUR2 ;
  97. EV3 = EVOL MANU 'TIME' LT 'UR ' LUR3 ;
  98. *
  99. * --- TRACER LES COURBES
  100. *
  101. TABDD = TABLE ;
  102. TABDD.1 = 'MARQ ETOI NOLI ' ;
  103. TABDD.2 = 'MARQ CARR NOLI ' ;
  104. TABDD.3 = 'MARQ TRIA NOLI ' ;
  105. TABDD.'TITRE' = TABLE ;
  106. TABDD.'TITRE'. 1 = MOT 'N 25' ;
  107. TABDD.'TITRE'. 2 = MOT 'N 62' ;
  108. TABDD.'TITRE'. 3 = MOT 'N 36' ;
  109. 'TITR' 'RADIAL DISPLACEMENT THICK-WALLED CYLINDER ' ;
  110. SI (NEG GRAPH 'N');
  111. DESS (EV1 ET EV2 ET EV3) TABDD 'LEGE' ;
  112. FINSI ;
  113. *
  114. * COMPARAISON AVEC LA THEORIE
  115. *
  116. UR1 = (MAXI LUR1) / 3.326E-6 ; UR1TH = 2. ;
  117. ERR1 = ((UR1 - UR1TH) / (ABS UR1TH)) * 100. ;
  118. UR2 = (MAXI LUR2) / 2.773E-6 ; UR2TH = 2. ;
  119. ERR2 = ((UR2 - UR2TH) / (ABS UR2TH)) * 100. ;
  120. UR3 = (MAXI LUR3) / 2.673E-6 ; UR3TH = 2. ;
  121. ERR3 = ((UR3 - UR3TH) / (ABS UR3TH)) * 100. ;
  122. 'RAPPORT TRANSITOIRE / STATIQUE ' ;
  123. 'NODE 25 THEORIE 'UR1TH 'CASTEM ' UR1 '%ERREUR 'ERR1;
  124. 'NODE 62 THEORIE 'UR2TH 'CASTEM ' UR2 '%ERREUR 'ERR2;
  125. 'NODE 36 THEORIE 'UR3TH 'CASTEM ' UR3 '%ERREUR 'ERR2;
  126. *
  127. SI ((MAXI ABS (PROG ERR1 ERR2 ERR3 )) < 5);
  128. ERRE 0;
  129. SINON;
  130. ERRE 5;
  131. FINSI;
  132. FIN ;
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  

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