Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге
12.2. Задача 9
restart:read «ris.m»;
Подписи шарниров на рисунке
nam:=array(0..1,[O,A]):
Начальные значения и размеры ( м)
R:=1: R2:=0.7: r:=0.1:a:=4:b:=3: c:=0.2:
Угловая скорость OA и период обращения
omegaOA:=2: T:=2*pi/omegaOA:
with(plots):with(plottools):
Количество кадров K
K:=36:
Создаем все кадры
for i from 0 to K do
t:=sin(pi*i/K)*T/3:
Закон движения АO
phi:=pi/6+t*omegaOA:
Координаты шарнира А N1
x[1]:= 2*R:
Координаты центра цилиндра О N0
x[0]:=x[1]-R2*cos(phi): y[0]:= R:
Координаты шарнира А N1
y[1]:=y[0]+R2*sin(phi):
Горизонтальный шток
x[3]:=x[0]-a: y[3]:=0:
x[4]:=x[3]+2*a: y[4]:=0:
Вертикальный шток
x[2]:=x[1]: y[2]:=y[1]+b:
P[i]:=display(Cir(0,R),
cir4(0,R,phi),
Lin(3,4),Lin(2,1),
Cir(1,r),Cir(0,r),
seq(TEXT([x[j]+0.2,y[j]+0.2],nam[j]),j=0..1)):od:
Изображение механизма в движении
PP:=display(seq(P[i],i=0..K),insequence=true,
thickness=2,
scaling=constrained,
axes=none):
display(Поверхность(-1.5*R,-r,R/2,0.1),
Поверхность(-1.5*R,r,R/2,-0.1),
Поверхность(4*R,r,R/2,-0.1),
Поверхность(4*R,-r,R/2,0.1),
Стенка(2*R-r,3*R,R/2,0.1),
Стенка(2*R+r,3*R,R/2,-0.1),PP);
