Jump to content

SPAM !!!!!!!!!!!!!!!! [THREAD UFFICIALE!!!]


principino1984

Recommended Posts

se volete anzi senza il se... beccatevi sta roba...

 

kill(all);

 

U0():=1/(2*E*A)*(1/rho^2*Xi_Mf+Xi_N)*l;

 

Xi_Mf:1/3*(2*(C_0^2+C_1^2+C_2^2)+C_0*C_1+C_1*C_2+C_2*C_0);

 

Xi_N: (+X +1/sqrt(3)*Y +1/l/sqrt(3)*(2*C_2-C_0-C_1))^2+ (-X +1/sqrt(3)*Y +1/l/sqrt(3)*(2*C_0-C_1-C_2))^2+ (+1/2*X -2/sqrt(3)*Y +1/l/sqrt(3)*(2*C_1-C_2-C_0))^2;

 

geom:[A=2*%pi*r*delta,rho=r/sqrt(2),delta=xi*r,l=eta*r];

 

U1():=ratsimp(ev(U0(),geom,infeval));

 

eqns:[diff(U1(),C_0)=0,diff(U1(),C_1)=0,diff(U1(),C_2)=0];

 

iper:solve(eqns,[C_0,C_1,C_2]);

iso:[C_0=0,C_1=0,C_2=0];

 

U_iper() :=ev(U1(),iper);

U_iso() :=ev(U1(),iso);

 

KX_iper() :=X/ratsimp(ev(diff(U_iper(),X),Y=0));

KY_iper() :=Y/ratsimp(ev(diff(U_iper(),Y),X=0));

KX_iso() :=X/ratsimp(ev(diff(U_iso(),X),Y=0));

KY_iso() :=Y/ratsimp(ev(diff(U_iso(),Y),X=0));

 

ratio_X() :=KX_iper()/KX_iso();

ratio_Y() :=KY_iper()/KY_iso();

 

N_0 :X+1/(l*sqrt(3))*(2*C_2-C_0-C_1);

Mf_0 :C_0;

 

sigma:N_0/A+Mf_0/(%pi*r^3*delta/64)*r;

 

sigma_ratio:ev(ratsimp(ev(ev(sigma,iper)/ev(sigma,iso),geom,infeval)),X=1,Y=0);

 

plot2d(ratio_X(),[eta,1,10]); /*,["xlabel,l/r ratio"],["ylabel,KX_iper/KX_iso ratio"]);*/

 

plot2d(sigma_ratio,[eta,0,250]); /*,["xlabel,l/r ratio"],["ylabel,sigma_iper/sigma_iso ratio"]);*/

 

se qualcuno usa maxima capirá...:fiufiu::fiufiu:

Link to comment
Share on other sites

se volete anzi senza il se... beccatevi sta roba...

 

kill(all);

 

U0():=1/(2*E*A)*(1/rho^2*Xi_Mf+Xi_N)*l;

 

Xi_Mf:1/3*(2*(C_0^2+C_1^2+C_2^2)+C_0*C_1+C_1*C_2+C_2*C_0);

 

Xi_N: (+X +1/sqrt(3)*Y +1/l/sqrt(3)*(2*C_2-C_0-C_1))^2+ (-X +1/sqrt(3)*Y +1/l/sqrt(3)*(2*C_0-C_1-C_2))^2+ (+1/2*X -2/sqrt(3)*Y +1/l/sqrt(3)*(2*C_1-C_2-C_0))^2;

 

geom:[A=2*%pi*r*delta,rho=r/sqrt(2),delta=xi*r,l=eta*r];

 

U1():=ratsimp(ev(U0(),geom,infeval));

 

eqns:[diff(U1(),C_0)=0,diff(U1(),C_1)=0,diff(U1(),C_2)=0];

 

iper:solve(eqns,[C_0,C_1,C_2]);

iso:[C_0=0,C_1=0,C_2=0];

 

U_iper() :=ev(U1(),iper);

U_iso() :=ev(U1(),iso);

 

KX_iper() :=X/ratsimp(ev(diff(U_iper(),X),Y=0));

KY_iper() :=Y/ratsimp(ev(diff(U_iper(),Y),X=0));

KX_iso() :=X/ratsimp(ev(diff(U_iso(),X),Y=0));

KY_iso() :=Y/ratsimp(ev(diff(U_iso(),Y),X=0));

 

