Ana Sayfa / ETABS & SAP2000 / Sonlu Eleman Temelleri
Sonlu Eleman Temelleri — Eleman Tipi Seçimi ve Mesh Yakınsaması
Frame/Shell/Solid eleman ayrımı, Membrane vs Shell vs Thick/Thin Shell seçimi, mesh yoğunluğunun sonuç doğruluğuna etkisi — CST elemanının kayma kilitlenmesi (shear locking) davranışıyla somutlaştırılmış.
transpose() yardımcı fonksiyonu eklendi; regresyon testi üç bağımsız kontrolle (boyut, simetriklik, rijit gövde hareketi altında sıfır kuvvet — temel bir FE sağlık kontrolü) kilitlendi. Bu, Modal Analiz sayfasındaki eigenSolve hatasıyla aynı dersi tekrarlıyor: gerçekçi bir örnek denenmeden hata gizli kalabiliyor.
İçindekiler
1. Teori — Eleman Tipleri ve Yakınsama 2. ETABS'ta Eleman Tipi ve Mesh Ayarları 3. Örnek 1 (Kolay) — Kiriş Elemanı, Tek Elemanla Tam Çözüm 4. Örnek 2 (Orta) — CST Elemanı, Rijitlik Matrisi Montajı 5. Örnek 3 (Zor) — Rijit Gövde Hareketi Testi (FE Sağlık Kontrolü) 6. Örnek 4 (Çok Zor) — CST'nin Kayma Kilitlenmesi (Shear Locking)1. Teori — Eleman Tipleri ve Yakınsama
ETABS/SAP2000 üç temel eleman ailesi sunar: Frame (çubuk/kiriş — kolon, kiriş, çapraz), Shell/Area (kabuk/alan — döşeme, perde), Solid (katı — kalın temel, özel detaylar). Doğru eleman seçimi, hem sonuç doğruluğu hem hesaplama süresi açısından kritiktir.
Membrane vs Shell vs Thick/Thin Shell
Shell: her ikisi birden (in-plane + eğilme) Thick Shell: enine kayma deformasyonu dahil (Mindlin) Thin Shell: ihmal edilir (Kirchhoff) CSI resmi dokümantasyonu; sheerforceeng.com karşılaştırmalı analiz
Kritik pratik sonuç: Bir döşemeyi membrane olarak modellerseniz (yalnızca düzlem-içi rijitlik), ETABS bu döşemeyi otomatik olarak ağa örmez ("cookie-cut" mesh) ve yük dağılımı gerçekçi olmayan sonuçlar üretebilir — bir mühendisin belgelediği örnekte, membrane döşeme kolonlara tamamen yanlış bir yük dağılımı vermiştir. Kalın (500mm+) döşeme/perdelerde Thick Shell tercih edilir; ince döşemelerde (150-200mm) Thin Shell genelde yeterlidir.
Mesh Yakınsaması — Genel İlke
Eleman boyutu küçüldükçe (mesh sıklaştıkça) sonuç, gerçek çözüme yakınsar — ama her eleman tipi her problem için eşit hızda yakınsamaz. Bu, Örnek 4'te somut olarak gösterilecektir.
2. ETABS'ta Eleman Tipi ve Mesh Ayarları
Duvar/döşeme kesiti tanımlanırken tip seçimi (Membrane / Shell-Thin / Shell-Thick / Layered) burada yapılır.
"Auto Mesh Objects into Structural Elements" — üç temel seçenek: kirişlerde/duvar kenarlarında ağ örme, görünür grid çizgilerinde ağ örme, maksimum eleman boyutuyla ek bölme. Eğrisel geometrilerde (kavisli teras vb.) dördüncü seçenek ("further subdivide... maximum element size") ile ağ sıklaştırılır.
Düz, açıklıksız duvarlar varsayılan olarak ağ örülmez (tek eleman olarak kalır) — kavisli veya açıklıklı duvarlarda otomatik dikdörtgen ağ devreye girer.
ETABS mesh yakınsamasını otomatik kontrol etmez — mühendis, kritik bir bölgede (örn. büyük açıklıklı, keskin köşeli bir döşeme) mesh boyutunu kademeli küçültüp sonucun (moment, gerilme) belirli bir toleransın (örn. %1-2) altında değiştiğini görene kadar tekrarlamalıdır.
3. Örnek 1 — Kiriş Elemanı, Tek Elemanla Tam Çözüm
Sistem: Konsol kiriş, L=3m, E=200 GPa, I=8.36×10⁻⁶ m⁴ (IPE300 benzeri), uçta P=10 000 N.
| Eleman sayısı | δ (FEM, m) | Fark |
|---|---|---|
| 1 (tek Euler-Bernoulli elemanı) | 0.053828 | %0.0000 |
4. Örnek 2 — CST Elemanı, Rijitlik Matrisi Montajı
Sistem: 2m×1m dikdörtgen bir plaka, köşegen boyunca iki CST (Constant Strain Triangle) üçgenine bölünmüş. E=200 GPa, ν=0.3, t=10mm.
| Eleman | Düğümler | Alan (m²) | Ke[0][0] (N/m) |
|---|---|---|---|
| 1 | (0,0)-(2,0)-(2,1) | 1.0 | 5.495×10⁸ |
| 2 | (0,0)-(2,1)-(0,1) | 1.0 | 7.692×10⁸ |
Toplam alan kontrolü: 1.0+1.0 = 2.0 m² = 2m×1m ✓ — iki üçgenin toplamı orijinal dikdörtgeni tam olarak kaplıyor, montajın geometrik tutarlılığının ilk kontrolü.
5. Örnek 3 — Rijit Gövde Hareketi Testi (FE Sağlık Kontrolü)
Amaç: Her sonlu eleman rijitlik matrisinin sağlaması gereken temel bir fiziksel kural: eğer eleman, deformasyon olmadan (rijit cisim gibi) tamamen bir yöne kayarsa, iç kuvvet sıfır olmalıdır (eleman "gerinmiyor").
Sistem: Örnek 2'deki Eleman 1'in (Ke, 6×6) düzeltilmiş rijitlik matrisi.
| Test | Sonuç kuvvet vektörü (N) | Durum |
|---|---|---|
| Rijit x-ötelemesi | [0, 0, ~10⁻⁷, ~10⁻⁷, 0, 0] | ✓ Makine hassasiyetinde sıfır |
| Rijit y-ötelemesi | [0, 0, ~10⁻⁸, 0, 0, 0] | ✓ Makine hassasiyetinde sıfır |
6. Örnek 4 — CST'nin Kayma Kilitlenmesi (Shear Locking)
Amaç: Örnek 1'de kiriş elemanının tek elemanla tam sonuç verdiğini gördük. CST elemanı için durum tam tersidir — eğilme baskın problemlerde kaba mesh kullanmak ciddi hatalara yol açar.
Fiziksel Mekanizma
Saf eğilmedeki bir kesitte gerçek yerdeğiştirme alanı doğrusal olmayan bir fonksiyondur (Euler-Bernoulli teorisine göre kesitler düzlem kalır ama dönerek eğrilir). CST elemanının şekil fonksiyonları ise doğrusaldır ve gerinim alanı sabittir (adının kaynağı budur: Constant Strain Triangle). Bir CST üçgeni saf eğilmeye zorlandığında, gerçekte var olmayan yapay bir kayma gerinimi üretir — bu, elemanı fiziksel olarak olması gerekenden çok daha rijit gösterir.
| Yükleme tipi | CST davranışı | Kaynak |
|---|---|---|
| Saf eksenel (çekme/basınç) | Tam doğru — gerinim gerçekten sabit | Teorik olarak kesin |
| Saf eğilme, kaba mesh (1-2 eleman/derinlik) | Aşırı rijit — gerçek deplasmanın önemli ölçüde altında sonuç | ScienceDirect (Cook) ders kitabı özeti |
| Saf eğilme, ince mesh (çok sayıda eleman/derinlik) | Gerçek çözüme yakınsar (O(h²) oranıyla) | Grokipedia CST yakınsama özeti |