Télécharger la plaquette Cast3M


fr en

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. tab1 . HYPOTHESE_DEFORMATIONS = 'QUADRATIQUE' ;
  64.  
  65. PASAPAS tab1 ;
  66.  
  67. CONTR = tab1 . CONTRAINTES . 13 ;
  68. SIXX1 = EXTR CONTR SMXX 1 1 1 ;
  69. CONTR = tab1 . CONTRAINTES . 20 ;
  70. SIXX2 = EXTR CONTR SMXX 1 1 1 ;
  71.  
  72. S1 = 3.43818E6 ;
  73. S2 = 1.04167E6 ;
  74.  
  75. err1 = abs ((SIXX1 - S1)/S1 * 100) ; 
  76. err2 = abs ((SIXX2 - S2)/S2 * 100) ; 
  77.  
  78. si ((err1 > .1) ou (err2 > .1)); 
  79.  erre 5; 
  80. finsi ; 
  81.  
  82. fin ;
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  

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