Bu çalışma, hava destekli ve su destekli levha benzeri yapıların küresel sualtı geçici basınç darbelerine verdiği tepkiyi kapalı biçimli analitik çözümlerle incelemiştir. Araştırmacılar, yapıyı tek serbestlik dereceli yığılmış kütle ve yay-kütle sistemleriyle temsil etmiş; çözümleri 4,2 gram küresel TNT yükünün, suyla dolu konik bir şok tüpündeki yayılımını modelleyen tam bağlı DYSMAS akışkan-yapı etkileşimi simülasyonlarıyla karşılaştırmıştır. Yay-kütle modeli, su destekli sistemlerde iletilen basıncın tepe değerini bütün doğrulama örneklerinde %5,5’in altında; hava destekli sistemlerde tepe levha hızını %16’nın altında hatayla tahmin etmiştir. Bununla birlikte model, gerçek levhanın uzaysal eğilmesini, yüksek titreşim modlarını, kavitasyonu, akışkan-yapı ayrılmasını ve tam akustik yayılım alanını çözmeyen indirgenmiş mertebeli bir yaklaşımdır.

Araştırmanın temel sonucu, küresel dalga geometrisinin yapısal tepki üzerinde ikincil bir ayrıntı olmadığıdır. Küresel dalgalarda basınç ile akışkan parçacık hızı aynı fazda değildir; “ardıl akış” olarak adlandırılan bileşen, basınç cephesi geçtikten sonra da akışkan hareketini ve yapı-akışkan arasındaki enerji alışverişini sürdürebilir. Çalışma, akışkan yüklemesini temsil eden α, küresel yayılma ve ardıl akışı temsil eden β ile darbe süresini temsil eden γ parametrelerinin birlikte titreşim, sönüm ve rezonans rejimlerini belirlediğini göstermektedir. Düzlemsel dalga kabulü β = 0 sınırına karşılık gelir; β büyüdükçe düzlemsel ve küresel model tahminleri birbirinden belirgin biçimde ayrılmaktadır.

Türkiye açısından değerlendirme: Çalışmanın yaklaşımı; gemi ve denizaltı yapıları, sualtı araçlarının basınç gövdeleri, açık deniz enerji tesisleri, sualtı boru hatları, kapalı akışkan sistemlerindeki basınç darbeleri ve sualtı altyapısının geçici hidrodinamik yüklere karşı ön tasarımı açısından Türkiye’de uyarlanabilir bir analiz çerçevesi sunabilir. Ancak Türkiye’de kullanılacak malzeme, geometri, sınır koşulu, su derinliği ve yükleme senaryoları için modelin yerel deneylerle, daha ayrıntılı sonlu eleman-akışkan modelleriyle ve ölçülmüş basınç verileriyle doğrulanması gerekir. Araştırma Türkiye’de üretilmiş bir deniz platformunu, yerli bir malzemeyi, gerçek bir saha olayını veya Türkiye’ye özgü güvenlik ve maliyet sonuçlarını incelemediğinden ülkede doğrudan performans ya da dayanım sonucu çıkarılamaz.

Araştırmanın temel sorusu nedir?

Sualtı geçici basınç dalgaları; sualtı patlamaları, yapıların içe çökmesi, su darbesi, yapısal titreşimler ve gaz boşluklarının çökmesi gibi farklı olaylardan kaynaklanabilir. Bu dalgalar bir levhaya veya basınç sınırına ulaştığında yapı hareket eder, dalganın bir bölümü yansır ve yapı suyla çevriliyse bir basınç bileşeni karşı tarafa iletilir.

Önceki analitik çalışmaların büyük bölümü dalga cephesini düzlemsel kabul etmiştir. Bu kabul uzak alan veya çok büyük eğrilik yarıçapı için kullanışlı olabilir; ancak kaynağa görece yakın sonlu yarıçaplı dalgalarda basınç cephesi küreseldir. Küresel dalgada basınç ile parçacık hızı arasında faz farkı oluşur ve akışkan içinde geçici olarak depolanan enerji ardıl akış meydana getirir.

Araştırmanın ana sorusu şöyledir: Küresel yayılma ve ardıl akış, hava veya suyla desteklenen yapıların titreşimini, sönümünü, rezonansını, levha hızını ve karşı tarafa iletilen arayüz basıncını nasıl değiştirir?

Çalışmanın literatüre getirdiği yaklaşım

Çalışma dört kapalı biçimli model geliştirmiştir:

• Su destekli yay-kütle sistemi,
• Hava destekli yay-kütle sistemi,
• Su destekli yığılmış kütle sistemi,
• Hava destekli yığılmış kütle sistemi.

Yay-kütle modeli hem yapının ataleti hem de eşdeğer rijitliği içerir. Yığılmış kütle modeli ise rijitliği ihmal ederek yalnızca atalet tepkisini ele alır. Araştırma, bu iki yaklaşımın yalnızca belirli yükleme ve doğal frekans koşullarında birbirine yakın sonuç verdiğini göstermiştir.

Şekil 1’deki hava ve su destekli yapılar neyi temsil etmektedir?

