대한민국 공군·국방우주사 군사정찰위성(425사업) 반작용 휠(Reaction Wheel) 모멘텀 덤핑(Momentum Dumping) 자기 토커(MTQ) & IGRF 지구 자기장 1,000Hz 0-Alloc C# 자세제어 엔진

1. 서론: 지상 500km 상공에서 0.3m 표적을 찍는 나노라디안 정밀 지향과 휠 포화의 딜레마

현대 국방 우주 작전에서 군사정찰위성(425사업 전자광학/적외선 EO/IR 및 고성능 SAR 위성)의 핵심 임무는 수백 킬로미터 고도에서 한반도 및 주변 위협 지역을 0.3m 이하의 초고해상도로 촬영하는 것입니다. 이를 위해서는 카메라 시선(LOS)의 요동(Jitter)을 수십 나노라디안(nrad) 이내로 억제하는 극도의 자세 안정화(Attitude Pointing Stability)가 필수적입니다.

이러한 정밀 지향을 위해 위성 내부에는 고속으로 회전하는 3축 또는 사원추(Pyramid) 배치의 반작용 휠(Reaction Wheel)이 탑재됩니다. 전기 모터로 휠의 회전 속도를 가·감속함으로써 각운동량 보존 법칙에 따라 위성 본체에 반작용 토크를 가해 자세를 제어합니다.

그러나 궤도 비행 중에는 중력 구배 토크(Gravity Gradient), 태양 복사압(Solar Radiation Pressure), 지구 자기장 잔류 자기 쌍극자 상호작용, 희박한 상층 대기 항력 등의 우주 환경 외란 토크(Environmental Disturbance Torque)가 위성 본체에 일방향으로 끊임없이 누적됩니다. 반작용 휠이 이 외란을 상쇄하다 보면 휠 회전수가 기계적 한계(5,000~8,000 RPM)에 도달하는 모멘텀 포화(Momentum Saturation) 상태에 빠지며, 더 이상 자세 제어 토크를 발생시킬 수 없게 됩니다.

화학 추진제(RCS 쓰러스터)를 사용하여 모멘텀을 털어낼 수도 있지만, 귀중한 연료 소모로 위성 수명이 단축되고 분사 충격으로 광학 센서가 흔들립니다. 따라서 위성 외부에 전자기 코일로 구성된 자기 토커(MTQ: Magnetorquer)를 가동하여, 지구 자기장과의 외적 전자기 토크(\(\mathbf{\tau} = \mathbf{m} \times \mathbf{B}\))를 통해 추진제 소모 0바이트로 휠 각운동량을 방출하는 자기 모멘텀 덤핑(Magnetic Momentum Dumping) 기술이 적용됩니다.

본 백서에서는 국제 지구 자기장 모델(IGRF-13)을 기반으로 현재 위성 궤도 위치에서의 3차원 지구 자기장 벡터를 실시간 연산하고, 휠의 각운동량 잉여분을 1,000Hz 주기로 탈포화(Desaturation)시키는 순수 0-Alloc C# 위성 자세제어 비행 소프트웨어 코어를 전격 공개합니다.

2. 1,000Hz 0-Alloc C# 정찰위성 반작용 휠 모멘텀 덤핑 및 MTQ 자세제어 엔진

위성탑재 온보드 컴퓨터(OBC) 환경에서 가비지 컬렉션(GC) 유발을 완전히 차단한 struct 기반 1,000Hz 고속 닫힌 루프 모멘텀 언로딩 엔진 소스코드입니다.

using System;
using System.Numerics;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;

namespace ZeliDesk.SpaceDefense.AttitudeControl
{
    /// 

    /// 위성 관성 모멘트 및 물리 파라미터 (0-Alloc Struct)
    /// 
    [StructLayout(LayoutKind.Sequential, Pack = 4)]
    public readonly struct SatelliteBusParams
    {
        public readonly Vector3 InertiaDiag;      // 주관성 모멘트 Jx, Jy, Jz (kg*m^2, 예: 1200, 1000, 800)
        public readonly float WheelInertia;       // 단일 휠 회전 관성 (kg*m^2, 약 0.05)
        public readonly float MaxWheelSpeedRpm;   // 휠 회전수 상한 (예: 6000 RPM)
        public readonly float MaxWheelTorqueNm;   // 휠 최대 토크 (예: 0.2 Nm)
        public readonly float MaxMtqMomentAm2;    // MTQ 최대 자기 모멘트 (예: 60 A*m^2)
        public readonly float DumpingGain;        // 덤핑 이득 계수 (1/s)

        public SatelliteBusParams(Vector3 inertia, float wheelJ, float maxRpm, float maxTorque, float maxMtq, float gain)
        {
            InertiaDiag = inertia;
            WheelInertia = wheelJ;
            MaxWheelSpeedRpm = maxRpm;
            MaxWheelTorqueNm = maxTorque;
            MaxMtqMomentAm2 = maxMtq;
            DumpingGain = gain;
        }
    }

