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 / Pnömatik Olarak Tahrik Edilen Antagonistik Pouch Motor Çiftinin Model Tabanlı Kontrolü
Mühendislik

Pnömatik Olarak Tahrik Edilen Antagonistik Pouch Motor Çiftinin Model Tabanlı Kontrolü

Araştırma, şişirildiğinde kısalan iki yumuşak pnömatik pouch motorun birbirine antagonistik biçimde bağlandığı bir sistemde hassas konum takibini incelemektedir.

14/08/2026  Veri Anla 25 görüntüleme
Pnömatik Olarak Tahrik Edilen Antagonistik Pouch Motor Çiftinin Model Tabanlı Kontrolü

Araştırma, şişirildiğinde kısalan iki yumuşak pnömatik pouch motorun birbirine antagonistik biçimde bağlandığı bir sistemde hassas konum takibini incelemektedir. Bir motorun basıncı artırılırken diğerinin basıncı azaltıldığında iki aktüatör karşıt yönde kuvvet üreterek aralarındaki hareketli yükü sağa veya sola konumlandırmaktadır. Pouch motorların mekanik davranışı ve içlerindeki sıkıştırılabilir havanın basınç dinamiği doğrusal olmadığı için araştırmacılar yalnız konum ile basınç arasındaki statik ilişkiyi değil, hava debisinin pouch içi basınca ve basıncın mekanik kuvvete dönüşmesini de kontrol modeline dahil etmiştir.

Çalışmanın ana katkısı, basınç dinamiğini içeren analitik sistem modeli üzerine kurulan yeni bir doğrusal olmayan takip kontrolcüsüdür. Bilinmeyen model hataları, sürtünme ve diğer bozucular tek bir toplu bozucu \(F\) altında ele alınmış; bu büyüklüğü çevrim içinde tahmin etmek amacıyla Immersion and Invariance (I&I) yaklaşımına dayanan doğrusal olmayan bir gözlemci geliştirilmiştir. Basınç dinamiği kontrol girişinin mekanik yük denklemine doğrudan girmemesine yol açtığı için kontrol yasası iki kademeli, iç içe bir kayma değişkeni üzerinden tasarlanmıştır.

Araştırmacılar yeni kontrolcüyü iki ayrı referans yaklaşımıyla karşılaştırmıştır. Birincisi, aynı mekanik sistemi kullanmasına rağmen pouch içindeki basınç dinamiğini ihmal eden indirgenmiş model tabanlı kontrolcüdür. İkincisi ise fiziksel deneylerde kullanılan klasik PID kontrolcüdür. Böylece yalnız “model tabanlı kontrol PID'den iyi mi?” sorusu değil, pnömatik basınç dinamiğini modele eklemek gerçekten ölçülebilir bir avantaj sağlıyor mu? sorusu da doğrudan sınanmıştır.

Fiziksel sistemde iki pouch motor doğrusal kızak üzerindeki 60 gramlık bir bloğu hareket ettirmiştir. Basınçlar dijital regülatörlerle ayarlanmış, her pouch içindeki basınç ayrı sensörlerle ölçülmüş, konum doğrusal enkoderle izlenmiş ve kontrol döngüsü Arduino Uno R4 WiFi üzerinde 100 Hz'de çalıştırılmıştır. Deneylerde sinüzoidal ve kübik iki konum yörüngesi kullanılmış ve her kontrolcü–yörünge kombinasyonu beş kez tekrarlanmıştır.

Kübik yörüngede yeni tam model kontrolcüsünün ortalama mutlak hatası 0,04 mm ve RMS hatası 0,04 mm'dir. Basınç dinamiğini ihmal eden kontrolcüde bu değerler sırasıyla 0,05 ve 0,06 mm, PID'de 0,05 ve 0,07 mm olarak ölçülmüştür. Araştırmacılar tam model kontrolcüsünün indirgenmiş model kontrolcüsüne kıyasla kübik yörüngede RMSE'yi %33,3 ve MAE'yi %20 azalttığını bildirmektedir.

Sinüzoidal yörüngede yeni kontrolcünün MAE değeri indirgenmiş kontrolcüyle aynı, 0,04 mm'dir; RMSE ise sırasıyla 0,05 ve 0,06 mm'dir. Bu nedenle basınç dinamiğini modele dahil etmenin yararı bütün hata ölçütlerinde ve bütün yörüngelerde aynı büyüklükte değildir. Ayrıca maksimum tekil hata açısından indirgenmiş kontrolcü iki deney profilinde de yeni kontrolcüden az miktarda daha düşük değer vermiştir.

Çalışmanın sonuçları, sıkıştırılabilir havanın dinamiğini kontrol modeline dahil etmenin antagonistik pouch motorların takip doğruluğunu artırabileceğini ve basınç sinyallerini daha düzgün hâle getirebileceğini göstermektedir. Bununla birlikte yeni kontrolcü referans yörüngesinin üçüncü türevine ihtiyaç duymakta, daha karmaşık bir fiziksel modele dayanmakta ve deneylerde kullanılan sistem dışındaki yükler, keskin bozucular, belirgin stiction, değişken sızıntı veya çok serbestlik dereceli yumuşak robotlarda henüz doğrulanmış değildir.

Türkiye açısından değerlendirme: Çalışmanın yöntemi Türkiye'de yumuşak robotik, pnömatik yapay kaslar, hafif robotik mekanizmalar, medikal robotik ve güvenli insan–robot etkileşimi üzerine çalışan araştırma ekipleri açısından değerli bir kontrol örneği sunmaktadır. Özellikle düşük maliyetli ve kolay üretilebilen yumuşak aktüatörlerde yalnız mekanik konum modelinin değil pnömatik besleme dinamiğinin de kontrol tasarımına katılması, hassasiyet gerektiren uygulamalarda dikkate alınabilecek bir yaklaşımdır. Buna karşılık çalışma cerrahi işlem, hasta veya klinik görev üzerinde test edilmemiştir; “robotik cerrahiye uygunluk” makalenin motivasyonlarından biridir, klinik başarı sonucu değildir.

Pouch motor nedir?

Pouch motor, birbirine bağlanmış ince ve esnek iki tabakanın hava ile şişirilmesi sonucu eksenel olarak kısalan yumuşak pnömatik aktüatördür. Araştırmadaki motorlar, uzamayan termoplastik tabakaların özel bir lazer kaynak sistemiyle birleştirilmesiyle üretilmiş pouch'lardan oluşmaktadır.