Şekil 1’de küresel bir basınç dalgası, su içinden dairesel bir levhaya doğru ilerlemektedir. Hava destekli durumda levhanın arka yüzünde hava bulunduğu için karşı tarafa iletilen su basıncı ihmal edilir. Su destekli durumda ise levhanın iki tarafında da su vardır ve hareket eden levha ile arka taraftaki su arasında ek bir akışkan-yapı etkileşimi oluşur.

Yapı, birim alan başına kütlesi bulunan rijit bir eleman ve eşdeğer bir yayla temsil edilmiştir:

\[ m = \rho_s h \]

Burada m birim alan başına yapısal kütleyi kg/m², ρ_s yapı malzemesinin yoğunluğunu kg/m³, h ise levha kalınlığını metre cinsinden gösterir.

Yay-kütle sisteminin hareket denklemi şöyledir:

\[ \frac{d^2w}{dt^2}+\mu^2w=\frac{P_{total}}{m} \]

w levhanın merkezindeki tek boyutlu yer değiştirmeyi, μ yapının havadaki doğal frekans parametresini, P_total levhaya etkiyen net basıncı ve m birim alan başına kütleyi temsil eder. Model, gerçek levhanın yüzey boyunca değişen eğilmesini tek bir hareket koordinatına indirger. Eğilme rijitliği, geometri ve kenar koşulları tek bir eşdeğer yay teriminin içine gömülmüştür.

Model hangi fiziksel ayrıntıları içermez?

Tek serbestlik dereceli indirgeme nedeniyle model:

• Levha yüzeyindeki uzaysal deformasyon dağılımını,
• Birden fazla eğilme ve titreşim modunu,
• Levha kalınlığından geçen elastik dalgaları,
• Kenar bölgelerindeki yerel gerilme ve deformasyonları,
• Plastik şekil değiştirme veya hasarı,
• Kavitasyon ve akışkan-yapı ayrılmasını,
• Yapının oluşturduğu tam üç boyutlu akustik yayılım alanını

doğrudan çözmez. Bu nedenle sonuçlar ayrıntılı yerel gerilme veya hasar tahmini değil, baskın zamansal sistem tepkisinin indirgenmiş bir temsilidir.

Küresel dalgada parçacık hızı neden farklıdır?

Küresel bir basınç dalgasındaki akışkan parçacık hızı iki bileşenle tanımlanmıştır:

\[ v=\frac{P}{\rho_f c_f}+\frac{1}{\rho_f R}\int_0^t P\,dt \]

Burada:

• v, akışkan parçacık hızı, m/s;
• P, akustik basınç, Pa;
• ρ_f, akışkan yoğunluğu, kg/m³;
• c_f, akışkandaki ses hızı, m/s;
• R, basınç kaynağı ile yapı arasındaki dalga yarıçapı veya mesafe, m;
• t, zaman, s

olarak tanımlanır.

İlk terim kaynaktan dışarı yayılan akustik enerjiyi temsil eder. İntegral içeren ikinci terim ardıl akıştır. Düzlemsel dalga sınırında R çok büyüdüğünden bu ikinci terim küçülür. Sonlu yarıçaplı küresel dalgada ise ardıl akış, basınç azalmış olsa bile akışkan hareketini sürdürerek yapı ile akışkan arasındaki enerji alışverişini değiştirebilir.

Yansıyan ve iletilen basınç ne anlama gelmektedir?

Su destekli sistemde levhaya etkiyen net basınç:

\[ P_{total}=P_i+P_r-P_t \]

olarak yazılmıştır. P_i gelen, P_r yansıyan ve P_t iletilen basınç bileşenidir.

Çalışmadaki “yansıyan” ve “iletilen” basınçlar, tam uzaysal akustik dalga alanları değildir. Bunlar, levha merkezindeki akışkan-yapı arayüzünde hız ve basınç sürekliliğinden elde edilen indirgenmiş arayüz bileşenleridir. Model, yapının hareketini yeni bir küresel akustik kaynak olarak çözmemektedir.

Arayüzdeki hız sürekliliği su destekli durumda şu şekilde verilmiştir:

\[ \frac{dw}{dt}=v_i-v_r=v_t \]

Bu ilişkilerin sadeleştirilmesiyle çalışmada su destekli yapı için net basınç:

\[ P_{total}=2P_r \]

biçiminde elde edilmiştir.

Gelen basınç darbesi nasıl modellenmiştir?

Kaynakta üretilen küresel basınç dalgası, yapı konumuna ulaştığı anda sabit bir R mesafesinde değerlendirilmiştir:

\[ P_i=\frac{P_{mag}}{R}\exp\left(-\frac{t}{n}\right) \]

P_mag/R terimi arayüzdeki başlangıç basıncı P₀ olarak birleştirildiğinde:

\[ P_i=P_0\exp\left(-\frac{t}{n}\right) \]

elde edilir. P₀ basınç darbesinin yapıdaki başlangıç büyüklüğünü, n ise saniye cinsinden üstel azalma zaman sabitini gösterir. Bu yaklaşım basıncın zamanla üstel azaldığını varsayar.