ratio_X() :=KX_iper()/KX_iso();

ratio_Y() :=KY_iper()/KY_iso();

 

N_0 :X+1/(l*sqrt(3))*(2*C_2-C_0-C_1);

Mf_0 :C_0;

 

sigma:N_0/A+Mf_0/(%pi*r^3*delta/64)*r;

 

sigma_ratio:ev(ratsimp(ev(ev(sigma,iper)/ev(sigma,iso),geom,infeval)),X=1,Y=0);

 

plot2d(ratio_X(),[eta,1,10]); /*,["xlabel,l/r ratio"],["ylabel,KX_iper/KX_iso ratio"]);*/

 

plot2d(sigma_ratio,[eta,0,250]); /*,["xlabel,l/r ratio"],["ylabel,sigma_iper/sigma_iso ratio"]);*/

 

se qualcuno usa maxima capirá...:fiufiu::fiufiu:

nn stai bene...:hypo:

Link to comment
Share on other sites

se volete anzi senza il se... beccatevi sta roba...

 

kill(all);

 

U0():=1/(2*E*A)*(1/rho^2*Xi_Mf+Xi_N)*l;

 

Xi_Mf:1/3*(2*(C_0^2+C_1^2+C_2^2)+C_0*C_1+C_1*C_2+C_2*C_0);

 

Xi_N: (+X +1/sqrt(3)*Y +1/l/sqrt(3)*(2*C_2-C_0-C_1))^2+ (-X +1/sqrt(3)*Y +1/l/sqrt(3)*(2*C_0-C_1-C_2))^2+ (+1/2*X -2/sqrt(3)*Y +1/l/sqrt(3)*(2*C_1-C_2-C_0))^2;

 

geom:[A=2*%pi*r*delta,rho=r/sqrt(2),delta=xi*r,l=eta*r];

 

U1():=ratsimp(ev(U0(),geom,infeval));

 

eqns:[diff(U1(),C_0)=0,diff(U1(),C_1)=0,diff(U1(),C_2)=0];

 

iper:solve(eqns,[C_0,C_1,C_2]);

iso:[C_0=0,C_1=0,C_2=0];

 

U_iper() :=ev(U1(),iper);

U_iso() :=ev(U1(),iso);

 

KX_iper() :=X/ratsimp(ev(diff(U_iper(),X),Y=0));

KY_iper() :=Y/ratsimp(ev(diff(U_iper(),Y),X=0));

KX_iso() :=X/ratsimp(ev(diff(U_iso(),X),Y=0));

KY_iso() :=Y/ratsimp(ev(diff(U_iso(),Y),X=0));

 

ratio_X() :=KX_iper()/KX_iso();

ratio_Y() :=KY_iper()/KY_iso();

 

N_0 :X+1/(l*sqrt(3))*(2*C_2-C_0-C_1);

Mf_0 :C_0;

 

sigma:N_0/A+Mf_0/(%pi*r^3*delta/64)*r;

 

sigma_ratio:ev(ratsimp(ev(ev(sigma,iper)/ev(sigma,iso),geom,infeval)),X=1,Y=0);

 

plot2d(ratio_X(),[eta,1,10]); /*,["xlabel,l/r ratio"],["ylabel,KX_iper/KX_iso ratio"]);*/

 

plot2d(sigma_ratio,[eta,0,250]); /*,["xlabel,l/r ratio"],["ylabel,sigma_iper/sigma_iso ratio"]);*/

 

se qualcuno usa maxima capirá...:fiufiu::fiufiu:

 

stai + che bene :D:ROTFL:

Link to comment
Share on other sites

Povero apix che è sempre all'uni..come direbbe lui..CHE PATACCA!!!

:uhahah::uhahah::uhahah:

 

apix ... for President """""""""""""""":fiufiu::nutkick::nonno::uhahah:

 

 

lasciamo stare va che il patacca del mio prof oggi ha voluto tirare fino alle 18:15... 4ore e 15 di lezione. per la miseria almeno è bella (progetto motore) ma i neuroni dopo un pò mi mollano in pausa eterna e punto stop... e questo che andava per la sua strada... sono arrivato in quel di fidenza alle 20:00 panicoooooooo meno male che domani è la festa dei lavoratori:nutkick::nutkick:

Link to comment
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

×
×
  • Create New...