Hidrolik darbe veya hidrolik impuls sistemleri, basınçlandırılmış bir akışkanın enerjisini çok kısa sürelerde mekanik harekete dönüştürür. Bu mekanizmalar darbeli delme, titreşimli makineler, hidrolik aktüatörler ve yüksek dinamik yük altında çalışan çeşitli endüstriyel sistemlerde kullanılabilir. Ancak sürekli tekrarlanan basınç darbeleri yalnız akışkanı değil; pistonu, mili, yatakları ve sürtünme yüzeylerini de çevrimsel mekanik yüklere maruz bırakır.

N. Smakova ve çalışma arkadaşları bu problemi yalnız hidrolik basınç kontrolü açısından değil, yüzey bütünlüğü ile dinamik sistem davranışını aynı modelde birleştirerek ele alıyor. Araştırmanın merkezinde, düşük frekanslı hidrolik impuls modülünün tekrarlanan darbeleri altında çalışan bir mil bulunuyor. Mil yüzeyinin pürüzlülüğü, işleme koşulları ve mekanik gerilme durumu; sistemin titreşim genliği, sürtünme kayıpları ve çalışma kararlılığıyla birlikte değerlendiriliyor.

Deneysel çalışmada çapı 60 mm, uzunluğu 120 mm olan 30KhGSA alaşımlı çelik silindirik numuneler rotary multi-blade turning, yani döner çok ağızlı tornalama yöntemiyle işlenmiştir. İşleme aralıkları 600–900 rpm dönme hızı, 0,10–0,30 mm/dev ilerleme ve 0,50–1,00 mm kesme derinliği olarak belirlenmiştir. Ölçülen ortalama yüzey pürüzlülüğü \(R_a\), 0,63–1,25 µm aralığında bulunmuştur.

Aynı sistem için piston hareketi, çalışma ve dönüş odası basınçları, sönüm, yay rijitliği ve sürtünme etkilerini içeren dinamik bir matematiksel model geliştirilmiştir. Model MATLAB/Simulink ve Python SciPy ortamlarında dördüncü dereceden Runge–Kutta yöntemi kullanılarak çözülmüş; fiziksel sistem ile sanal model arasındaki veri alışverişi için MQTT, görselleştirme için Node-RED Dashboard ve MATLAB Live Script kullanıldığı belirtilmiştir.

Kaynağın sonuçlarına göre dijital ikiz ve parametre optimizasyonu sonrasında normalize titreşim genliği 1,00'dan 0,78'e düşerek yaklaşık %22 azalmış; tepe basınçlarında %15–18 azalma ve yüzey kalitesinde yaklaşık %20 iyileşme raporlanmıştır. Yazarlar deneysel ve hesaplanan dinamik veriler arasındaki farkın %5'in altında kaldığını belirtmektedir.

Ancak sonuçlar dikkatli yorumlanmalıdır. Makale iki farklı optimum kesme parametresi seti vermekte, hidrolik hareket denkleminde basınç ile kuvvet arasındaki etkin piston alanını açıkça göstermemekte ve Random Forest tabanlı optimizasyon katmanı için ayrıntılı makine öğrenmesi doğrulama metrikleri sunmamaktadır. Dolayısıyla çalışma güçlü bir bütünleşik mühendislik yaklaşımı ortaya koysa da bazı modelleme ayrıntılarının daha açık doğrulanması gerekmektedir.

Hidrolik impuls sistemi ne yapar?

Bir hidrolik impuls mekanizmasında çalışan akışkan kısa süreli basınç değişimleriyle mekanik elemanlara momentum aktarır. Basıncın çok kısa zamanda değişmesi

\[ \frac{dp}{dt} \]

büyüklüğünün yüksek olmasına neden olabilir. Bu basınç dalgalarının yayılma hızı kaynakta

\[ c=\sqrt{\frac{K}{\rho}} \]

ile ifade edilmektedir.

Burada bu denklem özelinde \(K\), akışkanın hacim modülü; \(\rho\) ise yoğunluğudur. Kaynak daha sonraki hareket denkleminde aynı \(K\) sembolünü yay rijitliği için de kullanmaktadır. İki \(K\) fiziksel olarak aynı büyüklük değildir.

