From 8265f7f7afaf7e3ed9d9f3947f04e59f38261248 Mon Sep 17 00:00:00 2001 From: rio <> Date: Tue, 10 Jan 2017 10:59:05 +0000 Subject: [PATCH] mise a jour du test DFC cf. Readme --- .../Test_R_dynamique_explicite/README | 25 +++++++++++ .../acceleration.pdf | Bin 0 -> 19705 bytes .../cube_dynamique_explicite.info | 9 +++- .../cube_dynamique_explicite_princ.maple.ref1 | 26 +++++++++-- .../Test_R_dynamique_explicite/plotPdf.gnu | 41 ++++++++++++++++++ .../Test_R_dynamique_explicite/vitesse.pdf | Bin 0 -> 19171 bytes 6 files changed, 96 insertions(+), 5 deletions(-) create mode 100644 Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/acceleration.pdf create mode 100644 Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/plotPdf.gnu create mode 100644 Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/vitesse.pdf diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/README b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/README index 8e54e6f..2e5b85f 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/README +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/README @@ -2,6 +2,10 @@ Auteur ------------------------------------------------------ Julien Troufflard (troufflard@univ-ubs.fr) +Gérard Rio (gerard.rio@univ-ubs.fr) +NB: GR: modification du test janvier 2017: +- remplacement du maple.ref1 complet: les données initiales étaient fausses +- inclusion de commentaires, figures dans le Readme ------------------------------------------------------ Mots-cles @@ -19,12 +23,33 @@ test simple du mode de calcul : dynamique explicite (DFC) TYPE_DE_CALCUL dynamique_explicite "" +L'objectif est ici d'observer le fonctionnement de l'algorithme classique DFC. +On observe le déplacement, vitesse, accélération transversalement au sens imposé. +Ces mouvements dépendent donc ici du coefficient de poisson, c-a-d -0.3 ce qui +se passe dans le sens imposé. + +Remarques: On observe l'apparition de fréquences dites numériques. +Pour supprimer ou atténuer ces fréquences numériques on pourra utiliser: + 1) soit des schémas numériques adaptés, type : + - en explicite Tchamwa et Chung lee, + - en implicite HHT, + 2) ou une loi de comportement modifiée par l'ajout d'un bulk viscosity. + +(cf. documentions théoriques et utilisateur d'herezh) + ------------------------------------------------------ Description du calcul ------------------------------------------------------ traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE) - loi elastique ISOELAS + - chargement suivant une vitesse constante (une rampe de déplacement) + +\figures: acceleration.pdf vitesse.pdf \legende: La première figure montre l'évolution de +l'accélération gamma2 du noeud 8. On observe que globalement on a bien une vitesse +transversale globalement = -0.3 la vitesse imposée. On observe également des oscillations +numériques, ce qui est classique en DFC.. La figure 2 montre le même type +d'évolution pour la vitesse \fin_legende ------------------------------------------------------------- Grandeurs de comparaison diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/acceleration.pdf b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/acceleration.pdf new file mode 100644 index 0000000000000000000000000000000000000000..e274cd2eca8f83bd38cf4015451ebd543cab10ff GIT binary patch literal 19705 zcmd_SWmsKX&_6~X0fGer1c!sWy9W>M?(Xic!QI_0IKkaL!QI_GXpr1YiJ? zWS{^50JWg8nVr53@YhV+PM=?2*HTX(ii-=%*3L#>+Z@U%wIz8de2EFh|Jnnpt(_$M z7&Mp)3h5P;U6g|Z4~lX0Yk0kP!YJ+!GH(e*cS3vf)seid@xp{fGR)5F4lP5}CfZg! zl~`=~ON;WXGM5K{42>!FI9!jX&f%HrX)nWT!hns@Hf zc}j4qcjhwU@(6VI1j_|}pfA0~rZb6#&1oP%#0S?NX@1#{67lZQ`C-qRr(n3J&_>bJ z3xTY0^$3cG1Wcv(^32)!+Xa9YsU$*j=#7j>JKw|2jsESeb&fm@tdoH(w2BgR983+W zQL?)}eSKo*T&Dt?FBU{ZFdn*bT*lfM%Bd6l3n3w=!k+etL-O>F9;A_Xwd-hyEKV?x z1v$(X{(Wlg8@V8<$ZJ{ez$B~6oxo9KIaFtZ4J$TQD^)UwR<-^6G}?82XowUzCUsVg zuaHMpb%G6{RnK4?Fgb>yj3+rR>!Izc5nH!5$F(XxO(Tn+Xo=*1EYWfyheL^QDBsre zd=bjKt16NI4otkN+zi-1i|VY%Q~$5`TwE!M1E8JW?~EZ}o6w#z+yxt2(e_k*EOi2B zR}e@@EYlDlJOjmf?X1!rpaR47we=iLP3+c8n>?fA7#?RQ)GmlZp^1yU%eTZ05Y-@X z^~g~cxyko+Y?1}8SQ)XnqXTB@UGxW71qvxBdD(=FYZmqT61Q0-yW1buB2)2z%WZl= zN=)2+b0xoOtf_C<+0aP?4HBZV8!Wj7;<4Zm^kS+X))|y4-=i__?{Z3^kfk!F73KF@ zOMOj#4-03ISsVG%q(aC0v$ou%YbKQuK@PjSc-Tf}GNU&-=9l;8xUVO*dO)|f3QMow zt4GeGv*rD}6nQf9GmtH3oe%*K-;SCf?|J?HyQ@M7S({H1gG>;zb$!ycePBhL%Yj$@ zlZVK*lQ~8!8*bVg06QoynoC#IDOZ-vrLHiWH)=}~fw9#&Rx5mV`UjtmiqM#{m)~W- zm1vFs;hHbJ0+o8H%>wmOa3WIfQ^lDleBN|0KsyD$Mmdy->q;56gNYrRzL-B$La@ug zK_L@7V`pkx)cjplP=vbyU&Yt(J_^*ey}ihR!xEiWV!3(Xqkg6kA@h?%Ty?Y>h6DHz zUdhb+h?~PAN)@D}-4?@>-g{eGl;%`}YRo7cH&x!4UfkUW6iKJz?Tqw^fVqp2eI3Zg z6dESf$rxk&a9Hf4g?0#%I#j)mT#U`ReOnc=9CTC>p(7TFhK9miQ;6~}Bj7fft8dQO zL@o9Ls_YM&`5ko;im=2!GMw4)lQ(1~GKvzhxtVY%4@(qV;Pmp!u{-S=XuxG+tlP(n z*)6A27#jA}gR?3`&3E82Uw)Rf2%S+r?>F8Xna*w-&c;~2u*JDLE_su=pK%3q*E)%g zuqI$E@!?wnwy}o$;&@~$#HllBn`a}2HgidUMKzDtW#;;R3k~62{EMQd6YSkl(zWq& z@lHJ9CW6pA=^QY3Za#ZQ1Wc9_Fl-#dzV5!I z2MzCw@hFM-PS=Al2^fgZ=@d=~ERPgNk z2)^8c)nIcnCT(`rXl~9Oz=&3wARSr9cPTqIuq?=RBLEqmWR{oMr9d29K*>1=$KG5b zoUCbpYfSBG%YRtEXn5bc^v;YO)07JOvI}O%*2&HP^<4WwbQi-Y*8Tze&IgR57l;o& zU0IL_yU@uWRBy(1oYxp-<^lHDpWeh00Ja9?QpXsRsTzjz>|``ayk+vh=L~Xq=elDU zVT@vg3K`EtQ+!#7F;y5Yz~6*QdQyYEQmyrif&N@U8`HfqX~X0j%5Na7S-zRbrGVmO zTaY;WLSKnzr^@>CD>(in>CGOmT^UDzdGOL)G6Z{)0utp`Ymm}@eHnLwi&0P#!ifgs zfo64^fT0k1j`4fzo&hdOxi2$(k3HDB#;&EK@;6^%6-NU(UShhrYpayK*s{~$Rc~vP zDk~AgQ~f}qxxV&!w;R*dM?^GH;qw=Y$nc6HiI=7Hk^o{Sh;44d;5-OwhXilD@1IZ6 zFZN-^rlhENA;bNP;s`U5_!Az?Qn7~{DfEi4lM3moVQ!~_*IU3P9cARkkoCQOxSx>+ z{~%~u{9e78vtZE4>JpS-4w^VFCBj(|<09u;ku{lfTznd|e&cY9(uvQL1*LDH_s8x8 z{Of67d)myNeu3K=EfWjIo`d$U ziTu}1$J0PRYtjN}X_#pk{(c*x;_N0j)$${)=}Sx1$(@*i1yPm0W@2RI2mi>hss3cZ zSj9`wK!869Bx)=pnS2zyH<^W3Cf#cu@WObqcXw*V_A?6PN{6N=2lhkTUthcp(rJM`n4wIw{bu9gO9yd)qr24PCNOlax#Z#g;yp<5h3R`9LXdB4 zcQD7{zLnOA$;c2~)Nn9F8YmD>STRaM$&c{H6s#1+6k9l*ukWF>1gBr82KCUjWqzRE z^-Uy8GCvFPG|0GhFiTCkSVF)Q&|@>dy)CeZ}6@E4+oU#-Bz>ZhEe~ zWEjc(sSU6~lK63Nv)IpP-dSytOpa(p8}krb$_E*0Ezjr(hEA4CbipW9Qe{)LIwX|& zc5B0~Zg#<){GJ(n2C0N}QD^0F33D;JDKB%OtZWO~xy}C^g$&=toe4u&w=~L(Q68^* zAYV5*`t&nLYj@Cvw5;%Q;5?rCo$&(Bcj}9S)@~TG@8sWM^PC^@^7t)D7-QdUMyMl) z>1{VQgF~0rZYJ_(4vzZ;cD4t+-|h`Plfu;*uJM1L7KFB=!{>LU1Fn1nPdlU%@y6JP zh$K7|@@8D$ASSP$Ce&XGK$l%#II6>BO*gF56Mdd%sI~B6e@q;!Ze@D9Mr%kHg1O97 zuP?oO>DY#VN!6PnflRX^@N{O#1ZGx?pW7|9(0RWuywL4Krnz-QH%V`ON^;_~kWm1gOYV z?j}TNggZzoykB2V9y4oT-1RpeS6_ETGa}VZHOg{$A71_(+A+iWz!+FTF%veiQc3bL z>MF8A8v#LYMm3CRCMKyygU;9##htmQw&65ACT4pa}@1 zXrs8Dqu6QicpjoOFD<|fVAgow(g&5@+oY>}X?+jRDm>PM9~BVV%|3b{i0-$>k#Cx@ zH+9=i#vf$DKd!=46GC3$rR-OwDtq(cL@%r#ZxZH3fF!n;{{Wq*^l@@x|QS z4x(i==;uHZ$;jIkkBftaGVUtV>-ddm4|yVPJ+*#r>eNt!DYROrLK*om>wN31D72p9 zNz6hdl`w<64kHGmY7O5gupD{qFp_lw&V8wsqOp3|oN5@vvvlj?OSK?fv{#pgaP|A;=%gOMG&LFK%OL3B_*GZH4TnyXiQuhsmn=BND?Z%4CP(FO?dm`vU3P9 zI0-2{d=1F+pBuAd&0C%F%T8zP*TW>lr|6sD?7lhtRf>qFV*K%vcgo6d zYTATAthHH+%j(@n%n%#*L6%)vCmP7WT0q9`7C$xrhp5*6!wK04LcGl>Z}v}3O)G9O zZbIf*^`srqraBk3Av2G3XEYx;C_c0XI9Q+{weJ-grka>Yh&pKu51tNUDl2-4tyUq` zQpBEELWv>NFiW^e$8cpg$(K@D3DPN1$z;~rzbd}4rsFD-kKLro;m|&DVf%nQ^$w2S zxXWz*qxXp2L!;Rt5i~xYmb_8*dl8tK<%NsR<0?+g93%5$^3Bc+kB<9UcrZqzWovm_B!uyxs?BwD zZ3?;KiTdOEw?1YuQX5wnHVt6{nfh!=z42L1Ube#!2`U4V#H7Qw*O$h{6HJj`kunBG zK2ukgSC(62)l@}FRELV{qOc^S;kLU}&)enT>FBkj$MH6b)W?xmQi#sHtU2H$+(WBr z)$Y$`$fS@9blOHM@jcF#i$G%#UTs8e8Z3W90*n4~)NfzQaFVS%p1Ht|ry%C=wT_}$ zA)u$ri|m~@UEq7;AR_^3%3jO+ub(dQ^X*DGH88Zo$KuF^cZ($zr!be;>hxNDr`2F^ zi2b*Qi1_oVl$#!gUt>ttwx!?fM-eN`)QVA1B{tk-`VAS&1=h{FdP8WwkzCIuDbhie z=&@3#ATilmpOa`~1SFqBmuISmhZ5n3TyQ!h#}q*5Ls|tbVefT=^`apIR*9h&KCvAh zPXe1X4TiLM5ED-jioU568I;p^?mFBKo~A%PtS>*0d3f?uLmi6L%&UsMR&W;dwW z<$}K&OWfprkbDMsT2x;knsp}YWwB2}R`j_5L7uxHEDs4DK8%1yP7g(AEOsDF1V>=F zCf#S3hY9cV^jGB%BP$~@<0)k2&79WTVb`BJ(;+P9DqM%)=uCx-=h^iZYhoZT$G;K8 z-k4qQ1X*f`a9pl&y{#v|YClfDJ9ZO`B{e;V78xD}#SmjH8kSET&s3A=e+bQ+6d7WG zYUwhu&-@7MtE@Xcj-8xMt#atlC!oEiLJQeUW1RWv!Va4(+QP>R*+Bb3usBwPtL_0Io zxf-)nJ<;|9*Ju#*aCr=1_x#A}H8E`~l`DKz=psYOU(kK0&4P=viMS2{h_R9rhW(a% z^uvy~*PE8E%a>_&>DQzqyCbNDa~$M)MvlshXnAzuRge}a1y-uD9u@_rZ3l@5SSTBH zMh{feWm48I)hNe+14x5S8=fXk_YFHR|V zAV%giS2)4S1iEQSaSiEJQ(wzen=4tcwMf^?TJmExnqs4t!MDwVenyQ1O zMT@AA8^0ilOcG>CNyE_x;#(U13}%r2dofB4F($zwmh_4KYK7&ZN)ul<_^jfp5h*IU zA~%nZT}M^qB)$76iwborz1kFO=M#bgkRL^G1m*Et!YeSKn{^`9JXt`lv`8c-h5Q=g z-Y%1f7}}$Jr-sS(EACXG3RdLSF~bTz$r@17i}!)^(`7(_rdMf+tqtV2iku(Cz#4Pn z>)SjGJr7V_PD^#R8OZIq#cBUOulJh=GUmK{;4gf zHM`QKgi2+neB)gj(zB?S;G83QcR@qzFQ*Z-Kv#1?{c5xX50gPg6_mYT&|i-U8D=6t zig5K(>?5IZ)V(k#J?BJ)+;^X}yuLG)NG-Y2P4DMa!)*B;Ux3&*?@*c?tvOw{|CO1| z;5s$6&GhRRce{~V>ns}>d5E`f9Mv&|F_=%IoCKSj#w6(x0VwAmEG|$8&n;?)tlB;; zET)k{>)-hH@PsHmY>!p_#YYH-uqs3@^pyec9MV$LirpkK8$(rKrucw9rUF)Zn;U@Ro{!&~_KkAZa7_PeVP;6eMQchBm zhFLmLjj`lBH+B86etCT*vyy2+HfH3}xv9P9i~ULyAqM%7;-lA%KiHNOKfpX6d{AK30UMg zgHZGj;8=a$i!*V8B<4VKF4Z=Tj3-_4d*rqLC}WS{Z^q4{^}bw9dGZDb36qc z@Y=#-sOq~PJfdz5AtZ7-#FaJn3`owKU-tm&xa%~7k4s6Hpt3fEA#|=|bV@;L845t& zl?wxS*+O$)TGf|=(n3!oIS%N*u=W>VC{cumjp}6ircb*}D60E?xh@MdUky{^E{wlQ z=e#CNo!4}2yy2Cs@QP8@{?(LY-0NVl!ylaIVs|3NnsNCfU*H(@RW|W8)vGp=5ZT!^ zQ{p%+RyYk3#BvkZWfnQ6P`P0N-q|i)A4& zYT}r6yA=*e4&KrbNQpi>KK_Vn7d5IL{BYQ;sV`WNOS=!+ojJ45!al8(`HGjsUqNN^ z#V9&W%)+P*EkBW3qfG&3;>&VBN;YauD&Zu$*B#k%lfGpkIm?feE=qMeONbla>K2>s z9q!huJyYE_RCR8ALy|CE4VRA~;*Z^TyGof`E z5y0Uqq#1Sa=VbYk5&Nbw8j6(T%>3v8(eozy9(|Z+Da1qkjiMkJ4_hX-@C;q3|BN(x zgGL$DIfV;Ni;p9!a2ZtzL{kSU@zid~H0<2<0Gsu`FmeHeM2J7By>tvmED80yVr7o1 z-OyTWI_^flng)A`rFhq|2TbZ^@v)Kd$u(rKB|OQozUa>#E8l~2sE#pxF=EG{Sh1nq zZ4z^$(9X3X*g2&nmlxF~I!R=x$de6`zlf)Ji$x-T(MZS+pfGjq{0@g4~RL&%W zJ2hiWYXYjU(&vuBm$8DZf;27{hN;M;3SqIGxlkZu+P$h$AB@5_H~c+}`KS#8VZNK4 zDKdo`(17^IYPmJjdaeicH8x^9r{t=O-LD7*g>;{yC}WHrA2e>ahO`6>kDm}f#@Bms z9g#-cBOo&m=2zt2jTy{I^OXg?!;DFbk`Dbs0G&LqZ;ojyBTbYtOi4GBFUP_X04As1 z&?R7FuFXE1W>ln?AVAP1p-R5?d$=IDqXeKV~CA&IVFv;9}t`2JTU0~deWhz?mAM~ ztoy~g6NXCa4qHf;mYOHwSY}&+b4ZsY0{PrEN9%4bBeh^4hoqxfYW}luR+$-RAOtI2 z7p+P}`NZ71BVJUM47c2NXoJZ&sx#JP($a7Rt5mbF0X*Okk#HT1qp)rGA!8LI{#@@J zbh^i0z;tWlHH$?LtaO1oA@$XD*j9kK$3wL#Lm3<{3EmG(DMjFMqA(HQ^rrB`r3es?mgxR#1hv2#5TG|;8xWmf3tVzvY_ z@|OF-<1=ziHjHB(1h!wmD)T9`4aS~mEzv?uIWaen@n>Ak{Pm{zizvz6t}pqU>6@A{ zuley_Ce@qJn*|72J$Bus5{$CE}!EUg&csvC18wU z>irV`l9;=oWWM1<@h)=`Mv}&rM6GJRm|GG@5=V?1G~ocOVi%NBK!GcqMga6FWyb=R zH3By&e%vOy8+*USK8tf8X2Jdd={cKyb64`)a@us2dixXe3b!N-fdA=PasTonKfyVSQeUVl zW)95^qA9K=CcW=5-Cea-qE>vG{EqkP4OO@JX1ui{o})(q-5mz%BMkwVqbB4`7l(|% zHZO|J1}+^0n{?xL3}t=a7Q6`4*uky}8BX4sa{a43dz6xor59#SB)3%V@XI{H4RAaA zUJqh~^>biiK{F@j%=*l{~qAXNg`q1_40R zyyC=R*sEk+v35np1pH}t+Y0UB3jggyloS&x9K>%-qFv1i5uxk-uYS zEIfK8Y)U@EDmv?2NVj_@R?%(itz&l%O`US?_7{u5FoQ`JR&Ek-64# zjqIWm%xScodq!IoHdA;X_I;TtgG#0yi4yU|{T|~_VPUYepRKJB*RfF)^J*8?i{GUr zqqyqYXy%pI*jm&3K+B@LU`Y$rSW3w2!t7(X97#9H4GeHus45bz?nUf`4xr)3>-(6i z7u^pXeY#?nAkDKb;O@2QjKu{BtdF+DJf`1HBc@#M`WhdoEc)7Rn7l*#hT$#8Hs0dM zEb2!(zB}C#&&5=wzDZzAwN%8s%-LGLeDN55WK?Y~a3h7?KF10Tv;3Pgd)rIBtnwe0 z33Xr4%%X1XHGMt=hhj$~AGYz=iM9t0&AFB0c`g(~9B`>#R)Cj;j0U=fC+U1d+{f@s ztV6(9kBsBQg**T5kIh@l*eh(Ki&?bT#{+j!Go$5LOJWjIWaVD6si(`A?m{ z=2jyg>vw#_2!ds*#Clk*f_sOD<-QVkOV~TGu?t{{PHnt=tiffQg@-T(NYaf1qf{S|kB)sjAB z#GSQY-y0z(&|`;YztirF*K5Muuhs^*U$2pSqOQ3tev8BFRD9ThPq#tql=Z}1Gwwo2 z4{u5%!sXa7k%dtbt&f{a1)g?%p|#a<<5CrZ4D-zb!6K3##d|Lj}wV#op9io)yd9;s4% z6nzEq&O&hBtfs#)n8yaXnq;W&nS<~rnxme(??b3Mn(YhQOZ5}qK-L>>f}H%|djlqD zv@=Ygy&A#Z)8o8P@P2)A3y$H1u*1Nl`yIpU6I{ER&mPYDiji=EU3*Jl_uJJ|sI0EK z*Iw-9RvSw(xsUlKnEAk1pqR3 zh;%@4J5R4rT&Q94z;PvMY&!Eg|Ma!j4(CDvaS+zz>FdhtA|5V5o~vjGZ7!#N-Bse^ zGeK0mrf|RVZw$IK8x@$c5JW12ubrEjOGm#^41ja_Dv}?Ok%9 zgwiYB6TzBFm+r*gy%;H9g6162VT;{oLk$mE;W%y`uHvMWaJY0iM*e}AKVvIdq_o#7 znQ=t{>*&0V^5FroX?-$FLz47rB#?}=htP+3o$Fu#y4Ts~67!m0(x8Jp9j1=qOGwQ3 z9G8fpmL+xbdrsH} zk_#}+nCY5<|Gi3j<@8iiuo(YVa@rDm~eS(XXRC3>D zc^vWm8$S>%h4H0vpKI?Ms|G2H>MKn9krsznQ)@QAhgx!v?xq)XaGV?UCu~J-c*Vyzl*LY69Km9lvUk%Rfa4#)Ff$zxQ2YaOrTq&T{slP# zX#YY@{|lT10-k?nM5#xQNOS|>drvreFlV~*^BVVw#1X4BTuvoRxIjK;MJoYfvA)P4km1`^I1ij`kLL#iaQ1?^@HupFzb&or&GX=A4#%Shq;R z3oo4tpJGJ9`NYyw(S9u4fNhBif@H~~P2=FTH1lbJ%f9NGnT~{!k8JAKc=yGq z2*DJ~b!1cGYFEMzM>lF5ym3jEZ3JnZ1JFS{voDm@^QYdIPMDyoT$XEKF)5dC)j>9- zSsoP)$RRWm+6nEx-%TF{%M#)h;`Nb|{P0%BFEa8`w3&ea-g-R((KGdCqZ^0#k>#u} zh}+{bA^pewuqph(q?$~;EhOnGhYIT)qpa>3o%`~weB3oPE@98ce`w)X`WX=5rDI|K zQwhIN`zHYYCoug4xAW`U>e?7v*;(2=gW69d32k#ADlWszrzopH{z2c&LEp|;SNqRm zL2YwmGiLzFpCzO}i_EkQZ2=5_R^$DNn^Uqf(f}wKSZDyW3=E6_T3R~VKk##5J8d&# zT^hJ+PZ+lB}J+xdMRY zXM?AU55F#+W&&JU`gA~W`=4O>Q%b=vu>3Dk{Or(P$L={RfKQnJsf(v#y5|yJZCm{( z%KwgNfc0gaZGjvJTNqeAHF%P5!M4v z(b&%U&mviSD=Ra7^QVbF&*~58|7qmEbOjXV^PHagQUs{d0236L=;;CUzz-uM)Bos0 z$ISd(kL6bz8X5rO)90x@@RyE`?zz9G-(P(`wNY39E$!fcBFHBR{?yO^LV^Moz+}iL zp-}%43gw4q$^TpkdfGsR0;OOJoGk5gQpnFd3^_}A3*)CWO2D(9`CC@RQ=-T(MgPi* z_%l=FKT|QEB>F2Wf|-$lJFpc}NXg}Z$`8c%&k4uIlrj$f9#!Etc3r&zgM)4c|r^^o{{Zt-?E zAChOC8a`vH#bxOz4}TF6!~x?Bx!6K^`h^u=;tkPyv7+N*t6;Ob=KFi5^%m1*Gg?(Z8OJRCv&j>x|bo|rSgs5zg?jQjPa$?nqa)gK}^P`={Sv`eNrU#?Ucq z#rH~gOUsAlB$SE^ce_;>6&&=pmez0+E?{@eq|4LyZ1 z7ZSeYa`S-F$*T{%c}8?b$~$tnq6hIX$3-0ZGzI%fWcem07VQMDn|GnaD;srYsqacH zBH(RM_R$=ZT2%_m9Fjh!E`lkCCk?ps>*-CzOP7<1+c8tJvw#YY@d zrWb+JXHX_zFjD=P^D$@eDs(#6R6jwoPV7ZBLZStT z;V@baIFhiP$ZXJ)^I6B6(jZSnHuv~!O$-g?zl!rhzK62dlqpm8^jcF>ltnPIM_{K-aid7jp z7;UU$vETHZ(X;pMlZ{{Yo6Q6}E_F9o#M~!%i|gr4-JIw)9KId4)vNZOG9I7`PQvug z%PmSng2C6F!cXYTYRtE0KzkdEm6DznF!Zoa@#9c*|3q7UU3^%`dF6dmT^E?m9=7#+ z7m|(F$6ocCY>5bnnn)|F!v&9D%@T1<7MVI$nf8u+mMpY+jo}KI^eWj?X(wf_Ebr*< zm_3qi<#7ff$WbN>Oi8;?lK@Z(**8epb_RU#QrUuEGWACCWpcJNo168bN9f|jWn)63 z*TiGZQO}$7d^`7)yfN3 z=)g)li?R>Sf%s#ZfsxEKE3YVMSq(J=ylqi4HiNGy%snntxv^EdIhK`MI0$oU52xdFsbeZL7gxOYq}>6MwNOSq zOZ|9V^`IO0G_eRFX60Z>(qD@c{1MFf8egmyyRLt9jUIs4_Pa=3!+D8t-^YLc zO=hm$+Q%Ca+-dadm`DqLZ4|w3Ty+|B8(|%3O3(!I4BBmZfbW82WyV#3(3R&S)|KoO zwf;F4N?U+U60}8J>Gio1waVbQlVVb&PLY;nrTi}aPCas5c|>?vrysi*n4t(fSa>E= zI|RLq`J5KGEEAhRHUX@(f-fP)1e(bqn<5h`t4XU}i1iGJslW7F)C^`M28k4l^Q6JB z3DobmihN;lCuTpUdv7$?&%($}T`EgFe;f&6FDXx+BDL+@SKf%8)p;WPIB}HFnfjQc z^yBt+q~aD-p=YG~$3@@Qwye^#g74G+cH4QC)&Q~i=OrNfc zN6XO>%#CqGNy6c=IV!xg6r;`SpjSaTT0U-#y3A2_h#Ubgu1fit1zeGpukoFbt~#jK z@S8CL5SN0MyUBwb_w>x|vMJ@5R%u#?u(v!SZ=~0|_){bM{nb=zupBx!a(|G-$wDdx zFmJFe1~s83HKd?zGku%BL0yXkvn^z3r>yu0&mysW-F7N@pIBM>wkg~lS33TNyYb^l zfjSxco{t9EwalD^E1U4Ow@zx9c>N2{z@}41qJ)NYWBt~?#o^`=D!qh?gSnad6C2-q zCXWa~oF4?%bk=Ov;3cd6mR%~dADUdTHN2ruSw9AYtE0A)mGtDveYc{v#k^j#&VMbr zM*ao&p1w5zX9ERh1QiuJrmsZTfi56vCuyaUd`|t4?=WLdvOdYKW4-g*Q!qdHS9^m~Y-b;!6JI}`j<(0diilPqY#HG19`qx>i% zpJqSk<7$METM=^k5EpV3%%79;L39|Q;p^(Fi9t^clo)(Zq!+_E(6gO$KGeA)i9?yI zIz%}{!5qo6eRpY48K(M`Q%s-R(F7MryyftEF&{8+ytZiq+gvq4V6r2&q z^pQ#r8ro0Q&)kl!jVR-RP|Snn3*JFNr=u2{eF~PO>MUC!`{d0X z{_<|22y@W|7+D%`^@tF)(}CHNwk6U+i~ZR5<>SZ3E1FO2REU_3K5m58{h5kLLfQ^} z{uUVagQ1OuT{vYS^AFngp|F$o6cJ36CZ*8q?yZmmA0`(?Op!Lr>N%|T!O=5w~%v)JWoC?G|Gz$XfU7Mv#+oo5RRwg7?d2J6e`5Y#1**l|17;50= z^_vEfHgW3AOU}oRBT&AOCFWhob1n4h{g`o$K3A07n{viKk>YN11_hB*I3FQo&so;r z$Mc}8oBW@ ziTB&WT^?FB%-4kb3)|AM5`O zSy2@&|1yZQpO~jCqOjA=l9RMd(XD6UH0@H(Joi!`)omJ)FYp3$KYL5DT z2Lsz|uhv+9aQR_S@i`sO;B#2gb$%O|cP8tEZpHP(0M`7d`0=zxQLsM9tW|->X$;us z_7J?qt&3lIO6cjO-jO|7C$k$Lr00PNyhRI{AHJfX6YU(PH2mY$%g)dAP%#L{cTkz9 zVSKH6Z#M^ifXQV%-jr!BCMrCyd^F*v;vZZ`5yUN-Kc${B98T?YLJBiS%=6yQFdIVcsJN3zTDQ? z^`1sZ$}XQ1bXa|1t?Waf!VwsJv)su`c;j%UkfQOv*tUW~h1^bGy}hr1QDn{RaltC? z(T2SR_!!S0iJ~?BCTLU5&B{!>T@R-qQ+E8d>?Udq+UxCfFk^I+w=V4h%X^54EJ zgMb*0xlG+*t6I{X*yAkkm7>27+!0@JB_4w#J*Aq5kbTcC8W6FTQ-R5G^YY~>+;^R! z-hcob1?PF=T|i0JyilCHYPYX%k}T(lHrGp*Tt^>20N2s2sj&Bsn{a3c&&andRG}qm z`kUUB)nxX(aM3fnITjK@#V z`})%@*=M22Nv`8ALO0Yq($I|9`h1u_9$icg_D1#aEW~XhfbU}8Joi8JPW6`p7@NGs0}$R>^S5%tloqBn;%C?LyUi|H?)zq;vK)@>m#qiN5+P>CpwZ0huFu_^y9PdlmHRD1NfrdFb) zruGVI&;(iEG_Hnip^RoRAI6c)w1)KQZ9kZ{uQJGWpin#M;fUo-OeiaoI7d0w5b~~W zT+x64v+F^`zXkRYXZNO!?LVTxPhqYjOn)BIqJ26Q_0N#l)0wEBCz$@1vG^!) z%Ps(n|H(ZR&H|r6%yF@dAA@}Y=2e#d-HVs0FrQxV_v8NX@RLj85Cw3~Qh(O`Ci)KNVkyj`yzoR@g$j~oiAQ40|x%tFj!=e_pBhU9$f zf@eoW@Vdqi3AaTW{b2eB-RkGCubhk@BSVMChr%f}&_WgzTE2ylWoW+kP8t3WIsZCO z1`Hh1G5#5Cr~Cap)?Xp!zaqHLfz79I=Fe!mkfI8|Am9HdqU}Eq=K;gkzm5ey9Yg!I z;{FoFvjAlO2hsNb6wUv;(e`I;{xp!^qU}sfKO^N&(eYo=^5?&SQSyIA&3{GE{}w&} zWukvaWd1Ig|B9afr$qiwM9+Z`4A9>H6+QoXUh+4^{ReDE_dkd2=%r}9x{&zFcwd2W z;){xkLdyoC(h;*J3=PDHu1qLqQAxE{@Ya-#tAVcl8ac(4` zp|NUkv3X_l0>rD}ZOp$7@0pW792_(2pNs8xlTt~ws z^UOLG2$en9!G8?>3-)-jdAdJ`{vB-ixz2y-N{;$${UUnA8>P32U?ax0fUn$Z*xklT(Ahxd>L8KQ06m!hGP$^FEIellKq z-|ADs9XS#DK{MFDG>RUu>UxHsFurBTE5H6g1o$eAp>rhk)?A*{s8rsE<8dN64YkDC zh%5D-eAz+H@ZwOV_-v~SKRvXqg#4BCx*1Pufv!GncNp0Y%Dsk7>whThcPCBvCocVO zg*`dsf7>7ZGwG@9C#d=4+kyXn8qgD2T-(mZ`014JlVJgo&gbu+#lMySfbOSXe8QLb zjI?c}Y%Fzw$9R7Sqik&Lo=W}zL&dcJxcmupDjMtA89hN?zyYX$I2;Q--OmC3JQECj z+{hDH1sni~@BWF2(f!k^>7O|K<&OWb5c+3_@N<^`i^r4vW!36noL~%ECQxq~KUY^* z)5v?#SCGF0OtSn_T_(O-_ITXG4x;dP5--V`M+S_ z^uN5u({Xs9eLpJ#xTyf``*~Ks1DEu_FU?;{`hBhZA3sAz{r_k7P5;YvK0D+;?fZ8p z`M-DL@kIZ>#lD}H$Ip9&zwG;2=fCc({~xjM=ausN(s?rFU*`RD(ofGfuzfaeQDeO) z7xKJO{P`3Edtjyk?QfVq{qJJq13r}qm>U8h(bYD#u>{ak(K1jm0Z5GO?5x~!X|yy9EI?NOR*#mM{z>J3r!g@+P3iA6p!)`{*}u_fXjz$nF6QqvI$Gv` z@yE)*{LlVqfLS&`FZlPitSn6buE)aqPZ`hxr~A)(v@}c%{~VW=k&)(~JOgP@e)@0Y z(lW9z{kt9uD{$ZWdp#yv#(&jgV*HmrXqlLSH`xDd3)}(!MPvDwTxgkCY5v6@6D!ld zbi_KWWVWmK8HFiRoYUSeXCC3o{GrKlKJ&lK++q3vf~X zvp*Jknt!erS{9}!U-);PS(tyHubquHFe}RD`Qc=;#xD9#HuH4LAY*B12Y9xO<^TWy literal 0 HcmV?d00001 diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info index ca7a623..a90c94b 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info @@ -1,7 +1,7 @@ dimension 3 -niveau_commentaire 3 +niveau_commentaire 0#3 TYPE_DE_CALCUL ---------------------------------------------- @@ -66,7 +66,12 @@ para_pilotage_equi_global ---------------------------------------------- para_syteme_lineaire ---------------------------------------------- -para_affichage ---------------------------------------------- +para_affichage #------------ +#----------------------------------# +# PARAMETRE | VALEUR # +#----------------------------------# +FREQUENCE_AFFICHAGE_INCREMENT 100000 +FREQUENCE_AFFICHAGE_ITERATION 100000 resultats ---------------------------------------------- COPIE 0 diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 index 969b3c3..63fadb1 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 @@ -1,7 +1,7 @@ #fichier au format maple6 ############################################################################################### -# Visualisation elements finis : Herezh++ V6.687 # -# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# Visualisation elements finis : Herezh++ V6.785 # +# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # # http://www-lg2m.univ-ubs.fr # ############################################################################################### @@ -15,6 +15,20 @@ # >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit # correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions # des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds + # >> le nombre de ref de noeuds, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference de noeuds + # le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti + # >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference d'element et de pti + # le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # # puis pour chaque maillage # >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie , # puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees @@ -47,6 +61,12 @@ #---------------------------------- torseurs de reactions -------------------------------- #0 0 (nombre de torseurs et nombre total de grandeurs associees) # +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref --------------- +#0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# #---------------------------------- grandeurs aux noeuds -------------------------------- #1 2 (nombre de noeuds, nombre total de grandeurs associees) # noeud_8 [2]X [3]Y [4]Z [5]V2 [6]GAMMA2 ; @@ -77,4 +97,4 @@ # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs # etc. pour tous les points d'integration - element - 1.000000000000e-03 1.010000000000e+00 9.970574516706e-01 9.970574516706e-01 -3.777977376754e+00 -4.888446928246e+05 1.000000000000e-03 2.134381140592e-01 7.863544198951e-01 7.863544198951e-01 9.851975296540e-03 -2.955587360158e-03 -2.955587360161e-03 2.068916079252e+03 + 1.000073916208e-03 1.010000739162e+00 9.970574516706e-01 9.970574516706e-01 -5.723375540897e+00 -2.061055598212e+07 1.000073916208e-03 2.134382702626e-01 7.863544198951e-01 7.863544198951e-01 9.852692719186e-03 -2.955587360158e-03 -2.955587360161e-03 2.069118889116e+03 diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/plotPdf.gnu b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/plotPdf.gnu new file mode 100644 index 0000000..a45fbcd --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/plotPdf.gnu @@ -0,0 +1,41 @@ +set term X11 font "Bold,20" +show output +# set key 0.02,1000. + set xtics 1. + set ytics 1. +set size 0.98, 1. +set terminal pdf lw 2 font "Times-Roman,10" #"Helvetica" 14 + set key center top +set ylabel "Acceleration / 10000000 " +set xlabel "temps x 1000000" +set xrange [0.:10] +#set yrange [-6.:6.] +##set y2range [0:9.] +##set y2label "Axial strain [%]" +set grid + set output "acceleration.pdf" +# set title " Strain = 2 x (strain_max in log)-0.0035 " font "Times-Roman,10" + show title + set style data line +plot 'cube_dynamique_explicite_princ.maple' u ($1*1000000.):($6/10000000.) \ +t 'acceleration ' lt 1 lw 2, \ +'cube_dynamique_explicite_princ.maple.ref1' u ($1*1000000.):($6/10000000.) \ + t 'acceleration ref' lt 6 lw 1\ + +#--------------- + set output "vitesse.pdf" +# set title " Strain = 2 x (strain_max in log)-0.0035 " font "Times-Roman,10" + show title +set autoscale +set xrange [0.:10] +set yrange [-14.:8.] + set style data line +plot 'cube_dynamique_explicite_princ.maple' u ($1*1000000.):($5/1.) \ +t 'vitesse ' lt 1 lw 2, \ +'cube_dynamique_explicite_princ.maple.ref1' u ($1*1000000.):($5/1.) \ + t 'vitesse ref' lt 6 lw 1\ + + + +# pause -1 "Hit return to continue" + diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/vitesse.pdf b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/vitesse.pdf new file mode 100644 index 0000000000000000000000000000000000000000..21ae96a0422beb3f6356f80d9b02808908c772c8 GIT binary patch literal 19171 zcmd_Sby%EB6E7H?5Flu97~ExMa1HM65Zv88A;I0<%Mzb+uJ@RrgcBu9s3lM2r#0#D+v!yIXdJ1O~7G>0n-eep)WN}_RdJ1L$)D_%tTj#B<%I?uYQhI%PHYJLH zK3@`xIa{iyZsfw){QmWOI~=&!{pTWCgXmDSGNHNMeAQ*7EFei-w{b^@{i}_W!i}yc z`>&+9)otKR;55Xn)UgBV`g?2mU+XNs^o$*i*UA)N6p*0Beu2+)ymp-QZDXUm^bx^x zz5^D|Yf_Q?kHR&N(Qdsq^^o8NW3Cxk4V3>6$G>Ah8ESY z0ziAuu|8M@zr+}|G!daB&Z82Hhfkm_FjeD1aV1ok437R`BzH+hO{vmT9G{%d9@1Z| zp5in+$S&7H6cjfuwC`{<#VLrX!>bt5Q!I&AC}`4OoT4h{w*V!lD6viM+&<8Ah*bux zs*9qeJH%ysY{*(9gU43gZ(d;pl!`K<)ir{3(!&u~Sp`~Rd@@%^0a=_ulMr|0Vkbk?sn@575e8GWG$9?Jrmd>N$6ZIA4v(L^te6dVHFxS+>k zmH4Sh$w;J$1blfnNNQtBkhd5KBPmHfC@h78XyK9gzl(?Gc}AJ< zYPh8}oFRiFrJ@hCgoYc-0+N3TM;z7fTFw0?hovLkkh|*}jKsh%n{|N4YWgx2{SfcH z`^#0i>le~BemGK=tHN;4GN}~Kk@3O&L&4|BYC-23Ciy+aq==Zx6>rars5LKWD4fD} zv|4xl4hWd&D(u;Oe62>4S1iT+ttt#xs^mma%dEz}DdFg-=&wL}FqH-DoCrsp((mf{ z#=zWICaWgrIwlH7%)})I$?BrC%egOhl!~J>wfIWS4U+Y&j7q^C1t0ukDA*oFjp2ph zVa~RekcrGcYwS#LA9G=po|2O&M8v_?Xh?nD&SpMn&qm(Crcr4xdS0F&7v$H#RCe8L z|6#2^;YA9aRL_|x%Qv4nj+82#DM`Y`(|6D=Q|05yc9*KTHY4TLhKlmK6sRXwUKo{O zW4%wuIK|9)(XtMy_}BUlzI;CHI9`QUznT2PgWewT{WV!uW0U~`1q@y+e1u32ARlr*-NWku#tzG@9w9VY6iONstyIR?S2lI5{X zm#>d*vRj5@B?UU7CD%W_ci;vGH4paMHcy4F9dLvowW%E_`D%)z=D`VjWl)0I6j%&D z7l&*(=h)d{kb03NHI_R(N+TMl^4~Y9Wojb=z zrPVBH@9Z%4gja25gnpyHJc}y7oN?JJ-gd(3hs(z#X3LtXcUE5Fg9-3}H}(9nT~U6n z;h?gdu0(Bw#!`OsrCc`!Veg8d>}j=C4p6#GDUhOhi+$m1OE-VV$>W;>cvs?7rJggA%qAj~fZI$G4Az7@nH;*X+@HQl#D04;gME;^fmJihOAoihmBUz@>Ks z7kwR|N?KT4t!wVuT_g@g2a+sDbFBod#5(RJVC-^2ZEY?y2ITgPrHWix!(vwse#0Qp z&g#(CA`g0hR0XejnT?#X+#I1d?fVfN6J7ODT>Ck4%6P#?T)5$VIp>V+l(4l{46JV=n}-jlXr zom%v%hVTyH58IHS9*g)D{s=WIA&gcf=9BBP!R{|R#@J2*&&x6mU(8bJi*G1#@s&K1 zIhCa%D^~IeU;3deSdib>Ize8DPu04bO5gHv2lZ9>$)Tl&nPHTHH)DLQnD(Orlz>R~ zy30gL>&Mvn&}@c>Pu*|Hp24XXY$TJ2?InhNeUY_)gbxi@C?NY(R5%@LA9?6Y>RUQ? zztT2_cY$-&Zg1bF!*hgs!Amd0Tn3gIE7})XFHLGwjUXfwTjM`AO~`-u+wJ`xd;bgBWr3Vv z)_aEgBE34 zHvup!7?_zn^jN{b(FD=~;IDGbvL?nB210i30Bs0~1psDe0f2z)0G$VJVMu2o>`s9D zQ^?OFhq$Aii~Vl{_*wr0^_K*c0n93n2DVQ2_iY(@{31&NoE=>rdT~EWVF;IqiK~T? ziITX`JaRS`STIsh@0q^_wYgm6Ag@FGV$iJP1?i>2a2?TI)0@?p^ zIBL=Ia#xwC)$_`9T2D{uYX+x?n}Wb8gY^3dm>9vvh#nZeF=oCZm;kirus{*NH-a8w zc{Ffj3SHugg2ouWGFe(B3Of^Z{3R1ocD24qo_60`T6znGO!v6WgdUZj96xP)8`p9+ z8&!Qe+txBUYc+xccmxPy(tplVpSy5zu++_ms_%~Z=7Z5~6`V9vMbjc23pNzZ&Xc{( z^w)Qy+I-%XI6YtMYU>8)4li0G9+kmc(_44Ki_&FAoY-i{TWNpxw(BIFxzX6~@$4j$ z-l%)j?Dh=62DM$CR5RZsv#Dfi(1Nl(Y%-vobhFydb6nB-c;k*wfd)o@X;oKj%RL#* zrx*_%&5`V*Aw~`Qn|hWIm%Nl}38OxvK7Cq&Cq0MsZ}9MXVlabDY;}Ex^N$VajOxY8 zt`od2(jKGN?H!>Lit2M;xM@kdt-KhHhIYg?{PwZK5najklS}h#dJOO){S@3=OfWdW z)(~E{^;Jnl&<#am=`D2DzC)Dnqdi?85eY#WACn8%#HBu;9I3KhP80rljBlEJ#!7ox7q^aMi-ajwNU+{-4!u4xa@vGbbWvV45-D(T!o_U}m)9 z-i30QnpM23EoAffAVN*etsafFi4<;2bmPG);muJKV7u*crsJ4mOWa7^A2Zi6I7wfB z$8d=yu{{{}s#oj7LNr35P-Loqv1X~y2Tgit`YL0dmOQmyg=8>>>fn?|h+qO(Xg*s+ zN5QL2Dj|tbX6;=6ZGZE)dRbU>2x&wPjTfy8Y`PU531d-z5(RNHXVHx#niCBX1JM_p zjTeNMGCV~@B9)?B$Bxv_fk;k{jJnZIx=qv*^f28|8OXoVU4L44J#YUGym~C{F3{`5 zBJDtV719-kXS^)TVj(SJxjS1sDw3*!&maM!#Qur{FYd8W4B}4~LOK%fv1$%?i^A90 z?3V{ok!1*lhb&V`ZoH0JF6XT~3CNZI&Ig~yQS5Sk^mg(ViFpdAF%DPYa=&D-ziBV( zc8#6{*mA3|rNei2TjLi^jr*Npwyxo3u`1yzvE_F5V2kE~(^Jt?Gem=Xn8DQ4(=!Ef=0<^as8yYF9T6_s7f zh>MP@dXlTKKfT@>AMVdh+rS#C27RaZfy&XVt6}GKIsCq>Z8?>>=}zytwS}XmRx3w0 znD$}wjqwxGFB}9@F!&YiI;7)bkt#^tS%KAX9to|%nwCsT1u8RZ(uAS!B5I-)qu(E4 zL`~42*)LMm&O8mtI*T)gqAf+HRrVc|N5%q&b*v*0h0Um9hX^!CuuyK*@Vaol(p@MT z^Gs@0;0SG0E@nb^34B~?`St_V<&?SeP{(IfjGe{T``OfW8!BzRC+YhUt$2J_1Di_% z{DR%OyfSocTzA+hjxpg8tvENw;@eA)ENs{+pK1`zDuK?I^=alhOzdy-EHz=Kx^3#L zt4=k(B|z^VUHBpe%wT_Uv}U;o*X$a{;HP8;3Kay!(;-rn%hD^sm0VA`dK9_4?8bYk zuP38NG+?9X)@CT8mUoq~qKWNiFRI9Wgv9|Z!jGl9wG@RLeD^~Q+Y~hC%J#UbT_iZ9KGt* zZ|lwCxt<{~cD+cwr(JJ0jta68cjzaP9NKm9Z#~366c?B~uHT%p6Vc;s%5*0|6=*Fy zE7bOm6o-oMx()U=W2R|`y*LmnvV>2e&CIu89V9P3>Pig=^x{3j-+h&cdgeH?(5mY3 z?ubK^otxv5-NdecyRxXn|7Ck_!P}Q63WWn-eQDEalZ$gi%MfSmWpuDKQpvF>$zQQn z(wM2Jw!2U{4dI|XdrA-vYr?`exBi71Sd#Z%M1aD$=td+P=NeFL)XiO#iCZn}t&nfj zJ;w<<)2Fi$P-Ymy?Mh^Lxov*n;g^f_BrxNtkgbf!gmh=!*pCzxou~YxtxmTEnIsx* zr&~l`x#Ht?tp|(y{LE*Q9M9&7TMk@*Fuzj+WuH`Ub66)&G`=+7=|*{KGk5#3s&K0{ z{Icv-=BxAGz+x^2_qCXrcxDaiX8+Sx`b#+0S)_Ek|-1St4Wu*#Ix-8DnRhznH<4 z5!GApXVcj8Y)2dUxNzI@r*5&JO`Io3nil)d3lm=U*7BK_=}6ymzgg^S)$ zkX~lCg}fzdOO?^OI!L|Q4a#&+_PXA78Fw-hy?fT;Xq9+Z?NlzhMCoBnsz^+I+XxW|{a3Y5^?VYuvPvhN5{5Ucf$^Ls|Ef^Uld40jqR5qlXrNaiQh2TH6fC@iYb81%rb%A7FiGVZR~X)8no-)>ICKo zPGWg|PN+^4j}N9IpY2@HxK2Q{n;M9m>FFpjv1veZ)E|(&XAoBUS_XYrP@9a5=53#B zAi=OvCa#4J%veI|I>G01RN1-uD-HGW%k&oM&46rb+y0RfJ@?d2W^sh4;&-QCHQTgfu%D(x5TNIjO}**5l+?}QSKC9(O-$bGNjkl3HBaALHHNhdill{Z`W;EBr)>kL?{|}T zqu3|o@(!mt!^nK544U=N9Mtr**!Qb7np*`LPB~D@Pv?zy^BVOO)lPqdKnm9+VM4of1WZj+Mz zB1^{J3rnqOE+^}K-&$4L4jo=pqwO&Zd;m#;6wXJZv>%!nG_%o;Cm5fzCX zhG^ps#8kNpTR7bvXv$vX7hov1z@raB`#D>}2-XkB-&L(yJii zZ0f|XQ~bAA*;Ar-bLbhvo(eb%Y>I8mHCtamO)*owUtH5bsaEE)D{BFd7|#GYZQlSmq{95s7-&&J03RC{VNQz$d9#Ile-tHr^M8N}dTx}NS%GziYE()F>abzA(>y_^` z45KQkb=KvN&yjQz%1UbJmcnaNS}!b~G-o`{D3NmsYT_$h7~fmSJbF>!Bj~bTHEj)} z5z$ALHj|u?ygNmMgQP2;9ptH`ip1Wje_K0byR7MTx`%9n%B*fh`;@DU@LgxR@o7r; zMN9LD-HpL#fn9LJG8Q^p_>xXEL5LZfwUA3D11#qlk?{@}sKaj|#X zw|=-wf^uzi28Z{2&zJ|Ki2N{6}l!S2yqFdFs^1s>0}m#jf=ccvNobG+Q>HS1u#H)bEndY`^h2o*FfM1YUb~Lhd-oJ`Ce;OsrCi6{J~>FPwfn$ z3Vew=_AHpNMKn_EO^vDW^PD~{hqsNNm4ypM$V2an03OO?VN@yMGZcHA93-s#)yd$w z+P+!YsGN6%H1onK>bRcn7Yb1{VJId*hy7)nJf`^eITIuvVb+i?NxBqWgZA1I)y6lH zBin42L8&F>6k3Cc^bVrucZb5fFKLLbVVHE8x$DunK1p;i!ar@#-TUmVliFvhDCcl` za<)6&o)D2Ly9JJydF`;B;4!#ZIs#=+Ef_X^1%$6zwdH(dNMd*yry*9@!@oBChP{#E zRA5l!lX5s7f_Y61tXxND1p=xc@>}8Ar}ctQuLBu;p+E$eE3aaG-&Q`+2>wXn9{0)6 zuJffkd{2Gny9yW^dZ;!{XtA9IDc;Hiji;o12$NlJI}i&H**^{_gukfo(BHY3M&HSR z9$$Gx_-TI$U^_d2Z^+y|?DMD|Fyb}svF!?H<=*`GzMJ=ISCsq>n4G7vbY5pT=P z3nkxCzTDOelP56wY&*CG`wF&pB>W15y)==R`Lcu1l;hKBH-D5U`eT8-wf936L~q4z zgAVy5->5`)vQf-5j17+!`&%W5p`qfl9{p)CHs%X!u)a&>i!2V}v?hN)(ug#VpV2A@HZOeGk zWZ(3vf$>F@4BidMmCTjul-M$s$Tr|x1-$whT{9fc3W2X8O(ptMJn8;kn{(F~KIb>Q zlJMqyT#}192)%l}gTIwoTxPu%^AwT{zwc4h>qngHEBG9^D3~IwBl>_cfaVRX8=-yAi*6@g--4 z154uR`72%dRGx`%&Uw=|eQIHvyqn0lsYMv^NFP4trkT+A3?jLiRD*YHP^)#_TwD@0 z-0Jm(2V78w45-;v7dfeRtP=yd>oAYX0;*vyJqr12GHF5myw9O+q8d zVWnamaH(qE+k%ZqBa0BV7rtQ|gf+eRZp4etc-(PP&#Dqn-Si@2I#5;_-#+sFU_eg# zfkFrGS#uB>Y4SLO5@R(%hm5pQGJ9CFGoC?Q6R)+F$KuOU&xDd`M;)_dPdaU_CVF(= zppvx-)5^&zIICLajNZgDJ8Q(alTM*v5EETg*xdIrw6s(lu1-g{l%`TIMX>J^W&F$c zq*eh<`OQG&6O>np(V`{BCaB3%;VvuR5vsb)B5(#OnT}K4FG!wt&n<_ALdipa>jJJ?LSn)uE$Y9y+E* zj8o>J=#e>Zh&?s>VlsFwt*O??9-pw>g{P?vc3$Xg`Srr|XV5&^v3VP83>Jr{h6s z;WoK~Ey*J;minc@g_N8UMuwh?%GiVulr5&bDsVd!NQ{hvJZ zFnF|@2%ON`i&26xWjgVMBl^gVt&G-izs8Xt8Bp>sMD!Yru>^@?AQ-dw`#FNG98wl9 z1GVuxbiVbiD0qTGn- z?fc>l2cG`oxitY89ijpl9bE#jD?kCnmCB_E`^WvfY9i+auZcR=Y@av$xI#zOi8^qP zXQv>2aY|0gyAlL=aqEXxc5drWuegqczET3=oZtYwIKDxYZzf*Rig(TNKgc6_ywXVG zd+?RS|6mAGDPRA~>fi&$J8u8N`8Nk&_Ql<2PY~tBEea(ON6ur=-d+f!*`5m`a(Hue zR4c5Dp|5C*qoGm;t$6x3wZhR{Z$^)_q8^>w1%9qJn(L9aor-_?m|L)9N%su1_~ltU z*ha5b_mGwZy*2R`#t-?m&j$uHlt&b3J0vRnjtS$vai6xb!-?N)?JU?Mjnl|@Bnx#G z!z<$zXvw&c$%eK2o7C<+df!>)2e*Tu5E_#1q$BTo2OQH! z!W2GlxKH31RGzCMO$y=nNM5Z)>doSC>b5(xEV=;D(Azu~y@QOrr2KdGoOtPOapyPO zDRkdgtQ9hICBnLLuYZ%giE|a(Z08pTo{k{W5(S7diH~{@$E<}R!ftDxN=E)@UgG|| z^O=PjbmH%fn)a&62-!7%n*OB`O9lX;mso<@>yj6fs3+#rH{#;QSi9lTSCJ}L$GqrO zAQg?>CA^|j59m0;x6kuOf_{WgUASDAUc2*Wj(n6(N^(2wdzq7kQ{$^8frjM^I6yxV z(km8r^(FLI^8dDsE=XmMeniGzMcqqUd^ig9!7h3vWJq!9&{p&<1l_Shk6uZlnIeux z)1QzV!;Txvv4X5iib#uYoz<2&9ubSmkm1T+1b#b?TWTo((xGv+|LAoh+fuaI^>L5P z&Ps-P_X^CZ7_v!B@C<=CciTqi`FJaI^(1JkrV7?1H}D8 zcmTlvfK>m61wrtoKkvop#17d80MG)?c)J+va5#c*$q?og-;^T^B6p=3Q(XySOQAJr$=Tz|me{J=d0B2U1LPaBd1eMMi+08pb zH1}sD1R3AinR+OioL@2#XAd&e=_%!Q_elo`HWF*9=BWysd@jM`zA4w;s~Rq{!*gt! zo^{D6rf`;<5<`FvcCr0GMi1p_d3XP`n_aD^5fFoU%M|PUHSl(*`)nj~yBkxh$UwGtyMDepi7enU<;8j>s`EO8VJ9#hZCd+ zlD4LH_XY0dtYl(l0a?2cX$eA5Xp;v?-9wx5j>aaA7Pe*p8cAcw5G|ZN{>)N#vA4H2 zvAG}k!>Im1SMOW?&s;%-`6sTc4$x+SV8HBPFaQksu(7fKS04}u#{(bduQDtw0Jc9q zAQ0%Gu6rIGo!^k$|3+~468kwH{}YLc+Cr|q-3x*Fmk<~w9why9;*AX;69^f#6J&6} z2b}#Uw5?*NYHM+SVGi(MOa3yk_t^WdiTwpM{|UbT>!rKjfo6_-JN3J@Y0>-%G;i~+ zr+ZLvf$yh4s_ z>Xh`&qpSdaeL=0sreUwUB~yOrEM|Cp^->L|@$XL?7RsUWZV{hvv|3DDVh14!G$$@9 zoRhe^xiBo9tQeg^PkYN%9ob@R79a|)GNWZqw0JJ;=9AB3K5`{E#C~V1ISFMioP18T z`cB<#zE!MQNALMH`)Z5Tq-*KL1>X`a>1|=n?ZC(9jE~Nc?ie(wKknexeC>M95DF!( z8gKZhTtU@6XQn=+KsW1czU~!!dRD!8doyNUNMUP3-j_{qV9})GoDt4lhrrdv4(`-> z7ut{^DdwBDz)8dB2%hfJww52&cEnttG^Neu@)U(x7;($ip$d+gN2AwjSvms_8c;Nc zttPjisbW(+`1Nxhd-9D~x}*n+^Z7Bbr1;`?2EYHNA(lH)$?Je%tSl8e@qL-3UEJAK zxSxNiF2Ams_p0MbU3Vq$s-v;&;OfMH;ETpPs}h2qlAe-?1U7S4RkhCoMkZq|0UiW- z7hzht^R#C11R?X<<^(tu6#C+19&1kf2_^5~N+`RgC*f_1dVVPVsOF!H+Bp@7vh_}? zY_7H1`I^fNe997*Qfjnhl=O<9z~~SjJEG)M(+>T-ko_CXxA9?bYu1S8pbn!WGN|#- zQZ*;vJMj2p)obo!rOUm;0+A9JhNIRq8P_(=1F-^MT5*iFKGQWP)S-E$wZN&R$z@r> z;G3E{h93|uCLL;4o>b_ggN)122*@jt7G%X6r)ZZWeER-r%!eysR8f= z+B|w{{~3uR@0P)Z3RMJ$wvw#Hvuf+;A;Uyaq)3G#i-se8=Zr&GStM zlk5$(vOcMyT}dJN%7yQqKfxDV%Y+t#Ch*J7dL1hE-1OX9^n+yMIq>e|-c>n0I|*hP zqlK=(yNj;#AoHsyoDh*%B! zE*5oN&=wVUMT-m9(2sU}hKn$XPD-X^vWU%JME6mIA~pBtN`-!)??OCrXCtZJMz?hT z76sNs1RJ>#dB0>95^qg@l9nk38J=3$Y7qXDQm^TwxHOM6ZAU3>V#4*%$fSLf=?#tb z&Dlymsap%V=Aj`{5@|oWqCkN@>}w?3uZrEqU$x?UQQ#oAXc^Mf`1y#Gr;U!my+>` zqqcW#eYU~KwCxmD;)3<)m&LXPU4(-3(tet4uQGfaQskGUneRR@g#%--^akb z5X(=P$on3S>d5tsjZ0s$Hm{O6(SX@>dN(qxHLV-GsGk!P-^PTSVNw4)(sN2G)=8CxM zBp+F0p+sqKvx5e1 zRvGSgq`&VQsIJNmi+e0RZ>sNvaoi$ycI;I}h((-WO?T~dcspUL- z`!f61uJ3}afsh4iA-i!UPdf06(y84Q=!(NTEZY%t4sq(CWmJ*&CdICcb-Agde~7P}%FZ-aZ(WbFIe%<0eotKZ ztfsQP^Y1jN4@_Z<>LhI`NvD5s2QAybAr;@YH^^q}L3O1L$S^fmn&jfJsGq1YApaI< zZ23{TPtd`!O>qJ4;uRf6P#P*T=a*Myw0C%Jhhki2x}rW#xS3x=PZ?~ye`pI3Yk%Th z)NJ9U%&ps*OfX`OtIV2Tdb6eY3LsyHWZt#VOV-&S!=S^X(OAwBj~-@S4wIrd^e!n7 z-CDR2YWbblYJgX4ADV&xkMtE1SoG^2k)va!nU4;>KA75mx4y}VEtj!4+IzMJem*IPY|pttg@eE<4x zSDbi3^UDEN+s1X}x%6TEYxKK)H#x(JfLx8ci;Izp3uv{jkxzF&dWPDvOAia5lfAFL z;J3FIU51iqmONq1kJ8I@&lEsWdg-$IdYUfv<%?pO<^p`V_L~ZvNMCqVvJ_Z1e0HH=dCp15}0}LvRwjk^S%+xm2w?R7g;f5>(YJLG?^+o^(ghS9+Mg{Gm zQDcnC$EVAhMHHoT2CNP&K#|-bdp<(n&L6jSqa%cCpE2brhsS1^$THH*zn+D_h2-k{ zzHBt&h;hc`4T3r?72y!|!qz+^cgH$yXI>$1CJ4e@2wD6@AL6!UY~!56sKUO?(mFu= z%{%&BakWzW}m1r2JBG!+L z6#;0Ra*Jne2lCg+m6aGxk*`P<6VC-214asU=yp3)@X%vA5-?BEcN;YN_AODCw|mW5kl(vd#*gjl}Af$@s|=e{FvoXYMwt%o|D>Z_Yuo?ZOT> zeZ0U*O!!UKQv-w_ac+z~jK*vDy>c~`1T$PRK#xwRgI7nBK?X3^CP^P#KaoreAHaL` z(|w*mkaw(IUUwVcnn&sSD!L{*+&9#R3J0TkOUqwQJ-+M6y~Ez0vhEzGCtHDMRQ*81W7h+FoZvu^vt9JQSNr&r{Qn?+I_B|l*3GK6YI#F^c(&AAG443k-Go$ltM36e;yeNMZp z*cOoqn>1Ej6X+8notqQb_E#v9bvi5Db?6*zUDB2*Lot&yw!|Dlm6g=#tn~wWg~794 z&5EUMlS>Or<8sSFPTN_+u3s*AI$~&8Yf)!Sn)$N7FY%aBL%}|2z$PSG# zeVoyDG)Z`eYKHWcw{s@V^``L>I?=h|9G8@3OUIO*?PCy!=u#@;FHz`i*4X7K)1)0k zXk;xdNQv^uJ}+f{y>8tszV$Y1!!VNM4?Ixy;bjW;DS3kyVW1bs{bl2fOLg#-H%SV% zB1?uC)2qH*Lxs~5#3Xq!cWzDXGO`zkD5Nhu*bqk===|T6h1;{uyvPT<9BO;Il9ELF z;`jt{AP)Vl=wv@#N9__?M|C;;%7H5B3pW%0?Y?zxxD~EDet!ZcVfVc$a$Y-7$8q1R z2VKXLR#JZg2xowLz(K9LSVM?)tD2zu7^Kei0W{kSPyh0?@PiZMT^G6zj+erSueUAH zY+LC3$<@W&fir}PXt$@Z9aUhYIP|?Mq^yGo;Z|deuYGqgDw((EWqR|;YB!H@MHAX8 z4+SP}YrW*$UQfCJ^Oi_&m2c<~BHIK@si@z==s(9Nnb>d0vsjh3 z`N#2!1TI+Z5L;8acByNsoPM>N>1EyI7<>x%oK~XEY6@=14b?(KZ+l4sTYcUERvlwPr@ftGQwUm$3~B^ysgewB3o$%h6v7M2kwUy0eh z(fRaHA6)O(MWjs>Xq`M?j{$|OkWuR6mzs8;4W7z|3H#YtoZ+o?I> zs?9MdfvT**+=^HMxVO#~-ML&LLD_-&T)Y&PBK(2!G0{9B;wmdKW+HONB4|UEpLjY$ zlTwY`-8GbG$lgvC#Zv`}2gc@*R|P54$Pq>3=9|Rs4Q2Iy9=?dL9br*W zS;c={o_abk+x+@wdX$JLizS@cQGAn5N~D7*1@}4?CvK#1u4t=dXdrE=*}KS1b_D6f za<_atXO_67Jc-}slDlQmh|j9Kjlg#|W1umiHJ84zdkH)8^t^M?sD-|arCA7X?<6Ry zJ+yIKE2=cQsl%_uvHVgM_(ij+ZnvtoPVRke-6>MPCANuGLM``P8OwYDvKyUME$#i= zr9b1lRzLr?THP0KHzE&eN@c0!8OE^>VbAJ6tD8~~cHT$@wxHaSY@WAq|BKi6-p%!g zznztX9sH-;_dd|%?=G?XFq8i-4pk>c)~*wP9JqgtL^3BV8n^dO$)D9FiSRVrQuWOhL=4 zlwN9)x}!pXecKhJRoEPIOxkx#*J>K2V9GR^7(dSGSx+ToT(@t*ZO0TwkYCnPu_u+t0IQ@*aD)%2j~YP8Flf( z=G0n_!{{>g9{Z#X|BIY|h3MRS*V#CJi5c{Jh|7Om%K!1$J~%e--I+hV>*DHKB4YpF z@UH)icY?T0e??l|NA~<$asLs;g8-EO3-9{Bist{Ucl}{De=5jt-gS0%h-dxYLk{t< z|MHMO{15Sx-=}|h-hX-4@6&(rp8rzOzalVy70kbR&;L~-{|nx82;>6M_kVcLe}0hwNz0N;`9NzZ>^gGz_bDjU1 zlbOCAOQ3e-T7^5{>6%e{uiux zbi3_?73U3zTF55r;}S}pl}6nyF!6fiD;xUdkz+5Zoq$k@%q*1g@m=PLk;hh9hTG$> z)Cl9s(A&JAAp-+pmeI}TB7)~;pWo5fr#NGQI&ZkCQ@QQ@%(bKSG}AH5?LB)~M$;ZW zx@ec!xkfC^lx`C^kq(zn?WHdAS9)W7+>jRO#)sYuoyjxUD80>P^g42%4TCkm{7ixh zFr-LexGRofqe^RDs_M&oJD!?>``*KxKmD0%*>>*m{6MAbbgQQb_=%I8>Z#(YwP1Ro zkqPiq1l%Q~$sZr2D*|Am+Vohy49ipnI~cfwQB< zectVR`(tA;RVcjtY);Mv%zJ-@zzHC+GW|Kfq5JgFli#flYM_ zV`uYw=nK*S6U3r$fB5qS`EJM;V*OihlYS348qS* z{!bo6@gJ?$f#HK;)we{#V5_RBsbNul1FoQlO^Hh%`1Z9ge}C(pf|5#Ve?Cc>{-^wX zD6(|ubZq1gKXkBiP_U5=Kk({c;o>&D;QIjm;n9^SzSKYG?4eJ8=r`w|`VIaUn*;$Y z|J&osf2ZFM75>rmfAkytyJ`77 z=U$b6srS!8--nDlJt()dh4H-!dDtla3=DUH+$;e8hUtTU7n?BT`7@BaAOIR80}Dqx z0FVjD%ES(!F?V*h=Vpck#Y2c@jt2JT7Di4?c8+GWKkt4(jGM8G5#)I`|9!Fl>}}v* z#`-_UB4Go_-4HvohcP^i=8v%V{}gg}G%-bD0kEji(M0lDt&_g`sTTx@^m1G00nKqlZXd_XoXmVfYZaBxAa z)?fHoK)`>t3j*J7>wo3r~DPRtbgYNvHCJc$okKAA+&$A3uL>$cJjBfkcs%a{D5ri_h-1jw!zNA{*QLq?+?>|=L50+lMf8J zQ1o{`F#F%-0%T{s-&_9HE-S}B%R+Y9zw>dhu>NBl9H4*7iUV}7n15}91AK3%{!U~2 z`}_hqSpO*(4%UB;g9CCg>hJY&0$KhcQ%=zD<8^j4fZT#|e0U6%vW2I~z0TYp8ID literal 0 HcmV?d00001