opti echo 1; ******************* cette ligne fait 72 caractères ********************* * modele2.dgibi * - sans frottement * - sans la prise en compte du contact poincon/poudre * - uniaxial * - profil droit de la matricice * - sans aiguille * - profil droit des poincons * - modele Cam Clay * titr 'cas1'; opti dime 2 elem qua4 mode plan; *Figures 'O' ou 'N' graph='O'; trait='O'; force1=1.; press1=100.; * p1=-.005 -0.01; p2=0.005 -0.01; p3=-0.005 0.01; p4=0.005 0.01; l12=p1 droi 40 p2; l24=p2 droi 40 p4; l43=p4 droi 40 p3; l31=p3 droi 40 p1; surf1=dall l12 l24 l43 l31; * mo1=mode surf1 mecanique elastique plastique cam_clay qua4; * mat1=mate mo1 youn 1000. nu 0.2 e0 2. cohe 0.5 m 1.2 p0 5. lamd 0.1 kapa 0.01 g1 416.666667; * cl1=bloq l12 uy; cl2=bloq l24 ux; cl3=bloq l31 ux; cltot1=cl1 et cl2 et cl3; * *force2=(-1)*force1; *fo1=fsur mass (0. force2) mo1 l43; fo1=press mass press1 mo1 l43; listt1=prog 0 pas 1. 20.; listf1=prog 0 pas 0.1 1. pas -0.1 0.; evf1=evol manu 'temps' listt1 'force' listf1; cha1=char meca fo1 evf1; * si(neg graph 'N'); TRACE qual surf1; DESS (EVF1*force1) ; vectf1=vect fo1 1. fx fy jaun; trac vectf1 (l12 et l24 et l43 et l31); finsi; * tcalc1=prog 0. pas 1. 20.; * tab1=table; tab1.blocages_mecaniques=cltot1; tab1.modele=mo1; *tab1.pas_ajuste=vrai; tab1.chargement=cha1; tab1.caracteristiques=mat1; tab1.temps_calcules=tcalc1; tab1.grands_deplacements=vrai; *tab1.grandes_deformations=vrai; pasapas tab1; * opti sauv '/home/gv219722_vrh012860/Desktop/Simulations/save1.sortgibi'; sauv tab1; * list TAB1.ERREUR; list TAB1.CONV; * si(neg trait 'N'); geotot=extr tab1.modele mail; d=dime tab1.deplacements; i=0; repeter fig d; dep1=tab1.deplacements.i; def0=defo geotot dep1 0. blan; def1=defo geotot dep1 1. rouge; titr 'Graph' i; trac (def0 et def1); i=i+1; fin fig; finsi; fin;