Cilt: 1 Sayı: 3 sh. 81-88 Ekim 1999

 

TEK KRİSTAL SİSTEMLERDE ESNEKLİK ÖZELLİKLERİN İNCELENMESİ

(AN INVESTIGATION ON ELASTIC PROPERTIES OF SINGLE CRYSTAL SYSTEMS)

 

Coşkun İŞÇİ, E. Mine İLK

Dokuz Eylül Üniversitesi, Fen-Edebiyat Fakültesi, Fizik Bölümü, Alsancak, İZMİR

 

ÖZET / ABSTRACT

    Esneklik teorisinin tek kristal sistemlere uygulanmasıyla bulunan esneklik sabitleri ve ultrases hızları arasındaki bağıntılar kullanılarak ses hızının, Young modülünün ve sıkışabilirlik katsayısının, kristal doğrultularına bağlılıkları verilmiştir. Tek kristallerin esneklik özellikleri anizotropi gösterir. Örnek olarak, kübik simetriye sahip Ni tek kristali, nümerik hesaplamalarda kullanılmıştır. İkinci mertebeden gerilim ve deformasyon tensörleri ve dördüncü mertebeden esneklik sabitleri tensörü ile ilgili temel bilgi verildikten sonra, farklı doğrultulardaki hızlarla Cij sabitleri arasındaki bağıntılar çıkarılmıştır. Nümerik hesaplamalarda Mathematica 2.2 paket programı kullanılarak, hız ve Young modülü yüzeyleri bulunmuş ve sonuçlar grafik olarak çizilmiştir.

    By applying elasticity theory to single crystal systems, the relations between elastic constants and ultrasonic velocities, and the dependence of velocities, Young modulus and compressibility coefficient on crystallographic directions have been given. Elastic properties of single crystals show anisotropic character. As an example, Ni single crystal, which has a cubic symmetry, has been chosen for numerical calculations. After giving the theoretical background of second order stress , strain tensors and fourth order elastic constant tensor, the relations between velocities at different directions and Cij constants have been derived. For numerical calculations Mathematica 2.2 packet program has been used to calculate velocity and Young modulus surfaces and the results have been plotted.

ANAHTAR KELİMELER / KEY WORDS

Tek kristaller, Esneklik, Ultrasonik ses hızı, Young modülü

Single crysytals , Elasticity, Ultrasonic velocities, Young modulus

   

1. GİRİŞ

    Katı cisimlerin esneklik özelliklerinin incelenmesi ile atomlar ve moleküller arasındaki bağlayıcı kuvvetler hakkında bilgi elde edilebilir. Özel koşullarda büyütülen tek-kristaller (mono-crystal), esneklik açısından anizotropik özellik gösterirler. Farklı kristal doğrultularında ultrases hızı ve esneklik sabitlerinin ölçülmesi ile, esneklik özellikler ve bağlar hakkında yorumlar yapılabilir.

    Daha önce yapılan deneyler sonucunda elde edilen veriler ( Truell vd...,1969 ) analiz edilerek, farklı kristal yapıları için nümerik hesaplar yapılır ve akustik hız, Young modülü, sıkışabilirlik gibi fiziksel parametrelerin farklı kristal doğrultuları için grafikleri oluşturulabilir.

2. ESNEKLİK TEORİSİ

    Üzerine kuvvetler uygulanan bir katı cisim şekil değişikliğine uğrar. Bu kuvvetler kaldırıldığında bu şekil değişikliği yok olursa bu tür değişikliklere, esnek (elastik) şekildeğişimi denir. Birim yüzeye etki eden kuvvete gerilim (stress) denir. Anizotropik ortamlarda bu ikinci dereceden bir gerilim tensörü olan Tij ile gösterilir. Şekil değişikliği veya deformasyon (strain) de yine ikinci mertebeden bir tensörü olan Eij ile gösterilir ve aşağıdaki gibi ifade edilir.

