Badanie wzrostu pęknięć zmęczeniowych i prawo Parisa dla stali

Co mierzy badanie zmęczeniowego wzrostu pęknięcia

Badanie zmęczeniowego wzrostu pęknięcia mierzy, jak istniejące pęknięcie wydłuża się pod wpływem powtarzających się obciążeń. Próbka jest najpierw poddawana wstępnemu pękaniu zmęczeniowemu, tak aby badanie rozpoczynało się od kontrolowanego, ostrego pęknięcia, a nie od nieznanej wady powstałej podczas obróbki. W trakcie dalszego obciążania długość pęknięcia mierzy się przy wybranych liczbach cykli, a wynik podaje się jako przyrost długości pęknięcia przypadający na cykl:

$$\frac{\mathrm{da}}{\mathrm{dN}}$$

$a$ oznacza długość pęknięcia, a $N$ — liczbę cykli obciążenia. Mierzona wielkość nie jest zatem wytrzymałością na rozciąganie, granicą plastyczności ani twardością. Jest to szybkość akumulacji uszkodzeń dla konkretnego pęknięcia, historii obciążenia, geometrii próbki, środowiska oraz procedury pomiarowej.

ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego — od zakresu bliskiego progowi aż do niestabilności sterowanej wartością Kmax. Strong evidence

[1][1] ISO 12108:2018: Metallic materials — Fatigue testing — Fatigue crack growth method. International Organization for Standardization. ISO standard, 2018.

W standardowym badaniu stosuje się zwykle przede wszystkim obciążenia otwierające w trybie I do próbki zwartego rozciągania, próbki rozciąganej z pęknięciem centralnym lub próbki o podobnej geometrii. Siła obciążająca może mieć stałą amplitudę, przy zdefiniowanym współczynniku obciążenia $R=K_{\min }/K_{\max }$, podczas gdy pęknięcie rozwija się w szeregu znanych warunków intensywności naprężeń. Norma ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego „od zakresu bliskiego progowi do niestateczności sterowanej przez $K_{\max }$” (ASTM International, 2023). Norma ISO 12108:2018 podobnie określa badania od progowego zakresu współczynnika intensywności naprężeń, $\Delta K_{\mathrm{th}}$, do początku szybkiego niestatecznego pękania, głównie dla izotropowych materiałów metalicznych poddawanych przeważnie liniowo-sprężystemu obciążeniu w trybie I przy stałym współczynniku obciążenia $R$.

Zakres ten ma znaczenie. Wynikiem jest zależność między szybkością wzrostu pęknięcia a siłą napędową opisaną mechaniką pękania, a nie jedna stała liczba opisująca każdy proces zniszczenia konstrukcyjnego w stali.

Szybkość wzrostu pęknięcia a trwałość do jego inicjacji

Inicjacja pęknięcia i jego propagacja to różne pomiary. Badanie trwałości do inicjacji może rozpoczynać się od gładkiej, polerowanej lub karbowanej próbki i rejestrować liczbę cykli potrzebnych do powstania pęknięcia zgodnie z przyjętym kryterium wykrywania. Kryterium to może obejmować pęknięcie powierzchniowe o określonej długości, zmianę podatności, sygnał związany ze spadkiem potencjału lub inną reakcję przyrządu pomiarowego. Wyniki silnie zależą od stanu powierzchni, wtrąceń, naprężeń resztkowych, kształtu karbu, mikrostruktury oraz czułości metody kontroli.

Badanie zmęczeniowego wzrostu pęknięcia ogranicza znaczną część tej niepewności, ponieważ rozpoczyna się od celowo utworzonego pęknięcia. Badanie odpowiada na węższe pytanie: gdy istnieje pęknięcie o znanym rozmiarze, jak daleko rozwinie się ono podczas każdego kolejnego cyklu? Procedura zwykle wyróżnia trzy obszary fizyczne.

Progowy zakres współczynnika intensywności naprężeńProgowy zakres współczynnika intensywności naprężeń Zdefiniowana proceduralnie granica ΔK, poniżej której w określonych warunkach badania nie wykrywa się trwałego przyrostu pęknięcia.

Przy małej sile napędowej pęknięcie może pozostać nieaktywne albo rozwijać się jedynie okresowo. Jest to obszar bliski progowi, w którym zmierzona szybkość zbliża się do wartości progowej związanej z wybraną procedurą. Powyżej tego obszaru zachodzi stabilna propagacja podkrytyczna. Pęknięcie rozwija się przy każdej grupie cykli, chociaż szybkość może się zmieniać wraz ze wzrostem intensywności naprężeń. Przy dostatecznie wysokim obciążeniu pozostałe przewężenie nie jest już w stanie utrzymać stabilnego wzrostu i rozpoczyna się szybkie pękanie. To końcowe zdarzenie nie jest kolejnym zwykłym punktem na długotrwałej krzywej wzrostu zmęczeniowego; jest to niestateczność, na którą silnie wpływają $K_{\max }$, skrępowanie próbki i odporność na pękanie.

To rozróżnienie pozwala uniknąć częstego błędu. Stal może charakteryzować się długą trwałością do inicjacji, ale stosunkowo dużą szybkością propagacji po pojawieniu się pęknięcia, albo może wykazywać wczesne powstawanie pęknięć, po którym następuje powolny stabilny wzrost. Sama krzywa szybkości wzrostu pęknięcia nie pozwala przewidzieć liczby cykli potrzebnych do zainicjowania pęknięcia w elemencie bez pęknięć.[2][2] Fatigue crack growth testing. TWI. TWI technical knowledge, 2024.

Ten sam problem pojawia się przy porównywaniu norm. TWI poinformował w 2024 roku, że ASTM E647 i ISO 12108 stosują różne kryteria dotyczące progu oraz zmniejszania $\Delta K$. W rezultacie dwa laboratoria mogą badać nominalnie podobną stal i uzyskać różne wyniki w obszarze bliskim progowi, przy czym żaden z tych zbiorów danych nie jest automatycznie nieważny. Na zgłaszaną krzywą wpływają geometria próbki, kalibracja długości pęknięcia, sposób sterowania obciążeniem, środowisko, współczynnik obciążenia oraz zasady redukcji danych.

Rola $\mathrm{da}/\mathrm{dN}$, $\Delta K$ i $K_{\max }$

Podstawowym parametrem siły napędowej w liniowo-sprężystym badaniu wzrostu pęknięcia zmęczeniowego jest zakres współczynnika intensywności naprężeń:

$$\Delta K=K_{\max }-K_{\min }$$

Dla danej próbki $K$ zależy od przyłożonego naprężenia, długości pęknięcia i współczynnika korekcji geometrii. Często stosuje się uproszczoną postać:

$$K=Y\sigma \sqrt{\pi a}$$

gdzie $Y$ oznacza geometrię próbki i pęknięcia, $\sigma$ jest przyłożonym naprężeniem, a $a$ — długością pęknięcia. Zatem to samo naprężenie nominalne nie wywołuje takiej samej siły napędowej pęknięcia w próbkach o różnych rozmiarach pęknięcia lub różnej geometrii. Dlatego wyniku dotyczącego wzrostu pęknięcia nie można przenosić między konfiguracjami bez sprawdzenia właściwego rozwiązania dla intensywności naprężeń oraz granic jego stosowalności.

Zależność Parisa i Erdoğana, przedstawiona w ich publikacji z 1963 roku, opisuje środkową część krzywej, odpowiadającą umiarkowanym szybkościom wzrostu, w następujący sposób:

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

Wykładnik prawa ParisaWykładnik prawa Parisa Dopasowane nachylenie wykresu logarytmu szybkości wzrostu pęknięcia względem logarytmu zakresu współczynnika intensywności naprężeń w wybranym przedziale.

W niektórych sprawozdaniach zamiast $m$ stosuje się $n$. Na wykresie logarytmicznym $\mathrm{da}/\mathrm{dN}$ w funkcji $\Delta K$ $m$, lub $n$, jest nachyleniem, a $C$ — wyrazem wolnym przy przyjętych jednostkach i konwencji dopasowania. Stałe te są zatem warunkowe. Zmiana jednostek zmienia wartość liczbową $C$; zmiana współczynnika obciążenia, środowiska, temperatury, mikrostruktury stali lub przedziału dopasowania może zmienić obie stałe.

Paris i Erdoğan opowiadali się za oceną w szerokich zakresach szybkości wzrostu pęknięcia i z wykorzystaniem wielu próbek, zamiast traktowania wąskiego dopasowanego odcinka jako uniwersalnego prawa materiałowego. Ostrzeżenie to pozostaje istotne, ponieważ równanie Parisa nie opisuje dobrze obszaru progowego i nie uwzględnia końcowego przyspieszenia sterowanego przez $K_{\max }$. W pobliżu progu dominujące znaczenie mogą mieć zamykanie pęknięcia i bariery mikrostrukturalne. W pobliżu niestateczności decydujące stają się maksymalna intensywność naprężeń i odporność na pękanie.

Standardowe analizy całkowania uszkodzeń zasadniczo wykorzystują dane z badań próbek z pęknięciami przy stałej amplitudzie. Strong evidence

[3][3] Damage Tolerance Design Handbook. AFGROW. AFGROW Damage Tolerance Design Handbook, 2024.

AFGROW stwierdza, że dane wykorzystywane w standardowych analizach całkowania uszkodzeń pochodzą na ogół z badań przy stałej amplitudzie na próbkach zawierających pęknięcia, przy czym większość takich badań obejmuje norma ASTM E647. Dane te można następnie całkować wzdłuż historii intensywności naprężeń w konstrukcji, aby oszacować trwałość do wzrostu pęknięcia. To późniejsze obliczenie jest założeniem przyjętym w ocenie, a nie częścią podstawowego pomiaru.

W pracy NIST z 2001 roku dotyczącej dwóch ferrytyczno-perlitycznych stali rurociągowych wykorzystano zakrzywione próbki M(T) i analizowano rozrzut szybkości, przypisując zmierzoną zmienność stałej $C$ prawa Parisa. W sprawozdaniu NIST z 1976 roku dotyczącym spawanych elementów ze stali zawierającej 5% Ni również dopasowano dane dotyczące wzrostu pęknięć zmęczeniowych za pomocą potęgowego prawa Parisa oraz określono $C$ jako wyraz wolny, a $n$ jako nachylenie w układzie logarytmiczno-logarytmicznym. Żaden z tych przykładów nie przekształca $C$ w niezależną od gatunku stałą stali.

Przepisy dotyczące konstrukcji odzwierciedlają to ograniczenie. FAA AC 25.571-1D dotyczy analizy wzrostu pęknięć, prawdopodobnych lokalizacji uszkodzeń zmęczeniowych oraz planowania kontroli konstrukcji samolotów transportowych, natomiast 14 CFR 25.571 wymaga, w stosownych przypadkach, przeprowadzania ocen tolerancji uszkodzeń popartych wynikami badań przy obciążeniach powtarzalnych. Badanie dostarcza zmierzonych danych dotyczących zachowania pęknięcia podczas wzrostu; geometria, widmo obciążenia, założenia dotyczące początkowych wad oraz niezawodność kontroli decydują o znaczeniu tego zachowania w konstrukcji.

Prawo Parisa: równanie, parametry i ograniczenia

Zależność Parisa–Erdoğana

Prawo Parisa zapisuje się zwykle jako

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

gdzie $a$ oznacza długość pęknięcia, $N$ — liczbę cykli obciążenia, a $\mathrm{da}/\mathrm{dN}$ — prędkość wzrostu pęknięcia zmęczeniowego. Wyrażenie $\Delta K=K_{\max }-K_{\min }$ oznacza zakres współczynnika intensywności naprężeń, zwykle wyrażany w MPa√m lub ksi√in. $C$ i $m$ są parametrami dopasowania, a nie uniwersalnymi właściwościami, które można odczytać z oznaczenia gatunku stali.

P. C. Paris i F. Erdoğan wprowadzili tę zależność empiryczną w artykule z 1963 roku pt. „A Critical Analysis of Crack Propagation Laws”. Ich praca połączyła zmierzone prędkości przyrostu pęknięcia z zakresem sprężystego współczynnika intensywności naprężeń i wykazała, że prawo potęgowe może opisywać znaczną część stabilnego wzrostu pęknięcia zmęczeniowego. Równanie opisuje pewien zakres zachowania, a nie każdy etap propagacji pęknięcia zmęczeniowego. W pobliżu progu zamykanie pęknięcia, mikrostruktura i wpływ historii obciążenia mogą powodować spadek prędkości poniżej prostej Parisa. Gdy $K_{\max }$ zbliża się do granicy kontrolowanej przez odporność materiału na pękanie, przyspieszenie prowadzące do niestabilnego pękania powoduje kolejne odstępstwo od tej zależności.

Ograniczenie zakresu ma istotne znaczenie. Norma ASTM E647-23 obejmuje wyznaczanie prędkości wzrostu pęknięcia zmęczeniowego „od zakresu bliskiego progowi do niestabilności kontrolowanej przez Kmax”, natomiast norma ISO 12108:2018 określa badania od progowego zakresu współczynnika intensywności naprężeń $\Delta K_{\mathrm{th}}$ do początku szybkiego niestabilnego pękania. Żadna z tych norm nie sprawia, że każdy zmierzony punkt w tym przedziale można uwzględnić w jednym dopasowaniu prawa Parisa. Metoda badania określa sposób wyznaczania prędkości, długości pęknięcia, siły, podatności, intensywności naprężeń i kryteriów ważności; analityk nadal musi zidentyfikować przedział, w którym prawo potęgowe stanowi akceptowalny opis.

Oryginalne rozważania Parisa–Erdoğana również przemawiały za oceną zachowania związanego ze wzrostem pęknięcia w szerokim zakresie prędkości oraz za wykorzystaniem wielu próbek, zamiast przypisywania dużego znaczenia wąskiemu fragmentowi dopasowania. Jest to cenna przestroga przy badaniu stali. Prosta, która wygląda na liniową w ciągu jednej dekady $\Delta K$, może uzasadniać lokalne przybliżenie inżynierskie, ale nie może potwierdzać zachowania w pobliżu $\Delta K_{\mathrm{th}}$, w pobliżu niestabilności pękania ani przy innej historii obciążenia.

Znaczenie $C$ oraz $m$ lub $n$

Zastosowanie logarytmów dziesiętnych daje

$$\log \left(\frac{\mathrm{da}}{\mathrm{dN}}\right)=\log C+m\log \left(\Delta K\right).$$

Na wykresie $\mathrm{da}/\mathrm{dN}$ w funkcji $\Delta K$ w skali logarytmicznej $m$ jest nachyleniem, a $\log C$ — wyrazem wolnym na osi pionowej dla wybranej jednostki $\Delta K$. W niektórych sprawozdaniach zamiast $m$ stosuje się $n$, dlatego

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^n$$

ma takie samo znaczenie, gdy $n=m$. Zmienia się oznaczenie, ale nie zmienia się dopasowana zależność.

Wykładnik opisuje, jak gwałtownie zmierzona prędkość wzrostu pęknięcia zmienia się wraz z zakresem współczynnika intensywności naprężeń w wybranym przedziale. Większy wykładnik oznacza bardziej strome nachylenie prostej. Nie oznacza to, że dana stal jest po prostu „bardziej odporna na zmęczenie” w każdych warunkach. Wyraz wolny określa poziom prędkości, ale $C$ ma wymiar. Jeżeli $\mathrm{da}/\mathrm{dN}$ wyraża się w m/cykl, a $\Delta K$ w MPa√m, wówczas $C$ ma jednostkę

\[ \mathrm{m/cykl}\,(\mathrm{MPa}\sqrt{\mathrm m})^{-m}. \]

Zmiana jednostki $\Delta K$ z MPa√m na ksi√in zmienia wartość liczbową $C$, mimo że krzywa fizyczna pozostaje niezmieniona. Podana wartość $C$ bez jednostek prędkości wzrostu pęknięcia, jednostek współczynnika intensywności naprężeń, wykładnika oraz przedziału dopasowania jest informacją niepełną.

Wyraz wolny nie musi również odpowiadać bezpośrednio zaobserwowanemu punktowi. Większość dopasowanych danych nie obejmuje $\Delta K=1$ w wybranym układzie jednostek, a wyraz wolny uzyskuje się przez ekstrapolację w przestrzeni logarytmicznej. Sprawia to, że $C$ jest wrażliwe na dopasowane nachylenie oraz na zakres uwzględnionych danych. W sprawozdaniu NIST dotyczącym spawanych elementów ze stali zawierającej 5% Ni wskazano na $C$ jako wyraz wolny, a na $n$ jako nachylenie dopasowania logarytmicznego $\mathrm{da}/\mathrm{dN}$ w funkcji $\Delta K$. Opis ten jest matematycznie precyzyjny, ale nie należy go mylić ze stwierdzeniem, że $C$ jest stałą metalurgiczną.

Zmienność zmierzonych szybkości wzrostu pęknięć w stali rurociągowej można przedstawić jako zmienność współczynnika C prawa Parisa dla badanego zbioru danych. Limited evidence

W odrębnym badaniu NIST dotyczącym dwóch rurociągowych stali ferrytyczno-perlitycznych zastosowano zakrzywione próbki M(T) i przeanalizowano „rozrzut zmierzonych prędkości”, przypisując go stałej prawa Parisa $C$. Jest to model niepewności dotyczący tych pomiarów i warunków. Nie eliminuje on rozrzutu powodowanego przez ustawienie próbki, kształt frontu pęknięcia, naprężenia własne, pasmowość mikrostruktury, decyzje operatora ani statystyczny rozkład przyrostów pęknięcia.

Dlaczego dopasowanie logarytmiczne nie jest uniwersalną stałą materiałową

Dopasowanie prawa Parisa odnosi się do konfiguracji badania i określonej populacji danych. Szczególnie istotny jest współczynnik obciążenia $R=K_{\min }/K_{\max }$, ponieważ zamykanie pęknięcia i maksymalny poziom intensywności naprężeń mogą zmieniać zmierzoną prędkość przy tej samej wartości $\Delta K$. Znaczenie ma również środowisko: powietrze, próżnia, wilgotność, roztwory wodne i warunki zawierające wodór mogą prowadzić do różnych prędkości propagacji w tej samej stali. Temperatury, częstotliwość, produkty korozji i czasy utrzymania obciążenia mogą dodatkowo zmieniać mechanizm.

Geometria wpływa na sposób obliczania $K$ oraz na dokładność odwzorowania siły napędowej pękania. Próbki zwartościowe, próbki rozciągane z pęknięciem centralnym oraz zakrzywione próbki M(T) mogą powodować różne warunki więzów, frontu pęknięcia i naprężeń własnych. Prawidłowe rozwiązanie dla współczynnika intensywności naprężeń ogranicza jedno ze źródeł błędu, ale nie sprawia, że próbki stają się mechanicznie identyczne. Należy również określić stan stali: przykłady obejmują rurociągowe stale ferrytyczno-perlityczne, spawane elementy ze stali zawierającej 5% Ni, strefy wpływu ciepła złączy spawanych, stale ulepszane cieplnie oraz odpowiednią specyfikację wyrobu lub spawania. Samo określenie „stal” nie jest wystarczającym opisem materiału.

ASTM E647 i ISO 12108 stosują różne kryteria progowe i kryteria zmniejszania ΔK. Strong evidence