Basınç dalgası pistona ulaştığında akışkanın enerjisi pistonun kinetik enerjisine ve bağlı mekanik elemanların deformasyonuna dönüşür. Sistem tekrar tekrar çalıştığında mil, bağlantılar ve sürtünme çiftleri çevrimsel yük altında kalır. Bu nedenle yalnız ortalama basınç değil; darbe genliği, titreşim, sönüm ve malzemenin yorulma davranışı önem kazanır.

Yüzey bütünlüğü neden dinamik davranışı etkileyebilir?

Bir makine parçasının yüzey bütünlüğü yalnız geometrik pürüzlülükten ibaret değildir. Kavram genel olarak yüzey pürüzlülüğü, yüzeye yakın mikroyapı, sertlik veya mikrosertlik, plastik deformasyon ve artık gerilmeler gibi özellikleri kapsar.

Hidrolik impuls modülündeki mil hareket ederken yüzeyler arasındaki gerçek temas koşulları sürtünme kuvvetini etkileyebilir. Daha kaba bir yüzey, temas tepelerinde daha yüksek yerel gerilmeler oluşturabilir; yüzey işleme sırasında meydana gelen mikroyapısal değişiklikler de tekrarlanan yük altında dayanımı etkileyebilir.

Kaynak bu nedenle yüzey işleme teknolojisiyle dinamik sistem davranışını birbirinden bağımsız iki problem olarak görmemektedir. Amaç, mekanik işleme parametrelerini değiştirerek daha uygun bir yüzey oluşturmak ve bunun titreşim ile enerji aktarımına etkisini dijital model üzerinden birlikte değerlendirmektir.

Düşük Frekanslı Hidrolik Darbe Sistemi Nasıl Modelleniyor?

Çalışmada hidrolik impuls modülü, basınç farkının pistonu harekete geçirdiği kütle–yay–sönümleyici tipinde bir mekanik sistem olarak modelleniyor; oda basınçları piston hızına bağlanıyor ve modelin zaman cevabı MATLAB/Simulink ile Python SciPy ortamında sayısal olarak çözümleniyor.

Piston hareket denklemi

Kaynak sistem dinamiğini

\[ M\frac{d^2x}{dt^2} + C\frac{dx}{dt} + Kx = P_A(t)-P_B(t)-F_c \]

denklemiyle ifade etmektedir.

Burada:

• \(M\): hareketli elemanların eşdeğer kütlesi,
• \(C\): viskoz sönüm katsayısı,
• \(K\): bu denklemde geri dönüş yayının rijitliği,
• \(P_A(t)\): çalışma odası basıncı,
• \(P_B(t)\): dönüş odası basıncı,
• \(F_c\): sürtünme ve sızdırmazlık kayıplarını temsil eden terimdir.

Bu denklem kaynakta bu biçimde verilmiştir. Bununla birlikte \(P_A-P_B\) bir basınç farkını, diğer terimler ise kuvvet boyutunu taşımaktadır. Hidrolik silindir mekaniğinde basınç farkının kuvvete dönüşmesi normalde etkin piston alanıyla ilişkilidir. Kaynak denklemde bu alanı açıkça göstermediği ve basınçların önceden kuvvete normalize edilip edilmediğini açıklamadığı için denklemin boyutsal tanımı tam olarak açık değildir.

Basınç ile piston hızı arasındaki bağlantı

Çalışma ve dönüş odası basınçları kaynakta

\[ P_A(t)=P_0+\rho c\frac{dx}{dt} \]

ve

\[ P_B(t)=P_0-\rho c\frac{dx}{dt} \]

olarak tanımlanmıştır.

Bu modelde piston hızı arttıkça iki oda arasındaki basınç farkı da büyür. Başlangıç koşulları ise

\[ x(0)=0, \qquad v(0)=0 \]

ve

\[ P_A(0)=P_B(0)=P_0=10\ \mathrm{MPa} \]

şeklinde verilmiştir.

Tek darbenin enerjisi

Pistonun çarpma anındaki kinetik enerjisi

\[ E=\frac12Mv^2 \]

olarak hesaplanmaktadır. Burada

\[ v=\frac{dx}{dt} \]

