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

UptoLike

32
4.2. Сложное движение точки
> restart:
Точка
M
движется относительно тела D .
По заданным уравнениям относительного движения точки
M
() ( )()
tπ2sin110tSOM +==
и движения тела
()
2
t6.1t4tφφ +==
определить для заданного момента времени
8
1
t =
абсолютную скорость и
абсолютное ускорение точки .
M
Схема механизма показана на рисунке.
Решение
> S := t -> 10 * (1 + sin(2*Pi*t)):
phi := t -> 4*t + 1.6*t^2:
> with(plots):
with(plottools):
> pos := proc(psi, s, t1)
local
point_x, point_y;
point_x := s*cos(psi);
point_y := s*sin(psi);
PLOT(plottools[point]([point_x, point_y], symbol=circle,
color=red, symbolsize=15),
op(plot([S(t), phi(t), t=0 .. t1], coords=polar,
color=green)),
line([0, 0], [25*cos(psi), 25*sin(psi)], color=red,
linestyle=4));
end proc:
> pos2 := proc(psi, s, t1)
local
телоD, телоB,
pointB, point1, point2, point3, point4,
h, OA;
OA := 30;
h := 0.25;
point1 := [h*cos(psi+Pi/2),h*sin(psi+Pi/2)];
point2 := [h*cos(psi+Pi/2) + OA*cos(psi),
h*sin(psi+Pi/2) + OA*sin(psi)];
point3 := [point2[1] + 2*h*cos(psi+3*Pi/2),
point2[2] + 2*h*sin(psi+3*Pi/2)];
point4 := [h*cos(psi+3*Pi/2), h*sin(psi+3*Pi/2)];
     4.2. Сложное движение точки
> restart:
Точка M движется относительно тела D .
По заданным            уравнениям относительного движения точки M
     OM = S (t ) = 10 (1 + sin(2πt ))
и движения тела
     φ = φ(t ) = 4 t + 1.6 t 2
                                               1
определить для заданного момента времени t = абсолютную скорость и
                                               8
абсолютное ускорение точки M .
Схема механизма показана на рисунке.




Решение
> S := t -> 10 * (1 + sin(2*Pi*t)):
  phi := t -> 4*t + 1.6*t^2:
> with(plots):
  with(plottools):
> pos := proc(psi, s, t1)
    local
      point_x, point_y;
    point_x := s*cos(psi);
    point_y := s*sin(psi);
    PLOT(plottools[point]([point_x, point_y], symbol=circle,
                            color=red, symbolsize=15),
          op(plot([S(t), phi(t), t=0 .. t1], coords=polar,
             color=green)),
          line([0, 0], [25*cos(psi), 25*sin(psi)], color=red,
               linestyle=4));
  end proc:
> pos2 := proc(psi, s, t1)
    local
      телоD, телоB,
      pointB, point1, point2, point3, point4,
      h, OA;
    OA := 30;
    h := 0.25;
    point1 := [h*cos(psi+Pi/2),h*sin(psi+Pi/2)];
    point2 := [h*cos(psi+Pi/2) + OA*cos(psi),
                h*sin(psi+Pi/2) + OA*sin(psi)];
    point3 := [point2[1] + 2*h*cos(psi+3*Pi/2),
                point2[2] + 2*h*sin(psi+3*Pi/2)];
    point4 := [h*cos(psi+3*Pi/2), h*sin(psi+3*Pi/2)];
                                  32