13-5S3A.html | AcademiaXXI

Кирсанов М.Н.

Решебник.Теоретическая механика.

 М.:ФИЗМАТЛИТ, 2002.

Иллюстрации и программы к книге

13.5. Пример 3

restart:

with(plottools): with(plots):

Количество кадров и размеры
K:=12:a:=50:R:=10:h:=30: L:=35:
Диск1:=cylinder([0,0,0],R,h):
Диск2:=cylinder([0,0,-5],R,3):
Ось:=cylinder([0,0,-6],1,7):
Стержень:=cuboid([-2,0,h/2-1],[2,a+2,h/2+1]):
Поршень:=cuboid([R,0,0],[-h,h,h]):
Опора:=cuboid([R,50,-6],[1+R,150,h+6]):
for i to K do
A:=sin(2*Pi*i/K)*Pi/14+6:
P1:=display(rotate(PLOT3D(Диск1,ORIENTATION(0,-30)),0,0,A)):
P2:=display(rotate(PLOT3D(Диск2),0,0,-A*a/R*cos(A)+A),Ось):

      P3:=display(rotate(PLOT3D(Стержень),0,0,A)):
   F:=arrow([0,-L,0],color=red):
      P[i]:=display(translate(P1,0,R*A,0),
   translate(P2,a*sin(A),a*cos(A)+R*A,h/2),
   translate(P3,0,R*A,0),
   translate(F,-10,R+h+L+a*cos(A)+R*A,h/2),
   translate(Поршень,0,R+a*cos(A)+R*A,0));

od:

Изображение механизма в движении

     PP:=display(seq(P[i],i=1..K),insequence=true,
                               thickness=2,
                               scaling=constrained,
                               axes=none):

display(PP,Опора);