Akademik araştırmalar, anlaşılır dil

Verianla | Akademik Araştırmalardan Türkçe Ekonomi ve Bilim İçerikleri

27 Eylül 2026, Pazar
VERİANLABağımsız bilim yayıncılığı
Menüyü aç veya kapat
...
Home / Uygulamalı Bilimler / Mühendislik / Sınır Elemanları Yöntemiyle Üç Boyutlu Kuvvet Alanı Yeniden Yapılandıran Çok Modlu Görsel-Dokunsal Algılama
Mühendislik

Sınır Elemanları Yöntemiyle Üç Boyutlu Kuvvet Alanı Yeniden Yapılandıran Çok Modlu Görsel-Dokunsal Algılama

Bu çalışma, robotik sistemlerin bir nesneye temas ederken kuvveti, yüzey dokusunu, yer değiştirmeyi ve sıcaklığı aynı sensör üzerinden algılayabilmesi için BEM-Tactile adlı fizik tabanlı bir yöntem sunmaktadır.

03/08/2026  Veri Anla 23 görüntüleme
Sınır Elemanları Yöntemiyle Üç Boyutlu Kuvvet Alanı Yeniden Yapılandıran Çok Modlu Görsel-Dokunsal Algılama

Bu çalışma, robotik sistemlerin bir nesneye temas ederken kuvveti, yüzey dokusunu, yer değiştirmeyi ve sıcaklığı aynı sensör üzerinden algılayabilmesi için BEM-Tactile adlı fizik tabanlı bir yöntem sunmaktadır. Araştırmacılar, görsel-dokunsal bir sensörün elastomer yüzeyindeki işaretçilerin üç boyutlu hareketini kullanarak temas kuvvetlerini Sınır Elemanları Yöntemiyle (Boundary Element Method, BEM) yeniden yapılandırmıştır. Yöntem, sonlu eleman yaklaşımındaki gibi sensörün bütün hacmini ağ elemanlarına ayırmak yerine yalnızca yüzey sınırlarını kullanmaktadır. Böylece hacim içindeki ölçülemeyen yer değiştirmelerin tahmin edilmesine gerek kalmadan, kamera tarafından izlenen yüzey hareketleri doğrudan üç boyutlu kuvvet alanına dönüştürülmektedir.

Araştırma, ticari bir GelSight mini sensörün yalnızca algılama kapağı değiştirilerek yürütülmüştür. Standart 0,8 mm işaretçiler 0,5 mm’ye küçültülmüş, sıcaklıkla renk değiştiren termokromik mürekkep kullanılmış ve sensör malzemesi Ecoflex 00-45 ile değiştirilmiştir. Kırmızı görüntü kanalı sıcaklık değişimini, yeşil kanal işaretçi hareketini, bütün RGB kanalları ise yüzey dokusunu hesaplamak için kullanılmıştır. BEM ağı 440 yüzey elemanı ve 442 düğümden oluşmuş; bunların 63’ü kamera tarafından takip edilen işaretçilere karşılık gelmiştir. Yer değiştirme hataları x, y ve z yönlerinde sırasıyla 0,10 mm, 0,07 mm ve 0,20 mm olarak bildirilmiştir. Kesme kuvveti hataları Fx için 0,240 N, Fy için 0,202 N; 0–7 N aralığındaki normal kuvvet hatası ise 0,66 N olmuştur.

Türkiye açısından değerlendirme: Yaklaşım; Türkiye’de robotik kavrayıcılar, insan-robot iş birliği, savunma ve havacılık sistemleri, medikal robotik, protezler, hassas montaj, gıda işleme ve tekstil yüzeyi analizi alanlarında araştırma potansiyeli taşımaktadır. Özellikle tek bir optik sensörle kuvvet, kayma, sıcaklık ve yüzey dokusunun birlikte ölçülmesi; değişken biçimli veya sıcak nesnelerin robotlarla güvenli şekilde tutulmasına yönelik çalışmalar için kullanılabilir. Türkiye’de uygulamaya geçilmeden önce sensörün yerel olarak üretilebilen elastomer ve mürekkeplerle yeniden kalibre edilmesi, farklı aydınlatma koşullarında sınanması, endüstriyel kirlenme ve aşınmaya karşı dayanıklılığının ölçülmesi ve belirli robotik görevlerde gerçek zamanlı gecikmesinin doğrulanması gerekir. Araştırma, Türkiye’deki fabrikalarda doğrudan kullanıma hazır bir ürün, güvenlik sertifikası, maliyet analizi veya seri üretim performansı sunmamaktadır.

Çalışmanın en önemli sonucu, üç boyutlu kuvvet dağılımının yalnızca toplam kuvvet olarak değil, sensör yüzeyindeki 63 farklı noktaya karşılık gelen vektör alanı biçiminde hesaplanabilmesidir. Bununla birlikte sistemin doğruluğu sıcaklık, aydınlatma dağılımı, temas derinliği ve işaretçilerin yüzey dokusunu örtmesi gibi etkenlerden etkilenmektedir. Özellikle 50 °C’de termokromik renk değişimi, hesaplanan derinliğin ve kuvvetin zamanla yapay biçimde yükselmesine yol açmıştır. Bu nedenle sistemin sıcaklıkla değişen ortamlarda nicel kuvvet ölçümü için yerinde kalibrasyona veya renkten bağımsız derinlik hesaplama yöntemlerine ihtiyaç duyduğu gösterilmiştir.

Araştırmanın temel sorusu nedir?

İnsan dokunma duyusu; basınç, sürtünme, yüzey biçimi, hareket ve sıcaklık gibi çok sayıda bilgiyi aynı anda birleştirebilmektedir. Robotik dokunma sensörlerinde ise bu büyüklüklerin aynı donanımla eş zamanlı ölçülmesi güçtür. Görsel-dokunsal sensörler, deforme olabilen şeffaf veya yarı şeffaf bir elastomer yüzeyin kamera ile görüntülenmesine dayanır. Kamera, yüzeydeki işaretçilerin hareketini ve temas yüzeyinin görüntüsünü izleyerek biçim ile yer değiştirme bilgisi üretir.

