10.07.2015 Views

L'approche ONERA pour le cumul des dommages en fatigue ... - SF2M

L'approche ONERA pour le cumul des dommages en fatigue ... - SF2M

L'approche ONERA pour le cumul des dommages en fatigue ... - SF2M

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

L’approche Onera <strong>pour</strong> <strong>le</strong> <strong>cumul</strong> <strong>des</strong><strong>dommages</strong> <strong>en</strong> <strong>fatigue</strong> sous chargem<strong>en</strong>tscomp<strong>le</strong>xes et température variab<strong>le</strong>J.L. ChabocheO.N.E.R.A., Châtillon, France(*) Jean-Louis.Chaboche@onera.fr2 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Pièces soumises à <strong>des</strong> sollicitationsthermomécaniques sévères3 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Analyses de durée de vie <strong>des</strong> structuresLOIS DECOMPORTEMENTStructureCALCULS DEσεLOISLOIS DESTRUCTURESD’ENDOMMAGEMENTPROPAGATIONchargem<strong>en</strong>t<strong>en</strong> serviceDcouplagedurée de vieà amorçagepropagationet ruptureMECANIQUE DE L’ENDOMMAGEMENT4 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Analyses de durée de vie <strong>des</strong> structuresLOIS DECOMPORTEMENTStructureCALCULS DEσεLOISSTRUCTURESD’ENDOMMAGEMENTchargem<strong>en</strong>t<strong>en</strong> servicedurée de vieà amorçagePost-traitem<strong>en</strong>t5 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Fatigue thermomécanique Lois de comportem<strong>en</strong>t : élasto-viscoplasticité Analyse inélastique cyclique Loi d’<strong>en</strong>dommagem<strong>en</strong>t interaction <strong>fatigue</strong>-fluage interaction <strong>fatigue</strong>-<strong>en</strong>vironnem<strong>en</strong>t Problème de la dép<strong>en</strong>dance <strong>en</strong> températureExposé ess<strong>en</strong>tiel<strong>le</strong>m<strong>en</strong>t basésur <strong>le</strong>s travaux et approchesdéveloppées par l’<strong>ONERA</strong><strong>en</strong> coopération avec Snecma / Turboméca et EMP retour sur 40 années d’activité quelques nouvel<strong>le</strong>s propositions6 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Plan Plan de la l'exposéprés<strong>en</strong>tation Loi de comportem<strong>en</strong>t Modè<strong>le</strong> de <strong>fatigue</strong> pure à <strong>cumul</strong> non linéaire Modè<strong>le</strong> d'interaction Fatigue-Fluage Modè<strong>le</strong> d'interaction Fatigue-Oxydation-Fluage Applications à <strong>des</strong> élém<strong>en</strong>ts structurauxsous chargem<strong>en</strong>ts thermomécaniques cycliques Prise <strong>en</strong> compte de la température variab<strong>le</strong>7 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Loi de comportem<strong>en</strong>t élasto-viscoplastique unifiéepot<strong>en</strong>tiel :Ω =Kn +1J( σ − X )−KR−kn+1Chaboche, 1977contrainte seuilnormalité :multiplicateur :écrouissagecinématique :p ∂Ω 3 σ'−X'& ε = = p&=∂σ2 J8 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012vp&n( σ − X )( σ − X ) − R − kn2 p p∂ΩJp& = ==∂σKsuperposition de l’écrouissage isotrope :• sur <strong>le</strong> seuil :X= ∑iX iR• sur <strong>le</strong> modu<strong>le</strong> d’écrouissage : γ = γ (R)• sur la contrainte de « trainée » : K = K (R)contrainte de trainée& ε3: & εArmstrong-Frederick, 1966 :X&2 pi= C3i& ε − γ iXip&restauration dynamiquemodu<strong>le</strong> d’écrouissage


