Pensez à lire la Charte avant de poster !

$\newcommand{\K}{\mathbf K}$


Les-Mathematiques.net - Cours de mathématiques supérieures
 Les-Mathematiques.net - Cours de mathématiques universitaires - Forum - Cours à télécharger

A lire
Deug/Prépa
Licence
Agrégation
A télécharger
Télécharger
127 personne(s) sur le site en ce moment
E. Cartan
A lire
Articles
Math/Infos
Récréation
A télécharger
Télécharger
Théorème de Cantor-Bernstein
Théo. Sylow
Théo. Ascoli
Théo. Baire
Loi forte grd nbre
Nains magiques
 
 
 
 
 

MEF : Poutre soumise à charge triangulaire

Envoyé par huffman 
MEF : Poutre soumise à charge triangulaire
l’an passé
Bonjour à tous,

Je me trouve actuellement confronté à un problème d'élements finis que je n'arrive pas à résoudre...

Le problème est une poutre de Berbouilli avec un charge triangulaire


Traduit en équation de cette manière :

edit : petite erreur sur la 4eme conditions, c'est v'''(L) et non v'''(x)



Ma solution analytique (vérifiée avec wolframalpha) est la suivante :


Après quelques petits calculs, j'obtiens la formulation faible suivante :



Suite à cela je calcul ma matrice de rigidité (qui semble bonne vis-à-vis de ce que j'ai pu lire sur internet)


Le calcul des forces nodales me parait plus bancale, après être passé dans un repère local, j'obtient la formule suivante :


Et après intégration :


Ensuite j'ajoute mes conditions initiales.
Pour bloquer le v(0)=0, je supprime la première ligne et colonne de ma matrice de rigidité
et pour rajouter le -kv(L)Ni(L) du vecteur de force, je rajoute un -k sur l'avant dernière ligne-avant dernière colonne de ma matrice de rigidité.

Après calcul de matlab, j'obtient qqch du genre :


La forme est pas mal, mais visiblement une petite erreur s'est glissée dans mes calculs....
Je pense que mon erreur vient de mon vecteur de force, et/ou de mes conditions aux limites que je ne dois pas correctement mettre en place...


Un expert en élements finis voit-il une erreur dans mon résonnement ?

Merci d'avance à tous



Edité 2 fois. La dernière correction date de l’an passé et a été effectuée par huffman.
Re: MEF : Poutre soumise à charge triangulaire
l’an passé
Après analyse détaillée de mon raisonnement, je viens de me rendre compte que l'erreur vient simplement d'un mauvais signe au moment d'appliquer la condition aux limites n°4... Merci à tous.



Edité 1 fois. La dernière correction date de l’an passé et a été effectuée par AD.
Seuls les utilisateurs enregistrés peuvent poster des messages dans ce forum.

Cliquer ici pour vous connecter

Liste des forums - Statistiques du forum

Total
Discussions: 134 237, Messages: 1 292 707, Utilisateurs: 23 272.
Notre dernier utilisateur inscrit Toutous.


Ce forum
Discussions: 2 219, Messages: 16 430.

 

 
©Emmanuel Vieillard Baron 01-01-2001
Adresse Mail:

Inscription
Désinscription

Actuellement 16057 abonnés
Qu'est-ce que c'est ?
Taper le mot à rechercher

Mode d'emploi
En vrac

Faites connaître Les-Mathematiques.net à un ami
Curiosités
Participer
Latex et autres....
Collaborateurs
Forum

Nous contacter

Le vote Linux

WWW IMS
Cut the knot
Mac Tutor History...
Number, constant,...
Plouffe's inverter
The Prime page