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/CQCTam İkinci Dereceden Kombinasyon YöntemiTam İkinci Dereceden Kombinasyon anlamına gelen, SRSS'in bir genellemesi olarak birbirine yakın periyotlu modlar arasındaki istatistiksel korelasyonu da hesaba katan modal kombinasyon kuralı. 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 DenklemiHareket DenklemiBir dinamik sistemin atalet, sönüm ve rijitlik kuvvetleri ile dışarıdan uygulanan yükler arasındaki dengeyi ifade eden, sistemin zamana bağlı tepkisini tanımlayan temel diferansiyel denklem. 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ü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ı. 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 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. matrisidir. Her biri (N×N) boyutundadır; N yapının toplam 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ı. 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 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. 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\} = \omega^2 \cdot [M] \cdot \{\phi\}$$

ya da standart özdeğer formunda:

$$\bigl([K] - \omega^2 \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 FrekansAçısal FrekansBir titreşim hareketinin, radyan/saniye biriminde ifade edilen dönme hızı; doğal frekansın 2π ile çarpılmasıyla elde edilen ve hareket denkleminde doğrudan kullanılan büyüklük. ($\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 - Table 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^2 \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^2 / M_n$$

ya da daha pratik ifadeyle:

$$\alpha_n = L_n^2 / (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 kuvvetiKesme KuvvetiBir kesitin bir tarafındaki kuvvetlerin, kesite paralel bileşkesi; kesitte kayma gerilmesi doğuran iç kuvvet. 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İvme SpektrumuTek serbestlik dereceli bir sistemin farklı doğal periyotlarda belirli bir deprem kaydına verdiği azami ivme tepkisini gösteren eğri. kullanılacaksa:

$$a_{n,max} = \omega_n^2 \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_1^2 + R_2^2 + R_3^2 + \cdots + R_n^2}$$

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^2 + \zeta_j^2)\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ılarda ve burulmaBurulmaBir elemanın ekseni etrafında dönme eğilimi yaratan, kesitte kayma gerilmesi doğuran dış moment etkisi. etkisi belirgin sistemlerde kullanılır.

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 - Table 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ümlemeRayleigh SönümlemeBir yapının sönüm matrisini, kütle ve rijitlik matrislerinin belirli katsayılarla ağırlıklı toplamı olarak tanımlayan, zaman tanım alanında analizlerde yaygın kullanılan matematiksel sönüm modeli. (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ıSönümleme OranıBir sistemin sahip olduğu sönümlemenin, kritik sönümleme değerine oranı; yapının dinamik tepkisinin genliğini belirleyen boyutsuz parametre (ξ). (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.

Sıkça Sorulan Sorular

Modal analiz nedir ve statik analizden farkı nedir?

Modal analiz, yapının doğal frekanslarını ve mod şekillerini özdeğer çözümüyle bulup dinamik tepkisini bu modlarda inceler. Statik analiz yükleri doğrudan dengeye uygularken modal analiz kütle ve rijitliğin dinamik etkisini dikkate alır.

Modal analizde kaç mod dikkate alınmalıdır?

Mod sayısı, her yatay doğrultudaki etkin kütle katılım oranı en az yüzde 95 olacak şekilde seçilmelidir; ayrıca tek başına yüzde 3'ten fazla katılım sağlayan modlar da hesaba katılır. Bu nedenle sabit bir sayı yerine yapının kütle ve rijitlik dağılımına göre yeterli mod alınır.

Modal kütle katılım oranı nedir ve neden yüzde 95 aranır?

Bu oran, seçilen modların toplam kütlenin ne kadarını temsil ettiğini gösterir. Yüzde 95 eşiği, deprem tepkisinin büyük bölümünün hesaba katılmasını sağlar.

SRSS ve CQC mod birleştirme yöntemleri ne zaman kullanılır?

SRSS, birbirinden yeterince ayrık modların tepkilerini istatistiksel olarak birleştirebilir. Yakın frekanslı ve korelasyonlu modlarda CQC tercih edilir.

Sönüm oranı modal analiz sonuçlarını nasıl etkiler?

Sönüm oranı her moddaki titreşim tepkisinin ne kadar azaltılacağını belirler. Rayleigh sönümünde bu oran frekansa bağlı değişebildiği için birleşik sonuç da etkilenir.

Burulma ve yüksek mod etkileri modal analizde nasıl anlaşılır?

Mod şekillerinde dönme bileşenlerinin belirginleşmesi burulma etkisine işaret eder. Katılım oranına ulaşmak için çok sayıda yüksek mod gerekiyorsa bu modların tepkileri ayrıca incelenmelidir.