Bu sensörler yüzey geometrisi ve doku görüntülemesinde başarılı olsa da yüzeyde ölçülen yer değiştirmeyi fiziksel kuvvet dağılımına dönüştürmek ayrı bir ters problem oluşturmaktadır. Makine öğrenmesine dayalı yaklaşımlar etiketlenmiş kuvvet verilerine ihtiyaç duyabilir ve eğitim aralığı dışındaki temaslarda sınırlı genelleme gösterebilir. Basit polinom uyarlamaları ise elastomer içindeki fiziksel dengeyi doğrudan temsil etmez.

Araştırmanın ana sorusu, yüzeydeki işaretçilerin ölçülen üç boyutlu hareketlerinden fiziksel temele sahip, hesaplama yükü düşük ve gerçek zamanlı kullanılabilecek bir üç boyutlu kuvvet alanının elde edilip edilemeyeceğidir. İkinci soru ise kuvvet hesabı için gereken optik işaretçilerin sıcaklık ve doku ölçümleriyle birbirine karışmasının nasıl azaltılabileceğidir.

BEM-Tactile ile BEMT sensörü arasındaki fark nedir?

BEM-Tactile, yüzey yer değiştirmelerinden temas kuvvetlerini hesaplayan genel yöntemdir. BEMT sensörü ise bu yöntemin doğrulanması için algılama kapağı değiştirilmiş ticari GelSight mini sensöre araştırmacıların verdiği addır.

Çalışma tamamen yeni bir kamera veya optik kutu geliştirmemiştir. Ticari sensörün görüntüleme kutusu, kamera ve aydınlatma mimarisi korunmuş; yalnızca elastomer algılama kapağı değiştirilmiştir. Araştırmacılar bu tercihin amacını yeni bir sensör ürünü geliştirmekten çok, BEM tabanlı kuvvet hesaplamasının mevcut bir donanım platformunda uygulanabilirliğini göstermek olarak açıklamaktadır.

Sınır Elemanları Yöntemi neden kullanılmıştır?

Sonlu Elemanlar Yöntemi (Finite Element Method, FEM), elastomerin bütün hacmini küçük üç boyutlu elemanlara ayırarak gerilme ve yer değiştirmeleri hesaplayabilir. Ancak görsel-dokunsal sensör yalnızca dış yüzeydeki işaretçilerin hareketini ölçmektedir. Elastomerin iç kısmındaki yer değiştirmeler kamera tarafından gözlenemediği için FEM uygulamasında hacim içi durumların ayrıca tahmin edilmesi veya ara değerlerle doldurulması gerekir.

BEM ise fiziksel problemi yalnızca cismin sınırı üzerinde ayrıklaştırır. Çalışmada sensörün üst, alt ve yan yüzeyleri yüzey elemanlarına bölünmüş; hacim içinde düğüm oluşturulmamıştır. Bu yaklaşım, ölçümün zaten yüzeyde yapılması nedeniyle görsel-dokunsal sensörün veri yapısıyla doğrudan uyumludur.

Her düğümün üç boyutlu yer değiştirme ve yüzey çekme vektörleri şu şekilde tanımlanmaktadır:

\[ \mathbf{u}_i = \begin{bmatrix} u_x^{(i)} \\ u_y^{(i)} \\ u_z^{(i)} \end{bmatrix}, \qquad \mathbf{t}_i = \begin{bmatrix} t_x^{(i)} \\ t_y^{(i)} \\ t_z^{(i)} \end{bmatrix} \]

Burada \(\mathbf{u}_i\), i numaralı düğümün x, y ve z yönlerindeki yer değiştirmesini; \(\mathbf{t}_i\) ise aynı noktadaki birim yüzey alanına etkiyen çekme veya temas gerilmesi vektörünü göstermektedir.

BEM’in temel matris denklemi çalışmada şu biçimde verilmiştir:

\[ \mathbf{H}\mathbf{u}=\mathbf{G}\mathbf{t} \]

\(\mathbf{H}\) yer değiştirme katsayılarını, \(\mathbf{G}\) yüzey çekme katsayılarını içeren matristir. Sensör geometrisi, elastik özellikler ve yüzey ağı değişmediği sürece bu matrisler önceden hesaplanabilmektedir.

Sınır koşulları nasıl tanımlanmıştır?

BEM ağı üç düğüm grubuna ayrılmıştır:

  • Grup A: Algılama bölgesindeki 63 optik işaretçinin bulunduğu üst yüzey düğümleridir. Bu noktalardaki yer değiştirmeler kamera tarafından ölçülür, kuvvetler bilinmemektedir.
  • Grup B: Sensör kapağının rijit pencereye bağlı alt yüzeyindeki 195 düğümdür. Bu yüzey sabit kabul edildiği için yer değiştirme sıfırdır.
  • Grup C: Algılama bölgesi dışındaki diğer sınır düğümleridir. Bu bölgeler dış kuvvetten arındırılmış kabul edildiği için yüzey çekmesi sıfırdır.

Bu sınır koşulları matematiksel olarak:

\[ \mathbf{u}_B=0, \qquad \mathbf{t}_C=0 \]

şeklinde tanımlanmıştır. Bilinen \(\mathbf{u}_A\) yer değiştirmeleri kullanılarak bilinmeyen \(\mathbf{t}_A\) temas çekmeleri çözülmektedir. Her düğüm çevresindeki alan \(A\) ile çarpıldığında düğüm kuvvetleri elde edilmektedir:

\[ \mathbf{f}_A=A\mathbf{t}_A \]

Araştırmacılar sabit matris işlemlerini önceden bir \(\mathbf{W}\) araç matrisinde birleştirmiştir. Böylece çalışma sırasında kuvvet hesabı yaklaşık olarak aşağıdaki doğrudan dönüşüme indirgenmektedir:

\[ \mathbf{f}_A=\mathbf{W}\mathbf{u}_A \]

Bu dönüşüm sayesinde her kamera karesinde ölçülen işaretçi yer değiştirmelerinden yeni bir üç boyutlu kuvvet alanı hesaplanabilmektedir. Bununla birlikte ana metinde sistemin saniyedeki kare sayısı, toplam hesaplama gecikmesi veya kullanılan işlemci donanımı için ayrıntılı bir performans tablosu verilmemiştir. Bu nedenle “gerçek zamanlı” nitelemesi gösterimlerdeki çevrim içi çalışmaya dayanmakta, nicel gecikme değeriyle desteklenmemektedir.

