1 ISMANS Calcul de structures en bureau d’études Calcul de structures en bureau d’études École thématique IN2P3 École thématique IN2P3 Porquerolles, 15 au 19 mai Porquerolles, 15 au 19 mai Calcul de structures en Calcul de structures en bureau d’études bureau d’études Jean Jean - - Charles CRAVEUR, ISMANS Charles CRAVEUR, ISMANS Arnaud CADIOU, SUBATECH Arnaud CADIOU, SUBATECH
53
Embed
École thématique IN2P3 Porquerolles, 15 au 19 mai Calcul ... · ISMANS 4 Calcul de structures en bureau d’études « Cahier des Charges » - calculer la réponse de la bielle
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
1ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
École thématique IN2P3École thématique IN2P3
Porquerolles, 15 au 19 maiPorquerolles, 15 au 19 mai
Calcul de structures en Calcul de structures en bureau d’étudesbureau d’études
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
La structure est extérieurement hypostatique : elle manque de liaison. La somme des forces extérieures
est nulle mais cela ne suffit pas à garantir la position d’équilibre : il y en a une infinité.
Quand le problème est écrit sous forme matricielle F = Kq et que l’on veut l’inverser pour calculer les
déplacements, l’hypostaticité se traduit par la non-inversibilité de la matrice de raideur, représentant
ainsi l’infinité de solutions possibles. Les termes diagonaux nuls du système à inverser sont les trop
célèbres « pivots nuls »…
11ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Pour une structure 3D, il faut bloquer les 6 mouvements de corps solide (modes rigides). Pour une structure 2D, il faut bloquer les 3
modes rigides.
Comment faire ici, sachant qu’il ne faut pas introduire de liaison qui modifierait la réponse
attendue du système ? Et comment contrôler le bien-fondé des liaisons introduites ?
12ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
13ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Il faut vérifier après le calcul qu’il n’y a pas de réaction sur un degré de liberté arbitrairement fixé pour éliminer un mode rigide. Si la réaction
est nulle, le « nœud » (en fait son ddl) n’avait pas tendance à se déplacer et l’avoir bloqué n’a rien
changé. Si la réaction est non nulle, le nœud avait tendance à se déplacer sous charge et sa
fixation a changé le problème : les conditions supplémentaires ne sont pas adaptées.
14ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Bielle2.sfieldAmplitude des déplacements
15ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
16ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
17ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Les déplacements sont petits, compatibles avec les hypothèse de linéarité. Mais les contraintes équivalentes dépassent la limite élastique : on ne peut pas prédire ce qui se passe dans la bielle.
D’où vient le problème et comment vérifier la tenue de la structure sous charge ?
18ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
On répartit la charge sur une surface au lieu de la concentrer sur une ligne : simule l’action de
l’axe qui n’est pas modélisé.
19ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
On bloque le déplacement radial à l’intérieur de l’alésage sur une ouverture angulaire
représentative du contact avec l’axe : simule l’action de l’axe qui n’est pas modélisé.
20ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Bielle3.sfieldAmplitude des déplacements
21ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
22ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
23ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
La contrainte est quasiment une fois et demie plus petite (1,9 GPa au lieu de 3 GPa) alors que le
modèle est beaucoup plus « physique » que précédemment.
Qu’est-ce qui peut provoquer ce problème puisque la vraie bielle tient sans problème ?
24ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Pour réduire les temps de calcul : exploitation de la symétrie, quart de bielle.
Les perçages sont de petite dimension : on suppose qu’il est possible de les éliminer du modèle.
quart1.sfield
25ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Conditions de symétrie à ne pas oublier, et 4500 N de résultante sur ¼ de la structure.
26ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Amplitude des déplacements
27ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
28ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Zoom sur les contraintes équivalentes extrapolées
29ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Les résultats ont beaucoup évolué entre ces deux modèles. Est-ce normal, est-ce
physique ?
QUESTION OUVERTE
30ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
La limite élastique du matériau est dépassée mais bien que les frontières des éléments ne soient pas affichées, on devine le maillage.
Et si le problème venait du choix des éléments, de leur
degré, de la façon de construire le maillage ?
31ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Maillage structuré (« réglé »)Eléments hexaédriques de degré 1
bielle-extrud-deg1.sfield
32ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Amplitude des déplacementsbielle-extrud-deg1.sfield
33ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
34ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
35ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
« Simplification »
- La bielle est relativement élancée,- La matière est symétriquement répartie par rapport à un plan moyen,- Les conditions aux limites sont symétriques par rapport à ce plan moyen,- Le chargement a une résultante qui est contenue dans le plan moyen.
36ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
On peut commencer par un modèle plan :
- Pas de dépouilles,
- Pas de congés hors plan.
MEMBRANES
quart3.sfield
37ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Amplitude des déplacements
38ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
39ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées (zoom)
40ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
La contrainte équivalente dépasse la limite élastique, même pour un modèle 2D.
N’a-t-on rien oublié ?
Et les coussinets ?
41ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
On modélise les deux coussinets, réalisés en alliage de cuivre CUSn9P.
Module d’élasticité : 115 000 MPa
Coefficient de Poisson : 0.33
Limite en rupture : 820 MPa
Limite élastique : 800 MPa
42ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’étudesbielle-coussinets1.sfield
43ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’étudesbielle-coussinets1.sfield
44ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Au montage, le diamètre extérieur des coussinets est de 0.05 mm supérieur au diamètre des alésages. Après le frettage, il existe des contraintes internes
dues à l’assemblage.
Mise en interférence initiale et calcul de
l’équilibre pour retrouver des positions compatibles
et l’état de contrainte initial.
45ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
46ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
47ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
La contrainte équivalente est de l’ordre de 350 MPa après montage, dans les coussinets. La contrainte est plus faible dans la bielle.
Se superpose à ce champ de contrainte celui qui est induit par la
traction exercée sur la bielle.
Peut-on prévoir ce qui se passe en contraintes en utilisant tous les
résultats précédents ?
48ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
NON
La traction crée des contraintes plutôt « positives ».
Car en ayant les contraintes équivalentes, on n’a pas le signe de la contrainte.
La somme des deux n’est pas la somme des contraintes équivalentes…
Le frettage crée des contraintes plutôt « négatives ».
49ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’étudesbielle-coussinets2.sfield
Charges et conditions aux limites (hors symétrie)
50ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Amplitude des déplacements
51ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
52ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études
Contraintes équivalentes extrapolées
53ISMANS
Calcul de structures en bureau d’étudesCalcul de structures en bureau d’études