Instabiele taaie scheurgroei in een gasturbineschijf
De vraag was hoe diep een scheur in de naaf van een gasturbineschijf kan groeien voordat die bij vol toerental instabiel gaat scheuren: circa 12–14 mm bij 3600 tpm, ruim boven wat ultrasoon onderzoek betrouwbaar detecteert. Een breuktaaiheidsproef volgens ASTM E1820 op rotorstaal, gekoppeld aan een elastisch-plastisch scheurmodel in Abaqus via het raakcriterium voor de J-integraal.
- Project
- Casestudy, Politecnico di Milano
- Mijn rol
- Taaiheidsproef, FEA en verificatie
- Software
- Abaqus/Standard
- Methode
- Raakcriterium voor de J-integraal (Paris–Tada–Zahoor–Ernst)
- ≈14,0mmkritieke beginscheurdiepte in de naaf vóór instabiele taaie scheurgroei, 3600 tpm
- 12,0mmmet de experimentele R-curve (−14%)
- 167,2kJ/m²geldige initiatiebreuktaaiheid, JIc (ASTM E1820)
- ±2%J in het diepste punt, contouren 4–10
Aangetoonde vaardigheden
Aanpak
Voor een halfronde fout die wordt aangenomen in de boring van een roterende gasturbineschijfnaaf de beginscheurdiepte bepalen waarboven stabiele taaie scheurgroei bij nominaal toerental instabiel wordt, met een breuktaaiheidscurve die met een fysieke proef is vastgesteld in plaats van een aangenomen materiaaleigenschap.
- De elastisch-plastische breuktaaiheid (JIc) en de scheurgroeiweerstandscurve (J-R) van het schijfmateriaal experimenteel bepalen volgens ASTM E1820, met de methode van elastisch ontlasten en herbelasten op één CT-proefstuk.
- Een elastisch-plastisch FE-model in twee stappen opbouwen in Abaqus van een halfronde scheur in de naafboring onder centrifugale belasting bij nominaal toerental, met een globaal schijfmodel en een lokaal verfijnd submodel, en de J-integraal bepalen met de domein- (contour)integraalmethode.
- Het raakcriterium van Paris–Tada–Zahoor–Ernst toepassen op de resulterende curve van de scheurdrijvende kracht om de kritieke beginscheurgrootte te berekenen.
- Het FE-resultaat verifiëren met een onafhankelijke analytische handberekening en een controle op padonafhankelijkheid van de J-integraal, en de gemeten taaiheid valideren met gepubliceerde gegevens voor vergelijkbare rotorstalen.
- De twee delen van de studie verbinden met een gevoeligheidscontrole: de berekening van de kritieke scheurgrootte herhalen met de eigen R-curve uit één proefstuk in plaats van de referentiecurve.
-
Deel 1: de breuktaaiheid beproeven
ASTM E1820
Een standaard CT-proefstuk (compact tension) van martensitisch roestvast rotorstaal met 9–12% chroom, zoals gebruikt in de energietechniek, is met vermoeiing voorgescheurd, van zijgroeven voorzien en beproefd volgens de methode van ASTM E1820 met elastisch ontlasten en herbelasten op één proefstuk, zodat één proef een volledige J-R-curve oplevert.
Het proefstuk is quasistatisch belast met periodieke gedeeltelijke ontlast- en herbelastcycli; de elastische compliantie van elke cyclus wordt omgerekend naar een momentane scheurlengte, zodat één proefstuk de volledige reeks J–Δa-punten levert waarvoor anders zes proefstukken nodig zouden zijn.
- J berekend als som van de elastische en plastische component in elk ontlastpunt (ASTM E1820, bijlage A1).
- Machtsfunctie als weerstandscurve gefit op het gekwalificeerde databereik: J = C1ΔaC2.
- JIc bepaald op het snijpunt van de fit met de constructielijn met 0,2 mm offset.
- Fractografie in 9 punten na de proef om de op compliantie gebaseerde schattingen van de scheurlengte te verifiëren.
Gefitte weerstandscurve: J = 29,49 Δa0,642 (kJ/m², mm), gekwalificeerd over Δa = 0,239–1,721 mm. Met deze taaiheid ligt het materiaal midden in het gepubliceerde bereik voor vergelijkbare rotorstalen; de waarde wordt meegenomen als onafhankelijke, op beproeving gebaseerde controle van de kritieke scheurgrootte van de schijf.
- Geldige initiatietaaiheid, JIc
- 167,2 kJ/m²
- Vloeigrens / treksterkte
- 703 / 836 MPa
- Proefstuk W, B / Bn
- 50,8 mm, 25,4 / 20,4 mm
- Beginscheur, a0
- 33,78 mm (a0/W = 0,665)
Technische tekening van het CT-proefstuk (ASTM E1820): hoofdmaten, GD&T-toleranties, detail van de chevronkerf en isometrische weergave van de bewerkte helften. Selecteer om te vergroten. 
Representatieve elektromechanische beproevingsmachine (klasse Instron Model 640), ingericht voor beproeving volgens ASTM E1820, met een clip-gauge-extensometer over de scheurmond. 
Weerstandscurve J–Δa: proefgegevens, constructielijnen volgens ASTM E1820 en regressie als machtsfunctie. Selecteer om te vergroten. Materiaaleigenschappen en afmetingen van het proefstuk Eigenschap / afmeting Symbool Waarde Elasticiteitsmodulus E 206.817 MPa Vloeigrens (0,2%-rekgrens) σYS 703 MPa Treksterkte σUTS 836 MPa Breedte proefstuk W 50,8 mm Dikte / netto dikte (met zijgroeven) B / Bn 25,4 / 20,4 mm Beginlengte (vermoeiingsvoorscheur) a₀ 33,78 mm (a₀/W = 0,665) -
Het schijfmateriaal kiezen
Materiaalkeuze
Voordat het beproefde staal als schijfmateriaal wordt gekozen, wordt het vergeleken met de twee andere realistische materiaalklassen voor een roterende turbineschijf, op sterkte, taaiheid, bedrijfstemperatuur en kosten.
De martensitische roestvaste kwaliteit met 9–12% Cr is gekozen: het is hier de enige kandidaat met een direct beproevingsresultaat in plaats van een aangenomen eigenschap, ze dekt ruim de verhoogde temperatuur in het FE-model en ze biedt een duidelijk betere corrosiemarge dan het laaggelegeerde staal, tegen een fractie van de meerprijs van een nikkelsuperlegering. Die laatste blijft de aangewezen upgrade als een hetere of snellere toepassing dat vraagt.
Vergeleken schijfmaterialen Eigenschap Laaggelegeerd NiCrMoV Roestvast staal 9–12% Cr (beproefd) Nikkelsuperlegering (klasse IN718) Vloeigrens (kamertemperatuur) ≈ 600–800 MPa 703 MPa (gemeten) ≈ 1030–1170 MPa Breuktaaiheid, JIc ≈ 150–300 kJ/m² 167 kJ/m² (gemeten) ≈ 80–150 kJ/m² Max. bedrijfstemperatuur ≈ 450–480 °C ≈ 580–600 °C ≈ 650–700 °C Bestendigheid tegen corrosie / oxidatie Matig Goed Zeer goed Relatieve materiaal- en bewerkingskosten Referentie (1,0×) ≈ 1,2–1,4× ≈ 4–6× -
Deel 2: de gescheurde schijf modelleren
Abaqus/Standard
Een vereenvoudigd maar maatvoeringstechnisch representatief axisymmetrisch schijfprofiel (verdikte naaf, taps toelopend lijf, verdikte rand) wordt bij nominaal toerental geanalyseerd met een halfronde scheur die in de naafboring wordt aangenomen, met een globaal model dat lokaal wordt verfijnd door een gestuurd submodel.
Beide modellen zijn opgebouwd als 2D-axisymmetrische meshes en daarna geroteerd tot een 3D-sector met losse mesh (orphan mesh) om de plaatsafhankelijke J-integraal langs het halfronde scheurfront te bepalen. De snijvlakken van het submodel worden aangestuurd met verplaatsingen uit de geconvergeerde globale oplossing: de gangbare werkwijze om een singulariteit aan de scheurtip op te lossen zonder het hele componentmodel te verfijnen.
- Globale mesh: elementen van 8 mm, verfijnd tot 4 mm bij de scheur, met een gefocuste ring van 24 elementen rond de scheurtip.
- Elastisch-plastisch materiaal: E = 170,5 GPa, σy = 550 MPa (verstevigingstabel met ware spanningen), representatief voor een metaaltemperatuur van de schijf van ≈520–540 °C.
- Belasting: alleen een rotatievolumekracht bij ω = 377 rad/s; scheurvlakken vrij gelaten, symmetrie op de snijvlakken van de geroteerde sector.
- J-integraal opgevraagd als contour- (domein)integraal over 10 opeenvolgende ringen, zowel in het diepste punt als in het oppervlaktepunt van het scheurfront.
- Nominaal toerental
- 3600 tpm (377 rad/s)
- Straal van boring tot rand
- 252–960 mm
- Aangenomen scheur
- 10 mm, a = c
- Model
- 2 stappen: globaal model + submodel
Geometrie van de schijf: isometrische weergave, doorsnede van de aangenomen scheur in de naafboring (Ø10 mm) en gemaatvoerde halve doorsnede (naafhoogte 250 mm, afrondingen R50/R60, randhoogte 170 mm). Selecteer om te vergroten. 
Globaal model: 3D-orphan-mesh met aangebrachte belasting en symmetrierandvoorwaarden op de snijvlakken. 
Submodel: lokaal verfijnde mesh met belasting op het scheurfront en door het globale model aangestuurde randvoorwaarden. 
Gefocuste ringmesh rond de scheurtip, 24 elementen rond de tip, die de singulariteit oplost. Selecteer om te vergroten. -
De kritieke scheurgrootte bepalen
Raakcriterium
De curve van de scheurdrijvende kracht, uit drie submodeloplossingen bij beginscheurdieptes van 10, 15 en 20 mm, wordt vergeleken met de weerstandscurve van het materiaal volgens het instabiliteitscriterium van de raaklijn (Paris–Tada–Zahoor–Ernst).
De referentie-R-curve wordt langs de as van de scheurdiepte verschoven (stippellijn in de grafiek hieronder) tot ze precies raakt aan de curve van de drijvende kracht. Die verschuiving, en niet het raakpunt zelf, is de kritieke beginscheurgrootte: a0,cr = 13,99 mm.
Het raakpunt zelf ligt op een grotere diepte, a* = 16,21 mm (J = 75,6 kJ/m²), wat binnen 0,1% overeenkomt met de oorspronkelijke grafische constructie van verschuiven en vergelijken. Onder a0,cr blijft de drijvende kracht over het hele groeipad onder de weerstandscurve, zodat elke kleine scheurgroei vanzelf tot stilstand komt.
- Kritieke beginscheurdiepte
- ≈14,0 mm (13,99 mm)
- Raakpunt
- a* = 16,21 mm, J = 75,6 kJ/m²
Scheurdrijvende kracht uit het submodel, diepste punt Beginscheurdiepte, a (mm) J-integraal, diepste punt (kJ/m²) 10 23,79 15 60,0 20 142,0
Constructie volgens het raakcriterium
Scheurdrijvende kracht tegen de scheurdiepte, uitgezet tegen de referentie- en de eigen R-curve. Beweeg over een markering voor details.
Controle van de resultaten
Verificatie bevestigt dat het model correct is doorgerekend; validatie bevestigt dat de materiaalinvoer overeenkomt met onafhankelijke, gepubliceerde gegevens. Beide zijn nodig voordat de kritieke scheurgrootte een technische conclusie kan dragen.
- Domeinonafhankelijkheid: ±2% in het diepste puntJ stabiliseert vanaf contour 4 (23,64–23,87 kJ/m²), wat bevestigt dat de mesh de scheurtip oplost zonder contourafhankelijke drift.
- Globaal model tegenover submodel: 3,2% overeenstemmingDe maximale spanning stijgt bij lokale verfijning van 698 naar 721 MPa: de verwachte grootte en richting voor een goed opgezet submodel.
- Handberekening: juiste orde en richtingDe analytische berekening geeft σθ(boring) ≈ 850 MPa tegenover een FEM-waarde van 340–450 MPa: een overschatting met de factor die te verwachten is als de verdikking van de naaf wordt verwaarloosd, precies zoals bedoeld.
- Validatie met literatuur: 167 kJ/m² ligt midden in het bereikDe gepubliceerde JIc voor vergelijkbare rotorstalen ligt ruwweg tussen 100 en 300 kJ/m²; de gemeten waarde ligt in het midden, passend bij een normaal taai materiaal.
| Kwalificatiecontrole volgens ASTM E1820 | Verkregen waarde | Resultaat |
|---|---|---|
| Regressie-exponent, C₂ < 1 | C₂ = 0,642 | Vervuld |
| Dikte, B > 10 JQ/σY | 25,4 mm > 1,97 mm | Vervuld |
| Ligament, b₀ > 10 JQ/σY | 17,02 mm > 1,97 mm | Vervuld |
| Verhouding scheurlengte, 0,45 ≤ a₀/W ≤ 0,70 | a₀/W = 0,665 | Vervuld |
Technisch oordeel en aanbevelingen
Naast het getal voldoet/voldoet niet laat de analyse zien waar het ontwerp marge heeft en waar een werkelijke beslissing over geschiktheid voor gebruik meer nodig heeft dan deze studie met één belastinggeval biedt.
-
Ruime marge boven de detectiegrens van NDO
Een kritieke fout van 12–14 mm ligt ruim boven de detectiegrens van enkele millimeters waarmee de ultrasone en oppervlaktemethoden voor niet-destructief onderzoek die gewoonlijk op rotorboringen worden toegepast betrouwbaar detecteren. Dat ondersteunt verdere inzet onder een periodiek inspectieprogramma dat op die marge is afgestemd.
-
Aandachtspunt
Overtoeren verkleinen de marge; nog niet gekwantificeerd
De centrifugale spanning schaalt met ω² en de elastische drijvende kracht J met het kwadraat van de spanning, dus J schaalt met ω4: 10% overtoeren bij een proef vergroot J met circa 1,46× en verkleint de kritieke scheurgrootte. Een volledige kwantificering vraagt om het herhalen van de raakconstructie bij overtoeren; dat wordt aanbevolen als vervolgwerk.
-
Aandachtspunt
Baseer een werkelijke beslissing op een R-curve uit meerdere proefstukken
De gevoeligheid van 14% (analyse 3 in de technische details) tussen twee afzonderlijk verdedigbare weerstandscurves betekent dat een proef met één proefstuk, hoe goed ook gekwalificeerd, op zichzelf geen beslissing over de levensduur of het inspectie-interval van een werkelijke rotor mag dragen.
-
Scope: één statisch belastinggeval
De analyse betreft één constante toestand van 3600 tpm met alleen centrifugale belasting: geen thermische of montagespanningen, geen wisselende start-stopcycli, en een geïdealiseerde in plaats van een werkelijk aangetroffen foutgeometrie.
Technische details
De drie analysestappen achter de gerapporteerde kritieke scheurgrootte.
Details van de analyse Spanningen in globaal model en submodel, bepaling van de J-integraal en gevoeligheid voor de bron van de R-curve
Elke stap van de FE-werkwijze, van globaal model tot gevoeligheidsstudie, is op zijn beurt geverifieerd voordat de kritieke scheurgrootte als gerapporteerd resultaat wordt beschouwd.
01. FE-spanningsresultaten: globaal model en submodel
Doel
De volledige schijf onder belasting bij nominaal toerental doorrekenen voor het algemene spanningsveld en een eerste J-integraal, en daarna bij de scheurtip lokaal verfijnen met een submodel dat door de globale oplossing wordt aangestuurd, zonder de kosten van het verfijnen van de hele schijf.
De stijging van 3,2% bij verfijning, in de verwachte richting, is wat een goed opgezet submodel hoort te laten zien, en geen inconsistentie tussen beide modellen.
02. Bepaling van de J-integraal en de curve van de drijvende kracht
Doel
Bevestigen dat de J uit de domeinintegraal over 10 contourringen mesh-geconvergeerd is, en daarna de drie submodeloplossingen (a = 10, 15, 20 mm) fitten tot een doorlopende curve van de scheurdrijvende kracht, CDF(a), voor de raakconstructie in stap 4.
CDF(a) = 0,9157 a² − 15,65 a + 88,72 (a in mm, J in kJ/m²)
Het diepste punt bepaalt de curve, omdat het zowel numeriek stabiel is als de plaats met de grootste verhindering langs het scheurfront.
03. Gevoeligheid voor de bron van de R-curve
Doel
De twee delen van de studie rechtstreeks verbinden: de raakconstructie uit stap 4 herhalen met de experimentele R-curve uit één proefstuk uit deel 1 (C₁ = 29,49, C₂ = 0,642) in plaats van de referentiecurve die representatief is voor meerdere proefstukken (C₁ = 54,45, C₂ = 0,412).
De verschuiving van 14% tussen twee afzonderlijk verdedigbare weerstandscurves is de praktische onderbouwing van de aanbeveling om een werkelijke beslissing over geschiktheid voor gebruik op een curve uit meerdere proefstukken te baseren.
