Iyi geceler. Doktora tezim icin 9 katli moment-resisting frame (MRF) uzerinde Incremental Dynamic Analysis (IDA) calisiyorum, Vamvatsikos-Cornell 2002 metodolojisi. OpenSees'te TCL script'imi kurdum ama runtime cok uzun — 30 kayit × ortalama 12 IM basamagi = 360 NLTH kosum, tek kayit 22 dakika. Toplam 132 saat. Dizustumu ofiste birakip bahçede oturma konforuna erisemeyecegim galiba.

Script'i paylasiyorum bir parcasini, elestirilmesini istiyorum — belki bir seyleri yanlis yapiyorum, belki optimizasyon mumkun:

# Hunt & Fill: hunt factor 0.2, fill step 0.05 (Vamvatsikos-Cornell oneriler)\nset hunt_factor 0.2\nset fill_step 0.05\n# Kayit olcekleme: Sa(T1,5%) icin target, T1=1.45 s\n# Recorder: max interstory drift (9 drift, her kata)\n# Convergence: Krylov-Newton, tol 1e-6, maxIter 100

Ozellikle "flatline" tespit algoritmasinda sikintidayim — IDA egrisi dikey asimptota yaklasinca convergence kaybediyor, bu collapse mi yoksa numerik sorun mu? Vamvatsikos scale-back veriyor ama bazen 10 iterasyon doner duruyor. OpenSees'te bu tecrubesi olan var mi?

bim_berk yazdi:

bazen 10 iterasyon doner duruyor

Berk, klasik IDA problem. Convergence kaybi gercekten collapse olabilir de, rigid-body mode olabilir, veya P-Delta instability. Ayirt etmek icin:

(1) Maksimum drift > %10 (tipik collapse proxy) + ardindan convergence kaybi = gercek collapse, flatline dogru. (2) Maksimum drift < %5 ama convergence kaybi = numerik sorun, hunt factor 0.2'yi 0.1'e indir, fill step 0.05'i 0.025. (3) Maksimum drift %5-10 arasi = belirsiz, scale-back yapip tekrar kosma, hala yakinsamiyorsa "collapse" ilan ediyoruz.

Runtime kismi: OpenSees TCL single-thread. Python wrapper (OpenSeesPy) + multiprocessing ile 30 kaydi paralel koşabilirsin (her CPU core bir kayit), 8-core makinede 132 saat / 8 = 16.5 saate duser. Bahçede oturursun yani.

Berk, Furkan'in paralel onerisi dogru, ben de ekleyeyim: Krylov-Newton yerine "Newton -initial" bazen yakinsamayan noktalarda daha saglam. Ben OpenSeesPy'de su algorithm fallback zincirini kullaniyorum: Krylov-Newton → Newton → ModifiedNewton → NewtonLineSearch → BroydenFirstFormula. Hiçbiri yakinsamazsa "non-converged, IM=X, record=Y" logluyorum ve bir sonraki IM seviyesine geciyorum.

Bir de kayit olceklemesi: Sa(T1,5%) mi Sa_avg (Eads et al. 2015) mi? Sa_avg daha dusuk dispersion veriyor, ama akademik tartismasi hala suruyor.

Berk, IDA konusuna vallahi duskunum. Vamvatsikos-Cornell 2002 klasik ama bu is 20 yilda evrildi, FEMA P-58 (2018) aslinda Multiple Stripe Analysis (MSA)'ye cok kullanisli bir alternatif sunuyor. IDA her kayda ayri basamaklar ile zengin egri veriyor, MSA hedef IM seviyelerinde toplu stripe — hesap maliyeti genelde yari yariya, dispersion degerleri yakin. Tezin zorunlu olarak IDA icin mi yaziliyor, yoksa MSA da kabul gorur mu?

Yine Baker-Cornell (2005) conditional mean spectrum yaklasimi kayit secimi acisindan Sa(T1) standart olceklemeyi guzel tamamliyor, spektral bias azaliyor. Tezin danismaniyla tartismadiysan aklina dusur.

ank_ezgi yazdi:

Multiple Stripe Analysis (MSA)'ye cok kullanisli bir alternatif sunuyor

