Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге
12.2. Задача 2
restart:
Подписи шарниров на рисунке
nam:=array(0..4,[O,A,B,C,D]):
Считываем библиотеку стандартных программ для рисования
read «ris.m»;
Начальные значения и размеры ( м)
OA:=2: a:=0.15: b:=3: BC:=10:R:=1:r:=0.1:
Угловая скорость OA и период обращения
omegaOA:=2: T:=2*pi/omegaOA:
Координаты опоры O
x[0]:=0:y[0]:=0:
with(plots):with(plottools):
Warning, the name changecoords has been redefined
Warning, the name arrow has been redefined
Количество кадров K
K:=24:
Создаем все кадры
for i from 0 to K do
t:=i/K*T:
Закон движения АO
phi:=pi/4+t*omegaOA:
Координаты шарнира А N1
x[1]:= OA*cos(phi): y[1]:= OA*sin(phi):
Координаты ползуна C N3
x[3]:=x[2]+BC: y[3]:=0:
Координаты B N2
x[2]:=x[1]+a/2: y[2]:=0:
Координаты цилиндра D N4
x[4]:=x[1]/2+BC-3*R: y[4]:=R:
Координаты центра кулисы N5
x[5]:=x[1]: y[5]:=0:P[i]:=display(box(5,a,b),box(3,3*a,2*a),Cir(4,R),Line(0,1,8),
Line(3,2,0),Cir(1,r),Cir(4,r),cir4(4,R,x[1]/R/2),
seq(TEXT([x[j]-0.3,y[j]+0.7],nam[j]),j=0..4)): od:
Изображение механизма в движении
PP:=display(seq(P[i],i=0..K),insequence=true,
thickness=2,
scaling=constrained,
axes=none):
display(Поверхность(x[0]+5*R,2*R,5*R,-0.3,1),
Поверхность(x[0]+7*R,-2*a,5*R,0.3,1),PP,Опора(0,0.1));