* fichier : plas7.dgibi ************************************************************************ * Section : Mecanique Plastique ************************************************************************ * Test Plas7.dgibi: Jeux de données * * --------------------------------- * * * * CAS TEST DU 91/06/13 PROVENANCE : TEST ******************************************************* * Test plas7.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 X ; SINO ; OPTI ECHO 0 ; FINSI ; SAUT PAGE; * ******************************************************* * * * TEST PLAS7 * * * * RECTANGLE EN FLEXION PURE (DEFORMATIONS 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 ECHO 1; * OPTI DIME 2; OPTI MODE PLAN DEFO; OPTI ELEM QUA4; * 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 'PLAS7 : MAILLAGE'; TRAC 'QUAL' SUR1; FINSI; * NBNO1 = NBNO SUR1; NBEL1 = NBEL SUR1; * *+++++++++++++++++++++++++++++++++++++++++++++++++++++* * definition du modele, des caracteristiques * *+++++++++++++++++++++++++++++++++++++++++++++++++++++* * MOD1 = MODE SUR1 MECANIQUE ELASTIQUE PLASTIQUE PARFAIT QUA4; 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 = TABLE; * * On prend comme liste des pseudo-temps * le deplacement u. * LREEL1 = PROG 0.029435 0.1 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 = 90.57; MOM2.2 = 135.2; MOM2.3 = 138.4; MOM2.4 = 139.0; MOM2.5 = 139.2; MOM2.6 = 139.3; MOM2.7 = 139.4; * MOM2.8 = 543.4; MOM2.9 = 557.7; *Boucle pour chaque pas de chargement : fonction de UA. I = 0; REPETER BLOC1 7; I = I + 1; CONT1 = TAB1.CONTRAINTES.I; DEPL1 = TAB1.DEPLACEMENTS.I; CHPO10 = CHAN CHPO MOD1 CONT1; * Trace facultatif de la deformee * (verification du chargement) SI (I EGA 7); SI (NEG GRAPH 'N'); TITR 'PLAS7 : DEFORMEE Increment :' I; 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 = EXTRA CHPO10 SMXX POI1; LIST2 = LIST2 ET (PROG SM); FIN BLOC4; LIST3 = LIST2 * LISD2; TITR 'CONTRAINTE SMXX Increment :' I; ESMXX = EVOL MANU Y LISD2 SMXX LIST2; TITR 'PRODUIT Y * SMXX Increment :' I; EYSMXX = EVOL MANU Y LISD2 YSMXX LIST3; SI (I EGA 7); SI (NEG GRAPH 'N'); DESS ESMXX; DESS EYSMXX; FINSI; FINSI; MOM1.I = INTG EYSMXX; MOM1.I = -1.E-4 * MOM1.I; * Calcul de l'ecart entre theorie et calcul. ERG.I = 100. * (ABS ((MOM1.I - MOM2.I) / MOM2.I)); * Extraction de la contrainte maximale SMXX * pour les pas 1 et 7. SI (I EGA 1); MOM1.8 = EXTR CHPO10 SMXX PA; MOM1.8 = MOM1.8 / 1.E6; ERG.8 = 100. * (ABS ((MOM1.8 - MOM2.8) / MOM2.8)); FINSI; SI (I EGA 7); MOM1.9 = EXTR CHPO10 SMXX (D1 POIN 2); MOM1.9 = MOM1.9 / 1.E6; ERG.9 = 100. * (ABS ((MOM1.9 - MOM2.9) / MOM2.9)); FINSI; SI (NEG GRAPH 'N'); TITR 'SMXX Increment :' I; TRAC MOD1 (EXCO CONT1 SMXX); FINSI; FIN BLOC1; *+++++++++++++++++++++++++++++++++++++++++++++++++++++* * affichage des resultats * *+++++++++++++++++++++++++++++++++++++++++++++++++++++* MESS 'RESULTATS : PLAS7'; MESS '-------------------'; SAUT LIGN; MESS 'Moment de flexion du aux contraintes SMXX le long de AB.'; MESS '--------------------------------------------------------'; SAUT LIGN; MESS ' UAmax = 0.029435 mm :'; MESS ' ------------------- '; MESS ' Moment theorique :' MOM2.1 'kN.m'; MESS ' Moment calcule :' MOM1.1 'kN.m'; MESS ' Soit un ecart de :' ERG.1 '%'; SAUT LIGN; MESS ' UAmax = 0.1 mm :'; MESS ' -------------- '; MESS ' Moment theorique :' MOM2.2 'kN.m'; MESS ' Moment calcule :' MOM1.2 'kN.m'; MESS ' Soit un ecart de :' ERG.2 '%'; SAUT LIGN; MESS ' UAmax = 0.2 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM2.3 'kN.m'; MESS ' Moment calcule :' MOM1.3 'kN.m'; MESS ' Soit un ecart de :' ERG.3 '%'; SAUT LIGN; MESS ' UAmax = 0.3 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM2.4 'kN.m'; MESS ' Moment calcule :' MOM1.4 'kN.m'; MESS ' Soit un ecart de :' ERG.4 '%'; SAUT LIGN; MESS ' UAmax = 0.4 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM2.5 'kN.m'; MESS ' Moment calcule :' MOM1.5 'kN.m'; MESS ' Soit un ecart de :' ERG.5 '%'; SAUT LIGN; MESS ' UAmax = 0.5 mm :'; MESS ' --------------'; MESS ' Moment theorique :' MOM2.6 'kN.m'; MESS ' Moment calcule :' MOM1.6 'kN.m'; MESS ' Soit un ecart de :' ERG.6 '%'; SAUT LIGN; MESS ' UAmax = 1.0 mm :'; MESS ' ---------------'; MESS ' Moment theorique :' MOM2.7 'kN.m'; MESS ' Moment calcule :' MOM1.7 'kN.m'; MESS ' Soit un ecart de :' ERG.7 '%'; SAUT LIGN; MESS 'Contrainte maximale smxx :'; MESS '------------------------ :'; SAUT LIGN; MESS ' UAmax = 0.029435 mm'; MESS ' -------------------'; MESS ' Valeur theorique :' MOM2.8 'MPa'; MESS ' Valeur calculee :' MOM1.8 'MPa'; MESS ' Soit un ecart de :' ERG.8 '%'; SAUT LIGN; MESS ' UAmax = 1.0 mm'; MESS ' --------------'; MESS ' Valeur theorique :' MOM2.9 'MPa'; MESS ' Valeur calculee :' MOM1.9 'MPa'; 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