Sonlu Eleman Yöntemine Giriş: Mühendis İçin

Sonlu eleman yöntemi (SEY / FEM — Finite Element Method), karmaşık geometri ve yük koşullarındaki yapıların sayısal analizi için kullanılan en kapsamlı yaklaşımdır. Yapı sonsuz serbestlik dereceli sürekli bir ortam olarak değil, sonlu sayıda elemana bölünmüş ayrık bir sistem olarak çözülür. Her eleman basit şekil fonksiyonlarıyla (shape functions) tanımlanır; elemanların bir araya getirilmesiyle global denklem sistemi oluşturulur.

TBDY 2018, perde duvarların modellemesini çubuk elemandan kabuk (shell) sonlu elemana zorunlu olarak geçirmiş; etkin kesit rijitliği çarpanlarını Tablo 4.2 ile standartlaştırmış ve ağ bağımsızlık kontrolünü mühendis sorumluluğu kapsamına almıştır. Türkiye'de SAP2000 ve ETABS en yaygın kullanılan SEY yazılımlarıdır; her ikisi de TBDY 2018 uyumlu analiz modülleri sunmaktadır.

1. Yapısal Analiz Yöntemlerinin Karşılaştırması

Tablo 1: Yapısal Analiz Yöntemlerinin Karşılaştırması

Dikkat: SEY sonuçlarının güvenilirliği, eleman tipi seçimine ve sınır koşullarının doğruluğuna kritik biçimde bağlıdır. Yanlış eleman tipi seçimi %30–60 oranında hatalı iç kuvvet üretebilir.

2. Temel Kavramlar

Tablo 2: Temel Kavramlar

3. SEY Formülasyonu

3.1 Temel Denklem

Yapısal SEY, sanal iş ilkesine (virtual work principle) dayanır:

$$\int_V [\delta \varepsilon]^T [\sigma] \, dV = \int_V [\delta u]^T \{b\} \, dV + \int_S [\delta u]^T \{t\} \, dS$$

Matris formunda:

$$[K] \{u\} = \{F\}$$

$[K]$ global rijitlik matrisi, $\{u\}$ düğüm deplasman vektörü, $\{F\}$ dış yük vektörüdür.

3.2 Şekil Fonksiyonları

1B lineer eleman (çubuk/kiriş) şekil fonksiyonları:

$$N_1(\xi) = \frac{1 - \xi}{2}, \quad N_2(\xi) = \frac{1 + \xi}{2}, \quad \xi \in [-1, +1]$$

2B dörtgen eleman (Q4) şekil fonksiyonları — izoparametrik formülasyon:

$$N_i(\xi, \eta) = \frac{1}{4}(1 + \xi_i \xi)(1 + \eta_i \eta), \quad i = 1, 2, 3, 4$$

İzoparametrik yaklaşımda geometri ve deplasman aynı şekil fonksiyonları ile tanımlanır: $x = \sum N_i x_i$, $y = \sum N_i y_i$.

3.3 Birim Şekil Değiştirme–Deplasman İlişkisi

$$\{\varepsilon\} = [B] \{u_e\}$$

$[B]$ matrisi şekil fonksiyonlarının türevlerinden oluşur:

$$[B] = \begin{bmatrix} \dfrac{\partial N_1}{\partial x} & 0 & \dfrac{\partial N_2}{\partial x} & 0 \\[8pt] 0 & \dfrac{\partial N_1}{\partial y} & 0 & \dfrac{\partial N_2}{\partial y} \\[8pt] \dfrac{\partial N_1}{\partial y} & \dfrac{\partial N_1}{\partial x} & \dfrac{\partial N_2}{\partial y} & \dfrac{\partial N_2}{\partial x} \end{bmatrix}$$

3.4 Eleman Rijitlik Matrisi

$$[K_e] = \int_V [B]^T [D] [B] \, dV$$

Nümerik entegrasyon (Gauss-Legendre kuralı):

