Введение в Maple и рассмотрение задач теоретической механики. Эксаревская М.Е - 30 стр.

UptoLike

30
> W['B'] := W[`A_вр`] + W[`B_вр`]*sin(psi);
3
3560
280:W
B
+=
> w2 := vector([W['B']/W_размер, 0]):
w3 := vector([W[`B_вр`]*sin(psi)/W_размер,
W[`B_вр`]*cos(psi)/W_размер]):
> Wb := arrow(pointB, w2, .35, 1, .20, color=blue):
Wb_вр := arrow(pointB, w3, .35, 1, .20, color=blue):
> display(draw(Pi/2),
Wa, Wb, Wb_вр,
textplot([w1['x'] + dxy, 35 + w1['y'] + dxy, 'Wa'],
color=blue),
textplot([pointB[1] + w2[1] + dxy, dxy, 'Wb'],
color=blue),
textplot([pointB[1] + w3[1] + dxy,
pointB[2] + w3[2] + dxy, 'Wb_вр'], color=blue),
scaling=CONSTRAINED);
> display(draw(Pi/2),
Wa,
arrow(pointA, w2, .35, 1, .20, color=blue),
arrow([pointA[1] + w1['x'], pointA[2] + w1['y']],
w3, .35, 1, .20, color=blue),
textplot([w1['x'], 35 + w1['y'] - dxy, 'Wa'], color=blue),
textplot([pointA[1] + w2[1] + dxy, pointA[2], 'Wb'],
color=blue),
textplot([25, 20, 'Wb_вр'], color=blue),
scaling=CONSTRAINED);
Найдём ускорение точки .
C Воспользуемся формулой
() ()
ц
AC
вр
ACAC
wwww ++=
и запишем её в проекциях на координатные оси.
> W['C_x'] := W[`A_вр`] + epsilon['AB']*AC*sin(psi):
> W['C_y'] := -W[`A_ц`] + epsilon['AB']*AC*cos(psi):
> W['C'] := sqrt(W['C_x']^2 + W['C_y']^2);
()
400003120280:W
2
C
++=
> w4 := vector([W['C_x']/W_размер, W['C_y']/W_размер]):
> W['B'] := W[`A_вр`] + W[`B_вр`]*sin(psi);
                  560 3
    WB := 280 +
                    3
> w2 := vector([W['B']/W_размер, 0]):
  w3 := vector([W[`B_вр`]*sin(psi)/W_размер,
                W[`B_вр`]*cos(psi)/W_размер]):
> Wb := arrow(pointB, w2, .35, 1, .20, color=blue):
  Wb_вр := arrow(pointB, w3, .35, 1, .20, color=blue):
> display(draw(Pi/2),
          Wa, Wb, Wb_вр,
          textplot([w1['x'] + dxy, 35 + w1['y'] + dxy, 'Wa'],
                   color=blue),
          textplot([pointB[1] + w2[1] + dxy, dxy, 'Wb'],
                   color=blue),
          textplot([pointB[1] + w3[1] + dxy,
                    pointB[2] + w3[2] + dxy, 'Wb_вр'], color=blue),
          scaling=CONSTRAINED);




> display(draw(Pi/2),
          Wa,
          arrow(pointA, w2, .35, 1, .20, color=blue),
          arrow([pointA[1] + w1['x'], pointA[2] + w1['y']],
                w3, .35, 1, .20, color=blue),
          textplot([w1['x'], 35 + w1['y'] - dxy, 'Wa'], color=blue),
          textplot([pointA[1] + w2[1] + dxy, pointA[2], 'Wb'],
                   color=blue),
          textplot([25, 20, 'Wb_вр'], color=blue),
          scaling=CONSTRAINED);




Найдём ускорение точки C. Воспользуемся формулой
                   вр     ц
     wC = w A + wC ( A ) + wC ( A )
и запишем её в проекциях на координатные оси.
> W['C_x'] := W[`A_вр`] + epsilon['AB']*AC*sin(psi):
> W['C_y'] := -W[`A_ц`] + epsilon['AB']*AC*cos(psi):
> W['C'] := sqrt(W['C_x']^2 + W['C_y']^2);

    WC :=   (280 + 120 3 )2
                              + 40000
> w4 := vector([W['C_x']/W_размер, W['C_y']/W_размер]):

                                        30