Needs["dsvsolve`"];

Input

γ=Pi/4;
example={$$domains->{{-2*b, -2*b, 2*b, b}, {-b, b, b, 2*b}},
$$coefficients->{1, 1},
$$traction->0,$$wheretraction->{0, 0},
$$bending->{-MM*Cos[γ], -MM*Sin[γ]}}; dSVcShowInput[example]

Risultati simbolici

A=14b ^2;xG=0;yG=((3 b)/14);Iη=((50 (b^4))/3);Iξ=((673 (b^4))/42);a0=-((9 MM)/(673 (Sqrt[2]) (b^3)));a1=((3 MM)/(50 (Sqrt[2]) (b^4)));a2=-((21 (Sqrt[2]) MM)/(673 (b^4)));σ[x_,y_]=a0-a1 x-a2 y;tanβ=-a1/a2;β=ArcTan[tanβ];σmax=σ[2*b,-b];σmin=σ[-2 *b,2*b];

Risultati numerici

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

Output

sol = {AB, rG, JB, sigma, gobj, details} = dSVcSolve[example];
gsol=dSVcShowOutput[example, sol]

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