D:\j_Iam_TEMP\_jh\_jh_ex2.wxmx

# Умова

Визначити силу гідростатичного тиску на бокову грань "MN" параболоїдального каналу
висотою Hmn = 1 метр та шириною B = 3 метри. Рівень води у каналі H = 2 метри.

kill(all);

ro:1000;g:9.81;H:2;B:3;Hmn:1;

fh(x):=x**2;fHab(x):=Hmn;fH(x):=H;Fx(h):=sqrt(h);diff(Fx(h),h);

dh(h):=1/(2*sqrt(h));

plot2d([fh(x),fHab(x),fH(x)],[x,0,2]);

D:\j_Iam_TEMP\_jh\_jh_ex2.png

plot2d(Fx(h),[h,0,H]);

D:\j_Iam_TEMP\_jh\_jh_ex2a.png

fp(h):=ro*g*(H-h);

# P_x - горизонтальна проекція сили тиску

P_x:integrate(fp(h)*B,h,0,Hmn);

mP_x:integrate(fp(h)*(H-h)*B*h,h,0,Hmn);

mP_x_down:integrate(fp(h)*B*h,h,0,Hmn);

h_D:mP_x/P_x;

h_D_:mP_x_down/P_x;

h_D+h_D_;

D:\j_Iam_TEMP\_jh\_jh_ex2_P_x.jpg

# P_z - вертикальна проекція сили тиску

P_z:integrate(fp(h)*B*dh(h),h,0,Hmn);

mP_z:integrate(fp(h)*B*dh(h)*(Fx(h)),h,0,Hmn);

x_C:mP_z/P_z;

D:\j_Iam_TEMP\_jh\_jh_ex2_P_z.jpg

# P - Сила гідростатичного тиску

P:sqrt(P_x**2+P_z**2);

phi_rad:atan(P_z/P_x),numer;

phi_grad:atan(P_z/P_x)*(180/%pi),numer;

d:tan(P_z/P_x);

Plot2d

fh_D(x):=h_D_;

plot2d([fh(x),fHab(x),fH(x),fh_D(x),[discrete,[0.44,0.44],[0,1]],
        [discrete,[0.44],[0.45]]],[x,0,1.5],
    [legend, "parabola","Top_box", "Water", "P_x", "P_z", "point"], 
    [style, [lines, 5,5], lines, [lines, 3,1], lines, lines, [points, 3,2]], 
    [point_type, circle]);

D:\j_Iam_TEMP\_jh\_jh_ex2c.png

eq:y=x_C;plot2d([fh(x),fHab(x),fH(x),fh_D(x),eq],[x,0,2],[y,0,4]);

# Answer

P_x;P_z;P;h_D_;x_C;phi_grad;

D:\j_Iam_TEMP\_jh\_jh_ex2.jpg

Scale screen 300 px

P_x/300;P_z/300;P/300;