Kenta Endo, Manami Kanamaru ve Keita Tanaka tarafından geliştirilen çalışma, 3U sınıfı nanosatellitlerde hızlı yönelme ile hassas hedefleme arasındaki temel değiş tokuşu azaltmak amacıyla değişken hızlı kontrol moment jiroskopları (Variable-Speed Control Moment Gyroscope, VSCMG) için iki aşamalı bir momentum yönetimi stratejisi öneriyor. Sistem büyük açılı manevranın ilk bölümünde yüksek tekerlek açısal momentumu kullanarak yüksek jiroskopik tork üretirken, tutum hatası 5°'nin altına düştüğünde tekerlek momentumu önceden azaltılıyor. Böylece uydu hedef yönelime yaklaşırken düşük ataletli 3U gövdede aşırı salınıma yol açan jiroskopik tork küçültülüyor. Deneysel olarak tanımlanmış gimbal modeli ve tezgâhta ölçülmüş tekerlek dinamiğiyle desteklenen kapalı çevrim simülasyonda 90°'lik manevra 25,30 saniyede tamamlanmış, ortalama yönelme hızı 3,55°/s ve nihai tutum doğruluğu 0,42° olmuştur. En önemli sınırlılık, bu uzay aracı düzeyi sonuçların fiziksel uydu üzerinde değil simülasyonda elde edilmesidir.

Çalışmanın karşılaştırmasında klasik üç reaksiyon tekerlekli sistem aynı 90° manevrada 116 saniyede yakınsamış, maksimum gövde hızı 0,78°/s ve doğruluğu 0,08° olmuştur. İki aşamalı CMG çözümü bu referansa göre yaklaşık 4,5 kat daha yüksek yönelme hızı sağlamaktadır. Geleneksel, tek momentum seviyeli CMG simülasyonunda ise yaklaşık 5,2° nihai hata oluşurken iki aşamalı sistem 0,42°'ye ulaşmıştır. Yazarların genel karşılaştırmasındaki yaklaşık 12 katlık doğruluk iyileşmesi bu iki farklı ayarlı durum arasındadır; makalenin ayrıca yaptığı sabit geri besleme kazançlı kontrollü karşılaştırmada momentum anahtarının tek başına doğruluğu yaklaşık 0,72°'den 0,36°'ye düşürdüğü, yani yaklaşık iki kat iyileştirdiği gösterilmektedir.

Çalışmanın mühendislik yaklaşımı, daha karmaşık bir sürekli kontrol yasası geliştirmekten ziyade aktüatörün fiziksel momentum durumunu yönetmeye dayanıyor. Büyük hata bölgesinde her CMG için yüksek momentum seviyesi \(h_1=1.41\times10^{-3}\,\mathrm{N\,m\,s}\), hedef yakınında ise bunun yarısı olan \(h_2=7.05\times10^{-4}\,\mathrm{N\,m\,s}\) kullanılıyor. Böylece yüksek manevra kabiliyeti manevranın büyük bölümünde korunurken hassas yönelme evresine girilmeden önce tork otoritesi azaltılıyor.

Sonuçlar ayrıca yöntemin tüm çalışma koşullarında otomatik olarak yüksek hassasiyet garanti etmediğini gösteriyor. ±10% atalet, ±20% kontrol kazancı ve ±50% aktüatör gürültüsü değişimleriyle yapılan 60 koşuluk indirgenmiş-model Monte Carlo taramasında manevraların tamamı yakınsamış olsa da nihai doğruluk 0,51° ile 1,57° arasında değişmiştir. Dolayısıyla özellikle yüksek aktüatör gürültüsünde 1° hedefinin aşılabildiği görülmektedir.

3U CubeSat'lerde neden hız ve hassasiyet aynı anda zorlaşıyor?

Reaksiyon tekerlekleri küçük uydularda yüksek yönelim hassasiyeti sağlayabilir; ancak ürettikleri tork sınırlı olduğu için büyük açılı manevralar yavaş gerçekleşir. Kontrol moment jiroskopları (CMG) ise yüksek hızda dönen bir tekerleğin açısal momentumunu bir gimbal ekseni etrafında yönlendirerek çok daha yüksek tork üretebilir. Bunun bedeli, özellikle düşük ataletli nanosatellitlerde, yüksek depolanmış momentumun hedefe yaklaşırken aşırı tork, overshoot ve eksenler arası salınım oluşturabilmesidir.

