Bu çalışma, buz sınıfı bir gemi pervanesinin buz, su ve pervane arasındaki kısa süreli fakat yüksek maliyetli sayısal etkileşiminden elde edilen gerilme kayıtlarını uzun süreli yorulma çatlağı ilerleme analizinde kullanabilmek için yeni bir gerilme spektrumu yeniden yapılandırma yöntemi geliştirmiştir. Araştırmacılar; LS-DYNA içindeki yapılandırılmış keyfî Lagrange–Euler yöntemi, kohezyonlu buz kırılması modeli, yağmur akışı çevrim sayımı, Paris yasası ve FRANC3D çatlak ilerleme yazılımını tek bir hesaplama zincirinde birleştirmiştir. Yeniden yapılandırılmış spektrumla 51 ilerleme adımı sonunda hesaplanan çatlak uzunluğu 79,15 mm, özgün kısa spektrumun sonsuz tekrarıyla hesaplanan uzunluk ise 77,18 mm olmuş ve aradaki fark yüzde 2,55 olarak bildirilmiştir. Temel sınırlılık, sonuçların tek bir sayısal çalışma durumuna dayanması ve gerçek pervane, deniz buzu veya kutup seferi verileriyle doğrudan doğrulanmamış olmasıdır.

Sayısal modelde 4,115 metre çapında, dört kanatlı, nikel-alüminyum bronz R sınıfı bir pervane; üç metre uzunluğunda bir buz bloğu, su ve hava alanıyla tamamen eşleştirilmiştir. Kararlı olarak tanımlanan 1–2 saniyelik kesitte pervane itmesi yaklaşık −150 kN ile +230 kN arasında yüksek frekanslı ve darbeli biçimde değişmiştir. Çatlak analizi için seçilen 1,16 saniyelik anda kanat kökündeki en yüksek von Mises gerilmesi 254 MPa, en yüksek asal gerilme ise 322 MPa olarak hesaplanmıştır. Çekme yönündeki asal gerilmenin kanat kökünün arka yüzünde yoğunlaşması, başlangıç çatlağının bu bölgeye yerleştirilmesinin gerekçesi olmuştur.

Yeniden yapılandırma yöntemi kısa gerilme kaydını nokta nokta rastgele çoğaltmamaktadır. Özgün kayıttan çıkarılan gerilme blokları; Paris yasasına göre birim zamanda oluşturdukları çatlak sürücüsü, yağmur akışı gerilme aralığı dağılımı, blok birleşimindeki süreklilik ve yakın zamanda kullanılan blokların tekrarından kaçınma ölçütleriyle seçilmektedir. Bu yöntemle 25, 100, 200, 500 ve 1.000 saniyelik spektrumlar üretilmiştir. Bunların süreye göre normalize edilmiş Paris yasası sürücü oranları sırasıyla 0,982, 0,983, 1,002, 1,005 ve 1,020’dir. Çatlak ilerleme hesabında, özgün kayda en yakın oranlardan birini sağlayan 200 saniyelik spektrum kullanılmıştır.

Başlangıç çatlak uzunluğu, incelenen değişkenler arasında en güçlü etkiyi göstermiştir. Çatlak uzunluğu 2 mm’den 10 mm’ye çıkarıldığında, çalışmada tanımlanan referans kırılma sürücü düzeyine ulaşmak için gereken çevrim sayısı 4,05 × 10⁸ çevrimden 0,38 × 10⁸ çevrime düşmüş; azalma yüzde 90,7 olmuştur. Başlangıç çatlak şekil oranı a/c = 0,5’ten 1,0’a yükseltildiğinde ise kendi karşılaştırma grubundaki referans düzeye ulaşmak için gereken çevrim sayısı yüzde 42,9 azalmıştır. Bu değerler nihai kopma ömrü değil, önceden belirlenen kırılma sürücüsü düzeylerine erişme süreleridir.

Türkiye açısından: Yöntem; Türkiye’deki gemi inşa, deniz mühendisliği, pervane üretimi, klaslama araştırmaları ve kutup koşullarına uygun gemi ekipmanı geliştirme çalışmalarında kısa süreli yüksek doğruluklu simülasyonlarla uzun süreli çatlak değerlendirmesi arasında bağlantı kurmak için incelenebilir. Türkiye’de uygulanmadan önce yerel üretim yöntemleriyle oluşturulan pervane malzemelerinin Paris yasası katsayıları, döküm kusurları, artık gerilmeler ve gerçek çatlak dağılımları deneysel olarak belirlenmelidir. Model farklı buz kalınlıkları, ilerleme hızları, pervane devirleri, kavite koşulları ve görev profilleriyle yeniden çalıştırılmalı; sonuçlar tam ölçekli veya model ölçekli ölçümlerle doğrulanmalıdır. Çalışmadan Türkiye’deki belirli bir pervanenin güvenli kullanım süresi, bakım aralığı veya izin verilebilir çatlak boyutu doğrudan çıkarılamaz.

Araştırmanın çözmeye çalıştığı temel problem nedir?

Buz sınıfı gemilerin pervaneleri, normal hidrodinamik kuvvetlerin yanında kanatların buzu kesmesi, ezmesi ve parçaları su içinde hızlandırması nedeniyle kısa süreli ve yüksek genlikli yüklerle karşılaşabilir. Buzla temas sürekli ve düzgün değildir. Bir kanat buz bloğuna yaklaşır, temasa girer, buzu kırar, parçaları ayırır ve yeniden sudaki hidrodinamik çalışma durumuna döner. Bu olayların tekrarlanması, pervane kanat kökünde değişken genlikli gerilme çevrimleri oluşturur.

Bu etkileşimin yüksek doğruluklu sayısal çözümü önemli hesaplama kaynağı gerektirir. Saniyeler düzeyindeki bir buz-su-pervane simülasyonu bile saatlerce sürebilir. Buna karşılık yorulma çatlağının ilerlemesi milyonlarca veya yüz milyonlarca yük çevriminin değerlendirilmesini gerektirir. Çalışmanın hedeflediği boşluk, kısa süreli sayısal gerilme kaydıyla uzun süreli kırılma mekaniği hesabı arasındaki zaman ölçeği farkıdır.

