Ana Sayfa / Yazılar / Sonlu Elemanlar FEM-05
Dört Eleman İki Yüz Elli Altıyı Yenerse
FEM-02'de Q4 elemanı kaba ağda on kat yanılıyordu. İkinci dereceden Q8 elemanı aynı ağda %2,5 hata veriyor — ve dört elemanı, iki yüz elli altı Q4 elemanından daha doğru. Ama bu üstünlüğün bir bedeli var ve bedel çok keskin.
1. Kilitlenmenin kaynağına dönelim
FEM-02'de kilitlenmenin sebebini şöyle koymuştuk: bir eleman, temsil edemediği bir şekil değiştirme alanına zorlandığında sahte kayma şekil değiştirmesi üretir ve gerçekte olduğundan çok daha rijit davranır.
Eğilmede şekil değiştirme kesit yüksekliği boyunca doğrusal değişir. Q4 elemanının yerdeğiştirme alanı iki doğrusaldır; türevi, yani şekil değiştirme, eleman içinde neredeyse sabit kalır. Doğrusal bir dağılımı üretemez.
2. Q8: sekiz düğüm, doğrusal şekil değiştirme
Q8 elemanı köşelerde dört, kenar ortalarında dört düğüm taşır. Şekil fonksiyonları ikinci dereceden (serendipity ailesi), dolayısıyla:
İntegrasyon için tam kural 3×3 Gauss, indirgenmiş kural 2×2'dir. Q8'de indirgenmiş integrasyon yaygın bir tercihtir ve tek bir sahte mod üretir — Q4'ün iki modu kadar yıkıcı değildir.
3. Kenar ortası düğümü paylaşılır
4. Yama sınaması
Yeni bir eleman formülasyonunun ilk sınavı, sabit şekil değiştirme alanını tam olarak üretip üretemediğidir.
| Eleman | İntegrasyon | Yama sınaması hatası | Sonuç |
|---|---|---|---|
| hesaplanıyor… | |||
5. Kilitlenme karşılaştırması
FEM-02'deki aynı konsol: L = 1000, H = 100, t = 10 mm, uç yükü 1000 N. Referans —.
| Ağ | Q4 tam | Q8 tam | Q8 indirgenmiş |
|---|---|---|---|
| hesaplanıyor… | |||
6. Serbestlik derecesi başına verim
"Dört eleman iki yüz elli altıyı yendi" ifadesi tek başına adil değil, çünkü Q8 elemanı daha çok düğüm taşır. Doğru karşılaştırma serbestlik derecesi üzerinden yapılır.
| Model | Eleman | Serbestlik derecesi | Hata |
|---|---|---|---|
| hesaplanıyor… | |||
Dörtgen sayısı üzerinden yapılan karşılaştırmalar yanıltıcıdır. Q8 daha çok düğüm istiyor — ama aynı doğruluğa çok daha az serbestlik derecesiyle ulaşıyor, ve çözüm maliyeti serbestlik derecesiyle belirlenir.
7. Hata neden sıfıra gitmiyor?
Tabloya dikkatlice bakınca bir şey göze çarpıyor: Q8'in hatası ağ sıklaştıkça düşüyor ama sıfıra gitmiyor, yaklaşık %0,2 dolayında duruyor.
Bunun pratik sonucu şu: yakınsama incelemesi, referansınızdan daha doğru olamaz. Q4 ile çalışırken bu sorun görünmezdi, çünkü elemanın kendi hatası referansın hatasını fazlasıyla gölgeliyordu. Yüksek dereceli elemana geçtiğinizde referansı da sorgulamanız gerekir.
8. Bedel: kenar ortası düğümü
Şimdiye kadar Q8 her ölçütte üstün çıktı. Bedeli burada.
Q8'in kenar ortası düğümü kenarın tam ortasında olmak zorunda değil — izoparametrik dönüşüm onu başka yere de taşıyabilir. Ama taşındığında dönüşüm hızla bozulur.
| Kenar ortası kayması | Uç sehimi | Hata | Düzgün ağa göre |
|---|---|---|---|
| hesaplanıyor… | |||
9. İki bozulmayı karşılaştırmak
FEM-04'te Q4'ün köşe düğümleri kaydırıldığında hatanın %13,95'ten %18,76'ya çıktığını ölçmüştük — bağıl olarak yaklaşık %35 kötüleşme. Burada Q8'in kenar ortası düğümü kaydırıldığında bağıl kötüleşme kırk kattan fazla.
Bunun sebebi anlaşılır. Q8'in üstünlüğü ikinci dereceden yerdeğiştirme alanını doğru temsil etmesinden geliyor. Kenar ortası düğümü kaydığında dönüşüm o ikinci dereceden alanı çarpıtıyor — yani elemanı üstün kılan şeyin kendisi bozuluyor.
10. Pratikte ne yapmalı
- Eğilmenin baskın olduğu yerde ikinci dereceden eleman kullanın. Q4 ile kaba ağda çalışmak, FEM-02'de görüldüğü gibi %70 hata üretebilir.
- Karşılaştırmayı serbestlik derecesi üzerinden yapın. Eleman sayısı yanıltıcıdır; maliyet serbestlik derecesiyle belirlenir.
- Kenar ortası düğümlerini ortada tutun. Eğrisel sınırlarda kaydırma gerekebilir, ama kaymayı gerektiği kadarla sınırlayın.
- Kenar ortası düğümünün paylaşıldığını doğrulayın. Paylaşılmazsa elemanlar kopar ve bu, sonucu bakınca hemen belli olmaz.
- Referansınızı sorgulayın. Yüksek dereceli elemanla çalışırken kalan hata elemandan değil referanstan geliyor olabilir.
- Ağ kalitesi gereksinimini gevşetmeyin. Yüksek derece, kötü ağa toleransı artırmaz; azaltır.
11. Kendini sınama
- Q8 neden eğilmede kilitlenmiyor? Cevabınızda "şekil değiştirmenin derecesi" nerede geçiyor?
- Dört Q8 elemanı 256 Q4 elemanından daha doğru. Bu karşılaştırma neden tek başına adil değil, ve adil olanı nedir?
- Q8'in hatası neden sıfıra gitmiyor? Bu elemanın mı yoksa referansın mı kusuru?
- Kenar ortası düğümü kaydığında Q8'i üstün kılan hangi özellik bozuluyor?
- Bir kenarın orta düğümü iki elemanca paylaşılmazsa ne olur? Sonuca bakınca fark eder misiniz?
- Q8'de indirgenmiş integrasyon kaç sahte mod üretir? Q4'e göre neden daha az tehlikeli?
Kaynaklar
- Serendipity ailesi ikinci dereceden şekil fonksiyonları ve izoparametrik formülasyon — standart sonlu elemanlar sonuçları; ülke yönetmeliğinden bağımsızdır.
- Yama sınaması — eleman geçerliliğinin standart ölçütü; bkz. FEM-02.
- Referans çözüm: kesme deformasyonu dâhil Timoshenko kirişi. Bu yazının 7. bölümü tam olarak bu referansın sınırını konu ediyor.
- Ders modülleri ve MATLAB uygulamaları için: Kwon, Y.W. & Bang, H., The Finite Element Method using MATLAB. Anlatım, türetme ve bütün sayılar kendi motorumuzdandır.
Yazıdaki bütün sayılar engine/fem2d.js tarafından üretilir ve tests/fem2d.test.js içinde ayrıca sabitlenmiştir (61/61). Q8 elemanı yama sınamasını hem tam hem indirgenmiş integrasyonda makine hassasiyetinde geçer; kenar ortası düğümlerinin komşu elemanlarca paylaşıldığı ayrı bir vakayla denetlenir.
Yüksek derece, kötü ağa toleransı artırmaz; azaltır.