Przedział dopasowania należy podawać wraz z $C$ i $m$ lub $n$: należy podać zakres $\Delta K$, zakres $\mathrm{da}/\mathrm{dN}$, współczynnik obciążenia, środowisko, temperaturę, geometrię próbki, orientację próbki, metodę sterowania siłą oraz procedurę redukcji danych. Normy ASTM E647-23 i ISO 12108:2018 nie ustanawiają identycznych kryteriów progowych ani kryteriów dotyczących malejącego $\Delta K$; TWI wyraźnie wskazuje, że ich procedury różnią się w tych obszarach. Wartości progowej uzyskanej zgodnie z jedną normą nie należy zatem porównywać tak, jakby była automatycznie równoważna wartości uzyskanej zgodnie z drugą.

AFGROW stwierdza, że dane stosowane w standardowych analizach całkowania uszkodzeń pochodzą na ogół z badań przy stałej amplitudzie na próbkach z pęknięciem, najczęściej objętych normą ASTM E647. Ocena konstrukcji polega następnie na całkowaniu wybranej zależności wzrostu pęknięcia w całym widmie obciążenia, często z uwzględnieniem poprawek geometrycznych i założeń dotyczących inspekcji. Jest to odrębny etap od pomiaru $\mathrm{da}/\mathrm{dN}$. FAA AC 25.571-1D i 14 CFR 25.571 umieszczają analizę wzrostu pęknięcia, prawdopodobne lokalizacje uszkodzeń zmęczeniowych, wspartą badaniami ocenę powtarzalnego obciążenia oraz planowanie inspekcji w ramach zgodności z wymaganiami dotyczącymi tolerowania uszkodzeń w statkach powietrznych transportowych.

Praktyczna zasada jest prosta: nie należy przenosić wartości $C$ ani $m$ między różnymi współczynnikami obciążenia, środowiskami, geometriami próbek ani stanami stali bez odpowiednich dowodów. Należy ponownie dopasować właściwe dane lub uzasadnić przeniesienie za pomocą zweryfikowanych korekt i granic niepewności. Równanie Parisa jest wartościowe właśnie dlatego, że stanowi kontrolowane przybliżenie — a nie dlatego, że jest kompletnym opisem zmęczenia stali.

Intensywność naprężeń, ΔK, Kmax i współczynnik obciążenia R

Badania wzrostu pęknięć zmęczeniowych nie korelują przyrostu pęknięcia wyłącznie z nominalnym naprężeniem. Korelują zmierzoną szybkość, zazwyczaj $\mathrm{da}/\mathrm{dN}$, z elastyczną intensywnością naprężeń na wierzchołku pęknięcia. W przypadku obciążenia dominująco w trybie I istotnymi wielkościami są maksymalny współczynnik intensywności naprężeń $K_{\max }$, wartość minimalna $K_{\min }$, ich zakres

$$\Delta K=K_{\max }-K_{\min },$$

oraz współczynnik obciążenia

$$R=\frac{K_{\min }}{K_{\max }}.$$

Dla ustalonej próbki i długości pęknięcia $R$ jest również często zapisywany jako $R=F_{\min }/F_{\max }$, ponieważ obie wartości intensywności naprężeń są proporcjonalne do przyłożonej siły. Definicje te są ze sobą powiązane, ale nie opisują tej samej cechy cyklu. $\Delta K$ mierzy amplitudę elastycznej siły napędzającej pęknięcie na jego wierzchołku; $R$ opisuje położenie tego cyklu, w tym jego rozciągającą składową średnią. Dwa badania mogą mieć taką samą wartość $\Delta K$ i różne wartości $R$, albo taki sam nominalny zakres naprężeń i różne wartości $\Delta K$.

Obciążenie w trybie I w zakresie sprężysto-liniowym

Tryb I jest trybem otwierania: powierzchnie pęknięcia oddzielają się normalnie do nich w miarę zwiększania obciążenia rozciągającego. ASTM E647-23 i ISO 12108:2018 dotyczą przede wszystkim pomiarów wzrostu pęknięć zmęczeniowych w warunkach dominującego obciążenia w trybie I, opisywanego mechaniką pękania w zakresie sprężysto-liniowym, zazwyczaj przy stałym współczynniku obciążenia $R$. „Sprężysto-liniowy” nie oznacza, że stal nie zawiera obszarów plastycznych. Oznacza, że odpowiedź próbki i pole na wierzchołku pęknięcia można w wystarczającym stopniu opisać za pomocą mechaniki pękania w zakresie sprężysto-liniowym, a ewentualna niewielka strefa plastyczna jest traktowana jako lokalna poprawka, a nie jako dominująca skala konstrukcyjna.

Dla idealizowanego pęknięcia w elemencie poddanym naprężeniu rozciągającemu współczynnik intensywności naprężeń w trybie I często wyraża się jako

$$K_I=Y(a/W,\dots )\sigma \sqrt{\pi a},$$

gdzie $\sigma$ jest przyłożonym nominalnym naprężeniem, $a$ jest długością pęknięcia, $Y$ jest bezwymiarowym współczynnikiem geometrycznym, a $W$ jest charakterystycznym wymiarem próbki. Współczynnik ten może zależeć od długości pęknięcia, szerokości i grubości próbki, układu obciążenia oraz odległości od swobodnych krawędzi lub uchwytów. W badaniu sterowanym siłą równoważne wyrażenie wykorzystuje przyłożoną siłę zamiast nominalnego naprężenia.

Wartości cykliczne wynikają bezpośrednio z zależności:

\[ K_{\max}=Y\sigma_{\max}\sqrt{\pi a}, \qquad K_{\min}=Y\sigma_{\min}\sqrt{\pi a}. \]

W konsekwencji,

$$\Delta K=Y\Delta \sigma \sqrt{\pi a}$$

jeżeli ten sam współczynnik geometryczny obowiązuje w całym cyklu. W miarę wzrostu pęknięcia $a$ zwiększa się, a $\Delta K$ zazwyczaj rośnie, nawet jeśli zakres przyłożonej siły pozostaje stały. Ten rosnący czynnik napędzający wyjaśnia, dlaczego badanie przy stałej amplitudzie nie daje stałej szybkości wzrostu pęknięcia na całej jego długości.

Podejście oparte na intensywności naprężeń ma swoje ograniczenia. Jeżeli strefa plastyczna jest zbyt duża w stosunku do wymiarów pozostałego przekroju, jeżeli więzy na froncie pęknięcia ulegają znaczącym zmianom lub jeżeli sprężyste odciążenie nie jest rozsądnym przybliżeniem, nominalna wartość $\Delta K$ może niewystarczająco charakteryzować cykl na wierzchołku pęknięcia. ASTM E647-23 obejmuje pomiary od zakresu bliskiego progowi aż do niestabilności sterowanej przez $K_{\max }$, jednak obowiązujący zakres nadal zależy od wielkości próbki, wymagań dotyczących pozostałego przekroju, jakości pęknięcia oraz procedury badania.

Wpływ geometrii na współczynnik intensywności naprężeń

Geometria nie jest poprawką dodawaną po przeprowadzeniu badania; stanowi część mierzonej wielkości z zakresu mechaniki pękania. Próbka kompaktowa, próbka z pęknięciem środkowym poddana rozciąganiu oraz próbka z jednostronnym karbem zginana mogą doświadczać tego samego nominalnego zakresu naprężeń, a mimo to dawać różne wartości $K$, ponieważ różnią się funkcjami $Y$. Długość pęknięcia również występuje bezpośrednio w wyrażeniu \(\sqrt a\). Dlatego dłuższe pęknięcie doświadcza większej intensywności naprężeń przy tym samym nominalnym naprężeniu.

Rozróżnienie to ma znaczenie przy porównywaniu danych dotyczących stali. Stal rurociągowa o strukturze ferrytyczno-perlitycznej badana w zakrzywionej próbce M(T) nie jest poddawana działaniu takich samych więzów na wierzchołku pęknięcia jak stal badana w próbce kompaktowej, nawet jeśli w obu sprawozdaniach podano ten sam nominalny zakres siły. W badaniu NIST z 2001 r. dotyczącym dwóch stali rurociągowych o strukturze ferrytyczno-perlitycznej wykorzystano zakrzywione próbki M(T), a rozrzut zmierzonych szybkości potraktowano poprzez przypisanie go stałej $C$ w prawie Parisa. Wynik ten nie oznacza, że $C$ jest właściwością niezależną od geometrii; pokazuje natomiast, w jaki sposób niepewność dopasowanego prawa szybkości może uwzględniać zmienność doświadczalną i materiałową w ramach określonego układu badawczego.

Ta sama ostrożność dotyczy złączy spawanych. W sprawozdaniu NIST z 1976 r. dotyczącym złączy spawanych ze stali zawierającej 5% Ni dopasowano dane za pomocą postaci prawa Parisa, identyfikując $C$ jako wyraz wolny, a $n$ jako nachylenie wykresu logarytmicznego $\mathrm{da}/\mathrm{dN}$ w funkcji $\Delta K$. Dopasowany wykładnik i wyraz wolny opisują ten zbiór danych oraz warunki badania. Nie należy przenosić ich do innej geometrii spoiny, innego stanu naprężeń własnych, innej temperatury, środowiska lub współczynnika obciążenia bez odpowiednich dowodów.

Paris i Erdoğan wprowadzili powszechnie zapisywaną zależność

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

w swojej pracy z 1963 r. Na osiach logarytmicznych $m$ — nazywane w niektórych sprawozdaniach $n$ — jest nachyleniem, a $C$ jest wyrazem wolnym. Równanie to jest użytecznym opisem środkowego zakresu, ale nie stanowi kompletnego prawa zmęczeniowego. Paris i Erdoğan opowiadali się za oceną obejmującą szerokie zakresy szybkości wzrostu pęknięcia i wiele próbek, zamiast przedstawiania wąskiego odcinka dopasowania jako uniwersalnej stałej materiałowej.

Współczynnik R i zjawiska zamykania pęknięcia

$\Delta K$ i $R$ należy podawać łącznie, gdy wymagają tego warunki badania. Na przykład przy $R=0.1$ minimalna rozciągająca intensywność naprężeń stanowi 10% wartości maksymalnej. Przy ujemnym $R$ część cyklu może mieć charakter ściskający. Badanie przy $R=0.5$ może mieć taką samą wartość $\Delta K$ jak badanie przy $R=0.1$, a jednocześnie wyższą wartość $K_{\max }$ i większe średnie obciążenie rozciągające. Ich szybkości wzrostu pęknięcia nie muszą być takie same.

Jedną z przyczyn jest zamykanie pęknięcia. Chropowate powierzchnie przełomu, pozostałości tlenków, odkształcenie plastyczne i naprężenia własne mogą utrzymywać powierzchnie pęknięcia w kontakcie podczas części nominalnie rozciągającego cyklu. Pęknięcie doświadcza wówczas mniejszego efektywnego zakresu otwierania, niż wynikałoby to z obliczonej zewnętrznie wartości $\Delta K$. Zwiększanie $R$ na ogół skraca część cyklu, w której pęknięcie pozostaje zamknięte, chociaż odpowiedź zależy od mikrostruktury stali, ścieżki pęknięcia, środowiska, więzów w próbce i historii obciążenia. Znaczenie ma również $K_{\max }$: zachowania w pobliżu progu, uszkodzenie materiału na wierzchołku pęknięcia oraz zbliżanie się do niestabilnego pękania nie mogą być opisane wyłącznie za pomocą $\Delta K$.

Dlatego identyczne nominalne zakresy naprężeń mogą prowadzić do różnych szybkości wzrostu pęknięcia. Geometria może zmieniać $Y$; długość pęknięcia może zmieniać $a$; inna wartość $R$ może wpływać na zamykanie pęknięcia i $K_{\max }$; a naprężenia własne mogą przesuwać lokalny cykl. ASTM E647-23 i ISO 12108:2018 stosują również różne kryteria wyznaczania progu i zmniejszania $\Delta K$, jak zauważa TWI, dlatego podana wartość $\Delta K_{\mathrm{th}}$ nie jest automatycznie porównywalna między normami. W dokumentacji badania należy określić geometrię próbki, zakres długości pęknięcia, współczynnik obciążenia, środowisko, procedurę wyznaczania progu oraz ograniczenia dotyczące ważności wyników.

Szczegóły te mają znaczenie, gdy dane są przenoszone z krzywej laboratoryjnej do oceny tolerancji uszkodzeń. AFGROW stwierdza, że standardowe dane do całkowania uszkodzeń pochodzą na ogół z badań przy stałej amplitudzie prowadzonych na próbkach z pęknięciami. Oceny statków powietrznych zgodnie z FAA AC 25.571-1D i 14 CFR 25.571 łączą następnie prawa wzrostu pęknięcia z dowodami dotyczącymi obciążeń powtarzalnych, prawdopodobnych miejsc uszkodzeń zmęczeniowych, krytycznych długości pęknięć oraz progów inspekcji. Stałe Parisa są danymi wejściowymi w tym łańcuchu, a nie jego zamiennikiem.

ASTM E647: Zakres, próbki i redukcja danych

ASTM E647-23 jest głównym amerykańskim standardem dotyczącym pomiaru prędkości wzrostu pęknięć zmęczeniowych w materiałach metalicznych. Zakres tego standardu obejmuje obszar od reżimu bliskiego progowi aż po niestabilność kontrolowaną przez Kmax, a wyniki są wyrażane za pomocą zakresu współczynnika intensywności naprężeń w ujęciu liniowo-sprężystym, ΔK. Zakres ten ma znaczenie, ponieważ krzywa wzrostu pęknięcia zmęczeniowego nie jest pojedynczą właściwością odczytywaną bezpośrednio ze świadectwa dotyczącego stali. Jest wynikiem uzyskanym przy użyciu określonej próbki, historii obciążenia, metody pomiaru pęknięcia, środowiska oraz procedury redukcji danych.

Badanie zazwyczaj prowadzi do wyznaczenia zależności między prędkością wzrostu pęknięcia, $\mathrm{da}/\mathrm{dN}$, a przyłożonym zakresem współczynnika intensywności naprężeń:

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

P. C. Paris i F. Erdoğan przedstawili tę zależność empiryczną w swojej publikacji z 1963 roku. Na wykresie logarytmicznym $\mathrm{da}/\mathrm{dN}$ w funkcji ΔK $m$ jest nachyleniem, a $C$ wyrazem wolnym, gdy równanie zapisuje się w standardowej postaci. W niektórych sprawozdaniach zamiast $m$ jako wykładnika stosuje się $n$. Symbole nie stanowią głównego problemu. Najważniejsze są zakres badania oraz dane wykorzystane do ich wyznaczenia.

Paris i Erdoğan opowiadali się za oceną wzrostu pęknięcia w szerokim zakresie prędkości oraz z wykorzystaniem wielu próbek, zamiast przedstawiania wąskiego dopasowanego odcinka jako prawa uniwersalnego. Ta ostrożność pozostaje istotna w przypadku stali. Wartość $m$ dopasowana w centralnym, w przybliżeniu liniowym obszarze krzywej logarytmiczno-logarytmicznej nie opisuje zachowania w pobliżu progu, szybkiego wzrostu w pobliżu pęknięcia, opóźnienia wzrostu po przeciążeniu ani każdego możliwego stosunku naprężeń.

Od reżimu bliskiego progowi do niestabilności kontrolowanej przez Kmax

ASTM International stwierdza, że ASTM E647 obejmuje wyznaczanie prędkości wzrostu pęknięć zmęczeniowych „od reżimu bliskiego progowi aż po niestabilność kontrolowaną przez Kmax” (2023). W obszarze małych prędkości wzrostu podstawową wielkością jest progowy zakres współczynnika intensywności naprężeń, ΔKth: zakres, poniżej którego pęknięcie w określonych warunkach badania albo nie propaguje, albo rozwija się zbyt wolno dla zdefiniowanej metody pomiarowej. W pobliżu progu niewielkie zmiany w zamykaniu pęknięcia, stanie powierzchni, naprężeniach szczątkowych, środowisku lub historii obciążenia mogą powodować znaczne zmiany pozornej prędkości wzrostu.

Na drugim końcu zakresu decydujące znaczenie ma Kmax. Kmax jest maksymalnym współczynnikiem intensywności naprężeń w cyklu, a jego wzrost może przesunąć próbkę w kierunku pęknięcia statycznego lub szybkiego, niestabilnego wydłużenia pęknięcia, nawet gdy samo ΔK nie opisuje w pełni danego stanu. Obszar środkowy jest tym, w którym najczęściej stosuje się zależność Parisa, jednak tej wygody nie należy utożsamiać z pełnym pokryciem krzywej. Wzrost może odbiegać od linii prostej zarówno przy małych, jak i dużych wartościach ΔK.

ISO 12108:2018 opisuje powiązany zakres badania — od ΔKth do początku szybkiego, niestabilnego pęknięcia — przede wszystkim dla izotropowych materiałów metalicznych poddanych dominującemu liniowo-sprężystemu obciążeniu w trybie I przy stałym stosunku sił $R$, gdzie $R=K_{\min }/K_{\max }$. ASTM E647 i ISO 12108 nie są zamiennymi oznaczeniami jednej identycznej procedury. TWI poinformował w 2024 roku, że standardy te stosują różne kryteria dotyczące progu oraz zmniejszającego się ΔK. Wynik oznaczony jako „progowy” wymaga zatem wskazania standardu, sekwencji obciążenia oraz podstawy akceptacji.

Stosunek naprężeń ma szczególne znaczenie, ponieważ zamykanie pęknięcia i efektywna część cyklu mogą zmieniać się wraz z $R$. Dwa badania tej samej stali, przy takim samym nominalnym ΔK, lecz przy różnych stosunkach naprężeń, mogą dać różne prędkości wzrostu. Znaczenie może mieć również częstotliwość, między innymi wskutek nagrzewania, pękania wspomaganego przez środowisko oraz procesów zależnych od czasu. Z tego powodu użyteczny wynik badania ASTM E647 powinien obejmować co najmniej typ próbki, wszystkie główne wymiary, grubość, orientację materiału, częstotliwość obciążenia, środowisko, stosunek naprężeń oraz zastosowaną historię obciążenia. Należy także określić stan stali: postać wyrobu, obróbkę cieplną, spoinę lub materiał rodzimy, a także istotny stan mikrostrukturalny, ponieważ czynniki te mogą wpływać zarówno na zmierzoną krzywą, jak i na jej rozrzut.

Wynik badania jest pomiarem wydłużenia pęknięcia w określonych warunkach. Nie jest automatycznie dopuszczalnym parametrem konstrukcyjnym. Analiza tolerancji uszkodzeń może wymagać modelu progowego, korekty dla obszaru bliskiego progowi, ograniczenia Kmax, uwzględnienia zamykania pęknięcia lub zastosowania innej reguły wzrostu pęknięcia. AFGROW’s 2024 Damage Tolerance Design Handbook stwierdza, że dane do standardowych analiz całkowania uszkodzeń pochodzą na ogół z badań przy stałej amplitudzie na próbkach z pęknięciami, z których większość jest objęta ASTM E647. W analizach konstrukcji następnie całkuje się wybrane prawo wzrostu w warunkach zmiennej historii współczynnika intensywności naprężeń; to późniejsze obliczenie jest odrębne od samego badania.

Typowe konfiguracje próbek z pęknięciem

