pulisco la memoria

kill(all);

inserisco i diagrammi di sollecitazione
in funzione di P, CO, CA, CB

definisco delle funzioni ausiliarie lineari

una prima funzione vale 1 sul primo estremo del tratto di
trave, e 0 sul secondo

variabile indipendente è un'ascissa curvilinea
adimensionale xi che vale 0 al primo estremo e 1 al secondo

N1 : 1 - xi ;

la seconda funzione vale 0 sul primo estremo 
e 1 sul secondo estremo. Tali funzioni
hanno andamento lineare

N2 : xi ;

ora posso definire una quantità linearmente variabile
da un f1 a un f2 in xi da 0 a 1 come

f1*N1 + f2*N2;

vado a definire le caratteristiche di sollecitazione

             |-----P-----|----CB----|----CA----|----CO----|

N_OA :    +P        +CB/b        +0       -CO/b       $
N_AB :    -P*l/a    -CB*a/b/l  -CA*b/a/l  +CO*l/a/b   $
N_BO :    +P*b/a    +0         +CA/a      -CO/a       $

Mf_OA:    +0       +0          +0*N1+CA*N2 +CO*N1+0*N2$
Mf_AB:    +0       +0*N1+CB*N2 +CA*N1+0*N2 +0         $
Mf_BO:    +0       +CB*N1+0*N2    +0       +0*N1+CO*N2$

manca taglio

definisco l'energia potenziale elastica sulla struttura
in funzione di P, CO, CB, CA

U:
    integrate(
        (
            +  N_OA^2/2/E/A
            + Mf_OA^2/2/E/J
            +  T_OA^2/2/G/A * eta
        ) *a 
        ,xi, 0, 1
    )
;

mancano contributi altri tratti