13.11.2012 Views

Veza naprezanja i deformacija

Veza naprezanja i deformacija

Veza naprezanja i deformacija

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 1 -<br />

4. VEZA IZMEðU NAPREZANJA I DEFORMACIJA<br />

S a d r ž a j:<br />

4.1 Linearizacija općeg zakona ponašanja materijala<br />

4.2 Tenzor koeficijenata elastičnosti<br />

4.3 Heksagonalna simetrija (uslojeni materijal)<br />

4.4 Laméove konstante elastičnosti i opći Hookov zakon<br />

4.5 Fizikalno značenje modula elastičnosti<br />

4.5.1 Jednoosno naprezanje<br />

4.5.2 Čisti posmik<br />

4.5.3 Hidrostatički tlak<br />

4.5.4 Zakon ponašanja u termoelastičnosti<br />

4.6 Potencijalna energija<br />

4.7 Duhamel – Neumanove jednadžbe<br />

4.8 Primjeri


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 2 -<br />

4.1 Linearizacija općeg zakona ponašanja materijala<br />

Općeniti zakon ponašanja materijala se izražava kao funkcionalna veza izmeñu tenzora<br />

<strong>deformacija</strong> i tenzora <strong>naprezanja</strong>:<br />

σ = S( ε)<br />

(4.1)<br />

Svaka funkcijska ovisnost se dade izraziti beskonačnim redom potencija (Mac Laurin,<br />

Taylor), te nakon odbacivanja članova s višim redom derivacija, preostaju prva dva člana<br />

beskonačnog reda kod čega je prvi konstantan i opisuje referentno naprezanje, a drugi daje<br />

linearnu vezu izmeñu <strong>deformacija</strong> i pripadnog <strong>naprezanja</strong>.<br />

Na ovaj način općeniti zakon ponašanja je lineariziran, a materijal koji se ponaša po ovom<br />

zakonu često se naziva Hookeov materijal.<br />

Teorija temeljena na Hookeovom zakonu je linearna teorija elastičnosti, a meñuovisnost<br />

<strong>deformacija</strong> i <strong>naprezanja</strong> se izražava izrazom:<br />

gdje je:<br />

σij – tenzor <strong>naprezanja</strong><br />

0<br />

ij = σij<br />

+ Cijkl εkl<br />

σij 0 – tenzor referentnog <strong>naprezanja</strong><br />

Cijkl – tenzor koeficijenata elastičnosti<br />

εkl – tenzor <strong>deformacija</strong><br />

i,j,k,l – indeksi koordinatnih smjerova<br />

σ (4.2)<br />

U daljnjem tekstu su nabrojena osnovna svojstva tenzora Cijkl, kao i tenzori koeficijenata<br />

elastičnosti za neke vrste materijala od kojih je posebno istaknut materijal koji zadovoljava<br />

heksagonalni uvjet simetrije.<br />

4.2 Tenzor koeficijenata elastičnosti<br />

Linerano elastično ponašanje materijala opisano je općenito tenzorom koeficijenata<br />

elastičnosti Cijkl i u trodimenzionalnom prostoru ima 81 komponentu.<br />

Meñutim, dokazano je da uvjeti ravnoteže dovode do simetričnosti tenzora <strong>naprezanja</strong> što ima<br />

za posljedicu da tenzor koeficijenata elastičnosti mora biti simetričan po prva dva indeksa:<br />

Cijkl = Cjikl (4.3)<br />

Takoñer, simetričnost tenzora <strong>deformacija</strong> εij = εji povlači simetričnost tenzora koeficijenata<br />

elastičnosti po preostala dva indeksa:<br />

Cijkl = Cijlk (4.4)<br />

Činjenica da rotacija krutog tijela ne dovodi do promjene <strong>naprezanja</strong> postavlja uvjet<br />

simetričnosti po indeksima k,l tenzora koeficijenata elastičnosti, no kako je taj uvjet već prije<br />

zadovoljen vrijedi i inverzan zaključak.


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 3 -<br />

Simetričan tenzor koeficijenata elastičnosti svaki antisimetrični tenzor <strong>deformacija</strong> će prevesti<br />

u nultenzor, a nesimetrični tenzor će “simetrirati”.<br />

Npr:<br />

σ = ε ; k,<br />

l = 1, 2, 3<br />

(4.5)<br />

13<br />

C13kl kl<br />

ili uz izjednačavanje komponenata simetrično po indeksima k,l:<br />

13<br />

Slijedi:<br />

1311<br />

11<br />

1322<br />

C1312 = C1321 ; C1313 = C1331 ; C1323 = C1332 (4.6)<br />

22<br />

1333<br />

33<br />

1312<br />

( ε + ε ) + C ( ε + ε ) + C ( ε + )<br />

σ = C ε + C ε + C ε + C<br />

ε (4.7)<br />

12<br />

Isti rezultat se dobiva računajući sa simetričnim tenzorom <strong>deformacija</strong> εij iz čega se može<br />

zaključiti da se može računati s gradijentom polja pomaka bez prethodnog simetriranja, a<br />

saznanje da je tenzor <strong>deformacija</strong> simetričan koristi samo utoliko da se postavi uvjet na tenzor<br />

koeficijenata elastičnosti Cijkl, da mora biti simetričan po drugom paru indeksa (k, l).<br />

Tenzor Cijkl koji opisuje elastična svojstva realnog materijala mora zadovoljiti uvjet pozitivne<br />

definitnosti:<br />

21<br />

1313<br />

Cijkl ij kl<br />

ij<br />

kl<br />

13<br />

31<br />

1323<br />

ε ε > 0;<br />

∀ε<br />

≠ 0 , ∀ε<br />

≠ 0<br />

(4.8)<br />

gdje su εij, εkl proizvoljni tenzori drugog reda različiti od nultenzora, što podrazumijeva da je<br />

determinanta njihove matrice različita od nule.<br />

Navedena svojstva simetričnosti reduciraju broj različitih komponenata tenzora Cijkl s 81 na<br />

36 komponenata.<br />

Do daljnjeg pojednostavljenja tenzora Cijkl dolazi se na temelju analize termodinamičkih<br />

procesa na elastičnom kontinuumu.<br />

Slobodna energija u jedinici termoelastičnog volumena uz prethodnu linearizaciju je opisana<br />

funkcijom:<br />

gdje je:<br />

1 ( , T)<br />

= C ε ε − β ε Θ + F ( Θ)<br />

εij ijkl ij kl ij ij 0<br />

F (4.9)<br />

2<br />

Θ = T−T0 - razlika temperature izmeñu početnog i promatranog trenutka<br />

F0 (Θ) - suma svih članova koji ovise samo o temperaturi<br />

βij - tenzor temperaturnih svojstava materijala<br />

εij, εkl<br />

Cijkl<br />

- tenzor <strong>deformacija</strong><br />

- tenzor koeficijenata elastičnosti<br />

F (εij, T) - termodinamička funkcija - slobodna energija<br />

Za slučaj da se tijelo nalazi u sredini s konstantnom temperaturom, termodinamička funkcija<br />

postaje:<br />

1<br />

( ij,<br />

T = const)<br />

= Cijkl<br />

εij<br />

kl<br />

F ε ε<br />

(4.10)<br />

2<br />

23<br />

32


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 4 -<br />

i predstavlja izvršeni rad uslijed deformacije elastičnog tijela pri izotermičkim uvjetima.<br />

⎛ F ⎞<br />

dF ⎜<br />

∂<br />

= ⎟<br />

⎜ ⎟<br />

= Cijkl<br />

εkl<br />

dε<br />

⎝<br />

∂εij<br />

⎠<br />

T=<br />

const<br />

Očigledno da je dF totalni diferencijal funkcije F pa mora biti zadovoljena relacija:<br />

ij<br />

2<br />

∂ F<br />

∂ε<br />

∂ε<br />

kl<br />

= C<br />

ijkl<br />

kl<br />

2<br />

∂ F<br />

=<br />

∂ε<br />

∂ε<br />

ij<br />

= C<br />

ij<br />

klij<br />

(4.11)<br />

(4.12)<br />

Posljednja jednadžba daje uvjet simetrije po parovima indeksa (i,j) i (k,l) što konačno dovodi<br />

do tenzora koeficijenata elastičnosti s 21 različitom komponentom kojima je opisano elastično<br />

linearno ponašanje općenito anizotropnog materijala.<br />

Promatra se potpuni tenzor koeficijenata elatičnosti zbog praćenja redoslijeda indeksa pri<br />

transformaciji koordinatnog sustava.<br />

C ε = σ (4.13)<br />

ij<br />

ijkl<br />

Tenzor koeficijenata elastičnosti u ortonormiranom koordinatnom sustavu, može se zapisati u<br />

