numero i punti come 12345 a partire dal vincolo in alto in figura

pulisco memoria

kill(all);

provo a definire le reazioni vincolari per SOLO equilibrio
(nel caso mi servano per ricavare i grafici delle caratteristiche di sollecitazione, 
o siano richieste dal testo)
(ovvero, nel caso no mi servano, non le ricavo...)

tz : -F+Z1+Z3+Z5 = 0;

r2y : -Z3*c-Z5*2*c = 0;

r2x : (Z1+Z5)*b=0;

linsolve( [ tz,r2y,r2x ] , [Z1,Z3,Z5] ), globalsolve=true;

diagrammi di momento flettente, tratto 12

Mf1 : - Z1 * t;

Mf2 : 0 ;

Mt  : 0 ;

U12 : integrate(Mf1^2/2/E/J+Mf2^2/2/E/J+Mt^2/2/G/Kt, t , 0 , b);

diagrammi di momento flettente, tratto 23

Mf1 : (F-Z1) * t;

Mf2 : 0 ;

Mt  : Z1*b ;

U23 : integrate(Mf1^2/2/E/J+Mf2^2/2/E/J+Mt^2/2/G/Kt, t , 0 , c);

diagrammi di momento flettente, tratto 54

Mf1 : Z5 * t;

Mf2 : 0 ;

Mt  : 0 ;

U54 : integrate(Mf1^2/2/E/J+Mf2^2/2/E/J+Mt^2/2/G/Kt, t , 0 , b);

diagrammi di momento flettente, tratto 43

Mf1 : Z5 * t;

Mf2 : 0 ;

Mt  : Z5 * b ;

U43 : integrate(Mf1^2/2/E/J+Mf2^2/2/E/J+Mt^2/2/G/Kt, t , 0 , c);

energia potenziale elastica totale

U : U12 + U23 + U54 + U43;

calcolo freccia

d : diff(U,F);

definisco geometria per valutazione numerica

dim:	[a=800,b=1000,c=600,F=1000,C=500000,
	 E=210000,G=210000/2/(1+3/10),
	 J=(40^4-36^4)*%pi/64,Kt=(40^4-36^4)*%pi/32];

ev(d,dim);
%,numer;