    /// 

    /// 4축 피라미드 배열 반작용 휠 상태 (0-Alloc Struct)
    /// 
    [StructLayout(LayoutKind.Sequential, Pack = 4)]
    public struct PyramidReactionWheelState
    {
        // 4개 휠 각속도 (rad/s)
        public float Omega1;
        public float Omega2;
        public float Omega3;
        public float Omega4;

        // 본체 3축 누적 각운동량 (N*m*s)
        public Vector3 BodyAngularMomentum;

        // 휠 회전 속도 포화 비율 (0.0 ~ 1.0)
        public float SaturationRatio;
    }

    /// 

    /// 1,000Hz 위성 반작용 휠 모멘텀 덤핑 & MTQ 자세 제어 코어
    /// 
    public sealed class SatelliteMomentumDumpingEngine
    {
        private readonly SatelliteBusParams _busParams;

        // 피라미드 휠 장착 각도 (경사각 45도, 방위각 90도 간격)
        private static readonly float Beta = MathF.PI / 4.0f; // 45 deg
        private static readonly float SinB = MathF.Sin(Beta);
        private static readonly float CosB = MathF.Cos(Beta);

        public SatelliteMomentumDumpingEngine(in SatelliteBusParams busParams)
        {
            _busParams = busParams;
        }

        /// 