shemi:<br />

kl<br />

⎡C1111<br />

C1211<br />

C1311<br />

C2111<br />

C2211<br />

C2311<br />

C3111<br />

C3211<br />

C3311<br />

⎤<br />

⎢C<br />

⎥<br />

1112 C1212<br />

C1312<br />

C2112<br />

C2212<br />

C2312<br />

C3112<br />

C3212<br />

C3312<br />

⎢<br />

⎥<br />

⎢C1113<br />

C1213<br />

C1313<br />

C2113<br />

C2213<br />

C2313<br />

C3113<br />

C3213<br />

C3313<br />

⎥<br />

⎢C<br />

⎥<br />

1121 C1221<br />

C1321<br />

C2121<br />

C2221<br />

C2321<br />

C3121<br />

C3221<br />

C3321<br />

⎢<br />

⎥<br />

C ijkl = ⎢C1122<br />

C1222<br />

C1322<br />

C2122<br />

C2222<br />

C2322<br />

C3122<br />

C3222<br />

C3322<br />

⎥ (4.14)<br />

⎢C<br />

⎥<br />

1123 C1223<br />

C1323<br />

C2123<br />

C2223<br />

C2323<br />

C3123<br />

C3223<br />

C3323<br />

⎢<br />

⎥<br />

⎢C1131<br />

C1231<br />

C1331<br />

C2131<br />

C2231<br />

C2331<br />

C3131<br />

C3231<br />

C3331<br />

⎥<br />

⎢C<br />

⎥<br />

1132 C1232<br />

C1332<br />

C2132<br />

C2232<br />

C2332<br />

C3132<br />

C3232<br />

C3332<br />

⎢<br />

⎥<br />

⎣C1133<br />

C1233<br />

C1333<br />

C2133<br />

C2233<br />

C2333<br />

C3133<br />

C3233<br />

C3333<br />

⎦<br />

Tenzori se iz jednog u drugi koordinatni sustav transformiraju po zakonu:<br />

gdje su:<br />

= n<br />

= n<br />

Cijkl, Cij – komponente transformiranog tenzora<br />

C<br />

C<br />

ijkl<br />

ij<br />

iα<br />

n<br />

iα<br />

n<br />

jβ<br />

n<br />

jβ<br />

kγ<br />

C<br />

n<br />

αβ<br />

lδ<br />

C<br />

αβγδ<br />

Cαβγδ, Cαβ – komponente tenzora koje se transformiraju<br />

(4.15)<br />

niα, njβ, nkγ, nlδ – kosinusi kutova što ih zatvaraju koordinatne osi osnovnog i transformiranog<br />

koordinatnog sustava<br />

i, j, k, l – indeksi transformiranog tenzora<br />

α, β, γ, δ – indeksi osnovnog tenzora


4.3 Heksagonalna simetrija (uslojeni materijal)<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 5 -<br />

Crtež 4.1 Model uslojenog materijala<br />

Odabrani model materijala ima istaknuta svojstva u dva meñusobno okomita smjera i to u<br />

ravnini x1,x3 je izotropan što ima za posljedicu da se komponente tenzora ne mijenjaju pri<br />

zaokretanju koordinatnog sustava za proizvoljni kut ϑ oko koordinatne osi x2.<br />

Po pravcu koordinatne osi x2 materijal ima različita svojstva od onih u ravninama uslojenosti,<br />

meñutim komponente tenzora se ne mijenjaju ukoliko se koordinatni sustav zarotira oko osi x1 ili<br />

osi x3 za kut ϑ = nπ, n=1, 2, …, N.<br />

Ova konstatacija se temelji na činjenici da svojstva materijala ostaju ista bez obzira u kojem<br />

koordinatnom sustavu se promatrali, ali tenzor koeficijenata elastičnosti se mijenja promjenom<br />

koordinatnog sustava. Meñutim, u konkretnom slučaju pravac na kojem leži koordinatna os<br />

odreñuje jedno istaknuto svojstvo materijala, te da bi se dobilo to svojstvo i u zarotiranom<br />

sustavu potrebno je da se zarotirana os x2 poklopi s prvobitnim pravcem osi x2, a to ima za<br />

posljedicu da tenzor koeficijenata elastičnosti ostaje nepromijenjen. Koeficijenti transformacije<br />

koordinatnih sustava, a takoñer i komponenata tenzora izraženih u tom koordinatnom sustavu za<br />

slučaj rotacije oko ishodišta ortonormiranog koordinatnog sustava, su kosinusi kutova što ih<br />

zatvaraju koordinatne osi osnovnog i transformiranog koordinatnog sustava.<br />

xi xα x1 x2 x3<br />

x 1 n11 n12 n13<br />

x 2 n21 n22 n23<br />

x 3 n31 n32 n33<br />

(4.16)<br />

Ako se koordinatnim osima lokalnog sustava pridruže indeksi α, β, γ, δ = 1, 2, 3, a osima<br />

transformiranog sustava i, j, k, l = 1, 2, 3, tada se veza izmeñu indeksa i-α, j-β, γ-k, δ-l, itd.<br />

izražava vezama u (4.16).


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 6 -<br />

Za ortotropan materijal u koordinatnom sustavu čije se osi poklapaju s istaknutim svojstvima<br />

materijala, rotacijom koordinatnog sustava oko dvije koordinatne osi za kut ϑ = ± nπ; n=1,2,…,N<br />

tenzor koeficijenata elastičnosti se pojednostavljuje:<br />

⎡C1111<br />

0 0 0 C2211<br />

0 0 0 C3311<br />

⎤<br />

⎢ 0 C<br />

⎥<br />

1212 0 C2112<br />

0 0 0 0 0<br />

⎢<br />

⎥<br />

⎢ 0 0 C1313<br />

0 0 0 C3113<br />

0 0 ⎥<br />

⎢ 0 C<br />

⎥<br />

1221 0 C2121<br />

0 0 0 0 0<br />

⎢<br />

⎥<br />

C ijkl = ⎢C1122<br />

0 0 0 C2222<br />

0 0 0 C3322<br />

⎥ (4.17)<br />

⎢ 0 0 0 0 0 C<br />

⎥<br />

2323 0 C3223<br />

0<br />

⎢<br />

⎥<br />

⎢ 0 0 C1331<br />

0 0 0 C3131<br />

0 0 ⎥<br />

⎢ 0 0 0 0 0 C<br />

⎥<br />

2332 0 C3232<br />

0<br />

⎢<br />

⎥<br />

⎣C1133<br />

0 0 0 C2233<br />

0 0 0 C3333<br />

⎦<br />

Na crtežu 4.1 shematski je prikazan materijal koji zadovoljava heksagonalni uvjet simetrije<br />

oko koordinatne osi x2. Takav materijal je izotropan u ravnini x1,x3, a u smjeru koordinatne osi<br />

x2 ima različita svojstva od onih u preostala dva smjera. Rotacijom koordinatnog sustava oko osi<br />

x2 za proizvoljni kut ϑ komponente tenzora koeficijenata elastičnosti moraju ostati<br />

nepromijenjene.<br />

Dovoljno je poći od tenzora koeficijenata elastičnosti za ortotropan materijal u lokalnom<br />

koordinatnom sustavu (4.17) (istaknuta svojstva materijala se podudaraju s koordinatnim osima).<br />

Za heksagonalnu simetriju oko osi x2, koeficijenti transformacije su:<br />

xi xα x1 x2 x3<br />

x 1 n11 = cosϑ n12 = 0 n13 = sinϑ<br />

x 2 n21 = 0 n22 = 1 n23 = 0<br />

x 3 n31 = −sinϑ n32 = 0 n33 = cosϑ<br />

i,α=1,2,3 (4.18)<br />

Tablica (4.18) odreñuje vrijednost projekcija ortova lokalnog koordinatnog sustava x1,x2,x3 na<br />

koordinatne osi globalnog sustava x 1 ,x 2 ,x 3 (crtež 4.2).


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 7 -<br />

Crtež 4.2.<br />

Obzirom da su koordinatni sustavi ortonormirani, komponente tenzora su jednake i<br />

transformiraju se iz jednog u drugi koordinatni sustav kao i koordinatni vektori (crtež 4.2).<br />

Koristeći zakon transformiranja tenzora (4.15) i koeficijente transformacije (4.18), mogu se<br />

izračunati komponente tenzora (4.17) u koordinatnom sustavu zarotiranom za proizvoljni kut ϑ<br />

oko koordinatne osi x2.<br />

C<br />

Npr. za komponentu C 1111 :<br />

1111<br />

= n<br />

4<br />

11<br />

+ n<br />

+ n<br />

C<br />

2<br />

11<br />

2<br />

12<br />

1111<br />

n<br />

n<br />

2<br />

13<br />