Şişme sırasında kesenin kesiti genişlerken iki ucu birbirine yaklaşır. Bu nedenle basınç arttıkça aktüatörün toplam uzunluğu azalır ve mekanik çekme hareketi elde edilir.

Antagonistik çift neden kullanılıyor?

Tek bir pouch motor yalnız şişirilerek kısalma eğilimindedir. İki motor karşılıklı bağlandığında ise biri kısalırken diğeri uzayabilecek şekilde iki yönlü doğrusal hareket oluşturulabilir.

Çalışmadaki düzende hareketli blok iki motor arasında bulunmaktadır. Pouch motor 2 şişirilip pouch motor 1 boşaltıldığında blok pozitif \(x\) yönüne; tersi yapıldığında negatif \(x\) yönüne hareket etmektedir.

Bu düzen insan kaslarındaki antagonistik çiftlere benzer bir mekanik işlev görmektedir: iki aktüatör aynı hareket ekseni üzerinde karşıt yönlerde çalışır.

Deney düzeneği nasıl görünüyor?

Çalışmanın Şekil 1'i dört önemli bileşeni birlikte göstermektedir: şişirilmiş tek bir pouch'un geometrik kesiti, seri bağlı pouch motorun profili, iki motor arasındaki hareketli yükün şematik antagonistik bağlantısı ve gerçek laboratuvar prototipi.

Gerçek düzende pouch motorlar yatay doğrusal kızak üzerindeki 60 gramlık yükün iki tarafına bağlanmıştır. Motorların hava hatları ayrı regülatörlere gider; pouch basınçları sensörlerle ve yükün konumu doğrusal enkoderle ölçülür.

Geometrik model nasıl oluşturuldu?

Araştırmacılar şişmiş pouch'un yaklaşık silindirik bir profile sahip olduğunu varsaymaktadır. Pouch yüzey malzemesi uzamaz kabul edildiği için sabit pouch uzunluğu:

\[ L_0=2R\theta \]

ve şişmiş durumdaki eksenel uzunluk:

\[ L=2R\sin\theta \]

olarak yazılmaktadır.

Bunların birleştirilmesiyle:

\[ L=L_0\frac{\sin\theta}{\theta} \]

elde edilir. Burada \(L_0\) pouch'un sabit yüzey uzunluğunu, \(L\) şişmiş pouch'un eksenel uzunluğunu, \(R\) geometrik yarıçapı ve \(\theta\) merkez açısını temsil etmektedir.

Kısalma neden Taylor açılımıyla sadeleştirildi?

Bir pouch motor \(n_L\) adet seri pouch içerdiğinde toplam kısalma:

\[ \delta L=n_LL_0 \left( 1-\frac{\sin\theta}{\theta} \right) \]

olmaktadır.

Çalışma bunu ikinci dereceden Taylor yaklaşımıyla:

\[ \delta L \approx n_LL_0\frac{\theta^2}{6} \]

biçiminde sadeleştirmiştir.

Araştırmacılar pouch şişmesinin pratikte 45°'yi aşmadığını ve bu açıda söz konusu yaklaşımın truncation hatasının yaklaşık %3 olduğunu belirtmektedir.

Düzeltme katsayısı \(\nu\) neden gerekli?

Gerçek pouch motor ideal geometriden sapmaktadır. Termoplastik malzemenin elastikliği, kaynak çizgisinin kalınlığı ve iki tabakanın birleştirilme biçimi şişmeyi sınırlayabilir. Araştırmacılar bu etkileri doğrudan ayrı ayrı modellemek yerine deneysel olarak belirlenen bir düzeltme katsayısı \(\nu\) kullanmıştır.

İki pouch için açılar:

\[ \theta_1= \sqrt{ \frac{6}{n_LL_0} \left( x_M-x_0-\frac{x}{\nu} \right) } \]

ve:

\[ \theta_2= \sqrt{ \frac{6}{n_LL_0} \left( x_0+\frac{x}{\nu} \right) } \]

olarak ifade edilmektedir.

Burada \(x\) hareketli bloğun konumunu, \(x_M\) toplam hareket geometrisiyle ilişkili mesafeyi, \(x_0\) yarı hareket mesafesini ve \(\nu\) deneysel düzeltme faktörünü temsil etmektedir.

\(\nu\) deneysel olarak nasıl belirlendi?

Açık çevrim deneyinde pouch motorlara tamamlayıcı bir basınç adımı uygulanmıştır:

  • \(P_1=30\) mbar,
  • \(P_2=350\) mbar.

Model ile deney arasındaki RMSE, \(\nu=1\) ile 5 arasında 0,1 adımlarla taranmış ve en düşük hata:

\[ \nu=3.5 \]

için elde edilmiştir.

Çalışmanın Şekil 2'sinde \(\nu=3,5\) kullanılan simülasyon eğrisinin gerçek açık çevrim cevabını \(\nu=2\) durumuna göre belirgin biçimde daha iyi izlediği görülmektedir.

Araştırmacıların duyarlılık analizine göre \(\nu\)'daki ±%10 değişim RMSE'yi yaklaşık %20 artırmış, ancak mutlak RMSE yine yaklaşık \(5\times10^{-3}\) mm sınırı içinde kalmıştır.

Pouch hacmi nasıl modelleniyor?

Geometrik hacim denklemi trigonometrik terim için Taylor açılımıyla sadeleştirildiğinde:

\[ V(\theta)= \frac{n_LL_0^2}{3} \left( D+\frac{2}{3}d \right)\theta \]

yaklaşımı elde edilmektedir.

\(D\) ve \(d\), pouch motor kesit geometrisinin boyutlarıdır. Araştırmacılar bu hacim yaklaşımındaki truncation hatasının \(\theta=45^\circ\) için yaklaşık %13 olduğunu belirtmektedir.

Bu hata ayrıca ayrı bir düzeltme denklemiyle giderilmemekte; diğer modelleme hatalarıyla birlikte toplu bilinmeyen bozucu \(F\)'nin içinde ele alınmaktadır.

Pouch basıncı nasıl mekanik kuvvete dönüşüyor?

Çalışmanın kullandığı temel ilişki:

\[ F_{\mathrm{pouch}} = -P_{\mathrm{pouch}} \frac{\partial V}{\partial x} \]

