* fichier : plas6.dgibi ************************************************************************ * Section : Mecanique Plastique ************************************************************************ * Test Plas6.dgibi: Jeux de données * * --------------------------------- * * * * CAS TEST DU 91/07/23 PROVENANCE : BIRET ******************************************************* * Test plas6.dgibi: Jeux de données * * --------------------------------- * ******************************************************* * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE * SINON SI GRAPH DIFFERENT DE N TOUS * LES GRAPHIQUES SONT AFFICHES GRAPH = 'N' ; SAUT PAGE; SI (NEG GRAPH 'N') ; OPTI ECHO 1 ; OPTI TRAC PSC ; SINO ; OPTI ECHO 0 ; FINSI ; SAUT PAGE; ******************************************************* * * * TEST PLAS6 * * * * RECTANGLE EN FLEXION PURE (CONTRAINTES PLANES, * * PLASTICITE PARFAITE) * * * * Test V.P.C.S. numero SSNP11/89 * * * * References : R. Hill, Mathematical theory of * * plasticity, Oxford, Oxford Clavendon Press, 1950. * * * * Description : * * ----------- * * * * axe y * * D | A * * |----- ---- UA epaisseur a = 12.5 mm * * | | / hauteur b = 100.0 mm * * | | / * * | | / * * | |/ * * O |-----/-----> axe x * * | /| * * | / | * * | / | * * | / | * * ----- * * C UB B * * * * Chargement : Deplacement impose variant * * lineairement de UA a -UA entre A et B, * * UA croissant jusqu'a : UAmax = 0.5 mm . * * * * Conditions aux limites : * * sur CD : u = 0. * * au point O : v = 0. * * * * Proprietes des materiaux : * * Materiau elastique, parfaitement plastique * * E = 210 000 MPa , NU = 0.3 , SIGY = 483 MPa * * * ******************************************************* OPTI DIME 2; OPTI MODE PLAN CONT; OPTI ELEM QUA4; GRAPH2 = GRAPH; TEMPS; ******************************************************* * geometrie : maillage * ******************************************************* PA = 12.5E-3 50.E-3; PB = 12.5E-3 -50.E-3; PC = 0. -50.E-3; PD = 0. 50.E-3; NB1 = 32; NB2 = 4; D1 = PA DROI NB1 PB; D2 = PB DROI NB2 PC; D3 = PC DROI NB1 PD; D4 = PD DROI NB2 PA; SUR1 = DALL D1 D2 D3 D4 PLAN; SI (NEG GRAPH 'N'); TITR 'PLAS6 : MAILLAGE'; TRAC QUAL SUR1; FINSI; NBNO1 = NBNO SUR1; NBEL1 = NBEL SUR1; ******************************************************* * calcul du modele,caracteristiques * ******************************************************* MOD1 = MODE SUR1 MECANIQUE ELASTIQUE PLASTIQUE PARFAIT; MAT1 = MATE MOD1 YOUN 2.1E11 NU 0.3 SIGY 483.E6; ******************************************************* * conditions aux limites * ******************************************************* CL1 = BLOQ UX D3; CL2 = BLOQ UY (SUR1 POIN PROC (0. 0.)); CL3 = BLOQ UX D1; CL = CL1 ET CL2 ET CL3; * Creation du chpoint de deplacement impose variable : * deplacement impose variant lineairement de u a -u * entre A et B, u croissant jusqu'a 0.5 mm . I = 0.; LISD1 = PROG 1.; REPETER BLOC10 NB1; I = I + 1.; YY = 1. - (2. * (I / NB1)); LISD1 = LISD1 ET (PROG YY); FIN BLOC10; LISD1 = LISD1 * 1.E-3; CHPO1 = MANU CHPO (CHAN 'POI1' D1) 1 UX LISD1; CHA1 = DEPI CL3 CHPO1; ******************************************************* * objets relatifs a la procedure PASAPAS * ******************************************************* TAB1 = TABL ; *On prend comme liste des pseudo-temps le deplacement u LREEL1= PROG 0.02875 0.05 0.1 0.15 0.2 0.3 0.4 0.5 1.0; * Chargement EVOL1 = EVOL MANU UA LREEL1 COEFF LREEL1; CHAR1 = CHAR DIMP CHA1 EVOL1; ******************************************************* * resolution * ******************************************************* TAB1.'CARACTERISTIQUES' = MAT1; TAB1.'CHARGEMENT' = CHAR1; TAB1.'MODELE' = MOD1; TAB1.'TEMPS_CALCULES' = LREEL1; TAB1.'BLOCAGES_MECANIQUES' = CL; PASAPAS TAB1; ******************************************************* * extraction des resultats * ******************************************************* * Initialisation MOM1 = TABL ; LISD2 = LISD1 * 50.; MOM2 = TABL ; ERG = TABL ; MOM2.1 = 805.0; MOM2.2 = 1074.; MOM2.3 = 1174.; MOM2.4 = 1193.; MOM2.5 = 1199.; MOM2.6 = 1204.; MOM2.7 = 1205.; MOM2.8 = 1206.; MOM2.9 = 1207.5; *Boucle pour chaque pas de chargement : fonction de UA. I = 0; REPETER BLOC1 9; I = I + 1; CONT1 = TAB1.CONTRAINTES.I; DEPL1 = TAB1.DEPLACEMENTS.I; * Passage du champs de contrainte en chpoints * (utilisation des nouveaux chamelems) CHPO10 = CHAN CHPO MOD1 CONT1; * Trace facultatif de la repartition * de la contrainte SMXX SI (I EGA 9); SI (NEG GRAPH2 'N'); TITR 'PLAS6 : CONTRAINTE SMXX'; BB = EXCO 'SMXX' CONT1; TRAC BB MOD1; FINSI; FINSI; * Trace facultatif de la deformee * (verification du chargement) SI (I EGA 9); SI (NEG GRAPH2 'N'); TITR 'PLAS6 : DEFORMEE'; DEF0 = DEFO SUR1 DEPL1 0.; DEF1 = DEFO SUR1 DEPL1 2. VERT; TRAC SUR1 (DEF0 ET DEF1); FINSI; FINSI; * Calcul du moment de flexion du aux contraintes * smxx s'exercant le long de AB. LIST2 = PROG; J = 0; NB = NB1 + 1; REPETER BLOC4 NB; J = J + 1; D = EXTR LISD2 J; POI1 = SUR1 POIN PROC (12.5E-3 D); SM = EXTR CHPO10 SMXX POI1; LIST2 = LIST2 ET (PROG SM); FIN BLOC4; LIST3 = LIST2 * LISD2; TITR 'CONTRAINTE SMXX'; ESMXX = EVOL MANU Y LISD2 SMXX LIST2; TITR 'PRODUIT Y * SMXX'; EYSMXX = EVOL MANU Y LISD2 YSMXX LIST3; SI (I EGA 9); SI (NEG GRAPH2 'N'); DESS ESMXX; DESS EYSMXX; FINSI; FINSI; MOM1.I = INTG EYSMXX; MOM1.I = -1.E-3 * MOM1.I; ERG.I = 100. * (ABS ((MOM1.I - MOM2.I) / MOM2.I)); FIN BLOC1; ******************************************************* * affichage des resultats * ******************************************************* MESS 'RESULTATS : PLAS6'; MESS '-------------------'; SAUT LIGN; MESS 'Moment de flexion du aux contraintes SMXX le long de AB.'; MESS '--------------------------------------------------------'; SAUT LIGN; MESS ' UAmax = 0.02875 mm :'; MESS ' ------------------'; MESS ' Moment theorique :' MOM1.1 'N.m'; MESS ' Moment calcule :' MOM2.1 'N.m'; MESS ' Soit un ecart de :' ERG.1 '%'; SAUT LIGN; MESS ' UAmax = 0.05 mm :'; MESS ' --------------- '; MESS ' Moment theorique :' MOM1.2 'N.m'; MESS ' Moment calcule :' MOM2.2 'N.m'; MESS ' Soit un ecart de :' ERG.2 '%'; SAUT LIGN; MESS ' UAmax = 0.1 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM1.3 'N.m'; MESS ' Moment calcule :' MOM2.3 'N.m'; MESS ' Soit un ecart de :' ERG.3 '%'; SAUT LIGN; MESS ' UAmax = 0.15 mm :'; MESS ' ---------------'; MESS ' Moment theorique :' MOM1.4 'N.m'; MESS ' Moment calcule :' MOM2.4 'N.m'; MESS ' Soit un ecart de :' ERG.4 '%'; SAUT LIGN; MESS ' UAmax = 0.2 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM1.5 'N.m'; MESS ' Moment calcule :' MOM2.5 'N.m'; MESS ' Soit un ecart de :' ERG.5 '%'; SAUT LIGN; MESS ' UAmax = 0.3 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM1.6 'N.m'; MESS ' Moment calcule :' MOM2.6 'N.m'; MESS ' Soit un ecart de :' ERG.6 '%'; SAUT LIGN; MESS ' UAmax = 0.4 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM1.7 'N.m'; MESS ' Moment calcule :' MOM2.7 'N.m'; MESS ' Soit un ecart de :' ERG.7 '%'; SAUT LIGN; MESS ' UAmax = 0.5 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM1.8 'N.m'; MESS ' Moment calcule :' MOM2.8 'N.m'; MESS ' Soit un ecart de :' ERG.8 '%'; SAUT LIGN; MESS ' UAmax = 1.0 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM1.9 'N.m'; MESS ' Moment calcule :' MOM2.9 'N.m'; MESS ' Soit un ecart de :' ERG.9 '%'; SAUT LIGN; MESS 'Donnees sur le maillage'; MESS '-----------------------'; MESS ' Nombre d elements :' NBEL1 'QUA4'; MESS ' Nombre de noeuds :' NBNO1 ; SAUT LIGN; ******************************************************* * code fonctionnement * ******************************************************* ERGMAX = MAXI (PROG ERG.1 ERG.2 ERG.3 ERG.4 ERG.5 ERG.6 ERG.7 ERG.8 ERG.9); SI (ERGMAX