kill(all);

calcolo in ipotesi di parete sottile

A_ps : integrate(                        t * rm ,theta,0,2*%pi);

J_ps : integrate(  (rm * cos(theta))^2 * t * rm ,theta,0,2*%pi);

calcolo senza ipotesi di parete sottile

A : 
integrate(
    integrate(                        
        r 
        ,theta,0,2*%pi
    )
,r,rm-t/2,rm+t/2)$
fullratsimp(%);

J : 
integrate(
    integrate(                        
        (r * cos(theta))^2 * r 
        ,theta,0,2*%pi
    )
,r,rm-t/2,rm+t/2)$
fullratsimp(%);

grafico

plotme : fullratsimp(ev(J_ps/J,t=eta,rm=1));
wxplot2d(plotme,[eta,0,1]);

calcolo instabilità

dati sezione

rm : 10;
t  : 1 ;

l:1000;

E:210000;

Fcr : %pi^2*E*J_ps/(2*l)^2,infeval,numer;