BEM ile FEM karşılaştırması ne göstermiştir?

Araştırmacılar, ek bilgi dosyasındaki hesaplama deneylerinde 442 düğümlü BEM çözümünün 25.137 düğümlü yüksek çözünürlüklü FEM çözümüne çok yakın yer değiştirme ve yüzey çekme profilleri verdiğini bildirmektedir. Sayısal olarak bu karşılaştırma yaklaşık 57 kat daha az düğüm anlamına gelmektedir.

Buna karşılık sonuç bölümünde aynı avantaj “yaklaşık sekiz kat daha az düğüm” şeklinde ifade edilmiştir. Ana metindeki 442 ve 25.137 düğüm değerleriyle sonuç bölümündeki sekiz kat ifadesi birbiriyle uyuşmamaktadır. Farklı FEM ağlarının karşılaştırılmış olması mümkün olsa da yüklenen çalışmada bu ayrım açıklanmamıştır.

İlgili FEM karşılaştırmaları Note S1’e yönlendirilmiştir. Ek bilgi dosyası yüklenen çalışma içinde bulunmadığı için ağ yakınsaması, kullanılan elastik malzeme parametreleri, hata ölçütleri ve hesaplama süreleri ayrıntılı olarak incelenememiştir. Dolayısıyla BEM’in hesaplama üstünlüğü çalışmanın desteklediği güçlü bir sonuç olmakla birlikte “eşdeğer doğruluk” iddiasının bütün teknik koşulları ana metinden doğrulanamamaktadır.

Ticari sensör üzerinde hangi değişiklikler yapılmıştır?

DeğişiklikUygulamaAmaç
İşaretçi küçültmeÇap 0,8 mm’den 0,5 mm’ye düşürülmüştür.İşaretçilerin yüzey dokusunu örtmesini azaltmak.
Termokromik mürekkepSıcaklıkla maviden açık griye dönüşen üst katman kullanılmıştır.Sıcaklık dağılımını optik olarak ölçmek.
Kırmızı arka plan mürekkebiTermokromik katmanın altına sabit kırmızı mürekkep yerleştirilmiştir.Renk değişimi sırasında işaretçi kontrastını korumak.
Şeffaf silikon yapıştırıcıÜç katmanlı işaretçinin tabanında kullanılmıştır.İşaretçiyi elastomer yüzeye sabitlemek.
Elastomer değişimiSolaris yerine Ecoflex 00-45 kullanılmıştır.Yoğun doğrulama deneylerinde mekanik dayanıklılığı artırmak.

İşaretçi çapının 0,8 mm’den 0,5 mm’ye düşürülmesi, kaplanan yüzey alanını yaklaşık %61 azaltmıştır. Küçük işaretçiler, kamera görüntüsünde yeterli piksel kontrastı sağlamayı sürdürürken doku görüntüsünde daha az örtülmeye neden olmuştur.

Optik kanallar birbirinden nasıl ayrılmıştır?

Üç katmanlı işaretçinin en üstünde mavi termokromik mürekkep, ortasında kırmızı ve sıcaklığa duyarsız mürekkep, altında ise şeffaf silikon taban bulunmaktadır. Düşük sıcaklıkta işaretçi ağırlıklı olarak mavi görünür. Sıcaklık 40 °C’nin üzerine çıktığında üst katman açık griye dönüşmeye başlar ve alttaki kırmızı katman daha görünür hâle gelir.

Tek bir kamera görüntüsü üç farklı amaçla işlenmiştir:

  • Kırmızı kanal: İşaretçi rengindeki değişimden sıcaklık haritası çıkarılmıştır.
  • Yeşil kanal veya gri yoğunluk: İşaretçilerin yatay düzlemdeki konumları ve hareketleri takip edilmiştir.
  • RGB kanallarının tamamı: Fotometrik stereo yöntemiyle yüzey eğimleri ve derinlik haritası hesaplanmıştır.

Bu spektral ayrım, farklı ölçümlerin ayrı kameralarla yapılması yerine aynı optik kare içinde doğal olarak hizalanmasını sağlamaktadır. Ancak sıcaklıkla oluşan renk değişiminin fotometrik stereo derinliğini etkilemesi tamamen ortadan kaldırılamamıştır.

Üç boyutlu yer değiştirme nasıl ölçülmüştür?

Sensör, düzlem içindeki hareketi işaretçi takibiyle; yüzeye dik hareketi ise fotometrik stereo derinlik haritasıyla ölçmektedir. Derinlik haritasının ham değeri milimetre değil piksel tabanlı olduğu için fiziksel kalibrasyon yapılmıştır.

Çapı 2 mm olan yuvarlak uçlu bir prob, algılama yüzeyindeki dokuz farklı konuma bastırılmıştır. Bu konumlar merkez, köşeler ve kenar orta noktaları çevresine dağıtılmıştır. Motorlu doğrusal sahneden alınan gerçek girinti derinliği ile kamera algoritmasının verdiği ham derinlik karşılaştırılmıştır.

Gerçek derinlik ile ham derinlik arasında yaklaşık doğrusal bir ilişki bulunmuştur:

\[ z(x,y)=k(X,Y)\,\tilde{z}(\tilde{x},\tilde{y}) \]

Burada \(\tilde{z}\) piksel tabanlı ham derinlik, \(z\) milimetre cinsinden gerçek derinlik ve \(k(X,Y)\) sensör yüzeyindeki konuma bağlı kalibrasyon katsayısıdır. Katsayı ikinci dereceden bir yüzeyle modellenmiştir:

\[ k(X,Y)=\sum_{i=0}^{2}\sum_{j=0}^{2}a_{ij}X^iY^j \]

X ve Y, 0 ile 1 arasında normalize edilmiş görüntü koordinatlarıdır. Katsayıların yüzey boyunca önemli ölçüde değişmesi, RGB aydınlatmanın sensör üzerinde homojen olmadığını göstermiştir. Merkez ve üst-sol bölgelerdeki hassasiyet ile alt-sağ bölgedeki hassasiyet aynı değildir.

Piksel koordinatlarını milimetreye dönüştürmek için yaklaşık 0,065 ölçek katsayısı kullanılmıştır. Ancak \(a_{ij}\) katsayılarının bulunduğu Table S3 yüklenen çalışmada yer almadığından kalibrasyon yüzeyi yeniden oluşturulamamaktadır.

