Ana Sayfa / Yazılar / Sonlu Elemanlar FEM-01
Ağı Sıklaştırmak Her Zaman Doğru Cevap Değildir
Sonlu elemanlarda bir sonuç yanlışsa refleks şudur: daha çok eleman. Bu yazıda, tek elemanın tam sonucu verdiği ve otuz iki elemanın hatalı kaldığı bir durum var. Fark ağda değil, yayılı yükün düğümlere nasıl aktarıldığında.
1. Sonlu elemanlar neyi çözer?
Doğrudan rijitlik yönteminde yapı, uçlarından birbirine bağlanmış elemanlara ayrılır. Her eleman için, uç yerdeğiştirmeleriyle uç kuvvetlerini ilişkilendiren bir rijitlik matrisi yazılır; bunlar ortak düğümlerde toplanarak sistem matrisi kurulur ve
ETABS, SAP2000 ve ABAQUS gibi programlar arka planda bunu yapar. "Kara kutu" hissi çoğunlukla bu taşımanın görünmez olmasından gelir.
2. Kiriş elemanının gizli üstünlüğü
2 boyutlu çerçeve elemanının eğilme kısmı, kübik Hermite şekil fonksiyonlarına dayanır. Bu keyfi bir seçim değildir: eleman içinde yayılı yük yokken Euler–Bernoulli kirişinin denklemi
Bu, sürekli ortam elemanlarından (düzlem gerilme, katı) belirgin bir farktır. Orada şekil fonksiyonları gerçek çözümü üretmez ve ağ sıklaştırma gerçekten gereklidir. Kiriş elemanı ayrıcalıklıdır — ve bu ayrıcalık bir sonraki bölümde kaybedilecek.
3. Sorun: yayılı yük düğümde durmaz
Bir kirişe düzgün yayılı yük etkidiğinde, eleman içinde artık d⁴v/dx⁴ = 0 değildir:
Ama tam olarak bitmiyor. Çünkü {F} vektörü doğru kurulursa, kübik yaklaşım düğümlerde hâlâ tam sonucu verebilir. Mesele, w'yi düğümlere nasıl taşıdığımızdır.
4. İki yol: yalın ve tutarlı
Yalın (lumped) yaklaşım
Sezgisel olan yol: yükü düğümlere pay et. Uzunluğu ℓ olan bir eleman için toplam yük wℓ'dir; yarısı bir uca, yarısı diğerine verilir.
Tutarlı (consistent) yük vektörü
Doğru yol, eşdeğerliği kuvvet dengesinden değil virtüel işten kurmaktır: düğüm yükleri, gerçek yayılı yükle aynı virtüel işi yapacak şekilde seçilir. Şekil fonksiyonları üzerinden integral alındığında ortaya bir moment terimi de çıkar:
5. Sayısal karşılaştırma
Ankastre konsol, düzgün yayılı yük. Kapalı-form çözüm δ = wL⁴/(8EI) = —. Aşağıdaki tablo sayfa açılırken çözücüden üretilir.
| Eleman | Yalın yük | Hata | Tutarlı yük | Hata |
|---|---|---|---|---|
| hesaplanıyor… | ||||
6. Neden tutarlı yük tam sonuç veriyor?
Tutarlı yük vektörü, eleman içindeki yayılı yükün düğümlerdeki etkisini tam olarak temsil eder. Bir kirişte ankastrelik momenti wℓ²/12'dir; tutarlı vektör bunu aynen düğüme aktarır. Sonuçta sistem denklemi, gerçek yükleme altındaki düğüm yerdeğiştirmelerini veren denklemle özdeş hâle gelir.
Yalın yaklaşım bu moment çiftini atladığı için, konsolda elemanın uçlarını gerçekte olduğundan daha serbest bırakır. Serbest kalan uç daha çok döner, daha çok sehim üretir — tabloda yalın sütununun büyük çıkmasının sebebi budur.
7. Yakınsama mertebesi: O(h²)
Yalın yaklaşım yanlış olsa da tutarsız değildir: ağ sıklaştıkça doğru cevaba yaklaşır. Tabloda hata sütununa dikkat edin — her incelemede hata tam olarak dörde bölünüyor.
Buradan pratik bir araç çıkar: bir modelin yakınsayıp yakınsamadığını bilmek için tam çözümü bilmeniz gerekmez. Ardışık iki ağdaki farkın nasıl küçüldüğüne bakmak yeter.
8. Sınırı da söyleyelim: eleman İÇİ tam değildir
Bunun doğrudan sonucu şudur: eleman içindeki maksimum momenti veya sehimi merak ediyorsanız, ağ sıklaştırmak yine gerekir. Tek elemanlı bir konsolun uç sehimi tamdır, ama açıklık ortasındaki sehim değildir. Aynı şekilde iç kuvvet diyagramları da eleman sayısı arttıkça düzelir.
Yani "ağ sıklaştırmak gereksiz" demiyoruz. Söylediğimiz şey, ağ sıklaştırmanın yanlış yük vektörünü düzeltmediği; onu ancak yavaşça örttüğüdür.
9. Motorun kendi denetimi
Bu yazının sayıları sitedeki engine/frame-solver.js çözücüsünden gelir. Yazıyı hazırlarken çözücü altı bağımsız kapalı-form çözüme karşı tarandı:
| Kontrol | Kapalı form | Çözücü | Durum |
|---|---|---|---|
| hesaplanıyor… | |||
10. Pratikte ne yapmalı
- Yayılı yükü programın kendisine tanımlayın, elle düğümlere bölmeyin. Ticari programlar tutarlı yük vektörünü kendileri kurar; eliyle "yükü düğümlere paylaştıran" kullanıcı bu üstünlüğü kaybeder.
- Bir sonuç şüpheliyse önce yük tanımına bakın, ağa değil. Ağ sıklaştırmak pahalıdır ve yanlış yük vektörünü düzeltmez.
- Yakınsamayı hata oranıyla sınayın. İki ağ çözün; farkın dörde bölünmesi beklenen davranıştır. Bölünmüyorsa modelde başka bir sorun vardır.
- Eleman içi büyüklükler için ağ sıklaştırın. Düğüm sonuçlarının tam olması, diyagramların tam olduğu anlamına gelmez.
- Basit bir kapalı-form vakayla programı sınayın. Yeni bir yazılıma geçtiğinizde ilk yapılacak iş, δ = PL³/3EI gibi bilinen bir çözümü koşturmaktır.
11. Kendini sınama
- Ucundan tekil yükle yüklenmiş bir konsolu ikiye bölerseniz sonuç değişir mi? Neden?
- Yayılı yükte tutarlı yük vektörünün moment terimi neden −wℓ²/12 ve +wℓ²/12 olarak zıt işaretlidir?
- Bir modelde ağı iki katına çıkardınız ve sonuç %8 değişti; tekrar iki katına çıkardınız, %2 değişti. Bu davranış beklenen midir? Nihai değeri nasıl kestirirsiniz?
- Tek elemanlı konsolun uç sehimi tam çıkıyorsa, açıklık ortasındaki sehim neden tam değildir?
- Düzlem gerilme elemanlarında da tek eleman tam sonuç verir mi? Kiriş elemanından farkı nedir?
Kaynaklar
- Doğrudan rijitlik yöntemi ve Hermite şekil fonksiyonları — standart sonlu elemanlar sonuçları; ülke yönetmeliğinden bağımsızdır.
- Tutarlı yük vektörünün virtüel iş üzerinden türetilmesi; ankastrelik momenti wℓ²/12 ile özdeşliği.
- Ders modülleri ve MATLAB uygulamaları için: Kwon, Y.W. & Bang, H., The Finite Element Method using MATLAB. Bu yazıdaki anlatım ve türetme kendi sözlerimizledir; kaynak, konu sıralaması ve kapsam için atıf olarak verilmiştir.
Yazıdaki bütün sayılar engine/frame-solver.js tarafından üretilir ve tests/frame-solver.test.js içinde ayrıca sabitlenmiştir (17/17). Doğrulama, üçüncü taraf bir yazılımın çıktısıyla karşılaştırma yoluyla değil, kapalı-form elastisite çözümleriyle yapılmıştır.
Ağ sıklaştırmak yanlış yük vektörünü düzeltmez; yalnızca yavaşça örter.