Initiation au calcul des structures dans le domaine

  • Slides: 33
Download presentation
Initiation au calcul des structures dans le domaine plastique Elasto - plasticité en petite

Initiation au calcul des structures dans le domaine plastique Elasto - plasticité en petite transformation Cours Passer à la première page

Objectifs • Introduction à la plasticité « classique » • Aborder la résolution de

Objectifs • Introduction à la plasticité « classique » • Aborder la résolution de problème non linéaire • Utiliser un code de calcul en non linéaire Passer à la première page

Plan • Aspect physique - modèles analogiques • Elasto-plasticité des barres (1 D) étude

Plan • Aspect physique - modèles analogiques • Elasto-plasticité des barres (1 D) étude des treillis • Elasto-plasticité des poutres (3 D ==> 1 D) étude des portiques (rotule plastique) • Critères 3 D & règles d ’écoulement plastique Passer à la première page

Supports : • Polycopié de cours • Sur le site intranet de l ’école

Supports : • Polycopié de cours • Sur le site intranet de l ’école • Cours, • Exercices, projets • Doc d’utilisation et tutorials des logiciels • ALGOR • CASTEM Evaluation : • 2 notes : • TD - TA • Projet sur 10 Passer à la première page

Aspects physiques Elastiques (instantanées - réversibles) Déformations Visqueuses (fct du temps) Plastiques (irréversibles -

Aspects physiques Elastiques (instantanées - réversibles) Déformations Visqueuses (fct du temps) Plastiques (irréversibles - non linéaire) Passer à la première page

Ecrouissage imposé mesuré les Essais Fluage - recouvrance imposé mesuré Relaxation imposé mesuré Passer

Ecrouissage imposé mesuré les Essais Fluage - recouvrance imposé mesuré Relaxation imposé mesuré Passer à la première page

Modèles rhéologiques Modèles « linéaires » ==> solides visco - élastiques Le ressort ou

Modèles rhéologiques Modèles « linéaires » ==> solides visco - élastiques Le ressort ou de l ’amortisseur peuvent être non linéaire Passer à la première page

 Modèles non linéaires Modèles de base de la plasticité Passer à la première

Modèles non linéaires Modèles de base de la plasticité Passer à la première page

Essai de traction Impossible à mesurer ==> Norme observation Ecrouissage monotone s A so

Essai de traction Impossible à mesurer ==> Norme observation Ecrouissage monotone s A so Domaine élastique Cycles charges- décharges ep(A) ee (A) e Passer à la première page

Ce n ’est pas si simple !! 2ème passage 1 er passage Il faut

Ce n ’est pas si simple !! 2ème passage 1 er passage Il faut connaître l ’historique du chargement Caractère incrémental des lois de comportement en plasticité On utilise un temps cinématique = temps réel / Passer à la première page

Ecrouissage ISOTROPE Ecrouissage monotone s Même augmentation en traction et compression 2 so Wdef

Ecrouissage ISOTROPE Ecrouissage monotone s Même augmentation en traction et compression 2 so Wdef élastique e Ecrouissage CINEMATIQUE Effet Baushinger durcissement dans un sens adoucissement dans l ’autre Passer à la première page

Modélisation s ET Etat présumé Etat réel Ds ss(E) Etat actuel (s, E) E

Modélisation s ET Etat présumé Etat réel Ds ss(E) Etat actuel (s, E) E Critère Loi d ’écoulement plastique e Passer à la première page

Evolution élasto-plastique des Treillis Résolution Explicite • Analytique (RDM) • par la MEF •

Evolution élasto-plastique des Treillis Résolution Explicite • Analytique (RDM) • par la MEF • Modèles EPP et EPE • Chargements cycliques Résolution Numérique (MEF) • Algorithmes de résolution des Pb non linéaire • Newton-Raphson et N-R modifié • Algorithme de projection sur le critère • Calcul du résidus {Fint}-{Fext} Passer à la première page

Application Structure à étudier Modèle EPP (pas d ’écrouissage) Passer à la première page

Application Structure à étudier Modèle EPP (pas d ’écrouissage) Passer à la première page

N 2 A' so S 7 A 2 B 3 Décharge élastique 5 Phase

N 2 A' so S 7 A 2 B 3 Décharge élastique 5 Phase élasto-plastique en compression 1 Phase élastique 2 Phase élasto-plastique En E ruine de la structure // à 1 En B ruine de la structure En A : 6 Phase élastique En C (F = 0) Contraintes résiduelles : 7 Phase élasto-plastique Domaine d ’élasticité 1 C' 3 so S 6 4 En D (F = - so. S) 4 5 D ne vérifient pas les conditions de compatibilité Le domaine d'élasticité n'est pas modifié au cours du chargement car le matériau est EPP C E N 1 On ne peut pas calculer directement les déformations plastiques de la barre 2 Les déformations plastiques D Il y a plastification de la barre 2 en compression Passer à la première page

N 2 Ecrouissage ISOTROPE N 1 Passer à la première page

N 2 Ecrouissage ISOTROPE N 1 Passer à la première page

N 2 Ecrouissage CINEMATIQUE N 1 Passer à la première page

N 2 Ecrouissage CINEMATIQUE N 1 Passer à la première page

Résolution numérique Principe Incrément de charge donné Calcul élastique ==> Pour chaque élément Élastique

Résolution numérique Principe Incrément de charge donné Calcul élastique ==> Pour chaque élément Élastique ou non ? ALGORITHME DE PROJECTION ==> Assemblage et calcul du résidu {R} = {DFext}-{DFint} Si Il faut itérer Passer à la première page

Algorithme de résolution • Maillage éléments finis • Définition de l ’historique du chargement

Algorithme de résolution • Maillage éléments finis • Définition de l ’historique du chargement (incréments DF) F • Définition des lois de comportement so , ET • Calcul de [K] F = K(u) U DF Pour chaque incrément {DF} Initialisation du résidu {R} <== {DF} Tant que || {R} || > e Solution cherchée {R 1} Calcul de {DU} = [K]-1 {R} Calcul de [K(u)] Pour chaque élément Convergence la plus rapide déformation ==> projection {DFint}e Convergence lente fin pour Nouveau résidu {R} = {R} - {DFint} fin tant que Impression des résultats de l ’incrément fin pour U Passer à la première page

PROJECTION sur le critère de plasticité Etat présumé si s+Dse > ss(E) élément en

PROJECTION sur le critère de plasticité Etat présumé si s+Dse > ss(E) élément en cours de plastification Projection sur le critère ET Dse ss(E) Etat actuel (s, E) s De Cas particulier s = ss(E) incrément plastique Etat réel si non s+Dse < ss(E) incrément élastique Dep Dee E e Passer à la première page

A vous de jouer. . . Passer à la première page

A vous de jouer. . . Passer à la première page

Plasticité des poutres Rappels Hypothèses de Bernoulli 1 2 L. C. Elastique Equations «

Plasticité des poutres Rappels Hypothèses de Bernoulli 1 2 L. C. Elastique Equations « PFD » 3 Intégrée sur la section 4 4 équations pour 4 champs ==> résolution En plasticité ce qui change c ’est la LC intégrée on pose Passer à la première page

Objectif : obtenir la LC élasto-plastique Essai de flexion pure ==> Mf = cte

Objectif : obtenir la LC élasto-plastique Essai de flexion pure ==> Mf = cte Début de plastification des fibres pour Pour LC élasto-plastique Passer à la première page

Facteur de forme plastique d ’une section Représente la réserve vis à vis de

Facteur de forme plastique d ’une section Représente la réserve vis à vis de la plastification totale : 1 pas de réserve Passer à la première page

Essai de flexion simple Négligeons T et utilisons la LC ==> Passer à la

Essai de flexion simple Négligeons T et utilisons la LC ==> Passer à la première page

Etude de Comportement asymptotique ==> notion de rotule plastique Modèle simplifié - rotule plastique

Etude de Comportement asymptotique ==> notion de rotule plastique Modèle simplifié - rotule plastique Application Etude des portiques charges limites Passer à la première page

Phase élastique ==> critère ==> Rotule plastique en A pour Phase élasto-plastique critère ==>

Phase élastique ==> critère ==> Rotule plastique en A pour Phase élasto-plastique critère ==> F 2 charge de ruine Passer à la première page

A vous de jouer. . . Passer à la première page

A vous de jouer. . . Passer à la première page

Plasticité 3 D Critère (fonction de charge) État actuel est intérieur au domaine d

Plasticité 3 D Critère (fonction de charge) État actuel est intérieur au domaine d ’élasticité État actuel se situe sur la frontière du domaine (état plastique) Un état extérieur au domaine d’élasticité est physiquement impossible à obtenir Le domaine d’élasticité représente l’ensemble des états s admissibles. Ecrouissage « E » Modèle EPE Doit respecter les symétries matérielles Modèle EPP Passer à la première page

Critère de Von Mises (1910) Pour les métaux ==> utilisation du déviateur des s

Critère de Von Mises (1910) Pour les métaux ==> utilisation du déviateur des s Fonction de charge Partie sphérique { ==> essai de traction : s. VM Contrainte de Von Mises Passer à la première page

Critère de Tresca (1870) Contraintes principales s. T Contrainte de Tresca Critère dit de

Critère de Tresca (1870) Contraintes principales s. T Contrainte de Tresca Critère dit de cisaillement maximal Tresca plus sévère de 13% Passer à la première page

Loi d ’écoulement plastique Principe de Hill ==> De est convexe dep est normal

Loi d ’écoulement plastique Principe de Hill ==> De est convexe dep est normal à la surface de charge Dissipation plastique associée à Loi de normalité Von Mises ==> Passer à la première page

Principe de la résolution numérique Soit la solution de Pour chaque élément Aux points

Principe de la résolution numérique Soit la solution de Pour chaque élément Aux points de Gauss La sol. Supposée élastique Si OK Si non il faut projeter sur De Schéma idéal schéma numérique Passer à la première page