kill(all)$globalsolve:true$ /*vyčistíme si paměť a nastavíme řešič*/

rovnice rovnováhy tuhého tělesa

RR:[
R_x+N_2*√(2)/2=0,           /*směr x*/
R_y+N_1+N_2*√(2)/2-F=0,     /*směr y*/
R_y*2*b+N_1*b+F*b=0         /*kolem z*/
];

fyzikální rovnice pro jednotlivé deformovatelné pruty

Δ_1:N_1*b/E/S;
Δ_2:N_2*b*√(2)/E/S;

deformační podmínka

DP:'((Δ_1-δ)/b=(2*Δ_2/√(2))/(2*b));

dosazením z FR do DP za jednotlivá Δ dostaneme:

DP:[ratsimp(''DP)];

řešení soustavy RR (3 rovnice) a DP (1 rovnice) pro neznámé reakce Rx a Ry a síly v prutech N1 a N2:

linsolve(append(RR,DP),[R_x,R_y,N_1,N_2]);

vyjádření řešení v jednotlivých prutech

'(N_1)=factorsum(N_1);             /*zápis z pohledu Maximy čuňárna, ale vypadá to hezky...*/
'(N_2)=factorsum(N_2);