Bu yazıda, MR sönümleyiciler tarafından kontrol edilen yapıların güvenilirliğini korumak için tasarlanmış, zaman gecikmesi kompanzasyonlu bir Stokastik Optimal Yarı Aktif Kontrol yöntemi (SOSC-PSO) önerilmiştir.
Research Article
Bu yazıda, MR sönümleyiciler tarafından kontrol edilen yapıların güvenilirliğini korumak için tasarlanmış, zaman gecikmesi kompanzasyonlu bir Stokastik Optimal Yarı Aktif Kontrol yöntemi (SOSC-PSO) önerilmiştir.
Manyetorolojik (MR) sönümleyicilerin yarı aktif kontrol sistemlerinde kullanımı önemli bir zorlukla karşı karşıyadır: stokastik uyarımlar altında inşaat mühendisliği yapılarının güvenilirliğini azaltan geri besleme süreçlerinin neden olduğu zaman gecikmesi. Bu makale, bu sorunu ele almak ve yapısal güvenilirliği korumak için Fiziksel Stokastik Optimal kontrol teorisinden (PSO) yararlanarak, zaman gecikmesi telafisi (SOSC-PSO) ile bir Stokastik Optimal Yarı Aktif Kontrol yöntemi önermektedir. Önerilen yöntem, yarı aktif kontrol kuvvetini hem mevcut hem de önceki durumların bir fonksiyonu olarak türetir ve kontrol sürecindeki zaman gecikmelerini telafi eder. Kontrol etkinliğini optimize etmek için, temel parametreler sistem için bir güvenilirlik kriterine göre ayarlanır. Stokastik sismik uyarımlar altında tek serbestlik dereceli ve çok serbestlik dereceli yapılar üzerinde yapılan doğrulama analizleri, zaman gecikmelerinin MR sönümleyicilerin performansını önemli ölçüde bozduğunu göstermektedir. Bununla birlikte, zaman gecikmesi kompanzasyonlu SOSC-PSO yöntemi, kontrol etkinliğini önemli ölçüde artırır ve optimize edilmiş parametrelerle, yapısal kontrol sisteminin güvenilirliğini parametre optimizasyonu olmayan yöntemlerin ötesinde artırır.
Depremler ve aşırı rüzgarlar gibi felaket olayları karşısında mühendislik yapılarının performansının iyileştirilmesi, inşaat mühendisliği topluluğu içinde birincil endişe kaynağı olmaya devam etmektedir. Titreşimleri azaltmak için etkili bir teknik olan yapısal kontrolün, bu tür yapıların hem güvenliğini hem de işlevselliğini arttırdığı kanıtlanmıştır 1,2,3. Son birkaç on yılda, bu amaç için çeşitli ileri yöntemler ve teknolojiler geliştirilmiştir. Bu yöntemler, kontrol cihazlarını çalıştırmak için kullanılan enerji türüne göre genel olarak dört kategoriye ayrılabilir: aktif, yarı aktif, pasif ve hibrit kontrol sistemleri 4,5,6,7.
Aktif kontrolde, gerekli kontrol kuvveti doğrudan kontrol cihazları aracılığıyla uygulanır ve bu da önemli miktarda enerji gerektirir 8,9,10. Yarı aktif kontrol ise, aktif sistemlere kıyasla çok daha az enerji gerektiren, kontrol sinyallerine dayalı olarak kontrol cihazlarının özelliklerinin (sönümleme veya sertlik gibi) ayarlanmasını içerir11. Pasif kontrol, aksine, sisteme herhangi bir harici enerji girişi olmadan enerji dağılımına dayanır 12,13,14. Hibrit sistemler, daha etkili performans elde etmek için aktif/yarı aktif ve pasif kontrol stratejilerinin özelliklerini birleştirir15. Bu yaklaşımlar arasında yarı aktif kontrol, düşük enerji tüketimi ve yüksek verimlilik dengesi nedeniyle özellikle umut verici olarak kabul edilmektedir 16,17,18. Manyetolojik (MR) sönümleyici, üstün dinamik sönümleme özellikleri ile en etkili yarı aktif kontrol cihazlarından biri olarak kabul edilmektedir 19,20,21,22.
Bununla birlikte, MR damperlerini kullanan yarı aktif kontrol sistemlerinde bir zorluk ortaya çıkar, çünkü geri besleme mantığı kaçınılmaz olarak zaman gecikmelerine neden olur. Bu gecikmelere tipik olarak aşağıdakiler dahil olmak üzere çeşitli faktörler 23,24,25 neden olur: (i) sensör verilerinin alınması ve işlenmesi, (ii) kontrolör26 tarafından gerekli kontrol kuvvetinin hesaplanması, (iii) transdüserler aracılığıyla MR damperlerine27 sinyal iletimi ve (iv) MR damperlerinin kendileri tarafından gerçek kuvvet üretimi28. Bu tür zaman gecikmeleri, üretilen kontrol kuvveti ile yapı tarafından beklenen kuvvet arasında tutarsızlıklara yol açabilir ve kontrol etkinliğini önemli ölçüde tehlikeye atabilir29. Ne yazık ki, MR damperleri için mevcut kontrol algoritmalarının çoğu bu gecikmeleri hesaba katmamaktadır.
Ek olarak, katastrofik olayların doğal rastgeleliği nedeniyle, herhangi bir etkili yarı aktif kontrol algoritması, stokastik uyarımlar altında performansı koruyabilmelidir. Klasik bir stokastik optimizasyon tekniği olan Doğrusal Kuadratik Gauss (LQG) kontrol yöntemi, MR damperleri ile donatılmış yapılarda titreşimleri azaltmak için araştırılmıştır. Örneğin, Dyke ve ark. hızlanma geri beslemesi30 kullanarak MR damperleri ile donatılmış yapıların sismik tepkisini iyileştirmek için bir LQG kırpılmış-optimal kontrol stratejisi önerdi. Ying ve ark. MR sönümleyicili doğrusal olmayan yapılar için stokastik ortalama alma ve dinamik programlama tekniklerini içeren, kırpılmamış yarı aktif stokastik optimal kontrol stratejisini tanıttı31. Diğer çalışmalar, taban yalıtımlı binaların19 sismik performansını artırmak için modal tabanlı LQG kontrolü uygulamış ve rüzgar türbini kulelerinde32 titreşim azaltma için etkinliğini analiz etmiştir. Bununla birlikte, beyaz Gauss gürültüsünü varsayan klasik LQG yöntemleri, depremler veya kuvvetli rüzgarlar gibi gerçek dünya uyarımlarının durağan olmayan, Gauss olmayan doğasını modellemek için pek uygun değildir. Bunu ele almak için, geleneksel LQG yöntemlerinin sınırlamalarının üstesinden gelen ve çeşitli stokastik uyarımları ele almak için daha doğru bir çerçeve sunan fiziksel tabanlı stokastik optimal (PSO) kontrol kavramıgeliştirilmiştir 33,34 35. Çalışmalar, PSO tabanlı yarı aktif stokastik optimal kontrolün, stokastik sismik kuvvetlere maruz kalan hem doğrusal hem de doğrusal olmayan yapılar için titreşim azaltmayı önemli ölçüde iyileştirdiğini göstermiştir36.
Bu nedenle, MR damperleri25,37 tarafından kontrol edilen yapıların güvenilirliğini artıran, zaman gecikmelerini telafi eden sağlam bir yarı aktif kontrol yöntemine acil ihtiyaç vardır. Ayrıca, stokastik uyarımlar altında optimum kontrol performansını sağlamak için, önerilen kontrol yönteminin kritik parametrelerini güvenilirliğe dayalı bir yaklaşım kullanarak optimize etmek esastır. Bu nedenle, MR sönümleyicilere sahip yapıların güvenilirliğini artırmak için bu yazıda zaman gecikmesi kompanzasyonlu bir Stokastik Optimal Yarı Aktif Kontrol yöntemi (SOSC-PSO) önerilmiştir.
Access restricted. Please log in or start a trial to view this content.
1. Stokastik optimal yarı aktif kontrol yöntemi
Bir MR damperinin kontrol etkisi, kaçınılmaz zaman gecikmesinden önemli ölçüde etkilendiğinden, kontrol edilen yapının performansını artırmak için zaman gecikmesi kompanzasyon algoritmasına sahip yarı aktif bir kontrol yöntemi geliştirilmiştir. Ayrıca, dış uyarımların doğasında bulunan rastgelelik, dinamik tepkilerin bariz belirsizliklerine neden olur. PSO kontrolü, yapı kontrol sisteminin güvenilirliğini sağlamak için önerilen yöntemin kritik parametrelerini optimize etmek için tanıtıldı.
1.1 Yarı aktif kontrolün kompanzasyon algoritması
Genellik kaybı olmadan, stokastik uyarımlara tabi hareket denkleminin şu şekilde verildiği, zaman gecikmeli MR sönümleyicileri tarafından kontrol edilen bir n-Serbestlik Derecesi (n-DOF) yapısını düşünün:
(1)
Burada
yapının sırasıyla n x n kütle, sönümleme ve sertlik matrislerini temsil eder; n, kontrol edilen yapının serbestlik derecesi sayısını temsil eder. Yapının yer değiştirme, hız ve ivme vektörleri sırasıyla ile temsil
edilir. Sembollerin üzerindeki küçük bir ve iki nokta, birinci ve ikinci kez türevleri temsil eder. Us (t - tD), MR damperlerinin zaman gecikmeli kontrol kuvvetinin r boyutlu vektörünü temsil eder ve r, damper sayısını temsil eder; t zamanı temsil eder; tD , bu çalışmadaki tüm sönümleyiciler arasında aynı olduğu varsayılan zaman gecikmesini ifade eder, burada
, burada integral sayısı
ile ve
örnekleme periyodunu belirtir.
stokastik uyarımların p boyutlu vektörünü temsil eder ve
dış uyarımlarla ilişkili rastgeleliği karakterize eden stokastik parametre vektörünü temsil eder. Boyutu,
dış uyarımlar için kullanılan modele bağlıdır, ancak yapının mekanik serbestlik dereceleri ile ilgili değildir. Bs ve D , sırasıyla damperlerin ve harici uyarımların konumlarını belirten n x r ve n x p matrisleridir. Durum uzayı gösteriminde, Eşitlik (1) şu şekilde yazılır:
(2)
Burada 2n boyutlu durum vektörünü temsil eder
;
2n x 2n sistem matrisini temsil eder;
MR sönümleyicilerin 2n x r konum matrisini temsil eder;
Dış uyarımların 2N x P konum matrisini temsil eder. Bu parametreler şu şekilde ifade edilir:
,
,
, (
3)
Nerede
, aynı sıraya sahip bir özdeşlik matrisini
belirtir. Hesaplama kolaylığı için, sürekli durum uzayı denklemi Eşitlik (2) ayrık biçimde şu şekilde ifade edilebilir:
(4)
Zaman noktasının
olarak basitleştirildiği
yer. Ve
,
,
2n x 2n, 2n x r ve 2n x p matrislerini gösterir, bunlar şu şekilde ifade edilir:
,
, (
5)
Nerede
, örnekleme dönemini belirtir.
Aktif kontrol ile benzer etkinliği elde etmek için, zaman gecikmeli MR damper tabanlı kontrol için Hrovat algoritmasına38 dayalı basit ve verimli bir kontrol yöntemi önerilmiştir:
(6)
burada
MR damperi için \ zaman noktasında\
yarı aktif kontrol kuvveti sinyalini temsil eder;
zaman gecikmesi telafisi ile referans aktif kontrol kuvvetini temsil eder;
MR damperinin değiştirilebilir olan maksimum sönümleme kuvvetini temsil eder;
mutlak değer işlem sembolünü temsil eder;
MR damperinin maksimum ve minimum Coulombic kuvvetlerini temsil eder;
viskoz sönümleme katsayısını temsil eder;
pistonun damper silindirine göre hareket hızı olan MR damperindeki hız girişini temsil eder. Eşitlik (6)'da,
MR damperinin tasarlanmış parametrelerini temsil edin.
Denklem (6), MR damperinin yarı aktif kontrol kuvvetinin zaman gecikmeli olarak hesaplanmasını gösterir. Denklem (4)'teki zaman adımındaki
yarı aktif kontrol kuvvetinin
, zaman adımındaki
aktif kontrol kuvveti
ve zaman adımındaki
durum değişkenine
dayalı olarak hesaplandığı görülmektedir. Aktif kontrol kuvvetini
elde etmek için, geleneksel maliyet fonksiyonu şu şekilde belirlenir:
(7)
Burada
sistem durumunun 2n x 2n simetrik pozitif yarı kesin ağırlıklandırma matrisini temsil eder;
kontrol kuvvetinin r x r simetrik pozitif belirli ağırlıklandırma matrisini temsil eder. Dış uyarımın yok denecek kadar az olduğu görülmektedir. Aslında, Eşitlik (7)'ye dayanarak, yapı kontrol sistemi, harici uyarma türü ne olursa olsun optimum kontrol etkisini elde edebilir. Ayrık form olarak, Eşitlik (7)39 olarak ifade edilir:
(8)
Maliyet fonksiyonunun
en aza indirilmesi, koşullu bir aşırı değer problemine yol açar ve aktif kontrol kuvveti
25 ile hesaplanır:
(9)
Burada
, ağırlıklandırma matrisleri ve 37 tarafından belirlenen, zaman adımındaki
durum değişkeni
ve
zaman adımındaki
aktif kontrol kuvveti
için kontrol kazançlarını
belirtin. Kaçınılmaz zaman gecikmesi nedeniyle, elde edilen kontrol kuvveti
yapıya zaman noktalarında
etki eder. Bu nedenle, Eşitlik (9)'daki aktif kontrol kuvveti şu şekilde hesaplanır:
(10) buyurmuştur.
Kavramın netliği için, Denklem (10)'daki yarı aktif kontrol kuvveti, Denklem (6)'ya dayalı olarak aşağıdaki gibi ifade edilir:
(11)
Nerede
, zaman noktasındaki
aktif kontrol kuvvetini belirtir. Eşitlik (11)'i Eşitlik (10) ile değiştirmek,
(12) buyurmuştur.
NOT: Hesaplama
için zaman noktasından
itibaren
tüm durum
değerleri gereklidir. Bu nedenle, aşağıdaki durum tahmin yöntemitanıtılmıştır 40.
(13)
Nerede
, tahmin edilen değişkeni belirtir. Eşitlik (13) yinelemesi ile, 'deki
durum
hesaplanır. Aktif kontrol kuvveti Ua(k) şu şekilde ifade edilir:
(14) buyurmuştur.
Yarı aktif kontrol kuvveti, Eqs'in birleştirilmesiyle elde edilebilir. (6), (12) ve (13).
1.2 Yapı kontrol sisteminin güvenilirlik analizi
Olasılık koruma ilkesine göre, artırılmış sistem
ve
MR sönümleyicili yapı olasılık korunur ve aşağıdaki genelleştirilmiş olasılık yoğunluk evrimi denklemleri (GDEE'ler)38 tarafından yönetilir:
(15) buyurmuştur.
(16) buyurmuştur.
Burada
,
sırasıyla sistem durumunun ve yarı aktif kontrol kuvvetinin karşılık gelen bileşenlerini ve zaman gecikmeli olarak belirtir;
ve
artırılmış örnek sistemlerinin
ortak olasılık yoğunluk fonksiyonlarını belirtir ve
sırasıyla;
'nin
örnek uzayı içindeki örneği belirtir;
ve
) sırasıyla sistem durumu ve yarı aktif kontrol kuvveti bileşenlerinin birinci dereceden zaman kısmi türevlerini belirtir. Ortak olasılık yoğunluk fonksiyonları
ve
Eqs çözülerek elde edilebilir. (15) ve (16) aşağıdaki başlangıç koşullarıyla:
(17) buyurmuştur.
(18) buyurmuştur.
Burada
Dirac delta fonksiyonunu gösterir;
ve
'nin
deterministik başlangıç değerlerini sırasıyla belirtir
;
örneğin
olasılık yoğunluk fonksiyonunu belirtir.
Anlık olasılık yoğunluk fonksiyonları
ve
örnek uzayın alanının tek boyutlu integrali
ile elde edilebilir:
(19) buyurmuştur.
(20) buyurmuştur.
Nerede
, 'nin
dağıtım etki alanını belirtir.
Eqs'ye dayalıdır. (19) ve (20), ilgili fiziksel büyüklüklerin tam olasılıksal bilgilerinin, algoritma parametreleriyle ilişkili ilişkileri tanımlandığı takdirde, kolayca türetildiği görülmektedir. İlgili fiziksel büyüklüklerin güvenilirliği, aşağıdaki ilgili performans fonksiyonu ile hesaplanabilir:
(21) buyurmuştur.
(22) buyurmuştur.
Burada
ve
sırasıyla ilgili durum miktarının ve yarı aktif kontrol kuvvetinin hesaplanan güvenilirliğini belirtir;
kontrol edilen yapının i'inci ilgili fiziksel miktarının eşdeğer uç değer vektörünü belirtir;
j'inci kontrol kuvvetinin eşdeğer uç değer vektörünü belirtir; sembollerdeki '-' şapkası, eşdeğer aşırı değer vektörünü41 gösterir;
kontrol edilen yapının k'inci DOF'unun i'inci ilgili fiziksel miktarını belirtir;
J-inci kontrol kuvvetini belirtir;
,
,
,
,
ve
sırasıyla yapının ve yapıya monte edilen MR damperlerinin ilgili fiziksel büyüklüklerinin ve DOF'unun sayısını belirtir.
ve
ve
eşiklerini
gösterir;
rastgele olayın olasılığını belirtir.
1.3 Parametrelerin optimizasyonu
Genliklerinin
oranı ve
geri besleme kontrol modalitesinde kontrol etkinliğini38 yüksek oranda etkiler. Bu nedenle, en iyi kontrol etkinliğini
elde etmek ve
optimize etmek gerekir. Ayrıca, MR damper tabanlı bir kontrol yöntemi olarak, kontrol etkinliği de Eşitlik (6)'daki cihaz parametrelerinden
etkilenir. Yapı kontrol sisteminin optimal güvenilirliğini elde etmek için, güvenilirliğe dayalı kriter Eşitlik (25)'te önerilmiştir.