Üç boyutsuz parametre neyi ifade etmektedir?

Su destekli yay-kütle sistemi için üç temel boyutsuz parametre tanımlanmıştır:

\[ \alpha=\frac{2\rho_f c_f}{\mu m} \]

\[ \beta=\frac{c_f}{\mu R} \]

\[ \gamma=\mu n \]

Hava destekli durumda levhanın yalnızca bir tarafı suyla temas ettiği için α parametresindeki 2 katsayısı bulunmaz:

\[ \alpha_{air}=\frac{\rho_f c_f}{\mu m} \]

Bu fark, aynı fiziksel yapı ve akışkan için su destekli sistemin daha güçlü akışkan yüklemesine maruz kalması anlamına gelir.

Doğal frekans nasıl hesaplanmıştır?

Tam kenetlenmiş dairesel levhanın doğal frekans parametresi klasik levha teorisiyle hesaplanmıştır:

\[ \mu=\frac{\lambda_{mn}^2}{a^2}\sqrt{\frac{D}{m}} \]

\[ D=\frac{Eh^3}{12(1-\nu^2)} \]

Burada λ_mn titreşim moduna bağlı boyutsuz özdeğer, a levha yarıçapı, D eğilme rijitliği, E elastisite modülü, h kalınlık ve ν Poisson oranıdır.

Gösterim uyarısı: Denklem, μ’yü açısal doğal frekans biçiminde vermektedir. Metin daha sonra f = μ/(2π) dönüşümünden söz etmekte, ancak Tablo 1’de sütun “μ (Hz)” olarak adlandırılmakta ve bu değerlerin doğrudan modelde kullanıldığı belirtilmektedir. Hz ile rad/s arasındaki bu gösterim ayrımı, model yeniden uygulanırken açıklığa kavuşturulması gereken bir kaynak tutarsızlığıdır.

Salınımlı ve salınımsız tepki nasıl ayrılmıştır?

Su destekli yay-kütle sisteminin basınç denklemi üçüncü dereceden doğrusal bir diferansiyel denkleme dönüşmektedir. Homojen çözümün karakteristik kökleri α ve β’ye bağlıdır:

• Bir gerçek kök ve iki karmaşık eşlenik kök oluşursa sistem az sönümlü veya salınımlı kabul edilmiştir.
• Üç kök de gerçekse sistem aşırı sönümlü veya salınımsız kabul edilmiştir.

Çalışmanın hareket denkleminde ayrı bir yapısal sönüm terimi bulunmamaktadır. Grafiklerde görülen “sönüm”, akışkan-yapı eşleşmesinin karakteristik kökler üzerindeki etkisinden doğmaktadır. Enerji akışkana aktarılırken salınım genliği azalabilir; bu durum malzeme içi viskoz veya yapısal sönümün modele dâhil edildiği anlamına gelmez.

Şekil 5’te α-β düzleminin yalnızca sınırlı bir bölümü aşırı sönümlü çıkmıştır. Çalışmaya göre aşırı sönümlü tepkinin oluşabildiği en yüksek β değeri yaklaşık 0,19’dur. β büyüdüğünde ardıl akış, basınç darbesinin ardından akışkan hareketini sürdürerek yapısal salınımı destekleyebilir.

Metinde α = 5 ve β = 0,04 sistemi aşırı sönümlü; α = 5 ve β = 0,06 sistemi az sönümlü olarak tanımlanmıştır. Ancak Şekil 5(b)’de salınımlı eğrinin rengi ile lejanttaki “salınımlı/salınımsız” renk tanımları ters görünmektedir. Metin, parametreler ve eğrilerin davranışı birlikte değerlendirildiğinde lejantta renk eşleşmesi hatası bulunması olasıdır.

α ve β titreşim frekansını nasıl değiştirmiştir?

Şekil 6’daki kontur haritalarında az sönümlü bölge için hem titreşim frekansının hem de salınım azalım hızının α ve β arttıkça yükseldiği görülmektedir. α’nın büyümesi, yapıya göre daha güçlü akışkan etkisini ve ek akışkan kütlesini temsil eder. β’nın büyümesi ise ardıl akışı ve yapı çevresindeki sürdürülen akışkan hareketini güçlendirir.

Çalışmada iki örnek verilmiştir:

• α = 0,1 ve β = 1,6 için boyutsuz sönüm oranı yaklaşık 0,013’tür. Tepki uzun süre salınmaya devam etmektedir.
• α = 1,4 ve β = 0,27 için sönüm oranı yaklaşık 1’e yaklaşmaktadır. İletilen basınç yaklaşık bir salınımdan sonra sıfıra yaklaşmaktadır.

İkinci sistem matematiksel olarak karmaşık köklere sahip az sönümlü bölgede bulunsa da tepki biçimi güçlü akışkan sönümü nedeniyle aşırı sönümlü sisteme benzemektedir.

Rezonans çalışmada nasıl tanımlanmıştır?

