Veza naprezanja i deformacija
Veza naprezanja i deformacija
Veza naprezanja i deformacija
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⎥⎦