Integralność konstrukcji
Wykład Nr 4
Metoda naprężenia nominalnego
Wydział Inżynierii Mechanicznej i Robotyki
Katedra Wytrzymałości, Zmęczenia Materiałów i Konstrukcji
http://zwmik.imir.agh.edu.pl/dydaktyka/dla_studentow/imir/imir.html
4.1. NAPRĘŻENIA NOMINALNE (lub średnie) - S i NAPRĘŻENIA LOKALNE -
a) rozciąganie pręta pryz- matycznego: y = S;
b) zginanie pręta pryzma- tycznego:
y = S, gdy S < Re,
y max < S, gdy S >Re; c) rozciąganie elementu z karbem: y S
y max = kt S, gdy kt S Re;
y < kt S, gdy kt S >Re.
Rys 4.1. Przykłady rozkładu naprężeń nominalnych S i lokalnych y w przekrojach wzdłuż osi x.
4.2. WYKRES WÖHLERA (tzw. KRZYWA S-N) MATERIAŁU
Rys 4.2. Krzywa S-N dla gładkich próbek ze stali A517 przy zginaniu obrotowym, z naprężeniem średnim m = 0.
4.2. WYKRES WÖHLERA (tzw. KRZYWA S-N) MATERIAŁU
Pojęcia podstawowe:
a) Wytrzymałość zmęczeniowa trwała materiału największa amplituda naprężenia a przy której nie dochodzi do zniszczenia próbki.
Wytrzymałość zmęczeniową trwałą wyznacza się ją z krzywej S - N dla próbek gładkich, jako:
asymptotę Z = a, przy Nf (stale zwykłej jakości i niskostopowe)
W tym przypadku jest to największa amplituda naprężenia, przy której nie nastąpi zniszczenie zmęczeniowe próbki.
wartość Z = a przy Nf = 107 lub 108 cykli, gdy brak asymptoty (np. stopy Al, Cu)
Wytrzymałość zmęczeniowa trwała jest stałą materiałową, ale zależy od sposobu obciążenia, np. przy zginaniu jest o 10-15 % wyższa niż przy rozciąganiu.
Stale: rozciąganie przy R = -1 Z 0.5 Rm (wartość niższa w stalach o wysokiej wytrzymałości)
4.2. WYKRES WÖHLERA (tzw. KRZYWA S-N) MATERIAŁU
Pojęcia podstawowe:
b) Wytrzymałość zmęczeniowa ograniczona największa amplituda naprężenia a, przy której nie nastąpi zniszczenie próbki przed upływem danej liczbie cykli Nf (np. Nf =105).
c) Zmęczenie wysokocyklowe naprężenia są na tyle niskie ze można pominąć odkształcenia plastyczne
d) Zmęczenie niskocyklowe typowo w zakresie 102-104 cykli, znaczne odkształcenia plastyczne.
Czynniki wpływające na wytrzymałość zmęczeniową:
obecność karbu,
naprężenia średnie m,
środowisko ,
Mikrostruktura,
naprężenia resztkowe (w związku z wpływem naprężenia średniego cyklu m).
4.3. MATEMATYCZNY OPIS KRZYWEJ S - N MATERIAŁU
Jeżeli krzywa Wöhlera (S-N) może być we współrzędnych podwójnie logarytmicznych aproksymowana linią prostą, to do jej opisu używa się zależności
a versus Nf w formie:a) równania:
a= A N
f B(4.1 a)
Rys. 4.3a Ilustracja opisu matematycznego krzywej Wöhlera wg równania (4.1a)
4.3. MATEMATYCZNY OPIS KRZYWEJ S - N MATERIAŁU
Jeżeli krzywa Wöhlera (S-N) może być we współrzędnych podwójnie logarytmicznych aproksymowana linią prostą, to do jej opisu używa się zależności
a versus Nf w formie:a) równania:
a= A N
f B(4.1 a) b) równania Basquina:
a=
f’ (2N
f)
b (4.1 b)Rys. 4.3b Ilustracja opisu matematycznego krzywej Wöhlera wg równania (4.1b)
4.3. MATEMATYCZNY OPIS KRZYWEJ S - N MATERIAŁU
(4.1 a) (4.1 b)
Stałe materiałowe A, B lub f’, b wyznacza się z dopasowania do równania (4.1 a) lub (4.1 b) danych z badań na próbkach gładkich.
Przy dużych odkształceniach plastycznych należy używać naprężenia rzeczywistego
Ponieważ 2Nf jest liczbą nawrotów obciążenia (1 cykl=2 nawroty), to f’ można interpretować jako wartość a, przy której następuje zniszczenie próbki po jednym nawrocie (półcyklu), tj. przy 2Nf = 1 (Nf = 0.5).
~
a4.3. MATEMATYCZNY OPIS KRZYWEJ S - N MATERIAŁU
(4.1 b)
Komentarz do równania Basquina (4.1b):
Gdyby własności materiału przy obciążeniu cyklicznie zmiennym były takie, jak przy obciążeniu monotonicznym, to naprężenie
f’
byłoby równe rzeczywistemu naprężeniu niszczącemu ( ), - por. rys. 2.4 i rów. (2.10) - gdyż próbę monotonicznego rozciągania można traktować jako jeden nawrót obciążenia zmęczeniowego. Jednak f’ różni się nieco od , gdyż
f’
wyznacza się przez ekstrapolację do Nf = 0.5 prostej dopasowanej do punktów (
a, Nf) otrzymanych z badań zmęczeniowych, gdy wartości materiału uległy zmianie na skutek cyklicznego umocnienia lub osłabienia (por. p. 3.3). Podobnie jak , naprężenie
f’
jest zawsze wyższe od niszczącego naprężenia inżynierskiego i od Rm , przy czym różnica ta jest mniejsza dla metali o wyższej wytrzymałości, które wykazują małe odkształcenia plastyczne. Wartości b dla różnych metali są na ogół zbliżone. ~
f ~
f ~
f4.3. MATEMATYCZNY OPIS KRZYWEJ S - N MATERIAŁU
(4.1 a)
(4.1 b)
a = A Nf B
a = f’ (2Nf)b Tabela 4.1 Parametry materiałowe występujące w równaniach (4.1a) i (4.1b)
Materiał Re Rm a = f’(2Nf)b=ANfB
f’ A b=B MPa MPa MPa MPa ---
stale
AISI 1015
normalizowana 227 415 976 886 -0.14 Man-Ten
walcowana na gorąco 322 557 1089 1006 -0.115 RQC-100
hart. i odpuszczana 683 758 938 897 -0.0648 AISI 4142
hart. i odpuszczana 1584 1757 1937 1837 -0.0762 AISI 4340
lotnicza 1103 1172 1758 1643 -0.0977 metale
nieżelazne
Al 2024-T4 303 476 900 839 -0.102 Ti-6Al-4V
przesycony i starzony 1185 1233 2030 1889 -0.104
4.3. MATEMATYCZNY OPIS KRZYWEJ S - N MATERIAŁU
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝑹𝒎 = 𝟏 𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 + 𝒌𝒇𝒎𝝈𝒎 𝑹𝒎
𝟐
= 𝟏
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝝈𝒇′ = 𝟏
ai
mii=1 n
RainFlow
= f(t) t
=
zrReguła P-M 𝑩𝒇 𝑵𝒊
𝑵𝒇𝒊 = 𝟏
B
fk
fi, k
fmikt, , , 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏 𝟏 + 𝜶 𝝆 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏
𝟏 + 𝜷 𝝆
Ni
ai ,mi
ariN
fi
arNf R=-1
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.1. Prezentacja wyników badań zmęczeniowych materiału (tzn. na próbkach gładkich) przy niezerowych naprężeniach średnich (
m = Sm)Gdy
m
0 to wyniki badań zmęczeniowych materiału przedstawia się według jednej z poniższych trzech koncepcji.a) R = const
Rys. 4.4 Krzywe S-N materiału przy stałym współczynniku asymetrii cyklu ( R = const.) Gdyby prezentowane tu wyniki przedstawiać jako dane a vs Nf, to najwyżej leżałaby krzywa R=-1 a najniżej krzywa R=0.
Np. dla Nf =104:
R 0 -0.5 -1
a (MPa) 410 530 5704.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.1. Prezentacja wyników badań zmęczeniowych materiału (tzn. na próbkach gładkich) przy niezerowych naprężeniach średnich (
m = Sm)Gdy
m
0 to wyniki badań zmęczeniowych materiału przedstawia się według jednej z poniższych trzech koncepcji.b) m = const
Rys. 4.5 Krzywe S-N materiału przy stałym naprężeniu średnim m = const)
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.1. Prezentacja wyników badań zmęczeniowych materiału (tzn. na próbkach gładkich) przy niezerowych naprężeniach średnich (
m = Sm)Gdy
m
0 to wyniki badań zmęczeniowych materiału przedstawia się według jednej z poniższych trzech koncepcji.c) Nf = const
Rys. 4.6 Wykresy stałej wartości (Nf=const) Uwaga: wykresy Nf=const na rys. 4.6 otrzymano z wykresów m = const z rys. 4.5 (por. te same oznaczenia punktów na obu rysunkach).
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.2. Znormalizowany wykres a /ar
Jeżeli każdą z krzywych Nf=const (rys. 4.6) przedstawi się w formie znormalizowanego wykresu a/ar versus m, gdzie ar - wytrzymałość zmęczeniowa przy m = 0 (R = -1) dla danego Nf, to wszystkie takie wykresy mają następujące dwa wspólne punkty:
(a/ar = 1; m = 0) oraz
(a/ar = 0; m = Rm).
Rys. 4.7 wskazuje, że występuje tendencja do konsolidacji punktów (a/ar; m) dla różnych Nf w pojedynczą krzywą.
Rys. 4.7 Znormalizowany wykres amplitudy w funkcji naprężenia średniego otrzymany z wykresów na rys. 4.5
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.3. Matematyczny opis zależności a/ar versus m Aproksymacja linii a/ar versus m:
a) równanie Goodmana (prosta): (4.2)
b) Równanie Gerbera (parabola): (4.3)
c) Równanie Morrowa (prosta): (4.4)
1
m m ar
a
R
0 przy
, 1
2
mm m ar
a
R
1
f m ar
a
f’ amplituda niszcząca po 1 nawrocie obciążenia (2Nf = 1), por. równanie (4.1b) i rys. 4.3b
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.3. Matematyczny opis zależności a/ar versus m
a) równanie Goodmana (prosta): (4.2)
b) Równanie Gerbera (parabola): (4.3)
c) Równanie Morrowa (prosta): (4.4)
1
m m ar
a
R
0 przy
, 1
2
m
m m ar
a
R
1
f m ar
a
Równanie (4.2) - najlepsze wyniki dla materiałów o niskiej ciągliwości.
Równanie (4.3) - najlepsze wyniki dla materiałów o wysokiej ciągliwości (wydłużenie procentowe w próbie rozciągania > 5 %, por p. 2.1). Przewiduje ono, niezgodnie z doświadczeniami, niekorzystny wpływ m<0 na wytrzymałość zmęczeniową. Założenie zachowawcze: przy m 0 - linia punktowana pozioma.
Równanie (4.4) - lepsza zgodność z eksperymentem w porównaniu z (4.2). Dobra aproksymecja wyników dla wszystkich materiałów ciągliwych.
Metale kruche (żeliwo): równanie (4.2) prowadzi do wyników niezachowawczych (punkty doświadczalne leżą pod prostą Goodmana). Stosuje się do nich specjalne równania.
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.4. Wyznaczenie trwałości przy niezerowym naprężeniu średnim m Podstawowa idea:
Dla danego materiału (scharakteryzowanego przez Rm lub f) trwałość zmęczeniowa przy dowolnej kombinacji amplitudy a i niezerowego naprężenia średniego m jest taka sama, jak przy amplitudzie ar i m=0.
Takie podejście jest dogodne, gdy dysponujemy tylko krzywą Wöhlera dla
m= 0, a chcemy wyznaczyć trwałość N
f(lub wytrzymałość zmęczeniową
a) przy
m 0.
Wtedy:
N
f(
a,
m0) = N
f(
ar,
m=0)
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.4. Wyznaczenie trwałości przy niezerowym naprężeniu średnim m
N
f(
a,
m0) = N
f(
ar,
m=0)
Z równania Goodmana (4.2) można wyznaczyć ar jako: (4.5)
Trwałość przy (a, m0) można wyznaczyć podstawiając do równania Basquina (4.1b):
prawą stronę równania (4.5) zamiast ar, otrzymując:
(4.6)
Z równania Morrowa (4.4) mamy: (4.7)
Uwzględniając (4.7) i równanie Basquina (4.1b) otrzymamy zależność:
(4.8) również określaną jako: równanie Morrowa
m m a ar
R
1
ar=
f (2N
f)
b
f b fm m
a
N
R 2
1
f m a ar
1
a= (
f -
m) (2N
f)
b4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
4.4.4. Wyznaczenie trwałości przy niezerowym naprężeniu średnim m Np. przy tej samej amplitudzie a
m/Rm ar/a (wg. 4.5)
0.2 1.25
0.5 2
m a m
ar
R
1
1
(4.5)
4.4. WPŁYW NAPRĘŻEŃ ŚREDNICH (
mlub S
m)
i=1 n
RainFlow
= f(t) t
=
zrReguła P-M 𝑩𝒇 𝑵𝒊
𝑵𝒇𝒊 = 𝟏
B
fk
fi, k
fmikt, , , 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏 𝟏 + 𝜶 𝝆 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏
𝟏 + 𝜷 𝝆
Ni
ai ,mi
ariN
fi
arNf R=-1
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝑹𝒎 = 𝟏 𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 + 𝒌𝒇𝒎𝝈𝒎 𝑹𝒎
𝟐
= 𝟏
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝝈𝒇′ = 𝟏
ai
mi4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
Rzeczywiste przebiegi obciążeń w czasie (tzw. historie obciążenie - czas) spotykane w warunkach eksploatacyjnych mają zazwyczaj charakter zmiennoamplitudowy.
Rys. 4.8 Siła w lewym kulistym przegubie zawieszenia samochodu w czasie przejazdu przez tory kolejowe
Przykłady :
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
Rys. 4.9 Maksymalne naprężenia zginające w połączeniu skrzydła z kadłubem w czasie jednego lotu samolotu o nieruchomych skrzydłach; a ) historia rzeczywista, b ) historia uproszczona
Przykłady obciążeń eksploatacyjnych:
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
Rys. 4.10 Zapis naprężeń w drążku kierowniczym samochodu: a) rzeczywista historia obciążenia; b) fragment historii obciążenia w czasie jazdy po nierównościach; c) obciążenie w czasie manewrowania
Przykłady obciążeń eksploatacyjnych:
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.1. Reguła Palmgrena - Minera
Założenia: Jeżeli amplituda
a,ipowtarza się przez N
icykli, a liczba cykli do zniszczenia określona z krzywej S-N przy tej amplitudzie wynosi N
f,i, to część trwałości zużytej przy
a,iwynosi N
i/N
f,i. Zniszczenie nastąpi, gdy:
1
,
i f
i
N
N
(4.9a)tzn. trwałość przewidywana wynosi:
N f , P M N i
(4.9b)4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.1. Reguła Palmgrena - Minera
1
,
i f
i
N
N
(4.9a)N f , P M N i
(4.9b)Rys. 4.10 Schemat objaśniający wykorzystanie reguły P - M do przewidywania trwałości materiału przy zmiennych amplitudach naprężeń dla przypadku: m = 0 (R = -1)
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.1. Reguła Palmgrena - Minera
1
,
i f
i
N
N
(4.9a)N f , P M N i
(4.9b)Jeżeli jedna i ta sama sekwencja obciążenia, którą wtedy można nazwać okresem, jest powtarzana wiele razy, np. lot samolotu, to:
1
, 1
okres i
f i
f
N
B N
gdzie: Bf - liczba powtórzeń okresu (4.10)okres i
f i
N N
,
1
uszkodzenie zmęczeniowe w 1 okresieJeżeli w jakichś cyklach historii obciążenie czas występują niezerowe naprężenia średnie, to Nf,i trzeba wyznaczyć np. z równań (4.6) lub (4.8).
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.2. Efekty interakcji obciążeń
Zjawisko to polega na tym, że w zmiennoamplitudowej historii obciążenia uszkodzenie zmęczeniowe Di spowodowane danym cyklem i (a,i, m,i) może być inne, niż przy obciążeniu stałoamplitudowym, tzn.:
i f
i
N
D
,
1
(4.11) gdzie: Nf, i trwałość przy obciążeniu stałoamplitudowym o
parametrach a,i, m,i
W zależności od historii obciążenia (spektrum obciążenia), materiału, poziomu średniego naprężenia spektrum i geometrii elementu może być:
i f
i
N
D
,
1
i f
i
N
D
,
1
niekorzystny efekt interakcji
korzystny efekt interakcji
(4.12a)
(4.12b)
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.2. Efekty interakcji obciążeń
Ponieważ reguła Palmgrena - Minera nie uwzględnia efektu interakcji obciążeń, w bardzo wielu przypadkach może dawać wyniki wysoce niezgodne z doświadczeniem, zarówno nadmiernie zachowawcze, jak i niezachowawcze. Może być:
1
100 N
100
N rzeczywiste
f P M f
,
Sposoby uwzględniania efektu interakcji obciążeń:
1) nieliniowe reguły kumulacji uszkodzeń 2) względna reguła P M
3) uwzględnienie amplitud poniżej trwałej wytrzymałości zmęczeniowej
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
Ad. 2 Względna reguła Palmgrena-Minera (Schütz, 1972)
Założenie: jeżeli dwie historie obciążenia są dostatecznie podobne, to odchylenia od reguły P - M mają te same kierunki i względne wartości.
Jeżeli dla jednego spektrum znamy gdzie to odpowiednio trwałości rzeczywiste i obliczone z reguły P-M, to dla drugiego spektrum które jest
„podobne” będzie:
obl f eksp
f
N
N N
f eksp, N
f oblobl obl
"
f eksp"
ff eksp
f
N N N
N
a stąd:
obl
"
obl"
f ekspN
fN
f ekspN
fN
(4.13)Wada: brak ogólnego kryterium „podobieństwa” spektrum.
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
Ad. 2 Względna reguła Palmgrena-Minera (Schütz, 1972) Praktyczne zastosowanie:
a) historia eksploatacyjna inna niż projektowa (zmiana zadań urządzenia, inne niż przewidziano warunki eksploatacji),
b) spektrum eksploatacyjne nie zostało ocenione prawidłowo,
c) nie jest możliwe przeprowadzenie w laboratorium badań symulujących pełną historię obciążenia w eksploatacji, np.: w przypadku spektrum obciążenia o długim „ogonie”
małych amplitud ze względów czasowych trzeba pominąć znaczną liczbę „małych” cykli
Rys. 4.11 Ilustracja konieczności pominięcia „małych” cykli w badaniach laboratoryjnych
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
Rys. 4.11 Ilustracja konieczności pominięcia „małych” cykli w badaniach laboratoryjnych
N - liczba przekroczeń danego poziomu amplitudy a
liczba cykli uwzględniona w badaniach laboratoryjnych: 107 liczba cykli przewidywana w eksploatacji: 109
liczba cykli pominiętych w badaniach laboratoryjnych Npom = 109-107 cykli = 9.9x108 cykli Zysk na czasie badań przy założeniu częstości obciążenia 20 Hz:
109 cykli = 578 dni; 108 cykli = 58 dni; 107 cykli = 6 dni
Widma lotnicze: pominięcie cykli o amplitudach poniżej 0.5Z - wzrost trwałości o 10 - 30 %.
„Małe” cykle w realistycznych, nieregularnych historiach obciążenia mogą się okazać szkodliwe, gdy w materiale istnieją już mikrouszkodzenia zmęczeniowe (także pasma poślizgów) spowodowane przez poprzedzające cykle.
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
Ad. 3 uwzględnienie amplitud poniżej trwałej wytrzymałości zmęczeniowej
1 3 2
10
7 a
Z
(log)N f
Rys. 4.12 Różne propozycje modyfikacji krzywej S-N do obliczeń trwałości przy obciążeniach zmienno- amplitudowych, 1 - obciążenie stałoamplitudowe (krzywe Wöhlera)
Modyfikacja krzywej S-N wg linii 2 lub 3.
Poprawa ocen trwałości przy zmiennych amplitudach przy użyciu linii 2 lub 3
jest możliwe tylko przy niekorzystnych efektach interakcji (por. równanie
4.12a). przy korzystnych efektach interakcji (por. równanie 4.12b) użycie linii
2 lub 3 spowoduje pogorszenie ocen N
f.
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
i=1 n
RainFlow
= f(t) t
=
zrk
fi, k
fmikt, , , 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏 𝟏 + 𝜶 𝝆 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏
𝟏 + 𝜷 𝝆
Ni
ai ,mi
ariN
fi
arNf R=-1
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝑹𝒎 = 𝟏 𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 + 𝒌𝒇𝒎𝝈𝒎 𝑹𝒎
𝟐
= 𝟏
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝝈𝒇′ = 𝟏
ai
miReguła P-M 𝑩𝒇 𝑵𝒊
𝑵𝒇𝒊 = 𝟏
B
f4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.3. Zliczanie cykli – metoda Rainflow
W przypadku nieregularnych historii obciążenia (por. np. rys. 4.8- 4.10) nie jest jasne, jakie wydarzenie uznać za cykl obciążenia. W licznych metodach liczenia cykli, które zaproponowano, wysunięto rozmaite propozycje. Obecnie za najbardziej racjonalne metody liczenia cykli uważa się techniki typu Rainflow (pierwsza propozycja - T. Endo, Japonia, 1968). W metodzie Rainflow zawsze uwzględnia się zakres między najwyższym maksimum i najniższym minimum.
Rys.4.13 Podstawowe wydarzenia obciążenia nieregularnego
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.3. Zliczanie cykli – metoda Rainflow
Rys. 4.13 Podstawowe wydarzenia obciążenia nieregularnego
Rys. 4.14 Warunek naliczania cyklu metodą Rainflow
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
4.5.3. Zliczanie cykli – metoda Rainflow
Rys. 4.15 Przykład naliczania cykli metodą Rainflow
4.5. OBCIĄŻENIA ZMIENNOAMPLITUDOWE
k
fi, k
fmikt, , , 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏 𝟏 + 𝜶 𝝆 𝒌𝒇= 𝟏 + 𝒌𝒕− 𝟏
𝟏 + 𝜷 𝝆
ai ,mi
ariN
fi
arNf R=-1
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝑹𝒎 = 𝟏 𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 + 𝒌𝒇𝒎𝝈𝒎 𝑹𝒎
𝟐
= 𝟏
𝒌𝒇𝝈𝒂
𝝈𝒂𝒓 +𝒌𝒇𝒎𝝈𝒎 𝝈𝒇′ = 𝟏
ai
miReguła P-M 𝑩𝒇 𝑵𝒊
𝑵𝒇𝒊 = 𝟏
B
fi=1 n
RainFlow
= f(t) t
=
zrNi
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
Pęknięcia zmęczeniowe i w rezultacie zniszczenie elementów konstrukcyjnych zostają z reguły zainicjowane w karbach (nieciągłości geometryczne, jak otwory, odsadzenia, rowki itp.). Przyczyna - spiętrzenie naprężeń spowodowane karbem, którego miarą jest współczynnik koncentracji naprężeń k
t(por p.6
„Przypomnienie” i rys.4.1c)
k
tzależy od: geometrii elementu, sposobu obciążenia
k
tnie zależy od: wielkości obciążenia, materiału, wielkości elementu
Uwaga: definicja naprężenia nominalnego S może się opierać na przekroju netto lub brutto, a jej wybór wpływa na wartość k
t.
W przykładzie z rys. 4.1c może więc być:
w d t S P
lubt w S P
Wartości kt można znaleźć w różnych poradnikach.
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
Rys. 4.15. Przykłady zmienności kt dla różnych karbów w zależności od geometrii
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.1. Wpływ karbu przy obciążeniach statycznych.
Rys.4.16 Element z karbem a) i rozkład naprężeń dla różnych przypadków: b) odkształcenie liniowo - sprężyste; c) lokalne płynięcie w materiale ciągliwym; d) płynięcie całego przekroju w materiale ciągliwym; e) naprężenie niszczące dla próbki z materiału kruchego
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.1. Wpływ karbu przy obciążeniach statycznych.
Rys.4.16
Materiały ciągliwe (rys. 4.16 b-d): stan naprężenia w przekroju karbu przed zniszczeniem (rys. 4.16d) jest taki, jak w próbce gładkiej o przekroju An. Stąd zniszczenie próbki z karbem, gdy:
S = naprężenie niszczące w próbce gładkiej o przekroju An, tj.:
S = Re (płynięcie przekroju netto), S = Rm (utrata spójności) Materiały kruche (rys.4.16e): utrata spójności, gdy:
max Rm czyli
t m
k
S R
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.2. Wpływ karbu przy obciążeniach zmęczeniowych - współczynnik działania karbu.
Gdyby o wytrzymałości zmęczeniowej decydowało naprężenie na dnie karbu,
to byłoby:
f t af
a
k
N S
N
(4.14) gdzie:
a(Nf) - wytrzymałość zmęczeniowa próbki gładkiej
Sa(Nf) - wytrzymałość zmęczeniowa próbki z karbem wyrażona w naprężeniach nominalnych
a i Sa - przy tej samej trwałości Nf Doświadczenie wskazuje, że:
f ta f
a
k
N S
N
(4.15) Współczynnik działania karbu kf (polskie oznaczenie k), definicja:
ar ar
f
S
k
(4.16) gdzie i S odnoszą się do R = - 1 i długiej trwałości (N = 106 107 cykli)
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.2. Wpływ karbu przy obciążeniach zmęczeniowych - współczynnik działania karbu.
Rys. 4.17 Wpływ karbu przy zginaniu obrotowym na krzywą S - N, stopu aluminium oraz porównanie wytrzymałości zmęczeniowej zredukowanej przy użyciu kt i kf
Wnioski z rys. 4.17:
krzywa - - - wg równania (4.14) leży pod eksperymentalną krzywą S-N próbki z karbem dla wszystkich trwałości
krzywa wg równania (4.16) leży pod eksperymentalną krzywą S-N próbki z karbem dla niskich trwałości. Oznacza to, że stosunek wytrzymałości zmęczeniowej próbki gładkiej do wytrzymałości zmęczeniowej próbki z karbem zależy od trwałości:
(4.17)
f f far f ar
f
f N k
N S
k N
) (
)
(
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.2. Wpływ karbu przy obciążeniach zmęczeniowych - współczynnik działania karbu.
Współczynnik wrażliwości na karb (definicja):
1 0
1
1
q
k q k
t f
(4.18)
Wartości graniczne q:
q=1, kf = kt (najwyższy możliwy wpływ karbu na wytrzymałość zmęczeniową)
q=0, kf =1 (karb nie wpływa na wytrzymałość zmęczeniową)
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.3. Przyczyny efektu kf < kt - interpretacja jakościowa
(4.19) 1) Gradient naprężeń w karbie (por. rys. 4.18)
Rozkład naprężeń w przekroju karbu y(x) przy założeniu materiału idealnie liniowo - sprężystego;
gradient naprężeń - miara spadku naprężeń ze wzrostem odległości x punktu od karbu
a) uszkodzenie zmęczeniowe w pewnej małej, skończonej objętości materiału
Rys. 4.18 Interpretacja wytrzymałości zmęczeniowej jako średniego naprężenia w skończonej odległości od wierzchołka karbu
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.3. Przyczyny efektu kf < kt - interpretacja jakościowa 1) Gradient naprężeń w karbie
a) uszkodzenie zmęczeniowe w pewnej małej, skończonej objętości materiału (4.19)
kf wg (4.19) będzie tym bardziej różnić się od kt im większy gradient naprężeń, a więc im mniejszy promień karbu . Trend zgodny z doświadczeniem, jak pokazuje rys. 4.19.
k średnia amplituda naprężenia między x i x
f
S
a
t0
k
Rys. 4.19 Współczynniki działania karbu dla różnych promieni karbu wyznaczone doświadczalnie z równania (4.16) dla stali miękkiej przy zginaniu obrotowym
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.3. Przyczyny efektu kf < kt - interpretacja jakościowa 1) Gradient naprężeń w karbie
b) teoria najsłabszego ogniwa
Przy ustalonej wartości max region wysokich naprężeń koniecznych do inicjacji uszkodzenia w miejscu defektu mikrostrukturalnego jest tym mniejszy, im wyższy gradient dy /dx.
Argument statystyczny: im mniejsza objętość materiału poddanego działaniu wysokich naprężeń, tym niższe prawdopodobieństwo, że znajdzie się tam defekt mikrostruktury, w którym nastąpi inicjacja pęknięcia (por. p. 1.2).
Stąd współczynnik kf będzie niższy przy większym gradiencie naprężeń, a więc mniejszym promieniu .
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.3. Przyczyny efektu kf < kt - interpretacja jakościowa 1) Gradient naprężeń w karbie
c) obecność pęknięcia (por. rys. 4.20)
Wierzchołek pęknięcia o długości l w próbce gładkiej znajduje się w strefie wyższych naprężeń, niż wierzchołek takiego samego pęknięcia w próbce z karbem.
Potwierdzenie: obecność tzw. pęknięć niepropagujących w próbkach z ostrymi karbami poddanych zmęczeniu wysokocyklowemu (Nf=106-107 cykli) przy amplitudach poniżej wytrzymałości zmęczeniowej.
Rys. 4.20 Próbka gładka i próbka z karbem przy tych samych naprężeniach lokalnych w miejscu zainicjowania pęknięcia (l = 0)
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.3. Przyczyny efektu kf < kt - interpretacja jakościowa 2) Odkształcenia plastyczne w karbie
Dotyczy zmęczenia niskocyklowego we wszystkich materiałach i zmęczenia wysokocyklowego w materiałach o bardzo wysokiej ciągliwości:
W strefie plastycznej karbu a< kt Sa, stąd musi być kf < kt
Rys. 4.21 Efekt odwróconego płynięcia w niewielkim obszarze w pobliżu karbu przy amplitudzie naprężeń S
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.4. Empiryczne oszacowanie kf
Najczęściej używane równania empiryczne:
a) Równanie Petersona:
k k
f
t
1 1
1
(4.20)
- stała materiałowa (zależna od sposobu obciążenia):zginanie, rozciąganie:
mm - stopy Al
mm - stale niskowęglowe normalizowane
mm - stale hartowane i temperowane
0 51 0 25
0 064 . .
.
wyżarzane lub
skręcanie: skr 0.6
Stale o podwyższonej i wysokiej wytrzymałości:
R 550 M Pa
2070 mm
025 .
0
m8 . 1
R
m MPa
(4.21)
- promień dna karbu4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.4. Empiryczne oszacowanie kf a) Równanie Petersona:
k k
f
t
1 1
1
(4.20)
2070 mm 025
. 0
8 . 1
Rm
MPa (4.21)dla Rm 550MPa
Rys. 4.22 Współczynnik wrażliwości na karb q (a) i wartości stałej
(b) dla stali wg równania Petersona (4.20).4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.4. Empiryczne oszacowanie kf b) Równanie Neubera:
(4.22)
(4.23)
Rys. 4.23 Współczynnik wrażliwości na karb q (a) i wartości stałej dla stali (b) wg równania Neubera
1
1
t1
f
k k
- stała materiałowa (zależna od sposobu obciążenia)
- promień dna karbuMPa mm R
m586
log 134
dla Rm 1520 MPa (stale)
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.4. Empiryczne oszacowanie kf
b) Równanie Neubera:
(4.22)
1
1
t1
f
k k
a) Równanie Petersona:
1 1
t1
f
k k
(4.20)Równania (4.20) i (4.22) nadają się do przybliżonego oszacowania kf dla karbów konstrukcyjnych (stosunkowo łagodnych).
Jeżeli karb jest głęboki i ostry, to lepszym podejściem jest mechanika pękania.
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.4. Empiryczne oszacowanie kf
Wnioski z rys. 4.22a) i 4.23a):
dla danego materiału: q rośnie z ;
dla danej klasy materiałów: q rośnie z Rm;
rozbieżność między kf i kt jest największa dla materiałów o dużej ciągliwości i ostrym karbie (por. też rys. 4.19).Rys. 4.22 a) Rys. 4.23 a)
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.5. Wpływ karbu przy średnich i krótkich trwałościach
Generalnie, stosunek wytrzymałości próbki gładkiej ar do wytrzymałości próbki z karbem Sar zależy od trwałości, por. równanie (4.17) i rys. 4.17:
f f far f ar
f
f N k
N S
k N
) (
)
(
(4.17)
Rys. 4.17 Wpływ karbu przy zginaniu obrotowym na krzywą S-N, stopu aluminium oraz porównanie wytrzymałości zmęczeniowej zredukowanej przy użyciu kt i kf.
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.5. Wpływ karbu przy średnich i krótkich trwałościach A) Metale o dużej ciągliwości:
Wpływ odwróconego płynięcia (por. rys. 4.21) jest tym większy, im wyższe naprężenia, a więc im niższa trwałość.
Stąd kf zmienia się od kf = kf (duże trwałości) do kf 1 (małe trwałości).
Rys. 4.24 Wyniki badań metalu ciągliwego ilustrujące zależność wpływu karbu od trwałości.
Punkty z wykresu S - N (rys. a) zostały użyte do otrzymania kf = a /Sa (rys. b)
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.5. Wpływ karbu przy średnich i krótkich trwałościach A) Metale o dużej ciągliwości:
Gdyby o wytrzymałości zmęczeniowej elementu konstrukcyjnego decydowała tylko amplituda naprężenia na dnie karbu a, to:
Rys. 4.25. Wyjaśnienie trendów widocznych na rys. 4.24 przy pomocy koncepcji odwróconego płynięcia dla materiału sprężysto - idealnie plastycznego.
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.5. Wpływ karbu przy średnich i krótkich trwałościach A) Metale o dużej ciągliwości:
Rys. 4.25a):
brak uplastycznienia (kt Sa Re), a,A=kt Sa, stąd kf = kt (4.24) por. zakres (a) wykresu na rys. 4.25d
Rys. 4.25d):
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.5. Wpływ karbu przy średnich i krótkich trwałościach A) Metale o dużej ciągliwości:
Rys. 4.25b)
odwrócone płynięcie tylko w otoczeniu karbu (kt Sa Re)
a,A=Re, stąd kf =Re/ Sa (4.25) por. zakres (b) wykresu na rys. 4.25d
Rys. 4.25d):
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.5. Wpływ karbu przy średnich i krótkich trwałościach A) Metale o dużej ciągliwości:
Rys. 4.25c):
odwrócone płynięcie w całym przekroju netto (Sa Re) – por. rys. 4.16d.
Jednorodny stan naprężenia w przekroju karbu, podobnie jak w próbce gładkiej
a=Sa, stąd kf 1 (4.26) por. zakres (c) wykresu na rys. 4.25d)
Rys. 4.25d):
4.6. OBLICZENIA TRWAŁOŚCI ELEMENTÓW
KONSTRUKCYJNYCH ZAWIERAJĄCYCH KARBY
4.6.5. Wpływ karbu przy średnich i krótkich trwałościach A) Metale o dużej ciągliwości:
Rys. 4.25 d)
Wniosek:
Linia kf (Sa) z rys. 4.25d) dobrze przybliża w sensie jakościowym trend w wartościach kf [Nf(Sa)] z rys. 4.24b). Różnica między poziomem kf =kf i wartością kt (przy długich trwałościach) wskazuje jednak na dodatkowy wpływ innych niż odwrócone płynięcie czynników na kf (por. p. 4.6.3).
Rys. 4.24 b)