Колебание цилиндра
В вышеприведенной команде опция points дает
количество кадров в единичном интервале времени,
а опция frames - полное число кадров.
Вы можете, изменяя углы, начальную угловую
скорость, точку подвеса и даже сам объект
получать различные движения, например
> Object:=cylinder;
Axes:=[3,0,10];
Center:=[1,1,-3];
![[Maple Math]](Images/Phys_pend_cyl1.gif)
![[Maple Math]](Images/Phys_pend_cyl2.gif)
![[Maple Math]](Images/Phys_pend_cyl3.gif)
> Mass:=rho*physics[body_inert](m,Axes,object=Object,'volume');
![[Maple Math]](Images/Phys_pend_cyl4.gif)
> motion(Axes,[op(Center),op(Angles)],Velocity,type=Type,
object=Object,static,
style=contour,contours=20,color=wheat,axes=normal,scaling=constrained,
abels=[r1,r2,r3],labelfont=[TIMES,BOLD,14]);
![[Maple Plot]](Images/Phys_pend_cyl5.gif)
> n0:=physics[conv_coord]((r,v)->r,[1,0,0],move_c(v),frame=inv_zveer,angles=Angles);
l0:=physics[conv_coord]((r,v)->r,Center,move_c(v),frame=zveer,angles=Angles);
d0:=evalf(sqrt(evalm([0,l[2],l[3]]&*[0,l[2],l[3]])));
![[Maple Math]](Images/Phys_pend_cyl6.gif)
![[Maple Math]](Images/Phys_pend_cyl7.gif)
![[Maple Math]](Images/Phys_pend_cyl8.gif)
> motion([a,0,h],[x,y,z,op(Angles)],[omega[theta]],density=rho,intencity=[0,0,g],
type=Type,object=Object,equations);
![[Maple Math]](Images/Phys_pend_cyl9.gif)

![[Maple Math]](Images/Phys_pend_cyl10.gif)
![[Maple Math]](Images/Phys_pend_cyl11.gif)
![[Maple Math]](Images/Phys_pend_cyl12.gif)
![[Maple Math]](Images/Phys_pend_cyl13.gif)
![[Maple Math]](Images/Phys_pend_cyl14.gif)
![[Maple Math]](Images/Phys_pend_cyl15.gif)
![[Maple Math]](Images/Phys_pend_cyl16.gif)
![[Maple Math]](Images/Phys_pend_cyl17.gif)
![[Maple Math]](Images/Phys_pend_cyl18.gif)
![[Maple Math]](Images/Phys_pend_cyl19.gif)
![[Maple Math]](Images/Phys_pend_cyl20.gif)
![[Maple Math]](Images/Phys_pend_cyl21.gif)

> motion(Axes,[op(Center),op(Angles)],Velocity,type=Type,object=Object,equations);
![[Maple Math]](Images/Phys_pend_cyl22.gif)
![[Maple Math]](Images/Phys_pend_cyl23.gif)

![[Maple Math]](Images/Phys_pend_cyl24.gif)

![[Maple Math]](Images/Phys_pend_cyl25.gif)
![[Maple Math]](Images/Phys_pend_cyl26.gif)

![[Maple Math]](Images/Phys_pend_cyl27.gif)
Численное решение задачи
> Q:=motion(Axes,[op(Center),op(Angles)],Velocity,type=Type,object=Object,solution);
![[Maple Math]](Images/Phys_pend_cyl28.gif)
Строим график зависимости угла от времени
> plot(Q[1],0..3);
![[Maple Plot]](Images/Phys_pend_cyl29.gif)
Генерируем анимацию движения физического
маятника
> motion(Axes,[op(Center),op(Angles)],Velocity,type=Type,object=Object,style=patch,
lightmodel=`light4`,shading=`zhue`,frames=35,points=70,scaling=constrained);
 
|