12-2-10.html

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

12.2. Задача 10

  restart:

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

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

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

  read «ris.m»;

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

  OA:=2: L:=3: a:=0.15: 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/3:

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

  phi:=3*pi/4+t*omegaOA:    

                   Координаты  O

  x[0]:=-phi*R1:y[0]:=R1:     

                 Координаты  А                         N1

  x[1]:=x[0]+OA*cos(phi): y[1]:=y[0]+OA*sin(phi):

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

  x[3]:=x[2]+4: y[3]:=y[2]+R2:

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

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

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

  display(Поверхность(-7,0,5,0.3,1),
Поверхность(-1,R1-2*R2,4,0.3,1),
Поверхность(3,R1-2.5*a,3,0.3,1),PP);