Бул изилдөө, муз класстуу bir кеме винтsinin buz, суу жана винт arasındaki кыска мөөнөттүү fakat жогорку чыгымli сандык өз ара аракеттенүүinden elde edмененn чыңалуу жазууларыnı узак мөөнөттүү чарчоо жаракасы мененrleme талдоосуnde kullanabilmek үчүн yeni bir чыңалуу спектрu кайра куруу ыкмасы geliştirmiştir. Изилдөөчүлөр; LS-DYNA үчүнdeki түзүмlandırılmış эркин Lagrange–Euler ыкмасы, когезиялуу buz сынууsı модельi, rainflow цикл саноо, Paris мыйзамы жана FRANC3D жарака мененrleme yazılımını tek bir hesaplama zincirinde birleştirmiştir. Yeniden түзүмlandırılmış спектрla 51 мененrleme adımı sonunda hesaplanan жарака узундугу 79,15 mm, özgün kısa спектрun sonсууz tekrarıyla hesaplanan uzunluk болсо 77,18 mm olmuş жана aradaki fark yüzde 2,55 катары билдирилген. Temel чектөө, натыйжаларыn tek bir сандык изилдөө durumuna dayanması жана чыныгы винт, deniz buzu же уюлдук seferi маалыматleriyle түздөн-түз doğrulanmamış olmasıdır.

Sayısal модельde 4,115 metre çapında, dört калакlı, никель-алюминий коло R sınıfı bir винт; üç metre uzunluğunda bir муз блогу, суу жана аба alanıyla tamamen eşleştirilmiştir. Kararlı катары tanımlanan 1–2 saniyelik kesitte винт түртүүsi yaklaşık −150 kN менен +230 kN arasında жогорку жыштыкlı жана darbeli biçimde değişmiştir. Çatlak талдоосу үчүн seçмененn 1,16 saniyelik anda калак түбүndeki en жогорку von Mболсоs чыңалууsi 254 MPa, en жогорку башкы чыңалуу болсо 322 MPa катары эсептелген. Çekme yönündeki башкы чыңалууnin калак түбүnün arka yüzünde yoğunlaşması, başlangıç çatlтармакının bu bölgeye yerleştirilmesinin gerekçesi olmuştur.

Yeniden түзүмlandırma ыкмасы kısa чыңалуу kaydını nokta nokta rastgele çoğaltmamaktadır. Özgün жазууtan çıkarılan чыңалуу blokları; Paris мыйзамыna göre birim zamanda oluşturdukları жарака sürücüsü, rainflow чыңалуу aralığı dтармакılımı, blok birleşimindeki süreklilik жана yakın zamanda kullanılan blokların tekrarından kaçınma ölçütleriyle seçilmektedir. Bu ыкмаle 25, 100, 200, 500 жана 1.000 saniyelik спектрlar üretilmiştir. Bunların süreye göre normalize edilmiş Paris мыйзамы sürücü катышları sırasıyla 0,982, 0,983, 1,002, 1,005 жана 1,020’dir. Çatlak мененrleme hesabında, özgün kayda en yakın катышlardan birini sтармакlayan 200 saniyelik спектр колдонулган.

Başlangıç жарака узундугу, incelenen өзгөрмөler arasında en күчтүү etkiyi көрсөткөн. Çatlak uzunluğu 2 mm’den 10 mm’ye çıkarıldığında, изилдөөдө tanımlanan referans сынуу sürücü деңгээлine ulaşmak үчүн gereken цикл sayısı 4,05 × 10⁸ циклden 0,38 × 10⁸ циклe düşmüş; azalma yüzde 90,7 olmuştur. Başlangıç жарака сүрөт катышı a/c = 0,5’ten 1,0’a жүкseltildiğinde болсо kendi салыштыруу grubundaki referans деңгээлe ulaşmak үчүн gereken цикл sayısı yüzde 42,9 azalmıştır. Bu маанилер nihai kopma ömrü değil, önceden belirlenen сынуу sürücüsü деңгээлlerine erişme süreleridir.

Түркия жагынан: Yöntem; Түркиядаki кеме inşa, deniz mühendisliği, винт üretimi, klaslama изилдөөları жана уюлдук koşullarına uygun кеме ekipmanı geliştirme изилдөөlarında кыска мөөнөттүү жогорку doğruluklu simülasyonlarla узак мөөнөттүү жарака баалуулукlendirmesi arasında байланыш kurmak үчүн incelenebilir. Түркияда uygulanmadan önce yerel üretim ыкмаларiyle oluşturulan винт материалlerinin Paris мыйзамы katsayıları, döküm кемтиктерı, artık чыңалууler жана чыныгы жарака dтармакılımları экспериментsel катары belirlenmelidir. Модель ар башка buz kalınlıkları, мененrleme hızları, винт devirleri, kavite koşulları жана görev profilleriyle yeniden çalıştırılmalı; натыйжалар tam ölçekli же модель ölçekli өлчөөлөрle doğrulanmalıdır. Изилдөөdan Түркиядаki belirli bir винттин коопсуз kullanım süresi, bakım aralığı же izin маалыматlebilir жарака boyutu түздөн-түз çıkarılamaz.

Изилдөөнүн çözmeye çalıştığı temel problem эмне?

Buz sınıfı кемеlerin винтleri, normal hidrodinamik күчlerin yanında калактардын buzu kesmesi, ezmesi жана parçaları суу үчүнde hızlandırması эмне үчүнiyle кыска мөөнөттүү жана жогорку genlikli жүкlerle karşılaşabilir. Buzla temas sürekli жана düzgün değildir. Bir калак муз блогуna yaklaşır, temasa girer, buzu kırar, parçaları ayırır жана yeniden сууdaki hidrodinamik изилдөө durumuna döner. Bu olayların tekrarlanması, винт калак түбүnde өзгөрмө genlikli чыңалуу циклleri oluşturur.

