12-2-4.html | AcademiaXXI

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

12.2. Задача 4

restart:

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

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

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

read «ris.m»;

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

pi:=evalf(Pi):

OA:=2: 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):

Warning, the name changecoords has been redefined

Warning, the name arrow has been redefined

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

K:=36:

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

for i from 0 to K do

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

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

phi:=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: y[2]:=y[0]-R2:
  x[5]:=L: y[5]:=y[0]:
  x[6]:=L+R2: y[6]:=y[2]:
      Координаты груза C                         N3

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

x[4]:=x[6]: y[4]:=y[3]-b:

P[i]:=display(box(4,b,b),
Cir(0,r),Cir(1,r),Cir(0,R1),Cir(2,R2),
Line(0,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]-0.3,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(Поверхность(-9,0,7*R1,0.4,1),PP,Опора(2,0.1));