pulisco la memoria

kill(all);

qui dovrei definire schema statico e simboli...
TODO...

definisco le equazioni di equilibrio

eqtx : -NA = 0 $

eqty : +P/2 -NB = 0 $

eqrOz : +CA +NA*(l+r)+CB-NB*r $

nell'ultima equazione ho lasciato implicito
lo "=0"

utilizzo la funzione di soluzione sistemi 
lineari di equazioni

linsolve (      
    [ eqtx , eqty , eqrOz ]
, 
    [ NA , CB , NB ] 
);

inserisco le caratteristiche di azione interna
(= di sollecitazione per la sezione di trave)
lungo i tratti di struttura

tratto 1, curvilineo

N1 : P/2 * sin(phi) $

T1 : P/2 * cos(phi) $

Mf1 : -CA + P/2*r*sin(phi) $

tratto 2 , rettilineo, (NON) funzione di
una ascissa curvilinea "s"

N2 : P/2 $

T2 : 0 $

Mf2 : -CA + P/2 * r $

calcolo energia pot. elastica sui tratti

tratto 1

U1 : integrate(
    (
              Mf1^2 /2/E/A/delta/r
        +      N1 ^2 /2/E/A
        + xi * T1 ^2 /2/G/A
        -      N1*Mf1  /E/A      /r
    )*r
,
    phi,0,%pi/2
) ;

manipolazioni di base, per vedere l'effetto
che fanno

U1 : expand(U1);

procedo col tratto 2

U2 : integrate(
    +     Mf2^2/2/E/J
    +      N2^2/2/E/A
    + xi * T2^2/2/G/A
,
    s, 0, l
);

U2 : expand(U2);

calcolo l'energia potenziale elastica totale 
(sul quarto di maglia di catena...),
e procedo con Castigliano

U : U1 + U2$

utilizzo Castigliano per trovare la rotazione
della sezione al punto di applicazione di CA 
(punto A)

rotA : diff( U, CA , 1 );

impongo rotazione nulla in A in coerenza
con l'originale vincolo doppiopendolo

sols : linsolve(  rotA=0  , CA  );

sols;

ridefinisco CA sulla base del valore dello stesso
che soddisfa la condizione di compatibilità
col vincolo

CA : ev( CA , sols );

energia potenziale elastica sul quarto di
maglia di catena, in funzione del carico applicato
P, e dei parametri geometrici

U : fullratsimp (ev( U ) );