Télécharger la plaquette Cast3M


fr en

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. ;
  123. 'RAPPORT  TRANSITOIRE / STATIQUE       '   ;
  125. 'NODE 25   THEORIE 'UR1TH 'CASTEM  ' UR1 '%ERREUR 'ERR1;
  127. 'NODE 62   THEORIE 'UR2TH 'CASTEM  ' UR2 '%ERREUR 'ERR2;
  129. 'NODE 36   THEORIE 'UR3TH 'CASTEM  ' UR3 '%ERREUR 'ERR2;    
  130. *
  131. SI ((MAXI ABS (PROG ERR1 ERR2 ERR3 )) < 5);
  132. ERRE 0;
  133. SINON;
  134. ERRE 5;
  135. FINSI;
  136. FIN ;
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  

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