* fichier : beton.dgibi ************************************************************************ * Section : Mecanique Elastique ************************************************************************ * Test Beton.dgibi: Jeux de données * * --------------------------------- * * * * CAS TEST DU 93/01/19 PROVENANCE : NAH OPTION ECHO 1; GRAPH='N'; SAUT PAGE ; * *-----------------------------------------------------* * * * Exemple simple de calcul mécanique élastique * * plastique beton fonctionnant avec les nouvelles * * structures MMODEL et MCHAMELEM * * * * Un cube en béton soumis à un chargement de * * traction uniaxial suivant l'axe OZ * * * * La limite en traction simple du béton LTR1 est * * égale à 4 MPa * * * *-----------------------------------------------------* * Définition des options * OPTI DIME 3 ELEM CUB8 MODEL TRIDIM ; *opti impi 9 ; N = 1 ; * * Maillage d'un cube de N*N*N éléments CUB8 * DENS N ; * A1 = 0 0 0 ; A2 = N 0 0 ; A3 = N N 0 ; A4 = 0 N 0 ; B1 = 0 0 N ; B2 = N 0 N ; B3 = N N N ; B4 = 0 N N ; * A1A2 = A1 D N A2 ; A2A3 = A2 D N A3 ; A3A4 = A3 D N A4 ; A4A1 = A4 D N A1 ; B1B2 = B1 D N B2 ; B2B3 = B2 D N B3 ; B3B4 = B3 D N B4 ; B4B1 = B4 D N B1 ; A1B1 = A1 D N B1 ; A2B2 = A2 D N B2 ; A3B3 = A3 D N B3 ; A4B4 = A4 D N B4 ; * A1A2A3A4 = A1A2 A2A3 A3A4 A4A1 DALLE PLAN ; B1B2B3B4 = B1B2 B2B3 B3B4 B4B1 DALLE PLAN ; * A1A2B2B1 = A1A2 A2B2 (INVE B1B2) (INVE A1B1) DALLE PLAN ; A4A3B3B4 = (INVE A3A4) A3B3 B3B4 (INVE A4B4) DALLE PLAN ; * A2A3B3B2 = A2A3 A3B3 B2B3 (INVE A2B2) DALLE PLAN ; A1A4B4B1 = (INVE A4A1) A4B4 B4B1 (INVE A1B1) DALLE PLAN ; * TOTAL = PAVE A1A2A3A4 B1B2B3B4 A1A2B2B1 A4A3B3B4 A2A3B3B2 A1A4B4B1 ; * *TRAC TOTAL CACH (1000 -2000 1000) ; * * * Formulation, matériau défini à l'aide des * opérateurs MODE et MATR * MO = MODE TOTAL MECANIQUE ELASTIQUE PLASTIQUE BETON CUB8 ; MA = MATE MO YOUN 40000 NU 0.20 ALPH 1E-5 LTR1 4 ; * * Conditions de blocages * ENC1 = BLOQ UX A1A4B4B1 ; ENC2 = BLOQ UY A1A2B2B1 ; ENC3 = BLOQ UZ A1A2A3A4 ; ENC = ENC1 ET ENC2 ET ENC3 ; * * Blocage pour déplacements imposés * * END1 = BLOQ UX A2A3B3B2 ; * END2 = BLOQ UY A4A3B3B4 ; END3 = BLOQ UZ B1B2B3B4 ; * END = END1 ET END2 ET END3 ; END = END3 ; EN = ENC ET END ; * * * Valeur des déplacements imposés * * FEXT1 = DEPI END1 1E-4 ; * FEXT2 = DEPI END2 1E-4 ; FEXT3 = DEPI END3 1E-4 ; * FEXT = FEXT1 ET FEXT2 ET FEXT3 ; FEXT = FEXT3 ; * * Force en fonction du temps * lt1 = prog 0. 10. ; lt2 = prog 0. 10. ; ev = evol rouge manu temps lt1 force lt2 ; chat = char DIMP fext ev ; * * Préparation à nonlin * ltt = prog 0 1 2 ; tab1 = table; tab1.blocages_mecaniques = EN; tab1.caracteristiques = ma; tab1.modele = mo; tab1.chargement = chat; tab1.temps_calcules = ltt; tab1.hypothese_deformations = lineaire; pasapas tab1; * * Traitement des résultats * ptt = prog 0 ; psz = prog 0 ; ic = 0 ; repeter mabou 2 ; ic = ic + 1 ; dd = tab1 . deplacements . ic; ptt = ptt et ( prog ( extr dd uz b3 ) ) ; ss = tab1.contraintes.ic ; sz = extr ( chan chpo mo ss ) smzz b3 ; psz = psz et ( prog sz ) ; fin mabou ; * evsz = evol vert manu deplacement ptt contrainte psz ; *dess evsz ; * * Code de fonctionnement * si ( (abs ( sz - 2. )) < 1.E-6 ) ; erre 0 ; sino ; list sz; erre 5 ; finsi ; * fin ;