4. Synthèse des résultats

On observe que les résultats sont très proches de la solution analytique, aussi bien pour le schéma d’Euler que pour le schéma des différences centrées. Les instants de changement de raideur obtenus par ces deux schémas sont identiques.

Les bilans énergétiques présentent une erreur inférieure à \(0,1J\) dans les deux cas. Cette valeur est à relativiser compte-tenu du pas de temps très faible qui a tendance à réduire les écarts, mais un pas de temps aussi faible est nécessaire pour traiter ce problème raide.

Enfin, il y a une adéquation parfaite (à la précision numérique près) entre les forces de contact et la cinématique, ce qui signifie que le traitement du contact par pénalisation est correctement effectué.