Badania według ASTM E647 wykorzystują próbki, w których celowo wprowadzone pęknięcie rozwija się pod wpływem kontrolowanego obciążenia cyklicznego. Próbka kompaktowa rozciągana, powszechnie oznaczana jako C(T), jest szeroko stosowana, ponieważ wymaga stosunkowo niewielkiej ilości materiału i zapewnia silne, mierzalne pole napędzające pęknięcie. Jej geometria z obciążeniem przykładanym przez sworznie sprawia również, że istotne stają się ustawienie oraz efekty związane z linią obciążenia. Wymiary i grubość nie są drugorzędnymi szczegółami sprawozdawczymi: wpływają na skrępowanie, plastyczność, kształt frontu pęknięcia oraz zasadność interpretacji liniowo-sprężystej.

Próbki rozciągane ze środkowym pęknięciem, oznaczane jako M(T), są przydatne, gdy pożądane jest szerokie przewężenie oraz centralne pęknięcie, szczególnie w badaniach blach, płyt i spoin. W pracy NIST z 2001 roku dotyczącej dwóch ferrytyczno-perlitycznych stali rurociągowych wykorzystano zakrzywione próbki M(T) i przeanalizowano rozrzut zmierzonych prędkości, przypisując go zmienności stałej $C$ w prawie Parisa. Badanie to pokazuje, dlaczego dopasowana stała może odzwierciedlać zmienność eksperymentalną i materiałową w takim samym stopniu jak ustaloną właściwość stali.

W badaniach wzrostu pęknięć zmęczeniowych spotyka się również układy z próbkami rozciąganymi i zginanymi z jednostronnym karbem, często zapisywane jako SEN(T) i SEN(B), jednak dokładne oznaczenie i konfiguracja muszą być podane, a nie wywnioskowane ze zdjęcia. W przypadku złączy spawanych może być konieczna orientacja próbki, która umieszcza pęknięcie względem metalu spoiny, linii wtopienia lub strefy wpływu ciepła. W przypadku płyt i blach znaczenie ma kierunek walcowania. Znaczenie mają również orientacje poprzeczne i wzdłużne, zwłaszcza gdy wtrącenia, pasmowość lub anizotropowa odporność na pękanie wpływają na rozwój pęknięcia.

Grubość próbki wpływa na skrępowanie w kierunku grubości oraz na stopień, w jakim front pęknięcia pozostaje odpowiedni dla przybliżenia stanu odkształcenia płaskiego. Cienki materiał może wykazywać większą plastyczność i większe zróżnicowanie frontu pęknięcia; gruba próbka może zapewniać większe skrępowanie, nie czyniąc jednak wszystkich pozostałych założeń prawidłowymi. Zakrzywione próbki wprowadzają własne poprawki geometryczne. Efekty te uwzględnia się za pomocą właściwego rozwiązania dla współczynnika intensywności naprężeń, a nie przez samo porównywanie nominalnych zakresów siły.

Oznaczenie stali powinno pozostać powiązane z wynikiem. Na przykład „złącze spawane ze stali zawierającej 5% Ni” nie dostarcza informacji równoważnej gatunkowi stali, normie wyrobu ani dokumentacji procedury spawania. W 1976 roku NIST opublikował dopasowanie prawa potęgowego Parisa do danych dotyczących wzrostu pęknięć zmęczeniowych w złączach spawanych ze stali zawierającej 5% Ni, określając $C$ jako wyraz wolny, a $n$ jako nachylenie wykresu logarytmiczno-logarytmicznego. Taki wynik wymaga podania lokalizacji spoiny, orientacji, środowiska oraz stosunku obciążenia, zanim będzie można porównać go z danymi dotyczącymi stali płytowej lub rurowej.

Przyrostowy pomiar pęknięcia i obliczanie prędkości wzrostu

Podstawowym pomiarem jest przyrost długości pęknięcia w czasie określonej liczby zmierzonych cykli. Długość pęknięcia można śledzić w odstępach czasu podczas cyklicznego obciążania za pomocą odpowiedniej metody bezpośredniej lub pośredniej, przy czym należy odnotować metodę i jej kalibrację. Surowy zapis nie jest jeszcze prędkością wzrostu pęknięcia. Zawiera informacje o sile, przemieszczeniu lub podatności, liczbie cykli, oszacowaniach długości pęknięcia oraz możliwym szumie pomiarowym.

Redukcja danych przekształca kolejne obserwacje długości pęknięcia w przyrostową prędkość wzrostu, w ujęciu koncepcyjnym:

$$\frac{\mathrm{da}}{\mathrm{dN}}\approx \frac{\Delta a}{\Delta N}$$

Odpowiadające ΔK oblicza się na podstawie zakresu przyłożonego obciążenia, długości pęknięcia, geometrii próbki oraz jej wymiarów, wykorzystując zależność współczynnika intensywności naprężeń właściwą dla danej konfiguracji. Każdy zredukowany punkt danych łączy zatem zmierzoną długość pęknięcia ze skorygowaną geometrycznie siłą napędzającą pęknięcie. Niedopuszczalne jest przepisanie zakresu siły z jednej próbki i bezpośrednie porównanie go z ΔK z innej próbki.

Następnie analityk analizuje dane $\mathrm{da}/\mathrm{dN}$ w funkcji ΔK, sprawdza występowanie nieważnych lub niestabilnych odcinków i dopasowuje prawo wyłącznie do określonego obszaru. Dopasowanie parysowskie wykonuje się zwykle dla w przybliżeniu liniowego odcinka środkowego na osiach logarytmicznych. Jednostki $C$ zależą od jednostek wybranych dla długości pęknięcia, liczby cykli i współczynnika intensywności naprężeń, dlatego wartości $C$ nie można porównywać bez przeliczenia jednostek. Rozrzut należy przedstawiać lub kwantyfikować, a nie ukrywać za pomocą pojedynczej dopasowanej linii.[4][4] 14 CFR 25.571: Damage-tolerance and fatigue evaluation of structure. Electronic Code of Federal Regulations. 14 CFR, 2024.

To rozróżnienie warunkuje późniejsze decyzje inżynierskie. FAA AC 25.571-1D obejmuje analizę wzrostu pęknięć, wspieraną badaniami ocenę działania pod powtarzalnym obciążeniem oraz planowanie inspekcji w ramach wytycznych dotyczących tolerancji uszkodzeń samolotów transportowych. Zgodnie z 14 CFR 25.571 oceny tolerancji uszkodzeń obejmują prawdopodobne lokalizacje i mechanizmy uszkodzeń zmęczeniowych, analizy obciążeń powtarzalnych poparte dowodami z badań oraz, w stosownych przypadkach, progi inspekcji oparte na wzroście pęknięcia. ASTM E647 zapewnia jedną z podstawowych podstaw doświadczalnych, jednak stałe uzyskują znaczenie inżynierskie dopiero po określeniu próbki, warunków, niepewności oraz zamierzonego zastosowania w analizie tolerancji uszkodzeń.

ISO 12108:2018 i jej relacja do ASTM E647

Badanie wzrostu pęknięcia zmęczeniowego jest procedurą pomiarową, zanim stanie się modelem materiałowym. Mierzoną wielkością jest zazwyczaj $\mathrm{da}/\mathrm{dN}$, czyli przyrost pęknięcia na cykl obciążenia, przedstawiany w funkcji zakresu współczynnika intensywności naprężeń w ujęciu liniowo-sprężystym, $\Delta K$. Znane prawo Parisa,

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m,$$

zostało wprowadzone przez P. C. Parisa i F. Erdoğana w 1963 roku. Na wykresie logarytmicznym $m$ jest nachyleniem, a $C$ wyrazem wolnym, choć w niektórych opracowaniach zamiast $m$ stosuje się $n$. Stałe te opisują wybrany zbiór danych z badań przeprowadzonych w określonych warunkach; nie są automatycznie uniwersalnymi właściwościami materiału.

To rozróżnienie ma znaczenie przy porównywaniu ISO 12108:2018 z ASTM E647-23. Obie normy dotyczą szybkości wzrostu pęknięć zmęczeniowych, ale nie określają identycznych sposobów uzyskania tego samego zbioru danych. Ich definicje, metody wyznaczania progu, procedury zmniejszania $\Delta K$, wymagania dotyczące próbek oraz kryteria akceptacji mogą wpływać na dopasowaną krzywą, zwłaszcza w pobliżu obszaru progowego.

Zakres ISO i określone założenia dotyczące obciążenia

ISO 12108:2018 określa badania służące do wyznaczania szybkości wzrostu pęknięcia zmęczeniowego — od progowego zakresu współczynnika intensywności naprężeń, $\Delta K_{\mathrm{th}}$, do początku szybkiego niestatecznego pękania. Norma odnosi się przede wszystkim do izotropowych materiałów metalicznych poddanych obciążeniu w przeważająco liniowo-sprężystym zakresie, w trybie I, przy stałym współczynniku obciążenia $R$, gdzie

$$R=\frac{K_{\min }}{K_{\max }}.$$

Zakres ten ma istotne znaczenie w badaniach stali. Stal rurociągowa ferrytyczno-perlityczna, stal konstrukcyjna ulepszana cieplnie oraz metal spoiny mogą być określane wspólnym mianem „stali”, jednak ich charakterystyka wzrostu pęknięć może się różnić w zależności od mikrostruktury, obszaru spoiny, naprężeń własnych, środowiska oraz współczynnika obciążenia. ISO 12108 nie eliminuje tych zmiennych. Ustanawia kontrolowany sposób pomiaru określonej odpowiedzi materiału w zdefiniowanych warunkach.

Założenie liniowej sprężystości oznacza, że przyłożone współczynniki intensywności naprężeń są odpowiednimi parametrami opisującymi siłę napędową na wierzchołku pęknięcia. Założenie to staje się mniej wiarygodne w przypadku rozwinięcia się znacznej plastyczności, gdy próbka jest zbyt mała w stosunku do zamierzonej długości pęknięcia lub gdy pękanie zachodzi w wyniku mechanizmów, których nie opisuje obciążenie w przeważająco trybie I. Obciążenie w mieszanym trybie, historie obciążenia o zmiennej amplitudzie oraz znaczny wpływ naprężeń własnych wymagają dodatkowego uwzględnienia, a nie automatycznego zastępowania ich wynikiem badania ISO przy stałej amplitudzie.

ASTM E647-23 również obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego — od obszaru bliskiego progowi do niestateczności sterowanej przez $K_{\max }$ — przy czym wyniki wyraża się za pomocą $\Delta K$. ASTM opisuje badanie w sposób umożliwiający pomiary w szerokim zakresie wzrostu pęknięcia, w tym w obszarze, w którym stabilna propagacja przechodzi w szybkie pękanie. Wydany przez AFGROW w 2024 roku Damage Tolerance Design Handbook stwierdza, że dane wykorzystywane w standardowych analizach całkowania uszkodzeń są na ogół uzyskiwane w badaniach przy stałej amplitudzie na próbkach zawierających pęknięcia, z których większość jest objęta ASTM E647.

Żadne z tych stwierdzeń nie oznacza, że standardowe badanie bezpośrednio dostarcza prawa wzrostu pęknięcia dla kadłuba samolotu, zbiornika ciśnieniowego czy elementu mostu. Ocena konstrukcji musi uwzględniać współczynniki korekcyjne geometrii, widma obciążeń, skutki zamykania pęknięcia, naprężenia własne, możliwości inspekcji oraz przyjęte kryterium zniszczenia. FAA AC 25.571-1D umieszcza analizę wzrostu pęknięć i planowanie inspekcji w ramach oceny tolerancji uszkodzeń samolotów transportowych, natomiast 14 CFR 25.571 wymaga uwzględnienia prawdopodobnych lokalizacji uszkodzeń zmęczeniowych, analiz obciążeń cyklicznych popartych wynikami badań oraz mających zastosowanie progów inspekcji związanych ze wzrostem pęknięć.

Procedury wyznaczania progu i zmniejszania $\Delta K$

Obszar progowy jest miejscem, w którym pozornie niewielkie wybory proceduralne mogą mieć istotne konsekwencje interpretacyjne. $\Delta K_{\mathrm{th}}$ nie jest po prostu stałą liczbą, która pojawia się po określeniu gatunku stali. Zależy on od współczynnika obciążenia, środowiska, częstotliwości, historii próbki, długości pęknięcia, rozdzielczości pomiaru oraz kryterium stosowanego do stwierdzenia, że pęknięcie zatrzymało się lub rośnie z szybkością mniejszą od określonej wartości.

Badanie ze zmniejszaniem $\Delta K$ rozpoczyna się od większej siły napędowej, a następnie stopniowo zmniejsza się $\Delta K$, zazwyczaj przez zmianę zakresu obciążenia w miarę wydłużania się pęknięcia. Celem badania jest ustalenie punktu, w którym zmierzona szybkość wzrostu osiąga przyjęte kryterium progowe. Podejście to jest efektywne, jednak na wynik mogą wpływać rozwój zamykania pęknięcia, skutki historii obciążenia oraz niewystarczający czas na ustalenie przez pęknięcie ustalonej odpowiedzi po każdym zmniejszeniu obciążenia. Dlatego progu uzyskanego podczas zmniejszania obciążenia nie można interpretować bez uwzględnienia procedury, w której został wyznaczony.

TWI zwraca szczególną uwagę, że ASTM E647 i ISO 12108 stosują odmienne kryteria dotyczące progu i zmniejszania $\Delta K$. Różnica ta nie ma charakteru wyłącznie redakcyjnego. Jeżeli dwa laboratoria zbadają nominalnie identyczną stal zgodnie z tymi dwiema normami, zgłoszone przez nie wartości $\Delta K_{\mathrm{th}}$ mogą się różnić, ponieważ „próg” został operacyjnie zdefiniowany poprzez odmienne wymagania dotyczące przyrostu pęknięcia, pomiaru szybkości oraz redukcji obciążenia. Rozbieżność ta może następnie wpływać na dolny zakres dopasowania prawa Parisa oraz na obliczenia trwałości, w których krzywa jest całkowana w kierunku progu.

Ta sama ostrożność dotyczy górnego zakresu. ASTM E647 odnosi się do zakresu aż do niestateczności sterowanej przez $K_{\max }$, podczas gdy ISO 12108 wskazuje początek szybkiego niestatecznego pękania. Badanie może wykazywać pozornie gładką zależność $\mathrm{da}/\mathrm{dN}$ od $\Delta K$, dopóki pęknięcie nie zbliży się do warunków, w których dane dotyczące wzrostu w zakresie liniowo-sprężystym nie opisują już mechanizmu pękania będącego czynnikiem kontrolującym. Danych w pobliżu tego przejścia nie należy włączać do dopasowania prawa Parisa tylko dlatego, że osie wykresu pozostają logarytmiczne.

W artykule z 1963 roku Paris i Erdoğan opowiadali się za oceną zachowania podczas wzrostu pęknięcia w szerokich zakresach szybkości wzrostu oraz z wykorzystaniem wielu próbek. Stanowisko to jest bardziej rygorystyczne niż powszechna praktyka dopasowywania wąskiego odcinka centralnego i przedstawiania $C$ oraz $m$ jako uniwersalnych stałych dla stali. Raport NIST z 1976 roku dotyczący złączy spawanych ze stali zawierającej 5% Ni przedstawia konwencjonalne nazewnictwo stosowane w raportach: $C$ jest wyrazem wolnym, a $n$ nachyleniem wykresu logarytmicznego. W badaniu NIST z 2001 roku, obejmującym dwie ferrytyczno-perlityczne stale rurociągowe badane z użyciem zakrzywionych próbek M(T), rozrzut zmierzonych szybkości przypisano zmienności stałej $C$ w prawie Parisa. Analiza statystyczna jest użyteczna, ale pokazuje również, dlaczego pojedynczy dopasowany wyraz wolny może ukrywać znaczną zmienność między poszczególnymi próbkami.

Dlaczego norm nie można traktować jako wzajemnie wymiennych receptur

ISO 12108:2018 i ASTM E647-23 dotyczą ściśle powiązanych zagadnień inżynierskich, ale nie są identycznymi procedurami laboratoryjnymi. Wynik oznaczony jako „ASTM E647” zawiera informacje o procedurze ASTM, konfiguracji próbki oraz kryteriach progowych lub kryteriach niestateczności stosowanych w tej procedurze. Wynik oznaczony jako „ISO 12108” ma inną historię proceduralną. Usunięcie tych oznaczeń pozostawia liczbę, która może wyglądać na porównywalną, lecz ukrywa warunki, w których ją uzyskano.

Ma to szczególne znaczenie przy przenoszeniu danych do oprogramowania do oceny tolerancji uszkodzeń. AFGROW i podobne narzędzia całkują zależność wzrostu pęknięcia w funkcji historii naprężeń konstrukcji, często wykorzystując skorygowany geometrycznie $\Delta K$ oraz wybrane prawo $\mathrm{da}/\mathrm{dN}$. Całkowanie wzmacnia różnice systematyczne: niewielka zmiana w obszarze małych szybkości wzrostu może zmienić przewidywaną liczbę cykli do momentu inspekcji, natomiast odmienne ujęcie obszaru dużych wartości $K_{\max }$ może zmienić przewidywane zbliżanie się do pękania.

Normy należy zatem porównać przed połączeniem zbiorów danych. Należy sprawdzić współczynnik obciążenia, środowisko, częstotliwość, typ próbki, metodę pomiaru długości pęknięcia, historię zmniejszania $\Delta K$, definicję progu oraz sposób traktowania niestatecznego wzrostu. Jeżeli zamierzone zastosowanie obejmuje spoinę, płytę anizotropową, obciążenie o zmiennej amplitudzie lub znaczną plastyczność, wynik laboratoryjny wymaga również wyraźnego uzasadnienia jego przeniesienia na konstrukcję.

Nie ma podstaw naukowych do stwierdzenia, że ISO 12108 powszechnie zastępuje ASTM E647 ani że zachodzi sytuacja odwrotna. Uzasadnionym wyborem jest norma, której zakres i procedury odpowiadają materiałowi, obciążeniu oraz decyzji dotyczącej konstrukcji. Jeżeli porównuje się wyniki uzyskane według obu norm, zależność od metody należy zgłosić jako część materiału dowodowego, a nie traktować jej jako szumu eksperymentalnego.

Materiały do badań: blacha, pręt, stal rurociągowa i złącza spawane

Stale rurociągowe ferrytyczno-perlityczne[5][5] Fatigue Crack Growth Rates in Pipeline Steels Using Curved M(T) Specimens. National Institute of Standards and Technology. NIST publication, 2001.

Stal rurociągowa nie jest adekwatnie opisywana za pomocą jednej krzywej prawa Parisa. Przydatnym przykładem jest opublikowane w 2001 roku badanie National Institute of Standards and Technology (NIST), Fatigue Crack Growth Rates in Pipeline Steels Using Curved M(T) Specimens. W ramach tych prac zbadano dwie ferrytyczno-perlityczne stale rurociągowe, stosując zakrzywione próbki do rozciągania z pęknięciem centralnym, czyli próbki M(T). Zakrzywioną geometrię wybrano w celu odtworzenia pewnych aspektów badań wzrostu pęknięć w rurze przy jednoczesnym zachowaniu kontrolowanej zależności między przyłożonym obciążeniem a współczynnikiem intensywności naprężeń.

Znaczenie tego wyniku wynika nie tyle z wartości liczbowych stałych, ile ze sposobu potraktowania rozrzutu wyników. NIST przeanalizował zmienność zmierzonych szybkości wzrostu pęknięć, przypisując ją stałej $C$ prawa Parisa. W standardowym równaniu

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m,$$

