Needs["dsvsolve`"];

Input

example={$$nodes->{{-L, 0}, {0, 0}},
$$edges->{{1, 2}},$$absname->s,
$$constraints->{{{"hinge", 1}}, {{"hinge", 1}}},
$$bodyloads->{{1, {0, -(q*s)/L, 0}}},
$$nodalloads->{{1, {0, 0, 0}}},
$$predeformations->{{1, {0, 0, 0}}},
$$cdisplacements->{{1, {0, 0, 0}}},
$$EIoverEAL2->0};
(*All constraints are initialized to hinge constraints; use "BFConstraintsPalette[BFConstraintTypes];" to input different constraints.*)
BFShowInput[example]

Output

(*BFClassify[example]; sol = {{Ns, Qs, Ms}, details} = BFForcesSolve[example];*)
(*A=Transpose[First[BFStaticProblem[example]]];MatrixForm[A]
BFShowRigidMotions[example]*)
BFEquations[example,"Nodal"]; sol = {{us, vs}, {Ns, Qs, Ms}, details} = BFDisplacementsSolve[example];
gsol=BFShowOutput[example, sol]

CellPrint/@Rest/@Last[sol];
Print/@Cases[gsol,Graphics[___],Infinity];

-D[(((L^3) q (7 s-10 (s^3)+3 (s^5)))/(360 EI)),s]/.s->1