Şekil 1: Zaman gecikmesi kompanzasyonu ile stokastik optimal yarı aktif yapı kontrolü için parametre optimizasyonunun akış şeması. Kısaltmalar: GDEE = Olayın genelleştirilmiş diferansiyel denklemi ; PDF = Olasılık yoğunluk fonksiyonu. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.
Özetle, zaman gecikmesi kompanzasyonu ile MR sönümleyicilere dayalı stokastik optimal yarı aktif yapı kontrolünü gerçekleştirmek için iki aşama gereklidir:
1.3.1 Eşitlik (7)'de gösterilen maliyet fonksiyonunun
en aza indirilmesiyle, numuneler anlamında durum geri besleme mantığı elde edilir, yani Eşitlik 12)
1.3.2 Güvenilirliğe dayalı performans fonksiyonunun
en aza indirilmesiyle, istatistik anlamında en uygun parametreler elde edilir. Şekil 1, adım 1.3.2'deki kontrolörün parametrelerinin optimizasyon akış şemasını sunar.
İki döngü katmanı içeren adım 1.3.2'deki aşağıdaki adımlar aşağıdaki gibidir:
1.3.2.1 Stokastik parametre vektörü
ile karakterize edilen stokastik uyarımların olasılık atanmış uzayının bölümlenmesi, ilişkili atanmış olasılıklarla birlikte bir dizi temsili noktanın
tanımlanmasına izin verirPq'lar. Bu yaklaşım, olarak gösterilen stokastik uyarma için numune süreçlerinin verimli bir şekilde oluşturulmasını sağlar, kolayca
elde edilebilir.
1.3.2.2 Yarı aktif yapı kontrol sisteminin parametre optimizasyonu için, maliyet-fonksiyon ağırlıklarının
değerlerini başlatın veya güncelleyin ve
. Zaman gecikmeli ilişkili durum geri besleme kontrol kuvveti, bakınız Eşitlik (12), daha sonra hesaplanır.
1.3.2.3 Eqs'de gösterilen GDEE'leri çözerek yapısal tepkilerin ve kontrol kuvvetinin olasılık yoğunluk fonksiyonlarının (PDF'ler) elde edilmesi. (17)–(20):
Optimizasyon çevrimdışıdır, tamamlanır ancak çevrimiçi değildir. Önerilen yöntemin gerçek yapı kontrol sistemlerinde uygulanabilirliğinde, optimal
ve
ulaşılabilir hale getirilmiştir ve gerçek uygulanabilir süreçte optimizasyona ihtiyaç duyulmamaktadır.
Numune uyarımına tabi tutulan yarı aktif yapı kontrol sisteminin deterministik dinamik analizi, bu sayede ilgili fiziksel büyüklükler
ve türevleri
hesaplanmıştır.
Toplam Varyasyon Azaltma (TVD) özelliklerine sahip modifiye edilmiş Lax-Wendroff şeması gibi bir sonlu fark yöntemi kullanılarak, Genelleştirilmiş Diferansiyel Denklemler (GDEE'ler) çözülebilir ve ortak olasılık yoğunluk fonksiyonları
için sayısal çözümler türetilebilir.
Deterministik dinamik analiz ve sonlu farklar yönteminin yukarıdaki iki adımını tekrarlayarak ve tüm temsili noktaların
üzerinden geçerek, olasılık yoğunluk fonksiyonları toplama ile elde edilebilir:
(23) buyurmuştur.
(24)
Nerede
, örnek noktasıyla
ilişkili alt etki alanının alan ölçüsünü temsil eder. Olayın Genelleştirilmiş Diferansiyel Denklemini (GDEE) çözmek için kullanılan yöntem, olasılık yoğunluk evrimi yöntemi (PDEM)42 olarak bilinir.
1.3.2.4 İlgili fiziksel büyüklüklerin PDF'leri, güvenilirliklerini değerlendirmek için kullanılır ve bunlar daha sonra olasılık kriterinin performans fonksiyonuna
dahil edilir.
1.3.2.5 Parametre optimizasyon süreci için sonlandırma kriterlerinin karşılanıp karşılanmadığını değerlendirin. Koşullar yerine getirilmezse, Adım 1.3.2.2'ye dönün; Eğer öyleyse, zaman gecikmeli kontrol sistemi için en uygun parametreler belirlenebilir. Bu çalışmada, parametre optimizasyonu için parametre güncellemelerini kolaylaştıran ve sonlandırma koşullarını tanımlayan MATLAB'daki genetik algoritma (GA) araç kutusu kullanılmıştır. Genetik algoritma ile, optimal değerler on yineleme içinde elde edilebilir ve yakınsama, yerel minimumma sorunları olmadan stabildir. GA'nın hesaplama maliyeti, parçacık sürüsü optimizasyonundan veya gradyan tabanlı yöntemlerden daha büyük olmasına rağmen, GA, bu makaledeki gibi karmaşık ve farklılaştırılamayan problemlerle başa çıkmada iyidir. Bu nedenle, GA uygulanır. Ve hızlı yakınsama özelliğinden dolayı, GA'nın hesaplama maliyeti kabul edilebilir.
1.3.2.6 Optimum güvenilirlik kontrolünü elde etmek için aşağıdaki performans fonksiyonu
formüle edilmiştir. Parametrelerin
optimal değerlerini belirlemek için simge durumuna küçültün
ve
:
} (25)
(26) buyurmuştur.
Daha önce bahsedildiği gibi, kontrolörün ağırlıklandırma matrisleri simetriktir, burada matristeki
elemanlar yer değiştirmeye, hıza ve bunların etkileşimine atanan ağırlıklara karşılık gelirken, matristeki
elemanlar kontrol kuvveti ile ilgili ağırlıkları temsil eder. Daha önceki çalışmalardan elde edilen bulgulara38 dayanarak, maliyet-fonksiyon ağırlıklarının
konfigürasyonlarının değiştirilmesinin ve
kontrol etkinliği üzerinde minimum etkiye sahip olduğu görülmüştür. Sonuç olarak, bu çalışmada aşağıda özetlendiği gibi basitleştirilmiş bir konfigürasyon benimsenmiştir:
,
(27)
Nerede
, tanımlanacak durum ağırlıklandırma matrisinin katsayısını belirtir;
tanımlanacak kontrol ağırlıklandırma matrisinin katsayısını belirtir;
özdeşlik matrisini belirtir. Kontrol etkisi, değerlerin kendilerinden
değil, oranından
etkilenir. Ayrıca,
Denklem (6), MR damperinin minimum Coulombic kuvvetini belirtir, bu da damperdeki akım girişinin sıfır olduğu anlamına gelir
. Daha sonra, Eşitlik (26) şu şekilde daha da ifade edilebilir: T
(28) buyurmuştur.
En uygun endişe katsayıları nerede
. Mühendislikte makul hale getirmek için, optimum aralıkları MR damper fabrikalarının üretim kapasitesine dayanmaktadır.
Açıkçası, Eşitlik (28) ile elde edilen optimal parametreler, önerilen yarı aktif kontrol algoritmasının zaman gecikmeli optimal kontrol etkinliğini korur, bu da yapı kontrol sisteminin stokastik uyarımlara maruz kalan dengeli optimal güvenilirliği elde etmesini sağlar.
2. Vaka çalışması
Bir MR damperi için zaman gecikmesi telafisi ile önerilen yarı aktif kontrol yönteminin etkinliğini analiz etmek için, bir MR damperine bağlı düzlemsel tek katlı bir kesme çerçevesi yapısı, Şekil 2'de gösterildiği gibi yatay stokastik sismik yer hareketine tabi tutulmuştur. Yarı aktif kontrollü yapı sisteminin parametreleri aşağıdaki gibidir: yapı kütlesi
, doğal dairesel frekans
, sönümleme oranı
. Geçmiş deneyimlere ve piyasadaki mevcut damper kantarlarına göre, yapısal yer değiştirme, hız, ivme ve kontrol kuvveti eşik değerleri sırasıyla 10 mm, 100 mm/s, 1.500 mm/s2 ve 150 kN'dir. Stokastik sismik yer hareketini temsil etmek için, fiziksel olarak motive edilmiş stokastik yer hareketi modeli kullanılmıştır43:
(29) buyurmuştur.
Nerede
, mühendislik sahasındaki zaman alanındaki yer hareketini temsil eder ve
ana kayadaki frekans alanındaki yer hareketini belirtir. Vektör
, mühendislik sahasının yüzeyindeki yer hareketinin stokastik doğasını karakterize eder. Parametreler
, baskın frekans
ve eşdeğer sönümleme
dahil olmak üzere saha toprağının özelliklerini tanımlayan stokastik değişkenlerdir. Vektör
, kaynak özelliklerinden ve yayılma yolundan
etkilenen ana kayadaki yer hareketinin stokastik doğasını temsil eder ve bu aşamada yer alan stokastik değişkenlerin sayısını gösterir.
dairesel frekansı ifade eder ve i hayali birimdir.
Mühendislik sahasının baskın frekansı
ve eşdeğer sönümleme oranı
, saha toprağının dinamik özelliklerini karakterize eden temel parametrelerdir. Bu büyüklüklerin olasılık dağılımı ve istatistiksel parametreleri, belirli bir mühendislik alanı sınıfından toplanan sismik ivme kayıtlarına dayalı olarak belirlenebilir. Açıklama amacıyla, [150, 250] m/s kesme dalgası hız aralığına sahip bir saha sınıfı düşünülmüş ve yer hareketinin tasarım karakteristik periyodu 0,45 s olarak ayarlanmıştır.
Kritik parametreler
,
her ikisi de parametre tanımlaması için log-normal bir dağılımı takip eden karşılıklı bağımsız stokastik değişkenler olarak ele alındı. İstatistiksel parametreler aşağıdaki gibidir: 'nin 
ortalaması sırasıyla 12 rad/s ve 0.1 idi. ve varyasyon
katsayıları sırasıyla 0.42 ve 0.35 idi. Ana kayadaki yer hareketi, 0.20 g'lık bir tepe yer ivmesine karşılık gelen 2 m/s2'lik bir Fourier genliğine sahip bir Gauss beyaz gürültü süreci olarak modellenmiştir. Ana kaya hareketini oluşturmak için kullanılan faz açısı şu şekilde
gösterildi: . Bu nedenle, Sb =1 ve
. Bu yaklaşım, belirli bir aşma olasılığı44 ile koşullu stokastik yer hareketlerinin modellenmesi olarak düşünülebilir. Log-normal dağılım varsayımı, gözlemlenen yer hareketi parametrelerinin çarpık doğasını yakalama kabiliyeti nedeniyle yer hareketi modellemesinde yaygın olarak kullanılmaktadır.
Bu varsayımın duyarlılığını değerlendirmek için, normal ve gama dağılımları da dahil olmak üzere alternatif olasılık dağılımları kullanılarak ek analizler yapılmıştır45,46. Bulgular, genel eğilimler tutarlı kalırken, log-normal dağılımın, özellikle yer hareketi yoğunluğu ölçümlerinin kuyruk davranışını yakalamak için gözlemlenen verilere en iyi uyumu sağladığını göstermektedir
Teğet küreler yöntemi47 sayesinde, ilişkili atanmış olasılıklara sahip 221 temsili noktadan oluşan bir koleksiyon seçildi ve temsili yer ivmeleri sentezlendi. Örnekleme frekansı 50 Hz ve yer hareketlerinin süresi 20.48 s idi. Simüle edilmiş yer hareketine durağan olmayan bir yoğunluk atamak için, düzgün bir modülasyon fonksiyonu kullanıldı ve formülasyonu aşağıdaki gibiydi42:
(30) buyurmuştur.
Burada
sırasıyla 2 ve 16 s alın.
Stokastik sismik yer hareketinin ortalama ve standart sapması ve temsili sismik yer hareketinin zaman geçmişi Şekil 3'te gösterilmektedir. Ortalamanın genliği (0.06 m/s2), standart sapmanın (0.8 m/s2) genliğinin ~%8'i idi, bu da fiziksel olarak motive edilmiş stokastik yer hareketi modelinin sıfır ortalama özelliğini sergilediğini gösteriyor. Bu arada, sismik yer hareketi hem zamansal hem de frekans alanlarında dikkate değer durağan olmayan davranışlar sergiledi.

Şekil 2: Manyetolojik sönümleyicili tek katlı bir kesme çerçevesinin çizimi.
yapı kütlesini, sönümleme oranını ve sertliği belirtir;
stokastik sismik uyarımı belirtir. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 3: Seçilen sismik yer hareketlerinin istatistikleri ve temsili örnekleri. (A) Stokastik sismik yer hareketinin ortalama ve standart sapması; (B) Temsili sismik yer hareketinin zaman geçmişi. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.
2.1 Zaman gecikmesinin etki analizi
Zaman gecikmesinin bir MR damperinin yarı aktif kontrol etkinliği üzerindeki etkisini analiz etmek için, Şekil 4 , kontrol edilen yapının yer değiştirmesini, hızını ve ivme RMS'sini (kök-ortalama-kare yanıtı) göstermektedir. Ağırlıklandırma matrisi oranı ve damper parametreleri
, zaman gecikmesi olmadan optimize edilmiş sonuçlar olan ,
, ve
idi38. Kontrol edilen yapının zaman gecikmeli tepkileri, zaman gecikmeli olmayan kontrollü yapıdan daha büyüktü ve zaman gecikmesinin kontrol etkisi üzerindeki etkisi, zaman gecikmesinin artmasıyla birlikte periyodikliği temsil ediyordu. MR sönümleyici ile yarı aktif kontrol yöntemi bir tür geri besleme yöntemi olduğundan, etkinin periyodikliğinin kontrol edilen yapının doğal periyodu (T = 0.56 s) ile ilişkili olduğu düşünülmüştür.
Maksimum MRS yer değiştirmesi, hızı ve ivmesi 24.6 mm, 270.0 mm/s ve 3111.3 mm/s2 olan kontrolsüz yapının tepkileri ile karşılaştırıldığında, MR sönümleyici kontrollü yapı kayda değer bir etki elde etti. Zaman gecikmesi25 ile aktif kontrollü yapının kararsızlığının aksine, MR damper tabanlı yarı aktif kontrol yöntemi, zaman gecikmeli olsa bile hala bir miktar etki elde etti.

Şekil 4: Kontrol edilen yapının zamanla verdiği yanıtların maksimum kök-ortalama-kareleri. (A) Yer değiştirme; (b) hız; (C) Hızlanma. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.
Zaman gecikmesinin yarı aktif olarak kontrol edilen bir yapının güvenilirliği üzerindeki etkisini analiz etmek için, Şekil 5 , artan zaman gecikmesi ile yer değiştirme, hız ve ivmenin güvenilirliğini göstermektedir. Herhangi bir zaman gecikmesi ile kontrol edilen yapı yanıtlarının güvenilirliği, zaman gecikmesi olmayan değerlere göre daha düşüktü, bu da kontrol edilen yapının güvenilirliğinin zaman gecikmesi ile azaldığı anlamına geliyor. Bu arada, maksimum RMS'ye benzer şekilde, güvenilirlik kontrol etkisi, artan zaman gecikmesi ile periyodiklik gösterdi. Ayrıca, zaman gecikmeli veya gecikmesiz kontrollü yapının güvenilirliği, yer değiştirme, hız ve ivme güvenilirliği 0.0954, 0.1058 ve 0.1111 olan kontrolsüz yapıların değerlerinden daha yüksekti.

Şekil 5: Artan zaman gecikmesi ile kontrollü yapının yanıtlarının güvenilirliği. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.
Özetle, yapının tepkileri yarı aktif kontrol yöntemi ile hafifletilmiştir. Bununla birlikte, RMS veya güvenilirliğinden bağımsız olarak, yarı aktif kontrol yönteminin kontrol etkinliği, zaman gecikmesi ile azalmıştır. Bu nedenle, zaman gecikmesini telafi etmek gerekir.
2.1.1 Zaman gecikmesi kompanzasyon yönteminin analizi
Önerilen zaman gecikmesi kompanzasyon yönteminin etkinliğini analiz etmek için, kontrolsüz (Unc) için yer değiştirme, hız ve ivmenin RMS zaman geçmişleri, kompanzasyon kontrollü olmayan zaman gecikmeli sistem (TDN-SAC-PSO) ve kompanzasyon kontrollü zaman gecikmeli sistem (TDC-SAC-PSO) zaman gecikmesinin 0,1 s olduğu Şekil 6'da karşılaştırılmıştır. Ağırlıklandırma matrisinin ve yarı aktif kontrol yönteminin parametreleri, adım 2.1'deki ile aynıdır.

Şekil 6: Farklı yöntemlerle kontrol edilen yapıların kök-ortalama-kare zaman geçmişi karşılaştırması. (A) Yer değiştirme; (b) hız; (c) hızlanma; (D) Kontrol kuvveti. Kısaltmalar: Unc = kontrolsüz; TDC-SAC-PSO = kompanzasyon kontrollü zaman gecikmeli sistem; TDN -SAC-PSO = kompanzasyon kontrolsüz zaman gecikmeli sistem. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.






































Access restricted. Please log in or start a trial to view this content.
Bu yazıda, zaman gecikmesinin MR damper tabanlı yöntemin kontrol etkinliği üzerindeki etkisini amaçlayan, zaman gecikmesi kompanzasyonlu yarı aktif bir kontrol yöntemi önerilmiştir. Önerilen yöntemde, kritik parametreler güvenilirlik kriterine göre optimize edilmiştir. Kontrol etkinliğini karşılaştırarak, aşağıdaki sonuçlar ele alınmaktadır:
(1) Bir MR damperine dayalı yarı aktif kontrol yöntemi, aktif kontrol yönteminden daha sağlamdır. Zaman gecikmeli olsa b...
Access restricted. Please log in or start a trial to view this content.
Fiziksel Stokastik Optimal kontrol teorisinin (PSO) tanıtılmasıyla birlikte, bu yazıda, MR damperleri tarafından kontrol edilen yapıların güvenilirliğini korumak için tasarlanmış, zaman gecikmesi kompanzasyonlu bir Stokastik Optimal Yarı Aktif Kontrol yöntemi (SOSC-PSO) önerilmiştir. Önerilen yöntemdeki zaman gecikmesini telafi etmek için, yarı aktif kontrol kuvveti sadece mevcut durumların değil, aynı zamanda ayrık durum uzayında önceki zaman adımlarındaki durumların ve kontrol kuvvetle...
Access restricted. Please log in or start a trial to view this content.
Tüm yazarların beyan edeceği herhangi bir çıkar çatışması yoktur.
Yazarlar, Hebei Eyaleti Doğa Bilimleri Vakfı'nın (Hibe No. E2023210007) desteğini minnetle kabul ederler.
Access restricted. Please log in or start a trial to view this content.
| Name | Company | Catalog Number | Comments |
|---|---|---|---|
| MATLAB | MathWorks | 2016 | Önerilen yöntemin hesaplanması için kullanılması |
Access restricted. Please log in or start a trial to view this content.
Request permission to reuse the text or figures of this JoVE article
Request Permission