Low-Thrust Orbital Trajectory Guidance and Control Using Mean Orbital Elements
تتقصى هذه الورقة تطبيق العناصر المدارية المتوسطة في صياغة قانون توجيه تكيفي اتجاهي لعمليات النقل المداري ذات الدفع المنخفض، مما يثبت متانتها وفعاليتها مقارنة بالعناصر المتذبذبة في التعامل مع مشكلات التحسين عالية الأبعاد مع الاضطرابات الطبيعية عبر أنواع مختلفة من المناورات.
المؤلفون الأصليون:Houman Hakima, Michael C.F. Bazzocchi
لطالما اعتمد السفر عبر الفضاء على دفعات قوية وقصيرة من الوقود لتحريك المركبات الفضائية من مدار إلى آخر، تماماً مثل العداء الذي ينطلق بقوة من مكعبات البداية. ومع ذلك، برزت في العقود الأخيرة طريقة أكثر هدوءاً وكفاءة: الدفع الكهربائي. هذه الأنظمة، التي تُسمى غالباً محركات الدفع المنخفض، لا تنفجر بضربة مفاجئة، بل توفر دفعاً لطيفاً ومستمراً لشهور أو حتى لسنوات. وبين أنها تستهلك وقوداً أقل بكست كثير من الصواريخ الكيميائية التقليدية، إلا أنها تتطلب أوقاتاً أطول بكثير لإتمام الرحلة. ويكمن التحدي أمام المهندسين في أن هذا الدفع البطيء والمستقر يجب أن يكون موقوتاً وموجهاً بدقة متناهية عبر آلاف المدارات. فإذا كان نظام التوجيه حساساً للغاية تجاه التمايلات الصغيرة والسريعة في مسار المركبة الفضائية الناتجة عن شكل الأرض غير المنتظم أو السحب الجوي، فقد تهدر المركبة الفضائية وقودها في تصحيح مسارها باستمرار. ويصبح السؤال: كيف يمكن توجيه مركبة يتم دفعها بهمس بدلاً من صرخة، لضمان وصولها إلى وجهتها دون استنفاد إمداداتها المحدودة من الوقود؟
في دراسة حديثة، بحث الباحثون في طريقة محددة لحل مشكلة التوجيه هذه من خلال تغيير كيفية رؤيتهم لموقع المركبة الفضائية. فبدلاً من تتبع الموقع الدقيق واللحظي للمركبة الفضائية، والذي يتذبذب باستمرار بسبب الجذب الثقالي والمقاومة الهوائية، اقترحوا استخدام متوسط سلس للمدار. وقد اختبروا طريقة توجيه تُعرف باسم "التوجيه التكيفي الاتجاهي"، والتي تحسب أفضل زاوية لتوجيه المحرك من خلال وزن أهمية الأهداف المدارية المختلفة، مثل تغيير حجم المدار أو ميله. وقارن الفريق بين نهجين: أحدهما يستخدم البيانات الخام والمتذبذبة للمسار الفعلي للمركبة الفضائية، والمعروفة بالعناصر المماسية (osculating elements)، والآخر يستخدم المسار المتوسط المحسوب، والمعروف بالعناصر المتوسطة (mean orbital elements). ومن خلال إجراء عمليات محاكاة حاسوبية مفصلة لثلاثة أنواع مختلفة من الرحلات الفضائية، سعوا لمعرفة أي طريقة ستسمح للمركبة الفضائية بالوصول إلى هدفها بكفاءة أكبر وبجهد أقل على آليات التوجيه الخاصة بها.
قام الباحثون بمحاكاة ثلاث مهام متميزة لاختبار نظريتهم. تضمنت الأولى رفع قمر صناعي من ارتفاع منخفض يبلغ 400 كيلومتر إلى ارتفاع أعلى يبلغ 1,500 كيلومتر. وحاكت الثانية رحلة معقدة من مدار انتقالي إلى مدار ثابت بالنسبة للأرض (geostationary orbit)، مما يتطلب تغييرات في حجم المدار وشكله وميله. أما الحالة الثالثة، فكانت انتقالاً إلى مدار بيضاوي عالي يُستخدم للاتصالات، حيث كان على المركبة الفضائية تعديل سرعتها وشكلها بعناية لتجنب الهبوط منخفضاً جداً في الغلاف الجوي. وفي كل سيناريو، تم نمذجة المركبة الفضائية بكتلة قدرها 500 كيلوغرام وتزويدها بمحرك دفع كهربائي عالي الكفاءة قادر على إنتاج قوة ضئيلة ولكن مستقرة. وتتبع الفريق كمية الوقود المستهلكة، والوقت الذي استغرقته المناورات، وعدد المرات التي اضطرت فيها المركبة الفضائية لتدوير جسمها لاتباع أوامر التوجيه.
أظهرت النتائج أن كلا الطريقتين يمكنهما توجيه المركبة الفضائية بنجاح إلى وجهتها، لكنهما قدمتا مقايضات مختلفة. فعندما استخدم نظام التوجيه البيانات الخام اللحظية، وصلت المركبة الفضائية إلى هدفها بشكل أسرع قليلاً واستهلكت جزءاً ضئيلاً أقل من الوقود. وهذا أمر متوقع، حيث كان النظام يتفاعل مع الحالة الحقيقية والفورية للمدار. ومع ذلك، تطلب هذا النهج من نظام التحكم في وضعية المركبة الفضائية إجراء تعديلات أكثر تكراراً وسرعة على زاوية المحرك. في المقابل، تطلب نظام التوجيه القائم على المدار المتوسط أن تغير المركبة الفضائية اتجاه توجيهها بشكل أبطأ وأكثر استقراراً. وكان الحد الأقصى لمعدل تغيير اتجاه دوران المركبة الفضائية أقل بكثير عند استخدام بيانات المدار المتوسط. على سبيل المثال، في محاكاة الرحلة إلى المدار الثابت بالنسبة للأرض، كان متوسط معدل تغير زاوية التوجيه حوالي 0.2 درجة في الثانية باستخدام الطريقة الملساء، مقارنة بقيم أعلى قليلاً باستخدام طريقة البيانات الخام.
خلصت الدراسة إلى أنه بينما يعد نهج البيانات الخام أكثر دقة من الناحية التقنية فيما يتعلق بالوقود والوقت، فإن طريقة المدار المتوسط الملساء غالباً ما تكون أكثر عملية للمهام الواقعية. فالمركبات الفضائية لديها قدرة ميكانيكية محدودة على تدوير أجسامها بسرعة، ومحاربة أوامر التوجيه السريعة والمتذبذبة باستمرار يمكن أن يؤدي إلى تآكل المكونات أو تجاوز الحدود الفيزيائية للمركبة. ومن خلال استخدام المدار المتوسط، يوفر نظام التوجيه مساراً أكثر لطفاً وقابلية للتنبؤ لتتبعه المركبة الفضائية، مما يقلل الطلب على أجهزة التوجيه الخاصة بها. ووجد الباحثون أنه بالنسبة للمدارات القريبة من خط الاستواء، حيث يسبب شكل الأرض أكبر التذبذبات، كان هذا الفرق في جهد التحكم ملحوظاً بشكل خاص. وتشير الدراسة إلى أنه بالنسبة للعديد من المهام المستقبلية، وخاصة تلك التي تشمل أقماراً صناعية صغيرة ذات قدرة محدودة على المناورة، فإن تصميم نظام التوجيه حول السلوك المتوسط للمدار بدلاً من تقلباته اللحظية يوفر حلاً قوياً وفعالاً. ويتيح هذا النهج للمهندسين التخطيط لمهام الدفع الكهربائي طويلة الأمد بثقة أكبر في أن المركبة الفضائية لن تُثقل كاهلها بالحاجة إلى تصحيح مسارها باستمرار.
ملخص تقني: توجيه والتحكم في المسارات المدارية ذات الدفع المنخفض باستخدام العناصر المدارية المتوسطة
بيان المشكلة تعد أنظمة الدفع ذات الدفع المنخفض، التي تتميز بنبض نوعي عالٍ، حيوية بشكل متزايد للمهمات التي تتطلب تغيرات كبيرة في السرعة (delta-v)، مثل نقل المدارات، والحفاظ على الموقع المداري، والسفر بين الكواكب. ومع ذلك، فإن تحسين مسارات الدفع المنخفض يمثل تحدياً عالي الأبعاد يتفاقم بسبب التسارعات الاضطرابية (مثل استدارة الأرض، والسحب الجوي) التي تمنع وجود حلول تحليلية دقيقة. تعتمد قوانين التوجيه التقليدية غالباً على العناصر المدارية المتذبذبة (osculating orbital elements)، والتي تمثل المدار "الحقيقي" اللحظي. ويجادل المؤلفون بأن العناصر المتذبذبة تكون دون المستوى الأمثل للتحكم طويل الأمد بسبب التقلبات الدورية القصيرة (مثل تأثير J2 والسحب الجوي). هذه التقلبات السريعة تستلزم تعديلات متكررة وعالية المعدل في اتجاه وضعية المركبة الفضائية للحفاظ على اتجاه الدفع المطلوب، مما قد يتجاوز قدرات نظام تحديد وتحكم الوضعية (ADCS) ويزيد من جهد التحكم.
المنهجية يبحث البحث في تطبيق العناصر المدارية المتوسطة (mean orbital elements)—وهي المعاملات التي تعمل على ترشيح التباينات الدورية القصيرة لتمثيل متوسط المدار بطيء التغير—في صياغة قانون التوجيه التكيفي الاتجاهي (DAG).
الديناميكا والاضطرابات: تستخدم الدراسة معادلات غاوس المتغيرة (GVEs) لنمذجة الحركة المدارية. وهي تأخذ في الاعتبار صراحةً اضطرابات J2 (استدارة الأرض) في اشتقاق العناصر المتوسطة، وتتضمن كلاً من J2 والسحب الجوي في المحاكاة العددية لضمان الدقة.
اشتقاق العناصر المتوسطة: يعرّف المؤلفون العناصر المتوسطة (ϵˉ) عن طريق طرح الحدود الدورية القصيرة (Δϵsp) المستمدة من اضطرابات J2 من العناصر المتذبذبة (ϵosc). وقد تم توفير تعبيرات تحليلية لهذه التباينات الدورية القصيرة لجميع العناصر المدارية الستة الكلاسيكية (نصف المحور الأكبر، الانحراف المركزي، الميل، خط طول العقدة الصاعدة، حجة الحضيض، والشذوذ الحقيقي).
صياغة قانون التوجيه: يتم استخدام قانون DAG، الذي يدمج اتجاه دفع صافٍ كجموع مرجح لاتجاهات الدفع المثلى لعناصر مدارية محددة. يتم تعريف اتجاه الدفع من خلال زوايا التوجيه داخل المستوى (α) وخارج المستوى (β).
زوايا التوجيه المثلى: تم اشتقاق الحلول التحليلية للزوايا التي تعظم معدل التغير لكل عنصر. ويتناول البحث قضايا التنفيذ العملي، مثل تصحيح قيود النطاق لدالة الظل العكسي (inverse-tangent) واختيار المسار الأقصر للعناصر الزاوية (مثل خط طول العقدة الصاعدة، وحجة الحضيض) لتقليل استهلاك الوقود.
التكيف: يتضمن القانون معاملات تكيف (Rϵ) تقوم ديناميكياً بتعديل وزن اتجاهات الدفع بناءً على الخطأ بين قيم العناصر الحالية والمستهدفة.
المحاكاة العددية: تمت محاكاة ثلاث دراسات حالة باستخدام مركبة فضائية كتلتها 500 كجم ومحرك أيوني شبكي بقوة 235 ميلي نيوتن (Isp=4,220 ثانية):
الحالة 1: من مدار أرضي منخفض (LEO) إلى مدار أرضي منخفض مرتفع (زيادة نصف المحور الأكبر).
الحالة 2: من مدار انتقال جيومداري (GTO) إلى مدار جيومداري (GEO) (تقليل متزامن لنصف المحور الأكبر، والانحراف المركزي، والميل).
الحالة 3: من مدار انتظار إلى مدار مولنيا (تغييرات متزامنة في نصف المحور الأكبر، والانحراف المركزي، وحجة الحضيض).
المساهمات الرئيسية
مقارنة العناصر المتوسطة والمتذبذبة: يقدم البحث تقييماً مقارناً لقوانين التوجيه القائمة على العناصر المتذبذبة مقابل تلك القائمة على العناصر المتوسطة في وجود الاضطرابات الطبيعية.
منطق التوجيه المنقح: يفصل البحث اشتقاق زوايا التوجيه المثلى للعناصر المتوسطة ويقدم منطقاً محدداً للتعامل مع حالات عدم الاستمرارية في زوايا التوجيه (مثل مشكلات الالتفاف في دوال الظل العكسي وتقليل الخطأ الزاوي للعناصر الدورية).
استراتيجية الوزن: توضح الدراسة مخطط وزن ديناميكي حيث يتم تأخير أولوية تصحيح عناصر معينة (مثل الانحراف المركزي) حتى يزد সাই المدار، وذلك لتجنب المرور عبر مناطق الحضيض المنخفضة التي قد تسبب سحباً جوياً مفرطاً.
النتائج تشير النتالئد العددية عبر دراسات الحالة الثلاث إلى ما يلي:
الوقود والوقت: حققت قوانين التوجيه القائمة على العناصر المتذبذبة أداءً أفضل قليلاً من حيث إجمالي استهلاك الوقود ومدة المناورة (على سبيل المثال، في الحالة 2: 43.70 كجم مقابل 43.80 كجم؛ و65.30 يوماً مقابل 65.46 يوماً). ويعزى ذلك إلى أن العناصر المتذبذبة تمثل الحالة الحقيقية اللحظية، مما يسمح بالدفع الأمثل نظرياً في كل لحظة.
جهد التحكم: أظهرت قوانين التوجيه القائمة على العناصر المتوسطة ملفات تعريف أكثر سلاسة لزوايا التوجيه. كانت معدلات التغير القصوى والمتوسطة لزوايا التوجيه أقل عند استخدام العناصر المتوسطة.
في الحالة 2 (التحول من GTO إلى GEO)، قلل نهج العنصر المتوسط من أقصى معدل تغير لزاوية التوجيه داخل المستوى للانحراف المركزي (α˙e) من 4.4525 درجة/ثانية إلى 4.4471 درجة/ثانية، ومن المعدل المتوسط من 0.6466 درجة/ثانية إلى 0.6438 درجة/ثانية.
في الحالة 3 (تحول مولنيا)، قلل نهج العنصر المتوسط من أقصى معدل تغير لزاوية حجة الحضيض داخل المستوى (α˙ω) من 13.0630 درجة/ثانية إلى 12.9823 درجة/ثانية.
التقارب: كان الفرق بين العناصر المتذبذبة والمتوسطة أكثر وضوحاً عند الارتفاعات المنخفضة حيث تكون اضطرابات J2 أقوى، ولكن الاتجاهات تقاربت عموماً مع زيادة حجم المدار، باستثناء التحولات ذات الانحراف المركزي العالي حيث تسبب سحب الحضيض تباعداً.
الأهمية والادعاءات يخلص المؤلفون إلى أنه بينما يعد التوجيه القائم على العناصر المتذبذبة أكثر كفاءة بشكل طفيف من حيث الوقود والوقت، فإن التوجيه القائم على العناصر المتوسطة هو الأكثر قابلية للتطبيق في المهمات الفضائية العملية، خاصة للمركبات الفضائية ذات سلطة التحكم في الوضعية المحدودة. إن التطور الأكثر سلاسة لزوايا التوجيه المطلوبة بواسطة توجيه العناصر المتوسطة يقلل من الطلب على نظام التحكم في وضعية المركبة الفضائية، مما يتجنب الحاجة إلى تعديلات سريعة وعالية التردد في الوضعية التي قد يتطلبها توجيه العناصر المتذبذبة بسبب التذبذبات الدورية القصيرة.
يدعي البحث بتواضع أنه بالنسبة للمدارات منخفضة الميل، فإن الفرق في جهد التحكم ليس كبيراً، ولكن بالنسبة للمدارات الاستوائية وشبه الاستوائية، فإن الفرق جوهري. ويقترح أن العمل المستقبلي يجب أن يدمج مصادر اضطراب إضافية (مثل السحب، والهارمونيات الزونية من الرتب الأعلى J3,J4، وتأثيرات الأجرام الثالثة) لزيادة دقة الحدود الدورية القصيرة المستخدمة في حسابات العناصر المتوسطة.