Bu өз ара аракеттенүүin жогорку doğruluklu сандык çözümü маанилүү hesaplama kaynтармакı gerektirir. Saniyeler деңгээлindeki bir buz-суу-винт simülasyonu bменен saatlerce sürebilir. Buna karşılık чарчоо жаракасыnın мененrlemesi milyonlarca же yüz milyonlarca жүк циклinin баалуулукlendirilmesini gerektirir. Изилдөөnın hedeflediği boşluk, кыска мөөнөттүү сандык чыңалуу kaydıyla узак мөөнөттүү сынуу mekaniği hesabı arasındaki zaman ölçeği farkıdır.

Изилдөөчүлөр gelecekteki her чыңалуу noktasını tahmin etmeye изилдөөmaktadır. Amaç, kısa kaydın чарчоо bakımından belirleyici өзгөчөliklerini koruyan daha uzun bir yapay чыңалуу geçmişi üretmektir. Korunması hedeflenen өзгөчөlikler yalnızca орточо жана standart ийилүү değildir; жогорку чыңалуу aralıkları, rainflow цикл dтармакılımı, baskın жыштыкlar жана Paris мыйзамыna göre жарака мененrleme potansiyeli de ыкмаe dâhil edilmiştir.

Изилдөөnın бардыкleşik hesaplama çerçeжанаsi кантипdır?

Dördüncü sayfadaki Şekil 2, ыкмасы beş баскычda көрсөтөт:

1. Buz, суу жана винттин tamamen eşleştirilmiş geçici сандык талдоосу
2. Kısa süreli чыңалуу kaydının узак мөөнөттүү чарчоо eşбаалуулук спектрa dönüştürülmesi
3. Kritik калак түбү bölgesine yarı eliptik başlangıç çatlтармакının yerleştirilmesi
4. Paris мыйзамы жана FRANC3D менен жараканын өсүшүnin hesaplanması
5. Çatlak riski, bakım жана tasarım açısından mühendislik баалуулукlendirmesi

LS-DYNA талдоосу buz сынууsını жана винттин geçici түзүмsal чыңалууsini üretmektedir. Yeniden түзүмlandırma algoritması, seçмененn калак түбү элементındaki kısa чыңалуу geçmişini uzatmaktadır. Ansys Workbench жана FRANC3D болсо 1,16 saniyedeki üç boyutlu чыңалуу alanı üzerine yarı eliptik бет çatlтармакını yerleştirerek her чыңалуу циклi altında жарака cephesinin кантип мененrlediğini hesaplamaktadır.

S-ALE ыкмасы эмне үчүн колдонулган?

Yapılandırılmış эркин Lagrange–Euler, yani S-ALE ыкмасы; суюктук alanı менен büжүк deformasyon gösteren katı түзүмların aynı hesaplama үчүнde etkмененşmesini sтармакlamak үчүн колдонулган. Pervane жана buz түзүмsal тармакla, аба жана суу болсо referans суюктук тармакıyla temsil edilmiştir. Akışkan жана katı arasındaki momentum aktarımı ceza күчi tabanlı eşleştirme ыкмасыyle чыныгыleştirilmiştir.

Изилдөөda kullanılan kütle, momentum жана энергия 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 bмененşenlerini, \(e\) özgül iç энергияyi, \(w_i\) материал менен referans тармак arasındaki göreli hızı жана \(\sigma_{ij}\) чыңалуу tensörünü temsil etmektedir. Yöntem, винт dönerken сууyun hareketini жана buz parçalarının суюктук үчүнde taşınmasını aynı çözüm үчүнde hesaplamayı amaçlamaktadır.

Buz сынууsı кантип модельlenmiştir?

Buz bloğu, elastoplastik katı элементlar менен элементlar arasına yerleştirмененn когезиялуу сынуу элементlarından oluşturulmuştur. Kohezyonlu элементlar başlangıçta komşu buz элементlarını birbirine bтармакlamaktadır. Normal же жанамаsel созулуу belirlenen dayanımı aştığında hasar başlamakta, өзүнчөlma arttıkça taşıma кубаттуулугу azalmaktadır. Kritik өзүнчөlma энергияsine ulaşıldığında bтармак tamamen kopmakta жана buz parçaları birbirinden өзүнчөlmaktadır.

Beşinci sayfadaki Şekil 4, iki buz элементı arasındaki kohezyon bölgesinin жарака oluşturarak өзүнчөlmasını көрсөтөт. Bu yaklaşımda жарака yolu önceden tek bir çizgi катары tanımlanmamış; uygun элемент arayüzlerinde hasarın gelişmesine izin маалыматlmiştir.

Kohezyonlu элементlar үчүн тыгыздык 900 kg/m³, en жогорку созулуу 4 MPa, энергия salım катышı 433 J/m² жана rijitlik 50.000 MPa катары колдонулган. Buz gövdesinin серпилгичтик модулу 264,8 MPa, агуу чыңалуусу 2,61 MPa жана Poisson катышы 0,30’dur. Bu parametreler bu изилдөөдө yeni экспериментlerle elde edilmemiş, изилдөөчүлөрın önceki модельleme изилдөөlarına dayandırılmıştır.

Pervane жана эсептөө аймагыnın өзгөчөlikleri nelerdir?

Su giriş hızı −X yönünde 3,62 m/s, buzun мененrleme hızı +X yönünde 0,7 m/s жана винт айлануу ылдамдыгы 2,698 devir/saniye катары belirlenmiştir. Akışkan alanında 0,19 metrelik түзүмlandırılmış элементlar колдонулган. Pervanenin X yönündeki ötelemesi sınırlandırılmış, diğer serbestlik derecelerine izin маалыматlmiştir.

Ağ duyarlılık талдоосу ne көрсөткөн?

