09.12.2012 Views

STAVEBNÍ OBZOR č. 06/2004

STAVEBNÍ OBZOR č. 06/2004

STAVEBNÍ OBZOR č. 06/2004

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>STAVEBNÍ</strong> <strong>OBZOR</strong> 6/<strong>2004</strong> 173<br />

finální podoby<br />

Porovnání řešení Boussinesqova<br />

s pružnou vrstvou<br />

Obě řešení porovnáme na příkladu, kdy svisle na podloží<br />

působí konstantní plošné zatížení v kruhu. Pro srovnání<br />

použijeme výpo<strong>č</strong>et průběhu napětí na ose rota<strong>č</strong>ní symetrie.<br />

V případě Boussinesqova řešení pomocí rovností (3), (5)<br />

integrací získáme<br />

3 ⎡ z ⎤<br />

σ z = f0<br />

⎢1−<br />

.<br />

2 2 3/<br />

2 ⎥ (19)<br />

⎣ ( r + z ) ⎦<br />

Správnost řešení můžeme ověřit v [3].<br />

V případě pružné vrstvy (zadání přibližuje obr. 3) řešíme<br />

diferenciální rovnice, které dostaneme z (18) transformací<br />

do polárních souřadnic. V případě polárních souřadnic je<br />

rovnováha vyžadována soustavou rovnic<br />

w<br />

j<br />

j−1<br />

2 ( j ) w = ( −1)<br />

2 f .<br />

∆w − α (18)<br />

Obr. 3. Rovnoměrné spojité zatížení působící v kruhu<br />

1<br />

j−1<br />

2<br />

+ w ( ) ( 1)<br />

2<br />

j,<br />

x − jα<br />

w j = − f . (20)<br />

x<br />

j,<br />

xx<br />

z<br />

Ty jsou známy jako modifikované diferenciální rovnice<br />

Besselovy. Vlivem rota<strong>č</strong>ní symetrie došlo k další redukci<br />

dimenze. Znamená to, že nyní již řešíme úlohu jednodimenzionální.<br />

Té přísluší okrajové podmínky<br />

w ( ∞)<br />

= 0 ,<br />

j<br />

V jednotlivých intervalech využijeme známé řešení<br />

kde<br />

j<br />

dw j ( x)<br />

= 0 .<br />

dx<br />

x = 0<br />

w<br />

x ∈ 0; r w<br />

2<br />

j<br />

0 α 2<br />

j<br />

x ∈ r<br />

j<br />

j−1<br />

p1<br />

() x = K 1 I ( j x)<br />

+ ( −1)<br />

,<br />

j<br />

) w () x = K 2 K ( jαx),<br />

; ∞ j<br />

0<br />

α =<br />

C1 f z f z<br />

, f z = − , wp1<br />

=<br />

C2<br />

C2<br />

C1<br />

z<br />

,<br />

(21)<br />

(22)<br />

Dvě zbývající integra<strong>č</strong>ní konstanty 1, K 2 ur<strong>č</strong>íme z podmínek<br />

spojitosti průhybu a první derivace v bodě x = r.<br />

Dosazením do (22) získáme soustavu dvou lineárních algebraických<br />

rovnic<br />

odkud pomocí Kramerova pravidla s využitím identity<br />

vypo<strong>č</strong>teme<br />

j<br />

j<br />

K<br />

K<br />

I<br />

0<br />

2<br />

H.<br />

π<br />

C 1 = Eoed<br />

, C<br />

2<br />

2 =<br />

8H<br />

1<br />

1<br />

I<br />

I<br />

0<br />

( jαr)<br />

K ( jαr)<br />

+ I ( jαr)<br />

K ( jαr)<br />

j<br />

j<br />

1<br />

1<br />

Před zpětným dosazením výsledků (22) a (25) do (7) si<br />

připomeneme identitu<br />

∑ ∞<br />

n=<br />

0<br />

[ ( 2 + 1)<br />

z]<br />

K j j<br />

n sin n π . z<br />

( −1)<br />

= .<br />

2<br />

( 2n<br />

+ 1)<br />

4<br />

Zavedeme-li j = 2n+1, n = 0, 1, 2, ..., pak v jednotlivých<br />

intervalech můžeme psát<br />

−∑ ∞<br />

n=<br />

0<br />

⎧<br />

⎨<br />

⎩<br />

⋅ ∑ ∞<br />

n=<br />

0<br />

Odtud pomocí (9) dopo<strong>č</strong>teme průběh svislého napětí σ z.<br />

