* fichier : compar_syco_plast.dgibi ************************************************************************ * Section : Mecanique Viscoplastique ************************************************************************ ******************************************************* **------------------------------------------------------------------- MESS ' COMPARAISON ' ; MESS ' CALCUL DYNAMIQUE COMPORTEMENT VISCOPLASTIQUE SYCO1 '; MESS ' avec parametre p de la loi SYCO1 valant 1 '; MESS ' CALCUL DYNAMIQUE COMPORTEMENT ELASTOPLASTIQUE '; MESS ' en 2D defo plane ' ; MESS ' Unites: N Pa m kg s' ; MESS '' ; **-------------------------------------------------------------------- *** Options de calcul opti dime 2 elem QUA4 mode plan defo ; ******************************************************* *** maillage pa1 = 0. 0.; pa2 = 5. 0.; pa3 = 5. 3.; pa4 = 0. 3.; li12 = pa1 droi 1 pa2; li23 = pa2 droi 1 pa3; li34 = pa3 droi 1 pa4; li41 = pa4 droi 1 pa1; su1 = dall li12 li23 li34 li41; elim 0.001 su1; * ********************************************************** *** MODELE & MATERIAU *** ********************************************************** Epsprog0 = ('PROG' 0. 0.01 0.1 0.5 ); Sigprog0 = ('PROG' 0. 500. 1000. 1200.) * 1.e6; lcxx = 'EVOL' 'MANU' 'Epsilon' Epsprog0 'Sigma' Sigprog0 ; * ab1 = extr lcxx 'ABSC'; ord1 = extr lcxx 'ORDO'; eey1 = (extr ord1 2)/(extr ab1 2); lsm1 = ord1 enle 1 ; lep1 = (ab1 enle 1) - (lsm1 / eey1) ; evec = evol vert manu eps lep1 sig lsm1 ; *dess (lcxx et evec) titr ' Courbes de traction et d ecrouissage (vert)' ; * nu0 = 0.3 ; rho0= 7800.; ********************************************************** *** Modele & materiau *** * * elements standards mod1 = 'MODELISER' su1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'VISCOPLASTIQUE' 'SYCO1' 'QUA4'; mat1 = 'MATERIAU' mod1 'YOUN' eey1 'NU' nu0 'RHO' rho0 'ECRO' evec 'PSYC' 1. 'DSYC' 1.d8 ; * ********************************************************** *** CL et DEPLACEMENTS IMPOSES *** ********************************************************** cl1 = bloq (pa1 et pa4) 'UX'; cl2 = bloq (pa1 et pa2) 'UY'; cl3 = bloq (pa2 et pa3) 'UX'; cltot = (cl1 'ET' cl2 'ET' cl3 ) ; * *... chargement en deplacement impose umax = 1.; uimp = depi cl3 umax; * xx1 = prog 0. 1.; yy1 = prog 0. 50.; evochar = evol manu 'temps' xx1 'charge' yy1; char1 = char meca uimp evochar; * * *=======================================================* * * * RESOLUTION * * * *=======================================================* xx0 = prog 0. 1.15e-3 1.25e-2; * tab1 = table; tab1.dynamique = vrai; tab1.modele = mod1; tab1.caracteristiques = mat1; tab1.blocages_mecaniques = cltot; tab1.chargement = char1; tab1.temps_calcules = xx0; tab1.precision = 1.e-5; pasapas tab1; * *post traitement sig1 = tab1.contraintes. 2; ut1 = tab1.deplacements. 2; vt1 = tab1.vitesses. 2; tt1 = tab1.temps. 2; var1 = tab1.variables_internes. 2; * * val1 = extr var1 'EPSE' 1 1 3; val2 = extr var1 'VP' 1 1 3; val3 = extr ut1 'UY' pa3; ******************************************************* * ******************************************************* **------------------------------------------------------------------- MESS '' ; MESS ' CALCUL DYNAMIQUE COMPORTEMENT ELASTOPLASTIQUE '; MESS ' en 2D defo plane ' ; MESS ' Unites: N Pa m kg s' ; MESS '' ; **-------------------------------------------------------------------- *** Options de calcul opti dime 2 elem QUA4 mode plan defo; ******************************************************* *** maillage pa1 = 0. 0.; pa2 = 5. 0.; pa3 = 5. 3.; pa4 = 0. 3.; li12 = pa1 droi 1 pa2; li23 = pa2 droi 1 pa3; li34 = pa3 droi 1 pa4; li41 = pa4 droi 1 pa1; su1 = dall li12 li23 li34 li41; elim 0.001 su1; * ********************************************************** *** MODELE & MATERIAU *** ********************************************************** nu0 = 0.3 ; rho0= 7800.; ********************************************************** *** Modele & materiau *** * * elements standards * mod1 = 'MODELISER' su1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'PLASTIQUE' ; mat1 = 'MATERIAU' mod1 'YOUN' eey1 'NU' nu0 'RHO' rho0 'ECRO' evec ; ********************************************************** *** CL et DEPLACEMENTS IMPOSES *** ********************************************************** cl1 = bloq (pa1 et pa4) 'UX'; cl2 = bloq (pa1 et pa2) 'UY'; cl3 = bloq (pa2 et pa3) 'UX'; cltot = (cl1 'ET' cl2 'ET' cl3 ) ; * *... chargement en deplacement impose umax = 1.; uimp = depi cl3 umax; * xx1 = prog 0. 1.; yy1 = prog 0. 50.; evochar = evol manu 'temps' xx1 'charge' yy1; char1 = char meca uimp evochar; * *=======================================================* * * * RESOLUTION * * * *=======================================================* xx0 = prog 0. 1.15e-3 1.25e-2; * tab2 = table; tab2.dynamique = vrai; tab2.modele = mod1; tab2.caracteristiques = mat1; tab2.blocages_mecaniques = cltot; tab2.chargement = char1; tab2.temps_calcules = xx0; tab2.precision = 1.e-5; pasapas tab2; * *post traitement sig1 = tab2.contraintes. 2; ut1 = tab2.deplacements. 2; vt1 = tab2.vitesses. 2; tt1 = tab2.temps. 2; var1 = tab2.variables_internes. 2; * * xal1 = extr var1 'EPSE' 1 1 3; xal3 = extr ut1 'UY' pa3; aa1 = ABS(val1 - xal1); aa3 = ABS(val3 - xal3); * mess 'aa1, aa3 =' aa1 aa3 ; si (aa1 >EG 1.e-8); erre 5; finsi; si (aa3 >EG 1.e-8); erre 5; finsi; fin;