Bu araştırma, metal sıvama işleminde sonlu eleman simülasyonundan elde edilen merdane hareket yörüngesinin gerçek sıvama makinesine aktarılırken oluşan sistematik sapmayı ölçmek ve azaltmak için bir yöntem geliştirmektedir. ABAQUS simülasyonunda analitik rijit gövde olarak tanımlanan merdaneye 1 × 1 × 1 mm boyutunda ayrık bir kabuk bağlanmış, bu kabuğun düğüm koordinatları Python ile kare kare çıkarılarak merdane merkez yörüngesi hesaplanmıştır. Simülasyon yörüngesi gerçek makine koordinat sistemine dönüştürülmüş; tek merdaneli sıvama makinesinin gerçek hareketi çekme telli deplasman sensörüyle ölçülmüştür. Üç telafisiz deneyde ortalama projeksiyon hatası 0,503 mm iken projeksiyon yöntemi ve Gauss Süreci Regresyonu ile oluşturulan normal-yön telafisinden sonra yeni makine yürütmesindeki ortalama hata 0,229 mm olarak raporlanmıştır. Kaynak bu değişimi %54,47 hata azalması olarak vermektedir. Bununla birlikte çalışma, telafi yönteminin bitmiş sıvama parçasının cidar kalınlığı, yüzey kalitesi veya boyutsal doğruluğunu aynı oranda geliştirdiğini göstermemektedir; deneysel doğrulama esas olarak merdane hareket yörüngesine yöneliktir.

Araştırmanın odak noktası klasik anlamda “en iyi sıvama yörüngesini” bulmak değildir. Çalışma, sonlu eleman modelinde tanımlanan bir yörüngenin gerçek CNC/makine hareketine geçirildiğinde neden tam olarak aynı geometrik yolu izlemediği ve bu simülasyon–makine farkının nasıl telafi edilebileceği sorusunu ele almaktadır. Yazarlar koordinat bileşenlerindeki X veya Z farkını ayrı ayrı modellemek yerine, ölçülen her noktanın ideal eğri üzerindeki en yakın projeksiyonunu bulup hatayı eğrinin yerel normal doğrultusunda tanımlamaktadır.

Gauss Süreci Regresyonu, bu normal-yön hatasının yörünge boyunca nasıl değiştiğini modellemek için kullanılmıştır. Üç telafisiz makine yürütmesinden elde edilen 370 projeksiyon örneği birleştirilmiş, yay uzunluğu boyunca 1,0 mm genişlikli kutular içinde ortalama alınarak 102 eğitim örneğine düşürülmüştür. GPR modelinde Matérn çekirdeği (\(\nu=1{,}5\)), sabit çekirdek ve beyaz gürültü terimi birlikte kullanılmış; modelin tahmin ettiği sistematik hata, simülasyon yörüngesine ters normal yönde uygulanarak telafi edilmiş yeni yörünge oluşturulmuştur.

Telafi sonucunun istatistiksel analizi de yalnız ortalama değere dayanmamaktadır. Kaynak, telafi öncesi ve sonrası mutlak projeksiyon hata dağılımlarını Welch t-testi ve Mann–Whitney U testiyle karşılaştırmıştır. Bildirilen p değerleri sırasıyla \(2{,}442\times10^{-4}\) ve \(2{,}040\times10^{-5}\) olup çalışmada kullanılan \(\alpha=0{,}05\) anlamlılık eşiğinin altındadır. Bu sonuç, incelenen makine ve yörünge koşulunda ölçülen hata azalmasının deneysel veri içinde rastlantısal değişimle açıklanmasının güç olduğunu desteklemektedir; farklı makineler veya farklı sıvama koşulları için aynı hata azalması henüz gösterilmemiştir.

Türkiye açısından çalışma, sonlu eleman simülasyonunun yalnız tasarım aracı olarak kalmayıp gerçek CNC/metal şekillendirme ekipmanına aktarılmasında ölçüm-temelli bir kapalı doğrulama zinciri kurulabileceğini göstermesi bakımından önemlidir. Türkiye'deki metal sıvama ve hassas şekillendirme uygulamalarında benzer bir yaklaşım, simülasyon yörüngesi ile gerçek eksen hareketi arasındaki sistematik sapmayı belirlemek için uyarlanabilir. Ancak mevcut çalışma tek bir kendi geliştirdikleri sıvama makinesi, tek bir çift konik yörünge ailesi ve sınırlı deney tekrarıyla yürütülmüştür; ayrıca yörünge uydurma deneyinde iş parçası kullanılmadığı açıkça belirtilmektedir. Bu nedenle 0,229 mm sonucu farklı makinelere veya nihai parça toleranslarına doğrudan aktarılmamalıdır.