çarpma anındaki piston hızıdır.

Bu basit ifade önemli bir tasarım gerilimini gösterir: hız arttığında darbe enerjisi \(v^2\) ile yükselir; fakat daha yüksek hız aynı zamanda basınç darbeleri, titreşim ve mekanik yükler üzerinde olumsuz sonuçlar yaratabilir. Sistem bu nedenle yalnız maksimum darbe enerjisine göre değil, enerji aktarımı ile kararlılık arasında denge gözetilerek optimize edilmektedir.

Dinamik kayıp fonksiyonu

Kaynak dinamik optimizasyon için

\[ J= \min \int_0^T \left[ (P_A-P_B)^2+\alpha v^2 \right]dt \]

ifadesini kullanmaktadır.

\(\alpha\), hızla ilişkili enerji kayıplarının göreli önemini belirleyen ağırlık katsayısı olarak tanımlanır. Basınç farkı ve hız farklı fiziksel birimlere sahip olduğundan, \(\alpha\)'nın aynı zamanda terimlerin ölçeklendirilmesi için gerekli birimsel ağırlığı taşıması gerekir; kaynak bu katsayının sayısal değerini veya birimini vermemektedir.

Dijital İkiz Hidrolik Sistemde Ne İşe Yarıyor?

Kaynakta dijital ikiz, fiziksel hidrolik impuls mekanizmasındaki basınç, sıcaklık ve titreşim sensörlerinden alınan verileri sanal dinamik modele taşıyan; hesaplanan sistem davranışını kontrol ve optimizasyon katmanına geri besleyen gerçek-zamanlı bir fiziksel–dijital eşleştirme yapısı olarak tasarlanmıştır.

Fiziksel sistem ile sanal sistem nasıl haberleşiyor?

Kaynak, fiziksel sistem ile sanal model arasındaki gerçek zamanlı veri değişimi için MQTT protokolünü kullanmaktadır. Basınç, sıcaklık ve titreşim sensörlerinden gelen veriler sanal modele aktarılır.

Veri akışının görselleştirilmesi için Node-RED Dashboard, hesaplanan sonuçların güncellenmesi için ise MATLAB Live Script kullanıldığı belirtilmektedir.

Bu yapı klasik bir çevrimdışı simülasyondan farklıdır. Amaç yalnız geçmiş veriyi analiz etmek değil, ölçülen fiziksel durum ile hesaplanan sanal durum arasında sürekli bağlantı kurmaktır.

Optimizasyon katmanında hangi değişkenler kullanılıyor?

Kaynak, dijital ikizin bilişsel katmanında Gradient Descent ile Random Forest Regressor kombinasyonu kullandığını belirtmektedir. Amaç fonksiyonu

\[ \Phi= w_1\frac{R_a}{R_{a,\mathrm{ref}}} + w_2\frac{P_{\max}}{P_{\mathrm{ref}}} + w_3\frac{A_{\mathrm{vib}}}{A_{\mathrm{ref}}} \rightarrow\min \]

olarak tanımlanmıştır.

Burada:

• \(R_a\): yüzey pürüzlülüğü,
• \(P_{\max}\): maksimum basınç,
• \(A_{\mathrm{vib}}\): titreşim genliği,
• \(w_1,w_2,w_3\): üç hedef arasındaki göreli önemi belirleyen ağırlıklardır.

Bu çok-ölçütlü yapı yalnız en düzgün yüzeyi bulmaya çalışmaz. Çok düşük pürüzlülük elde edilirken basınç tepe değerlerinin veya titreşimin kabul edilemez hale gelmesini önlemek için üç fiziksel kriter aynı amaç fonksiyonunda bir araya getirilmektedir.

Random Forest bölümü ne kadar ayrıntılı?

Makale Random Forest Regressor kullanımını bildiriyor; ancak kullanılan eğitim örneği sayısı, giriş özelliklerinin tam listesi, ağaç sayısı, maksimum derinlik, eğitim/test ayrımı, çapraz doğrulama yöntemi veya \(R^2\), RMSE gibi model performans ölçütleri verilmemektedir.