2<br />

13<br />

+ n<br />

C<br />

C<br />

2<br />

11<br />

1313<br />

3232<br />

n<br />

2<br />

12<br />

+ n<br />

+ n<br />

C<br />

2<br />

11<br />

2<br />

11<br />

1122<br />

n<br />

n<br />

2<br />

13<br />

2<br />

13<br />

+ n<br />

C<br />

C<br />

1331<br />

2<br />

11<br />

3311<br />

n<br />

2<br />

13<br />

+ n<br />

+ n<br />

C<br />

2<br />

11<br />

2<br />

13<br />

1133<br />

n<br />

n<br />

2<br />

12<br />

2<br />

12<br />

+ n<br />

C<br />

C<br />

2<br />

11<br />

2112<br />

3322<br />

n<br />

2<br />

12<br />

+ n<br />

+ n<br />

C<br />

2<br />

11<br />

4<br />

13<br />

1212<br />

n<br />

2<br />

12<br />

C<br />

C<br />

3333<br />

+ n<br />

Uvodeći u (4.19) koeficijente transformacije iz (4.18), dobiva se:<br />

1111<br />

4<br />

1111<br />

2<br />

2<br />

2<br />

11<br />

2121<br />

n<br />

2<br />

12<br />

+ n<br />

C<br />

2<br />

12<br />

1221<br />

n<br />

2<br />

13<br />

+<br />

C<br />

3223<br />

4<br />

( C1133<br />

+ C1313<br />

+ C1331<br />

+ C3113<br />

+ C3131<br />

+ C3311)<br />

+ sin C3333<br />

+<br />

(4.19)<br />

C = cos ϑ C + sin ϑcos<br />

ϑ<br />

ϑ (4.20)<br />

Iz uvjeta simetrije tenzora koeficijenata elastičnosti (4.3), (4.4) i (4.12) proizlazi da sljedeći<br />

koeficijenti moraju meñusobno biti jednaki:<br />

C<br />

C<br />

1133<br />

1313<br />

= C<br />

= C<br />

3311<br />

1331<br />

= C<br />

3113<br />

= C<br />

3131<br />

(4.21)<br />

Osim toga, obzirom da je materijal izotropan u ravnini x1,x3, takoñer moraju biti meñusobno<br />

jednaki koeficijenti:<br />

C C = (4.22)<br />

1111<br />

3333<br />

Jednadžba (4.20) uz relacije (4.21) i (4.22) postaje:


1111<br />

ili napisana na drugi način:<br />

1111<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 8 -<br />

4 4<br />

2 2<br />

( sin ϑ + cos ϑ)<br />

C + 2sin<br />

ϑcos<br />

ϑ ( C 2C<br />

)<br />

C =<br />

+<br />

1111<br />

1133<br />

1313<br />

2 2 2<br />

2 2<br />

( sin ϑ + cos ϑ)<br />

C + 2sin<br />

ϑcos<br />

ϑ ( − C + C C )<br />

C = +<br />

(4.23)<br />

1111<br />

Imajući u vidu fizikalni smisao transformacije, nije teško zaključiti iz (4.23) da je C 1111 =<br />

C1111 što ima za posljedicu da se drugi pribrojnik u jednadžbi (4.23) mora poništiti.<br />

Kut rotacije ϑ je proizvoljan pa izraz u zagradama mora biti jednak nuli što daje sljedeću vezu<br />

meñu komponentama:<br />

1111<br />

1133<br />

1313<br />

1313 3131 1<br />

C = C = ( C1111<br />

− C1133<br />

)<br />

(4.24)<br />

2<br />

Ispunjavajući uvjet heksagonalne simetrije oko koordinatne osi x1 i osi x3, tenzor koeficijenata<br />

elastičnosti se pojednostavljuje i opisuje svojstva izotropnog materijala:<br />

1<br />

2<br />

3<br />

1 2 3<br />

1 2 3 1 2 3 1 2 3<br />

1 C 1111 0 0 0 C 1122 0 0 0 C 1122<br />

2 0 (C 1111−C 1122) /2 0 (C 1111−C 1122) /2 0 0 0 0 0<br />

3 0 0 (C 1111−C 1122) /2 0 0 0 (C 1111−C 1122) /2 0 0<br />

1 0 (C 1111−C 1122) /2 0 (C 1111−C 1122) /2 0 0 0 0 0<br />

2 C 1122 0 0 0 C 1111 0 0 0 C 1122<br />

3 0 0 0 0 0 (C 1111−C 1122) /2 0 (C 1111−C 1122) /2 0<br />

1 0 0 (C 1111−C 1122) /2 0 0 0 (C 1111−C 1122) /2 0 0<br />

2 0 0 0 0 0 (C 1111−C 1122) /2 0 (C 1111−C 1122) /2 0<br />

3 C 1122 0 0 0 C 1122 0 0 0 C 1111<br />

ili u indeksnom zapisu:<br />

C = α δ δ + β δ δ + γ δ δ<br />

gdje su α, β i γ poznate konstante za odreñeni materijal.<br />

ijkl<br />

4.4 Laméove konstante elastičnosti i opći Hookov zakon<br />

ij<br />

kl<br />

ik<br />

jl<br />

il<br />

jk<br />

(4.25)<br />

Uvodeći Lamé-ove konstante λ = C1122 i μ = 1/2(C1111−C1122) u (4.25), dobiva se tenzor<br />

koeficijenata elastičnosti sa samo dvije neovisne komponente jer je i C1111 = λ + 2μ , pa veza<br />

izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> ima sljedeći oblik:


⎧σ<br />

⎪<br />

⎪σ<br />

⎪<br />

σ<br />

⎪<br />

⎪σ<br />

⎪<br />

⎨σ<br />

⎪<br />

⎪σ<br />

⎪σ<br />

⎪<br />

⎪σ<br />

⎪<br />

⎪⎩<br />

σ<br />

11<br />

21<br />

31<br />

12<br />

22<br />

32<br />

13<br />

23<br />

33<br />

⎫<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪⎭<br />

=<br />

⎡λ<br />

+ 2μ<br />

⎢<br />

⎢ 0<br />

⎢<br />

0<br />

⎢<br />

⎢ 0<br />

⎢<br />

⎢ λ<br />

⎢<br />

⎢ 0<br />

⎢<br />

⎢<br />

0<br />

⎢ 0<br />

⎢<br />

⎢⎣<br />

λ<br />

0<br />

μ<br />

0<br />

μ<br />

0<br />

0<br />

0<br />

0<br />

0<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 9 -<br />

0<br />

0<br />

μ<br />

0<br />

0<br />

0<br />

μ<br />

0<br />

0<br />

0<br />

μ<br />

0<br />

μ<br />

0<br />

0<br />

0<br />

0<br />

0<br />

λ<br />

0<br />

0<br />

0<br />

λ +<br />

0<br />

0<br />

0<br />

λ<br />

2μ<br />

0<br />

0<br />

0<br />

0<br />

0<br />

μ<br />

0<br />

μ<br />

0<br />

0<br />

0<br />

μ<br />

0<br />

0<br />

0<br />

μ<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

μ<br />

0<br />

μ<br />

0<br />

λ ⎤ ⎧ε<br />

⎥ ⎪<br />

0 ⎥ ⎪ε<br />

0<br />

⎥ ⎪<br />

ε<br />

⎥ ⎪<br />

0 ⎥ ⎪ε<br />

⎥ ⎪<br />

λ ⎥ ⋅ ⎨ε<br />

⎥ ⎪<br />

0 ⎥ ⎪ε<br />

⎥<br />

0 ⎪<br />

⎥<br />

ε<br />

⎪<br />

0 ⎥ ⎪ε<br />

⎥ ⎪<br />

λ + 2μ⎥⎦<br />

⎪⎩<br />

ε<br />

11<br />

21<br />

31<br />

12<br />

22<br />

32<br />

13<br />

23<br />

33<br />

⎫<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪⎭<br />

(4.26)<br />

Analizirajući jednoosno stanje <strong>naprezanja</strong>, npr. u koordinatnom smjeru x1, sve komponente<br />

tenzora <strong>naprezanja</strong> se poništavaju osim σ11, a ono je uz (4.25) i (4.26):<br />

( λ + μ)<br />

ε11<br />

+ λε22<br />

+ 33<br />

σ 11 = 2 λε<br />

(4.27a)<br />

Takoñer su mješovite komponente tenzora deformacije jednake nuli, pa se mogu napisati još<br />

samo dvije jednadžbe za smjerove x2 i x3:<br />

σ<br />

σ<br />

22<br />

33<br />

= 0 =<br />

= 0 =<br />

( λ + 2μ)<br />

ε22<br />

+ λε11<br />

+ λε33<br />

( λ + 2μ)<br />