$a$ oznacza długość pęknięcia, $N$ — liczbę cykli, $\Delta K$ — zakres współczynnika intensywności naprężeń, $C$ — wyraz wolny, a $m$ — nachylenie wykresu log–log zależności $\mathrm{da}/\mathrm{dN}$ od $\Delta K$. Zmiana $C$ przesuwa dopasowaną krzywą w pionie; nie oznacza to, że stal ma jedną dokładną, niezależną od geometrii szybkość wzrostu pęknięcia.

Rozróżnienie to ma istotne znaczenie w przypadku ferrytyczno-perlitycznych wyrobów w postaci blach i rur. Ferryt stanowi stosunkowo miękką osnowę, natomiast kolonie perlitu, pasmowość i wtrącenia mogą wpływać na odchylanie pęknięcia oraz lokalną odporność. Obróbka termomechaniczna może prowadzić do powstania różnych struktur ziarnistych wzdłuż grubości blachy oraz różnych właściwości w kierunku walcowania i w kierunku poprzecznym. Formowanie rur i spawanie szwu wprowadzają dodatkowe zmienne. Dwie stale o podobnej wytrzymałości na rozciąganie i nominalnym składzie chemicznym mogą zatem wykazywać różne wartości $C$, $m$, zachowanie progowe i rozrzut wyników.

Badanie NIST z zastosowaniem próbek M(T) pokazuje również, dlaczego zakrzywionej próbki reprezentatywnej dla rury nie należy traktować jako zamiennej z próbką kompaktową do rozciągania, wyciętą z płaskiej blachy. Geometria wpływa na więzy na froncie pęknięcia, podatność, kalibrację współczynnika intensywności naprężeń oraz część frontu pęknięcia mającą kontakt z lokalną strukturą materiału. Uzyskane dane mogą nadal służyć do dopasowania równania Parisa, jednak opis badania musi wskazywać formę próbki i metodę kalibracji.

ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego — od obszaru w pobliżu progu aż do niestabilności sterowanej przez $K_{\max }$. ISO 12108:2018 określa badania od $\Delta K_{\mathrm{th}}$ do początku gwałtownego niestabilnego pękania, przede wszystkim dla izotropowych materiałów metalicznych poddanych głównie liniowo-sprężystemu obciążeniu w trybie I, przy stałym współczynniku asymetrii cyklu $R$. Nie są to określenia zamienne. TWI wskazuje, że ASTM E647 i ISO 12108 stosują różne kryteria wyznaczania progu oraz kryteria dla zmniejszającego się $\Delta K$, dlatego przy podawaniu wartości progowej należy podać obok niej również normę i procedurę.

Złącza spawane ze stali zawierającej 5% Ni[6][6] Fatigue Crack Growth in 5% Ni Steel Weldments. National Institute of Standards and Technology. NBSIR 76-843, 1976.

Raport NIST Fatigue Crack Growth in 5% Ni Steel Weldments, opublikowany w 1976 roku jako NBSIR 76-843, pokazuje ten sam problem w przypadku materiału spawanego. Jego dane dopasowano za pomocą potęgowego prawa Parisa, przy czym $C$ opisano jako wyraz wolny, a $n$ jako nachylenie wykresu log–log zależności $\mathrm{da}/\mathrm{dN}$ od $\Delta K$. Oznaczenie różni się od częściej stosowanego $m$, lecz dopasowana zależność jest tego samego rodzaju empirycznym prawem.

Złącze spawane jest układem materiałowym, a nie jednorodną stalą. W stali zawierającej 5% Ni metal rodzimy zachowuje historię obróbki odziedziczoną po macierzystej blasze; metal spoiny krzepnie z ciekłego jeziorka i rozwija własny skład, układ segregacji oraz strukturę krzepnięcia; strefa wpływu ciepła (HAZ) ulega przemianom cieplnym, lecz nie zostaje ponownie stopiona. Odporność na wzrost pęknięcia może się gwałtownie zmieniać w poprzek tych obszarów. Próbka, w której pęknięcie przechodzi przez granicę wtopienia, może wykazywać szybkość wzrostu kontrolowaną przez najmniej odporny odcinek, przez odchylanie pęknięcia na granicy lub przez wpływ naprężeń szczątkowych, których nie ma w odizolowanym badaniu metalu rodzimego.

W związku z tym w przypadku wyniku uzyskanego dla złącza spawanego należy określić, czy pęknięcie propagowało w metalu rodzimym, metalu spoiny, HAZ, czy też w więcej niż jednym z tych obszarów. Należy również podać proces spawania, klasyfikację materiału dodatkowego, obróbkę cieplną po spawaniu, orientację spoiny oraz miejsce pobrania próbki. Określenie wyniku jako „krzywej dla stali zawierającej 5% Ni” ukrywa właśnie te cechy, które decydują o sposobie jego wykorzystania.

Dopasowanie krzywej nie zastępuje opisu materiału. Zależność Parisa i Erdoğana z 1963 roku miała stanowić empiryczny opis w użytecznym zakresie szybkości wzrostu pęknięcia, a ocena miała opierać się na szerokich zbiorach danych i wielu próbkach, a nie na wąskim odcinku przedstawianym jako uniwersalna stała materiałowa. W pobliżu progu dominować mogą zamykanie pęknięcia, chropowatość, naprężenia szczątkowe i bariery mikrostrukturalne. Przy wysokich wartościach $\Delta K$ wpływ $K_{\max }$, plastyczności i niestabilnego pękania sprawia, że proste dopasowanie według prawa Parisa jest nieodpowiednie.

Wpływ mikrostruktury, orientacji i obszaru spoiny

Blachy, pręty i rury należy identyfikować według postaci wyrobu, gatunku, obróbki cieplnej i orientacji próbki. W przypadku blach powszechnie stosowane oznaczenia, takie jak ASTM A516/A516M Grade 70 lub ASTM A572/A572M Grade 50, same w sobie nie definiują zachowania podczas wzrostu pęknięcia zmęczeniowego; nadal znaczenie mają grubość, sposób przetwarzania, populacja wtrąceń i orientacja. Próbka z pręta może obejmować strukturę poddaną obróbce w kierunku wzdłużnym, podczas gdy próbka z blachy w kierunku poprzecznym może w odmienny sposób przecinać wydłużone wtrącenia i pasma ferrytyczno-perlityczne. Sposób raportowania powinien jednoznacznie określać kierunki, na przykład L-T, T-L lub S-L, zamiast opierać się wyłącznie na określeniu „wzdłużny”.

Anizotropia zmienia zarówno szybkość wzrostu pęknięcia, jak i poprawność kalibracji współczynnika intensywności naprężeń. Kształt frontu pęknięcia może stać się nieregularny, gdy przemieszczający się front napotyka pasma, struktury ziaren byłego austenitu, wtrącenia lub dendryty metalu spoiny. Może to utrudniać pomiar długości pęknięcia i zwiększać pozorny rozrzut wyników bez zmiany podstawowego obciążenia cyklicznego.

Naprężenia szczątkowe mają równie istotne znaczenie w złączach spawanych i rurach formowanych. Nominalny współczynnik asymetrii cyklu $R$ nie opisuje w pełni lokalnego efektywnego współczynnika, gdy występują rozciągające lub ściskające naprężenia szczątkowe. Obróbka skrawaniem próbki może spowodować odprężenie; zachowanie stanu po spawaniu może te naprężenia utrzymać. Wybór ten zmienia zamykanie pęknięcia, $\Delta K_{\mathrm{eff}}$, pomiary wartości progowej, a niekiedy także pozorne nachylenie krzywej Parisa.

W ocenie konstrukcji badana krzywa jest danymi wejściowymi, a nie kompletną prognozą. AFGROW stwierdza, że standardowe dane do całkowania uszkodzeń pochodzą na ogół z badań przy stałej amplitudzie na próbkach z pęknięciami, najczęściej prowadzonych zgodnie z procedurami ASTM E647. Metodyka oceny tolerowania uszkodzeń opracowana przez FAA/Volpe National Transportation Systems Center oraz dokument FAA AC 25.571-1D wiążą takie dane z analizą wzrostu pęknięć, oceną obciążeń powtarzalnych i planowaniem inspekcji. Zgodnie z 14 CFR 25.571 oceny te muszą uwzględniać prawdopodobne lokalizacje i mechanizmy uszkodzeń zmęczeniowych, poparte badaniami dowody dotyczące obciążeń oraz obowiązujące progi inspekcji związane ze wzrostem pęknięć. Wiarygodna krzywa dla stali zawiera zatem informacje o geometrii próbki, środowisku, współczynniku asymetrii cyklu, orientacji, obszarze spoiny, stanie naprężeń szczątkowych i zakresie dopasowania. Bez tych zmiennych $C$ i $m$ są współrzędnymi na wykresie, a nie przenośnymi właściwościami „stali”.

Geometria próbki i praktyka obciążania

Szybkość wzrostu pęknięcia zmęczeniowego nie jest niezależną właściwością mierzoną niezależnie od próbki badawczej. Oblicza się ją na podstawie zaobserwowanego przyrostu pęknięcia oraz zakresu współczynnika intensywności naprężeń, ΔK, uzyskanego z rozwiązania właściwego dla danej geometrii. Zmiana kształtu próbki, położenia pęknięcia, grubości, układu obciążenia lub współczynnika sił może zmienić zgłaszaną wartość $\mathrm{da}/\mathrm{dN}$, nawet gdy stal pozostaje niezmieniona. Dopasowana zależność Parisa,

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m,$$

opisuje zatem konfigurację badania i określony zakres warunków, a nie uniwersalną stałą dla danego gatunku stali. Paris i Erdoğan wprowadzili tę zależność w 1963 roku, przy czym $C$ jest wyrazem wolnym, a $m$ nachyleniem wykresu logarytmicznego szybkości wzrostu pęknięcia w funkcji ΔK. W niektórych sprawozdaniach zamiast $m$ stosuje się $n$.

ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego — od zakresu bliskiego progowi aż do niestabilności sterowanej przez $K_{\max }$ — natomiast ISO 12108:2018 obejmuje badania od zakresu progowego, ΔKth, do szybkiego niestabilnego pękania. Obie normy dotyczą przede wszystkim obciążenia w dominująco liniowo-sprężystym zakresie, w trybie I, jednak stosowane w nich procedury wyznaczania progu oraz zmniejszania ΔK różnią się. TWI wskazuje w szczególności na te różnice między ASTM E647 a ISO 12108. W konsekwencji wyniku nie można oddzielić od normy i procedury, które zastosowano do jego uzyskania.

Koncepcje próbek kompaktowych i próbek z centralnym pęknięciem rozciąganych

Próbki kompaktowe, czyli CT, zawierają wykonaną mechanicznie karbę inicjującą oraz wstępnie wytworzone zmęczeniowo pęknięcie między dwoma otworami obciążającymi. Przyłożona siła otwiera pęknięcie w trybie I, a długość pęknięcia mierzy się w miarę jego przemieszczania przez przewężenie. ASTM E647 podaje rozwiązanie dla współczynnika intensywności naprężeń w geometrii CT, zwykle zapisywane w postaci

$$K=\frac{P}{B\sqrt{W}}f(a/W),$$

gdzie $P$ oznacza siłę, $B$ — grubość, $W$ — charakterystyczną szerokość, $a$ — długość pęknięcia, a $f(a/W)$ — bezwymiarową funkcję geometrii. Dokładne wyrażenie oraz definicje wymiarów muszą być zgodne z normą; podstawienie równania dla próbki CT do innego typu próbki prowadzi do formalnie precyzyjnej, lecz fizycznie błędnej wartości ΔK.

Próbki z centralnym pęknięciem rozciągane, czyli M(T), mają postać szerokiego pasa obciążanego rozciąganiem, z centralnym pęknięciem, zwykle rozchodzącym się symetrycznie od szczeliny lub karby inicjującej. Rozwiązanie dla współczynnika intensywności naprężeń uwzględnia stosunek długości pęknięcia do szerokości próbki oraz przyłożone naprężenie nominalne. W przeciwieństwie do próbki CT układ M(T) może odwzorowywać szeroki arkusz lub płytę z centralnym pęknięciem na wskroś, co jest użyteczne, gdy analizowany problem konstrukcyjny obejmuje obciążenie rozciągające działające w poprzek dużego przewężenia.

Jak stwierdza AFGROW w Damage Tolerance Design Handbook, większość standardowych danych do całkowania uszkodzeń pochodzi z badań próbek z pęknięciami przy stałej amplitudzie. Maszyna przykłada powtarzającą się maksymalną i minimalną siłę przy wybranym współczynniku sił $R=P_{\min }/P_{\max }$, a długość pęknięcia rejestruje się w funkcji liczby cykli. Stała amplituda nie eliminuje potrzeby kontrolowania $R$, częstotliwości, przebiegu czasowego, środowiska ani zakresu długości pęknięcia. Jedynie ogranicza w bardziej precyzyjny sposób historię obciążenia.

Pęknięcie należy najpierw wytworzyć zmęczeniowo w kontrolowanych warunkach. Karb wykonany piłowaniem nie jest równoważny naturalnie zaostrzonemu pęknięciu zmęczeniowemu: jego promień wierzchołka, pole naprężeń resztkowych oraz lokalne odkształcenie plastyczne mogą zmieniać wyniki uzyskiwane we wczesnym etapie. Dane pochodzące z obszaru blisko karbu są zazwyczaj odrzucane lub traktowane ostrożnie do czasu, aż pęknięcie uzyska wymagany kształt frontu i zostaną spełnione określone warunki ważności wyników.

Zakrzywione próbki M(T)

Zakrzywiona próbka M(T) stanowi ważny przykład tego, dlaczego geometria zasługuje na większą uwagę, niż często poświęca jej równanie Parisa. W badaniu NIST z 2001 roku mierzono szybkość wzrostu pęknięć zmęczeniowych w dwóch stalach rurociągowych o strukturze ferrytyczno-perlitycznej, stosując zakrzywione próbki M(T). Krzywiznę dobrano tak, aby dokładniej niż płaski pas odwzorowywała warunki obciążenia ściany rury lub powłoki cylindrycznej. Próbka wymagała zatem rozwiązania dla współczynnika intensywności naprężeń odpowiedniego dla zakrzywionej konfiguracji, a nie wyrażenia dla płaskiej próbki M(T) skopiowanego z podręcznika.

NIST przeanalizował rozrzut zmierzonych szybkości, przypisując go zmienności stałej $C$ w prawie Parisa. Takie ujęcie jest użyteczne w analizie niepewności, lecz nie należy go interpretować jako dowodu, że $C$ jest stałą, niezmienną właściwością materiału. Część rozrzutu może wynikać z kształtu frontu pęknięcia, lokalnej mikrostruktury, pomiaru podatności, niewspółosiowości, naprężeń resztkowych oraz dokładności kalibracji dla próbki zakrzywionej. Dopasowanie może wchłonąć te efekty do $C$, a niekiedy także do nachylenia $m$, nie ujawniając ich odrębnych przyczyn.

Badania zakrzywionych próbek M(T) pokazują również różnicę między pomiarem a zastosowaniem konstrukcyjnym. Prawo wzrostu pęknięcia dopasowane na podstawie próbki o geometrii przypominającej rurę może wspierać ocenę rurociągu, lecz dopiero po sprawdzeniu przez analityka, czy stopień skrępowania, krzywizna, współczynnik sił, stan spoiny i środowisko elementu odpowiadają warunkom badania. Zastosowanie przez NIST stali rurociągowych o strukturze ferrytyczno-perlitycznej nie uprawnia do przenoszenia uzyskanych stałych na austenityczną stal nierdzewną, niskostopową stal ulepszaną cieplnie ani na złącze spawane.

Kontrola współosiowości, grubości i plastyczności

Niewspółosiowość wprowadza zginanie do nominalnie rozciąganego układu badawczego. Jedna strona pęknięcia może wówczas doświadczać innej intensywności naprężeń niż druga, co prowadzi do zakrzywienia lub nierównomierności frontu pęknięcia i zmienia lokalny współczynnik $R$. Należy sprawdzić uchwyty, sworznie, powierzchnie czołowe próbki i osie obciążenia, aby siła przechodziła przez zamierzoną linię działania. Podczas badania należy kontrolować prostoliniowość frontu pęknięcia; wyraźnie skośny front świadczy o tym, że przyjęte rozwiązanie dwuwymiarowe może już nie odwzorowywać próbki.

Grubość kontroluje stopień skrępowania. Cienka próbka umożliwia większy przepływ plastyczny w kierunku przez grubość i może wykazywać zachowanie w stanie płaskiego naprężenia, natomiast próbka o wystarczającej grubości w większym stopniu zbliża się do stanu płaskiego odkształcenia. Stan płaskiego odkształcenia zwiększa skrępowanie w wierzchołku pęknięcia i może ograniczać wpływ powierzchni swobodnych, jednak jest to założenie wymagające sprawdzenia wymiarów i obciążenia, a nie automatyczna konsekwencja zastosowania stali. Rowki boczne mogą pomagać w utrzymaniu prostego frontu pęknięcia, lecz jednocześnie zmieniają efektywną geometrię i muszą zostać uwzględnione w odpowiednim obliczeniu.

Liniowo-sprężysta mechanika pękania wymaga występowania małoskalowego uplastycznienia: strefa plastyczna w wierzchołku pęknięcia musi pozostać mała w porównaniu z przewężeniem i innymi istotnymi wymiarami. Jeżeli naprężenie w przekroju netto zbliża się do granicy plastyczności albo strefa plastyczna stanowi znaczącą część pozostałego przewężenia, samo ΔK nie jest już wystarczającym opisem. Duża plastyczność może powodować zamykanie pęknięcia, stępienie, rozdzieranie lub wzrost w mieszanym zakresie sprężysto-plastycznym. Badanie może nadal mieć wartość, jednak jego wyniku nie należy przedstawiać jako konwencjonalnego pomiaru zgodnego z prawem Parisa bez odpowiednich zastrzeżeń.

Granica ta ma znaczenie w projektowaniu. FAA AC 25.571-1D oraz 14 CFR 25.571 włączają analizę wzrostu pęknięć, dowody dotyczące obciążeń powtarzalnych, prawdopodobne miejsca uszkodzeń zmęczeniowych oraz planowanie inspekcji do oceny tolerancji uszkodzeń statków powietrznych transportowych. Takie analizy zależą od prawa wzrostu pęknięcia, którego geometria próbki, praktyka obciążania i granice ważności są znane, zanim $C$ i $m$ zostaną wprowadzone do obliczeń konstrukcyjnych.

Badanie progu i zakres bliski progowi

Znaczenie operacyjne ΔKth

Zakres progowy współczynnika intensywności naprężeń, ΔKth, nie jest pojedynczą, uniwersalnie przenośną stałą materiałową. Jest wynikiem operacyjnym uzyskanym w ramach określonego badania: pęknięcie poddaje się obciążeniu cyklicznemu, zakres przyłożonego współczynnika intensywności naprężeń zmniejsza się lub poddaje innemu sposobowi kontroli, a zmierzona szybkość wzrostu zbliża się do określonego kryterium niepropagacji. Podana wartość progowa opisuje zatem granicę zaobserwowaną dla danej geometrii próbki, współczynnika asymetrii cyklu naprężenia, środowiska, stanu powierzchni, częstotliwości, historii obciążenia i procedury redukcji danych.