Araştırmanın temel problemi nedir?

Sonlu eleman analizinde merdane hareketi matematiksel olarak hassas bir eğri üzerinden tanımlanabilir. Gerçek makine ise bu yörüngeyi CNC interpolasyonu, koordinat dönüşümleri, aktarma mekanizmaları ve fiziksel eksen hareketleri üzerinden uygular. Bu nedenle simülasyonda oluşturulan geometrik yol ile gerçek merdane hareketi arasında sistematik farklar oluşabilir.

Araştırmacıların hedeflediği literatür boşluğu, metal sıvamada yörüngenin şeklinin optimize edilmesinden farklıdır. Çalışmanın temel sorusu, sonlu eleman simülasyonundaki merdane merkez yörüngesinin gerçek makineye nasıl güvenilir biçimde taşınacağı ve makinenin uyguladığı yörünge ile simülasyon yörüngesi arasındaki yerel geometrik hatanın nasıl telafi edileceğidir.

Kaynak özellikle klasik koordinat-bileşeni hata yöntemlerinin sürekli eğrisel yörüngelerde yerel geometrik sapmayı tam olarak temsil etmeyebileceğini savunmaktadır. Bu nedenle hata X ve Z eksenlerinde iki ayrı koordinat farkı yerine eğrinin yerel normal doğrultusunda tek bir işaretli geometrik hata olarak tanımlanmıştır.

Sonlu eleman modeli hangi parçayı temsil ediyor?

ABAQUS 2022 ortamında çift konik bir parçanın tek merdaneli sıvama modeli oluşturulmuştur. Model üç temel bileşenden oluşmaktadır:

• mandrel/çekirdek kalıp,
• 6061 alüminyum alaşımlı halka biçimli taslak,
• α tipi sıvama merdanesi.

Çekirdek kalıbın küçük uç çapı 30 mm, büyük uç çapı 160 mm, ilk yarı koni açısı 30°, ikinci yarı koni açısı 45°, toplam yüksekliği 78,39 mm ve iki konik bölüm arasındaki geçiş yayı yarıçapı 50 mm'dir.

Halka biçimli taslağın iç çapı 15 mm, dış çapı 150 mm ve kalınlığı 2,2 mm'dir.

Merdane geometrisi için ilgili model bölümünde çalışma açısı 30°, köşe yarıçapı 5 mm ve merdane çapı 200 mm olarak verilmektedir.

Sonlu eleman ağ yapısı nasıl kurulmuştur?

6061 alüminyum taslak C3D8R sekiz düğümlü azaltılmış integrasyonlu hegzagonal elemanlarla ağlanmıştır. Nihai model:

• 5850 eleman,
• 11.960 düğüm

içermektedir. Taslağın dış yüzeyindeki ağ yoğunluğunu geliştirmek amacıyla adaptif ağ yaklaşımı uygulanmıştır. Merdane ve mandrel analitik rijit gövde kabul edildiğinden bu iki bileşen için sonlu eleman ağı oluşturulmamıştır.

Temas koşulları nasıl tanımlanmıştır?

İki temel yüzey-yüzey temas çifti oluşturulmuştur: merdane–taslak ve mandrel–taslak. Normal doğrultuda “hard contact”, teğetsel yönde penalty sürtünme modeli kullanılmıştır.

Mandrel sabit açısal hızla dönerken merdane önceden tanımlanmış yörünge boyunca radyal dalma ve eksenel ilerleme hareketi gerçekleştirmektedir.

6061 alüminyum için hangi malzeme değerleri kullanılmıştır?

Simülasyon parametrelerinde neden bir kaynak tutarsızlığı var?

Simülasyon parametrelerini açıklayan metin, merdane kurulum açısını 30°, köşe yarıçapını 5 mm, mil hızını 120 r/min ve ilerleme oranını 0,4 mm/r olarak tanımlamaktadır. Aynı bölümdeki Tablo 3 ise ilk iki parametreyi sırasıyla 25° ve 6 mm vermektedir.

Kaynak bu iki farklı açı/yarıçap setinden hangisinin nihai simülasyonda kullanıldığını ayrıca açıklamamaktadır. Bu nedenle Verianla makalesinde herhangi bir set “doğru değer” olarak sessizce seçilmemiştir.

Analitik rijit merdanenin merkezi nasıl çıkarılmıştır?

