Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге
12.5. Задача 6
| > | restart: |
| > | read «ris.m»: |
Подписи шарниров на рисунке
nam:=[B,A,C,D]:Nnam:=[10,9,5,6]:
Начальные значения и размеры ( cм) (радиус ra уменьшен по сравнению с условием для того,
чтобы колесо А не катилось слишком далеко)
Rb:=25: rc:=24: Rc:=42: ra:=12: Ra:=25: a:=pi/4:h:=150:
sna:=sin(a): csa:=cos(a):
Угловое ускорение
| > | eps:=2.5: |
Координаты опор
| > | x[10]:=0: y[10]:=0: x[1]:=Rb*sna: y[1]:=-Rb*csa: x[3]:=0: y[3]:=Rb: y[4]:=y[3]: y[5]:=y[4]+Rc: x[11]:=60: y[11]:=y[5]: x[7]:=x[10]+(Rb+Ra-ra)*sna-h*csa:y[7]:=y[10]-(Rb+Ra-ra)*csa-h*sna: x[8]:=x[7]+0.7*h*csa:y[8]:=y[7]+0.7*h*sna: Количество кадров K |
| > | K:=12: with(plots):with(plottools): |
Создаем все кадры
| > | for i from 0 to K do |
| > | t:=i/K: |
| > | phi:=t^2/2*eps: S12:=55+phi*Rb/(Ra-ra)*Ra:S34:=40+phi*Rb/(Rc-rc)*rc: x[2]:=x[1]-S12*csa: y[2]:=y[1]-S12*sna: x[9]:=x[2]-ra*sna: y[9]:=y[2]+ra*csa: x[4]:=x[3]-S34: x[5]:=x[4]: x[6]:=x[5]+1.4*Rc: y[6]:=y[5]: P[i]:=display(Cir(5,1),Cir(9,1),pruzh(x[6],x[11],y[6],3,24), box(6,8,rc), Cir(9,Ra),Cir(9,ra), Cir(5,Rc),Cir(5,rc), Cir(10,Rb), cir4(10,Rb,-phi), cir4(5,Rc,phi*Rb/(Rc-rc)),cir4(5,rc,phi*Rb/(Rc-rc)), cir4(9,Ra,phi*Rb/(Ra-ra)),cir4(9,ra,phi*Rb/(Ra-ra)), seq(TEXT([x[Nnam[j]]+2,y[Nnam[j]]+6],nam[j]),j=1..4), seq(Line(2*i-1,2*i,9),i=1..3),Lin(7,8)): 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,Опора(10,2),Поверхность(-120,Rb+Rc-rc,180,5), Поверхность(x[7],y[7],h,5), Стенка(x[11],Rb+Rc-rc,50,-5)); |
| > |