To rozróżnienie ma znaczenie, ponieważ „brak zmierzonego wzrostu” nie musi oznaczać, że pęknięcie stało się fizycznie niezdolne do dalszej propagacji. Pęknięcie może zwiększyć swoją długość na krótkim odcinku, a następnie ulec zatrzymaniu, propagować się okresowo poniżej rozdzielczości układu pomiarowego albo wykazywać szybkość mniejszą od tej, którą laboratorium może odróżnić od niepewności pomiaru długości pęknięcia. Wartość progowa jest więc związana zarówno z czasem trwania i rozdzielczością obserwacji, jak i z kryterium przyjętym do uznania pęknięcia za nieaktywne.

ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego „od zakresu bliskiego progowi do niestabilności sterowanej przez Kmax”, przy czym wyniki wyraża się za pomocą zakresu współczynnika intensywności naprężeń w ujęciu liniowo-sprężystym, ΔK. ISO 12108:2018 określa badania od progowego zakresu współczynnika intensywności naprężeń ΔKth do początku szybkiego, niestabilnego pękania, przede wszystkim dla izotropowych materiałów metalicznych poddanych głównie liniowo-sprężystemu obciążeniu w trybie I przy stałym współczynniku asymetrii cyklu R. Zakresy tych norm częściowo się pokrywają, lecz nie oznacza to, że normy te można stosować zamiennie.

TWI wyraźnie wskazuje na różnice między kryteriami wyznaczania progu oraz kryteriami zmniejszania ΔK określonymi w ASTM E647 i ISO 12108. Różnica ta może zmienić punkt, w którym w sprawozdaniu z badania szybkość wzrostu pęknięcia zostaje zaklasyfikowana jako zachowanie progowe. Wartość uzyskana zgodnie z jedną normą powinna zatem zachować oznaczenie tej normy, konfigurację próbki i warunki badania podczas porównywania jej z inną wartością.

Zakres bliski progowi znajduje się również poza głównymi założeniami zależności Parisa. Paris i Erdoğan wprowadzili w swojej pracy z 1963 roku postać empiryczną

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m,$$

Na wykresie logarytmicznym zależności szybkości wzrostu pęknięcia od ΔK, C jest wyrazem wolnym, a m — w niektórych sprawozdaniach stosuje się również symbol n — jest nachyleniem. W środkowej, w przybliżeniu liniowej części krzywej szybkości wzrostu zależność ta może dobrze opisywać dane. W pobliżu progu na pozorny trend mogą jednak dominująco wpływać zamykanie się pęknięcia, bariery mikrostrukturalne i ograniczenia pomiarowe. Wąskiego odcinka użytego do dopasowania nie należy przedstawiać jako pełnego opisu zachowania stali.

Zmniejszanie obciążenia i metody zmniejszania ΔK

W badaniach progowych ΔK jest zwykle zmniejszane w miarę zwiększania długości pęknięcia; procedurę tę nazywa się zazwyczaj zmniejszaniem obciążenia lub metodą zmniejszania ΔK. Maszyna zmniejsza zakres przyłożonej siły stopniowo albo zgodnie z określonym harmonogramem, tak aby pęknięcie doświadczało coraz mniejszej siły napędowej. Po każdym przyroście monitoruje się długość pęknięcia, a badanie trwa do momentu, gdy zmierzona szybkość spełni odpowiedni warunek progowy.

Metoda ta stwarza praktyczny problem: pęknięcie zachowuje pamięć wcześniejszego obciążenia. Jeżeli ΔK jest zmniejszane zbyt szybko, pęknięcie może zachować za sobą strefę plastycznie odkształconego materiału i pozostawać zamknięte przez część nominalnego cyklu rozciągania. Jeżeli zmniejszanie jest zbyt wolne, próbka pochłania nadmierną liczbę cykli i może być narażona na oddziaływanie środowiska lub zmiany mikrostrukturalne niezwiązane z zamierzonym warunkiem progowym. Uzyskana szybkość nie jest zatem po prostu odpowiedzią na końcowy poziom obciążenia; odzwierciedla ona drogę obciążenia.

Dodatkową komplikację stanowią naprężenia szczątkowe. Próbka stali spawanej może zawierać rozciągające lub ściskające naprężenia szczątkowe, które zmieniają lokalną siłę napędową pęknięcia, nawet gdy obliczony na podstawie obciążenia zewnętrznego ΔK pozostaje niezmieniony. Ślady obróbki, wżery korozyjne, nieciągłości spoiny i chropowate powierzchnie przełomu mogą również tworzyć lokalne bariery lub sprzyjać przedwczesnej propagacji. W badaniu typu M(T) o zakrzywionej geometrii, obejmującym dwie stale rurociągowe o strukturze ferrytyczno-perlitycznej, NIST odnotował rozrzut zmierzonych szybkości i przeanalizował go, przypisując go zmienności stałej C w prawie Parisa. Takie ujęcie jest użyteczne w obliczeniach uszkodzeń, lecz nie przekształca zaobserwowanego rozrzutu w uniwersalną wartość progową.

ASTM E647 i ISO 12108 określają odmienne wymagania dotyczące zmniejszania ΔK i wyznaczania progu. Różnica ta wpływa na dopuszczalne historie obciążenia, akceptację danych oraz interpretację pęknięcia, które wydaje się zatrzymywać. Sprawozdanie laboratoryjne zawierające jedynie zapis „ΔKth =”, bez wskazania normy, współczynnika R, środowiska, częstotliwości, rodzaju próbki i kryterium progowego, pomija informacje niezbędne do oceny znaczenia tej liczby.

Z tego powodu wyniki progowe należy przedstawiać jako zdefiniowane badaniem zachowanie wzrostu pęknięcia, a nie jako właściwość całej klasy „stali”. Ten sam nominalny gatunek może wykazywać odmienne wyniki w materiale podstawowym, materiale w strefie wpływu ciepła i materiale spoiny. W sprawozdaniu NIST z 1976 roku dotyczącym złączy spawanych ze stali zawierającej 5% Ni dane dotyczące wzrostu pęknięcia dopasowano za pomocą prawa potęgowego Parisa, opisując C jako wyraz wolny, a n jako nachylenie; oznaczenie to pokazuje, że podczas porównywania zbiorów danych należy rejestrować również same konwencje dopasowania.

Zamykanie się pęknięcia, środowisko i pozorny brak wzrostu

Zamykanie się pęknięcia ma zasadnicze znaczenie dla interpretacji zakresu bliskiego progowi. Obliczony na podstawie obciążenia zewnętrznego zakres ΔK wykorzystuje minimalne i maksymalne obciążenia, lecz pęknięcie może nie otwierać się przy nominalnym obciążeniu minimalnym. Chropowate powierzchnie przełomu, osady tlenkowe, efekty śladu plastycznego wywołanego propagacją, produkty przemian oraz rozciągające lub ściskające naprężenia szczątkowe mogą zmieniać obciążenie otwarcia. Pęknięcie doświadcza wówczas mniejszego zakresu efektywnego, niż sugeruje nominalne ΔK. Niska pozorna wartość progowa może zatem odzwierciedlać zamykanie się pęknięcia, a nie samoistną niezdolność mikrostruktury stali do propagacji pęknięcia.

Współczynnik asymetrii cyklu R, definiowany na podstawie naprężenia minimalnego i maksymalnego, zmienia zarówno nominalny sposób obciążenia, jak i prawdopodobieństwo zamykania się pęknięcia. Znaczenie może mieć również częstotliwość. Przy niskiej częstotliwości wilgoć, wodór, tlen lub inne reaktywne składniki mają więcej czasu na dotarcie do wierzchołka pęknięcia. Przy wyższej częstotliwości reakcje środowiskowe mogą być ograniczone, choć nagrzewanie i dynamiczne oddziaływania maszyny mogą wprowadzać inne błędy. Powietrze laboratoryjne, woda o wysokiej czystości, roztwór soli i środowiska charakterystyczne dla warunków eksploatacji mogą prowadzić do odmiennego zachowania w zakresie bliskim progowi. Temperaturę i wilgotność należy uwzględnić w dokumentacji badania.

Stan powierzchni nie jest szczegółem drugorzędnym. Próbka polerowana, próbka po obróbce mechanicznej oraz powierzchnia skorodowana lub spawana mogą powodować odmienne warunki inicjacji i zamykania się pęknięcia. Zakrzywienie czoła pęknięcia i niejednorodna propagacja mogą również sprawić, że zmierzona długość pęknięcia na powierzchni będzie słabym przybliżeniem efektywnego trójwymiarowego rozmiaru pęknięcia.

Pozorny brak wzrostu stanowi praktyczne zagrożenie. Pęknięcie może propagować się poniżej granicy rozdzielczości optycznej lub podatnościowej; chropowata trajektoria pęknięcia może sprawiać, że jego rzutowana długość wydaje się niezmieniona; zanieczyszczenia mogą mostkować powierzchnie pęknięcia; albo pole naprężeń szczątkowych może tymczasowo zatrzymać jedną część czoła pęknięcia, podczas gdy inna część nadal się propaguje. Procedura progowa, w której taki stan rejestruje się jako zerowy wzrost, może być odpowiednia dla określonego kryterium inżynierskiego, lecz nie należy jej traktować jako dowodu trwałego zatrzymania pęknięcia.

AFGROW wskazuje, że dane dotyczące wzrostu pęknięć zmęczeniowych wykorzystywane w standardowych analizach całkowania uszkodzeń pochodzą na ogół z badań przy stałej amplitudzie obciążenia próbek zawierających pęknięcia, najczęściej zgodnych z ASTM E647. W ocenie konstrukcji dane te stosuje się następnie w ramach wybranej reguły wzrostu pęknięcia, widma obciążenia i modelu inspekcji. FAA AC 25.571-1D oraz 14 CFR 25.571 umieszczają analizę wzrostu pęknięć, prawdopodobne lokalizacje uszkodzeń zmęczeniowych, wspieraną badaniami ocenę oddziaływania powtarzalnych obciążeń oraz planowanie inspekcji w ramach praktyki oceny tolerancji uszkodzeń statków powietrznych transportowych. Ten łańcuch decyzji wyjaśnia, dlaczego ΔKth należy dokumentować jako wynik pomiaru, a nie przenosić do oceny jako bezwarunkową stałą materiałową.

Trzy obszary krzywej wzrostu pęknięcia zmęczeniowego

Krzywa wzrostu pęknięcia zmęczeniowego przedstawia przyrost pęknięcia na cykl obciążenia, $\mathrm{da}/\mathrm{dN}$, w funkcji zakresu współczynnika intensywności naprężeń, $\Delta K$, zwykle na osiach logarytmicznych. W przypadku pęknięcia rozwijającego się głównie w trybie I, w warunkach liniowo-sprężystych, krzywa zazwyczaj dzieli się na trzy praktyczne obszary: wzrost w pobliżu progu, w przybliżeniu prostoliniowy odcinek opisany prawem potęgowym oraz gwałtowne przyspieszenie prowadzące do niestateczności pękania. Granice tych obszarów nie są stałymi materiałowymi. Zmieniają się w zależności od współczynnika asymetrii cyklu naprężenia $R$, środowiska, geometrii próbki, efektów zamykania pęknięcia, mikrostruktury, naprężeń własnych oraz procedury zastosowanej do uzyskania danych.

Ma to znaczenie, ponieważ dopasowane równanie Parisa opisuje wyłącznie tę część krzywej, z której uzyskano jego stałe. Nie zastępuje ono zmierzonej krzywej.

Obszar I: wzrost w pobliżu progu

Na końcu krzywej odpowiadającym małym wartościom $\Delta K$ pęknięcie rośnie bardzo powoli, a niewielka zmiana warunków badania może spowodować dużą zmianę zmierzonej szybkości wzrostu. Progowy zakres współczynnika intensywności naprężeń, powszechnie zapisywany jako $\Delta K_{\mathrm{th}}$, jest granicą operacyjną: poniżej tej wartości, przy określonym czasie trwania badania i zastosowanej procedurze, nie wykrywa się trwałego przyrostu pęknięcia. Nie jest to uniwersalny punkt, w którym każde pęknięcie w stalowym elemencie przestaje się rozwijać.

Kilka efektów komplikuje ten obszar. Zamykanie pęknięcia może zmniejszać efektywną część cyklu napędzającą pęknięcie. Chropowate powierzchnie przełomu, osady tlenkowe, naprężenia własne oraz plastyczność w wierzchołku pęknięcia mogą zmieniać pozorną wartość progu. Progu zmierzonego przy współczynniku asymetrii cyklu naprężenia $R=0.1$ nie można automatycznie zastąpić wartością zmierzoną przy $R=0.7$. Wilgotne powietrze, próżnia, słona woda oraz środowiska zawierające wodór mogą również prowadzić do uzyskania innych szybkości wzrostu, zwłaszcza w podatnych na takie oddziaływania stalach o wysokiej wytrzymałości i złączach spawanych.

Sama metoda badania wpływa na raportowaną wartość progu. ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego „od obszaru w pobliżu progu do niestateczności kontrolowanej przez Kmax”, natomiast ISO 12108:2018 określa badania od $\Delta K_{\mathrm{th}}$ do początku gwałtownego niestatecznego pękania, przede wszystkim dla izotropowych materiałów metalicznych poddanych obciążeniu głównie w liniowo-sprężystym trybie I, przy stałym współczynniku asymetrii siły $R$. ASTM E647 i ISO 12108 nie stosują identycznych kryteriów wyznaczania progu ani kryteriów zmniejszania $\Delta K$, na co zwraca uwagę TWI. Porównywanie dwóch opublikowanych wartości $\Delta K_{\mathrm{th}}$ bez sprawdzenia normy, historii obciążenia, środowiska i reguły zakończenia badania może zatem prowadzić do pozornego porównania materiałów.

Ekstrapolacja dopasowania według prawa Parisa w kierunku tego obszaru jest szczególnie myląca. Ponieważ prawo potęgowe nie zawiera wbudowanego progu, przewiduje niezerową szybkość wzrostu przy każdej dodatniej wartości $\Delta K$. Przewidywana szybkość może być znacznie większa od doświadczalnie obserwowanej szybkości w pobliżu progu, prowadząc do niepotrzebnie krótkiego obliczonego odstępu między kontrolami. Z drugiej strony dopasowanie wykonane na podstawie zanieczyszczonych lub słabo rozdzielonych danych z obszaru w pobliżu progu może zaniżać szybkość wzrostu.

Obszar II: zachowanie zgodne z prawem Parisa

W środkowej części krzywej $\mathrm{da}/\mathrm{dN}$ często spełnia w przybliżeniu prawo potęgowe:

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

P. C. Paris i F. Erdoğan przedstawili tę zależność empiryczną w swojej publikacji z 1963 roku. Na wykresie logarytm–logarytm $m$ jest nachyleniem, a $C$ wyrazem wolnym, chociaż w niektórych opracowaniach zamiast $m$ jako wykładnika stosuje się $n$. Raport NIST z 1976 roku dotyczący danych o wzroście pęknięć zmęczeniowych w złączach spawanych ze stali zawierającej 5% Ni opisuje $C$ jako wyraz wolny, a $n$ jako nachylenie dopasowania logarytm–logarytm.

Jest to obszar, w którym równanie Parisa jest najbardziej użyteczne. Krzywa jest wystarczająco oddalona zarówno od efektów progowych, jak i od przyspieszenia poprzedzającego końcowe pękanie, dzięki czemu przybliżenie prostą może opisywać istotny przedział wzrostu pęknięcia. Inżynierowie mogą wprowadzić zmierzoną zależność do obliczeń całkowania uszkodzeń, zwykle całkując

$$N={\int }_{a_i}^{a_f}\frac{\mathrm{da}}{C[\Delta K(a){]}^m},$$

przy uwzględnieniu zmiany strukturalnego współczynnika intensywności naprężeń wraz z rozmiarem pęknięcia. AFGROW stwierdza, że dane dotyczące szybkości wzrostu pęknięć zmęczeniowych stosowane w standardowych analizach całkowania uszkodzeń pochodzą zazwyczaj z badań przy stałej amplitudzie, prowadzonych na próbkach z pęknięciem, z których większość jest objęta normą ASTM E647.

Nawet w tym obszarze $C$ i $m$ nie są uniwersalnymi oznaczeniami odnoszącymi się wyłącznie do gatunku stali. Zależą od konfiguracji próbki, grubości, orientacji pęknięcia, współczynnika asymetrii cyklu naprężenia, sposobu sterowania obciążeniem, środowiska, metody redukcji danych oraz wybranego przedziału dopasowania. W badaniu NIST z 2001 roku, dotyczącym dwóch ferrytyczno-perlitycznych stali rurociągowych badanych z użyciem zakrzywionych próbek M(T), rozrzut zmierzonych szybkości analizowano, przypisując go zmianom stałej $C$ prawa Parisa. Takie ujęcie przypomina, że rozrzut jest częścią wyniku, a nie jedynie niedogodnością, którą należy usunąć.

Paris i Erdoğan opowiadali się za oceną zachowania podczas wzrostu pęknięcia w szerokim zakresie szybkości wzrostu oraz z wykorzystaniem wielu próbek. Wąski odcinek prostoliniowy może dawać pozornie dokładne nachylenie, jednocześnie ukrywając znaczne różnice między próbkami lub zmianę mechanizmu. Równanie to jest użytecznym modelem określonego obszaru, a nie pełnym opisem pękania zmęczeniowego.

Obszar III: przyspieszenie prowadzące do niestateczności

Przy dużych wartościach $\Delta K$ szybkość wzrostu pęknięcia przekracza wartość wynikającą z prawa potęgowego dla obszaru środkowego i gwałtownie rośnie, gdy maksymalny współczynnik intensywności naprężeń zbliża się do warunku odporności materiału na pękanie. Istotną zmienną sterującą jest często $K_{\max }$, a nie samo $\Delta K$. Plastyczność w wierzchołku pęknięcia zwiększa się, pozostały przekrój materiału ulega osłabieniu, a stabilny przyrost cykliczny przechodzi w kierunku rozdzierania lub niestatecznego pękania.

Prosta odpowiadająca prawu Parisa, przedłużona do tego obszaru, zwykle zaniża szybkość wzrostu. Błąd ten może być poważny, ponieważ obliczenia uszkodzeń są wrażliwe na końcową część historii wzrostu pęknięcia: pęknięcie może rozwijać się przez wiele cykli w obszarze II, lecz po rozpoczęciu przyspieszenia szybko przebić pozostały przekrój materiału. Modele uwzględniające $K_{\max }$, odporność na pękanie, zachowanie związane z zamykaniem pęknięcia lub odrębną regułę wzrostu w pobliżu pękania są bardziej odpowiednie, gdy ocena obejmuje ten zakres.

Pełna krzywa jest zatem zapisem pomiarów z określonymi obszarami granicznymi, a nie pojedynczą prostą. ASTM E647-23 i ISO 12108:2018 opisują badania obejmujące zakres od zachowania progowego do niestateczności, mimo że ich procedury się różnią. W ocenach konstrukcji należy następnie wybrać tę część krzywej i to prawo wzrostu pęknięcia, które są odpowiednie dla rzeczywistego elementu. W przypadku statków powietrznych FAA AC 25.571-1D i 14 CFR 25.571 łączą analizę wzrostu pęknięć z prawdopodobnymi miejscami uszkodzeń zmęczeniowych, oceną powtarzalnych obciążeń popartą badaniami oraz planowaniem kontroli. Dopasowane stałe Parisa są danymi wejściowymi tego procesu decyzyjnego, a nie substytutem obszarów progowego wzrostu i niestateczności krzywej.

Analiza danych, rozrzut i niepewność