ABAQUS'ta analitik rijit gövde olarak tanımlanan merdanenin geometrik merkezinin sürekli hareket yörüngesi standart çıktıdan doğrudan alınamadığı için araştırmacılar yardımcı bir 1 × 1 × 1 mm ayrık kabuk kullanmıştır.

Bu küçük kabuk merdaneye tie constraint ile bağlanmış ve merdaneyle rijit biçimde birlikte hareket ettirilmiştir. Kabuk üzerindeki karakteristik bir köşe düğümünün koordinatları her analiz karesinde kaydedilmiştir.

COORD alan çıktısı ABAQUS ODB veritabanından Python 3.13 betiğiyle kare kare okunmuş; düğüm numarası ve üç boyutlu koordinatlar metin dosyasına yazılmıştır.

Karakteristik düğüm:

\[ A_i=(x_{A_i},y_{A_i},z_{A_i}) \]

ve merdane merkezi:

\[ P_i=(x_{P_i},y_{P_i},z_{P_i}) \]

olarak tanımlandığında kaynak temel ilişkiyi:

\[ P_i=A_i+(\Delta x,\Delta y,\Delta z) \]

biçiminde vermektedir.

Basılı geometrik ofsetler:

\[ \Delta x=\frac{a}{2}+\frac{d}{2}\cos\theta \]

\[ \Delta y=\frac{a}{2} \]

\[ \Delta z=\frac{d}{2}\sin\theta+h \]

şeklindedir. Burada \(a\) ayrık kabuğun kenar uzunluğu, \(d\) merdane kalınlığı, \(\theta\) merdane kurulum açısı ve \(h\) ek geometrik ofsettir.

Kaynak \(a=1\) mm, \(d=30\) mm ve \(\theta=30^\circ\) değerlerini verdikten sonra merdane merkezi için \(P_1=(x_A+13{,}433,\ y_A-22{,}267,\ z_A+0{,}5)\) biçiminde bir dönüşüm yazmaktadır. Bu sayısal ofset basılı Denklem (2)–(7)'den tüm eksenlerde doğrudan yeniden üretilemediği için kaynak içi koordinat dönüşümü belirsizliği olarak değerlendirilmelidir.

Simülasyon yörüngesi gerçek makineye nasıl aktarılmıştır?

Çıkarılan merdane merkez koordinatları önce ABAQUS koordinat sisteminden makine koordinat sistemine dönüştürülmüştür. Araştırmacılar bunun koordinat başlangıçlarının hizalanması ve makine geometrisine uygun uzaysal konum düzenlemeleriyle yapıldığını belirtmektedir.

Dönüştürülmüş yörüngenin Tablo 5'te 401 veri noktası bulunduğu gösterilmektedir. Ardından CNC makinesinin artımsal interpolasyon gereksinimi nedeniyle komşu noktalar arasındaki hareketler birinci dereceden fark alma yöntemiyle işlenmiştir.

Burada birim gösterimine dikkat edilmelidir. Şekil 10 eksenleri “m” olarak etiketlenirken sonraki fiziksel yörünge tablolarında ve grafiklerde X–Z koordinatları mm cinsinden gösterilmektedir. Kaynak iki gösterim arasındaki ölçek dönüşümünü ayrıntılı biçimde açıklamamaktadır.

Gerçek merdane yörüngesi nasıl ölçülmüştür?

Araştırmacılar makinenin kendi dahili sensörlerine bağımlı kalmamak amacıyla harici ve müdahalesiz bir ölçüm çözümü kullanmıştır. Seçilen mutlak çekme telli deplasman sensörü RS485 çıkışına sahiptir.

Kaynak “1024” değerinin tam fiziksel veya sayısal tanımını vermediğinden bu değer 0,098 mm çözünürlükle kaynak dışı bir formülle uzlaştırılmamıştır.

Sanal-gerçek eşleme sistemi ne yapıyor?

Gerçek zamanlı yörünge izleme için Unity3D 2022.3019f1c1 tabanlı bir sanal-gerçek eşleme sistemi geliştirilmiştir. Sensör verileri Modbus-RTU protokolü üzerinden alınmış ve XChart ile görüntülenmiştir.

Arayüzde sonlu eleman modelinden gelen hedef yörünge ile gerçek makine merdanesinin ölçülen yörüngesi aynı grafik üzerinde gösterilmektedir. 13–16. sayfalardaki görseller, bilgisayar arayüzüyle gerçek makine eksen konumlarının birlikte kontrol edildiğini göstermektedir.

Yörünge deneyi gerçek bir parça üzerinde mi yapılmıştır?

