대한민국 육군 K9 자주포 155mm 52구경장 사격 반작용력 주퇴복좌기(Hydro-Pneumatic Recoil) 오리피스 감쇠 및 10,000Hz 0-Alloc C# 완충 해석 엔진

제1장: 강내 사격 충격량 및 주퇴복좌계 비선형 동역학

포탄 질량 $m_p$, 추진 장약 질량 $m_c$, 포탄 총구 이탈 속도 $v_0$, 연소가스 평균 배출 속도 $u_g$, 및 주퇴부(포신+약실) 총 질량 $M_r$인 K9 자주포의 총 사격 충격량 $I_{\text{total}}$과 자유 주퇴 속도 $V_{r0}$는 다음과 같이 유도됩니다:

제2장: 가변 오리피스 유압 감쇠 및 질소 가스 폴리트로픽 단열 압축

주퇴 피스톤 면적 $A_p$, 변위 $x_r$에 따라 가변 설계된 유효 오리피스 면적 $a_o(x_r)$, 작동유 밀도 $\rho_o$, 복좌기 초기 체적 $V_0$ 및 압력 $P_0$에 대한 유체 감쇠력과 복좌력 수식은 다음과 같습니다:

제3장: 10,000Hz 0-Alloc C# K9 주퇴복좌 완충 해석 엔진