İşaretçinin başlangıç konumu \((x_0,y_0)\), temas sırasındaki konumu \((x_c,y_c)\) olarak tanımlanmıştır. Yer değiştirme bileşenleri:

\[ u_x=x_c-x_0,\qquad u_y=y_c-y_0,\qquad u_z=z(x_c,y_c) \]

şeklinde hesaplanmıştır. Bildirilen toplam kök ortalama kare hataları x yönünde 0,10 mm, y yönünde 0,07 mm ve z yönünde 0,20 mm’dir. Yüzeye dik yöndeki hata, iki düzlem içi bileşenden daha yüksektir.

Üç boyutlu kuvvet alanı nasıl doğrulanmıştır?

Algılama yüzeyi 440 dört düğümlü yüzey elemanına ve toplam 442 düğüme ayrılmıştır. Kamera tarafından izlenen 63 işaretçinin kuvvet vektörleri, algılama bölgesindeki üç boyutlu kuvvet alanını oluşturmuştur.

Doğrulama deneyinde düz disk biçimli bir prob, motorlu sahneyle sensöre bastırılmış ve yatay yönde hareket ettirilmiştir. Prob üzerine bağlı üç eksenli kuvvet/tork sensörü referans ölçüm sağlamıştır. Kesme kuvvetleri 0,3 mm, 0,4 mm ve 0,5 mm normal girinti altında test edilmiştir.

ÖlçümAralıkGenel RMSE
Fx kesme kuvveti−1,4 ile 2,3 N0,240 N
Fy kesme kuvveti−1,8 ile 2,1 N0,202 N
Fz normal kuvvet0 ile 7,0 N0,66 N

Kesme kuvveti doğruluğu normal girinti derinliği arttıkça kötüleşmiştir. 0,3 mm girintide Fx ve Fy hataları sırasıyla 0,169 N ve 0,158 N iken 0,5 mm girintide 0,310 N ve 0,248 N olmuştur.

Normal kuvvet yaklaşık 7 N’yi aştığında hesaplanan değerler referanstan belirgin biçimde ayrılmıştır. Araştırmacılar bunu disk probun algılama bölgesinin büyük kısmını kaplamasına, yüzeydeki ışık gradyanını azaltmasına ve fotometrik stereo derinlik hesabını bozmasına bağlamıştır.

Bu sonuçlar, kuvvet rekonstrüksiyonunun yalnızca elastik BEM modeline değil, önceki aşamada üretilen optik derinlik ve işaretçi hareketlerinin doğruluğuna da bağlı olduğunu göstermektedir.

Sıcaklık nasıl ölçülmüştür?

Termokromik işaretçiler yaklaşık 40–50 °C arasında belirgin renk değişimi göstermiştir. Araştırmacılar 20–55 °C aralığında bütün işaretçilerin kırmızı kanal yoğunluklarını ölçmüştür. En güçlü optik ayrım, kırmızı arka plan mürekkebinin kırmızı görüntü kanalıyla birlikte kullanıldığı durumda elde edilmiştir.

Bu düzenlemede kırmızı kanal yoğunluk aralığı 20 °C’de 26–75 iken 55 °C’de 44–167 aralığına yükselmiştir. Her işaretçi için kırmızı yoğunluk ve sıcaklık arasında ayrı doğrusal kalibrasyon uygulanmıştır. 40–50 °C ölçüm aralığında sıcaklık RMSE değerleri işaretçinin konumuna göre 0,26–1,40 °C arasında değişmiştir.

Hatalar sensörün orta bölgesinde daha düşük, sol, üst ve alt kenarlarda daha yüksek bulunmuştur. Bu dağılım, sensör kutusundaki aydınlatmanın uzamsal olarak eşit olmadığını göstermektedir. Dolayısıyla tek bir genel sıcaklık kalibrasyonu yerine her işaretçi için konuma özgü kalibrasyon gerekmiştir.

Sıcaklık kuvvet hesabını nasıl bozmuştur?

Termokromik mürekkebin renk değiştirmesi sıcaklık ölçümünü mümkün kılarken fotometrik stereo algoritmasına yeni bir hata kaynağı eklemiştir. Aynı 10 mm çaplı düz probla aynı normal girinti 40 °C ve 50 °C’de uygulanmıştır.

40 °C’de ölçülen yüzeye dik yer değiştirme ve kuvvet zaman boyunca kararlı kalmıştır:

  • Ortalama Uz: 0,38 ± 0,01 mm
  • Ortalama Fz: 1,70 ± 0,03 N

50 °C’de fiziksel girinti sabit tutulmasına rağmen hesaplanan derinlik yaklaşık 150 saniyede 0,38 mm’den 0,44 mm’ye yükselmiş, hesaplanan kuvvet ise 1,7 N’den 3,5 N’ye çıkmıştır. Derinlik ve kuvvetteki yapay artış logaritmik eğrilerle yüksek uyum göstermiştir:

\[ U_z=0{,}03\ln(t-48{,}6)+0{,}29,\qquad R^2=0{,}98 \]

\[ F_z=0{,}99\ln(t-39{,}5)-1{,}66,\qquad R^2=0{,}99 \]

Bu artış gerçek mekanik kuvvet artışı değildir. İşaretçilerin maviden kırmızıya yönelen renk değişimi, fotometrik stereo algoritmasının yüzeyi olduğundan daha derin algılamasına yol açmıştır. BEM de bu hatalı yer değiştirmeyi fiziksel giriş olarak kabul ettiği için kuvvet değerini yükseltmiştir.

Sonuç, sıcaklık ile kuvvetin aynı sensörde ölçülebilmesine rağmen iki ölçümün bütünüyle bağımsız olmadığını göstermektedir. Değişken sıcaklıklarda nicel kuvvet ölçümü için yerinde yeniden kalibrasyon, ışık dağılımı düzeltmesi veya renk değişimine duyarsız fotometrik stereo yöntemleri gereklidir.

Yüzey dokusu hangi çözünürlükte ölçülmüştür?

Sensör, 300, 500, 800 ve 1.000 µm dalga boyuna sahip paralel testere dişi dokulu metal yüzeylerle denenmiştir. Sensör yüzeye 500 gramlık ağırlıkla bastırılmıştır.

