Modal Analiz Nedir? Yapı Dinamiğinde Kullanımı

Modal analiz; çok serbestlik dereceli yapıların deprem davranışını, her birini bağımsız birer basit sistem gibi ele alarak analiz eden ve ardından istatistiksel mod birleştirme kurallarıyla birleştiren güçlü bir hesap yöntemidir. Bu yazıda hareket denkleminden eigenvalue analizine, modal dönüşümden %95 katılım kuralına ve SRSS/CQC seçimine kadar tüm süreç teknik olarak ele alınmaktadır.

---

Neden Modal Analiz?

Günümüzde birçok inşaat mühendisi, yapısal analiz programlarında günlük olarak Modal Analiz yapıyor. Peki bilgisayarın o kara kutusunun içinde neler yaşandığını biliyor muyuz?

Deprem mühendisliğinde modal analiz neden bu kadar merkezi bir yere sahip? Çünkü gerçek yapılar, binlerce serbestlik derecesine sahip son derece karmaşık sistemlerdir. Bu karmaşıklığı doğrudan çözmek yerine modal analiz, zekice bir matematiksel dönüşümle problemi N adet bağımsız basit probleme parçalar. Her parçayı ayrı ayrı çözer, ardından sonuçları istatistiksel kurallarla birleştirir.

Alternatif yöntem olan Zaman Tanım Alanı Analizi her zaman adımını tek tek hesaplar; bu son derece güçlü ama hesap açısından çok daha maliyetli bir yöntemdir. Modal analiz ise büyük çoğunluğu oluşturan rutin tasarım projelerinde hız, pratiklik ve TBDY 2018 uyumu açısından tercih edilen standarttır.

---

Yapı Dinamiğinin Kalbi: Hareket Denklemi

Yapı dinamiğinin en temel denklemi olan Hareket Denklemi ile başlayalım. Tek serbestlik dereceli bir sistem için bu denklem şu şekildedir:

$$m \cdot \ddot{u}(t) + c \cdot \dot{u}(t) + k \cdot u(t) = P(t)$$

Burada;

• m: Yapının kütlesi,

• c: Yapının sönümleme katsayısı,

• k: Yapının elastik (yay) rijitliği,

• $u(t)$, $\dot{u}(t)$, $\ddot{u}(t)$: Sırasıyla yapının yer değiştirmesi, hızı ve ivmesi,

• P(t): Yapıya dışarıdan etkiyen zamanla değişen dinamik kuvvettir.

Not: Deprem durumu için P(t) yerine eylemsizlik kuvveti olan $-m \cdot \ddot{u}_g(t)$ yazılır.

Gerçek hayatta tasarladığımız yapılar tek bir kütleden ibaret değildir; binlerce serbestlik derecesine sahip çoklu sistemlerdir (MDOF). Bu durumda yukarıdaki terimler tek bir sayı olmaktan çıkar ve devasa matrislere dönüşür:

$$[M] \cdot \ddot{u}(t) + [C] \cdot \dot{u}(t) + [K] \cdot u(t) = P(t)$$

Burada [M] kütle matrisi, [C] sönümleme matrisi, [K] ise rijitlik matrisidir. Her biri (N×N) boyutundadır; N yapının toplam serbestlik derecesi sayısıdır.

Bu devasa matris denklemini Zaman Tanım Alanında saniye saniye çözmek mümkündür. Ancak pratik tasarımda, yapının depremin 17. saniyesinde tam olarak nerede olduğundan ziyade, deprem boyunca yapacağı maksimum deplasman ve kuvvetler bizi ilgilendirir. İşte bu maksimum değerleri çok daha hızlı ve pratik bir şekilde bulmak için denklemi frekans tanım alanına taşırız.

---

Özdeğer (Eigenvalue) Analizi ve Mod Şekilleri

Yukarıdaki [M], [C] ve [K] matrisleri birbirine bağlı (coupled) denklemler üretir. Yani yapının çatısındaki bir hareket, temelindeki bir kuvveti doğrudan etkiler. Bu karmaşık ve birbirine dolanmış N adet denklemi çözmek zordur.

İşte burada imdadımıza Özdeğer (Eigenvalue) Analizi yetişir. Serbest titreşim ve sönümsüz (c=0) durumu için hareket denklemi:

$$[K] \cdot \{\phi\} = \omega2 \cdot [M] \cdot \{\phi\}$$