Araştırmacılar gelecekteki her gerilme noktasını tahmin etmeye çalışmamaktadır. Amaç, kısa kaydın yorulma bakımından belirleyici özelliklerini koruyan daha uzun bir yapay gerilme geçmişi üretmektir. Korunması hedeflenen özellikler yalnızca ortalama ve standart sapma değildir; yüksek gerilme aralıkları, yağmur akışı çevrim dağılımı, baskın frekanslar ve Paris yasasına göre çatlak ilerleme potansiyeli de yönteme dâhil edilmiştir.

Çalışmanın bütünleşik hesaplama çerçevesi nasıldır?

Dördüncü sayfadaki Şekil 2, yöntemi beş aşamada göstermektedir:

1. Buz, su ve pervanenin tamamen eşleştirilmiş geçici sayısal analizi
2. Kısa süreli gerilme kaydının uzun süreli yorulma eşdeğer spektruma dönüştürülmesi
3. Kritik kanat kökü bölgesine yarı eliptik başlangıç çatlağının yerleştirilmesi
4. Paris yasası ve FRANC3D ile çatlak ilerlemesinin hesaplanması
5. Çatlak riski, bakım ve tasarım açısından mühendislik değerlendirmesi

LS-DYNA analizi buz kırılmasını ve pervanenin geçici yapısal gerilmesini üretmektedir. Yeniden yapılandırma algoritması, seçilen kanat kökü elemanındaki kısa gerilme geçmişini uzatmaktadır. Ansys Workbench ve FRANC3D ise 1,16 saniyedeki üç boyutlu gerilme alanı üzerine yarı eliptik yüzey çatlağını yerleştirerek her gerilme çevrimi altında çatlak cephesinin nasıl ilerlediğini hesaplamaktadır.

S-ALE yöntemi neden kullanılmıştır?

Yapılandırılmış keyfî Lagrange–Euler, yani S-ALE yöntemi; akışkan alanı ile büyük deformasyon gösteren katı yapıların aynı hesaplama içinde etkileşmesini sağlamak için kullanılmıştır. Pervane ve buz yapısal ağla, hava ve su ise referans akışkan ağıyla temsil edilmiştir. Akışkan ve katı arasındaki momentum aktarımı ceza kuvveti tabanlı eşleştirme yöntemiyle gerçekleştirilmiştir.

Çalışmada kullanılan kütle, momentum ve enerji korunumu denklemleri özetle şöyledir:

\[ \frac{\partial \rho^r}{\partial t} = -\rho\frac{\partial u_i}{\partial x_i} -w_i\frac{\partial\rho}{\partial x_i} \]

\[ \rho\frac{\partial u_i^r}{\partial t} = \sigma_{ij,j} +\rho b_i -\rho w_j\frac{\partial u_i}{\partial x_j} \]

\[ \rho\frac{\partial e^r}{\partial t} = \sigma_{ij}u_{i,j} +\rho b_i u_i -\rho w_j\frac{\partial e}{\partial x_j} \]

Burada \(\rho\) yoğunluğu, \(u_i\) hız bileşenlerini, \(e\) özgül iç enerjiyi, \(w_i\) malzeme ile referans ağ arasındaki göreli hızı ve \(\sigma_{ij}\) gerilme tensörünü temsil etmektedir. Yöntem, pervane dönerken suyun hareketini ve buz parçalarının akışkan içinde taşınmasını aynı çözüm içinde hesaplamayı amaçlamaktadır.

Buz kırılması nasıl modellenmiştir?

Buz bloğu, elastoplastik katı elemanlar ile elemanlar arasına yerleştirilen kohezyonlu kırılma elemanlarından oluşturulmuştur. Kohezyonlu elemanlar başlangıçta komşu buz elemanlarını birbirine bağlamaktadır. Normal veya teğetsel çekme belirlenen dayanımı aştığında hasar başlamakta, ayrılma arttıkça taşıma kapasitesi azalmaktadır. Kritik ayrılma enerjisine ulaşıldığında bağ tamamen kopmakta ve buz parçaları birbirinden ayrılmaktadır.

Beşinci sayfadaki Şekil 4, iki buz elemanı arasındaki kohezyon bölgesinin çatlak oluşturarak ayrılmasını göstermektedir. Bu yaklaşımda çatlak yolu önceden tek bir çizgi olarak tanımlanmamış; uygun eleman arayüzlerinde hasarın gelişmesine izin verilmiştir.

Kohezyonlu elemanlar için yoğunluk 900 kg/m³, en yüksek çekme 4 MPa, enerji salım oranı 433 J/m² ve rijitlik 50.000 MPa olarak kullanılmıştır. Buz gövdesinin elastisite modülü 264,8 MPa, akma gerilmesi 2,61 MPa ve Poisson oranı 0,30’dur. Bu parametreler bu çalışmada yeni deneylerle elde edilmemiş, araştırmacıların önceki modelleme çalışmalarına dayandırılmıştır.

Pervane ve hesaplama alanının özellikleri nelerdir?

Su giriş hızı −X yönünde 3,62 m/s, buzun ilerleme hızı +X yönünde 0,7 m/s ve pervane dönüş hızı 2,698 devir/saniye olarak belirlenmiştir. Akışkan alanında 0,19 metrelik yapılandırılmış elemanlar kullanılmıştır. Pervanenin X yönündeki ötelemesi sınırlandırılmış, diğer serbestlik derecelerine izin verilmiştir.

Ağ duyarlılık analizi ne göstermiştir?

Pervane için 40, 50, 60 ve 70 mm olmak üzere dört ağ büyüklüğü karşılaştırılmıştır. İlk buz kesme aşamasındaki 0,5–1,25 saniye aralığında hesaplanan ortalama kuvvetler şöyledir:

