From 2ce0ee80433d1861ead6f7b36eedc684ba46e290 Mon Sep 17 00:00:00 2001 From: yize Date: Tue, 3 Jun 2025 11:00:42 +0800 Subject: [PATCH] 1 --- .../images/Multibody/坐标系示意图.png | Bin 0 -> 89845 bytes 多体动力学模块理论手册/多体动力学.md | 560 ++++++++++++++++++ 2 files changed, 560 insertions(+) create mode 100644 多体动力学模块理论手册/images/Multibody/坐标系示意图.png create mode 100644 多体动力学模块理论手册/多体动力学.md diff --git a/多体动力学模块理论手册/images/Multibody/坐标系示意图.png b/多体动力学模块理论手册/images/Multibody/坐标系示意图.png new file mode 100644 index 0000000000000000000000000000000000000000..a5a1c5b0f7cc6adc60a8a011560b61c15d82c722 GIT binary patch literal 89845 zcmdSBbzGBc{5KAW3W}g|RHPLNr39qehEggZC7qJO7~Lrc1(6cz8VJ%Q9g+g0H*%v> z7>te)gN&bt~Deiq^8y~GOEvHWLQ<5>@BTrEl5bXqFg^VDYmL!?PwZ!{+R58 z*!w?vS?em)O`GSDZmI)!UUHA{f@}(^R@lPdCkJkvH=L~CE0!2;Z#k8I>~Qs(!Znv` zq3P_?;E{?C2Q}_2jXs-V)URoUVNO62&z@3RT7mnOU%8H0STlmBNqmyj%+)kg?_6#K zgI|{1+48zEUp}oB$*KL~L6@Jm&*ALR%63T;yF{PMKsxSm9$l{+$sMJ_V#%;w4#yH@ z#cP{R@JEO5%cp0{>#X>$aW>xD5GhM~=A*(moV=$*i*^RXP}Dx&0uXkqynVF`4%u-pIUk}~|h zmM7&RIrrCnG7^$tYZCIGZC(&R&wf4<-)DXP{yg_7h~y&iuj|COM>^TB*5t+M=YCz2 zRS=(%Jl2qT_Kf(cVdi9EVdwnX-lfTx&y9G4(&4GDGYQE}_OmzXGxd90#PLzqFLhjW zl$9jR>>+%w%yLSUn)NcFqzWfV+RSkRV>4-R8f``d1Sd zFyO9^vMQ^Ly^{s22;W1#hj)ROSXo)6oXjmH)Me#=b|?M?xcl10#X*9f-`(Av&s~Vm z-pPtzKwMm$|Dhnipdc@?1+TNGoy#i^UOQ*DzX$m>j;w{VnUl4Hi?zKS>)E)kOzmA= z0C(@6P4xTcZ#peJtpAzG&iUuEhzsODd%`ck_mKbh*u<_$Y|gtDgNI=176yDB(Fx4Bcrvf%1h!WQOnMLNcD)nxc>e;yRKll z`gYcZghZO;ne5}29;9ojSKOxZ%ig9&En1p;-cD4JcY}6|FAm#@QapS5M7GZe|cQE4*>f!rOO-X2GI%_YdAONhvq0{`Q8~soH241gK z{V$5eP85J~qCo}+u;%-OHX*h}kt#Ly(SNc0^_i@~)8%Q2x;h}=M)=R^_c(}+8bX+J zHz%s?(~u{y5zb4u{!{vPPl%eqC-o6Gi}61_o|!m4K7N;3y2@cf_8+gU=fuWFo(7ob zECwx(r-(n}>4B2q|7!VAhS<`QBz$u1qBdJ(gs;$lEoJm zZk)uWd}EtN5EXpA1DW)njo@dE(2*@1G<$Qgr?3%@m<`+iB#Vw{MfA4Padit9wG~5j zr~k`ZW*?BB%wis}X>Mt-7(8yXH@qQTJN};yoz5CM$?RjcN{pMP37e;c&CEfXBBQax z|1`td@w4{h3kTUW#aBPD7#aFVuQ&cnA(|b?_j>UN&}qHzDZy8aW&PvX!^v3U$_;(Y zU1{aG{2W}jNh^>4Et=MATWM%IyQgc!#?66eO4ojIf+$kp8@HK%5pCje1}`U9#~u?c zmVtIFuppGi&d1?hlRI`7jR;k+IQSooHDM&f&32>L9bxTIA`Mc3f6#>U;u8T`B&(Jeh0#yCRH4UH=@^0YBqL3Xzg2lQ*tx*uVk?s$7|WBO=|0 zH_&Q)!Jjf$!czrEX+21Y(ynV^U6udfaB`*p%7u|T|F#w?8EM1>A;Tl#RVf>BN}4-F z9TL6$dP<47N{&Fx!++aL5No;>gs3{EPuJf5;={NIn8Onua=d>wu!D-&blz?J+n ztrZ=8!wf60mLZ-09{DO65!8~f)b65onwEycso6{q78clLN+-dKKfIqgP zHT!@~Bl`=)n)}Z7FXaoPXLNr2A9SWaWrcc%{wZp=H8WbbJuC5uDB>8RAe42scwwLH z_2hk$S2`kD><=#lTg$SnBs>3$AEdOb0W(ehO(|x&O+S3#K(KS;?JN;Ri_!k7q&}_M zqzHD1%_bU=rdX$kJX=7}85R!p^Wc{(+k&5L&rQ4sLLF~8*p`OGi`bVHbfwIFgCzaN z=;kSkE!!$91ARoxmH_l}=8QgMcs{bqrc80%R7nx{_mJ>2%%sI5`(7cg*sR{0ZM}u@ z#KB_1t*uX9R^5{+{)?Z#z$cFOz17GfJz-yLp*|Ze=iwAYbMLC}EBxR|FI^nY#*?^o zD_W0}KyJGX55OTFkH1p#U?K9(POOC#$rF0~+?( zqklulEDdQ&!=vq!LRVf3LLIZYqJ4Mu@|cV~A~Z9#oXx+j$&#uJm>J&6Yz9sxs^-F( zy`B;pHTRqUWO#`b8Sl5G0T<66Of!lHRYVJ@x#>zCONjmpVh?3lp`~$7wd0A33ry=R zn%v&upDQD~4I#b5?(szBDOR4{rYU~j+5vL&iwzGUz&xOhEz2_Zzuw=Tl{d19H1LJ| z1D`6a&DEY|)s8mXRy*2$-oaLKZIifN($P7P%J-Z~m2hi>NCda*_a^p`9c|_z_U$vC z&l(VwGj7{SWGbZW#jjBre3p@|Z4;2W@aGy+?F6a=JZT+RF4O;Za& z@8RvDeB1IbDv48)wycb|U?_SGMNnB!&EWL3_-WPnv|Hi7;4p>+{`{%f@w|p2PBfSG z;7c{c!$@|r$Vz@a(60!s-5C)hrZ7uuYOh?BqKDOm6W3Eqgzg7n01-Od2RiAB$pWIE zGjEvz--aYJ<0(Ja8^D14;ep2H+m$dGODcg?1-xymaL-z`74!P0bx+-IXoV`XLc{f5 zGkWsbY@?}eLs2H)T?mr3>Mw}f+eM{aqj7s!1HC6bWki{usqAUg?1U1py_#c`T@?Fr z!_Hw7C2}pS;?;crTk@c>WQcQcK2n(<(;^yDq2R`-U=;nT;7WH%pLkkCSv6AS;5Tb) zZjz=z4sYIo2k?eCW$`=cyG*JDp@rM3ekg?t%#Q6w7+j@AV8bGv-=lORc11z%=;-!M z38JD~5Laz8pT!Ubl#&J%C(`mR#@p(xO-M%4SlZwewr&xkSh@#t;yg@QmhWffsvHkS z7ESb{^O}7g+2}hUS+_@Ity3nNZ*RD!N$q398Ce-`v33 z?ytr7diI5!uo`9~xFqB89F7i=k!G$WZjh}!BUVu8`{i1`?L*}yl9;Ui3U2W+im6Dg5th-5?j z7=!ycQh#P+Po_@zXDlaXAKNQkL1fiCRc<8Y_b?l2+Kr;yrdx$(wvcMcYK{(7jf z)E~VNgYjZt-L#PwzdWel;~)uDc~kwu(;*z%m?#(y#l%c!Y}Cn0pwlbE-gwxT{qzDP zs9Szwyh*MK$pXpGT++_&aQ*FF9B`zg->QEBf?_%=Pzvzd8nGd0l$)Y4yXGbe*(98( z;;+AHm}ItRE^@%^sbEF+#H4N@_2%^ikph3tx2G>ua7}x%0MP3%JeWtZm@(1rc#I}z zcfKg*eVLVk%1nySPJsVb{CynUE!%tjN_z? zK7^WAJuJH-{>kp?rQ(+?$K@({(DhzTew;Z_eLt-JhjAe%2scU(S`;=OkKf5ZXMtq1 zGEH@|OzpNly{U_F|E)O={$yj~PZ6{{Gihx#*PN(6uIkzqbz}+T04qBr(EgNMY3fHM z6^caWX>)EJaBdYSIHjK&FfLTHIG+1cQ7?xt4wJ{=sbniNqpsfYv7{@`Wz7&F_DKHg zFz&GE*1QhwWh|^Tzx^pFgKBle9+}6pn$e^SSp6vZJKT~AC4pb^y7_*w5}jtHCdHx0 zl+2FqPIB*0T$#7$GKUEJs({aNPP0h>^)ASCwx#hxVcpwF%_Edj`}PLPm&Y+*nydiI zZ*Bofl}|0ByW=os{+dSM%oL2=RD~+dY4BxXKG^3%o$5#04T^iE#%rx!hb8@W6)T|M?!;2n5woRBi;zns73VpR-h(T)9IqWr9g#oQPIrCsFmS;u0 z{@#KY5hb_J9M}&6bAO5fTKv8J{ZEIRPuZ!lhi!o&c0k;lXRCP(&nuZxl_C{Fd6Ubb z@s+c#<6*HZ5(ZeDM6?sh&YVgIY|v;ZwS%zUC=i zb!YpR&zZAmN>yY1mHQ@;&o5UdF*Myy=XEQ`tyP`0m2jFS*3d5^gy9psHWu}R==+p> zY5X1|6yoP#^Jz2^vd;Gl4pWUhD-(pS@b?%p~-wy z%3F{#&6PH?pr0@(%3Jra+UHy5{yra=r<=+RWk*NPk2hwh#{5a`x7By^q+NK1V(4P{ z$2lKOg-g|?i3 zUMtYSr!O*_>kO`sO2a*$8sISRkHhw1$|W$W-h9cwyi$W;X*cTh11XuU@grdhb1NX$ zXSF7uP`mrk0{vqc6yHN`bHMEdVVWqjs3vX{JZB!SVu+!v(XIKYn}*m>E`=!{M=_HJ zzPo?g&po$q{M@u3iaQK1?CyvD_N)WIS7-eqc9oP0ca2!O4A4pT;;~ZGUJ>nMHIG#0 zbw8nBkN+4O+%0${0m0683iLDkol~@~u_Z%aC;F1kN1K+Z3=C9Zgjt(=r)~^axkY|? zk2OrB0Q>e?pT-c`Vs}FZ<8aG9oz>P)g&+u+VRcxf=K&E8dMIKJcOK6PLgI+smY)%$ zc&|amr}2FCY!GwIEeB6Fy4QoFH-5zc0hVOjO05yViH~jA-9fsJPGp#8fgPX&g)7{? z1Djg~k=vxWO7->-wN#9#rLK&tGQ92Q-QUV-HR+O3*?#pnU%gl07ED$y@M}a;+Jrk- z5rPjkE?#`TKyLxyUr_2keRFXdMd0!dv+Ek=A7Jp!XwpN7=eZ7mjAW2aT2JY*jD;T` zRtbEjuEz^0-CF?Djz$s#nDvhRmwl6?%WV`i#qxV0J}JLtm>?ZO!a{AnuE?xsjyH5? z4j3PMRMqD>=7lWo%X-ZQ=iMLiLVhm6NuQ!I8Op}~?9y?ia$_(NW8w&VX9D&0q|Q_p z6%cb>`kv%lKD<-Z|M21Z``?(qKrzenZub&UcQE*(?zdGKF^s)|kxj*{iEQcb*s4g?<;5sBR z{%T6fdg8Lwvwyik;nGqHdHz<=ykVO7gFC7W$Ks*hY4$709-;@JyK=snOvrYj zDi(vj@Ru`gJMaovX_vWy6&Z)>SrW(<7!T>2cALN8E~`2zI}&#bk~Q}aP$iO=m0!Yd zGicqLOCg^R3_&^z!+$go4rHXHQ>9s>_1ul+jUQ)xn1@uixJN#$#5f2A$?U$r}hfO_>#*80U;l0yc0B zRpAmxuw_<&67o)5pl(OFUssrv0Hi`@&CgZIi0#@|B5M?{bV)tOG% z$B=x`b_>Ig3d7^t-?_7gz<@Zs{|VRHTq5~Aqm4L2@yZal{5UL77s&d`Fm$uw#Sbj||4!fpaj4?hkPcp( zyZ!wpYm+u)+XlY-x^Awl8y*ixPpH19@ef}+q1qu3iZ#zboPvn_*Xi7OFdW~|kCoRy zX4);68cuNeyFdZ;f6{n6Nl3^~{EV#vIIq7)Qw!3~ojf*9F%FeZ(SRzV=$o!JggajX z&ZDpTNHxhb3LLv|oNXr21kx(k$8C2zXB;=JC4eN|jRLA56Ysa}Vz?_!`=;bmFMjX% z$HZZ;=(N_7rvc#7$MRyB34doAwRY@!jB8T(iS{{xc}qNpTI+^Qr&Fmg2~0v`2PWpc|5_D+<=Qaf{U>RzqJQN;PB2k zt=&+9_|IsI$OkSwFZX7uO3D?8Ku5XbszC@|Uc9GD4wA=xwRwyM1fI>zbB(ZJiC3Pr zJ4=JwvC+@ZwbE!EOqwrTpk_!U0DMghz`Jy&jc?sFdr=mF-F%}5i+0rK+BUh<*7vi_awzCzkrb2~$44q2aqPf4)+w*OUf zc^}D-c)lt>h7U~wj>8O!yfO@;?`&mMPd^JqTbEYKUsymJ>l?_xeg9UNji)5uf^dez zIo;^fh}b~~&j=FqwIyXYw?15T%I$JLHl%+$Wz6utd?y#)&(y{TxD-7eVi4&KHpMZtu#r_+EJbA zGGrpTnEHE{TT@_#T9{CJivlNe4$AHNIJSb{K+mfQ;XFTra+_?B$#1F zw6%{a+Phrd%QG3(ve-MhV04xpSnz-;Q(@%w3d`k+tM2#BbHO9b@ zU^75G37U0R%iaTS8*9w1Ah~C6J#AY3r97*|0(Xro+!b^gX->9qZFK(+vwDvg+gyCi zIbv!8^rDF$URPQdL02>hh>IW-SC1o-4Gxm3tTf3ztAw%()d|g&6o=a9b4J|3gWAqHb2ItH-lvDpSLj(XBdkkp2~)1? z-@X{DT*aK*v!r^-Q58R4R?6R1W_0^2A1-EPy{onECqSCWlGX7Kf~G`uiFzBsaTRST z#VzljAvKQ1$ufb9g}ym^#T{nh6!CkEmnYb3#hWIl<)u5`0(w*KrOX61KgvrLxt_q$ ziWX&L&|Duh#LIYECsBu`ZGQ$8$KR`Ls@S|PDXiZDEfns40y(wv0#QC zUcngG##9pjpb)b#FHYK;%cMpHwFnk71OjTVcs!;8OMn~f*Z2U4&f4UzX!BA-YjY`3iU z#Od+8jt2A!!}^#KC-%foD#gYMton+dyG~HsE+wUBXan2sm(yi=YWs4c(%Yad?l5n= zC2yq*KXI@*bL%KaB5u^k+}!F__#Hk%X*g;~6l}k_xtRtnB?^Q3NxpW6(kdvn_E*hK zlr{aBsa4k*X=04#_91v$7$PVMUK+^@G60oTCIr~?Sd=on>$`) zAitBtgmd|K2nHDtzU1Z`e-WiND0lm^d)2eBuR01Kf5M44L+yURR+5<4a)Rp~328W+ z(|bsL1cRp*qOMb$W@liUbZ(qj-}ozzS$i&hATIm^pkiGkb^VNn5;e5V%5=95tlWm9 zwT=$78fC@vd?&FG<5WvSy1Usrpm*W!Ltf?=+(eTE_*oUNoEsC0#)$(Q(_AfZRlDrP zall0N%22_*h)-gE?;&DuKoYoRg7bD5avp$wup0+d=$=(dn$q7B?E#~mmgjw+HEbtIWJ22ds(5msY)4W9 zH#>ME0**kAA_gGKGu7w<4}G$KWMY=Y2Udn`+)ArM9F@In{9(cmV8`oY6c^%4@3*($ ztj`@D9$TvAU@~}oEPb_t-HNn5V^V+G$!tFPJRwM|hSNxxsZf~|`4_FEE zKN2BDB`u9>x@|h`Wm)of!)k1N%@}Yy;n}=i8Pm)t0r*PLi6HdjT(#N*g#B8-!or|e zRcR{jh?iy0hN)Xg4?(bX{*}SeW=r?xIuP5m*Jfx3ANP26_*g$$BmO#V_<@yrHuG_} zp7L7_Y965T?ery0j}(nigvZjmlzTbY7Mh^dloHO5eY6iz0{SIW0A&+@ zGQR~IJ!KS!uNUBoX1=@ytLi63N%!KH#?Sq7sd`l|kAFPXFr|64y`C$c%5Le5pZrCQ z3ddE%O%q4-ZhLfx*WRf3t#~+Sf9OzN{FWsVqkBQIedEH7;<^9@E2(lq7BD6p_!)@I zi{#_yOgBq!=+ha^{^wn?fpE!~-KS;TmS?_}BN+koq<+^!-QGc{2xyWJ&vnpo&MKmK*soPXTa zR^+3G{F8kGu@n=2{Bsj+HlsKpR?;B28!^>2G?^R62b+;|v(lT7StU6;%J-t8<_X`A z9%VDgV+-Plf&&&c`%h_(IW%q@V{ud##cBVE^f^QKgLi2+c@4rt#FSBJvg7o`EHcHb#{yR?Z4Sd`2ku<2F<&y4CU@cKu z6H(`kD!WLL9J`?DIi9mx$K`;>2Dj7=oCD-^qwO3V1GH*0yV;vYR1VcYSmmX$dk@pD zmo;inUsqUar7#-_4Z!})us5rb*XJ*iu=0jB9TvZW)&OlX zFFsSn5~1^0aZs?dNR^nRsn3GABI4zi+uL$`6&uBcQ1k}xfvgDi`zvOS!OjI`^YvXA zuZSF1U-tCg%8{mb8IKHrXQ5^hYYD@xTd3S@Igb6Y%F0JEmUr~n%OLi`n!pA$yp64+& z3to#GGa!ON>Mr*$Km6jN*W0;WxL!qhJfackskFQ1`2x`}@rHb~A}% znr5L^or8ho&}hYV?Qj-9ArV;rj!Zo?bHNkp^;{em2cOEd2b{+5MZPRaIou^#^#FWo zQrS9x;>ZH7TktK~u9XrYB7vR&u+DKM`-Z&(B=Ginb(^ihy_YJs$wLLcXmeBrxYtH; z6}9a83kV(t2v95&UPRU-qT2ww9Xnx7hY}+O7a)TkWQ& z*pk@P6!+O2-+kCFi9(ce3w$%bGk~VIbQ|pS9G!cAWb%&*A6r12T?M9V?WKpc_Wa>- z6Tjnr7+bm(t(SXsHHa>Lm#6{7p=h0ya7TjN2M9|n=V_p4=p%QxqT0M|8M~5lNX0pQ zlhx|mkH4?*kU1Jfc%-GZ!`vdz4j*(l5GuYW1E<#tQOc?`!Fn_nF_B)tE9~b-;veN` z4%j_-nfoR3^>ju~`t-=}Nbw>uQhd*D_3_+PMP4e!axGs(LshujlE%y%%Od}Vz2P9v z@$@ylw6U35%01SPeJY4PX7i%VcSH@VRK;+$G8l>a04R2r6qp(fTP}uRS#Ax15}B(k zYRCKC4g)F{8FFdPcY%C_t_V!79Bf6a8e06UHB}Q4{+d@rt?MlHMl=Zmn%lo^)qwPM zqfeJUU4Aa5QSQgLnZ;R7yU9;1zpaZFJ#9T2BR?!>+HmDSsC{y16g`9YWQTr)GQZtM zA*%1Cs6bKh^KJz|%Osk85yVa4mTf4VLj7JRf<4{pn)+*ix9ao}fOx{1A}@1#IKIY8i&z1)tJPgk#fdN{tdQ$0DcVruSY~AJW%x%0 zGoKm|(>0l8i@b@o6#SjVRrBf19C1~txqZ?Hpo|-e=~P?$@X?E zo8G2eM&x3O(1n8zwP%aKD7(nBi00IIT43wA+TO87Ba{jv3IUit^J1eMDtX+wvaL|- zNaK+&cIcmr`5Nqgq zPl?oD@5Sb@e4G%P76+Lg=laQkepd9iq{!+N1aMLNs<MLg20?EKxK9{cB&{>X#S1?~qCUkB#343xUQB5(PAoYz;)jdT?TI71{F zy~A5mXzYHKYS}-gveYZw|BQ~q+A3UI~Ki zWEvf!xdLYD_P$r?__9}*F-nFbPy*{-Y%o|_(Ku3_<;6=ZOZ?@0gDKkp4n@od-&jaf z87j+}7HtF-S5!H7ZZ>?B8T2diRcw?qWEU@r^nH%fWhn|>^a$$Rg8_+&F92nm_?T1L z#!95?Ok^`qg?FQ4Mg3YIFJq~1F3tGXoOsc_9>VL6MB5h; zJG~3y5m@k!aGhz|8og9p_8Bo{)L|{kAh0bTc3=m=)_fe}Y`=oXNts4?P@aM4lU1@Z z$Cl>r*EC{6*7Zv=3F8^gnSa&J9A&~`rl@gPPpRPsK)&7R}i@dnW z13z+)s)E!UI@f1$*c7v*>g6l~d+oTir=ngycX!9CkpQYjK6%^gGR=W-K$x@$i;y@f zv00sLLTf7v$gU`Ze*ajWosX+QV~D zXbfh82-&Y&iT~VHmx*ar2INL9z^62#lwxbGRp?R)x<1VL1HPQGtbgLzM~)>^8ClYW zTK?ecG1cu9p0+k5NFC#9N7g{;P%N}rG=qbl?QnLFmJH|~A-=a-<_R{4BDC;(uizbp z4OtPmD3tnrKbV?L$%?@s5`(^;8YldBIV3=wtnRdTJ=VEy(vOBrg7fON5`FlSnA9G9 z7i@asQK_8$uxJi<_bx8?Q0UYo*^kv-Fra~=t=4U+eui6q6%hAEyU@&Y(~_5u+Ic1c zlRe^nC!{f^Q7+H`nT@X|8KR5FyqcM|6$HJ$TpRr>?oyUM(A{`E=2u_Mf&|?Rur)IA z8-_ng5KWQMjc8OWkq;-RfsS)0lbrUf!l~DX@@WY1({3M1^9cjpTd@MhU7G`9z%iCo z8K|NPl+xGqK&9JP44VJ$a=?-^b)M~J9~YmxbK~o8`y)7=?9);xD*s}c;`QUZn^8-9 z0&!V>`&)eye+N8J1?dB(O#p)7RIkm^t~d3a2P_2dz;vuGzgjMFm?J3c4~7i))pJOq zf|S~ZRTnSI%04FQe~fFos(6>-BZ1&!)Tyn@0EcE3`rT4Z(Tbc69^;z^thBXSo% zo)9JUHP_cnE*4nbsq5c_bupGk(W&4}>b`Q|I5f5#7^T7CEYEZ@NC)3)5!rk!jlR}% z@V-}#@mV^a6UWw2RzxgFY5u4+zV^kGkG@xThTFJ|Qbx;I*nUM4eUQ`_2dofGF_U9S zjGs4^-1*d-!dv~zGQ&ys<2@|5pc-uZlOoVy2fjo}+Q%Xsy#;GZt-BF$N%*rO8Sx>8 z^sB2GcxOJjJvaH) zq%8Pz@w)aKGIP$kQWU*{M4pRjIdj=Deh zRu2XCOeRaG@atAQpq=n64sqrw$YGAnG8q2$BgrdI5z9Z2SvK5WkCabL0JFsDn{!Fp z6xSH7U&hxOu!eBSK zMaMWavTTh3Z`pw9JRD+)l}W@!Sv~9cd@ma7ocdFHd#^P&ey`T6neiw|%v^Hojq)h& zTDop-YHADsId0rRIp4I15>kspft@f}d?Pb5EPHq5BrAeP%d1{92)fsRgpqWaCCbv2 z%La4!jCH5yjBRQA0U+-Hg{d8-IH6KhlBAAo4Sr_NsYzqR@AtDTyOcMqSK%zPr& zNpNE{PP&XSedGP{xTs$z3y5$!>0gi7yX)l(NV%-Y_Ns{|vaZZC*l$eyec8A?gWBdN z@EI*lvXR~D+KO;LwK!~&f^lV^1#JymfEe0)GR^(Rj|3yiL(WBQ=J=Wq$+P3hhSumJ0>gjlWpIN;YILI zM0?HN-0H(Kg2MTpM=$57(lR zOuU|Je#*X8EfOv=wWZU+=%phMWrvi+j!d+zG5pv~sdF19ECea9iSwbvm{=XqL3Y+JW|*SJACg z=u|qJS6@1h-n9_@rrp#=GSZ_>!IJk?m@_oq%j&$2}DJt!#fx_;>PQR>*(Mav9(Bt=Y z@U!`;Bi5h10pSrBV~av%4O}Rcuy=3zBoBWo6HOexfJM%y`|S7Uwg-Qrg6vG(Is71j zk9%Lb^HoMS_Ot)jR*&Z*XD2NFa%F_c{D$)7(SDHalmH4;iu)dq6X^K_8<~cWX_92M zT=CfGmejMUft&C!VkXkI(SrKI)7xf93lQz(b(l;brZNuj$Bfzc{lJ3DRIZbUvpbWP zW7oANd8(IteZJ6)x_CjvY=Q|kj00L6Dos^CY=khMjcp{OLnKD|yK%pcx-L<)0Xn?9 z8B4*{9ov?Vwi(-Swe?@@i)U!ZTmDE(SlTGyKD4PC(H(c&B%IOG5rZY*)r%GGzx!_%Ex^d8gx zmY-}f&Om$JfU|yWSywdy)4sKXm!J zS6mg-xHRgt8dSLA#nPu&8p-#p*x&)j@~fS3*SDscVnzdYhDwGTq;$StnJiPQs-we} zi~t&7C+#YHBAo-%+$|l;obx7t?j+N&ab#v;Rc7O|qOG=_KsMLvFQf0IFgfi9+*+)~ za}F)kl~S6dU?-b*Zy-Q`2E@|K6f>`>AvYSsfm3FQ&*!`n^)>`7-j?m>0AoLI)K+m_ zY&fhq{Nh`pYH7>
  • FSkdN(+0BU>V$oYq$3J-vu^fZ6g#VRame)ZE(_~g3w23`5M zPJywjs(qYFG^k*@;RRl@Czz*G4d3nK`ZrcLClD)5y5AxrZznd~C1o;0-T|4`FQf^p z9`i)HRAs1i7smVNm%nT1ZY>#A7?WqIyy9yniOC=+U2t-L7Iv`VIBij1)piXHa`Bw3 zzS{6?#`?P0m8T>JIeyuFcDMMwM^%|Bclbyjah;cxDw+v>6?i$1;)Jwh;N*F#Yu4Y{ zBUSSG3DHrNjYC+wYCnxjkX5`#@e};*g>(vv&`mGe?W!vwJ*;29uXXLjOJ9u%UtxdV*i@j2);Pu4f!*c2S7cjY% z@#8UxJsc13lEUqAOSMYZY5H;~XCKi8SU+g2hm7nL!N;!-Bat(Apd!;MO=@nxzdp6B zD1D%}Hm85E(KxL;p+)BC+{!8t==D*L6^>!Xq(CITPLOb*>Flj|6xTt6f>S- zg?dx;t8~GqM4a~R9De_MU}8Rg3tJ*<)jK6KWo!GfDY1_DD6cMrr!ovw87a76X7)VQ7@oLSS>;l^0ms z5lt85m7f!CNs~>Rv4|~GdkO@2q_qW0gx(;aUWHGqv-o}`$9nA}MT!t7*#x?U>r9D8 z#K_}4xr@3ku{(E6;(`Akr1=eZqegkha>ySkc7-WFqc#BkHj~1)MgIGba%rDfY+Gl$dv1+2- z?=~Ox5#j7cw|EoqEF~Znl%+{@KscVW6-Fwm6aadDo^))UqV%wDKv1me?A(F* zW=Mjr`zFZPBcaw~X|)kd7(HS@>DVjeukQ+;S-!k3?;W0X25Xtryn+@+ai zPBWbaGa4rFTEsNx4tq`%x@}o$eFjJ_CJH{%>*)`XodRRSrX#lnN~X}ff|(uQ(|xw) zqn3J+`2F4QNy@Vl-i+i5S!4S&Qs*FI_ccDMVCS)TGCwb{cpa57 zA&4N>H($Fo#9no^8xb<~hh=-N_(7s>3&INAUuG6lq zVBy&ofYBZ2mk>{a`cv8_GWb{4xQLLZZgL&ebigCJ_It>YCCe76k=Ao1J|lL2G<=SG zGZ1dQnXC$FnGq(|qYKQ!rv`A#Desl$g^lMTPVQ*aJCCS+$`kn9Cw{xK)hUmEBC16B zy@j$Ac)iy61&JA5$Da!mW4#Z- z$KLY#u+bFfE(p~~P z{v&_8Dn=%?Q&Dh$ra{YFZ*%8}TucNE6ZR_EeY==s)GZ6xPI9WOx8%gEx)u>)T&g-_ zaCL-6q7kOAxiFC#2$d{sR|44Wi?{0|4*I!Jd(p!y+J#o$ar{wzxNmUsji+x8t~<|g zx-%~NpJGoRi1lq{-6x1v4T-Y1pALF-oqCZp^b1ef^7ARHe45!hI@(1u5IDp1p3IDb zx@96c@H4_XAZ7QP$xh+4`dStm@8LZA zEOBMn?O6P=2JUjf7iYpLQx$z*_ ziC1o{q2$^{^)M?qf=*^lr9o}u8LG0hV0*xi-ROtm5bbdp4Yp%{PAt+BA6NcI+ECyMue(1=0p> zDLWKhIlB?=g~5|mkZ!~MrMd#!D(e?BWj;U>jj{3`XGd#qgUNC+O#*x|Zo6=Mz>(iE zsZs7QLFns9l1rH0vGfbP^|p72U99G*jax-z&+*k+o!;iivee|RzPYEi4HY=fdNHF|G?ck9qL{6?i_X8C(RY}D7ln!po5jC{kF`=W*u@6H2Ak;~z z)VIgge#*Ht573#vsff#mO2)*>v1uT=I1ELWDgo+&Y>`sK7jeNvdZne&Z_)42ZK=ds zgV8sm?#Hf`Yi5dj1$HqCN{%xjz+Erd83k-e@1spXwR|`}-)}{y)MrFz?SMN#Kzmx7 zCAfC}hZWrZK#^KK!ZQXuW%3-K*@y(=vqS63>otwT7D6t2`Fy8=`gE5-oNZr<>}bPF zi`6!4fH)Twq(SDZlO(4j_H`W8bBeXM*)g_GDrIK_6dm4@x zX|v8oJ1=lm^5N>mO&iSv19^Z7KwfBQNzncMcC3avBb!0Ide4xbMV{Al>;sJls|L2) zcQ6A4TMd7&@M_Fk>LZUrmZ%ObE;y@*poOrAH`wSG;}WFJb$oy*Yh%C)_=u3-LlwJ5LcQ&oDC$jwCGeZy5_WSXm@Z`X;L;eEF%$(dm3AMnx< zDF&-5;(71q(0tz50WVCa(|V>ww~Zl>#)YG6d4ELjjAMdfUlySvRk&q>cz_s4%2GV% znhkdP$(Wh&_guzJnUsqi-e1$&B|EBEDpA$l2Nh;htCSukh6uYhXh+JdLYU8ImAYI6 zlis}YpN-m2>CLqG+OIC~c1`8^QKqpwTJczOV?OijsGAxuvDfMWj_)dD7#g8)e2~zp zxc%xAM{+SOiBlczoeYZ&RfV|ty(eybd5z*@%@dWz<9f|Z z_cvNl>Ynpu4IjMQ3qQ{ITLpb5$Qm`cQShRr!dQC{|L`5yxa3l9n%SqmJKdHL+2^MT z<$llLjh)ufi|k%{<|A66^%n3u>AP>3mrI6Xmdufn_{3F4u9c9ou@aGVdxMhPWc)#| zu-15LXwDK1iVTy#C=9TZOp{Jw&|hjl8kigk6qB!+nfG!+bgrtc)a_Py7Ls9T`u3ng6c+c+aqE)zZ0RA;ezqNNLHEv+w8&ma%wjR*F++DvbSOcyxqo*Y3?k zCHL%->`+Ljuv#^8T{j!>V<0}uO;fw5{$dk~$| zXBf`3h&-$`W>_rE+PJAgGw#^8{CP+Mcm2KDThBEY02y;-kid`}{o;7F-8Ss?)F}Pr zsNp@_i$DQCq!4PKyIt`Vy|cz}115!XxAtZXq*rMSczk=wCfZLBRfQBh7A)!WA9?24 zR|(+$%X!$KBTE2p^Kwpej!Avd=3S8FJSX`eP;m47lC8F}PVhaDdm9ppdydM_tZdp>YjDO_&ZHD#T#x9<^@j<%&Dp?WC4z`wjezd(5fcj>okH+W!0eu9dKaHM8 zZbP|CYW>a9;_xPLnzK}OA?MZ2XC>~dy3j#q`iDePt)o3XB_yw?tLZg8lWQu-;?WAH zEJgSh-*oPq8DyaLj^Kq<`CL|IxD30|5rb898Xu@pi>SME@zl)J|A(`6$XW_tb)1eBRsp{eE4Z@9I2r@neZX^b{UD;F`E$a9+kUP2Yu= z+UYc%ZN*9UM0r*Oa6zjbCwHoyw0C;1OvGQ1538u)cFQ9T3?Hk$wF?-_2!l)IJ*~ET zftstd_e>x==To~!coRjlOy1by9UU`OUR^dRr}l1qEr!|~Ma;d8bHb*6YzGpfq^9}NLa3@YDv z(Y%W+U2fw@??HyPK@0BqU^y@LAhk}i}#B{b*j--;uVnVxp6w6H|dLXs;t(m?J zUF-{2b;T+gt~rng?oC2MTTlgdT{hSd*UxmPl1Q6Gml3zJ=bRtvD?gZ_<2nM;fW`5O zns9Z2NDyN{zd*LTp-yqi9nn^x!>WL>OCy7$3}UjnAAcAL1x$~`&!+R*g(jkPyXM?n zx`bRisi;$ZQ1*R$rfQgPjPQy`C(L3&8(V{ttDKW^NI14&{G@YLr8Ea-p>GRiNRZ=| zE59~Y+6n73JS`TS@tOBB(@aSH(d~wl$xRb3eJnq0z}IiJRZIPN1MmFkHIHUcP`XMQm;=#0)F&Mr7*12|U9s_Zgzx@NiNW=1quz*a z=a6qA3z`Sn*6!}fhX~^EyJN3|H{f45*=wrgp}IDTtvdj@ z#F8RiKU#^fcRs_l=e!xG1nwaB?XLfq`fK@(W2uql8enHH(D+09hSDOIma|Qkf!6pH z%bIhUW5fm8scCU<*p2af>&7DY4=%SX*~NckekMoR*VF`Cp~pF=`+{QU#))sW)0!Lv z4A>QW54KIT<6%TZg3c#O%v4R%ySzneqz4To_bl`#m6Ank9~@yNhVO!U8jK^JnKWF$ zZu&00PY-{y`t)XI!^O7|ELpiR*Sf1KQMS3#aoy2t$)lz+zydAhn?U~y-oVAxX4>5p z&F%*ZbWCtB-)ti2n|<@Y4Bltl>HC|>9%wNAZz+ytpF?^fl!&un_a20%>h3Hj8pt}e zu&&)x{sbgYU9$D2Q8E@nq4U~l#tCJ10D({0Nvc@maPB59MMoe=wvt&k)Gt`SEqmQX zZ2m~Mzs}cvl8*4lx>*HL9)bDy88Nb+j`qQ5j0Q9zIkHT;HoM`n*ct@kM`t3B54Ab+ zUsjqP;HzDiu|K2lZ;2){sipdxaD&^ol40xu{!Cx#&e)2i=#u$$>gOpMPvxPg&}H`z znfxU-3O-4Gi%YwncHITk^*%s9u1Qxs>)f$C8KEdHYt621cArTjk24ZR{);99$IUPQ zXdJZKWAPy{L#jLIu-wloo@VXY6owDSKB~I(3)Fa<6y%TlBT-{Hs&`DO1a^locLf$3 zq~~sSo%#j)s_&bBb>iBOItYiyWfR;^%|Se$A&%3;OmIU5CArQ1p#rs8HwLwd7`Uwe0^X_cQX6f0v*nQ~Mn`8g>&$fPY@!a>DIbh%; zH(K?7_)ni5-{@9r!!z7sd0O_M*|gT^`epZ$d}WaGz1y4@#k9~%6Q1rOwiPWt7 zk)2iK8d|%zky?5+87B0ocgkmK+il2bYMgqeu_|;SF&HvyR6pqAmw+sEsWkuB_#r=3 zNf@ivCq%iI$dBY=h#>!^~BCwb6u{ zxveM2v7L;&(oSht|Fm%&K;z3OAT`du-KxN7em7t>B&Y4%htLTI=Cz`#L1(gz|9e^4CVo zD1K?8qa$|9OL}yXudGjh77QuAC8)}Vug-;5P)3rA(^z0%E55W!csL$Rv?J}B{1=M? zU{lEjx{ykx_cmEHs>M8oC^qCo*S*vWSGEMPe<@Vt2eK2z#jG)7kQnA=qB%@JB$8{dtX&} z-XBepNgl(y?fx3ZN2wPY@4q{%yHvx{I)aWlu`>%$=#@S_uJary$e0?^_MB&aCy%s> zf2|62wGCk8;nWu_YdAUJf8b-iW0&dk8BtAN$e)bTsP9Mg!;N)Q@71NMEH>GQykNK; zGbPp`dgRT)0uqNQY>AKb6@GEch;`F>@%cOJwvz^xH`O|SZY04M$(Qjm9lt8?@|z@) zhx(Wor>e$uo=8KY21Z63<$H?&oGM06rW1BG+j=aKxda;$iPNT-Kb4c)ndPS*6OWHl zo*RM>|B7gov1B@t+Sg+BuRC9w54tEsH9x&Mg*dzHU>Lc<#YcO>c80GVvR2r9)f93( zVi(+DbG~wpQ}hyOWHhO_OugIKJNvj!;p&yc4WMDZyZUr|xNDR%RHk-AVgIz#<_mn5 z#b}`$sjwtxYY+vxqAcQSl4AUAUq|JGMSNbJbf432`VtAb2Z2`Wh(BigHW3?ao>)zu9?D;Y*B~ zt5$^xR@>7bT;)Ntv-IoyM{WXQ6$jB@eN=K|K;t)#zk#j*Dg`x~bx$+*_xHbV$`lIy zb>kWPao1Sw@JzK15P5Q<6SlEm>%NAIUYEH!MBMG=$SoTC8yln92d;}l!)^Rg9hEGN zuu};tW*;p$VX_w+2WH@^cqt>i;Vff4O^o!?6K`lWo8i7x}NIh}*f`)Yn2j2zdKAM4t|Z&(d+NT*`ySH}zaxnB%ipmS=1 z3SYZ0bIbv8&o$o;8vF4IKBQUq*G|hdc1-BpxWLxY1iUOf5RJc}vlH|3Aoz<%Ec5}S(&W_H?KOKqtePUjS`otd=u+^ivs?xeldFY z5ab0>Fn#pwpAAVgE@h!>{BQW(~*ixj4;Xr_XoQVs2vKMlY zFD>3k*AXy)8hMX>N>IoFAEy{4W+(KO@4Z1(QRg(~pzDM{HOH+Ng=f-U(42K{O{nC& z1SGPoeQ$RCjH5MB!rnHHTh-xoC0JfhczhdeO)w;L@Ma1LtXeZ`+`zCZjM%x*sx_;4 zkP@0vW`JVw)@u?aScV6#g*@{N?*3^gU|r(mhSZ4OoGeAFGH;2bajP5Ta%4k~^pIQ< zvDox50%U+PYb(zp)4>G3#g@3>T2V=U{Ugq)f9LQ?XWgClV`_*}y+_l2zXc8Wdg4Q} zAm0K-7Hb?Db&&Odom6jB$Zv7Vz^qy}Jg91g_(N~qy3dpM+b z!-V1M*JdQGZ&T$m=z7Asx8h_5G#XYSf?XePIwZWE^p$hQv!vI5eN*a@YXf zzj22vHCGSAf+L4Z-#`yNu6J5o%p@>lM;LdfjaQD(M*H!nLr_s+O32pB>C17698*jM zR&ppYf#GQ`YR|`P^`$+-%GAf;gcq?Phhl*rhW3aaZ&FnvkRHPOmaDM(6o$TLR)e2L zz*4l+Dndhb`2lIs_1iorTc{{4wUqdgRh7M)AhZ<2!{yYp3HYNnzmLP~dqo9$Am77^s}E6$mPEDyve*dG7MprG1?Z zu-7a8=B^n|dTYI(q@%+;apO^T?~2!vjMK3O-YvS?Rnj^-5`o_TBh5`w>qT(FnGZ2n z#J_9q2cE3DLkm+lu`4&HvioyLM5a)V7lZN6XYBJkDePWAwG=J?C13zDrm_Hs=~I|e{*s`a6ZFJUasT!)0P{q|&*NU94SlAQ4U7U= z_4zx^o}zk*&!#BzFtX+8oU`+9vWm zEqchy3zA!xZZ&;MVAN>XKu=?PFM)PpPqmnww(7WL>Q2HBpmItDXtEYyS$F^|mU&l( z6B95)c<|8MJGF2JwHN`#^^^lhW!b2Y_sRb%_w*MRVWt2uf0_jE>FZ4iU^;GJ;#oGY zF1exlnBQ*Njz30b{l}+uNA#agsn_~QflO?d+XzXSDSi;&^ z%5vebC-{DHc)C!A>2<0xwVw1#orGBccLXp*gL7i+Ej=t*k>tTV)u{Xzc>jJ=89AyX z)ThXtb#Mo-L{si5tUid$S_RTVu+kmQYZp6}dj@&t3t8tvx{asU(sQyl%SQUaRx2~y zg^;_t-v&Dbb1mp4DC;n4>t4ta{Em*X1*+yX#o2a91o?xO**h)ecp+D1s#pU=(L*I# zKm&{!7Z#1%Dq^6dAeHA4fHv;qA^118IJiB^)b&FvZjf*A)x2R&ZGXdGMcnC7GGYuuzRkQy3am zLyn;-3Mif1kBb{M0)Au!?7O;kJaC5aKE6wj^K~%%b>XJY5l5_YSn+?&hW#ogSW(;x8gZAG8h zw&Zy^;@`gPE3nY4T=~W05fjwoH8O%mAmPu)1~2?t_7E7n7me3=v3*+!uYDHjN`poz zACAmD(jM@&kWwdMaaK+Nmk>;Zu0LRBxb2iF&RyD5r`-0|#3J|tT1l+I6d*h4scYU7 z>DWO^Z=%j;|9ubeC~45vQfMhMFo$bO{^J9B${{rv&xcE<%Z(Pge#p$aah|Kc)DCC= z)k#Kc_9oNJuXzdbD}`goUh8ZiKj>=eHaL{{V3>B7Vk`6Dh(-3;#JK&n#?-Z)VyKcD z4#Y>ovb)m#uQ&{}?ia%g(Zw8dW)ea1N<+X(s!Rhz2QpC*8fXUoeTOqP4B4r6viQ=+(Ck@ytBO3`1^7uEt}eY zcrVY2ynIVwOK9!dykcZ>`fE@QXe!&KkaGVYY6?`DI(7A%SbsPP@NsZIq-n+X496+7 zP^G3&E6C;ipmM+Ye1imIyFW7V!YS%M8L5N_o2vK49izNcNN>Mln(gF?D*k(rLZH15 zl43JTeyX={d4HO)0W})eRvt>jqk|K?O^%*K1f4yc8J=;(yth#{Xl0D1U39F!#6VB? zox49OD#5G%GYXiQC^jx{)OqfAzce}hYKk*5ffq43=Dv`knXS55K#SE}tMoe)Uu>EZ znN*c)3n=!kc|gsOW@=!P!?{SIK;iS|vr*o><4NY2r9zHbG#IO?=8nED&83+>LbK1G zqW*!+SNqfvpSCSG(&pRDoxyd(hNZ*lgmA#)-&DO)5M7Az6&m%KbM?OQQWZ+i&eLkV-z^6y%LMS+cf}B^U^WnQXdI2D3uZ zme;9=JFSWGP32nDmItipoeu4_RV}isx5)03fEKNJSFnku+Cf60k5^{e>;?_=B0ASF ze{ILG?#eYPl;1x&#U5ys`p0$f+vaF$D8pj>VV|6ZuBLzfJ*}&0~4q`gB#!(8a+)MVdj%vqWf@WK>ajW!San)$Wnj z(#8Vd(Lv_ntXAIF`C=Ai&2C_U6Lfs3`nNKAO%DvKR5Z{vQa3=i2@dfm+Te|)h4Fe( z*Zj|IpEHoYnVRsHI)!|-HTMGbjv`zM$Ej;Ystblx(xoacd?GG;A~K+`un{Q48koec z7_}pl#3!S(p7p+Xn6=S>xyAT{{CO666HWx?gsf>2KGYCcV{@a;+z5(Y9`&e&kur6x zsa`&^Q? zK>qz;A+$O7uaye7p{R4acLJgo-ZOvN4BBkZ3^ob>G4|4Z#k=RM$0!;epT-+atiN>A zu_yba|9E-PJ}WZkvjGH5{B$jVqz}wGL{RBMvtHBhZGd)OK{D+O!k5_Z5$j9Ueo(%7 zc@{3P6Cxc1bKM3PAbFg5kbbF19T%@C%1 z*-dw5dhI7j`pRR3Ql38x4YpgU!ez1Esc%(zIKV%8RJg4{8oCH)_;Tq2*m#!%lx&0E znQIt;tjo&etqlvc1Ys@fLQfvPx?Z|oQo3fM&==PGxR@OI5*X5A*a@o^FV>KNSiypq}Q??qY6awm?W7sDMO|5 zCcn>`#af_-oUdB``xa`roSTXB1RlCtDXvwL4Kxne`5{MF9-Ew#4ErSnUpyk=D^+tD z!pZd?i&@j4R^1pASxCX#h++Fh0K=&rrPd3T-5XEsm}MEu-crGNbpGI;O3)8}v>*~# z!(l!sPuR-kXn${@Kv2DLFVgXa;)d3|wL+k$fI+G;ro`wy%lJ$mPLlhLFY$K%4S7uz zAFttDD$SfxI2k0T;tJIq_zE+cPTjl3oA~Cil3tTv=4iO$%lbxZRSu^Dd!G`1S z$g}9}B!Y|k*jL%bBDDTE7#m#bzl>Yy-1Ckf%S8^!(r~>CyyNe}Q#QPrC5zo^R14B; zgFNn;X{TSMjnMCj-F)<|@?orKchw3kT2%biGxl$u*}t?ug0pnl59rxSJhco=YGigr!aFM{dC`hJN+B)C(@WJa=gAgM$w=Pw`J#- z(8N8WNXb_@nk8Z=N7LXa<9TGZ^f)+cIabpW`vB1SS={nNzH90tX2ny#?e0_ohVDdj zK_TGW-$p6)Y3Y!um}kAeOeBjJ?g)nLm0}HDV-B%avk7;pCp!7>QBj~8lZek-Zz(x< zf-KC_uT&7z0)+~@qKP4nZTeB}Zs#2fr)FF!f~ zjmZqZ@-&_S;(57VBcftK5@l`+&K~DQoPE^{h=P`5}7tHX_9b4QJU4twG`cPHg8|dP$*@h$kun z59@Q@M70P4OY7(tQOYy0SABKyB{3mdh;~hf5<^@%8Vyft*9*9B1^kXe%?s` zF=y1VS}4sWF3dOMdk?UW!=F5dw-7IbO>}i9r_2sP{_P!3LH@8R&%xfh&i<4591@?A zD6#(71@JMX=lW=epp($1O3-|iKeH#8D;i59?VsC#{8HX{1Fj1d?bhl2L64`2CZRNt z)Yg9pEX(PG2OhAOWqkK)u?RZ(NQt%ZPgpo~Be&tw4Rs#_l>ZL!a(3|Vg}UuI^*1XcOp+39jm*Uo>s z@N&`jM*5$C+!=$GkNn1Jy;VtTB@E8g5(=G;Wxc(~3Ws(EAVa0Nc~#HXyqf0}?JEZGjI7dZ}f0u9jR<6!M>RO>%_3 zf>TwFG$FS3YTjJ!YuqLu1#wp53c5ZLHt;A&9o>0@H{oQTJA#;^P^t<-4&EGE6j_(Q ziu`S~{RN2Ts{pa6uE=igtF!QpX(zVQhyF#x+wa))IhTL0g=6snFIn|w=d&!UiU8V% zRnYV%{mfgBKD($lmIZ&8@3BB4LS%HZ!i%0z;h z>a>zFNas)qG8+_WwfhNjP!l|)tGK~x6ucPjcJOh@?%PyNJsQWWxOG5wvNe{M`7kb?hz{P;LV=Yb+EwM^2H z|A>u`R|7h&x_;hBY!5Zc7O5uT`^F9XU$1z|_+=N78HyyIgBs1ytN-n0-32_((HqZ< z-wATbmT8y*n*`AJAyB5+LXHxPifhBNua5oN#47co>)POEE*Z}JHJT3{n;m-&PsEj{ z{Ot7wx;oF*7P8SPI8aZ@_}xeRZf7)e;%7{;;*#LahQ^+d-5=bSr5TQ?vbExYXTohJ zzFZm$SlSNGURzLSW~EjA9C1Fng?=}yw|VGqD(J(DjkIoAjeP$lGK3>VTaM{0!^)FU zmkn{Ev8xf2*sm%UuCSWkxUxzOrza|qhPP}viCesCuW=`5Id+`18H@QA|9C(C)4EBC zxd%!YQ-&+^#1l|1?v9wC0-WbIR*PRM4OG*lboC}P-pdDLkAdDsLyGpL90t4fP^ZqpE!7Wtx{5p4e;D-**91R zrO`mrRAVqn|HG(HnnRz4qPd|PMzfP$+?-|$c$9c9pIFZ`f5|quMl3$TD|UCz=HFCu z4dso;{gQeOirgwbg+Kf+X!PLct=+$d+jTS4D&s$Dd@260nILY0P>{Y`cE4=xMK{ga zP8<}~ZbQ95a7|m5&h@iTO|U2$L_N*MU#8MZ1TPxi_M0J;?)xr7Pc9x_?EL${A-5fJ zLiyAi4nua=ay3gnfoBm!fn7^AZ-w=4XRU8}2#h-0CGmZWmq5{>kbHG#_P^pUN{AP} zaL*Wt2f2h!{W9On2yUgHjRMFokK)3J=-G|iD83MBksU;Ok`o3mUpje9x;DZWP^smh zPfB=SNQn=PgBC#JvQ&}giw*Kv5`X8)={ZZq20xYwDRU!W?a8w-E4EUqNQ|23sQ*mf zoXc!21HZBC-jOHI><9Et${XCbPxBwni8NPdFX~0sv-l12FCx#Y!#-35ID&k@)-hPh zCfFK3Uq~ZdXe!M!mqX`s49YMb`oZf6IB$tRy2!Bf$Q6M@rsk>wzGA>3K8}|o3dLTw zta%9h9UA`t5W!%G#6P(bln@LlfpAJVVI7j^q`{lhJH0)l-M#POchiT;zpHbNw9%_g zk52;QgO=BBl@06vs&(px)(0J2ugcXn@jbxtr*9!(YwBq~foTPQVZ#brAJb#6lhfy> z_MIC4x5=N4)Ti2XMM{PSQhsp^*yrLn~3n?=oR4h zLtNXb1-&-O^zC7et@{vp-j|ocN~p~W>)9HpPW@>PlyM)*6AbrO?cN$$PS)F74?MOyYBvc3K0Lk793NZ9=wJx_fUO+Nnv#D@oQcx7p5CR zgo@#g4XGE_xbUpem=?j#WkG&afusNXm5&{=)Pno7M>E*U$^A;<-fd5&SgHGke4=?u z$#2EP3lcsJy|7}mC(9In$)_;~nL1LF?dMVT+CBFsEhQL|scy++y(3qXe!ewP;cv7d zu0z0CTq0EDwOm_71;`=&h&jEq0en1AcHJ|9<4l*QUdCEFVw}=BWeyW0h=9C+!9iJ# z^9q%d@oc+W-bG4ihktk{WhC`R3>1{WFM1)h9eHUA9mn49gdm;1o2;ZoyWBmvmqnJF zPuq}d{MXVeDE=YZC>ZK*SA$!4zg8Ztnnbs!o~+;p!;;LJnV(esDpobFWmZk z*^B$iz(fdk7qm9<$r}kLChX9lNe;?|;@%0aR*GwL`&8rK1|vV12__l*{(`Tx_bx27 zE=bMwv9L>K4e#L2Gn$TbB(eDiCTwxLM*ApRgr(GeA)jKD{3R2>WB!W34!c9lB>>n# zKqn6^w0iwI`+jDl;iuY$JD&PgVH2`ukY|J6VFu5eu9;w<`ZAt)KB{r$1npCK0=e4} zXltT+2op7S#v0)ljG}Jn047)#=ARb z$+A8d+6_j~w8;XCLopdjZ6#$NZuH`vbc56OFQPRPSJ^-CO2q!Cmr!fs`7iW9bP$@8 z>cS7Yk%#p&)3UUYn#TO41;;15x)l<5`&oGr?&%#yhc~wA@nK2RmrPN&Y!QQ1n(7Ln zrw5L%WgkNo13ROcg*inUzX;3Yg((4fU61^xOSY#5l}2V16KS&KiU*zScwFjwdrs<< z;?kSdt0IeM)+a|mC11fkyabaaOIKBV_wHC-CSC(NjDE;y-k)!Q%s+&+-{3}haQQr|LrmqS{P%B6@$><10xT*m(^N-p?qJ2 zw%4>}13eg@4Qw2Iv*Tf+mruviihEma0c5IRll2D0r4x5&`cs~QZ~uabfVuU3+^gWb ze&klgsk>L2HDPOFp_ot^9R_{h+s%pH2yjW3T z>KD2pZ`NqHqLLo?5kr!bS)RaR2y^+VJ_m$BS4Qh1=G&t_>@11UZT89l?2D(!{NtLI ztxr~Y7Nne0h_bbY`G>>#76eMX@N{zTh2FZ;Jn2cr1-+Gr=$63ly2qFPG2Hzo8mcG? z@SRKEdTPn*d!?xM`HEQAnEDEkX26L~0e~u{2+yElUt|Hiwrjols^os*wF?_b-zdf* zCZ^`n|*%|g2<0Yag;i6q(Plvt0FSuP-rdQTmqX^ zKv>W?-m-g+I(OWy0Bx{gwz+o4N}xRg*&D@@G~ zxQH3I$UAz>toN1#$y4*gTGF7kh%0F5*VaEK7cF}+HoUrl3<+eB`4jl6Ma1`Z5M!ZT zs1thsk(Np8uxX4*P|$p6cX=2oGnNszyu=ne?67KmzPCj7=+|<_8}Jm(U~CZar@!@@ z;Or-ly`R>m`YtqG!p~0~rfHwOeX|w9z0!NNcMu)|6u+plS~;F?=^p@nyyfDpxOK(O zckf|=g{iCm_t;_;PQHQ(ey)Zw=IAQ2K^3ommTcvfV2L2X?hPibLgq73T$R5l2m$d zyoM7~Hdrcn@j$EXJbmUzE!fAM;53gZma%zF*66H8*$DHQEmtVt%B$k}cSATp0oWg) z8c*2jTaC)D7N(?gin-Sa-x13M=FKz~}vD~MT zORD%xIQ)|R!SoRm#%cZeest)P(b%;C@g=(DJN!e_LvvUkPUYH3Ov<$Z>0KT-e0W1X z2-hO($(VZMzlDllbHJ0hXMzH)9~u4Z^$=q6`|IwXzp4HG#jSX-4C^jeBBIe&Q+=&1 zFD2CP*CF=y6nDwhk-!c9fRIlt-H5B}Z3>!w3lq-5^%^ovS8A&4ye&#NmcSZJpw;)g zG-R<3^RJt{A@=E_y_*Qzg;@asY26T|%$RdO%9p+f@+q5JL2h3W^(7H&dq4qE@w^v8 zcbB$)2GYGnssS-f9=evn`P96tK^Wg>2&?x%~?0SOg8Ehcd__aPBd z>8y;i2)^Ui0yI%-M{otsEC_4X0a3wHgAdiX|$x?_eI0ssx z;^B#SE73J8+83@b)aAm_yfA@1`;(OSGlQQ)raX}wn!-6y)3T~+P2pmr)b^RFHbqP0 z#)Te+dj23zkk^zquRG>Tdmw#`^Go|J*%6EC9_`=&9^|2ZFS@ROHQrUyqK}z3|sN9PEL-jXxp zdcZGvmkhnZ6=9jq&U0^TS&7V5ji+x?A$L#2{V0*g>*+dkDt5|ZxCHXHO&Y$9PiRu4U(oFGNKyM{Lv?MNiLC*Q zdd`)4Pr8AJ&yd|JJ1hxMsB=Y)4qpV|7IIRjC=l>iIc0KoG1SVjPyZguijdm}y7-xUF43N?l$hC2H+{J&<9o;-uFn(BRdr9`jKVY_M8teyf+H zr~+5wP1u%i@R~yui{z41E09vLtKJ-d`D(u|{-#lr?B8Om+f$wE zh^tvrpmsltpJ=DnvRIc(INV>Y9++rFNqk8!FO zr`_~m04FKJao)wQDSCiad_8nD)#R~iQL+Aqa8$HFknx<;>)rZL)mW~{(_mq@Og|$~ z#S|!sEiK}yr-%s6MeLSzn;SI;Zix3+HMD)Na+{=8zh+hj`5*^pXBdYB1}`CZ^{K?-bI2q@l zhSC>uOtMK}BW?bQP~lEHzBJ-CGvz9lgYdub1tz$VMZbhq?~7)BXsx(}&8Z)r8XnDD zXw-Dkv^VGHi)1|M0bw6zV!=^5P=3e0=WU*%$pr`D2&Y7mfe=yTy!l@8=)V9c5YZCd9*%PUBeqN zRMmFGvFzzu_STo^6EWJzYgzKr;-Y@Du=@i{FHpKcBf6|cUY02qImy={cPwHuIlip-F!bN?DC@8z*k*}Q_Pnf7J z^iHxjOhbdPcCY8Jp0V5^)GQc(lB-ZbgN+7u)Z#o8mxW@Jh3cd1jID5>i-Ts2QI@0r z-xadk-pC=v583Y3vZ}bBwvEp7O4+4l0cyYk;^$8!#$%%;Lovf~^2lTN5infM9c<23 zDGj|AQ!_G&nlc!DVN4zTQPZYZfaaH=@A84MwdcM_v-`AVubU6mnGtP}n$g(@TaiL} zBVuG#-sr|<;c(o|#NvsK&KYbFJgQ5Mri;xTonbJ0&_10}kKRxUe#$0n`6u}da`sZ3 z8z*h>5od^bTTYD*j>xgDHrobyodQLa-E(OA*6Z7Tvo;9AFHxpF}%zaIvoY7&*cG8O({ge%o!u?d7f>hoUC+CKIXgypB&e%0Zm$ z41zjK!iL`LN6NbKbNlf1QOW2+)zxKHh>>I8QCLyK&NFRZ>_m-yxd^Oeyf~X-o(O)X zd^J8J)u){CR4@Y5q)w4~9FqHrygb7P)|Kp7+4{x($X#ZPij_IPUs-uDwKA9zVYVVg z=yxHZ6U$sm7SfF8keN2In%*wns5KjhI+5~fXJk9$A58AuX!|^*nQ9@$aW$Lc-^j;z z@hkPAUt|85WFuDAX3|&M_r5p6YF7G#Nobm%7drGq+F%an0cpYDa4^=ghFKsO1v4ziBR^>I%%hs~>rWBRc&P@_sw%G# zBhaxnYgi?McztE@*Ci9a`mQYyl=*yhBC@+GvnJ~z1!0vC4ahr}T=_AoYnC(uVO<0a z=deOBgG{OI+Y(xpAIE?O9Dr!`)x$hTxhYFAVMrZK^sl@a(NA9eRA;->ro|`A#!UQg zm2!Q^n_bK5hh-GYce^~Tgy=4CY`=yRpv`VOZcRhIBXtEAZINq^l$OKQh1bL!44p=L zB1?=`2+UzyXKuhg#W> z_Yphz9rh-=5#Ygpq>e}sh~rPVX5|(*GS+MrbmN^ARDpIF|hzM$XY^z^3H^ddXy z`gb-w(y0IfV)37^QPl-o-L2o~PuJP&yZ%Yo6Jyc6BsXb&b~Y+%FplUZRKNKqKTPEP z1unxYcaEQD9-R6+*z-GV?N2;alR`E(!m4-|hacS;!yetI)< zp#FUP-zr4oLq)gnb{qd(9BBFyBt@|eD;R+ zDP%AwV2xNi0&RgLEZjYa)T&#FSQ4MOa5RA$ADlsHppFiFx|-|F{W;{_MM=M2Z+^<3 zCcPjze1-DDV-aexG2!O6QdQxGe5=xa8~Hf1v#mCwVMwtny{Vla^`6O|Iwr@huY-yn z&pRgC)bIs!;JtD8-F!7%UIlY^CRu+!St+WVtkeMNeSZ5KM6s|bR-^in#ibkeKDF}Z zSDYYAtfepM+pK z&ES?g0@<~4K2Q)4z7QDh5_AmxK3tdZo*T&y@QBhlMQlG#{f^T@y-wk%0?o1dG5zY| zym~Ea?PuqKAl@pa>k#5{n!oXpUJRKR2_lxZLtT&sv&E{Yk2|o-`#Xgwl39R(TC>p-$ZgH_=?nF@#^yaf*z+l7DItXy*Y(4#xR^gG{y96vJf< zP*033IJ@-v-hJgay`=&R-_p;C17Gd{?Q zsrq9p>h7Tf=~&D^oh;jCir*g=GZwO$>zR6hD1URAw3vmN$^}*&H+iAN?k`mT_w_Ih0fzV_TZ`!h?=d{!&DMd4HCR}t^zHsf8OfoQa z@XKqXlb^iaf{^ISlJMY9lxoiXW@qQUl(xK6x;8V7r>9n$kIZerX0(C3z**hh?HS6m z%QNJ>TLU(8gHK0ycc_ov9khUOPCrxhQv8(;QiypaHkf(G0tr&$2ia8Q++SK@GOYoR zZ&u|TK>;C_sL~t3_AmRLK`vx|@zVGryZ!5;C2VQq<9d1>w1k+wI%k|`=yVVWt@JoA zarV4#X0Nsc$!SC%(E{JvLL1>e^v23gP5ZCAf7nbNBlxkJ&X z;Xn-<5?>Zysj-=2LIByahNQ#(jR$@UO-N}Q+x3!H`3@xuHz62c25<4O>nVo{-HINx&CnYSP79U&FiNJ z?s>CjZfcIh#M_a)+Em9$_uYt|%O)?`mryAsvCH%~oKz)XH$>JYplZnBIuRruwf9Hg zqa%;kI!Xk`VfWpE7!VG{9o#8a4UYUa2vGbkId-f;Yt2ILd2e|s5zuor&iHR%Mk|N9#Y@c5jK>WdD!95 zA+14+TQ~0K8lN#iNh`u`yCV?p8GiX2Rdk);{S30rdPZGMC(wAJ;fQYQ#_kadbmBDa zTgYUwQjG{_+tmI2B8=hatXvv?u-;4u%SyG$S^iA%5UQd|=?4Mi$922)k0oKsx6v9B zgoy5{z}|C?hXi;jN1`zvVM5#GE}U zqMJWWBst9ARr$^`L*P82QE0?ZQM>n_mr7wW=r=%jA}Z~0^*<@N%AOR~Ef^>5>#sK+ z|36oq{@Ov7F7#fWOf8P^#_k#hc6U*PdRyM%lb?G`OCFttH!7@D+iWL5CmQrmIz}BL z3+x^#WyYi}mE<_b;cFO=8y4T($KrY@JMW8GB6E1zpjrq8civ$R;Qf8Ap{ zXMyF6o7?3&>mqupJFaU6svSX*hzZy5?7hZ6Q3W&pJqJY5PM%$E^u1FDB3Cdmtc`cmNf%=7h5@A2uy0898LF+;K-;o-^^s?P3t$8;30ARTTESz?(wd@IdJR&HreQ_1k@tjAmM)oKW0)?gHLblt0IBUzQY8L$`L zN)3?TKTh=j>7oYg@J2`+!Cd5SuYI=VV8q%@WvXD9vdMwX3(lPPQ=i!S8|D+*h`pjD z**p&_O^*F*T|6z>Zw*7sKWIA#S@jc+)N_v$!dQ-wImNqis3|8yZImy4RaA3D1F-Nu ziv?a>-0jWklHpG8O`fKJ4bCyk6PmA2NVVa#2(=zo-WjW;3pOTM7>;p|e;_d)-=K(T8G$meu{$QHnU4}wUiB$US!pI$5|x+s44q{ zdRYExkof3Xt*MDzA+7A0ZijAFxc9;DOti>=+gEu}#e3{LM_Ff{p1n#B z5-(}B1H82Voq95(TO%=2)-LE8NkjuF?5=fZ;qb9`d)cM)GufW>3T@vL!oA%@+P$qi zQRXsj?{Tjaqbeo!)s}HlPkZKOU>yqNMif27M`Zeet#*M~@!RNwv@~ZT>`ihSJxY|7 zH?(NSR>3PiZM>7fHKUg0+Lfh0w%#pjKva(0Qk!pv|Is|SlK`hi0FR&o%=5)o9Lg9b zbwC<5GH~5&lxR2-aFD+sd`p6U|7he1_n68EE*L4qq8kq?~4fmN_7OQ+{q!-U!v$; z?w1KIO->-a+gTV?BmB!pDvY33j+R?qtM;~gpYEfg{j?v+d6TBzYdx~GnM?nWwSEo%W z^g}M(Kr{DXXtP5HzYq!cafMl)1h#)p;1axj5^-slYM@D43-C*%_TeOawmSN6a#Xxq z@oQHycVt%@NxDSlw8mUQ)tA)xb^X;O6Wk$S?BPSH(1Fkq;fqFI3&5bp^uy~P$40@t zH!7`ld2PQ=*Q7pCa$7b5UUcK!K8p2thM4Wv3h^{?Xno86{UFFZPHRJI_32&fXVlei zz!lBlj1_vbqJ2w$_4o1ezXM``7CnU*hS?z1Yxd1kOjL3AmOIe`PD>UE4mR=SInqs@ zgkMzrKd#<15b8Dj8&5hqouiyqWUC`dNXTwXr3j%?kv-d#v5lQ!kcz^{zRwgP$1?VP zELp}nlWptL|BL7Uc3#YVfA9OcuKQX)AKksFweX5EA|L$L3jFd~ zI6d9CtMshRtQVDN7o5|yoqPDq9$Z$q(tx7D4MP$Gw&E>KqH<=9g}Yn&c=uRQ0i~BH zXoX#k?IO#IITs?B75N@%0gV}4_*jQ(1ke2RFL#zqe^EMQ1EQEwAs3Z%&>%)j$RrFv zigoQOuKDQw`sIj}U#sZ$=ZZI#<=c8l+?lz}cd~2TYuEe=KU&_zfkBnt`iGxJ z%6qFEY;|~kxOTg6c(1vmn>z-KO`Y-hmWg*Z`{&Au;Gmmlz-)(Z&Ks0F^sCrt8ZBu{ zN-XH}qI~OFvR#FWuUuyZLld(2{lI)fDxk)5-70E|F%eUwI~^{+^eCpiA`*ewQD3bJ zs6cJ!A_`e5Q>BGj&}M?Oda;>bk($$@%HO~qZ+1Knj;z2!Y8NypGeBbcGwC7X7P{{p zXfffS#c z!nH185W^U}?4@Kuw0ak2Nhb(cEyeW_Zgg{v_i!P8`BkpubV8B>7<=~w7>XNFU|2gjAI{aPWdS?x{Vm7i~HM`_vyGCer*QWAPO7@OX ztG3Oy<`DW?zvPPLC9|H9O(64+I^9z?rBXuvnoV7>P*+q@eNgTBsjk7za%gYmH2xYo z6XUt9=+WM!?x(v4ujiu}aLmcV&tlzhTKeU;KArJWBIWdjn|faA%khvD#*?ONTjKZq znMOrXVx*)oTn_sD+LNTJ*6V^9q!4pjLE9UXEg^DcXJeGQ5xb(k(Xop*Iw<$EMD&B?uY*20v z_1ew#R*6lI`&+Do;&vvb5CMFAW&@&->8x&3WxfsP&zmyxdWypp!}&t|3~(0YIk)H> zT&xm|-x?>cy)oS}6bU1%3A1DZvmV~6St#qc-pME(xL+|ReSw@FS@F)% za8=P)~w+p;z7rY$8Hsd| zH64A+3hYV$f44Pk-WdK1mYfa{<>M)Hg5w5ckHV!|7k@M|*SwDHKXAlh<@ zNBQ*tQ*79`_i{6mAc35i19w;3D9^&?N&}uiAPxz1ZW^V$1c6N~f*VKNz6Uo3@>P=> zO~w2DGXyrH!tUF6zDF&j_sVuoY+uy5omg7lX|rNtdGy*!LAeRF5fG+<0A&OmynSp$ zm)`Ef`HLg%t(BsCK21p&sMQrLLHc&}ckvM-t#vgQjWx{?=3Y%eS%?3LeVb!zHVm)Fp`XFq{7tg1Xd8 zu856yhthiY95vSYxj!^L;0la~-bxUF#<;eh_^673MTfiTC9jh+_wDx##7G~Vq5>9!~WDh8$dLxfKC zW@+|a&fu)oDXLb0>UFp8xvXZbz*H`x= z#I&b)$=|$%^M83XQjPSZ_Os=4<)8dS7vyO+no{ z@xsG4&bDTQLax(bM%389e6qK)Fw2wz`&IG9y5{eCh~_OJ@YRhvv@_%vVM7$Gbfimw z%L;GouN($pI>=f5WmBKT?J>5^1wSo=+)&^vKm@$nM$+edhTdACZMN-*&j>QrN zzGKWzHx>afwWqm@{f?9n$U#PhM^Z_3YdPPhvtLJ@ zKHMFf7hh!YCR*o3uEM)44YMaLER?_%{n9i=p~E=c&WUo<;oeIYuW_irdu|NBYhvy@ zi~I6zqEF-?nGyN4W*TGqU{MSI%GFvMD_A-@CAVigwFW5`g#n`}=r*ynM@!Jb;PtZB zT^CX3fQ5qw!APGV6`{1d3|@@3lYFSpQIrAmr<3<2C|azmB+WxAIEc z_v{_Br`^GzQCPnW`HIKS#pUl}A2~I|a)I|(V`s|Us%igBS#kPPBK@n=qxMJ&L@yio zB6#@MS z$HL6zmkE1OZ+5BSk51)td z?br41NO5+1ZSL1BKRH_Rlsk2f2{q^I`k#W9%_}LY6JLt@6~$B%nqq^_FP$+`aLV8! zNwO_=|F!^ER0S?roR&v6jmokewT=)L!TB;5QZ(&V?J!%Xp`zEy=0l9mf~izn&NHXH5Lz* zx}`pr(`%+s9e;)mfW<7T3f87`t(HH>nc542P+kv@-1_JxbwY$B`cY#5CK^>{WFKl= zN>z1e$CDP&)z})1ypn4b65Lo`hv^*!9~8>kf1VI9unXT!?*<e8RvbyjSj}3CS$ICkStPOY2714ig#vLKj=;N8xi1#yD05y!CGcvX4>`cJ0{lw zU-fRVzLSMh3+wU+Mcj9dWm#G2-ESHje_)rs)JKMeUql=DrHnGWH2K#5U6B#lH4= zlGD6}$g41E4>I4b|o+CqJoy{@OV&i$K~Id_^f zbXpTtJArh{a5`kXR?sTK<$DEBavU{#&39fcZVY$BVE;0qWb%Q|uMT5Z88Ka@%?BZU zy1aeNu4P4Bcgq{zC#2!+KlfUGf^Gvi8xNY^Gu$#am5yKXufH$Kyij}GwR73-e|GRv zJr!GBz`Vutqs-578c!q;F7j2PB1YcqpAhYK>25j4jmRr^N*-g^Q>TKh>U2KM$1VAC?muXDUJ6#tOer%cr;pGOIuGO;UAY(cQbp^V(-d>DQ@maTc- zU=2?0mIPCyAWCJyhU}Snn~O5zwTRu{vazWs0^N7c{LF`XN;tb`Y2BZI^TJbSz&IaY zaEKyIxs89Hd7n6EK9s#-zOPjFOhwTDGB*NY=hn2uhaHRct&%^y20bleH@s7lc`^z} zHQ9dR!puobTN?AkhXwQQaT=?Q9o7YO1!;(#(XM4LS}g?y6CvCSE(UC*lYI-plITIA zszDHZJFmnrf1p%0NH*E@1pn%Y+nb{{A;H}a{kw;>kfZ%wsYB8()Fnj7QBu%f3fRK@ zI9&;BpZ&RP0eo~Y zH23BNQyOIUC1w2PL!+b}4(C_@N|uW*%a0_71@&I7^?)l3hYN0QxBV@DeAd<3=uMjSaZW|1E)T6Q)Kp))J<87bEZePhtw5vDd`RFjaU9e_s ziQow~iYiLNOUU)*R;Avq-@<4i0^T}V+uyF`Teo7f$2kr2s4P=97{!H2LRG*t-|eLG zHl}ceRPVi#Vd!e1e-e8hLUSRGZz*(%@u7v9ix!rx>c2CRTHN%%7CeCm8+J0(x-oE2S?x>S$S|B!@BxB+sUO-O>K> zug{k3pVbN|7Ukj+(Ftjd}S2L3f9L_$W1Nmg^8Bm*u z=l5MseiJXE^_?JD7e8iDe@6M|y!ws>U0(p1Dqaof=hu$AME$6&t#tHJ13_S4^6~f* z(^b*qR@IOaE3tKfEM z^9T3^^ZnA1*;Ib0a<2FLSi`WP8-o7Fej5jMbi}3#;*do_%k_>`;bqlrc-xM_y+82} zBhC!S7H6rQRPLz#yDg2COS2c)4Dom_TJe${AhH-8^DYNnH@xoH!hES5NNlH=1N}yf z_l!S-Ak%tr1{84`2qM1kZreg=J4@M>!3#^uQ`IrZ#ap}{9-j`brMY_YXniQ(Px!h_ z_Imvoum97Sr^NoNHzEiCmnRR_-N#hw@GM`UaaMZs}@!;SOr=sQgxyza; z_t9@U?xGlw(J-2O5K&Wp+h0MAsY@`gl8xzn`==&W7Y45gs>X3~Xzf1K#^75N_rp`g`14?T z!dv{zV4Aiu>>yJSlj!rtoD+4(&GUhq_VAlx-vi#g@w4kiNuQMEzMW|c#Sc6i zs#pM~jtrbfd^#ctcci>VNAV@740Q7nz+bWsH8W@mxh%2?z`BH>h*Q~XYd@PKG?=*v zDV-psZicWVx9;stCMo1dGZqAyi}Y6#L-1t=>_X|IeO|cNCJk$5U;}le=@h?Mc7X`ee#Q6pI7wl)};kXy~MZz zp^$~ZtneeJHTxDdjcqNn`(#yoeM0hW%R1D>K#57o_RoB^1j)M@+5U5dXj$No*#dh( z_ZizMx?O{;xy)qIo0-_$H_LkG{$7zEh$?g5Hh{%io(w0021D3^JlwA&31hgvc`cp; z-YRzB$971^%VAB*1q|cX#Lx4iNI>I1QQk;jQb!tlL6ly9lDEi2eGum~g>#JvfoWa- zA<6T-H0~Ih^K9+;wik62_Zua4?C+`v0yG?^X(%$~>xyhyO$gm%a@AVH<0x z7op_+zw)K03za?EQ#E+kPe)yp=HPqp61yAWDP1(kS8wLuNoZ-dt`^_sxZT(HF>2Qk z$aZ2>!G-ghO5YY(=Y%+OdxORJ%9HkY(TLTe&4NO7&YwxfMh56q&yQts=!sWLnZ10011S00$bL;# zsC*0T1F$?DJfJ%Jw%;n&s2*g`sBI>%x?^aS zA^hCK=a{WoN`LTl=hqzI#=#AmUfUsw8zw~=&!jlvUn~IX$}c+IE$?|f+_dCggO6b5 zo%#|eG5p)i6gSH#`;}jnmx3Pzch1JkdW|`+ukC>?-f6>HJbn9iz+&9P-6BfrViuoN z1()r=GNMiP%OmEpNK`W~E<$JrOyG=H}1H2ge@4U;=o7(}cq|Ap0#dIg# zx4?rD$)zQD%`xj9UfgT>O_Q%q&S}JESae`jz>YKB;Up^g%wcSzZ9-o%Z=%5F<+^$; z1Aot3g8s+3HBFr%pCL&^sJfEqCYB;}GIPJGO#J%Rmu735@~6_TBT2sZC3_;t(XjU8 zZn29jjn0WiZ+*lRq=0r9=%1$}e7>yerG;FFmK`1@JSP|2m=I=G<5#ws7l;Rno-Z6V z$enH3^qIV$gk<39m<0T~PbZhmDG4WyDqqWdGV11oE-T(z(Glr1)hpF=GT_nY^^d7^ zu65Was?ir3r@OVO<5}w|_9kM**n1J`H{mx4c0*e>)dxxqqt@2{B3M)JPRh(WPKlSP zndk^tnZ@ymwM7NFx>DU^Vk1dq?Nv&M(N%TC?IK1TyE{~%w%iRGt?@A(Z?aTFGB<@spU`>i-V5q`bhDz4)GbNu9vP;|vkTEL z)&IZSw6sD?Xe#!7lB8@H^&!LeSB*5?V+~zpvLrSx!fTQh*OToYG~1G2T^y~6W@?c% z`Po=c*F%e~gc$PM4sQ@!mZRjqS+f|uFFL2|Wm=KQ3y>&Z>&pww#VZVbb({QHUAC>fW;NK^+ru{S_ zIQ=6pg;}L}h%MfwM0=yril0GvgcN7ej1eW7bjkx`%iR*)QJ+$a3KG&aCJR^0sF`i$ zueOjtAQNrFBq_)$b4{W|tjpf?e8qe8RQ2AZeT8H!2sTuWzgrYrza$_%{LdEB{Cf4m z(B!sa(DIqUxTbGU#OhbPtWtQ4L*cj4wb9_$4KM03*x)%oVN@)k z&<0v!I`mr9s#9sa(6nMyQ|pWD^o$ZEfAr&C{CKDB6(UFA$Lfv7kHe&A5k9nFY2-J8 zoprHh!Fj5cCJ;TIvN_vKF04O4i2jlg_Aj@(u6op=$@3ukLLZ+g$&`^~u!dj2Xjgh$ zTmCgNjMZuP`dykf!`&rlWKnbdn^fUb?$IfS>H>kF8?Dzf%*j-TX;-s^;PzN8f9JL3 zq-)i1hi=+nrsoA~jz1D9P7PlxuyHOiF^3qzK6O~?v8@w{a7yHd^(yHVHD8v+B(kMc zc`1?&574|Ri7Aua54eD|2FIzyVBCf`d4o0rl^crjXOHp%mZ9 z*iKIGzI%}O4xRf9tnd+uZWy%ezlGoT=Tp923&t&PujY^jKPg9-jc{b)hyyL%h^ia( z7xDA=*b*Uovfc4={oN^R4c3-JwIiEDgIvsbOfp~b*}v<)H&;^YZmo3s+y-XL*;z1| zGdi$bs5Njd`0oU9xF=lFQ?NdVwLDlkLGqckkVQ3^@3Dg8n^LC?$BE^ncyXB9g7wfnXmMQLM2?3LUN~^5d+1=qf2FluYZg$ z9t7}~O2+5ZJOkOpjLTb!?63N>O#cAq+a>P;F(hI?XoCI+HQ%-=x`x-9aj{MIM6RL2 zrXA0Y(sIZl3!=g;9+Va6ckYNb`@a$A!a+^LAZVTG^OjQs^KH64OUNPgj$N%{O58mY z9~PnTRRc%MeVHDZ@_Wvl{~{hmGR>`%(?N4FU#~T~b?5|jcw1f`MH!A&Zgw>dm&$g; zgWaFvM(xCq^FahNI0N9(%DY*cTPRlgmpIh(YXFGQk(Cafw4zbwDd4TevVm5Gw^=m)uY5ybV zw-CJ#Bjo-^dKSzA2yyo{f7$@W$oMcIA6d~nOv-rLePoUH8^<-&zULPCU1Iaky5`4e z%ym*IP2s|$_NDq+^=t^qM~OFGZu;zU3i01v)lOx6g7{iIG-I-4g}Lx_OEO4es_$Vl zRl6FtewbZM)Hl@nCEF1r|Le7=1dC&? z%3TXY3*0JW)SJ= z5N7%w)cNme#q!RY-ebwwYq5ZtOZ;h|{h)AJ{Fa-gJqQd2&>m zVWq;p*Ze)a#KC>iF5KNk98#3uh~nBk5deC}INSQ!lJ>gcp{DL6D8_v}L^>MXn{&1^ zTXXB(;+K`7g5Q~^%`-fz%v@vnI(KC4`OSxoH)0!mOD5``8MR+R_rOm3mB$eVt!J61 z6}H$PGh@332*XyCPxbuCF5B`gO| zn=E%I#%vfPpnQ!blg~_zb>=*9MW01DJdnlan4{1&pwVjH)9!c zi+eb9fw!37y)mPD$#jb>U>71AK@BOEEt@X!!@DAvQ|5axOHe<{{=5$-g}PnHw`cbi zdS|9SjT1OBL9WHsngteV4P#V1Hw^Zl;4R>laB~Q4aFTn!ccsH!k7)gXX{S1*we-!b zFKPZKF1jvxOU?ca(p(dD*3TEcDG^{A%$n)PQ!?mcPr=r)IbNrT5}DQ{tg-}pY`04m zKBsCxxA3y{%pztI@y@xLsnYc0c%{_pvCJf37JM$Zefx!KYiQ$5DoWM%A(dX$kz5;Z z>_)&-uWH~NImy{UUDb}4yIZ_@_iWO}j$8X(l+(zC_tI|Dq^U0b^g=$rT+7U;LG=(T zyT4D4{lz`yPFH*vv;94&oFPA3s?DB3!aZ@l`kxM1pMF3~{vU2?d_UsC8M z9W)ChPb-&*xMAR*jqI$f?L3H=0r5}K;yICP^7&`wksqaqv8s(ZteNSvV0(Oj=qNoh z%^x6tJ8k5W-LBhxM66hs-)z1%{jpYyvYT#dKRrraZktXY8~f6Gc@lsElN)ZeVm-5| zSda7aQ(6MNMk}!+EXYPgjqIc7L3D?aB`;5hK9dJ-SF|mZ#Eq4^DGMqdAwuA1-@2Al zvYV(Y=I6zh4G64HP$jVm$(4F7tvG7U5@dmuI@x8W+uAuJPOdpy@^c05_o6R)TbSpE zP50T1Q57?|NAkLZO`!EG|2P@YU+h8}GmMgVD06x$UU9DJolWi2jYQ|2zj)g>)t%fN zR;I!@Hyh=5U%+UQic9C`&s^t0c5{8vtj}86e^ay~2kZ|CvE2{}Jxu1;a*=?z1MJ~i z%PIbAue9?eAe0?LIUC)L7ex^{2PF(LyoRXEAOUBETn%WyB)>Z=GwxvxRV$w%H{&91 zXO3&)6v@IX7E`GH-Z*~@P_C=X=ZS3zYJa7ATo)7DHJx!Q0(r}gOIXY)Sdw*^%ieVrI=y>!(`%`L^`2xCWv~S*JVZ7Q}jx zHfp{9YNs>ke&5fTKUA^*#A=iqD2icWY$Orn6LaK5xKt5F63irWQ6m zZ@XuhJXYP$yrzK&xh_5gT{t%(<_ex?depA6cLfJ(r(}HbR$p~r=xAH{ZouVX59*{G zGK+8bkylwje+m%I3&S&OJICxpbE}P;tcO}@FHec;DRsin6R>uJ<~v#HD16U0F`5oA^#`5AUmiaW=s%Jb9pR z65O#4S7s#X2PvI&JS>Z7{^@Sbi;&JG;gSyH!#EYz+v2k)v(g{9JFOhtBmk{c{va7B zI)4uCXPwB>+&>(4p70?$MY_g{)OW?l&~kS0&hkco51r$2i8EkZ5Dolu#O8Gc86^T5 zGt%>p(UjY>jQz}_lJ538lcS!cGM_{AM2RD#Kn9hV0`|wN9BNm)$Ms{s7cm zV?Q1`hoADI^{i{iO%y)U)-A7L1rhqiFfINEoq{@=~RM}nIg+Gu}V zXVyiE(*K;va&oW~BBC$zrsn=tlBIKtIkOX^bDHjcLzckl9}&BQ#B6t^AgC6DhG z*E{vES$Sz8_euwX4CE%>u2)~yg%in*-tdf|Wi>hheB_<>PlZ!kU`oZucu(~nQ}I5r z%xZ6f`La~>dDV)zyNz(1a&(Y<6ZBbOGL$!{n$46n zK>rl`=kkxw6 z^uk5U_bFUmcL68y3-1Mi%?Nm7%E8Rd)I<7~>0rRVNZUTId z4o}#ohA_sU-+Xs}6%^Mjl|9Vv9ex%$@7-xgr729d-oO|F&Rg=Ti zr3Ayf(|ns!6N>tC+CaV7-XM(k(JDEDWervW9n99`;WXJ&`>T=x`lxcAavq)kxTsr& z4gly9dRnZwqgi}k>Z8e(bCm6O1046Wna8hXT?ZDuz*YM( zjF0T3rPy%%&zrIbf;1@z<7Db1q9ZAsKpwN{VYF(7Hu%o#5a|u#0wdD9G-11EyWlv;Yj9_1>#?^ig8FBAY zFamRlic*h;9=^h`q{kT(Uu&+}^S+iEe}K)8 z0a5Z7lUbR)q0W$gK`s7%%Nfp;Rrj0~3^S7j05n7Lb)iA*YnJGDs(+wxzuc_|doi0* zEkkdAt9)-LNdq$Nsu-6*s$|z&yb`lvEepo{6hybXu`bj|64HEvV!LvDj_M96d|ITP zoa;i=Xb{bnYJ*s8-dRx#Y~M0wK=;ZDSkD{BO0ODOGUa`EB%5@=c3x9TB1`Ch?eTF5 z`3q?Q*SK@!TY|4!LaD9TfT@dqa@3Jsbg#_hUEUz&O{Sza`W^W6_zRek!4;9;qfpOz z_p)*hCoLAl)m|MkgJ-N2a}Tgc{WUc|Fn5#gsvLIm8=t3GFA=v~6gOr-EcMS!HXj{Z7UrsJ3JTD-A)96Fzp!II$4xjoD5OV@CfJ8YPzc3_EG%`nFJ%2V@$ z_X{SQ9mcK5mi=^6pXj}p{}b#_1M7=J&=DhUXhro%W~1jaE3Y-WuB2mKX^RN3r>D6j z)Ku}J3NhOyA?F{^yqrP%ZEAN@%GMVRm>({$flo>wk`JA=4oJ=gD~b1a_mI98E}L_z zX@zqV_%s^Pw$?&jPDzW}nMzLxlT>H{n zC}7ed>1aO$T|c`ARJk~sn1D^d90z@ds^KYDv~F{q_XlxV`OVZbZgBIoTfF*(SnBC~ z{KyfkQW2~rK0K(KJ8JGiP;2oTd#8Y^IE;>2H+WUo6bK4*%6!IBrzUBw^MPPlLEkaR zr>iHC^7e>!d!3A;ic@`eXA7#93SbTU4wGLpd9Iea=GA`MW4+CNDSVP*^qb%OSNvT| zKXar~_WJMl1`X91WgYD|9|=7j!Fm77zhJa@qA?<+Vs-kQmoVLORoqIl0P1XA8ttm@ zx%#fjEOs}3^{fwR%>ZkM#+HqpWKAay8RArPEcC_7>!c+mG?Jgv%n;sqgi@+#%wiyb z-UvS7O=Xygt8ST)JrAr|V1YM^#J9(e2i9Czt@Mo1`qZjA+3vjS#J}*r@T?ey>+`_@ z#PYavj95}@oOMr^3-V1J!uz=CyxNH5i_W#$!L^#iXmLcH)(t{G!4f@XHi2^WEPoZE zEGfO)Aa3J09j9%j@2PN8Bz8lpp;BxTm^VH9r)^0Ud*rFXzV~5L2|;Nmz*)ogz57vnU}U4LM=X53~Ky7+toxAEBa!trRE z<($!1kewV}EAeNJiRcEzkXB%j4u9oHQ($I%4nxzZ%3g4?J(E%8Sw64oxA0$)v}YxF zr<-N$SV2guuDn>6aH{W-Em7bl)?mp+mZYiVCd)Bc^nWNVF;bhwBaF0=e0&{Y9_z^Z zpUZ;B2z{tM)I7<1u-(!KNqo_i1bu9)dAt_&87hqT0$cUI`l>MYu{;iHH))37ErK9k z@*_VWzH}SwpINg*IsaXEBj0+t&3>WE#7JvnzCPSnE%@8{#Vg>1u~_B#O#AGj2cr9; zF@=@}j$+;3oQ6LW+wGzM!tj-%n7knT$jmma00u3Wj(ZUjaYJy#_3+=#-?|UAa&n6NXS<7>U6@K33+$7 zU!w`MlNaczd0JVR$-9&AUwz#vV0B9X@~F<>{`J-G#>1~Q@L9X3qP+(iF$GX@t{A2n z^P^W2He7%ja5~?8?2DIv3QGNHtR)_M*T+uKQ>6iP@EtUy+&B+hRp!-y$h(hiYB+{>DK;yCPs6j2@JgLw0h*NHA$dI` zA8o`wIr6pl60dRL0v(v&e*E=jm`btNz{}8`dW5i%@xdI~zIZD8eQ2=h{chL0$JWTZ zAE+2`W1o-W3%`iWo)0e9OS+461KeDyUyJp0ExHicW5#RtBKz%`UNZUO6BVcTlY5-t zlHRdvx-ougb3d6aajK z5S*0!A0c`d>h|9kQeYPJwf3jSB^Z||#`K!Z+jR;%7do;xCVAVvI9LTkb9P3m{_J;yQ%T}-4Ym;VTWEFp^N@NRlcq4jq~^78ez>&4tu*;Ydo z9yl(Pqu1=rybo_Cmt~$n<-3?7itQBvYl_cCDmocyUwvP;>Mr+}P1$8zjG1 zmQYcY&Gjln=f0f49!EqqXkTQpM>1EEol^OU{DjmcQN0ZaNQ*odUb1psE#G@;itSOv zIU%b|&J?{?3gTk#6>~1#6bs;cp?b+v@;qROm+##j&AK>>s!b+Otq6y0{7dr!x0emc zhdo_NVpGkaHD+<4oyW*q(OG3jX-x&yPpzE|^H6;`nH$AFOrYl$Z^v#?QDWVWAQw5q z_3(qo&CghFkmK%G{{c|mx22E`3INteRq?R#b*y&%AJ!?RK0bZ2wdgw}dw@T+@qK8C z#_fRID#?lJ^kD4$=cG%w+dNkz;qD6Y2`3oTrrrh6|dA&o7APSwR}s$h4uhZRaR(807zxW{A2sy++)Ai z-kQq_>7<+&j4i=1_5pCGS`4(QJZ(b;-yum{-_!kRV1IKJWruOyIi2_1BC#x3&!Sov zz1Q6>Cq+51vsg6wkAJ_qt2;7qzVXgm8-;LT(Ucxo1N=mn5Xi(d9#T1M87)0ER#X%g(Vz&cMz0vHO z-{14gyW5;XRrkfni$#_6-_6IL_^0IZ5B0xRlaG=55y3HC(x{mLKtf3QWFfM6^*2+@10G=$3f$ewe7tB)XLYZy zwx<`)5cIUK%smn~Y5&Hj=;5#so2mbMS@h+6v5|ZkVGm0qorOQEQF=J4S_w`8a7In!^+ zx0h&X)3JXD1{frr2%1*=p{Q$>xY)k z0L%s4ZbGE>0tk63^_Yt1QqK8rsJ=2pV807ZX-mg#X0fq*V@2iJ?{B3xSj$3C#+OZ1 z5$bq2_2iC_Cl-k-8!|mEZ`8LejTXe;>4f}BQU7lyTANsf_gN>}gJek|Wap?7lwL?p z6s+dhvU1}g&!up;u){IZU&~IFCxNdOc-N8uCA~WK)D)8hIc4pRS?}klqA#NO6LvLJqNVt(PKtyjziF%($Lp&0ZS=Fs5R# zUZmR9#|PbiyxCKo`}QKK^8Wsa+_BwM%@_HwdAT>n0j=Z`08CzZ2yXY!>^_V!s^6#Ld!{8Gm^()7CN`+Qd7-6HwDmIRW_qpeHjR|YIA`O(<$|?v| z$W-Jl)D|Rv)h8$w8|dNZT^s?dta|147W=$UxRkXRpYX&F^I&?r)_9jploHqZACK>j z{`J|~aoT!ZneCjY$*5-Bggx&u6V!Coun(!`paT;EjFG}{WhV!p>&%!!#Y(rzk%wG% zs6TJBWlC%ilJ zv6DJ80!OZ+$9iVpPC}-S3RK87{JJ839~iL*o^|e)i(SuejXvbcM*Cl6xVBo2f4Bqz zaO1;6OThmOK}|>|Y9+m6J9Mw;8qQ`m)!{YFibF8PWS006h-49(chK`@B6<>|7y=*JnK~)hgyKskUC9>Ht1ngyIT* zW~c%j%mJlI(LQkT{eK~-=!C!Wy!r3>Mp6JD)K?76J+E`z z1S3fhrkz&f@LE97;?BKUa@KVjrpAO1SpVn{2T|8=a$2TKF^Avi8Zu08-)VsDZMhoul-N3) zB82IVpsr2a4R#=9YGbz7^?Y-dVaOTH;v~6B#NepkkB1u0?Y$mFwsYMw=_JAJ-axVDu8`)06FhV{+Uv}~4P`z}?zbWIyavaNyn zLyF~w8aDr-@c!w;?U+Jrq1;_#J!n8&)-tN1*u+B3%t%oORL z^hTfce0KG3QsG%F)Z7Ky`j3(PuTF}F2uXD+bW$0~{Yv^P*7Y`DA7WUkL?`@US1F}Q z2!{vU&0hl~@m&vy__Ma){EOH!wU;rP+vD#0HL(J8q@1K|8~h3L8d>}cj(|W;WVWOdR@0p%ugmtN}cY zNs>cNAE8iZPNRB#kNJ3tDQlRD!ycFaz;|xV|WzGA~KK?pc zHVSne;vCBUzMqdz#ll$yUq1)#8Qj&*x2QaodKXCUM3`eDV462Z0{gztCHQ&Ibf5^dXW7JNVB_ugd(bwfvz|v<4%V{a(2hVqT(E z1S=v}iWYusda8r8aiMH$aydr6lqJhLf3s%_tE>)e1A2IQkQn;&w(%b?hV4813%+H@ zOaKuucUbM2DoZ|c$F)AFn)S&bwDxxw?M3Ic`p#bG>R|PXfv|E5398Vy%lBO}hmg#x zQ(GzLQcVl{C9R}wH(hq09{s4n7nted%4r!CA;IWb_-gz%KS$Zz2wPxk%mJ|9$b&@P zt^OzZJSVC@Y=&fLgTbmHHw#2n5zkWc1^g`m+XOk)?(-u_8Me_;J6Skc{1N`Z&EZsr zY4re{Sa*A9?{R#j10IG@f=;spqR3Y$ zusTPaY{AKGvM9OotK2iK91uKcJnxTD*o{S{u<4xJUq!nh`aNHicZSjS_|A6RzTeb|?$ z9a*ABdXid)stax@va}gXFymxluj@{JtOW?es#}gmu(dn9c4oEZK+j6E8W3BJQy!ex z))~i@-Flo*@CkGH`pniav4LM^ImA>)dVd8se4g>rzT{uO3%K>Y2%ed0T>LsuFx*Amf2OM|A(h@k7xRi{y0)nsZhC7Nysha z&ge$&A@^&R2y>hJ%;nov*xc_-A@|8`?w68HbI)a%sD{a0GUl?m{62lZzu)hFd-#0b zpZDdwU$5smr#ZdfY_d$4*wb(n(ra+m3hjTmK5oe$DJkC)743#VrE}{%8KSm3Eclst zWG(CLRg(Ae_s1#kJ3Xq?>*6EA(2&`|WclODzV4rPRdb}eX9a^+c**XondJwUDoqkj zddJL*$0tnMRRV+Mz=}zbXQXZ+#BAIPS8VeIof|Z!VU}ls-aDpcd4T-@o=ciE%NQpY zWBJIIQaN8QkiJ})n~kM443YR9>?0vz8xZbVqlGQh^htPV{x8%h4`lDvMckdN^g1~T zlq?-(f7a^BvGnX)YYjq(-$2Kvv{5Q=%b(Ez!9E{T9q<|dmz(*qHTz#CwN~T7dkGZ- zq<{5W7ZwQn(&&}=dWdc#TfBnT*Y+x#yUg@2Hm|2h>?SbU4fi52Z9Sf@{Z#=QyGPr6 zg_Y6Oavj>%TOIqS9kk3)RIohE0xy3CX<_GWLZIVtQ(p|RBR8@`DfCqx=KE-tx?kJO zUVeNUrT3+>E=aXZ;f&3ynt1_hn(JipUjd;H zIF(kiX17}DW`L;x1mzVHoNxZr#A%f}OBx^69!yH*n7*!bo#BD*_(rY6i3OOeH%MCu zf7A=|R+D`KT~YlIFz~ypWBL-VJq01Mx>Oa1Yt_?%PidsPYR=U?LC+olMRX%7tNwem z<18&-=bByO1eE7}Oj~YvsHUHsY}T4OY*Xt2yOkK(A{8(r;lku)Z0%vA@^D(XXo$M0 zOHZ?Z{+%P?H;xeoyI#2^x5nvQVTqr@Qj*EFbF58IuZn{OMPn54^#kipQAd3pPY}>+ zy)JwUp7v=)3oG340p&53iHbHW+}M>8aE1TQeEQO7;bdnSDO?Ry#G7oe{RS47mBY-VADK)Y)Z$V{;+D7|N6ay zQoZ#Jl;99>Tik3(`ENAsS?D`V7~`F9bHNY!hJ%Tk#2W_%TOMBc*+OM{$a(Y1`l#~5 z+nd_2^Kau!-QX@PDkpDFsPTlIO>!p~Z5KNZqNfuY)jKz8C!+ z=!zkO8v$VkNC%=y21WSfSxF!7lyZw;Z(MERwl`VRGnAr$FqB-`au!ThxBOAX)Ss>WLH^4Gj^lI@C{rTC6otmKrI(6_-bB1J+HZF*R3#b zV{MCYNZIb*a5CpB1FxDkZ)pwMVmn~q=CyOogJnyt+A?_LBP-9&m)Gd*dv9WEQb#)* z2jj{7hh=3N2c`6+JS_9F|ExxtQ#&p)ysf(8&**e$#DC8$#=fKr+kS7t0E+O(M>|Cn zhf5G19h`J4Cjm**x1%V3(wO$uazvC)blO=&eozS5rbM_gsGJ>`_*_Zk9$^9`nvv>j z^Ey|*uRfPnrn&to%u{pAc)TVGZR*+InBn+1{g5Vy5>!bBDL?NxbL> zF$K@Yn4{j`5>cNCxn??Ap`@C%_ce^W zxkm;4D@@XGRd0$n<-h`jtkD&+fOOg^PVj9N5A?!*n>r zkv75Nuceop+;MsGR=+DMd!;kL1Es%9xL%FDq*&^sETE{e!o6f4+JO_B5JaT$zE7p+c`kXv>%jN}l~OX0NL z5Dcf2az{;g%QB2t1rb8HzB4C&ey?KP z5V2D{S6zJ8H~zHq@mZvZ;lDF*f!(jL%~OIQsy*#(!NcJjQ5%@1#E{fWe})E)P)88% zBhXtu<>%(YCBld?5Kisgv1}N9ktIUi@G9GUND5YA{%5Y~O5w_+@&tLX0fm>_XvY_Q zPYuCz;`vuqr*#FA3QZv6>N473ZnzMxc-5$bPM^izIol4cdzd$w zW1pvwe>BFuTxks+8~$yqa5aj)<2~y-SPndLTo_t#_7>Q?;_s?Hb5BQQqrx+9%DX`Z zS&d*UMlP%Pjl~HDxPcVJm*33%=8sW41D+Y1M~8c8I@Ng&I^k8^V*);J@L5f(h<6q_ zEaZ;AyLV^-NsX%KB+Ea`vB7OV(w-;|z!y1(2@`7t)reSS#oP*^Dq5AAm5G(?E6=G* z3qvamJEQPEv$;7n0CHG1x8FWPJ#>2KK0K|S{qL!AOZ>cC{+g!Nw}IgGt7pI-L=CVq z+pSMbcrM%%l$6NGwFLHM?f3bG%SD?_fuLuo?0);NtDl5@t#yuCHB)bo2?X)4hdZdt zlGWx;#DY`vhw61xr(XpBxI^qsC711eQlBro7xiJ;=9MztEK3#JROv1<&J9_krvz9J z{yt4@2UiTX=jLXUdgOV}_2$vJUCw;}>cx;<{)WPDiXI-jsC@o^qmAo-^naj$JOfP= zA#VX{Bj9DLSf)u^gq3yJRkO6+mrMu|J|b%$?Y+u;9AK3ArCAVKnVL6<8;j6Vp%Nx; zv_3*rbLFd5@V%BiAsyFdd$t0$@fNc!$Y5Vzlb@{sVDL0*!nX0$MR*m&9Bu*jKsAa; z)oD~dI-8YFUmk1{?AzO9HJA6(t+?cWKFCuZ%O@c6td(yS=*S0{*)A+7);~`|S|2c| z^HIbC5b9*%+RCW&7FirLJ|Ujdtc=3%{Rm?HmBUtj>`(Nk+}FKc;*}z)NK(E&A-fJt zi0|%Y(k!qy5vQuRt~t5&ZIlu$)#|n^V1d zn;UU^BZTyvE3Lt?J|fKCkRq{&=Ut&h0O6Fbpa=wB+2&QZC~%{`026^XUho9!cg9yK zpto>MH;$1Xz-sFcoLc9S*N!SkJv4g9*2cBc_N3|wp=ZL3#c7MDKy&ylV#VD9eB_AI z)qv<#8E*Lby>b&??GheF1)>=wL9a;iCGIGeDpVoJLC{(MhA5#E zsI4{l_3>U7W(6v4Eim`-dPTdWj~z+_)E%p!noj2lD!V|;<&zWhomZefp1o_0NF$jF z)e9kQ0qE)NLjM#}t?2;Xu6AZi$bOC+7V5Tu9Dyg{pH=DUrivOUG#tib+ zvJ&Nf+w;x0S?gp6%Vv9#|3|{vpf`8@?gLQ{SX5}^gY#hg&Xv&jQN1s2tPhA{)N)cp z=}B5t(L;}FC=c0v5nRu4kaRQH7^U#xO6zg$Uag!*;(mvFn`B$g@6#1%vFG`*Fy-?D<7O`Re{Al~j0o1YdO0_ck zQe#O*zm>zF@;P+odv}>T_>^>)Z zQ?~OnI-mpbhv7|}H}8?9=(r52h|X8Fk_ zt1#h;JT;B_0m}nR-Liuvi)oCwBGL`8Z40~z-t1Kd$bgL zTf!}9?cQ&~xA;}uHc4(d8~^{o^j#S5fl?p-#G>H2H#fq{+4v6oaZLUoBd#{n=EU9&e@X*wf|y;}%G#^1gDyYX@`9{aDl7ItkIa+@Fv!cUE4dH>wQi=zCpb{C`{9x+ciUqAEe${>;WSj^5U4QIi z*jPqTjRQ{VS>Y-z6W&lslb>o3#!-z6*q8^eEW)6)a!{@!F9{I=G)O21)rc&l&DZ!+ zpl*cl%517=Uv{mbo216l=Jpd7UnCRJ_wR$;x$HAGbF26(?h0VIXo1>MckaU@A`wgn z{09joU8Mp__uGX+Xqr7X58A5=W(-Nssijwj3>kIK$cQ1(L~+wY6--QmJ2!pn$sO&> zb9KKfg?xNjOu2@eeB@D&4)fq_aIdddebZipCx=WcyTjJQOohUch9GO`&4_<8mQO#O zFIn)_0(EwGq@}JZ|Gj|lod0oxZ_`6vB@olSXXZ*Q3`YB!Mt2Oi9h3&?KFl-xg&r)W z7+xXt-WsHFaitbRB*pOrWm6%G#?;y2*o9tbS&?!1dVc>SMyp%7`+Cj#^4j6`KNpzu zR6>Ft+J+Fin}U1P^2*7(A*&L75VyG z6kTyc5i_P%z3Le+f_qC=pLMh^Tv5rD>QtG8Jt1v0CWx}5|C%Dq+%|8SRPcpyXE1jo zQ|lr64AC%jlnQe~Nk~emd%J?MP)2-_bD}E@AHNVjsN9P(7`uXws+UQx$QOguw<7^} zVd82R4eZM?stlui<>X5%#5G4cpz381>r06Orx*O&+A(NGPBgtN4S6UWw2~ZEFna>~ zFR=zp{kQC2j>7jmq(+VgLx|Cx2LK@418#oDevgIURIi4ILDw|ZX|5g5+QT18qE^Ox zXfMDe^B>Wc=$573=790Sq6Ey&Uwzw-wQ2Cv4M*)U0|@jf+>%e#O07~OIb*rILk$jq zb_CjaxL2Qc#gO;vxV2tn-Tdh77Ky$OMRX+~4(B{WN+UYqg|8%|40 zUqD#C!u!B;zrAkx=wtdiSs8ZV2sn7d=}cq^=-@hk;hjqhYYQP*{8uf}x>xb-8@=zY z6Wcee%0)Hk9%96`$zP~bKAb=Q8wWoG699TB6W*f^YP57ESD_}Njgj0jyWB5UpUwYO z#X8w+>f5P&0wl-|%7R81f0sjih z7@7uPqs%T*)~)x|PwtLycU%tt;?~EpoLNr#roq!&Q|(zq`G;;mT#_f2OxKza8$njY zLBvo;`-Hk%uoYoKP-c!*SL`OG&8+?F_Mlxi?QQL?o$&X!=t6p^$#b1k46J1*@79CA zeS;Lonl{ON@FmTU- zVO}X)a;fi-L(rEd)0*XHbvT8XMy@xTd4g+b;)X8!f-3|FoWxD@Lar9Tew4yY#VCRKb6<7cQpkVeUCd9@FLYt8tVL_^!Kc4S%;}( z+w%bnaLTw7y5nxE78(}d0DRt?$}wXt4-&j8&0y4abYdXNfk5N$SS<`O zJ)>H@Am^i8w3Eo$aV8IkyR<&)67edaFiGNNof#lQhC7|g!kQIL?lGb;9!U}-i z?H!^Fp|QZQm$Gij2#(h*wvQLRxZLpQpru@3wPWuZDAu=S+Qq5cylZe1>L zzSGXB!nALm#BZ(haFgW_@*r)&5^+PcXT5SJwXIIUJiOB8gw@l!_wJBU3D!WZWi3LW z_rl0>xPfo}7OyETlu<)(NsDrKTM@9X3m**s_tW2FZR0eZuZ(?^z(Z*7G%|Xz{BT@D z+VBuDo+eH+6X*STdhf1vOAMfPy0MfsloBHA8asYfO<4KbeDJ--l}+2nZQ}&k+E|z!wOE>Jn&5O*d>rXk&W&FI zC;R_h&2c)ORjsXyM9JoEj@2NL>+^O()-e%bR_vucuQUvC;FVG%xuJGT7mOA<8 zpA%U6`^Ag0>raco6|$aG&Fm8I41t5vD73TmPy*>R{Q>@OBq9Sia%0eg5&!Xyk-&#q zu;mX%PqnRHee4nJ3?EZLcToRi{U`AqV`qp_C`Lo1IcYcP{X`DJlW zH94_EJb(1X`hAO}v^n;dGI#ZYT!Rc$CI&J|;#cPlFbFcA*0Vh8>(R7xo3-OZn^#+R zPn@N|ob~+Szvw1fi55JE8==+q2IQHVwV;YsfAWJOdthy3#~J?hjk{QHa19S+@~B<#8cdBu2!NOyJ`VJr?yXXnERXR-O1ePc)0)Osf-Y%WYci@ zQW8a&zKeVOi92==r@zbQdoryfrVD=Is=>~;nbE^E*_EJ#YfhcbI$VRsKj5t9DzUng zbLQcnG@Po}dP>)rMR)h|l_253I|!iDTO^oeZ#X@?6!dYxQ#TsQb}aPpR`V(zO`<|+jcOm7_x>FvDXy-RZ|&kWHv50-FNG}cH?a%GEjNmb^i_V)> z3X|<216#~JM}Gw&)5~Z}cOTs^9XErE8qRvDZ^u0P%)i*Flyx;lE>3gki-FT}xKf9I zL%wON%>{9rs@prVuPowrr~KEYS=A4P)ZBxDVcbo3n4q^Jf3{vNu@cGChHbV@gr;wT zSP(TWAnW|ph@iC(Ps|p-@0Xi|;AB?M5#a+54f+6poP(uIsX9=o7mmBJI``bgL{A8I z4Q^S@iBxjxms!68AiZ_|180GfJ~Pj#8&8-_*%JL8+a_=4-wH4<A9-so&fRe7zdCV~yc>LRF5Lz7P35D+Ig)(jlr>IQqZh&Mw(DVE zZ`NE!m%x#)`@M`F!MK+O*=_Gy z;(r|w33)>|zlbe2RIZ>5VSA}0 zz*Pg4)8)t#Z;Z3cE}}|4F7i6hK-_$75pLF57nVK7Fi3j_p^S?wxvhWwUeQ|q9RH;z zg?;Smd2oepguhASx0>{sMtm$NO6s+F>w(258Q0wD^IcLuL(}X=FJ*>@FK1Q{`NaLn zKrvwnH=wHq*r(1Cmlo1*81jwP)uV`6ItV844eZB3`LUlDeXi~IKZg%uuP|<(#euW% zu6)D-t-g98e3H#;(jzT=o>9RG<^#!N+XYjR<(AJs<)6um^9Z>m5|n4wtiIFGM)t2i z4(=xJUjBrslO(VySRb*oMY*1}9~u-`(^i|sTOy~N`urT|(EAB7yeY6CK|qhbvuGnU-@@}Nmx|thqSA0KH>-qq?zdFT5~?W@4q2@`$0Be|oZCa!<#g?k zisIHTK6<57CKdJF7Rxf$QyV&nDa(Z0rg9VCzcTDY;9uH< zE}b!tmENI$pubx=OPEwq(Jgmi=(aZ?;Lk*t(M75K5+y~(E2Pa3?%~W>oIKi4|Gt`P zC(zCH;hLrG+D;CW+SSfwQe({)pmx@)He4wPNp57{+t`XX-BkUR)tJ3|>z~@ocX7ZZ zn|^hG`C4V1jPs5&zz`HPtGZsZlP(uyNht})6Nflf6RzGK%-h@X?}fO)NwiG+-BeJ+ zQ|*KuRklzB#yjLudUD0)OWL4u0kD~P3Hk>0ewZjygoO&p9Hs}*jG54RiM)iNPQ59) zW@3~xeZA2GIsGLZJyPQ^6m+1g29F3DEQs;}xABKU3vCQ+GzK zuikXchRj8Lc0?CY|+GIv-312aEG|Jf}-4hbaYihb6G)?7bUk z-;@e|sORWw;dUeEOapQvLomg%KWGK8Rk-$Q67UH@iG6gN5OUkgehqrldNwPsdzoO3 zWAa+ROfNz0o|XajyBD|HT7zQ)LOQG3(;d><5y@R9^PTEs>f(E+_H5W!mi}!{PY0$# z&F#BR$&>jIOHf~Coe9l7N(Kxh()gt_fu%2Yx>GKL zO7WZ-ynM)sZ{bRaLd*o~tyyziJ8mwg)jjZ?Ov8L0AM6FDUZmMVjd}Ej`$D=n`x5uM zgLS>6W4-Fv2Bg!EADLv2ihwDVlGo`tOax18fN|nb%b1Pbhl2=j7Wb2v5xgn<7Ip|V zpq0D!Dl+oM>jMd7kfyuyyaZTlXbPo92^-YJKW>I{d+58+nk*Or&|2lWvFSv+=blLW zt8&!KDK<`5y|$I>Z3$1(e@-7?*r>A}YFYA1&zLWpoT6ASA8U}`ew;VWPb?a2Xcw;N zup0zV=$VwsJ8b)~yy3&1MnKtXiXc;LYUBX)dLDHw6G73q7qye%ux=QEd{I63We(gZ zv#O&b{xsod)5E*6<|ZKG6?gmev?WgBvT9|V8ZIVc;+gH&CP0;z1Q^pwf&!K}m0O!w z=NoT}D1=N}F){?}M5wAJ0Jbe3BgJZk)M%(C&%7o5iU0%wDNE&l&{aVN-&<=;{5SfHgVgGJa2U&S zYlIY`=;UtH#hLfGG7jbgX9>)YR;4$GL#I=FYFF2eoR9iW9O;#<>VLrPL`iDz|1=eY$cECuDEU-_e3 z-+_fEA~DdBE@87Jf^J^Uc$AerL=jpoq!r`jL0jq^@_q0#)Oz$i5Ez89Jfg@Ge`gj` zn8L0E2ON>EEQq?p1TzpZ7?bsD1vT+O>}38{RRrZC!$VW5{hQ4=s0@rgZb#0v;wK>6 zO&r`q*3591YJg+kz{TEAdv0E0HN~s(O1hlG!UmDPI~M&ypc>Bnykti1*J`WK8pYjnN*z{9SyTMKv2nD+z5I`hJA=I)N(z< zD!o24w*4q8S=<764>6oz;Vro&+i$N(3)A0}>!plSZSl`nhvLo3$yuF_nKRUFa}#c-**WIHK@DJ6*Y#m@21i%KA?O$s6u}z8R8!#1QOL?1334M2x)suU8vb)AY zizNwC-7JKY(W-Xe+fDAh2T#6XVK{sV0xazeL1=to+r5b_WkERt`7I=hehe#@_! zeTCCy|FT@H2DV{)l+<1m5HyX#iM7nnH}9Wgrrz`U5Y?0_ymtznno7*OUJ&>D_u#wd zCR?jYlMi7LE$!wE0084TJug2#oYmqqqPR=~G9ELZc+0;63k@ZS9@ZaNCf|6DOO4XK z&bX==Hp zTk3MawIc(qz|UlTA@mZEc)0$)W|oX&#apj~ZG9fI0Lf6lc8Y(^voBokObKLO^iRbF z%)2LVlQtaTL-FevFD@79p4(OYXaGZkdZp;TC&lg18My(p%vM`a`zGhI!Bxy*-fp;VE#qo8zar=hJ@f|13KO zNn7c~K_vg#y>}JA@cGTYk-5F3C>|s_1er6wjO8A**!U7$rhfj(!BnWOGgdwTZt$hX zHEYmz3#q?grt1xS6A%4p`AKS3K`Y%tZ@#Xd+tI!B$nkqke$l)+2wLB)G2w_xm*SKkSO_(C7wej4i`USx6v&R4U4BdYAq#klgFf% z882i*uns*4jWcW8G<uO4CT_>^TuA519>uFK?T-wN|h;#C)F%A*s({B{N zpX7r{jEZZh%ITZsFL(LBmtHM|JVZaMHX#C=r5&e^rk9{+L?b`Y=PR$#odKEx)gC7| zotWRc+d8BE{?N5GR}e@}(g_|pr#jRFIyldWNXykGW%f(zk&d;ifr4%$CfynPC-<_) zG81v8l|zPfM~iZ9$U990PoM4jz9Z|cZO-SgSc*281$;^E*RsJK?KEyAUv>(#jwVx7re_}#7B4gXgQ0McX<8Fj)jgOU#@P8{s<_3v)o8hwB640t}rna7dRTkUJ7bFz3`$at< zd@!6Muklxb79L2)k!x~k?oj+~w0@tzgQvZ>z>wOBg-yY*aw7_SUzjD=LQMXXa&z0Q zV2ePg0u&IU4HAwPYZc}HJOz>gqa42*hBqYuawCZ~kn>B8N+!;c`i3R``OK)m1nKj& zjNKRKAlttl$oT}+4x><$x{UP%ZU%5w^yJc6c_uxP3bBpq_?Y7FiH0V){~?846W9eeOpLUPk9tV~T4vbX^)ht`j^woM zuQ@+fLzpZ?+!b+LcxUujbm5`wvV}fl0pO0GW#NaZQn=Ln2a?YhH1+c95+F%>YSVK# z3g>X9i=d$G*Wj|yk(OLRAw`02pXJkK8zOyIn0WzanhYxmIS&O~Al~d~*?*eB&JrA} zp9zvB{>zgQd5qTzHgbPBGi122;TCOr81?lX1ael_C*!VLMSbMm@UejlEFSfqCVlWUJ&n5zgt7xw{5;;SnKryu{`xQVw&#emwuEoPM;tKiuAORLM@8`L-=JI-KW+xbzQCcJE zsR!*4QQze;*#1QvgH0ErlXpq%4~0D1vZj@Y~pB`{U{g)kD z`G+@fAO7plr(>*#nEV`G&(gc@%cmcJaa&~&sx(%Wb_5lSK{6@>SGztroc zyGsUujQi9{_aT%wRM%(Z$2-o&nS9}Q&8Ef${4pAob|;bIj`X^YC+FS?mtitrmxVbh zXkXzeqX~A9&Ii1YzhF2)>-~K?K8>n(PRL52TLo>k>G4tA!e% z)DyFUjch<@AwT6RG!|!7I_etFqS$)NEWKW7dWb4}D(w8fS>pFDql4Vh{`J0S5dy)X zJ^L%BVMhQjGjf~uJIqK(?J49$s(p3^$dqDUzPYpt2;za~z7Ly<_jAav(YaWh9y$M7 zdfh)p^O!#k^dbnIWIO3$oV5k$_}$b>)quHIkz&Pyh{jY%Qsdd7ju?oTJ;yf z)f_v~AP%i*!$~V*S;SIdkQn_P$3pGFMZOiHy0s%}BBiBhe>%InoC*i>{DypjdwM({ zTR^5cSZp39?S~(*oj4kO!=&Q(X6O3s$=aU$tI;UEAMWRjc5UIy?iVUU^2NJQ*#-A! z47c20w3Q$Z7~QS%$)dRY3*zUXaqmq!(s}aKO$>M@x#ZrYypnbb#g?dDknCv2 z$g0Hfvj#>sI_P;WdGf7a`O)OPPHA!vJ8X1c;{+t~Y3|{JV}*{*7%RRSZ3<$-2vgW& zz3?&W2ANl{rHI~96<_HQvc=YyhZ6rGU`@_ye$B3Fa#{AR-l?BokAXTiLhj*dRVO|A zuU;AgrFj=;EM!MxKgiwegp&lAl*nbXZRL`UZd=L6uzn8p*C12ZZ;CUo3587Da&*7U zCerrq`8o6730L8Drw*E){+_nKarz(+&Dp#6GyB%gTfpYH!L`xMZln{PvyfCi(}b9A zp8Z=XF~s_Sr9k2K ze}spHXcy#`PT4l!cZz$fyj!z#p1jPYjY|}d49f# z@EkCe&a%uD27Z55*yz_U2zRL@Z@8qIO) z+4^ZJrC{>M&yzPz{F71*LMz-BJA-;u6wJ4a^Xq@V6bmhyE>dl-JO?<6SuIWOMkTYo z%Y_TgW(DgWjC&Es6npk^(|jO*ZQk{(=c~m%;BDw19au`LFyrd}9Ce&ok43R-RL^TG+&WdkY>JF-Z7we2s2f%}5S2E{=yNVY-*J zsyC!GtDECqPF%^=Lms05Ug5{<%#M0c*f=*VV#W~KqP2Q)fz5GkC1qq$qc(q0fx6q_ zuT&dyX~cNWpzBwV7%lXU-~H+|y35G@8*T#w=hfywrF$A1X_Flz=}ocg@^4H}?o1;b zN_cKJpH_l@&hQ$=>!N#pB%Di;`xu$RBvR>?9nN@w>wpSE_y!P7plC+Fuk-9i8Ux@T z;`Ny*7c71^$P+uY$o$B5HVUk&^NMAPR5h+Z1QZhC5O82^vD1jBCR8+rcNBhiXj?n; zTh1zENoYrA!)z-?_M86;fWu`(z*%4$oK+&LhCyT1Tbk3k*=*_rb6PRj@fi=TFn%`N zZ-#TT*JEOGO!XLMW+IS_E7fjs-(M4I%FMNH%2Mc)3JSLbTj?KHS~b95X~*old0Eq$GxW9!Eu zT1kh>)4@iDQGCoM9!86+WXaYETmlOmNnbe?=qAtn(R5>uyI|Gd+j(^)++FAwODu&c zAiYZW9{JdckywDo*UQVF@U26>2`ADaUw6l8iqw^_4Z7bMn$4Qal-U ztbk+fA$XADUX7$oa~2FO#f?4i=nn1$u<&8Z@AD+~J-mEju;(Mcgzg+oEiQ7&GB*;o zBv0Yp;>n_Ai%UAuw+VQtxEte)+n-qaN%_quZT(&Se?jFj@vehYd)}~@%T6dGtyug; z+*m82!BNX_G4)03bCc_TRGCd87F%vQSe~u*>}^&*tUL(3a-#dudsG0e3?wU%qO;XV zx1y|hiY&~gsdz>;SOi@&l<|UBwiF@kk=m^hZYPek(jJ`?)JV{a-n54+<)=^d0$WxP zi^KE$nL6~%53_~q{tlGKh0B@(>rx>8cgX5!-9MlEf__}>I$nUCFEN^9?!4r%bv%In z+%)m)r6B26<+m^`W1WMDz<$-H@;ZE{6s7#E7QHM_5D9%n)(RLG4Qi`B=Sj5;R$9nU zDpsHLb?>D?`ohqMRN*NCoA3bZ;Q6jxT6Xv=n^ug2ir&Gj?m;F)h8SG)QGDueX#t%Q z^2t9c17j{ko$~#vYe$Qj_o^50a`YEL#Lof7O)!qzmj~Omie1mFGz^T(`#f~M>*SML z(l>6>&lxaVmjx*RO7s=0&xeljii^XFslDH3=;1t25>O>uliu%XAJDy^0{ z&s#Rmvo0~RJ3bikHc=zuy&eu)E!-ZQ&T!LXjsb;=-e0IO_3@t2kP|P6Fobg#SDYDh zOygcahojq`x)m!Ao7PR=Imr7L4zOom3gYJt)mj&17w`kL5(5{80*#P!Lg_*my=JFO z9{WJyR}K>lLt_6lbz$L?uZb78i;-h&hNuF+sT3|5?W?8%h1BoA5vMeq<72QZ{2uEg zSxiVJ7gnQt{ZM-+?{$*?D+Pw;>f)y^2E8DZw=I4OLtY?Ig~8$R)vQ2yze_ZegOgu0 z#(Bd|f4y<4DBoN6Qv>#G07_l1hOfo3Mg#xqOsZoN?ScprhG2rve#WJ8Bwz+Q@dD}O zO3G-q*Gzr_{*G;3M73}IKk#}7H_v9|I_JX)_lTlG1<;CUF%<#cjr`C_pq;IrGynp# z=(W4U`w1}i9kcMY<;ea+v9S4Ba6ec3>;zT?J-fK^dKIFyQ8<&ERUz|7X4Pt{!oRNv zm(%xq-Kyh5#7!svZgPhCbhEoCG5RI->7JAG)?-&#zNmGW-5QC(deyhCV`z*Eh)|3Q zE;3h6NmY$*yy9r8dgA;{@ZQ;Cf>ANG!Et8Cc1L%u{$M2H*!mNe0e0^=ymOQ)ubUFr zymw0Jl7$yA#|X*(PmU~e(w@mU%D8K?(D6Ziges+cG(HY}aA&2jVKmI)*X!NU zhD+M(``|p}1051C?vH(J=j7bK^}RB6>{(h+u6g-dkCCy;?O{U%N@@CGLVxp7;mmJ7 zmPlY7-YMB;{UtOL$s!AM(`8nsdOb|w{thjvL*&n}-1DX6ag!v682 zHs`{;N9#rK+ijt#OJu3gVtlV+8dSQf&!3`;viW1J#ah)rLB9JfFC_}!M_)vK!U zfw&`9U|MPJ|4;1m@T0gYSa&#=boCfkS zjXwMaaFoKg*_SlGQL~Vbxj{BLf7$mz&S9k(FX-3&aG%t}szOKAkY&B-nlCdHsP=YQ z(%(54YXV{a(mHN`86)JEf(CE%RyW^{WggRgex`ANvRQR*@(uL*!eFgpSn6M~7GccPk#Oj?59}JIQb9iN z&ZTU-XDv2IPCtDGZb(@fT&fZr+=}zHD~^Eb<`ZW{&y#-S-`eT8!vyMbPy5pFSS-hC zBVm;{1xS2v7SD(v@Q%LMwkv~7uD8qXKZY*6nadmeI7X4(?7Y(+De|`xSAYrRq(ie| zl_?L#ud^}}0oNOJJI){chACmlaLshZt7Eho0~ zAj5Jw?KzcLXUh2K_{W6OBSUz7!Q2;nEefnu6Jb2`QMjX;IG%ecXX@>W zEZa9%q4CB{@UF}36~YUJ$d0HZv?WSKc|?ya;YM~%rkzjs+H@r~gL0%s??R$CcRK)RnTUPbcrTTm6&V@z@+^nMC2V zhpe}^BhHmSn7aNVh5eD(M`ORX8$lGP^%bEbbR{RfB=db`u(^*akS8Km;Ll4!)^PVt zp=iBO$mGPHFAF{%_f%1snSKJOmG_qg(}A;uAJi;o3qZCS%-dDpdRr<`dhJPljF#l^ z!%1t~A%=?@M-!8`sSe1;>1VD;de8maL(P6O%!SmoUcPs7Vdch3yd%oJ)0>6^HL4#r z)!4t-81>`k2>32iS0Rm^KKo#VdV;(}7 z-{SztH{W@i1ehwgKazLI3kc9IrI`TIkF{pwGnJ~CjWiG5TF6sN=U6K%wTFs|{ZsQ) zAkAo_)O`EbJyK=xnToSgf>Kh`b+~j$uJ+x_gQsyjijl)+YOk%bRqEM;)v;qlkW*eX z7Ub>t%xN^lvyiM(w)Hiw9pRxZksyx%F5_pkZQ}5N5|V$5oxW7S4&BzPz+nODWSQ>_C(lya$>={e0wY4&IyQjdKjTWh)AJ0A2NaZ@{oF$op!^kA_w(Gcx-%cX8ek}wUPUZp!W(=vCp|oxiB;O7g1a0 z{uAU0{{XQ^|1Cckm-VoXl z<6&^MxlWi8^g~bUJi>P&eoSdXFpeWq$ zaU12%gmC_PrTgcw_Ab|;^_2qXw~f`#XAK4pB`bfz^WRe=ve383o{q~&O4?gDW^Jnh zB~1ah%^~PJUlJq9WSPe~M;PeCJ*$Gw--%K3(9z*}KcezxEF;(Qkgezhx}jGe`0I$x zp4-c-c+0uN%v+BtN~2lmIGUA$7Zr$SfyY}x%-YIg01I6iDsUC?&74a8g-^FXTe ze$Xp^g77+OdGK`#Hr$W2Xfctz>-CN0HH330<1{ zB=R6H3D_XdYV8Q6N6S7al+_I{5xHE`Owx}(5|4Y_jeq<7asiZ@z*8GBU4HUl-tt*x zh@h{-r_W|Pg|wx0`!9`2A1}#Qb;yo9oy-Y8-AA|_&q`K0={K|kvVQ)-ZS30loyBp& zVfIRX`htM9DXfd6a;1xd`?KEu`vN*=*PC#26#x#%h{o;<$Z?HD)^6pB#|CXK?WkZ6 zduJ1zU1pWV!4aRWeZe;aBJ5Ovt6&T3uk3Tboi>*hd?bOAS$5wVR&XbmdL}(in@5bE zO8m1m>T9QUNmjV&7%=aA>D2s|2jM9PLU+7&ul&l^p3c9e*@Zpjif}!L&AH;c57Yc8 z!cO5zCT$*jel&F@I`eTI-rM zGi%nA??f=1sF#+zZ*r#-Z1yMwX6IgaSulSD2mBaQi^VmTdX&!)GcceYebkjwU{=u2N(2i>v&w88ZjTO zTILLX#1Hj%=%TsPb^Z>*O8zH`ul37}R2}zd1?oifH&WmJS)E>I#IzWPVcS|9r<(V< zymU2E_gwrC5)~hdb}o6&><*CbAl%ft+@{YB`w|!xXeJ3@MO`7nUY)uybo#x+h3%E* zhb@G7ZD8+$IO=mB7e3P0FKmbpy#Edj9^E4=zL5AlT{UbUaTXUe`pYw~4q+qAeUew( zsK5Piw#Q^kUoU)v7(E*sS#*Pb4!#=hQ*!cqV9Zh7_a=G8yWc-=f5o`9q zu!1~MNYjyE zxt8EfAVGMFi>xL}{On$9-;1hyaM{@;sa>zT)l%}OK9=?3iVuFXJ0rh^R~6V4z0G_1 z@qt5#+E2V~kaJT`V}Vqh@*^pV#%o4yt>u@~yt1r=jojJ;t4cm!4&)z->LMm#?ruAK z2wO{z4L=!Zp?@`h)Far$)y!&;;)QS5JXn^q!|LJUBU7vw0-0XwjP${`Waq5x-NV`6 z{}iNm-9L7AfW$WtAL{DSWK*wf=Q7ZBCMXyv`lwK1@Iv*=-+;a0 zIP?F>by}4cHc$pPT)5R8W|S`|JcMjPt`HBtwvK~cg27!E$R8Qz(a@r-tqw@E{)<4z z>pKsK%Z|1&n{NXOeo51yJbKi1)i)G({sEnm!k!sGW|-kv4gj@uC3<{s99;ad_4Jd8 zk|H->WIguDfag`iTxU?s%hPiB zTlzv5PT70tQ&J91xgVmB#v56Db+ULgI}$B3Yyz$ZaE(wCBLQB*^$$2BkH_ymHI?S5 z#Y9PeKG2$0G|;a*Su)pqM8PD{?i|QYSlSR3WntZYyz38>6kFn^o|$4B@!ZM@PfVwk zqM7O19>3@MaTAcY;`?6heAmurauz-ZC2N~I6QPOn^fyZj=DZBwC}ei;PhE7+&!Ay% z(dlcqF^P_ZlSM9W!Hp^r1DjrMWUiSXW|{$l9J;USKOi6J`*P!!b2j3VVs!6nhQUkM z(rd25#hjs@&-FU3nR!`=&+?XEQZ;qPPa@!U)QvV!&PmndNB{Oba^#n2)?t ztU4sRIx-P*{u|P-JV7UKm9qc1t@UtrrfGa^+hEtxXFjinX7oR6MMb5WCa4hl?Vk~A z!3sk$ep{at7@}rf5&Qdtc0Nb{e8IS8e3xlvV57#9r>4JOsBjbRkSaq&1)c9zjXd;a zC$cZv|2(}e5{@(TH#a@{B{egO^(}8QyzJds^@OB~IKRg9+U0LULB~wK1d*j~bcYQl zXQ-b6!uxjsFd2_*rF<@W@D@c6BQ|GT+Nd#fFt|oN-O<f+3l3s{(>X^HU6Dq9 zgh!u~8s6Lj3C+%!zusd71yMpl@~>ih?7c#?2!!5|jI2{Y0K)SPcKdI@%$flW`UQzdnH3>5Qt!nKHgCkaDN7bzYYwz~n>YEfSzFj( z4}u)wXUM>ZCjQ*Q4FpuR@s-|Sm!~Z^1J!o2No0UC4+U?Y6`j5u`t(;|t0GYLrDF8? zjn2a*!X7a!QZ)M0qRCO9l%V@>>*CA-m9Dk7{!eS$08^IyNF3t(vei>`QB{1Q@R+HZ zPTWh8J8YGgE7HRfkym{%lXF(@~>{Tk?zb_X`KAz`X1tM9+>n!Tzm5F zj&hKD0fYE7`4We_bYDaZ@A0O1AZicEc53VX!AE}Prp&*9$}9B&dPGhvj~v|{9F!V+Cjr)oZzTBLFg;r7u-ATM@jmpmCFu%Z^URPpc5U=_?KWYw;A~H z0-Yb;$pS)%;Hy~TV?hPV#aH&D7lnDJux*JdDQ5e9e&mQE4evw|jiTqPe zA^Aw7qvqct)B!v=Go&uAt_yUFRZvGsqvCwFfyip)!Pn4S3ejZ(sQ`qu{@a?WZrJP4 zSScO?hFvZhdBAiIeFt1O6hD<~FCSeTNb9g=HYf#032Vrkd`;KVsG7g&Gjr+Rx>0KVwFbRhWjWWtB0)}Ew~O%_KJ zF*Z$6DEMl3d~kM+KJt02OVJgE4U;IVMd4^uwoV+y$5@WNf_SIBo$pVia7B(on-Y|B1bU=rI?^|5z7kQKg*Foe= zReoFVB~-QJPm3{us6q1vW8=N2mYda<%cGs1=~e>+pvFWFo_bMA$lv8Qdp3<&XX=uq zAq6irsIt`+N-giVkDg@;N(Fi>5_SDn2Fh92D{A4^9nRZ--gHl! zKluJ}T=w=pe#9qW$zA~91k?i|(=E?VR^H`GK#C7hB3#h+nY*T1ZZw78A9f5WeY=^B zA+nxi+s6kubcaL|9gD~&F{kGM|Kl8eOGfWOSo9( zYf_@;kii51P->dxX6+Jw-MMJ3H|l{Q^Yn$k;_jU{gqtnhiZQP9G~CUrTzfqyriUNa zly?{5W$?nSl*E|$z)7k5?_<>>p?h6e-%Qqy^p$>kpmfOp$OXAkJ@Y_)))i`xmExh- zORM+JUo-wxNI#Q6Ti&=S@)Yx&24|)>f~~=c&kM#nfFsUP4DfC+Wbm(!KUBT4 zOxM|L9*fZt1}3=xOycr=(%&f`|9xtxNA#>aaww_Nmg)_aXY$HvSg|}&FDK)r)*A%y zwk#32Gin|(FXO=3$BZMKJ)P;Zn{3(9N_8QEpDQgFQLwc5J!E`3suUUEq`&*fsZVcR zStry&J~s_z+C5R7?aCjKOu&KOrLfav2$KL3zkxE;^ zq=~ptX)jAT$3oWId1kM0#HKUtzJ#_eJZ3}7aWC%T0`py?;79bTjD7m$F=pZM%`T;0 za-87OcmDw9yS1WA5w|fgC1c^s&tNtyLo)d@!4@ z^~9*4sIl;T+sn1GwRL`Jehaki>5jnz^*5Vux!$mwuq2ah zlMI&hUkzW`fR06qN~{f}y{IQwbVaxd2O{1^vr&m=GH*;1_7f>djg>iSjXE83#pl-? zpRqUi3D=ek4adJ}1il<=7&-uB+|Z>;{czZmmap>ww+QDryyDw{I^I2%?f7}Is-1H! zyI)geplLNT)kcN9G=de&I`u9dRe5oQCi0qTF^-?oa^05nxfTEa#a8(CK$4F7yrMY4pG@ROfsuDTqsFH5ELaCyD84_WQ{1ziX(#cxrzoKb=tQT?k4ur-fU0a z(37mJC;)`w{Hs|fS2(-TM2KT8Ys|Yy#)P|EJZe|fP;I5rLz8mMPe`%t)V@oI^RHaL z%JFY69B}nLN;!~9{FcJQYg6ik|Ja->OM22XFQq#>w|-oVdGcNSVS?Cc*z#Od#I@>c*;Pb&_{{;xhUy5gLF3KF>g#L}GN_gem?C26;+_n*tPhI0O` z{~BF-3R_)!SfLYlU9rd5{Fo{pZ9nCrD)zlv3;sf???e4=N~8kO%oNjr*m_%o&&eVZ zaN4)#%hmY=fSGHeKlecI{2RmN zeXoHAflG$ngbnE(lW2qzIJisF)=rgtGsRlfWBC<}{!BdGsohq075RpJ7z&!+_<+XdwQURGJS6k&(9kdy+vZLFU2-FlU; zSy$gAUu7#gvTf3!k1*Rk(?GzBDfa#C$;E5qk_zAdoc~!V8=~idwzEk`eq(NaMc;-)3vWHt)PEGMVW)CTbz|VCJxAqDa4Exl@lP`3G`im?jP}Zu z73tF!(f2CVto`6v;}qvITCSsax5iYh=GvJ=$3h5*{A4rG$+i1adYj}bhp+jZvU?z~ zVORDEziECN@57lUzZ@j28g@f|=tQ|Vv04FM9e@^83b_eS+`>bRA33gmRv)BCMABvf zRA6H+decfVNGlD%khL{zyQc&9yylYowz*l0M_dmd@8Z5;epw z!0%w0o8ce^&zGLS%(T^h)HNerUQ$+rQjPlzRQU;SS$<2FAK?(p`gM+j=rz4T@&w&* zoje8Gy3MzS_;u6Whuq7>+in=U|1>1bt{iXP0XoLe|9xhCgRhgw6Z)f40r{899gnZK z5*7~i2WhzWt@-&q+6wfGka=|onqr%%t%*MNGj65nDl%#4d+kG?#bR$H??4b@wO@B^ zOfu>VdhGc`#B$uPq!_>S#jmJUxuYpE$H5T3 zl~!?W(uvf%%{p}3B=JaCh-!6}jhNRycnRZXWE|6G+##w8uFxKc6O%pob7RhU>UApH zQdAc6Ux^9)=Uca&)5l(mm(NrEu$lj$w$|M$d}Khf{E3!Qf7OLZw>| zIY*+@IS-+u>4zd{?b)uDB|t5P$iM>vCq@U4s}nc6Y^uSNG$Fk7hLRdU#tzZQ#^=(&M(s z9MrySpe50{`_FIiqmNJO2Usgk4o_3K4&fq_gUZd?cPJ@?N{2grp95}e+|#Lmg$HFS z-iXkfW^>+3!wZwWIZ9M$1%P#+E&6@ql@ur1^?9lLyXHk*bdi^T?)VAUprqCs{d9+z zlbEWqnog=IjFze*GZg*8Mpb{{D^~0Su+Xqm(Uyb^Y&NYcR?dtK=^+1*_Olx=<2wIJ zeiT}BGvvn-N$0F_bI_xDs8JM7#; zvtU{;YdGaVZm(h+hQU2+U1~TFxZ$2Y`*0Ip_qi2xZlP$UMg)GoA(%(zy=XU#)lG-Yx?sTD3i{;O@d9ptL;rsQ~*;C0O zJxJwG{1-Z0$_l7_H&gVx6o|4J1}U|N23(ktf(Bn5zOQio8Ul+Gnu7hTB z?6x#R-VNUhk=pckim|U;5s_e+%}h(*ahI)bz3=PV$2Hdl>TAtKB!=redjBz=m^Q=; z5L!MQk^)Yxea(-@P9QU+yjs-bhmK5QzrXe2zMZ>`M7cS3G+#|jPVsvbeNwZ=b=!z1 zDLccyGn~7Rs|k7kgqE?3A`B!fD!su(WTEm`jF(>w^x<8#rFL%mWmB&Y{lZ>|?wa#9 zH}R3#yyp#6%m5b%npN}W_lE7Ae>hiwtGw`Q!t51AwNUVU4?(@<=ZN*o&2?h#)@cdG zpzCa3(m|$Adt&;!2b&e_=vMrXb;fFC{^|+~ADJ4|)XR|l{sYWqMN30Z3L--n;^ajd zO;$h6YMcWXwrGh3;}^}2z0+hL1x?RuDa1)V2akt?P((W)vvC70AkqS9ZzU(=<}l z4p#ez@$HSeWq98L-B!h3q4uPN+u{w`+JS(R)Q)puZht`jIK)c7P`ISOdG`?)c}o44 z=5y7><6Dy<;kyG;W89~Vd)F>hJy}u<3}<^?K>Y1~xM{I% z$Nhs{Yv_4b4KN6Ze*OE3^5rbZZ!HrHStr^34#~nzw`eayT{eTOoM`W31DvY!AXIz6 zPqII;C+ZTeUl;XG<>r{M5aE%&miu951KDPL*q?Uhc`fx8(h#sO+szgxBX0XKT0GEw!rv8?PssIt=(+AV7o=HgZogjb8{|tsc{}RlbJ@@N|WSxa+ zfR5?Cnf-&AEw?DfJc1 z1$6sJY1B_wjqG)-Q1pRDFJ*TvfTuPqDU*2y^uLu^~c_{Uz;;OnyL&a@HZ7cT+F zMA4!N05cRs}+}K?AyiUiTMN-uz81?fq$E4>Va+^#%WTT%94BH4mdKj;9 zs8GhktDQik_$sFHw(@pgswQ3c&!ERn{J-O*YT><&$|qF(y!|1gbvHJKK;+l1D|$Q_ zSLO>OvfWMq71K-?8RvB$p*2{1+YlAFzeLpUzY3`9B-WOkIABtSxEcleE3U2=`rI0N zn)A1~>h{S6T?{GOKOfL*uNS)$%pYg3tiRM!ZK3hEH^XuSSR;Z!?Q2sXBDdc7-Sg6p z-#T#S=r>K*-Be7T`0%scLo1DW>b&pjGH2dexZ8k+#>VB4AB*+`59B(}PwjVa(>ZSz z!}{cYM&azGn5e9q4?=|v!aQqRd6}jv!hu8|ZTKWt`CUDEEW|v<&*;qu%7R3^~N8mtl?E7{FAk$<}i8z?09#KsSFL zMcUu3QXcvm@lEeF2ju$Jl(kT}{+qYV&rPF( zCUwQy7CRFhyv}#?^IiPj;1^}zT%`reGahmEVh~JUVBO7KBI84rj*YvM|5>pY1v0?T z=TrzPjam-AKKcj$4u8|6EH-!5PkIX{b6VAtnjEpEAih-4PehtDl>8m?ar@Cl5b1^S z7;DEK4Y|Kx_kKd~89ivJ_#WzCzoCUBUTjh85J%aMn&T?pK2d9^mqa5x+88L=R{9GRkhfDtUuRRbB5uL8y ze~KrhVxfmj0yIJQv^LupATa;{WwSW7_t`*Nud#AsTRc8$eDBl8YyTX8;td?lpx0NH z>K75l#Td2P&6ywR8h z?sgxTZ1$*CM%9MyB9zZC$-koR2(*!U7mMC=e#QE}`M2M>V_C{TF1Vr)<)L}*cJWK` ze4tu4OA*k`$`7!=x#LV+c&p#GCh#NW<8watCX@Bensw2D!2bUG9+B~D-!hwgTTL3b z=s;&xj7eLWt{1xE&atP~KIh4DlQuM49K4N+LKv)H4o0N1?M}Xn-27{0SHOb}1P?PE zI1I4hZ6J&33QW!>8SS2U85sQf(mxzPC_lmh>TUh8M4Sk&F;Wo^ z`dRq$*1#_xw1MeMbv?}cr%hC`?;CU23k_yvO@I(s3NP`5+_G`@zzT%A_Ld%A_gdZpLXki!RQTSw@UmSEfmY8EJA~%it6v-D zRufTSn{+ae_&1{XNf)GeI&5X4P2%oT+XKbwDUm?tRua6H;7qWpzN>#02mOkO1QlDV zX98^Z;n2^B5}1)S;ClkjGP_gn9RGM2kO(QU#Zu860SiZhs;>booGtqO?VDN`(=&Uv zFTRnvXny9(w@gXnomUs|*)Ix{fYh(8pO&s6K&j zj0QXZ8vNaN=+mjVaKGmA#(BW^RC{m{NbHXo;R1FZ|Ma-on3XXJG+xPdp{>zXm$mY0b3j%Z31zEWl|ra(pag9>}_IP7#or zWk2i=kUGn>@k==4XyJg>QPBSHnG)vKjR!yDEjF&Pk$9UAkE8;sV6m$vz`H4kT<~1n zeFSv#zUSi7$nk>$EUwCfi$gk@gwgeyvuVL;xgk=YZHQtwT(I`3WViAXj4&15!W;l?r~O)}(fU?{+bD}4p1 z$)|v*waGM{Kl?msmN;`b;L!2u9I_z`=tj}6`lsFcCE-B+?UX{MHb)Q!cpz*hvYPS2 z52pGX;z}1VtAb7n%%bLN!b_JKiLM`hJgC}xp7kfC#I;H`WGmv}av{Gtcx&8P+VV-z z#keG~2M1b_2QIr3^OJBMqc_J!MS|Ez&}MgbK6f>EQzIqUxj^!V?4Q6yr=KL5Hj4eL zYpxcbRgyVmGIaENk<99_r8QfD!=cYmS5ko@V4#m|QcZYdG)hgtju+zAo8zz2Lgg$P zfRL;hu()%dCv*g3y;GEEI*FTkrWDev-Ivz#A^OR2ddF+siX3isTc5eC$Cg=8Drh$pdsHSD3#uLFxB?hRP2g3N);2x2e%UZj0wNyF z>vjOdWS^-bt#K*svryLi+d@YZmI}k7cdjyvsEWdyU4viCkz}EI16xFe3g( zI29X0u!SG6!JWdR+%TXMBcbKS(?k^Eo?IchD9|9k&320i&t1Hb`W!CKa%iN=`X<0@S~sA)w0~&A)A<@Y(2X zsG{WhNeRjVF@COD{P=$s4i}v*Q=Vq;xnm~c93V^po@ISZauo0n&}ymi_SY9;hjIoP z9^Xm#Dx2`e5+QvCtP;QR49qXC5VVmrY8#OWpa*>y`_29 zLGMZKiW{k)kC<5pPLy47xZ2XK&=j}tB0~a266&i;1s)0%#=cE!?XHPP6kFdL&Gq0Z|2{br@w&WaLB_@q6>Ye=E**RP`*?seevHmmq4r@hy(>ZxMUUc zP_E0~jvyi}@uN&AXJ`Kd$2_rF31BKYJ@kLzn+2D@-<>D%aew}`WOCQ^82#Jj`5l)0 zRlsrt34dg+7b`3e7!eZBo!&Y54v0#cnh5j>qT1|#KYu^gbY6Ae&hfQBvmlZ#mWsNV z?NAcr%D0!|R{_@A9rV%`sA&ApRKU-V&v)(2chiQCMQgepaCPVmJ2WOst}Qlw7GM*} z$7`zjay|-7>6=aWRenDdw?}l$8%6&ZzN_(-R$RN-1a07xy$e`iXe3;E5YQ8ehyPN7 zJQThG%=I6Jx%uFbaO0T43NMniP&HOO4p?WX`&QYC2?2_xCbDbm4eDRS8NY%j&la=9 zWukul^_Kuc+Eet2Ts&-j_nr<+skKcG5jakIB!6^Un}4hRW6B%BPFcxl05l)Kt@^tG zt`<~>Huk=6uxa5bt@w<7>RnR2OZK+*`P`ePkd^{}#TA?4^h2Wj*9Ha;ivD1q?)nPh=i6aJ*-KFQ#QzqTrS29-% zXuEcueq(X<(!Gj`HEh-g>;0g2-TeL^23N&V=5-VzeKUys4i&bw$#8xOrvQg`r}@vuT)SNzdh@UwY|g^O z{%+x${lXHl&F6%W>Pj(4WR{`QEtIBv`m%T&(K>by7`Y=|@rS7C z*bCA1OLstZy?J94tO7~Hlbp~R%`K!F2m#bZ_ng0*sR2wGrzsL;XL*I09PQIj%|%ipY~4Oy8Ph!U z6J~F|Y|j?VzsQbV;iceBajcZG2tPL;&n9Cm6 z*D34{h@FzFCOTpU&~&Z*+A#tB{Y17t-)j1@{_rh6ol*;jW^8f54!Ki%3>=6}B8x5Z zJJ({LuI7p-`+&DjE$sm#~DMWi#{h)2Q9H-Bp_-rJ%0ki@qXuF2i_T6wssBb6A{Xx=VOXV(cr)*vaMv zu;j%H!n~h_uJ^ZM-BrpJXr=4JVV0F>jF;h**8u_|Qlq%3F+IhNJ1`>!dCKE{ekFKI z5s1!zOp=H8uf$Wz6u2CQf&05M_^}lIpBr6LOMbiMNq+ zxzgaI2N;Jwv*r<-zNP%}8-z6j(wS@|Bq~DmTU0F8f z{cJPfJ)pvhKA_3h*`8|Y!Io~2CdUXTW~ys#v`*N6$-rR|ZkAtIRzG-FC=qi-abdX9 zInXt`_aPu45&%0a833zLTmo3IlBy(D5r^d+0hoz00kg$-34xrU!<#g{PPy&yDc zm3d)bl4#h;815A``$H<#E6Yf0AtP!%!1w9nFPC-+6R(?kW;^5rJ11gku(Q3b51z-0 ztv3l_gTy%yE2LyyT85&rh~k%P0E0OGb6-H=#XOPaOFg!y{Ea1oZow93qp|IKE;~9v zA))@Q*`(0yT-d-6y34sbU9jWfC*zEW9udTxN=-j9NLdIMx>dgjd+M$Xw_j@A1F)N^ zfc4}K`~IfuN7L`#H_Ld@)v&WWcm*`pT%=A)7X-UDNbb>_zHfeem zOZ=KLg7Tb07wbtPudXtR@G zPT0wa<*&N#l*BF#u@4pPmHe(RO59OUlS;AbtGBFxEwYQ4|;Vi9FRo z&DhU9$b9>fwV`DOL^6) zsc^1L>F?D_7|nn7U`JCy3&^VB^C|blZqjSm zgl_THtq10zn_k@^rDQ*64xF@Z-Z(zuL_dd632xXODIgv4MeNQAOAc$q$uAo-FSccd z1&igZbwbj7MyuOH5HRtzk+lfojDU5Gp6wDx%U?>D6)aEf^w^%WDDfw=81WL6f`Odo zIhaaVBG^CXHZ6WGsDPHVKCt>+g<4b97uIUBa^W}eQp6^(Q65O>yp`v2C4jI-f=x1_Ji7s?VA1mjf5wDerpJ~fg68t(CQI0S@3`SSW_h|QFec+PU>RwccDQ@w8$yi0|yX1GReL=@h-F4zvDb+c!w ziZKuD0C%0^d;yT@HB_SeqVDx?JDX7Jz^C@N0Xzj(6db%y4paXS99V+K1XJNJ#5M3w z%xE6(|3gFm2CP8kdKvItTY#=1hb^A|uP*{G>U#YXf0Ar>JZ#&2;Gna-xu3Ht`S`B5 z-zPd4&zlvnF~IKsHMp19I;Ii6#hxou{JlsCy{mwnzWQ=dWMDEougJ# zu>W4!|NiH$xGRDciOO91-!uHrWMPtk#vlEsSD;+}=iUDMe8?E}v+!(`c+LCSWm3(*-Dv-NcoC_kVB`LcqLcw@(w)gn2Ip{3eg~2hj1< zZQPE))z}{u@+ez+2f*Cf;Q}2WcDN}D?a z8>}uK17RZ1KNZ_^Ym#JN?1C)kRvAH8@D%TLmp~4-vG64tt)$~hFEdu?&1pTepcKq6 zf%Qm*cUOByPgFzX0M$9RE3Q_2TT9evKk)lVBx`l4^3Hu8B&y2Z^Y4WaPRWXI!)1o= zez?iE(Q+y{`*CQiUT_FQ?bUXDbyBwm3P!R|G-QJ#3SGz_xF|yP?rw$N#-^?it=cM| zQVxP^tImGB576|x<8qm1ZDsHK(kqqy5eq*1QPJs+44mp({b9Bf7}*Z>87v8OH(rr@ zSo1!2MZ<&x)|_2#pUOgX=$0tWp^LdQEmFb~b!od|0y&|=y*aO9@$KOaHn6`}xyw6J zh~0GH#8|i7vr+o}URf7vYgE^gPrg;oOAT7k{YTrAMXEMlu~*f8g_rBIu9?UV-LM!| zWhz~7?+s*4s+R%1r9&VKoLNsw=m$dvm{4g%{#jeLA^n5$nY+v&h&M#9qO{@b(W2Xi zKe8^mk78>|z?|Bt5R-JmhPx!BSJYDd{7P|V(9t+dX zqORnOtOKVu(PAvV@|t00wktMx=Ww!Lgvpm~7>5MZQXexmNGro+c!qfLd$;=>oRLeF zMkRNl{lp|AXaC`f-P8|m;#b(ETnx^(=RUN|RDEvaeTDWAZ#>HoGG?ZkU0{D1MD&EgRC+;bfo^y=-6n zb@R(qY0qJtTC_*7%i2uC9Y-Fb&6a?ziJZwNnm6&>yS$A09|5eo&Jo$7joxif`lt$n zl(-r{*KA{y8M)Lp@$A_>r1Ti)PNT{tU*a(^2RwyZ!cq>88iVXxd@)chGdvd2mFY%# zSCd4G$4ym5S7W?FK(@;4oWy)vH)qE*LXEU$r-no25x|&R3>V&NB@MSKH1)Iv=#Xt0quDt6JyxwX=Q623<-s3X0DSciBvG& zQG3APIv=Y$d52#~|I0EMMT^X~#RP@LywP!+;&r#<)9DcJN@OIpmRlNW28zI?prDB! zAi)8MIX<>>q7u>5w%R+mu%g7gM_L}q$iR@xrHAK8<<3Y?FgXQ}(hgk$y&AO)D@c>p zY${bJpqAQa%yYWz2`KG>I;EAaYF+*hyrwHY^g}~K6K2E{G8*B>8lr97h$g8v*>qCS zQSc`+-fqB2q^7vYEy%kxrWIG4HcoKW%{o}wGu+L%a@}nh=tlh1_D64{LS}fyCHJ9} zG|(vTM{217(g~TUNsf-#4ua{pyfLnzBm~0k7o__MI-KG(u74K3uRm&_1hYk~D&3h5 zw_EV>!(yb9LrB%Ql>uo|HVDdD0wzYEdNQ>ge!6D3^%$;f(8=}K=GG0m&r0L4<&eRC z43P|!wr}_04O|i_IkOEVIF_d~C#ww9CPHa_=F3NTO?G|CPWaWwvxnsrc z>Ge8)Ojvn{$xJe1wGu~Gf~Hw=%RIUVo^<(kgIN9Qb0`KWJUC^Ci(TDvWg zoDmR%r*{+}0?w7WX`Cy3)DS*A>KJzKdA`}^B2Q`4-<3un!RA$}`(3E7^EXN;Bjp*M zGSq%gtc_SWO)kT1j-XnZnh_HrKDqnCbpI5p~&xC4D9=lE7gaTd2 zegxYWZh|r>|#noV9z5 z^F1nVHnd?2)}5?zPACl$rpGg5xh%nFt0OA41!7<8$Eu{};ITPwn|6t-?d97Sdl)HE z#iYcte12)p_}v9~VC^1iO;Dnb%*;c8-Q}%6;jIz(-g$ma*HjI;XV!y4Q+(H)wr^C% z%$qP8(KP)aXU6JM9q(RK$)p}9H&wf5>-tjl%dUvWMxC0DolYq&=l$KwE?zA3sUpvh=qE8ZX#e3ExGm@J83RXdhY@5g+IOE6wgmC%O z3V(>6HKidM6M>m`r-XCe6d(g2c8&i~(5#U+Dl!8)8kueF8JxwN=bXW<&a`Fw4=qO) z7*}5p-6%?o(!9m!%KvC=)VSkMRLbtrh2VTf_#(m#`+l>+6w7CNJ zSc|w|s9950_SuDbjhtFfvaa{eoV2iNofru$c2~f-W2C`R{bm~#lMXn~$sv|rb?`_l z>d4M)yYt$+2oDDO<>BF*{g9Hg)yvVz5NW2%+H5|{%O$2__hd0Ur(rj76l%h{%kBxB znr%aeSnsDn3ttILr+slT(U8cSY86ozg-7GiGQD=eBk`DO*?x{~;fj5^KA7UUbp|=! zXU!P;I};S*yyDd4HZoZ?7ax^^s2=q^usP(Rz^QRZvh`~aoFG4)&@7&sUCnAS;}%q< zz^7}f4-@Hpi|SzjE+dUp*Pxc~WoLxB(?~UkYA17t21X!G!QMIDn-(%|0E<#2iWb1O zHhvuA5nI+;NY0cI|C1uQ7Q5-PT($ZHd& zxd~uaC;KBGFOQucak1Q?4hDvCwRf_ydAL4<(8_u{t))0X7L^yR&eCY5TB1=BKX-Dt zbD7y|$l3s@=zikx@ttz?(rnD0g`F4ZaRF$>!CKD;`&)}rx99^XyprF&;I$@BN;Eiz ztUdk?G&&5a=66Je5OZPYKoFT3Z~7yUPcCfOB>3s(+74WAt?>RN$r{|GT^rP}dkqe> zUtu_W07JRNAFf=Lowe?gLKaCH(=&s6hO;BYTn+bw?vd!M!Ku;z|z#E&Bv{3taf zOG8Xb5UW(bM7M)UF`F;T@#bh3)Wc`J@BI<_ThB&vr?=T_PtLWcGHuY_Sq6@HZSsnd z!m_V*1B~REv)Vqmx&x%oY}$KeZuOnU`RXfpchEhNZF0G#s*uze^x0Q>CHc21NDidw zM-yq?5k9`W5#!MP9l4wDqU{_i?k`y@RQBQ_!TV7>a*089=+H8seV@r4+4fQ8f-w}= zihKUun8sS41+=2VTU&){b9oEJf4Txh%Wj^PErMhq9DviMdp1JtUXxRQ9r-;QAlxcdTIDpXGzqun z1hG&^yj#eda!BmlC~<9@7+ms(t`2G!#oSxb?id+eziDh6ke3pzgv4@&lJkb*9%(z=xejZEd%g?2+JFWc5E#Dt2%nlzJm z9zmHvZqBuFyQ@p@hIWtk5o?sTZ4J^5wUbi_4r`pd>J3bIiXVBl(8Y@H-C15;6D8}m z0}w4y%Tu~1*O*Sqx&kA==i;!z7N(2F zSNjRlEOrKv4;NB}@WrNb4Sw&Sq02^t!bZD4nicGyt%-x&jsdUGE~d5Bs9@G~>aC}Y zyNrBr*Hs=qADwW zYSGU;ux9inxD39wP9E&CA@tdzcqB}wq0UO;(^-fdAskKX)*<|$QD;qn4&tmlyEeJr zutqoHYV~zUgS8(uXg;voux4lySyMKmBuDMOpOFoR$dnaw(*vh0vP-Q&P=_}7c%LUC znw;szEUxD!8E4K_TO<9Qwi~dJC$2DqQ6H2BxLf5_Mj8+~c%m8?dEg$0M?pj<~T(;f6rJ-L9~c__y{{$D?UA z>DO^C@J^gpRvCISvNn}ne3wvZU{|y^>O+ z_dFcY8k(*cQh-st%T)UYO5 zhXZw(vcY9wwgw=%AMJtk)|>=XeSBrj_I7ugHZI2?x7;ntj}>yWtX`+@X2W>48!Bvn zyfvWTAT%gvW~#lyF|CYjxNHecY&XIpWE7%a zN=L9<1-MCtb>$CXiTg_L7A_?=q%QROJJ$>(?U+OUZc0%i%5M1N_^E(k!C4`N&WwTC zHWak$KC(N$dt^Didf?+m2t?c3{p?&%=xI+dffU)1lrcJ(<%z?*8J|o$zZ~tT?}7KX zGVo7#t?5cKNVks8MsHW!!D#!B@>{U*RV-6_VMzV_;j^NqpCmrgZR0ijDaS7~luv)j zsg^ce;W*%go(69g;QexG2Ar3V8lbI?VVx9>>ahJaU^xth#%EQQ4MJ|VQ36@n#{%Y5 zsv-MFjkdMtGhj7AhGaeqa#i#+=Kyj!i*sS$LK=i7@N7^V_*Mz7|CVF7 z{mP(j2+A5wrNancA@iF2oZ`|9RRmX z`z{w_5T!P3s(}@TqD{oj)3_NFwr3=yVIe5o3RXM6rFr0Yhy<;)xZ)k0t#$*M6)Qex z^u=#D_#NhvxnH8EYQW!`CS~@zYf0RZ&Ug&0o>m=RHnkeb<%XCzSVHGr163kq7mPf4 zWl_a*jnoEyNp77@rAH!lc|og8n`M?knupylL__AaIsVyYj*g?KN1hqkp=1Wfc>9d& zP;jHCx}y_2a_U&(+N|GHCKiPZ>KC|;tQF~nxohq@hjZgm@G0of8MDISjnJBuEnB)l zJZEs5W~|@DsEynVgSCaSH67;zr_vzDa!0^bn$Y{*uI!#H6$Z(O+ZiQIAy3{!CMZzL z%Kt@5BcvNZ9Opzt$#Wy%{zw!->tkcckGXHeCia`{3FW^Dke{ejEohV56U`#Ber7=C zS?^Tmlw5~L4R68t-c?iMDbba@i2zqlbvN0kbV|i4x>A9ewJ`ztx6?b(-+GzI>FW2U zPDyP9R7YAF?5vj8@XB!csXl?eL%|~(dTrZ3T-hq!2CQm~jb)SVTvF=nqZA5;1yzO3 zX`+Ht??T!KA|1DL*clETDXPK4IXbQL?e6;pZYBd^jxqHtbF&c;IF=W_FhoBuX-xs* zp@ux3T{e|wz1`lZf6JJvZ8CHRyxA@tN!23Cky&59Ou7MR4EuK&rxYi zRa9cxIx*@)zX{fr8bs_#uPK4Zah#KlgG&=q*1Up^wL$D34t^%dl~a!Zd$Uw>-EZ_T z5v+JS;1HJS)cM>9rp5j0zf&I_n!|H+KtHs@XkW+oq(#d^t;&7$=^ZmsEJM0e7=BhW zZH39B1x{to=r5R?gcv9g3SHq1Qkbx;k|r35=)j9GTERfDi2A7rCc|l7a4S}UmX%k2 z93Zcoqc%zkn6|^*Srjp4r7L&I2Njs&IZw~f_OJ&14jobs{#cut-u}2SY%^7Im-yT! zYfhVx>cUBb`{j zYfNiMFh;4LY{w6o$^Cflx#ymH=bkzD|NqUoH@TT~1nMfPjn0ZTOw{RPW}{Q^ zX!)WS){*QWF-jd_r}L%MnJRM zEiH7TipI>m!?26xa|q}1l!B+ht9eCOQc$+S}r#rNGLW=d937Hhf# z-dov9(?EwGAVx6-r<}FxnCaD#0^oj!%DgP-{I;5r-4JD;gdpL-e+HV9(Jbd(y5=0N zvKy9hPRkeRz)U2TMQmdL&1G_2s9H8YbER&9*9FsZ7@DvLIT=kW;BpfSK~)q2z~>~; zh8<9F7@a9yXOB9Gg6z~7t!kotyhJ|5Scv-pSs|L0+M@X**~x?4MsbJ?zmmC>bvijy zC7y)oP^*`S@&N?J%}dQ)e$3GDslleWB%G*=t%Cz8sXHN!5HbX%6PRsFh8Ex>uP=J0 zZ|fO1r*ABs%>P4MmRdjtsPje~Q@wb=IjTn#wCxbyE-~^|ZG(V41|W&i=aaa0bkNFC z`(xy&M}cO&pBJmt)+nCvn*&Z0pnNdut|Ovv6Hll)HDoGZdfV(jgc2JF2?@kBCsMwt zeAm-0(iv)($AtEYNVh;9kJlS#TFLpxX`7eA+?AGRcpOR4Pa|0Q6FiJXj7gBvnj$kF z;Jo$%NO?av$nnvC)e7m_G_4dMI#ex0T%H~~1Q{Rp=e5!c1y|$L7q+_+8<^r;yXX-$ zV!JTcu0hMD3@(ZPd$O77<0m|A@(6kR*`C6VKx2TZnsM=l;FczE(pK`6D5~u?f6v8r z?4`MQFZi0WAkdDZKKJB^9^KmShsZfLd;ap*gJIX6l}wqsS1%8Sbyj1SyB``$+G?0u z&gBsMVX18r5~#jt;#!|JgdV1><3j_J234q#aG{PTK7`(^Z)QkQU5ljsb&k4SVy7>EM2Mi!z|-He}_C`huX=^nXoHqV2-sq#_qLyNT^Tt_bcI zcn!Ll24UV+X0{ESQfPi_S7|}htc!&37s~6FiWCQW=jboFAFFat8%J|KOXScYoC&xE zSF(ulp!AjK#)8z@HXt_?GD-~oOCrk_AoGVzE91en@wt9Z)jP=xiWzZ~>&&eUiDn=o z%5zLGUZntjqWHUp%O+ zb)po+r(dwK98pT%3M52d=qbl#AH3&+q2?q&qKQKvg&Pvb+3z;QGnbRt(JN#n{nwj) z+}&%Bw#%gP8N1G$pBKu!>n77YW7B0-Mk0;U$9UMF@a#%Rd-Ffha4Dg3}757?3Q#5wzJ2@nUE^a?s~+ZS;Mh; z?Nj}+q3dzI;3VQYrf1_X9QvhGa?+26GU%u}X$Ps6%%G$?p7e4OlpTKkgO2 z(qCUIZ9W3_ySeB1p5MK)Lf?&G!W*yVS7ae!;4Pm8clP{#3uCtmf8%F% z73kjV)wWOfY+dooXV@ty-tN0^FCXP=3wE@0^9|!Sg6ccZge|XuTkOXdJ`()n^XvZu DJ99K8 literal 0 HcmV?d00001 diff --git a/多体动力学模块理论手册/多体动力学.md b/多体动力学模块理论手册/多体动力学.md new file mode 100644 index 0000000..3d6a90a --- /dev/null +++ b/多体动力学模块理论手册/多体动力学.md @@ -0,0 +1,560 @@ + + +[TOC] +
    + +# 1. 基于模态叠加法的叶片变形计算 +## 1.1. 叶片模态计算原理 +在本节会出现时间导数和距离导数函数,约定时间导数用变量正上方加点,距离导数变量右上方加撇,如$u(z,t)$为时间t和距离z的函数,$\dot{u}$表示时间导数,$u^{\prime}$表示距离导数。 +### 1.1.1. 叶片模态计算 +把叶片看成是一端固定于轮毂另一端自由的柔性悬臂梁, 基于模态叠加法,叶片侧向变形$u(z,t)$用前N阶模态表示为 +$$ +u(z,t)=\sum_{a=1}^{N}\phi_{a}(z)q_{a}(t) \tag{1.1} +$$ +其中,$\phi_{a}(z)$为第a阶模态的振型函数,是只关于z的函数,z是柔性叶片径向距离,$z=o$表示悬臂固定端, $z=z$表示悬臂自由端。$q_{a}(t)$为第a阶模态广义坐标, 是只关于时间的函数。对于悬臂梁,正交模态的广义坐标习惯上取自由端的变形量,每个模态可用柔性叶片长度无量纲化和归一化,故自由端的模态值为1。 +Rayleigh-Ritz法,$\phi_{a}(z)$可表示为 +$$ +\phi_a(z)=\sum_{b=P}^{N+p-1} D_{a b} \varphi_b(z) \quad(a=1,2, \cdots, N ; b=P, P+1, \cdots, N+P-1)\tag{1.2} +$$ +式中$\varphi_b(z)$为形函数, $D_{a b}$为与第a个模态第b个形函数相关的系数。$u(z,t)$可表示为 +$$ +\begin{gathered} +u(z, t)=\sum_{a=1}^N \phi_a(z) q_a(t)=\sum_{a=1}^N\left[\sum_{b=P}^{N+P-1} D_{a b} \varphi_b(\mathrm{z})\right] q_a(t) \\ +=\sum_{b=P}^{N+P-1} \varphi_b(\mathrm{z}) \sum_{a=1}^N D_{a b} q_a(t)=\sum_{b=P}^{N+P-1} \varphi_b(z) d_b(t) +\end{gathered}\tag{1.3} +$$ +有 +$$ +d_b(t)=\sum_{a=1}^N D_{a b} q_a(t)\tag{1.4} +$$ +式中,$d_b(t)$为形函数$\varphi_b(z)$的广义坐标。 +假定悬臂梁的每个模态可用一个多项式表示,$\varphi_b(z)$可定义为 +$$ +\varphi_b(z)=\left(\frac{z}{Z}\right)^b \tag{1.5} +$$ +在悬臂梁的固定端,形函数的导数必为0, 故有$P \geq 2$ 。本例在柔性叶片计算中取$P = 2$,$N = 5$。叶片模态可表示为 +$$ +\phi_a(z)=\sum_{b=P}^{N+P-1} D_{a b} \left(\frac{z}{Z}\right)^b \tag{1.6} +$$ + +### 1.1.2. 相关系数$D_{ab}$计算 + +对于有N个自由度的保守定常系统,运动方程为 + +$$ +\sum_{b=P}^{N+P-1} m_{a b} \ddot{d}_b(t)+\sum_{b=P}^{N+P-1} k_{a b} d_b(t)=0 \tag{1.7} +$$ +其中,$m_{a b}$、$k_{a b}$与广义坐标$d_b(t)$对应的广义质量、广义刚度。 +当柔性梁以第$m( a = m +) $阶模态振动,有 +$$ +q_a(t)= \begin{cases}Q_a \sin \left(\omega_a t+\beta_a\right) & (a=m) \\ 0 & (a \neq m) +\end{cases} \tag{1.8} +$$ +所以 +$$ +d_b(t)=\sum_{a=1}^N D_{a, b} q_a(t)=D_{m, b} q_m(t)=D_{m, b} Q_m \sin \left(\omega_m t+\beta_m\right) \tag{1.9} +$$ +将(1.9)代入(1.7),有 +$$ +\sum_{b=P}^{N+P-1}\left(-\omega^2 m_{m b}+k_{m b}\right) D_{m, b}=0 \tag{1.10} +$$ +写出矩阵形式为 +$$ +\left(-\omega^2 M_d+K_d\right) D=0 \tag{1.11} +$$ +其中,$M_d$,$K_d$为与广义坐标$d_b(t)$对应的广义质量矩阵、广义刚度矩阵,$D$为模态系数矩阵$N*1$列向量。把(1.11)改写为 +$$ +\left(\frac{1}{\omega^2} \boldsymbol{I}-\boldsymbol{K}_d^{-1} \boldsymbol{M}_d\right) \boldsymbol{D}=\mathbf{0} \tag{1.12} +$$ +由方程式(4.12)可知有N个特征值$\frac{1}{\omega_a^2}$,对应的N个特征向量$D_a$,即 +$$ +\boldsymbol{D}_a=\left(D_{a, p} D_{a, p+1} \cdots D_{a, N+p-1}\right)^T \tag{1.13} +$$ + +### 1.1.3. 广义质量矩阵$M_d$和广义刚度矩阵$K_d$计算 +假设轮毂的旋转角速度为$\Omega$,且柔性叶片在挥舞方向和摆振方向的变形是相互独立的。不考虑叶片的锥角、叶尖质量。基于一维弹性体系统的模态叠加法,质量矩阵元素$m_{ij}$计算式为 +$$ +m_{i j}=\int_0^{R-R_H} \mu_B(z) \varphi_i(z) \varphi_j(z) d z \quad(i=1,2, \cdots, N ; j=1,2, \cdots, N) \tag{1.14} +$$ +式中$\mu_B(z)$为叶片分布线密度,$R$表示风轮半径,$R_H$为轮毂半径。式(4.14)对挥舞和摆振方向都是有效的。 +刚度矩阵元素$k_{ij}$计算式为 +$$ +k_{i j}=\int_0^{R-R_H} E I_B(z){\varphi}_i^{\prime \prime}(z) \varphi_j^{\prime \prime}(z) d z+\Omega^2 \int_0^{R-R_H}\left(\int_r^{R-R_H} \mu_B(\gamma)\left(R_H+\gamma\right) d \gamma\right) {\varphi}_i^{\prime}(z) \varphi_j^{\prime}(z) d z \tag{1.15} +$$ +其中$E I_B(z)$为叶片的分布刚度,$r$为叶片长度的虚变量。式(4.15)中的第二项为离心刚度项。 +至此, 通过式(4.14-4.15)求得广义质量矩阵$M_d$和广义刚度矩阵$K_d$后, 再通过式(4.12)求出模态系数$D_{ab}$,再根据式(4.6)可求出叶片的各阶模态。 + + +## 1.2. 叶片弹性恢复力计算 +柔性或弹性体的弯曲会产生弹性恢复力, 使其恢复到未变形位置。对于柔性叶片,与广义坐标$d_b(t)$对应的广义质量矩阵$M_q$和广义刚度矩阵$K_q$的计算式为 +$$ +m_{i j}=\int_0^{R-R_H} \mu_B(z) \phi_i(z) \phi_j(z) d z \tag{1.16} +$$ +$$ +k_{i j}=\int_0^{R-R_H} E I_B(z) \phi_i(z) \phi_j(z) d z+\Omega^2 \int_0^{R-R_H}\left(\int_r^{R-R_H} \mu_B(\gamma)\left(R_H+\gamma\right) d \gamma\right) \phi_i(z) \phi_j(z) d z \tag{1.17} +$$ +式中,$i=1,2, \cdots, N$,$j=1,2, \cdots, N$,N等于某方向表示叶片变形的模态个数。在接下来的式(4.12)-(4.20)中i,j也如此。 +对于大型风机叶片,由于转速相对较小, 通常可不考虑离心钢化作用。 +叶片固有频率$f_j$公式为 +$$ +f_j=\frac{1}{2 \pi} \sqrt{\frac{k_{i j}}{m_{i j}}} \tag{1.18} +$$ +其中下标j表示第j阶模态频率。 +叶片阻尼矩阵$C_q$的元素$C_{ij}$计算为 +$$ +c_{i j}=\frac{\zeta_j k_{i j}}{\pi f_j} \tag{1.19} +$$ +其中$\zeta_j$为第j阶模态阻尼 +柔性叶片在变形方向的弹性恢复力计算表达式为 +$$ +{ }^{\text {Ela }} F_i^B=-\sum_j k_{i j} q_j(t)-\sum_j c_{i j} \dot{q}_j(t) \tag{1.20} +$$ + +## 1.3. 有结构扭角时叶片变形和转角计算 +叶片在挥舞方向的变形用两个振动模态计算, 而叶片在摆振方向的柔性较小,故其变形用一个模态计算。设$\phi_{1F}(z)$、$\phi_{2F}(z)$为无扭角叶片的一 +、二阶挥舞模态,$\phi_{1E}(z)$为无扭角叶片的一阶摆振模态,$q_{1}(t)$、$q_{2}(t)$为叶尖一、二阶挥舞变形, $q_{3}(t)$为叶尖一阶摆振变形, 为表示方便, 在后面的表达中省略模态函数中的z 和叶尖变形中的t。在叶片局部坐标系中, 局部挥舞曲率可表示为$q_1 {\phi_{1\mathrm{~F}}}^{\prime \prime}+q_2 {\phi_{2 \mathrm{~F}}}^{\prime \prime}$,局部摆振曲率可表示为$q_3 {\phi_{1{\mathrm{E}}}}^{\prime \prime}$。假定局部扭角为$\theta_{\mathrm{tw}}(z)$, 为叶片结构预扭角和叶片桨矩角的和,则面外方向曲率为 +$$ +{u}^{\prime \prime}=\cos \theta_{\mathrm{tw}}\left[q_1 {\phi_{1\mathrm{~F}}}^{\prime \prime}+q_2 {\phi_{2 \mathrm{~F}}}^{\prime \prime}\right]+\sin \theta_{\mathrm{tw}} q_3 {\phi_{1{\mathrm{E}}}}^{\prime \prime} \tag{1.21} +$$ +面内方向的曲率为 +$$ +v^{\prime \prime}=-\sin \theta_{\mathrm{tw}}\left[q_1 {\phi_{1\mathrm{~F}}}^{\prime \prime}+q_2 {\phi_{2 \mathrm{~F}}}^{\prime \prime}\right]+\cos \theta_{\mathrm{tw}} q_3 {\phi_{1{\mathrm{E}}}}^{\prime \prime} \tag{1.22} +$$ +将式(1.21)和(1.22)表示为 +$$ +{u}^{\prime \prime}=q_1 {\varphi_1}^{\prime \prime}+q_2 {\varphi_2}^{\prime \prime}+q_3 {\varphi_3}^{\prime \prime} \tag{1.23} +$$ +$$ +v^{\prime \prime}=q_1 {\varphi_1}^{\prime \prime}+q_2 {\varphi_2}^{\prime \prime}+q_3 +{\varphi_3}^{\prime \prime} \tag{1.24} +$$ +定义${\varphi_1}^{\prime \prime} = \cos \theta_{\mathrm{tw}}{\phi_{1\mathrm{~F}}}^{\prime \prime}$;${\varphi_2}^{\prime \prime} = \cos \theta_{\mathrm{tw}}{\phi_{2\mathrm{~F}}}^{\prime \prime}$;${\varphi_3}^{\prime \prime} = \sin \theta_{\mathrm{tw}}{\phi_{1\mathrm{~E}}}^{\prime \prime}$为面外扭转形函数,$\psi_1^{\prime \prime}=-\sin \theta_{\mathrm{tw}}{\phi_{1\mathrm{~F}}}^{\prime \prime}$;$\psi_2^{\prime \prime}=-\sin \theta_{\mathrm{tw}}{\phi_{2\mathrm{~F}}}^{\prime \prime}$;$\psi_3^{\prime \prime}=\cos \theta_{\mathrm{tw}}{\phi_{1\mathrm{~E}}}^{\prime \prime}$为面内扭转形函数。对面外面内扭转形函数关于z进行两次积分得面外模态$\varphi_1$,$\varphi_2$,$\varphi_3$,和面内模态$\psi_1$、$\psi_2$、$\psi_3$。叶片变形统一在叶根坐标系$b_{m 0}\left(x_{b m 0} y_{b m 0} z_{b m 0}\right)$中表示, 在$x_{b m 0}$方向的变形(挥舞变形)$u$和$y_{b m 0}$方向的变形(摆振变形)$v$可表示为 +$$ +u=\sum_{i=1}^3 q_i \varphi_i \tag{1.25} +$$ +$$ +v=\sum_{i=1}^3 q_i \psi_i \tag{1.26} +$$ +面外转角$\theta_{u}$和面内转角$\theta_{v}$可表示为 +$$ +\theta_u=\frac{\partial u}{\partial z}=\sum_{i=1}^3 q_i \varphi_i^{\prime} \tag{1.27} +$$ +$$ +\theta_v=\frac{\partial v}{\partial z}=\sum_{i=1}^3 q_i \psi_i^{\prime} \tag{1.28} +$$ +$\theta_{u}$沿$y_{b m 0}$轴正向为正,$\theta_{v}$沿$x_{b m 0}$轴负向为正。 +假定叶片弯曲变形后总的叶片长度不变,则在叶片锥坐标系$z_{mc}$轴方向的投影必减少, 减少量表达式为 +$$ +w=-0.5 \sum_{i=1}^3 \sum_{j=1}^3 s_{i j} q_i q_j \tag{1.29} +$$ +式中有$s_{i j}=\int_0^z\left(\dot{\varphi}_i \dot{\varphi}_j^{\prime}+\psi_i^{\prime} \psi_j^{\prime}\right) d z$,为轴向缩减因子。 +根据式(1.25)、(1.26)、(1.29),变形后的叶片节点当前位置向量为 +$$ +\boldsymbol{u}(z, t)=u i_{m c}+v j_{m c}+[r(z)+w] k_{m c} \tag{1.30} +$$ +式中$r(z)$为叶片变形前从叶根到叶素节点的局部半径;$i_{m c}$、$j_{m c}$、$k_{m c}$为第个叶片锥坐标系各轴的基矢量。 + +# 2. 基于模态叠加法的塔架变形计算 +## 2.1. 塔架变形计算 +把塔架看成是倒立的悬臂梁, 且在自由端存在一个点质量$m_{Top}$,这个点质量即为塔顶质量,等于机舱、轮毂及风轮的质量和。假定塔架在前后向(fore-aft)和侧向(side-side)的变形是独立的,且塔架在这两个方向的刚度相等,故在两方向的固有频率和模态也是相同的。只要计算塔架一个方向的模态即可。 +塔架与叶片的模态计算原理及弹性恢复力计算公式是一样的,只是塔架为自由端有质量的倒立悬臂梁,其广义质量矩阵和阻尼矩阵的元素计算稍有区别,为 +$$ +m_{i j}=m_{\mathrm{Top}}+\int_0^H \mu_{\mathrm{T}}(z) \varphi_i(z) \varphi_j(z) d z \tag{2.1} +$$ +$$ +k_{i j}=\int_0^H E I_{\mathrm{T}}(z) \varphi_i^{\prime \prime}(z) \varphi_j^{\prime \prime}(z) d z-g \int_0^H\left[m_{\mathrm{Top}}+\int_z^H \mu_{\mathrm{T}}(\gamma) d \gamma\right] \varphi_i^{\prime}(z) \varphi_j^{\prime}(z) d z \tag{2.2} +$$ +上两式中的$\mu_{\mathrm{T}}(z)$为塔架的分布线密度,$E I_{\mathrm{T}}(z)$为塔架的分布刚度。从式(2.2)可知塔顶质量和塔架自身的重量减少了塔架的刚度。 +塔架在前后方向和侧向的变形各用两个振动模态计算。设$\phi_{\mathrm{1FA}}(z)$、$\phi_{\mathrm{2FA}}(z)$为塔架前后方向的一、二阶模态,$\phi_{\mathrm{1SS}}(z)$、$\phi_{\mathrm{2SS}}(z)$为塔架在侧向的一、二阶模态,$q_{\mathrm{1FA}}(t)$、$q_{\mathrm{2FA}}(t)$为塔顶一、二阶前后向变形,$q_{\mathrm{1SS}}(t)$、$q_{\mathrm{2SS}}(t)$为塔顶一、二阶侧向变形,为表示方便, 在后面的表达中省略模态函数中的z和塔顶变形中的t。在塔基坐标系$t_o(x_{t0}, y_{t0}, z_{t0})$中,塔架在任意高度位置的前后向变形$u_{\mathrm{FA}}$和侧向变形$v_{\mathrm{SS}}$表示为 +$$ +\left\{\begin{array}{c} +u_{\mathrm{FA}}=\phi_{\mathrm{1FA}} q_{\mathrm{1FA}}+\phi_{ \mathrm{2FA}} q_{ \mathrm{2FA}} \\ +v_{\mathrm{SS}}=\phi_{\mathrm{1SS}} q_{\mathrm{1SS}}+\phi_{\mathrm{2SS}} q_{\mathrm{2SS}} +\end{array}\right. \tag{2.3} +$$ +假定塔架在$z_{t0}$轴向无转角,有$y_{t0}$轴向转角$\theta_{\mathrm{FA}}$和$x_{t0}$轴向转角$\theta_{\mathrm{SS}}$为 +$$ +\left\{\begin{array}{l} +\theta_{\mathrm{FA}}=\frac{\partial u_{\mathrm{FA}}}{\partial z} \\ +\theta_{\mathrm{SS}}=\frac{\partial v_{\mathrm{SS}}}{\partial z} +\end{array}\right. \tag{2.4} +$$ +塔架的轴向缩减为 +$$ +\left\{\begin{array}{l} +w_{\mathrm{T}}=-0.5\left(\sum_{i=1}^2 \sum_{j=1}^2 s_{i j \mathrm{FA}} q_{\mathrm{iFA}} q_{j \mathrm{FA}}+\sum_{i=1}^2 \sum_{j=1}^2 s_{i j \mathrm{SS}} q_{i \mathrm{SS}} q_{j \mathrm{SS}}\right) \\ +s_{i j \mathrm{FA}}=\int_0^z \phi_{\mathrm{iFA}}^{\prime} \phi_{j \mathrm{FA}}^{\prime} d z \quad s_{i j \mathrm{SS}}=\int_0^z \phi_{i \mathrm{SS}}^{\prime} \phi_{j \mathrm{SS}}^{\prime} d z +\end{array}\right. \tag{2.5} +$$ +根据(2.3)和(2.5),变形后塔架任意节点当前位置向量为 +$$ +\boldsymbol{u}_{\mathrm{T}}(z, t)=u_{\mathrm{FA}} \boldsymbol{i}_{t 0}+v_{\mathrm{SS}} \boldsymbol{j}_{t 0}+\left[Z(z)+w_T\right] \boldsymbol{k}_{t 0} \tag{2.6} +$$ +式中$Z(z)$为塔架变形前从塔基到塔架节点的距离;$i_{t0}$、$j_{t0}$、$k_{t0}$为塔基坐标系各轴的基矢量。 + +# 3. 整机刚柔混合多体动力学建模 +在多体动力学建模过程中,把水平轴三叶片上风向海上浮式风电机抽象成由浮式平台、塔架、机舱、低速轴、轮毂、叶片、高速轴及电机构成的多体系统。其中浮式平台、机舱、高速轴及电机看成是刚性的,塔架、低速轴及叶片看成柔性的。塔架和叶片的变形用模态叠加法表示,低速轴的柔性用旋转惯量- 刚度-阻尼器表示。用Kane方法建立整机结构动力学模型时,把叶片离散成17个叶素单元体, 把塔架离散成11个单元体。 +## 3.1. 定义参考系 +对于结构复杂的多体系统,为方便描述建立各组成件的连体参考系,并确定各连体参考系间的转换关系。海上浮式风电机坐标系如图3.1所示, 图中$O(XYZ)$为惯性参考系,$p(x_{p}, y_{p}, z_{p})$、$t_i(x_{ti}, y_{ti}, z_{ti})$、$t_p(x_{top}, y_{top}, z_{top})$、$n(x_{n}, y_{n}, z_{n})$、$l(x_{l}, y_{l}, z_{l})$、$h(x_{h}, y_{h}, z_{h})$、$c(x_{mc}, y_{mc}, z_{mc})$、$b_{m0}(x_{bm0}, y_{bm0}, z_{bm0})$及$b_{mi}(x_{bmi}, y_{bmi}, z_{bmi})$分别表示浮式平台参考点$p$、塔架单元节点$t_i$、塔架顶点$t_p$、机舱参考点$n$、低速轴参考点$l$、轮毂参考点$h$、风轮锥顶点$c$、第$m$个叶片叶根节点$b_{m0}$及叶片离散单元(叶素)节点$b_{mi}$处点连体参考系。图中的$d$、$f$、$e$点分别表示平台、机舱及轮毂的质心。$\theta_l$为风轮轴倾斜角,$\theta_c$为叶片锥角。假定惯性参考系坐标轴的基矢量分别为$X=(1,0,0)^T$、$Y=(0,1,0)^T$、$Z=(0,0,1)^T$,根据各参考系间的变换矩阵可求得任意时刻各连体系各坐标轴的单位矢量。 + +
    + 图片描述 +
    图 1: 坐标系示意图
    +
    + +假定浮式平台相对固定参考系的转动为小角转动,塔架和叶片因变形产生的转动也为小角转动,故从惯性系到平台连体系、从塔基系到塔架单元局部系及从叶根系到叶素局部系间的变换可以不考虑转动角顺序,采用小角变换矩阵: +$$ +T_{S A T}\left(\theta_1, \theta_2, \theta_3\right)=\left(\begin{array}{ccc} +\frac{\theta_1^2 b+\theta_2^2+\theta_3^2}{a b} & \frac{\theta_3 a+\theta_1 \theta_2(b-1)}{a b} & \frac{-\theta_2 a+\theta_1 \theta_3(b-1)}{a b} \\ +\frac{-\theta_3 a+\theta_1 \theta_2(b-1)}{a b} & \frac{\theta_1^2+\theta_2^2 b+\theta_3^2}{a b} & \frac{\theta_1 a+\theta_2 \theta_3(b-1)}{\left(\theta_1^2+\theta_2^2+\theta_3^2\right) \sqrt{1+\theta_1^2+\theta_2^2+\theta_3^2}} \\ +\frac{\theta_2 a+\theta_1 \theta_3(b-1)}{a b} & \frac{\theta_1 a+\theta_2 \theta_3(b-1)}{a b} & \frac{\theta_1^2+\theta_2^2+\theta_3^2 b}{a b} +\end{array}\right) \tag{3.1} +$$ +式中$a = \theta_1^2+\theta_2^2+\theta_3^2$,$b = \sqrt{1+\theta_1^2+\theta_2^2+\theta_3^2}$,$\theta_1$、$\theta_2$、$\theta_3$为绕三个坐标轴的转动角。 +根据惯性系到平台参考系的小角变换矩阵有 +$$ +\left(\begin{array}{l} +x_{\mathrm{p}} \\ +y_{\mathrm{p}} \\ +z_{\mathrm{p}} +\end{array}\right)=T_{\mathrm{SAT}}\left(\theta_X, \theta_Y, \theta_Z\right)\left(\begin{array}{l} +X\\ +Y \\ +Z +\end{array}\right) \tag{3.2} +$$ +式中$\theta_X, \theta_Y, \theta_Z$为浮式平台横摇角、纵摇角及艏摇角。$x_p, y_p, z_p$为平台系三个坐标轴的基矢量。 +从平台参考系$p(x_{p}, y_{p}, z_{p})$通过小角变换到塔顶系$t_p(x_{top}, y_{top}, z_{top})$矢量变换式为 +$$ +\left(\begin{array}{l} +x_{\text {top }} \\ +y_{\text {top }} \\ +z_{\text {top }} +\end{array}\right)=T_{\mathrm{SAT}}\left(\theta_{x_{\mathrm{p}}}, \theta_{y_{\mathrm{p}}}, \theta_{z_{\mathrm{p}}}\right)\left(\begin{array}{l} +x_{\mathrm{p}} \\ +y_{\mathrm{p}} \\ +z_{\mathrm{p}} +\end{array}\right) \tag{3.3} +$$ +式中$\theta_{x_{\mathrm{p}}} = -\theta_{SS}$,$\theta_{y_{\mathrm{p}}} = \theta_{FA}$,$\theta_{z_{\mathrm{p}}} = 0$。 +从塔顶系$t_p(x_{top}, y_{top}, z_{top})$的$z_{top}$正向旋转偏航角$\theta_y$到机舱系$n(x_{n}, y_{n}, z_{n})$,机舱系的基矢量 +$$ +\left(\begin{array}{l} +x_{\mathrm{n}} \\ +y_{\mathrm{n}} \\ +z_{\mathrm{n}} +\end{array}\right)=\left(\begin{array}{ccc} +\cos \theta_y & \sin \theta_y & 0 \\ +-\sin \theta_y & \cos \theta_y & 0 \\ +0 & 0 & 1 +\end{array}\right)\left(\begin{array}{l} +x_{\text {top }} \\ +y_{\text {top }} \\ +z_{\text {top }} +\end{array}\right) \tag{3.4} +$$ +从机舱系$n(x_{n}, y_{n}, z_{n})$的$y_{n}$轴正向旋转轴倾斜角$\theta_l$到低速轴系$l(x_{l}, y_{l}, z_{l})$,低速轴系的基矢量为 +$$ +\left(\begin{array}{l} +x_1 \\ +y_1 \\ +z_l +\end{array}\right)=\left(\begin{array}{ccc} +\cos \theta_l & 0 & -\sin \theta_l \\ +0 & 1 & 0 \\ +\sin \theta_l & 0 & \cos \theta_l +\end{array}\right)\left(\begin{array}{l} +x_n \\ +y_{\mathrm{n}} \\ +z_{\mathrm{n}} +\end{array}\right)\tag{3.5} +$$ +从低速轴系$l(x_{l}, y_{l}, z_{l})$的$x_l$轴正向旋转风轮方位角$\theta_a$到轮毂系$h(x_{h}, y_{h}, z_{h})$,轮毂系随风轮转动,其$z_{h}$轴与叶片1方位角一致, +$$ +\left(\begin{array}{l} +x_{\mathrm{h}} \\ +y_{\mathrm{h}} \\ +z_{\mathrm{h}} +\end{array}\right)=\left(\begin{array}{ccc} +1 & 0 & 0 \\ +\cos \theta_a & \sin \theta_a & 0 \\ +-\sin \theta_a & \cos \theta_a & 0 +\end{array}\right)\left(\begin{array}{l} +x_1 \\ +y_1 \\ +z_{\mathrm{l}} +\end{array}\right) \tag{3.6} +$$ +绕轮毂系$h(x_{h}, y_{h}, z_{h})$的$x_{h}$轴正向旋转角度$2 \pi(m-1) / 3$后,再绕$y_{h}$轴负向旋转叶片锥角$\theta_c$得到第$m$个叶片的锥坐标系$c(x_{mc}, y_{mc}, z_{mc})$,通常三个叶片的锥顶点$c$、与轮毂参考点重合,有 +$$ +\left(\begin{array}{l} +x_{\mathrm{cm}} \\ +\boldsymbol{y}_{\mathrm{cm}} \\ +z_{\mathrm{cm}} +\end{array}\right)=\left(\begin{array}{ccc} +\cos \theta_c & 0 & \sin \theta_c \\ +0 & 1 & 0 \\ +-\sin \theta_c & 0 & \cos \theta_c +\end{array}\right)\left(\begin{array}{ccc} +1 & 0 & 0 \\ +\cos [2 \pi(m-1) / 3] & \sin [2 \pi(m-1) / 3] & 0 \\ +-\sin [2 \pi(m-1) / 3] & \cos [2 \pi(m-1) / 3] & 0 +\end{array}\right)\left(\begin{array}{l} +x_{\mathrm{h}} \\ +y_{\mathrm{h}} \\ +z_{\mathrm{h}} +\end{array}\right) \tag{3.7} +$$ +式中的$m=1,2,3$,指代三个叶片。绕锥坐标系$c(x_{mc}, y_{mc}, z_{mc})$的$z_{mc}$轴负向旋转桨矩角$\theta_p$得到叶根系$b_{m0}(x_{m0}, y_{m0}, z_{m0})$,有 +$$ +\left(\begin{array}{l} +x_{\mathrm{bm} 0} \\ +y_{\mathrm{bm} 0} \\ +z_{\mathrm{bm} 0} +\end{array}\right)=\left(\begin{array}{ccc} +\cos \theta_p & -\sin \theta_p & 0 \\ +\sin \theta_p & \cos \theta_p & 0 \\ +0 & 0 & 1 +\end{array}\right)\left(\begin{array}{l} +x_{\mathrm{cm}} \\ +y_{\mathrm{cm}} \\ +z_{\mathrm{cm}} +\end{array}\right) \tag{3.8} +$$ +从叶根系$b_{m0}(x_{m0}, y_{m0}, z_{m0})$的$z_{m0}$轴负向旋转叶素局部扭角$\theta_t$得到叶素局部系$b_{mi}(x_{mbi}, y_{mbi}, z_{mbi})$,有 +$$ +\left(\begin{array}{c} +x_{\mathrm{bmi}} \\ +y_{\mathrm{bmi}} \\ +z_{\mathrm{bmi}} +\end{array}\right)=\left(\begin{array}{ccc} +\cos \theta_t & -\sin \theta_t & 0 \\ +\sin \theta_t & \cos \theta_t & 0 \\ +0 & 0 & 1 +\end{array}\right)\left(\begin{array}{l} +x_{\mathrm{bm} 0} \\ +y_{\mathrm{bm} 0} \\ +z_{\mathrm{bm} 0} +\end{array}\right) \tag{3.9} +$$ + +## 3.2. 定义整机运动自由度 +对三叶片水平轴上风向海上浮式风电机共设立22个自由度,浮式平台设置6个自由度,三个平移自由度$q_1$,$q_2$,$q_3$,三个旋转自由度$q_4$,$q_5$,$q_6$。塔架4个自由度,分别是塔顶一阶、二阶前后向(fore-aft)变形,用$q_7$,$q_9$表示,塔顶一阶、二阶侧向(side-side)变形,用$q_8$,$q_{10}$表示。一个机舱偏航自由度$q_{11}$,一个低速轴转动自由度$q_{12}$,一个低速轴扭转柔性自由度$q_{13}$。每一叶片三个自由度,叶尖一阶、二阶挥舞(flapwise)变形$q_{14+3(m-1)}$、$q_{16+3(m-1)}$和一阶摆振(edgewise)变形$q_{15+3(m-1)}$下标中,${m =1,2,3}$表示不同的叶片。 +## 3.3. 基于Kane法的整机多体结构动力学建模 + +对同一多体系统广义速率有不同的选取方案, 系统中同一质量和同一刚体也将有不同的偏速度和偏角速度。故在选取广义速率时应尽量使对应偏速度和偏角速度的表达式简单,以便于求导运算,易于建立动力学方程和使方程间的耦合程度减弱。选取广义速率$u_r=\dot{q}_r, r=1,2, \cdots, 22$。$\dot{q}_r$为各自由度速度。在论文中用$\omega_X^{(r)}$表示刚体$X$对应第$r$个广义速率的偏角速度,在海上浮式风电机中刚体$X$指代平台$P$、塔架单元$T_i$,塔顶基板$T_{tp}$、机舱$N$、低速轴$L$,电机$G$、轮毂$H$,叶素单元$B_{mi}$。用$v_X^{(r)}$表示质点或参考点$x$的第$r$个广义速率的偏线性速度,同样$x$表示平台参考点$p$;平台质心点$d$、塔架单元节点$t_i$、塔顶点$t_p$、机舱参考点$n$、机舱质心点$f$、轮毂参考点$h$、轮毂质心点$e$及叶素单元节点$b_{mi}$。符号$\overrightarrow{r x y}$指在惯性系中,从点$x$到点$y$的空间位置向量。 +### 3.3.1. 各刚体偏角速度 +浮式$p$的偏角速度 +$$ +\boldsymbol{\omega}_{\mathrm{P}}^{(1)}=0 \quad \boldsymbol{\omega}_{\mathrm{P}}^{(2)}=0 \quad \boldsymbol{\omega}_{\mathrm{P}}^{(3)}=0 \quad \boldsymbol{\omega}_{\mathrm{P}}^{(4)}=\boldsymbol{X} \quad \boldsymbol{\omega}_{\mathrm{P}}^{(5)}=\boldsymbol{Y} \quad \boldsymbol{\omega}_{\mathrm{P}}^{(6)}=\boldsymbol{Z} \quad \boldsymbol{\omega}_{\mathrm{P}}^{(r)}=0(r=7,8, \cdots, 22) \tag{3.10} +$$ +塔顶$T_{tp}$的偏角速度 +$$ +\begin{aligned} +& \boldsymbol{\omega}_{\mathrm{Ttp}}^{(1)}=0 \quad \boldsymbol{\omega}_{\mathrm{Ttp}}^{(2)}=0 \quad \omega_{\mathrm{Ttp}}^{(3)}=0 \quad \omega_{\mathrm{Ttp}}^{(4)}=X \quad \omega_{\mathrm{Ttp}}^{(5)}=Y \quad \omega_{\mathrm{Ttp}}^{(6)}=Z \quad \omega_{\mathrm{Ttp}}^{(7)}=\phi_{1 F A}^{tt p} y_p \\ +& \boldsymbol{\omega}_{\mathrm{Ttp}}^{(8)}=-\phi_{1 S S}^{ttp} \boldsymbol{x}_p \quad \boldsymbol{\omega}_{\mathrm{Ttp}}^{(9)}=\phi_{2 \mathrm{FA}}^{\mathrm{ttp}} \boldsymbol{y}_p \quad \boldsymbol{\omega}_{\mathrm{Ttp}}^{(10)}=-\phi_{2 S S}^{\mathrm{ttp}} \boldsymbol{x}_p \quad \boldsymbol{\omega}_{\mathrm{Ttp}}^{(r)}=0(r=11,12, \cdots, 22) +\end{aligned} \tag{3.11} +$$ +其中$\phi_{1 F A}^{tt p}$、$\phi_{2 F A}^{tt p}$为塔顶前后向一、二阶模态,$\phi_{1 SS}^{tt p}$、$\phi_{2SS}^{tt p}$为塔顶侧向一、二阶模态。 +机舱$N$的偏角速度 +$$ +\omega_{\mathrm{N}}^{(r)}=\omega_{\mathrm{T}_{\mathrm{tp}}}^{(r)}(r=1,2, \cdots, 10) \quad \omega_{\mathrm{N}}^{(11)}=z_n \quad \omega_{\mathrm{N}}^{(r)}=0(r=11,12, \cdots, 22) \tag{3.12} +$$ +低速轴$L$的偏角速度 +$$ +\omega_{\mathrm{N}}^{(r)}=\omega_{\mathrm{N}}^{(r)}(r=1,2, \cdots, 11) \quad \omega_{\mathrm{H}}^{(12)}=\boldsymbol{x}_l \quad \omega_{\mathrm{H}}^{(13)}=\boldsymbol{x}_l \quad \omega_{\mathrm{H}}^{(r)}=0(r=14,15, \cdots, 22) \tag{3.13} +$$ +电机$G$的偏角速度 +$$ +\omega_{\mathrm{G}}^{(r)}=\omega_{\mathrm{N}}^{(r)}(r=1,2, \cdots, 11) \quad \omega_{\mathrm{G}}^{(12)}=\lambda x_l \quad \omega_{\mathrm{G}}^{(r)}=0(r=13,14, \cdots, 22) \tag{3.14} +$$ +式中$\lambda$为传动比。 +叶素单元${B_{mi}}$的偏角速度 +$$ +\begin{aligned} +& \boldsymbol{\omega}_{\mathrm{B} m i}^{(r)}=\omega_{\mathrm{H}}^{(r)}(r=1,2, \cdots, 13) \quad \omega_{\mathrm{B} m i}^{(14+3 m-3)}=\left(\varphi_l^{b m i}\right)^{\prime} y_{\mathrm{b} m 0}-\left(\psi_l^{b m i}\right)^{\prime} x_{\mathrm{b} m 0} \\ +& \boldsymbol{\omega}_{\mathrm{B} m i}^{(15+3 m-3)}=\left(\varphi_2^{\mathrm{bmi}}\right)^{\prime} y_{\mathrm{b} m 0}-\left(\psi_2^{\mathrm{bmi}}\right)^{\prime} x_{\mathrm{b} m 0} \quad \boldsymbol{\sigma}_{\mathrm{B} m i}^{(16+3 m-3)}=\left(\varphi_3^{\mathrm{bmi}}\right)^{\prime} \boldsymbol{y}_{\mathrm{b} m 0}-\left(\psi_3^{\mathrm{b} m i}\right)^{\prime} x_{\mathrm{b} m 0} +\end{aligned} \tag{3.15} +$$ +### 3.3.2. 质心或参考点偏线速度 +平台参考点$p$的偏线速度为 +$$ +\boldsymbol{v}_{\mathrm{p}}^{(1)}=\boldsymbol{X} \quad \boldsymbol{v}_{\mathrm{p}}^{(2)}=\boldsymbol{Y} \quad \boldsymbol{v}_{\mathrm{p}}^{(3)}=\boldsymbol{Z} \quad \boldsymbol{v}_{\mathrm{p}}^{(r)}=0(r=4,5, \cdots, 22) \tag{3.16} +$$ +平台质心$d$的偏线速度为 +$$ +\begin{gathered} +v_{\mathrm{d}}^{(1)}=v_{\mathrm{p}}^{(1)}(r=1,2,3) \quad v_{\mathrm{d}}^{(4)}=\omega_{\mathrm{p}}^{(4)} \times \overrightarrow{r p d} \quad v_{\mathrm{d}}^{(5)}=\omega_{\mathrm{P}}^{(5)} \times \overrightarrow{r p d} \\ +v_{\mathrm{d}}^{(6)}=\omega_{\mathrm{P}}^{(6)} \times \overrightarrow{r p d} \quad v_{\mathrm{d}}^{(r)}=0(r=7,8, \cdots, 22) ; +\end{gathered} \tag{3.17} +$$ +式中$\overrightarrow{r p d}$为平台质心$d$到参考点$p$的向量。 +塔架顶点$t_{tp}$的偏线速度为 +$$ +\begin{aligned} +& v_{\mathrm{ttp}}^{(r)}=v_{\mathrm{p}}^{(r)}(r=1,2,3) \quad v_{\mathrm{tpp}}^{(4)}=\omega_{\mathrm{P}}^{(4)} \times \overrightarrow{r p t_{t p}} \quad v_{\mathrm{ttp}}^{(5)}=\omega_{\mathrm{P}}^{(5)} \times \overrightarrow{r p t_{t p}} \quad v_{\mathrm{tpp}}^{(6)}=\omega_{\mathrm{P}}^{(6)} \times \overrightarrow{r p t_{t p}} \\ +& v_{\mathrm{ttp}}^{(7)}=x_{\mathrm{p}}-\left(S_{\mathrm{FA} 11}^{\mathrm{ttp}} q_7+S_{\mathrm{FA} 12}^{\mathrm{ttp}} q_9\right) z_{\mathrm{p}} \quad v_{\mathrm{ttp}}^{(8)}=x_{\mathrm{p}}-\left(S_{\mathrm{SS} 11}^{\mathrm{ttp}} q_8+S_{\mathrm{SSt} 12}^{\mathrm{ttp}} q_9\right) z_{\mathrm{p}} \\ +& v_{\mathrm{ttp}}^{(9)}=x_{\mathrm{p}}-\left(S_{\mathrm{FA} 12}^{\mathrm{ttp}} q_7+S_{\mathrm{FA} 22}^{\mathrm{ttp}} q_9\right) z_{\mathrm{p}} \quad v_{\mathrm{ttp}}^{(10)}=x_{\mathrm{p}}-\left(S_{\mathrm{SS} 12}^{\operatorname{ttp}} q_8+S_{\mathrm{SS} 22}^{\operatorname{ttp}}q_{10}\right) z_{\mathrm{p}} \\ +& v_{\text {ttp }}^{(r)}=0(r=11,12, \cdots, 22) +\end{aligned} \tag{3.18} +$$ +式中$S_{\mathrm{FA} 11}^{\mathrm{ttp}}$、$S_{\mathrm{FA} 12}^{\mathrm{ttp}}$、$S_{\mathrm{FA} 22}^{\mathrm{ttp}}$为塔顶前后向缩减因子,$S_{\mathrm{SS} 11}^{\mathrm{ttp}}$、$S_{\mathrm{SS} 12}^{\mathrm{ttp}}$、$S_{\mathrm{SS} 22}^{\mathrm{ttp}}$为塔顶侧向缩减因子;$\overrightarrow{r p t_{t p}}$为参考点$p$到塔架顶点$t_{tp}$的位置向量。塔架单元节点$t_i$的偏线速度计算与塔顶点的相似,只是位置向量和缩减因子不同。 +机舱质心$f$的偏线速度 +$$ +\begin{aligned} +& v_f^{(r)}=v_{\mathrm{tpp}}^{(r)}(r=1,2,3) \quad v_f^{(r)}=v_{\mathrm{ttp}}^{(r)}+\omega_{\mathrm{N}}^{(r)} \times \overrightarrow{r o f}(r=4,5, \cdots, 10) \\ +& \boldsymbol{v}_f^{(11)}=\omega_{\mathrm{N}}^{(11)} \times \overrightarrow{r o f} \quad v_f^{(r)}=0(r=12,13, \cdots, 22) +\end{aligned} \tag{3.19} +$$ +式中$\overrightarrow{r o f}$为参考点$o$到机舱质心$f$的位置向量。 +轮毂参考点$h$的偏线速度 +$$ +\begin{aligned} +v_{\mathrm{h}}^{(r)}=v_{\text {ttp }}^{(r)}(r & =1,2,3) \quad v_{\mathrm{h}}^{(r)}=v_{\mathrm{ttp}}^{(r)}+\omega_{\mathrm{N}}^{(r)} \times \overrightarrow{r o h}(r=4,5, \cdots, 10) \\ +v_{\mathrm{h}}^{(11)} & =\boldsymbol{\omega}_{\mathrm{N}}^{(11)} \times \overrightarrow{r o h} \quad v_{\mathrm{h}}^{(r)}=0(r=12,13, \cdots, 22) +\end{aligned} \tag{3.20} +$$ +式中$\overrightarrow{r o h}$为参考点$o$到轮毂参考点$h$的位置向量。 +轮毂质心$e$的偏线速度 +$$ +\begin{aligned} +v_{\mathrm{e}}^{(r)}=v_{\mathrm{h}}^{(r)}(r=1,2,3) \quad \nu_{\mathrm{e}}^{(r)}=\boldsymbol{v}_{\mathrm{h}}^{(r)}+\omega_{\mathrm{H}}^{(r)} \times \overline{r h e}(r=4,5, \cdots, 11) \\ +\boldsymbol{v}_{\mathrm{e}}^{(12)}=\boldsymbol{\omega}_{\mathrm{H}}^{(12)} \times \overrightarrow{r h e} \quad \boldsymbol{v}_{\mathrm{e}}^{(13)}=\boldsymbol{\omega}_{\mathrm{H}}^{(13)} \times \overrightarrow{r h e} \quad \boldsymbol{v}_{\mathrm{h}}^{(r)}=0(r=14,15, \cdots, 22) +\end{aligned} \tag{3.21} +$$ +式中$\overrightarrow{r h e}$为参考点$h$到轮毂质心$e$的位置向量。 +叶素节点$b_{mi}$的偏线速度 +$$ +\begin{aligned} +& v_{\mathrm{b} m i}^{(r)}=v_{\mathrm{h}}^{(r)}(r=1,2,3) \quad v_{\mathrm{bmi}}^{(r)}=v_{\mathrm{h}}^{(r)}+\omega_{\mathrm{H}}^{(r)} \times \overrightarrow{r h b_{m i}}(r=4,5, \cdots, 11) \\ +& \boldsymbol{v}_{\mathrm{bmi}}^{(\mathrm{I2)}}=\boldsymbol{\omega}_{\mathrm{H}}^{(12)} \times \overrightarrow{r h b_{m i}} \quad \boldsymbol{v}_{\mathrm{b} m i}^{(13)}=\boldsymbol{\omega}_{\mathrm{H}}^{(13)} \times \overline{r h b_{m i}} \\ +& \boldsymbol{v}_{\mathrm{b} m i}^{(14+3 m-3)}=\varphi_1^{\mathrm{b} m i} \boldsymbol{x}_{\mathrm{b} m 0}+\psi_1^{\mathrm{b} m i} \boldsymbol{y}_{\mathrm{b} m 0}-\left(S_{11}^{\mathrm{bmi}} q_{14+3 m-3}+S_{12}^{\mathrm{b} m i} q_{15+3 m-3}+S_{13}^{\mathrm{b} m i} q_{16+3 m-3}\right) z_{\mathrm{b} m 0} \\ +& \boldsymbol{v}_{\mathrm{b} m i}^{(15+3 m-3)}=\varphi_3^{\mathrm{b} m i} \boldsymbol{x}_{\mathrm{b} m 0}+\psi_3^{\mathrm{bmi}} \boldsymbol{y}_{\mathrm{b} m 0}-\left(S_{31}^{\mathrm{bmi}} q_{14+3 m-3}+S_{32}^{\mathrm{bmi}} q_{15+3 m-3}+S_{33}^{\mathrm{bmi}} q_{16+3 m-3}\right) z_{\mathrm{b} m 0} \\ +& v_{\mathrm{b} m i}^{(16+3 m-3)}=\varphi_2^{\mathrm{b} m i} x_{\mathrm{b} m 0}+\psi_2^{\mathrm{bmi}} y_{\mathrm{b} m 0}-\left(S_{21}^{\mathrm{b} m i} q_{14+3 m-3}+S_{22}^{\mathrm{bmi}} q_{15+3 m-3}+S_{23}^{\mathrm{b} m i} q_{16+3 m-3}\right) z_{\mathrm{b} m 0} +\end{aligned} \tag{3.22} +$$ +其中$\varphi_1^{\mathrm{b} m i} $、$\varphi_2^{\mathrm{b} m i} $、$\varphi_3^{\mathrm{b} m i} $表示第$m$个叶片第$i$个叶素节点处挥舞方向的扭转形函数,$\psi_1^{\mathrm{bmi}}$、$\psi_2^{\mathrm{bmi}}$、$\psi_3^{\mathrm{bmi}}$表示摆振方向的扭转形函数,$S_{11}^{\mathrm{b} m i}$、$S_{12}^{\mathrm{b} m i}$ 等为轴向缩减因子 +己知偏角速度后,可用它们表示刚体的角速度和角加速度为 +$$ +\omega_{\mathrm{X}}=\sum_{r=1}^{22} \omega_{\mathrm{X}}^{(r)} \dot{q}_r \tag{3.23} +$$ +$$ +\boldsymbol{\beta}_{\mathrm{X}}=\dot{\boldsymbol{\omega}}_{\mathrm{X}}=\sum_{r=1}^{22} \boldsymbol{\omega}_{\mathrm{X}}^{(r)} \ddot{q}_r+\sum_{r=1}^{22} \dot{\omega}_{\mathrm{X}}^{(r)} \dot{q}_r \tag{3.24} +$$ +己知偏线速度后,各质心或参考点的线速度和线加速度为 +$$ +v_x=\sum_{r=1}^{22} v_x^{(r)} \dot{q}_r \tag{3.25} +$$ +$$ +\alpha_x=\dot{v}_x=\sum_{r=1}^{22} v_x^{(r)} \ddot{q}_r+\sum_{r=1}^{22} \dot{v}_x^{(r)} \dot{q}_r \tag{3.26} +$$ +由式(3.24)和(3.26)可知角加速度和线加速度由两部分组成,一部分与自由度的加速度有关,另一部分与自由度的加速度无关。 + +### 3.3.3. 整机动力学方程 +海上浮式风电机各组成件受到的惯性力和力矩主要包括:浮式平台受到的惯性力$F_d$和惯性力矩$L_p$,塔架离散单元惯性力$F_{ti}$;机舱平动惯性力$F_f$和和偏航惯性力矩$L_{N}$,电机转动惯性力矩$L_{G}$;轮毂惯性力$F_H$和惯性力矩$L_{H}$,叶片离散单元惯性力$F_{Bmi}$。各刚体的惯性力统一表示为 +$$ +F_x=m_x a_x=m_x\left(\sum_{r=1}^{22} v_x^{(r)} \ddot{q}_r+\sum_{r=1}^{22} \dot{v}_x^{(r)} \dot{q}_r\right) \tag{3.27} +$$ +式中的下标$x$指代浮式平台质心$d$、塔架单元节点$t_i$及机舱质心$f$、轮毂质心$e$及叶素节点$b_{mi}$,$m_x$为质心或节点处的质量。根据动量矩定理,各惯性力矩统一表示为 +$$ +\begin{gathered} +\boldsymbol{L}_{\mathrm{Xx}}=\boldsymbol{\omega}_{\mathrm{X}} \times\left[I_{\mathrm{Xx}} \boldsymbol{x}\left(\boldsymbol{x} \cdot \boldsymbol{\omega}_{\mathrm{X}}\right)\right]+I_{\mathrm{Xx}} x\left(x \cdot \boldsymbol{\beta}_{\mathrm{X}}\right) \\ +=\boldsymbol{\omega}_{\mathrm{X}} \times\left[I_{\mathrm{Xx}} x\left(x \cdot \omega_{\mathrm{X}}\right)\right]+I_{\mathrm{Xx}} x\left[x \cdot\left(\sum_{r=1}^{22} \boldsymbol{\omega}_{\mathrm{X}}^{(r)} \ddot{q}_r+\sum_{r=1}^{22} \dot{\omega}_{\mathrm{X}}^{(r)} \dot{q}_r\right)\right] +\end{gathered} \tag{3.28} +$$ +式中的$X$指代平台$P$、机舱$N$、轮毂$H$及电机$G$,$x$指代过质心的转动轴的单位矢量,$I_{\mathrm{Xx}}$刚体$X$在$x$方向的转动惯性,根据Kane方法, 广义惯性力$F_r^*$为 +$$ +F_r^*=\sum_x F_x \cdot v_x^{(r)}+\sum_{\mathrm{X}} L_{\mathrm{X}} \cdot \boldsymbol{\omega}_{\mathrm{X}}^{(r)}(r=1,2, \cdots, 22) \tag{3.29} +$$ + +多体系统的主动力包括系统外部对系统内物体的作用力,也包括系统内物体的相互作用力。对于海上浮式风电机,广义主动力包括:广义重力 $ F^g_r $ 为 +$$ +F^g_r = m_d g Z \cdot v_d^{(r)} + \sum_i m_{ti} g Z \cdot v_{ti}^{(r)} + m_f g Z \cdot v_f^{(r)} + m_e g Z \cdot v_e^{(r)} + \sum_m \sum_i m_{bmi} g Z \cdot v_{bmi}^{(r)} \tag{3.30} +$$ +叶素节点处的气动力 $ F_{bmi}^A $ 和气动力矩 $ L_{Bmi}^A $,广义气动力为 +$$ +F_r^{Aero} = \sum_m \sum_i F_{bmi}^A \cdot v_{bmi}^{(r)} + \sum_m \sum_i L_{Bmi}^A \cdot \omega_{Bmi}^{(r)} \tag{3.31} +$$ + +电机传动反力矩 $ L^G $,由变速变桨控制策略确定,对应广义力为 +$$ +\begin{cases} +F_{13}^G = -\lambda L^G \\ +F_r^G = 0 \, (r \neq 13) +\end{cases} \tag{3.32} +$$ + +偏航刚度和阻尼产生的广义偏航力矩 $ L^{Yaw} $ 为 +$$ +\begin{cases} +L_{11}^{Yaw} = -K_{Yaw} q_{11} - C_{Yaw} \dot{q}_{11} \\ +L_r^{Yaw} = 0 \, (r \neq 11) +\end{cases} \tag{3.33} +$$ + +柔性传动轴刚度和阻尼产生的广义力矩 $ L^{Dri} $ 为 +$$ +\begin{cases} +L_{13}^{Dri} = -K_{Dri} q_{13} - C_{Dri} \dot{q}_{13} \\ +L_r^{Dri} = 0 \, (r \neq 13) +\end{cases} \tag{3.34} +$$ + +柔性塔架广义刚度和阻尼产生的广义恢复力 $ { }^{\mathrm{Ela}} \boldsymbol{F}_r^{\mathrm{T}} $ ($ r = 7, 8, 9, 10 $),和柔性叶片广义刚度和阻尼产生的广义恢复力 $ { }^{\mathrm{Ela}} \boldsymbol{B}_r^{\mathrm{T}} $ ($ r = 14, 15, \dots, 22 $) 都采用式 (1.20) 计算,$ r $ 为其它自由度时,广义恢复力为 0。 + +浮式平台参考点处受到的载荷分为两部分,即与平台运动加速度有关的力和无关的力,表示为 +$$ +F_i^{Platform} = F_i^{Platform1} + F_i^{Platform2} \\ +\left\{\begin{array}{l} +F_i^{\text {Platform1 }}=-A_{i j} \ddot{q}_j \\ +F_i^{\text {Platform2 }=}=F_i^{\text {Hydrostatic }}+F_i^{\text {Wave }}-\int_0^t K_{i j}(t-\tau) \dot{q}_j(\tau) d \tau+F_i^{\mathrm{Lines}}+{ }^{(2)} F_i+F_i^{\text {Monison }} +\end{array}\right. \tag{3.35} +$$ + +浮式平台参考点处受到的与平台运动加速度无关的力及力矩 $F_p^{Hydro}$、$L_p^{Hydro}$,与平台运动加速度有关的附加惯性力及力矩 ${ }^{\mathrm{A}} F_{\mathrm{p}}^{\text {Hydro }}$、${ }^{\mathrm{A}} L_p^{Hydro}$ 可表示为 + +$$ +\left\{\begin{array}{l} +\boldsymbol{F}_p^{\text {Hydro }}=F_1^{\text {Plafform2 }} \boldsymbol{X}+F_2^{\text {Platform2 } } \boldsymbol{Y}+F_3^{\text {Platform2 }} \boldsymbol{Z} \\ +\boldsymbol{L}_{\mathrm{p}}^{\text {Hydro }}=F_4^{\text {Platform2 } } \boldsymbol{X}+F_5^{\text {Platform2 }} \boldsymbol{Y}+F_6^{\text {platform2 }} \boldsymbol{Z} +\end{array}\right. \tag{3.36} +$$ +$$ +\left\{\begin{array}{l} +{ }^{\mathrm{A}} \boldsymbol{F}_{\mathrm{p}}^{\text {Hydro }}=\sum_{i=1}^6\left(-A_{1 i} X-A_{2 i} \boldsymbol{Y}-A_{3 i} Z\right) \ddot{q}_i \\ +{ }^{\mathrm{A}} \boldsymbol{L}_{\mathrm{p}}^{\text {Hydro }}=\sum_{i=1}^6\left(-A_{3 i} X-A_{4 i} \boldsymbol{Y}-A_{5 i} Z\right) \ddot{q}_i +\end{array}\right. \tag{3.37} +$$ + +广义水动力 $F_r^{Hydro}$ 为 +$$ +F_t^{Hydro} = F_p^{Hydro}\cdot v_p^{(r)} + L_p^{Hydro}\cdot \omega_p^{(r)} + A_p^{Hydro}\cdot v_p^{(r)} + A_p^{Hydro}\cdot \omega_p^{(r)} \tag{3.38} +$$ + +总广义主动力 $F_r$ 为 +$$ +F_r=F_r^{\mathrm{g}}+F_r^{\text {Aero }}+F_r^{\mathrm{G}}+F_r^{\mathrm{Yaw}}+F_r^{\mathrm{Dri}}+{ }^{\text {Ela }} F_r^{\mathrm{T}}+{ }^{\text {Ela }} F_r^{\mathrm{B}}+F_r^{\text {Hydro }} \tag{3.39} +$$ + +对于有 22 个自由度的海上浮式风电机系统,凯恩动力学方程表示为 +$$ +F_r + F_r^* = 0 \, (r=1,2,...,22) \tag{3.40} +$$ + +通过整理动力学方程表示为 +$$ +M_{22\times 22}({q}_r, t)\ddot{q}_{22\times 1} = f_{22\times 1}(\dot{q}_r,q_r,t) \tag{3.41} +$$ + +式中 $M_{22\times 22}$为对称质量矩阵: $\ddot{q}_{22\times 1}$ 为 22 个自由度加速度构成的列向量, 有 $\ddot{q}_{22\times 1} = (\ddot{q}_1 \ddot{q}_2 ... \ddot{q}_{22})^T$;$f_{22\times 1}(\dot{q}_r,q_r,t)$ 为与自由度加速度$\ddot{q}_r$无关的力向量。 +