Rezonans, periyodik ve sürekli bir harmonik yükün klasik kararlı durum rezonansı olarak değil, analitik çözümdeki özel çözüm katsayısının paydasının sıfıra yaklaşması olarak tanımlanmıştır. Temel koşullar şu ilişkilerle ifade edilmiştir:

\[ \alpha=\frac{1+\gamma^2}{\gamma}-\beta(1+\gamma^2) \]

\[ \beta=\frac{1}{\gamma}-\frac{\alpha}{\gamma^2+1} \]

\[ \beta\gamma^3-\gamma^2+(\alpha+\beta)\gamma-1=0 \]

Belirli α, β ve γ birleşimleri, basınç darbesinin zaman ölçeğini indirgenmiş sistemin doğal tepki ölçeğiyle eşleştirerek analitik çözümün çok büyük veya sınırsız büyümesine yol açmaktadır. Ancak bu “sınırsız tepki”, doğrusal tek modlu ideal modelin matematiksel tekilliğidir. Gerçek yapılarda plastikleşme, geometrik doğrusal olmama, hasar, kavitasyon, malzeme sönümü ve sınırlı enerji gibi etkiler büyümeyi sınırlar.

Çalışmadaki örneklerde:

• β = 1 ve γ = 0,1 için α yaklaşık 9,09 olduğunda belirgin bir rezonans tepesi oluşmuştur.
• α = 2 ve γ = 0,1 için β yaklaşık 8,02 olduğunda rezonans elde edilmiştir.
• α = 1 ve β = 1 için γ yaklaşık 0,56’da analitik rezonans görülmüştür.
• α = 1 ve β = 0 olan düzlemsel dalga örneğinde hiçbir pozitif γ değeri rezonans üretmemiştir.

γ denklemi üçüncü dereceden olduğundan aynı α-β birleşimi için birden fazla pozitif matematiksel rezonans kökü oluşabilir. Bunlar farklı basınç azalma sürelerinin aynı indirgenmiş sistemi etkili biçimde uyarmasına karşılık gelir.

Düzlemsel ve küresel dalga modelleri neden ayrışmıştır?

β parametresi dalga yarıçapını içerir:

\[ \beta=\frac{c_f}{\mu R} \]

R büyüdükçe β sıfıra yaklaşır ve küresel model düzlemsel dalga sınırına iner. Düzlemsel dalgada α > 2 olduğunda çözüm her zaman aşırı sönümlüdür. Küresel modelde ise α > 2 olsa bile yeterince büyük β, sistemi salınımlı bölgeye taşıyabilir.

Şekil 9, küresel ve düzlemsel modellerin maksimum iletilen basınç oranını göstermektedir. β küçükken iki yaklaşım birbirine yakındır. β yükseldikçe ardıl akış güçlenir ve iki modelin basınç tahminleri belirgin biçimde ayrılır. Bu sonuç, sonlu yarıçaplı veya yakın alan basınç darbelerinde dalga cephesini otomatik olarak düzlemsel kabul etmenin titreşim ve iletim tahminlerini değiştirebileceğini göstermektedir.

Sayısal doğrulama nasıl kurulmuştur?

Analitik çözümler, Dynamic System Mechanics Advanced Simulation (DYSMAS) yazılımındaki tam bağlı akışkan-yapı etkileşimi modelleriyle karşılaştırılmıştır. DYSMAS:

• Akışkan için Euler tipi Gemini çözücüsünü,
• Yapı için Lagrange tipi ParaDyn çözücüsünü,
• İki alan arasında standart bir bağlayıcı arayüzü

kullanmaktadır.

Gemini, sıkıştırılabilir ve viskozitesiz akışkanlarda şok dalgası yayılımı için yüksek dereceli Godunov tabanlı bir yöntem kullanır. Sayısal model; patlayıcının tepkisini, basınç dalgasının yayılmasını ve dalga cephesi ile hedef levha arasındaki etkileşimi birlikte simüle etmiştir.

Konik şok tüpü modeli

Şekil 2’de Naval Undersea Warfare Center Division Newport’taki konik şok tüpünün sayısal bir kopyası gösterilmektedir. Tüp yatay, suyla dolu ve iç koni açısı 2,6° olan dört bölümlü bir düzeneğe sahiptir:

1. Patlayıcı bölümü,
2. Konik bölüm,
3. Hedef levha,
4. Su veya hava destekli arka bölüm.

Model hesap yükünü azaltmak için düşey ve düzlem dışı yönlerde dörtte bir simetri kullanmıştır. Konik tüpün duvarları rijit, yalnızca hedef levha şekil değiştirebilir olarak modellenmiştir. Böylece çevre yapısının deformasyonu hariç tutularak hedef tepkisi izole edilmiştir.

Gelen basınç darbesi nasıl çıkarılmıştır?

Önce hedef levha bulunmayan bir simülasyon gerçekleştirilmiştir. Hedef konumundaki serbest alan basıncından:

• Başlangıç basıncı 24,1 MPa,
• Üstel azalma süresi 0,5 ms

olarak belirlenmiştir. Bu değerler analitik modeldeki:

\[ P_i=24{,}1\,\text{MPa}\,\exp\left(-\frac{t}{0{,}5\,\text{ms}}\right) \]

biçimindeki giriş darbesini oluşturmak için kullanılmıştır. Şekil 3(a), bu üstel eğrinin sayısal basınç geçmişinin genel düşüşünü temsil ettiğini, ancak simülasyondaki yerel dalgalanmaları içermediğini göstermektedir.

Ağ yakınsaması yeterli miydi?

Akışkan hücresi boyutunun sonucu belirlemediğini kontrol etmek için 0,30; 0,15 ve 0,05 cm hücrelerle üç simülasyon yapılmıştır. Karşılaştırma, 0,3175 cm kalınlığındaki çelik su destekli levha için gerçekleştirilmiştir.

Bu sonuçlara dayanarak 0,15 cm hücre boyutu ana simülasyonlar için yeterli kabul edilmiştir.

Geers hata ölçütleri neyi göstermektedir?

Zamana bağlı iki sinyalin yalnızca tepe değerini değil, bütün dalga biçimini karşılaştırmak için Geers hata ölçütleri kullanılmıştır:

\[ G_M=\sqrt{\frac{\sum c_i^2}{\sum m_i^2}}-1 \]

\[ G_P=1-\frac{\sum m_i c_i}{\sqrt{\sum c_i^2}\sqrt{\sum m_i^2}} \]

\[ G_C=\sqrt{G_M^2+G_P^2} \]

G_M büyüklük farkını, G_P dalga biçimi ve faz farkını, G_C ise ikisinin birleşik değerini gösterir. Sıfıra yakın değerler iki sinyal arasında daha güçlü uyuma karşılık gelir.

Su destekli model doğrulaması

Yedi su destekli poli-karbonat ve çelik levha örneği, geniş bir α-β parametre aralığını kapsayacak şekilde modellenmiştir. Analitik model parametreleri sayısal sonuçlara göre ayarlanmamış; malzeme, geometri ve giriş darbesinden bağımsız olarak hesaplanmıştır.

En düşük tepe basınç hatası, 0,3175 cm kalınlığındaki poli-karbonat levhada %0,17; en yüksek hata ise 0,9525 cm kalınlığındaki çelik levhada %5,32 olmuştur.

Şekil 4(a), en iyi eşleşen poli-karbonat örneğinde analitik ve sayısal basınçların benzer azalma eğilimi gösterdiğini ortaya koymaktadır. Şekil 4(b)’de kalın çelik levha için sayısal eğride daha belirgin salınımlar ve geç zaman farkları oluşsa da analitik model genel büyüklük ve azalma eğilimini korumaktadır.

Hava destekli model doğrulaması

Hava destekli sistemde karşı tarafa iletilen su basıncı yerine levhanın tepe hızı ve son yer değiştirmesi karşılaştırılmıştır.

Üç örnekte faz hatası 0,005 civarında veya altında; birleşik hata yaklaşık 0,08 veya altında kalmıştır. Analitik model tepe büyüklüğünü her durumda aynı kesinlikle yakalamamış olsa da hızın yükselme zamanını, genel büyüklüğünü ve azalma eğilimini nitel olarak temsil etmiştir.

Analitik ve sayısal sonuçlar neden tam olarak örtüşmemiştir?

DYSMAS modeli ile indirgenmiş çözüm arasındaki farkların başlıca nedenleri şunlardır:

• Analitik modelde levha tek bir rijit kütle ve eşdeğer yaydır.
• Uzaysal eğilme ve yüksek titreşim modları bulunmaz.
• Sayısal modelde görülebilen yoğunluk değişimleri ve kavitasyon analitik modelde yoktur.
• Akışkan ile levha arasındaki ayrılma hesaba katılmaz.
• Parçacık hızı doğrusal küresel akustik dalga biçiminde idealize edilmiştir.
• Yansıyan ve iletilen basınç, tam radyasyon alanı yerine levha merkezindeki arayüz değeridir.
• Sayısal model şok dalgasının yerel dalgalanmalarını içerirken analitik giriş tek bir üstel eğridir.

Su destekli yapılarda iletilen basınç ne zaman azalabilir?

Şekil 10, γ = 0,5 için α-β düzleminde yüksek ve düşük iletim bantları göstermektedir. Küçük α ve büyük β değerlerinde yapı parametreleri akışkan yüklemesine baskın gelmekte ve güçlü ardıl akış nedeniyle yansıyan basınç gelen basınca yaklaşabilmektedir. Bu durumda karşı tarafa iletilen arayüz basıncı küçülür.

Örnek olarak α = 0,1, β = 10 ve γ = 0,5 durumunda gelen ve yansıyan basınçlar birbirine yakın, iletilen basınç ise sıfıra yakın hesaplanmıştır. Buna karşılık α = 41,5, β = 2,4 ve γ = 0,5 durumunda yansıyan basınç küçülmüş ve iletilen basınç büyümüştür.

