Discussion :

[grand's]

Salut,

J'ai une fonction qui mets a jour la vitesse d'un objet.
Pour sa vitesse lineaire, pas de souci
j'ai la somme des forces exterieures et l'intervalle de temps, donc :
acceleration = forces / masse;
newvitesse = oldvitesse + acceleration * dt;
newposition = oldposition + vitesse * dt;

Mais cette objet peut egalement tourner.
Et je sais pas comment faire pour calculer la nouvelle vitesse angulaire.

Si quelqu'un peut eclairer ma lanterne,

[Jedai]

Si ton objet tourne selon un axe fixe (par rapport à l'objet), tu dois calculer le moment d'inertie J de ton objet selon cet axe, puis utiliser :
J * (dérivée de la vitesse de rotation) = moment des forces selon l'axe de rotation .
Si ton axe n'est pas fixe, ça se complique un peu (encore que certains autres cas soient relativement simple).
Je te conseille de consulter un petit cours de mécanique du solide, même superficiellement, ça t'éclairera un peu les idées !

--
Jedaï

[grand's]

bon j'ai repeché quelque chose dans un vieux cours de meca.

vAng = vAng + (~I * (M - (vAng ^ ( I * vAng))) * dt

avec :
I = matrice d'inertie
~I = matrice d'inertie inverse
M = moment

ca pourrait coller ?

[Matthieu Brucher]

Peut-être...

[grand's]

Peut-être...

j'aime bien

l'autre probleme c'est que je comprend pas cette equation ... et ca m'enerve
Si quelqu'un peut m'expliquer rapidement

[Matthieu Brucher]

Le nouvel angle est égal à l'ancien + une fraction dépendante de l'angle.
Logiquement, ça dépend du moment et de la matrice inverse de l'inertie, mais s'il y a de l'inertie, ça compense aussi le moment, donc on va retirer au moment une partie qui est fonction de l'angle et de l'inertie - enfin, ça me paraît louche, c'est quoi l'opérateur ^ ? -

[grand's]

le produit en croix de deux vecteur (cf code)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
// Vector cross product (u cross v)
inline	Vector operator^(Vector u, Vector v)
{
	return Vector(	u.y*v.z - u.z*v.y,
					-u.x*v.z + u.z*v.x,
					u.x*v.y - u.y*v.x );
}

[Matthieu Brucher]

C'est ce que je pensais. et * ?
Ca m'a l'air correct : s'il n'y a pas de vitesse angulaire au départ, elle sera fonction du moment et de l'inertie de l'objet.

[grand's]

le * c'est la multiplication d'une matrice par un vecteur (matrice 1x3)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
 
inline	Vector operator*(Matrix3x3 m, Vector u)
{
	return Vector(	m.e11*u.x + m.e12*u.y + m.e13*u.z,
					m.e21*u.x + m.e22*u.y + m.e23*u.z,
					m.e31*u.x + m.e32*u.y + m.e33*u.z);					
}

[Matthieu Brucher]

C'est logique en plus...
A priori, ça m'a l'air correct comme équation

[grand's]

a moi aussi maintenant,

surtout apres quelques recherches sur le net.
Et surtout apres avoir decouvert que l'auteur s'etait planté en mettant * au lieu de ^ (alors qu'il le mettait partout jusque la)

Merci bcp.