Bu nedenle Random Forest'ın sistem üzerinde ne kadar ek tahmin gücü sağladığı makaleden nicel olarak belirlenememektedir. Çalışmanın güçlü tarafı bütünleşik dijital ikiz mimarisidir; makine öğrenmesi bileşeninin bağımsız performansının kapsamlı biçimde doğrulandığı söylenemez.

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

Deney numuneleri

Deneysel çalışma 30KhGSA alaşımlı çelikten üretilmiş silindirik mil tipi numuneler üzerinde gerçekleştirilmiştir.

Yüzey pürüzlülüğü

\[ R_a= \frac1L \int_0^L |y(x)|\,dx \]

ifadesiyle tanımlanmıştır. Burada \(y(x)\) yüzey profilinin ortalama çizgiden sapmasını ve \(L\) örnekleme uzunluğunu temsil eder.

Kaynak ölçümlerin beş kez tekrarlandığını, ortalama değerlerin kullanıldığını ve standart sapma ile güven aralıkları üzerinden istatistiksel değerlendirme yapıldığını belirtmektedir. Ancak makalenin ana sonuç sayfalarında bu güven aralıklarının ayrıntılı tabloları sunulmamaktadır.

Optimum İşleme Parametreleri Hangileri?

Makale tek bir tutarlı optimum parametre üçlüsü vermiyor: yöntem bölümündeki matematiksel optimizasyon \(860\) rpm, \(0,21\) mm/dev ve \(0,70\) mm değerlerini bildirirken sonuçlar bölümü deneysel optimumu \(870\) rpm, \(0,23\) mm/dev ve \(0,75\) mm olarak yazıyor; iki set arasındaki ilişkinin ne olduğu kaynakta açıkça açıklanmıyor.

Bu değerlerin biri model tahmini, diğeri deneysel optimum olabilir; ancak makale bunu açık biçimde söylemediği için iki sonucu birbiriyle eşdeğer kabul etmek veya bunlardan birini hatalı ilan etmek doğru değildir.

Yüzey pürüzlülüğü nasıl değişiyor?

Kaynağın 5. sayfasındaki Şekil 2, \(R_a\)'nın kesme hızı ve ilerlemeyle değişimini üç boyutlu yüzey olarak göstermektedir. Grafiğin genel eğilimi daha yüksek dönme hızları ve daha düşük ilerleme değerlerinde daha küçük pürüzlülük değerleri elde edildiğini göstermektedir. :contentReference[oaicite:9]{index=9}

Sonuç bölümünde en iyi koşullarda yaklaşık

\[ R_a=0.63\ \mu\mathrm{m} \]

yüzey pürüzlülüğüne ulaşıldığı bildirilmektedir.

Basınç salınımları nasıl değişiyor?

Kaynağın 6. sayfasındaki Şekil 3, klasik kontrol ve dijital ikiz kontrollü sistemin impuls odası basıncını zamana göre karşılaştırmaktadır. Klasik sistemde başlangıçtaki salınımlar daha büyükken dijital ikiz kontrollü eğride genliğin daha hızlı azaldığı görülmektedir. :contentReference[oaicite:10]{index=10}

Yazarların genel değerlendirmesine göre dijital ikiz kullanımı tepe basınçlarını yaklaşık

\[ 15\%-18\% \]

azaltmaktadır.

Titreşimdeki %22 azalma

Şekil 4'te normalize edilmiş titreşim genliği optimizasyon öncesi

\[ A_{\mathrm{vib}}=1.00 \]

iken optimizasyon sonrasında

\[ A_{\mathrm{vib}}=0.78 \]

olarak gösterilmektedir.

Göreli değişim

\[ \frac{1.00-0.78}{1.00}\times100 = 22\% \]

olduğundan kaynakta bildirilen yaklaşık %22 azalma grafikle doğrudan uyumludur. :contentReference[oaicite:11]{index=11}

Yazarlar bu azalmanın dinamik yükleri düşürerek mekanik elemanların kullanım ömrünü iyileştirebileceğini savunmaktadır. Ancak makale uzun süreli yorulma ömrünü saat veya çevrim sayısı cinsinden ölçen bir dayanım deneyi sunmamaktadır; “ömrün uzaması” doğrudan yaşam süresi testi değil, azaltılmış titreşim ve yüklerden yapılan mühendislik yorumudur.