Dopasowanie C oraz m lub n

Zależność Paris–Erdoğana jest zwykle zapisywana jako

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

gdzie $\mathrm{da}/\mathrm{dN}$ oznacza szybkość wzrostu pęknięcia zmęczeniowego, $\Delta K$ jest zakresem współczynnika intensywności naprężeń, $C$ jest współczynnikiem, a $m$ — wykładnikiem. W niektórych opracowaniach zamiast $m$ stosuje się $n$. Zapis się zmienia, lecz problem dopasowania pozostaje taki sam. Na osiach logarytmicznych

$$\log (\mathrm{da}/\mathrm{dN})=\log C+m\log \left(\Delta K\right),$$

zatem $m$, lub $n$, jest nachyleniem, a $\log C$ — wyrazem wolnym dla wybranych jednostek.

To ostatnie zastrzeżenie decyduje o wyniku. $C$ nie jest niezależną od jednostek stałą materiałową. Zmiana jednostek $\Delta K$, szybkości wzrostu pęknięcia lub cyklu obciążenia zmienia jego wartość liczbową, nawet gdy dane fizyczne pozostają niezmienione. Dopasowana wartość $m$ również może się zmienić po zmianie wybranego przedziału $\Delta K$, szczególnie jeśli przedział obejmuje przejście w zakresie dolnego progu lub górny zakres, w którym na wzrost pęknięcia wpływają $K_{\max }$, zamykanie pęknięcia, plastyczność i niestateczność. Równanie Parisa jest zatem lokalnym opisem empirycznym części krzywej wzrostu pęknięcia, a nie pełnym opisem — od progu aż do pęknięcia.

Artykuł P. C. Paris i F. Erdoğana z 1963 roku ustanowił tę zależność jako empiryczne prawo propagacji pęknięcia, jednak ich argumentacja nie uzasadnia wyznaczania dwóch „uniwersalnych” stałych na podstawie kilku punktów. Autorzy postulowali ocenę wzrostu pęknięcia w szerokich zakresach szybkości oraz z wykorzystaniem wielu próbek. Regresja przeprowadzona dla jednej próbki i w wąskim przedziale może dać atrakcyjną linię prostą, a jednocześnie prowadzić do nieprawidłowej oceny zachowania potrzebnego do prognozowania trwałości.

Redukcja danych rozpoczyna się przed regresją. Długość pęknięcia mierzy się przy powtarzających się cyklach obciążenia, a szybkość uzyskuje się metodą siecznej, metodą wielomianową przyrostową lub inną procedurą różniczkowania. Niewielka zmiana zmierzonej długości pęknięcia podzielona przez dużą liczbę cykli ogranicza wpływ losowego szumu, ale może ukryć krzywiznę i lokalne zmiany szybkości wzrostu. Krótki przedział zachowuje lokalne zachowanie, lecz wzmacnia szum odczytu. Wygładzanie wielomianowe wprowadza kolejny wybór: rząd wielomianu, szerokość okna oraz sposób postępowania na końcach zapisu. Różne metody różniczkowania zastosowane do tych samych pomiarów surowych mogą dawać różne wartości $\mathrm{da}/\mathrm{dN}$, zwłaszcza w pobliżu progu, gdzie wzrost jest powolny.

Obliczenie $\Delta K$ ma własne dane wejściowe. Znaczenie mają geometria próbki, długość pęknięcia, zakres przyłożonej siły, grubość oraz właściwe wyrażenie kalibracyjne. Norma ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego — od zakresu w pobliżu progu do niestateczności sterowanej przez $K_{\max }$ — przy czym wyniki wyraża się za pomocą liniowo-sprężystego $\Delta K$. Norma ISO 12108:2018 obejmuje badania od $\Delta K_{\mathrm{th}}$ do początku gwałtownego niestatecznego pękania, przede wszystkim dla izotropowych materiałów metalicznych poddanych głównie liniowo-sprężystemu obciążeniu w trybie I przy stałym współczynniku asymetrii cyklu $R$. TWI zwraca uwagę, że ASTM E647 i ISO 12108 stosują różne kryteria dotyczące progu oraz zmniejszania $\Delta K$. Regresji wykonanej na podstawie danych uzyskanych według jednej normy nie należy porównywać z wynikami uzyskanymi według innej tak, jakby procedury te były wzajemnie zamienne.

Wybór przedziału dopasowania wymaga wyraźnego udokumentowania. Dopasowanie ograniczone do centralnego zakresu Parisa może być właściwe dla obliczenia, które pozostaje w tym zakresie, lecz nie powinno być ekstrapolowane na zachowanie progowe ani na zachowanie w pobliżu niestateczności. Z drugiej strony uwzględnienie wszystkich punktów w jednej linii najmniejszych kwadratów może zmusić równanie do reprezentowania zakresów, dla których nie było przeznaczone. Wartości odstające wymagają czegoś więcej niż usuwania ich według wizualnej oceny. Punkt może odzwierciedlać błąd przepisywania danych, niestabilny rozwój frontu pęknięcia, chwilowe odchylenie sterowania obciążeniem lub rzeczywiste zachowanie materiału. Usunięcie go bez zapisanej reguły prowadzi do zaniżenia niepewności.

Rozrzut między próbkami

Dane dotyczące wzrostu pęknięć zmęczeniowych różnią się między próbkami, nawet gdy nominalny materiał, wymiary i ustawienia badania są takie same. Stale rurociągowe ferrytyczno-perlityczne dobrze ilustrują ten problem. W swoim badaniu z 2001 roku, obejmującym dwie stale rurociągowe badane za pomocą zakrzywionych próbek M(T), NIST przeanalizował „rozrzut zmierzonych szybkości”, przypisując go stałej $C$ prawa Parisa. Takie postępowanie jest użyteczne, ponieważ przekształca rodzinę krzywych szybkości w rozkład parametrów, lecz nie należy go interpretować jako dowodu, że każde źródło zmienności jest fizycznie fluktuacją $C$. Przypisana zmienność może również obejmować błąd pomiaru długości pęknięcia, zmienność sterowania siłą, szum różniczkowania oraz niedoskonałość modelu.

Metalurgia wnosi rzeczywistą zmienność. Pęknięcie może napotykać ziarna ferrytu, kolonie perlitu, wtrącenia, pasmowość, strefy wpływu ciepła spoiny, naprężenia resztkowe lub różne lokalne warunki skrępowania frontu pęknięcia. Tekstura i orientacja mikrostruktury zmieniają odchylanie oraz zamykanie pęknięcia. W przypadku spawanych elementów ze stali 5% Ni raport NIST z 1976 roku dopasowywał dane dotyczące wzrostu pęknięć zmęczeniowych za pomocą potęgowego prawa Parisa i opisywał $C$ jako wyraz wolny, a $n$ jako nachylenie wykresu log–log. Metalu spoiny, granic wtopienia i materiału strefy wpływu ciepła nie należy automatycznie łączyć w jedną populację tylko dlatego, że mają to samo nominalne oznaczenie gatunku stali.

Rozrzut statystyczny oznacza, że powtarzane obserwacje różnią się wokół wybranego modelu populacji. Zmienność fizyczna oznacza, że próbki lub obszary mikrostruktury rzeczywiście mają różne charakterystyki wzrostu pęknięcia. Błąd pomiaru laboratoryjnego jest jeszcze czymś innym: jest to niepewność wprowadzana przez ekstensometr, kalibrację podatności, optyczną lub elektryczną metodę pomiaru długości pęknięcia, pomiar siły, zliczanie cykli, ustawienie próbki lub układ sterowania. Kategorie te mogą się nakładać. Na przykład nieregularność frontu pęknięcia może zmniejszać powtarzalność sygnału potencjału elektrycznego, powodując pozorny rozrzut pomiarowy, a jednocześnie odzwierciedlając rzeczywiste trójwymiarowe zachowanie pęknięcia.

Liczba próbek wpływa zarówno na dopasowaną średnią, jak i na poziom ufności pokładany w ogonach rozkładu. Dwie próbki nie pozwalają z dużą pewnością ustalić częstości występowania powolnej lub szybkiej odpowiedzi wzrostu. Wiele próbek ujawnia, czy krzywe są w przybliżeniu równoległe, czy $C$ się zmienia, podczas gdy $m$ pozostaje podobne, oraz czy zmienia się samo nachylenie. Regresja łączna może ukryć te wzorce. Analiza hierarchiczna lub przynajmniej oddzielne dopasowanie dla każdej próbki, a następnie udokumentowane podsumowanie, pozwala odróżnić zmienność resztową w obrębie próbki od zmienności między próbkami.

Sterowanie obciążeniem jest częścią tych danych. Zamierzony współczynnik asymetrii cyklu $R$, zakres siły, przebieg czasowy, częstotliwość i liczba cykli muszą pozostawać w granicach określonych w procedurze. Krótkotrwałe przeciążenie może zmienić zamykanie pęknięcia i opóźnić jego późniejszy wzrost; niedostatecznie kontrolowana siła maksymalna może przesunąć pęknięcie w kierunku innego warunku $K_{\max }$. Taka odpowiedź nie jest po prostu losowym szumem i nie należy jej uśredniać bez odnotowania tego zdarzenia.

Niepewność pomiaru i krzywe konserwatywne

Rozdzielczość pomiaru długości pęknięcia jest często głównym ograniczeniem pomiarowym. Jeżeli przyrząd rozróżnia tylko niewielką liczbę przyrostów wzrostu pęknięcia w bloku cykli, obliczona pochodna może być zdominowana przez kwantyzację i ocenę operatora. W pobliżu $\Delta K_{\mathrm{th}}$ zgłoszona szybkość może zatem charakteryzować się dużą niepewnością względną, nawet gdy bezwzględny błąd długości pęknięcia wydaje się niewielki. Rozdzielczość, dryf kalibracji, krzywizna frontu pęknięcia, zmiany podatności oraz wybrany algorytm różniczkowania powinny być podawane wraz z danymi dotyczącymi szybkości.

Krzywa konserwatywna musi określać, przed czym ma chronić konserwatywność. Krzywa dolnej granicy może chronić przed zaniżeniem prognozowanego wzrostu pęknięcia dla określonej stali, orientacji, $R$, środowiska, rodzaju próbki i przedziału szybkości. Nie pozostaje ona automatycznie konserwatywna po zmianie grubości, stanu spoiny, temperatury, narażenia korozyjnego lub widma obciążenia. Jedną z metod jest dopasowanie zależności średniej i dodanie marginesu statystycznego opartego na rozrzucie reszt. Inną jest zastosowanie górnej granicy predykcji dla $\mathrm{da}/\mathrm{dN}$, z uwzględnieniem niepewności zarówno regresji, jak i populacji próbek. Wybrany poziom ufności lub tolerancji ma znaczenie; określenie „dolna granica” bez tej informacji nie ma odtwarzalnego znaczenia statystycznego.

W ocenie konstrukcji konserwatywność zależy również od sposobu całkowania krzywej. AFGROW stwierdza, że dane wykorzystywane w standardowych analizach całkowania uszkodzeń na ogół pochodzą z badań przy stałej amplitudzie na próbkach z pęknięciem, z których większość jest objęta normą ASTM E647. Obciążenie eksploatacyjne rzadko jest sekwencją o stałej amplitudzie, dlatego analityk musi zdecydować, jak uwzględniać przeciążenia, interakcję obciążeń, zmienne $R$, naprężenia resztkowe oraz przejścia między zakresami wzrostu pęknięcia. Konserwatywna krzywa Parisa nie naprawi nieodpowiedniego modelu uszkodzeń.

Ramy oceny tolerancji uszkodzeń opracowane przez FAA/Volpe National Transportation Systems Center dokonują tego samego rozróżnienia między zmierzonymi danymi propagacji a decyzjami dotyczącymi konstrukcji. FAA AC 25.571-1D dotyczy analizy wzrostu pęknięć, planowania inspekcji oraz wspieranej badaniami oceny powtarzalnego obciążenia, natomiast 14 CFR 25.571 wymaga, aby oceny tolerancji uszkodzeń uwzględniały prawdopodobne miejsca i mechanizmy uszkodzeń zmęczeniowych. Dopasowane $C$ i $m$ są danymi wejściowymi do tych decyzji, a nie samymi decyzjami. Ich niepewność należy uwzględnić w obliczeniach wzrostu pęknięcia, odstępie między inspekcjami oraz ocenie wytrzymałości resztkowej, zamiast ukrywać ją za dodatkowymi miejscami dziesiętnymi.

Środowisko, częstotliwość, naprężenie średnie i zmienna amplituda

Powietrze, wilgoć i środowiska korozyjne

Szybkość wzrostu pęknięcia zmęczeniowego jest właściwością określonych warunków badania, a nie po prostu określonego gatunku stali. Pęknięcie propagujące w powietrzu laboratoryjnym może zachowywać się inaczej niż pęknięcie wystawione na działanie wilgotnego powietrza, wody morskiej, kondensatu, gazu zawierającego wodór lub elektrolitu z ochroną katodową. Wilgoć może ograniczać efekty zamykania pęknięcia, sprzyjać rozpuszczaniu anodowemu na wierzchołku pęknięcia oraz ułatwiać wnikanie wodoru do stali podatnych na jego działanie. Zmierzona szybkość może zatem wzrosnąć, szczególnie w pobliżu wartości progowej, gdzie niewielkie zmiany chemii wierzchołka pęknięcia i jego zamykania mogą decydować o tym, czy pęknięcie wydłuży się podczas cyklu.

Dlatego TWI podkreśla znaczenie wyznaczania szybkości wzrostu pęknięcia w warunkach środowiskowych istotnych dla danego elementu. Badanie w suchym powietrzu laboratoryjnym nie może automatycznie reprezentować krzyżulca konstrukcji morskiej, zakopanego rurociągu, zbiornika magazynowego zawierającego wilgotny gaz ani konstrukcji statku powietrznego narażonej na wilgotność atmosferyczną. Środowisko oddziałuje również z częstotliwością: powolny cykl daje procesom korozyjnym i procesom związanym z wodorem więcej czasu na działanie podczas każdego przebiegu obciążenia. Częstotliwość, która wydaje się nieszkodliwa w powietrzu, może prowadzić do niekonserwatywnych wyników w środowisku korozyjnym.

Znaczenie mają również próbka i metoda badania. ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego — od zakresu bliskiego wartości progowej do niestabilności sterowanej przez Kmax — przy czym wyniki wyraża się za pomocą zakresu współczynnika intensywności naprężeń w ujęciu liniowo-sprężystym, ΔK. ISO 12108:2018 określa badania od progowego zakresu współczynnika intensywności naprężeń, ΔKth, do początku szybkiego niestabilnego pękania, przede wszystkim dla izotropowych materiałów metalicznych poddawanych przeważnie liniowo-sprężystemu obciążeniu w trybie I, przy stałym współczynniku obciążenia R. Normy te nie są wzajemnie wymiennymi recepturami badań środowiskowych. TWI wskazuje wprost na odmienne kryteria wartości progowej i zmniejszania ΔK w ASTM E647 oraz ISO 12108. Zgłoszona wartość ΔKth powinna zatem zawsze być opatrzona informacją o zastosowanej normie, historii obciążenia, środowisku i procedurze.

Zależność Parisa i Erdoğana, wprowadzona w 1963 roku, jest powszechnie zapisywana jako

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

gdzie $\mathrm{da}/\mathrm{dN}$ oznacza przyrost pęknięcia na cykl, $C$ jest wyrazem wolnym, a $m$, oznaczane w niektórych opracowaniach jako $n$, jest nachyleniem wykresu logarytmicznego szybkości wzrostu pęknięcia w funkcji ΔK. W raporcie NIST z 1976 roku dotyczącym spoin ze stali zawierającej 5% Ni parametr $C$ opisano jako wyraz wolny, a $n$ jako nachylenie. Żadna z tych stałych nie jest niezależna od środowiska badania. Badanie w środowisku korozyjnym może przesunąć krzywą, zmienić jej nachylenie i zmienić pozorną wartość progową; nie powoduje jedynie dodania stałej poprawki do krzywej wyznaczonej w powietrzu.

Wpływ częstotliwości i historii obciążenia

Współczynnik obciążenia $R=K_{\min }/K_{\max }$ opisuje naprężenie średnie w ujęciu mechaniki pękania. Przy stałym ΔK zmiana R powoduje zmianę $K_{\max }$, $K_{\min }$, plastyczności w obszarze wierzchołka pęknięcia oraz zamykania pęknięcia. Wyższe rozciągające naprężenie średnie zwykle ogranicza zamykanie pęknięcia i może zwiększać szybkość wzrostu, chociaż odpowiedź zależy od mikrostruktury stali, naprężeń własnych w spoinie, środowiska i zakresu szybkości wzrostu. W pobliżu wartości progowej wpływ R może być szczególnie duży. W obszarze Parisa może być mniejszy, lecz pozostaje istotny, gdy maksymalny współczynnik intensywności naprężeń zbliża się do granicy pękania lub przyspieszenia wzrostu.

Częstotliwość stanowi dodatkową, niezależną zmienną. W suchym powietrzu badanie stali przy 10 Hz może w części obszaru Parisa dawać szybkość podobną do uzyskiwanej przy 1 Hz, lecz podobieństwo to należy wykazać, a nie zakładać. W wodzie lub reaktywnym elektrolicie wolniejsze badanie może prowadzić do większego wspomagania środowiskowego wzrostu pęknięcia. Jeszcze bardziej odmienny jest test z okresem utrzymania maksymalnego obciążenia: podczas takiego okresu mogą wystąpić pełzanie w obszarze wierzchołka pęknięcia, transport wodoru, korozja naprężeniowa lub wzrost pęknięcia pod obciążeniem stałym, dlatego zliczanie wyłącznie cykli może ukrywać uszkodzenia zależne od czasu.

Historia obciążenia również zmienia pozorną szybkość wzrostu. Przeciążenie może wytworzyć za wierzchołkiem pęknięcia plastyczny ślad i zwiększyć zamykanie pęknięcia, powodując przejściowe opóźnienie wzrostu. To samo przeciążenie może także wywołać lokalne stępienie lub rozgałęzienie pęknięcia. Niedociążenie może usunąć efekt zamykania i przyspieszyć późniejszy wzrost. Skutki te zależą od stosunku przeciążenia, R, długości pęknięcia, stanu ograniczenia materiału oraz odstępów między zdarzeniami. Kolejność ma znaczenie. Blok dużych obciążeń, po którym następują małe obciążenia, nie jest na ogół równoważny tym samym blokom zastosowanym w odwrotnej kolejności.

Efekty te ujawniają ograniczenia związane z traktowaniem $C$ i $m$ jako uniwersalnych stałych stali. W badaniu NIST z 2001 roku, obejmującym dwie ferrytyczno-perlityczne stale rurociągowe badane z użyciem zakrzywionych próbek M(T), analizowano „rozrzut zmierzonych szybkości”, przypisując go stałej $C$ prawa Parisa. Takie podejście opisuje niepewność w obrębie określonego zbioru danych; nie dowodzi jednak, że rozrzut, wpływ częstotliwości, R, środowiska i widma obciążenia można łącznie przedstawić za pomocą jednego regulowanego wyrazu wolnego. Naprężenia własne w spoinie mogą również zmieniać efektywne naprężenie średnie, nie ujawniając się we współczynniku obciążenia przyłożonego z zewnątrz.

Dane o stałej amplitudzie w warunkach widm eksploatacyjnych