40 ve 50 mm ağların sonuçları birbirine yakınken daha kaba ağlarda belirgin ayrışma görülmüştür. Araştırmacılar hesaplama süresi ile çözüm kararlılığı arasında denge kurmak için 50 mm ağı seçmiştir. Bu karşılaştırma sayısal ağ yakınsamasını desteklemektedir; gerçek ölçüme karşı mutlak doğruluk kanıtı değildir.

Buz kesme sırasında pervane yükleri nasıl değişmiştir?

Sekizinci sayfadaki Şekil 6, pervane itmesinin düzenli bir sinüs dalgası olmadığını göstermektedir. Bir saniyelik kararlı kesitte değerler yaklaşık −150 kN ile +230 kN arasında değişmiştir. Negatif tepeler, buz kesme direncinin bazı anlarda pervanenin pozitif hidrodinamik itmesinden daha büyük hâle geldiğini göstermektedir.

Aynı sayfadaki Şekil 7’de buz, kanat temas bölgesinde çok sayıda küçük parçaya ayrılmaktadır. Parçalar kanat geçişi çevresinde kümelenmekte, ardından su akışı ve pervane hareketiyle taşınmaktadır. Görsel sonuçlar yükün sürekli temas yerine kesme, ezilme, kopma ve parça ayrılması evrelerinden oluştuğunu göstermektedir.

Çatlak için kritik bölge nasıl seçilmiştir?

Kanat kökündeki elemanlar arasında en ağır en yüksek asal gerilme geçmişine sahip 505354 numaralı eleman seçilmiştir. Gerilme kaydı 0,482–2,730 saniye aralığından alınmıştır. Kayıttaki en yüksek çekme durumu 1,16 saniyede oluştuğu için bu andaki üç boyutlu gerilme alanı daha sonraki FRANC3D analizinin temel alanı olarak kullanılmıştır.

Dokuzuncu sayfadaki Şekil 8 ve Şekil 9’a göre:

• En yüksek von Mises gerilmesi 254 MPa’dır.
• Pervane malzemesinin bildirilen akma gerilmesi 268,3 MPa’dır.
• En yüksek asal gerilme 322 MPa’dır.
• Çekme gerilmesi kanat kökünün arka yüzünde yoğunlaşmaktadır.

Von Mises değeri bildirilen akma sınırının yaklaşık yüzde 94,7’sine karşılık gelmektedir. En yüksek asal gerilme ise doğrudan akma ölçütü olarak değil, Mode I açılma çatlağının yönünü ve çekme sürücüsünü belirlemek için kullanılmıştır.

Neden doğrudan kısa gerilme kaydı tekrar edilmemiştir?

Kısa bir kaydın aynen art arda tekrarlanması aynı olay dizisini düzenli aralıklarla yeniden üretir. Böyle bir işlem özgün kayıtta bulunmayan yapay periyodiklik oluşturabilir. Gerilme noktalarının bağımsız biçimde rastgele seçilmesi ise zaman sırasını, frekans yapısını ve çevrimlerin ortalama gerilmelerini bozabilir.

Önerilen yöntem bu iki uç yaklaşım arasında yer almaktadır. Gerilme geçmişinden kayan pencereyle kısa bloklar çıkarılmakta ve blok içindeki zaman sırası korunmaktadır. Yeni uzun spektrum, özgün kaydı başlangıç çapası olarak kullanmakta ve seçilen blokları doğrusal örtüşmeyle birbirine eklemektedir.

Yeniden yapılandırma maliyet fonksiyonu nasıl tanımlanmıştır?

Her aday gerilme bloğu için toplam maliyet şu biçimde hesaplanmıştır:

\[ C_0= \alpha_d C_d+ \alpha_h C_h+ \alpha_b C_b+ C_r+ \varepsilon \]

Ana ağırlıklar \(\alpha_d=1{,}8\), \(\alpha_h=0{,}9\) ve \(\alpha_b=0{,}12\) olarak belirlenmiştir. Her hedef süre için birden fazla rastgele başlangıçla tam spektrum oluşturulmuş ve özgün kayda en iyi uyumu sağlayan sonuç seçilmiştir.

Paris yasasına göre eşdeğer çatlak sürücüsü nasıl hesaplanmıştır?

Paris yasası, kararlı yorulma çatlağı ilerleme hızını gerilme şiddeti faktörü aralığıyla ilişkilendirir:

\[ \frac{da}{dN}=C(\Delta K)^m \]

Gerilme şiddeti faktörü aralığı yaklaşık olarak:

\[ \Delta K=Y\,\Delta\sigma\sqrt{\pi a} \]

biçiminde yazılmıştır. Aynı malzeme, aynı çatlak konumu, aynı geometri ve aynı referans çatlak durumu karşılaştırıldığında \(C(Y\sqrt{\pi a})^m\) ortak çarpandır. Bu nedenle göreli çatlak sürücüsü, yağmur akışıyla bulunan gerilme aralıklarının \(m\) kuvveti üzerinden karşılaştırılabilir:

\[ G= \frac{1}{T} \sum_{i=1}^{N_r} n_i(\Delta\sigma_i)^m \]

• \(G\), birim zaman başına Paris tipi çatlak sürücüsüdür.
• \(T\), gerilme geçmişinin süresidir.
• \(\Delta\sigma_i\), i numaralı yağmur akışı gerilme aralığıdır.
• \(n_i\), bu aralığın çevrim sayısıdır.
• \(m\), Paris yasası üssüdür.

Paris sürücüsü maliyeti:

\[ C_d= \left| \ln\left(\frac{G_{\text{next}}}{G_{\text{raw}}}\right) \right| \]

olarak tanımlanmıştır. Oranın bire yaklaşması, aday blok eklendiğinde birim zaman başına çatlak ilerleme potansiyelinin özgün kayıtla benzer kaldığını göstermektedir.