$$[K_e] = \sum_{i} \sum_{j} w_i w_j [B(\xi_i, \eta_j)]^T [D] [B(\xi_i, \eta_j)] |J(\xi_i, \eta_j)|$$

Tablo 3: Eleman Rijitlik Matrisi

Saha Notu: Q4 eleman için 2×2 Gauss entegrasyon (4 nokta) standarttır. Azaltılmış entegrasyon (1×1) hesabı hızlandırsa da "shear locking" ve hourglass kararsızlıklarına yol açabilir; yapısal mühendislik uygulamalarında önerilmez.

3.5 Malzeme Matrisi [D]

Düzlem gerilme (plane stress) — ince plak, döşeme:

$$[D] = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \dfrac{1-\nu}{2} \end{bmatrix}$$

Düzlem şekil değiştirme (plane strain) — zemin, baraj, uzun perdeler:

$$[D] = \frac{E(1-\nu)}{(1+\nu)(1-2\nu)} \begin{bmatrix} 1 & \dfrac{\nu}{1-\nu} & 0 \\[6pt] \dfrac{\nu}{1-\nu} & 1 & 0 \\[6pt] 0 & 0 & \dfrac{1-2\nu}{2(1-\nu)} \end{bmatrix}$$

Türkiye betonarme standart değerleri:

• $E_c = 3250\sqrt{f_{ck}} + 14\,000$ (MPa) — TS 500:2000 Denklem 3.1
• C30 beton: $E_c \approx 32\,000$ MPa
• Poisson oranı: $\nu = 0.2$ (TS 500:2000 §3.1)
• C50 ve üzeri: TS EN 1992-1-1:2004 §3.1.3 (TBDY 2018 §5.4.1)

4. Eleman Türleri

Tablo 4: Eleman Türleri

5. SEY Analiz Adımları

Adım 1 — Geometri ve Ağ Oluşturma

Yapı geometrisi tanımlanır. Eleman boyutu ve tipi seçilir (TBDY 2018 §4.5 model kuralları kapsamında). Düğüm koordinatları ve eleman bağlantı tablosu (connectivity) oluşturulur.

Adım 2 — Malzeme ve Kesit Özellikleri

• Beton: $E_c$, $\nu$, $\rho$ — TS 500:2000 §3.1
• Beton sınıfı C25–C80; C50 üzeri için TS EN 1992-1-1:2004 (TBDY 2018 §5.4.1)
• Donatı çeliği: B420C veya B500C nervürlü (TS 708:2016); çekme/akma oranı ≥ 1.35 (TBDY 2018 §5.3.1)
• Doğrusal olmayan analiz için akma gerilmesi, pekleşme katsayısı, kırılma enerjisi

Adım 3 — Sınır Koşulları ve Yükler

• Mesnet koşulları: ankastre, mafsallı, yay mesnet
• Yükler: TS 498:1997 (karakteristik yükler), deprem TBDY 2018 §4.4
• Deprem spektrumu: TBDY 2018 §2.3 ve TDTH 2018 haritası

Adım 4 — Çözüm

Lineer sistem çözüm yöntemleri:

• Doğrudan çözüm: Gauss eliminasyon, Cholesky ($[K]$ simetrik pozitif tanımlı için)
• Yinelemeli çözüm: Conjugate Gradient (büyük sistem, seyrek matris)
• Modal analiz: $\det([K] - \omega^2 [M]) = 0$ — TBDY 2018 §4.8

Adım 5 — Post-İşlem

Deplasman alanı $\{u\}$ çözüldükten sonra:

$$\{\varepsilon\} = [B]\{u_e\}, \quad \{\sigma\} = [D][B]\{u_e\}$$

Von Mises gerilmesi:

$$\sigma_{VM} = \sqrt{\sigma_x^2 + \sigma_y^2 - \sigma_x \sigma_y + 3\tau_{xy}^2}$$

6. TBDY 2018 Kapsamında SEY Modelleme Gereksinimleri

6.1 Genel Modelleme Kuralları (TBDY 2018 §4.5)