Opět po intervalech píšeme<br />

⋅∑ ∞<br />

n=<br />

0<br />

⎧<br />

⎨<br />

⎩<br />

⋅ ∑ ∞<br />

Jistě nás bude zajímat extrémní hodnota normálového<br />

0<br />

H<br />

G .<br />

2<br />

j−1<br />

p1<br />

j<br />

( jαr)<br />

+ ( −1)<br />

2 = K 2 K ( jαr)<br />

,<br />

j<br />

( jαr)<br />

= − K 2 K ( jαr),<br />

1<br />

w<br />

K 2 =<br />

j<br />

p1<br />

w<br />

j<br />

w<br />

K 1 = −<br />

j<br />

w<br />

2<br />

p1<br />

1<br />

j−1<br />

( −1)<br />

2 αrK<br />

( jαr)<br />

j−1<br />

( −1)<br />

2 αrI<br />

( jαr)<br />

.<br />

x ∈<br />

0; r<br />

1<br />

,<br />

2 w ⎧ π<br />

⎨<br />

H ⎩8H<br />

1<br />

0<br />

1<br />

=<br />

jαr<br />

p1<br />

n<br />

( x,<br />

z)<br />

z − − ( −1)<br />

= ∑ ∞<br />

n ( −1)<br />

K [ ( 2n<br />

+ 1)<br />

αr]<br />

⋅ I [ ( 2n<br />

+ 1)<br />

αx]<br />

,<br />

n=<br />

0<br />

( 2n<br />

1)<br />

(23)<br />

(24)<br />

(25)<br />

αr<br />

K<br />

2n<br />

+ 1<br />

αr<br />

⎡ + π ⎤⎫<br />

1<br />

0 sin<br />

⎬ ,<br />

2 1<br />

⎢<br />

z<br />

n +<br />

⎣ 2H<br />

⎥<br />

⎦⎭<br />

w<br />

x ∈ r;<br />

∞),<br />

w<br />

H<br />

⎧<br />

⎨<br />

⎩<br />

αr<br />

2n<br />

+ 1<br />

p1<br />

n<br />

( x,<br />

z)<br />

⋅ ⋅ ( −1)<br />

I (<br />

= ∑ ∞<br />

n=<br />

0<br />

n ( −1)<br />

I [ ( 2n<br />

+ 1)<br />

αr]<br />

⋅ K [ ( 2n<br />

+ 1)<br />

αx]<br />

( 2n<br />

1)<br />

(26)<br />

αr<br />

⎡ + π ⎤⎫<br />

1<br />

0 sin<br />

⎬ .<br />

2 1<br />

⎢<br />

z<br />

n +<br />

⎣ 2H<br />

⎥<br />

⎦⎭<br />

σ<br />

z<br />

x ∈<br />

⎧<br />

⎨<br />

⎩<br />

0; r<br />

,<br />

E<br />

2r<br />

H<br />

oed<br />

n<br />

( x,<br />

z)<br />

f 1−<br />

⋅ ⋅ ( −1)<br />

K (<br />

= ∑<br />

0<br />

∞<br />

z<br />

G<br />

n=<br />

n ( −1)<br />

K [ ( 2n<br />

+ 1)<br />

αr]<br />

⋅ I [ ( 2n<br />

+ 1)<br />

αx]<br />

n=<br />

0<br />

1<br />

σ<br />

z<br />

0<br />

x ∈ r;<br />

∞),<br />

E<br />

( 2n<br />

1)<br />

1<br />

1<br />

1<br />

[2<br />

[2<br />

⎡ + π ⎤⎫<br />

cos<br />

⎢<br />

z ⎬ ,<br />

⎣ 2H<br />

⎥<br />

⎦⎭<br />

2r<br />

H<br />

⎧<br />

⎨<br />

⎩<br />

oed<br />

n<br />

( x,<br />

z)<br />

f ⋅ ⋅ ( −1)<br />

I (<br />

= ∑<br />

0<br />

∞<br />

z<br />

G<br />

n=<br />

n ( −1)<br />

I [ ( 2n<br />

+ 1)<br />

αr]<br />

⋅ K [ ( 2n<br />

+ 1)<br />

αx]<br />

1<br />

0<br />

( 2n<br />

1)<br />

(27)<br />

[2n<br />

⎡ + π ⎤⎫<br />

cos<br />

⎢<br />

z ⎬ .<br />

⎣ 2H<br />

⎥<br />

⎦⎭<br />

1

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!