🩹 Mikropęknięcia (prawo Parisa)
Naprężenia cykliczne powoli powiększają mikropęknięcia w tempie da/dN = C·(ΔK)^m. Patrz, jak pęknięcie pełznie krętą ścieżką — a potem pęka katastrofalnie, gdy K_max osiąga odporność na pękanie K_IC.
O tej symulacji
Ta symulacja modeluje wzrost pęknięcia zmęczeniowego w metalowej próbce poddanej obciążeniu cyklicznemu, wykorzystując prawo Parisa da/dN = C·(ΔK)^m, standardowe równanie inżynierskie opisujące stabilną propagację pęknięć. W każdym cyklu zakres współczynnika intensywności naprężeń na wierzchołku pęknięcia ΔK = Y·Δσ·√(πa) jest przeliczany na podstawie bieżącej długości pęknięcia a, a pęknięcie postępuje o da/dN milimetrów na cykl. Wzrost trwa, aż szczytowa intensywność naprężeń K_max = Y·σ_max·√(πa) osiągnie odporność na pękanie materiału K_IC — wtedy próbka pęka katastroficznie w pojedynczej klatce, zamiast dalej wolno się rozwijać.
🔬 Co pokazuje
Prostokątną próbkę rozciąganą między σ_max a σ_min w każdym cyklu, z pęknięciem rosnącym od początkowej długości a₀ w głąb materiału. Wykres na żywo śledzi długość pęknięcia a względem liczby cykli N w skali logarytmicznej, zaznaczając krytyczną długość pęknięcia a_c, przy której K_max = K_IC, a widok 3D rozdziela próbkę i pokazuje powierzchnię przełomu po wystąpieniu zniszczenia.
🎮 Jak korzystać
Wybierz gotowy materiał z listy rozwijanej (stal A36, aluminium 7075, tytan Ti-6Al-4V, żeliwo) lub wprowadź własne wartości C, m i K_IC. Dostosuj suwaki początkowej długości pęknięcia a₀, σ_max, σ_min, szerokości próbki W oraz współczynnika geometrycznego Y, ustaw liczbę cykli na klatkę, by przyspieszyć lub spowolnić symulację, i użyj przycisków Pauza/Odtwórz oraz Reset, by sterować przebiegiem.
💡 Czy wiesz, że?
Ponieważ da/dN skaluje się z ΔK podniesionym do potęgi m, podwojenie przyłożonego zakresu naprężeń Δσ nie tylko podwaja szybkość wzrostu pęknięcia — dla typowej stali z m≈3 mnoży ją mniej więcej przez 2³ = 8, dlatego niewielkie zwiększenie obciążenia cyklicznego może drastycznie skrócić żywotność zmęczeniową elementu.
Często zadawane pytania
Czym jest prawo Parisa i co oznaczają C i m?
Prawo Parisa to empiryczne równanie, da/dN = C·(ΔK)^m, opisujące szybkość wzrostu pęknięcia zmęczeniowego na cykl obciążenia w reżimie stabilnej propagacji. C i m to stałe materiałowe wyznaczone z danych doświadczalnych: C ustala ogólną szybkość wzrostu, a m kontroluje, jak bardzo ta szybkość jest wrażliwa na zakres współczynnika intensywności naprężeń ΔK. W tym symulatorze C podaje się w jednostkach ×10⁻¹¹, a m mieści się zwykle w zakresie od około 2,9 dla tytanu do 4,5 dla żeliwa.
Jak obliczany jest współczynnik intensywności naprężeń?
Symulacja wykorzystuje standardowy wzór liniowo-sprężystej mechaniki pękania K = Y·σ·√(πa), gdzie σ to przyłożone naprężenie, a to bieżąca długość pęknięcia, a Y to bezwymiarowy współczynnik geometryczny uwzględniający kształt próbki (regulowany suwakiem współczynnika geometrycznego Y, domyślnie 1,12 dla pęknięcia krawędziowego). ΔK wykorzystuje zakres naprężeń Δσ = σ_max − σ_min, natomiast K_max wykorzystuje bezpośrednio σ_max; zniszczenie następuje, gdy tylko K_max osiągnie K_IC.
Co decyduje o momencie pęknięcia próbki?
Do zniszczenia dochodzi, gdy szczytowa intensywność naprężeń K_max = Y·σ_max·√(πa) osiąga odporność na pękanie materiału K_IC, podawaną w MPa·√m za pomocą suwaka K_IC. Wraz ze wzrostem pęknięcia a rośnie K_max, mimo że przyłożone naprężenia pozostają stałe, więc każda kombinacja materiału i obciążenia ma krytyczną długość pęknięcia a_c (pokazaną jako przerywana linia na wykresie), powyżej której pęknięcie staje się niestabilne, a próbka pęka nagle zamiast wolno się rozwijać.
Czy to pełna symulacja fizyczna pęknięcia, czy uproszczony model?
To uproszczony model inżynierski, a nie symulacja metodą elementów skończonych ani pełna symulacja fizyczna. Sama ścieżka pęknięcia jest losowo drżącą linią dla efektu wizualnego, a geometria próbki to podstawowy prostokątny pręt, a nie prawdziwa próbka do badań pękania. Fizyka wzrostu jest jednak prawdziwym prawem Parisa całkowanym numerycznie cykl po cyklu za pomocą rzeczywistego wzoru na intensywność naprężeń, więc trendy w długości pęknięcia, szybkości wzrostu i liczbie cykli do zniszczenia odzwierciedlają rzeczywiste zachowanie zmęczeniowe.
Dlaczego różne materiały pękają po tak różnej liczbie cykli?
Każde ustawienie (stal A36, aluminium 7075, tytan Ti-6Al-4V, żeliwo) ma własne wartości C, m i K_IC zaczerpnięte z typowych publikowanych danych zmęczeniowych, a szacowana statystyka liczby cykli do zniszczenia jest obliczana przez całkowanie prawa Parisa od bieżącej długości pęknięcia do długości krytycznej a_c. Materiały o niskim K_IC (jak żeliwo przy 20 MPa·√m) osiągają swoją krytyczną długość pęknięcia szybciej, a wysoki wykładnik Parisa m sprawia, że wzrost przyspiesza znacznie szybciej po rozpoczęciu pękania, więc niewielkie różnice we właściwościach materiału przekładają się na duże różnice w żywotności zmęczeniowej.