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

rovnice rovnováhy tuhého tělesa

RR:[
-N_1+R_x+N_3*sin(β)=0,      /*směr x*/
N_2+R_y+N_3*cos(β)=0,       /*směr y*/
N_2*4*a+R_y*2*a-R_x*a=0     /*kolem z*/
];

fyzikální rovnice pro jednotlivé deformovatelné pruty

Δ_1:N_1*l/E/S+α*ΔT*l;
Δ_2:N_2*l/E/S;
Δ_3:N_3*l/E/S;

deformační podmínky

DP:[
'(Δ_1=δ+a*φ),
'(Δ_2=-2*a*φ),
'(Δ_3=a*φ*(2*cos(β)-sin(β)))
];

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

DP:ratsimp(''DP);

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

linsolve(append(RR,DP),[φ,R_x,R_y,N_1,N_2,N_3])$

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);
'(N_3)=factorsum(N_3);