WYZNACZANIE KRZYWEJ εa−Nf
Badania zmęczeniowe próbek gładkich przy Epsilona=const., Epsilonmin/Epsilonmax = - 1. Rejestruje się pętlę histerezy, której kształt stabilizuje się najpóźniej w połowie trwałości. Z każdego badania otrzymuje się punkt (Epsilona,Nf) lub parę punktów (Sigmaa/E,Nf) i (Epsilonap,Nf). Liczba cykli Nf odpowiada momentowi zniszczenia próbki lub oznacza liczbę cykli
po której amplituda naprężeń Sigmaa spadła do 50% wartości początkowej. $\varepsilon_{a} = \varepsilon_{\text{ae}} + \varepsilon_{\text{ap}} = \frac{\sigma_{a}}{E} + \varepsilon_{\text{ap}}$ (7.1) gdzie εap εae σaodczytuje się z ustabilizowanej pętli histerezy.Pkty wykresów (εae,Nf)(εap,Nf)układają się na krzywych, które można opisać równaniem: $\varepsilon_{\text{ae}} = \frac{\sigma_{a}}{E} = \frac{\sigma_{f}^{'}}{E}\left( 2N_{f} \right)^{b}$ (7.2a) εap = εf′(2Nf)c(7.2b) to jest równanie krzywej S-N rów. Basquina σa = σf′(2Nf)b(7.3) Z podstawienia (7.2) do (7.1) otrzymujemy
równanie Coffina – Mansona:$\varepsilon_{a} = \frac{\sigma_{f}^{'}}{E}\left( 2N_{f} \right)^{b} + \varepsilon_{f}^{'}\left( 2N_{f} \right)^{c}$ (7.4). Procedurę wyznaczania krzywejεa − Nf ilustruje rys. 7.2↓
znaczenie stałych materiałowych w
równaniu 7.4 objaśnia rys. 7.3.↓
UWAGI O STAŁYCH MATERIAŁOWYCH W RÓWNANIU
COFFINA- MANSONA (7.4) Duże trwałości:εap ≪ εae
krzywa εa(Nf) ≅ εae(Nf) Krótkie trwałości:εae ≪ εap
krzywa εa(Nf) ≅ εap(Nf) Przecięcie wykresówεae(Nf) i εap(Nf)
w pkt (εap=εae,Ntr) Ntr - tzw. trwałość przejściowa.Z (7.4) i z definicji Ntr mamy:$\frac{\sigma_{f}^{'}}{E} = \left( 2N_{\text{tr}} \right)^{b} = \varepsilon_{f}^{'}\left( 2N_{\text{tr}} \right)^{c}$ stąd $N_{\text{tr}} = \frac{1}{2}\left( \frac{\sigma_{f}^{'}}{\varepsilon_{f}^{'}E} \right)^{\frac{1}{c - b}}$ Zmęczenie niskocyklowe:Nf < Ntr (εap>εae) Zmęczenie wysokocyklowe: Nf > Ntr (εap>εae)
Z równania (7.2b) $2N_{f} = \left( \frac{\varepsilon_{\text{ap}}}{\varepsilon_{f}^{'}} \right)^{\frac{1}{c}}$ i (7.3) wynika $\sigma_{a} = \frac{\sigma_{f}^{'}}{\left( \varepsilon_{f}^{'} \right)^{\frac{b}{c}}}\varepsilon_{\text{ap}}^{\frac{b}{c}}$ (7.6) Z równania Ramberga - Ozgooda (3.5) cyklicznej krzywej σ − ε: σa = H′εapn′ (7.7) Z porównania (7.6) i (7.7) $n^{'} = \frac{b}{c}$ (7.8a) i $H^{'} = \frac{\sigma_{f}^{'}}{\left( \varepsilon_{f}^{'} \right)^{\frac{b}{c}}}$(7.8b)
Związki (7.8) wskazują, że z sześciu parametrów materiałowych n′, b, c, H′, εf′, σf′ tylko cztery są niezależne. Jednak zależności (7.8) są spełnione tylko w sposób przybliżony, gdyż H′ i n′ są wyznaczane z innych eksperymentów niż pozostałe cztery stałe materiałowe. Materiały ciągliwe: znaczna poprawa korelacji równań (7.2), (7.4), (7.7) z danymi
doświadczalnymi jeżeli użyje się $\tilde{\sigma_{a}}\ ,\ \tilde{\varepsilon_{a}}\ ,\ \tilde{\varepsilon_{\text{ap}}}\text{\ \ }$zamiast σa , εa , εap
Wówczas: $\sigma_{f}^{'} \cong \tilde{\sigma_{f}},\ \ \varepsilon_{f}^{'} \cong \tilde{\varepsilon_{f}}$
PRZYBLIŻONE OSZACOWANIE STAŁYCH MATERIAŁOWYCH
STOSOWANYCH W METODZIE ODKSZTAŁCENIA LOKALNEGO
Jeżeli używa się naprężeń i odkształceń rzeczywistych, to: $\sigma_{f}^{'} \cong \tilde{\sigma_{f}},\ \ \varepsilon_{f}^{'} \cong \tilde{\varepsilon_{f}}$ (7.9) Metale o wysokiej wytrzymałości (względnie kruche): (niska wytrzymałość ale duże odkształcenia plastyczne, duża energia do zniszczenia)
niskie wartości $\tilde{\sigma_{f}},\ \ $ i wysokie $\tilde{\varepsilon_{f}}$, wysoka trwałość Ntr, stroma krzywaεa − Nf rys. 7.4 linia 3. Metale o dużej udarności: (duże pole pod wykresem rozciągania, stosunkowo duża wytrzymałość i dobre własności plastyczne)
Zachowanie pośrednie - rys. 7.4 linia 2.Rys. 7.4 Tendencje obserwowane w krzywej εa − Nf(a) i pętli histerezy (b) dla
metali: (1) o wysokiej wytrzymałości, (2) o wysokiej udarności, (3) ciągliwych
UWZGLĘDNIENIE WPŁYWU NAPRĘŻEŃ ŚREDNICH
Podejście doświadczalne na podstawie badań zmęczeniowych
próbek gładkich.a) Badania pod kontrolą odkształcenia: εa = const i εm = const gdzie: εm ≠ 0 – odkształcenie średnie.
Cykliczna relaksacja naprężenia średniego: σm ≠ const rys. 7.6a. Wysokie
wartości εap − σm osiąga wartość zerową po niewielkiej liczbie cykli. Niskie
wartości εap− po pewnej liczbie cykli szybkość relaksacji dσm(N)/dN maleje
i σm ustala się na poziomie niższym od początkowego. Rys. 7.6
Schemat cyklicznej relaksacji naprężenia średniego (a) i cyklicznego pełzania (b)
b) Badania pod kontrolą naprężenia: σa = const i σm = const
gdzie:σm ≠ 0 – naprężenie średnie.
Przy εap ≠ 0 cykliczne pełzanie, tzn w kolejnych cyklach może rosnąć
odkształcenie średnie εm, rys. 7.6b. Jeżeli uplastycznienie występuje tylko
lokalnie (np. w obszarze karbu), to proces pełzania jest ograniczony przez
kontrakcję sprężystego materiału otaczającego strefę plastyczną.
Uwaga: Cykliczne pełzanie i relaksacja mogą występować równocześnie i
towarzyszyć cyklicznemu umocnieniu (osłabieniu) materiału. Zależą od
czasu, a więc i częstotliwości obciążenia.
Wzory empiryczne - gdy brak krzywej εa−Nfdla danego poziomu naprężenia średniego.
równania Ramberga - Osgooda (3.5): $\varepsilon_{a} = \frac{\sigma_{a}}{E} + \left( \frac{\sigma_{a}}{H^{'}} \right)^{\frac{1}{n^{'}}}\ $wzór Morrowa (4.8):σa = (σf′−σm)(2Nf)b (7.8a): $H^{'} = \sigma_{f}^{'}/{\varepsilon_{f}^{'}}^{\frac{b}{c}}\text{\ \ }$ (7.8b): n′ = b/c
Doświadczenia dowodzą, że zależność między σa i εanie zależy od poziomu σm. Równania (4.8) i (7.18) - to samo rozwiązanie Nf (w przybliżeniu, bo związki (7.8) są spełnione tylko w sposób przybliżony). Z tych 4-ech równań wynika: $\varepsilon_{a} = \frac{\sigma_{f}^{'}}{E}\left( 1 - \frac{\sigma_{m}}{\sigma_{f}^{'}} \right)\left( 2N_{f} \right)^{b} + \varepsilon_{f}^{'}\left( 1 - \frac{\sigma_{m}}{\sigma_{f}^{'}} \right)^{\frac{c}{b}}\left( 2N_{f} \right)^{c}$ (7.18) Przy wysokich wartościach εa (zakres krótkich trwałości) równanie (7.18) prowadzi do przesadnie zachowawczych (tzn. zbyt niskich) ocen Nf. Powód: intensywna relaksacja naprężeń średnich. Dlatego częściej stosuje się zmodyfikowaną formę (7.17), pomijając wpływ σm w drugim członie (tj. na εap): $\varepsilon_{a} = \frac{\sigma_{f}^{'}}{E}\left( 1 - \frac{\sigma_{m}}{\sigma_{f}^{'}} \right)\left( 2N_{f} \right)^{b} + \varepsilon_{f}^{'}\left( 2N_{f} \right)^{c}$ (7.19).
OKREŚLENIE TRWAŁOŚCI PRZY ZMĘCZENIU
ZMIENNOAMPLITUDOWYM W celu określenia liczby powtórzeń Bf
danej historii obciążenia aż do powstania pęknięcia (zniszczenia) w elemencie z karbem należy: 1) Dla każdego naliczonego metodą Rainflow cyklu „i” o amplitudzie lokalnego odkształcenia εa i lokalnym naprężeniu średnim σm , i określić trwałość Nf , i z (7.18) lub (7.19).
2) Wyznaczyć liczbę powtórzeń bloku do zniszczenia Bf na podstawie odpowiedniej hipotezy kumulacji uszkodzeń zmęczeniowych – np. reguły Palrgena-Minera.
PODSUMOWANIE METODY ODKSZTAŁCENIA LOKALNEGO
i PORÓWNANIE Z METODĄ NAPRĘŻENIA NOMINALNEGO
Metoda odkształcenia lokalnego powstała na przełomie lat 1950 i 1960, jako sposób analizy zmęczenia niskocyklowego (części reaktorów atomowych, silniki odrzutowe). Później zaczęto ją stosować do analizy zmęczenia zmiennoamplitudowego, gdy w historii obciążenia występują cykle o dużych amplitudach, powodujące lokalne odkształcenia plastyczne w karbach
elementów konstrukcyjnych (np. elementy samochodów i innych pojazdów i maszyn, turbiny i generatory, elementy aparatów latających i in.).
Sposób uwzględniania uplastycznienia jest racjonalniejszy niż w metodzie naprężenia nominalnego. Metoda odkształcenia lokalnego jest zalecana głównie do analizy zmęczenia niskocyklowego (Nf < Ntr). Przy wysokich trwałościach obie metody dają tożsame wyniki. Zalecane użycie metody naprężenia nominalnego: gdy dysponujemy krzywą S-N elementu
konstrukcyjnego z danego materiału bardzo podobnego do elementu właśnie rozważanego (np. połączenia spawane, nitowane lub sworzniowe, osie o specjalnej obróbce powierzchni itp.). Taka krzywa automatycznie ujmuje już wpływ różnych złożonych czynników (metalurgicznych, technologicznych, geometrii, frettingu i in.), trudnych do uwzględnienia w analizie metodą
odkształcenia lokalnego. Może to równoważyć brak szczegółowej analizy lokalnego odkształcenia.Uwaga: Wartości Nf otrzymane metodą odkształcenia lokalnego odnoszą się do zniszczenia lub powstania znacznych pęknięć w małych (5-10 mm średnicy) próbkach gładkich. Dla elementu z
karbem przewidujemy więc, że po „liczbie” cykli Ni = Nf powstanie widoczne gołym okiem pęknięcie (tzw. pęknięcie o wymiarach inżynierskich ai) o długości rzędu 1-5 mm. Jeżeli powstanie takiego pęknięcia nie oznacza definitywnego zniszczenia elementu, można obliczyć pozostającą jeszcze trwałość Np korzystając z mechaniki pękania. Jest to tzw. podejście
dwustopniowe. Nf=Ni+Np gdzie: Nf - całkowita trwałość elementu konstrukcyjnego; Ni - tzw. okres inicjacji pęknięcia - równy liczbie cykli do powstania pęknięcia o wymiarze ai obliczonej metodą odkształcenia lub naprężenia nominalnego, Np - okres propagacji pęknięcia - równy liczbie cykli po których pęknięcie o wymiarze ai osiągnie wymiar krytyczny akr
ZGINANIE SPRĘŻYSTO-PLASTYCZNE PRĘTA
PRYZMATYCZNEGO O PRZEKROJU PROSTOKĄTNYM PRZY
OBCIĄŻENIACH CYKLICZNYCH (σa=const)
Rys. 6.4. Belka o przekroju prostokątnym przy
cyklicznym zginaniu (a) powodującym
uplastycznienie przy obciążeniu i odciążeniu.
Rozkład cyklicznych naprężeń (b) i odkształceń (c).
Zmienność naprężeń w funkcji odkształceń we
włóknach skrajnych (d) oraz zależność między
poziomem odkształceń we włóknach skrajnych a
obciążeniem cyklicznym Mg (e).
Niech Mg zmienia się cyklicznie między Mmax i
Mmin (rys.6.4a). Równania równowagi na poziomach Mg
= Mmax i Mg = Mmin Mmax=2t∫0hσmaxydy Mmin=2t∫0hσminydy (6.12) M=Mmax−Mmin=2t∫0hσydy (6.13) gdzie σ = σmax(y) − σmin(y) ; σmax(σmin) – rozkład naprężeń σ(y) gdy Mg = Mmax(Mg = Mmin) Na każdym poziomie momentu Mg obowiązuje
prawo płaskich przekrojów rys. 6.4c $\frac{\mathbf{\varepsilon}}{\mathbf{y}}\mathbf{=}\frac{\mathbf{}\mathbf{\varepsilon}_{\mathbf{h}}}{\mathbf{h}}$ (6.14) Równania (6.13) i (6.14) są formalnie identyczne z (6.2) i (6.1): M = 2t∫0hσgydy (6.2); M = 2t∫0hσydy (6.13) ; $\frac{\varepsilon}{y} = \frac{\varepsilon_{h}}{h}$ (6.1); $\frac{\varepsilon}{y} = \frac{\varepsilon_{h}}{h}$ (6.14); Stąd naprężenia i odkształcenia przy odciążaniu od Mmax do Mmin można obliczyć jak przy monotonicznym obciążaniu od O do Mg = M, przy czym: I. Mg , σg , ε zastępujemy przez M, σ, ε II. posługujemy się dwukrotnie rozszerzoną krzywą ε = f(σ): $\frac{\varepsilon}{2} = f\left( \frac{\sigma}{2} \right)\text{\ tj.\ }\varepsilon_{a} = f(\sigma_{a})$ (6.15); III. obliczamy: εmin = εmax − ε = εmax − 2εa σmin = σmax − σ = σ max − 2σa (6.16) ; IV. przy ponownym obciążeniu stosujemy tę samą
procedurę, ale początek układu (σ,ε) jest w
punkcie (σmin, εmin). Gdy Mg po raz drugi osiąga poziom Mmax, wierzchołek pętli histerezy znajduje się ponownie w punkcie (σmax,εmax)rys.6.4d. Przy wyidealizowanym założeniu, że własności materiału nie ulegają zmianie przy obciążeniach cyklicznych, punkt (,) poruszałby się cały czas po
jednej i tej samej zamkniętej pętli histerezy. Jeżeli operujemy amplitudami, a nie zakresami to, równania monotonicznej i cyklicznej krzywej odkształcenia są formalnie identyczne: εMAX = f(σMAX) εa = f(σa) (6.17) ; Stąd zależności między εh i Mmax oraz między εah i Masą też formalnie identyczne: εh = g(MMAX) εah = g(Ma) (6.18)
PRZYBLIŻONE WYZNACZANIE NAPRĘŻEŃ I ODKSZTAŁCEŃ
W STREFIE PLASTYCZNEJ KARBU - REGUŁA NEUBERA
Rys.6.5. Ilustracja reguły Neubera: a) element z karbem i krzywa odkształcenia wraz z hiperbolą Neubera; b) zmienność współczynnika koncentracji naprężeń ksigma i odkształceń kepsilon w karbie; c), d) zmienność odpowiednio odkształceń i naprężeń w karbie w funkcji naprężenia nominalnego (linia
przerywana - rozwiązanie liniowo-sprężyste). Jeżeli materiał w strefie karbu uplastycznia się, to lokalne odkształcenia są większe, niż ktS/E(rys.6.5c), a lokalne naprężenia niższe niż kt=S (rys.6.5d). Należy, więc zdefiniować oddzielnie współczynniki koncentracji dla naprężeń i
odkształceń: kσ=σ/S kε=ε/e (6.20) gdzie: e - odkształcenia nominalne związane z naprężeniem nominalnym S w myśl równania krzywej odkształcenia materiału S=f(e). Jeżeli S≤Re (co zazwyczaj ma miejsce), to: e = S/E (6.21) Reguła Neubera: $\sqrt{\mathbf{k}_{\mathbf{\sigma}}\mathbf{k}_{\mathbf{\varepsilon}}}\mathbf{=}\mathbf{k}_{\mathbf{t}}$ (6.22). Jeżeli przy osiowym rozciąganiu nastąpiłoby pełne uplastycznienie przekroju, (S = Re) to kσ → 1
a więc w myśl (6.22), kε → kt2 co przedstawia rys. 6.5b. Uwzględniając w (6.22) zależności (6.20) oraz (6.21). mamy: $\mathbf{\sigma\varepsilon =}\frac{\left( \mathbf{k}_{\mathbf{t}}\mathbf{S} \right)^{\mathbf{2}}}{\mathbf{E}}$ (6.23) Jeżeli ε = f(σ) ma postać krzywej Ramberga-Osgooda $\mathbf{\varepsilon =}\frac{\mathbf{\sigma}}{\mathbf{E}}\mathbf{+}\left( \frac{\mathbf{\sigma}}{\mathbf{H}} \right)^{\frac{\mathbf{1}}{\mathbf{n}}}$ (6.24) to po podstawieniu za epsilon w (6.23) prawej strony (6.24) dostajemy: $\mathbf{\sigma}^{\mathbf{2}}\mathbf{+}\left( \frac{\mathbf{E}}{\mathbf{H}^{\frac{\mathbf{1}}{\mathbf{n}}}} \right)\mathbf{\sigma}^{\frac{\mathbf{1}}{\mathbf{n}}\mathbf{+ 1}}\mathbf{-}\left( \mathbf{k}_{\mathbf{t}}\mathbf{S} \right)\mathbf{= 0}$ (6.25) Uwaga: W praktyce w celu wyznaczenia lokalnych odkształceń i naprężeń przy użyciu reguły Neubera
najwygodniej jest korzystać z układu równań (6.23) i (6.25).
W celu wyznaczenia amplitud, zastępujemy w (6.23) i (6.25) σ, ε, i S przez odpowiednio σa, εa i Sa
)gma awieniu za epsilon w (6.23) prawej strony (6.24) dosatjemy : dddddddddddddddddddddddddddddddddddddddsadadadsadasdgafgdf