Hayır; kaynak yörünge uydurma deneyinin iş parçası olmadan gerçekleştirildiğini açıkça belirtmektedir. Çalışmanın amacı bu aşamada merdanenin simülasyon yörüngesini makinede takip edip edemediğini ölçmektir. İlgili görsellerde bulunan kırmızı çerçeve de parçayı değil ölçek referansını temsil etmektedir.

Bu nedenle deney sonucunun bilimsel anlamı “gerçek makinedeki eksen/merdane hareketi simülasyondaki geometrik yola ne kadar yaklaşıyor?” sorusuyla sınırlıdır. Nihai şekillendirilmiş parçanın duvar kalınlığı, yüzey pürüzlülüğü veya geometrik toleransı bu deneyde sonuç değişkeni değildir.

Telafi öncesi üç deney ne gösterdi?

Makine yörüngesi üç kez yürütülmüş ve gerçek yörünge ile simülasyon yörüngesinin özellikle konik bölgede küçük fakat sistematik farklar gösterdiği görülmüştür. Gerçek hareket grafiklerde basamaklı, simülasyon yörüngesi ise düzgün bir eğri biçimindedir.

Üç deneyin ortalama değerlerinin birbirine yakın fakat aynı olmaması, yörünge sapmasının tekrarlanabilir bir sistematik bileşen yanında deneyden deneye değişen bileşenler de içerdiğini göstermektedir.

Projeksiyon yöntemi neden kullanılmıştır?

Gerçek makine yörüngesi basamaklı olduğu için aynı X koordinatındaki basit düşey farkı almak veya X ve Z hatalarını bağımsız değerlendirmek eğrinin yerel geometrisini tam olarak temsil etmeyebilir. Bu nedenle her ölçülen nokta için simülasyon eğrisi üzerindeki en yakın nokta hesaplanmıştır.

Ölçülen nokta:

\[ P_i=(x_i,y_i) \]

ve simüle edilen parametrik eğri:

\[ C(u)=(x(u),y(u)) \]

olarak tanımlandığında kareli uzaklık:

\[ D_i(u)=\|P_i-C(u)\|^2 \]

şeklindedir.

En yakın nokta koşulu:

\[ (x_i-x(u))x'(u)+(y_i-y(u))y'(u)=0 \]

eşitliğiyle elde edilmiş ve doğrusal olmayan denklem Newton yöntemiyle iteratif olarak çözülmüştür:

\[ u_{k+1}=u_k-\frac{f(u_k)}{f'(u_k)} \]

Bu işlem her ölçüm noktasının ideal eğri üzerindeki projeksiyon noktasını belirlemektedir.

Normal-yön hatası nasıl tanımlanmıştır?

Projeksiyon noktasına göre hata vektörü:

\[ \mathbf e=P-C(u^*) \]

olarak tanımlanmaktadır.

Eğrinin birim teğet vektörü:

\[ \mathbf t(u)= \frac{C'(u)}{\|C'(u)\|} \]

ve iki boyutlu birim normal vektörü:

\[ \mathbf n(u)=(-t_y,t_x) \]

olduğunda işaretli normal hata:

\[ e(u)=\left(P-C(u)\right)\cdot\mathbf n(u) \]

şeklinde hesaplanmaktadır.

Bu tanım, hatanın eğriye teğet bileşeninden ziyade ideal yörüngeden dışarı veya içeri doğru geometrik sapmasını modellemeyi amaçlamaktadır.

GPR veri kümesi nasıl oluşturuldu?

Üç telafisiz deneydeki bütün projeksiyon hata noktaları birleştirilmiştir. Kaynakta başlangıç veri kümesi:

370 projeksiyon örneği

olarak verilmektedir.

Gerçek ölçümdeki basamaklı yapı ve tekrarlı konumların GPR üzerindeki etkisini azaltmak amacıyla yay uzunluğu 1,0 mm genişlikli aralıklara bölünmüş ve her aralıktaki hata değerlerinin ortalaması alınmıştır. Bu işlemden sonra:

102 GPR eğitim örneği

kalmıştır.

Kaynak, veri sayısı sınırlı olduğu gerekçesiyle eğitim/test ayrımı yapmadığını açıkça belirtmektedir. Üç deneyin tamamı tek birleşik GPR hata modelinin eğitiminde kullanılmıştır.

GPR hangi çekirdeği kullanıyor?

Yörünge hatasının genel olarak sürekli fakat yerel düzensizlikler içerebildiği varsayımı nedeniyle Matérn çekirdeği kullanılmıştır. Düzgünlük parametresi:

\[ \nu=1{,}5 \]

olarak seçilmiştir.

Kaynakta kullanılan Matérn çekirdeği:

\[ k(u,u')= \sigma^2 \left( 1+\frac{\sqrt3 r}{l} \right) \exp \left( -\frac{\sqrt3 r}{l} \right) \]

ve:

\[ r=|u-u'| \]

şeklindedir.

Burada \(\sigma^2\) sinyal varyansını, \(l\) karakteristik uzunluk ölçeğini temsil eder.

Ölçüm gürültüsünü hesaba katmak amacıyla ayrıca beyaz gürültü terimi eklenmiştir:

\[ k_{total}(u,u')= k_{matern}(u,u') + \sigma_n^2\delta(u,u') \]

Uygulama, sabit çekirdek + Matérn çekirdeği + beyaz gürültü çekirdeği bileşiminden oluşmaktadır. Sinyal varyansı, uzunluk ölçeği ve gürültü seviyesi log marjinal olabilirlik maksimize edilerek optimize edilmiştir. Bununla birlikte kaynak nihai optimize edilmiş hiperparametrelerin sayısal değerlerini ana metinde raporlamamaktadır.

Telafi yörüngesi nasıl üretilmiştir?

GPR'nin belirli bir yay uzunluğundaki tahmini ortalama normal hatası \(\hat e(u)\) olarak alındığında telafi edilmiş yörünge:

\[ C_{comp}(u) = C(u) - \hat e(u)\mathbf n(u) \]

şeklinde oluşturulmuştur.

Yani simülasyon yörüngesi, makinenin sistematik olarak sapma eğilimi gösterdiği yönün tersine ve eğri normaline paralel biçimde kaydırılmıştır.

Şekil 24'te yörünge boyunca telafi miktarı değişken bir eğri olarak gösterilmekte ve grafikte ortalama telafi büyüklüğü 0,3873 mm olarak işaretlenmektedir. Bu değer yörüngenin her noktasına sabit bir 0,3873 mm ofset uygulandığı anlamına gelmez; telafi konuma göre değişmektedir.

Telafi sonrası hata ne kadar azaldı?

Telafi edilmiş yörünge makineye tekrar uygulanmış ve yeni gerçek hareket yörüngesi ölçülmüştür. Kaynak, telafi sonrası ortalama sapmayı 0,229 mm olarak bildirmektedir. Telafi öncesi üç deneyin birleşik ortalaması 0,503 mm'dir.

İstatistiksel test ne gösterdi?

Telafi öncesi veri kümesi üç telafisiz deneydeki bütün mutlak projeksiyon hatalarından, telafi sonrası veri kümesi ise telafi edilmiş yürütmedeki projeksiyon hatalarından oluşturulmuştur.

Kaynak iki ayrı test kullanmıştır:

Her iki p değeri de 0,05'in altındadır. Araştırmacılar bu sonucu telafi sonrasında ölçülen hata azalmasının kendi deney veri kümesinde istatistiksel olarak anlamlı olması şeklinde yorumlamaktadır.

Bu analiz yöntemin farklı makine veya farklı parça üzerinde genellenebilirliğini test etmez. İstatistiksel anlamlılık, mevcut deney düzeninde telafi öncesi ve sonrası hata dağılımlarının farklı olduğunu desteklemektedir.

İstatistiksel karşılaştırmada hangi referans kullanılmıştır?

Kaynak burada dikkat gerektiren bir ifade kullanmaktadır. Bölüm 4.2'nin ilk kısmı 0,229 mm değerini telafi edilmiş gerçek yörünge ile telafi edilmiş simülasyon yörüngesi arasındaki ortalama sapma olarak tanımlar. Devamındaki istatistik bölümünde ise hem telafisiz hem telafi edilmiş ölçümlerin “original finite element simulated trajectory” üzerine projekte edildiği yazmaktadır.

Kaynak bu iki değerlendirme tanımının nasıl uzlaştırıldığını ayrıntılı olarak açıklamamaktadır. Bu nedenle 0,229 mm ve p-değerleri kaynaktaki kendi tanımlarıyla korunmuş, referans yörüngesinin sessizce yeniden yorumlanmasından kaçınılmıştır.

GPR modelinin güçlü yanı nedir?

Çalışmanın önerdiği yaklaşım, hata için önceden sabit bir polinom biçimi dayatmamaktadır. GPR, yörünge konumuna göre sistematik normal hatanın ortalama değerini ve tahmin belirsizliğini birlikte üretebilmektedir.

Şekil 22'de kırmızı noktalar eğitim örneklerini, mavi eğri GPR tahmin ortalamasını, gölgeli alan ise %95 güven aralığı olarak etiketlenen tahmin bandını göstermektedir. Bu görsel, hata modelinin yörünge boyunca sabit bir ofset değil, yerel olarak değişen bir fonksiyon olduğunu açık biçimde göstermektedir.

GPR modelinin sınırı nedir?

Modelin eğitimi için yalnız üç telafisiz makine yürütmesinden elde edilen veriler kullanılmıştır. Kaynakta bağımsız eğitim/test ayrımı veya farklı çalışma koşullarına yönelik cross-validation bulunmamaktadır.

370 ham projeksiyon örneğinin aynı üç yürütmeden gelmesi, örnek sayısının 370 olmasının 370 bağımsız deney anlamına gelmediği konusunda önemlidir. Bu örnekler aynı fiziksel sistem üzerindeki yörünge boyunca toplanmış ilişkili noktalardır.

Telafi sonrası yeni yürütme yöntemin aynı makine ve aynı genel yörünge koşulundaki pratik etkisini göstermektedir; farklı parça geometrisi, farklı devir, farklı besleme oranı veya farklı makineye genelleme mevcut veriyle doğrulanmamıştır.

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

• Analitik rijit merdanenin hareket yörüngesi, ona bağlanan küçük ayrık kabuk üzerinden ABAQUS ODB verisinden çıkarılabilmiştir.
• Simülasyon yörüngesi makine koordinat sistemine dönüştürülerek gerçek tek merdaneli sıvama makinesinde yürütülmüştür.
• Harici deplasman sensörüyle gerçek makine hareketi ölçülmüş ve simülasyon yörüngesiyle karşılaştırılmıştır.
• Üç telafisiz yürütmenin ortalama projeksiyon hataları 0,458, 0,504 ve 0,546 mm; genel ortalama 0,503 mm olarak bulunmuştur.
• Projeksiyon yöntemi, yörünge hatasını yerel normal doğrultusunda tanımlamak için kullanılmıştır.
• GPR modeli üç deneyden alınan normal-yön hata örnekleriyle oluşturulmuştur.
• GPR tahmini kullanılarak simülasyon yörüngesine konuma bağlı normal-yön telafisi uygulanmıştır.
• Telafi sonrası yeni makine yürütmesinde ortalama hata 0,229 mm olarak ölçülmüştür.
• Kaynak hata azalmasını %54,47 olarak raporlamıştır.
• Welch t-testi ve Mann–Whitney U testi mevcut veri kümesinde telafi sonrası hata azalmasını istatistiksel olarak anlamlı bulmuştur.

Çalışmanın desteklemediği veya test etmediği sonuçlar

• Çalışma, bitmiş sıvama parçasının boyutsal hatasının %54,47 azaldığını göstermemektedir.
• Yörünge uydurma deneyinin ana bölümünde fiziksel iş parçası kullanılmamıştır.
• Cidar kalınlığı dağılımı, yüzey pürüzlülüğü veya şekillendirilmiş parça geometrisi telafi öncesi/sonrası karşılaştırılmamıştır.
• GPR farklı makine modellerinde bağımsız olarak doğrulanmamıştır.
• Farklı devir, ilerleme oranı veya farklı yörünge geometrilerinde aynı performans gösterilmemiştir.
• Makine dinamiği, titreşim ve kontrol cevabından kaynaklanan dinamik yörünge hataları mevcut modelin ana kapsamına dahil edilmemiştir.
• 102 GPR eğitim örneği için bağımsız bir test kümesi kullanılmamıştır.
• Çalışma çevrimiçi gerçek zamanlı adaptif GPR telafisi göstermemektedir; telafi modeli ölçülen verilerden oluşturulup yeni yörünge olarak tekrar yürütülmektedir.
• 0,229 mm değerinin bütün metal sıvama makinelerinde ulaşılabilir doğruluk sınırı olduğu gösterilmemiştir.

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

Çift konik modelin temel geometrik parametreleri

Sonlu eleman çözümünün teknik özeti

Merdane açısı ve köşe yarıçapı için kaynak metni ile Tablo 3 arasında 30°/5 mm ve 25°/6 mm farklılığı bulunduğundan bu iki parametre yukarıdaki teknik tabloda tek bir nihai proses değeri olarak birleştirilmemiştir.

Yörünge çıkarımının veri zinciri

Gerçek makine yörünge ölçümü

Projeksiyon analizi için deney sayıları

118 + 127 + 125 = 370 projeksiyon örneği, sonraki GPR modelinin başlangıç hata veri kümesini oluşturmuştur.

GPR veri hazırlama ve model parametreleri

Telafi öncesi ve sonrası sonuçlar

Telafi miktarının 0,3873 mm olması ile telafi sonrasındaki 0,229 mm artık hatanın birbirinden farklı kavramlar olduğu unutulmamalıdır. Birincisi yörüngeye uygulanan ortalama düzeltme büyüklüğünü, ikincisi telafi sonrasında kalan ortalama geometrik sapmayı temsil etmektedir.

İstatistiksel doğrulama

İki farklı test kullanılması, kaynağın yalnız normal dağılım varsayımına dayalı tek bir istatistiksel sonuca bağlı kalmamasını sağlamaktadır. Bununla birlikte telafisiz grup üç yürütmeden, telafi edilmiş grup tek yürütmeden oluşmaktadır; yörünge boyunca çok sayıda nokta bulunması bağımsız fiziksel deney tekrarlarının sayısını artırmamaktadır.

Şekiller sonuçları nasıl destekliyor?

Şekil 17, üç telafisiz yürütmede hedef ve gerçek yörüngelerin genel biçimde birbirine yakın olduğunu ancak büyütülmüş konik bölgede fark bulunduğunu göstermektedir. Şekil 18 bu farkın önemli nedenlerinden birini görünür hâle getirir: gerçek makine yörüngesi basamaklı, simülasyon eğrisi süreklidir.

Şekil 19–21'de siyah simülasyon eğrisi, yeşil gerçek makine yörüngesi ve kırmızı projeksiyon noktaları birlikte gösterilmektedir. Bu görseller projeksiyon yönteminin neden koordinat bileşeni farkından daha uygun görüldüğünü açıklamaktadır.

Şekil 22 GPR'nin yörünge boyunca değişen normal hatayı ve tahmin belirsizliğini modellediğini; Şekil 23 telafi edilmiş yörüngenin üç ölçüm setinin genel yönüne doğru kaydırıldığını; Şekil 24 ise telafi büyüklüğünün yay uzunluğuna göre sabit olmadığını göstermektedir.

Şekil 25 telafi edilmiş yörünge ile yeni ölçülen makine yörüngesinin yakınlığını; Şekil 26 telafi sonrası işaretli ve mutlak hataların yay uzunluğuna göre dağılımını; Şekil 27 ise üç telafisiz veri grubu ile telafi edilmiş veri grubunun kutu grafiklerini ve iki istatistiksel testin p değerlerini göstermektedir.

Kaynak içindeki teknik belirsizliklerin özeti

Araştırmacıların belirttiği gelecek çalışma sınırları

Çalışmanın sonuç bölümü, mevcut modelin esas olarak statik yörünge sapmalarına odaklandığını açıkça belirtmektedir. Makine dinamikleri, titreşim ve kontrol cevabının oluşturduğu dinamik hatalar mevcut çalışma kapsamına dahil edilmemiştir.

Ayrıca hata modeli sınırlı deneysel veri üzerine kurulmuştur. Araştırmacılar gelecekte dinamik hata modelleme yöntemlerinin eklenmesini ve farklı çalışma koşullarında daha geniş deney veri kümeleriyle yöntemin genellenebilirliğinin sınanmasını önermektedir.

Kaynak ve Yöntem Notu

Tam özgün çalışma adı: Roller Trajectory Extraction and Error Compensation Based on Projection Method and Gaussian Process Regression

Yazarlar: Wei Liang, Shi-Yuan Kong, Xin Zhao, Pi-Yao Liu.

Yazar sırası: Kaynaktaki sıra aynen korunmuştur.

Sorumlu yazar: Shi-Yuan Kong.

Eş katkı/eş birinci yazar: Kaynakta böyle bir beyan bulunmamaktadır.

Kurumlar: School of Mechanical and Electrical Engineering, Changchun University of Science and Technology, Changchun, China; Zhongshan Institute, Changchun University of Science and Technology, Zhongshan, China.

Dergi: Machines.

Yayınevi: MDPI.

Bibliyografik kayıt: Machines 2026, 14(6), 676.

DOI: 10.3390/machines14060676.

Resmî yayın bağlantısı:https://doi.org/10.3390/machines14060676

Yayın tarihi: 10 Haziran 2026.

Kaynak türü ve hakemlik durumu: Hakemli özgün araştırma makalesidir. Çalışma sonlu eleman simülasyonu, fiziksel makine yörünge ölçümü ve veri tabanlı hata telafisini birlikte kullanmaktadır.

Lisans: Creative Commons Attribution (CC BY).

Finansman: Çalışma Zhongshan Research Institute of Changchun University of Science and Technology tarafından yürütülen introduced innovative research team project kapsamında, CXTD2023006 proje numarasıyla desteklenmiştir.

Veri erişilebilirliği: Yazarlar çalışmanın sonuçlarını destekleyen ham verilerin talep üzerine sağlanacağını belirtmiştir. Ayrıca dönüştürülmüş yörünge koordinatları, işlenmiş koordinatlar, üç projeksiyon veri seti, telafi edilmiş yörünge ve telafi sonrası ölçülen koordinatlar ek materyalde Tablo S1–S7 olarak listelenmiştir.

Çıkar çatışması: Yazarlar çıkar çatışması olmadığını beyan etmiştir.

Yazar katkıları: Wei Liang kavramsallaştırma, yöntem, kaynak sağlama, inceleme/düzenleme, gözetim ve finansman ediniminde; Shi-Yuan Kong yazılım ve doğrulamada; Xin Zhao ilk taslak yazımında; Pi-Yao Liu veri düzenlemede görev almıştır. Kaynakta bütün yazarların yayımlanan sürümü okuyup onayladığı belirtilmiştir.

Bu Verianla makalesindeki sonlu eleman modeli, geometriler, malzeme ve temas değerleri, yörünge çıkarım yöntemi, sensör sistemi, projeksiyon denklemleri, GPR yapısı, eğitim veri sayıları, telafi değerleri, deney sonuçları ve istatistiksel analizler incelenen kaynak çalışmaya dayanmaktadır. Bibliyografik kimlik doğrulaması dışında dış kaynaklardan yeni bilimsel deney veya performans sonucu eklenmemiştir.

Temel deneysel sınır: Çalışmanın merdane yörünge uydurma deneyinde iş parçası kullanılmadığı kaynakta açıkça belirtilmektedir. Telafi sonrası deney de makinenin yörünge hareketini doğrulamaktadır; kaynak bitmiş çift konik parçanın telafi öncesi/sonrası boyutsal hassasiyeti, cidar kalınlığı veya yüzey kalitesi karşılaştırmasını sunmamaktadır. Bu nedenle %54,47 hata azalması merdane yörüngesi eşleşmesine aittir ve aynı oranda parça kalite artışı olarak yorumlanamaz.

GPR doğrulama sınırı: Üç telafisiz deneyden gelen 370 projeksiyon örneğinin tamamı model geliştirmede kullanılmıştır. 1,0 mm yay uzunluğu kutulaması sonrası 102 örnekle GPR eğitilmiş ve bağımsız rastgele eğitim/test ayrımı yapılmamıştır. Telafi edilmiş yörüngenin yeniden makinede yürütülmesi yöntemin aynı sistemdeki etkisini doğrulamakla birlikte farklı makineler ve proses koşulları için genellenebilirlik kanıtı değildir.

Kaynak içi proses parametresi notu: Merdane kurulum açısı ve köşe yarıçapı açıklama metninde 30° ve 5 mm; Tablo 3'te 25° ve 6 mm olarak verilmektedir. Kaynak hangi değer çiftinin nihai simülasyonda kullanıldığını tek anlamlı biçimde açıklamamaktadır.

Kaynak içi koordinat notu: Ayrık kabuk–merdane merkezi dönüşümünde verilen Denklem (2)–(7) ile kaynakta daha sonra yazılan \(P_1=(x_A+13{,}433,\ y_A-22{,}267,\ z_A+0{,}5)\) sayısal ofsetinin bütün bileşenleri doğrudan uyuşmamaktadır. Ayrıca Şekil 10'da makine koordinat eksenleri metre olarak etiketlenirken sonraki deneysel yörünge verileri mm ölçeğinde sunulmaktadır. Bu bilgiler kaynakta bulunduğu biçimde korunmuştur.

İstatistiksel referans notu: Bölüm 4.2'nin ilk açıklaması 0,229 mm'yi telafi edilmiş gerçek ve telafi edilmiş simülasyon yörüngesi arasındaki sapma olarak tanımlarken sonraki istatistik açıklaması telafi edilmiş ölçüm dahil bütün veri kümelerinin özgün sonlu eleman yörüngesine projekte edildiğini belirtmektedir. Kaynak bu iki referans tanımını ayrıca uzlaştırmamaktadır.

Model kapsamı: Araştırmacılar mevcut çalışmanın statik yörünge sapmalarına odaklandığını, makine dinamiği, titreşim ve kontrol cevabından kaynaklanan dinamik hataların incelenmediğini ve daha geniş deney veri kümeleriyle ek doğrulama gerektiğini belirtmektedir.