Eij = ui / xj ( i, j = 1,2,3)                                                                                                                                   (1)

Burada ui, yer değiştirme vektörü u nun bileşenleridir. xj, yer vektörü r'nin bileşenleridir. Hooke kanuna göre gerilim (zor), deformasyon (zorlanma) ile lineer olarak orantılıdır.

Tij = Cijkl Ekl                                                                                                                                                            (2)

veya

Eij = Sijkl Tkl                                                                                                                                                           (3)

ile ifade edilir. Burada Cijkl (elastic stiffness constants) ve Sijkl (elastic compliance constants) esneklik sabitleri olup 4. mertebeden simetrik tensörlerdir.

Cijkl = Cklij = Clkij                                                                                                                                                 (4)

Bu sabitler matris gösteriminde

Ti = Cij Ej ve Ei =Sij Tj                                                                                                                                           (5)

olarak yazılabilir (İşçi, 1977 ve Nye, 1957). Simetrik özelliklerden (Cij = Cji) dolayı, en düşük simetrili kristallerde 21 bağımsız, Cij esneklik sabiti vardır. Kübik kristallerde bu sayı 3'e, hexagonal simetrilerde 6'ya inmektedir. Cij ve Sij ler arasında

Sij =(-1)i+j D cij /D c                                                                                                                                               (6)

bağıntısı vardır (Truell vd...,1969). Burada D c, Cij matrisinin determinantı, D cij de, Cij elemanın minörüdür.

Akustik dalga denklemi,

r ( 2u i / t2) = ( Tij/ Xj) dir.                                                                                                                              (7)

Burada r cismin yoğunluğudur. Hooke kanununun ifadesi (7)'de yerine konulur ve bazı ara işlemler yapılırsa, düzlem dalga için çözüm;

ui = u0i exp [i(wt- k.r)]                                                                                                                                           (8)

olarak bulunur. Burada ui, yer değiştirme vektörünün bileşenleri; u0i ,vektörün maksimum değeri; w, açısal frekans; k yayılma vektörü ( k=(2p /l )n ) ; r, yer vektörü; n, doğrultu kosinüsleri (n1, n2, n3) dir. Dalga hızları için Christoffel denklemleri en uygun ifadelerdir (Akgöz vd..., 1976):

(Lik -r v2 d ik) uok = 0 (i, k = 1,2,3)                                                                                                                       (9)

Lik = Cijkl nj nl                                                                                                                                                      (10)

Burada Lik'lar, Christoffel sabitleridir. (9) denklemi açık olarak yazılırsa;

(L11-r v2) u01 + L12 u02 + L13 u03 = 0

L12 u01 + (L22-r v2)u02 + L23 u03 = 0                                                                                                                   (11)

L13 uo1 + L23 u02 + (L33-r v2)u03 = 0

denklemleri elde edilir. Bu denklemlerin çözümü için uok'ların katsayılarının determinantı sıfır olmalıdır. Bu ifade, r v2 ye göre kübik bir denklemdir. Bir yayılma doğrultusu için üç farklı yayılma hızı vardır. Hızlar yerine konulduğunda 3 tane birbirine dik yayılma vektörü bulunur. Bir boyuna, iki enine dalga mevcuttur. Bazı özel doğrultularda iki enine dalga hızı birbirine eşit olabilir.

u x n =0 boyuna, u . n =0 enine dalgaları verir.

2.1. Kübik simetriye sahip kristallerde (Ni, Cu gibi) Christoffel sabitleri ve hız ifadeleri

    Kübik kristallerde Christoffel sabitleri:

L11= n12 C11 + n22 C44 + n32 C44

L22= n12 C44 + n22 C11 + n32 C44

L33= n12 C44 + n22 C44 + n32 C11

L23= n2 n3 ( C12 + C44)                                                                                                                                        (12)

L31= n3 n1 ( C12 + C44)