Mil üzerindeki gerilme dağılımı

Kaynağın Şekil 5'inde mil boyunca hesaplanan gerilme dağılımı renk haritası olarak gösterilmektedir. En yüksek gerilmeler bağlantı ve geçiş bölgelerinde yoğunlaşmaktadır. Kaynak bunun ANSYS Workbench hesaplarıyla uyumlu olduğunu belirtmektedir. :contentReference[oaicite:12]{index=12}

Bununla birlikte makale bu sonlu elemanlar analizinin ağ yoğunluğu, sınır koşulları, temas modeli veya mesh-independence testi gibi ayrıntılarını vermemektedir. Bu nedenle şekil gerilmenin mekânsal eğilimini destekler; ancak bağımsız bir ayrıntılı FEA doğrulama çalışması değildir.

Deney ve model uyumu

Çalışmanın sonuç bölümünde hesaplanan ve deneysel dinamik sonuçlar arasındaki sapmanın

\[ <5\% \]

olduğu ifade edilmektedir.

Bu değer modelin deneysel davranışı iyi temsil ettiğine ilişkin önemli bir nicel iddiadır. Ancak ana metinde bu %5'in hangi değişkenlerin, kaç veri noktasının ve hangi hata metriğinin ortalaması olduğu ayrıntılı olarak tanımlanmamıştır. Bu nedenle sonuç “kaynakta bildirilen %5'ten düşük uyum hatası” olarak aktarılmalı, ayrıntılı RMSE veya güven aralığı gibi yorumlara genişletilmemelidir.

Kaynağın raporladığı ana sonuçlar

Çalışmanın desteklediği sonuçlar

• 30KhGSA mil numunelerinde rotary multi-blade turning parametreleri yüzey pürüzlülüğünü belirgin biçimde etkilemektedir.
• İncelenen aralıkta \(R_a\) değerleri 0,63–1,25 µm arasında raporlanmıştır.
• Hidrolik impuls sisteminin basınç ve piston dinamiği MATLAB/Simulink ve Python SciPy ile modellenmiştir.
• Dijital ikizli kontrol senaryosu kaynak simülasyonlarında basınç salınımlarını klasik kontrole göre daha hızlı bastırmaktadır.
• Normalize titreşim genliği 1,00'dan 0,78'e düşmüştür.
• Kaynak tepe basıncında %15–18 ve yüzey kalitesinde yaklaşık %20 iyileşme raporlamaktadır.
• MQTT, Node-RED ve MATLAB Live Script fiziksel–dijital veri akışı için bir mimari içinde kullanılmıştır.
• Yüzey işleme optimizasyonu ile hidrolik sistem kontrolünün aynı dijital ikiz çerçevesinde ele alınabileceği gösterilmektedir.

Çalışmanın kanıtlamadığı veya tam açıklamadığı noktalar

• İki farklı optimum işleme parametresi setinden hangisinin nihai optimum olduğu kaynakta net değildir.
• Hareket denklemindeki basınç farkının kuvvete dönüşümü için etkin piston alanı açıkça tanımlanmamıştır.
• Random Forest Regressor için eğitim büyüklüğü, hiperparametreler ve bağımsız doğrulama performansı verilmemiştir.
• “Self-learning” dijital ikizin uzun süreli çevrimiçi öğrenme başarısı ayrı bir kontrollü deneyle gösterilmemiştir.
• %20 yüzey iyileşmesinin başlangıç ve son değerlerini veren ayrıntılı veri tablosu sunulmamıştır.
• %5'ten düşük model sapmasının hangi hata metriğiyle hesaplandığı açıklanmamıştır.
• Titreşim azalmasının gerçek hizmet ömrünü ne kadar uzattığı uzun süreli yorulma deneyiyle ölçülmemiştir.
• Mil gerilme analizinin FEA mesh, sınır koşulu ve yakınsama ayrıntıları verilmemiştir.
• Sonuçlar bütün hidrolik darbe sistemlerine veya farklı malzemelere doğrudan genellenemez.

Kaynak ve Yöntem Notu

Tam özgün başlık:Dynamic Modeling and Surface Integrity Optimization of Low-Frequency Hydraulic Impulse Systems.