500 µm ve üzerindeki dokular kamera görüntülerinde ve derinlik haritalarında belirgin biçimde yeniden oluşturulmuştur. 300 µm doku ise algılanabilmiş olmakla birlikte daha belirsiz görünmüş ve işaretçi örtüşmesinden daha fazla etkilenmiştir. Tablo 1’de doku algılama sınırı 300 µm olarak verilse de ana metin, güvenilir ve açık yeniden yapılandırmanın daha çok 500 µm üzerindeki örneklerde sağlandığını göstermektedir.

ISO 25178-2:2012 yaklaşımına göre dokuz yüzey parametresi hesaplanmıştır:

SembolYüzey parametresiTemel anlamı
SqKök ortalama kare yükseklikYüzey yüksekliğinin genel değişkenliği
SskÇarpıklıkYükseklik dağılımının simetrisi
SkuBasıklıkYükseklik dağılımının sivriliği veya kuyruk yapısı
SpEn yüksek tepeOrta düzlemin üzerindeki en yüksek nokta
SvEn derin çukurOrta düzlemin altındaki en derin nokta
SzToplam maksimum yükseklikEn yüksek tepe ile en derin çukur arasındaki fark
SaAritmetik ortalama yükseklikYüzey sapmalarının ortalama büyüklüğü
SdqKök ortalama kare gradyanYüzey eğimlerinin büyüklüğü
SdrGeliştirilmiş yüzey alanı oranıDoku nedeniyle düz yüzeye eklenen alan oranı

Doku dalga boyu arttıkça Sq, Sp, Sv, Sz, Sa, Sdq ve özellikle Sdr yükselmiştir. Test yüzeylerinin diş geometrileri benzer olduğu için Ssk ve Sku değerlerinde düzenli bir dalga boyu ilişkisi görülmemiştir.

İşaretçilerin küçültülmesi doku ölçümünü ne kadar iyileştirmiştir?

Optik işaretçiler yer değiştirme ölçümü için gereklidir; ancak yüzey görüntüsünün bir bölümünü kapatarak doku ölçümünü bozar. Araştırmacılar 500 µm adımlı testere dişi yüzeyi dört sensör kapağıyla karşılaştırmıştır:

  • İşaretçisiz özgün GelSight kapağı
  • İşaretçisiz yeni kapak
  • 0,8 mm işaretçili özgün kapak
  • 0,5 mm işaretçili yeni kapak

Özgün 0,8 mm işaretçiler kullanıldığında gradyana duyarlı Sdq parametresi işaretçisiz duruma göre yaklaşık %55, Sdr ise %160’tan fazla yüksek hesaplanmıştır. Bunun nedeni işaretçilerin yüzey ayrıntısını örtmesi ve yerel aydınlatma-yansıma geometrisini değiştirmesidir.

0,5 mm işaretçiler Sdq ve Sdr sapmalarını belirgin biçimde azaltmıştır. Bununla birlikte Sa, Sq, Sp, Sv ve Sz gibi genliğe dayalı büyüklüklerde çoğu durumda %20’nin üzerinde sapma devam etmiştir.

Dolayısıyla işaretçilerin küçültülmesi doku ve kuvvet ölçümü arasındaki dengeyi iyileştirmiş, fakat yüksek doğrulukta yüzey metrolojisi sorununu tamamen çözmemiştir. Malzeme kalite kontrolü gibi hassas doku ölçümü gereken uygulamalarda işaretçisiz bölgeler, ek görüntüleme yöntemleri veya sonradan düzeltme algoritmaları gerekebilir.

Sıcak para deneyi ne göstermiştir?

İlk robotik uygulamada ısıtılmış bir madeni para sensör üzerine bastırılmış ve döndürülmüştür. Kamera, paranın üzerindeki “1” biçimini yüzey dokusu olarak görüntülemiştir. Aynı anda:

  • Sıcaklık artışı alanı ΔT,
  • Düzlem içi yer değiştirme alanı Uxy,
  • Yüzeye dik yer değiştirme Uz,
  • Kesme kuvveti alanı Fxy,
  • Normal kuvvet alanı Fz

hesaplanmıştır. Para döndürüldüğünde yer değiştirme ve kesme kuvveti okları dönme yönündeki temas hareketini göstermiştir. Her işaretçinin kuvvet vektörü kullanılarak toplam düzlem içi tork da hesaplanmıştır.

Bu uygulama, aynı temas olayındaki biçim, sıcaklık, kuvvet ve dönme bilgisinin tek bir kamera koordinat sisteminde birlikte izlenebildiğini göstermektedir.

Sıcak su dolu kap nasıl kavranmıştır?

İkinci uygulamada BEMT sensörü bir robot kavrayıcının parmaklarından biri olarak kullanılmıştır. Robot önce yumuşak ve boş bir kabı tutmuş, ardından kaba sıcak su dökülmüştür.

Su eklenirken sensör, kabın aşağı doğru hareketini yer değiştirme alanında gözlemlemiştir. Araştırmacılar bu hareketi başlangıç kaymasıyla ilişkili bir gösterge olarak değerlendirmiştir. Kuvvet alanı ise aşağı yönde dağılan sürtünme kuvvetlerini göstermiştir.

Yer değiştirme alanında elastomerin sıkışmasına bağlı olarak sağa ve sola açılan hareketler görülürken kuvvet alanında aynı yatay açılma bulunmamıştır. Bu fark, yüzey hareketinin kuvvetle birebir aynı olmadığı ve gerçek temas kuvvetini hesaplamanın yalnızca deformasyonu izlemekten farklı bilgi sağladığı sonucunu desteklemektedir.

Boş kap tutulduğunda düşey kesme kuvveti yaklaşık 0,11 N’ye, sıcak su eklendikten sonra 0,73 N’ye yükselmiştir. Araştırmacılar bu sonuçlardan kabın yaklaşık 22 gram, suyun ise 127 gram olduğunu tahmin etmiştir. Ortalama sıcaklık artışı 24 saniye içinde 8,2 °C olmuştur. Kuvvet değişimi hızlı gerçekleşirken termokromik mürekkep ve silikonun ısıl tepkisi daha yavaş olmuştur.

Kumaş deneyleri ne göstermiştir?