L12= n1 n2 ( C12 + C44)

ve [100] doğrultusu için ( n1 = 1, n2 = 0, n3 = 0 ) Christoffell denklemleri:

(C11 - r v2) u01 = 0

(C44 - r v2) u02 = 0                                                                                                                                              (13)

(C44 - r v2) u03 = 0

olarak verilir. Buradan vB boyuna dalga ve vE enine dalga hızları olmak üzere,

r v2B = C11 ve r v2E = C44                                                                                                                                   (14)

bulunur. Kübik kristallerde bazı özel doğrultularda ultrases hızı ile esneklik sabitleri arasında Çizelge 1'de verilen bağıntılardan ve ultrases hızı ölçümlerinden yararlanarak esneklik sabitleri hesaplanır.

Çizelge. 1. Kübik kristallerde ultrases hızı ile esneklik sabitleri arasındaki bağıntılar.

Yayılma

Doğrultusu

Polarizasyon

Doğrultusu

Bağıntı

[ 100]

[ 100]

r v2 = C11

[ 100]

(001) düzleminde

r v2 = C44

[ 110]

[ 110]

r v2 =1/2(C11 + C12 + 2C44 )

[ 110]

[ 001]

r v2 = C44

[ 110]

[ 110]

r v2 =1/2 (C11 - C12 )

[ 111]

[ 111]

r v2 =1/3 (C11 + 2C12 + 4C44)

[ 111]

(111) düzleminde

r v2 =1/3 (C11 + C44 - C12)

 

2.2. İzotropik Sistemlerde İlgili Bağıntılar

    İzotropik sistemlerde ( tüm polikristaller ) iki bağımsız esneklik sabiti ve iki hız vardır. Aşağıdaki ifadelerde verilen vB boyuna dalga hızı, vE enine dalga hızıdır:

v B = (C11 / r ) 1/2 , vE = (C44 / r ) 1/2

C44 = ½ (C11 - C12)                                                                                                                        (15)

Bu sistemlerde iki Lame sabitinden l ve m söz edilir (Truell vd..., 1969).

C11 = l + 2 m , C44 = m , C12 = l                                                                                                     (16)

2.3. Kübik ve İzotropik Sistemlerde Young Modülü, Hacim Modülü ve Sıkışabilirlik Katsayısı

    Kübik kristallerde; Young Modülünün tersi ni doğrultusunda,

(1/Y) = S11 - 2 (S11 - S12 - 0.5 S44) ( n12 n22 + n22n32 + n32n12)                                                                           (17)

ile verilir ( Nye, 1957). Burada,

S11 = (C11 + C12) / ((C11 - C12) (C11 + 2C12))

S12 = -C12 / ((C11 - C12) (C11 + 2C12))                                                                                                               (18)

S44 = 1 / C44

Sıkışabilirlik (linear compresibility), b = S11 + 2S12 ve hacimsel sıkışabilirlik, 3b 'dır.

İzotropik kristallerde bulk modülü, K = (3r vB2 - 4r vE2 ) / 3                                                                               (19)

Enine (shear) modülü, m = r vE2                                                                                                                            (20)

Young modülü, Y = r vE2 (3r vB2 - 4r vE2 ) / (r vB2 - r vE2 )                                                                               (21)

Poisson oranı, s = ½ (r vB2 - 2r vE2 ) / (r vB2 - r vE2 )                                                                                       (22)

bağıntıları ile verilir ( İşçi, 1983).

3. YÖNTEM

    Daha önce yapılan deneylerden elde edilen veriler değerlendirilerek ( İşçi vd...,1978), farklı yapıya sahip kristaller için nümerik hesaplar yapıldı. Akustik hız, Young Modülü, sıkışabilirlik gibi fiziksel parametrelerin farklı kristal doğrultuları için haritaları oluşturulabilir.

    Bazı özel doğrultularda nikel tek kristalinde yapılan ultrases hızı ölçümleri sonucu, kübik kristallerdeki 3 bağımsız esneklik sabiti, C11, C12 ve C44 hesaplanmıştır (İşçi ,1983).

