kill(all);

cons. l'equil. del quarto di telaietto

eqtx: XB=0 $

eqty: YA=0 $

eqtz: -P+ZA+ZB=0 $

eqrOx: -P*b+ZA*b+MxA=0 $

eqrOy: P*a-ZB*a+MyB=0 $

eqrOz: 0=0 $

risolvo il sistema lineare, versione naive

linsolve ([eqtx,eqty,eqtz,eqrOx,eqrOy],[ZA,YA,XB,MyB,MxA]) , globalsolve=True;

caratteristiche di azione interna sui tratti di struttura

tratto BC, ascissa curvilinea s in [0,b]

Mfxi_BC : ZB*s;

Mt_BC   : MyB;

Teta_BC : ZB ;

tratto AC, ascissa curvilinea s in [0,a]

Mfxi_AC : ZA * s;

Mt_AC   : MxA ;

Teta_AC : ZA  ;

calcolo energia potenziale su quarto di struttura

U_AC : integrate(

                   Mfxi_AC^2 / 2 / E / Jxixi
+                    Mt_AC^2 / 2 / G / Kt
+ shearcoeff_eta * Teta_AC^2 / 2 / G / A

, s , 0 , a

);

U_BC : integrate(

                   Mfxi_BC^2 / 2 / E / Jxixi
+                    Mt_BC^2 / 2 / G / Kt
+ shearcoeff_eta * Teta_BC^2 / 2 / G / A

, s , 0 , b

);

U:U_BC + U_AC;