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->0,$$wheretraction->{0, 0},
$$shear->{0, 0}, $$whereshear->{0, 0},
$$bending->{-MM*Cos[γ], -MM*Sin[γ]}, $$torsion->0}; dSVtShowInput[example]

Risultati

A=14b s;Iy=((34 (b^3) s)/3);Ix=((131 (b^3) s)/6);Mx=MM Cos[γ];My=MM Sin[γ];σ[x_,y_]=Mx/Ix*y-My/Iy*x;tanβ=Tan[γ]*Ix/Iy;β=ArcTan[tanβ];σmax=σ[-b,5/2b];σmin=σ[b,-5/2 b];

Versione 1

valori1={γ->20*Pi/180,MM->Quantity[2*10^4,"Kilonewtons"*"Centimeters"],b->Quantity[12,"Centimeters"],s->Quantity[1,"Centimeters"],σ0->Quantity[20,"Kilonewtons"/"Centimeters"^2]};
risultati1={A,Ix,Iy,β*180/Pi,σmax,σmin,σmax/σ0,σmin/σ0}/.valori1//N//NumberForm[#,3]&

Versione 2

valori2={γ->20*Pi/180,MM->Quantity[10^4,"Kilonewtons"*"Centimeters"],b->Quantity[12,"Centimeters"],s->Quantity[1,"Centimeters"],σ0->Quantity[20,"Kilonewtons"/"Centimeters"^2]};
risultati2={A,Ix,Iy,β*180/Pi,σmax,σmin,σmax/σ0,σmin/σ0}/.valori2//N//NumberForm[#,3]&

Versione 3

valori3={γ->80*Pi/180,MM->Quantity[5*10^4,"KiMlonewtons"*"Centimeters"],b->Quantity[14,"Centimeters"],s->Quantity[1,"Centimeters"],σ0->Quantity[20,"Kilonewtons"/"Centimeters"^2]};
risultati3={A,Ix,Iy,β*180/Pi,σmax,σmin,σmax/σ0,σmin/σ0}/.valori2//N//NumberForm[#,3]&

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 numerico versione 1

γ=20*Pi/180;MM=2*10^4;b=12;s=1;
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 numerico versione 2

γ=20*Pi/180;MM=10^4;b=12;s=1;
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 numerico versione 3

γ=80*Pi/180;MM=3*10^4;b=14;s=1;
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