ε33<br />

+ λε11<br />

+ λε22<br />

Iz (4.27a) i (4.27b) se mogu izraziti komponente <strong>deformacija</strong>:<br />

( ) ( ) 11<br />

33 22<br />

11<br />

3λ<br />

+ 2μ<br />

2 λ + μ<br />

(4.27b)<br />

11<br />

λ + μ<br />

λ<br />

ε = σ ; ε = ε = − ε<br />

(4.28)<br />

μ<br />

Ako se u izrazima (4.28) koeficijent uz σ11 zamijeni s 1/E, a uz ε11 s ν, dobiva se veza modula<br />

elastičnosti E i Poissonovog koeficijenta ν s Laméovim konstantama μ i λ:<br />

ili:<br />

( 3λ<br />

+ 2μ)<br />

μ<br />

λ<br />

E =<br />

; ν =<br />

(4.29)<br />

λ + μ<br />

2<br />

λ =<br />

ν E<br />

;<br />

μ =<br />

( λ + μ)<br />

( 1+<br />

ν)<br />

( 1−<br />

2ν)<br />

2 ( 1+<br />

ν)<br />

E<br />

(4.30)<br />

Za realan materijal Laméove konstante su pozitivne iz čega slijede uvjeti za veličinu<br />

Poissonovog koeficijenta:<br />

iz:<br />

λ<br />

ν = → ν > 0 , a iz<br />

2<br />

( λ + μ)<br />

μ =<br />

2<br />

Da bi λ bio pozitivan u (4.30) mora biti zadovoljena relacija:<br />

E<br />

( 1+<br />

ν)<br />

→<br />

ν < 1


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 10 -<br />

1<br />

1 − 2ν<br />

> 0 → 0 < ν <<br />

(4.31)<br />

2<br />

Poissonov koeficijent veličine 0.5 odgovara idealno nestišljivoj tekućini.<br />

Uvažavajući svojstvo simetričnosti tenzora <strong>naprezanja</strong> i tenzora <strong>deformacija</strong>, šest je<br />

komponenata tenzora <strong>naprezanja</strong> izraženo preko tenzora <strong>deformacija</strong> kod čega je γ ij = εij<br />

+ ε ji .<br />

Iz izraza (4.26), uz zamjenu Lamé-ovih konstanti elastičnim konstantama materijala E i ν,<br />

dobiva se opći Hookov zakon:<br />

⎧σ<br />

⎪<br />

⎪σ<br />

⎪<br />

⎪<br />

⎪σ<br />

⎨<br />

⎪σ<br />

⎪<br />

⎪σ<br />

⎪<br />

⎪<br />

⎩σ<br />

11<br />

22<br />

33<br />

12<br />

23<br />

31<br />

⎫<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎭<br />

=<br />

ili u indeksnom zapisu:<br />

⎧ε<br />

⎪<br />

⎪ε<br />

⎪<br />

⎪<br />

⎪ε<br />

⎨<br />

⎪γ<br />

⎪<br />

⎪<br />

γ<br />

⎪<br />

⎪<br />

⎩γ<br />

11<br />

22<br />

33<br />

12<br />

23<br />

31<br />

⎫<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎭<br />

=<br />

⎡ 1−<br />

ν<br />

⎢<br />

⎢ ν<br />

⎢<br />

⎢<br />

⎢ ν<br />

E<br />

⎢<br />

( 1+<br />

ν)<br />

( 1−<br />

2ν)<br />

⎢ 0<br />

⎢<br />

⎢<br />

⎢<br />

0<br />

⎢<br />

⎢ 0<br />

⎣<br />

1<br />

E<br />

⎡<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢<br />

⎢⎣<br />

ili u indeksnom zapisu:<br />

1<br />

− ν<br />

− ν<br />

0<br />

0<br />

0<br />

− ν<br />

1<br />

− ν<br />

0<br />

0<br />

0<br />

ij<br />

ν<br />

1−<br />

ν<br />

ν<br />

0<br />

0<br />

0<br />

kk<br />

ij<br />

ν<br />

ν<br />

1−<br />

ν<br />

0<br />

0<br />

0<br />

σ = λ ε δ + 2 μ e ;<br />

− ν<br />

− ν<br />

1<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

1−<br />

2ν<br />

2<br />

ij<br />

2(<br />

1+<br />

ν)<br />

1 ⎛ λ ⎞<br />

εij = ⎜σij<br />

− σkk<br />

δij<br />

⎟<br />

2μ<br />

⎝ 3λ<br />

+ 2μ<br />

⎠<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

0<br />

1−<br />

2ν<br />

2<br />

0<br />

2(<br />

1+<br />

ν)<br />

0 ⎤ ⎧ε<br />

⎥ ⎪<br />

0 ⎥ ⎪<br />

⎥<br />

ε<br />

⎪<br />

⎥<br />

0<br />

⎪<br />

⎥ ⎪ε<br />

⎥ ⋅ ⎨<br />

0 ⎥ ⎪γ<br />

⎥ ⎪<br />

⎥<br />

0 ⎪<br />

⎥ γ<br />

⎪<br />

1−<br />

2ν<br />

⎥ ⎪<br />

⎥<br />

2<br />

⎪<br />

⎦ ⎩γ<br />

11<br />

22<br />

33<br />

12<br />

23<br />

31<br />

⎫<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎭<br />

0 ⎤ ⎧σ11<br />

⎫<br />

⎥ ⎪ ⎪<br />

0 ⎥ ⎪σ<br />

⎪<br />

22<br />

⎥ ⎪ ⎪<br />

⎥ ⎪ ⎪<br />

0 ⎥ ⎪σ33<br />

⎪<br />

⎥ ⋅ ⎨ ⎬<br />

0 ⎥ ⎪σ12<br />

⎪<br />

⎥ ⎪ ⎪<br />

⎥<br />

0<br />

⎪<br />

σ<br />

⎪<br />

⎥ ⎪ 23⎪<br />

⎥ ⎪ ⎪<br />

2(<br />

1+<br />

ν)<br />

⎥⎦<br />

⎪⎩<br />

σ31⎪⎭<br />

(4.32)<br />

(4.33)<br />

Budući da jednadžbe (4.32), odnosno jednadžbe (4.33), ne sadrže nikakve parcijalne<br />

derivacije, oblik tih jednadžbi se pri transformaciji iz Descartesovog koordinatnog sustava u<br />

proizvoljni ortogonalni koordinatni sustav ne mijenja.<br />

U skladu s tim, za cilindrični koordinatni sustav vrijedi:


σ<br />

σ<br />

σ<br />

τ<br />

r<br />

ϕ<br />

z<br />

rϕ<br />

=<br />

=<br />

=<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 11 -<br />

( 1+<br />

ν)(<br />

1−<br />

2ν)<br />

[ ( 1−<br />

ν)<br />

ε + ν ( ε + ε ) ]<br />

E [ ( 1−<br />

ν)<br />

εϕ<br />

+ ν ( εz<br />

+ εr<br />

) ]<br />

( 1+<br />

ν)(<br />

1−<br />

2ν)<br />

( 1+<br />

ν)(<br />

1−<br />

2ν)<br />

= G γ<br />

rϕ<br />

E<br />

E<br />

,<br />

τ<br />

ϕz<br />

[ ( 1−<br />

ν)<br />

ε + ν ( ε + ε ) ]<br />

= G γ<br />

ϕz<br />

r<br />

z<br />

,<br />

τ<br />

ϕ<br />

r<br />

zr<br />

z<br />

ϕ<br />

= G γ<br />

zr<br />

⎫<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎭<br />

(4.34)<br />

dok je eksplicitni oblik jednadžbi (4.33) analogan jednadžbama u pravokutnom koordinatnom<br />

susatvu i za cilindrični koordinatni sustav glasi:<br />

ε<br />

ε<br />

ε<br />

γ<br />

ϕ<br />

=<br />

rϕ<br />

=<br />

=<br />

=<br />

[ σ − ν(<br />

σ + σ ) ]<br />

[ σ − ν(<br />

σ + σ ) ]<br />

[ σ − ν(<br />

σ + σ ) ]<br />

τ<br />

ϕ<br />

rϕ<br />

γ<br />

ϕ<br />

ϕz<br />

=<br />

ϕ<br />

τ<br />

ϕz<br />

4.5 Fizikalno značenje modula elastičnosti<br />

r<br />

z<br />

1<br />

E<br />

1<br />

E<br />

1<br />

E<br />

1<br />

G<br />

r<br />

z<br />

,<br />

r<br />

z<br />

z<br />

1<br />

G<br />

r<br />

,<br />

γ<br />

zr<br />

=<br />

1<br />

G<br />

τ<br />

zr<br />

⎫<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎬<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎭<br />

(4.35)<br />

