Równanie wzrostu pęknięć

Rysunek 1: Typowy wykres szybkości wzrostu pęknięć w funkcji zakresu intensywności naprężeń. Równanie Paryża -Erdogana pasuje do centralnego regionu liniowego reżimu B.

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
Rysunek 2: Reprezentacja geometryczna próbki do badania rozciągania w środku pękniętego

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

Rysunek 3: Schematyczne przedstawienie procesu przewidywania trwałości zmęczeniowej

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

Rysunek 4: Reprezentacja geometryczna próbki do badania rozciągania z karbem pojedynczej krawędzi

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