Çalışmada kullanılan 3U uydu için atalet matrisi şu şekilde alınmıştır:

\[ I = \operatorname{diag}(0.034,\;0.034,\;0.0067) \quad [\mathrm{kg\,m^2}] \]

X ve Y eksenlerindeki asal atalet momentleri 0,034 kg·m² iken Z eksenindeki değer yalnız 0,0067 kg·m²'dir. Z ekseni ataleti diğer iki eksenden yaklaşık beş kat daha küçüktür. Bu asimetri nedeniyle küçük bir çapraz eksen CMG torku bile yaw kanalında belirgin salınım oluşturabilmektedir.

Geliştirilen CMG donanımı nasıl yapılandırılmış?

Çalışmada dört CMG, üç gövde eksenine dengeli tork aktarabilmek için 54,73° eğimli piramit düzeninde yerleştirilmiştir. Fiziksel CMG modülünün boyutları 90 × 90 × 105 mm, toplam kütlesi 584 gramdır. Her bir prototip tekerleği 13,4 mm çapındadır ve pirinçten üretilmiştir.

Her tekerleğin açısal momentumu:

\[ h_i = I_w\Omega_i \]

ile tanımlanır. Burada \(I_w\) tekerleğin dönme ekseni etrafındaki atalet momenti, \(\Omega_i\) ise i'nci tekerleğin açısal hızıdır. Nominal durumda dört tekerlek aynı momentum büyüklüğünde çalıştırılmıştır.

Burada önemli bir ölçek ayrımı vardır. Fiziksel prototip tekerleğinin atalet momenti yaklaşık \(1.35\times10^{-7}\,\mathrm{kg\,m^2}\)'dir ve yaklaşık 10.000 rpm'de \(1.41\times10^{-4}\,\mathrm{N\,m\,s}\) momentum üretmektedir. Uzay aracı simülasyonunda kullanılan Faz 1 momentumu ise bunun on katı olan \(1.41\times10^{-3}\,\mathrm{N\,m\,s}\)'dir. Yazarlar prototipi uçuş ölçekli momentum üretiminin doğrudan gösterimi olarak değil, gerekli hız geçişlerinin ve kontrol zarfının gerçekleştirilebilirliğini sınamak amacıyla kullanmıştır.

Stratejinin ayırt edici özelliği, kontrol kazançlarını sürekli değiştirmek yerine aktüatörde depolanmış fiziksel açısal momentum miktarını değiştirmesidir. Dolayısıyla yöntem klasik gain scheduling, momentum dumping veya yalnız komut torkunu şekillendiren bir referans yöneticisiyle aynı değildir.

CMG torku neden momentum azaltılınca küçülüyor?

Dört-CMG piramit grubunun ürettiği tork çalışmada:

\[ T_C=C_j\dot{\delta} \]

şeklinde verilmektedir. \(T_C\) uydu gövdesine uygulanan CMG torkunu, \(C_j\) gimbal geometrisine ve tekerlek momentumuna bağlı Jacobian matrisini, \(\dot{\delta}\) ise dört gimbalin açısal hız vektörünü ifade eder.

Jacobian matrisinin tamamı tekerlek momentum büyüklüğü \(h\) ile çarpılmaktadır. Bu nedenle aynı gimbal hareketinde \(h\)'nin azaltılması mümkün tork büyüklüğünü de doğrudan azaltır. Önerilen yöntemin temel fiziksel mantığı budur: yüksek \(h\) hızlı manevra sağlar; hedef yakınında daha düşük \(h\) ise düşük ataletli gövdenin aşırı uyarılmasını sınırlar.

5° anahtarlama eşiği neden seçilmiş?

Yazarların amacı tekerlek hızındaki düşüşün hassas yönelme bölgesine girmeden tamamlanmasını sağlamak, ancak yüksek momentumlu evreyi gereğinden erken sonlandırmamaktır. Bu nedenle çalışma 5° hata eşiğini nominal değer olarak kullanmıştır.