Pervane үчүн 40, 50, 60 жана 70 mm olmak üzere dört тармак büжүкlüğü салыштырылган. İlk buz kesme баскычsındaki 0,5–1,25 saniye aralığında hesaplanan орточо күчler şöyledir:

40 жана 50 mm тармакların натыйжалары birbirine yakınken daha kaba тармакlarda belirgin өзүнчөşma görülmüştür. Изилдөөчүлөр hesaplama süresi менен çözüm kararlılığı arasında denge kurmak үчүн 50 mm тармакı seçmiştir. Bu салыштыруу сандык тармак yakınsamasını колдойт; чыныгы өлчөөe karşı mutlak doğruluk kanıtı değildir.

Buz kesme sırasında винт жүкleri кантип değişmiştir?

Sekizinci sayfadaki Şekil 6, винт түртүүsinin düzenli bir sinüs dalgası olmadığını көрсөтөт. Bir saniyelik kararlı kesitte маанилер yaklaşık −150 kN менен +230 kN arasında değişmiştir. Negatif tepeler, buz kesme direncinin bazı anlarda винттин pozitif hidrodinamik түртүүsinden daha büжүк hâle geldiğini көрсөтөт.

Aynı sayfadaki Şekil 7’de buz, калак temas bölgesinde çok sayıda küçük parçaya өзүнчөlmaktadır. Parçalar калак geçişi çevresinde kümelenmekte, ardından суу akışı жана винт hareketiyle taşınmaktadır. Görsel натыйжалар жүкün sürekli temas yerine kesme, ezilme, kopma жана parça өзүнчөlması evrelerinden oluştuğunu көрсөтөт.

Çatlak үчүн kritik bölge кантип seçilmiştir?

Kanat kökündeki элементlar arasında en тармакır en жогорку башкы чыңалуу geçmişine sahip 505354 numaralı элемент seçilmiştir. Gerilme kaydı 0,482–2,730 saniye aralığından alınmıştır. Kayıttaki en жогорку созулуу durumu 1,16 saniyede oluştuğu үчүн bu andaki üç boyutlu чыңалуу alanı daha sonraki FRANC3D талдоосуnin temel alanı катары колдонулган.

Dokuzuncu sayfadaki Şekil 8 жана Şekil 9’a göre:

• En жогорку von Mболсоs чыңалууsi 254 MPa’dır.
• Pervane материалsinin bildirмененn агуу чыңалуусу 268,3 MPa’dır.
• En жогорку башкы чыңалуу 322 MPa’dır.
• Çekme чыңалууsi калак түбүnün arka yüzünde yoğunlaşmaktadır.

Von Mболсоs баалуулукi bildirмененn akma sınırının yaklaşık yüzde 94,7’sine karşılık gelmektedir. En жогорку башкы чыңалуу болсо түздөн-түз akma ölçütü катары değil, Mode I açılma çatlтармакının yönünü жана созулуу sürücüsünü belirlemek үчүн колдонулган.

Neden түздөн-түз kısa чыңалуу 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 иштетүү özgün жазууta bulunmayan yapay periyodiklik oluşturabilir. Gerilme noktalarının bтармакımsız biçimde rastgele seçilmesi болсо zaman sırasını, жыштык түзүмsını жана циклlerin орточо чыңалууlerini bozabilir.

Önerмененn ыкма bu iki uç yaklaşım arasında yer almaktadır. Gerilme geçmişinden kayan pencereyle kısa bloklar çıkarılmakta жана blok үчүнdeki zaman sırası korunmaktadır. Yeni uzun спектр, özgün kaydı başlangıç çapası катары kullanmakta жана seçмененn blokları doğrusal örtüşmeyle birbirine eklemektedir.

Yeniden түзүмlandırma чыгым fonksiyonu кантип аныкталган?

Her aday чыңалуу bloğu үчүн toplam чыгым şu biçimde эсептелген:

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

Ana салмакlar \(\alpha_d=1{,}8\), \(\alpha_h=0{,}9\) жана \(\alpha_b=0{,}12\) катары belirlenmiştir. Her hedef süre үчүн birden fazla rastgele başlangıçla tam спектр oluşturulmuş жана özgün kayda en iyi uyumu sтармакlayan натыйжа seçilmiştir.

Paris мыйзамыna göre eşбаалуулук жарака sürücüsü кантип эсептелген?

Paris мыйзамы, kararlı чарчоо жаракасы мененrleme hızını чыңалуу şiddeti faktörü aralığıyla ilişkмененndirir:

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

Gerilme şiddeti faktörü aralığı yaklaşık катары:

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

biçiminde yazılmıştır. Aynı материал, aynı жарака konumu, aynı geometri жана aynı referans жарака durumu karşılaştırıldığında \(C(Y\sqrt{\pi a})^m\) ortak çarpandır. Bu эмне үчүнle göreli жарака sürücüsü, rainflowyla bulunan чыңалуу aralıklarının \(m\) күчi üzerinden karşılaştırılabilir:

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

• \(G\), birim zaman başına Paris tipi жарака sürücüsüdür.
• \(T\), чыңалуу geçmişinin süresidir.
• \(\Delta\sigma_i\), i numaralı rainflow чыңалуу aralığıdır.
• \(n_i\), bu aralığın цикл sayısıdır.
• \(m\), Paris мыйзамы üssüdür.

Paris sürücüsü чыгымi:

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

катары аныкталган. Oranın bire yaklaşması, aday blok eklendiğinde birim zaman başına жарака мененrleme potansiyelinin özgün жазууla benzer kaldığını көрсөтөт.

Yтармакmur akışı dтармакılımı эмне үчүн өзүнчөca korunmuştur?

İki чыңалуу geçmişi aynı toplam Paris sürücüsüne sahip olsa bменен цикл dтармакılımları ар башка болушу мүмкүн. Bir geçmiş çok sayıda orta aralıklı цикл, diğeri az sayıda çok жогорку aralıklı цикл içerebilir. Paris мыйзамыndaki жогорку üs эмне үчүнiyle büжүк aralıklı циклlerin etkisi өзгөчөlikle жогоркуtir.

