kill(all);

The nodes are numbered from the load application point "1", to the intermediate supported point "2", and finally to the clamped end "3".

considero i tratti in senso orario, a partire dal punto di applicazione del carico (punto 1), 
passando per il punto di appoggio intermedio (punto 2) fino ad arrivare all'incastro (punto 3).

Mf_12 :    0*(1-xi) +       F*a*xi$
Mf_23 :    0*(1-xi) + (F*b-R*b)*xi$

Mt_12 :   0$
Mt_23 : F*a$

U: integrate( (Mf_12^2/2/E/J + Mt_12^2/2/G/Jp) * a , xi, 0,1 )
  +integrate( (Mf_23^2/2/E/J + Mt_23^2/2/G/Jp) * b , xi, 0,1 );

delta_R : diff(U,R);

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

The following passages are NOT required by the exercise.

U : ev(U);

delta_F : diff(U,F);

dim:[a=800,b=1000,c=600,F=1000,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);