Yağmur akışı dağılımı neden ayrıca korunmuştur?

İki gerilme geçmişi aynı toplam Paris sürücüsüne sahip olsa bile çevrim dağılımları farklı olabilir. Bir geçmiş çok sayıda orta aralıklı çevrim, diğeri az sayıda çok yüksek aralıklı çevrim içerebilir. Paris yasasındaki yüksek üs nedeniyle büyük aralıklı çevrimlerin etkisi özellikle yüksektir.

Yağmur akışı gerilme aralıkları \(B\) sınıfa ayrılmış ve olasılık dağılımına dönüştürülmüştür. Güncellenen dağılım:

\[ P_{\text{next}} = \beta P_{\text{run}} + (1-\beta)P_{\text{cand}} \]

eşitliğiyle hesaplanmış ve \(\beta=0{,}92\) kullanılmıştır. Özgün ve güncellenmiş dağılım arasındaki fark toplam değişim uzaklığıyla ölçülmüştür:

\[ C_h= \frac{1}{2} \sum_{b=1}^{B} \left| P_{\text{next},b} - P_{\text{raw},b} \right| \]

On birinci sayfadaki Şekil 10, çalışan dağılımla aday bloğun dağılımının nasıl birleştirildiğini şematik olarak göstermektedir.

Blok birleşimindeki yapay çevrimler nasıl sınırlandırılmıştır?

Bir bloğun son gerilmesi ile yeni bloğun ilk gerilmesi arasında büyük sıçrama oluşursa yağmur akışı sayımı özgün kayıtta bulunmayan yüksek aralıklı bir çevrim oluşturabilir. Bu nedenle sınır maliyeti:

\[ C_b= \frac{ |\sigma_{\text{tail}}-\sigma_{\text{cand}}| }{ \operatorname{std}(\sigma_{\text{raw}}) } \]

şeklinde tanımlanmıştır. Bloklar arasında doğrusal örtüşme de uygulanmıştır.

Son beş seçim içinde kullanılan adaylara \(0{,}35\) tekrar cezası verilmiştir. Rastgele pertürbasyonun standart sapması \(0{,}015\)’tir. Bu iki bileşen, algoritmanın sürekli aynı birkaç bloğu seçmesini ve tek bir yerel çözüme kilitlenmesini azaltmayı amaçlamaktadır.

Hiperparametre çıkarım analizi ne göstermiştir?

On ikinci sayfadaki Şekil 11’de ana ağırlıklar, dağılım hafıza katsayısı, rastgele terim ve tekrar cezası değiştirilerek sonuç ölçütlerinin nominal durumdan yüzde farkı gösterilmiştir.

En duyarlı değişken Paris sürücüsü ağırlığı \(\alpha_d\) olmuştur. Bu terim tamamen çıkarıldığında:

• Genel Paris sürücüsü hatası yüzde 12,0 artmıştır.
• Yüzde 95 kuyruk sürücüsü hatası yüzde 12,3 artmıştır.
• Yüzde 99 kuyruk ölçütü yüzde 20,3 azalmıştır.

Son sonuç, tek bir ölçütün iyileşmesinin bütün spektrumun daha iyi olduğu anlamına gelmediğini göstermektedir. Yüksek kuyruk uyumu iyileşirken genel sürücü uyumu bozulabilmektedir.

Histogram ağırlığı çıkarıldığında yüzde 99 kuyruk ölçütü yüzde 17,8 kötüleşmiştir. Sınır sürekliliği ağırlığı çıkarıldığında genel Paris sürücüsü ölçütleri yaklaşık yüzde 10–11 kötüleşmiştir. Araştırmacılar nominal ayarların hiçbir ölçütte tek başına en iyi sonucu vermediğini, ancak ölçütler arasında dengeli bir çözüm oluşturduğunu değerlendirmiştir.

Yeniden yapılandırılmış spektrumlar özgün kayda ne kadar benzemektedir?

Yöntem 25, 100, 200, 500 ve 1.000 saniyelik beş spektrum üretmiştir. Yağmur akışı histogramları, çevrimlerin büyük bölümünün yaklaşık 50 MPa altındaki gerilme aralıklarında toplandığını ve yüksek aralıklarda uzun fakat seyrek bir kuyruk bulunduğunu göstermiştir.

On üçüncü sayfadaki aralık–ortalama matrislerinde yoğunluk merkezi düşük gerilme aralıkları ile yaklaşık 100–150 MPa ortalama gerilme çevresinde kalmıştır. Süre uzadıkça toplam çevrim sayısı artmış, ancak ana dağılım bölgesi büyük ölçüde korunmuştur.

On dördüncü sayfadaki hızlı Fourier dönüşümü sonuçlarında ana frekans tepeleri özellikle yaklaşık 30 Hz altındaki bölgede benzer konumlarda kalmıştır. Yeniden yapılandırma, özgün kayıtta bulunmayan baskın ve yüksek genlikli yeni bir frekans üretmemiştir.

200 saniyelik spektrum, birim zaman başına Paris sürücüsünde yalnızca yüzde 0,2 fark göstermesi nedeniyle çatlak ilerleme analizinde seçilmiştir. Seçilen spektrumda 43.525 yağmur akışı çevrimi bulunmaktadır.

Başlangıç çatlağı nasıl oluşturulmuştur?

FRANC3D içinde kanat köküne yarı eliptik yüzey çatlağı yerleştirilmiştir. Çalışmanın “2 mm başlangıç çatlak uzunluğu” olarak adlandırdığı ölçü yüzey boyunca toplam \(2c=2\) mm’dir. Şekil oranı \(a/c=0{,}6\) olduğundan yarı uzunluk \(c=1\) mm ve başlangıç derinliği \(a=0{,}6\) mm’dir. On altıncı sayfadaki Şekil 17 bu geometrinin hem yüzeyden hem kesitten görünümünü göstermektedir.