İndirgenmiş modelde eşik 2°, 5°, 8° ve 12° değerleri için ayrıca taranmıştır. Bu taramada tepe gövde hızı dört durumda da 3,55°/s ve nihai doğruluk yaklaşık 0,30° olarak kalmış; belirgin biçimde değişen değer momentum geçişinin başlama zamanı olmuştur. Eşik 2° iken geçiş 28,4 saniyede, 12° iken 23,9 saniyede başlamıştır.

Ancak bu eşik duyarlılığı sonuçları tam gimbal dinamiğinin yer aldığı ana simülasyonla aynı modelden üretilmemiştir. Yazarlar Bölüm 3.5'te bu çalışmanın hızlı parametre taraması için kullanılan indirgenmiş model olduğunu özellikle açıklamaktadır. Örneğin indirgenmiş modelde minimum tekillik göstergesi \(m_{\min}=1.09\) iken tam modelde yaklaşık 0,07'dir. Bu nedenle eşik tablosu mutlak performans tahmini değil, eşik değişiminin göreli etkisinin göstergesi olarak değerlendirilmelidir.

Tekillik sorunu neden önemli?

CMG sistemlerinde bazı gimbal geometrileri belirli bir yönde istenen torkun üretilemediği kinematik tekilliklere yol açabilir. Çalışma bu durumu aşağıdaki manipulability göstergesiyle izlemektedir:

\[ m= \sqrt{\det(A_gA_g^T)} \]

\(A_g=C_j/h\) boyutsuz geometrik Jacobian'dır. \(m=0\) tekil bir konfigürasyonu ifade eder. Çalışma, literatürde yaygın kullanılan yaklaşık \(m_0=0.1\) uyarı eşiğini temel almaktadır.

Ana 90° manevrada sistem iki kez bu eşik bölgesine yaklaşmış ve yaklaşık \(m_{\min}=0.07\) görülmüştür. Yazarlar yalnız pseudo-inverse yönlendirmeye güvenmenin bu bölge için sağlam olmadığını belirterek Generalized Singularity-Robust (GSR) inverse kullanmıştır.

GSR yönlendirme yasası:

\[ \dot{\delta} = C_j^T \left( C_jC_j^T+\lambda E \right)^{-1} T_r \]

biçimindedir. Burada \(T_r\) referans torku, \(\lambda E\) ise tekilliğe yaklaşıldığında matris terslemesini düzenleyen regularization terimidir. Operability yüksek olduğunda düzenleme etkisi çok küçülür; sistem tekilliğe yaklaştığında büyüyerek gimbal hız komutlarının kontrolsüz artmasını sınırlar.

VSCMG aynı zamanda reaksiyon tekerleği gibi kullanılabiliyor mu?

Evet. Çalışma ayrıca dört gimbali belirli açılarda kilitleyerek VSCMG tekerleklerinin yalnız hızlanma ve yavaşlamasından tork üretildiği bir RW modu simüle etmiştir. Bu mod yüksek hız için değil hassas tutum koruması için değerlendirilmiştir.

15°'lik roll manevrasında uydu 0,1° hata bandına 351 saniyede ulaşmış, son 10 saniyedeki en kötü durum kararlı hata 0,013° olmuştur. Dolayısıyla aynı VSCMG kümesinin yüksek torklu CMG modu ile oldukça yavaş ancak daha hassas RW modu arasında farklı görev fazlarına hizmet edebileceği gösterilmektedir.

Tekerlek hızını yarıya indirmek mikro-titreşimi neden azaltıyor?

Prototip tekerlekler ayrıca dengesizlik kaynaklı mikro-titreşim açısından değerlendirilmiştir. Kalıntı kütle eksantrisitesinin oluşturduğu senkron radyal kuvvet:

\[ F(\Omega)=U\Omega^2, \qquad U=m_we \]