Çizelge 2 Nikel tek kristalinde farklı doğrultularda ölçülen ultrases hızları

Yayılma doğrultusu

Polarizasyon doğrultusu

Oda sıcaklığında hızın değeri ( 103 m/s )

[110]

[110]

5.97

[110]

[001]

3.73

[110]

[110]

2.32

[100]

[100]

5.24

[100]

(100) düzleminde

3.69

[111]

[111]

6.17

[111]

(111) düzleminde

2.84

    Bu değerlerden hareket edilerek (001) ve (011) düzlemlerinde 5 derece aralıklarla, farklı doğrultuda 1 boyuna ve 2 enine dalga hızları hesaplanarak, grafikleri çizilmiştir. Bu grafikler Şekil 1 ve 2'de verildiği gibi (001) düzleminde 4-katlı simetri, (011) düzleminde 2-katlı simetri göstermiştir. Fiziksel bir özellik olan esneklik, kristal simetrisi ile aynı simetriyi göstermiştir( Nye, 1957). Nümerik hesaplamalar Mathematica 2.2 da yapılmıştır (Wolfram, 1991). Sonuçlar ve programlar aşağıda verilmiştir. Grafikler ayrıca Mathematica'da PolarPlot ile çizilmiştir. Bu çalışmada iki ayrı program kullanılmıştır. İlk programla, kübik kristallerde (001) ve (011) düzlemlerinde ultrases hızı yüzeyleri oluşturulmuş; ikinci program ile, Young modülü, sıkışabilirlik katsayısı, bulk (hacim) modülü gibi esneklik sabitleri hesaplanmıştır.

4. SONUÇLAR VE ÖNERİLER

    Bazı özel doğrultularda ( Çizelge 1 ) yapılan ultrases hızı ölçümleri sonucu, kübik kristallerdeki 3 bağımsız esneklik sabiti, C11, C12 ve C44 hesaplanmıştır. Ni için bu değerler aşağıya çıkarılmıştır:

Yoğunluk (r ) = 8.90x103 kg/m3

C11 = 24.65x1010 N/m2

C12 = 14.73x1010 N/m2

C44 = 12.47x1010 N/m2

    İlk programın sonuçları aşağıda verilmiştir. Grafik 5'er derece aralıklarla hesaplanarak çizdirilmiştir. Nümerik sonuçlardan 30 derece aralıklı olanlar aşağıda listelenmiştir.

                 Çizelge 3. Nikel tek kristalinde farklı doğrultularda hesaplanan ultrases hızları (103 m/s)

Açı

v1

v2

v3

va1

va2

va3

0

5.26

3.74

3.74

5.26

3.74

3.74

30

5.87

3.74

2.68

5.94

3.45

2.91

60

5.87

3.74

2.68

6.23

3.04

2.77

90

5.26

3.74

3.74

6.01

3.74

2.36

120

5.87

3.74

2.68

6.23

3.04

2.77

150

5.87

3.74

2.68

5.94

3.45

2.91

180

5.26

3.74

3.74

5.26

3.74

3.74

210

5.87

3.74

2.68

5.94

3.45

2.91

240

5.87

3.74

2.68

6.23

3.04

2.77

270

5.26

3.74

3.74

6.01

3.74

2.36

300

5.87

3.74

2.68

6.23

3.04

2.77

330

5.87

3.74

2.68

5.94

3.45

2.91

360

5.26

3.74

3.74

5.26

3.74

