Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге
8.3. Задача 2
> | restart: |
> | read «ris.m»; |
Подписи шарниров на рисунке
> | nam:=array(1..4,[‘O’,’A’,’B’,’C’]): Начальные значения и размеры ( cм) Замечание: OC уменьшено с 87см для того, чтобы механизм делал полный оборот |
> | OA:=29: BC:=47: AB:=41: |
Угловая скорость OA и период обращения
> | omegaOA:=4: T:=2*pi/omegaOA: |
Координаты опоры O
> | x[1]:=0:y[1]:=0: |
Координаты опоры C
> | x[4]:=59: y[4]:=0: |
Количество кадров K
> | K:=12:with(plots):with(plottools): |
Создаем все кадры
> | for i from 0 to K do |
> | t:=i/K*T: |
Закон движения OA
> | phi:=t*omegaOA+pi/3: Координаты шарнира A |
> | x[2]:=OA*cos(phi):y[2]:=OA*sin(phi): koord(4,2,3,BC,AB): P[i]:=display(Cir(3,1),Cir(2,1), seq(Line(i,i+1,6+i),i=1..3), seq(TEXT([x[j]-2,y[j]+3],nam[j]),j=1..4)): 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,1),Опора(4,1)); |
> |