Ana Sayfa / Yazılar / Sonlu Elemanlar FEM-06
Aynı İşlem, Zıt Anlam
FEM-05'te kenar ortası düğümünü kaydırmak hatayı kırk kat büyütüyordu. Eğrisel sınırda aynı işlem gereklidir. Ve bu yazıyı yazarken daha rahatsız edici bir şey çıktı: FEM-03'ün referansı olan 3,000, bu geometri için doğru değer değil.
1. Düz kenarda zarar, eğri kenarda düzeltme
FEM-05'in sonucu netti: Q8 elemanının kenar ortası düğümü kenarın ortasından kaydırıldığında hata %0,67'den %30,99'a çıkıyordu — kırk altı kat.
Şimdi delikli levhaya bakalım. Delik çeperi bir daire. Köşe düğümleri daire üzerinde, ama kenar ortası düğümü varsayılan olarak iki köşenin aritmetik ortasına konuyor — yani kirişin ortasına, dairenin değil.
İzdüşüm yapılmazsa delik, çokgenle temsil edilir. Ağdaki sapmayı ölçtük: —.
2. İzdüşüm tutarsız olamaz
Burada gizli bir tuzak var. Bir kenar iki elemanca paylaşılır. İzdüşüm her iki elemanda ayrı ayrı yapılırsa, ikisi farklı sonuç verirse elemanlar birbirinden kopar.
Motorumuz bunu yapısal olarak imkânsız kılıyor: kenar ortası düğümü bir kez üretilir, kenar anahtarıyla saklanır ve iki elemana da aynı düğüm verilir. İzdüşüm o tek üretim anında uygulanır. Tutarsızlık için yer yok — test paketinde ayrı bir vaka olarak korunuyor.
3. Üç yolun karşılaştırması
Aynı problem üç ayrı yolla çözüldü: Q4, izdüşümsüz Q8, izdüşürülmüş Q8. Delik yarıçapı a = 10, çeyrek model yarı genişliği W = 200, kademelendirme 3.
| Ağ | Q4 | Q8 düz kenar | Q8 izdüşürülmüş | Q8 DOF |
|---|---|---|---|---|
| hesaplanıyor… | ||||
4. İzdüşümün tek başına katkısı
Üçüncü ve dördüncü sütunlar aynı ağı, aynı eleman sayısını, aynı serbestlik derecesini kullanıyor. Aralarındaki tek fark dört düğümün konumu.
FEM-05'in aynadaki görüntüsü: orada dört düğümü yerinden oynatmak hatayı kırk kat büyütüyordu, burada yerine oturtmak birkaç kat küçültüyor.
5. Bir şey tutmuyordu
Tabloya bakınca şu göze çarpıyor: izdüşürülmüş Q8 ağ sıklaştıkça 3,000'e değil, 3,02 dolayına yakınsıyor. Ve orada duruyor.
FEM-05'te benzer bir durumla karşılaşmış ve sebebin referans olduğunu bulmuştuk. Aynı şüpheyle başladım: levhayı genişletirsem 3,000'e iner mi?
6. Kendi ölçüm hatam
Sebep şuydu: kademelendirme oransaldır. Levha genişledikçe delik çevresindeki ilk eleman halkası fiziksel olarak kalınlaşıyor, yani delik yakınındaki çözünürlük düşüyor. Sonlu genişlik etkisiyle ağ seyrelmesini aynı tabloda karıştırmışım.
Düzeltme basitti: her a/W oranında ağı ayrı ayrı yakınsatmak. Bu yapıldığında sonuç tekdüze ve fiziksel çıkıyor.
7. Referansın kendisi yanlıştı
| a/W | SCF | 3,000'den sapma | Heywood kestirimi |
|---|---|---|---|
| hesaplanıyor… | |||
Kirsch çözümü sonsuz levha varsayar. Bizim modelimiz sonlu: a/W = 0,05. Bu geometride doğru yığılma katsayısı 3,000 değil 3,021. Ve bu iki bağımsız yolla doğrulanıyor — a/W küçüldükçe 3,000'e yaklaşması, ve sonlu genişlik için bilinen bir yaklaşık bağıntıyla örtüşmesi.
8. FEM-03'e düzeltme
Yönü önemli: FEM-03 Q4'ün hatasını olduğundan küçük göstermiş. Yani yazının asıl tezi — gerilme yavaş yakınsar, kaba ağ tepe gerilmeyi düşük okur — zayıflamıyor, güçleniyor. Yanlış olan bildirilen sayının kendisiydi.
Bu, sitenin kendini altıncı kez düzeltmesi. Ve deseni yine aynı: hata, metni tekrar okurken değil yeni bir şey inşa ederken çıktı. Q8 olmasaydı bulunamazdı, çünkü Q4'ün kendi hatası referansın hatasını fazlasıyla gölgeliyordu.
9. Motor geçersiz ağı reddediyor
Ölçümler sırasında bir ağ reddedildi. Kademelendirme artırıldığında, radyal kenarlarda kenar ortası düğümü kenarın ortasında kalıyor ama eleman kuvvetle kademelendiği için orta düğüm elemanın "doğal" ortasından uzağa düşüyor — ve Jacobian determinantı negatife geçiyor.
Klasik kural burada kendini gösteriyor: Q8'in kenar ortası düğümü kenarın orta yarısında kalmalıdır. Ağır kademelendirme bunu radyal kenarlarda ihlal eder.
10. Pratikte ne yapmalı
- Eğrisel sınırda kenar ortası düğümünü gerçek eğriye oturtun. Düz bırakmak daireyi çokgenle temsil eder.
- Düz kenarda oturtmayın. Aynı işlem orada zarar verir — bkz. FEM-05.
- İzdüşümün paylaşıldığını doğrulayın. İki elemanda ayrı hesaplanırsa elemanlar kopar.
- Ağ parametrelerini değiştirirken tek değişkeni değiştirin. Bu yazıdaki ilk ölçümüm, genişlik ile ağ çözünürlüğünü aynı tabloda karıştırdığı için yanlıştı.
- Referansınızın varsayımlarını okuyun. Kirsch sonsuz levha varsayar; sizin levhanız sonlu.
- Kademelendirmeyi Q8'de ölçülü kullanın. Orta düğüm kenarın orta yarısında kalmalı.
11. Kendini sınama
- Kenar ortası düğümünü kaydırmak ne zaman hata, ne zaman düzeltmedir? Ayrımı ne belirler?
- İzdüşüm iki elemanda ayrı ayrı yapılsaydı ne olurdu? Sonuca bakınca fark eder miydiniz?
- Levha genişledikçe SCF neden 3,000'e yaklaşır?
- İlk ölçümüm neden tekdüze çıkmadı? Hangi iki etkiyi karıştırmıştım?
- FEM-03'ün hatası düzeltilince tezi zayıfladı mı güçlendi mi? Neden?
- Ağır kademelendirme Q8'de neden geçersiz eleman üretir? Q4'te aynı sorun neden yok?
Kaynaklar
- Kirsch çözümü — sonsuz levhada dairesel delik çevresindeki gerilme alanı; klasik elastisite sonucu.
- Sonlu genişlik düzeltmesi için bilinen yaklaşık bağıntılar (Heywood tipi). Bu yazıda kestirim olarak, bağımsız bir kontrol amacıyla kullanılmıştır.
- İzoparametrik formülasyon ve serendipity ailesi ikinci dereceden elemanlar — standart sonlu elemanlar sonuçları.
- Ders modülleri 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 (81/81). Sayfadaki canlı tablolar tarayıcıda makul sürede çalışacak ağlarla sınırlıdır; daha ince ağlarla elde edilen yakınsak değer (3,0209) test paketinde sabitlenmiştir. 6. bölümdeki ölçüm hatası düzeltilmiş hâliyle değil, hatasıyla birlikte anlatılmıştır — çünkü hatanın kendisi öğreticidir.
Referansınızın varsayımlarını okumazsanız, elemanın hatasıyla referansın hatasını karıştırırsınız.