* Exercice de non-linearite en contact : opti echo 1; opti dime 2 elem tri3 mode plan deformation; ** Donnees : e = 20000e6; poi = 0.2; ** Geometrie du cylindre : r = 5.0e-3; ** Geometrie du plan : lx = 15.0e-3; ly = 10.0e-3; ** Parametres du chargement : w = 70000; q = 0.; ** Calcul de la longueur theorique de contact : v = (((1-(0.2*0.2))/e)+((1-(0.2*0.2))/e))**-1; contact = 2*((4*(w/(2*r))*r)/(3.14*v))**0.5; a = contact/2; ** creation des points du systeme (pour creer le maillage du systeme) : *b est un angle b = atg( (1.2*a) / (((r**2)-((1.2*a)**2))**0.5) ); r_1 = 1.1*(r-((cos b)*r)); r_2 = 1.1*(r-((cos (2*b))*r)); r_3 = 1.1*(r-((cos (3*b))*r)); ** Densite du maillage : message 'Donner la densite. Attention a donner un entier!'; message '(Ne pas choisir une densite inferieure a 16)'; obtenir ds_gen5*entier; ds_gen4 = ds_gen5/2; ds_gen3 = ds_gen5/4; ds_gen2 = ds_gen5/6; ds_gen1 = ds_gen5/10; ds_cer = ds_gen5/16; ** Calcul de la pression theorique de contact : p = 0.000005; pr0 = ((2*(w/(2*r))))/(3.14*a); list_r = prog (-1.0*a) pas p a; list_n = prog (dime list_r)*1.; list_r2 = (pr0/a)*((((a**2)*list_n)-(list_r**2))**0.5); list_r3 = prog 0 pas p (2*a); ev_f = evol vert 'MANU' distance list_r3 pression list_r2; *dess ev_f; *points du maillage du demi-cercle : p1 = ((r*(sin b))+(lx/2)) (lx-((cos b)*r)); p2 = (coor 1 p1) (ly + r_1); p3 = ((lx/2)-(r*(sin b))) (coor 2 p2); p4 = (coor 1 p3) (coor 2 p1); p5 = ((lx/2)-(r*(sin (1.5*b)))) (lx-((cos (1.5*b))*r)); p6 = (coor 1 p5) (ly + r_2); p7 = ((lx/2)+(r*(sin (1.5*b)))) (coor 2 p6); p8 = (coor 1 p7) (coor 2 p5); p9 = ((r*(sin (2.5*b)))+(lx/2)) (lx-((cos (2.5*b))*r)); p10 = (coor 1 p9) (ly + r_3); p11 = ((lx/2)-(r*(sin (2.5*b)))) (coor 2 p10); p12 = (coor 1 p11) (coor 2 p9); p13 = ((lx/2)-r) lx; p14 = (lx/2) lx; p15 = ((lx/2)+r) lx; *points du maillage du plan : p16 = (coor 1 p2) ly; p17 = (coor 1 p3) ly; p18 = (coor 1 p17) (ly-(r_1)); p19 = (coor 1 p16) (coor 2 p18); p20 = (coor 1 p7) ly; p21 = (coor 1 p7) (ly-r_2); p22 = (coor 1 p6) (coor 2 p21); p23 = (coor 1 p6) ly; p24 = (coor 1 p11) ly; p25 = (coor 1 p11) (ly-r_3); p26 = (coor 1 p10) (ly-r_3); p27 = (coor 1 p10) ly; p28 = lx ly; p29 = lx 0.; p30 = 0. 0.; p31 = 0. ly; p32 = (lx/2) ly; p33 = (lx/2) ly; ** creation des maillages : *Maillage du demi-cercle : l12 = droi ds_cer p1 p2; l23 = droi ds_gen5 p2 p3; l34 = droi ds_cer p3 p4; l432 = cerc ds_gen5 p4 p14 p32; l321 = cerc ds_gen5 p32 p14 p1; ct1 = l12 et l23 et l34 et l432 et l321; mail1 = surf ct1 plan; *trac mail1; l45 = cerc ds_gen4 p4 p14 p5; l56 = droi ds_cer p5 p6; l67 = droi ds_gen4 p6 p7; l78 = droi ds_cer p7 p8; l81 = cerc ds_gen4 p8 p14 p1; ct2 = l45 et l56 et l67 et l78 et l81 et l12 et l23 et l34; mail2 = surf ct2 plan; *trac mail2; l89 = cerc ds_gen3 p8 p14 p9; l910 = droi ds_cer p9 p10; l1011 = droi ds_gen3 p10 p11; l1112 = droi ds_cer p11 p12; l125 = cerc ds_gen3 p12 p14 p5; ct3 = l89 et l910 et l1011 et l1112 et l125 et l56 et l67 et l78; mail3 = surf ct3 plan; *trac mail3; l1213 = cerc ds_gen2 p12 p14 p13; l1314 = droi ds_gen2 p13 p14; l1415 = droi ds_gen2 p14 p15; l159 = cerc ds_gen2 p15 p14 p9; ct4 = l1213 et l1314 et l1415 et l159 et l910 et l1011 et l1112; mail4 = surf ct4 plan; *trac mail4; mail_ar = (mail1 et mail2 et mail3 et mail4) coul vert; trac mail_ar; *Maillage du plan : l1633 = droi ds_gen5 p16 p33; l3317 = droi ds_gen5 p33 p17; l1718 = droi ds_gen1 p17 p18; l1819 = droi ds_gen4 p18 p19; l1916 = droi ds_gen1 p19 p16; ct5 = l1633 et l3317 et l1718 et l1819 et l1916; mail5 = surf ct5 plan; l2021 = droi ds_gen1 p20 p21; l2122 = droi ds_gen3 p21 p22; l2223 = droi ds_gen1 p22 p23; l2317 = droi ds_gen3 p23 p17; l1620 = droi ds_gen3 p16 p20; ct6_2 = l2021 et l2122 et l2223 et l2317 et l1620 et l1718; ct6 = ct6_2 et l1819 et l1916; mail6 = surf ct6 plan; l2324 = droi ds_gen2 p23 p24; l2425 = droi ds_gen1 p24 p25; l2526 = droi ds_gen3 p25 p26; l2627 = droi ds_gen1 p26 p27; l2720 = droi ds_gen2 p27 p20; ct7_2 = l2324 et l2425 et l2526 et l2627 et l2720 et l2021; ct7 = ct7_2 et l2122 et l2223; mail7 = surf ct7 plan; l2728 = droi ds_gen1 p27 p28; l2829 = droi ds_gen1 p28 p29; l2930 = droi ds_gen1 p29 p30; l3031 = droi ds_gen1 p30 p31; l3124 = droi ds_gen1 p31 p24; ct8_2 = l2728 et l2829 et l2930 et l3031 et l3124 et l2425; ct8 = ct8_2 et l2526 et l2627; mail8 = surf ct8 plan; mail_pl = (mail5 et mail6 et mail7 et mail8) coul jaune; trac mail_pl; mail_cnt = (impo mail (l432 et l321) (l1633 et l3317)) coul rouge; *mail_fro = IMPF mail_cnt; ** Construction du modele : *mode_fro = model mail_fro frottement coulomb; *mate_fro = materiau mode_fro mu 0; modpar = model mail_ar mecanique elastique; matpar = materiau modpar youn e nu poi; modppla = model mail_pl mecanique elastique; matppla = materiau modppla youn e nu poi; ** Application du chargement : l_tmp = prog 0. 1.; FORY1 = FORC (0. (-1*W.)) (l1314 et l1415); VECT1 = VECT FORY1 FX FY bleu; *TRACE VECT1 (MAIL_AR ET l1314 et l1415) TITRE 'CHARGE'; LFORC = PROG 0. 1.; EVF1 = EVOL rouge MANU 'temps' L_TMP 'force' LFORC ; CHA1 = CHAR MECA FORY1 EVF1; ** Conditions aux limites : * cl_rig = rela 'ENSE' uy (l1314 et l1415); cl_pl = bloq UX UY l2930; cl_p2 = bloq UX p32; clt = cl_pl et cl_p2 et cl_rig; ** Resolution : p = 0.05; pas_cal =prog 0. pas p 1; tab1 = table; tab1 . 'TEMPS_CALCULES' = pas_cal; tab1 . 'MODELE' = modpar et modppla; tab1 . 'CARACTERISTIQUES' = matpar et matppla; tab1 . 'BLOCAGES_MECANIQUES' = clt; tab1 . 'CHARGEMENT' = CHA1; tab1 . 'CONTACT' = mail_cnt; *tab1 . 'FROTTEMENT' = mail_fro; tab1 . 'GRANDS_DEPLACEMENTS' = VRAI; * pasapas tab1; *OPTI SAUV JELE; *SAUV ; *OPTI rest JELE; *rest ; ** Animation : * @cartoon tab1 (mail_ar et mail_pl) 1. noscil; * ** Resultats : rf = peche tab1 'REACTIONS'; rl1 = redu rf (l1633 et l3317); rl2 = exco fy rl1; ev = evol rouge 'CHPO' rl2 scal (l1633 et l3317); dess ev XBOR 0 0.0006; dess (ev_f) XBOR 0 0.005; s = max rl2; list s; *('XBOR' -1.0 1.0) *AAA = tab1.'DEPLACEMENTS'. 1; *DEF1 = defo (l1633 et l3317) AAA 100. ; *TRAC (DEF1) ; *AAA = tab1.'DEPLACEMENTS'. 1; *DEF1 = defo (mail_ar et mail_pl) AAA 2 ; *TRAC (DEF1); *AAA2 = EXTR AAA UY P33; *LIST AAA2; rf = peche tab1 'DEPLACEMENTS'; rl1 = redu rf (l1633 et l3317 et l1620 et l2317); rl2 = exco UY rl1; ev = evol bleu 'CHPO' rl2 scal (l1633 et l3317 et l1620 et l2317); dess ev XBOR (-1.0*0.01) 0.01;