Równanie wzrostu pęknięć
Równanie wzrostu pęknięcia służy do obliczania rozmiaru pęknięcia zmęczeniowego powstającego pod wpływem obciążeń cyklicznych. Rozwój pęknięć zmęczeniowych może doprowadzić do katastrofalnej w skutkach awarii, szczególnie w przypadku samolotów. Równanie wzrostu pęknięć można wykorzystać do zapewnienia bezpieczeństwa, zarówno w fazie projektowania, jak i podczas eksploatacji, poprzez przewidywanie wielkości pęknięć. W konstrukcji krytycznej obciążenia można rejestrować i wykorzystywać do przewidywania wielkości pęknięć, aby zapewnić konserwację lub wycofanie z eksploatacji przed uszkodzeniem któregokolwiek z pęknięć.
Trwałość zmęczeniową można podzielić na okres inicjacji i okres wzrostu pęknięć. Równania wzrostu pęknięć są wykorzystywane do przewidywania wielkości pęknięć począwszy od danej początkowej skazy i zazwyczaj opierają się na danych doświadczalnych uzyskanych z testów zmęczeniowych przy stałej amplitudzie .
z najwcześniejszych równań wzrostu pęknięć opartych na zakresie intensywności naprężeń w cyklu obciążenia ( ) jest Parisa-Erdogana
gdzie pęknięcia i wzrost pęknięcia zmęczeniowego dla pojedynczego cyklu obciążenia . Opracowano różne równania wzrostu pęknięć podobne do równania Parisa-Erdogana, aby uwzględnić czynniki wpływające na tempo wzrostu pęknięć, takie jak współczynnik naprężeń, przeciążenia i skutki historii obciążeń.
Zakres intensywności stresu można obliczyć z maksymalnego i minimalnego natężenia stresu dla cyklu
Współczynnik geometrii jest używany do powiązania naprężenia pola dalekiego z intensywnością naprężenia wierzchołka pęknięcia za pomocą
- .
Istnieją standardowe odniesienia zawierające współczynniki geometrii dla wielu różnych konfiguracji.
Historia równań propagacji pęknięć
Na przestrzeni lat zaproponowano wiele równań propagacji pęknięć, aby poprawić dokładność prognoz i uwzględnić różnorodne efekty. Prace Heada, Frosta i Dugdale'a, McEvily'ego i Illga oraz Liu dotyczące zachowania pęknięć zmęczeniowych położyły podwaliny pod ten temat. Ogólną postać tych równań propagacji pęknięć można wyrazić jako
długość pęknięcia jest oznaczona przez zastosowanych cykli obciążenia jest określona przez , zakres naprężeń przez , a materiał parametry przez do . symetrycznych długość pęknięcia od linii symetrii jest definiowana jako całkowitej długości pęknięcia .
w postaci są prawdziwym równaniem różniczkowym nie modelują procesu wzrostu pęknięć w sposób ciągły przez cały cykl Jako takie, osobne algorytmy zliczania cykli lub identyfikacji, takie jak powszechnie stosowany algorytm zliczania opadów deszczu , są wymagane do określenia wartości maksymalnych i minimalnych w cyklu. Chociaż wykazano, że zliczanie opadów deszczu zostało opracowane z myślą o metodach naprężeń / odkształceń, które również działają w przypadku wzrostu pęknięć. Opracowano również niewielką liczbę rzeczywistych równań wzrostu pęknięć zmęczeniowych pochodnych.
Czynniki wpływające na tempo wzrostu pęknięć
reżimy
typowy wykres tempa wzrostu pęknięć w funkcji zmiennego natężenia naprężeń lub siły napędowej wierzchołka pęknięcia na skali logarytmicznej Zachowanie szybkości wzrostu pęknięć w odniesieniu do zmiennego natężenia naprężeń można wyjaśnić w różnych reżimach (patrz rysunek 1) w następujący sposób
Reżim A: Przy niskich szybkościach wzrostu zmiany w mikrostrukturze , średnim naprężeniu (lub stosunku obciążenia) i środowisku mają znaczący wpływ na szybkość propagacji pęknięć. Obserwuje się, że przy niskich współczynnikach obciążenia szybkość wzrostu jest najbardziej wrażliwa na mikrostrukturę, aw materiałach o niskiej wytrzymałości jest najbardziej wrażliwa na współczynnik obciążenia.
Reżim B: W średnim zakresie tempa wzrostu zmiany mikrostruktury, średniego naprężenia (lub współczynnika obciążenia), grubości i środowiska nie mają znaczącego wpływu na szybkość propagacji pęknięć.
Reżim C: Przy dużych szybkościach wzrostu propagacja pęknięć jest bardzo wrażliwa na zmiany mikrostruktury, średniego naprężenia (lub współczynnika obciążenia) i grubości. Skutki środowiskowe mają stosunkowo mniejszy wpływ.
Efekt współczynnika stresu
Cykle o wyższym współczynniku naprężeń mają zwiększone tempo wzrostu pęknięć. Efekt ten jest często wyjaśniany za pomocą zamykania pęknięć , która opisuje obserwację, że powierzchnie pęknięć mogą pozostawać w kontakcie ze sobą przy obciążeniach powyżej zera. Zmniejsza to efektywny zakres współczynnika intensywności naprężeń i tempo wzrostu pęknięć zmęczeniowych.
Efekty sekwencyjne
ZA Równanie podaje tempo wzrostu dla pojedynczego cyklu, ale gdy obciążenie nie ma stałej amplitudy, zmiany obciążenia mogą prowadzić do tymczasowych wzrostów lub spadków tempa wzrostu. Aby rozwiązać niektóre z tych przypadków, opracowano dodatkowe równania. Tempo wzrostu jest spowolnione, gdy w sekwencji ładowania wystąpi przeciążenie. Obciążenia te generują strefy plastyczne, które mogą opóźniać tempo wzrostu. Dwa godne uwagi równania do modelowania opóźnień występujących podczas narastania pęknięcia w obszarze przeciążenia to:
- Model Wheelera (1972)
- =
gdzie jest strefą plastyczną odpowiadającą i-temu cyklowi, który występuje po przeciążeniu, a jest odległością między pęknięciem oraz zasięg strefy plastycznej przy przeciążeniu.
- modelu Willenborga
Równania wzrostu pęknięć
Równanie progowe
Do przewidywania tempa wzrostu pęknięć w rejonie progowym wykorzystano następującą zależność
Równanie Paryża-Erdoğana
Do przewidywania tempa wzrostu pęknięć w reżimie pośrednim wykorzystuje się równanie Parisa-Erdoğana
Równanie Formana
W 1967 roku Forman zaproponował następującą zależność, aby uwzględnić zwiększone tempo wzrostu spowodowane współczynnikiem naprężeń i zbliżaniem się do odporności na K
Równanie McEvily'ego-Groegera
McEvily i Groeger zaproponowali następujący związek mocy, który uwzględnia skutki zarówno wysokich, jak i niskich wartości
- .
równanie NASGRO
Równanie NASGRO jest używane w programach wzrostu pęknięć AFGROW, FASTRAN i NASGRO. Jest to ogólne równanie, które obejmuje niższą szybkość wzrostu w pobliżu progu wzrostu zbliżającą się do , a także uwzględnienie średniego efektu naprężenia poprzez włączenie współczynnika naprężeń. . Równanie NASGRO jest
gdzie do fa , , , , i to współczynniki równania.
Równanie McClintocka
W 1967 roku McClintock opracował równanie na górną granicę wzrostu pęknięcia w oparciu o cykliczne przemieszczenie otwarcia wierzchołka pęknięcia.
gdzie jest naprężeniem płynięcia, i w zakresie
Równanie Walkera
Aby uwzględnić efekt współczynnika stresu, Walker zasugerował zmodyfikowaną postać równania Parisa-Erdogana
gdzie jest parametrem który reprezentuje wpływ stosunku naprężeń na szybkość wzrostu pęknięć zmęczeniowych. Zazwyczaj przyjmuje wartość około ale może wahać się między . Ogólnie przyjmuje się, że ściskająca część cyklu obciążenia nie ma wpływu na wzrost co daje Można to fizycznie wytłumaczyć, biorąc pod uwagę, że pęknięcie zamyka się przy zerowym obciążeniu i nie zachowuje się jak pęknięcie pod obciążeniem ściskającym. W bardzo plastycznych materiałach ściskające przyczynia się do wzrostu pęknięć zgodnie .
Równanie Elbera
poprzez wprowadzenie poziomu intensywności naprężenia , przy którym następuje kontakt. Poniżej tego poziomu nie ma ruchu na wierzchołku pęknięcia, a zatem nie ma wzrostu. Efekt ten został wykorzystany do wyjaśnienia efektu współczynnika naprężeń i zwiększonego tempa wzrostu obserwowanego przy krótkich pęknięciach. Równanie Elbera jest
Równanie materiałów ciągliwych i kruchych
Ogólną postać szybkości wzrostu pęknięć zmęczeniowych w materiałach ciągliwych i kruchych podaje wzór
gdzie i . Opierając się na różnych mechanizmach posuwania się pęknięć i osłaniania wierzchołków pęknięć w metalach, ceramice i materiałach międzymetalicznych , obserwuje się, że tempo wzrostu pęknięć zmęczeniowych w metalach jest istotnie zależne od terminu, w ceramice od K. , a związki międzymetaliczne mają prawie podobną zależność od i K warunki.
Prognozowanie życia zmęczeniowego
Programy komputerowe
Istnieje wiele programów komputerowych, które implementują równania wzrostu pęknięć, takie jak Nasgro , AFGROW i Fastran . Ponadto istnieją również programy realizujące probabilistyczne podejście do wzrostu pęknięć, które obliczają prawdopodobieństwo awarii w całym okresie eksploatacji elementu.
Programy wzrostu pęknięć zwiększają pęknięcie od początkowego rozmiaru wady, aż przekroczy odporność na pękanie materiału i zawiedzie. Ponieważ odporność na pękanie zależy od warunków brzegowych, odporność na pękanie może zmieniać się od płaskich odkształcenia dla półkolistego pęknięcia powierzchniowego do płaskich warunków naprężeń dla pęknięcia przelotowego. Odporność na kruche pękanie w płaskich warunkach naprężeń jest zwykle dwa razy większa niż w przypadku płaskich odkształceń. Jednak ze względu na szybkie tempo wzrostu pęknięcia pod koniec jego żywotności, wahania odporności na kruche pękanie nie wpływają znacząco na żywotność elementu.
Programy wzrostu pęknięć zazwyczaj zapewniają wybór:
- metody liczenia cykli w celu wyodrębnienia ekstremów cykli
- współczynniki geometrii, które dobierają kształt pęknięcia i przyłożone obciążenie
- równanie wzrostu pęknięć
- modele akceleracji/opóźnienia
- właściwości materiału, takie jak granica plastyczności i odporność na pękanie
Rozwiązanie analityczne
Współczynnik intensywności stresu jest określony przez
gdzie jest przyłożonym równomiernym naprężeniem rozciągającym działającym na próbkę w kierunku prostopadłym do płaszczyzny pęknięcia, i jest bezwymiarowym parametrem, który jest zależy od geometrii próbki. Zmienia się intensywność naprężeń
gdzie zakresem cyklicznej amplitudy
Zakładając, że początkowy rozmiar pęknięcia wynosi pęknięcia próbki można obliczyć za pomocą jak
Powyższe równanie potrzeby można je rozwiązać numerycznie.
Przypadek I
Dla zamknięcie pęknięcia ma znikomy wpływ na tempo wzrostu pęknięcia, a równanie Paris-Erdogan może być użyte do obliczenia trwałości zmęczeniowej próbki, zanim osiągnie ona krytyczny rozmiar pęknięcia jak
Model wzrostu pęknięć o stałej wartości 𝛽 i R = 0
W przypadku modelu wzrostu pęknięć Griffitha-Irwina lub pęknięcia środkowego o długości rysunku 2, mamy i jest niezależny od pęknięcia długość. Można również Zakładając, że powyższa całka upraszcza się do
całkując powyższe wyrażenie dla przypadków i całkowita liczba cykli obciążenia jest dana przez m
Teraz, aby krytyczny rozmiar był bardzo duży w porównaniu z początkowym rozmiarem pęknięcia da
Powyższe wyrażenia analityczne dla całkowitej liczby cykli obciążenia do złamania Y . W przypadkach, w których środkowego pękniętego (CCT), do obliczenia można zastosować całkowanie numeryczne. .
Przypadek II
Ponieważ zjawisko , aby obliczyć trwałość zmęczeniową próbki, zanim osiągnie ona krytyczny rozmiar pęknięcia za do jako
Obliczenia numeryczne
przydatny, gdy od . Początkowy rozmiar pęknięcia jest uważany za . Współczynnik intensywności naprężenia przy obecnym rozmiarze pęknięcia przy użyciu maksymalnego przyłożonego naprężenia as
-
. jest mniejsza niż odporność na pękanie pęknięcie nie osiągnęło rozmiaru krytycznego i symulacja jest kontynuowana z bieżącym rozmiarem pęknięcia, aby obliczyć intensywność naprężeń przemiennych jako
Teraz, zastępując współczynnik intensywności naprężeń w równaniu Parisa-Erdogana, przyrost wielkości pęknięcia się jako
gdzie jest rozmiar kroku cyklu Nowy rozmiar pęknięcia staje się
gdzie indeks do bieżącego kroku iteracji. Nowy rozmiar pęknięcia jest używany do obliczenia intensywności naprężenia przy maksymalnym zastosowanym naprężeniu dla następnej iteracji. Ten iteracyjny proces jest kontynuowany aż
Po spełnieniu tego kryterium zniszczenia symulacja jest zatrzymywana.
Schematyczne przedstawienie procesu prognozowania trwałości zmęczeniowej przedstawiono na rysunku 3.
Przykład
Współczynnik intensywności naprężeń w próbce SENT (patrz rysunek 4) pod wpływem wzrostu pęknięć zmęczeniowych jest określony wzorem
Do obliczeń brane są pod uwagę następujące parametry
- mm, mm, mm, = ,
MPa, , .
Krytyczną długość pęknięcia obliczyć, gdy Jak
Rozwiązując powyższe równanie, długość
Teraz, powołując się na równanie Paris-Erdogan daje
Poprzez numeryczne całkowanie powyższego wyrażenia, całkowitą liczbę cykli obciążenia do zniszczenia otrzymuje się jako .
Linki zewnętrzne
- Forman, RG; Shivakumar, V.; Kardynał ŚJ; Williams, Kalifornia; McKeighan, PC (2005). „Baza danych wzrostu pęknięć zmęczeniowych do analizy tolerancji uszkodzeń” (PDF) . FAA . Źródło 6 lipca 2019 r .
- Gallagher, JP; Giessler, FJ; Berens, AP; Engle, Jr, JM „Podręcznik projektowania odpornego na uszkodzenia USAF: Wytyczne dotyczące analizy i projektowania konstrukcji samolotów odpornych na uszkodzenia. Wersja B” . Zarchiwizowane od oryginału w dniu 9 lipca 2019 r . Źródło 9 lipca 2019 r .
- „Podręcznik oceny tolerancji uszkodzeń, tom I: Wprowadzenie, mechanika pękania, propagacja pęknięć zmęczeniowych” (PDF) . Federalna Administracja Lotnictwa. 1993 . Źródło 16 lipca 2019 r .