Większość standardowych analiz całkowania uszkodzeń rozpoczyna się od danych uzyskanych przy stałej amplitudzie na próbkach z pęknięciami. AFGROW stwierdza, że dane dotyczące szybkości wzrostu pęknięcia zmęczeniowego wykorzystywane w takich analizach są na ogół uzyskiwane w badaniach ze stałą amplitudą, przy czym większość badań obejmuje ASTM E647. Takie rozwiązanie jest praktyczne: oddziela pomiar $\mathrm{da}/\mathrm{dN}$ od późniejszego całkowania prawa wzrostu pęknięcia w odniesieniu do historii naprężeń konstrukcji. Nie stanowi ono dowodu, że obciążenie eksploatacyjne ma stałą amplitudę.

Krzywa Parisa dopasowana do wyników uzyskanych przy stałej amplitudzie odpowiada na ograniczone pytanie: jak szybko pęknięcie rosło przy zastosowanym ΔK, R, częstotliwości, środowisku, geometrii próbki i procedurze redukcji danych? Prognoza eksploatacyjna dotyczy innego zagadnienia. Musi uwzględniać sekwencję obciążeń, zmienny współczynnik geometryczny, naprężenia własne, opóźnienie, przyspieszenie, sposób traktowania wartości progowej oraz możliwość niestabilności przy wysokim $K_{\max }$. Stosowanie

$$\mathrm{da}/\mathrm{dN}=C(\Delta K{)}^m$$

cykl po cyklu może być uzasadnione w obszarze Parisa, gdy efekty wzajemnego oddziaływania obciążeń są niewielkie, a kalibracja jest odpowiednia. Staje się jednak wątpliwe w pobliżu ΔKth, podczas sekwencji przeciążeń oraz w pobliżu szybkiego pękania, gdzie potęgowe prawo Parisa nie opisuje pełnej krzywej wzrostu pęknięcia.

Paris i Erdoğan postulowali ocenę wzrostu pęknięcia w szerokich zakresach szybkości oraz z wykorzystaniem wielu próbek, zamiast przedstawiania wąskiego dopasowanego odcinka jako uniwersalnego prawa. Ta ostrożność pozostaje istotna w przypadku konstrukcji stalowych. Krzywa zmierzona w powietrzu przy jednej wartości R i jednej częstotliwości może służyć do interpolacji w porównywalnych warunkach; nie uzasadnia automatycznie ekstrapolacji do wilgotnego, powolnego widma eksploatacyjnego o zmiennej amplitudzie. Odrębne przedstawienie danych dotyczących wzrostu pęknięcia zmęczeniowego w spoinach ze stali zawierającej 5% Ni w opracowaniu NIST ilustruje, dlaczego należy zachować widoczność postaci materiału i stanu spoiny, w jakich zgłoszono wyniki.

Praktyka dotycząca statków powietrznych transportowych wyraźnie pokazuje to rozróżnienie. FAA AC 25.571-1D umieszcza analizę wzrostu pęknięć, planowanie inspekcji i ocenę tolerancji uszkodzeń w ramach opartego na badaniach systemu wykazywania zgodności, natomiast 14 CFR 25.571 wymaga uwzględnienia prawdopodobnych lokalizacji i postaci uszkodzeń zmęczeniowych, analiz obciążeń powtarzalnych popartych dowodami z badań oraz, w stosownych przypadkach, progów inspekcji wyznaczanych na podstawie wzrostu pęknięć. Ocena konstrukcji wykorzystuje zatem prawo wybrane do określonego celu i zwalidowane względem odpowiednich obciążeń, a nie dopasowanie Parisa oderwane od warunków badania.

Poza prawem Parisa: prawa wzrostu pęknięć stosowane w analizie tolerancji uszkodzeń

Zależność Parisa–Erdoğana jest ważna, ale nie stanowi pełnego opisu wzrostu pęknięć zmęczeniowych. W ich publikacji z 1963 r. wprowadzono postać empiryczną

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

gdzie $\mathrm{da}/\mathrm{dN}$ oznacza szybkość wzrostu pęknięcia, $\Delta K$ — zakres współczynnika intensywności naprężeń, $C$ — współczynnik związany z wyrazem wolnym, a $m$ — nachylenie wykresu logarytmicznego w układzie log–log. W niektórych opracowaniach zamiast $m$ stosuje się $n$. Równanie to opisuje przede wszystkim zakres pośredni: przybliża w przybliżeniu prostoliniową część krzywej wzrostu pęknięcia w układzie log–log, leżącą między zachowaniem w pobliżu progu a gwałtownym przyspieszeniem poprzedzającym niestabilność pękania.

To ograniczenie ma znaczenie w analizie tolerancji uszkodzeń. ASTM E647-23 obejmuje wyznaczanie szybkości wzrostu pęknięcia zmęczeniowego „od zakresu w pobliżu progu do niestabilności sterowanej przez Kmax”, natomiast ISO 12108:2018 określa badania od $\Delta K_{\mathrm{th}}$ do początku gwałtownego niestabilnego pękania. Dopasowana prosta Parisa przedstawia zatem tylko część zmierzonego zachowania, chyba że określone zostaną jej zakres i zamierzone zastosowanie. P. C. Paris i F. Erdoğan opowiadali się również za oceną w szerokich zakresach szybkości wzrostu pęknięcia oraz z użyciem wielu próbek, a nie za traktowaniem wąskiego dopasowanego odcinka jako uniwersalnej stałej materiałowej.

Późniejsze obliczenia są decyzją dotyczącą modelowania konstrukcji. Wybrane prawo musi odpowiadać danym z badań, historii naprężeń, geometrii pęknięcia, środowisku, współczynnikowi obciążenia oraz kryterium zniszczenia przyjętemu dla danego elementu.

Korekty uwzględniające próg i niestabilność pękania

Przy małych wartościach $\Delta K$ wzrost pęknięcia może zbliżać się do wartości progowej, zamiast nadal przebiegać zgodnie z prostą Parisa. Prosta ekstrapolacja według prawa Parisa może przewidywać wzrost pęknięcia w sytuacji, gdy badanie wykazuje znikomy wzrost lub brak propagacji. Jest to niebezpieczne, gdy prowadzi do zawyżenia trwałości, ale może również być nadmiernie konserwatywne, jeśli zostanie zastosowane poza zakresem pomiarowym. Wartość progowa nie jest jedną liczbą obowiązującą dla całego gatunku materiału: zależy od współczynnika obciążenia $R=K_{\min }/K_{\max }$, efektów zamykania pęknięcia, stanu powierzchni, historii obciążenia, naprężeń własnych, temperatury, korozji oraz procedury badawczej.

Pomiar wartości progowej zależy również od normy. TWI wskazuje, że ASTM E647 i ISO 12108 stosują różne kryteria wyznaczania progu oraz zmniejszania $\Delta K$. W konsekwencji wartości oznaczonej jako $\Delta K_{\mathrm{th}}$ nie należy przenosić między bazami danych bez sprawdzenia, w jaki sposób w badaniu zbliżano się do progu, jak potwierdzano wzrost pęknięcia oraz czy współczynnik siły pozostawał stały. ISO 12108:2018 jest ukierunkowana przede wszystkim na izotropowe materiały metalowe poddawane zasadniczo liniowo-sprężystemu obciążeniu w trybie I przy stałym współczynniku siły $R$; warunki te określają znaczenie uzyskanego wyniku.

W górnym zakresie krzywej wzrost pęknięcia przyspiesza, gdy $K_{\max }$ zbliża się do warunku odporności materiału na pękanie. Dopasowanie według prawa Parisa może nie uwzględniać tego zakrzywienia i może prowadzić do nierealistycznie długiej przewidywanej trwałości resztkowej, jeśli całkowanie jest kontynuowane aż do krytycznej długości pęknięcia. Dlatego w obliczeniach tolerancji uszkodzeń zazwyczaj osobno wprowadza się ograniczenie związane z niestabilnością pękania, stosując dopuszczalną wartość $K_{\max }$, krytyczny współczynnik intensywności naprężeń, całkę J lub kryterium pękania sprężysto-plastycznego, jeśli jest to właściwe, a także określone wymaganie dotyczące wytrzymałości resztkowej.

Niektóre prawa wzrostu pęknięcia uwzględniają zarówno składnik progowy, jak i składnik związany z niestabilnością. Równanie NASGRO jest szeroko stosowanym przykładem w oprogramowaniu do analizy wzrostu pęknięć w lotnictwie; jego składniki modyfikują licznik o charakterze podobnym do prawa Parisa w celu uwzględnienia zachowania progowego, a mianownik — w celu uwzględnienia zbliżania się do niestabilności pękania. Równania typu Formana również wprowadzają $K_{\max }$ lub składnik związany z odpornością na pękanie, tak aby wzrost przyspieszał w miarę zbliżania się do warunku krytycznego. Nie są one jednak automatycznym zamiennikiem badań. Ich parametry wymagają określonego zbioru danych, stanu materiału, środowiska oraz metody kalibracji. Prawo, które odtwarza dane z próbek kompaktowych rozciąganych wykonanych z jednego gatunku stali, może nie odtwarzać danych dotyczących metalu spoiny, strefy wpływu ciepła lub konstrukcji w pełnej skali.

Modele wrażliwe na współczynnik $R$ i $K_{\max }$

Postać prawa Parisa zawiera $\Delta K$, ale nie zawiera jawnie $R$. Dwa badania o tym samym $\Delta K$, lecz różnym $R$, mogą wykazywać różne szybkości wzrostu, ponieważ różnią się maksymalną intensywnością naprężeń, zamykaniem pęknięcia, naprężeniem średnim i historią strefy plastycznej. Efekt ten jest szczególnie istotny w przypadku stali wykazujących znaczne zamykanie pęknięcia, konstrukcji spawanych z naprężeniami własnymi oraz wtedy, gdy obciążenie eksploatacyjne zawiera wysokie wartości szczytowe naprężeń rozciągających.

Pierwszym rozszerzeniem jest korekta zależna od $R$, zastosowana do współczynnika lub wykładnika prawa Parisa. Zależności typu Walkera modyfikują na przykład efektywny zakres za pomocą funkcji $R$, często zapisywanej w postaci zawierającej $(1-R{)}^{\gamma }$. Wykładnik $\gamma$ jest dopasowywany, a nie wynika z oznaczenia gatunku stali. Taka zależność może opisywać serię badań przy stałej amplitudzie, ale nie należy traktować jej jako dowodu, że wszystkie historie obciążenia o tym samym nominalnym $R$ będą miały taki sam przebieg.

Modele typu Formana przywiązują większą wagę do $K_{\max }$, często łącząc składnik $(\Delta K{)}^m$ podobny do składnika prawa Parisa z mianownikiem, który zmniejsza się, gdy $K_{\max }$ zbliża się do odporności na pękanie. Modele typu NASGRO dodatkowo jawnie uwzględniają próg, współczynnik naprężeń oraz zakrzywienie w pobliżu niestabilności. Inne reguły inżynierskie, w tym korekty oparte na naprężeniu średnim lub zamykaniu pęknięcia, mogą być wybierane, gdy dostępne dane je uzasadniają. Ich nazwy nie eliminują potrzeby określenia jednostek, geometrii pęknięcia, współczynnika siły, środowiska oraz sposobu uwzględnienia górnej granicy skutków przeciążeń.

Różnica ta ma znaczenie praktyczne. Badanie przy stałej amplitudzie na próbce z pęknięciem dostarcza zależności wzrostu dla kontrolowanej ścieżki obciążenia. Rzeczywisty statek powietrzny lub konstrukcja ciśnieniowa podlega obciążeniom widmowym, opóźnieniu wzrostu po przeciążeniach, naprężeniom własnym, a czasami także oddziaływaniu wielu pęknięć. Zastosowanie korekty współczynnika $R$ poza zakresem jej kalibracji może dać precyzyjnie wyglądającą, lecz niepopartą dowodami prognozę trwałości.

Reguły omawiane w praktyce FAA i AFGROW

Damage Tolerance Assessment Handbook opracowany przez FAA/Volpe National Transportation Systems Center przedstawia propagację pęknięć zmęczeniowych jako zbiór możliwych do zastosowania reguł, a nie jedno obowiązkowe równanie. Jego metodyka obejmuje mechanikę pękania, dane dotyczące wzrostu pęknięć, całkowanie widma obciążeń, wytrzymałość resztkową oraz odstępy między przeglądami, przy czym wybór prawa wzrostu zależy od celu oceny i dostępnych dowodów. AFGROW podobnie stwierdza, że dane wykorzystywane do standardowych analiz całkowania uszkodzeń są na ogół uzyskiwane w badaniach przy stałej amplitudzie na próbkach z pęknięciami, a większość takich badań obejmuje ASTM E647.

W praktyce analityk może zastosować zależność Parisa dla zmierzonego zakresu pośredniego, regułę skorygowaną pod kątem progu dla wzrostu przy małych wartościach $\Delta K$, równanie wrażliwe na współczynnik $R$ lub $K_{\max }$ dla zmiennego naprężenia średniego oraz kryterium niestabilności pękania do zakończenia całkowania. Oprogramowanie może oferować modele Formana, Walkera, NASGRO, interpolację tabelaryczną lub inne opcje, ale lista dostępnych modeli nie stanowi specyfikacji materiału. Analityk musi udokumentować, którą regułę wybrano, w jaki sposób dopasowano jej współczynniki, jakie dane wyłączono oraz czy stosowane jest zachowanie dolnej granicy, czy wartości średnie.

Takie podejście jest zgodne z FAA AC 25.571-1D i 14 CFR 25.571. Oceny tolerancji uszkodzeń samolotów komunikacyjnych obejmują prawdopodobne lokalizacje i postacie uszkodzeń zmęczeniowych, analizy obciążeń powtarzalnych poparte dowodami z badań, obliczenia wzrostu pęknięć, wytrzymałość resztkową oraz planowanie przeglądów. Zasadnicze pytanie nie brzmi: „Jaka jest stała Parisa dla tej stali?”. Chodzi o to, czy wybrany opis wzrostu pęknięcia, połączony z geometrią, widmem obciążenia, sposobem uwzględnienia niepewności oraz możliwościami kontroli, zapewnia podstawę do wymaganej decyzji dotyczącej konstrukcji.

Przykłady NIST pokazują, dlaczego rozróżnienie to pozostaje istotne. W badaniu z 2001 r. dwóch ferrytyczno-perlitycznych stali rurociągowych zastosowano zakrzywione próbki M(T) i przeanalizowano rozrzut zmierzonych szybkości, przypisując go stałej $C$ prawa Parisa. W raporcie NIST z 1976 r. dotyczącym spoin ze stali zawierającej 5% Ni opisano $C$ jako wyraz wolny, a $n$ jako nachylenie dopasowania w układzie log–log. Dopasowane wielkości opisują konkretne doświadczenia i stany materiału. Nie czynią one z prawa Parisa niezależnej właściwości każdego elementu wykonanego z tego samego nominalnego gatunku stali.

Od danych dotyczących przyrostu pęknięcia do oceny tolerancji uszkodzeń

Całkowanie przyrostu pęknięcia od rozmiaru początkowego do krytycznego

Badanie przyrostu pęknięcia zmęczeniowego daje zależność opisującą szybkość wzrostu pęknięcia; ocena tolerancji uszkodzeń wykorzystuje tę zależność do prognozowania trwałości konstrukcji. Są to różne operacje. Mierzoną wielkością jest zwykle przyrost pęknięcia na cykl, $\mathrm{da}/\mathrm{dN}$, natomiast w ocenie należy określić, ile cykli eksploatacyjnych potrzeba, aby pęknięcie zwiększyło się od założonego lub wykrytego rozmiaru początkowego $a_0$ do rozmiaru granicznego $a_c$.

W zakresie Paris powszechnie stosuje się zależność

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m$$

gdzie $C$ jest wyrazem wolnym, a $m$, oznaczane w niektórych raportach jako $n$, jest nachyleniem wykresu log–log zależności $\mathrm{da}/\mathrm{dN}$ od $\Delta K$. P. C. Paris i F. Erdoğan wprowadzili tę zależność empiryczną w swojej publikacji z 1963 roku. Jest ona użyteczna, ale wyłącznie w zakresie reprezentowanym przez dane. Sama nie opisuje zachowania progowego, efektów zamknięcia pęknięcia, wzrostu pęknięcia w pobliżu niestateczności ani każdej sekwencji obciążenia.

Etap całkowania łączy zmierzoną zależność z rozwiązaniem dla współczynnika intensywności naprężeń właściwym dla rzeczywistego detalu konstrukcyjnego:

$$\Delta K=Y\left(a\right)\Delta \sigma \sqrt{\pi a}$$

W tym równaniu $Y\left(a\right)$ uwzględnia geometrię, kształt pęknięcia, powierzchnie swobodne, otwory pod łączniki, stopki spoin, usztywnienia i inne efekty brzegowe; $\Delta \sigma$ oznacza zakres przyłożonych naprężeń. Po podstawieniu otrzymujemy

$$N={\int }_{a_0}^{a_c}\frac{\mathrm{da}}{C[\Delta K(a){]}^m}$$

dla obciążenia o stałej amplitudzie w zakresie Paris. Wynik jest obliczoną liczbą cykli, a nie stałą materiałową. Zmiana geometrii pęknięcia, współczynnika naprężeń $R$, pola naprężeń resztkowych, warunków środowiskowych lub założonej wielkości wady początkowej zmienia prognozowaną trwałość, nawet gdy $C$ i $m$ pozostają niezmienione.

Rzeczywiste widma eksploatacyjne wymagają prowadzenia obliczeń blok po bloku obciążenia. Dla każdego zakresu naprężeń i naprężenia maksymalnego analityk oblicza $K_{\max }$, $K_{\min }$ oraz $\Delta K$, wybiera regułę przyrostu pęknięcia, wyznacza nowy rozmiar pęknięcia, aktualizuje współczynnik geometryczny i powtarza ten proces. AFGROW stwierdza, że dane dotyczące szybkości przyrostu pęknięć zmęczeniowych wykorzystywane w standardowych analizach całkowania uszkodzeń uzyskuje się zazwyczaj w badaniach próbek z pęknięciami poddawanych obciążeniu o stałej amplitudzie, przy czym większość takich badań obejmuje norma ASTM E647. Model trwałości eksploatacyjnej rozszerza zatem wyniki laboratoryjne na problem obciążenia o zmiennej amplitudzie; nie odtwarza po prostu badania.

Wybrane prawo szybkości wzrostu musi również odpowiadać części krzywej, przez którą przebiega analizowane pęknięcie. ASTM E647-23 obejmuje wyznaczanie szybkości przyrostu pęknięć zmęczeniowych — od zakresu bliskiego progowi aż po niestateczność sterowaną przez $K_{\max }$. ISO 12108:2018 określa badania od progowego zakresu współczynnika intensywności naprężeń, $\Delta K_{\mathrm{th}}$, do początku gwałtownego pękania niestatecznego, przede wszystkim dla izotropowych materiałów metalicznych poddawanych obciążeniu w dominująco liniowo-sprężystym trybie I, przy stałym współczynniku obciążenia $R$. TWI podaje, że ASTM E647 i ISO 12108 stosują różne kryteria wyznaczania progu oraz kryteria dla malejącego $\Delta K$. W związku z tym wartości progowej nie można przenosić między normami bez sprawdzenia sposobu jej wyznaczenia.