Kod slučaja izotropije, prethodno definiran modul elastičnosti ima jednostavno fizikalno<br />

značenje. Na crtežu 4.3 prikazan je ''radni'' dijagram jednoosnog vlačnog <strong>naprezanja</strong> za tri<br />

različite vrste materijala. Fizikalne karakteristike materijala E, G ili ν odreñuju se praćenjem<br />

pojedinih stanja <strong>naprezanja</strong>, crtež 4.4, korištenjem standardne laboratorijske tehnike.<br />

σ<br />

Čelik<br />

Lijevano željezo<br />

Aluminij<br />

Crtež 4.3 σ−ε dijagrami za čelik, lijevano željezo i aluminij<br />

ε


4.5.1 Jednoosno naprezanje<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 12 -<br />

Razmotrimo test jednoosnog <strong>naprezanja</strong> na uzorku podvrgnutom vlaku u smjeru osi x (vidi<br />

crtež 4.4). Stanje <strong>naprezanja</strong> potpuno je opisano sljedećim tenzorom <strong>naprezanja</strong>:<br />

σ<br />

ij<br />

⎡σ<br />

=<br />

⎢<br />

⎢<br />

0<br />

⎢⎣<br />

0<br />

1+<br />

ν ν<br />

Koristeći relaciju ε ij = σij<br />

− σkkδij<br />

, dobiva se pripadajući tenzor <strong>deformacija</strong>:<br />

E E<br />

σ<br />

σ<br />

F<br />

F<br />

ε<br />

ij<br />

⎡σ<br />

⎢E<br />

⎢<br />

= ⎢ 0<br />

⎢<br />

⎢ 0<br />

⎣<br />

τ<br />

0<br />

0<br />

0<br />

0<br />

ν<br />

− σ<br />

E<br />

0<br />

τ<br />

τ<br />

Mt<br />

τ<br />

Mt<br />

0⎤<br />

0<br />

⎥<br />

⎥<br />

0⎥⎦<br />

0<br />

⎤<br />

⎥<br />

⎥<br />

0 ⎥<br />

ν ⎥<br />

− σ⎥<br />

E ⎦<br />

Crtež 4.4 Posebni slučajevi stanja <strong>naprezanja</strong><br />

Zbog toga je E = σ / εx<br />

i predstavlja tangens kuta odnosno nagib tangente na krivulju<br />

naprezanje-<strong>deformacija</strong>, dok je ν = −εy<br />

/ εx<br />

= −εz<br />

/ εx<br />

odnos poprečne i uzdužne deformacije.<br />

Standardni sustavi mjerenja lako mogu sakupiti podatke o aksijalnom naprezanju i poprečnoj ili<br />

uzdužnoj deformaciji, i preko ovog jednog testa odrediti obje konstante za materijal koji nas<br />

zanima.<br />

4.5.2 Čisti posmik<br />

Ako je tanki puni cilindar podvrgnut torzijskom opterećenju (kako je prikazano na crtežu 4.4),<br />

stanje <strong>naprezanja</strong> na površini cilindričnog uzorka je odreñeno tenzorom <strong>naprezanja</strong>:<br />

σ<br />

ij<br />

⎡0<br />

=<br />

⎢<br />

⎢<br />

τ<br />

⎢⎣<br />

0<br />

τ<br />

0<br />

0<br />

0⎤<br />

0<br />

⎥<br />

⎥<br />

0⎥⎦<br />

p<br />

p<br />

p


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 13 -<br />

Koristeći Hookeov zakon, izračuna se pripadajući tenzor <strong>deformacija</strong>:<br />

ε<br />

ij<br />

⎡ 0<br />

=<br />

⎢<br />

⎢<br />

τ / 2μ<br />

⎢⎣<br />

0<br />

τ / 2μ<br />

Dakle modul posmika je odreñen s μ = τ / 2εxy<br />

= τ / γxy<br />

, i ovaj je modul jednostavno nagib<br />

krivulje posmično naprezanje – posmična <strong>deformacija</strong>.<br />

4.5.3 Hidrostatički tlak<br />

Posljednji primjer je vezan s jednolikim tlakom (ili naprezanjem) kockastog uzorka, kako je<br />

prikazano na crtežu 4.4. Ovaj se tip pokusa izvodi postavljanjem uzorka u visoko tlačnu komoru.<br />

Stanje <strong>naprezanja</strong> u ovom slučaju je odreñeno sa:<br />

σ<br />

ij<br />

⎡−<br />

p<br />

=<br />

⎢<br />

⎢<br />

0<br />

⎢⎣<br />

0<br />

0<br />

− p<br />

0<br />

0<br />

0<br />

0⎤<br />

0<br />

⎥<br />

⎥<br />

0⎥⎦<br />

0 ⎤<br />

0<br />

⎥<br />

⎥<br />

= −p<br />

δ<br />

− p⎥⎦<br />

Ovo je izotropno stanje <strong>naprezanja</strong> i deformacije slijede iz Hookeovog zakona:<br />

ε<br />

ij<br />

⎡ 1−<br />

2ν<br />

− p<br />

⎢ E<br />

⎢<br />

= ⎢ 0<br />

⎢<br />

⎢<br />

0<br />

⎢⎣<br />

0<br />

1−<br />

2ν<br />

− p<br />

E<br />

0<br />

ij<br />

0<br />

⎤<br />

⎥<br />

⎥<br />

0 ⎥<br />

⎥<br />

1−<br />

2ν<br />

− p<br />

⎥<br />

E ⎥⎦<br />

Deformacija koja predstavlja promjenu volumena materijala je definirana sa<br />

ϑ = ε = −3<br />

1−<br />

2ν<br />

p / što se može zapisati kao:<br />

kk<br />

( ) E<br />

gdje se k E / [ 3(<br />

1−<br />

2ν)<br />

]<br />

p = −k<br />

ϑ<br />

(4.36)<br />

= naziva zapreminski modul elastičnosti. Ova dodatna elastična konstanta<br />

predstavlja omjer tlaka i deformacije i može se odnositi kao volumna krutost materijala.<br />

Primijetimo da kada se Poissonov koeficijent bliži vrijednosti 0.5, zapreminski modul postaje<br />

beskonačan a materijal nema nikakvu volumnu deformaciju i praktično postaje nestlačiv.<br />

Ova diskusija o elastičnim modulima za izotropne materijale dovela nas je do definicije pet<br />

konstanti λ, µ ( ili G)<br />

, E, ν i k . Meñutim, treba imati na umu da su samo dvije od njih potrebne za<br />

opis karakteristika materijala. Iako smo već razvili neke veze izmeñu pojedinih modula, može se<br />

pronaći još mnogo drugih relacija. U stvari, može se pokazati da je svih pet konstanti u<br />

meñusobnom odnosu, i ako su dane bilo koje dvije od njih preostale tri se mogu odrediti<br />

koristeći jednostavne formule. Rezultati ovih veza zgodno su prikazani u tablici 4.1.<br />

Tipične nominalne vrijednosti elastičnih konstanti za pojedinu vrstu materijala dane su u<br />

tablici 4.2. Ovi moduli predstavljaju srednje vrijednosti, i za pojedine materijale može doći do<br />

stanovitih varijacija.<br />

U prethodnim poglavljima je naglašavano kako cilindrične koordinate imaju široku primjenu<br />

kod mnogih problema.


E,ν<br />

E,<br />

k<br />

E,<br />

E,<br />

μ<br />

λ<br />

ν,<br />

k<br />

ν,<br />

μ<br />

ν,<br />

λ<br />

k,<br />

k,<br />

μ<br />

λ<br />

ν,<br />

λ<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 14 -<br />

Tablica 4.1 <strong>Veza</strong> izmeñu konstanti koje definiraju karakteristike materijala<br />

Ε<br />

E<br />

E<br />

E<br />

E<br />

3k(<br />

1−<br />

2ν)<br />

2μ(<br />

1+<br />

ν)<br />

λ(<br />

1+<br />

ν)(<br />

1−<br />

2ν)<br />

ν<br />

9kμ<br />

6k<br />

+ μ<br />

9k(<br />

k − λ)<br />

3k<br />

− λ<br />

μ(<br />

3λ<br />

+ 2μ)<br />

k − μ<br />

ν<br />

ν<br />

3k<br />

− E<br />

6k<br />

E − 2μ<br />

2μ<br />

2λ<br />

E + μ + R<br />

ν<br />

ν<br />

ν<br />

3k<br />

− 2μ<br />

6k<br />

+ 2μ<br />

λ<br />

3k<br />

− λ<br />

λ<br />

2(<br />

λ + μ)<br />

k<br />

E<br />

3(<br />

1−<br />

2ν)<br />

k<br />

μE<br />