On farklı örme yapısı, lif türü ve renge sahip kumaş, sensör üzerine 500 gramlık sabit yük uygulanarak ölçülmüştür. Her kumaş için dokuz yüzey parametresi hesaplanmıştır.

Dokuzuncu kumaş, Ssk ve Sku dışındaki parametrelerin çoğunda en yüksek değerlere ulaşarak en belirgin dokuya sahip örnek olarak ayrılmıştır. Üçüncü, yedinci ve onuncu kumaşların yüzey parametreleri daha düşük bulunmuş ve daha az belirgin dokular gösterdiği değerlendirilmiştir.

Ssk ve Sku değerleri kumaşlar arasında gürültülü ve geniş hata aralıklı olduğundan tek başlarına güvenilir doku göstergeleri olarak görülmemiştir. Çalışma, kumaş sınıflandırmasının tek bir yüzey büyüklüğü yerine birden fazla parametrenin birlikte değerlendirilmesini gerektirdiğini göstermektedir.

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

  • Yüzeyde ölçülen üç boyutlu işaretçi hareketleri BEM kullanılarak noktasal üç boyutlu kuvvet vektörlerine dönüştürülebilmektedir.
  • Sınır tabanlı çözüm, sensör hacminin tamamının elemanlara ayrılmasını gerektirmemektedir.
  • 63 işaretçi, algılama bölgesindeki üç boyutlu kuvvet dağılımını temsil edebilmektedir.
  • 0,5 mm işaretçiler, 0,8 mm işaretçilere göre doku örtüşmesini ve özellikle gradyana dayalı yüzey parametrelerindeki hataları azaltmıştır.
  • Kırmızı ve yeşil kanalların ayrılması sıcaklık ile işaretçi konumunun aynı görüntüden ölçülmesini sağlamıştır.
  • Sensör sıcaklık, yer değiştirme, kuvvet ve doku bilgisini robotik temas ve kavrama görevlerinde eş zamanlı olarak sunabilmiştir.
  • Normal sıkıştırma arttıkça kesme kuvveti doğruluğu azalmıştır.
  • Sıcaklığa bağlı renk değişimi düzeltilmezse derinlik ve kuvvet değerlerinde önemli sürüklenmeye neden olabilmektedir.

Çalışmanın kanıtlamadığı sonuçlar

  • BEMT yönteminin bütün elastomerler, sensör geometrileri veya ticari görsel-dokunsal sensörlerde aynı doğrulukla çalışacağı gösterilmemiştir.
  • Yöntemin gerçek zamanlı gecikmesi, saniyedeki kare sayısı ve farklı işlemcilerdeki hesaplama maliyeti ana metinde nicel olarak verilmemiştir.
  • BEM ile FEM karşılaştırmasının bütün hata ve yakınsama ayrıntıları ek dosya bulunmadığı için doğrulanamamıştır.
  • 50 °C çevresindeki kuvvet sürüklenmesi çözülmüş değildir; yalnızca tanımlanmış ve kalibrasyon gereksinimi belirtilmiştir.
  • 300 µm dokunun yüksek doğrulukla yeniden oluşturulduğu gösterilmemiştir; bu doku algılanmış ancak belirsiz kalmıştır.
  • İşaretçi küçültmenin bütün yüzey parametrelerini yüksek doğrulukla düzelttiği gösterilmemiştir.
  • Uzun süreli aşınma, çevrim ömrü, kimyasal dayanıklılık, toz, yağ ve endüstriyel kirlenme etkileri incelenmemiştir.
  • Robotik uygulamalarda otonom kapalı çevrim kavrama kontrolü veya hata durumunda otomatik tepki gösterilmemiştir.
  • Yöntemin seri üretim maliyeti, ticari uygunluğu veya güvenlik sertifikasyonu değerlendirilmemiştir.

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

Çalışmanın temel güçlü yönü, kuvvet hesaplamasını yalnızca veri odaklı bir tahmin problemi olarak değil, elastik cisim sınır koşullarına dayanan fiziksel bir ters problem olarak ele almasıdır. Bu yapı, hesaplanan kuvvetlerin hangi mekanik varsayımlardan üretildiğini açıklanabilir hâle getirmektedir.

İkinci güçlü yön, aynı sensörde dört ölçüm türünün birlikte değerlendirilmesidir. Yüzey şekli, işaretçi hareketi, kuvvet alanı ve sıcaklık aynı kamera görüntüsünde uzamsal olarak hizalıdır. Ayrı kameralar veya zaman içinde birleştirilen görüntüler gerekmemektedir.

Üçüncü güçlü yön, yöntem doğrulamasının yalnızca basit basma deneyleriyle sınırlı kalmamasıdır. Kesme, sıcak temas, döndürme, yumuşak kap kavrama ve farklı kumaş yüzeyleri gibi çeşitli temas biçimleri denenmiştir.

Başlıca sınırlılıklar nelerdir?

Önemli deney ayrıntılarının büyük bölümü ek bilgi dosyasına bırakılmıştır. BEM-FEM karşılaştırması, sensör kapağının üretimi, deney protokolleri, kalibrasyon katsayıları, sıcaklık etkisi, cilt kalınlığı, işaretçi sürüklenmesi ve karmaşık çok noktalı kuvvet alanları Note S1–S10 içinde açıklanmıştır. Bu ekler yüklenen çalışma içinde bulunmamaktadır.

Elastomerin mekanik davranışı BEM hesabında doğrusal elastik sınır modeliyle temsil edilmektedir. Büyük deformasyonlar, viskoelastik gecikme, histerezis, sıcaklığa bağlı malzeme değişimi ve yerel sürtünme davranışının modele nasıl dâhil edildiği ana metinde ayrıntılı biçimde açıklanmamıştır.

Doğrulama çoğunlukla toplam Fx, Fy ve Fz kuvvetleri üzerinden yapılmıştır. Her bir işaretçi noktasında hesaplanan yerel kuvvet vektörlerinin bağımsız referans ölçümü yapılmamıştır. Bu nedenle toplam kuvvet doğruluğu, yerel kuvvet dağılımının her noktada aynı doğruluğa sahip olduğunu tek başına kanıtlamaz.

