Needs["dsvsolve`"];

Input

example={$$nodes->{{-b, b/2}, {-b, -1/2*b}, {-b, (-5*b)/2}, {-b, (5*b)/2}, {b, b/2}, {b, (5*b)/2}, {b, -1/2*b}, {b, (-5*b)/2}},
$$edges->{{4, 1}, {1, 2}, {2, 3}, {2, 7}, {7, 8}, {7, 5}, {5, 6}, {1, 5}},
$$thicknesses->s,$$absname->t,$$thname->z,
$$traction->F,$$wheretraction->{b, -b/2},
$$shear->{0, 0}, $$whereshear->{0, 0},
$$bending->{0,0}, $$torsion->0}; dSVtShowInput[example]

Output simbolico

sol = {AB, rG, JB, CS, FS, KT, sigma, tau, gobj, details} = dSVtSolve[example];
gsol=dSVtShowOutput[example, sol]

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

Output

Output simbolico

sol = {AB, rG, JB, CS, FS, KT, sigma, tau, gobj, details} = dSVtSolve[example];
gsol=dSVtShowOutput[example, sol]

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

Caratteristiche geometriche della sezione

valori={F->Quantity[2*10^4,"Kilonewtons"],b->Quantity[14,"Centimeters"],s->Quantity[2,"Centimeters"],σ0->Quantity[20,"Kilonewtons"/"Centimeters"^2]};
A=Last[Last[sol][[4]]];
{{Iy,Iyx},{Ixy,Ix}}=Last[Last[sol][[10]]];
Print["Area = ",A,"=",A/.valori//N//(NumberForm[#1,3]&)];
Print["Ix = ",Ix,"=",Ix/.valori//N//(NumberForm[#1,3]&)];
Print["Iy = ",Iy,"=",Iy/.valori//N//(NumberForm[#1,3]&)];
Print["Ixy = ",Ixy,"=",Ixy/.valori//N//(NumberForm[#1,3]&)];

Tensione normale

{a0,a1,a2}=Last[Last[sol][[16]]];
σ[x_,y_]=a0+-a1*x-a2*y;

Print["(Subscript[σ, z])(x,y)=",σ[x,y],"=",σ[x,y]/.valori//N//(NumberForm[#1,3]&)]

Punti e valori di massima sollecitazione normale

xmax=b;ymax=-5 b/2;
xmin=-b;ymin=5/2b;

σmax=σ[xmax,ymax];
σmin=σ[xmin,ymin];

Print["(Subscript[σ, max])=",σmax,"=",σmax/.valori//N//(NumberForm[#1,3]&)]
Print["(Subscript[σ, min])=",σmin,"=",σmin/.valori//N//(NumberForm[#1,3]&)]