kill(all);

load(scifac);

"-----------------------------------Expression de la charge repartie"

define(q(x),q0);

"-----------------------------------Effort total du a  la charge repartie, selon y"

define(r(x),integrate(-q(t),t,0,x));

"------------------------------------Au total"

r(l);

"-------------------------------------Moment en A du à la charge repartie, selon z"

combine(define(mr(x),integrate(-t*q(t),t,0,x)));

"------------------------------------pour x = l on aura"

combine(mr(l));

"--------------------------Mise en place et resolution du systeme lineaire a  2 equations 2 incoonnues"

S1:[Fya+Fyb+r(l)=0, Fya*l+2*l*Fyb+mr(l)=0];

solve(S1,[Fya, Fyb]);

Fyb:-l*q0/2;

"-------------------------------------Calcul du moment du torseur de cohesion sur [OA] en G"
"Moment du a  la charge repartie"

define(m(x),integrate(t*q(t),t,0,x));

"Moment total du torseur de cohesion"

define(mt(x),-m(x));

"--------------------------------------Calcul du moment du torseur de cohesion sur [AB] en G"
"On somme tous les moments a  droite de G. Il ne reste que Fyb pour S1 !"

define(mts(x),(2*l-x)*Fyb);

"------------------------------------On precise les constantes : q0=1N/m et l=8m"

q0:1;

l:8;

"-------------------------Verification de la continuite du moment en B"

mt(l);

mts(l);

"----------------------------On dessine le diagramme du moment flechissant :"

H(x) := if x > l then 0 else mt(x) $

G(x) := if x < l then 0 else mts(x)$

 plot2d ([H(x), G(x)], [x,0,2*l])$

kill(all);