Çatlak, üretim sırasında oluşabilecek küçük döküm kusurunu veya işletme sırasında buz darbesiyle oluşabilecek erken hasarı temsil etmektedir. Çatlak yönü, 1,16 saniyedeki çekme gerilmesi alanına göre seçilmiştir.

Paris yasası malzeme parametreleri nelerdir?

Bu değerler çalışma içinde yapılan yeni malzeme deneylerinden değil, kaynak gösterilen önceki nikel-alüminyum bronz çatlak ilerleme çalışmasından alınmıştır.

Yeniden yapılandırılmış spektrumun çatlak sonucu nasıl doğrulanmıştır?

Aynı çatlak geometrisi ve aynı hesaplama ayarlarıyla iki ayrı 51 adımlı ilerleme hesabı yapılmıştır:

• Özgün kısa gerilme kaydının sürekli tekrarı
• 200 saniyelik yeniden yapılandırılmış spektrumun sürekli tekrarı

Bu karşılaştırma, yeniden yapılandırılmış spektrumun aynı model ve aynı başlangıç koşulları altında özgün kısa kaydın çatlak ilerleme sonucuna yakın kaldığını göstermektedir. Bağımsız bir ölçümle doğrulama değildir; iki sayısal yükleme yaklaşımı arasındaki iç tutarlılık sınamasıdır.

Çatlak ilerleme süreci nasıl gelişmiştir?

On yedinci sayfadaki Şekil 18, yarı eliptik çatlak cephesinin 51 adım boyunca büyümesini göstermektedir. Çatlak cephesi genel olarak düzgün ve yarım daireye yakın biçimini korumuştur. Toplam çevrim sayısı ilk adımlarda hızla yükselmiş ve yaklaşık 25. adımdan sonra 4,0 × 10⁸ çevrim çevresinde yataylaşmıştır.

Bu görünüm, her hesap adımının eşit sayıda çevrime karşılık gelmediğini göstermektedir. İlk geometrik ilerleme adımları için daha fazla yük çevrimi gerekirken çatlak büyüdükçe aynı ilerleme adımı daha az çevrimde tamamlanmaktadır.

En yüksek Mode I gerilme şiddeti faktörü:

• İlk aşamada yaklaşık 340–400 MPa√mm düzeyindedir.
• 25. adım çevresinde yaklaşık 1.100 MPa√mm’ye yükselmiştir.
• 51. adımda yaklaşık 1.850 MPa√mm’ye ulaşmıştır.

Son değer bildirilen 3.952,85 MPa√mm kırılma tokluğuna ulaşmamaktadır. Bu nedenle 51. adım nihai kırılma veya kopma noktası değildir.

Çatlak ucundaki gerilme ve ağ yapısı ne göstermektedir?

On sekizinci sayfadaki Şekil 19, 1, 25 ve 51. adımlardaki çatlak açıklığını ve çatlak ucu gerilme alanını göstermektedir. Ağ yalnızca çatlak cephesi çevresinde yoğunlaştırılmış, uzak bölgelerde daha kaba bırakılmıştır.

İlk adımda yüksek gerilme küçük bir bölgeye yoğunlaşmaktadır. İlerleyen adımlarda çatlak açılması görünür hâle gelmekte, yüksek gerilme bölgesi genişlemekte ve iki çatlak ucundaki yoğunlaşma artmaktadır. Görüntülerde yer değiştirmeler 30 kat büyütülmüştür; gösterilen açıklık gerçek ölçekli deformasyon değildir.

Hangi kırılma modu baskındır?

On sekizinci ve on dokuzuncu sayfalardaki sonuçlar, üç gerilme şiddeti faktörü arasında \(K_I\)’in açık biçimde baskın olduğunu göstermektedir. \(K_I\) çatlak yüzeylerinin birbirinden açıldığı Mode I yüklemeyi temsil eder.

• \(K_I\), çatlak büyüdükçe düzenli ve güçlü biçimde artmıştır.
• \(K_{II}\), daha küçük genlikli ve dalgalı kalmıştır.
• \(K_{III}\), işaret değiştirmiş ancak \(K_I\)’den çok daha küçük olmuştur.
• J-integrali \(K_I\)’e benzer biçimde artmıştır.
• T-gerilmesi çatlak cephesinin büyük bölümünde negatif kalmıştır.

Bu sonuç, incelenen yükleme durumunda çatlağın esas olarak açılma modu tarafından ilerletildiğini, kayma bileşenlerinin ise yerel ve ikincil katkı yaptığını göstermektedir. Farklı buz temas açıları veya yükleme yönlerinde aynı mod sıralaması garanti edilmemektedir.

Başlangıç çatlak uzunluğu nasıl etkili olmuştur?

Şekil oranı \(a/c=0{,}6\) sabit tutularak 2, 3, 5 ve 10 mm başlangıç yüzey uzunlukları karşılaştırılmıştır. Bu grup için ortak karşılaştırma düzeyi, en yüksek \(K_I\)’in yaklaşık 1.200 MPa√mm ve en yüksek J-integralinin yaklaşık 1,2 × 10⁻⁵ MJ/mm² düzeyine ulaşmasıdır.

Yirminci sayfadaki Şekil 21’de başlangıç çatlağı uzadıkça bütün \(K_I\) ve J eğrileri sola kaymaktadır. Bu, aynı kırılma sürücü düzeyine daha az çevrimde ulaşıldığını göstermektedir. Etki doğrusal değildir: başlangıç uzunluğunun beş katına çıkması, referans çevrim sayısını yaklaşık onda bire indirmiştir.

Başlangıç çatlak şekil oranı nasıl etkili olmuştur?

Başlangıç yüzey uzunluğu 2 mm sabit tutulmuş ve \(a/c\) oranı 0,5, 0,6, 0,7, 0,8 ve 1,0 olarak değiştirilmiştir. Oran büyüdükçe çatlak cephesi daha derin ve daha dairesel hâle gelmektedir.