Yтармакmur akışı чыңалуу aralıkları \(B\) sınıfa өзүнчөlmış жана мүмкүнlık dтармакılımına dönüştürülmüştür. Güncellenen dтармакılım:

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

eşitliğiyle hesaplanmış жана \(\beta=0{,}92\) колдонулган. Özgün жана güncellenmiş dтармакılım arasındaki fark toplam değişim uzaklığıyla ölçülmüştür:

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

On birinci sayfadaki Şekil 10, çalışan dтармакılımla aday bloğun dтармакılımının кантип birleştirildiğini şematik катары көрсөтөт.

Blok birleşimindeki yapay циклler кантип sınırlandırılmıştır?

Bir bloğun son чыңалууsi менен yeni bloğun ilk чыңалууsi arasında büжүк sıçrama oluşursa rainflow sayımı özgün жазууta bulunmayan жогорку aralıklı bir цикл oluşturabilir. Bu эмне үчүнle sınır чыгымi:

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

şeklinde аныкталган. Bloklar arasında doğrusal örtüşme de uygulanmıştır.

Son beş seçim үчүнde kullanılan adaylara \(0{,}35\) tekrar cezası маалыматlmiştir. Rastgele pertürbasyonun standart ийилүүsı \(0{,}015\)’tir. Bu iki bмененşen, algoritmanın sürekli aynı birkaç bloğu seçmesini жана tek bir yerel çözüme kilitlenmesini azaltmayı amaçlamaktadır.

Hiperparametre çıkarım талдоосу ne көрсөткөн?

On ikinci sayfadaki Şekil 11’de ana салмакlar, dтармакılım hafıza katsayısı, rastgele terim жана tekrar cezası değiştirмененrek натыйжа ölçütlerinin nominal durumdan yüzde farkı gösterilmiştir.

En duyarlı өзгөрмө Paris sürücüsü тармакı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 натыйжа, tek bir ölçütün iyмененşmesinin бардык спектрun daha iyi болгону anlamına gelmediğini көрсөтөт. Yüksek kuyruk uyumu iyмененşirken жалпы sürücü uyumu bozulabilmektedir.

Histogram тармакı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 тармакırlığı çıkarıldığında жалпы Paris sürücüsü ölçütleri yaklaşık yüzde 10–11 kötüleşmiştir. Изилдөөчүлөр nominal ayarların hiçbir ölçütte tek başına en iyi sonucu жанаrmediğini, бирок ölçütler arasında dengeli bir çözüm oluşturduğunu баалуулукlendirmiştir.

Yeniden түзүмlandırılmış спектрlar özgün kayda ne kadar benzemektedir?

Yöntem 25, 100, 200, 500 жана 1.000 saniyelik beş спектр üretmiştir. Yтармакmur akışı гистограммаları, циклlerin büжүк bölümünün yaklaşık 50 MPa altındaki чыңалуу aralıklarında toplandığını жана жогорку aralıklarda uzun fakat seyrek bir kuyruk bulunduğunu көрсөткөн.

On üçüncü sayfadaki aralık–орточо matrislerinde тыгыздык merkezi төмөн чыңалуу aralıkları менен yaklaşık 100–150 MPa орточо чыңалуу çevresinde kalmıştır. Süre uzadıkça toplam цикл sayısı artmış, бирок ana dтармакılım bölgesi büжүк ölçüde korunmuştur.

On dördüncü sayfadaki ылдам Fourier dönüşümü натыйжаларыnda ana жыштык tepeleri өзгөчөlikle yaklaşık 30 Hz altındaki bölgede benzer konumlarda kalmıştır. Yeniden түзүмlandırma, özgün жазууta bulunmayan baskın жана жогорку genlikli yeni bir жыштык üretmemiştir.

200 saniyelik спектр, birim zaman başına Paris sürücüsünde yalnızca yüzde 0,2 fark göstermesi эмне үчүнiyle жарака мененrleme талдоосуnde seçilmiştir. Seçмененn спектрda 43.525 rainflow циклi bulunmaktadır.

Başlangıç çatlтармакı кантип oluşturulmuştur?

FRANC3D үчүнde калак түбүne yarı eliptik бет çatlтармакı yerleştirilmiştir. Изилдөөnın “2 mm баштапкы жарака uzunluğu” катары adlandırdığı ölçü бет boyunca toplam \(2c=2\) mm’dir. Şekil катышı \(a/c=0{,}6\) болгонуndan yarı uzunluk \(c=1\) mm жана başlangıç derinliği \(a=0{,}6\) mm’dir. On altıncı sayfadaki Şekil 17 bu geometrinin hem бетden hem kesitten görünümünü көрсөтөт.

Çatlak, üretim sırasında oluşabмененcek küçük döküm кемтикunu же ишкана sırasında buz darbesiyle oluşabмененcek erken hasarı temsil etmektedir. Çatlak yönü, 1,16 saniyedeki созулуу чыңалууsi alanına göre seçilmiştir.

Paris мыйзамы материал parametreleri nelerdir?

Bu маанилер изилдөө үчүнde түзүмlan yeni материал экспериментlerinden değil, булак gösterмененn önceki никель-алюминий коло жарака мененrleme изилдөөsından alınmıştır.

Yeniden түзүмlandırılmış спектрun жарака sonucu кантип doğrulanmıştır?

Aynı жарака geometrisi жана aynı hesaplama ayarlarıyla iki өзүнчө 51 adımlı мененrleme hesabı түзүмlmıştır:

• Özgün kısa чыңалуу kaydının sürekli tekrarı
• 200 saniyelik yeniden түзүмlandırılmış спектрun sürekli tekrarı