şeklindedir.

Bu denklem, pouch hacminin eksenel konuma göre değişimi ile iç basınç çarpımının mekanik kuvvet oluşturduğunu ifade eder.

İki pouch motorun kuvvetleri zıt işaretli olduğundan hareketli yük üzerindeki net kuvvet, pouch 1 ve pouch 2 basınçlarının ve anlık geometrik konumun birlikte belirlediği antagonistik farktan oluşmaktadır.

Hareketli yükün dinamik denklemi

Modelde 60 gramlık yük için temel hareket denklemi:

\[ M\ddot{x} + c\dot{x} + \frac{\kappa_1P_1} {2\nu \sqrt{x_M-x_0-x/\nu}} - \frac{\kappa_2P_2} {2\nu \sqrt{x_0+x/\nu}} + F =0 \]

biçimindedir.

Burada \(M\) hareketli kütleyi, \(c\) eşdeğer sönümü, \(P_1\) ve \(P_2\) pouch basınçlarını, \(\kappa_1\) ve \(\kappa_2\) geometrik katsayıları, \(F\) ise modelleme hataları ve bozucuların toplu etkisini temsil etmektedir.

Basınç dinamiğini eklemek neden önemli?

Bir basınç regülatörüne belirli bir basınç komutu vermek pouch içindeki basıncın anında aynı değere ulaşacağı anlamına gelmez. Hava sıkıştırılabilir; boruların, bağlantıların ve pouch'un hacmi vardır ve valf üzerinden geçen debi sınırlıdır.

Bu nedenle çalışma genel basınç değişimini:

\[ \dot{P} = \Gamma \frac{ \dot{m} - \frac{\partial V}{\partial x}\dot{x} }{V} \]

ilişkisiyle modellemektedir.

İzentropik ideal gaz varsayımında:

\[ \Gamma=\gamma P_{\mathrm{pouch}} \]

kullanılmaktadır. \(\dot{m}\) pouch'a giren veya çıkan net hava kütle debisini, \(V\) anlık pouch hacmini ve \(\gamma\) politropik indeksi temsil eder.

Kontrolcü neden doğrudan basıncı değil debiyi hesaplıyor?

Basınç dinamikleri hesaba katıldığında mekanik sistemin doğal kontrol girdisi pouch'a giren veya çıkan hava kütle debisidir. Fakat deney düzeneğindeki SMC regülatörleri doğrudan debi komutu değil hedef basınç komutu kabul etmektedir.

Araştırmacılar bu iki alan arasındaki dönüşümü ISO 6358 pnömatik orifis modeliyle kurmuştur. Model subsonik ve choked akış bölgelerini ayrı olarak ele almakta ve ideal iki atomlu gaz için kritik basınç oranını:

\[ b=0.528 \]

olarak kullanmaktadır.

Kontrol yasasının istediği debiden upstream regülatör basıncına geçebilmek amacıyla orifis modeli ters çevrilmiştir.

ISO 6358 modeli gerçek regülatörü tamamen temsil ediyor mu?

Hayır. Kaynak açık biçimde modelin sabit keskin kenarlı orifis varsayımına dayandığını ve proportional valve davranışını tüm ayrıntılarıyla temsil etmediğini belirtmektedir.

Orifis çapı \(D_c\) ve uzunluğu \(L_c\) gerçek regülatör iç geometrisi olarak kesin biçimde bilinmemektedir. Araştırmacılar bunun özellikle deneysel basınç sinyallerindeki pürüzlülüğün kaynaklarından biri olabileceğini değerlendirmektedir.

Sistemdeki belirsizlikler nasıl temsil edildi?

Modelde sürtünme, parametre hataları, Taylor yaklaşımı hataları ve diğer modellenmemiş etkiler:

\[ F \]

adlı tek bir toplu bozucu altında birleştirilmiştir.

Temel varsayım:

\[ |\dot{F}|\leq\epsilon \]

şeklindedir.

\(\epsilon\)'nun kontrol tasarımında bilinmesi gerekmemektedir; ancak bozucunun değişim hızının sonlu olması gerekir.

Bu bozucu varsayımı neleri dışarıda bırakıyor?

Çalışmanın kendi açıklamasına göre keskin süreksizlikler bu varsayımın kapsamına girmez. Örneğin:

  • belirgin deadband,
  • stiction,
  • ani sert çevresel temas,
  • çok keskin ve süreksiz bozucular

mevcut teorik çerçevenin doğrudan kapsadığı durumlar değildir.

Araştırmacılar pouch motorları deney başında önceden şişirerek deadband oluşma olasılığını azaltmıştır.

Immersion and Invariance gözlemcisi ne yapıyor?

Toplu bozucunun gerçek değeri ölçülmediği için kontrolcü doğrudan \(F\)'yi kullanamaz. Araştırmacılar bunu tahmin etmek amacıyla I&I tabanlı doğrusal olmayan bir gözlemci geliştirmiştir.

Tahmin hatası:

\[ \Delta= \hat{F}+\eta(\dot{x})-F \]

ve yardımcı fonksiyon:

\[ \eta(\dot{x})=-M\alpha\dot{x} \]

olarak tanımlanmıştır.

Buradaki \(\alpha\), gözlemcinin yakınsama hızı ile gürültü hassasiyeti arasındaki temel ayar parametresidir.

Gözlemci için hangi kararlılık sonucu elde edildi?

Kaynakta Lyapunov analizi sonucunda tahmin hatasının:

\[ |\Delta| \rightarrow \frac{\epsilon}{\alpha} \]

mertebesinde sonlu bir bölgeye yakınsadığı gösterilmektedir.

Teorik olarak daha büyük \(\alpha\) daha küçük hata sınırı ve daha hızlı yakınsama sağlayabilir. Ancak gerçek sistemde yüksek \(\alpha\), konumdan sayısal türevle elde edilen hız sinyalindeki gürültüyü büyütmektedir.

Neden deneyde gözlemci kazancı düşürüldü?

Simülasyonda:

\[ \alpha=40 \]

kullanılırken fiziksel deneylerde:

\[ \alpha=10 \]

seçilmiştir.

Araştırmacılar deneysel sistemde \(\alpha>10\) olduğunda hız kestirimindeki yüksek frekanslı gürültünün aktüatör gecikmesiyle etkileşerek sürekli salınımlara yol açtığını gözlemlemiştir.