TBDY 2018, DBYBHY 2007 (mülga, yerine TBDY 2018)'ye kıyasla önemli yenilikler getirmiştir:

• Perde duvarlar: Boşluksuz tüm betonarme perdeler kabuk (shell) sonlu elemanlarla modellenmelidir (TBDY 2018 §4.5.3.4). Gövde kayma çerçevesi (bant model) hiçbir koşulda kullanılamaz.
• Döşemeler: Tüm plak döşemeler 2B kabuk sonlu elemanlarla modellenir (TS 500 §11.3.2, TBDY 2018 §4.5.5).
• Kolon ve kirişler: Çubuk (frame) eleman olarak tanımlanmaya devam edilir (TBDY 2018 §4.5.2).
• Poligon perdeler: Çekirdek perdeler dahil her geometri kabuk elemanlarla modellenmelidir (TBDY 2018 §4.5.3.5).

6.2 Etkin Kesit Rijitliği Çarpanları (TBDY 2018 Tablo 4.2)

Tablo 5: Etkin Kesit Rijitliği Çarpanları (TBDY 2018 Tablo 4.2)

Kritik Not: TBDY 2018 Madde 4.5.8.3 uyarınca etkin kesit rijitliği çarpanları yalnızca deprem içeren yük kombinasyonlarında uygulanır. G+Q kombinasyonlarında brüt kesit rijitlikleri kullanılır. Bu durum SAP2000/ETABS'ta iki ayrı analiz seti gerektirmektedir; yanlış uygulama iç kuvvetlerin hatalı hesabına yol açar.

6.3 Diyafram Modelleme (TBDY 2018 §4.5.6)

Tablo 6: Diyafram Modelleme (TBDY 2018 §4.5.6)

7. Ağ Bağımsızlık Testi (Mesh Convergence)

SEY sonuçları ağ yoğunluğuna bağlıdır. Ağ yeterince rafine edildiğinde sonuçlar yakınsar.

Tablo 7: Ağ Bağımsızlık Testi (Mesh Convergence)

h-refinement vs. p-refinement:

• h-refinement: Eleman boyutu küçültülür → daha fazla düğüm, aynı eleman tipi
• p-refinement: Eleman boyutu sabit, şekil fonksiyonu derecesi artırılır (Q4 → Q8)

Ağ bağımsızlık testi yalnızca gerilme yoğunlaşması olan bölgelerde gereklidir. Saint-Venant ilkesi gereği yerel gerilme lokal etkidir; tüm modelin ağını inceltmek verimsizdir.

8. Türkiye Saha Koşulları

8.1 Deprem Etkisi ve TBDY 2018

Türkiye'nin %96'sı deprem kuşağındadır; SEY analizlerinde deprem yükü çoğunlukla boyutlandırmayı yöneten yük olmaktadır. TBDY 2018 ile SEY'in üç ana uygulama alanı:

1. Modal Hesap Yöntemi (Mod Birleştirme) — TBDY 2018 §4.8: CQC veya SRSS ile modal katkıların birleştirilmesi
2. Zaman Tanım Alanında Doğrusal Hesap — TBDY 2018 §4.8.2: Gerçek deprem kayıtları
3. Doğrusal Olmayan Analiz (İtme/Zaman Tanım) — TBDY 2018 §5: Plastik mafsal, hasar değerlendirmesi

8.2 Zemin Koşulları ve Eleman Seçimi

Tablo 8: Zemin Koşulları ve Eleman Seçimi

9. Doğrusal Olmayan SEY

Deprem etkisi altında betonarme yapılarda doğrusal olmayan davranış iki kaynaktan gelir:

1. Malzeme Doğrusal Olmayanlığı: Beton çatlama, donatı akması (TBDY 2018 §5.4)
2. Geometrik Doğrusal Olmayanlık: P-Δ etkisi (TBDY 2018 §4.9)

Newton-Raphson iterasyon şeması:

$$\{R\}^{(k)} = \{F_{ext}\} - [K_T]^{(k-1)}\{u\}^{(k-1)}$$

Tablo 9: Doğrusal Olmayan SEY

Dikkat: TBDY 2018 §5.6 uyarınca performans değerlendirmesinde doğrusal olmayan analiz kullanılıyorsa, en az 11 yer hareketi kaydı (DD-2 düzeyinde ölçeklendirilmiş) ile zaman tanım alanında analiz yapılmalıdır.

10. Çözümlü Problemler

Problem 1: Tek Elemanlı Çelik Çubuk (Temel Formülasyon)

Veriler:

• $L = 1.0$ m, $A = 4.0 \times 10^{-4}$ m², $E = 200$ GPa (çelik)
• Sol uç ankastre ($u_1 = 0$), sağ uçta $F = 20$ kN aksiyel çekme

Çözüm:

Eleman rijitlik matrisi:

$$[K_e] = \frac{EA}{L} \begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix} = 8 \times 10^7 \begin{bmatrix} 1 & -1 \\ -1 & 1 \end{bmatrix} \text{ N/m}$$

Sınır koşulları uygulandıktan sonra ($u_1 = 0$):

$$8 \times 10^7 \times u_2 = 20 \times 10^3$$ $$u_2 = \frac{20 \times 10^3}{8 \times 10^7} = 0.25 \text{ mm}, \quad \sigma = \frac{F}{A} = 50 \text{ MPa}$$

Kontrol: $\sigma = 50$ MPa $< f_y = 355$ MPa (S355 çelik, TS EN 10025-2) $\checkmark$

Problem 2: Betonarme Kiriş — SEY ve Analitik Karşılaştırma

Veriler:

• Basit mesnetli kiriş: $L = 6$ m, $b = 300$ mm, $h = 500$ mm, C30 beton
• $E_c = 32\,000$ MPa (TS 500:2000 Denklem 3.1)
• Yayılı yük $q = 20$ kN/m

Çözüm:

Atalet momenti:

$$I = \frac{bh^3}{12} = \frac{300 \times 500^3}{12} = 3.125 \times 10^9 \text{ mm}^4$$

Eğilme rijitliği (brüt kesit — düşey yük kombinasyonu):

$$EI = 32\,000 \times 3.125 \times 10^9 = 1.0 \times 10^{14} \text{ N\cdot mm}^2 = 1.0 \times 10^5 \text{ kN\cdot m}^2$$

Analitik çözüm:

$$v_{maks} = \frac{5qL^4}{384EI} = \frac{5 \times 20 \times 6^4}{384 \times 1.0 \times 10^5} = 3.375 \text{ mm}$$

2 elemanlı SEY çözümü ($L_e = 3$ m, kiriş eleman):

$$[K_e] = \frac{EI}{L_e^3} \begin{bmatrix} 12 & 6L_e & -12 & 6L_e \\ 6L_e & 4L_e^2 & -6L_e & 2L_e^2 \\ -12 & -6L_e & 12 & -6L_e \\ 6L_e & 2L_e^2 & -6L_e & 4L_e^2 \end{bmatrix}$$

SEY sonucu: $v_2 \approx 3.37$ mm (analitik ile fark < %0.2%)

Kontrol: TS 500:2000 §12.4: $v_{maks} = 3.375$ mm $\ll L/250 = 24$ mm $\checkmark$

Problem 3: 6 Katlı Karma Yapı — Etkin Rijitlik ve Modal Periyot

Senaryo: 6 katlı çerçeve+perde betonarme bina; $H_{kat} = 3.0$ m, $H = 18$ m; DD-2 deprem düzeyi; ZC zemin; $S_{DS} = 0.95g$; $S_{D1} = 0.40g$

Çözüm:

Adım 1 — Kesit özellikleri (C30, $E_c = 32\,000$ MPa):

Tablo 10: Problem 3: 6 Katlı Karma Yapı — Etkin Rijitlik ve Modal Periyot

