Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге
13.8. Пример 2
> | restart: |
> | read «ris.m»; |
Подписи шарниров на рисунке
> | nam:=[O,A,B]: Pазмеры ( cм) r1:=10: R1:=20: |
> | R2:=20: OA:=60: OB:=OA-r1-R2: |
Координаты опоры O
> | x[1]:=0:y[1]:=0: |
Количество кадров K
> | K:=12:with(plots):with(plottools): |
Создаем все кадры
> | for i from 0 to K do |
Закон движения АO
> | Амплитуда:=pi/10: |
> | phi:=3*pi/4+sin(2*pi*i/K)*Амплитуда: Координаты шарнира А |
> | x[2]:= OA*cos(phi): y[2]:= OA*sin(phi): Координаты шарнира B x[3]:= OB*cos(phi): y[3]:= OB*sin(phi): CLR:=red: cc1:=Cir(2,1),Cir(3,1): CLR:=green: cc2:=Cir(2,R1),Cir(2,r1): CLR:=magenta:cc3:=Cir(3,R2): CLR:=black: phi1:=-3*phi: phi2:=3*phi: P[i]:=display(cc1,Cir(1,OA+R1),cc2,cc3, cir4(2,R1,phi1),cir4(2,r1,phi1), cir4(3,R2,phi2), Lin(1,2),seq(TEXT([x[j]+3,y[j]+4],nam[j]),j=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,Опора(1,1)); |
> |