Bonjour
J'ai remarqué une divergence des toupies qui utilisent les formules d'évolution donné sur la page 6 (alors toupie generale) en combinaison avec Runge Kutta et l'integrateur Euler-cromer. Au début j'ai pensé que c'était une faute d'implementation des integrateurs comme ce problème n'aparaisse pas avec Newmark, mais je pense le problème est au niveau des formules données sur la page 6. On a la condition que si theta = 0, alors on calcule les accelerations d'une manière différent pour éviter la division par sin(theta). J'ai implementé la condition que abs(theta + precision) >0 avec precision = 10e-6 mais il me semble que ceci ne suffit pas comme condition, surtout pour l'intégrateur Runge Kutta qui fait plusieurs fois appel à la méthode f et donc les erreurs vont s'agrandir à chaque appel (division par sin(theta)).
Je pense que Newmark ne pose pas des problèmes comme il a un mécanisme de contrôle.
Est-ce que vous une idée pour éviter ce défaut? J'ai voulu juste augmenter la borne en posant comme condition theta > 0.05 mais ceci me semble assez arbitraire et on n'obtient plus les bonnes valeurs.
Merci pour votre réponse