Clear the memory, or, kill all the variable definitions

kill( all ) ;

equilibrium equations for the quarter
(simplified) ladder frame chassis

eqtx : UA = 0 ;

eqty : VD = 0 ;

eqtz : WA + WD - F = 0;

eqrOx : - F * c + WD*(b+e) + ThetaD = 0;

eqrOy : PhiA - WA*a + F*a = 0;

eqrOz : 0 = 0 ;

define a list of equations and a list of unknowns

eqlst : [ eqtx  , eqty  , eqtz ,
          eqrOx , eqrOy , eqrOz ] ;

unlst : [ UA , WA , PhiA , VD , WD, ThetaD ];

solve the linear system of equations

linsolve( eqlst , unlst ) ;

there is one dependent equation, so I'm unable to solve 
the problem of obtaining all the reaction forces/moments by equilibrium alone

unlst_red : [ UA , WA , PhiA , VD , WD];

sols : linsolve( eqlst , unlst_red ) ;

option A for assignment:

''unlst_red : ev( unlst_red , sols ) ;

massive assignment of the values that solve the equations to the 
unknowns. The ''unlst_red syntax forces an immediate evaluation of the
symbol, such that the actual assignment is

[UA,WA,PhiA,VD,WD] : [UA=0,WA=((e-c+b)*F+ThetaD)/(e+b),
                      PhiA=-(a*c*F-a*ThetaD)/(e+b),VD=0,
                      WD=(c*F-ThetaD)/(e+b)];

option B for assignment

linsolve( eqlst , unlst_red ) , globalsolve=True ;

we now define the cross section resultants for the D->B beam, 
as a function of the s= 0..a arclength coordinate

N__DB : 0 $

Mf_eta__DB : VD * s $

Mf_xi__DB : WD * s ;

Mt__DB : ThetaD;

S_eta__DB : WD;

S_xi__DB  : VD$

we now evaluate the internal strain energy for the DB beam

U__DB : integrate ( 
               ( 
                       J_etaeta__DB* M_xi__DB^2
                   +    J_xixi__DB*M_eta__DB^2
                   + 2*J_xieta__DB* M_xi__DB*M_eta__DB
                )/( 
                    2*E__DB*(J_xixi__DB*J_etaeta__DB-J_xieta__DB^2)
                )
            +
                N__DB^2/2/E__DB/A__DB
            +
                ( 
                      chi_xi__DB * S_xi__DB^2
                    + chi_eta__DB * S_eta__DB^2
                    + chi_etaxi__DB * S_xi__DB * S_eta__DB
                )/(
                    2*G__DB*A__DB
                )
            +
                Mt__DB^2 / 2 / G__DB/ Kt__DB
 , s , 0 , a );