Bon, j'ai enfin pu trouver le temps de reprendre personnellement en détail le cas de la toupie chinoise et vous avez raison, les nombres indiqués dans ce message https://moodlearchive.epfl.ch/2019-2020/mod/forum/discuss.php?d=32800#p69826 étaient erronés. Voici ce que j'obtiens :
Système au temps t = 0.01
Toupie 1: Ω=(0.5, 0.109555, 0) Ω'=(50, -0.04447, 0)
A : (0, 0, 0)
G : (-9.25034e-05, 0.000169326, 0.148246)
Energie totale: 0.0178868
Projection de LA sur k: 0.00063428
Projection de LA sur a: 0.000630888
Produit mixte [vitesse de rot,L,a]: 4.43369e-20
Système au temps t = 0.02
Toupie 1: Ω=(1.00199, 0.108668, 6.28118) Ω'=(50.1989, -0.088762, -0.200133)
A : (-3.19801e-05, 5.85392e-05, 0)
G : (-0.000193234, 0.000161627, 0.148246)
Energie totale: 0.0178869
Projection de LA sur k: 0.00063428
Projection de LA sur a: 0.000630842
Produit mixte [vitesse de rot,L,a]: -3.03608e-20
Système au temps t = 0.03
Toupie 1: Ω=(1.50801, 0.107448, 6.27512) Ω'=(50.6023, -0.121974, -0.606102)
A : (-0.000126622, 0.000157685, 0)
G : (-0.000315499, 0.000169559, 0.148245)
Energie totale: 0.0178869
Projection de LA sur k: 0.000634279
Projection de LA sur a: 0.000630716
Produit mixte [vitesse de rot,L,a]: 3.944e-21
Système au temps t = 0.04
Toupie 1: Ω=(2.01974, 0.106117, 6.26333) Ω'=(51.1723, -0.133044, -1.17946)
A : (-0.000303106, 0.000266471, 0)
G : (-0.000471499, 0.000185348, 0.148245)
Energie totale: 0.017887
Projection de LA sur k: 0.000634276
Projection de LA sur a: 0.000630524
Produit mixte [vitesse de rot,L,a]: -4.49986e-20
qui correspond plus, je crois, à vos données. Non ?