Bu sonuç, gerçek bir koruyucu levhanın basıncı tamamen engellediği anlamına gelmez. Çalışmadaki değerler, indirgenmiş arayüz modelinin belirli boyutsuz parametreler altındaki matematiksel davranışıdır.

Hava ve su destekli yapılar arasındaki fark

Aynı boyutsuz α, β ve γ değerleri kullanıldığında hava ve su destekli sistemlerin titreşim frekansı, fazı ve sönüm rejimi aynı matematiksel yapıyı izler. Hava destekli boyutsuz hız genliği, çalışmadaki karşılaştırmada su destekli değerin iki katı olmuştur.

Ancak gerçek boyutlu sistemlerde aynı yapı ve aynı akışkan için α değerleri eşit değildir. Su destekli yapıda su iki taraftan etkili olduğu için akışkan yüklemesi daha güçlüdür. Bu durum:

• Genliği,
• Titreşim frekansını,
• Fazı,
• Sönümü,
• Rezonans bölgesini,
• Az veya aşırı sönümlü olma durumunu

değiştirebilir. Bu nedenle gerçek hava ve su destekli yapılar yalnızca eşit boyutsuz parametrelerle değil, kendi boyutlu fiziksel değerleriyle karşılaştırılmalıdır.

Yığılmış kütle modeli ne zaman yetersiz kalmıştır?

Yığılmış kütle yaklaşımı yapısal rijitliği sıfır kabul eder:

\[ \frac{d^2w}{dt^2}=\frac{P_{total}}{m} \]

Yay-kütle ve yığılmış kütle modelleri, doğal frekans çok küçük olduğunda benzer sonuç vermiştir. Bu durumda yükleme süresi içinde belirgin bir geri çağırıcı yay etkisi gelişmez ve tepki atalet ağırlıklı olur.

Doğal frekans arttıkça yay-kütle modeli daha düşük iletilen basınçlar üretmiş ve iki çözüm ayrışmıştır. Çalışmanın önemli sonuçlarından biri, yalnızca birim alan başına kütleyi artırmanın yığılmış kütle kabulünü geçerli kılmamasıdır. Kütle sonsuza yaklaşsa bile yay-kütle sisteminde doğal frekanstan kaynaklanan geri çağırıcı terim bulunmaya devam eder; yığılmış kütlede ise böyle bir mekanizma yoktur.

Çalışmanın güçlü yönleri nelerdir?

• Küresel dalga ve ardıl akış etkisi kapalı biçimli analitik çözümlere dâhil edilmiştir.
• Hava ve su destekli dört farklı yapısal idealizasyon aynı çerçevede ele alınmıştır.
• Model parametreleri sayısal sonuçlara göre sonradan uydurulmadan karşılaştırılmıştır.
• Yedi su destekli ve üç hava destekli örnek geniş bir α-β alanını kapsamıştır.
• 90 milyondan fazla hücreli tam bağlı sayısal model kullanılmıştır.
• Ağ yakınsaması tepe değer, darbe ve bütün zaman geçmişi üzerinden sınanmıştır.
• Hata değerlendirmesi yalnızca tepe değerlere değil, Geers büyüklük ve faz ölçütlerine de dayandırılmıştır.
• Düzlemsel/küresel ve yığılmış kütle/yay-kütle varsayımlarının geçerlilik bölgeleri karşılaştırılmıştır.

Çalışmanın sınırlılıkları nelerdir?

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Çalışmada yeni bir fiziksel deney gerçekleştirilmemiştir.
• Doğrulama, başka çalışmalarda deneysel olarak sınanmış bir sayısal yazılıma karşı yapılmıştır; analitik çözüm doğrudan ölçüm verisine karşı sınanmamıştır.
• Levha tek serbestlik derecesine indirgenmiştir.
• Plastik deformasyon, hasar ve kırılma modellenmemiştir.
• Kavitasyon ve akışkan-yapı ayrılması analitik çözümde yoktur.
• Yükleme tek üstel basınç darbesiyle temsil edilmiştir.
• Tam üç boyutlu yansıyan ve yayılan akustik alan çözülmemiştir.
• Rezonanstaki sınırsız büyüme fiziksel değil, doğrusal model tekilliğidir.
• Doğal frekansın Hz veya rad/s olarak kullanımı kaynakta tam açık değildir.
• Şekil 5(b)’de lejant ile eğri davranışları arasında olası renk eşleşmesi sorunu vardır.
• Aşırı sönümlü bölge, konik tüpün sabit mesafesinin erişilebilir β aralığını sınırlaması nedeniyle sayısal doğrulamada incelenmemiştir.
• Model yalnızca dairesel ve tam kenetlenmiş hedeflerin eşdeğer tek modlu tepkisini temsil etmektedir.

Çalışma neyi desteklemektedir?

