대한민국 해병대 고속상륙정 솔개(LSF-II) 155톤 차체 해상상태 3 파랑 충격 흡수 세그먼트 스커트(Skirt) 4분면 챔버 압력 능동 배분 및 1,000Hz 0-Alloc C# 고속 활주 복원 엔진

1. 에어쿠션 챔버 압축성 유체 역학과 압력 상태방정식

체적 $V_c(t)$를 갖는 에어쿠션 챔버 내부의 공기 질량 보존 및 이상기체 단열 압축/팽창 과정은 다음과 같은 비선형 미분방정식에 의해 지배됩니다:

여기서 $\gamma = 1.4$는 비열비, $\rho_a$는 공기 밀도, $\dot{m}_{in}$은 리프트 팬 유입 질량 유량입니다. 스커트 하단과 해수면 사이의 틈새 면적 $A_{gap}(t)$를 통해 외부로 방출되는 오리피스 유출 유량 $\dot{m}_{out}$은 베르누이 방정식에 의해 도출됩니다:

선수가 파도와 조우하여 틈새 높이가 0으로 수축($h_{cushion} \le \eta_{wave}$)하면 유출이 차단되면서 챔버 내부 압력 $P_c$가 $15\,\text{kPa}$ 이상으로 급상승하여 강력한 유체 슬래밍 반력이 형성됩니다.

2. 4분면 챔버 압력 차등 배분 및 선수 박힘(Bow Plow-in) 복원 모멘트

LSF-II 차체 무게중심(CG)을 기준으로 4분면 챔버(1: 좌선수, 2: 우선수, 3: 좌선미, 4: 우선미)에 작용하는 압력 적분에 의한 합성 피치($M_{\text{pitch}}$) 및 롤($M_{\text{roll}}$) 복원 모멘트는 다음과 같습니다:

선수 박힘(Bow Plow-in) 위험도 지수 $PHI$는 전방 침하량과 피치 각속도의 가중치 합으로 정의됩니다:

$PHI \ge 0.75$ 감지 시 선수측 셔터 밸브를 $100\%$ 급속 개방하여 선수 챔버 압력을 $30\%$ 가압함으로써 차체 피치업 복원 모멘트를 형성하여 전복을 선제 차단합니다.

3. 1,000Hz 0-Alloc C# Lsf2SkirtPressureEngine.cs 구현

아래의 C# 구현체는 4분면 스커트 내부 압력 센서, 초음파 해수면 틈새 센서, 함정 6자유도 IMU 자이로를 1ms 주기로 수신하여 가비지 컬렉션(GC) 0% 무할당으로 4채널 셔터 댐퍼 개도율을 산출하는 실시간 온디바이스 엔진입니다.

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

namespace ZeliDesk.Defense.MarinesHovercraft
{
    /// <summary>
    /// LSF-II 공기부양정 4분면 센서 텔레메트리
    /// </summary>
    [StructLayout(LayoutKind.Sequential)]
    public struct Lsf2HovercraftTelemetry
    {
        // 4분면(1:좌선수, 2:우선수, 3:좌선미, 4:우선미) 스커트 챔버 실측 압력 [kPa]
        public Vector4 ChamberPressuresKpa;
        // 4분면 해수면 틈새 높이 [m]
        public Vector4 GapHeightsM;
        public Vector3 HullAnglesDeg;       // Roll, Pitch, Yaw [deg]
        public Vector3 HullAngularRateDps;  // Roll, Pitch, Yaw 각속도 [deg/s]
        public float HovercraftSpeedKnots;  // 대지/대수 속력 [knots]
    }

    /// <summary>
    /// 1,000Hz 4분면 셔터 댐퍼 및 자세 복원 출력 지령
    /// </summary>
    [StructLayout(LayoutKind.Sequential)]
    public struct Lsf2DamperCommand
    {
        public Vector4 ShutterOpenings;     // 4분면 셔터 밸브 개도율 [0.0 ~ 1.0]
        public float PitchRestoringMomentKn;// 피치 복원 모멘트 [kN.m]
        public float RollRestoringMomentKn; // 롤 복원 모멘트 [kN.m]
        public float PlowInHazardIndex;     // 선수 박힘 위험도 PHI [0.0 ~ 1.0]
        public bool PlowInWarning;          // 선수 박힘 긴급 경보 플래그
    }