Bu karşılaştırma grubunda kullanılan ortak uç düzey yaklaşık 1.750 MPa√mm en yüksek \(K_I\)’dir. Bu referans, çatlak uzunluğu çalışmasındaki 1.200 MPa√mm düzeyinden farklıdır.

Şekil oranı da çatlak ilerlemesini hızlandırmıştır; ancak eğriler arasındaki ayrım başlangıç uzunluğu değişimindeki kadar büyük değildir. İncelenen parametre aralığında başlangıç çatlak uzunluğu daha baskın değişkendir.

Çalışmadaki “eşdeğer saat” değerleri nasıl yorumlanmalıdır?

Her 200 saniyelik yapay spektrumda 43.525 yağmur akışı çevrimi bulunmaktadır. Tablolardaki saat değerleri, gerekli çevrim sayısının bu sabit çevrim yoğunluğuna bölünmesiyle hesaplanmıştır:

\[ t_{\mathrm{eq}} = \frac{N_{\mathrm{ref}}}{43.525} \times \frac{200}{3.600} \]

Bu sürelerin önemli yorum sınırları vardır:

• Gerçek bir geminin kesintisiz olarak aynı buz kesme durumunda çalıştığını varsaymaz; yalnızca matematiksel tekrardır.
• Nihai kopma veya \(K_{\mathrm{IC}}\)’ye ulaşma süresi değildir.
• İki parametrik tabloda kullanılan referans \(K_I\) düzeyleri aynı değildir.
• Farklı görev profilleri, açık su seyri, değişken buz kalınlığı ve bakım duruşları hesaba katılmamıştır.
• Gerçek takvim saati veya güvenli bakım aralığı olarak kullanılamaz.

Yaklaşık 600 MPa√mm düzeyi ne anlama gelmektedir?

Araştırmacılar parametrik eğrilerde \(K_I\) yaklaşık 600 MPa√mm’yi geçtikten sonra eğimin belirgin biçimde arttığını, yaklaşık 1.100 MPa√mm sonrasında ise ilerlemenin çok hızlandığını gözlemlemiştir.

Bu değerler nikel-alüminyum bronz için evrensel bir “kararlı-kararsız geçiş sınırı” olarak yorumlanmamalıdır. Sonuç; kullanılan Paris katsayılarına, çatlak geometrisine, seçilen gerilme alanına, sayısal ilerleme adımlarına ve bu çalışmanın yük spektrumuna özgüdür. Deneysel eşik veya standardize edilmiş kabul sınırı değildir.

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

• Buz, su ve pervane yüklerini tek bir tam eşleştirilmiş modelde birleştirmektedir.
• Buz kırılması kohezyonlu elemanlarla ve parça ayrılmasıyla temsil edilmiştir.
• Kısa süreli gerilme verisi doğrudan tekrarlanmak yerine yorulma ölçütleriyle yeniden yapılandırılmıştır.
• Paris sürücüsü, yağmur akışı dağılımı, sınır sürekliliği ve tekrar çeşitliliği aynı amaç fonksiyonunda kullanılmıştır.
• Ana hiperparametrelerin ayrı ayrı çıkarıldığı duyarlılık analizi sunulmuştur.
• Yağmur akışı histogramı, aralık–ortalama matrisi, frekans spektrumu ve Paris sürücü oranı birlikte değerlendirilmiştir.
• Özgün ve yeniden yapılandırılmış yüklerle hesaplanan çatlak uzunlukları doğrudan karşılaştırılmıştır.
• Mode I, Mode II, Mode III, J-integrali ve T-gerilmesi çatlak cephesi boyunca incelenmiştir.
• Başlangıç çatlak uzunluğu ve şekil oranı ayrı parametrik gruplarda değerlendirilmiştir.

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

• Çalışmada yeni fiziksel buz kesme veya pervane gerilme deneyi yapılmamıştır.
• Modelin malzeme ve buz kırılma parametreleri önceki çalışmalardan alınmıştır.
• Tek bir R sınıfı pervane geometrisi incelenmiştir.
• Tek buz bloğu, tek giriş hızı, tek buz hızı ve tek pervane devri kullanılmıştır.
• Yalnızca bir kanat kökü elemanındaki gerilme geçmişi yeniden yapılandırılmıştır.
• Yeniden yapılandırma, kısa kayıtta hiç bulunmayan yeni deniz veya buz durumlarını üretemez.
• 25–1.000 saniyelik bütün spektrumlar aynı kısa sayısal kayıttan türetilmiştir.
• Spektrum doğrulaması bağımsız bir uzun süreli gerçek gerilme kaydıyla yapılmamıştır.
• Özgün ve yeniden yapılandırılmış spektrum karşılaştırması aynı sayısal model içinde yapılmıştır.
• 1,16 saniyedeki tek uzamsal gerilme alanının değişken genlikli yükleme boyunca temsil edici olduğu varsayılmıştır.
• Paris yasasının geçerli olduğu kararlı çatlak ilerleme bölgesi dışındaki davranış ayrıntılı modellenmemiştir.
• Çatlak kapanması, çevrim sırası etkisi, aşırı yük geciktirmesi ve artık gerilmeler açıkça ele alınmamıştır.
• Üretim kusurlarının gerçek dağılımı veya olasılıksal belirsizliği hesaba katılmamıştır.
• Hesaplanan referans süreler nihai kırılma veya gerçek görev ömrü değildir.
• Çatlak uzunluğu ve şekil oranı çalışmalarında aynı referans \(K_I\) düzeyi kullanılmamıştır.
• Çalışma hakem değerlendirmesinden geçmemiştir.

Çalışma neyi desteklemektedir?