• Küresel dalga eğriliği ve ardıl akış, yapısal titreşim ve basınç iletimi tahminlerini değiştirebilir.
• α ve β, az sönümlü veya aşırı sönümlü tepki rejimini belirleyen temel parametrelerdir.
• γ, darbe süresi ile yapısal tepki zamanını karşılaştırarak rezonans koşullarında önemli rol oynar.
• Yay-kütle modeli, incelenen sayısal örneklerde baskın basınç ve hız eğilimlerini kabul edilebilir hatayla temsil edebilir.
• Düzlemsel dalga yaklaşımı yalnızca β küçükken küresel modelle yakın sonuç verir.
• Yüksek kütle tek başına yapısal rijitliği ihmal etmek için yeterli değildir.

Çalışma neyi kanıtlamamaktadır?

• Gerçek bir gemi, denizaltı veya sualtı aracının hasar dayanımını kanıtlamamaktadır.
• Belirli bir levhanın sualtı patlamasına karşı güvenli olduğunu göstermemektedir.
• Malzeme akması, çatlama, kopma veya kalıcı deformasyon öngörüsü sunmamaktadır.
• Rezonans koşullarında gerçek basıncın fiziksel olarak sonsuza çıkacağını göstermemektedir.
• İletilen arayüz basıncını uzak alandaki tam akustik basınç alanı olarak vermemektedir.
• Bütün dalga biçimlerinin üstel azalım modeliyle temsil edilebileceğini kanıtlamamaktadır.
• Türkiye’deki deniz yapıları için doğrulanmış bir tasarım standardı oluşturmamaktadır.

Gündelik ve mühendislik uygulamaları açısından anlamı

Çalışmanın esas değeri doğrudan bir nihai yapı tasarlamasından çok, ayrıntılı simülasyon veya deney öncesinde baskın parametrelerin hızlı biçimde taranabilmesini sağlamasıdır. Bir mühendisin malzeme yoğunluğunu, kalınlığı, doğal frekansı, akışkan özelliklerini, dalga yarıçapını ve darbe süresini değiştirerek sistemin hangi titreşim rejimine yaklaşacağını önceden değerlendirmesine yardımcı olabilir.

Yaklaşım; sualtı basınç gövdeleri, kapalı su sistemlerindeki geçici basınç sınırları, gemi kaplamaları, açık deniz platformları, sualtı boruları ve şok tüpü deneylerinin ön tasarımında azaltılmış mertebeli bir tarama aracı olarak incelenebilir. Gerçek tasarım kararları için malzeme doğrusal olmaması, hasar, kaynak ve bağlantılar, üç boyutlu geometri, farklı basınç biçimleri ve deneysel doğrulama ayrıca gereklidir.

Çalışmanın Yöntemi ve Bulguları

Yöntemsel tasarım

Doğrulama için kullanılan malzemeler ve yapısal aralıklar

Su destekli örneklerden seçilmiş sayısal sonuçlar

Ana bulgular

• Yay-kütle analitik çözümü, incelenen su destekli örneklerde tepe iletilen basıncı %5,5’in altında hatayla tahmin etmiştir.
• Su destekli örneklerde toplam darbe hatası %8’in altında, ortalama darbe hatası %6,63 olmuştur.
• Hava destekli örneklerde tepe levha hızı hatası %1,79 ile %15,35 arasında kalmıştır.
• Hava destekli örneklerde son yer değiştirme hatası %1,62 ile %9,21 arasında olmuştur.
• 0,15 cm ağ ile 0,05 cm ağ arasındaki tepe basınç farkı yalnızca %1,1’dir.
• β yaklaşık sıfıra yaklaşırken küresel ve düzlemsel model sonuçları birbirine yaklaşmıştır.
• β arttığında ardıl akış, düzlemsel modelin öngörmediği salınımlı tepkileri sürdürebilmiştir.
• Aşırı sönümlü bölgenin oluşabilmesi için β’nın yaklaşık 0,19’dan küçük olması gerekmiştir.
• Belirli α-β-γ birleşimleri analitik çözümde rezonans tekilliği oluşturmuştur.
• Aynı boyutsuz parametrelerde hava destekli hız genliği, su destekli genliğin iki katı hesaplanmıştır.
• Gerçek boyutlu karşılaştırmada su destekli yapı daha güçlü akışkan yüklemesi, farklı frekans ve farklı sönüm göstermiştir.
• Doğal frekans yükseldiğinde yay-kütle ve yığılmış kütle modelleri belirgin biçimde ayrışmıştır.
• Birim alan kütlesinin büyümesi tek başına rijitliğin ihmal edilmesini haklı çıkarmamıştır.

Grafiklerin birlikte gösterdiği mekanizma

Şekil 3, giriş basıncının üstel yaklaşımını ve ağ yakınsamasını; Şekil 4, analitik-sayısal zaman geçmişi uyumunu; Şekil 5 ve 6, α-β düzlemindeki titreşim ile sönüm rejimlerini göstermektedir. Şekil 7 ve 8, rezonansın yalnızca tek bir parametreye değil, α, β ve γ’nın birlikte oluşturduğu dar bölgelere bağlı olduğunu ortaya koymaktadır.