ile modellenmiştir. \(m_w\) tekerlek kütlesi, \(e\) kalıntı eksantrisitesi ve \(\Omega\) tekerlek açısal hızıdır. Kuvvet hızın karesiyle arttığı için tekerlek hızı 10.000 rpm'den 5000 rpm'ye indirildiğinde bu senkron titreşim bileşeni yaklaşık dörtte bire düşmektedir.

Kaynak, 2,5–15 µm eksantrisite aralığı için 10.000 rpm'de yaklaşık 16–99 mN, 5000 rpm'de ise yaklaşık 4–25 mN senkron radyal kuvvet tahmin etmektedir. Bunlar yörüngede ölçülmüş payload jitter değerleri değil, imalat toleranslarından hareketle yapılan aktüatör düzeyi tahminlerdir.

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

Uzay aracı dinamiği nasıl modellenmiş?

Uydu dönme dinamiği Euler denklemiyle modellenmiştir:

\[ I\dot{\omega}_s+ \omega_s\times(I\omega_s) = T_C+T_d \]

Burada \(I\) uydu atalet matrisi, \(\omega_s\) gövde açısal hızı, \(T_C\) CMG kontrol torku ve \(T_d\) dış bozucu torktur.

Ana kapalı çevrim modelinde dış bozucu:

\[ T_d=T_a+T_g \]

olarak alınmış ve 560 km irtifadaki aerodinamik sürükleme ile yerçekimi gradyanı etkileri dahil edilmiştir. Kullanılan atmosfer yoğunluğu \(2.108\times10^{-13}\,\mathrm{kg/m^3}\), kesit alanı 0,03 m² ve sürükleme katsayısı 2,2'dir.

Kaynakta tahmin edilen karakteristik bozucu torklar; maksimum yerçekimi gradyanı için yaklaşık \(4.9\times10^{-8}\,\mathrm{N\,m}\), aerodinamik sürükleme için yaklaşık \(8.0\times10^{-9}\,\mathrm{N\,m}\), güneş radyasyon basıncı için yaklaşık \(4.4\times10^{-9}\,\mathrm{N\,m}\) ve kalıntı manyetik dipol için \(2.5\times10^{-8}\) ile \(2.3\times10^{-7}\,\mathrm{N\,m}\) aralığındadır.

Güneş radyasyon basıncı ve kalıntı manyetik tork ana simülasyona dahil edilmemiştir. Yazarlar yaklaşık 25 saniyelik hızlı manevrada bunların oluşturacağı açısal impulsun CMG'nin değiştirdiği gövde açısal momentumundan çok küçük olduğunu hesaplamaktadır. Bununla birlikte çalışma, daha uzun zaman ölçeklerinde kalıntı manyetik dipolün önemli bir bozucu olabileceğini de açıkça belirtmektedir.

Quaternion geri besleme denetleyicisi

Komut torku quaternion hata geri beslemesiyle oluşturulmaktadır:

\[ T_e=-K_{qp}q_v-K_{qd}\omega_s \]

\(q_v\), hedef ve mevcut tutum arasındaki hata quaternionunun vektör bölümünü; \(K_{qp}\) oransal tutum kazancını; \(K_{qd}\) açısal hız sönümleme kazancını ifade eder.

Nominal kararlılık analizi için kullanılan Lyapunov adayı:

\[ V= \frac{1}{2} \omega_s^TI\omega_s+ 2K_{qp}(1-q_4) \]

ve ideal koşullardaki türevi:

\[ \dot V= -K_{qd}\|\omega_s\|^2 \le 0 \]

olarak verilmiştir. Bu sonuç nominal quaternion geri besleme çevriminin asimptotik kararlılığını desteklemektedir. Ancak kaynak açıkça bunun tam iki-aşamalı anahtarlamalı VSCMG sisteminin biçimsel kararlılık kanıtı olmadığını belirtmektedir; gimbal saturasyonu, anahtarlama, aktüatör gecikmesi ve model belirsizliği tam kanıta dahil değildir.

Gimbal modeli deneysel olarak nasıl çıkarılmış?

Gimbal motoru, sürücüsü ve encoder'ı birlikte kara kutu sistem olarak ele alınmıştır. 6 V genlikli, 0,1–100 Hz frekans taraması 400 Hz örnekleme hızıyla uygulanmış ve ikinci dereceden transfer fonksiyonu elde edilmiştir:

\[ G(s)= \frac{108600} {s^2+210.9s+1460} \]

Model hold-out doğrulama segmentinde %93,2 zaman alanı uyumu sağlamıştır. Bununla birlikte çalışma farklı genliklerde tekrar deneyleri veya bağımsız step-input doğrulaması yapılmadığını da sınırlılık olarak kaydetmektedir.

Gimbal hız kontrolü için PID, LQR ve Integral-type Optimal Servo (IOS) karşılaştırılmıştır. 0,5 Hz sinüzoidal takipte IOS yaklaşık 0,02 saniyelik faz gecikmesi göstermiş; PID yaklaşık 0,31 s, LQR ise yaklaşık 0,14 s gecikmiştir. Step bozucuda IOS asimptotik olarak sıfır kalıcı hataya ulaşırken PID için −16,9°/s ve LQR için 260°/s kalıcı sapma raporlanmıştır. Bu nedenle kapalı çevrim simülasyonlarında IOS tercih edilmiştir.

Tekerlek dinamiği sentetik değil, ölçülen geçişe dayanıyor

Tekerlek hız denetleyicisi prototip üzerinde kalibre edilmiştir. 10.000 rpm'den 5000 rpm'ye ölçülen yavaşlama zaman serisi kapalı çevrim uydu simülasyonuna doğrudan dinamik giriş olarak verilmiştir. Ölçülen geçiş birinci dereceden yaklaşık modelle \(\tau\approx0.6\) s zaman sabiti, %94,7 uyum ve yaklaşık 70 rpm RMSE ile temsil edilebilmektedir; %95 yerleşme süresi yaklaşık 1,8 saniyedir.

90° manevrada ana sonuç

İki aşamalı kontrol olmadan geleneksel CMG durumu yaklaşık 27,61 saniyede yakınsamış ancak roll ekseninde 5,17° undershoot ve yaw ekseninde yaklaşık 3,3° salınım oluşmuştur.

İki aşamalı momentum kontrolü kullanıldığında aynı 90° roll hedefi 25,30 saniyede gerçekleştirilmiş, ortalama gövde yönelme hızı 3,55°/s ve nihai tutum hatası 0,42° olmuştur.

Bu karşılaştırmada iki aşamalı CMG, reaksiyon tekerlekli referansa göre yaklaşık 4,5 kat daha yüksek yönelme hızı sağlarken, reaksiyon tekerleğinin 0,08°'lik nihai doğruluğu daha iyidir. Dolayısıyla çalışma “CMG her açıdan reaksiyon tekerleğinden daha iyi” sonucunu desteklememektedir; önerilen yöntemin avantajı hız ve sub-degree doğruluk arasında daha dengeli bir çalışma noktası sağlamasıdır.

12 katlık doğruluk iyileşmesi yalnız momentum anahtarından mı geliyor?

Hayır; makalenin kendisi bu noktayı ayrıca kontrol etmiştir. Ana karşılaştırmada iki aşamalı ve iki aşamasız simülasyonlarda quaternion geri besleme kazançları farklıdır. Bu nedenle yaklaşık 5,2° → 0,42° değişiminin tamamını yalnız momentum anahtarına atfetmek doğru değildir.

Yazarların Şekil 13'te yaptığı sabit kazançlı kontrollü karşılaştırmada aynı P–D denetleyici üzerinde momentum anahtarı etkinleştirildiğinde nihai doğruluk yaklaşık 0,72°'den 0,36°'ye iyileşmiştir. Dolayısıyla kontrol kazançları sabit tutulduğunda momentum anahtarının kendi katkısı yaklaşık iki katlık doğruluk artışıdır.

Aynı karşılaştırmada yalnız türev kazancının iki katına çıkarılması veya jiroskopik feedforward eklenmesi, iki aşamalı momentum yönetiminin sağladığı son doğruluğu üretmemiştir. Yazarlar bu sonucu iyileşmenin esas bileşeninin aktüatör momentumunun fiziksel olarak azaltılması olduğu şeklinde yorumlamaktadır.

