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

> display(draw1(Pi/2),
          Wa, Wax, Way,
          Wb, Wb_вр, Wb_ц,
          textplot([w1['x'] + dxy, 35 + w1['y'] + dxy, 'Wa'],
                   color=blue),
          textplot([w1['x'] + 1.5*dxy, 35, 'Wa_вр'], color=blue),
          textplot([1.5*dxy, 35 + w1['y'] + dxy, 'Wa_ц'],
                   color=blue),
          textplot([pointB[1] + 20 + dxy, dxy, 'Wb'], color=blue),
          textplot([pointB[1] + 20*sin(psi) + dxy,
                    pointB[2] + 20*cos(psi) + dxy, 'Wb_вр'],
                    color=blue),
          textplot([pointB[1] - 20*cos(psi) + dxy,
                    pointB[2] + 20*sin(psi) + dxy, 'Wb_ц'],
                    color=blue),
          scaling=CONSTRAINED);




На рисунке показаны:
   •   направление ускорения точки B (в предположении замедленного
       движения ползуна);
   •   направление вращательного ускорения точки B вокруг точки A (в
       предположении о положительном значении ε AB );
   •   направление центростремительного ускорения точки B (в действи-
       тельности оно равно нулю, т.к. ω AB = 0 ).
      Величины указанных ускорений найдем из проекций уравнения (*)
на оси координат. Если в результате получим отрицательные значения, то
направления ускорений надо изменить на противоположные.
      Запишем уравнение (*) в проекции на ось Oy .
> -W[`A_ц`] + W[`B_вр`]*cos(psi) = 0:
  W[`B_вр`] := solve(%, W[`B_вр`]);
            1120 3
    WB _ вр :=
               3
     Так как вращательное ускорение получилось положительным, то на-
правление на рисунке указано верно. Поскольку WB _ вр = ε AB AB , то мы мо-
жем найти угловое ускорение шатуна АВ:
> epsilon['AB'] := W[`B_вр`] / AB;
         16 3
    ε AB :=
           3
Запишем уравнение (*) в проекции на ось Ox :
                                    29