Télécharger betdynlmt.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : betdynlmt.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *** essai de traction uniaxiale
  5.  
  6.  
  7. OPTI DIME 3 ELEM CUB8 ;
  8.  
  9. * Defintion du maillage *
  10. * cube à 8 noeuds *
  11.  
  12. P1 = 0. 0. 0. ;
  13. P2 = 0. 1. 0. ;
  14. L1 = DROI 1 P1 P2 ;
  15. COTEG = L1 TRAN (0. 0. 1) DINI 1. DFIN 1.;
  16. GEOM = VOLU COTEG TRAN (1. 0. 0.) DINI 1. DFIN 1.;
  17. COTED=FACE GEOM 2 ;
  18.  
  19. *** Modèle et matériau
  20.  
  21. MOD1= MODE GEOM MECANIQUE ELASTIQUE PLASTIQUE_ENDOM
  22. BETON_DYNAR_LMT ;
  23. MAT1 = MATE MOD1 'YOUN' 57.e9 'NU' 0.2 'RHO' 2500
  24. 'F0' 0.3 'Q1' 1.5 'Q2' 0.8 'Q3' 1.0
  25. 'SGM0' 70.e6 'XN' 15 'NVP' 1.5
  26. 'MVP' 1.e-4 'K' 15 'MDT' 0.5e-4
  27. 'NDT' 3. 'AT' 30000 'BT' 1.6
  28. 'MDC' 0.5e-3 'NDC' 20 'AC' 3000
  29. 'BC' 4. 'ED0' 9.5e-5;
  30.  
  31.  
  32. *** Conditions aux limites
  33.  
  34. CL1 = BLOQ P1 UZ UY;
  35. CL2 = BLOQ L1 UZ ;
  36. CL3 = BLOQ COTEG UX ;
  37. CLDEPL = BLOQ COTED UX ;
  38.  
  39. CLTOT = CL1 et CL2 et CL3 et CLDEPL ;
  40.  
  41. *** chargement
  42.  
  43. FO1 = DEPI CLDEPL 1. ;
  44. PROG1 = PROG 0. 20. ;
  45. PROG2 = PROG 0. 2.e-4 ;
  46. EVOL1 = EVOL MANU PROG1 PROG2 ;
  47. LIST1 = PROG 0. PAS 1. 20. ;
  48.  
  49. CHAR1 = CHARGEMENT FO1 EVOL1 DIMP ;
  50.  
  51. *** préparation du pasapas
  52.  
  53. tab1 = TABLE ;
  54. tab1 . PLASTIQUE = vrai ;
  55. tab1 . MODELE = mod1 ;
  56. tab1 . CARACTERISTIQUES=mat1 ;
  57. tab1 . CHARGEMENT = char1 ;
  58. tab1 . BLOCAGES_MECANIQUES = cltot ;
  59. tab1 . TEMPS_CALCULES = list1 ;
  60. tab1 . TEMPS_SAUVES = list1 ;
  61. tab1 . DYNAMIQUE = vrai ;
  62. tab1 . MOVA = 'D' ;
  63.  
  64. PASAPAS tab1 ;
  65.  
  66. CONTR = tab1 . CONTRAINTES . 13 ;
  67. SIXX1 = EXTR CONTR SMXX 1 1 1 ;
  68. CONTR = tab1 . CONTRAINTES . 20 ;
  69. SIXX2 = EXTR CONTR SMXX 1 1 1 ;
  70.  
  71. S1 = 3.43818E6 ;
  72. S2 = 1.04167E6 ;
  73.  
  74. err1 = abs ((SIXX1 - S1)/S1 * 100) ;
  75. err2 = abs ((SIXX2 - S2)/S2 * 100) ;
  76.  
  77. si ((err1 > .1) ou (err2 > .1));
  78. erre 5;
  79. finsi ;
  80.  
  81. fin ;
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  

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