Bu yazımızda, modern inşaat mühendisliği yazılımlarının kalbinde yatan Sonlu Elemanlar Metodunun ne olduğunu ve bu metodu kullanan programların (SAP2000, ETABS, ANSYS vb.) arka planda nasıl çalıştığını inceleyeceğiz.

Bu tarz programları kullanan her mühendisin, arka planda dönen matematiği ve seçtiği eleman tiplerinin (çubuk, kabuk vb.) ne anlama geldiğini detaylıca bilmesi gerekir. Üniversitelerde iki dönem boyunca okutulan bu devasa konunun tamamını bir yazıya sığdırmak imkânsız olsa da, bu yazıda en temel ve kritik noktalara değineceğiz.

Direkt Rijitlik Yöntemi

Direkt rijitlik yönetminde düğüm noktaları ve eleman bağlantı şemасı
Görsel 1 - Figure 1 - Sonlu elemanlar metodunda düğüm noktaları ve mesh ağ bağlantı şeması

Sonlu elemanlar metodunu kullanan yapısal analiz programlarının temeli Direkt RijitlikRijitlikRijitlik, bir elemanın veya yapının yük altında ne kadar az şekil değiştirdiğini gösteren özelliktir; aynı yük altında daha az yer değiştiren sistem daha rijittir. Yöntemi’ne dayanır. Bu yöntemde karmaşık ve sürekli bir sistem, sonlu sayıda idealize edilmiş küçük parçalara (elemanlara) bölünür. Bu elemanların birbirine bağlandığı noktalara ise düğüm noktasıDüğüm NoktasıYapısal analiz modelinde çubuk, kiriş veya kabuk elemanlarının birbirine bağlandığı, konumu ve serbestlik dereceleri (deplasman ve dönmeleri) tanımlanan referans noktası. denir.

Sistemin davranışı şu meşhur matris denklemi ile ifade edilir:

F=[K]⋅u{F}=[K]⋅{u}

Bu denklemde;

  • F: Dış kuvvetler (yük) vektörünü,

  • [K]: Sistemin rijitlik matrisini,

  • u: Düğüm noktalarındaki yer değiştirmeler (deplasmanDeplasmanBir yapı noktasının, yük etkisi altında ilk konumuna göre yaptığı yer değiştirme; yapısal analizin temel çıktılarından biri.) vektörünü ifade eder.

Program, yapıya etkiyen yükleri F ve yapının rijitliğini [K] tanımladıktan sonra, denklemi çözerek u yer değiştirmelerini bulur. Bu yöntemde bilinmeyen değerler sadece düğüm noktalarında hesaplanır.

Lokal ve Global (Sistem) Rijitlik Matrisleri

Yapımızı oluşturacak her bir elemanın kendi eksenine göre bir lokal rijitlik matrisiRijitlik MatrisiBir yapı sistemindeki düğüm noktası yer değiştirmeleri ile bu noktalara etkiyen kuvvetler arasındaki ilişkiyi gösteren matris. vardır. Örneğin, üç boyutlu uzayda bir çerçeve (kirişKirişAğırlıklı olarak eğilme momenti ve kesme kuvveti taşıyan, döşeme yüklerini kolon veya perdelere aktaran yatay taşıyıcı eleman./kolonKolonAğırlıklı olarak eksenel basınç kuvveti taşıyan, düşey ve yatay yükleri kirişten temele aktaran düşey taşıyıcı eleman.) elemanının her bir düğüm noktasında 6 adet serbestlik derecesiSerbestlik DerecesiBir yapısal sistemin konumunu tam olarak tanımlamak için gereken bağımsız yer değiştirme (öteleme/dönme) bileşeni sayısı. (3 öteleme, 3 dönme) bulunur. İki ucu olan bir çubuk elemanın lokal rijitlik matrisi 12×12 boyutlarındadır.