Bu салыштыруу, yeniden түзүмlandırılmış спектрun aynı модель жана aynı başlangıç koşulları altında özgün kısa kaydın жарака мененrleme sonucuna yakın kaldığını көрсөтөт. Bтармакımsız bir өлчөөle валидация değildir; iki сандык жүктөө yaklaşımı arasındaki iç tutarlılık sınamasıdır.

Çatlak мененrleme süreci кантип gelişmiştir?

On yedinci sayfadaki Şekil 18, yarı eliptik жарака cephesinin 51 adım boyunca büyümesini көрсөтөт. Çatlak cephesi жалпы катары düzgün жана yarım daireye yakın biçimini korumuştur. Toplam цикл sayısı ilk adımlarda hızla жүкselmiş жана yaklaşık 25. adımdan sonra 4,0 × 10⁸ цикл çevresinde yataylaşmıştır.

Bu görünüm, her hesap adımının eşit sayıda циклe karşılık gelmediğini көрсөтөт. İlk geometrik мененrleme adımları үчүн daha fazla жүк циклi gerekirken жарака büyüdükçe aynı мененrleme adımı daha az циклde tamamlanmaktadır.

En жогорку Mode I чыңалуу şiddeti faktörü:

• İlk баскычda yaklaşık 340–400 MPa√mm деңгээлindedir.
• 25. adım çevresinde yaklaşık 1.100 MPa√mm’ye жүкselmiştir.
• 51. adımda yaklaşık 1.850 MPa√mm’ye ulaşmıştır.

Son баалуулук bildirмененn 3.952,85 MPa√mm сынуу бышыктыгыna ulaşmamaktadır. Bu эмне үчүнle 51. adım nihai сынуу же kopma noktası değildir.

Çatlak ucundaki чыңалуу жана тармак түзүмsı эмнени көрсөтөт?

On sekizinci sayfadaki Şekil 19, 1, 25 жана 51. adımlardaki жарака ачыкlığını жана жарака ucu чыңалуу alanını көрсөтөт. Ağ yalnızca жарака cephesi çevresinde yoğunlaştırılmış, uzak bölgelerde daha kaba bırakılmıştır.

İlk adımda жогорку чыңалуу küçük bir bölgeye yoğunlaşmaktadır. İlerleyen adımlarda жарака açılması görünür hâle gelmekte, жогорку чыңалуу bölgesi genиштетүүkte жана iki жарака ucundaki yoğunlaşma artmaktadır. Görüntülerde yer değiştirmeler 30 kat büyütülmüştür; gösterмененn ачыкlık чыныгы ölçekli deformasyon değildir.

Hangi сынуу modu baskındır?

On sekizinci жана on dokuzuncu sayfalardaki натыйжалар, üç чыңалуу şiddeti faktörü arasında \(K_I\)’in ачык biçimde baskın болгонуnu көрсөтөт. \(K_I\) жарака бетlerinin birbirinden açıldığı Mode I жүктөөyi temsil eder.

• \(K_I\), жарака büyüdükçe düzenli жана күчтүү biçimde artmıştır.
• \(K_{II}\), daha küçük genlikli жана dalgalı kalmıştır.
• \(K_{III}\), işaret değiştirmiş бирок \(K_I\)’den çok daha küçük olmuştur.
• J-интегралы \(K_I\)’e benzer biçimde artmıştır.
• T-чыңалуусу жарака cephesinin büжүк bölümünde negatif kalmıştır.

Bu натыйжа, incelenen жүктөө durumunda çatlтармакın esas катары açılma modu tarafından мененrletildiğini, kayma bмененşenlerinin болсо yerel жана ikincil katkı yaptığını көрсөтөт. Farklı buz temas açıları же жүктөө yönlerinde aynı mod sıralaması garanti edilmemektedir.

Başlangıç жарака узундугу кантип etkili olmuştur?

Şekil катышı \(a/c=0{,}6\) sabit tutularak 2, 3, 5 жана 10 mm başlangıç бет uzunlukları салыштырылган. Bu grup үчүн ortak салыштыруу деңгээлi, en жогорку \(K_I\)’in yaklaşık 1.200 MPa√mm жана en жогорку J-интегралыnin yaklaşık 1,2 × 10⁻⁵ MJ/mm² деңгээлine ulaşmasıdır.

Yirminci sayfadaki Şekil 21’de başlangıç çatlтармакı uzadıkça бардык \(K_I\) жана J eğrмененri sola kaymaktadır. Bu, aynı сынуу sürücü деңгээлine daha az циклde ulaşıldığını көрсөтөт. Etki doğrusal değildir: başlangıç uzunluğunun beş katına çıkması, referans цикл sayısını yaklaşık onda bire indirmiştir.

Başlangıç жарака сүрөт катышı кантип etkili olmuştur?

Başlangıç бет uzunluğu 2 mm sabit tutulmuş жана \(a/c\) катышı 0,5, 0,6, 0,7, 0,8 жана 1,0 катары değiştirilmiştir. Oran büyüdükçe жарака cephesi daha derin жана daha dairesel hâle gelmektedir.

Bu салыштыруу grubunda kullanılan ortak uç деңгээл yaklaşık 1.750 MPa√mm en жогорку \(K_I\)’dir. Bu referans, жарака узундугу изилдөөsındaki 1.200 MPa√mm деңгээлinden ар башкаdır.

Şekil катышı da жараканын өсүшүni hızlandırmıştır; бирок eğrмененr arasındaki өзүнчөm başlangıç uzunluğu değişimindeki kadar büжүк değildir. İncelenen parametre aralığında баштапкы жарака uzunluğu daha baskın өзгөрмөdir.

Изилдөөdaki “эквиваленттүү саат” маанилерi кантип yorumlanmalıdır?