Termokromik mürekkebin kullanılabilir sıcaklık aralığı yaklaşık 40–50 °C ile sınırlıdır. Oda sıcaklığı çevresindeki küçük değişimler doğrudan hassas biçimde ölçülememektedir. Isıl tepki de mekanik kuvvet tepkisinden daha yavaştır.

Çalışma henüz hakem değerlendirmesinden geçmemiştir. Ayrıca ana metindeki BEM-FEM düğüm azaltma oranı gibi bazı sayısal ifadeler kendi içinde uyumlu değildir.

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

Deneysel sistemin özeti

BileşenKullanılan yapıÇalışmadaki görevi
Temel sensörTicari GelSight mini görüntüleme kutusuKamera, RGB aydınlatma ve görüntü işleme altyapısı
Algılama kapağıAraştırmacılar tarafından değiştirilen kapakElastomer deformasyonu, işaretçi hareketi ve sıcaklık ölçümü
ElastomerEcoflex 00-45Temas altında deforme olan algılama yüzeyi
İşaretçiler0,5 mm termokromik mikro işaretçilerYer değiştirme takibi ve sıcaklık ölçümü
Kuvvet hesabıSınır Elemanları YöntemiÖlçülen yer değiştirmeyi üç boyutlu kuvvet alanına dönüştürmek
Görüntü işlemeGelSight SDK ve OpenCVDerinlik haritası ve işaretçi konumu çıkarmak
Referans kuvvetÜç eksenli kuvvet/tork sensörüBEMT kuvvet ölçümlerini doğrulamak
Hareket sistemiMotorlu doğrusal sahneProb konumu ve girinti derinliğini kontrol etmek

BEM ağının yapısı

Ağ özelliğiDeğer
Yüzey elemanı sayısı440
Toplam düğüm sayısı442
İşaretçi düğümleri, Grup A63
Sabit alt yüzey düğümleri, Grup B195
Diğer serbest sınır düğümleri, Grup C184
Hesaplanan yerel kuvvet vektörü sayısı63 adet üç boyutlu vektör

Ölçüm performansı

Ölçüm türüÖlçüm aralığıRMSEÖnemli koşul
Ux yer değiştirmesi0–1,00 mm0,10 mmDüzlem içi işaretçi takibi
Uy yer değiştirmesi0–1,00 mm0,07 mmDüzlem içi işaretçi takibi
Uz yer değiştirmesi0–1,00 mm0,20 mmKonuma özgü derinlik kalibrasyonu
Fx kuvveti−1,4–2,3 N0,240 NGirinti arttıkça hata yükselmiştir.
Fy kuvveti−1,8–2,1 N0,202 NGirinti arttıkça hata yükselmiştir.
Fz kuvveti0–7,0 N0,66 N7 N üzerinde belirgin sapma oluşmuştur.
Sıcaklık40–50 °C0,26–1,40 °Cİşaretçi konumuna bağlı hata
Doku dalga boyu300 µm algılama sınırıTek bir RMSE verilmemiştir.500 µm üzeri dokular daha açık yeniden oluşturulmuştur.

Basınç derinliğine göre kesme kuvveti hatası

Normal girintiFx RMSEFy RMSE
0,3 mm0,169 N0,158 N
0,4 mmAra değer grafikte verilmiştir; metinde ayrı sayısal RMSE bulunmamaktadır.Ara değer grafikte verilmiştir; metinde ayrı sayısal RMSE bulunmamaktadır.
0,5 mm0,310 N0,248 N

Sıcaklık-kuvvet etkileşimi

KoşulUz davranışıFz davranışıYorum
40 °C0,38 ± 0,01 mm düzeyinde kararlı1,70 ± 0,03 N düzeyinde kararlıOptik-sıcaklık etkileşimi sınırlıdır.
50 °CYaklaşık 0,38 mm’den 0,44 mm’ye yükselmiştir.Yaklaşık 1,7 N’den 3,5 N’ye yükselmiştir.Renk değişimi derinlik ve kuvveti yapay biçimde büyütmüştür.

Şekillerin temel mesajları

  • Şekil 1: Dokuz konumlu derinlik kalibrasyonunu, konuma bağlı kalibrasyon yüzeyini ve işaretçi hareketlerinden oluşturulan üç boyutlu yer değiştirme alanını göstermektedir.
  • Şekil 2: 442 düğümlü BEM ağını, kuvvet doğrulama düzeneğini ve hesaplanan kuvvetlerle üç eksenli referans sensör değerlerinin karşılaştırmasını sunmaktadır.
  • Şekil 3: Termokromik işaretçilerin 20–55 °C arasındaki kırmızı kanal değişimini ve sıcaklık kalibrasyon hatalarının sensör yüzeyindeki dağılımını göstermektedir.
  • Şekil 4: Aynı fiziksel girintide 40 °C’de kararlı kalan, 50 °C’de ise renk değişimi nedeniyle yapay olarak yükselen derinlik ve kuvvet sonuçlarını göstermektedir.
  • Şekil 5: 300–1.000 µm dalga boylu testere dişi yüzeylerin görüntü ve derinlik haritalarını, ayrıca dokuz yüzey parametresinin değişimini göstermektedir.
  • Şekil 6: 0,8 mm ve 0,5 mm işaretçilerin doku görüntüsü ile yüzey parametrelerine etkisini karşılaştırmaktadır.
  • Şekil 7: Sıcak bir madeni paranın bastırılması ve döndürülmesi sırasında sıcaklık, yer değiştirme, kuvvet ve tork ölçümlerini göstermektedir.
  • Şekil 8: Robot kavrayıcının sıcak su doldurulan yumuşak kabı tutarken sıcaklık ve kesme kuvvetini eş zamanlı izlemesini göstermektedir.
  • Şekil 9: On farklı kumaşın yüzey görüntülerini ve hesaplanan dokuz mikro yapı parametresini karşılaştırmaktadır.

Yöntemsel değerlendirme

BEM-Tactile yöntemi, kamera tarafından ölçülebilen yüzey yer değiştirmelerini doğrudan mekanik sınır koşullarıyla ilişkilendirmesi bakımından güçlüdür. Önceden hesaplanan dönüşüm matrisinin her görüntü karesinde tekrar çözülmesi gereken işlemleri azaltması, çevrim içi kuvvet alanı hesaplamasını desteklemektedir.

