Animation des Rades

Hier wird eine Animation des Präzedierenden Rades dargestellt. Die einzelnen Vektoren und Frequenzen entsprechen betragsmässig nicht den im vorigen Abschnitt berechneten. Der Schrarze Vektor beschreibt die Gravitationskraft, der blaue Vektor die Rotations-und Drehimpulsachse und der grüne Vektor das (durch die Gravitationskraft) ausgeübte Drehmoment.

>

#Veranschaulichung des Rades, der Vektoren und weiterer 3D-Objekte Torus1:=torus([0,0,0],8,40): Torus2:=rotate(Torus1,0,Pi/2,0): ### WARNING: calls to `translate` from the `polytools` package should be replaced by polytools[translate] rad:=translate(Torus2,30,0,0): VOS:=vector([30,0,0]): VGr:=vector([00,0,-59]): POS:=arrow([0,0,0],VOS,2,4,0.4,color=red): PGr:=arrow(VOS,VGr,2,4,0.1,color=black): PRot:=arrow([0,0,0],scalarmul(VOS,2),2,4,0.2,color=blue): PD:=arrow(VOS,[30,-20,0],2,4,0.4,color=green): wiese:=plot3d(-60,x=-35..35,y=-35..35,color=green): haus:=cuboid([-5,-5,-60],[5,5,-0.05]):

>

frames:=30: for i from 0 by 1 to frames do Ani1[i] := rotate(rad,i*5*2*Pi/frames,0,i*2*Pi/frames): Ani2[i] := rotate(POS,i*5*2*Pi/frames,0,i*2*Pi/frames): Ani3[i] := rotate(PGr,0,0,i*2*Pi/frames): Ani4[i] := rotate(PRot,i*5*2*Pi/frames,0,i*2*Pi/frames): Ani5[i] := rotate(PD,0,0,i*2*Pi/frames): Ani[i] := display({Ani1[i],Ani2[i],Ani3[i],Ani4[i],Ani5[i],wiese,haus},axes=none,orientation=[25,75]); od:

>	display([seq(Ani[i],i=0..frames)],insequence=true,scaling=constrained,title="Animation des Rades");

>