• Kullanılan S-ALE ve kohezyonlu eleman düzeni, darbeli buz kesme kuvvetleri ve buz parçalanması üretebilmiştir.
• İncelenen durumda kanat kökü yüksek çekme ve eşdeğer gerilme bölgesi olarak belirlenmiştir.
• Önerilen blok tabanlı yöntem, özgün kısa kaydın yağmur akışı ve baskın frekans özelliklerini yaklaşık olarak koruyabilmiştir.
• 200 saniyelik spektrumun normalize Paris sürücüsü özgün kaydın yüzde 0,2 üzerindedir.
• Yeniden yapılandırılmış ve özgün yüklemelerin 51 adımlık çatlak uzunluğu sonuçları yüzde 2,55 farklıdır.
• İncelenen yükleme durumunda Mode I açılma bileşeni çatlak ilerlemesinde baskındır.
• Başlangıç çatlak uzunluğu, seçilen parametre aralığında şekil oranından daha güçlü etki göstermiştir.
• Daha büyük başlangıç kusurları referans kırılma sürücüsüne çok daha az çevrimde ulaşmıştır.

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

• Yöntemin bütün buz sınıfı pervanelerde aynı doğruluğu sağlayacağını kanıtlamaz.
• Modelin gerçek kutup koşullarındaki yükleri deneysel belirsizlik içinde doğru tahmin ettiğini göstermez.
• Yeniden yapılandırılmış spektrumun gerçek uzun süreli görev spektrumuyla aynı olduğunu göstermez.
• 79,15 mm çatlak uzunluğunun gerçek bir pervanede aynı çevrim sayısında oluşacağını kanıtlamaz.
• Tablolardaki saatlerin gerçek seyir veya takvim ömrü olduğunu göstermez.
• 2 mm, 5 mm veya 10 mm’nin izin verilebilir bakım sınırı olduğunu belirlemez.
• Yaklaşık 600 MPa√mm değerinin evrensel hızlanma eşiği olduğunu kanıtlamaz.
• Çatlağın 51. adımda nihai kırılmaya ulaştığını göstermez.
• Farklı buz kalınlıkları, kavite, korozyon veya değişken sıcaklık altında aynı sonuçların geçerli olduğunu göstermez.
• Gerçek bakım aralığını veya ekonomik olarak en uygun muayene planını hesaplamaz.

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

Teknik yöntem özeti

Gerilme spektrumu yeniden yapılandırma ayarları

Ana geçici yük bulguları

Spektrum doğrulama bulguları

Başlangıç çatlak geometrisi sonuçlarının karşılaştırması

İki yüzde farklı karşılaştırma gruplarına ve farklı referans kırılma sürücüsü düzeylerine aittir. Mutlak saat değerleri iki tablo arasında doğrudan karşılaştırılmamalıdır.

Şekillerin temel mesajları

• Şekil 1: Nikel-alüminyum bronz için önceki çalışmadan alınan düşük ve yüksek çevrimli S-N eğrisini göstermektedir.
• Şekil 2: Buz-su-pervane analizinden spektrum üretimine ve çatlak değerlendirmesine uzanan beş aşamalı yöntemi göstermektedir.
• Şekil 3: S-ALE akışkan ve yapı düğümleri arasındaki ceza kuvveti eşleştirmesini göstermektedir.
• Şekil 4: Buz elemanları arasındaki kohezyon bölgesinin hasar görerek çatlak oluşturmasını göstermektedir.
• Şekil 5: R sınıfı pervane geometrisini ve hava-su-buz hesaplama alanını göstermektedir.
• Şekil 6: Buz kesme sırasında −150 ile +230 kN arasında değişen darbeli pervane itmesini göstermektedir.
• Şekil 7: Buz bloğunun kanat çevresinde ezilmesi ve parçaların suya dağılmasını göstermektedir.
• Şekil 8: 1,16 saniyedeki von Mises ve en yüksek asal gerilme dağılımlarını göstermektedir.
• Şekil 9: Kritik kanat kökü elemanını ve 0,482–2,730 saniyelik asal gerilme geçmişini göstermektedir.
• Şekil 10: Yağmur akışı aralık dağılımlarının hafıza katsayısıyla birleştirilmesini göstermektedir.
• Şekil 11: Yeniden yapılandırma hiperparametrelerinin çıkarılması veya değiştirilmesinin ölçütlere etkisini göstermektedir.
• Şekil 12–13: Farklı sürelerdeki yağmur akışı histogramı ve aralık–ortalama matrislerinin ana yapısını koruduğunu göstermektedir.
• Şekil 14–15: Frekans yapısını ve beş sürenin normalize Paris sürücüsü oranlarını göstermektedir.
• Şekil 16: Çatlak hesabında kullanılan 200 saniyelik spektrumun çevrim sırasını göstermektedir.
• Şekil 17: Yüzey uzunluğu 2 mm ve derinliği 0,6 mm olan başlangıç çatlağını göstermektedir.
• Şekil 18–20: Çatlak cephesi, çevrim birikimi, \(K_I\), \(K_{II}\), \(K_{III}\), J-integrali ve T-gerilmesinin 51 adımlık gelişimini göstermektedir.
• Şekil 21: Başlangıç çatlak uzunluğu büyüdükçe kırılma sürücüsünün daha az çevrimde yükseldiğini göstermektedir.
• Şekil 22: Daha büyük \(a/c\) oranlarının ilerlemeyi hızlandırdığını, ancak etkinin uzunluk değişimine göre daha sınırlı olduğunu göstermektedir.

Kanıt düzeyine göre sonuçların durumu

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı: A Paris-law equivalent spectrum reconstruction method for fatigue crack propagation of ice-class propellers in polar environments

Yazarlar: Xiaoli Yuan, Xinyao Wang, Weipeng Xiong, Mi-An Xue, Yanqi Wang ve Chao Wang.

Yazar sıralaması: Xiaoli Yuan birinci, Xinyao Wang ikinci, Weipeng Xiong üçüncü, Mi-An Xue dördüncü, Yanqi Wang beşinci ve Chao Wang altıncı yazardır.

Eşit katkı bilgisi: Yüklenen sürümde eşit katkıda bulunan yazar bilgisi yer almamaktadır.

