Ana Sayfa / Yazılar / Sonlu Elemanlar FEM-03
Gerilme Yerdeğiştirmeden Yavaş Yakınsar
Delikli levhada gerilme yığılması, tam çözümü bilinen bir problem. Kaba ağ tepe gerilmeyi yarıya yakın düşük gösteriyor — yani emniyetsiz tarafta. Ve aynı eleman sayısıyla yalnızca yerleşimi değiştirmek hatayı %48'den %1'e indiriyor.
1. Neden bu problem?
Delikli levha, sonlu elemanlar öğretiminin klasik problemidir ve bunun iyi bir sebebi var: tam çözümü bilinir. Dolayısıyla sayısal sonucun ne kadar yanlış olduğu tahmin edilmez, ölçülür.
Ama seçilmesinin asıl sebebi başka. Önceki iki yazıda hep yerdeğiştirmeye baktık. Mühendislik kararları ise çoğunlukla gerilmeye dayanır — ve ikisi aynı hızda yakınsamaz.
2. Kirsch çözümü
Sonsuz bir levhada dairesel bir delik ve tek eksenli çekme için, delik çeperindeki teğetsel gerilme:
İkinci değer daha az bilinir ve bu yüzden değerlidir: bağımsız bir kontrol sağlar. Bir modelin tepe gerilmeyi doğru bulması tesadüf olabilir; aynı anda yandaki basıncı da doğru bulması olamaz.
3. Model ve simetri
Problem iki eksene göre simetrik olduğundan çeyrek model yeterlidir. Simetri sınır koşulları basittir: her simetri ekseninde, eksene dik yerdeğiştirme sıfırdır.
Delik yarıçapı a = 10, çeyrek modelin yarı genişliği W = 200 alındı; a/W = 0,05 olduğundan sonlu genişlik etkisi küçüktür ve Kirsch'in sonsuz levha varsayımına yakınız.
4. Yakınsama
| Ağ | Eleman | Serbestlik | Yığılma katsayısı | Hata |
|---|---|---|---|---|
| hesaplanıyor… | ||||
5. Neden gerilme daha yavaş?
Sonlu eleman yöntemi yerdeğiştirmeyi bilinmeyen olarak alır. Gerilme, yerdeğiştirmenin türevinden hesaplanır:
FEM-02'deki konsolla karşılaştırın: 256 elemanda yerdeğiştirme hatası %4 idi. Burada 512 elemanda gerilme hatası hâlâ —. Aynı mertebede bir ağ, gerilmede belirgin biçimde daha kötü sonuç veriyor.
Bunun pratik karşılığı şudur: yerdeğiştirmesi yakınsamış görünen bir model, gerilme bakımından hâlâ yakınsamamış olabilir. Sehim kontrolü için yeterli bir ağ, gerilme kontrolü için yeterli değildir.
6. Hangi yönden yakınsıyor — ve neden önemli
FEM-02'de yerdeğiştirmenin daima aşağıdan yakınsadığını, yani modelin gerçek yapıdan rijit olduğunu görmüştük. Burada tepe gerilme de aşağıdan yakınsıyor — ama sonucu çok daha ciddi.
Önemli bir uyarı: yerdeğiştirme için geçerli olan "daima aşağıdan" kuralı gerilme için genel bir teorem değildir. Bu problemde öyle davranıyor, çünkü tepe gerilme keskin bir gradyanın tepesindedir ve kaba ağ onu yumuşatır. Başka problemlerde — özellikle tekillik civarında — gerilme yukarıdan da yaklaşabilir, hatta sınırsız büyüyebilir.
7. İkinci, bağımsız kontrol: basınç
Yük doğrultusundaki delik çeperinde Kirsch −1 öngörüyor: çekme altındaki bir levhada, deliğin yanında basınç gerilmesi oluşur.
| Ağ | σ(θ=0) / σ∞ | Hata |
|---|---|---|
| hesaplanıyor… | ||
Bu değerin de doğru yakınsaması, modelin yalnızca bir sayıyı tutturmadığını, gerilme alanının bütününü doğru temsil ettiğini gösterir.
8. Elemanın yeri, sayısından önemli
Şimdi yazının en pratik bulgusu. Aşağıdaki tabloda ağ hiç değişmiyor — hep 32×16, yani 512 eleman. Değişen tek şey, elemanların delik çeperine doğru ne kadar sıklaştırıldığı.
| Kademelendirme | Eleman | Yığılma katsayısı | Hata |
|---|---|---|---|
| hesaplanıyor… | |||
Bu, serinin üç yazısında da tekrarlanan temanın en keskin hâli:
- FEM-01: doğru yük vektörü, daha çok elemandan iyidir.
- FEM-02: doğru eleman ve integrasyon seçimi, daha çok elemandan iyidir.
- FEM-03: doğru eleman yerleşimi, daha çok elemandan iyidir.
Üçünde de "ağı sıklaştır" refleksi en pahalı ve en az etkili çözüm oldu. Ağ sıklaştırmak işe yarar, ama önce doğru soruyu sormak gerekir: hata nereden geliyor?
9. Motorun kendi sınırı
Bu motor yoğun bir doğrusal denklem çözücüsü kullanır. Bellek serbestlik derecesinin karesiyle, süre küpüyle büyür. Öğretim amaçlı modeller için yeterlidir; daha büyük ağlarda motor çökmek yerine açıkça reddeder ve sebebini söyler.
10. Pratikte ne yapmalı
- Gerilme için ayrı bir yakınsama incelemesi yapın. Sehim yakınsadı diye gerilme yakınsamış sayılmaz.
- Elemanları gradyanın olduğu yere koyun. Delik, çentik, köşe, yük uygulama noktası. Düzgün bölünmüş bir ağ, elemanların çoğunu hiçbir şeyin olmadığı yere harcar.
- Tepe gerilmenin düşük çıktığını varsayın. Yakınsama incelemesi yapmadıysanız, bulduğunuz değer gerçeğin altındadır.
- Mümkünse kapalı-form bir kontrol bulun. Kirsch gibi klasik çözümler, yeni bir modelin ya da yeni bir yazılımın ilk sınavı olmalıdır.
- Birden fazla noktada kontrol edin. Tek bir sayıyı tutturmak tesadüf olabilir; gerilme alanının iki farklı noktasını tutturmak olamaz.
- Simetriyi bilinçli kullanın. Aradığınız davranış simetri kırıyorsa çeyrek model onu göstermez.
11. Kendini sınama
- Gerilme neden yerdeğiştirmeden bir mertebe daha az doğrudur?
- Kaba ağda tepe gerilmenin düşük çıkması neden emniyetsiz taraftır?
- Kirsch çözümünde θ = 0 noktasındaki gerilme neden basınçtır? Çekme uygulanan bir levhada bu nasıl olur?
- 512 elemanlı iki model kurdunuz; biri %48, diğeri %1 hata veriyor. Aradaki tek fark ne olabilir?
- Yerdeğiştirme için "daima aşağıdan yakınsar" diyebiliyoruz. Gerilme için neden aynı şeyi diyemiyoruz?
- Bir modelde sehim yakınsamış görünüyor. Gerilme sonuçlarına güvenmek için ne yaparsınız?
Kaynaklar
- Kirsch çözümü — sonsuz levhada dairesel delik çevresindeki gerilme alanı; klasik elastisite sonucu, ülke yönetmeliğinden bağımsızdır.
- Düzlem gerilme izoparametrik formülasyon ve yakınsama mertebesi — standart sonlu elemanlar sonuçları.
- Ders modülleri ve atölye problemleri için: Kwon, Y.W. & Bang, H., The Finite Element Method using MATLAB. Bu yazıdaki anlatım, türetme ve bütün sayılar kendi motorumuzdandır; kaynak konu kapsamı için atıf olarak verilmiştir.
Yazıdaki bütün sayılar engine/fem2d.js tarafından üretilir ve tests/fem2d.test.js içinde ayrıca sabitlenmiştir (35/35). Doğrulama iki bağımsız noktada yapılmıştır: delik tepesindeki yığılma katsayısı (3,000) ve yan noktadaki basınç (−1,000). Çeyrek model ağı, kullanılan bütün bölüntülerde pozitif Jacobian veriyor ve bu ayrıca test edilmektedir.
Sehim yakınsadı diye gerilme yakınsamış sayılmaz; kaba ağ tepe gerilmeyi düşük gösterir.