Needs["dsvsolve`"];

Input

example={$$nodes->{{-L, L}, {0, L}, {0, 0}, {L, 0}},
$$edges->{{1, 2}, {2, 3}, {2, 4}},$$absname->s,
$$constraints->{{{"rollv", 1}}, {{"clamp", 1, 2}, {"clamp", 2, 3}}, {{"rollh", 2}}, {{"rollv", 3}}},
$$bodyloads->{{1, {0, -q, 0}}},
$$nodalloads->{{3, {q*L, 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];