GSR tekillik önleme olmadan ne oluyor?

Çalışmanın tekillik testinde GSR kullanılmadığında sistem yaklaşık 6,04 saniyede tekil konfigürasyona ulaşmış ve gerekli tork yönünü üretemediği için hedef tutuma ulaşamamıştır.

GSR inverse etkin olduğunda uydu hedefe 26,44 saniyede ulaşmıştır. Operability göstergesi 6,97 saniyede \(m=0.0519\)'a kadar düşmüş, ancak regularization sayesinde sistem tekilliğe yakın bölgeden geçerek yeniden daha iyi koşullanmış bir gimbal konfigürasyonuna dönmüştür.

Bir CMG arızalanırsa sistem çalışmaya devam ediyor mu?

Çalışma CMG4'ün tamamen arızalandığı ve Jacobian'daki ilgili sütunun kaldırıldığı üç-CMG durumunu da simüle etmiştir. Bu senaryoda roll hata bileşeni yaklaşık 85,06 saniyede, yaw bileşeni ise yaklaşık 89,54 saniyede yakınsamıştır. Kararlı durum tutum doğruluğu 0,95° içinde kalmıştır.

Bununla birlikte arıza sırasında yaw hatası 23. saniyede 39,40°'ye kadar çıkmıştır. Bu büyük geçici sapma, üç aktif CMG'nin hedefe en kısa tutum yolunu sağlayacak optimum tork yönünü her anda üretemediğini göstermektedir. Dolayısıyla sonuç “tek CMG arızası performansı etkilemez” anlamına gelmez; yalnız sistemin ciddi performans kaybıyla da olsa hedefe ulaşabildiğini göstermektedir.

Yaklaşık 85 saniyelik yakınsama dört-CMG iki aşamalı durumdaki 25,3 saniyeden yaklaşık 3,4 kat daha yavaştır. Buna rağmen çalışmanın 116 saniyelik reaksiyon tekerleği referansından daha hızlı kalmaktadır.

Anahtarlama eşiği duyarlılığı

İndirgenmiş modelde kullanılan eşik taraması aşağıdaki sonuçları vermiştir:

Bu tarama, 2°–12° aralığında manevra hızının eşiğe çok hassas olmadığını; daha büyük eşiğin esas olarak momentum azaltımını daha erken başlattığını göstermektedir. Ancak tablo tam gimbal modeline ait değildir ve 0,30° mutlak doğruluk değeri ana tam-model sonucu olan 0,42° ile karıştırılmamalıdır.

Monte Carlo analizi yöntemin ne kadar sağlam olduğunu gösteriyor?

İndirgenmiş kapalı çevrim modelinde 60 Monte Carlo koşusu yapılmıştır. Uydu ataleti ±10%, geri besleme kazançları ±20% ve aktüatör gürültüsü ±50% değiştirilmiştir.

Tüm koşular hedef hata bandına yakınsamış ve tepe gövde hızı 3,55 ± 0,00°/s olarak kalmıştır. Ortalama yerleşme süresi 32,7 ± 3,8 saniyedir. Nihai doğruluk ise 0,87 ± 0,24° olmuş ve 0,51°–1,57° arasında değişmiştir.

Bu sonuç yöntemin manevra hızının incelenen belirsizliklere karşı sağlam olduğunu, ancak son yönelim doğruluğunun aktüatör gürültüsünden belirgin biçimde etkilendiğini göstermektedir. En yüksek gürültü koşullarında hata 1° hedefini aşabildiği için çalışma bütün belirsizlik kombinasyonlarında sub-degree doğruluk garantisi vermemektedir.

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

