Hodnocení únavové odolnosti svařovaných konstrukcímechanika.fs.cvut.cz/content/files/DPZ/weld_fatigue_2014_dpz.pdf · poškození v rámci standardních únavových postprocesorů(ANSYS

Hodnocení únavové odolnosti svařovaných konstrukcí

Jurenka Josef, Ph.D.Odbor pružnosti a pevnosti

Ústav mechaniky, biomechaniky a mechatronikyFakulta strojní, ČVUT v Praze

[email protected]

TechSoft Engineering, spol. s [email protected]

Obsah

1) Svařované konstrukce a metody predikce životnosti

2) Využití MKP při predikci únavového poškození

a) Predikce únavové životnosti na základě nominálních napětí

b) Predikce únavové životnosti na základě hot-spot napětí

c) Predikce únavové životnosti na základě vrubových elastických napětí v koření/patě svaru

d) Predikce únavové životnosti na základě lokálních elastických napětí v koření/patě svaru

3) Aplikace programů ANSYS při predikci únavového poškození svařovaných

konstrukcí - ANSYS nCode DesignLife – SEAM WELD, SPOT WELDS

2

1) Svařované konstrukce a metody predikce životnosti

3

Zvláštnosti svařovaných konstrukcí

• Nehomogenity materiálu:

– Mikrotrhliny, dutiny, neprůvary , přeložky, výpaly.

• Zbytková napětí v povrchové vrstvě:

– V oblasti svaru dosahují až meze kluzu a prudce klesají v jeho okolí.

• Geometrie svaru:

– Poloměry v patě a u kořene svaru, sklon obrysu svaru.

– Tyto vlivy se zohledňují statistickým vyhodnocením únavových zkoušek, zkoušky na reálných vzorcích předpokládají, že uvedené vlivy jsou „zahrnuty a pokryty“ křivkami životnosti.

• Specifika modelování MKP:

– 2D, 3D modely bez svaru, se svarovými detaily.

4

Postupy predikce únavového poškození

• Posouzení podle standardů:

– ČSN 05 0120 - Výpočet svarových spojů strojních konstrukcí.

– ČSN 73 1401 - Navrhování ocelových konstrukcí.

– ČSN EN 1993-1 – Navrhování ocelových konstrukcí.

– Eurocode 3.

– IIW (International Institue of Welding) - Recommendations for fatigue design of welded joints and components (IIW document XIII-1965-03 / XV-1127-03, 2005 Update)

• Posouzení na základě speciálních metodik:

– Metodika “Volvo“ implementovaná v ANSYS nCode DesignLife – SEAM WELDS

– Metodika Rupp, Störzel and Grubisic - SPOTWELDS

– Metodika Radaj, Seeger, Olivier - ANSYS nCode DesignLife nebo ANSYS Fatigue module

– Metodika FEMFAT

• Posouzení na základě vlastních únavových zkoušek.

5


• Rozsah použití:

– Konstrukce svařené z válcovaných nebo tažených polotovarů z refiticko-perlitické nebo bainitické oceli s maximální mezí kluzu fy = 960 MPa a austenitických korozivzdorných ocelí a slitin Al.

– Nelze použít v oblasti nízkocyklové únavy.

– Nelze použít v případě působení korozního prostředí a zvýšených teplot v oblasti creepu.

6

Definice používaných vztažných napětí

• Přístupy k predikci životnosti jsou poplatné uvažované (vypočtené) napjatosti:

7

Oblast nominálních napětí(Nominal stress approach)

Extrapolovaná (tvarová) napětí (Structural stresses)

Hot-spot napětí

Vrubová napětí(Notch stresses)

Lokální elastická napětí(FIN K)


8

Popis namáhání materiálu Popis únavové životnosti Stanovení únavové životnosti

Nominální napětí v daném průřezu

S-N křivky vyjádřené pomocí nominálních napětí

Kumulace únavového poškození (Palmgren-Miner)

Tvarové napětí (hot-spot) voblasti paty svaru

S-N křivky vyjádřené pomocí hot-spot napětí

Efektivní vrubové napětí v oblasti paty (kořene) svaru

S-N křivky vyjádřené pomocí vrubových napětí

Lokální napětí elastické napětí, resp. faktor intenzity napětí

