Voila je travail actuellement sur un solver élément finis, Le problème est actulement linéaire, l'équation d'équilibre est de la forme A.x=b, avec A matrice de coeffs réel, x le vecteur inconnues et b un vecteur de coeffs réels.
La résolution se fait en inversant la matrice A avec un pivot de gauss. (x=A^-1.F)
J'aimerais a présent intégrer un comportement non linéaire, qui se caractériserais par un changement de l'équation d'équilibre: A1.x^2+12.x=F (dans un premier temps, et si possible passer au degré 3).
Existe il une méthode de résolution directe?, ou un algorithme pas trop gourmant en temps d'execution pour ce genre de problème?
Merci d'avance pour votre aide.
Partager