********************************************* *Etude loi de comportement plastique cinematique* * Plaque de plâtre BA20 * * * ********************************************** *----------Donnée----------------- Lo= 0.500; TITR 'traction BA20 2D'; OPTI DIME 2 MODE PLAN DEFO ELEM QUA4; opti trac open; *----------Maillage----------------- * S1: plan inférieur A= 0. 0.; B= Lo 0.; C1 = Lo 0.0125; D1 = 0. 0.0125; dens1=300; AB = DROI dens1 A B; BC = DROI 8 B C1; CD = DROI dens1 C1 D1; DA = DROI 8 D1 A; S1 = DALL AB BC CD DA; F = point s1 proc (0. 0.0125/2); E = point s1 proc (Lo. 0.0125/2); trac s1; *-----------modèle de comportement du matériaux platre------ M1 = mode S1 MECANIQUE ELASTIQUE PLASTIQUE CINEMATIQUE ; MAT1 = mate M1 youn 2.82E9 nu 0.3 RHO 845. ALPH 1.75E-6 SIGY 4.4E6 H 2.82E8; *-----------------conditions limites-------------------- CL1 = BLOQ UX DA; CL2 = BLOQ UY F; CL3 = BLOQ UY E; CLTOT = CL1 et CL2 et CL3; *-------------imposer une forCe pr un de deplacement F0 = -6.E6; F1 = PRES MASS M1 F0 BC; LI1 = PROG 0. 20. 40.; LI2 = PROG 0. 1. 0.; EVF1 = EVOL MANU 'Temps' LI1 'P(t)' LI2; VEC1 = VECT F1 1.E-5 FX FY ROUG; * TRAC VEC1 S1; CHAMECA = CHAR 'MECA' F1 EVF1; * -------- RESOLUTION PASAPAS----------------* TAB1 = table; TAB1. 'MODELE' = M1; TAB1. 'CARACTERISTIQUES' = MAT1; TAB1. 'BLOCAGES_MECANIQUES' = CLTOT; TAB1. 'CHARGEMENT' = CHAMECA; TAB1. 'CONVERGENCE_FORCEE' = faux; TAB1. 'GRANDS_DEPLACEMENTS' = vrai; TAB1.'TEMPS_CALCULES' = PROG 0. PAS 0.05 40. ; TAB1. 'GRANDS_DEFORMATIONS' = vrai; * TAB1. 'PRECFLEX' = 100.; * TAB1. 'PRECISION' = 100.; * TAB1. 'MAXITERATION' = 100.; PASAPAS TAB1; toto0=defo (tab1.deplacements.0) S1; toto=defo (tab1.deplacements.20) S1 vert; toto1=defo (tab1.deplacements.40) S1 roug; trac (toto0 et toto et toto1); depref=prog; time=prog; pref= S1 POIN PROC (0.5 0.0125/2); toto=dime tab1.temps; i=0; j=0; repe bouc0 (toto); timi=tab1.temps.j; time=time inse (i + 1) timi; depi=extr (tab1.deplacements.i) UX pref; depref=depref inse (i + 1) depi; i=i+1; j=j+1; fin bouc0; evdep=evol manu time depref; * dess evdep; sigmref=prog; epsref=prog; toto=dime tab1.temps; i=0; repe bouc0 (toto); sigmi=exco (tab1.contraintes.i) SMXX; sigmi=chan chpo sigmi M1; sigmi=extr sigmi SMXX pref; sigmref=sigmref inse (i + 1) sigmi; epsii=EPSI M1 (tab1.deplacements.i); epsii=exco epsii EPXX; epsii=chan chpo epsii M1; epsii=extr epsii EPXX pref; epsref=epsref inse (i + 1) epsii; i=i+1; fin bouc0; ldc=evol manu 'eps' epsref 'sigm' sigmref; dess ldc; @EXCEL1 ldc 'loi de comportement.csv'; @EXCEL1 evdep 'déplacement centre bc.csv'; * sauver le résultat OPTI SAUV 'BA20 2D.resu'; SAUV; fin;