Yazarlar: N. Smakova, S. Pankov, B. Baisadykov, V. Zelenkov, D. Toibazarov ve A. A. Karypov.

Kurumlar: Kazakh University of Technology and Business named after K. Kulazhanov; National Defense University of the Republic of Kazakhstan; Scientific and Technical Center Almaty Special Design Bureau “Alatau”.

Dergi: Eurasian Physical Technical Journal.

Bölüm: Engineering.

Cilt / sayı / sayfalar: 23, No. 1(55), 113–120.

Yayın tarihi: 30 Mart 2026.

DOI: Yüklenen makale PDF'sinde 10.31489/2026N1/113-120 yazmaktadır. Derginin güncel web metadata'sında aynı makale için 10.31489/2025N2/113-120 bağlantısı da kullanılmaktadır. Bu Verianla içeriğinde bilimsel kaynağın kendi PDF'sinde basılı DOI korunmuştur ve metadata uyuşmazlığı ayrıca kaydedilmiştir.

Hakemlik: Hakemli araştırma makalesidir. Dergi normal araştırma makalelerinde en az iki bağımsız değerlendiricinin kullanıldığı single-blind peer-review süreci uyguladığını bildirmektedir.

Lisans: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International, CC BY-NC-ND 4.0.

Çalışma türü: Deneysel mekanik işleme + hidrolik dinamik modelleme + sayısal simülasyon + dijital ikiz ve optimizasyon.

Numune: 30KhGSA alaşımlı çelik, 60 mm çap ve 120 mm uzunlukta silindirik mil tipi numuneler.

Sayısal araçlar: MATLAB/Simulink, Python SciPy ve dördüncü dereceden Runge–Kutta integrasyonu. Kaynak ayrıca gerilme sonuçlarının ANSYS Workbench hesaplarıyla uyumlu olduğunu belirtmektedir.

Dijital ikiz altyapısı: MQTT veri aktarımı, Node-RED Dashboard görselleştirmesi ve MATLAB Live Script güncellemesi.

Optimizasyon: Kaynak Gradient Descent ve Random Forest Regressor kombinasyonu bildirmektedir. Ancak makine öğrenmesi eğitim ve doğrulama ayrıntıları ana metinde sınırlıdır.

Finansman: Araştırma, Kazakistan Cumhuriyeti Bilim ve Yükseköğretim Bakanlığı Bilim Komitesi tarafından finanse edilen IRN AR196079/0222 numaralı proje kapsamında gerçekleştirilmiştir.

Çıkar çatışması: Yazarlar finansal, kişisel, yazarlık veya araştırma sonuçlarını etkileyebilecek başka bir çıkar çatışması bulunmadığını beyan etmektedir.

CRediT katkıları: N. Smakova: kavramsallaştırma, metodoloji, ilk taslak ve gözetim; S. Pankov: araştırma ve veri kürasyonu; B. Baisadykov: deneysel doğrulama; V. Zelenkov: yazılım ve modelleme; D. Toibazarov: görselleştirme ve veri analizi; A. A. Karypov: kaynaklar ve teknik destek.

Kaynak içi tutarlılık notu: Bölüm 2.4 optimum işleme parametrelerini 860 rpm, 0,21 mm/dev ve 0,70 mm olarak verirken sonuçlar bölümü 870 rpm, 0,23 mm/dev ve 0,75 mm bildirmektedir. Ayrıca aynı K sembolü hem akışkan hacim modülü hem yay rijitliği için kullanılmıştır ve hareket denklemindeki basınç farkının etkin piston alanıyla kuvvete çevrilmesi açıkça gösterilmemiştir.

Verianla Live: Yüksek derecede uygundur. En anlamlı canlı görselleştirme; fiziksel hidrolik modül ile dijital ikiz arasındaki MQTT veri akışı, impuls odasındaki \(P_A\) ve \(P_B\) basınçlarının zamanla değişimi, piston \(x(t)\) ve \(v(t)\) hareketi, ardından klasik kontrol ile dijital ikiz kontrolünün titreşim ve basınç sönümlemesinin karşılaştırılmasıdır. Ham deney noktaları olmadan kaynak grafiklerinden yeni sayısal veri türetilmemelidir.