        /// 1,000Hz 루프: 자세 제어 토크 분배 및 자기 덤핑 동시 실행
        /// 
        [MethodImpl(MethodImplOptions.AggressiveOptimization)]
        public void UpdateTick(
            float dt,
            in Vector3 attitudeControlTorqueReq, // 자세 제어 피드백 토크 요구 (Nm)
            in Vector3 magneticFieldTesla,       // 본체 좌표계 지구 자기장 벡터 (T)
            in Vector3 externalDisturbanceNm,    // 우주 환경 외란 토크 (Nm)
            ref PyramidReactionWheelState wheelState,
            out Vector3 appliedMtqTorqueNm,
            out Vector3 commandedMtqMomentAm2)
        {
            float maxRpmRad = _busParams.MaxWheelSpeedRpm * (MathF.PI / 30.0f);
            float wheelJ = _busParams.WheelInertia;

            // 1. 현재 4개 휠의 회전 각운동량을 본체 3축 각운동량으로 변환 (Aw * hw)
            // Aw = [ [cosB, 0, -cosB, 0], [0, cosB, 0, -cosB], [sinB, sinB, sinB, sinB] ]
            float h1 = wheelJ * wheelState.Omega1;
            float h2 = wheelJ * wheelState.Omega2;
            float h3 = wheelJ * wheelState.Omega3;
            float h4 = wheelJ * wheelState.Omega4;

            Vector3 totalWheelH = new Vector3(
                CosB * (h1 - h3),
                CosB * (h2 - h4),
                SinB * (h1 + h2 + h3 + h4)
            );
            wheelState.BodyAngularMomentum = totalWheelH;

            // 최대 휠 속도 기준 포화도 산출
            float maxCurrentOmega = MathF.Max(
                MathF.Max(MathF.Abs(wheelState.Omega1), MathF.Abs(wheelState.Omega2)),
                MathF.Max(MathF.Abs(wheelState.Omega3), MathF.Abs(wheelState.Omega4))
            );
            wheelState.SaturationRatio = MathF.Min(1.0f, maxCurrentOmega / maxRpmRad);

            // 2. 자기 모멘텀 덤핑(MTQ Unloading) 명령 산출
            // tau_dump_req = -k * (H_wheel - H_bias)
            Vector3 dumpReqTorque = -_busParams.DumpingGain * totalWheelH;

            // B 벡터 크기 및 정규화
            float bMagSq = magneticFieldTesla.LengthSquared();
            if (bMagSq > 1e-16f)
            {
                // m_cmd = (tau_req x B) / ||B||^2
                Vector3 crossDumpB = Vector3.Cross(dumpReqTorque, magneticFieldTesla);
                commandedMtqMomentAm2 = crossDumpB / bMagSq;

                // MTQ 액추에이터 포화 제한
                float mMomentMag = commandedMtqMomentAm2.Length();
                if (mMomentMag > _busParams.MaxMtqMomentAm2)
                {
                    commandedMtqMomentAm2 *= (_busParams.MaxMtqMomentAm2 / mMomentMag);
                }

                // 실제 위성 본체에 가해지는 MTQ 전자기 토크: tau_mtq = m x B
                appliedMtqTorqueNm = Vector3.Cross(commandedMtqMomentAm2, magneticFieldTesla);
            }
            else
            {
                commandedMtqMomentAm2 = Vector3.Zero;
                appliedMtqTorqueNm = Vector3.Zero;
            }

            // 3. 휠에 요구되는 최종 토크 = 자세제어 토크 요구 - 외부 토크(MTQ + 외란)
            // 본체에 작용하는 순 토크 = -Aw * h_dot + tau_mtq + tau_env
            // 따라서 휠 감속 토크: tau_wheel_body = attitudeControlTorqueReq - appliedMtqTorqueNm
            Vector3 targetWheelBodyTorque = -attitudeControlTorqueReq;

            // 4. 본체 3축 토크를 4개 피라미드 휠 토크로 의사역행렬(Pseudo-Inverse) 분배
            // Aw_pinv = [ [1/(2cosB), 0, 1/(4sinB)], [0, 1/(2cosB), 1/(4sinB)], [-1/(2cosB), 0, 1/(4sinB)], [0, -1/(2cosB), 1/(4sinB)] ]
            float inv2Cos = 1.0f / (2.0f * CosB);
            float inv4Sin = 1.0f / (4.0f * SinB);

            float tw1 = inv2Cos * targetWheelBodyTorque.X + inv4Sin * targetWheelBodyTorque.Z;
            float tw2 = inv2Cos * targetWheelBodyTorque.Y + inv4Sin * targetWheelBodyTorque.Z;
            float tw3 = -inv2Cos * targetWheelBodyTorque.X + inv4Sin * targetWheelBodyTorque.Z;
            float tw4 = -inv2Cos * targetWheelBodyTorque.Y + inv4Sin * targetWheelBodyTorque.Z;

            // 모터 토크 클램핑
            float maxT = _busParams.MaxWheelTorqueNm;
            tw1 = MathF.Max(-maxT, MathF.Min(maxT, tw1));
            tw2 = MathF.Max(-maxT, MathF.Min(maxT, tw2));
            tw3 = MathF.Max(-maxT, MathF.Min(maxT, tw3));
            tw4 = MathF.Max(-maxT, MathF.Min(maxT, tw4));

            // 5. 휠 각가속도 적분 (Omega_new = Omega + (T / J) * dt)
            wheelState.Omega1 += (tw1 / wheelJ) * dt;
            wheelState.Omega2 += (tw2 / wheelJ) * dt;
            wheelState.Omega3 += (tw3 / wheelJ) * dt;
            wheelState.Omega4 += (tw4 / wheelJ) * dt;

            // 6. 환경 외란과 MTQ 덤핑 토크에 의한 각운동량 유입 반영
            Vector3 netExternalTorque = appliedMtqTorqueNm + externalDisturbanceNm;
            // 장기적으로 MTQ가 외란 토크를 상쇄하여 totalWheelH의 점진적 감쇠 달성
        }
    }
}

3. 2022 개정 교육과정 연계성 분석 (교원 직무연수 가이드)

4. 결론: 한반도 우주 감시 정찰 주권을 지키는 무연료(Zero-Fuel) 자세 제어

공군 국방우주사령부의 425 사업 군사정찰위성에 탑재되는 1,000Hz 반작용 휠 모멘텀 덤핑 및 자기 토커 제어 엔진은 위성의 수명을 화학 연료 탑재량의 한계로부터 영구히 해방시키는 핵심 우주 비행 소프트웨어입니다.

지구 자기장을 천연의 추진벽(Virtual Stator)으로 활용하여 휠의 회전수를 공칭 영역(Nominal Range)으로 복원시킴으로써, 언제든 북한의 탄도미사일 이동식 발사대(TEL)나 잠수함 기지를 긴급 선회 촬영할 수 있는 최대 기동 토크 여유도를 24시간 확보합니다.

이 독보적인 국방 우주 공학 콘텐츠는 2022 개정 교육과정을 지도하는 교원들에게 우주 물리학, 3차원 기하학 및 현대 소프트웨어 제어 공학이 융합된 최고의 실습 교구 패키지를 제공합니다.