SONLU ELEMANLAR DERS 2

1. SONLU ELEMANLAR DERS 2

2. Bir problemi sonlu elemanlarla formüle etmede üç yaklaşım vardır. • Direk formülasyon • Minimum toplam potansiyel enerji formülasyonu • Ağırlıklı hata formülasyonu

3. DİREK FORMÜLASYON Direk formülasyonu bir örnekle açıklayalım: ÖRNEK: w1 Yanda boyutları verilen değişken kesitli çubuk üst kenarından sabitlenmiş olup alt kenarından P kuvveti ile yüklenmektedir. Elastisite modülü E olan bir malzemeden üretilmiş olan bu çubuğun uzunluğu boyunca farklı noktalarda ne kadar deplasman yaptığını bulunuz. Çözümde çubuğun ağırlığını ihmal ediniz. L y w2 P

4. ÖN İŞLEM (PREPROCESSING) • İlk önce model elemanlara ve düğümlere bölünür. • Sonucun hassasiyetini arttırmak için eleman ve düğüm sayısı arttırılır. • Bu örnekte 5 düğüm ve 4 eleman kullanılacaktır.

5. 1.Çubuğun eleman ve düğümlere ayrılması 1 1 l1 1.eleman u1 A1 2 2 u2 2.eleman l2 A2 3 3 u3 3.eleman l3 A3 4 4 u4 4.eleman l4 A4 5 5 u5 P P P

6. 2. Bir eleman için çözüm geliştirilir. Bu eleman için ortalama gerilme: l Bu eleman için ortalama birim uzama: Dl Kesit alanı A F keş keş x Bu ifade F=k.x e benziyor F olarak bulunur. O halde

7. Yani modeli, kesit alanları farklı 4 yayın ucuca eklenmesi, bir ucunun sabitlenerek diğer ucundan kuvvetin uygulanması şeklinde düşünebiliriz. O halde

8. Tüm düğümler üzerine etkiyen kuvvetleri gösteren serbest cisim diyagramını çizersek R1 Statik dengeden her bir düğüme etkiyen kuvvetlerin toplamının sıfır olması gereklidir. Buna göre 1. düğüm k1(u2-u1) k1(u2-u1) 2. düğüm k2(u3-u2) k2(u3-u2) 3. düğüm k3(u4-u3) k3(u4-u3) 4. düğüm k4(u5-u4) k4(u5-u4) 5. düğüm P

9. Bu eşitliği kuvvet terimlerinin eşitliğin diğer tarafına atarak tekrar düzenlersek: Bu denklemleri matris formunda yazarsak:

10. Matris ifadesinde uygulanan kuvvet ve reaksiyon kuvveti ayırarak düzenleme yaparsak: Bu ifadeyi genel bir şekilde yazarsak:

11. Sınır şartları uygulanırsa 1 nolu düğümün deplasmanı sıfırdır. Bu nedenle global direngenlik matrisinde u1 e karşılık gelen satır ve sütun 0 olurken köşe eleman 1 olur. Yani:

12. Bu yeni matris formu Gauss eliminasyon yöntemi ile çözülerek düğümlerdeki deplasman değerleri bulunur. Deplasman değerleri bulunduktan sonra önceki bağıntılardan reaksiyon kuvvetleri bulunabilir. Şimdi bir eleman için direngenlik matrisi ve global direngenlik matrisinin nasıl oluşturulduğunu inceleyelim.

13. 3. Bir eleman için denklemler çıkarma fi=keş(ui+1-ui) • Problemdeki her bir elemanda iki düğüm vardır. • İki düğümle alakalı olarak tek yöne hareket var. • O halde deplasmanları bakımından iki bilinmeyen var. • İki bilinmeyeni bulmak için iki denkleme ihtiyaç var. i. düğüm ui i+1.düğüm Ui+1 fi+1=keş(ui+1-ui) veya fi=keş(ui-ui+1) • Yanda görüldüğü gibi iç kuvvetleri gösterebiliriz. • Statik denklemlerinden iç kuvvetlerin toplamının sıfır olduğunu söyleyebiliriz. i. düğüm ui i+1.düğüm Ui+1 fi+1=keş(ui+1-ui)

14. Deplasman matrisi Elemanın direngenlik matrisi İç kuvvetler matrisi

15. 4. Sistemin tümünü temsil eden global direngenlik matrisinin bulunması • Problemimizde 4 tane eleman ve 5 tane düğüm vardı. • Serbestlik derecesi düğümün kaç yönde hareket ettiğini gösterir. • Bizim problemimizde her düğüm sadece 1 yönde hareket etmektedir. O halde her düğüm bir serbestlik derecesine sahiptir. • Global K matrisinin boyutları bu tanıma göre düğüm sayısı.serbestik derecesi x düğüm sayısı.serbestlik derecesi • Sonuç olarak problemimizde K matrisi 5x5 boyutundadır.

16. 1. elemanın direngenlik matrisi

17. 2. elemanın direngenlik matrisi

18. 3. elemanın direngenlik matrisi

19. 4. elemanın direngenlik matrisi

20. Herbir elemanın glabal direngenlik matrisinin içinde bulunduğu matrisler toplanırsa

21. Sınır şartları uygulanırsa 1 nolu düğümün deplasmanı sıfırdır. Bu nedenle global direngenlik matrisinde u1 e karşılık gelen satır 0 olurken köşe eleman 1 olur. Yani:

22. Bu yeni matris formu Gauss eliminasyon yöntemi ile çözülerek düğümlerdeki deplasman değerleri bulunur. Deplasman değerleri bulunduktan sonra önceki bağıntılardan reaksiyon kuvvetleri bulunabilir.

23. Reaksiyon kuvvetlerinin hesabı 1. YOL: 2. YOL:

24. Gerilmelerin hesaplanması

25. ÖRNEK 2- ISI PROBLEMİ 6 5 4 3 1 2 Bir evin dış duvarı, yukarıdaki tabloda gösterilen malzemeleri içermektedir. Odadaki sıcaklık Tiç=70 ºF olup dış ortam sıcaklığı Tdış=20 ºF dir. Duvarın kesit alanı 150 ft2 olduğuna göre duvar boyuncaki sıcaklık dağılımını bulunuz.

26. ÖN İŞLEM (PREPROCESSING) • İlk önce model elemanlara ve düğümlere bölünür. • Sonucun hassasiyetini arttırmak için eleman ve düğüm sayısı arttırılır. • Bu örnekte 7 düğüm ve 6 eleman kullanılacaktır.

27. 1.Problemin eleman ve düğümlere ayrılması 3.eleman 4.eleman 5.eleman 2.eleman T7=70ºF T1=20ºF 2 6 7 4 1 5 3 (1) (2) (3) (4) (5) (6)

28. 2. Bir eleman için çözüm geliştirilir. Isı transferinin iletim ve taşınım olarak iki tipi vardır. (2), (3), (4) ve (5) nolu elemanlar iletim mevcut olup düzenli rejim halindeki termal davranışları Fourier Yasası kullanılarak modellenebilir. Enerji yüksek sıcaklıktan düşük sıcaklığa doğru moleküler aktivite tarafından taşınmaktadır. Isı transferi oranı olarak tanımlanabilecek ısı akısı, Fourier yasası ile Isı akışının sıcaklığın azalması yönünde olduğunu gösterir alan sıcaklık gradyanı ısı akısının x bileşeni ısı iletim katsayısı

29. Eşitlik farklı bir formda şu şekilde yazılabilir: l k U ısı iletim faktörü diye adlandırılırken birim alandan geçen ısı geçişini gösterir. Isı direncinin tersidir ve U=k/l dir. O halde: Ti+1 qx Ti x

30. (1) ve (6) nolu elemanların ısı davranışı Newton’un soğuma yasası kullanılarak modellenebilir. Taşınım ısı transferi, hareket halindeki bir akışkanın farklı sıcaklıktaki bir yüzeye temas etmesi sonucu olur. Newtonun soğuma yasası ile verilen ısı akısı; ısı taşınım katsayısı alan ısı akısı Akışkanın sıcaklığı yüzey sıcaklığı

31. U faktörü kullanılarak bu denklem farklı bir formda yazılabilir.U=h dır. O halde; Düzenli rejim koşullarında yüzeydeki enerji dengesi, yüzeyden, iletim ve taşınımla olan ısı akılarının birbirine eşit olması gerektiğini vurgular. l qtaşınım qiletim Ti+1 Ti=Ts k Tf x

32. Duvarın dış yüzeyinde iletimden dolayı olan ısı kaybı taşınımla olan ısı kaybına eşit olmalıdır. Duvarın içerisinde ise Duvarın iç yüzeyinde iletimden dolayı olan ısı kaybı taşınımla olan ısı kaybına eşit olmalıdır.

33. Bilinen sıcaklıklarla (T1=20 ºF ve T7=70 ºF) alakalı olan terimleri bilinmeyen sıcaklıklarla alakalı olan terimlerden ayırırsak:

34. Bunu matris formunda yazarsak Global direngenlik matrisi=[K]G Burada başlangıç koşulları uygulanmış durumdadır. Reaksiyon kuvvetleri de söz konusu olamaz. Gauss eliminasyon yöntemi ile bu matris çözülebilir.

35. 3. Bir eleman için denklemler çıkarma • Problemdeki her bir elemanda iki düğüm vardır. • Her düğümün bir serbestlik derecesi vardır. • Her eleman için iki denklem çıkarılabilir. • Bu denklemler düğümlerdeki sıcaklıkları ve eşdeğer direngenlik matrislerini içermelidir.

36. İletim için ısı akıları Matris formunda sıcaklık matrisi İletim için direngenlik matrisi Isı akısı matrisi Matris formunda Taşınım için ısı akıları sıcaklık matrisi taşınım için direngenlik matrisi Isı akısı matrisi

37. 1. elemanın direngenlik matrisi

38. 2. elemanın direngenlik matrisi

39. 3. elemanın direngenlik matrisi

40. 4. elemanın direngenlik matrisi

41. 5. elemanın direngenlik matrisi

42. 6. elemanın direngenlik matrisi

43. Herbir elemanın glabal direngenlik matrisinin içinde bulunduğu matrisler toplanırsa

44. Sınır şartları uygulanırsa 1 nolu düğümün ve 7 nolu düğümün sıcaklıkları biliniyor. Bunla alakalı düzenlemeler yapılırsa.

45. Matris düzenlemesi yapılırsa Bu matrisi Gauss eliminasyon metodu ile çözebilir ve bilinmeyen terimleri elde edebiliriz.