ya da standart özdeğer formunda:

$$\bigl([K] - \omega2 \cdot [M]\bigr) \cdot \{\phi\} = \{0\}$$

Bu matematiksel işlem sayesinde, N serbestlik derecesine sahip sistemimizin N adet çözümünü buluruz:

• N adet Doğal Açısal Frekans ($\omega$_n) → özdeğerler (eigenvalues): Yapının hangi ritimlerde sallanmak istediğini gösterir.

• N adet Mod Şekli ($\{\phi\}_n$) → özvektörler (eigenvectors): O ritimlerde sallanırken yapının nasıl bir şekil alacağını tanımlar.

Doğal periyot bu frekanslardan elde edilir: $T$_n = 2\pi / \omega_n

Tablo 1 — Mod Şkillerinin Temel Özellikleri

Özellik	Açıklama
Ortogonallik	Farklı modların şekilleri birbirine diktir: φᵢᵀ[M]φⱼ = 0 (i≠j)
Ölçek bağımsızlığı	Mod şekilleri yalnızca şekli tanımlar, genlik normalleşmeye bağlıdır
Sıralama	Birinci mod en uzun periyotlu, yüksek modlar giderek kısalan periyotlu
Fiziksel anlam	Her mod şekli o modda binanın aldığı “dans pozisyonunu” gösterir

---

Modal Dönüşüm: Bağlı Denklemlerden Bağımsız Denklemlere

Özdeğer analizinin güzelliği, bize bu dönüşümü mümkün kılan matematiksel araçları vermesidir. Modal dönüşüm şu adımlarla gerçekleşir:

Adım 1 — Modal Koordinatlara Geçiş:

Fiziksel yer değiştirme vektörü $\{u(t)\}$, mod şekilleri matrisi $[\Phi]$ ve modal koordinatlar $\{q(t)\}$ cinsinden yazılır:

$$\{u(t)\} = [\Phi] \cdot \{q(t)\}$$

Adım 2 — Denklemi Modal Koordinatlara Taşımak:

Bu dönüşüm ve ortogonallik özelliği kullanıldığında matris denklemi şuna dönüşür:

$$M$$_n \cdot \ddot{q}_n(t) + C_n \cdot \dot{q}_n(t) + K_n \cdot q_n(t) = L_n \cdot \ddot{u}_g(t)

Burada Mn, Cn, Kn n. modun modal kütlesi, modal sönümleme ve modal rijitliğidir. Ln ise modal katılım faktörüdür.

Adım 3 — Normalleştirme:

Mn = 1 (birim kütle normalleşmesi) yapıldığında denklem çok daha sade hale gelir:

$$\ddot{q}_n(t) + 2\zeta$$_n\omega_n \cdot \dot{q}_n(t) + \omega_n² \cdot q_n(t) = \Gamma_n \cdot \ddot{u}_g(t)

Burada $\Gamma$_n modal katılım faktörü, $\zeta$_n ise o modun sönümleme oranıdır.

Kritik Nokta: Artık elimizde birbirinden tamamen bağımsız N adet SDOF denklemi var. Her birini ayrı ayrı çözmek, dev matris sistemini çözmekten çok daha kolaydır.

Belirli mühendislik kabulleri yaptığımızda bu bağımsızlık koşulu sağlanır:

1. Davranışın lineer-elastik sınırda kalması

2. Klasik (orantısal) sönümleme matrisine sahip olunması

3. Doğal titreşim periyotlarının birbirine çok yakın olmaması

---

Mod Kütlesi Katılım Oranı ve %95 Kuralı

Modal analizin en kritik pratik sorusu şudur: “Kaç mod almak yeterlidir?”

Bu sorunun yanıtı Mod Kütlesi Katılım Oranı kavramında gizlidir.

Tanım:

n. modun kütlesel katılım oranı şu formülle hesaplanır:

$$\alpha$$_n = \Gamma_n² / M_n

ya da daha pratik ifadeyle:

$$\alpha$$_n = L_n² / (M_n \cdot M_{toplam})

Bu oran, ilgili modun toplam yapı kütlesinin yüzde kaçını “temsil ettiğini” gösterir. Tüm modların katılım oranları toplanırsa %100 elde edilir.

TBDY 2018 — %95 Kuralı (Madde 4.8.1.2):