Şekil 9, düzlemsel ve küresel dalga tahminlerinin β arttıkça ayrıştığını göstermektedir. Şekil 10, su destekli sistemlerde düşük ve yüksek iletim bantlarını; Şekil 11 ve 12 ise rezonans ile sönümün birlikte değerlendirilmesi gerektiğini göstermektedir. Şekil 13, hava ve su destekli koşullar arasındaki boyutsuz benzerliği fakat boyutlu fiziksel farkı; Şekil 14 ise yapısal rijitliğin ihmal edildiği yığılmış kütle yaklaşımının geçerlilik sınırını açıklamaktadır.

Sonuçların yorum sınırı

Modelin düşük hata değerleri, incelenen DYSMAS örneklerinin baskın zamansal davranışının tek serbestlik dereceli bir analitik sistemle temsil edilebildiğini göstermektedir. Bu sonuç, modelin gerçek yapılardaki yerel gerilmeyi, hasarı, kavitasyonu veya bütün akustik alanı aynı doğrulukla tahmin ettiğini göstermez. Analitik çözüm, ayrıntılı simülasyon ve deneylerin yerini alan nihai bir sertifikasyon aracı değil, parametre duyarlılığı ve ön değerlendirme yöntemidir.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı: Response of Structures Subjected to Spherical Underwater Transient Pressure Pulses

Yazarlar: Carlos Javier, Shyamal Kishore, Michael Galuska, Michael Papa, James LeBlanc, Helio Matos ve Arun Shukla.

Yazar sıralaması: Kaynakta verildiği biçimde korunmuştur.

Eş birinci yazar: Eşit katkı veya eş birinci yazarlık bilgisi yer almamaktadır.

Sorumlu yazar: Carlos Javier. Çalışmada verilen iletişim adresi carlos.r.javier.civ@us.navy.mil şeklindedir.

Kurum 1: Naval Undersea Warfare Center, Division Newport, 1176 Howell Street, Newport, Rhode Island 02841, ABD.

Kurum 2: Dynamic Photo Mechanics Laboratory, Department of Mechanical, Industrial and Systems Engineering, University of Rhode Island, Kingston, Rhode Island 02881, ABD.

DOI:10.2139/ssrn.7197531

Resmî kaynak bağlantısı:SSRN güncel kayıt sayfası

Önceki SSRN sürümü: Aynı başlık için SSRN 6571436 ve DOI 10.2139/ssrn.6571436 numaralı daha eski bir kayıt da bulunmaktadır. Bu içerikte yüklenen ve 7197531 numarasını taşıyan güncel sürüm esas alınmıştır.

Yayın platformu: SSRN.

Dergi: Belirli bir hakemli dergi adı veya dergi kabulü bu sürümde yer almamaktadır.

Özgün yayınevi: Hakemli dergi yayınevi bulunmamaktadır; çalışma SSRN preprint platformunda sunulmuştur.

Yayın yılı: 2026.

Kaynak türü: Kapalı biçimli analitik çözümler ve tam bağlı akışkan-yapı etkileşimi simülasyonları içeren modelleme çalışması.

Hakemlik durumu: Bu çalışma hakem değerlendirmesinden geçmemiş bir preprinttir; sonuçları bu sınırlılık dikkate alınarak okunmalıdır.

Finansman: Araştırma, National Institute for Undersea Vehicle Technology tarafından NIUVT Cooperative Research and Development Agreement NCRADA-NUWCDIVNPT-23-2214 kapsamında finanse edilmiştir.

Teşekkür: Araştırmacılar, sualtı basınç dalgası yayılımı konusundaki katkısı için Robert Koch’a ve çalışmanın konusunun geliştirilmesine katkısı için Elizabeth Magliula’ya teşekkür etmektedir.

Yazar katkıları: Çalışmada CRediT veya görev bazlı ayrı yazar katkı beyanı yer almamaktadır.

Çıkar çatışması: Ayrı bir çıkar çatışması veya rekabet eden menfaat beyanı bu sürümde yer almamaktadır.

Veri ve kod erişimi: Açık veri deposu, model kodu veya tekrar üretim paketi bağlantısı verilmemiştir.

Yöntemsel sınır: Analitik çözüm tek serbestlik dereceli indirgenmiş bir modeldir. Uzaysal deformasyon, yüksek modlar, plastik hasar, kavitasyon, akışkan-yapı ayrılması ve tam akustik radyasyon alanı çözülmemiştir. Doğrulama yeni fiziksel deneylere değil, DYSMAS simülasyonlarına dayanmaktadır.

Kaynak tutarlılığı notu: Doğal frekans parametresinin Hz ile rad/s gösterimi kaynakta tam tutarlı değildir. Şekil 5(b)’deki salınımlı ve salınımsız eğrilerin renkleri de lejantla ters eşleşiyor görünmektedir. Bu noktalar kaynak metni sessizce değiştirilmeden belirtilmiştir.

Bu Türkçe içerik, yüklenen çalışmanın bilimsel metni, denklemleri, tabloları, şekilleri, ekleri ve sonuçları temel alınarak hazırlanmıştır. Çalışmada bulunmayan bir deney sonucu, saha başarısı, hasar dayanımı, askerî performans, ticari uygulanabilirlik veya güvenlik garantisi eklenmemiştir.