Bonjour,
J'ai fais un programme matlab pour déterminer les fréquences propres d'une poutre 2D (quadrangle). J'utilise la méthode des sous-espaces.
Cela marche très bien pour une poutre encastrée-libre, je trouve les mêmes résultats que Abaqus (fréquences et déplacements des noeuds).
Le pb se corse lorsque je supprime les conditions limites (vibrations libres).
Les 1ere fréquences que je trouve sont soient réel soit imaginaires (mais très faibles) et les fréquences qui correspondent au mode non-rigide sont infinies....
D'après ce que j'ai lu (google books), si la matrice de masse est bien conditionnée, cela ne devrait pas poser de pb !!
J'ai fait classiquement :
Melementaire=N*Rho*N'*det(J)*w t(gauss)*Ep(def. plane)
-> Ca marche pas
et aussi, après assemblage (Mglobal)
Mglobal=diag(sum(M))
-> Idem
Quelqu'un a une idée pour m'éclairer ?
je peux faire autre chose ?
Il est a noter que c'est deux techniques donnent les bons résultats pour la poutre encastrée-libre...
Merci d'avance !!!
A+
-----