Ezgi ablam, danismanim aslinda ilk MSA oneriyordu, ben "IDA daha gorsel" dedim inat ettim. Simdi runtime duvari'na tosladim, geri donmeyi dusunmek lazim galiba. Tezimin amaci collapse fragility curve cikarmak, her ikisi de gider.

CMS onerin icin tesekkurler, PEER NGA-West2'den kayit sectim Sa(T1) olcekleme ile ama CMS'e gecip bias'i duzeltmek danisma toplantisinin gundemine giriyor.

Berk, pratik bir not — OpenSees'te 9 kat MRF tek kayit 22 dakika biraz yuksek goruniyor bence, element sayina bagli. Fiber section element sayisi cok mu? Nonlinear Beam-Column ile force-based formulation + 5 integration point kafi genelde, 8+ koyarsan gereksiz yuk. Bir de model damping Rayleigh mi modal mi? Modal damping mass-proportional + stiffness-proportional yerine modal-tangent kullaniyorsan runtime ciddi dusuyor.

Benim 8 katli modelde tek NLTH 7-8 dakika oluyor eski Xeon workstation'da. Kontrol et.

Pinar, integration point 7 koydum perdelerde, kirislerde 5 — perdelerde 5'e indirmem gerek belki. Damping Rayleigh (alpha-beta) yapiyorum, modal-tangent'e gecebilirim. Tesekkurler, gelecek hafta calistiririm bu optimizasyonlari, 22 dakika 12 dakikaya duserse paralel ile beraber cok rahatlarim.

Berk, doktora tezi derken bir nefes al ve basarilar. Ben tezim icin 2012'de SAP2000 ile IDA yapmaya calismistim, rezilce bir sey cikmisti — SAP2000 o zamanlar NLTH otomasyon yoktu, her kayit elle reset, dizustu haftasonu calisiyordu, dondurucu kahve. OpenSees dogdugunda "Allah seni korusun David McKenna" demistim.

Tezi yazmaya basladiginda fragility curve tarafi zevkli, analiz monoton. Geziye fragility yaparken konuyu takip ederim.

132 saat → 16.5 saat → belki 10 saat + 30 kayit fragility + danisman "MSA deseydin" + IDA asimptot hunt&fill = kara cuma. Berk cok çay iç.

Bilmedigim konu, sessiz okuyorum. Ama 9 kat MRF icin IDA yapmak derken cocuk tam 30 kayit kosuyor, tezi saglam olur. Tebrikler Berk, takip ediyorum.

Berk, Perform-3D'de IDA yari-otomatik macro var (NLTH batch, IM sweep), OpenSees kadar esnek degil ama hunt&fill benzeri bir algoritma ile cikti verir. Runtime daha yavas muhtemelen ama script/debug tarafi yok. Doktora isinde OpenSees kalacak herhalde, uygulama projesinde Perform-3D.

Yapisal tez, iyi is. Soru: IDA kayit setinde zemin sinifi uyumu nasil yonetiliyor? Kayitlar farkli zemin sinifindan geldigi zaman IM olcekleme tek basina yetmiyor. Deaggregation + site-specific scaling dusunuldu mu? Yoksa yine akademik olarak kayit seti "genel" mi?

geoteknik_selim yazdi:

Kayitlar farkli zemin sinifindan geldigi zaman IM olcekleme tek basina yetmiyor

Selim abi, 30 kayit hepsi Vs30 ZC-ZD araliginda (180-760 m/s) sectim, PEER NGA-West2 filtrelenmis. Site-specific scaling akademik tarafta tartismali — Baker-Lee (2018) "records as-is" yaklasimini oneriyor, spektral bias CMS ile duzelt. Tezimde CMS ekliyor olacagim.

Herkese tesekkurler — yapacaklarim: (1) OpenSeesPy + multiprocessing, (2) algorithm fallback zinciri (Begum'dan), (3) modal-tangent damping + perde integration point 5 (Pinar'dan), (4) CMS olcekleme danisman goruşu sonrasi, (5) IDA + MSA paralel kosum sansli olursam tezde ikisi birden. Sonuc ciktiginda forumda paylasirim, bir ay sonra tahmini.