Résolution de 2 equa diff d'ordre 2 à 4 inconnues
Bonjour à tous dans le cadre d'un projet je m'intéresse au mouvement du pendule double et j'ai réussis à établir les équations du mouvements qui
ne sont pas linéaires et assez compliquées.
Je compte donc utiliser python pour les résoudre avec n'importe quel module (odeint par exemple) mais je ne vois pas du tout comment faire :roll:
Pour vous donner une idée voici les équations :
Ӫ1=(-m2+m1)l1sin(θ1-θ2)cos(θ1-θ2)(θ͘1)2-m2l2sin(θ1-θ2)(θ͘2)2-m1gsinθ1-m2gsin(θ1- θ2)cosθ2)/(m1l1+m2l1sin2(θ1-θ2))
Ӫ2=((m1+m2)l1sin(θ1-θ2)(θ͘1)2+m2l2sin(θ1-θ2)cos(θ1-θ2)(θ͘2)2+(m1+m2)gsin(θ1- θ2)cosθ1)/(m1l2+m2l2sin2(θ1-θ2))
En gros j'aimerai afficher le mouvement pour un laps de temps donner ou faire un portrait de phase avec l1,l2,m1,m2 et g fixés ainsi que des conditions
initiales sur θ1 et θ2 ainsi que sur θ͘1 et θ͘2
Est ce que vous auriez des conseils svp ?
