Geçerli Olmak Yetmez

1. İzoparametrik dönüşüm

Dörtgen bir elemanın rijitliğini doğrudan gerçek geometri üzerinde integre etmek zordur: her eleman farklı şekildedir. İzoparametrik yaklaşım bunu tersine çevirir — integrasyon her zaman birim kare üzerinde yapılır ve gerçek eleman, bu kareden bir dönüşümle elde edilir.

2. Jacobian ne taşır?

Dönüşümün türevi Jacobian matrisidir. Determinantı, ana elemandaki bir alan öğesinin gerçek elemanda ne kadar alana karşılık geldiğini verir:

Şekil fonksiyonlarının gerçek koordinatlara göre türevleri de Jacobian'ın tersi üzerinden hesaplanır. Dolayısıyla Jacobian, dönüşümün kalitesini hem integrasyona hem şekil değiştirme hesabına taşır.

3. detJ yalnızca paralelkenarda sabittir

Bu, sezgiye aykırı ama kesin bir sonuçtur ve doğrudan sınanabilir:

Bunun anlamı şudur: çarpıklık tek başına detJ'yi bozmaz. Bozan şey, karşılıklı kenarların paralelliğini yitirmesidir — yani yamuklaşma ve konikleşme. Bu ayrım ağ üretiminde önemlidir: kayık ama paralelkenar bir ağ, yamuk bir ağdan daha iyidir.

4. Kalite ölçütleri

Ticari programların "mesh quality" raporlarında gördüğünüz ölçütler bu geometriden çıkar:

• Jacobian oranı — min(detJ)/max(detJ). Dönüşümün eleman içinde ne kadar değiştiği. 1'e yakın olması istenir.
• En-boy oranı — kenar uzunluklarının oranı. Yüksek olması, bir yöndeki gradyanın diğerinden çok daha kötü temsil edildiği anlamına gelir.
• Çarpıklık — köşe açılarının 90°'den sapması.
• Konikleşme — karşılıklı kenar uzunluklarının oranı; yamuklaşmanın ölçüsü.

Bu ölçütler göstergedir, kesin bir hata sınırı vermezler. Asıl soru şu: ne kadar kötüleşirse ne kadar hata olur?

5. Yama sınaması bozulmaya dayanıyor mu?

FEM-02'de yama sınamasının bir elemanın geçerliliğinin standart ölçütü olduğunu görmüştük. Şimdi düzgün bir ağın iç düğümlerini kademeli olarak kaydırıp sınamayı tekrarlıyoruz. Sınır düğümleri sabit tutuluyor, böylece problem aynı kalıyor.

6. Ama doğruluk bozuluyor

Aynı bozulmayı, sabit olmayan bir alana — eğilen bir konsola — uygulayalım. FEM-02'deki referans: Timoshenko çözümü.

7. Gerekli ama yeterli değil

Bu iki tablo birlikte okunduğunda yazının asıl noktası çıkıyor:

Pratik karşılığı şudur: bir programın "elemanlarım yama sınamasını geçiyor" demesi, sizin ağınızın iyi olduğu anlamına gelmez. İkisi farklı şeyleri ölçer.

8. Yön yine emniyetsiz

Bu serideki bir örüntü burada da tekrarlanıyor. FEM-02'de yerdeğiştirmenin daima aşağıdan yakınsadığını, FEM-03'te tepe gerilmenin düşük çıktığını görmüştük. Bozulma da aynı yönde çalışıyor: eleman daha rijit hâle geliyor, sehim daha da küçük çıkıyor.

Yani bozuk ağın hatası, zaten var olan kilitlenmenin üstüne biniyor. İkisi birbirini tersine çevirmiyor, topluyor.

9. Ölçüt ne kadar anlamlı?

Jacobian oranı bir gösterge olduğuna göre, gerçekten hatayla ilişkili mi? Tablodaki iki sütunu yan yana koyalım:

Düzeltilmiş ölçümde gerilme ve sehim hataları birlikte büyüyor; ağır bozulmada gerilme hafifçe öne geçiyor. Jacobian oranı düştükçe her ikisi de kötüleşiyor — yani ölçüt anlamlı, ama tek başına bir hata sınırı vermiyor.

10. Pratikte ne yapmalı

1. Negatif Jacobian bir uyarı değil, bir durdurma sebebidir. Dönüşüm tersinir değilse sonuç anlamsızdır; "uyarıyı geç" demeyin.
2. Yama sınamasını geçmeyi ağ kalitesiyle karıştırmayın. İkisi farklı şeyleri ölçer; birincisi eleman formülasyonunu, ikincisi sizin ağınızı.
3. Kritik bölgede düzgün ağ kurun. Gradyanın büyük olduğu yerde bozuk eleman, hem FEM-03'teki yavaş gerilme yakınsamasıyla hem de buradaki bozulma hatasıyla çarpışır.
4. Paralelkenar ağı yamuk ağa tercih edin. Kayıklık detJ'yi bozmaz; yamuklaşma bozar.
5. Ölçütleri mutlak eşik gibi kullanmayın. "J oranı 0,7 üzerindeyse sorun yok" diye bir kural yoktur; ölçüt göstergedir, doğrulama yakınsama incelemesiyle yapılır.
6. Bozuk ağın hatasının yönünü bilin. Daha rijit, yani daha küçük sehim, yani emniyetsiz taraf.

11. Kendini sınama

1. İzoparametrik yaklaşımın "izo" kısmı neyi anlatır?
2. detJ hangi eleman şekillerinde sabit kalır? Çarpık bir eleman bu kümeye girebilir mi?
3. Bir eleman yama sınamasını geçiyorsa, o ağdaki sonucun doğru olduğunu söyleyebilir misiniz?
4. Bozulma sehimi hangi yönde değiştirir? Bu, kilitlenmeyle nasıl etkileşir?
5. Negatif Jacobian ne anlama gelir ve neden hesabı durdurmak gerekir?
6. Jacobian oranı 0,6 olan bir ağda hatanın ne olduğunu söyleyebilir misiniz? Neden?
7. Bu yazıdaki gerilme ölçümü neden ilk kurulduğunda yanlıştı?