Türkiye Bina Deprem Yönetmeliği 2018, Madde 4.8.1.2 uyarınca iki koşul birlikte sağlanmalıdır:

(a) (X) ve (Y) deprem doğrultularında her bir mod için hesaplanan taban kesme kuvveti modal etkin kütleleri toplamının bina toplam kütlesinin %95’inden daha az olmaması:

$$\sum m_{txn} \geq 0.95 \cdot m$$_t

ve

$$\sum m_{tyn} \geq 0.95 \cdot m$$_t

(b) Katkısı %3’ten büyük olan bütün modlar ayrıca göz önüne alınacaktır.

Yani yazılım çıktısında kümülatif kütlesel katılım oranı her iki yatay doğrultuda (X ve Y) %95’i geçene kadar mod almaya devam etmelisiniz. Buna ek olarak, tek başına %3’ten fazla kütle katılımı olan herhangi bir mod, kümülatif toplam %95’i geçmiş olsa bile hesaba dahil edilmelidir.

---

Modal Analiz ve Süperpozisyon Mantığı

Bu kabullerin geçerli olduğu durumlarda, koskoca bir binayı birbirinden bağımsız çalışan basit sarkaçlar (SDOF sistemler) gibi düşünebiliriz.

Her bir modu tek tek, sanki diğerleri yokmuş gibi analiz eder, her modun kendi maksimum tepkisini (moment, kesme kuvveti, yer değiştirme vb.) Spektrum Eğrisinden okuyarak kolayca hesaplarız. Yapılması gereken zorlu işlemi parçalara bölerek inanılmaz derecede basitleştirmiş oluruz!

Modal analizin özü şudur: Birden çok serbestlik derecesine sahip sistemlerin karmaşık davranışını, bağımsız basit modların toplamı şeklinde analiz etmek.

Her mod için maksimum tepki hesabı şu formülle yapılır:

$$S_{n,max} = \Gamma$$_n \cdot S_d(T_n, \zeta_n)

Burada Sd(Tn, ζn) o modun periyoduna ve sönümleme oranına karşılık gelen Deplasman Spektrumu değeridir. İvme spektrumu kullanılacaksa:

$$a_{n,max} = \omega$$_n² \cdot S_{n,max} = \Gamma_n \cdot S_a(T_n, \zeta_n)

---

Mod Birleştirme Kuralları: SRSS ve CQC

Geldik en kritik noktaya: Yapının toplam (gerçek) davranışını bulmak için bu farklı modlardaki değerleri nasıl toplayacağız?

Dümdüz toplama (Mutlak Toplam) yapamayız! Çünkü birinci modun maksimum yer değiştirmesi depremin 4. saniyesinde olurken, ikinci modun maksimum tepkisi 12. saniyede gerçekleşiyor olabilir. Bütün modların maksimum değerlerini aynı anda oluyormuş gibi dümdüz toplarsanız, gerçekte asla oluşmayacak kadar abartılı (çok güvenli ama ekonomik olmayan) sonuçlar bulursunuz.

Bu eşzamanlılık sorununu çözmek için istatistiksel Mod Birleştirme Kuralları geliştirilmiştir:

SRSS — Karelerin Toplamının Karekökü

$$R_{toplam} = \sqrt{R$$₁² + R₂² + R₃² + \cdots + R_n²}

Ne zaman kullanılır? Doğal titreşim frekansları birbirinden yeterince uzak olan (iyi ayrışmış) sistemlerde tercih edilir. İki modun periyotları birbirinin %20’sinden fazla farklıysa SRSS yeterlidir. İstatistiksel olarak iyi sonuç verir.

Sınırlılığı: Modlar birbirine yaklaştığında modlar arasındaki çapraz korelasyonu ihmal eder ve hatalı sonuç verebilir.

CQC — Tam Karesel Birleştirme

$$R_{toplam} = \sqrt{\sum_{i} \sum_{j} \rho_{ij} \cdot R$$_i \cdot R_j}

Burada ρᵢⱼ modlar arası korelasyon katsayısıdır:

$$\rho_{ij} = \frac{8\sqrt{\zeta$$_i \zeta_j}\,(\zeta_i + \beta_{ij}\zeta_j)\,\beta_{ij}^{3/2}}{(1 - \beta_{ij}^2)^2 + 4\zeta_i\zeta_j\beta_{ij}(1 + \beta_{ij}^2) + 4(\zeta_i² + \zeta_j²)\beta_{ij}^2}

