/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 0.8.7 ] */

/* [wxMaxima: comment start ]
pulisco la memoria
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
kill ( all ) ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
definisco equazioni di equilibrio
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
equaz. equilibrio traslazione x,y,z
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
eqtx  :  XB = 0 ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eqty  : YA = 0 ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eqtz  : +ZB -P +ZA = 0 ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
equaz. equilibrio rot. x,y,z attorno ad O
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
eqrxO : +CA +ZA*b -P*b = 0 ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eqryO : +CB -ZB*a +P*a = 0 ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eqrzO : 0 = 0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
imposto sistema di equazioni di equilibrio,
trascuro la rotazione z in quanto non indipendente dalle altre
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
eqns : [ eqtx , eqty , eqtz , eqrxO , eqryO ] ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
definisco incognite (ZB è un parametro)
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
unks : [ XB, CB, YA, ZA, CA ] ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
risolvo il sistema lineare
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
sols : linsolve( eqns , unks ) , globalsolve=true  ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
controllo le assegnazioni
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
XB;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
YA;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
ok, le variabili sono definite
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
definisco i momenti flettente (componente non nulla 
a curvatura fuori piano) e torcente
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
              | contributo ZB  | contributo P |
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
Mf_BC :         +ZB * x             + 0  $
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Mt_BC :         +ZB * a         - P * a  $
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Mf_AC :         +ZB * y         - P * y  $
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Mt_AC :         +ZB * b             + 0  $
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
ricavo l'energia potenziale elastica della struttura
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
U : 
integrate(
  Mf_BC^2/2/E/J + Mt_BC^2/2/E/Jp ,
  x,
  0,
  b
)
+
integrate(
  Mf_AC^2/2/E/J + Mt_AC^2/2/E/Jp ,
  y,
  0,
  a
) ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
ridefinisco U come la riduzione a denominator comune 
(+ qualche semplificazione elementare) di se stesso
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
U : fullratsimp( U );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
utilizzo il teorema di Castigliano per ricavare loù
spostamento in z del punto B 
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
dispzB : diff ( U , ZB ) ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
risolvo equazione di compatibilita al vincolo
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
linsolve (   dispzB = 0  ,  ZB  ) , globalsolve=true ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
ridefinisco l'energia potenziale elastica esplicitando 
la dipendenza di ZB da P
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
U :

fullratsimp(

ev ( U )

);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
utilizzo il teorema di castigliano per trovare,
nel quarto di struttura modellato per antisimmetria,
la freccia sotto il carico P
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
dispzC : diff ( U , P );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
ROBA VARIA
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
assegnazione della lista di variabili sulla base della lista
di definizioni valutate secondo ipotesi sols

''unks : ev ( unks , sols ) ;

''unks forza la valutazione di unks, che evidentemente 
non avviene spontaneamente
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
[ pippo , pluto , paperino ] : [ 1 , ciao , miao ];
pippo;
pluto;
paperino;
kill(pippo);
pippo;
   [wxMaxima: comment end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$