kill(all);

Considero come punto 1 il punto di applicazione della forza F e sostituisco all'appoggio 
intermedio la reazione vincolare R: successivamento imporrò la condizione di compatibilità
con il vincolo sostituito.
Per definire le funzioni lineari, utilizzo un'ascissa curvilinea normalizzata xi, che 
vale 0 nel primo punto del tratto considerato e 1 nel punto finale.

errato: 
Mf_34: (R-F)*c*xi;
il momento flettente non si annulla in xi=0, ma eguaglia -F*c

Mf_34: (R-F)*c*xi;

Mf_12: -F*c*xi;
Mf_23: -F*c;
Mf_34: (R-F)*c*xi-F*c;

utilizzerei solo i contributi di momento flettente e torcente, trascurando 
sforzo normale e taglio

ometterei quindi 

N_12: 0;
N_23: -F;
N_34: 0;

altrimenti risulterebbe opportuno introdurre anche in contributi energetici 
(e quindi deformativi) associati al taglio

Espliciterei piuttosto l'assenza di momento torcente sulla struttura.

integro ora, sui vari tratti di trave, i contributi di energia potenziale elastica.

U: integrate( (Mf_12^2/2/E/J ) *c, xi, 0, 1)
  +integrate( (Mf_23^2/2/E/J ) *a, xi, 0, 1)
  +integrate( (Mf_34^2/2/E/J ) *c, xi, 0, 1);

calcolo cedimento al punto di applicazione di R, utilizzando Castigliano

delta_R: diff(U,R);

risolvo la reazione incognita R imponendo delta_R nullo

linsolve(delta_R=0,R),globalsolve=true;

ridefinisco U per valutazione di R secondo la nuova ipotesi di cui sopra
(saltando questo passo U potrebbe rimanere funzione di R e non della sua espressione risolta in funzione di P)

U: ev(U);

riduco a denominatore comune per snellire la forma

U:fullratsimp(U);

calcolo freccia sotto forza F

delta_F: diff(U,F);

valuto numericamente su dimensionamento specifico

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,Jp=(40^4-36^4)*%pi/32];

ev(delta_F,dim)$
%,numer;

controllo su FEM per avere una verifica indipendente