AIDT 온디바이스 STEM 광섬유 브래그 격자(FBG) 파장 복조 구조 건전성 감시 서보 엔진

코어 굴절률 주기 변조, 브래그 파장 \(\lambda_B = 2 n_{eff} \Lambda\), 광탄성-열광학 변형률 복조 10,000Hz C# 엔진

1. 개요

광섬유 브래그 격자(FBG, Fiber Bragg Grating) 센서는 광섬유 코어 내부에 자외선 간섭 패턴을 조사하여 주기적인 굴절률 변조 격자를 형성함으로써 특정 브래그 파장(\(\lambda_B\))의 빛만을 좁은 대역폭으로 반사시키는 광소자입니다. 외부 인장 응력이나 온도가 가해지면 격자 주기(\(\Lambda\))와 코어 유효 굴절률(\(n_{eff}\))이 변하여 반사 파장이 편이되므로, 전자기 간섭(EMI)이 극심한 항공우주 및 국방 복합소재 구조물의 변형률을 계측하는 데 최적의 신뢰성을 발휘합니다. 본 고에서는 광탄성 및 열광학 효과를 분리하여 마이크로스트레인(\(\mu\epsilon\))을 고속 복조하는 10,000Hz 0-Alloc C# 서보 엔진을 구현합니다.

2. 수학적 모델링 및 지배 방정식

기본 브래그 반사 파장 조건은 다음과 같습니다:

\[ \lambda_B = 2 n_{eff} \Lambda \]

종방향 기계적 변형률 \(\epsilon_z\)와 온도 변화 \(\Delta T\)가 동시에 작용할 때, 브래그 파장의 상대적 편이 \(\frac{\Delta \lambda_B}{\lambda_B}\)는 다음과 같이 정식화됩니다:

\[ \frac{\Delta \lambda_B}{\lambda_B} = (1 - p_e) \epsilon_z + (\alpha_\Lambda + \zeta_n) \Delta T \]

여기서 \(p_e\)는 유리카 코어의 유효 광탄성 계수(Photoelastic Coefficient)이며, 다음과 같이 포아송비 \(\nu\)와 광탄성 텐서 성분들로 결정됩니다:

\[ p_e = \frac{n_{eff}^2}{2} [p_{12} - \nu(p_{11} + p_{12})] \approx 0.22 \]

\(\alpha_\Lambda = \frac{1}{\Lambda}\frac{\partial\Lambda}{\partial T} \approx 0.55 \times 10^{-6}/\text{K}\)는 열팽창 계수, \(\zeta_n = \frac{1}{n_{eff}}\frac{\partial n_{eff}}{\partial T} \approx 6.7 \times 10^{-6}/\text{K}\)는 열광학 계수입니다. 기준 FBG 센서를 통한 온도 보상 후 순수 기계적 변형률 \(\epsilon_z\)를 정확히 추출합니다.

3. 실시간 피드백 제어 아키텍처

graph TD A[광대역 ASE 광원 광섬유 인입] --> B[항공기 날개 내장 FBG 격자 어레이 입사] B --> C[변형률 및 온도에 따른 특정 파장 브래그 협대역 반사] C --> D[파장 가변 파브리-페로 필터 F-P 광검출기 고속 스캔] D --> E[반사 피크 검출 및 가우시안 곡선 피팅 서브피코미터 파장 추출] E --> F[온도 센서 채널과 감산 연산 열팽창 간섭 성분 완전 배제] F --> G[복합재 날개 실시간 마이크로스트레인 및 비틀림 하중 산출] G --> H[10,000Hz 0-Alloc C# 구조 한계 하중 경보 및 능동 플랩 제동 서보] H --> B

4. 고성능 0-Alloc C# 서보 엔진 구현

GC 멈춤 없이 10,000Hz 주기로 다중 FBG 채널의 파장 이동을 추적하는 C# 구조입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.FiberBraggGrating
{
    public readonly record struct FbgConfig(
        double BaseBraggWavelengthNm,   // 기준 파장 lambda_B (예: 1550.0 nm)
        double PhotoelasticConst,       // 유효 광탄성 상수 pe (예: 0.22)
        double ThermalCoeffPerKelvin,   // 온도 민감도 (alpha + zeta) (예: 7.25e-6 / K)
        double MaxSafeStrainMicroStrain // 구조 안전 한계 변형률 (예: 3000.0 uE)
    );

    public ref struct FbgServoState
    {
        public double CurrentPeakWavelengthNm;
        public double MechanicalStrainMicroStrain;
        public double TemperatureCompensatedShiftNm;
        public double StressMegapascal;
        public bool IsStructuralYieldWarning;
    }

    public static class FbgInterrogationServo
    {
        private const double YoungsModulusGpa = 72.0; // 실리카 광섬유 영률 (약 72 GPa)

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void DemodulateStrain(
            in FbgConfig cfg,
            double measuredPeakWavelengthNm,
            double ambientTempDeltaK,
            ref FbgServoState state)
        {
            state.CurrentPeakWavelengthNm = measuredPeakWavelengthNm;

            // 1. 순수 총 파장 편이 계산
            double deltaLambda = measuredPeakWavelengthNm - cfg.BaseBraggWavelengthNm;

            // 2. 온도 변화에 따른 파장 이동량 계산 및 차감 보상
            // deltaLambda_T = lambda_B * (alpha + zeta) * deltaT
            double thermalShift = cfg.BaseBraggWavelengthNm * cfg.ThermalCoeffPerKelvin * ambientTempDeltaK;
            double mechanicalShift = deltaLambda - thermalShift;
            state.TemperatureCompensatedShiftNm = mechanicalShift;

            // 3. 기계적 변형률 복조 (Microstrain 단위: 1e-6)
            // deltaLambda_M / lambda_B = (1 - pe) * epsilon
            double strainFraction = mechanicalShift / (cfg.BaseBraggWavelengthNm * (1.0 - cfg.PhotoelasticConst));
            double microstrain = strainFraction * 1e6;
            state.MechanicalStrainMicroStrain = microstrain;

            // 4. 응력 계산 (Hooke 법칙 sigma = E * epsilon)
            state.StressMegapascal = (YoungsModulusGpa * 1000.0) * strainFraction;

            // 5. 구조 건전성 한계 경보 판정
            state.IsStructuralYieldWarning = Math.Abs(microstrain) > cfg.MaxSafeStrainMicroStrain;
        }
    }
}

5. 국방 및 교육 현장 응용

  • 국방 스텔스 무인전투기(UCAV) 및 잠수함 압력선체 구조 건전성 감시(SHM): 고G 기동 시 복합소재 주익 변형률 실시간 모니터링 및 심해 잠항 잠수함 선체 좌굴(Buckling) 징후 사전 경보.
  • 2022 개정 물리 및 융합과학 교과 교구: 광섬유 전반사, 브래그 회절 조건, 고체 재료역학의 응력-변형률 선도를 가상 센서 모듈로 직관적 실험 실습.

6. 결론

본 광섬유 브래그 격자(FBG) 파장 복조 서보 엔진은 10,000Hz 0-Alloc C# 연산을 통해 서브 피코미터의 파장 편이를 실시간 복조하여 항공 및 국방 플랫폼의 극한 구조 건전성 감시를 완벽히 수행합니다.