Ne zaman kullanılır? Titreşim periyotları birbirine çok yakın olan modların birbirini etkilediği yapılarda — özellikle düzensiz veya asimetrik yapılar, burulma etkisi belirgin sistemler.

Neden CQC daha genel? SRSS, CQC’nin özel bir halidir: Modlar birbirinden yeterince uzaksa korelasyon katsayıları ρij ≈ 0 olur ve CQC formülü otomatik olarak SRSS’e indirgenir.

Tablo 2 — SRSS ve CQC Yöntemlerinin Karşılaştırması

Kriter	SRSS	CQC
Periyot ayrışması	Geniş (Tᵢ/Tⱼ < 0.9)	Dar veya geniş (her durumda çalışır)
Burulma baskın yapılar	Yanlış sonuç riski var	Doğru sonuç verir
Hesap karmaşıklığı	Basit	Biraz daha karmaşık
TBDY 2018 tercihi	Uygun koşullarda kabul edilir	Genel olarak önerilen
Yazılım desteği	Evrensel	Evrensel

---

Rayleigh Sönümleme: Sönümleme Matrisi Nasıl Atanır?

Özdeğer analizinde sönümleme matrisini ([C]) ihmal ettiğimizi fark etmiş olabilirsiniz. Bunun nedeni, [C]‘nin modal dönüşümü basitleştiren ortogonallik özelliğini genellikle sağlamamasıdır. Peki sönümleme modal analize nasıl dahil edilir?

Rayleigh Sönümleme (Orantısal Sönümleme) en yaygın çözümdür:

$$[C] = \alpha [M] + \beta [K]$$

Burada α (alfa) ve β (beta) katsayıları, iki seçilen mod için hedef sönümleme oranı (genellikle %5) sağlayacak şekilde belirlenir:

$$\alpha = \frac{2\zeta \cdot \omega$$_i \omega_j}{\omega_i + \omega_j} $$\beta = \frac{2\zeta}{\omega$$_i + \omega_j}

Rayleigh sönümlemenin pratik sonucu:

Her modun sönümleme oranı bu katsayılarla belirli bir frekans bağımlılığı izler:

$$\zeta$$_n = \frac{\alpha}{2\omega_n} + \frac{\beta \cdot \omega_n}{2}

Bu formülden görüldüğü üzere, hedef modların dışındaki modlar için sönümleme oranı farklı değerler alır — bazı modlarda hedef değerden daha az, bazılarında daha fazla sönümleme oluşur. Bu durum, seçilen iki referans modun dikkatli belirlenmesini gerektirir.

Yazılımda pratik uygulama: SAP2000 ve ETABS’ta bu katsayılar “Proportional Damping” veya “Rayleigh Damping” seçeneğiyle tanımlanır. TBDY 2018 bağlamında betonarme yapılar için ζ = %5 standart değerdir.

---

Sonuç

Modal analiz, deprem mühendisliğinin en zarif çözümlerinden biridir. Binlerce serbestlik derecesine sahip dev bir matris problemini, fiziksel anlam taşıyan bağımsız modlara parçalayarak çözülebilir kılmak; hem matematiksel derinlik hem de mühendislik sezgisi gerektiren bir yaklaşımdır.

Bu yazıda öğrendiklerimizi özetlersek:

• Hareket denklemi SDOF’tan MDOF’a taşındığında matris formunu alır ve doğrudan çözümü zorlaşır.

• Özdeğer analizi yapının doğal frekanslarını ve mod şekillerini verir; bu iki bilgi modal analizin temelidir.

• Modal dönüşüm, bağlı denklemleri bağımsız SDOF denklemlere dönüştürür; bu sayede her mod ayrı ayrı analiz edilebilir.

• Mod kütlesi katılım oranı kaç mod almanın yeterli olduğunu belirler; TBDY 2018 bu değerin en az %95 olmasını zorunlu kılar.

• SRSS ve CQC modların maksimum tepkilerini istatistiksel olarak birleştirir; genel tercih CQC olmalıdır.

• Rayleigh sönümleme, sönümleme matrisini kütle ve rijitlik matrislerinin doğrusal kombinasyonu olarak tanımlar.

Kara kutunun içinde dönen bu mekanizmanın farkında olan bir mühendis, hem daha doğru modeller kurar hem de yazılım çıktısını sorgulayarak hataları erkenden yakalar.