Równie istotny jest górny punkt zakresu. W miarę wzrostu pęknięcia $K_{\max }$ może zbliżyć się do odporności materiału na pękanie lub do obniżonej wartości dopuszczalnej, właściwej dla danego stopnia skrępowania, temperatury, grubości i środowiska. Rozmiar krytyczny $a_c$ jest zatem wyznaczany przez wymaganie dotyczące wytrzymałości resztkowej, a nie tylko przez koniec dopasowania prawa Paris. W przypadku prostego sprawdzenia liniowo-sprężystego warunek można zapisać jako

$$K_{\max }\left(a_c\right)=K_{\mathrm{allow}}$$

jednak $K_{\mathrm{allow}}$ może być niższe niż laboratoryjna wartość $K_{\mathrm{IC}}$, gdy na konstrukcję wpływają plastyczność, geometria, wzajemne oddziaływanie innego pęknięcia, korozja lub niestateczne rozdzieranie. Analiza wytrzymałości resztkowej ma na celu ustalenie, czy element z pęknięciem może nadal przenosić określone obciążenie graniczne lub niszczące. Całkowanie przyrostu pęknięcia określa, ile czasu potrzeba, aby osiągnąć rozmiar, przy którym wymaganie to nie jest już spełnione. Potrzebne są oba rodzaje analizy.

Założenia dotyczące wady początkowej determinują pierwszą część całki. W przypadku stalowego detalu spawanego ocenę można prowadzić od stopy spoiny lub wskazania braku przetopu; blasze walcowanej można przypisać wadę produkcyjną, wżer korozyjny albo pęknięcie o rozmiarze odpowiadającym możliwościom inspekcji; próbka obrobiona mechanicznie może rozpoczynać analizę od celowo wytworzonej, ostrej szczeliny zmęczeniowej. Założona wada musi być powiązana z jakością wykonania, możliwościami inspekcji, koncentracją naprężeń i przyjętym kształtem pęknięcia. Traktowanie $a_0$ jako dowolnej, bardzo małej liczby może prowadzić do dużego i nieuzasadnionego zwiększenia prognozowanej trwałości, ponieważ początkowy etap wzrostu może obejmować wiele cykli, szczególnie w pobliżu progu.

Progi i odstępy inspekcji

Program inspekcji przekształca obliczenia przyrostu pęknięcia w decyzję dotyczącą utrzymania. Próg inspekcji jest rozmiarem pęknięcia lub stanem wykrywalnego uszkodzenia, przy którym oczekuje się, że dana metoda inspekcji wykryje istotną wadę z wymaganą skutecznością. Nie jest on automatycznie równy $\Delta K_{\mathrm{th}}$. Pęknięcie może znajdować się poniżej zakresu wzrostu progowego, a mimo to mieć znaczenie dla konstrukcji, jeżeli wytrzymałość resztkowa jest mała, jeżeli może połączyć się z sąsiednim pęknięciem lub jeżeli metoda inspekcji nie pozwala wiarygodnie określić jego rozmiaru.

Odstęp musi zapewniać wystarczający margines między wykrywalnym wzrostem a stanem krytycznym. Uproszczony przebieg można przedstawić następująco:

$$a_0\to a_{\det }\to a_c$$

Obliczona liczba cykli od $a_0$ do $a_{\det }$ opisuje czas dostępny na to, aby wada stała się wykrywalna. Liczba cykli od $a_{\det }$ do $a_c$ wyznacza dostępny odstęp między inspekcjami, z uwzględnieniem współczynników bezpieczeństwa, niepewności widma, błędu wymiarowania, ryzyka niewykrycia oraz konsekwencji pękania niestatecznego. Ocena może wymagać, aby odstęp był krótszy niż nominalnie obliczony okres wzrostu, tak aby pęknięcie wykryte w pobliżu granicy wykrywalności zostało naprawione przed osiągnięciem $a_c$.

Rozrzut wyników przyrostu pęknięcia należy uwzględniać w sposób jawny. W badaniu NIST z 2001 roku, dotyczącym dwóch ferrytyczno-perlitycznych stali rurociągowych badanych za pomocą zakrzywionych próbek M(T), analizowano „rozrzut zmierzonych szybkości”, przypisując go stałej $C$ prawa Paris. Podejście to ilustruje istotę problemu: dwa nominalnie identyczne elementy mogą nie podążać za średnią dopasowaną krzywą. W zależności od uzasadnienia bezpieczeństwa można zamiast niej wybrać konserwatywną krzywą górnej granicy, granicę ufności albo rozkład probabilistyczny. W raporcie NIST z 1976 roku dotyczącym spawanych elementów ze stali zawierającej 5% Ni opisano $C$ jako wyraz wolny, a $n$ jako nachylenie dopasowania log–log, jednak te dopasowane parametry nadal pozostają związane ze stanem elementu spawanego, procedurą badania i dopasowanym zakresem $\Delta K$.

Planowanie inspekcji zależy również od rodzaju uszkodzenia. Pęknięcie na wskroś w rozciąganym panelu, pęknięcie powierzchniowe przy stopie spoiny oraz liczne pęknięcia w pobliżu otworu pod łącznik nie mają wspólnego współczynnika geometrycznego ani tego samego problemu detekcji. Wzajemne oddziaływanie pęknięć, zmiana kierunku pęknięcia, naprężenia resztkowe, szczegóły naprawy i redystrybucja obciążenia mogą wymagać modelu wykraczającego poza obliczenia przyrostu pojedynczego, izolowanego pęknięcia z wykorzystaniem prawa Paris.

FAA AC 25.571-1D i 14 CFR 25.571

Tolerancja uszkodzeń konstrukcji samolotów transportowych sytuuje te obliczenia w ramach szerszego systemu wykazywania zgodności. FAA AC 25.571-1D zawiera wytyczne dotyczące oceny tolerancji uszkodzeń i zmęczenia konstrukcji samolotów transportowych, w tym analizy przyrostu pęknięć i planowania inspekcji. Okólnik doradczy stanowi wytyczne, a nie zastępstwo dla przepisu, lecz wskazuje dowody i powiązania analityczne oczekiwane w ocenie wykazującej zgodność.

Obowiązujący przepis, 14 CFR 25.571, wymaga uwzględnienia prawdopodobnych miejsc i mechanizmów uszkodzeń zmęczeniowych, w tym obszarów, w których cykliczne obciążenia, koncentracja naprężeń, połączenia, wycięcia, zamocowania lub szczegóły wykonania mogą inicjować pękanie. Wymaga również, w stosownych przypadkach, ocen powtarzalnego obciążenia popartych badaniami, zamiast polegania wyłącznie na abstrakcyjnej krzywej materiałowej. Wyniki badań mogą określać zachowanie reprezentatywnej konstrukcji, dróg przenoszenia obciążeń, napraw, połączeń i detali, których nie da się uchwycić w wyniku badania małej próbki.

W ramach tego systemu analityk określa prawdopodobne miejsca i rozmiary wad, stosuje reprezentatywne widma obciążeń w locie i na ziemi, oblicza przyrost pęknięcia, sprawdza wytrzymałość resztkową oraz ustala progi i odstępy inspekcji, gdy inspekcja jest głównym środkiem zabezpieczającym. Wynik oceny tolerancji uszkodzeń stanowi zatem łańcuch obejmujący: wyznaczanie przyrostu pęknięcia zgodnie z normami, całkowanie zależne od geometrii i widma obciążenia, ograniczenie wynikające z wytrzymałości resztkowej oraz poparty wynikami badań plan inspekcji. Stałe prawa Paris są elementem tego łańcucha, ale go nie definiują.

Lista kontrolna kontroli jakości dla wiarygodnego badania stali

Wynik badania szybkości wzrostu pęknięcia zmęczeniowego jest wiarygodny tylko wtedy, gdy inżynier może odtworzyć, co badano, w jaki sposób przeprowadzono badanie oraz w którym miejscu przedstawione równanie przestaje być poparte pomiarami. ASTM E647-23 obejmuje szybkość wzrostu pęknięcia zmęczeniowego od zakresu bliskiego progowi aż do niestabilności sterowanej przez $K_{\max }$, natomiast ISO 12108:2018 obejmuje badania od $\Delta K_{\mathrm{th}}$ do szybkiego niestabilnego pękania. Są to powiązane normy, ale nie są wzajemnie wymiennymi oznaczeniami: TWI wskazuje, że ASTM E647 i ISO 12108 stosują różne kryteria wyznaczania wartości progowej oraz zmniejszania $\Delta K$.

Dokumentacja materiału i próbki

Należy dokładnie zapisać oznaczenie stali zgodnie z certyfikatem. Określenia „stal węglowa” lub „stal wysokowytrzymała” nie są wystarczające. Należy podać, czy materiałem jest ASTM A36, ASTM A516 Grade 70, API 5L X65, EN 10025-2 S355J2, czy materiał o innym oznaczeniu, oraz dołączyć świadectwo materiałowe, numer wytopu lub odlewu, specyfikację wyrobu i postać wyrobu. Blacha, pręt walcowany na gorąco, rura bezszwowa, rura spawana i wyrób kuty mogą mieć różną mikrostrukturę, naprężenia własne i charakterystykę wzrostu pęknięć, nawet gdy nominalne wartości wytrzymałości wydają się podobne.

Dokumentacja powinna obejmować grubość wyrobu, szerokość, a w stosownych przypadkach także średnicę, kierunek walcowania, kierunek poprzeczny oraz orientację próbki, z zastosowaniem jednoznacznego oznaczenia, takiego jak L-T lub T-L. Należy określić miejsce pobrania próbki w obrębie blachy, rury lub złącza spawanego. W przypadku próbki spawanej należy wskazać proces spawania, klasyfikację materiału dodatkowego, metal spoiny, linię wtopienia, strefę wpływu ciepła oraz obszar metalu rodzimego. Wynik opisany jako dotyczący „stali spawanej”, bez wskazania lokalizacji pęknięcia, jest zbyt ogólny, aby można go było później wykorzystać w analizie konstrukcji.

Należy w pełni opisać obróbkę cieplną: normalizowanie, hartowanie i odpuszczanie, wyżarzanie, odprężanie lub stan dostawy; jeżeli są znane, należy zachować również temperaturę obróbki i czas wygrzewania. Należy uwzględnić granicę plastyczności, wytrzymałość na rozciąganie, wydłużenie, twardość, analizę chemiczną oraz, w stosownych przypadkach, wyniki próby udarności Charpy’ego. Obserwacje metalograficzne powinny identyfikować ferrytowo-perlityczne, bainityczne, martenzytyczne lub inne istotne składniki struktury, zamiast traktować samo oznaczenie gatunku jako pełny opis materiału.

Geometrię próbki należy opisać za pomocą oznaczenia normowego i wymiarów. Przykłady obejmują próbkę zwartego rozciągania, C(T), próbkę rozciąganą z pęknięciem centralnym, M(T), oraz próbkę do zginania z jednostronnym karbem, SE(B), wraz z podaniem grubości, szerokości, długości odcinka czynnego, promienia karbu i sposobu zamocowania. Należy zapisać metodę obróbki skrawaniem oraz końcowy stan powierzchni. Trzeba wskazać, czy zastosowano rowki boczne, a także podać ich głębokość, kąt i cel wykonania. W badaniu NIST dotyczącym dwóch rurociągowych stali ferrytyczno-perlitycznych zastosowano zakrzywione próbki M(T); szczegół ten ma znaczenie, ponieważ krzywizna zmienia sposób obliczania współczynnika intensywności naprężeń i nie można go usunąć z wyniku, podając wyłącznie $C$ i $m$.

Należy oddzielnie udokumentować wstępne wytwarzanie pęknięcia zmęczeniowego i badanie główne. Należy podać metodę wykonania karbu początkowego, współczynnik asymetrii obciążenia podczas wstępnego wytwarzania pęknięcia, częstotliwość, przebieg obciążenia, maksymalną siłę, końcową długość pęknięcia wstępnego oraz wszelkie wymagane kontrole prostoliniowości lub symetrii. Karb wykonany mechanicznie nie jest pęknięciem zmęczeniowym. Przejście od karbu do pęknięcia rozwiniętego w sposób naturalny może wpływać na pierwsze raportowane punkty danych.

Dokumentacja maszyny, środowiska i monitorowania pęknięcia

Należy zidentyfikować stanowisko badawcze, siłomierz, uchwyty, wersję oprogramowania oraz daty kalibracji. Należy zapisać kalibrację siły i przemieszczenia, kontrole osiowości, weryfikację podatności oraz metodę zastosowaną do obliczenia współczynnika intensywności naprężeń. Trzeba podać, czy współczynnik asymetrii obciążenia był stały:

$$R=\frac{K_{\min }}{K_{\max }}=\frac{P_{\min }}{P_{\max }}.$$

Należy również podać częstotliwość, przebieg obciążenia, maksymalną i minimalną siłę, tryb sterowania obciążeniem, metodę zliczania cykli oraz przerwy w badaniu. Wynik uzyskany przy stałej amplitudzie nie może być automatycznie przenoszony na widmo o zmiennej amplitudzie. AFGROW stwierdza, że dane wykorzystywane w standardowych analizach całkowania uszkodzeń pochodzą na ogół z badań przy stałej amplitudzie na próbkach z pęknięciem, co oznacza, że historia obciążenia jest częścią materiału dowodowego, a nie jedynie przypadkowym ustawieniem maszyny.

Środowisko badania wymaga równie precyzyjnego opisu. Należy zapisać, czy badanie przeprowadzono w powietrzu laboratoryjnym, czy w warunkach zanurzenia, a także temperaturę, wilgotność względną, ciśnienie — w stosownych przypadkach — tożsamość i stężenie roztworu, pH, zawartość rozpuszczonego tlenu, szybkość przepływu, wstępne kondycjonowanie próbki oraz czas ekspozycji. W przypadku badań w środowisku wodnym należy podać nazwę związku zamiast określenia „słona woda”, na przykład stężenie chlorku sodu w gramach na litr. Należy wskazać, czy występowały produkty korozji, ładowanie wodorem, ochrona katodowa lub dodatki inhibitorów. Należy również podać, czy pęknięcie badano w warunkach obwodu otwartego, czy przy kontrolowanych warunkach elektrochemicznych.

Należy wyjaśnić, w jaki sposób mierzono długość pęknięcia: metodą podatności, spadku potencjału, mikroskopem optycznym, mikroskopem pomiarowym, cyfrową korelacją obrazu lub inną metodą. Należy podać szczegóły kalibracji, rozdzielczość, zasady wygładzania danych, założenia dotyczące frontu pęknięcia oraz metodę zastosowaną do przeliczenia przyrostów długości pęknięcia na $\mathrm{da}/\mathrm{dN}$. Należy zapisać liczbę odczytów przypadających na przedział cykli oraz sposób postępowania z pominiętymi lub niestabilnymi odczytami. Jeżeli front pęknięcia był zakrzywiony, rozgałęziony lub asymetryczny, należy zachować surowe obserwacje zamiast raportować jedną, bezkrytycznie przyjętą długość.

Prezentacja krzywej bez nadmiernych twierdzeń

Przed przedstawieniem dopasowanego równania należy opublikować zmierzone wartości $\mathrm{da}/\mathrm{dN}$ oraz odpowiadające im wartości $\Delta K$, $K_{\max }$, długości pęknięcia i liczby cykli. Należy wskazać wykluczenia oraz ich przyczyny techniczne. Trzeba pokazać rozrzut wyników między próbkami, a nie tylko linię średnią. W badaniu NIST z 2001 r. dotyczącym stali rurociągowej przeanalizowano „rozrzut zmierzonych szybkości”, przypisując go stałej $C$ prawa Parisa; podejście to jest użyteczne tylko wtedy, gdy obserwacje leżące u jego podstaw pozostają widoczne.

Należy podać zakres dopasowania zarówno w odniesieniu do $\Delta K$, jak i $\mathrm{da}/\mathrm{dN}$, liczbę punktów, liczbę próbek, metodę regresji, sposób ważenia, podstawę logarytmu oraz przedział ufności lub predykcji. Zależność Parisa i Erdoğana z 1963 r. jest powszechnie zapisywana w postaci

$$\frac{\mathrm{da}}{\mathrm{dN}}=C(\Delta K{)}^m.$$

W tym równaniu $C$ jest wyrazem wolnym, a $m$, w niektórych opracowaniach oznaczane jako $n$, jest nachyleniem wykresu logarytmicznego. Należy podać jednostki $\Delta K$, $\mathrm{da}/\mathrm{dN}$ i $C$; zmiana MPa√m na ksi√in zmienia wartość liczbową wyrazu wolnego. Nie należy określać $C$ i $m$ jako uniwersalnych stałych stali. Opisują one określony stan materiału, orientację, geometrię, $R$, środowisko, procedurę pomiarową i przedział dopasowania.

Należy oznaczyć granice stosowalności: wymagania dotyczące stanu płaskiego odkształcenia, ograniczenia grubości i długości odcinka czynnego, założenia małoskalowego uplastycznienia, zakres długości pęknięcia, maksymalną wartość $K$, procedurę wyznaczania wartości progowej oraz początek niestabilności. ASTM E647-23 i ISO 12108:2018 różnią się pod względem procedur wyznaczania wartości progowej i zmniejszania $\Delta K$, dlatego przyjęta norma oraz odstępstwa od niej muszą być umieszczone obok krzywej.

Na koniec należy rozdzielić trzy elementy: zmierzone dane, dopasowane parametry prawa Parisa oraz dowolną ekstrapolowaną krzywą projektową. Krzywa projektowa może uwzględniać współczynniki rozrzutu, granice konserwatywne, założenia dotyczące wartości progowej lub inną regułę wzrostu pęknięcia na potrzeby oceny konstrukcji. Nie są to dodatkowe dane z badań. FAA AC 25.571-1D i 14 CFR 25.571 zaliczają analizę wzrostu pęknięć, prawdopodobne miejsca uszkodzeń, ocenę obciążeń powtarzalnych popartą badaniami oraz planowanie kontroli do elementów wykazania zgodności w ramach koncepcji tolerowania uszkodzeń; decyzje te wymagają założeń wykraczających poza dopasowanie laboratoryjne. Identyfikowalność jest częścią wyniku.

Źródła

1. [1]International Organization for Standardization. ISO 12108:2018: Metallic materials — Fatigue testing — Fatigue crack growth method. ISO standard, 2018. https://www.iso.org/standard/73809.html
2. [2]TWI. Fatigue crack growth testing. TWI technical knowledge, 2024. https://www.twi-global.com/technical-knowledge/job-knowledge/fatigue-crack-growth-testing
3. [3]AFGROW. Damage Tolerance Design Handbook. AFGROW Damage Tolerance Design Handbook, 2024. https://www.afgrow.net/applications/DTDHandbook/Sections/page7_2_2_1.aspx
4. [4]Electronic Code of Federal Regulations. 14 CFR 25.571: Damage-tolerance and fatigue evaluation of structure. 14 CFR, 2024. https://www.ecfr.gov/current/title-14/chapter-I/subchapter-C/part-25/subpart-C/section-25.571
5. [5]National Institute of Standards and Technology. Fatigue Crack Growth Rates in Pipeline Steels Using Curved M(T) Specimens. NIST publication, 2001. https://www.nist.gov/publications/fatigue-crack-growth-rates-pipeline-steels-using-curved-mt-specimens
6. [6]National Institute of Standards and Technology. Fatigue Crack Growth in 5% Ni Steel Weldments. NBSIR 76-843, 1976. https://nvlpubs.nist.gov/nistpubs/Legacy/IR/nbsir76-843.pdf