3(<br />

3μ<br />

− E)<br />

E + 3λ<br />

+ R<br />

6<br />

k<br />

2μ(<br />

1+<br />

ν)<br />

3(<br />

1−<br />

2ν)<br />

λ(<br />

1+<br />

ν)<br />

3ν<br />

k<br />

k<br />

3λ<br />

+ 2μ<br />

3<br />

µ<br />

E<br />

2(<br />

1+<br />

ν)<br />

3kE<br />

9k<br />

− E<br />

μ<br />

E − 3λ<br />

+ R<br />

4<br />

3k(<br />

1−<br />

2ν)<br />

2(<br />

1+<br />

ν)<br />

μ<br />

λ(<br />

1−<br />

2ν)<br />

2ν<br />

μ<br />

3<br />

( k − λ)<br />

2<br />

μ<br />

λ<br />

Eν<br />

( 1+<br />

ν)(<br />

1−<br />

2ν)<br />

3k(<br />

3k<br />

− E)<br />

9k<br />

− E<br />

μ(<br />

E − 2μ)<br />

3μ<br />

− E<br />

λ<br />

3kν<br />

1+<br />

ν<br />

2μν<br />

1−<br />

2ν<br />

λ<br />

2<br />

k − μ<br />

3<br />

Tablica 4.2 Tipične vrijednosti elastičnih konstanti za česte inženjerske materijale<br />

E(GPa) ν μ (GPa) λ (GPa) k (GPa) α (10 -6 /°C)<br />

Aluminij 68.9 0.34 25.7 54.6 71.8 25.5<br />

Beton 27.6 0.20 11.5 7.7 15.3 11.1<br />

Bakar 89.6 0.34 33.4 71 93.3 18<br />

Staklo 68.9 0.25 27.6 27.6 45.9 8.8<br />

Plastika 28.3 0.40 10.1 4.4 47.2 11.4<br />

Guma 0.0019 0.499 0.654*10 -3 0.326 0.326 18.7<br />

Čelik 207 0.29 80.2 20.4 164 13.5<br />

λ<br />

λ


4.5.4 Zakon ponašanja u termoelastičnosti<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 15 -<br />

Dobro je poznato da promjena temperature kod slobodnog elastičnog tijela izaziva<br />

deformiranje. Dakle, opće polje <strong>deformacija</strong> nastaje usljed mehaničkih i toplinskih djelovanja. U<br />

smislu teorije malih <strong>deformacija</strong>, ukupna <strong>deformacija</strong> može se rastaviti na mehaničku i toplinsku<br />

komponentu:<br />

ij<br />

( M)<br />

ij<br />

( T)<br />

ij<br />

ε = ε + ε<br />

(4.37)<br />

Ako je T0 uzeta kao referentna temperatura i T kao proizvoljna temperatura, toplinske<br />

deformacije kod slobodnog tijela mogu se napisati u linearnom obliku:<br />

( T)<br />

ij<br />

ε = α T − T )<br />

(4.38)<br />

ij(<br />

0<br />

gdje je α ij tenzor koeficijenata toplinskog širenja. Primijetimo da temperaturna razlika uzrokuje<br />

toplinske deformacije. Ako je materijal uzet kao izotropan, tada α ij mora biti dijagonalni tenzor<br />

drugog reda i izraz (4.38) postaje:<br />

( T)<br />

ij<br />

ε = α ( T − T ) δ<br />

(4.39)<br />

gdje je α konstanta materijala i naziva se koeficijent toplinskog širenja. Tablica 4.2 daje tipične<br />

vrijednosti za konstante nekih vrsta materijala. Primijetimo da kod izotropnih materijala<br />

temperaturna promjena neće uzrokovati posmična <strong>naprezanja</strong>. Koristeći (4.37) u kombinaciji s<br />

mehaničkim relacijama, dobiva se:<br />

0<br />

ij<br />

1+<br />

ν ν<br />

ε ij = σij<br />

− σkk<br />

δij<br />

+ α(<br />

T − T0<br />

) δij<br />

(4.40)<br />

E E<br />

Pripadajući rezultati za <strong>naprezanja</strong> izražena preko <strong>deformacija</strong> mogu se zapisati kao:<br />

σ = ε −β<br />

( T − T )<br />

(4.41)<br />

ij<br />

Cijkl kl ij 0<br />

gdje je β ij tenzor drugog reda koji sadrži toplinske koeficijente. Ovaj rezultat se ponekad naziva<br />

Duhamel-Neumannov zakon ponašanja termoelastičnosti. Izotropan slučaj se lako dobije samo<br />

invertirajući izraz (4.40) kako bi se dobilo:<br />

σ = λ ε δ + 2μ ε − ( 3λ<br />

+ 2μ)<br />

α ( T − T ) δ<br />

(4.42)<br />

ij<br />

kk<br />

ij<br />

ij<br />

Nakon što smo razvili potrebnih šest konstitutivnih relacija, sustav jednadžbi elastičnosti je<br />

sada upotpunjen s petnaest jednadžbi (<strong>deformacija</strong>-pomak, ravnoteža, Hookeov zakon) s petnaest<br />

nepoznanica (pomaci, deformacije, <strong>naprezanja</strong>). Kako bi se riješili pojedini problemi u<br />

inženjerstvu očito su potrebna još dodatna pojednostavljenja.<br />

0<br />

ij


4.6 Potencijalna energija<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 16 -<br />

Element elastične sredine deformira se pod djelovanjem normalnih i posmičnih <strong>naprezanja</strong>.<br />

Uz pretpostavku da su deformacije proporcionalne naprezanjima, rad koji se obavi prilikom<br />

deformiranja iznosi općenito 1 2 σ ⋅ u ⋅ dA .<br />

σx<br />

u2<br />

dx<br />

u1<br />

σx<br />

Crtež 4.5 Promjena volumena i oblika uslijed deformiranja<br />

gdje je σ naprezanje, u pomak u smjeru <strong>naprezanja</strong>, a dA površina na koju djeluje naprezanje σ.<br />

Ako se, nadalje, pretpostavi da se naprezanje mijenja od točke do točke, a isto tako da su pomaci<br />

dviju pobočaka površine dy⋅dx različiti, stanje <strong>naprezanja</strong> i <strong>deformacija</strong> može se označiti kako je<br />

prikazano na crtežu 4.6.<br />

dy<br />

σx<br />

u<br />

dx<br />

u+du<br />

σx<br />

σx+dσx<br />

Crtež 4.6 Stanje <strong>deformacija</strong> i <strong>naprezanja</strong> na diferencijalnom elementu<br />

Rad na tom elementu diferencijalne veličine zbog <strong>naprezanja</strong> σx može se ovako izraziti:<br />

1<br />

=<br />

2<br />

dz<br />

τzy<br />

τzx<br />

[ ( σ + dσ<br />

)( u + du)<br />

− σ u]<br />

dydz<br />

dWx x x<br />

x<br />

Negativan predznak drugog člana dolazi zbog toga što je vektor <strong>naprezanja</strong> suprotnog smjera<br />

od vektora pomaka. Ako se zanemari mala veličina višeg reda, taj se izraz može svesti na ovaj<br />

oblik:<br />

1<br />

=<br />

2<br />

pošto su σx i u funkcije koordinate x.<br />

1 d<br />

2 dx<br />

( σ du + dσ<br />

u)<br />

dydz<br />

= ( σ ⋅ u)<br />

dx dydz<br />

dWx x x<br />

x<br />

No iste pobočke zbog posmičnog <strong>naprezanja</strong> bočno se pomiču i ako su njihovi pomaci v i w u<br />

smjeru osi y i z, dobije se ukupni rad na pomacima pobočki površine dy⋅dz.<br />

( σ u + τ v + τ w)<br />

dx dydz<br />

1 d<br />

= x xy<br />

.<br />

2 dx<br />

dWx zx<br />

v<br />

σx<br />

w


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 17 -<br />

Uzme li se u obzir i rad na dva ostala para pobočaka elementarne prizme, dobije se:<br />

∂<br />

∂<br />

⎤<br />

( σ u + τ v + τ w)<br />

+ ( σ v + τ w + τ u)<br />

+ ( σ w + τ u + τ v)<br />

dx dydz<br />

1 ⎡ ∂<br />

x = x xy zx<br />

y yz xy<br />

z zx<br />

2 ⎢<br />

⎣∂x<br />

∂y<br />

∂z<br />

⎥<br />

.<br />

⎦<br />

dW yz<br />

Ukoliko postoje zapreminske sile, rad koji one obave iznosi:<br />

Ukupni rad je prema tome:<br />

dW =<br />

1<br />

2<br />

⎡<br />

⎢<br />

σ<br />

⎣<br />

x<br />

⎛ ∂σ<br />

+ u<br />

⎜<br />

⎝ ∂x<br />

⎛ ∂τ<br />

+ w<br />

⎜<br />

⎝ ∂x<br />

∂u<br />

+ τ<br />

∂x<br />

x<br />

zx<br />

xy<br />

∂τ<br />

+<br />

∂y<br />

xy<br />

∂τ<br />

+<br />

∂y<br />

zy<br />

∂v<br />

+ τ<br />

∂x<br />

∂τ<br />

+<br />

∂z<br />

xz<br />

∂σ<br />

+<br />

∂z<br />

xz<br />

z<br />

1<br />

2<br />

∂w<br />

+ σ<br />

∂x<br />

⎞⎤<br />

+ Z<br />

⎟<br />

⎟⎥<br />

⎠⎦<br />

( X u + Y ⋅ v + Z ⋅ w)<br />

dx dydz<br />

⎞ ⎛ ∂τ<br />

+ X<br />

⎟ + v<br />

⎜<br />

⎠ ⎝ ∂x<br />

⋅ .<br />

y<br />

∂v<br />

+ τ<br />

∂y<br />

xy<br />

dx dydz<br />

yz<br />

∂σ<br />

+<br />

∂y<br />

∂w<br />

+ τ<br />

∂y<br />

y<br />

yx<br />

∂τ<br />

+<br />

∂z<br />

yz<br />

∂u<br />

+ σ<br />

∂y<br />

z<br />

⎞<br />

+ Y<br />

⎟ +<br />

⎠<br />

∂w<br />

+ τ<br />

∂z<br />

zx<br />

∂u<br />

+ τ<br />

∂z<br />

zy<br />

∂v<br />

+<br />

∂z<br />

(4.43)<br />

Izrazi u okruglim zagradama jednaki su nuli jer predstavljaju uvjete ravnoteže, pa ostaje za<br />

čitavo tijelo:<br />

∫∫∫ [ ( σxε<br />

x + σyε<br />

y + σzεz<br />

) + 2(<br />

τxyεxy<br />

+ τyzε<br />

yz + τzxε<br />

) ] dx dydz<br />

1<br />

= (4.44)<br />

2<br />

W zx<br />

Rad vanjskih sila praktički potpuno prelazi u potencijalnu energiju elastičnog tijela. Promjena<br />

temperature je naznatna, a isto tako su neznatne ostale promjene pri deformiranju tijela, pa se<br />

njihov udio zanemaruje:<br />

Ep = W.<br />

Komponente <strong>naprezanja</strong>, meñutim, mogu se izraziti komponentama <strong>deformacija</strong> i obratno, pa<br />

se jedanput nakon uvrštenja prije izvedenih izraza za σx, ... τxy ... dobije:<br />

E<br />

ili, ako se λ izrazi sa μ,<br />

2<br />

2<br />

2<br />

[ λε<br />

+ 2μ<br />

( ε + ε + ε ) + 4μ(<br />

ε + ε + ) ] dx dydz<br />

1<br />

2 2 2 2<br />

= W =<br />

v x y z xy yz zx<br />

(4.45)<br />

2<br />

p ∫∫∫ ε<br />

E<br />

2ν<br />

λ =<br />

1−<br />

2ν<br />

⋅μ<br />

⎡ 1−<br />

ν 2 2⎤<br />

= 2μ<br />

⎢ e1<br />

− e2<br />

dx dydz<br />

(4.45a)<br />

⎣2<br />

( 1−<br />

2ν)<br />

p ∫∫∫ ⎥ ⎦<br />

Drugi put, zamjenom komponenata <strong>deformacija</strong> s komponentama <strong>naprezanja</strong>, dobije se:<br />

E<br />

p<br />

1<br />

W<br />

⎡ 2 2 2 ν 2 2 2 2<br />

= = ( σx<br />

+ σy<br />

+ σz<br />

) − I1<br />

+ 2(<br />

τxy<br />

+ τyz<br />

+ τzx<br />

)] dx dydz<br />

