Télécharger beton.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : beton.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Beton.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. * CAS TEST DU 93/01/19 PROVENANCE : NAH
  8. OPTION ECHO 1;
  9. GRAPH='N';
  10. SAUT PAGE ;
  11. *
  12. *-----------------------------------------------------*
  13. * *
  14. * Exemple simple de calcul mécanique élastique *
  15. * plastique beton fonctionnant avec les nouvelles *
  16. * structures MMODEL et MCHAMELEM *
  17. * *
  18. * Un cube en béton soumis à un chargement de *
  19. * traction uniaxial suivant l'axe OZ *
  20. * *
  21. * La limite en traction simple du béton LTR1 est *
  22. * égale à 4 MPa *
  23. * *
  24. *-----------------------------------------------------*
  25. * Définition des options
  26. *
  27. OPTI DIME 3 ELEM CUB8 MODEL TRIDIM ;
  28. *opti impi 9 ;
  29. N = 1 ;
  30. *
  31. * Maillage d'un cube de N*N*N éléments CUB8
  32. *
  33. DENS N ;
  34. *
  35. A1 = 0 0 0 ;
  36. A2 = N 0 0 ;
  37. A3 = N N 0 ;
  38. A4 = 0 N 0 ;
  39. B1 = 0 0 N ;
  40. B2 = N 0 N ;
  41. B3 = N N N ;
  42. B4 = 0 N N ;
  43. *
  44. A1A2 = A1 D N A2 ;
  45. A2A3 = A2 D N A3 ;
  46. A3A4 = A3 D N A4 ;
  47. A4A1 = A4 D N A1 ;
  48. B1B2 = B1 D N B2 ;
  49. B2B3 = B2 D N B3 ;
  50. B3B4 = B3 D N B4 ;
  51. B4B1 = B4 D N B1 ;
  52. A1B1 = A1 D N B1 ;
  53. A2B2 = A2 D N B2 ;
  54. A3B3 = A3 D N B3 ;
  55. A4B4 = A4 D N B4 ;
  56. *
  57. A1A2A3A4 = A1A2 A2A3 A3A4 A4A1 DALLE PLAN ;
  58. B1B2B3B4 = B1B2 B2B3 B3B4 B4B1 DALLE PLAN ;
  59. *
  60. A1A2B2B1 = A1A2 A2B2 (INVE B1B2) (INVE A1B1)
  61. DALLE PLAN ;
  62. A4A3B3B4 = (INVE A3A4) A3B3 B3B4 (INVE A4B4)
  63. DALLE PLAN ;
  64. *
  65. A2A3B3B2 = A2A3 A3B3 B2B3 (INVE A2B2)
  66. DALLE PLAN ;
  67. A1A4B4B1 = (INVE A4A1) A4B4 B4B1 (INVE A1B1)
  68. DALLE PLAN ;
  69. *
  70. TOTAL = PAVE A1A2A3A4 B1B2B3B4 A1A2B2B1 A4A3B3B4
  71. A2A3B3B2 A1A4B4B1 ;
  72. *
  73. *TRAC TOTAL CACH (1000 -2000 1000) ;
  74. *
  75. *
  76. * Formulation, matériau défini à l'aide des
  77. * opérateurs MODE et MATR
  78. *
  79. MO = MODE TOTAL MECANIQUE ELASTIQUE PLASTIQUE
  80. BETON CUB8 ;
  81. MA = MATE MO YOUN 40000 NU 0.20 ALPH 1E-5
  82. LTR1 4 ;
  83. *
  84. * Conditions de blocages
  85. *
  86. ENC1 = BLOQ UX A1A4B4B1 ;
  87. ENC2 = BLOQ UY A1A2B2B1 ;
  88. ENC3 = BLOQ UZ A1A2A3A4 ;
  89. ENC = ENC1 ET ENC2 ET ENC3 ;
  90. *
  91. * Blocage pour déplacements imposés
  92. *
  93. * END1 = BLOQ UX A2A3B3B2 ;
  94. * END2 = BLOQ UY A4A3B3B4 ;
  95. END3 = BLOQ UZ B1B2B3B4 ;
  96. * END = END1 ET END2 ET END3 ;
  97. END = END3 ;
  98. EN = ENC ET END ;
  99. *
  100. *
  101. * Valeur des déplacements imposés
  102. *
  103. * FEXT1 = DEPI END1 1E-4 ;
  104. * FEXT2 = DEPI END2 1E-4 ;
  105. FEXT3 = DEPI END3 1E-4 ;
  106. * FEXT = FEXT1 ET FEXT2 ET FEXT3 ;
  107. FEXT = FEXT3 ;
  108. *
  109. * Force en fonction du temps
  110. *
  111. lt1 = prog 0. 10. ;
  112. lt2 = prog 0. 10. ;
  113. ev = evol rouge manu temps lt1 force lt2 ;
  114. chat = char DIMP fext ev ;
  115.  
  116. *
  117. * Préparation à nonlin
  118. *
  119. ltt = prog 0 1 2 ;
  120. tab1 = table;
  121. tab1.blocages_mecaniques = EN;
  122. tab1.caracteristiques = ma;
  123. tab1.modele = mo;
  124. tab1.chargement = chat;
  125. tab1.temps_calcules = ltt;
  126. tab1.hypothese_deformations = lineaire;
  127. pasapas tab1;
  128. *
  129. * Traitement des résultats
  130. *
  131. ptt = prog 0 ;
  132. psz = prog 0 ;
  133. ic = 0 ;
  134. repeter mabou 2 ;
  135. ic = ic + 1 ;
  136. dd = tab1 . deplacements . ic;
  137. ptt = ptt et ( prog ( extr dd uz b3 ) ) ;
  138. ss = tab1.contraintes.ic ;
  139. sz = extr ( chan chpo mo ss ) smzz b3 ;
  140. psz = psz et ( prog sz ) ;
  141. fin mabou ;
  142. *
  143. evsz = evol vert manu deplacement ptt contrainte psz ;
  144. *dess evsz ;
  145. *
  146. * Code de fonctionnement
  147. *
  148. si ( (abs ( sz - 2. )) < 1.E-6 ) ;
  149. erre 0 ;
  150. sino ;
  151. list sz;
  152. erre 5 ;
  153. finsi ;
  154. *
  155. fin ;
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  

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