Bu nedenle \(\alpha=10\), incelenen donanım için pratikte kullanılabilir en yüksek gözlemci kazancı olarak seçilmiştir.

İç içe kayma değişkeni neden kullanılıyor?

Basınç dinamikleri dahil edildiğinde hava debisi mekanik yükün ivme denklemine doğrudan girmemektedir. Bu nedenle klasik tek katmanlı bir sliding-variable yaklaşımı yeterli değildir.

Önce konum hatası:

\[ e=x_d-x \]

olarak tanımlanır.

Birinci kayma değişkeni:

\[ \zeta=\dot{e}+K_1e \]

ve ikinci, iç içe değişken:

\[ \tau=\dot{\zeta}+K_2\zeta \]

olarak oluşturulur.

Kontrolcünün hedefi önce \(\tau\)'yı, ardından \(\zeta\)'yı ve sonuçta konum hatası \(e\)'yi sıfıra veya zamanla değişen bozucu altında küçük bir sınırlı bölgeye yaklaştırmaktır.

Antagonistik fazlalık nasıl çözülüyor?

İki pouch motor aynı doğrusal hareketi karşıt kuvvetlerle üretebildiğinden sistem aktüasyon bakımından redundant yapıdadır.

Tam model kontrolcüsünde araştırmacılar hava debilerine:

\[ \dot{m}_1+\dot{m}_2=0 \]

kısıtını uygulamaktadır.

Böylece bir motor hava alırken diğer motor karşılık gelen biçimde hava verir ve antagonistik davranış korunur.

Tam model kontrolcüsünün teorik takip sınırı

Zamanla değişen bozucular altında takip hatasının nihai sınırı kaynakta:

\[ \rho' = \frac{\epsilon}{K_1K_2} \sqrt{ \frac{2}{K_1} \left( \frac{1}{2K_2M^2} + \frac{K_i}{2\alpha} \right) } \]

olarak verilmiştir.

Kararlılık için \(K_2>K_1\) koşuluna ek olarak \(\alpha\) için kaynakta türetilen alt sınırın sağlanması gerekir.

Bozucu sabit olduğunda, yani:

\[ \dot{F}=0 \]

durumunda araştırmacılar referans hareketin üstel olarak stabilize edildiğini göstermektedir.

Basınç dinamiğini ihmal eden kontrolcü nasıl farklı?

İndirgenmiş modelde pouch içi basıncın zamana bağlı gelişimi ayrı bir durum olarak ele alınmaz. Regülatöre verilen basınç doğrudan kontrol girdisi kabul edilir.

Antagonistik basınçlar:

\[ P_2=2P_0+\beta-P_1 \]

ilişkisiyle birbirine bağlanmaktadır.

Böylece iki pouch basıncının toplamı yaklaşık sabit tutulur ve sistem sertliğinin büyük ölçüde korunması amaçlanır.

Basitleştirilmiş kontrolcünün avantajı nedir?

İndirgenmiş kontrolcü daha az model bilgisi gerektirir ve tam kontrolcüden farklı olarak referans yörüngesinin jerk değerine, yani:

\[ \dddot{x}_d \]

ihtiyaç duymaz.

Araştırmacılar bu nedenle ani konum değişikliklerinin istendiği veya basınç dinamiğinin yeterince küçük kabul edilebildiği uygulamalarda indirgenmiş yaklaşımın pratik bir alternatif olabileceğini belirtmektedir.

Basitleştirmenin bedeli nedir?

Basınç dinamiğinin kontrol modelinden çıkarılmasıyla bu dinamiğin etkisi fiilen modellenmemiş bozucuya aktarılmaktadır. Sonuçta:

  • gözlemcinin tahmin etmesi gereken toplu bozucu büyüyebilir,
  • daha agresif kazançlar basınç sisteminin modellenmemiş dinamiğini uyarabilir,
  • basınç sinyallerinde daha fazla salınım oluşabilir.

Deneysel Şekil 6–8 bu yorumla uyumlu biçimde, tam model kontrolcüsünde daha yavaş değişen bozucu kestirimi ve bazı koşullarda daha düzgün basınç profili göstermektedir.

Simülasyon modelinde hangi parametreler kullanıldı?

ParametreDeğerAnlam
\(M\)60 gHareketli yük
\(c\)100 kg/sEşdeğer sönüm
\(L_0\)13,5 mmPouch geometrik uzunluğu
\(d\)14 mmPouch geometrisi
\(D\)9 mmPouch geometrisi
\(n_L\)4Motor başına seri pouch sayısı
\(x_M\)11 mmGeometrik hareket parametresi
\(x_0\)5,5 mmYarı hareket parametresi
\(V_0\)3 × 10−10 m³Boru ve fitting dead volume modeli
\(T_0\)298 KAtmosfer sıcaklığı
\(P_0\)101325 PaAtmosfer basıncı
\(\gamma\)1,4Politropik indeks
\(R_s\)287 J/(kg·K)Kuru havanın özgül gaz sabiti
\(\nu\)3,5Deneysel geometri düzeltme katsayısı

Simülasyonda bozucu nasıl tanımlandı?

Bilinmeyen dış bozucu:

\[ F= 0.2\log\left[0.1(1+t)\right]\ \mathrm{N} \]

olarak verilmiştir.

Değişim hızı:

\[ \dot{F}= \frac{0.2}{1+t} \leq0.2 \]

olduğundan kontrol teorisindeki sınırlı-bozucu-türevi varsayımını sağlamaktadır.

%10 parametre belirsizliği testi ne gösterdi?

Araştırmacılar \(c\), \(M\), \(\kappa_1\) ve \(\kappa_2\) parametrelerini simülasyonda her zaman adımında nominal değerlerinin %90 ve %110'u arasında değiştirirken kontrolcüye nominal değerleri vermiştir.

Bu deney fiziksel bir yük değişimi değil, gözlemci ve kontrolcülerin hızlı değişen model uyumsuzluğuna tepkisini görmek için oluşturulmuş sayısal robustness testidir.

İki kontrolcü de sonunda referansı takip etmiştir; ancak tam model kontrolcüsünün başlangıç sapması daha küçük olmuş ve 0,1 mm referans bandına indirgenmiş kontrolcüden yaklaşık 0,1 saniye daha erken ulaştığı raporlanmıştır.

Parametre belirsizliği toplu bozucu büyüklüğünü yaklaşık 2 N artırmış ve pouch basınçlarının sınırlarına yaklaşmasına neden olmuştur. İndirgenmiş kontrolcüde \(P_2\)'nin alt basınç sınırına ulaşması pouch'un gevşemesine yol açabilecek bir durum olarak tartışılmıştır.

Fiziksel deney sistemi hangi donanımı kullandı?

BileşenKullanılan donanım
Pnömatik regülatörSMC ITV1030
DACMCP4922, 12 bit
Basınç sensörüMS4525DO, 14 bit
Doğrusal enkoderHEDS-9200
Encoder strip250 LPI
Quadrature counterLS7366R
Kontrol kartıArduino Uno R4 WiFi
Örnekleme frekansı100 Hz
PC seri iletişim115200 baud
Hareketli blok60 g

Deneylerde basınç sınırları neden uygulandı?

Pouch motorların hasar görmesini önlemek için:

\[ P_{\max}=400\ \text{mbar} \]

üst sınırı kullanılmıştır.

Motorların yetersiz basınç nedeniyle gevşemesini önlemek için:

\[ P_{\min}=30\ \text{mbar} \]

alt sınırı uygulanmıştır.

Her deney başlangıcında \(P_1=P_2=150\) mbar durumunda orta konum kalibre edilmiştir.

Kontrolcü kazançları

KontrolcüSimülasyonDeney
Tam model kontrolcüsü, Denklem (25)\(K_1=320,\ K_2=360,\ \alpha=40\)\(K_1=320,\ K_2=360,\ \alpha=10\)
İndirgenmiş model, Denklem (33)\(K_1=180,\ \alpha=40\)\(K_1=180,\ \alpha=10\)
PIDDeneysel karşılaştırmanın ana parçası değil\(K_p=0,1,\ K_i=0,05,\ K_d=0,01\)

Kazançlar nasıl seçildi?

Model tabanlı kontrolcülerin kazançları tamamen otomatik optimizasyonla bulunmamıştır. Araştırmacılar \(K_1\) ve \(K_2\)'yi teorik kararlılık koşullarını koruyarak artırmış ve maksimum yer değiştirmeye göre overshoot %10'u aşana kadar ayarlamıştır.

\(\alpha\) ise gözlemci tahmininde sürekli salınım görülene kadar yükseltilmiştir.

PID parametreleri, kübik takip testinde model tabanlı kontrolcülerin rise time'ına yaklaşacak ve overshoot'u azaltacak şekilde seçilmiştir.

Dolayısıyla sonuçlar bu belirli tuning prosedürüne aittir; bütün olası PID veya model tabanlı kontrol ayarlarının global optimum karşılaştırması değildir.

Deneylerde hangi referanslar kullanıldı?

İki referans tipi kullanılmıştır:

  • yumuşatılmış kübik geçişler içeren ±1 mm konum profili,
  • 0,05 Hz sinüzoidal konum profili.

Kübik profilin geçiş bölgelerinde birinci türev uç noktalarda sıfır olacak şekilde seçilmiştir. Bu özellik tam model kontrolcüsünün ihtiyaç duyduğu yüksek mertebeli referans türevlerinin daha düzenli olmasını sağlamaktadır.

Her deney kaç kez tekrarlandı?

Tablo 4'teki maksimum hata, MAE ve RMSE değerleri her kontrolcü ve yörünge için beş deneyin ortalama ve standart sapmalarıdır.

Hatalar konum enkoderinin 0,01 mm çözünürlüğüne yuvarlanmıştır.

Sinüzoidal deneyde sonuçlar

KontrolcüMaksimum hata (mm)MAE (mm)RMSE (mm)
Tam model, Denklem (25)0,19 ± 0,050,04 ± 0,000,05 ± 0,00
İndirgenmiş model, Denklem (33)0,16 ± 0,010,04 ± 0,000,06 ± 0,00
PID0,25 ± 0,010,07 ± 0,000,09 ± 0,00

Sinüzoidal görevde tam ve indirgenmiş model tabanlı kontrolcülerin MAE değerleri aynıdır. Tam model kontrolcüsü daha düşük RMSE üretmiş, fakat maksimum hata açısından indirgenmiş kontrolcü daha düşük değer vermiştir.

Kübik deneyde sonuçlar

KontrolcüMaksimum hata (mm)MAE (mm)RMSE (mm)
Tam model, Denklem (25)0,22 ± 0,060,04 ± 0,000,04 ± 0,00
İndirgenmiş model, Denklem (33)0,21 ± 0,000,05 ± 0,010,06 ± 0,00
PID0,24 ± 0,020,05 ± 0,000,07 ± 0,01

Kübik görev, basınç dinamiğinin modele eklenmesinin etkisinin en açık görüldüğü deneydir. Araştırmacılar tam model kontrolcüsünün indirgenmiş kontrolcüye göre RMSE'de %33,3 ve MAE'de %20 iyileşme sağladığını raporlamaktadır.

“%33,3 daha doğru” ifadesi nasıl okunmalı?

Bu oran bütün hata ölçütlerine uygulanmamalıdır. Kaynakta %33,3 iyileşme:

\[ 0,06\ \mathrm{mm} \rightarrow 0,04\ \mathrm{mm} \]

kübik-yörünge RMSE değişimine aittir.

Aynı deneyde maksimum hata:

\[ 0,21\ \mathrm{mm} \rightarrow 0,22\ \mathrm{mm} \]

şeklindedir; yani maksimum tekil hata yeni kontrolcüde az miktarda daha yüksek çıkmıştır.

Bu nedenle çalışmanın güçlü sonucu, tam model kontrolcüsünün hata dağılımının geneli bakımından daha iyi RMSE/MAE vermesidir; her anda ve her hata ölçütünde daha küçük değer ürettiği gösterilmemiştir.

Basınç profilleri ne gösterdi?

Şekil 6 ve 7'de model tabanlı kontrolcülerin genel olarak PID'den biraz daha yüksek pouch basıncı kullandığı görülmektedir. Araştırmacılar bunu daha küçük teorik takip hata sınırları için gereken yüksek kontrol kazançlarıyla ilişkilendirmektedir.

Kaynağa göre kontrolcüler arasındaki basınç farkı herhangi bir anda 50 mbar'ın altındadır.

Pouch motor 2'nin basınç profilinde daha yüksek tepe ve daha hızlı değişimler bulunmuştur. Yazarlar bunun küçük bir hava kaçağı nedeniyle regülatörün daha sık dolum yapmasından kaynaklanabileceğini belirtmektedir.

Bozucu gözlemcisi deneyde nasıl davrandı?

Şekil 8, tam model ve indirgenmiş model kontrolcüleri için gözlemcinin tahmin ettiği toplu bozucuyu göstermektedir.

Tam model kontrolcüsünde tahmin edilen bozucu daha düşük genlikli ve daha yavaş değişirken indirgenmiş kontrolde daha büyük tepe ve değişimler görülmektedir.

Araştırmacılar bunu, basınç dinamiği indirgenmiş modelden çıkarıldığında söz konusu dinamiğin gözlemci tarafından “bilinmeyen bozucu” olarak karşılanmak zorunda kalmasıyla açıklamaktadır.

Basınç dinamiğini modele katmanın temel mühendislik mesajı

Çalışmanın deneyleri yalnız modelin daha karmaşık olması nedeniyle daha iyi sonuç çıktığını göstermeyi amaçlamamaktadır. Ana mekanizma şudur:

Regülatör komutu → hava akışı → pouch basıncı → pouch hacmi ve geometrisi → mekanik kuvvet → yük konumu.

Bu zincirin orta kısmındaki basınç oluşumu ihmal edildiğinde kontrolcü mekanik sistemin komuta gerçek cevabındaki gecikme ve akış etkilerini toplu hata olarak görmek zorunda kalır. Bu çalışma mevcut prototipte bu ara dinamiğin açıkça modellenmesinin takip RMSE'sini azaltabildiğini göstermektedir.

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

  • İki pnömatik pouch motor fiziksel antagonistik çift hâlinde hassas iki yönlü doğrusal konum takibi için kullanılabilmiştir.
  • Pouch geometrisi, mekanik yük ve sıkıştırılabilir hava basınç dinamiğini birleştiren analitik model deneysel sistem için kontrol tasarımında kullanılabilmiştir.
  • I&I tabanlı doğrusal olmayan gözlemci toplu model belirsizliği ve bozucuyu kontrol içinde telafi edebilmiştir.
  • İç içe kayma değişkenli tam model kontrolcüsü için zamanla değişen bozucu altında nihai sınırlılık ve sabit bozucu altında üstel yakınsama koşulları analitik olarak türetilmiştir.
  • Basınç dinamiğini ihmal eden indirgenmiş model tabanlı kontrolcü de deneylerde PID'den daha düşük genel takip hataları üretmiştir.
  • Tam model kontrolcüsü kübik deneyde indirgenmiş kontrolcüye göre RMSE'yi %33,3 ve MAE'yi %20 azaltmıştır.
  • Tam model kontrolcüsü iki deney profilinde de PID'den daha düşük MAE ve RMSE üretmiştir.
  • Basınç dinamiğinin modele katılması deneylerde daha düzgün basınç davranışı ve daha küçük gözlemci bozucu tahminiyle ilişkilendirilmiştir.
  • Sistem yumuşatılmış yaklaşık 2 mm'lik konum değişimini iki saniyeden kısa sürede takip edebilmiştir.

Çalışmanın göstermediği sonuçlar

  • Yeni kontrolcünün bütün yumuşak pnömatik aktüatörlerde %33,3 daha doğru olduğu gösterilmemiştir.
  • %33,3 değeri bütün hata ölçütlerinin iyileşmesini ifade etmemektedir.
  • Yeni kontrolcü deneylerde her durumda en düşük maksimum anlık hatayı üretmemiştir.
  • Çalışma çok serbestlik dereceli yumuşak robot üzerinde deney yapmamıştır.
  • Robotik cerrahi görevi veya klinik uygulama gerçekleştirilmemiştir.
  • İnsan veya hayvan üzerinde deney yapılmamıştır.
  • Değişen payload kütleleri fiziksel deneyde sistematik olarak test edilmemiştir.
  • Dışarıdan uygulanan keskin darbe/temas bozucularına karşı deneysel robustness gösterilmemiştir.
  • State-dependent leakage veya pressure-dependent belirsizlikler ayrı dinamik modeller olarak gözlemciye dahil edilmemiştir.
  • Deadband ve stiction kontrol teorisinin mevcut bozucu varsayımı altında kapsanmamaktadır.
  • ISO 6358 orifis yaklaşımının gerçek dijital regülatör iç akışını eksiksiz temsil ettiği gösterilmemiştir.
  • Çok yüksek kontrolcü kazançları için ayrık-zaman kararlılık sınırları türetilmemiştir.
  • Uzun dönem pouch motor ömrü veya yüksek basıncın malzeme dayanıklılığı üzerindeki etkisi incelenmemiştir.

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

Modelden gerçek aktüasyona kontrol zinciri

Verianla Live: Pouch motor takip kontrolü nasıl çalışıyor?

Bu akış tam model kontrolcüsünün araştırmada kullanılan temel fiziksel ve matematiksel zincirini gösterir.

AşamaİşlemBilimsel anlamı
1Referans konum \(x_d\)Sinüzoidal veya yumuşatılmış kübik hedef hareket belirlenir.
2Konum hatası \(e=x_d-x\)Gerçek yük konumu ile hedef arasındaki fark hesaplanır.
3I&I bozucu gözlemcisiModelleme hatası, sürtünme ve diğer toplu bozucu etkileri tahmin edilir.
4İç içe kayma değişkenleri\(\zeta=\dot e+K_1e\) ve \(\tau=\dot\zeta+K_2\zeta\) oluşturulur.
5Gerekli hava kütle debisiBasınç dinamiğini içeren tam model kontrol yasası debi talebini hesaplar.
6ISO 6358 ters akış modeliİstenen debi regülatör upstream basınç komutuna dönüştürülür.
7Antagonistik basınçlandırmaBir pouch şişerken diğeri boşalır; karşıt kuvvetler yükü konumlandırır.
8Basınç ve konum geri beslemesiP1, P2 ve yük konumu 100 Hz kontrol döngüsünde tekrar ölçülür.
 

Deneysel kontrol karşılaştırması

Verianla Live: Basınç dinamiğini modele eklemek takip hatasını nasıl değiştirdi?

Veriler çalışmanın Tablo 4 sonuçlarıdır. Her değer beş fiziksel deneyden elde edilen ortalama ± standart sapmayı gösterir. Daha düşük MAE ve RMSE daha düşük genel takip hatası anlamına gelir; maksimum hata ayrı bir ölçüttür.

YörüngeKontrolcüMaksimum hata (mm)MAE (mm)RMSE (mm)Kaynak
SinüzoidalTam model – Denklem (25)0,19 ± 0,050,04 ± 0,000,05 ± 0,00Tablo 4
Sinüzoidalİndirgenmiş model – Denklem (33)0,16 ± 0,010,04 ± 0,000,06 ± 0,00Tablo 4
SinüzoidalPID0,25 ± 0,010,07 ± 0,000,09 ± 0,00Tablo 4
KübikTam model – Denklem (25)0,22 ± 0,060,04 ± 0,000,04 ± 0,00Tablo 4
Kübikİndirgenmiş model – Denklem (33)0,21 ± 0,000,05 ± 0,010,06 ± 0,00Tablo 4
KübikPID0,24 ± 0,020,05 ± 0,000,07 ± 0,01Tablo 4
 

Tablo 4'ün ana bilimsel mesajı

Deneysel sonuçlar üç ayrı noktayı birlikte göstermektedir.

  1. Her iki model tabanlı kontrol yaklaşımı da MAE ve RMSE açısından PID taban çizgisinden daha düşük takip hatası üretmiştir.
  2. Basınç dinamiğini içeren tam model özellikle kübik yörüngede indirgenmiş modele göre daha düşük MAE ve RMSE sağlamıştır.
  3. Maksimum tekil hata sıralaması MAE/RMSE sıralamasıyla aynı değildir; tam model kontrolcüsü maksimum hata bakımından indirgenmiş kontrolcüyü geçmemiştir.

Bu nedenle basınç dinamiği eklenmesinin deneysel avantajı “her hata değerini küçülttü” şeklinde değil, genel takip hata dağılımını ve özellikle RMSE'yi iyileştirdi şeklinde ifade edilmelidir.

Simülasyon ve deney kanıtlarının ayrılması

SonuçKanıt türüYorum sınırı
Nominal modelde iki kontrolcünün benzer davranmasıMATLAB simülasyonuModel ve kontrolcü aynı nominal parametreleri kullanıyor
%10 parametre belirsizliğinde tam modelin daha iyi transient davranışıMATLAB robustness simülasyonuParametreler yapay olarak %90–110 arasında değiştiriliyor
Sinüzoidal RMSE 0,05 mmFiziksel prototip, 5 tekrarMevcut 60 g yük ve donanım koşulları
Kübik RMSE 0,04 mmFiziksel prototip, 5 tekrarMevcut 60 g yük ve donanım koşulları
%33,3 RMSE iyileşmesiFiziksel kübik-yörünge karşılaştırmasıTam model C(25) ile indirgenmiş C(33) arasındaki RMSE farkı
Lyapunov kararlılık sonuçlarıAnalitik matematiksel ispatÇalışmada belirtilen model ve bozucu varsayımlarına bağlı

Araştırmanın güçlü yönleri

  • Yalnız kontrol simülasyonu değil gerçek yumuşak pnömatik aktüatör deneyleri içerir.
  • Basınç dinamiğinin modelde bulunmasının etkisini doğrudan indirgenmiş modelle karşılaştırır.
  • Model tabanlı yöntemlere karşı pratik bir PID taban çizgisi kullanır.
  • Bilinmeyen bozucular için doğrusal olmayan gözlemci içerir.
  • Kontrolcü ve gözlemcinin kararlılık özelliklerini Lyapunov analiziyle ele alır.
  • Deneysel sonuçları tek denemeyle değil beş tekrar üzerinden ortalama ve standart sapmayla raporlar.
  • Model parametresi \(\nu\) için açık çevrim deneysel kalibrasyon ve duyarlılık analizi yapar.
  • Gerçek donanımın basınç sınırlarını ve örnekleme hızını açık biçimde raporlar.
  • Model basitleştirmenin hangi koşullarda yine de tercih edilebileceğini tartışır.

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

  • Deneysel sistem tek serbestlik dereceli doğrusal mekanizmadır.
  • Yük yalnız 60 g'lık mevcut prototip koşulunda kapsamlı olarak raporlanmıştır.
  • Orifis geometrisi gerçek proportional regulator iç yapısını tam olarak temsil etmemektedir.
  • Orifis çapı ve etkili uzunluğunun kesin değeri bilinmemektedir ve akış yönüyle değişebilir.
  • Deneysel basınç sinyallerinde olası leakage etkisi ayrı parametre olarak modellenmemiştir.
  • State-dependent bozucular gözlemci tasarımında özel olarak ayrıştırılmamıştır.
  • Haricî yük değişimleri ve dışarıdan uygulanan bozucular fiziksel deneyde sistematik robustness testi olarak kullanılmamıştır.
  • Deadband, stiction ve keskin temas mevcut bozucu varsayımının kapsamı dışında kalmaktadır.
  • Hız konum ölçümünün ayrık türeviyle elde edildiğinden yüksek gözlemci kazançlarında gürültü problemi oluşmuştur.
  • Sürekli-zaman kararlılık analizi, çok yüksek kazançlı ayrık-dijital uygulamanın bütün sınırlarını tanımlamamaktadır.
  • Tam model kontrolcüsü jerk bilgisi gerektirir.
  • Çalışma çok-DoF yumuşak manipülatörde gerçekleştirilmemiştir.
  • Cerrahi veya biyomedikal uygulama doğrudan test edilmemiştir.

Kaynak ve Yöntem Notu

Tam özgün çalışma adı: Model-Based Control of Antagonistic Pair of Pneumatically Actuated Pouch Motors

Yazarlar: Syed Arshad Hussain ve Enrico Franco.

Sorumlu yazar: Syed Arshad Hussain.

Eş katkı/eş birinci yazarlık: Kaynakta eş katkı veya eş birinci yazarlık beyanı bulunmamaktadır.

Kurum: Mechanical Engineering Department, Imperial College London, Exhibition Road, London SW7 2AZ, United Kingdom.

Dergi: Actuators.

Yayınevi: MDPI, Basel, Switzerland.

Bibliyografik bilgi: Actuators 2026, 15(6), 332.

DOI:10.3390/act15060332

Resmî yayın bağlantısı:MDPI – Actuators 15(6), 332

Yayın tarihi: 11 Haziran 2026. Çalışma 23 Nisan 2026'da alınmış, 5 Haziran 2026'da revize edilmiş ve 9 Haziran 2026'da kabul edilmiştir.

Kaynak türü ve hakemlik durumu: Actuators dergisinde yayımlanmış hakemli araştırma makalesidir. Çalışma analitik doğrusal olmayan kontrol tasarımı, MATLAB simülasyonu ve fiziksel prototip deneylerini birlikte içermektedir.

Lisans: Creative Commons Attribution (CC BY).

Akademik editör: Zhuming Bi.

Finansman: Yazarlar araştırmanın dış finansman almadığını bildirmiştir.

Etik kurul: Uygulanabilir değildir.

Bilgilendirilmiş onam: Uygulanabilir değildir.

Veri erişilebilirliği: Yazarlar verilerin makul talep üzerine kendilerinden temin edilebileceğini belirtmektedir.

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

Teşekkür: Yazarlar deneylerde kullanılan pouch motorları sağladığı için Imperial College London Hamlyn Centre'dan Mark Runciman'a teşekkür etmektedir.

Yazar katkıları: Kavramsallaştırma ve araştırma Syed Arshad Hussain ile Enrico Franco; metodoloji ve analiz her iki yazar; yazılım ve ilk yazım Syed Arshad Hussain; inceleme ve danışmanlık Enrico Franco tarafından yürütülmüştür.

Deney kapsamı: Fiziksel deney düzeneği iki antagonistik pouch motorun doğrusal kızak üzerindeki 60 g yükü hareket ettirmesinden oluşmaktadır. Kontrol döngüsü 100 Hz'de çalıştırılmış ve her deneysel hata sonucu beş tekrar üzerinden raporlanmıştır.

Model varsayımı: Çalışma ideal gaz ve izentropik sıkıştırma/genleşme varsayımı kullanır. Şişmiş pouch'un yaklaşık silindirik profil aldığı varsayılır ve ideal geometriden sapmaların bir bölümü deneysel \(\nu=3,5\) düzeltme katsayısıyla karşılanır.

Geometrik yaklaşım sınırı: Pouch kısalma denklemindeki Taylor yaklaşımının 45°'de yaklaşık %3, hacim denklemindeki sadeleştirmenin ise yaklaşık %13 truncation hatası bulunduğu kaynakta belirtilmektedir. Hacim yaklaşımı hatası toplu bilinmeyen bozucu \(F\)'nin içinde ele alınmıştır.

Bozucu varsayımı: Kontrol teorisi bilinmeyen toplu bozucunun değişim hızının sınırlı olmasını, \(|\dot F|\leq\epsilon\), gerektirmektedir. Kaynak deadband, stiction ve keskin çevresel temas gibi süreksiz davranışların bu varsayım tarafından kapsanmadığını açıkça belirtmektedir.

Regülatör modeli sınırı: ISO 6358 sabit orifis modeli proportional digital regulator davranışının temsilî yaklaşımıdır. Etkili orifis çapı ve uzunluğu kesin olarak bilinmemekte ve akış yönüne göre değişebilmektedir.

Gözlemci sınırı: Daha yüksek \(\alpha\) teorik olarak daha hızlı bozucu tahmini sağlarken deneylerde hızın konum sinyalinden ayrık türevle hesaplanması nedeniyle gürültü büyümektedir. Mevcut donanımda \(\alpha>10\) sürekli salınımla ilişkilendirildiği için deneylerde \(\alpha=10\) seçilmiştir.

Kontrolcü karşılaştırma sınırı: PID ve model tabanlı kontrolcüler belirli heuristik tuning prosedürleriyle ayarlanmıştır. Çalışma tüm olası gain setleri üzerinde global controller optimization yapmamaktadır.

%33,3 sonucu için kapsam notu: Kaynakta bildirilen %33,3 doğruluk iyileşmesi, kübik deneyde tam model kontrolcüsünün RMSE değerinin 0,06 mm'den 0,04 mm'ye düşmesine dayanmaktadır. Maksimum hata, sinüzoidal ve kübik testlerin ikisinde de indirgenmiş model kontrolcüsünde az miktarda daha düşüktür. Bu nedenle oran bütün hata metriklerinin veya bütün görevlerin %33,3 iyileştiği anlamına gelmemektedir.

RMSE yuvarlama notu: Tablo 4'te sinüzoidal tam-model RMSE değeri enkoder çözünürlüğüne yuvarlanmış biçimde 0,05 mm olarak verilmiştir; sonuç bölümü ise iki deney profili için RMS hatalarının 0,05 mm'nin altında olduğunu ifade etmektedir. Yuvarlanmamış değerler kaynakta sunulmadığından Verianla içeriğinde kesin sayısal karşılaştırma için Tablo 4'teki değerler esas alınmıştır.

Basınç kullanımı sınırı: Model tabanlı kontrolcüler deneylerde PID'ye göre genel olarak biraz daha yüksek pouch basınçları kullanmıştır. Yazarlar bunun pouch ömrü ve cerrahi güvenlik gereksinimleri açısından teorik olarak istenmeyebileceğini belirtmekte, ancak mevcut sistemde kontrolcüler arası basınç farkının 50 mbar'dan düşük olduğunu bildirmektedir.

Gerçek dünya kapsamı: Araştırma gerçek fiziksel pouch motor prototipi üzerinde doğrulanmıştır; buna rağmen cerrahi robot, insan, hayvan, klinik işlem veya çok serbestlik dereceli yumuşak manipülatör üzerinde deney gerçekleştirilmemiştir.

Gelecek çalışma: Yazarlar modeli çok serbestlik dereceli yapılara genişletmeyi, ayrık-zaman kontrol tasarımını incelemeyi ve state-dependent bozucular ile leakage etkilerini ele alabilecek genişletilmiş gözlemciler geliştirmeyi önermektedir.

Bu Verianla içeriğindeki mekanik model, kontrol denklemleri, deney koşulları, hata değerleri ve bilimsel sonuçlar yalnız incelenen çalışmaya dayanmaktadır. Dış doğrulama yalnız resmî bibliyografik yayın kimliğini kontrol etmek amacıyla kullanılmış; dışarıdan yeni deneysel veya kontrol performansı sonucu eklenmemiştir.


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