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

13.8. Пример 2

restart:

with(plots):with(plottools):

    Количество кадров
  K:=12:pi:=evalf(Pi):
  r1:=10: R1:=20:

R2:=20: OA:=60: OB:=OA-r1-R2:

Опора:=cylinder([0,0,0],OA+R1,8,capped=false,style=PATCHNOGRID),        cuboid([-3,-80,-3],[3,80,0]),        cylinder([0,0,-4],1,10):

Водило:=cuboid([-3,-3,0],[65,3,2]):

>	Блок:=cylinder([0,0,2],R1,4),       cylinder([0,0,6],r1,6),       cylinder([0,0,-2],1,16): Колесо:=cylinder([0,0,7],R2,6),cylinder([0,0,-2],1,16): for i to K do Амплитуда:=pi/15:

>	phi:=sin(2*pi*i/K)*Амплитуда:

>	P1:=display(rotate(PLOT3D(Блок),0,0,3*phi)):

>	P2:=display(rotate(PLOT3D(Колесо),0,0,-3*phi)):

>	P3:=display(rotate(PLOT3D(Водило),0,0,-phi)):

>	x2:= OA*cos(phi): y2:= OA*sin(phi):

>	x3:= OB*cos(phi): y3:= OB*sin(phi): P[i]:=display(translate(P1,x2,y2,0), translate(P2,x3,y3,0),P3);

>

od:

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

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

>	display(PP,Опора);