12-2-5.html

Кирсанов М.Н. Решебник.Теоретическая механика.  М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге

12.2. Задача 5

  restart:

               Подписи шарниров на рисунке

  nam:=array(0..3,[O,A,D,C]):

              Считываем библиотеку стандартных программ для рисования

  read «ris.m»;

             Начальные значения и размеры         ( м)

   OA:=6: L:=3: a:=0.15: b:=0.5:
 BC:=10:R1:=1.2:R2:=1:r:=0.1:

    Угловая  скорость  OA и период обращения       

  omegaOA:=1: T:=2*pi/omegaOA:

   with(plots):with(plottools):

                 Количество кадров K

  K:=24:

                Создаем все кадры

  for i from 0 to K do

  t:=evalf(sin(i/K*Pi))*T/6:

          Закон движения  АO

  phi:=pi/4+t*omegaOA:    

                   Координаты  O

  x[0]:=-phi*R1:y[0]:=0:
x[7]:=-phi*R1:y[7]:=R1:         
            Координаты цилиндра D            N2     
  x[2]:=L:    y[2]:=y[0]-R2+R1:
  x[5]:=L:    y[5]:=y[0]+R1:
  x[6]:=L+R2: y[6]:=y[2]:
       Координаты груза C                         N3   

     x[3]:=x[6]:  y[3]:=-x[0]-5:

   x[4]:=x[6]:  y[4]:=y[3]-b:
 x[8]:=-5:  y[8]:=-R1:
 x[9]:=-5:y[9]:=-6:
 x[10]:=0:y[10]:=-R1:
 x[1]:=x[8]-2.4*a: y[1]:=y[0]-sqrt(OA^2-(x[0]-x[1])^2):

   P[i]:=display(box(4,b,b),box(1,2*a,3*a),
 Cir(0,r),Cir(1,r),Cir(0,R1),Cir(2,R2),
 Line(7,5,0),Line(6,3,0),Line(1,0,7),
 cir4(2,R2,-x[0]/R2),cir4(0,R1,-x[0]/R1),
 PLOT(seq(TEXT([x[j],y[j]+0.3],nam[j]),j=0..3))):
od: Конец цикла формирования кадров

                                 Изображение механизма в движении

   PP:=display(seq(P[i],i=0..K),insequence=true,
                              thickness=2,
                              scaling=constrained,
                              axes=none):

  display(Lin(8,9),Lin(10,8),Поверхность(-3,-R1,3*R1,0.3),PP,Опора(2,0.1));