Ulusal Mekanik Kongresi,03-07 Eylül 2007,ISPARTA
ZAMAN İNTEGRASYON METODLARININ DİNAMİK KOMPOZİT PROBLEMLERE UYGULANMASI
İbrahim H. Güzelbey, Ahmet Erkliğ
Gaziantep Üniversitesi, Makine Mühendisliği Böl., 27310 Şehitkamil, Gaziantep email: guzelbeyih@gantep.edu.tr
ÖZET
Dinamik kompozit problemlerinin, sonlu elemanlar analizi dolaysız integrasyon yöntemleri olarak bilinen Newmark, Wilson-θ, Merkezi Fark ve Houbolt yöntemleriyle gerçekleştirilmiştir. Dinamik kompozit problemleri çözebilecek şekilde iki boyutlu sonlu elemanlarla çözüm programı hazırlanmıştır. İki boyutlu, farklı fiber yönlenme açısına sahip dinamik kompozit problemler çözülmüştür. Sonuçlar ANSYS sonlu elemanlar paket programı ile karşılaştırılmıştır.
ABSTRACT
A finite element analysis of dynamic composite problems has been carried out. Newmark, Wilson-θ, Central Difference and Houbolt Methods have been used as direct integration methods. Two dimensional finite element program has been written for the solution of dynamic composite problems. Two cases have been carried out for symmetric layered composite. The results have been compared with ANSYS finite element package
1.GİRİŞ
Kompozit malzemeler, üretim teknolojilerinin gelişmesi ve geleneksel malzemeler kullanılarak ulaşılabilecek sınırların mühendisliğin birçok alanında sonlara dayanması sonucu giderek artan bir şekilde mühendislik uygulamalarında yapı malzemesi olarak yerlerini sağlamlaştırmaktadır. Birçok alanda geleneksel malzemelerin yerini almaya başlayan kompozit malzemeler özellikle yüksek mukavemet özellikleri ve hafiflikleri ile son dönemlerde tasarımcıların vazgeçilmez araçlarından biri olmuşlardır. Bu çalışmada iki boyutlu kompozit çubuğun dinamik analizi sonlu elemanlar metodu kullanılarak yapılmıştır.
Sonlu elemanlar dinamik problemlerinin çözümü dolaysız integrasyon metotları olarak anılan bazı sayısal yöntemler ile gerçekleştirilebilir [1,2]. Dolaysız integrasyon yöntemlerinde hareket denklemi sayısal integrasyon ile adım adım çözülür. Çeşitli dinamik temas, darbe ve elasto-dinamik problemlerin çözümünde bu yöntemler kullanılmıştır [1,3,4,5]. Bunların, en ok kullanılanları Newmark, Wilson-θ, Merkezi Farklar and Houbolt yöntemleridir.
Bu çalışmada dolaysız integrasyon metotları, iki boyutlu, dinamik, kompozit çubuk problemlerinin çözümü için iki boyutlu sonlu elemanlar programına uygulanmıştır. Farklı fiber yönlenme açısına sahip dinamik kompozit çubuk problemi çözülmüş ve sonuçlar ANSYS sonlu elemanlar paket programı ile karşılaştırılmıştır.
2.DİNAMİK ÇÖZÜM ALGORİTMASI
Sonlu elemanlar dinamik çözümünün sistematiği Şekil2.1’de verilmiştir. Bu sistematikte de görüldüğü üzere değişik zaman integrasyon yöntemlerinin içerisinde sonlu eleman çözümleri yapılmıştır.
Şekil2.1 Dinamik problemler için genel sayısal hesap algoritması
t = t + Δt
zaman < toplam zaman Dur
Kompozit Malzeme Özellikleri Girildi Kütle, Sönüm ve Rijitlik Matrislerini Hesapla Sınır durumlarını uygula ve matrisleri indirge
Etkili Yük ve Rijitlik Matrislerini Hesapla Çöz: KˆΔδt+Δt =Fˆ Başla H E Za man Art ım ı Döngüsü Tepki Kuvvetlerini ve Yerdeğiştirmeleri Değiştir İvme, Hız ve Yerdeğiştirmeleri Hesapla
İvme, Hız ve Yerdeğiştirmeleri Yeni Zaman Artımı İçin Güncelle
Tepki Kuvvetlerini ve Bilinmeyen Yer değiştirmeleri
Hesapla
Koordinatları Güncelle
3.İKİ BOYUTLU KOMPOZİT ÇUBUK 3.1.Basit mesnetli kompozit çubuk
T300/934 grafit/epoksi kompozit çubuğun boyutları ve yükleme durumu Şekil 3.1 ‘de verilmiştir. Çubuk, 55/-55/-55/55 simetrik oryantasyonlu dört katmandan oluşmaktadır. Çubuk ortotropik malzemeden yapılmıştır ve mekanik mukavemet değerleri;
Elastisite Modülü E1=54 Mpa, E2 =18 Mpa =E3
Rijitlik Modülü G12 =8,6 Mpa = G13 =G23
Poisson oranı υ12 =0,25 υ21 =0,25 υ13 = 0,25
Zaman artımı Δt = 0.001
A noktasının zamana bağlı olarak x eksinindeki gerilme grafiği Şekil 3.2’de verilmiştir. B noktalarının sehim y eksinindeki grafiği Şekil 3.3’te verilmiştir. Sonuçlar iki boyutlu FE programında farklı zaman integrasyon metotları ile ve ANSYS ile çözülmüştür. Zaman artımı, Δt = 0.001 saniye olarak seçilmiştir.
Şekil 3.1 Sonlu elemanlar modeli
0.1 m L=1 m 1000 N/m 0.1 m y x B A
-2.0E-03 -1.8E-03 -1.6E-03 -1.4E-03 -1.2E-03 -1.0E-03 -8.0E-04 -6.0E-04 -4.0E-04 -2.0E-04 0.0E+00 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 Zaman (sn) Y ek s ini nd ek i s e h im (m ) Ansys Merkezi Fark Houbolt Newmark Wilson
Şekil 3.2 B noktasının y yönündeki sehimin zamana göre değişimi
-7.00E+05 -6.00E+05 -5.00E+05 -4.00E+05 -3.00E+05 -2.00E+05 -1.00E+05 0.00E+00 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 Zaman (sn) X ek s in dek i ger il m e ( P a) Ansys Merkezi Fark Houbolt Newmark Wilson
3.2.Eksenel yönde kompozit çubuğa kuvvet uygulanması
Bu çalışmada, eksenel simetrik olan, T300/934 grafit/epoksi kompozit çubuğun x ekseninden uygulanan kuvvet altındaki durumu ele alınmıştır. Çubuk, 55/-55/-55/55 simetrik oryantasyonlu dört katmandan oluşmaktadır. Sonlu eleman modeli Şekil 3.4’te görülmektedir. Malzeme özellikleri bir önceki örnektekiyle aynı alınmıştır. A noktasının zamana göre x yönündeki uzama değişimi Şekil 3.5’da görülmektedir.
Şekil 3.4 Sonlu elemanlar modeli
0.0E+00 1.0E-06 2.0E-06 3.0E-06 4.0E-06 5.0E-06 6.0E-06 7.0E-06 8.0E-06 9.0E-06 1.0E-05
0.0E+00 1.0E-03 2.0E-03 3.0E-03 4.0E-03
Zaman (sn) X y önündek i uz am a ( m ) Ansys Merkezi Fark Houbolt Newmark Wilson
Şekil 3.5 A noktasının x yönündeki uzamasının zamana göre değişimi L = 1 m
F 0.1 m
F A
4.SONUÇLAR VE DEĞERLENDİRMELER
Kompozit çubuğun dinamik analizi iki boyutlu sonlu elemanlar programı (FEM) yardımıyla ve farklı doğrudan integrasyon yöntemleri ile gerçekleştirildi. Farklı doğrudan integrasyon yöntemleri ile yapılmış olan sonuçlar ANSYS sonuçları ile birebir örtüşmektedir. Bu durum hazırlanmış olan programın kompozit malzemelerin dinamik analizinde başarılı şekilde çalıştığını göstermektedir.
Bu çalışmanın amacı, daha sonra hazırlanacak 3 boyutlu dinamik kompozit ve temas darbe problemlerini için bir altyapı oluşturmaktır. Yazılmış olan iki boyutlu sonlu elemanlar programı, 3 boyutlu problemleri çözebilecek şekle çevrilecek ve temas darbe problemlerini çözebilecek duruma getirilmeye çalışılacaktır.
KAYNAKLAR
[1] Bathe K.J. and Wilson E.L., “Numerical Methods in Finite Element Analysis”, Prentice-Hall, 1976
[2] Guzelbey I.H., Erklig A., The Efficiency of Direct Integration Methods in Elastic Contact-Impact Problems, Acta Mechanica Sinica, Vol 4, 395-401, 2005
[3] Zienkiewicz O.C., The Finite Element Methods in Engineering. 2nd ed. Mc Graw-Hill, 1971
[4] Hughes T.J., “The Finite Element Methods:Linear Static and Dynamic Finite Element Analysis”, Prentice Hall, 1987
[5] Chen W.H., Yeh J.T. Three-Dimensional Finite Element Analysis of Static and Dynamic Contact Problems with Friction. Comp. & Struc, 35 (5), 541-552, 1990