DÜZGÜN YAYILI YÜKE MARUZ ORTOTROP İK KOMPOZİT ANKASTRE KİRİŞTE SEHİM HESABI
Alaattin AKTAŞ, Cemil TANIK
Kırıkkale Üniversitesi, Mühendislik Fakültesi, Makine Mühendisliği Bölümü, Kırıkkale
ÖZET
Bu çalışmada T300/976 Grafit-Epoksi ve Karbon-Epoksi malzemeden yapılan, yayılı yüke maruz bir ortotropik kompozit ankastre kiriş göz önüne alınmıştır. Bu kirişin sehiminin analitik ve sonlu elemanlar çözümü yapılmış, bu çözümler karşılaştırılmıştır. Sonlu elemanlar yönteminin uygulanmasında Ansys 50a programı kullanılmıştır.
Hesaplamalar fiber doğrultusu 0°, 15°, 30°, 45°, 60°, 75°, 90° için yapılmıştır.
Anahtar Kelimeler : Sehim, Ankastre kiriş, Anizotropik levha
CALCULATING DEFLECTION OF THE ORTHOTROPIC COMPOSITE CANTILEVER SUBJECTED TO A UNIFORMLY DISTRIBUTED LOAD
ABSTRACT
In this study, an orthotropic composite cantilever, which is made of T300/976 Grafit-Epoxy and Carbon-Epoxy materials subjected to a uniformly distributed load, is considered. The deflection of this cantilever is calculated by finite element method and calculated analytically. These results are compared. For finite element calculation, Ansys 50a program is used. The calculations are made for 0°, 15°, 30°, 45°, 60°, 75°, 90° of fiber directions.
Key Words : Deflection, Cantilever, Anisotropic plate
1. GİRİŞ
Günümüzde düşük özgül ağırlığa karşılık yüksek mukavemet ve korozyona dayanıklılığı nedeniyle kompozit malzemeler; uçaklarda, uzay araçlarında ve inşaat yapılarında metal malzemelere nazaran daha fazla tercih edilmektedir (Chang ve Richard 1982; Aktaş ve Karakuzu 1998).
Son zamanlarda yaygın olarak kullanılan kompozit kirişler kullanılma alanı nedeniyle büyük önem kazanmıştır. Bu nedenle birçok araştırmacı kompozit kirişler üzerinde çalışmıştır. Karakuzu ve Özcan (1996) bir ankastre kirişteki gerilmelerin anizotropik elastisite çözümünü farklı yükleme durumları için bulup, kirişin elasto-plastik gerilme analizini yapmıştır. Barbero ve ark., (1993) I profilli bir
kirişin kesme-düzeltme faktörünü ve sehimini sonlu elmanlar metoduyla kirişin kesit boyutlarını değiştirerek hesaplamıştır.
Genellikle kirişlerin kullanıldığı tüm yapılarda kirişlerin çökmesi istenmeyen bir durumdur.
Çökmenin önlenmesi için onun sehim değerinin bilinmesi gerekir. Bu öneminden dolayı, çalışmada bir ankastre kompozit kirişin sehim denkleminin anizotropik elastisite çözümü yapılıp, kirişin maksimum sehimi farklı iki kompozit malzeme kullanılarak bulunmuştur. Elde edilen çözüm aynı malzemeler kullanılarak sonlu elmanlar metoduyla da hesaplanıp, bulunan bu iki çözüm karşılaştırılmıştır. Maksimum sehim fiber doğrultusu 0°, 15°, 30°, 45°, 60°, 75°, 90° için hesaplanmıştır.
2. PROBLEMİN TANIMLANMASI VE KİRİŞİN SEHİMİNİN BULUNMASI
Bu çalışmada düzgün yayılı yüke maruz L uzunluğunda, b yüksekliğinde ve h genişliğinde bir kompozit ankastre kiriş gözönüne alınmıştır (Şekil 1). θ açısı x ekseniyle fiber doğrultusu arasındaki açı olup, gerilme-uzama bağıntısı aşağıdaki gibi verilir (Lekhnitskii, 1968; Jones, 1975).
1 q
h
O z
x θ O b
2 L
y y
Şekil 1. Yayılı yüke maruz kompozit ankastre kiriş
xy 16 y 12 x 11
x =a σ +a σ +a τ ε
xy 26 y 22 x 12
y =a σ +a σ +a τ
ε (1)
xy 66 y 26 x 16
xy =a σ +a σ +a τ γ
Burada;
y , v x u
y
x ∂
=∂
∂ ε
=∂
ε ve
x v y u
xy ∂
+∂
∂
=∂
γ (2)
olarak verilir.
Eğer cos(θ) = m ve sin(θ) = n ile gösterilirse aij
katsayıları aşağıdaki gibi tanımlanır.
4 22 2 2 66 12 4 11
11 s m (2s s )m n s n
a = + + +
2 2 66 22 11 4 4 12
12 s (m n ) (s s s )m n
a = + + + −
4 22 2 2 66 12 4 11
22 s n (2s s )m n s m
a = + + +
3 66 12 22 3 66 12 11
16 (2s 2s s )m n (2s s s )mn
a = − − − − − (3)
3 66 12 22 3 66 12 11
26 (2s 2s s )n m (2s s s )nm
a = − − − − −
) n m ( s n m ) s s 4 s 2 s 2 ( 2
a66= 11+ 22− 12− 66 2 2− 66 4 + 4
Burada;
12 66 2 22 1 12 12 1
11 G
s 1 E , s 1 E , s
E ,
s 1 ν = =
−
=
= (4)
d i r .
Yayılı yüke maruz kompozit ankastre kirişteki gerilmeler aşağıdaki gibi verilmiştir. (Lekhnitskii, 1968; Karakuzu, 1996)
−
+ −
+
− +
−
=
σ 5b
y 4 b
y . 4 a a a
4 a a 2 2 b
y 1 12 b .x a a h q I 2
y qx
3 3 2 11 2 16 11
66 12 2
2
11 16 2
x
− + −
−
=
σy 33
b y 4 b
y 1 3 h 2
q (5)
−
−
−
−
=
τ 3
3
11 16 2
2
xy b
y 4 b y ha y qa 4 b I 2 qx
(5) denklemleri (1) denkleminde yerine konup, bu denklemler
x u
x ∂
=∂
ε denklemiyle eşitlenir ve x’e
göre integre edilirse x yönündeki yerdeğiştirme fonksiyonu,
) y ( g h b
qxy a 2 h b a
qxy a 4
h b
qxy a 2 h b
y qx a 3 h b
y qx a 2 bh 10
qxy a 3 bh a 5
qxy a bh 10
qxy a 9 bh 4
qx a h 2
qx u a
3 3 66 3
11 3 2 16
3 3 12 3
2 2 16 3
3 11 66
11 2 16 12
2 16 12
+ +
−
+
−
−
− +
+
−
−
=
(6)
şeklinde elde edilir.
Benzer düşünceyle (5) denklemleri (1) denkleminde yerine konup, bu denklemler
y v
y ∂
=∂
ε denklemiyle
eşitlenir ve y’ye göre integre edilirse y yönündeki yer değiştirme fonksiyonu
h b
y qx a 3 bh a 20
qy a a 3 bh a 2
qy a a bh 4
qy a 3 bh a 5
qy a a 3 bh a 10
qy a 3 bh 2
qxy a 3 bh a
qxy a a h 2
qy v a
3 2 2 12 11
2 66 12 11
2 26 16 2 22 2
11 2 2 16 12 11
2 2 12 26
11 16 12
22 + − − + + − − −
−
=
f(x)
h b a 2
qy a a h b a
qy a a h b 2
qy a h b a
qy a a 2 h b a
qy a h b
qxy a 2 h b a
qxy a a 4
3 11
4 66 12 3
11 4 26 12 3
4 22 3
2 11
4 2 16 12 3
11 4 2 12 3
3 26 3
11 3 16
12 + + − − + + +
− (7)
olarak bulunur. v’yi bulmak için f(x)’in bulunması gereklidir. Bunun için (1) denklemindeki γxy ile eşitlenirse
bh a 10
qy a a 13 bh
qy a 3 bh a 5
qy a 6 bh a 5
qy a a 8 h b
qx a 2 bh 5
qx a 6 bh a 5
qx a 4 bh 10
qx a 9 h 2
q a
11 66 16 26
2 11 3 16 11
16 12 3
3 11 66
11 2 16 12
26 ) 2 ( xy ) 1 (
xy −γ = + − + − + − − −
γ
f (x) g(y) 0
h b a
qy a a 6 h b
qy a 4 h b a
qy a 8 h b a
qy a a 8
3 11
3 66 16 3
3 26 3
2 11
3 3 16 3
11 3 16
12 + + − + ′ + ′ =
− (8)
bulunur. Bulunan bu son denklem 0 (sıfır)’a eşit olduğundan x’e bağlı terimlerin toplamı ile y’e bağlı terimlerin toplamı ayrı ayrı birer sabit sayıya eşit olması gerekir. c ve d sabitler olmak üzere
c ) x ( f ) x (
A + ′ =
d ) y ( g ) y (
B + ′ = (9) Burada A(x) x’li terimlerin, B(y) y’li terimlerin toplamıdır.
( ) ( )
x f x cA + ′ = (10) olduğuna göre,
[
c A(x)]
dx )x (
f =
∫
− (11) dir. Buradan;d cx h b 2
qx a bh
a 20
qx ) a a 12 a 8 a a 9 ) ( x (
f 3
4 11 11
2 66 11 2 16 12
11 − + + − +
−
= (12)
bulunur. f(x) değeri (7) denkleminde yerine konursa ankastre kirişin denklemi c ve d sabitlerine göre bulunmuş olur.
x = L, y =0 noktasında 0
dxdv = sınır şartı sehim denklemine uygulanırsa c sabiti,
h b a 10
qx ) q L a 20 Lq b a a 12 Lq b a 8 Lq b a a 9
c ( 3
11
2 3 2 11 2
66 11 2
2 16 2
12
11 − − +
−
= (13)
elde edilir. (13) denklemi sehim denkleminde yerine konursa sadece d sabiti bilinmiyen olarak kalır.
Ankastre kirişte x = L, y = 0 noktasında v = 0
olduğundan bu değerler sehim denkleminde yerine konursa:
h b a 20
) q L a 20 q L b a a 12 q L b a 8 q L b a a 9
d ( 3
11
4 2 11 2
2 66 11 2
2 2 16 2 2 12
11 − + +
−
= (14)
bulunur. Bu değer de sehim denkleminde yerine konursa sehim denklemi;
bh 20
qx a 9 h b
qx L a 2 bh 5
Lqx a 6 bh a 5
Lqx a 4 bh 10
Lqx a 9 h b 2
q L a 3 bh 5
q L a 3 bh a 5
q L a 2 bh 20
q L a v 9
2 12 3
3 11 66
11 2 16 12
3 4 11 2
66 11
2 2 16 2
12 + − + + − + − −
−
=
bh a 5
qy a a 3 bh a 10
qy a 3 bh 2
qxy a 3 bh a
qxy a a h 2
qy a h b 2
qx a bh 5
qx a 3 bh a 5
qx a 2
2 11
2 2 16 12 11
2 2 12 26
11 16 12 22
3 4 11 2 66 11
2 2
16 − + − + − − +
+ (15)
h b
qy a h b
qxy a 2 h b a
qxy a a 4 h b
y qx a 3 bh a 20
qy a a 3 bh a 2
qy a a bh 4
qy a 3
3 4 2 12 3
3 , 26 3
11 3 16 12 3
2 2 12 11
2 66 12 11
2 26 16 2
22 − − − − + +
+
h b a 2
qy a a h b a
qy a a h b
qy a h b a
qy a a 2
3 11
4 66 12 3
11 4 26 16 3
4 22 3
2 11
4 2 16
12 − + +
−
olarak bulunur. Ankastre kirişin simetri eksenindeki sehim denklemini bulmak için bu son denklemde
y = 0 konup, gerekli sadeleştirmeler yapılırsa simetri eksenine ait sehim denklemi,
( )
h b 2
) x x L 4 L 3 ( q L a
x a a 4 a 3 a 8 bh 3 20
q v 3
3
4 3 4 2 11
66 11 2 16 12
+ + −
−
− + −
= (16)
olarak bulunur.
3. SONLU ELEMANLAR ANALİZİ
Sonlu elemanlar çözümü ile analitik çözümü karşılaştırmak için L = 150mm, b = 40 mm, h = 1 mm, q = 2N/mm ve fiber doğrultusu sırayla 0°, 15°, 30°, 45°, 60°, 75° ve 90° olarak alınmıştır.
Ankastre kirişin sehim değeri T300/976 Grafit-Epoksi ve Karbon-Epoksi malzemeleri için ayrı ayrı hesaplanmıştır. T300/976 Grafit-epoksi malzemesinin malzeme özellikleri Larry (1995)’den alınmış, CYTEC Aerospace (İngiltere)’den temin edilen karbon epoksi malzemesinin malzeme özellikleri ise çekme cihazı ve strain gauge yardımıyla bulunmuştur (Tablo 1).
Tablo 1. T300/976 Grafit-Epoksi ve Karbon-Epoksi Malzemelerinin Malzeme Özellikleri
E1 (GPa) E2 (GPa) G12 (GPa)ν12
T300/976 Grafit-Epoksi 156 13 7 0.23
Karbon-Epoksi 60.2 59.2 3.91 0.25
Problemin modellenmesi ve çözümü Ansys 50a’da yapılıp modelde 540 eleman 1735 düğüm, eleman olarak da 8 düğümlü Shell 91 elemanı kullanılmıştır (Ansys User’s Manual). Problemin sonlu elemanlar çözümü her iki malzeme için Şekil. 2-15.'te gösterilmiştir.
Sonlu eleman ve anizotropik elastisite çözümü her iki malzeme için x = 0, y = 0 noktasının sehimi fiber açısına göre Şekil 16 ve Şekil 17'de grafik olarak verilmiştir.
Şekil 2. T300/976 Grafit-Epoksi malzemesinden yapılan ankastre kirişin θ = 0° için sonlu eleman modeli ve çözümü
Şekil 3. Karbon-Epoksi malzemesinden yapılan ankastre kirişin θ = 0° için sonlu eleman modeli ve çözümü
Şekil 4. T300/976 Grafit-Epoksi malzemesinden yapılan ankastre kirişin θ = 15° için sonlu eleman modeli ve çözümü.
Şekil 5. Karbon-Epoksi malzemesinden yapılan ankastre kirişin θ=15° için sonlu eleman modeli ve çözümü.
Şekil 6. T300/976 Grafit-Epoksi malzemesinden yapılan ankastre kirişin θ=30° için sonlu eleman modeli ve çözümü
Şekil 7. Karbon-Epoksi malzemesinden yapılan ankastre kirişin θ = 30° için sonlu eleman modeli ve çözümü
Şekil 8. T300/976 Grafit-Epoksi malzemesinden yapılan ankastre kirişin θ = 45° için sonlu eleman modeli ve çözümü
Şekil 9. Karbon-Epoksi malzemesinden yapılan ankastre kirişin θ = 45° için sonlu eleman modeli ve çözümü
Şekil 10. T300/976 Grafit-Epoksi malzemesinden yapılan ankastre kirişin θ = 60° için sonlu eleman modeli ve çözümü
Şekil 11. Karbon-Epoksi malzemesinden yapılan ankastre kirişin θ = 60° için sonlu eleman modeli ve çözümü
Şekil 12. T300/976 Grafit-Epoksi malzemesinden yapılan ankastre kirişin θ = 75° için sonlu eleman modeli ve çözümü
Şekil 13. Karbon-Epoksi malzemesinden yapılan ankastre kirişin θ = 75° için sonlu eleman modeli ve çözümü
Şekil 14. T300/976 Grafit-Epoksi malzemesinden yapılan ankastre kirişin θ=90° için sonlu eleman modeli ve çözümü
Şekil 15. Karbon-Epoksi malzemesinden yapılan ankastre kirişin θ = 90° için sonlu eleman modeli ve çözümü
4. SONUÇLAR
Şekil 16. (E1>E2) ve Şekil 17. (E1 ≅ E2)’de kirişin ucundaki sehim fiber açısına bağlı olarak verilmiştir.
Her iki durumda sonlu eleman ve anizotropik elastisite çözümü birbirine yakın bulunmuştur.
0 0,5 1 1,5 2 2,5
0 15 30 45 60 75 90
Fiber Açısı
Maksimum Yerdeğiştirme (mm)
Analitik Ansys50a
Şekil 16. T300/976 Grafit-Epoksi kirişin analitik ve sonlu eleman çözümünün karşılaştırılması
0 0,5 1 1,5 2
0 15 30 45 60 75 90
Fiber Açısı
Maksimum Yerdeğiştirme (mm
Analitik Ansys50a
Şekil 17. Karbon-Epoksi kirişin analitik ve sonlu eleman çözümünün karşılaştırılması
E1 > E2 durumunda θ = 90°’de eğilme rijitliği düşük olduğundan bu açıda sehim maksimum olmaktadır.
E1 ≅ E2 olması durumunda θ = 0° ve θ = 90°’lerde eğilme rijitliği yaklaşık aynı olduğundan kirişin sehimi de yaklaşık aynıdır. Buna karşılık θ=45°’de kirişin eğilme rijitliği en düşük seviyesindedir. Bu yüzden sehim bu açıda maksimum olmaktadır.
Kompozit ankastre kiriş kullanılan yapılarda sehimin istenilmediği durumlarda, E1 > E2 için θ = 0° ve E1 = E2 olduğu durumlarda ise θ = 0°, θ = 90°
alınmalıdır.
5. KAYNAKLAR
Aktaş, A. and Karakuzu, R. 1998. “Mechanics of Composite Materials and Structures NATO Advanced Study ınstitute. Portugal” Failure Analysis of Two dimensional Corbon-Epoxy Composite Plate Pinned-Joint, 12-24 July 1998 Portugal, Vol. 3, 381-397
DeSalvo, G. J. And R. W. Gorman. 1989. ANSYS User’s Manual. Revision 50a Houston, PA: Swanson Analysis Systems, Inc.
Ever J. Barbero, Roberto Lopez-Anido, ve Julio F.
Davalos. 1993. On the Mechanics of Thin-walled Laminated Composite Beams. Journal of Composite
Materials. 27. 806-829.
Fu-Kuo Chang, Richard A. Scott ve George S.
Springer. 1982. Strength of Mechanically Fastened Composite joints Journal of Composite Materials.
16. 470-494.
Jones, R. M. 1975. Mechanics of Composite Materials. McGraw-Hill Kogakusha, Tokyo.
Karakuzu, R. Ve Özcan, R. 1996. Exact Solution of Elasto-Plastic Stresses in a Matel Matrix Composite Beam of Arbitrary Orientation Subjected to Transverse Loads. Composites Science and Technology. 56. 1383-1389
Larry, B. L., Mahmood, M. 1995. Two Dimensional Modeling of Composite Pinned-Joint Failure.
Journal of Composite Materials.
Lekhnitskii, S. G. 1968. Anisotropic Plates. Gordon and Breach Science Publishers. New York.