• İki seviyeli tekerlek momentumu, düşük ataletli bir 3U uydu modelinde CMG'nin yüksek manevra kabiliyeti ile hedef yakınındaki daha düşük tork ihtiyacını aynı aktüatör kümesinde birleştirebilir.
• Tam kapalı çevrim simülasyonunda 90° roll manevrası 25,30 saniyede, 3,55°/s karakteristik hız ve 0,42° doğrulukla tamamlanmıştır.
• Reaksiyon tekerlekli referansa göre manevra hızı yaklaşık 4,5 kat artmıştır; reaksiyon tekerleği ise daha yüksek nihai doğruluk sağlamıştır.
• Sabit geri besleme kazançlı karşılaştırma, momentum anahtarının kendi başına nihai doğruluğu yaklaşık iki kat iyileştirebildiğini göstermektedir.
• GSR yönlendirme, incelenen tekilliğe yakın manevrada pseudo-inverse yaklaşımına göre daha sağlam davranmıştır.
• Tek CMG arızası simülasyonunda sistem hedefe ulaşmış ancak manevra yaklaşık 3,4 kat yavaşlamış ve büyük geçici yaw sapması oluşmuştur.
• 5000 rpm'lik hassas yönelme evresi, 10.000 rpm evresine göre dengesizlik kaynaklı senkron mikro-titreşim kuvvetini model gereği yaklaşık dört kat azaltmaktadır.

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

• 25,3 saniyelik 90° manevra ve 0,42° doğruluk fiziksel bir 3U uyduda veya hava yataklı tam sistem testinde deneysel olarak gösterilmemiştir.
• Fiziksel prototip simülasyonda kullanılan uçuş ölçekli açısal momentumun tamamını üretmemektedir; prototip daha düşük momentum seviyesinde hız ve gimbal zarfını doğrulamaktadır.
• Lyapunov analizi tüm iki-aşamalı anahtarlamalı sistem için biçimsel kararlılık kanıtı değildir.
• Çalışma yalnız tek eksenli 90° büyük açılı ana manevrayı ayrıntılı olarak göstermektedir; keyfî eigenaxis manevraları ve ardışık çok hedefli görevler henüz doğrulanmamıştır.
• Birden fazla CMG arızası için kapsamlı fault detection, isolation and recovery sistemi test edilmemiştir.
• Prototipte kullanılan COTS motorlar uzay kalifikasyonlu değildir; radyasyon, termal vakum, uzun ömürlü rulman/lubrikasyon ve elektromanyetik uyumluluk doğrulaması yapılmamıştır.
• Mikro-titreşim değerleri entegre optik payload üzerinde ölçülmüş jitter değildir.
• Türkiye'ye özgü bir CubeSat, görev veya yörünge altyapısı üzerinde test yapılmadığından Türkiye için doğrudan performans sonucu çıkarılamaz.

Kaynak ve Yöntem Notu

Kaynak türü ve yayın durumu: Bu çalışma hakemli bir araştırma makalesidir ve preprint değildir. Aerospace dergisinde 2026 yılında cilt 13, sayı 7, makale 582 olarak yayımlanmıştır.

Özgün çalışma: Kenta Endo, Manami Kanamaru ve Keita Tanaka, “A Two-Step Variable-Speed Control Moment Gyroscope Control Strategy for 3U Nanosatellite Attitude Maneuvers”. DOI: 10.3390/aerospace13070582.

Sorumlu yazar: Keita Tanaka.

Kurumlar: Kenta Endo ve Keita Tanaka, Tokyo Denki University Graduate School of Science and Engineering; Manami Kanamaru, Shibaura Institute of Technology College of Engineering bünyesindedir.

Çalışmanın gerçek kanıt düzeyi: Uzay aracı düzeyindeki tutum kontrol sonuçları kapalı çevrim sayısal simülasyona dayanmaktadır. Fiziksel prototip deneyleri gimbal motor modelinin tanımlanması, tekerlek hız geçişinin ölçülmesi ve gerekli gimbal/tekerlek hız zarfının gerçekleştirilebilirliğinin sınanmasıyla sınırlıdır. Çalışma sistem düzeyinde deneysel tutum kontrol gösterimi değildir.

Donanım: Dört CMG'den oluşan piramit modülü 90 × 90 × 105 mm boyutlarında ve 584 g kütlededir. Fiziksel prototip tekerleği 6,06 g kütleli, 13,4 mm çaplı pirinç flywheel'dir. Tekerlek atalet momenti yaklaşık \(1.35\times10^{-7}\,\mathrm{kg\,m^2}\)'dir.

