Finite Element solution using BEAM2D elements

parameters={
l->1.3,
d1->0.04,
d2->0.08,
p->700,
M1-> 300,
F1->900,
E1->215
};

Kbeam=((2 Iz E1)/(L^3)){{6,3L,-6,3L},
{3L,2(L^2),-3L,(L^2)},
{-6,-3L,6,-3L},
{3L,(L^2),-3L,2(L^2)}};
Kbeam//MatrixForm//TraditionalForm

Fbeam={Vi,Mi,Vj,Mj}+((p L)/12){6,L,6,-L};
Fbeam//MatrixForm//TraditionalForm

Stiffness mx and load vector for each elements

K1=Kbeam/.{L->1.5l,Iz->(d2^4)π/64}/.parameters;
K1//MatrixForm

Global stiffness mx and load vector

...

Condensed stiffness mx and solution for unconstrained DoFs