kill(all)$globalsolve:true$

Cvičení 10, Příklad 2: nosník na 2 podporách, spojité zatížení A.sin(B.x)

assume(l>0,X>0);

q:q_0*sin(%pi*x/l);

wxplot2d(ev(q,q_0:1,l:1),[x,0,1])$

RR

RR:[
R_A+R_B-integrate(q,x,0,l)=0,
R_B*l-integrate(q*x,x,0,l)=0];

linsolve(RR,[R_A,R_B]);

M:expand(R_A*X-integrate(q*(X-x),x,0,X));

φ:expand(integrate(-M/EJ,X)+C_1);
w:expand(integrate(φ,X)+C_2);

OP:[
at(w,X=0)=0,
at(w,X=l)=0
];

linsolve(OP,[C_1,C_2]);

w:expand(''w);

wxplot2d(ev(diff(w,X,0),q_0:1,l:1,EJ:1,eval),[X,0,1],[ylabel,"w"])$
wxplot2d(ev(diff(w,X,1),q_0:1,l:1,EJ:1,eval),[X,0,1],[ylabel,"φ"])$
wxplot2d(ev(-diff(w,X,2)*EJ,q_0:1,l:1,EJ:1,eval),[X,0,1],[ylabel,"M"])$
wxplot2d(ev(-diff(w,X,3),q_0:1,l:1,EJ:1,eval),[X,0,1],[ylabel,"T"])$
wxplot2d(ev(diff(w,X,4),q_0:1,l:1,EJ:1,eval),[X,0,1],[ylabel,"q"])$