Loi de comportem<strong>en</strong>t viscoplastique du monocristalmodè<strong>le</strong> de Cail<strong>le</strong>taud (1987)critère de Schmid :<strong>en</strong> viscoplasticité :norma<strong>le</strong> ausn r plan defs= τs− xs− rsglissem<strong>en</strong>tsnms= 21rrrrl r(s s s s)n ⊗ l + l ⊗ndirection deglissem<strong>en</strong>tv & =ssfKsss& γ = v&Sg( τ − xss)contrainte résolue :r rs s sτ = n . σ . l = m : σsfs= τs− rs≤ 0x&matrice d’écrouissage(interactions <strong>en</strong>tre systèmes)rrs= c & γ −dssrss= Q ∑ h 1−& γs( exp( −bvs ) )9 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Essais de <strong>fatigue</strong> thermomécanique1200σ (MPa)εm750superalliagemonocristallin AM11100 T° C6000-800-1 0 ε (%) 1Essais : AM1 (essais ENSMP )Models :950macroscopiqueviscoplasticité cristalline1600 σ (MPa)1600σ (MPa)750110009500600-1200-0.8 0 ε (%) 0.8-800 -0.8 0 ε (%)0.610 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Plan Plan de la l'exposéprés<strong>en</strong>tation Loi de comportem<strong>en</strong>t Modè<strong>le</strong> de <strong>fatigue</strong> pure à <strong>cumul</strong> non linéaire Modè<strong>le</strong> d'interaction Fatigue-Fluage Modè<strong>le</strong> d'interaction Fatigue-Oxydation-Fluage Applications à <strong>des</strong> élém<strong>en</strong>ts structurauxsous chargem<strong>en</strong>ts thermomécaniques cycliques Prise <strong>en</strong> compte de la température variab<strong>le</strong>11 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


σσ 1Cumul <strong>des</strong> <strong>dommages</strong> – Effets de séqu<strong>en</strong>ce dans<strong>le</strong>s essais à 2 niveauxN 1N 2σ 2NAcier 4340Règ<strong>le</strong> de Minerniveaux décroissants :NNf11+NNf22≤ 1niveaux croissants :NNf11+NNf22≥ 112 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Bases du <strong>cumul</strong> non linéaire <strong>des</strong> <strong>dommages</strong>dDdNexemp<strong>le</strong> : <strong>fatigue</strong> (uniaxia<strong>le</strong>) = F ( σ , σ D)f a ,1Dbase de nombreuses théoriesdont <strong>le</strong> modè<strong>le</strong> différ<strong>en</strong>tiel Oneravariab<strong>le</strong>s non séparab<strong>le</strong>sσ a1,σ 1σ a2,σ 2NNN1+2≤f N1 f 21D 10N1N f 113 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012N2N f 21NN f


Cumul linéaire <strong>des</strong> <strong>dommages</strong>exemp<strong>le</strong> : <strong>fatigue</strong> (uniaxia<strong>le</strong>)= F ( σ , σ ) φ(D)dDdNfa1variab<strong>le</strong>s séparab<strong>le</strong>sDσ a1,σ 1σ a2,σ 2Règ<strong>le</strong> de Miner :NNN1+2=f N1 f 21D 10N1N f 114 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012N2N f 21NN f