Adım 2 — TBDY 2018 Tablo 4.2 etkin rijitlikler:

$$EI_{kolon,etkin} = 0.70 \times 166.7 = 116.7 \text{ MN\cdot m}^2$$ $$EI_{kiriş,etkin} = 0.35 \times 172.8 = 60.5 \text{ MN\cdot m}^2$$

Adım 3 — Periyot üzerindeki etki ($T \propto 1/\sqrt{K}$, kütle sabit):

$$\frac{T_{etkin}}{T_{brüt}} = \sqrt{\frac{k_{brüt}}{k_{etkin}}} \approx \sqrt{\frac{1}{0.50}} = 1.41$$

Adım 4 — TBDY 2018 ampirik periyot (Denklem 4.22):

$$T_A = C_t \cdot H^{0.75} = 0.1 \times 18^{0.75} = 0.849 \text{ s}$$

Etkin kesit ile beklenen modal periyot: $T_{1,etkin} \approx 1.05\text{–}1.30$ s

Adım 5 — Taban kesme kuvvetine etkisi:

$S_{ae}(T_{1,etkin})$ azalan koldan hesaplandığında ($T$ artışı → $S_{ae}$ düşer):

$$V_t = m \cdot \frac{S_{ae}(T)}{R}$$

Brüt kesit kullanımı bu durumda daha kısa periyot → daha yüksek $S_{ae}$ → aşırı kesme kuvveti anlamına gelir.

TBDY 2018 §4.7.3 kontrolü: $T_1 \leq 1.4 T_A = 1.4 \times 0.849 = 1.189$ s sınırı aşılırsa Mod Birleştirme Yöntemi zorunludur.

Sonuç: Etkin kesit rijitliği çarpanları hem modal periyotu (%41 uzama) hem de iç kuvvet dağılımını önemli ölçüde değiştirir. Brüt kesit kullanımı iç kuvvetleri yanıltıcı biçimde değerlendirmeye yol açar.

11. Sık Yapılan Hatalar

Tablo 11: Sık Yapılan Hatalar

12. Temel Formüller Özeti

Tablo 12: Temel Formüller Özeti

Referanslar

• TBDY 2018 (RG 18.03.2018) — §4.5.2–4.5.8, Tablo 4.2; §4.8, §5.4–5.6
• TS 500:2000 — §3.1, §7.3.2, §12.4
• TS EN 1992-1-1:2004 — §3.1.3, §5.4
• TS EN 1997-1:2012 — Bölüm 6, 9 (Zemin-yapı etkileşimi)
• TS 498:1997 — Karakteristik yük değerleri
• TS 708:2016 — B420C, B500C donatı özellikleri
• Zienkiewicz, O.C. & Taylor, R.L. — The Finite Element Method, Vol. 1–3
• Logan, D.L. — A First Course in the Finite Element Method, 6. Baskı
• Bathe, K.J. — Finite Element Procedures, MIT Press

Kaynaklar

1. BYKHY 2015.
2. TS EN 1991-1-2 — CEN — Avrupa Standardizasyon Komitesi (Eurocode). https://eurocodes.jrc.ec.europa.eu
3. Yapısal Analiz.

İlgili Hesaplama Araçları

Bu konuyla ilgili ücretsiz mühendislik hesaplama araçlarımızla ön tasarım ve kontrol yapabilirsiniz:

• Yük Kombinasyonu Hesaplama
• Deprem Yükü Hesaplama (TBDY 2018)

Önemli Mühendislik Uyarısı: Bu içerik yalnızca bilgilendirme amaçlıdır; nihai tasarım, hesap ve uygulama kararları, güncel yönetmelikler ile proje koşulları çerçevesinde yetkili bir inşaat mühendisinin denetiminde alınmalıdır. Sayısal örnekler ve formüller genel mühendislik pratiğini yansıtır; her projenin kendine özgü zemin, yük ve çevre koşulları proje müellifince ayrıca değerlendirilmelidir.