Кирсанов М.Н. 
Решебник.Теоретическая механика.
 М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге

12.7. Задача 4

Динамические реакции вращающегося вала

>    restart:

>    with(plots):with(plottools):

    Количество кадров

>     K:=12:
 
Перекос маховика (рад)
 alpha:=0.12:

>     Маховик:=cylinder([0,0,0],4,1,style=PATCHNOGRID):
 
Подвижные оси координат
cl:=COLOR(RGB,0,0,0):
c:=FONT(SYMBOL,16),cl:bl:=color=black:
        tx:=TEXT([5,0,-4],"x",cl),
            TEXT([0,5,-4],"y",cl),
            TEXT([0,0,14],"z",cl),
            line([0,0,-4],[0,5,-4],bl),
            line([0,0,-4],[5,0,-4],bl),
            line([0,0,-4],[0,0,13],bl):
 L1:=line([0,0,0],[0,5,0],bl):
 L2:=line([5,0,0],[0,0,0],bl):
 L3:=line([0,0,0],[0,0,5],bl):
 T1:=TEXT([0,0,5],"z",c):
 T2:=TEXT([0,5,0],"x",c):
 T3:=TEXT([5,0,0],"h",c):
 Стержень:= cuboid([0.2,0,8.8], [-0.2,6,9.2],color=red):
 ОпораА :=cuboid([-3,-1,-4.5],[1,1,-4]):
 ОпораВ :=cuboid([-3,-1,6.5],[1,1,6]):
 Груз:=cylinder([0,6,8.8],0.5,0.5,color=blue):
 Ось:=cylinder([0,0,-5],0.2,14):
for i to K do
 A:=2*Pi*i/K:
 P[i]:=display(rotate(PLOT3D(Маховик,L1,L2,L3,T1,T2,T3),alpha, 0,A),ОпораА,ОпораВ,
  rotate(PLOT3D(Стержень,Груз,tx),0, 0,A),
  rotate(PLOT3D(Ось,ORIENTATION(80,45)),0, 0,A)):
 od:

Warning, the name changecoords has been redefined

Warning, the name arrow has been redefined

>    display(seq(rotate(P[i],0,Pi/2,0),i=1..K),insequence=true,scaling=constrained,title="Динамические реакции вращающегося вала",lightmodel=light2);

[Maple Plot]

>