Her 200 saniyelik yapay спектрda 43.525 rainflow циклi bulunmaktadır. Tablolardaki saat маанилерi, gerekli цикл sayısının bu sabit цикл yoğunluğuna bölünmesiyle эсептелген:

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

Bu sürelerin маанилүү yorum sınırları vardır:

• Gerçek bir кемеnin kesintisiz катары aynı buz kesme durumunda çalıştığını varsaymaz; yalnızca matematiksel tekrardır.
• Nihai kopma же \(K_{\mathrm{IC}}\)’ye ulaşma süresi değildir.
• İki parametrik таблицаda kullanılan referans \(K_I\) деңгээлleri aynı değildir.
• Farklı görev profilleri, ачык суу seyri, өзгөрмө buz kalınlığı жана bakım duruşları hesaba katılmamıştır.
• Gerçek takvim saati же коопсуз bakım aralığı катары kullanılamaz.

Yaklaşık 600 MPa√mm деңгээлi эмне маанини билдирет?

Изилдөөчүлөр parametrik eğrмененrde \(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 болсо мененrlemenin çok hızlandığını gözlemlemiştir.

Bu маанилер никель-алюминий коло үчүн evrensel bir “kararlı-kararsız geçiş sınırı” катары yorumlanmamalıdır. Sonuç; kullanılan Paris katsayılarına, жарака geometrisine, seçмененn чыңалуу alanına, сандык мененrleme adımlarına жана bu изилдөөнүн жүк спектрuna özgüdür. Deneysel eşik же standardize edilmiş kabul sınırı değildir.

Изилдөөnın күчтүү yönleri nelerdir?

• Buz, суу жана винт жүкlerini tek bir tam eşleştirilmiş модельde birleştirmektedir.
• Buz сынууsı когезиялуу элементlarla жана parça өзүнчөlmasıyla temsil edilmiştir.
• Kısa süreli чыңалуу маалыматsi түздөн-түз tekrarlanmak yerine чарчоо ölçütleriyle yeniden түзүмlandırılmıştır.
• Paris sürücüsü, rainflow dтармакılımı, sınır sürekliliği жана tekrar çeşitliliği aynı amaç fonksiyonunda колдонулган.
• Ana hiperparametrelerin өзүнчө өзүнчө çıkarıldığı duyarlılık талдоосу сууnulmuştur.
• Yтармакmur akışı гистограммаı, aralık–орточо matrisi, жыштык спектрu жана Paris sürücü катышı бирге баалуулукlendirilmiştir.
• Özgün жана yeniden түзүмlandırılmış жүкlerle hesaplanan жарака uzunlukları түздөн-түз салыштырылган.
• Mode I, Mode II, Mode III, J-интегралы жана T-чыңалуусу жарака cephesi boyunca incelenmiştir.
• Başlangıç жарака узундугу жана сүрөт катышı өзүнчө parametrik gruplarda баалуулукlendirilmiştir.

Изилдөөnın başlıca чектөөлөрү nelerdir?

• Изилдөөda yeni физикалык buz kesme же винт чыңалуу экспериментi түзүмlmamıştır.
• Модельin материал жана buz сынуу parametreleri önceki изилдөөlardan alınmıştır.
• Tek bir R sınıfı винт geometrisi incelenmiştir.
• Tek муз блогу, tek giriş hızı, tek buz hızı жана tek винт devri колдонулган.
• Yalnızca bir калак түбү элементındaki чыңалуу geçmişi yeniden түзүмlandırılmıştır.
• Yeniden түзүмlandırma, kısa жазууta hiç bulunmayan yeni deniz же buz durumlarını üretemez.
• 25–1.000 saniyelik бардык спектрlar aynı kısa сандык жазууtan türetilmiştir.
• Spektrum валидацияsı bтармакımsız bir узак мөөнөттүү чыныгы чыңалуу kaydıyla түзүмlmamıştır.
• Özgün жана yeniden түзүмlandırılmış спектр салыштырууsı aynı сандык модель үчүнde түзүмlmıştır.
• 1,16 saniyedeki tek uzamsal чыңалуу alanının өзгөрмө genlikli жүктөө boyunca temsil edici болгону varsayılmıştır.
• Paris мыйзамыnın geçerli болгону kararlı жарака мененrleme bölgesi dışındaki davranış өзүнчөntılı модельlenmemiştir.
• Çatlak kapanması, цикл sırası etkisi, aşırı жүк geciktirmesi жана artık чыңалууler ачыкça ele alınmamıştır.
• Üretim кемтиктерının чыныгы dтармакılımı же мүмкүнlıksal belirsizliği hesaba katılmamıştır.
• Hesaplanan referans süreler nihai сынуу же чыныгы görev ömrü değildir.
• Çatlak uzunluğu жана сүрөт катышı изилдөөlarında aynı referans \(K_I\) деңгээлi kullanılmamıştır.
• Изилдөө рецензиядан өткөн эмес.

Изилдөө neyi колдойт?

• Kullanılan S-ALE жана когезиялуу элемент düzeni, darbeli buz kesme күчleri жана buz parçalanması üretebilmiştir.
• İncelenen durumda калак түбү жогорку созулуу жана eşбаалуулук чыңалуу bölgesi катары belirlenmiştir.
• Önerмененn blok tabanlı ыкма, özgün kısa kaydın rainflow жана baskın жыштык өзгөчөliklerini yaklaşık катары koruyabilmiştir.
• 200 saniyelik спектрun normalize Paris sürücüsü özgün kaydın yüzde 0,2 üzerindedir.
• Yeniden түзүмlandırılmış жана özgün жүктөөlerin 51 adımlık жарака узундугу натыйжалары yüzde 2,55 ар башкаdır.
• İncelenen жүктөө durumunda Mode I açılma bмененşeni жараканын өсүшүnde baskındır.
• Başlangıç жарака узундугу, seçмененn parametre aralığında сүрөт катышından daha күчтүү etki көрсөткөн.
• Daha büжүк başlangıç кемтиктерı referans сынуу sürücüsüne çok daha az циклde ulaşmıştır.

Изилдөө neyi далилдебейт?

• Yöntemin бардык муз класстуу винттерде aynı doğruluğu sтармакlayacтармакını далилдебейт.
• Модельin чыныгы уюлдук koşullarındaki жүкleri экспериментsel belirsizlik үчүнde doğru tahmin ettiğini göstermez.
• Yeniden түзүмlandırılmış спектрun чыныгы узак мөөнөттүү görev спектрuyla aynı болгонуnu göstermez.
• 79,15 mm жарака узундугуnun чыныгы bir винтde aynı цикл sayısında oluşacтармакını далилдебейт.
• Tablolardaki saatlerin чыныгы seyir же takvim ömrü болгонуnu göstermez.
• 2 mm, 5 mm же 10 mm’nin izin маалыматlebilir bakım sınırı болгонуnu belirlemez.
• Yaklaşık 600 MPa√mm баалуулукinin evrensel hızlanma eşiği болгонуnu далилдебейт.
• Çatlтармакın 51. adımda nihai сынууya ulaştığını göstermez.
• Farklı buz kalınlıkları, kavite, korozyon же өзгөрмө sıcaklık altında aynı натыйжаларыn geçerli болгонуnu göstermez.
• Gerçek bakım aralığını же ekonomik катары en uygun muayene planını hesaplamaz.

Изилдөөnın Yöntemi жана Bulguları

Teknik ыкма özeti

Gerilme спектрu кайра куруу ayarları

Ana geçici жүк bulguları

Spektrum валидация bulguları

Başlangıç жарака geometrisi натыйжаларыnın салыштырууsı

İki yüzde ар башка салыштыруу gruplarına жана ар башка referans сынуу sürücüsü деңгээлlerine aittir. Mutlak saat маанилерi iki таблица arasında түздөн-түз karşılaştırılmamalıdır.

Şekillerin temel mesajları

• Şekil 1: Nikel-alüminyum bronz үчүн önceki изилдөөдөn alınan төмөн жана жогорку циклli S-N eğrisini көрсөтөт.
• Şekil 2: Buz-суу-винт талдоосуnden спектр üretimine жана жарака баалуулукlendirmesine uzanan beş баскычlı ыкмасы көрсөтөт.
• Şekil 3: S-ALE суюктук жана түзүм düğümleri arasındaki ceza күчi eşleştirmesini көрсөтөт.
• Şekil 4: Buz элементları arasındaki kohezyon bölgesinin hasar görerek жарака oluşturmasını көрсөтөт.
• Şekil 5: R sınıfı винт geometrisini жана аба-суу-buz эсептөө аймагыnı көрсөтөт.
• Şekil 6: Buz kesme sırasında −150 менен +230 kN arasında değişen darbeli винт түртүүsini көрсөтөт.
• Şekil 7: Buz bloğunun калак çevresinde ezilmesi жана parçaların сууya dтармакılmasını көрсөтөт.
• Şekil 8: 1,16 saniyedeki von Mболсоs жана en жогорку башкы чыңалуу dтармакılımlarını көрсөтөт.
• Şekil 9: Kritik калак түбү элементını жана 0,482–2,730 saniyelik башкы чыңалуу geçmişini көрсөтөт.
• Şekil 10: Yтармакmur akışı aralık dтармакılımlarının hafıza katsayısıyla birleştirilmesini көрсөтөт.
• Şekil 11: Yeniden түзүмlandırma hiperparametrelerinin çıkarılması же değiştirilmesinin ölçütlere etkisini көрсөтөт.
• Şekil 12–13: Farklı sürelerdeki rainflow гистограммаı жана aralık–орточо matrislerinin ana түзүмsını koruduğunu көрсөтөт.
• Şekil 14–15: Frekans түзүмsını жана beş sürenin normalize Paris sürücüsü катышlarını көрсөтөт.
• Şekil 16: Çatlak hesabında kullanılan 200 saniyelik спектрun цикл sırasını көрсөтөт.
• Şekil 17: Yüzey uzunluğu 2 mm жана derinliği 0,6 mm olan başlangıç çatlтармакını көрсөтөт.
• Şekil 18–20: Çatlak cephesi, цикл birikimi, \(K_I\), \(K_{II}\), \(K_{III}\), J-интегралы жана T-чыңалуусуnin 51 adımlık gelişimini көрсөтөт.
• Şekil 21: Başlangıç жарака узундугу büyüdükçe сынуу sürücüsünün daha az циклde жүкseldiğini көрсөтөт.
• Şekil 22: Daha büжүк \(a/c\) катышlarının мененrlemeyi hızlandırdığını, бирок etkinin uzunluk değişimine göre daha sınırlı болгонуnu көрсөтөт.

Kanıt деңгээлine göre натыйжаларыn durumu

Kaynak жана Yöntem Notu

Изилдөөnı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 жана Chao Wang.

Yazar sıralaması: Xiaoli Yuan birinci, Xinyao Wang ikinci, Weipeng Xiong üçüncü, Mi-An Xue dördüncü, Yanqi Wang beşinci жана Chao Wang altıncı авторdır.

Eşit katkı bilgisi: Yüklenen sürümde eşit katkıda bulunan автор bilgisi орун алган эмес.

Sorumlu автор: Weipeng Xiong.

Sorumlu автор e-postası: xiongweipeng@hhu.edu.cn

Kurumsal байланышlar:

1. College of Harbour, Coastal and Offshore Engineering, Hohai Uniжанаrsity, Nanjing, Jiangсуу 210098, Çin.
2. College of Mechanics and Engineering Science, Hohai Uniжанаrsity, Nanjing, Jiangсуу 211100, Çin.
3. College of Shipbuilding Engineering, Harbin Engineering Uniжанаrsity, Harbin 150001, Çin.

Yazar–мекеме eşleştirmesi:

• Xiaoli Yuan: Hohai Uniжанаrsity bünyesindeki birinci жана ikinci мекемеlar.
• Xinyao Wang: Hohai Uniжанаrsity bünyesindeki birinci мекеме.
• Weipeng Xiong: Hohai Uniжанаrsity bünyesindeki birinci мекеме.
• Mi-An Xue: Hohai Uniжанаrsity bünyesindeki birinci мекеме.
• Yanqi Wang: Hohai Uniжанаrsity bünyesindeki birinci мекеме.
• Chao Wang: Harbin Engineering Uniжанаrsity bünyesindeki üçüncü мекеме.

DOI:10.2139/ssrn.7197210

Yayın platformu: SSRN.

Resmî булак байланышsı:SSRN изилдөө kaydı

Yayın yılı: 2026.

Kesin SSRN жүктөө tarihi: Yüklenen sürümde ачык bir SSRN gönderim tarihi орун алган эмес; kesin gün bilgisi bu текшерүүde bтармакımsız катары doğrulanamamıştır.

Dergi durumu: Sayfa üstbilgisinde “Ocean Engineering xxx (xxxx) xxxxxx” ifadesi bulunmaktadır. Cilt, sayı, makale numarası, kabul tarihi же nihai цикл içi yayın kaydı bulunmadığı үчүн hakemli dergi yayını doğrulanamamıştır.

Yayınevi: Nihai hakemli yayın үчүн kesin bir yayınevi kaydı жок. İncelenen sürüm SSRN üzerinde yayımlanan bir preprinttir.

Kaynak türü: Sayısal суюктук-түзүм өз ара аракеттенүүi, algoritmik чыңалуу спектрu кайра курууsı жана сынуу mekaniği tabanlı parametrik талдоо; preprint.

Hakemlik durumu: Бул изилдөө bir preprinttir жана рецензиядан өткөн эмес.

Yazar katkıları:

• Xiaoli Yuan: Kavramsallaştırma, ыкма жана gözden geçirme-düzenleme.
• Xinyao Wang: İlk taslak, валидация, маалымат düzenleme жана ыкма.
• Weipeng Xiong: Yöntem жана gözden geçirme-düzenleme.
• Mi-An Xue: Yöntem жана gözden geçirme-düzenleme.
• Yanqi Wang: Kavramsallaştırma жана gözden geçirme-düzenleme.
• Chao Wang: Yöntem жана gözden geçirme-düzenleme.

Finansman: Изилдөө; Çin Ulusal Doğa Bilimleri Vakfı 52501389 numaralı fonu, Jiangсуу Eyaleti Doğa Bilimleri Vakfı BK20241523 numaralı fonu, Shanghai Jiao Tong Uniжанаrsity State Key Laboratory of Ocean Engineering GKZD010091 numaralı fonu жана KY10100250185 numaralı temel изилдөө projesi tarafından desteklenmiştir.

Çıkar çatışması: Yazarlar изилдөөyı etkмененyebмененcek bilinen bir mali çıkar же kişболсоl ilişki bildirmemiştir.

Veri erişimi: Verмененrin talep üzerine paylaşılacтармакı белгиленген. Açık маалымат depoсуу байланышsı берилген эмес.

Kod erişimi: Gerilme спектрu кайра куруу algoritmasının kodu үчүн ачык depo байланышsı орун алган эмес.

Etik kurul: Изилдөө insan же hayvan katılımcı içermeyen сандык mühendislik изилдөөsıdır. Yüklenen sürümde etik kurul bilgisi орун алган эмес.

Bu Türkçe ачыкlama yalnızca жүкlenen 23 sayfalık изилдөөнүн metni, denklemleri, таблицаları, сандык модель өзүнчөntıları жана Şekil 1–22’deki görsel натыйжалар temel alınarak hazırlanmıştır. Изилдөөda bulunmayan чыныгы deniz экспериментi, tam ölçekli өлчөө, bakım standardı, izin маалыматlebilir жарака boyutu же коопсуз ишкана süresi eklenmemiştir. Dış булакlar yalnızca bibliyoграфик идентификациясыn, DOI’nin жана расмий SSRN kaydının doğrulanması үчүн колдонулган.

Изилдөөnın en күчтүү yönü, кыска мөөнөттүү buz-суу-винт simülasyonunu жарака sürücüsüne göre seçмененn узак мөөнөттүү bir спектрla жана FRANC3D жараканын өсүшүyle birleştirmesidir. En маанилүү валидация sonucu, özgün жана yeniden түзүмlandırılmış спектрların 51 adım sonunda yüzde 2,55 ар башка жарака узундугу üretmesidir. Муну менен бирге bu салыштыруу aynı сандык ortam үчүнdedir жана чыныгы узак мөөнөттүү чыңалуу өлчөөünün yerini tutmamaktadır.

Başlangıç жарака узундугуnun 2 mm’den 10 mm’ye çıkarılmasıyla bildirмененn yüzde 90,7 azalma, nihai түзүмsal сынуу ömrüne değil, uzunluk parametrik изилдөөsında seçмененn referans \(K_I\) жана J-интегралы деңгээлine ulaşmak үчүн gereken цикл sayısına aittir. Benzer biçimde таблицаlardaki эквиваленттүү саатler чыныгы takvim süresi же bakım aralığı değildir.

Yöntemin mühendislik uygulamasına geçebilmesi үчүн ар башка винт geometrмененri, материал partмененri, başlangıç кемтикu dтармакılımları, buz kalınlıkları, buz hızları, devir sayıları жана görev profilleri altında sınanması; bтармакımsız узак мөөнөттүү жүк жазууларыyla karşılaştırılması жана экспериментsel жарака мененrleme тестleriyle doğrulanması керек.