Bununla birlikte ölçüm zinciri ardışık hata kaynakları içermektedir. Kamera görüntüsü işaretçi konumunu ve derinliği üretmekte, bu yer değiştirmeler BEM’e girmekte ve BEM kuvveti hesaplamaktadır. Aydınlatma, renk, işaretçi örtüşmesi veya fotometrik stereo hatası yer değiştirme değerini bozduğunda fiziksel kuvvet çözümü de hatalı sonuç üretmektedir.

Bu nedenle sistemin başarısı yalnızca BEM denklemlerine değil; sensör malzemesinin mekanik kalibrasyonuna, aydınlatma homojenliğine, renk ayrımına, işaretçi tespitine ve temas geometrisine birlikte bağlıdır.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı:Multimodal Visuotactile Sensing with 3D Force Field Reconstruction based on Boundary Element Method

Yazarlar ve sıraları: Longteng Yu; Qi Wang; Haoyu Wu; Zonghao Gu; Dabiao Liu.

Eş birinci yazar veya eş katkı bilgisi: Çalışmada eş birinci yazarlık ya da eş katkı bildirimi yer almamaktadır.

Sorumlu yazarlar: Longteng Yu ve Dabiao Liu.

Kurumsal bağlantılar:

  1. Hangzhou International Innovation Institute, Beihang University, Hangzhou, Çin.
  2. Key Laboratory of Biomechanics and Mechanobiology, Ministry of Education; Beijing Advanced Innovation Center for Biomedical Engineering; School of Biological Science and Medical Engineering, Beihang University, Pekin, Çin.
  3. Zhejiang Lab, Hangzhou, Çin.
  4. Department of Engineering Mechanics, School of Aerospace Engineering, Huazhong University of Science and Technology, Wuhan, Çin.
  5. Hubei Key Laboratory of Engineering Structural Analysis and Safety Assessment, Huazhong University of Science and Technology, Wuhan, Çin.
  6. State Key Laboratory of Intelligent Manufacturing Equipment and Technology, Wuhan, Çin.

DOI: 10.2139/ssrn.7201959.

Dergi: Hakemli bir dergi adı veya kabul edilmiş dergi sürümü çalışmada yer almamaktadır.

Yayın platformu: SSRN.

Özgün yayınevi: Yer almamaktadır. SSRN bilimsel dergi veya özgün yayınevi değil, Elsevier bünyesindeki erken aşama araştırma ve preprint paylaşım platformudur.

Yayın yılı: 2026.

Kaynak türü: Sensör değişikliği, mekanik modelleme, kuvvet kalibrasyonu ve robotik uygulama deneyleri içeren araştırma preprinti.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Sonuçlar bağımsız hakem değerlendirmesi ve olası sonraki düzeltmeler öncesindeki erken araştırma bulguları olarak değerlendirilmelidir.

Resmî kaynak bağlantısı:https://papers.ssrn.com/sol3/papers.cfm?abstract_id=7201959

Yazar katkıları: Longteng Yu kavramsallaştırma, biçimsel analiz, fon sağlama, araştırma, görselleştirme, ilk taslak ve düzenleme; Qi Wang yazılım, araştırma, yöntem ve görselleştirme; Haoyu Wu araştırma, yöntem ve yazılım; Zonghao Gu araştırma ve görselleştirme; Dabiao Liu kavramsallaştırma, fon sağlama, danışmanlık, görselleştirme ve düzenleme görevlerini üstlenmiştir.

Fon desteği: Çalışma Çin Ulusal Doğa Bilimleri Vakfının 12302248, 12272146 ve 12572188 numaralı hibeleriyle desteklenmiştir. Dabiao Liu için ayrıca Hubei Provincial Natural Science Foundation 2026AFC0619 ve Central Universities Fundamental Research Funds 2024BRA009 destekleri belirtilmiştir.

Bu Türkçe açıklama, yüklenen çalışmanın 39 sayfası; denklemleri, Tablo 1’i ve Şekil 1–9’daki görsel ve sayısal sonuçlar incelenerek hazırlanmıştır. Çalışmada bulunmayan dış bilimsel sonuçlar eklenmemiş; dış kontrol yalnızca başlık, yazarlar, DOI, resmî platform ve yayın durumunun doğrulanması amacıyla kullanılmıştır.

Çalışmada çok sayıda sonuç ek bilgi dosyasına yönlendirilmiştir. Note S1–S10, Fig. S14, Table S3–S5 ve Movie S1–S2 yüklenen sürümde bulunmamaktadır. Bu nedenle BEM-FEM hesaplama karşılaştırmasının ayrıntıları, algılama kapağının bütün üretim parametreleri, deney protokolleri, kalibrasyon katsayıları, uzun dönem işaretçi sürüklenmesi, farklı elastomer kalınlıkları ve karmaşık çok noktalı kuvvet alanlarının sonuçları tam olarak değerlendirilememiştir.

Ana metinde 442 düğümlü BEM çözümünün 25.137 düğümlü FEM çözümüyle karşılaştırıldığı bildirilirken sonuç bölümünde BEM’in yaklaşık sekiz kat daha az düğüm kullandığı yazılmıştır. Verilen düğüm sayıları yaklaşık 57 kat fark göstermektedir. Bu tutarsızlık yazarlar tarafından açıklığa kavuşturulmadan kesin bir hesaplama tasarrufu oranı kullanılmamalıdır.

Yöntem yerel kuvvet alanını hesaplamakla birlikte doğrulama deneyleri esas olarak algılama bölgesine etkiyen toplam Fx, Fy ve Fz kuvvetlerine dayanmaktadır. Her işaretçi noktasındaki yerel kuvvet vektörü bağımsız bir referans kuvvet sensörüyle doğrulanmamıştır.

Termokromik renk geçişi, özellikle 50 °C’de fotometrik stereo derinliğini ve buna bağlı kuvvet hesabını sistematik biçimde yükseltmiştir. Bu nedenle bildirilen çok modlu algılama, sıcaklık ile kuvvet ölçümlerinin her koşulda tamamen ayrıştırıldığı anlamına gelmemektedir.


Paylaş:

Yorumlar incelendikten sonra yayımlanır.Gönderdiğiniz yorum onay sürecine alınır ve uygun bulunduğunda görünür hâle gelir.

Bir yorum bırakın

E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir

Your experience on this site will be improved by allowing cookies Cookie Policy