Odolnost proti šíření únavových trhlin vyjádřená pomocí materiálových charakteristik růstových zákonů

Predikce šíření únavových trhlin Sčítání inkrementálních přírůstků (Parisův zákon)

MKP modelování svařovaných konstrukcí

• Nosníkové modely: Vhodné pouze pro výpočty nominálních napětí, resp. silových účinků Nominální přístupy (příhradové konstrukce)

• Skořepinové modely: Vhodné pro výpočet hot-spot napětí + využití v rámci některých speciálních postprocesingových metodik (nCode)

– Bez modelování svarů (konstrukční uzel není namáhán ohybem)

– S modelování svarů:

• Modelování pouze konce svaru

• Modelování celého svaru

• Modelování svaru pomocí tuhých vazeb

• Modelování svaru zvýšením tuhosti příslušných skořepinových elementů

• Objemové modely: Vhodné pro výpočty lokálních napětí Přístupy přes vrubová napětí a lokální elastická napětí – FIN K

– Modely s hrubou výpočetní sítí je možné využít pro výpočet hot-spot napětí.

– Modely s jemnou výpočetní sítí a uvažováním efektivního zaoblení paty svaru je možné použít pro výpočet vrubových napětí – možnost hodnotit kumulaci únavového poškození v rámci standardních únavových postprocesorů (ANSYS Fatigue module, nCode, PragTic, apod.)

– Výpočet FIN K vyžaduje modelování trhlin – speciální metodiky a postupy.

9

2a) Predikce únavové životnosti na základě nominálních napětí

10

Predikce únavové životnosti na základě nominální napětí

• Výpočet pomocí analytických vztahů nebo pomocí MKP– Jmenovitá nominální napětí:

– Modifikovaná nominální napětí:

11

N


• Volba S-N křivek konstrukčním uzlům podle International Institue ofWelding FAT třídy S-N křivek:

12



13


• Každá S-N křivka je charakterizována únavovou životnosti detailu na bázi 2106 cyklů a danou hladinou rozkmitu napětí XX a odpovídá pravděpodobnosti porušení 2,3% křivka FATXX.

14

– V případě normálového namáhání svaru je sklon S-N křivek m = 3 a mez únavy je definována na bázi 1107 cyklů.

– V případě smykového namáhání svaru je sklon S-N křivek m = 5 a mez únavy odpovídá bázi 1108

cyklů.– V případě namáhání konstrukčního

uzlu vysokým počtem cyklů není definována mez únavy. Vodorovná část S-N křivky je nahrazena šikmou se sklonem m = 22.

• Kumulace únavového poškození je určena pomocí Palmgen-Minerovasumačního pravidla.

FAT56


• Predikce probíhá podle S-N křivek přiřazených jednotlivým třídám konstrukčních detailů FAT.

• S-N křivky zahrnují následující faktory typické pro svařované konstrukce:

– Koncentraci napětí od makrogeometrie konstrukčního uzlu.

– Koncentraci napětí od tvaru svaru a imperfekcí svarového spoje spojené s uvažovanou technologií svařování.

– Směr působení nominálních napětí.

– Zbytková napětí od svařování.

– Metalurgické podmínky.

– Parametry svařovacího procesu a finálních úprav.

– Třída kvality provedení svařování B podle ISO 5817.

15

V rámci S-N křivek jednotlivých tříd FAT jsou zahrnuty určité geometrické nepřesnosti. V případě významnějších geometrických úchylek je nominální napětí zvětšeno o faktor km,eff.


• MKP modely neobsahují svarové detaily (objemové, skořepinové a nosníkové modely):

16

Stanovení nominálních napětí

2b) Predikce únavové životnosti na základě (tvarových) hot-spot

napětí

17

Predikce únavové životnosti na základě hot-spot napětí

• Stanovení pomocí MKP výpočtu nebo experimentálně (tenzometricky):

18

Napětí podél hrany konstrukce (nezávisí na tloušťce desky t)

Napětí po ploše konstrukce (závisí na tloušťce desky t)b) a)

N


• Postup výpočtu hot-spot napětí na základě MKP analýz (IIW):

19

MKP výpočet

Velikost elementů