Ancak analiz lokal eksenlerde değil, tüm yapıyı kapsayan genel (global) eksenlerde yapılır. Bu nedenle her elemanın lokal rijitlik matrisi, bir rotasyon (dönüşüm) matrisi ile çarpılarak Global (Sistem) Rijitlik Matrisine dönüştürülür ve birleştirilir. Sisteme etkiyen yük vektörü ve global rijitlik matrisi oluşturulduktan sonra düğüm noktalarının yer değiştirmeleri topluca hesaplanır.

Galerkin Metodu ve Yaklaşık Hesap (Şekil Fonksiyonları)

Bernstein polinomu ile fonksiyon yaklaşımı animasyonu
Görsel 2 - Figure 2 - Bernstein polinomu ile fonksiyon yaklaşımı: polinom derecesi arttıkça gerçek çözüme yaklaşma animasyonuKaynak: Source: James pic, Wikimedia CommonsLisans: License: Public Domain

Sadece düğüm noktalarındaki değerleri bildiğimizi söylemiştik. Peki iki düğüm noktası arasındaki, yani elemanın ortasındaki bir noktada meydana gelen gerilmeyi veya yer değiştirmeyi nasıl buluyoruz?

İşte burada diferansiyel denklemleri basitleştiren Galerkin Metodu gibi sayısal yaklaşımlar devreye girer. Çözmemiz gereken gerçek problem 3. dereceden karmaşık bir eğri (polinom) olabilir. Ancak biz, iki düğüm noktası arasına 1. veya 2. dereceden daha basit şekil fonksiyonları (shape functions) atayarak yaklaşık bir hesap yaparız.

Sonlu elemanlar metodunda 1D şekil fonksiyonları ve yaklaşık çözüm gösterimi
Görsel 3 - Figure 3 - 1D şekil fonksiyonları (shape functions): iki düğüm noktası arasında polinom interpolasyonu ile yaklaşık çözüm

Kullanılan polinomun derecesi ve eleman sayısı (mesh sıklığı) artırıldıkça, yaklaşık sonuçlarımız analitik gerçek sonuçlarla birebir örtüşmeye başlar. FEM’in ana mantığı budur: Çözümü imkânsız veya çok zor olan büyük bir problemi, kolay çözülebilir binlerce küçük parçaya (elemana) bölmek ve yaklaşık ama çok hassas bir sonuç elde etmek.

FEM’de Kullanılan Temel Eleman Çeşitleri

Doğru analizin sırrı, doğru elemanı seçmekten geçer. Programı kullanan mühendis ne yaptığını biliyorsa gerçeğe çok yakın sonuçlar alır; aksi takdirde “çöp girer, çöp çıkar” (garbage in, garbage out) kuralı işler.

Kafes Eleman: Sadece eksenel (normal) kuvvet taşır. Düğüm noktalarında moment aktarmaz (mafsallıdır).

Kafes/Çubuk Eleman rijitlik matrisi şeması
Görsel 4 - Figure 4 - Kafes/çubuk eleman: lokal rijitlik matrisi ve düğüm serbestlik dereceleri

Çerçeve Elemanı: Kolon ve kirişleri modellemek için kullanılır. Düğüm noktası başına altı serbestlik derecesi bulunur. Eksenel kuvvetNormal KuvvetBir kesite, kesit düzlemine dik (eksen doğrultusunda) etkiyen, çekme veya basınç olarak ortaya çıkan iç kuvvet., kesme kuvvetiKesme KuvvetiBir kesitin bir tarafındaki kuvvetlerin, kesite paralel bileşkesi; kesitte kayma gerilmesi doğuran iç kuvvet., eğilme ve burulmaBurulmaBir elemanın ekseni etrafında dönme eğilimi yaratan, kesitte kayma gerilmesi doğuran dış moment etkisi. momentlerine karşı rijitliği vardır.

Çerçeve elemanı: global rijitlik matrisine dönüşüm ve lokal-global eksen bağlantısı şeması
Görsel 5 - Figure 5 - Çerçeve elemanı: global rijitlik matrisine dönüşüm ve lokal-global eksen bağlantısı

