pulisco la memoria

kill(all);

definisco (per assegnazione) il momento flettente sui tratti AB e BC

MfAB : R*s ;
MfBC : R*(b+s)-P*s ;

UAB : integrate( MfAB^2/2/E/J ,  s , 0 , b );

UBC : integrate ( MfBC^2/2/E/J , s, 0, a );

U : UAB + UBC ;

U : expand( U );

U;

calcolo frecce in A e B, positive se verso il basso; utilizzo il secondo teorema di Castigliano

fA : - diff( U , R );

fB : + diff( U , P );

essendo la risposta della struttura assunta lineare, dovrebbe risultare che
fA = alpha * P + beta * R; ricaviamo tali coefficienti

alpha : factor( coeff(fA,P) );
beta  : factor( coeff(fA,R) );

verifichiamo che le espressioni di fA tal quale e fA espresso come combinazione lineare siano equivalenti

ratsimp( fA - (alpha*P + beta*R) );

impongo nullo fA, ricavando il valore della reazione vincolare R

defR : linsolve( [fA=0] , [R] );

R : ev( R , defR );

R;

ratsimp( ev( fA ) );

fB : factor( ev( fB ) );

fB;