3.74

                                                                                  

                                                                          (a)                                                       (b)

                Şekil 1.  Ni tek kristalde ultrases hız yüzey kesitleri ( 103 m/s ).

                                (a) (100 ) düzleminde: 1, enine dalga (v3); 2, enine dalga (v2); 3, boyuna dalga (v1).

             (b) ( 110 ) düzleminde: 1, enine dalga (va3); 2, enine dalga (va2); 3, boyuna dalga (va1).

    İkinci programdan elde edilen sonuçlar aşağıda listelenmiş ve grafikleri çizilmiştir:

S44 = 8.02 10-13 m2 / N , S11 = 7.35 10-13 m2 / N , S12 = -2.75 10-13 m2 / N

Lineer Sıkışabilirlik = 1.85 10-13 (SI), Hacim Sıkışabilirliği = 5.55 10-13 (SI)

 

    Çizelge 4.  Nikel tek kristalinde farklı doğrultularda hesaplanan Young modülü değerleri (1012 N/m2)

Açı

YA

YB

0

1.36

1.36

30

1.97

2.05

60

1.97

2.98

90

1.36

2.32

120

1.97

2.98

150

1.97

2.05

180

1.36

1.36

210

1.97

2.05

240

1.97

2.98

270

1.36

2.32

300

1.97

2.98

330

1.97

2.05

360

1.36

1.36

   Şekil 2. Ni tek kristalinde Young modülü yüzeyleri ( 1012 N/m2 ). (a) (100) düzleminde, (b) (110) düzleminde.

    Young modülü Şekil.2'deki grafiğikte görüleceği üzere [111] doğrultusunda, Young modülü en büyük değere sahiptir. Yani Ni gibi kübik kristallerde, bu doğrultuda atomlar arası bağlar daha kuvvetlidir. Cisim bu doğrultuda gerilim ve basınca daha fazla dayanmaktadır. Bu sabitin küçük olduğu doğrultularda, maddenin deformasyonu daha kolaydır (Smile, 1995).

    Bu konuda en tipik örnek HCP yapıya sahip grafittir. Grafitin erime noktası ( 3000 oC civarı) çok yüksek olmasına rağmen, belirli düzlemler arası bağlar çok zayıf olması nedeni ile kolayca işlenip pota veya elektrotlar yapılabilir (İşçi, 1978).

    Katı cisimlerin esneklik özelliklerinin incelenmesi ile atomlar ve moleküller arasındaki bağlayıcı kuvvetler hakkında bilgi elde edilebilir. Tek-kristaller esneklik açısından anizotropik özellik gösterirler. Farklı doğrultularında ultrases hızı ve esneklik sabitlerinin ölçülmesi ile, esneklik özellikler ve bağlar hakkında yorumlar yapılabilir.

5. KAYNAKLAR

AKGÖZ, Y.C.; İŞÇİ, C; SAUNDERS, G.A. (1976): "The Elastic Constants of Sb-As Alloy Single Crystals", J. Material Science, V.11, p.297.

İŞÇİ, C; PALMER, S.B. (1978): "An Ultrasonic Study of the Magnetic Phases of Dy", J.Physics F: Metal Physics, V.8, p.247.

İŞÇİ, C.(1977): "Doktora Tezi", Hull University, England .

İŞÇİ, C. (1978): " Hexagonal Kristallerin Esneklik Özellikleri", E.Ü. Fen Fak. Dergisi,

V.3, p.279, İzmir.

İŞÇİ, C (1983) : " Ni ve Ni-V Alaşımlarının Magnetik ve Esneklik Özellikleri",

Doçentlik Tezi, Ege Üniv. Fizik Bölümü, İzmir.

NYE, J.F. (1957): "Physical Properties of Crystals: Their representation by Tensors and Matrices", Oxford Universty Press, London.

SMILE, R.W. (1995): "None Destructive Testing" PH Diversfield Inc., New York.

TRUELL, R.; ELBAUM, C; CHICK, B.B. (1969): "Ultrasonic Methods in Solid State Physics", Academic Press, New York.

WOLFRAM, S. (1991): "Mathematica" 2nd ed.; Addision-Wesley Com. Inc., New York.