Cumul non linéaire du dommage : modè<strong>le</strong> Onera classiqueLoi d’évolution différ<strong>en</strong>tiel<strong>le</strong> :dD=Dα (...)⎛⎜⎝σ aM(...)⎞⎟⎠βdN(1978)Point clé :α fonction du chargem<strong>en</strong>tforme de la fonction alpha (1978) :1 −α=aσσuσa − l(...)−σmaxα = 1 sous la limite de <strong>fatigue</strong>effets de contrainte moy<strong>en</strong>ne cont<strong>en</strong>us dans <strong>le</strong>s fonctions M(…) et σ l (…)σl= σ1 −bσD( moy)M =M (1−b'σomoy)intégration <strong>en</strong>treD = 0 et D = 1(chargem<strong>en</strong>t périodique)N−β1 ⎛ σ a⎞ σ u − σ maxf ==1 − α ⎜⎜⎟⎟⎝ M(...) ⎠ a σ a lcourbes de Wöh<strong>le</strong>r avec deux asymptotes− σ (...)⎛⎜⎜⎝σ a⎞M(...)⎟⎟⎠−β15 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Cumul non linéaire du dommage :généralisation de la formulation OneradD=α(...)D1 −α1 −α =dN=N (...)fNN*ff(...)(...)Dα(...)NdN*f(...)loi « cachée »courbe SN de référ<strong>en</strong>cel’expression de N fpeut être quelconqueCHOIX LIBRE(2006)16 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012- Réunion <strong>SF2M</strong>-GT4, 6.07.2006- Chapitre "Cumul <strong>des</strong> <strong>dommages</strong>"La Fatigue <strong>des</strong> Matériaux et <strong>des</strong>Structures (Tome 4),Pineau-Bathias Eds, Hermès, 2008


Chargem<strong>en</strong>ts spectraux ou par blocs (Gassner)spectre Gaussi<strong>en</strong>,par blocstiré <strong>des</strong> travaux deC.M. Sonsino (L.B.F.)J.I.P. 2006, Paris17 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


En moy<strong>en</strong>ne, la « réalité » est <strong>en</strong>core plus « sévère »Damage sum :Dsum=∑nNif iBasée sur courbes SN à deux p<strong>en</strong>tes : γ et1,0000γ∗ = 2γ - 20,7500facteur5 ou 6safeProbability0,50000,2500experim<strong>en</strong>ts0,0000Eulitz, K.G., Kotte, K.L.,Materials Week 2000, Frankfurt, Sept. 2000cité par C. M. Sonsino, ASTM STP 1439, 200518 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 20120,01 0,1 1 10Damage sum


Prediction de ces résultats620 conditions :différ<strong>en</strong>ts niveaux de maxima,spectres linéaires ou Gaussi<strong>en</strong>sdiffér<strong>en</strong>tes p<strong>en</strong>tesdispersion naturel<strong>le</strong> prise <strong>en</strong> compte :14 % variation <strong>en</strong> contrainte1γSa1009080706050403020100Sa1001 10 100 1000 10000 100000N <strong>cumul</strong>é90807060504030201001 10 100 1000 10000 100000N <strong>cumul</strong>éProbability0,750,50,250β * = 0. 45 γ320 pointssimulationsexperim<strong>en</strong>ts0,0100 0,1000 1,0000 10,0000Damage sumLes modè<strong>le</strong>s de <strong>cumul</strong> nonlinéaire sont capab<strong>le</strong>s detraduire l’écart à la règ<strong>le</strong> deMinerChaboche et al., Darmstadt, 200919 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


A rec<strong>en</strong>t application to gigacyc<strong>le</strong> <strong>fatigue</strong> regime250Huang Ph.D. Thesis, 2010Huang et al., 2011Stress Amplitude (MPa)200Giga Simp<strong>le</strong>,R=-1Prior LCF 0.62% 100Cyc<strong>le</strong>s,giga,R=-1Chaboche model,0.62%,100cyc<strong>le</strong>s+GigaChaboche model,0.34%,300cyc<strong>le</strong>s+GigaChaboche model,0.47%,200cyc<strong>le</strong>s+GigaPrior LCF 0.47% 200Cyc<strong>le</strong>s,giga,R=-1Young's modulusmeasurem<strong>en</strong>tFrequ<strong>en</strong>cy changemeasurem<strong>en</strong>tPrior LCF 0.34%,300Cyc<strong>le</strong>s,giga,R=-11501,0E+05 1,0E+06 1,0E+07 1,0E+08 1,0E+09 1,0E+10Cyc<strong>le</strong>sDA 42 C-Mn steelPrior LCF strain controlLCF0.34, 0.47, 0.62 % strain amplitudeN / N fVHCF20 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Plan Plan de la l'exposéprés<strong>en</strong>tation Loi de comportem<strong>en</strong>t Modè<strong>le</strong> de <strong>fatigue</strong> pure à <strong>cumul</strong> non linéaire Modè<strong>le</strong> d'interaction Fatigue-Fluage Modè<strong>le</strong> d'interaction Fatigue-Oxydation-Fluage Applications à <strong>des</strong> élém<strong>en</strong>ts structurauxsous chargem<strong>en</strong>ts thermomécaniques cycliques Prise <strong>en</strong> compte de la température variab<strong>le</strong>21 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Interaction non-linéaire <strong>fatigue</strong> fluage (modè<strong>le</strong> FatFlu)fluagedDdt=F( σ D)c ,<strong>fatigue</strong>dDdN=F( σ , σ D)f a ,1évolutions différ<strong>en</strong>tesDfluage<strong>fatigue</strong>D 10t1t c 1N 2Nf21NN ftt c22 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Interaction non-linéaire <strong>fatigue</strong> fluage (modè<strong>le</strong> FatFlu)fluagedDdtc=Fc( σ , D )c<strong>fatigue</strong>dDfdN=Ff( σ , σ , D )afac<strong>cumul</strong>ation <strong>des</strong> deux <strong>dommages</strong> à chaque cyc<strong>le</strong> :hypothèse : interaction additive :( σ , D , D ) dt + F ( σ , σ , D D )dNdD = dDc + dD f = Fcc f f a f ,effets d’interactioncdD( σ , D + D ) dt + F ( σ , σ D D )dN= dDc + dD f = Fcc f f a , f +cdD =intégration du fluage à chaque cyc<strong>le</strong> := prévision fluage purN cN f= prévision <strong>fatigue</strong> pureFD( σ , D) dt + F ( σ , σ D)dNc f a ,φcdD =N( D) φ f( D) dN + dNcNf23 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Méthode <strong>des</strong> « tranches planes » : depuis <strong>le</strong>s années 70 Contraintes d’origine thermique Hypothèse de Bernouilli (poutre)zM xFM y Découpage <strong>en</strong> élém<strong>en</strong>ts uniaxiauxσ zxε z= A + Bx+CyyTε zlibreε z=A+CycompressioncompressiontractionBord d’attaqueBord de fuite0Bord d’attaqueyBord de fuite24 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Validation de la loi de <strong>cumul</strong> non linéaire <strong>fatigue</strong>-fluageSNECMAsuperalliagebase NickelIN 100M53 1st stage turbine blade (15 channels)tested in <strong>ONERA</strong> lab (1973-74)25 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Validation de la loi de <strong>cumul</strong> non linéaire <strong>fatigue</strong>-fluageSNECMAN exp410superalliagebase NickelIN 1001973 - 1984310id<strong>en</strong>tificationprévisions(isothermes)superalliagebase NickelIN 1001 er étage aubes de turbine M53(15 canaux) essais <strong>ONERA</strong> (1973)10 22prévisions(anisothermes)composants2 410 10 3 10 calcN26 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Validation de la loi de <strong>cumul</strong> non linéaire <strong>fatigue</strong>-fluageN exp410superalliagebase NickelIN 1001973 - 1984310id<strong>en</strong>tification10 2prévisions(isothermes)prévisions(anisothermes)1 er étage aubes de turbine M53(15 canaux) essais <strong>ONERA</strong> (1973)10 210 3 410Ncalc27 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 20122composants


Validation de la loi de <strong>cumul</strong> non linéaire <strong>fatigue</strong>-fluageN exp410superalliagebase NickelIN 1001973 - 1984310id<strong>en</strong>tificationPlaque IN100 avec gradi<strong>en</strong>tsthermiques cycliques (Lesne, 1984)10 22prévisions(isothermes)prévisions(anisothermes)composants2 410 10 3 10 calcN28 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Validation de la loi de <strong>cumul</strong> non linéaire <strong>fatigue</strong>-fluagemesures DDV amorçage29 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Plan Plan de la l'exposéprés<strong>en</strong>tation Loi de comportem<strong>en</strong>t Modè<strong>le</strong> de <strong>fatigue</strong> pure à <strong>cumul</strong> non linéaire Modè<strong>le</strong> d'interaction Fatigue-Fluage Modè<strong>le</strong> d'interaction Fatigue-Oxydation-Fluage Applications à <strong>des</strong> élém<strong>en</strong>ts structurauxsous chargem<strong>en</strong>ts thermomécaniques cycliques Prise <strong>en</strong> compte de la température variab<strong>le</strong>30 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Modè<strong>le</strong> <strong>fatigue</strong> – oxidation – fluageFAT.OX.FLUF. Gal<strong>le</strong>rneau, 1995phénoménologie :micro-initiationFatiguemicro-propagationdD i= f i(...)dNdD f (...)dNp=psurfacedDOxydationox =fox(...)dtdDeffets du tempsFluagec =fc(...)dtN = N + NRipmicro-initiationinteractionFatigue-Oxydationmicro-propagationinteractionFatigue-FluageAM1 revêtu C1A31 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


2∆σ ∆σ ∆σFatOxFlu : AM1, 1100°C, R = -1500Courbes calculéesid<strong>en</strong>tification300/ 2 (MPa)∆σ / 250 Hz200 Essais :50 Hz fluageoxydation5 Hz5 Hz0.5 Hzmainti<strong>en</strong> t 0.5 Hzhold=90s10010 0 10 2 10 4 106Nombre de cyc<strong>le</strong>s à rupture32 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


2∆σ ∆σ ∆σFatOxFlu : AM1, 950°C, R = -1850prévisionscourbes calculées500tmns lqssllslfacteurs d’Arrhénius<strong>pour</strong> l’influ<strong>en</strong>ce dela température sur <strong>le</strong>fluage et l’oxydation/ 2 (MPa)∆σ / 2300Essais :50 Hz5 Hz0.5 Hz0.05 Hzmainti<strong>en</strong>tqt ttt hold=90s.sq0.05 Hzl50 Hzs l m5 Hz0.5 Hz20010 0 10 1 10 2 10 3 t10 4 10 5 10 6 10 7Nombre de cyc<strong>le</strong>s à rupture33 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


FatOxFlu : monocristal AM1, validation globa<strong>le</strong>10 8610N expessais <strong>pour</strong> id<strong>en</strong>tificationprévision essais isothermesprévision essais TMF410Essais TMFuniaxiaux102N100calc10010 2 10 4 10 6 10 8Nombre de cyc<strong>le</strong>s à ruptureEssais par Snecma,Onera,ENSMP,ENSMA Poitiers,CEAT34 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Plan Plan de la l'exposéprés<strong>en</strong>tation Loi de comportem<strong>en</strong>t Modè<strong>le</strong> de <strong>fatigue</strong> pure à <strong>cumul</strong> non linéaire Modè<strong>le</strong> d'interaction Fatigue-Fluage Modè<strong>le</strong> d'interaction Fatigue-Oxydation-Fluage Applications à <strong>des</strong> élém<strong>en</strong>ts structurauxsous chargem<strong>en</strong>ts thermomécaniques cycliques Prise <strong>en</strong> compte de la température variab<strong>le</strong>35 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Essais de <strong>fatigue</strong> thermomécaniqueavec gradi<strong>en</strong>t de paroirefroidissem<strong>en</strong>t par airtempératurechauffageinductioncontrainteTempératures stationnaires :paroi externe :paroi interne :T i nt = 900°Pacou, Poirier, Nouailhas, 199736 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012T ext = 1050°CC


Essais de <strong>fatigue</strong> thermomécanique avec gradi<strong>en</strong>t de paroitempérature :7206606001000800∆σ∆σ/ 2 (MPa)testspredictions540480600AM1 420IN 100360400300240180120200loi de comportem<strong>en</strong>tet d’<strong>en</strong>dommagem<strong>en</strong>tdéterminée <strong>en</strong> 197710 2 103Pacou, Poirier, Nouailhas, 1997N R41037 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Essais de <strong>fatigue</strong> thermomécanique avecgradi<strong>en</strong>t de paroi et trous de refroidissem<strong>en</strong>text<strong>en</strong>somètrethermocoup<strong>le</strong>sspiresinduction HFsimulation d’un cyc<strong>le</strong> stabilisécontrainte - déformationEssaisOneraPacou, Poirier, 2001OnzetrousSimulation detempérature(au maximum du cyc<strong>le</strong>)Thèse J.-M. Cardona, ENSMP, 20020 500 100038 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Essais de <strong>fatigue</strong> thermomécanique avecgradi<strong>en</strong>t de paroi et trous de refroidissem<strong>en</strong>t500∆σ/2 Mpa400300200100RuptureAmorçage InitiationEssais essaisRupture Initiation calculée Modè<strong>le</strong>0100 1000 Nombre de Ncyc<strong>le</strong>sR10000--> prévisions effet du gradi<strong>en</strong>t trop conservativeslocalPrévision de l’initiation à la position correcte (trou 2 à gauche)39 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Gradi<strong>en</strong>t effect for the crack initiation at ho<strong>le</strong>PhD Thesis Myriam KAMINSKI2007, Oneramaximal stress for crackinitiation analysisho<strong>le</strong>0.424 mmN R Ø0,4> N R Ø2> N R Ø4to attain a crack <strong>le</strong>ngth of300 µmdistance to ho<strong>le</strong> (mm)40 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Plate with ho<strong>le</strong>s of various sizes AM1 – 950°CΦ = 4, 2, 1, 0.4 mm and slits1000R = 0.1 - 5 Hzσ nom, maxsmooth5001 mm2 mm4 mm0.4 mmslits3 combined effects :gradi<strong>en</strong>t effectsca<strong>le</strong> effect (porosities)multiaxiality effect (thickness) 100100 1000 10000 100000 1000000 1000000041 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 201220010 2 10 3 10 4 10 5 10 6 N init 10 7


Plate with ho<strong>le</strong>s of various sizes AM1 – 950°CΦ = 4, 2, 1, 0.4 mm and slits1000R = 0.1 - 5 Hzσ nom, maxsmooth5001 mm2 mm4 mm0.4 mmslits20010 2 10 3 10 4 10 5 10 6 N init 10 742 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012Model predictions – local maximum100100 1000 10000 100000 1000000 10000000


Plate with ho<strong>le</strong>s of various sizes AM1 – 950°CΦ = 4, 2, 1, 0.4 mm and slits1000R = 0.1 - 5 Hzideal smooth (no porosity)σ nom, maxsmooth5001 mm2 mm4 mm0.4 mmslitsModel predictions – local average :φ( x ) =⎛exp⎜⎜⎜⎝⎛−⎜⎜⎝x ⎞d ⎟⎟⎠2⎞⎟⎟⎟⎠43 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 201220010 2 10 3 10 4 10 5 10 6 N init 10 7D(x ) =Ωφ(x ξ )dV100Vd = 120µm100 1000 10000 100000 1000000 10000000Ω( x ) = ∫ − ξ1(x )∫φ(x −ξ)D( ξ )dξ


Fatigue thermomécaniquemultiaxia<strong>le</strong> de l’AM1Ph.D. V. Bonnand, 2006réponsecisail<strong>le</strong>m<strong>en</strong>t / températuretraction-torsion sur tubesσ mises44 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Fatigue thermomécanique multiaxia<strong>le</strong> de l’AM1validation du modè<strong>le</strong> de viscoplasticité cristalline de l’Eco<strong>le</strong> <strong>des</strong> Minesessai de traction/torsion : ∆ε m /2=0,49 % - ∆γ/2=0,425 %plasticstraintotalstraintestsmodelplasticstraintotalstraintraction-compressiontorsioncyc<strong>le</strong> stabilisé: contraintes-déformations globa<strong>le</strong>s (moy<strong>en</strong>nes)45 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Plan Plan de la l'exposéprés<strong>en</strong>tation Loi de comportem<strong>en</strong>t Modè<strong>le</strong> de <strong>fatigue</strong> pure à <strong>cumul</strong> non linéaire Modè<strong>le</strong> d'interaction Fatigue-Fluage Modè<strong>le</strong> d'interaction Fatigue-Oxydation-Fluage Applications à <strong>des</strong> élém<strong>en</strong>ts structurauxsous chargem<strong>en</strong>ts thermomécaniques cycliques Prise <strong>en</strong> compte de la température variab<strong>le</strong>46 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Applications à température variab<strong>le</strong> et missions comp<strong>le</strong>xesInterpolation <strong>des</strong> coeffici<strong>en</strong>ts <strong>en</strong> fonction de la températureTempérature variab<strong>le</strong> – missions comp<strong>le</strong>xes :Fluage : pas de problème, il suffit d’intégrerFatigue : modè<strong>le</strong> non incrém<strong>en</strong>tal, “<strong>en</strong> cyc<strong>le</strong>s” Technique de sé<strong>le</strong>ction <strong>des</strong> cyc<strong>le</strong>s : Rain Flow (multiaxial) Prise <strong>en</strong> compte de la température variab<strong>le</strong>NormalisationMéthode de type Taira Cyc<strong>le</strong>s comp<strong>le</strong>xes non proportionnelsPrise <strong>en</strong> compte <strong>des</strong> effets de déphasageForme “incrém<strong>en</strong>ta<strong>le</strong> simplifiée”47 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Démarche basée sur la définition d’une température effective1. Normalisation, sans imposer <strong>des</strong> exposants id<strong>en</strong>tiques :S( t)=σ ( t )σ * (T( t))t =2. Technique du Rain-Flow sé<strong>le</strong>ction <strong>des</strong> cyc<strong>le</strong>s : { i} i 1, n3. Définition de l’amplitu<strong>des</strong> et de la va<strong>le</strong>ur moy<strong>en</strong>ne de S( t )amplitude :« c<strong>en</strong>tre » associé :amplitude effective :S aeqSSa effméthode du « range »4. Définition de la température effective*T5. Définition de l’amplitude effective <strong>pour</strong> chargem<strong>en</strong>t non proportionnel6. Utilisation dans la loi d’<strong>en</strong>dommagem<strong>en</strong>t à <strong>cumul</strong> linéaire ou non linéaire48 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Démarche basée sur la définition d’une température effective1. Normalisation, sans imposer <strong>des</strong> exposants id<strong>en</strong>tiques :S at =2. Technique du Rain-Flow sé<strong>le</strong>ction exposant <strong>des</strong> cyc<strong>le</strong>s variant : { i} i 1, nde façon monotone<strong>en</strong> fonction de TS( t3. Définition de l’amplitu<strong>des</strong> et de la va<strong>le</strong>ur moy<strong>en</strong>ne de S( t ))=σ ( t )σ * (T( t))amplitude : 1« c<strong>en</strong>tre » associé :S aeqSméthode du « range »amplitude effective :Sa eff4. Définition de la température effective*T10 4 log N f5. Définition de l’amplitude effective <strong>pour</strong> chargem<strong>en</strong>t non proportionnel6. Utilisation dans la loi d’<strong>en</strong>dommagem<strong>en</strong>t à <strong>cumul</strong> linéaire ou non linéaire49 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Définition de la température effective*Ts&( t)Simplification du <strong>cumul</strong> linéaire sur cette étapeIllustré sur <strong>le</strong> cas d’une loi puissance à seuil :S− S1 a D=eff*f B (T )N(T )β (T )SJS( t)En isotherme, proportionnel, cette loi <strong>en</strong> cyc<strong>le</strong>s peut résulter de l’intégration de :dDdtβ ( T(t))=*4 B ( T(t))J( S(t)− S )B*− SD( T(t))( T(t))β ( T(t))−1s&( t)avec :3( S( & t )) S' & ( t ) : S' & ( t )s &( t ) = J =déviateur de S &2forme incrém<strong>en</strong>ta<strong>le</strong> quiservira aussi <strong>pour</strong> <strong>le</strong> chargem<strong>en</strong>t non proportionnel50 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Cas du chargem<strong>en</strong>t anisotherme proportionnelIntégration numérique sur <strong>le</strong> cyc<strong>le</strong> :dDdN1= =∫ N*f cyc<strong>le</strong>β(T( t ))4 B*(T( t ))J( S( t ) − S )B*−SD(T( t ))(T( t))β (T( t )) −1s( & t)dt*Intégration analytique avec T = T = Cte :dDdN=Nf1( T*)=SaeffB−*S(TD*(T)*)*β ( T)Id<strong>en</strong>tification :donne*Téquation non linéaire à une inconnue :- Newton- point fixe51 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Cas du chargem<strong>en</strong>t isotherme non proportionnelIntégration numérique sur <strong>le</strong> cyc<strong>le</strong> :dDdN1= =∫ N*f cyc<strong>le</strong>β4 B*J( S( t ) − S )B*−SDβ −1Intégration analytique <strong>en</strong> uniaxial proportionnel :dDdN=Nf1( SId<strong>en</strong>tification :*donne S aeff :*aS eff− SD=⎛β⎜⎜4⎝∫cyc<strong>le</strong>J*aeff)=( S( t ) − S )S*aeff− Ss( & t )dtSβ− SD*B1ββ − 1D s( & t ) dt⎟ ⎟ ⎞⎠Js&( t)S( t)52 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Cas du chargem<strong>en</strong>t anisotherme non proportionnelIntégration numérique sur <strong>le</strong> cyc<strong>le</strong> :dDdN1= =∫*N f ( S( t ))cyc<strong>le</strong>β(T( t ))4 B*(T( t ))J( S( t ) − S )B*− SD(T( t ))(T( t ))β(T(t )) −1s( & t )dt*Intégration avec T = T = Cte :dDdN=*S aeffNf1( S( t), TId<strong>en</strong>tification :*1)=∑J( S − S )obt<strong>en</strong>u <strong>en</strong>suite comme dans <strong>le</strong> cas isotherme⎡⎢⎢⎢⎣=i+1 i+1 i+14 + −*iJi1 JiBdonne*T− S(TD*(T*β ( Téquation non linéaire à une inconnue :Jchargem<strong>en</strong>t linéaire par morceaux+petite approximation))*)−Ji− SB*D*(T(T)- Newton- point fixe*)β (T*⎤⎥⎥⎥⎦)53 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012


Amplitude effective avec chargem<strong>en</strong>t non proportionnel*aS eff− SD=⎛β⎜⎜4⎝∫cyc<strong>le</strong>J( S( t ) − S )SDβ − 1s( &t ) dtCas particulier de la traction-torsion déphasée de 90° :−⎟ ⎟ ⎞⎠1βavec S D = 0S*aeff=Saeff⎛ π ⎞β ⎟ ⎜⎝ 2 ⎠1ββparamètrefixé2,0Modè<strong>le</strong> Onera (1993) :1,5S1( r + ) r*a = Sa1eff effrparamètreajustab<strong>le</strong>S*a_eff / Sa_eff1,00,5Onera (1993)approche proposée54 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 20120,00 5 10 15 20 25rr = β − 1


Application à l’INCO 718 à 550°C1000- INCO 718 550°C -essais Onera 1987essais Onera 2005amplitude de contrainte équiva<strong>le</strong>nte MPatractiontorsiontract-tors. <strong>en</strong> phasetract.tors. hors-phasetorsion 1987traction-torsion1987TC 88 - Rsig = -1TC 88 Reps =-1 ou 0hors phase 1987traction-compr. p<strong>le</strong>inesOnera Ne Re=-1 tm=90Snecma NE Re=-1 0.5HzTOS (mars 2006) recaléTOS limite de <strong>fatigue</strong>100100 1000 10000Nf100000 1000000 10000000*Saeff−SD⎛ β=⎜⎜⎝4∫cyc<strong>le</strong>J( S( t ) − S )55 Séminaire Commission Fatigue <strong>SF2M</strong> – 6 décembre 2012−SDβ<strong>en</strong> phasetractiontorsiontraction-torsionhors phase (90°)traction-torsionprévision avec :S D = 0prévision avec :S D = 35011β ⎛ S ⎞ β− 1π aeffs( t ) dt =⎜⎜βSa− S⎟ ⎟ ⎞&eff2 Sa− S⎟⎟⎠ ⎝ eff D⎠( )D

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!