Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге
12.6. Задача 8
> | restart: |
> | read «ris.m»: |
Подписи шарниров на рисунке
nam:=[A,B,C,D]:Nnam:=[2,4,5,8]:
Начальные значения и размеры ( cм) (радиус ra уменьшен по сравнению с условием для того,
чтобы колесо А не катилось слишком далеко)
rc:=18: Rc:=32: rd:=18: Rd:=25:
Угловое ускорение
> | eps:=5: |
Координаты опор
> | x[7]:=0: y[7]:=0: x[1]:=-Rc: y[1]:=0: x[3]:=Rc: y[3]:=0: x[5]:=0: y[5]:=rc: x[2]:=x[1]: x[4]:=x[3]: y[6]:=y[5]: y[8]:=y[6]-rd: Количество кадров K |
> | K:=12: with(plots):with(plottools): |
Создаем все кадры
> | for i from 0 to K do |
> | t:=i/K: |
> | phi:=t^2/2*eps: S12:=40+phi*Rc:S34:=120-phi*Rc: S56:=133-phi*rc/(Rd+rd)*Rd: y[2]:=-S12: y[4]:=-S34: x[6]:=S56: x[8]:=x[6]: P[i]:=display(Cir(7,Rc),Cir(7,rc), Cir(8,1), box(2,8,8),box(4,8,8), Cir(8,Rd),Cir(8,rd), cir4(7,Rc,phi),cir4(7,rc,phi), cir4(8,Rd,phi*rc/(Rd+rd)),cir4(8,rd,phi*rc/(Rd+rd)), seq(TEXT([x[Nnam[j]]+14,y[Nnam[j]]+6],nam[j]),j=1..4), seq(Line(2*i-1,2*i,9),i=1..3)): 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,Опора(7,2),Поверхность(80,y[8]-Rd,60,5)); |
> |