Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 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)); |
| > |
| > |