Deneysel model tanımlama: Gimbal sistemi 0,1–100 Hz swept-sine girişle karakterize edilmiş ve ikinci dereceden dinamik model hold-out segmentinde %93,2 uyum sağlamıştır. Tekerlek için 10.000 → 5000 rpm hız geçişi doğrudan tezgâhta ölçülmüş ve simülasyon dinamiğine veri olarak aktarılmıştır.

Ana simülasyon: 560 km yörünge irtifasında aerodinamik ve yerçekimi gradyanı bozucuları altında 90° roll manevrası gerçekleştirilmiştir. İki aşamalı stratejide momentum seviyeleri \(h_1=1.41\times10^{-3}\) ve \(h_2=7.05\times10^{-4}\,\mathrm{N\,m\,s}\), anahtarlama eşiği 5°'dir.

Kararlılık sınırı: Kaynağın özetinde “Lyapunov-based stability guarantee” ifadesi yer almakla birlikte Bölüm 2.2.3, Lyapunov analizinin yalnız sabit momentumlu nominal çevrim, ideal tork takibi ve sıfır bozucu varsayımlarına uygulandığını açıkça belirtmektedir. Yazarlar aynı bölümde bu analizin anahtarlamalı iki-aşamalı sistem için biçimsel kararlılık kanıtı olmadığını ve ortak Lyapunov fonksiyonu veya average-dwell-time türü bir analizin gelecekte yapılması gerektiğini söylemektedir. Verianla makalesinde daha geniş bir kararlılık garantisi çıkarılmamıştır.

Kaynak içi sayısal ifade tutarsızlığı: Bölüm 2.3'te 5° anahtarlama eşiğinin yüksek momentumlu Faz 1'e “more than 94° of the 90° maneuver” bıraktığı yazılmıştır. 90° toplam manevrada 94°'den fazla bölüm fiziksel olarak mümkün değildir. İfadenin olası amaçlanan biçimi tahmin edilmemiş ve sessizce düzeltilmemiştir.

Model karşılaştırması uyarısı: Ana tam modelde nihai doğruluk 0,42° ve \(m_{\min}\approx0.07\) iken eşik duyarlılığı ile Monte Carlo taramaları daha hızlı bir indirgenmiş model kullanmaktadır. Bu modelde örneğin eşik taraması doğruluğu 0,30° ve \(m_{\min}=1.09\)'dur. Kaynak bu farkı açıkça belirtmekte ve indirgenmiş model sonuçlarının mutlak performans tahmini yerine göreli eğilim/sağlamlık göstergesi olarak okunmasını istemektedir.

12× karşılaştırmasına ilişkin yorum: Kaynağın 5,2° → 0,42° karşılaştırmasında geri besleme kazançları aynı değildir. Makale bu karıştırıcı etkiyi ayrıca inceleyerek sabit kazançlı deneyde iki-aşamalı momentum anahtarının 0,72° → 0,36° düzeyinde yaklaşık iki kat doğruluk iyileştirmesi sağladığını raporlamaktadır. Bu nedenle “12 kat” değeri momentum anahtarının izole etkisi olarak sunulmamıştır.

Finansman: Araştırma kısmen Tokyo Denki University ve Tokyo Metropolitan Government University Startup Support Program tarafından finanse edilmiştir.

Veri erişilebilirliği: Çalışmada kullanılan verilerin sorumlu yazardan talep üzerine edinilebileceği, devam eden proje geliştirmeleri nedeniyle kamuya açık olmadığı bildirilmiştir.

Proje: Çalışmanın HATOSAT projesinin bir parçası olarak yürütüldüğü belirtilmiştir.

Yapay zekâ kullanımı: Yazarlar makalenin dil düzenlemesinde yapay zekâ destekli araçlar kullandıklarını ve içeriğin tüm sorumluluğunu üstlendiklerini açıklamıştır.

Çıkar çatışması: Yazarlar çıkar çatışması bildirmemiştir.

Lisans: Makale Creative Commons Attribution (CC BY) lisansı altında açık erişim olarak yayımlanmıştır.

Resmî kaynak: https://doi.org/10.3390/aerospace13070582