4μ<br />

∫∫∫ (4.46)<br />

⎢⎣<br />

1+<br />

ν


E<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 18 -<br />

1 ⎡ 1 2 2⎤<br />

= W =<br />

I1<br />

I2<br />

dx dydz<br />

2 ⎢ −<br />

(4.46a)<br />

μ ⎣2<br />

( 1−<br />

ν)<br />

p ∫∫∫ ⎥ ⎦<br />

U gornjim izrazima, e1 i e2 su invarijante tenzora <strong>deformacija</strong>, a I1 i I2 invarijante tenzora<br />

<strong>naprezanja</strong>. Tim izrazima potvrñuje se tvrdnja da je potencijalna energija invarijantna veličina.<br />

Izjednače li se izrazi (4.45a) i (4.46a), mogu se i invarijante drugog stupnja meñusobno izraziti,<br />

naime:<br />

2 ⎡ 1−<br />

ν 2 2⎤<br />

1 2<br />

4μ ⎢ e1<br />

− e2<br />

= I1<br />

− I2<br />

2(<br />

1 2 ) ⎥<br />

⎣ − ν ⎦ 2(<br />

1+<br />

ν)<br />

ili<br />

Zbroje li se izrazi za normalne komponente <strong>naprezanja</strong>, dobije se:<br />

σx + σy + σz = I1 = (3λ + 2μ) εv = (3λ + 2μ) e1<br />

1−<br />

2ν<br />

1 = I . (4.47)<br />

2(<br />

1−<br />

2ν)μ<br />

e 1<br />

Uvrsti li se taj izraz u jednadžbu, dobije se:<br />

2<br />

2<br />

2<br />

ν(<br />

2 − ν)<br />

2 ( 1+<br />

ν)<br />

2<br />

1<br />

4μ<br />

e = I − I . (4.48)<br />

Izrazi za rad vanjskih sila, odnosno potencijalnu energiju deformiranog tijela (4.44), (4.45) ili<br />

(4.46), predstavljaju ukupni rad. Izdvoji li se, meñutim, iz deformiranog tijela element jedinične<br />

zapremine koji se pod djelovanjem komponenata <strong>naprezanja</strong> σx, ... , τzx deformira i komponente<br />

<strong>deformacija</strong> su εx, ... εxy, ... , onda je njegova potencijalna energija:<br />

[ σ ε + σ ε + σ ε + 2 ( τ ε + τ ε + τ ) ]<br />

1<br />

ep = w = x x y y z z xy xy yz yz zxεzx<br />

.<br />

2<br />

Ta se potencijalna energija ili rad naziva specifičnom potencijalnom energijom ili specifičnim<br />

radom. Ona se opet može izraziti ili samo komponentama <strong>deformacija</strong>:<br />

ili samo komponentama <strong>naprezanja</strong>:<br />

⎡ 1−<br />

ν 2 ⎤<br />

p = w = 2μ<br />

⎢ e1<br />

− e ,<br />

2(<br />

1 2 ) ⎥<br />

⎣ − ν ⎦<br />

e 2<br />

1 ⎡ 1 2 ⎤<br />

p = w =<br />

I1<br />

I .<br />

2 ⎢ −<br />

μ 2(<br />

1 ) ⎥<br />

⎣ + ν ⎦<br />

e 2<br />

Deriviranjem izraza za specifičnu potencijalnu energiju po komponentama <strong>naprezanja</strong> dobije<br />

se:<br />

∂ep<br />

1 ⎡ 2I1<br />

∂I1<br />

∂I2<br />

⎤ 1<br />

=<br />

[ x ( y z ) ] x<br />

x 2 ⎢ ⋅ − = σ − ν σ + σ = ε<br />

∂σ<br />

μ 2(<br />

1 )<br />

⎥<br />

⎣ + ν ∂σx<br />

∂σx<br />

⎦ E<br />

ili<br />

∂e<br />

∂τ<br />

p<br />

xy<br />

1 ⎡ ∂I<br />

=<br />

2<br />

⎢−<br />

μ<br />

⎣<br />

∂τ<br />

2<br />

xy<br />

⎤ τxy<br />

⎥ = = ε<br />

⎦<br />

2μ<br />

xy,


<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 19 -<br />

tj. odgovarajuće komponente <strong>deformacija</strong>. I obratno, deriviranjem istih izraza po komponentama<br />

<strong>deformacija</strong><br />

∂e<br />

∂ε<br />

p<br />

x<br />

⎡ 1−<br />

ν ∂e1<br />

∂e<br />

= 2μ<br />

⎢ e1<br />

⋅ −<br />

⎣1−<br />

2ν<br />

∂εx<br />

∂ε<br />

∂e<br />

∂ε<br />

p<br />

xy<br />

⎡ ∂e<br />

= 2μ<br />

⎢−<br />

⎣<br />

∂ε<br />

2<br />

xy<br />

