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

13.5. Пример 2

>

restart:

>	read "ris.m":

               Подписи точек на рисунке
nam:=["O1",A,B,D,"O3",C]:Nnam:=[1,2,3,4,11,7]:
R1:=40: r1:=20: R3:=30: r4:=20: R4:=41:a:=pi/9: h:=99:
stok:=95:
b:=1.5:sna:=sin(a): csa:=cos(a):

             Угловая  скорость       

>

omeg:=2:

                     Координаты

>

x[11]:=0:        y[11]:=0:   x[4]:=0:          y[4]:=-R3:   x[6]:=-R3*csa:    y[6]:=R3*sna:   x[5]:=x[6]-h*sna: y[5]:=y[6]-h*csa:   x[1]:=x[5]-R1*csa: y[1]:=y[5]+R1*sna:   y[12]:=y[1]:   y[7]:=y[12]:  y[8]:=y[12]:y[10]:=y[12]-r4:   y[9]:=y[10]:  x[3]:=x[1]: y[3]:=y[1]+R1:             Количество кадров K

>	K:=12: with(plots):with(plottools):

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

>	for i from 0 to K do

>	t:=2*pi*i/K:  phi:=t*omeg+3*pi/4: F:=arrow([x[9]-50,y[9]],evalm([40,0]),4,12,0.2,color=red): x[2]:=x[1]+r1*cos(phi):y[2]:=y[1]+r1*sin(phi):

>	x[12]:=x[2]:x[8]:=x[12]-b: x[7]:=x[8]-1.5*stok:x[9]:=x[2]*(R4-r4)/R4-3*stok: x[10]:=x[9]+2*stok:

>	P[i]:=display(Cir(1,R1),Cir(11,R3),F, Cir(7,R4),Cir(7,r4),Cir(7,1), Cir(2,1),cir4(7,R4,-x[2]/R4), cir4(7,r4,-x[2]/R4), cir4(1,R1,phi),cir4(11,R3,-phi*R1/R3), box(12,2*b,0.7*R1),box(9,10,R4-r4), seq(Line(2*i-1,2*i,9),i=1..4),Lin(9,10), seq(TEXT([x[Nnam[j]]+2,y[Nnam[j]]+8],nam[j]),j=1..6)): od:

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

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

Warning, the name changecoords has been redefined

Warning, the name arrow has been redefined

>	display(PP,Опора(1,2),Опора(11,2),Поверхность(-360,y[7]-R4-1,160,6), Поверхность(-180,y[7]-r4-3,20,3), Поверхность(-180,y[7]-r4+3,20,-3), Поверхность(-180,y[7]-3,20,3), Поверхность(-180,y[7]+3,20,-3));

>

>