5. Fasta ¨ amnens elastiska egenskaper (h˚ ardhet o. dyl)
[AM 22, Kittel s.80- ]
Med elasticitet menas enkelt sagt det hur fasta ¨ amnen ger efter n¨ ar man pressar p˚ a dem. Den klassiska elasticitetsteorin ¨ ar i sj¨ alva verket en kontinuitets-teori men utg˚ aende fr˚ an atomernas v¨ axelverkningsmodeller ¨ ar det m¨ ojligt att h¨ arleda denna teori.
Den grundl¨ aggande terminologin ¨ ar p˚ a engelska, svenska och finska:
- stress = tryck eller sp¨ anning = j¨ annitys. Har tryck-enheter, kraft/area.
- strain = sp¨ anning eller utt¨ ojning eller deformation = j¨ annitys tai venym¨ a. ¨ Ar enhetsl¨ os.
Den svenska terminologin ¨ ar allts˚ a lite problematisk d˚ a ordet sp¨ anning kan ha tv˚ a olika betydelser, men engelskan ¨ ar klar. Med “stress” menas allts˚ a trycket som ett material uts¨ atts f¨ or, och “strain”
hur mycket den dras ut (positiva v¨ arden) eller pressas ihop (negativa v¨ arden). P˚ a svenska ¨ ar det bra att anv¨ anda ordet “utt¨ ojning” f¨ or det kan inte missf¨ orst˚ as och har intuitivt r¨ att f¨ ortecken.
Den klassiska elasticitetsteorin ger sambandet mellan tryck och utt¨ ojning, och h¨ arleds till n¨ asta
utg˚ aende fr˚ an atom¨ ar struktur i material.
5.1. H¨ arledning av elasticitetstensorn
Betrakta en godtycklig interatom¨ ar potential V . Vi s˚ ag i f¨ orra kapitlet att de alltid har en t¨ amligen j¨ amn potentialgrop kring atomernas j¨ amviktsl¨ agen R . Allts˚ a kan vi utveckla potentialen med hj¨ alp av den tredimensionella versionen av Taylors serie
f (r + a) = f (r) + a∇f (r) + 1
2! (a · ∇)
2f (r) + 1
3! (a · ∇)
3f (r) + . . . Vi betraktar atomer som har f¨ orflyttats fr˚ an positionerna
R till r . Nu kan man skriva:
r(R) = R + u(R) (1)
d¨ ar u ¨ ar en m¨ angd vektorer som beskriver f¨ orflyttningarna.
Om vi nu har en potentialfunktion V (r), och utvidgar den kring j¨ amviktsavst˚ andet R f¨ or en atom som har f¨ orflyttats lite fr˚ an potentialgropens mitt till positionen r , f˚ ar vi f¨ or
den totala energin U i systemet
RInteratomär potential V(r)
r
Interatomärt avstånd r
u(R)
U ≈ N 2
X V (R − R
0) + 1 2
X
RR0
(u(R) − u(R
0)) · ∇V (R − R
0) (2)
+ 1 4
X
RR0
[(u(R) − u(R
0)) · ∇]
2V (R − R
0) + . . . (3)
d¨ ar summorna l¨ oper ¨ over atomens och dess grannars j¨ amviktspositioner R och R ’.
Den f¨ orsta termen i utvecklingen P V (R) ¨ar ju bara gittrets j¨amviktsenergi eller kohesions-energi.
Inom elasticitetsteori har nollpunktens l¨ age ingen skillnad, och man v¨ aljer f¨ or enkelhets skull att flytta nollpunkten s˚ a att U = 0 d˚ a alla atomer ¨ ar i sina j¨ amviktsl¨ agen.
Termen ∇V (R − R
0) f¨ orsvinner uppenbart, ty den ¨ ar ju helt enkelt kraften p˚ a atomer i deras j¨ amviktsl¨ agen, som m˚ aste vara ≡ 0 i en kristall i j¨amvikt.
Den andra derivatan ¨ ar den f¨ orsta termen i utvecklingen som blir kvar. I den harmoniska ap-
proximationen f¨ or fasta ¨ amnen arbetar man bara med denna term. De h¨ ogre termerna kallas
anharmoniska termer, och har betydelse t.ex. f¨ or v¨ arme-expansion.
I den klassiska f¨ orsta ordningens elasticitetsteori anv¨ ander man sig av den harmoniska approxi- mationen. Omr˚ adet d¨ ar denna approximation fungerar bra kallas linj¨ ar elasticitet.
Den harmoniska termen
[(u(R) − u(R
0)) · ∇]
2V (R − R
0) har inalles 9 termer, som ocks˚ a kan skrivas
X
µ,ν
[u
µ(R) − u
µ(R
0)]D
µν(R − R
0)[u
ν(R) − u
ν(R
0)]
d¨ ar µ och ν l¨ oper ¨ over koordinaterna x, y, z. Vi har allts˚ a ersatt vektornotationen med en index-notation. D inneh˚ aller nu derivatorna i formen:
D
µν= ∂
2V (r)
∂r
µ∂r
νoch vi kan skriva hela den harmoniska potentialenergin
U = − 1 4
X
R,R0 µ,ν
[u
µ(R
0) − u
µ(R)]D
µν(R
0− R)[u
ν(R
0) − u
ν(R)]
(Bytet av ordning ¨ ar m¨ ojligt p.g.a. gittrets inversions-symmetri och den harmoniska formen: ett teckenbyte i f¨ orflyttningen leder till samma energi.)
Hittills har f¨ orflyttningarna i atomernas platser u (R ) kunnat vara helt godtyckliga:
Men detta ¨ ar inte s¨ ardeles sannolikt f¨ or en elastisk f¨ orflyttning av atomer i j¨ amvikt. Mycket troligare
¨ ar att f¨ orflyttningsf¨ altet u(R) varierar j¨ amnt och kontinuerligt. Detta ¨ ar intuitivt enkelt och f¨ orst˚ a
om man t¨ anker p˚ a hur f¨ orflyttningarna kan uppst˚ a fr˚ an ett sk. dilatationscenter (t.ex. en st¨ orre atom i en kristall); det ¨ ar t¨ amligen klart att dessa ¨ ar kontinuerliga kring den st¨ orre atomen:
F¨ or att komma vidare m˚ aste antar vi dessutom att v˚ art gitter ¨ ar ett monatom¨ art Bravais-gitter,
s˚ a att alla punkter och “bindningar” ¨ ar ekvivalenta ((det kan uppst˚ a intressanta komplikationer i
icke-Bravais-gitter, som dock ¨ ar klart ¨ overkurs f¨ or denna kurs)).
Om vi nu definierar en kontinuerlig funktion u (r ), som ¨ ar = u (R ) i Bravais-gittrets punkter R kan vi g¨ ora Taylor-approximationen
u(R
0) = u(R) + (R
0− R) · ∇u(r)
r=Roch d¨ armed
U = − 1
4
X
R,R0 µ,ν
[(R
0− R) · ∇u
µ(R)]D
µν(R
0− R)[(R
0− R) · ∇u
ν(R)]
= − 1
4
X
R,R0 µ,ν,σ,τ
(R
0σ− R
σ) ∂u
µ(R)
∂x
σD
µν(R
0− R)
(R
0τ− R
τ) ∂u
ν(R)
∂x
τd¨ ar vi allts˚ a bara i senare steget skrivit om nablorna i summeringsnotation.
Nu kan man ytterligare ta ut ∇-termerna som ju beror bara p˚ a R , och s¨ atta in termerna med R
0i
en inre summa. D˚ a dessutom nu (R
0− R) ju m˚ aste ocks˚ a vara en vektor i Bragg-gittret, kan man
byta summerings-indexet R
0till R
00och f˚ a:
U = − 1 2
X
µ,ν,σ,τR
∂u
µ(R)
∂x
σ∂u
ν(R)
∂x
τ1 2
X
R00
R
00σD
µν(R
00)R
00τ(4)
Den senare termen ¨ ar en tensor av fj¨ arde rang, betecknas med
E
σµτ ν= − 1 2
X
R
R
σD
µν(R)R
τ(5)
och kallas elasticitetstensorn.
Det viktiga ¨ ar nu att E
σµτ νenbart beror p˚ a vektorerna R , dvs. atomernas perfekta positioner, inte p˚ a utt¨ ojningsf¨ altet. Allts˚ a ¨ ar E
σµτ νen konstant (med 81 komponenter) f¨ or ett specifikt material.
Nu vidare f¨ or att u ju varierar l˚ angsamt, kan man skriva om ekv. 4 som en integral,
U = 1 2
X
µ,ν,σ,τ
Z
dr ∂u
µ(R)
∂x
σ∂u
ν(R)
∂x
τE ¯
σµτ ν(6)
d¨ ar
E ¯
σµτ ν= E
σµτ νΩ
d¨ ar Ω ¨ ar volymen per en primitiv enhetscell i Bravais-gittret (och allts˚ a ocks˚ a volymen per atom).
Ekvation 6 ¨ ar utg˚ angspunkten f¨ or hela den klassiska elasticitetsteorin. Den ger energin f¨ or ett specifierat utt¨ ojningsf¨ alt u om man k¨ anner materialets elasticitetstensor E.
Trots att vi utgick fr˚ an atomer i v˚ ar h¨ arledning, ser man att ekvationen ¨ ar en ren kontinuitetsekvation.
F¨ or att sammanfatta detta ur fysikalisk synvinkel: med en h¨ arledningen som b¨ orjar fr˚ an interatom¨ ara potentialer, kommer man med ett par fysikaliskt mycket rimliga approxi- mationer fram till en helt kontinuerlig teori! D¨ armed ¨ ar klassisk (linj¨ ar) elasticitetsteori v¨ almotiverad, och anv¨ ands vitt uttryckligen som en kontinuum-teori! Moderna “Fini- te Element Method”-simuleringar av materials mekanisk egenskaper, som ¨ ar mycket anv¨ anda i industrin, baserar sig p˚ a linj¨ ar elasticitetsteori.
Tensorn E har hela 81 element(!), s˚ a detta verkar inte alltf¨ or trevligt att arbeta med. Det visar sig
dock att man snabbt kan drastiskt minska p˚ a antalet termer som m˚ aste beaktas genom att betrakta
olika symmetrier.
Den f¨ orsta symmetrin ¨ ar enkel. Det ¨ ar uppenbart att ekvation 5 ¨ ar symmetrisk med avseende p˚ a utbyte av indexena σ ↔ τ . Likas˚ a ¨ ar D
µνsymmetrisk med avseende p˚ a utbyte av µ ↔ ν p.g.a.
den harmoniska approximationen. Allts˚ a r¨ acker det med att best¨ amma E
σµτ νf¨ or de sex v¨ ardena
xx, yy, zz, yz, zx, xy
f¨ or indexparet µν och στ var f¨ or sig. Detta har redan s¨ ankt antalet element som b¨ or betraktas till 6 × 6 = 36.
N¨ asta steg f¨ or att reducera index f¨ oljer fr˚ an att betrakta en rotation av kristallen. Det ¨ ar uppenbart att en rotation kan inte ¨ andra p˚ a kristallens energi. S˚ a om vi nu g¨ or rotationen
u(R) = ∂ω × R,
s¨ atter in denna modifikation av gittret i ekv. 4 och kr¨ aver att U = 0, kan man l¨ att visa att U kan bero bara p˚ a en symmetrisk kombination av u-derivatorna,
ε
σµ= 1 2
∂u
µ∂x
σ+ ∂u
σ∂x
µ(7)
Nu kan man skriva om ekv. 4 i formen
U = 1 2
Z dr
X
µ,ν,σ,τ
ε
σµc
σµτ νε
τ ν
(8)
d¨ ar
c
σµτ ν= − 1 8Ω
X
R
[R
σD
µνR
τ+ R
µD
σνR
τ+ R
σD
µτR
ν+ R
µD
στR
ν] (9)
Betraktelse av dessa ekvationer visar att uppenbart ¨ ar ekvationerna symmetriska vid utbyte av σµ ↔ τ ν. Fr˚ an den senare ekvationen ser man ocks˚ a genast att c
σµτ ν¨ ar invariant vid transformationerna σ ↔ µ och τ ↔ ν. Detta reducerar antalet oberoende index i c
σµτ νtill 21 (antalet kommer av de 6 diagonal-elementen som bildas av matrisen som bildas av de 6 indexparena σµ plus de 6 indexparena τ ν = σµ + de 15 elementen i ena halvan av matrisen p˚ a n˚ agondera sidan diagonalelementen).
I deth¨ ar skedet ¨ ar det dags att introducerar en ny notation, med l˚ angt f¨ arre index ¨ an i notationen
ovan.
5.1.1. Interlud: samband med termiska f¨ orskjutningar
F¨ ore vi forts¨ atter med elasticitet, p˚ apekar vi en viktig fysikalisk till¨ aggsaspekt om f¨ orskjutningarna. Vi ser senare p˚ a kursen att atomernas termisk f¨ orskjutningar kring j¨ amviktsl¨ ahet ¨ ar av storleksordningen 0.1 ˚ A vid normala temperaturer.
F¨ orskjutningarna fr˚ an elasticitetsf¨ alt kan vara mycket mindre, ∼ 0.001 ˚ A. Betyder detta att de ¨ ar betydelsl¨ osa vid normala temperaturer?
Nej! Sambandet ¨ ar det att de elastiskaf¨ orskjutningarna f¨ orflyttar p˚ a j¨ amviktspunkten f¨ or de termiska.
I andra ord, om man tar atomernas momentana platsers medeltal ¨ over en l˚ ang tid, f˚ ar man den
elastiskt f¨ orskjutna positionen.
F¨ orskjutningar i kristaller kan m¨ atas med t.ex. r¨ otngen- eller elektrondiffraktion. D˚ a de m¨ ater ¨ over
makroskopiska tider, medan atomers vibrationer sker p˚ a ps-tidskalor, ger m¨ atningarna uttryckli-
gen medeltalsf¨ orskjutningarna av atomer (de termiska f¨ orskjutningarna p˚ averkar breddningen av
diffraktionspikarna).
5.1.2. Ingenj¨ orsnotation f¨ or elasticitet
Notation som anv¨ andes ovan ¨ ar den vetenskapligt b¨ ast motiverade och elegantaste, men samti- digt ocks˚ a n˚ agot otymplig med alla sina index. Den ¨ overl¨ agset vanligare notationen ¨ ar den s.k.
ingenj¨ orsnotationen (“engineering notation”) f¨ or elasticitet. Den utg˚ ar fr˚ an det att man ens i v¨ arsta fall inte beh¨ over mer ¨ an 36 olika element f¨ or att beskriva elasticitet.
I notationen ¨ overf¨ or man koordinatparena µν och στ till index fr˚ an 1-6 enligt f¨ oljande:
xx → 1, yy → 2, zz → 3, , yz → 4, zx → 5, xy → 6 (10) Nu introduceras nya elasticitetselement C
αβenligt transformationerna
C
αβ= c
σµτ νd¨ ar α ↔ σµ, β ↔ τ ν enligt samma regler som i (10). P˚ a liknande s¨ att omvandlas komponenterna ε
µνtill komponenter e, men s˚ a att
e
µν=
ε
µνd˚ a µ = ν
2ε
µνd˚ a µ 6= ν
och nu kan vidare e ocks˚ a betecknas e
αenligt transformationsreglerna 10.
Nu kan allts˚ a ekvationen f¨ or energit¨ atheten i kristallen skrivas i den betydligt enklare formen U =
12X
αβ
Z
dre
αC
αβe
β(11)
Konstanterna C
αβkalla elastiska styvhets-konstanter eller elastiska moduler och bildar allts˚ a en 6 × 6-matris. Elementen i den inversa matrisen S = C
−1kallas “elastic compliance constants”
eller helt enkelt elastiska konstanter.
5.1.3. Geometrisk tolkning av e och C
Efter all denna h¨ arledningsrumba har man l¨ att tappat taget till vad i all v¨ arlden detta nu hade och g¨ ora med materials h˚ ardhet, och hurdan modifiering av materialet e:s komponenter nu egentligen betyder ?
F¨ or att f¨ orst˚ a detta b¨ orjar vi med att betrakta f¨ oljande enkla, likformiga modifikation av ett r¨ atvinkligt koordinatsystem med enhetsvektorerna i , j och k
x
0= (1 +
xx)i +
xyj +
xzk
y
0=
yxi + (1 +
yy)j +
yzk (12)
z
0=
zxi +
zyj + (1 +
zz)k
vars geometriska tolkning ¨ ar en utt¨ ojning av koordinataxlarna som beskrivs av
xx,
yyoch
zz,
samt en vridning eller skjuvning som beskrivs av de ¨ ovriga indexena:
Notera: detta 6= ε som anv¨ andes ovan! Men vi ser snart att de nog har ett samband.
Variablerna << 1 d˚ a vi r¨ or oss inom det linj¨ ara elastiska omr˚ adet.
Om man g¨ or denna deformation av koordinataxlarna, hur p˚ averkas d˚ a positionen av en atom som ursprungligen ligger vid r = (xi + yj + zk) ? Efter deformationen kommer dess plats att vara r
0= xx
0+ yy
0+ zz
0, och f¨ orflyttningen (“displacement”) u (samma som i ekvation (1))
u(r) = r
0− r = x(x
0− i) + y(y
0− j) + z(z
0− k)
som genom ins¨ attning i 12 blir
u(r) = (x
xx+ y
yx+ z
zx)i + (x
xy+ y
yy+ z
zy)j + (x
xz+ y
yz+ z
zz)k som kan ocks˚ a skrivas i en mer allm¨ an form
u(r) = u
x(r)i + u
y(r)j + u
z(r)k
H¨ ar ser vi allts˚ a att deformationf¨ altet u(r) kan skrivas med hj¨ alp av enkla transformationer av koordinataxlarna som beskrivs av koefficienterna . Vidare ser man genast att
xx= ∂u
x∂x ,
xy= ∂u
y∂x , osv.
Dessa ekvationer g¨ aller f¨ or ¨ ovrigt som f¨ orsta approximation ocks˚ a f¨ or icke-likformiga distortioner av rymden.
Nu kan man igen definiera koefficienter e
αβp˚ a basen av koefficienterna . Koefficienterna med samma index ¨ ar helt enkelt
e
xx=
xx= ∂u
x∂x ; e
yy=
yy= ∂u
y∂y ; e
zz=
zz= ∂u
z∂z (13)
och koefficienter med olika index definieras s˚ a att de beror p˚ a ¨ andringen av vinkeln mellan koordinataxlarna
e
xy= x
0· y
0u
yx+
xy= ∂u
x∂y + ∂u
y∂x e
yz= y
0· z
0u
zy+
yz= ∂u
y∂z + ∂u
z∂y (14)
e
zx= z
0· x
0u
zx+
xz= ∂u
x∂z + ∂u
z∂x
J¨ amf¨ orelse med ekvationerna f¨ or e ovan visar att dessa definitioner faktiskt ¨ ar helt samma.
I.o.m. att ¨ andringarna i koordinataxlarna ¨ ar sm˚ a, g¨ aller e
xy= x
0· y
0= |x
0||y
0| cos θ ≈ √
1 + 2
xxp1 + 2
yycos θ
≈ (1 +
xx)(1 +
yy) cos θ ≈ (1 +
xx+
yy) cos θ d¨ ar alla termer av ordningen
2l¨ amnats bort omedelbart.
Nu ser vi vad sambandet med koefficienterna ε, och e ¨ ar. F¨ or koefficienterna med samma index
g¨ aller helt enkelt
αα= ε
αα= e
αα= e
nmedan f¨ or koefficienterna med olika index g¨ aller
αβ+
βα= 2ε
αβ= e
αβ= e
nd¨ ar indexet n best¨ ams av villkoren (10).
F¨ or att notationens-konfusionen skulle vara fullst¨ andig, anv¨ ander en del k¨ allor ε n¨ ar de menar e som ¨ ar definierat h¨ ar, osv. Om man arbetar med dessa borde man d¨ arf¨ or alltid definiera sin notation i b¨ orjan p˚ a n˚ agot entydigt s¨ att !
Alla dessa index kallas “strain” p˚ a engelska, p˚ a svenska kan man kalla dem sp¨ annings-komponenter eller utt¨ ojnings-komponenter.
Men h¨ ar kommer vi i forts¨ attningen att anv¨ anda fr¨ amst ingenj¨ orsnotationen med C
ijoch e
ij.
Nu kan vi ocks˚ a l¨ att ber¨ akna hur en transformation av gittret ¨ andrar p˚ a volymen. Efter deformationen blir volymen
V
0= V (x
0· y
0× z
0) u V (1 + e
xx+ e
yy+ e
zz)) (15)
d¨ ar vi igen l¨ amna bort termer O(e
2), och volymf¨ or¨ andringen, dilatationen, δ = V
0− V
V u e
xx+ e
yy+ e
zzSambandet med tryck ¨ ar nu l¨ att att h¨ arleda. Enligt termodynamikens andra grundlag ¨ ar dE = T dS − P dV
Vid T = 0 g¨ aller allts˚ a dE = −P dV . Detta ¨ ar allts˚ a ekvationen f¨ or gaser, och P ¨ ar en skal¨ ar som kallas hydrostatiskt tryck. F¨ or gitter arbetar man dock med ett generaliserat tryck i tensorform som ofta betecknas σ
ijoch kallas “stress” p˚ a engelska. Det kunde kallas p˚ afrestning eller belastning p˚ a svenska, men vi anv¨ ander helt enkelt tryck. Med “stress” menas det yttre arbete som g¨ ors p˚ a ett system, och dess komponenter har motsatt tecken till det inre trycket i systemet. Notera att detta leder ocks˚ a till att det ¨ ar fullt m¨ ojligt att ha negativa tryck inom materialfysiken, i motsats till klassisk termodynamik f¨ or gaser d¨ ar tryck ¨ ar alltid positiva.
Energiekvationens generalisering blir nu
dU/V = X
ij
σ
ijde
ijd¨ ar U/V ¨ ar energit¨ atheten i systemet. J¨ amf¨ orelse med ekv. (8) visar nu att
σ
ij= X
kl
c
ijklkl(16)
Denna lag som vi h¨ arlett s¨ ager allts˚ a att trycket har ett linj¨ art beroende p˚ a utt¨ ojningen, och kallas ocks˚ a Hookes lag.
(Lagen var ursprungligen var inte precis given i denna tensorform: Hooke sj¨ alv levde p˚ a 1600-talet och framf¨ orde lagen i formen “ut tensio sic vis”).
Komponenterna i tryck-tensorn σ har en klar geometrisk tolkning. De ¨ ar helt enkelt en kraft
per enhetsarea som p˚ averkar ett r¨ atvinkligt block av ett material som ¨ ar arrangerat l¨ angs med
koordinataxlarna
T.ex. σ
11= X
x¨ ar en kraft/area i X-riktningen som verkar p˚ a x-sidan, σ
12= X
yen kraft/area i X-riktningen som verkar p˚ a y-sidan osv. Den stora bokstaven betecknar alltid kraftens riktning, den lilla sidan den verkar p˚ a. Storheterna har enheter av kraft/area eller energi/volym, allts˚ a tryck-enheter.
Det existerar allts˚ a nio av dessa kraft-komponenter. Men f¨ or att kroppen skall h˚ allas i statiskt
j¨ amvikt, kan skjuvningskomponenterna inte vara oberoende. Betrakta t.ex. f¨ oljande fall:
Nu m˚ aste Y
x= X
yty annars hamnar blocket i en rotation. P˚ a motsvarande s¨ att g¨ aller Y
z= Z
yoch Z
x= X
z. De kvarvarande sex komponenterna kan v¨ aljas vara t.ex. X
x, Y
y, Z
z, Y
z, Z
xoch
X
y.
Nu kan man skriva om ekvationen (16) i ingenj¨ orsnotation s˚ a att
σ
11σ
22σ
33σ
23σ
31σ
12
=
X
xY
yZ
zY
zZ
xX
y
=
C
11C
12C
13C
14C
15C
16C
12C
22C
23C
24C
25C
26C
13C
23C
33C
34C
35C
36C
14C
24C
34C
44C
45C
46C
15C
25C
35C
45C
55C
56C
16C
26C
36C
46C
56C
66
e
xxe
yye
zze
yze
zxe
xy
(17)
H¨ ar har vi anv¨ ant oss av det faktumet att de sj¨ alvst¨ andiga elastiska konstanterna C ¨ ar h¨ ogst 21 till
antalet och skrivit t.ex. C
14i st¨ aller f¨ or C
41.
5.1.4. Ytterligare reducering av antalet elastiska moduler
Hittills har v˚ ar h¨ arledning g¨ allt f¨ or alla Bravais-gitter. Om man nu dessutom betraktar symmetrier hos ett specifikt gitter, kan man reducera antalet ytterligare.
T.ex. f¨ or kubiska system p˚ ast˚ ar vi att energins (11) t¨ athet U
tkan skrivas enbart med
U
t=
12C
11(e
2xx+ e
2yy+ e
2zz) +
12C
44(e
2yz+ e
2zx+ e
2xy) + C
12(e
yye
zz+ e
zze
xx+ e
xxe
yy) (18)
dvs. att termer som inneh˚ aller
(e
xxe
xy+ · · ·); (e
yze
zx+ · · ·); (e
xxe
yz+ · · ·);
inte f¨ orekommer. Detta kan visas genom att betrakta rotationer av gittret runt de fyra < 111 >-
axlarna s˚ a att koordinat-axlarna byts:
Lite omtanke visar att detta ¨ ar ekvivalent med f¨ oljande permuteringar av koordinat-axlarna:
x → y → z → x
−x → z → −y → −x x → z → −y → x
−x → y → z → −x
Dessa permuteringar ¨ andrar ocks˚ a tecken p˚ a utt¨ ojningsskomponenterna, t.ex. e
x,−y= −e
xy.
Uppenbart ¨ ar nu energit¨ atheten (18) invariant vid transformationerna. D¨ aremot kan tecknet p˚ a
termer av typen e
xxe
xy¨ andras vid transformationen, vilket kan komma att ¨ andra p˚ a energin.
D¨ armed kan dessa termer inte f¨ orekomma, och allts˚ a m˚ aste d¨ artillh¨ orande elastiska moduler vara
≡ 0.
Genom ins¨ attning i energit¨ athetens ekvation (11) kan man f¨ orvissa sig om att de numeriska faktorerna ocks˚ a ¨ ar r¨ att.
Allts˚ a beh¨ over vi bara tre (!) elastiska moduler C
11, C
12och C
44f¨ or att fullst¨ andigt beskriva de elastiska egenskaperna f¨ or kubiska system !
Detta ¨ ar en ganska dramatiskt framg˚ ang om man t¨ anker p˚ a att vi b¨ orjade med 81 m¨ ojliga index.
Genom att betrakta liknande transformationer kan man reducera antalet elastiska moduler f¨ or andra
Bravais-gitter. H¨ ar f¨ oljer en tabell ¨ over antalet moduler som beh¨ ovs:
De elastiska konstanterna kan i princip h¨ arledas om man k¨ anner till v¨ axelverkningarna mellan
atomer. Men de kan ocks˚ a m¨ atas direkt genom ultraljudsm¨ atningar (de har ett n¨ ara samband med
gittervibrationer, som behandlas senare p˚ a kursen).
5.2. Elastiska egenskaper f¨ or kubiska gitter
H¨ adanefter betraktar vi enbart kubiska system. F¨ or dessa kan allts˚ a samband mellan yttre tryck och utt¨ ojning skrivas p˚ a f¨ oljande s¨ att:
X
xY
yZ
zY
zZ
xX
y
=
C
11C
12C
120 0 0
C
12C
11C
120 0 0
C
12C
12C
110 0 0
0 0 0 C
440 0
0 0 0 0 C
440
0 0 0 0 0 C
44
e
xxe
yye
zze
yze
zxe
xy
(19)
Vi h¨ arleder nu n˚ agra makroskopiska h˚ ardhets-egenskaper utg˚ aende fr˚ an dessa mikroskopiska defini- tioner.
Det enklaste m˚ attet p˚ a ett materials elastiska h˚ ardhet ¨ ar bulkmodulen B, som helt enkelt ber¨ attar
hur mycket volymen p˚ a ett material ¨ andras om man trycker ihop det med ett hydrostatiskt yttre
tryck P . Den definieras som
B = −V dP dV
Inversen p˚ a bulkmodulen K = 1/B kallas kompressibilitet.
Ett konstant positivt yttre tryck motsvarar fallet X
x= Y
y= Z
z= −dP , som leder till en storleksf¨ or¨ andring e
xx= e
yy= e
zz. Alla skjuvningskomponenter X
y, e
xyetc. ¨ ar nu noll. Allts˚ a f˚ as
dP = −X
x= −(C
11e
xx+ C
12e
xx+ C
12e
xx) = −(C
11+ 2C
12)e
xxVolymf¨ or¨ andringen vid en transformation ¨ ar nu (se ekv. 15)
dV = V (e
xx+ e
yy+ e
zz) = 3e
xxV Allts˚ a f˚ as
B = −V −(C
11+ 2C
12)e
xx3e
xxV = 1
3 (C
11+ 2C
12)
Notera att B inte alls beror p˚ a C
44. Detta har sin f¨ orklaring i att C
44¨ ar ett m˚ att p˚ a skjuvning (vridning) av en kristall, och en j¨ amn ihoppressning leder ju inte till n˚ agon skjuvning..
Allts˚ a har vi p˚ a detta s¨ att visat att utg˚ aende fr˚ an rent mikroskopiska begrepp kan man h¨ arleda fasta
¨ amnens makroskopiska transformationsegenskaper.
H¨ ar ¨ ar en tabell p˚ a n˚ agra experimentella elastiska konstanter f¨ or grund¨ amnen:
F¨ or att g¨ ora en uppskattning p˚ a storleksordningen av den elastiska kompressionen, l˚ at oss ber¨ akna
hur mycket perfekt koppar skulle pressas ihop av ett tryck p˚ a 1 kbar. Med v¨ ardena i tabellen f˚ as B = 1.42×10
11N/m
2= 1420 kbar. Nu f˚ as dV /V = −dP /B = 0.0007 = 0.07%. Kompressibiliteten
¨ ar allts˚ a ganska s˚ a liten, vilket ju inte ¨ ar ¨ overraskande d˚ a den vardagliga erfarenheten ju klart s¨ ager att metaller inte g˚ ar att pressa ihop mycket.
Tv˚ a andra viktiga m˚ att f¨ or materials h˚ ardhet ¨ ar Youngs modul Y och Poissons kvot µ. Young’s modul kallas ocks˚ a ibland vilseledande den elastiska modulen E eller skjuvningsmodulen (“shear modulus”).
Dessa definieras p˚ a f¨ oljande s¨ att. Betrakta ett fritt block av material, som dras ut med ett tryck
X
x, medan alla andra tryck-komponenter ¨ ar = 0.
Nu kommer dess form att ¨ andras s˚ a att e
xx> 0 och e
yy= e
zz< 0. Youngs modul definieras som Y = X
xe
xxoch Poissons kvot som
µ = δw/w
δl/l = −e
yy/e
xxGenom ins¨ attning i ekv. 19 kan man visa att f¨ or kubiska kristaller g¨ aller (jfr. r¨ akne¨ ovningen) f¨ or utt¨ ojning l¨ angs med en h100i-riktning:
Y = (C
11+ 2C
12) C
11− C
12C
11+ C
12och
µ = C
12C
11+ C
12Notera att dessa ¨ ar de mikroskopiska definitionerna p˚ a Young’s modul. De ¨ ar inte n¨ odv¨ andigtvis
alls samma v¨ arden som kan erh˚ allas med en makroskopisk m¨ atning, och dessutom beror de p˚ a
kristallriktning (dessa ekvationer antar implicit att man betraktar en kubisk kristall och drar den ut
i 100-riktning).
F¨ or ber¨ akningar om Youngs modul i andra kristallriktingar, se
[Heino et al, Europhys. Lett. 41 (1998) 273].
Fr˚ an tabellen ovan kan man r¨ akna ut Young’s modul f¨ or n˚ agra vanliga metaller. H¨ ar ges ocks˚ a v¨ arden f¨ or j¨ amf¨ orelses skull f¨ or diamant, kolnanor¨ or och grafen (de tv˚ a sistn¨ amnda ¨ ar ¨ annu inte k¨ anda med h¨ og noggranhet, olika m¨ atningar ger avvikande v¨ arden, som dock ¨ ar av samma storleksordning som de i tabellen).
Amne ¨ Y (GPa)
W 408
Cu 72.5
Ni 150.8
Au 46.5
K 1.1
Diamant 1040
Grafen upp till ∼ 2400 har rapporterats
[Lee et al, Nano Lett. 12 (2012) 4444]Kolnanor¨ or upp till ∼ 1800 har rapporterats
[Treacy et al, Nature 381 (1996) 678]Diamant ¨ ar faktiskt ur denna synvinkels¨ att det h˚ ardaste materialet som finns – men samtidigt ¨ ar det mycket sk¨ ort, vilket begr¨ ansar dess nyttighet i praktiska till¨ ampningar!
D¨ aremot ¨ ar grafen och kolnanor¨ or inte sk¨ ora, vilket g¨ or dem mycket intressanta f¨ or mekaniska
till¨ ampningar – om man bara kunde ¨ overs¨ atta den h¨ oga mekaniska styrkan fr˚ an enskilda atomlager
till makroskopisk skala...
5.2.1. Isotropiska material
[Hirth and Lothe]
Fr˚ an beskrivningen ovan ¨ ar det klart att kristallens inre orientation kan p˚ averka dess elastiska egenskaper. Ett viktigt specialfall ¨ ar dock det d¨ ar s˚ a inte ¨ ar fallet, allts˚ a d˚ a de elastiska egenskaperna
¨ ar oberoende av riktning. S˚ adana material kallas isotropiska material.
F¨ or att se vad detta inneb¨ ar f¨ or de elastiska konstanterna, betrakta en expansion av x − axeln utan skjuvning eller modifikation av de ¨ ovriga axlarna. F¨ or enkelhets skull g¨ or vi detta f¨ or ett kubiskt system. Vi anv¨ ander nu indexena ε, som g˚ ar att transformera, men f¨ or vilka ekv. (19) g¨ aller p˚ a samma s¨ att som f¨ or e pga. den kubiska symmetrin. Allts˚ a ¨ ar ε
xx> 0 och alla andra ε = 0 . Ins¨ attning i ekv. 17 ger
σ
11= C
11ε
xx(20)
σ
22= σ
33= C
12ε
xx= C
12C σ
11(21)
Betrakta nu en rotation av 45
◦runt z-axeln, som har en transformations-matris
T
ij=
√1 2
√1
2
0
−
√12
√1
2
0
0 0 1
Man kan visa att en rotation omvandlar ε och σ enligt formeln ε
0ij= X
lm
T
ilT
jmε
lmoch motsvarande. f¨ or σ. D¨ armed f˚ as t.ex.
ε
0xy= T
11T
21ε
xx= −
12ε
xx(i.o.m. att alla andra termer utom ε
xx¨ ar 0 blir bara detta kvar) och
σ
120= T
11T
21σ
11+ T
12T
22σ
22= −
12σ
11+
12σ
22Men nu om vi betraktar den sista raden i ekv. (19) f¨ or det roterade koordinatsystemet, ser vi att
σ
120= 2C
44ε
012och genom sammanslagning av de tv˚ a senaste ekvationerna med ekv. (20) f˚ as
2C
44= C
11− C
12(22)
Detta ¨ ar allts˚ a villkoret f¨ or att ett kubiskt material ¨ ar isotropiskt. I.o.m. att ett isotropiskt material inte beror p˚ a gitterriktningar, kan man ocks˚ a f¨ or ¨ ovriga symmetrier visa att de oberoende gitterkonstanternas antal faktiskt reduceras till bara 2.
Allts˚ a har vi nu f¨ or isotropiska material reducerat antalet beh¨ ovliga elastiska konstanter ytterligare, till 2. L¨ angre ner ¨ an det kan man inte komma f¨ or fasta ¨ amnen.
F¨ or icke-isotropiska material beskrivs graden av anisotropi (“anisotropy ratio”) A som
A = 2C
44C
11− C
12som ¨ ar allts˚ a ≡ 1 f¨ or isotropiska material. H¨ ar ¨ ar n˚ agra v¨ arden p˚ a A
[Hirth-Lothe s. 837]:
Al 1.21
Cr 0.67
Cu 3.21
Au 2.9
W 1.00
C (diamant) 1.21
Si 1.56
Ge 6.35
F˚ a grund¨ amnen ¨ ar allts˚ a isotropiska: n˚ agra av de f˚ a exemplen ¨ ar W och (n¨ astan) Al. Men isotropiska
¨ amnen ¨ ar ¨ and˚ a extremt viktiga inom ingenj¨ orsvetenskapen i.o.m. att m˚ angkristallina ¨ amnen, vilket ju metaller oftast ¨ ar i praktiken, ¨ ar n¨ astan perfekt isotropiska p.g.a. den slumpm¨ assiga orientationen av kristallkornen. D¨ arf¨ or ¨ ar beskrivning av isotropiska material mycket viktig i praktiskt byggnaddsarbete och val av material. F¨ or deras beskrivning anv¨ ands ocks˚ a oftast i st¨ allet f¨ or de elastiska modulerna C helt enkelt Youngs modul Y och Poissons f¨ orh˚ allande µ, som allts˚ a nu r¨ acker f¨ or fullst¨ andig beskrivning av materialets linj¨ ara elastiska egenskaper.
Youngs modul Y och Poissons f¨ orh˚ allande µ ¨ ar ocks˚ a l¨ atta att m¨ ata genom att makroskopiskt dra
ut eller pressa ihop ett material och m¨ ata dess l¨ angd- resp. bredd-f¨ or¨ andring.
5.3. V¨ atskor
Nu har vi allts˚ a kommit ner fr˚ an 21 oberoende elastiska moduler till bara 2, genom att kontinuerligt f¨ orenkla strukturen i materialet. En logisk n¨ asta fr˚ aga ¨ ar att vilken f¨ orenkling skall man g¨ ora f¨ or att komma ner fr˚ an 2 till 1?
Svaret ¨ ar att man skall g¨ ora materialet till en v¨ atska! Definition p˚ a en v¨ atska ¨ ar i sj¨ alva verket att det ¨ ar ett kondenserat ¨ amne d¨ ar man inte kan m¨ ata en skjuvmodul, allts˚ a t.ex. Y . D˚ a ¨ ar den enda kvarvarande elastiska modulen hur mycket v¨ atskans volym kan kompresseras vid yttre tryck, dvs.
bulkmodulen B!
Den exakta definitionen f¨ or en v¨ atska inneh˚ aller i sj¨ alva verket ett numeriskt gr¨ ansv¨ arde f¨ or en
skjuvnings-effekt, som s¨ ager att ett fast ¨ amne ¨ ar ett material vars skjuvviskositet ¨ ar st¨ orre ¨ an 10
14.6poise.
5.4. Plasticitet och h˚ ardhet
5.4.1. Plasticitet
Hittills har vi hela tiden antagit att den harmoniska approximationen g¨ aller, och sett att det leder till ett enkelt linj¨ art f¨ orh˚ allande mellan en yttre kraft p˚ a ett material och dess deformation.
Fr˚ an vardaglig erfarenhet vet vi dock v¨ al att detta f¨ orh˚ allande inte kan h˚ alla f¨ or godtyckligt stora deformationer: f¨ or eller senare g˚ ar materialet s¨ onder.
Men redan f¨ ore ett material g˚ ar s¨ onder kan man skilja p˚ a olika beteende-omr˚ aden. Det elastiska omr˚ adet (av vilket den linj¨ ara delen som beskrevs ovan ¨ ar den f¨ orsta) k¨ annetecknas av att deformationer inte ¨ ar best˚ aende, utan att om man upph¨ or att p˚ averka p˚ a materialet med en kraft,
˚ aterst¨ aller den sig i sitt ursprungliga l¨ age.
Med plasticitet menas s˚ adana transformationer av material som leder till en best˚ aende ¨ andring i dess form efter att man slutat deformera det. H¨ ar ¨ ar det allts˚ a fr˚ aga om en sorts hysteresis-effekt.
Att man deformerar materialet s˚ a starkt att det ursprungliga l¨ aget inte ˚ aterst¨ alls betyder dock inte n¨ odv¨ andigtvis ¨ annu att det g˚ ar s¨ onder: materialets egenskaper p˚ a atomniv˚ a kan ¨ andras f¨ ore den g˚ ar slutgiltigt s¨ onder.
Skillnaden mellan elastiskt och plastiskt beteende syns tydligt i s.k. stress-strain-kurvor, t.ex.
Tv˚ a storheter som anv¨ ands f¨ or att beteckna denna praktiska h˚ ardhet ¨ ar flytgr¨ ansen “yield strength”.
Den best¨ ams genom som punkten d¨ ar trycket har l¨ att till en f¨ orflyttning av 0.002 i utt¨ ojningen fr˚ an
slutet p˚ a det elastiska omr˚ adet (se bilden ovan). En annan viktig storhet ¨ ar den dragh˚ allfastheten
(“tensile strength”), som ¨ ar helt enkelt punkten d¨ ar utt¨ ojningen b¨ orjar minska i ett utdragningstest:
H¨ ar ¨ ar en tabell ¨ over n˚ agra vanliga metallers flytgr¨ anser och tensila styrkor:
Det m˚ aste dock betonas att dessa v¨ arden varierar starkt inom samma material beroende p˚ a hur den
har tillverkats. En massa faktorer: kornstorlek, dislokationst¨ athet, orenhetsgrad, tom. hur metallen
har processerats, kan inverka med mer ¨ an en storleksordning p˚ a dess makroskopiska styrka! Dessa
beroenden ¨ ar i allm¨ anhet inte v¨ al k¨ anda eller f¨ orst˚ adda.
Ett undantag ¨ ar dock beroendet av flytgr¨ ansen σ
yp˚ a kornstorlek. F¨ or den g¨ aller, f¨ or en metall som annars har liknande egenskaper, den empiriskt erh˚ allna Hall-Petch-lagen som s¨ ager att
σ
y= σ
0+ k
y√ d (23)
d¨ ar d ¨ ar kornstorleken i medeltal, och σ
0och k
ymaterialkonstanter som b¨ or best¨ ammas.
Fr˚ an ekvationen ser man att desto mindre kornstorlek, desto h˚ ardare blir materialet. Detta ¨ ar orsaken till att nanokristallina material ¨ ar av stort intresse: bara med att minska p˚ a kornstorleken, kan man ¨ oka p˚ a h˚ ardheten. Experiment och datorsimuleringar visar att detta kan lyckas ner till storlekar d ∼ 10 nm.
Experimentellt har man kunnat visa att koppar med en kornstorlek d = 23 nm kan ha en flytgr¨ ans p˚ a ung. 770 MPa, allts˚ a faktiskt en storleksordning b¨ attre ¨ an v¨ ardet i tabellen ovan.
[Youssef, Appl. Phys.Lett. 85 (2004) 929]