12-2-1.html | AcademiaXXI

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

12.2. Задача 1

restart:read «ris.m»;

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

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

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

OA:=1.5: AB:=OA: BC:=7:R:=1:r:=0.1:a:=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:=12:

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

for i from 0 to K do

t:=sin(pi*i/K)*T/12:

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

phi:=pi/6+t*omegaOA:

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

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

Координаты шарнира B N2

x[2]:=2*x[1]: y[2]:=0:

Координаты шарнира C N3

x[3]:=x[2]+BC: y[3]:=0:

Цилиндр:=circle([x[3],y[3]],R,color=red):

P[i]:=display(box(2,4*a,2*a),
cir4(3,R,-x[2]/R),
Цилиндр,Line(0,1,7),Line(1,2,7),Line(2,3,8),Cir(2,r),Cir(1,r),Cir(3,r),seq(TEXT([x[j]+0.1,y[j]+0.4],nam[j]),j=0..3)):od:

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

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

display(Поверхность(x[0]+5*R,-R,6*R,0.3),
Поверхность(1.1,y[2]-2.5*a,2*R,0.25),
PP,Опора(0,0.1));

>