OBSAH ANOTACE .............................................................................................................................................2 1
ÚVOD ............................................................................................................................................3
2
ZADÁNÍ PŘÍKLADU .......................................................................................................................3
3
PODLOŽÍ ZATÍŽENÉ POMĚRNÝM VODOROVNÝM PŘETVOŘENÍM...............................................4
4
5
6
3.1
VÝPOČET PODLE ČSN 73 0039 ..................................................................................................4
3.2
VÝPOČET PODLE TEORIE PRUŽNOSTI ............................................................................................5
3.2.1
Výpočet programem NEXIS 32 ............................................................................................7
3.2.2
Výpočet programem ANSYS ................................................................................................7
3.2.3
Odchylky ve výpočtech programů NEXIS a ANSYS ...............................................................8
VÝMĚNA PODLOŽÍ ........................................................................................................................9 4.1
VÝPOČET PODLE ČSN 73 0039 ................................................................................................10
4.2
VÝPOČET PODLE TEORIE PRUŽNOSTI ..........................................................................................11
4.2.1
Výpočet programem NEXIS 32 ..........................................................................................11
4.2.2
Výpočet programem ANSYS ..............................................................................................12
APLIKACE REOLOGICKÉ KLUZNÉ SPÁRY.................................................................................12 5.1
ŘEŠENÍ PLOŠNÉ ......................................................................................................................12
5.2
VÝPOČET PODLE TEORIE PRUŽNOSTI ..........................................................................................13
5.2.1
Výpočet programem NEXIS 32 ..........................................................................................14
5.2.2
Výpočet programem ANSYS ..............................................................................................14
5.2.3
Výpočet programem NEXIS pro Ea = 0,001 MPa .................................................................18
5.2.4
Výpočet programem ANSYS Ea = 0,001 MPa .....................................................................18
ZÁVĚR.........................................................................................................................................18
SEZNAM LITERATURY: .......................................................................................................................20 PŘÍLOHY .............................................................................................................................................21 METODIKA VÝPOČTU SMYKOVÉHO NAPĚTÍ OD PŘETVOŘENÍ DLE ČSN 73 0039 METODIKA VÝPOČTU DLE ČSN 73 0039 – PODDAJNĚJŠÍ ZEMINA METODIKA VÝPOČTU PŘI APLIKACI KLUZNÉ SPÁRY PRŮBĚHY VNITŘNÍCH SIL SROVNÁNÍ ÚČINKU METOD SNIŽUJÍCÍ TŘENÍ V ZÁKLADOVÉ SPÁŘE VAN DER POELŮV NOMOGRAM
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
ŘEŠENÍ NAPJATOSTI ZÁKLADOVÝCH KONSTRUKCÍ NA ÚČINKY DEFORMAČNÍCH ZATÍŽENÍ Řešitel:
Petr Maňásek, VŠB – TU Ostrava, Fakulta stavební obor: 36 - 16 - 8 Průmyslové a pozemní stavitelství student V. ročníku
Vedoucí práce:
Doc. Ing. Radim Čajka, CSc. VŠB – TU Ostrava, Fakulta stavební
ANOTACE Práce se zabývá stanovením účinků poměrného vodorovného přetvoření podloží na základové konstrukce. Díky tomuto přetvoření dochází v základové spáře ke vzniku třecích sil, které mohou nabývat vysokých hodnot a tím podstatně ovlivnit celkový návrh základu. Stanovení těchto účinků lze v zásadě provést dvěma způsoby. První možností je použití normy ČSN 73 0039, jenž představuje řešení plošné. Druhé řešení je provedení výpočtu metodou konečných prvků (MKP), jenž představuje řešení podle teorie pružnosti. V práci je věnován prostor oběma těmto metodám. Jejich výsledky jsou porovnány a analyzovány. Následující část práce se zabývá možností snížení nežádoucích třecích sil v základové spáře. Pozornost je věnována aplikaci reologické kluzné spáry a výměně podloží za poddajnější zeminu bezprostředně pod základovou konstrukcí.
ANNOTATION The thesis is interested in a description of a effect of a rationed horizontal stress deformation of a foundation construction. This deformation leads to occurring of the friction force at the footing bottom of the construction foundation, which can reach of high values and involved the global design of the foundation. The first way, how to determinate the effect is a possibility to use ČSN 73 0039, which describes solution at the flat. The second way, how to determinate it, is to use the calculation by the Finite Element Method (FEM). This solution accords to the theory of elasticity. Both of this methods are at the work described and used. The results received by them are compared and analysed. The next part of the thesis is interested at the reduction of this friction forces. The focus is set to the application of a rheology shear gripe and a exchange of a subsoil by the flexible soil directly behind the foundation construction.
-2-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
1 ÚVOD Stanovení účinků na základové konstrukce od poměrného vodorovného přetvoření podloží lze provést v zásadě dvěma způsoby. První možností je využití postupu daného normou ČSN 73 0039 [1], jenž představuje řešení plošné. Druhou možnost jak stanovit tyto účinky je provést výpočet metodou konečných prvků (MKP), představující řešení založené na principech teorie pružnosti. Řešení podle teorie pružnosti je založeno na 2D modelu svislého řezu základové konstrukce včetně podloží základu. Pro výpočet MKP bylo použito výpočtových systémů NEXIS 32 a ANSYS. Účinky stanovené těmito metodami jsou následně porovnány a analyzovány. První fáze řešení se zabývá stavem, kdy se základová konstrukce nachází v přímém kontaktu s podložím, jenž je zatíženo poměrným vodorovným přetvořením terénu. Druhá fáze výpočtů je zaměřena na možnost snížení tření v základové spáře. Jedná se o aplikaci reologické kluzné spáry či výměny podloží bezprostředně pod základovou konstrukcí za poddajnější zeminu. Vzhledem k omezeným možnostem rozsahu práce jsou podrobné teoretické postupy jednotlivých metod uvedeny v příloze. Hlavní část je zaměřena na vlastní výpočty.
2 ZADÁNÍ PŘÍKLADU Schéma zadání, parametry základové desky a podloží jsou znázorněny na Obr. 1. Vnitřní síly a posuvy budou stanoveny na základovém pásu o délce l = 16 m, šířce b = 1 m a tloušťce h = 0,5 m. Velikost poměrného vodorovného přetvoření podloží uvažujeme εmax = 5·10-3. Základovou půdu tvoří ulehlý písek třídy S3 podle ČSN 73 1001 [2] s úhlem vnitřního tření ϕ = 32°, soudržností c = 0 kPa a modulem deformace Edef = 20 MPa při Poissonově čísle ν = 0,3. Provozní výpočtová hodnota rovnoměrného svislého napětí v kontaktní spáře je σv = 240 kPa. y
16.0 0.5
1. 0
x z
S3
+ε Obr. 1
+ε Schéma zadání
Plošné řešení dle ČSN 73 0039 je dáno přímo touto normou a je podrobně popsáno v příloze. V případě výpočtu MKP, jak již bylo uvedeno, je řešení založeno na 2D modelu svislého řezu. Z tohoto důvodu je nezbytné zvolit dostatečnou hloubku podloží, do které se mohou účinky přetvoření projevovat. Podobně je tomu i v případě přesahu podloží přes okraj základové konstrukce. V případě nedostatečných rozměrů podloží by došlo k negativnímu ovlivnění výstupů vnitřních sil. Rozměry podloží byly stanoveny experimentálně na základě několika výpočtů. Jako postačující se v tomto případě jeví přesah i hloubka cca 10,0 m. Vzhledem k symetrii základové desky (Obr. 1) jsou výpočty provedeny pouze pro pravou polovinu základové konstrukce. Levá polovina je samozřejmě zatížena shodně a tudíž i výsledné hodnoty musejí být rovněž naprosto stejné. -3-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
3 PODLOŽÍ ZATÍŽENÉ POMĚRNÝM VODOROVNÝM PŘETVOŘENÍM 3.1
Výpočet podle ČSN 73 0039
Nejprve je nutno stanovit hloubku tlumící vrstvy a, do které se projeví účinky působení vodorovného poměrného přetvoření terénu: 0 , 53 0 , 53 a = 0,75 ⋅ L0,56 1 − e −0,94⋅b = 0,75 ⋅ 16,0 0,56 1 − e −0,94⋅1,0
(
)
(
)
a = 2,16m
Podle ČSN 73 0039 [1] (tab. 4, příloha 1) určíme hodnotu součinitele η: σ v = 0,24MPa b = 1,0m ⇒ η = 1,19 a 2,16 = = 2,16 ≥ 1,5 b 1,0 Dále stanovíme korekční součinitel µε podle tab. 2 ČSN 73 0039 [1], uvažující délku základové konstrukce: 15m < L = 16m < 30m ⇒ µ ε = 0,85 Vlastní poměrné protažení základové konstrukce dle [1] uvažujeme ε eig = 1 ⋅ 10 −3 Oedometrický modul přetvárnosti Eoed základové půdy stanovíme pomocí vztahu dle ČSN 73 1001 [2]: 2ν 2 2 ⋅ 0,3 2 β = 1− = 1− = 0,743 1 −ν 1 − 0,3 1 1 E oed = ⋅ E def = ⋅ 20,0 = 27,0MPa β 0,743
Hodnota smykového napětí v základové spáře má tedy výsledný tvar: τ x = β x ⋅ µ ε (ε max − ε eig ) ⋅η ⋅ E oed = β x ⋅ 0,85(5 ⋅ 10 −3 − 1 ⋅ 10 −3 ) ⋅ 1,19 ⋅ 27 ⋅ 10 3
τ x = 109,242β x
Tabulka 1
Výsledné smykové napětí τx a posuv ∆x po délce základu
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
0,000
0,231
0,694
1,157
1,620
2,083
2,546
3,009
3,472
3,704
βx
0,000
0,001
0,031
0,100
0,192
0,295
0,404
0,515
0,628
0,685
τx [kPa] ∆x [mm]
0,000
0,163
3,388
10,944
21,012
32,276
44,120
56,267
68,584
74,784
0,0
1,7
5,1
8,5
11,9
15,3
18,7
22,1
25,5
27,2
ξ=
x a
Dále je nutno ověřit zda se smykové napětí τx vyvíjí po délce základu pouze v oblasti pružnoplastických poměrných přetvoření, tzn. zda dojde k překročení smykové únosnosti zeminy: -4-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
τ u = σ ⋅ tan ϕ + c = 240,0 ⋅ tan 32,0 o + 0,0 τ u = 150kPa ≥ τ max = 74,8kPa Nerovnost je splněna, smyková únosnost zeminy je větší a k prokluzu tedy nedojde. Na základě znalosti průběhu smykového napětí τ po délce základu stanovíme hodnotu smykové síly sumací tohoto napětí po dílčích úsecích délky Li, šířky bi a průměrné hodnoty smykového napětí τi od okraje základu do počítaného bodu x. Ohybový moment vyvolaný touto smykovou silou určíme jako součin poloviny tloušťky základu s hodnotou smykové síly Ti v příslušném místě x. Tabulka 2
Výsledné smykové síly Tx a ohybové momenty Mx po délce základu
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
Tx [kN]
237,3
237,3
235,9
229,0
213,2
186,6
Mx [kN·m]
59,3
59,3
59,0
57,3
53,3
46,6
6,5
7,5
8,0
148,4
98,3
35,8
0,0
37,1
24,6
9,0
0,0
Z výše uvedených hodnot je možno sestrojit pro větší představu grafy průběhů vnitřních sil po délce základové konstrukce. Tyto průběhy, včetně všech dalších zjištěných průběhů, jsou uvedeny rovněž v příloze. 3.2
Výpočet podle teorie pružnosti
Je-li podloží zatíženo poměrným vodorovným přetvořením +ε, pak toto zatížení vyvodí v konstrukci naprosto stejný stav, jako v případě zatížení samotné základové konstrukce poměrným vodorovným přetvořením –ε. Vzhledem k tomu, že zadání zatížení do základové konstrukce a následné řešení tohoto stavu je jednodušší, než by tomu bylo v případě přetváření podloží, s výhodou tuto možnost využijeme. Abychom mohli porovnávat takto stanovené účinky přetvoření na základovou konstrukci s plošným výpočtem dle ČSN 73 0039, musíme uvažovat shodné zatížení. Při výpočtu je tedy základová konstrukce zatížena poměrným vodorovným přetvořením o velikosti: ε = µ ε (ε max − ε eig )
ε = 0,85(5 ⋅ 10 −3 − 1 ⋅ 10 −3 ) = −3,4 ⋅ 10 −3
Toto poměrné vodorovné přetvoření můžeme zadat přímo, pokud to program umožňuje (NEXIS), nebo lze konstrukci zatížit změnou teploty (ANSYS), jenž toto přetvoření v konstrukci vyvodí podle vztahu ε = α ⋅ ∆T (1) kde
α ∆T
je součinitel teplotní délkové roztažnosti podle [4] je tato hodnota součinitele pro beton rovna α = 10·10-6°C-1 změna teploty
Z předchozího vyplývá, že požadovaného přetvoření ε = –3,4·10-3 v konstrukci docílíme − 3,4 ⋅ 10 −3 = 10 ⋅ 10 −6 ⋅ ∆T ∆T = −340°C jestliže ji zatížíme teplotní změnou o velikosti ∆T = −340°C . -5-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Jednotlivé výpočty byly provedeny pomocí programů NEXIS 32 a ANSYS, se vstupními materiálovými parametry: Vstupní materiálové parametry pro výpočet MKP
Tabulka 3
základová půda S3
základová konstrukce C 20/25
Edef,s = 20MPa
Edef,c = 29000MPa
-3
ρc = 2500 kg·m-3 νc = 0,2 αc = 10·10-6 [/°C]
ρs = 1750 kg·m νs = 0,3 cs = 0,0kPa
Je zřejmé, že dochází ke stěnovému namáhání prvků. Vzhledem k tomu, že tvar, zatížení i podepření tělesa je podél osy y neměnné, při výpočtu uvažujeme s rovinnou deformací, tzn. přetvoření εy = 0. Výpočtem MKP obdržíme výsledná normálová napětí v krajních vláknech prvku σmax a σmin (Obr. 2) a smykové napětí τx v základové spáře. y
σ z
=
σN =
N A
σM = ±
+
b = 1,0m
N ⋅e ⋅z Iy
σmin h = 0,5m
-
N
+
e
-
N
+
+
σmax
A = 0,5m2 Obr. 2
Průběh napětí po průřezu
Známe-li maximální a minimální hodnotu napětí v krajních vláknech a předpokládáme-li lineární průběh napětí po výšce prvku, lze ze vztahu pro napětí v prvku za předpokladu kombinace tahu (tlaku) a ohybu N M σ z = (− ) ± z (2) A Iy odvodit vztahy pro stanovení normálové síly v prvku Nx a její výstřednosti e: σ − σ min N = max ⋅A 2
e=
(σ max − σ min ) ⋅ I y
(3) (4)
N ⋅h
Platnost těchto vztahů byla ověřována na prostém betonovém nosníku, jenž byl vystaven známému zatížení. Výpočtem MKP lze určit napětí po průřezu a výslednými vztahy (3) a (4) stanovit vnesené zatížení. Vzhledem k tomu, že takto stanovené zatížení (na základě průběhu napětí po průřezu) bylo se zatížením vneseným totožné, pak lze prohlásit, že pro dané podmínky (lineární průběh napětí po průřezu, kombinace tahu a tlaku) jsou tyto vztahy správné. -6-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Pro ohybový moment pak platí M = N ⋅e
(5)
Smykovou sílu Tx získáme stejným způsobem jako v případě výpočtu podle ČSN 73 0039. 3.2.1
Výpočet programem NEXIS 32
Schéma zadání modelu řezu a vstupní materiálové vlastnosti jsou znázorněny na následujícím obrázku: x
ε, Edef,c, ρc, νc
u=0
u=0
z
Edef,s , ρs, νs w=0
Obr. 3
Tabulka 4
Schéma zadání vstupní geometrie a materiálových vlastností
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
-2,8
-8,6
-14,7
-21,5
-29,6
-40,7
-59,6
-116,2
-502,6
Tx [kPa]
-308,2
-307,5
-301,8
-290,2
-272,1
-246,5
-211,4
-161,2
-73,3
0,0
Nx [kN]
348,9
348,2
342,5
331,0
313,2
288,0
253,5
204,9
131,4
0,0
e [m]
0,143
0,144
0,153
0,171
0,197
0,228
0,263
0,293
0,283
-
Mx [kN·m]
50,0
50,3
52,5
56,6
61,6
65,8
66,7
60,0
37,2
0,0
∆x [mm]
0,0
1,7
5,1
8,6
12,0
15,4
18,8
22,3
25,7
27,4
Porovnáním maximální hodnoty smykového napětí τx s hodnotou mezní smykové únosnosti zeminy τu = 150kPa zjistíme, že smykové napětí na okraji základu je vyšší. K tomu je třeba přihlédnout při stanovení příslušné smykové síly Tx, kde mezní smyková únosnost τu omezí průběh napětí po délce základové konstrukce. 3.2.2
Výpočet programem ANSYS
Tabulka 5
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
-3,4
-10,4
-17,6
-25,5
-34,4
-45,9
-63,4
-107,6
-
Nx [kN]
424,8
424,0
417,5
404,1
383,4
354,4
315,2
260,4
165,0
0,0
e [m]
0,138
0,139
0,147
0,164
0,189
0,219
0,251
0,276
0,270
-
M [kN·m]
58,4
58,8
61,5
66,4
72,5
77,7
79,2
71,9
44,6
-
∆x [mm]
0,0
2,05
6,15
10,25
14,35
18,45
22,55
26,65
30,75
32,8
-7-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
3.2.3
Odchylky ve výpočtech programů NEXIS a ANSYS
Srovnáme-li hodnoty vnitřních sil získané MKP, pak je patrné, že v případě ANSYSu jsou nezanedbatelně vyšší. Tento rozdíl spočívá již v samotných rozdílných deformacích ∆x. Tabulka 6
Výsledné hodnoty posuvů ∆x po délce základové konstrukce dosažených různými modely NEXIS
ANSYS
x [m] od středu základu
0
8,0
0
8,0
Nx [kN]
348,9
0,0
424,8
0,0
∆x [mm]
0,0
27,4
0,0
32,8
Rozdíly ve vodorovných posunech spočívají v tom, že přetvoření působí pouze ve svislé rovině modelu, kdežto je-li přetvoření zadáno pomocí teploty, pak se zadaný součinitel teplotní délkové roztažnosti α pro případ rovinné deformace přepočítává. Potvrzení tohoto předpokladu bylo zjišťováno v ANSYSu 3D modelem jednotkové tloušťky (Obr. 4). Tento model byl zadán tak, aby představoval právě případ rovinné deformace, v němž změna teploty základové desky vyvolá přetvoření. Nejprve bylo řešení provedeno se stejnými vstupními materiálovými charakteristikami jako v případě 2D modelu, tzn. isotropními materiálovými konstantami. y
z
u=0 v=0 w=0
x
b = 1,0m
Obr. 4
Schéma zadání 3D modelu
Takto zadaný 3D model dosahuje naprosto stejných deformací jako 2D model, do něhož je přetvoření vneseno pomocí změny teploty. Z toho vyplývá, že vliv teploty se projevuje pravděpodobně i ve směru tloušťky prvku, což přesně neodpovídá pojetí zatížení vodorovným poměrným přetvořením a tím vznikají odchylky ve výpočtu. Vzhledem k tomu, že deformace u obou prvků jsou naprosto stejné, pak se vliv teploty v ostatních směrech musí projevit i u 2D modelu. Vlivu teploty v jiném než podélném směru je možno zabránit zadáním ortotropního součinitele teplotní délkové roztažnosti. S materiálovou hodnotu součinitele teplotní délkové roztažnosti budeme uvažovat pouze ve směru podélném, tedy ve směru přetváření, kdežto ve zbývajících směrech bude tato hodnota nulová, čímž se vliv teploty v těchto směrech vyloučí. Ostatní materiálové parametry zůstanou zachovány. Výpočtem 3D prvku se zadaným ortotropním součinitelem α bylo dosaženo naprosto stejných posuvů jako v případě 2D modelu řešeným programem NEXIS, který umožňuje zadat přetvoření přímo, tedy bez pomoci teploty.
-8-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Tabulka 7
Výsledné hodnoty posuvů ∆x po délce základové konstrukce dosažených různými modely NEXIS 2D model ε
ANSYS 2D model αiso
ANSYS 3D model αiso
ANSYS 3D model αortho
x [m] od středu základu
8,0
8,0
8,0
8,0
∆x [ mm]
27,4
32,8
32,8
27,4
Z výše uvedeného lze učinit závěr, že vodorovné přetvoření lze v konstrukci vyvolat teplotním zatížením prvku. Je však nutno zajistit, aby se teplotní změna projevila pouze ve směru poměrného vodorovného přetvoření, tedy materiál je nutno zadat s těmito charakteristikami: Konečné vstupní materiálové charakteristiky základové konstrukce
Tabulka 8
základová konstrukce C 20/25 αx = 10·10-6 [°C-1] αy = 0,0 [°C-1] αz = 0,0 [°C-1]
Edef,c = 29000MPa -3
ρc = 2500 kg·m νc = 0,2
Hodnoty vnitřních sil a posuvů získané řešením 3D modelu programem ANSYS při uvažování ortotropního součinitele α jsou uvedeny v následující tabulce. Tabulka 9
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce 3D modelu
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
-2,9
-8,7
-14,7
-21,2
-28,7
-38,2
-52,8
-89,7
-242,2
Tx [kPa]
-331,2
-330,5
-324,7
-313,0
-295,1
-270,1
-236,6
-191,1
-119,9
0,0
Nx [kN]
354,1
353,4
347,9
336,7
319,5
295,3
262,6
217,0
137,5
0,0
e [m]
0,138
0,139
0,148
0,165
0,189
0,219
0,251
0,276
0,270
-
Mx [kN·m]
48,8
49,1
51,4
55,4
60,4
64,8
66,0
59,9
37,2
0,0
∆x [mm]
0,0
1,7
5,1
8,6
12,0
15,4
18,8
22,3
25,7
27,4
Tyto vnitřní síly odpovídají přesně řešení 2D modelu, do něhož je přetvoření vneseno pomocí teplotní změny, při uvažování ortotropního součinitele α. Hodnoty smykových napětí jsou u konce základu vyšší než je příslušná smyková únosnost zeminy τu a proto je nutné jejich průběh pro výpočet smykové síly Tx omezit. Hodnoty vnitřních sil získané řešením 2D modelu programy ANSYS a NEXIS se sice liší, ale z praktického hlediska lze tyto rozdíly považovat za zanedbatelné.
4 VÝMĚNA PODLOŽÍ Snížení tření v základové spáře od vodorovného přetváření základové konstrukce či zeminy lze docílit výměnou vrstvy podloží bezprostředně pod základem za zeminu s nízkým deformačním modulem Edef. Jedná se především podle [2] o zeminy tříd F5 - F8, kde se tento modul pohybuje řádově v rozmezí 1,5 - 5,0 MPa. -9-
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Pod základovou konstrukci tedy použijeme například vrstvu zeminy třídy F5 o mocnosti 0,5 m, přesahující 1,0 m za okraj základu (Obr. 5). Okolní zemina zůstává nepozměněna. Rovněž přetvoření základové konstrukce je opět ε = –3,4·10-3. Zadané vstupní materiálové parametry: Vstupní materiálové parametry pro výpočet MKP
Tabulka 10
základová půda S3
4.1
základová půda F5
základová konstrukce C 25/30
Edef,s,1 = 20,0MPa
Edef,s,2 = 2,0MPa
Edef,c = 29000MPa
ρs,1 = 1750 kg·m-3 νs,1 = 0,3
ρs,2 = 2000 kg·m-3 νs,2 = 0,4
cs,1 = 0,0kPa
cs,2 = 30,0kPa
ρc = 2500 kg·m-3 νc = 0,2 αx = 10·10-6°C-1; αy = αz = 0,0°C-1
Výpočet podle ČSN 73 0039
Mocnost poddajnější vrstvy zeminy a1 = 0,5 m je menší než hloubka tlumící vrstvy a = 2,16 m a podle ČSN 73 0039 platí a = a1 = 0,5m Podle ČSN 73 0039 (tab. 4, příloha 1) určíme hodnotu součinitele η pro zeminu třídy F5: σ v = 0,24MPa b = 1,0m ⇒ η = 1,05 a 0,5 = = 0,5 b 1,0 Oedometrický modul přetvárnosti Eoed poddajnějsí základové půdy určíme dle [2]: 2ν 2 2 ⋅ 0,4 2 = 1− β = 1− = 0,467 1 −ν 1 − 0,4 1 1 E oed = E def = ⋅ 2,0 = 4,286 MPa β 0,467 Hodnota smykového napětí v základové spáře má výsledný tvar: τ x = β x ⋅ µ ε (ε max − ε eig ) ⋅η ⋅ E oed = β x ⋅ 0,85(5,0 ⋅ 10 −3 − 1,0 ⋅ 10 −3 ) ⋅ 1,05 ⋅ 4,286 ⋅ 10 3
τ x = 15,301 ⋅ β x Výsledné smykové napětí τx a posuvu ∆x po délce základu
Tabulka 11
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
0,0
1,0
3,0
5,0
7,0
9,0
11,0
13,0
15,0
16,0
βx
0,000
0,073
0,513
1,005
1,503
2,002
2,501
3,001
3,501
3,750
τx [kPa] ∆x [mm]
0,0
1,1
7,8
15,4
23,0
30,6
38,3
45,9
53,6
57,4
0,0
1,7
5,1
8,5
11,9
15,3
18,7
22,1
25,5
27,2
ξ=
x a
- 10 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Ověříme zda dojde k překročení smykové pevnosti zeminy: τ u = σ ⋅ tan ϕ + c = 240,0 ⋅ tan 32,0 o + 30,0
τ u = 180kPa ≥ τ max = 57,4kPa K prokluzu mezi základem a zeminou nedojde, průběh smykového napětí není třeba omezovat smykovou únosností zeminy. Stanovíme hodnotu podélné síly Tx a ohybového momentu Mx. Tabulka 12
Výsledné tahové síly Tx a ohybový moment Mx po délce základu
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
Tx [kN]
216,1
215,9
211,6
200,0
180,8
154,0
119,6
77,5
27,7
0,0
Mx [kN·m]
10,9
10,7
10,1
9,2
8,1
6,7
5,1
3,3
1,1
0,0
4.2
8,0
Výpočet podle teorie pružnosti
4.2.1
Výpočet programem NEXIS 32 x
ε, Edef ,c, ρc, νc z
Edef,s,1, ρs,1, νs,1
Edef,s,1, ρs,1, νs,1
Edef,s,1, ρs,1, νs,1
u=0
u=0
Edef,s,2, ρs,2, νs,2
w=0
Schéma zadání vstupní geometrie a materiálových vlastností
Obr. 5
Tabulka 13 Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
-2,4
-7,3
-12,3
-17,4
-22,8
-28,7
-35,6
-51,5
-152,4
Tx [kPa]
-202,6
-202,0
-197,2
-187,4
-172,5
-152,4
-126,7
-94,5
-51,0
0,0
Nx [kN]
123,8
123,2
120,0
114,5
105,3
91,5
75,2
55,4
29,2
0,0
e [m]
0,176
0,177
0,179
0,182
0,186
0,192
0,198
0,207
0,246
-
Mx [kN·m]
21,8
21,8
21,5
20,8
19,6
17,6
14,9
11,4
7,2
0,0
∆x [mm]
0,0
1,7
5,1
8,5
11,9
15,4
18,8
22,2
25,6
27,3
Maximální hodnota smykového napětí nepřesahuje mezní smykovou únosnost zeminy
τu = 180kPa.
- 11 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
4.2.2
Výpočet programem ANSYS
Tabulka 14
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
-1,4
-4,2
-7,1
-10,0
-13,1
-16,3
-19,7
-25,9
-58,1
Tx [kPa]
-105,3
-105,0
-102,2
-96,5
-88,0
-76,5
-61,8
-43,8
-21,0
0,0
Nx [kN]
120,9
120,5
117,6
111,8
103,0
91,1
76,0
57,4
31,2
0,0
e [m]
0,187
0,187
0,188
0,189
0,190
0,190
0,188
0,188
0,212
-
Mx [kN·m]
22,6
22,5
22,1
21,1
19,6
17,3
14,3
10,8
6,6
0,0
∆x [mm]
0,0
1,7
5,1
8,5
11,9
15,4
18,8
22,2
25,6
27,3
Ani v tomto případě není dosaženo mezní smykové únosnosti podloží τu = 180 kPa.
5 APLIKACE REOLOGICKÉ KLUZNÉ SPÁRY Nejprve je uvedeno řešení plošné na základě příspěvků [4] a [5] a dále opět řešení MKP podle teorie pružnosti. 5.1
Řešení plošné
Při zkouškách Technického a zkušebního ústavu stavebního (TZÚS) [4], [5] bylo zjištěno viskoelastické chování asfaltových izolačních pásů při dlouhodobém zatížení. Tohoto poznatku bylo využito a následně stanoven empirický vztah pro výpočet smykového napětí v základové spáře při aplikaci reologické kluzné spáry. Ovšem jak vyplývá z protokolu měření [5] tento vztah byl stanoven na základě čtyř hodnot měření za předpokladu lineární závislosti rychlosti posunu zkušebního bloku a hodnoty smykového napětí. Rovněž u teplotní závislosti je předpokládán lineární vztah mezi skutečným posunem a teplotou. Pro výpočet je nutno znát očekávanou rychlost vodorovného přetváření terénu vε,max. Jestliže příčinou vodorovného přetvoření terénu je hlubinné dobývání ložisek1 pak tuto maximální rychlost stanovíme pomocí vztahů uvedených v [5]. Budeme předpokládat hloubku těžené sloje pod povrchem terénu hs = 800 m a rychlost postupu porubní fronty vf = 40 m za měsíc. Nejprve určíme celkovou dobu poklesu bodu x 2 ⋅ hs 2 ⋅ 800 T= = = 40 měsíců vf 40 a maximální rychlost vodorovného přetváření terénu 8,925 ⋅ ε 8,925 ⋅ 0,005 vε ,max = = = 1,116 ⋅ 10 −3 m za měsíc t 40,0 =
1,116 ⋅ 10 −3 = 4,305 ⋅ 10 −10 m ⋅ s −1 30 ⋅ 24 ⋅ 3600
Smykový odpor τv proti vizkóznímu tečení ve vzdálenosti x od těžiště vodorovných sil v základu určíme podle vztahu uvedeného v [4]:
1
Jestliže vznik přetvoření má jiný původ, např. se může jednat o vlastní smršťování betonového základu, lze tuto rychlost samozřejmě stanovit odlišným způsobem.
- 12 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
τ vx = (1,5 − 0,1 ⋅ ∆T ) ⋅ 10 9 ⋅ vux + 1,5 kde
vux
je rychlost posuvu v kluzné spáře ve vzdálenosti x od těžiště vodorovných sil [m·s-1] vux = x ⋅ vε ,max
∆T
odchylka teploty v kluzné spáře od základní teploty T0 = 12 °C ∆T = T − 12
Předpokládáme-li průměrnou celoroční teplotu kluzné spáry T = 4 °C, pak ∆T = 4 − 12 = −8 o C
Výslednou sílu Tx určíme sumací po dílčích úsecích délky Li, šířky bi a průměrné hodnoty smykového napětí τi od okraje základu do počítaného bodu x. Ohybový moment vyvolaný touto smykovou silou stanovíme jako součin poloviny tloušťky základu s hodnotou smykové síly Ti v příslušném místě x. Tabulka 15
Výsledné hodnoty vnitřních sil a posuvu po délce základové konstrukce
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
vux [m·s-1]
0,0
2,153 E-10
6,458 E-10
1,076 E-09
1,507 E-09
1,937 E-09
2,368 E-09
2,798 E-09
3,229 E-09
3,444 E-09
5.2
τvx [kPa]
1,500
1,995
2,985
3,975
4,966
5,956
6,946
7,936
8,926
9,421
Nx [kN]
43,685
42,811
40,321
36,841
32,370
26,910
20,459
13,018
4,587
0,000
∆x [mm]
0,0
1,7
5,1
8,5
11,9
15,3
18,7
22,1
25,5
27,2
Mx [kN·m]
10,9
10,7
10,1
9,2
8,1
6,7
5,1
3,3
1,1
0,0
Výpočet podle teorie pružnosti
Pro výpočet základové konstrukce s aplikací reologické kluzné spáry MKP je použito výpočtového 2D modelu uvedeného již v předchozích případech, ovšem s tím rozdílem, že je upraven právě o vlastní konstrukci spáry. Výpočet základové konstrukce zatížené poměrným vodorovným přetvořením s aplikací reologické kluzné spáry podle teorie pružnosti je zaměřen na získání materiálových vlastností asfaltu použitého v kluzné spáře, za předpokladu, že výsledná smyková napětí budou srovnatelná s napětími získanými plošným řešením. Do výpočtu vstupují pro každý materiál tři základní parametry: Edef, ν a ρ. Poissonovo číslo, vzhledem k tomu, že v případě asfaltu se jedná o látku nestlačitelnou, můžeme předpokládat νa = 0,45. Rovněž objemovou hmotnost lze předem stanovit ρa = 1000 kg·m-3. Problém se tedy omezil na stanovení deformačního modulu Edef,a. Tento deformační modul byl stanoven experimentálně postupným snižováním hodnot modulu Edef,a a následným porovnáváním s plošným řešením. Schéma zadání vstupní geometrie, materiálových vlastností a vlastní konstrukce kluzné spáry je znázorněna na Obr. 6.
- 13 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
5.2.1
Výpočet programem NEXIS 32
ε, Edef ,c, ρc, νc Edef,a , ρa, νa
x
u=0
u=0
z
Edef,s , ρs, νs w=0
Obr. 6
Schéma zadání vstupní geometrie a materiálových vlastností včetně detailu konstrukce spáry
Tabulka 16
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
-0,5
-1,4
-2,4
-3,5
-4,2
-5,3
-6,3
-7,0
-7,7
Tx [kPa]
-30,7
-30,5
-29,6
-27,7
-24,7
-20,9
-16,1
-10,3
-3,7
0,0
Nx [kN]
37,8
37,6
36,4
34,1
30,6
25,8
19,9
12,8
4,7
0,0
e [m]
0,136
0,136
0,140
0,148
0,159
0,174
0,194
0,217
0,240
-
Mx [kN·m] ∆x [mm]
5,1 0,0
5,1 1,7
5,1 5,1
5,0 8,5
4,9 11,9
4,5 15,3
3,9 18,7
2,8 22,1
1,1 25,5
0,0 27,2
Abychom obdrželi výpočtem MKP vnitřní síly srovnatelné s hodnotami stanovenými plošným řešením je nutno zadat do výpočtu tyto vstupní materiálové charakteristiky asfaltu: Edef,a = 0,010 MPa νa = 0,45 ρa = 1000 kg·m-3 5.2.2
Výpočet programem ANSYS
V tomto výpočtu jsou již uvažovány získané materiálové vlastnosti asfaltu. Tabulka 17
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce
x [m] od středu základu
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
-0,5
-1,4
-2,4
-3,3
-4,3
-5,2
-6,1
-7,0
-7,2
Tx [kPa]
-30,1
-30,0
-29,0
-27,1
-24,3
-20,5
-15,7
-10,1
-3,6
0,0
Nx [kN]
31,2
31,0
30,1
28,1
25,2
21,3
16,5
10,6
4,0
0,0
e [m]
0,135
0,135
0,139
0,146
0,157
0,172
0,191
0,214
0,228
-
Mx [kN·m]
4,2
4,2
4,2
4,1
4,0
3,7
3,1
2,3
0,9
0,0
∆x [mm]
0,0
1,7
5,1
8,5
11,9
15,3
18,7
22,1
25,5
27,2
- 14 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Právě hodnota modulu pružnosti asfaltu je velmi důležitý vstupní parametr pro výpočet MKP. Vzhledem k tomu, že asfalt je viskoelastická látka, tento modul se velmi výrazně mění s časem a teplotou. Při dlouhé době zatížení a vysoké teplotě jeho hodnota klesá, kdežto při nízkých teplotách a krátkodobém zatížení nabývá hodnot vyšších. Z tohoto důvodu se oproti látkám pevným namísto označení modulu pružnosti E užívá název modul tuhosti S(T,t), který je závislý na teplotě a době zatížení. Vzhledem k tomu, že do výpočtů tento modul vstupuje jako materiálová konstanta zůstaneme u názvosloví pro látku pevnou, tedy modulu pružnosti E. Z protokolu měření [5] se podařilo sestavit závislost modulu pružnost na čase, a to tak, že příslušnými body měření byla proložena regresní mocninná křivka. Bohužel z protokolu není zcela zřejmá přesná skladba kluzné vrstvy. Rovněž je toto měření vázáno na teplotu při které bylo prováděno. Modul pružnosti asfaltu můžeme určit pomocí následujících vztahů. Známe-li celkovou deformaci pásu γ a hodnotu příslušného smykového napětí τ, pak modul pružnosti ve smyku stanovíme pomocí vztahu G=
τ γ
(6)
mezi modulem pružnosti a modulem pružnosti ve smyku platí závislost G= E 21 + ν
(7)
kde v je Poissonovo číslo; pro asfalty jako látku nestlačitelnou předpokládáme v = 0,5 a odtud E = 3⋅G Tabulka 18
(8)
Určení modulů pružnosti závislých na čase z měření TZÚSu
Průměrné smykové napětí
Celkový posun
Celkové poměrné zkosení
Modul pružnosti ve smyku
Modul pružnosti
Průměrná teplota
τ [Pa]
s [mm]
γ [–]
G(t) [Pa]
E(t) [Pa]
[°C]
1
2373
0,820
0,0683
34 728
104 183
18,0
2
2373
1,190
0,0992
23 930
71 790
18,0
3
2373
1,560
0,1300
18 254
54 763
18,0
4
2373
1,890
0,1575
15 067
45 201
18,0
5
2373
2,220
0,1850
12 827
38 482
18,0
6
2373
2,550
0,2125
11 167
33 502
18,0
7
4633
3,400
0,2833
16 351
49 052
21,0
8
4633
5,380
0,4483
10 333
31 000
21,0
9
4633
7,230
0,6025
7 689
23 067
21,0
Den
- 15 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Závislost modulu pružnosti na čase získaná proložením bodů mocninnou regresí
Graf 1
120 000 y = 105782x -0,593 R2 = 0,8742
100 000
Ea [Pa]
80 000 60 000 40 000 20 000 0 0
1
2
3
4
5
6
7
8
9
den
Ea [Pa]
Mocninný (Ea [Pa])
Závislost modulu pružnosti na čase má tedy výsledný tvar E a = 105782 ⋅ d −0,593 kde
d
(9)
představuje dobu ve dnech.
Několik příkladů hodnot modulů pružnosti pro různé časy stanovené tímto vztahem jsou uvedeny v následující tabulce: Tabulka 19
Moduly pružnosti stanovené pomocí vztahu získaného ze zkoušek TZÚS
počet dní Modul pružnosti Ea [Pa]
E a = 105782 ⋅ d −0,593 počet dní
Modul pružnosti Ea [Pa]
E a = 105782 ⋅ d −0,593
1
7
14
30
60
90
120
150
365
33363
22119
14076
9332
7337
6187
5420
4864
3199
730 (2 roky)
1825 (5 let)
3560 (10 let)
2121
1232
829
Modul pružnosti v závislosti na čase a teplotě můžeme stanovit také pomocí Van der Poelova nomogramu [7] uvedeného v příloze. Pro stanovení je ovšem nezbytné znát hodnotu penetrace a bodu měknutí KK použitého asfaltu. Při zkouškách TZÚSu při určování vhodného typu asfaltového pásu pro konstrukci kluzné vrstvy byla potřebná hodnota penetrace měřena, ovšem jehlou o hmotnosti 300 g po dobu 10 s. Podle [3] a potažmo i pro stanovení modulu pružnosti pomocí nomogramu se penetrace stanovuje jehlou o hmotnosti 100g po dobu 5 při 25 °C. Nebudeme-li tuto skutečnost uvažovat, pak na základě naměřené hodnoty penetrace asfaltového pásu pen = 15,0 uvedené v [6] a hodnoty bodu měknutí KK pro běžné asfaltové pásy TKK = 85,0 °C2 lze stanovit potřebnou hodnotu penetračního indexu PI: PI =
2
log 800 − log p 25 20(1 − 25 A) , kde A = 1 + 50 A TKK − 25
Tuto hodnotu lze vyhledat v nabídce asfaltových pásů různých výrobců.
- 16 -
(10), (11)
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
A=
log 800 − log 15,0 = 0,0288 85,0 − 25
PI =
20(1 − 25 ⋅ 0,0288) = 2,3 1 + 50 ⋅ 0,0288
Samozřejmě tyto hodnoty se týkají pouze asfaltu, ale přesto je využijeme pro srovnání. Asfaltový pás se bude chovat zřejmě odlišně než vzorek samotného asfaltu, neboť obsahuje navíc nosnou vložku. Na základě stanovených hodnot PI a bodu měknutí KK již můžeme odměřením z nomogramu určit hodnoty modulů pružnosti pro různé časy při teplotě T = 18 °C (teplota zkoušek TZÚS). Tyto hodnoty jsou uvedeny v následující tabulce společně s hodnotami stanovenými podle vztahu (9). V tabulce je ještě navíc uveden jeden případ. Vzhledem k rozdílnému způsobu stanovení penetrace bychom při menším zatížení jehly zcela jistě obdrželi hodnotu menší, odhadem pen = 8,0. Tabulka 20
Srovnávací tabulka modulů tuhosti asfaltu Modul pružnosti Ea [Pa] pro teplotu T = 18 °C Van der Poelův nomogram
Mocninná křivka získaná ze zkoušek TZÚS
pen 15,0 KK 85 °C PI +2,3
pen 8,0 KK 85 °C PI +1,3
E = 105782·d-0,593
1 den
65 000
180 000
105 782
30 dní
6 000
15 000
14 076
1 rok
1 000
1 500
3 199
Při srovnání hodnot je nutno uvážit několik faktorů. Jednak křivka a potažmo získaný vztah je stanoven na základě ne zcela přesných informací a také od hodnot modulů stanovených nomogramem nelze očekávat významnější přesnost. Obzvláště uvážíme-li, že nomogram je sestaven empiricky ze vzorků cca 40 druhů asfaltů [7]. Vzhledem k těmto okolnostem lze přesto konstatovat, že hodnoty se poměrně shodují. Zvláště pak vývoj hodnot modulů v čase je velmi podobný. Hodnota modulu pružnosti asfaltu Ea = 0,01 MPa, uvažovaná při prvním výpočtu s aplikací kluzné spáry by podle Tabulky 19 odpovídala přibližně jednomu měsíci zatížení. Ovšem v základové spáře při běžném provozu nelze uvažovat s teplotou T = 18,0 °C. Z tohoto důvodu by asfalt modulu pružnosti dosáhl hodnoty 0,01 MPa až za delší časový úsek. Výpočet plošným řešením vůbec tyto materiálové charakteristiky neuvažuje jako v případě řešení MKP. Celý výpočet je tak vlastně založen na rychlosti přetváření podloží a teplotě. Poslední uvedený výpočet je proveden pro případ, kdy hodnota modulu tuhosti asfaltu klesne na Ea = 0,001 MPa, což zhruba podle Tabulky 19 odpovídá době zatížení jednoho roku.
- 17 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
5.2.3
Výpočet programem NEXIS pro Ea = 0,001 MPa
Tabulka 21
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce
x [m] od středu základu
5.2.4
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
τx [kPa]
0,0
0,0
-0,1
-0,2
-0,3
-0,4
-0,5
-0,6
-0,7
-0,2
Tx [kPa]
-2,7
-2,7
-2,6
-2,5
-2,2
-1,9
-1,4
-0,9
-0,2
0,0
Nx [kN]
3,8
3,8
3,6
3,4
3,1
2,5
2,0
1,3
0,5
0,0
e [m]
0,228
0,229
0,232
0,233
0,235
0,243
0,248
0,253
0,250
-
Mx [kN·m]
0,867
0,865
0,840
0,792
0,717
0,615
0,483
0,317
0,119
0,000
∆x [mm]
0
1,7
5,1
8,5
11,9
15,3
18,7
22,1
25,5
27,2
Výpočet programem ANSYS Ea = 0,001 MPa
Tabulka 22
Výsledné hodnoty vnitřních sil a posuvů po délce základové konstrukce
x [m] od středu základu
τx [kPa] Tx [kPa] Nx [kN] e [m] Mx [kN·m] ∆x [mm]
0
0,5
1,5
2,5
3,5
4,5
5,5
6,5
7,5
8,0
0,0 -6,2 3,1 0,226 0,708 0,0
0,0 -6,2 3,1 0,226 0,705 1,7
-0,1 -6,1 3,0 0,227 0,686 5,1
-0,2 -5,9 2,8 0,229 0,647 8,5
-0,3 -5,6 2,5 0,232 0,588 11,9
-0,4 -5,2 2,1 0,236 0,505 15,3
-0,5 -4,8 1,7 0,241 0,397 18,7
-0,6 -4,2 1,1 0,246 0,262 22,1
-0,7 -3,6 0,4 0,239 0,096 25,5
-0,7 0,0 0,0 0,000 27,2
Z výsledných hodnot vnitřních sil je patrné, že snížením hodnoty modulu pružnosti asfaltu Ea o jeden řád dojde k velmi výraznému snížení těchto sil. Z praktického hlediska jsou tyto síly již zanedbatelné.
6 ZÁVĚR Srovnávací výpočty ukázaly, že účinky poměrného vodorovného přetvoření na konstrukce lze velmi dobře počítat MKP při použití 2D modelu svislého řezu základové konstrukce včetně podloží. Přetvoření je však nutno zadat přímo, nebo jej vnést pomocí teplotní změny, avšak s vyloučením jejího vlivu ve směrech odlišných od směru přetváření, tedy za použití ortotropního součinitele teplotní délkové roztažnosti. Srovnáme-li průběhy tahových sil získané řešením stavu, v němž je pouze podloží resp. základová deska zatížena poměrným vodorovným přetvořením, podle ČSN 73 0039 a MKP, pak lze shledat, že síly odpovídající plošnému řešení jsou výrazně nižší. Tento výrazný rozdíl vzniká díky odlišnému pojetí obou výpočtů. V případě řešení dle ČSN 73 0039 se hodnoty tahových sil získávají z průběhu napětí v základové spáře, a tudíž se jedná o hodnotu síly v základové spáře. Tato síla se prohlašuje za působící tahovou sílu v prvku. Řešením MKP se získávají normálové síly v základové konstrukci z průběhu normálových napětí po výšce prvku. Při řešení MKP není vyloučena tahová únosnost zeminy, což je zcela jistě také jeden z důvodů vyšších hodnot vnitřních sil. Rovněž do výpočtu podle ČSN 73 0039 nevstupují materiálové vlastnosti základové konstrukce jako v případě řešení podle teorie pružnosti. - 18 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
Velice zajímavé výsledky přináší výpočty, při nichž je uvažováno s výměnou základové půdy. Zatímco při výpočtu dle ČSN 73 0039 nebylo dosaženo výraznějšího snížení podélné tahové síly a potažmo smykového napětí v základové spáře, výpočtem dle teorie pružnosti byly zjištěny více než 60% poklesy těchto vnitřních sil. Při stanovení této síly plošným řešením se uvažuje pouze s vrstvou poddajnější zeminy. Řešení MKP však umožňuje využití základové půdy i pod touto poddajnou vrstvou. Při aplikaci reologické kluzné spáry dochází k velmi výraznému poklesu smykového napětí v základové spáře. Otázkou však zůstává, jak je vztah plošného řešení pro určení tohoto napětí přesný, neboť, jak již bylo uvedeno, je stanoven na základě čtyř měření proložených lineární regresí, což se jeví jako poměrně malý počet vstupních údajů pro určení přesnější závislosti. Rovněž teplotní závislost je určena za předpokladu lineárního vztahu mezi skutečným posunem a teplotou. V tomto případě se vůbec neuvažuje se změnou materiálových vlastností asfaltového pásu s rostoucí dobou působícího zatížení, a celý výpočet je tak založen na rychlosti posunu podloží a teplotě v základové spáře. Výpočet MKP může tuto skutečnost postihnout a navíc má tu výhodu, že umožňuje stanovit množství dalších výstupů. Ovšem při výpočtu MKP je nutno znát vstupní materiálové charakteristiky asfaltu, potažmo aplikovaného asfaltového pásu, což se vzhledem k malému množství zjištěných dat týkajících se chování asfaltu z hlediska dlouhodobého zatěžování jeví jako podstatný problém. Nejvýraznějšího snížení vnitřních sil bylo dosaženo aplikací reologické kluzné spáry, ovšem k zajímavému poklesu dojde i při výměně základové půdy řešením MKP. Míru účinnosti jednotlivých opatření lze porovnat pomocí vybraných průběhů vnitřních sil, jenž jsou umístěny v příloze. Modul tuhosti z hlediska krátkodobého zatížení je znám vzhledem k jeho zásadnímu použití v silničním stavitelství a stanovuje jej celá řada zkušeben. Ovšem z praktického hlediska pro případ zatížení poměrným vodorovným přetvořením podloží jsou tyto závislosti prakticky nevyužitelné, neboť toto zatížení je dlouhodobé. Viskoelastické vlastnosti asfaltových pásů by se při krátkodobém zatížení ani neprojevily. V případě krátkodobého přenosu sil lze v základové spáře uvažovat pouze s fyzikálním třením.
Asfalt také při zatížení mění svou viskozitu a to velmi výrazně. Avšak vzhledem k relacím, v jakých se nacházejí běžná napětí v základových spárách, lze tuto změnu bez problémů zanedbat. Z hlediska dlouhodobého zatěžování asfaltů je možnost získání vstupních parametrů velice skromná a v podstatě se je kromě zkoušek TZÚSu nepodařilo získat. Jistou možnost poskytuje Van der Poelův nomogram, avšak ten neumožňuje, alespoň ne z hlediska dlouhodobého, stanovení modulu tuhosti s větší přesností. Rovněž se týká pouze „čistých“ asfaltů a dá se zcela jistě předpokládat, že asfaltové pásy budou mít hodnoty jiné. Liší se například tím, že navíc obsahují nosnou vložku z jiného materiálu a v dnešní době se asfaltové pásy modifikují polymery, takže jeho využití je poměrně omezeno. Právě využití polymery modifikovaných asfaltových pásů, vzhledem k tomu, že se modifikace polymery využívá právě k vylepšení vlastností „klasických“ asfaltových pásů, by mohl přinést zajímavý posun v této problematice. Výpočetní systém ANSYS by zcela jistě umožnil simulaci stavu dlouhodobého přetváření či jeho předpokládané změny. Rovněž umožňuje zadání viskoelastického materiálu, tedy materiálu, který mění své vlastnosti v závislosti na teplotě a době zatížení. K tomu je však nezbytně nutné znát potřebné vstupní údaje.
- 19 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
SEZNAM LITERATURY: [1] [2] [3] [4] [5] [6] [7]
ČSN 73 0039 Navrhování objektů na poddolovaném území, Vydavatelství norem, Praha 1990 ČSN 73 1001 Základová půda pod plošnými základy, Vydavatelství Úřadu pro normalizaci a měření, Praha 1987 ČSN 65 7062 Penetrační zkouška ČSN P ENV 1992-1-1 Navrhování betonových konstrukcí, Český normalizační institut, Praha 1994 Balcárek V., Bradáč J.: Použití asfaltových izolačních pásů jako kluzné spáry staveb na poddolovaném území, Pozemní stavby 2 – 1982 Kraus L., Balcárek V., Kostka F.: Studie možnosti snížení tření v základové spáře staveb na poddolovaném území pomocí živičných izolačních pásů; protokoly měření TZÚS, Ostrava 1980 Van der Poel, C.: A general system describing the visco-elastic properties of bitumen and its relation to rountine test data, Journal of Applied Chemistry, vol 4, pp 221 - 236, 1954
- 20 -
Studentská vědecká odborná činnost 2003 VŠB – Technická univerzita Ostrava, Fakulta stavební Ludvíka Podéště 1875, 708 33 Ostrava-Poruba http://fast.vsb.cz
PŘÍLOHY Metodika výpočtu smykového napětí od přetvoření dle ČSN 73 0039 Metodika výpočtu dle ČSN 73 0039 – poddajnější zemina Metodika výpočtu při aplikaci kluzné spáry Průběhy vnitřních sil Srovnání účinku metod snižující tření v základové spáře Van der Poelův nomogram
- 21 -