주퇴-복좌 전 과정을 0.1ms(10,000Hz) 주기로 무할당 수치 적분하여 포탑 차체로 전달되는 충격 하중을 해석하는 순수 C# 실시간 엔진 구조체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Defense.K9Artillery
{
    /// <summary>
    /// K9 포신 사격 입력 파라미터 (0-Alloc stack-only)
    /// </summary>
    public readonly struct K9FiringEventInput
    {
        public readonly float GunElevationDeg; // 포신 사격 고각 theta (deg)
        public readonly float GasImpulseNs;    // 강내 연소가스 충격량 (N*s)
        public readonly float ShotDurationSec; // 탄자 강내 이동 시간 (s)
        public readonly float AmbientTempC;    // 외부 대기 온도 (degC)

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public K9FiringEventInput(float elevation, float impulse, float duration, float tempC)
        {
            GunElevationDeg = elevation;
            GasImpulseNs = impulse;
            ShotDurationSec = duration;
            AmbientTempC = tempC;
        }
    }

    /// <summary>
    /// 주퇴복좌 거동 시뮬레이션 출력 구조체
    /// </summary>
    public readonly struct K9RecoilStateOutput
    {
        public readonly float RecoilPositionM;     // 포신 후퇴 변위 x_r (m)
        public readonly float RecoilVelocityMps;   // 포신 후퇴 속도 v_r (m/s)
        public readonly float HydraulicBrakeForceN;// 유압 오리피스 제동력 (N)
        public readonly float RecuperatorForceN;   // 가스 복좌력 (N)
        public readonly float NetTrunnionLoadN;    // 차체/포이(Trunnion) 전달 총 하중 (N)
        public readonly bool IsBufferOverstroke;   // 완충 한계 초과 경보 플래그

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public K9RecoilStateOutput(
            float posX, float velV, float fHyd, float fRecup, float fTrunnion, bool overstroke)
        {
            RecoilPositionM = posX;
            RecoilVelocityMps = velV;
            HydraulicBrakeForceN = fHyd;
            RecuperatorForceN = fRecup;
            NetTrunnionLoadN = fTrunnion;
            IsBufferOverstroke = overstroke;
        }
    }

    /// <summary>
    /// 10,000Hz (0.1ms tick) 0-Alloc K9 주퇴복좌 해석 코어 엔진
    /// </summary>
    public ref struct K9RecoilDampingEngine
    {
        private const float RecoilMassKg = 2450.0f;     // 포신+약실 주퇴부 질량 (kg)
        private const float MaxRecoilStrokeM = 0.58f;   // 최대 허용 주퇴 한계 거리 (m)
        private const float PistonAreaM2 = 0.018f;      // 주퇴 피스톤 면적 (m^2)
        private const float FluidDensity = 870.0f;      // 유압 작동유 밀도 (kg/m^3)
        private const float CdCoeff = 0.62f;            // 오리피스 유량 계수
        private const float RecupAreaM2 = 0.009f;       // 복좌기 피스톤 면적 (m^2)
        private const float RecupP0 = 6.8e6f;           // 복좌기 질소 초기 압력 (Pa, ~68 bar)
        private const float RecupV0 = 0.012f;           // 복좌기 초기 체적 (m^3)
        private const float Gamma = 1.4f;               // 질소 가스 단열 지수

        private float _positionM;
        private float _velocityMps;
        private float _simTimeSec;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public K9RecoilStateOutput Step10000Hz(in K9FiringEventInput firing, float dt)
        {
            _simTimeSec += dt;

            // 1. 강내 사격 연소가스 추력 F_gas(t) 구형파/삼각파 모델
            float fGas = 0.0f;
            if (_simTimeSec <= firing.ShotDurationSec && firing.ShotDurationSec > 0.0f)
            {
                // 충격량 삼각형 근사 피크 추력
                fGas = (2.0f * firing.GasImpulseNs / firing.ShotDurationSec) *
                       (1.0f - (_simTimeSec / firing.ShotDurationSec));
            }

            // 2. 가변 오리피스 면적 산출 a_o(x) : 후퇴 변위에 따라 제동 압력 일정 유지
            // x가 증가할수록 오리피스 면적이 좁아지도록 설계
            float xClamped = Math.Clamp(_positionM, 0.0f, MaxRecoilStrokeM);
            float ao = MathF.Max(1.2e-4f, 8.5e-4f * (1.0f - (xClamped / MaxRecoilStrokeM) * 0.85f));

            // 3. 유압 오리피스 제동력 F_brake
            float vSign = MathF.Sign(_velocityMps);
            float vAbs = MathF.Abs(_velocityMps);
            float fHyd = 0.0f;
            if (ao > 1e-6f)
            {
                float pressureDrop = (FluidDensity / (2.0f * CdCoeff * CdCoeff)) *
                                     MathF.Pow(PistonAreaM2 / ao, 2.0f) * vAbs * vAbs;
                fHyd = PistonAreaM2 * pressureDrop * vSign;
            }

            // 4. 질소 가스 복좌력 F_recup (단열 압축)
            float vGasCurrent = MathF.Max(1e-4f, RecupV0 - RecupAreaM2 * _positionM);
            float pGas = RecupP0 * MathF.Pow(RecupV0 / vGasCurrent, Gamma);
            float fRecup = RecupAreaM2 * pGas;

            // 5. 중력 성분 및 마찰력
            float elevRad = firing.GunElevationDeg * (MathF.PI / 180.0f);
            float fGravity = RecoilMassKg * 9.80665f * MathF.Sin(elevRad);
            float fFriction = 3500.0f * vSign; // 레일 슬라이딩 마찰

            // 6. 차체 포이(Trunnion) 전달 총 하중
            float fTrunnion = fHyd + fRecup;

            // 7. 뉴턴 운동방정식 가속도 및 0.1ms 수치 적분 (Euler-Cromer)
            float netForce = fGas - fHyd - fRecup - fFriction - fGravity;
            float accel = netForce / RecoilMassKg;

            _velocityMps += accel * dt;
            _positionM += _velocityMps * dt;

            // 포신 원위치(x=0) 복귀 시 완충 정지 바운드
            if (_positionM < 0.0f)
            {
                _positionM = 0.0f;
                if (_velocityMps < 0.0f) _velocityMps = 0.0f;
            }

            bool isOverstroke = _positionM > (MaxRecoilStrokeM * 0.96f);

            return new K9RecoilStateOutput(
                _positionM,
                _velocityMps,
                fHyd,
                fRecup,
                fTrunnion,
                isOverstroke
            );
        }
    }
}

제4장: 2022 개정 교육과정 연계 교과 융합 지도안 (교사용 가이드)