Sorumlu yazar: Weipeng Xiong.

Sorumlu yazar e-postası: xiongweipeng@hhu.edu.cn

Kurumsal bağlantılar:

1. College of Harbour, Coastal and Offshore Engineering, Hohai University, Nanjing, Jiangsu 210098, Çin.
2. College of Mechanics and Engineering Science, Hohai University, Nanjing, Jiangsu 211100, Çin.
3. College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, Çin.

Yazar–kurum eşleştirmesi:

• Xiaoli Yuan: Hohai University bünyesindeki birinci ve ikinci kurumlar.
• Xinyao Wang: Hohai University bünyesindeki birinci kurum.
• Weipeng Xiong: Hohai University bünyesindeki birinci kurum.
• Mi-An Xue: Hohai University bünyesindeki birinci kurum.
• Yanqi Wang: Hohai University bünyesindeki birinci kurum.
• Chao Wang: Harbin Engineering University bünyesindeki üçüncü kurum.

DOI:10.2139/ssrn.7197210

Yayın platformu: SSRN.

Resmî kaynak bağlantısı:SSRN çalışma kaydı

Yayın yılı: 2026.

Kesin SSRN yükleme tarihi: Yüklenen sürümde açık bir SSRN gönderim tarihi yer almamaktadır; kesin gün bilgisi bu incelemede bağımsız olarak doğrulanamamıştır.

Dergi durumu: Sayfa üstbilgisinde “Ocean Engineering xxx (xxxx) xxxxxx” ifadesi bulunmaktadır. Cilt, sayı, makale numarası, kabul tarihi veya nihai çevrim içi yayın kaydı bulunmadığı için hakemli dergi yayını doğrulanamamıştır.

Yayınevi: Nihai hakemli yayın için kesin bir yayınevi kaydı bulunmamaktadır. İncelenen sürüm SSRN üzerinde yayımlanan bir preprinttir.

Kaynak türü: Sayısal akışkan-yapı etkileşimi, algoritmik gerilme spektrumu yeniden yapılandırması ve kırılma mekaniği tabanlı parametrik analiz; preprint.

Hakemlik durumu: Bu çalışma bir preprinttir ve hakem değerlendirmesinden geçmemiştir.

Yazar katkıları:

• Xiaoli Yuan: Kavramsallaştırma, yöntem ve gözden geçirme-düzenleme.
• Xinyao Wang: İlk taslak, doğrulama, veri düzenleme ve yöntem.
• Weipeng Xiong: Yöntem ve gözden geçirme-düzenleme.
• Mi-An Xue: Yöntem ve gözden geçirme-düzenleme.
• Yanqi Wang: Kavramsallaştırma ve gözden geçirme-düzenleme.
• Chao Wang: Yöntem ve gözden geçirme-düzenleme.

Finansman: Çalışma; Çin Ulusal Doğa Bilimleri Vakfı 52501389 numaralı fonu, Jiangsu Eyaleti Doğa Bilimleri Vakfı BK20241523 numaralı fonu, Shanghai Jiao Tong University State Key Laboratory of Ocean Engineering GKZD010091 numaralı fonu ve KY10100250185 numaralı temel araştırma projesi tarafından desteklenmiştir.

Çıkar çatışması: Yazarlar çalışmayı etkileyebilecek bilinen bir mali çıkar veya kişisel ilişki bildirmemiştir.

Veri erişimi: Verilerin talep üzerine paylaşılacağı belirtilmiştir. Açık veri deposu bağlantısı verilmemiştir.

Kod erişimi: Gerilme spektrumu yeniden yapılandırma algoritmasının kodu için açık depo bağlantısı yer almamaktadır.

Etik kurul: Çalışma insan veya hayvan katılımcı içermeyen sayısal mühendislik araştırmasıdır. Yüklenen sürümde etik kurul bilgisi yer almamaktadır.

Bu Türkçe açıklama yalnızca yüklenen 23 sayfalık çalışmanın metni, denklemleri, tabloları, sayısal model ayrıntıları ve Şekil 1–22’deki görsel sonuçlar temel alınarak hazırlanmıştır. Çalışmada bulunmayan gerçek deniz deneyi, tam ölçekli ölçüm, bakım standardı, izin verilebilir çatlak boyutu veya güvenli işletme süresi eklenmemiştir. Dış kaynaklar yalnızca bibliyografik kimliğin, DOI’nin ve resmî SSRN kaydının doğrulanması için kullanılmıştır.

Çalışmanın en güçlü yönü, kısa süreli buz-su-pervane simülasyonunu çatlak sürücüsüne göre seçilen uzun süreli bir spektrumla ve FRANC3D çatlak ilerlemesiyle birleştirmesidir. En önemli doğrulama sonucu, özgün ve yeniden yapılandırılmış spektrumların 51 adım sonunda yüzde 2,55 farklı çatlak uzunluğu üretmesidir. Bununla birlikte bu karşılaştırma aynı sayısal ortam içindedir ve gerçek uzun süreli gerilme ölçümünün yerini tutmamaktadır.

Başlangıç çatlak uzunluğunun 2 mm’den 10 mm’ye çıkarılmasıyla bildirilen yüzde 90,7 azalma, nihai yapısal kırılma ömrüne değil, uzunluk parametrik çalışmasında seçilen referans \(K_I\) ve J-integrali düzeyine ulaşmak için gereken çevrim sayısına aittir. Benzer biçimde tablolardaki eşdeğer saatler gerçek takvim süresi veya bakım aralığı değildir.

Yöntemin mühendislik uygulamasına geçebilmesi için farklı pervane geometrileri, malzeme partileri, başlangıç kusuru dağılımları, buz kalınlıkları, buz hızları, devir sayıları ve görev profilleri altında sınanması; bağımsız uzun süreli yük kayıtlarıyla karşılaştırılması ve deneysel çatlak ilerleme testleriyle doğrulanması gerekmektedir.