Typ a) – napětí závisí na tloušťce Typ b) - napětí nezávisí na tloušťce

a) 0,4t x tb) < 4 x 4 mm

Napětí v rohových uzlech Napětí v rohových uzlech

a) 0,4t x tNapětí v rohových uzlech

a) t x tb) 10 x 10 mm

Napětí v mid-side uzlech Napětí v mid-side uzlech

tths σσσ 0,14,0 67,067,1

ttths σσσσ 4,19,04,0 72,024,252,2

tths σσσ 5,15,0 50,050,1

mmmmmmhs σσσσ 1284 33

mmmmhs σσσ 155 50,050,1


• Postup výpočtu na základě experimentu (tenzometrického měření):

• Přepočet na napětí:

20

Experiment

Typ a) Typ b)

Dva tenzometry Tři tenzometry

Tři tenzometry

tths εεε 0,14,0 67,067,1

ttths εεεε 4,19,04,0 72,024,252,2

mmmmmmhs εεεε 1284 33

hshs εEσ



21

Stejná sada křivek FAT S-N křivek jako v případě nominálních napětí



22

Stejná sada křivek FAT S-N křivek jako v případě nominálních napětí


• MKP modely (ne)obsahují svarové detaily (objemové, skořepinové modely):

23

Pat

h-2



24

• Předpokládá se, že únavová trhlina bude iniciována právě na konci svaru.

[Niemi]

[Fayard]

M. Aygül: Fatigue Analysis of Welded Structures Using the Finite Element Method, Chalmers university oftechnology, Gothenburg, Sweden 2012



25

[Niemi]• Přístup umožňuje zahrnout jak

geometrii svaru, tak i jeho tuhost.• Pata svaru může být “provařena“.• Elementy by měly být druhého

řádu.• Tloušťka “svarových“ elementů by

měla odpovídat tloušťce svaru.[Fayard]• Doporučeny jsou lineární

elementy. • Kořen svaru není “provařen“.• Hot-spot napětí lze odečítat

přímo v těžišti přilehlých elementů.

[Niemi, Eriksson]• Tuhost styčníku má větší význam

než geometrický tvar svaru.• Pro aplikaci je nutné upravit MKP

síť ve styčníku podle obr.• U přeplátování je kolmé spojení

realizováno tuhou vazbou.


2c) Predikce únavové životnosti na základě vrubových elastických

napětí v kořeni/patě svaru

26

Predikce únavové životnosti na základě vrubových elastických napětí v kořeni/patě svaru

• Elastická (fiktivní) napětí v kořeni/patě svaru.

• Pro zahrnutí statistické povahy a rozptylu parametrů popisujících tvar svarů a chování materiálu byl navržen postup, na základě kterého je skutečný tvar svaru nahrazen efektivním. Pro konstrukční ocele a slitiny Al byl stanoven efektivní rádius zaoblení paty, resp. kořene svaru o velikosti R1 (platí pro konstrukce s

tloušťkou stěny větší než 5 mm). Vypočtená napětí jsou srovnávána s jedinou S-N křivkou FAT225 (IIW).

27


• Hypotéza průměrných napětí (the stress averaging approach) – Neuber.

28

Fricke W., "Guideline for Fatigue Assessment by Notch Stress

Analysis for Welded Structures", IIW Doc. XIII-2240r1-

08/XV-1289r1-08, 2008.

• Aplikace tohoto přístupu vyžaduje detailní modelování svarových spojů:

29


• Při predikci únavového poškození se předpokládá, že únavová trhlina je iniciována v patě nebo v kořeni svaru – v místě zaoblení.

• Navržené efektivní zaoblení paty, resp. kořene svaru R1 je použitelné pro konstrukce s tloušťkou stěny t 5 mm (t 5 mm R0,05; FAT630).

• Do výpočtu lze zahrnout další významné geometrické úchylky a imperfekce – pomocí MKP modelování.

• Efektivní napětí je třeba chápat

jako smluvní (modelové zaoblení

R1), nelze je tedy měřit

přímo na konstrukci.

30


N

• Pro velikost efektivního zaoblení R1 (konstrukční ocele a Al slitiny) a standardní kvalitu provedení svaru jsou doporučeny (IIW) S-N křivky třídy FAT225, resp. FAT71 pro ocele, resp. Al slitiny.