    /// <summary>
    /// 1,000Hz 무할당(Zero-Alloc) 솔개-II 스커트 압력 배분 및 자세 복원 엔진
    /// </summary>
    public ref struct Lsf2SkirtPressureEngine
    {
        private const float NominalCushionPressureKpa = 8.5f; // 정격 쿠션압 8.5 kPa
        private const float NominalHeightM = 1.2f;            // 정격 부양고 1.2m
        private const float ChamberAreaM2 = 45.0f;            // 4분면 개별 면적 약 45 m^2
        private const float ArmXMeters = 7.5f;                // 무게중심-전후 챔버 모멘트 암
        private const float ArmYMeters = 3.2f;                // 무게중심-좌우 챔버 모멘트 암

        private float _filteredPhi;

        public Lsf2SkirtPressureEngine(float initialPhi = 0.0f)
        {
            _filteredPhi = initialPhi;
        }

        /// <summary>
        /// 1,000Hz (dt = 0.001s) 주기 4분면 챔버 압력 능동 배분 루프
        /// </summary>
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void Update(in Lsf2HovercraftTelemetry telemetry, ref Lsf2DamperCommand cmd)
        {
            // 1. 선수 박힘(Bow Plow-in) 위험도 PHI 계산
            float minBowGap = MathF.Min(telemetry.GapHeightsM.X, telemetry.GapHeightsM.Y);
            float bowSinkage = MathF.Max(0.0f, NominalHeightM - minBowGap);
            float pitchRate = MathF.Abs(telemetry.HullAngularRateDps.Y);

            float rawPhi = (bowSinkage / NominalHeightM) + (pitchRate * 0.04f);
            _filteredPhi += (rawPhi - _filteredPhi) * 0.05f;
            cmd.PlowInHazardIndex = Math.Clamp(_filteredPhi, 0.0f, 1.0f);
            cmd.PlowInWarning = (cmd.PlowInHazardIndex >= 0.70f);

            // 2. 피치 및 롤 복원 요구 가중치 산출
            // 피치각(앞으로 기울어지면 음수 가정): 피치업을 위해 선수 챔버 가압 필요
            float pitchErr = -telemetry.HullAnglesDeg.Y; // 앞 숙임 시 양수 보정 요구
            float rollErr = telemetry.HullAnglesDeg.X;   // 우측 기울임 시 좌측 감압/우측 가압 요구

            // 3. 4분면 셔터 개도율 분배 (기본 0.70 기준)
            float baseOpening = 0.70f;
            float bowBoost = Math.Clamp(pitchErr * 0.05f + (cmd.PlowInWarning ? 0.25f : 0.0f), -0.3f, 0.3f);
            float sternBoost = -bowBoost * 0.7f;
            float rollDiff = Math.Clamp(rollErr * 0.04f, -0.2f, 0.2f);

            // 1: 좌선수, 2: 우선수, 3: 좌선미, 4: 우선미
            float v1 = Math.Clamp(baseOpening + bowBoost - rollDiff, 0.1f, 1.0f);
            float v2 = Math.Clamp(baseOpening + bowBoost + rollDiff, 0.1f, 1.0f);
            float v3 = Math.Clamp(baseOpening + sternBoost - rollDiff, 0.1f, 1.0f);
            float v4 = Math.Clamp(baseOpening + sternBoost + rollDiff, 0.1f, 1.0f);

            cmd.ShutterOpenings = new Vector4(v1, v2, v3, v4);

            // 4. 합성 복원 모멘트 계산
            // M_pitch = ( (P1+P2) - (P3+P4) ) * Area * ArmX
            float fBow = (telemetry.ChamberPressuresKpa.X + telemetry.ChamberPressuresKpa.Y) * ChamberAreaM2;
            float fStern = (telemetry.ChamberPressuresKpa.Z + telemetry.ChamberPressuresKpa.W) * ChamberAreaM2;
            cmd.PitchRestoringMomentKn = (fBow - fStern) * ArmXMeters;

            // M_roll = ( (P2+P4) - (P1+P3) ) * Area * ArmY
            float fRight = (telemetry.ChamberPressuresKpa.Y + telemetry.ChamberPressuresKpa.W) * ChamberAreaM2;
            float fLeft = (telemetry.ChamberPressuresKpa.X + telemetry.ChamberPressuresKpa.Z) * ChamberAreaM2;
            cmd.RollRestoringMomentKn = (fRight - fLeft) * ArmYMeters;
        }
    }
}

4. 2022 개정 물리·지구과학·정보 교육과정 연계표