Levha Eleman: Sadece düzlem içi rijitliği vardır, düzlem dışı rijitliği sıfırdır (kağıt gibi bükülebilir). Yalnızca membran kuvvetlerini taşır.

Levha (membrane) eleman: düzlem içi membran kuvvetleri ve yük vektörü dağılımı şeması
Görsel 6 - Figure 6 - Levha (membrane) eleman: düzlem içi membran kuvvetleri ve yük vektörü dağılımı

Plak (Plate) Eleman: Levhanın aksine sadece düzlem dışı bükülme (eğilme momentleri ve kesme) rijitliğine sahiptir.

Plak (plate) eleman: düzlem dışı bükülme rijitliği ve FEM çözüm adımları şeması
Görsel 7 - Figure 7 - Plak (plate) eleman: düzlem dışı bükülme ve FEM çözüm adımları

Kabuk Eleman: Levha ve Plak elemanların birleşimidir. Hem düzlem içi hem de düzlem dışı rijitlikleri bulunur.

Modellemede Ağ (Mesh) Yapısı ve Hassasiyeti

Sonlu elemanlar programlarında yapılan en yaygın ve tehlikeli hatalardan biri, sisteme “Otomatik Mesh At” deyip analiz sonuçlarını doğrudan kabul etmektir. Sisteminizi ne kadar küçük elemanlara (mesh) bölerseniz, Galerkin metodu gibi yaklaşımların ürettiği matematiksel hata payı o kadar düşer ve gerçek analitik sonuca o kadar yaklaşırsınız.

Ancak burada kritik bir denge vardır:

  • Kaba Ağ: Elemanların çok büyük seçilmesidir. Kaba bir mesh, matematiksel doğası gereği yapıyı olduğundan daha rijit gösterir. Bu da sistemde oluşacak yer değiştirmeleri ve gerilmeleri eksik (yanlış) hesaplamanıza neden olur.

  • İnce Ağ: Elemanların çok küçük seçilmesidir. Sonuçlar çok hassastır ancak bilgisayarın çözeceği denklem (matris) boyutu devasa boyutlara ulaşır ve analiz günlerce sürebilir.

Bu yüzden profesyonel mühendisler Mesh Bağımsızlık Testi (Mesh Convergence Study) yaparlar. Mesh boyutunu kademeli olarak küçültüp sonuçların (örneğin maksimum çökmenin) artık değişmediği o optimum noktayı (optimum mesh boyutunu) bulur ve analizi orada tamamlarlar.

Dinamik Analizlere Geçiş: Kütle ve Sönüm Matrisleri

Yazımızın başlarında sistemin davranışını F=[K]⋅u formülü ile, yani sadece Rijitlik Matrisi [K] üzerinden tanımlamıştık. Ancak bu formül sadece rüzgarsız, depremsiz, durağan (Statik) yükler için geçerlidir.

Eğer bir önceki yazımızda bahsettiğimiz gibi yapıya Zaman Tanım Alanında AnalizZaman Tanım Alanında AnalizYapının hareket denklemini gerçek veya yapay bir deprem ivme kaydına karşı her zaman adımında çözerek tepkisinin tam olarak hesaplandığı dinamik analiz yöntemi. yapacaksanız, yani deprem gibi dinamik ve değişken bir yük etki edecekse, bu statik denklem yetersiz kalır. Dinamik analizlerde program, rijitlik matrisinin yanına yapının ataletini temsil eden Kütle Matrisini [M] ve enerjiyi sönümlemeSönümlemeBir yapının titreşim enerjisini zamanla dağıtarak sönümleyen, salınım genliğini azaltan içsel veya eklenmiş direnç mekanizması. kapasitesini temsil eden Sönüm Matrisini [C] de eklemek zorundadır.

Denklem artık şu meşhur Dinamik Hareket Denklemine dönüşür:

[M]⋅{u¨}+[C]⋅{u˙}+[K]⋅{u}={F(t)}[M] \cdot \{\ddot{u}\} + [C] \cdot \{\dot{u}\} + [K] \cdot \{u\} = \{F(t)\}

Burada {u¨}\{\ddot{u}\} ivme vektörünü, {u˙}\{\dot{u}\} ise hız vektörünü temsil eder. Yani sonlu elemanlar programı, depremin her bir saniyesi için bu devasa matrisleri birbiriyle çarparak binanın tepkisini hesaplar.

Doğrusal Olmayan (Non-Lineer) Analizlere Giriş

Şu ana kadar anlattığımız mantıkta rijitlik matrisinin [K] analiz boyunca sabit kaldığını varsaydık. Yani çeliğin veya betonun üzerindeki yük kalktığında eski haline geri döneceği (Elastik/Lineer davranış) kabul edilir. Ancak gerçek bir depremde işler böyle yürümez!

Büyük bir depremde betonBetonÇimento, agrega (kum, çakıl) ve suyun belirli oranlarda karıştırılmasıyla elde edilen, zamanla sertleşerek yüksek basınç dayanımı kazanan kompozit yapı malzemesi. çatlar, donatılar akar (plastik mafsal oluşur) ve yapı kalıcı olarak şekil değiştirir. TBDY 2018’de Şekildeğiştirmeye Göre Değerlendirme (ŞGDT) olarak geçen bu durum, programda Doğrusal Olmayan (Non-Lineer) Analiz yapılmasını gerektirir.

Non-lineer bir analizde yapı hasar aldıkça elemanların rijitliği düşer. Sonlu elemanlar programı bu durumu simüle edebilmek için sabit bir [K] matrisi kullanamaz. Yükün her artış adımında, hasar gören elemanların lokal rijitliklerini düşürerek tüm sistemin global rijitlik matrisini [K] baştan hesaplamak zorundadır. Bu da lineer (doğrusal) analizlere kıyasla sonlu elemanlar programlarının işlem gücünü ve çözüm süresini katbekat artıran en üst düzey analiz aşamasıdır.

Sıkça Sorulan Sorular
Sonlu Elemanlar Metodu nedir?

Sonlu Elemanlar Metodu, karmaşık yapısal sistemleri küçük elemanlara ve düğüm noktalarına bölerek sayısal olarak analiz eder. Sistem, elemanların rijitlikleri birleştirilerek oluşturulan denklem takımıyla çözülür.

Sonlu Elemanlar Metodu nasıl çalışır?

Program önce elemanların rijitliklerini, yük vektörünü ve düğüm noktalarını tanımlar. Ardından rijitlik denklemini çözerek düğüm noktalarındaki yer değiştirmeleri ve bunlardan türetilen sonuçları hesaplar.

Kafes, çerçeve ve kabuk elemanları arasındaki fark nedir?

Kafes elemanı yalnızca eksenel kuvvet taşır, çerçeve elemanı ise eksenel kuvvet, kesme, eğilme ve burulma rijitliklerini temsil eder. Kabuk elemanı, levha ve plak davranışlarını birlikte modelleyerek düzlem içi ve düzlem dışı rijitlik sağlar.

Mesh boyutu sonlu elemanlar analizini nasıl etkiler?

Kaba ağ, yaklaşık çözüm hatasını artırabilir ve yapıyı olduğundan daha rijit gösterebilir. Daha ince ağ hassasiyeti artırır ancak çözülmesi gereken denklem sistemini büyütür, bu nedenle mesh bağımsızlık testi yapılır.

Sonlu elemanlar analizinde doğrusal olmayan çözüm neden daha uzun sürer?

Doğrusal olmayan analizde beton çatladıkça ve donatı aktıkça eleman rijitlikleri değişir. Program her yük artışında global rijitlik matrisini yeniden hesapladığı için doğrusal analize göre daha fazla işlem gücü ve süre gerekir.