• V těchto S-N křivkách je zahrnut vliv zbytkového napětí.

31


• Další modifikaci únavových křivek (FAT225) lze provést s ohledem na následující faktory:

– Střední napětí (zbytkového napětí)

– Tloušťka stěny (S-N křivky jsou definovány pro tloušťku 25 mm)

– Finální úpravy svarových spojů (broušení; otryskávání)

– Zvýšené teploty (vliv teploty 100 až 600°C)

– Korozní prostředí

• Konkrétní podobu uvedených modifikací je možné najít v příslušných normativech (IIW).

32


• MKP modelování – submodeling (požadavek jemné MKP sítě):

33



34



35



• Kumulace únavového poškození se provádí podle S-N křivky FAT225 (IIW recommendation).

36

FAT225

Aplikace programů ANSYS při predikci únavového poškození svařovaných konstrukcí – Fatigue module

• Experimentální validace relevance S-N křivky FAT225:

37


• Experimentální validace relevance S-N křivky FAT225:

38


Rekapitulace

39

Výhody Nevýhody

Přístup přes nominální napětí• Jednoduché výpočty• Poměrně značné historické zkušenosti• Často používané a aplikované na poměrně rozsáhlý

soubor konstrukčních detailů• Jsou dostupná potřebná experimentální data a analytické

vztahy• Vhodné pro predikci únavového poškození kumulovaného

v patě, resp. v kořeni svaru

• Predikce závislá na kategorizaci konstrukčních detailů• Nelze zahrnout libovolné geometrické parametry a

imperfekce• Méně přesné v případech komplexních konstrukčních uzlů• Vliv tloušťky nelze zahrnout

Přístup přes hot-spot napětí• Je potřeba méně S-N křivek• Přijatelně přesné• Jednoduché výpočty• Zahrnutí geometrických parametrů svarových spojů• Aplikace trubkové konstrukce

• Výsledky jsou závislé na hustotě a provedení MKP sítě • Různé způsoby vyhodnocení napětí• Vliv tloušťky nelze zahrnout• Použitelné pouze v případech, kdy dochází k iniciaci

únavových trhlin v patě svaru

Přístup přes vrubová napětí• Jediná S-N křivka• Možno zahrnout geometrické parametry svařovaného

spoje a imperfekce• Predikce únavového poškození kumulovaného v patě,

resp. kořeni svaru.

• Podkladem predikce jsou výsledky MKP analýz• Výsledky jsou závislé na hustotě MKP sítě a efektivním

zaoblení paty svaru• Poměrně náročné na modelování – rozsáhlé MKP modely

2d) Predikce únavové životnosti na základě lokálních elastických

napětí v kořeni/patě svaru

40

Predikce únavové životnosti na základě lokálních elastických napětí v kořeni/patě svaru

• V blízkosti svaru (kořen/pata) je uvažován defekt typu trhlina – ostrý vrub.

• Posouzení životnosti je založeno na výpočtu faktoru intenzity napětí K.

• Výpočet FIN K na základě analyticko-empirických vztahů (např. IIW) nebo pomocí MKP (ANSYS, MSC, ABAQUS apod.).

41

• Příklady tabelovaných analyticko-empirických vztahů pro výpočet FIN K (IIW recommendation):

42

benkbenbenmemkmemmem MYσMYσaπK ,,


• Predikce únavového poškození je založena na sčítání inkrementálních přírůstků délky trhliny, které jsou vypočteny pomocí tzv. Parisova vztahu:

• Doporučené hodnoty parametrů růstových zákonů:

43

fina

mth

m aKCNN

aKKpokusKC

N

a

inia00 dΔ0

d

dΔΔ,Δ

d

d

Doporučené parametry predikčních modelů (IIW recommendation)

Materiál Ocele Al slitiny

C0 1,5810-11[m, MPam] 1,2710-9 [m, MPam]

m 3 3

Kth 6,0-4,56R > 2 2,0-1,5R > 0,7


Hodnocení únavové odolnosti svařovaných konstrukcímechanika.fs.cvut.cz/content/files/DPZ/weld_fatigue_2014_dpz.pdf · poškození v rámci standardních únavových postprocesorů(ANSYS

Documents