Кирсанов М.Н. Решебник.Теоретическая механика. М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге
9.4. Пример 2
restart:
read «ris.m»:
Подписи шарниров на рисунке
nam:=array(1..6,[O,A,B,C,D,E]):
Начальные значения и размеры ( cм)
OA:=30: AB:=30: BC:=35: DE:=42:
Время, амплитуда качания
T:=pi/18: Амплитуда:=pi/14:
Координаты опор O, C, E
x[1]:=0:y[1]:=0: x[4]:=52:y[4]:=0: x[6]:=52:y[6]:=30:
Количество кадров K
K:=12:with(plots):with(plottools):
Создаем все кадры
for i from 0 to K do
t:=i/K*T:
phi:=sin(2*pi*t/T)*Амплитуда+pi/3:
Координаты шарнира A
x[2]:=OA*cos(phi):
y[2]:=OA*sin(phi):
koord(2,4,3,AB,BC):
y[5]:=’y[5]’: #’Очистка’ переменной
tga:=(x[2]-x[3])/(y[2]-y[3]);
x[5]:=x[3]+(y[5]-y[3])*tga;
y[5]:=solve((y[5]-y[6])^2+(x[5]-x[6])^2-DE^2)[2];
P[i]:=display(seq(Cir(i,0.5),i=2..3),Cir(5,1.2),
seq(TEXT([x[j]-2,y[j]+3],nam[j]),j=1..6),Line(5,6,6),seq(Line(i,i+1,6+i),i=1..3)):
od:
Изображение механизма в движении
PP:=display(seq(P[i],i=0..K),insequence=true,
thickness=2,
scaling=constrained,
axes=none):
display(PP,seq(Опора(i,1),i=[1,4,6]));
Warning, the name changecoords has been redefined
Warning, the name arrow has been redefined
> |