2<br />

x<br />

⎤<br />

⎥ = λe1<br />

⎦<br />

⎤<br />

⎥ = 2μ<br />

ε<br />

⎦<br />

xy<br />

+ 2με<br />

dobiju se komponente <strong>naprezanja</strong>. Isto, dakako, vrijedi i za ostale komponente tenzora<br />

<strong>naprezanja</strong> odnosno <strong>deformacija</strong>. Tako je:<br />

∂e<br />

∂ε<br />

∂e<br />

∂ε<br />

p<br />

xx<br />

p<br />

xy<br />

= σ<br />

= τ<br />

4.7 Duhamel – Neumanove jednadžbe<br />

xx<br />

xy<br />

Koristeći se principom superpozicije, ukupna se <strong>deformacija</strong> εij može napisati u obliku:<br />

ε<br />

ij<br />

= ε<br />

σ<br />

ij<br />

∂e<br />

∂σ<br />

∂e<br />

∂τ<br />

+ ε<br />

T<br />

ij<br />

xx<br />

p<br />

p<br />

xy<br />

= ε<br />

= ε<br />

xy<br />

= τ<br />

xx<br />

.<br />

xy<br />

x<br />

= σ<br />

x<br />

(4.49)<br />

(4.50)<br />

pri čemu je σ<br />

εij - <strong>deformacija</strong> nastala zbog <strong>naprezanja</strong> i T εij - <strong>deformacija</strong> nastala zbog promjene<br />

temperature.<br />

Promjena temperature je:<br />

a <strong>deformacija</strong> tijela zbog temperaturne razlike je:<br />

T T T − = Δ (4.51)<br />

T<br />

ij<br />

0<br />

ε = α ΔT<br />

δ<br />

(4.52)<br />

pri čemu je T temperatura tijela, T0 referentna temperatura, α koeficijent toplinskog rastezanja.<br />

Naprezanje tijela zbog temperaturne razlike je jednako nuli ako je deformiranje tijela<br />

slobodno, tj.:<br />

T =<br />

ij<br />

σ 0<br />

(a)<br />

Ako se u generalizirani Hookeov zakon (4.33) uključi izraz (4.52), dobiju se Duhamel –<br />

Neumanove jednadžbe:<br />

1 ⎛ λ ⎞<br />

ε ij = ⎜σij<br />

− σkkδij<br />

⎟ + α ΔTδij<br />

(4.53)<br />

2μ<br />

⎝ 2μ<br />

+ 3λ<br />


ili u obratnom obliku:<br />

T<br />

ij =<br />

Za i ≠ j je σ 0 i ε 0 .<br />

ij<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 20 -<br />

ij<br />

kk<br />

ij<br />

( 2μ<br />

+ 3λ)<br />

αΔT<br />

ij<br />

σ = 2με + λε δ −<br />

δ<br />

(4.54)<br />

T<br />

ij =<br />

Budući da promjena temperature utječe samo na dijagonalne članove tenzora deformacije, u<br />

eksplicitnom se obliku Duhmael – Neumanove jednadžbe dobiju dodavanjem vrijednosti (4.52)<br />

na vrijednosti εxx, εyy i εzz izraza (4.33).<br />

4.8 Primjeri<br />

PRIMJER 4.1<br />

Stanje deformacije u točki zadano je tenzorom deformacije:<br />

[ ε ]<br />

ij<br />

=<br />

−<br />

10 4<br />

⎡ 1<br />

⎢<br />

⎢<br />

0.<br />

3<br />

⎢⎣<br />

0.<br />

8<br />

0.<br />

3<br />

2<br />

− 0.<br />

5<br />

Zadano je takoñer: E = 210 GPa, ν = 0.3. Potrebno je:<br />

(a) napisati tenzor <strong>naprezanja</strong><br />

(b) odrediti modul kompresije<br />

(a) tenzor <strong>naprezanja</strong><br />

0.<br />

8⎤<br />

− 0.<br />

5<br />

⎥<br />

⎥<br />

0.<br />

5⎥⎦<br />

Komponente tenzora <strong>naprezanja</strong> će se odrediti pomoću jednadžbi (4.32).<br />

Prije toga teba odrediti Lameove konstante μ i λ:<br />

λ =<br />

Volumenska <strong>deformacija</strong> je:<br />

E 210<br />

μ = G = = =<br />

2<br />

νE<br />

( 1+<br />

ν)<br />

2(<br />

1+<br />

0.<br />

3)<br />

=<br />

0.<br />

3⋅<br />

210<br />

( 1+<br />

ν)(<br />

1−<br />

2ν)<br />

( 1+<br />

0.<br />

3)(<br />

1−<br />

2⋅<br />

0.<br />

3)<br />

ε<br />

v<br />

= ε<br />

x<br />

+ ε<br />

Pomoću jednadžbi (4.32) dobije se:<br />

σ<br />

σ<br />

σ<br />

x<br />

y<br />

z<br />

= 2με<br />

= 2με<br />

= 2με<br />

x<br />

y<br />

z<br />

+ λε<br />

+ λε<br />

+ λε<br />

v<br />

v<br />

v<br />

y<br />

+ ε<br />

= 2 ⋅80.<br />

77 ⋅10<br />

⋅1⋅10<br />

= 2 ⋅80.<br />

77 ⋅10<br />

⋅ 2 ⋅10<br />

z<br />

=<br />

= 2 ⋅80.<br />

77 ⋅10<br />

⋅ 0.<br />

5 ⋅10<br />

3<br />

3<br />

3<br />

−4<br />

−4<br />

80.<br />

77<br />

GPa<br />

= 121.<br />

15 GPa<br />

−4<br />

−4<br />

( 1+<br />

2 + 0.<br />

5)<br />

⋅10<br />

= 3.<br />

5⋅10<br />

−4<br />

3<br />

3<br />

3<br />

−4<br />

+ 121.<br />

15⋅10<br />

⋅3.<br />

5 ⋅10<br />

−4<br />

+ 121.<br />

15⋅10<br />

⋅ 3.<br />

5⋅10<br />

−4<br />

+ 121.<br />

15 ⋅10<br />

⋅ 3.<br />

5⋅10<br />

=<br />

=<br />

58.<br />

5<br />

77.<br />

7<br />

=<br />

MPa<br />

50.<br />

48<br />

MPa<br />

MPa


τ<br />

τ<br />

τ<br />

xy<br />

yz<br />

zx<br />

= Gγ<br />

= Gγ<br />

= Gγ<br />

xy<br />

yz<br />

zx<br />

Tenzor <strong>naprezanja</strong> je:<br />

(b) modul kompresije<br />

3<br />

3<br />

3<br />

<strong>Veza</strong> izmeñu <strong>naprezanja</strong> i <strong>deformacija</strong> - 21 -<br />

−4<br />

= 80.<br />

77 ⋅10<br />

⋅ 2⋅<br />

0.<br />

3⋅10<br />

−4<br />

= 80.<br />

77 ⋅10<br />

⋅ 2⋅<br />

0.<br />

8⋅10<br />

−4<br />

= 80.<br />

77 ⋅10<br />

⋅ 2⋅<br />

( −0.<br />

5)<br />

⋅10<br />

=<br />

4.<br />

86<br />

MPa<br />

= −8.<br />

46 MPa<br />

= 12.<br />

9 MPa<br />

[ σ ] =<br />

⎢<br />

4.<br />

86 77.<br />

7 −8.<br />

16<br />

⎥<br />

MPa<br />

ij<br />

⎡58.<br />

5<br />

⎢<br />

⎢⎣<br />

12.<br />

9<br />

4.<br />

86<br />

−8.<br />

16<br />

Modul kompresije odreñen je izrazom iz tablice 4.1:<br />

PRIMJER 4.2<br />

12.<br />

9<br />

58.<br />

48<br />

E 210<br />

K = =<br />

= 175 GPa<br />

3<br />

( 1−<br />

2ν)<br />

3(<br />

1−<br />

2⋅<br />

0.<br />

3)<br />

Stanje <strong>naprezanja</strong> u točki zadano je tenzorom <strong>naprezanja</strong>:<br />

Tijelo je zagrijano za ΔT = 20 o C.<br />

80<br />

⎤<br />

⎥<br />

⎥⎦<br />

[ σ ] =<br />

⎢<br />

− 60 30 20<br />

⎥<br />

MPa<br />

ij<br />

⎡<br />

⎢<br />

⎢⎣<br />

− 40<br />

− 60<br />

Zadano je: E = 210 GPa, ν = 0.3; α = 1.2 . 10 -5 1/ o C<br />

Napisati tenzor deformacije.<br />

20<br />

− 40⎤<br />

⎥<br />

50⎥⎦

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

Saved successfully!

Ooh no, something went wrong!