Determining stress-based bending mode limits for the Vera C. Rubin Observatory M1M3 active mirror system

기술 요약: Vera C. Rubin 천문대의 M1M3 능동 거울 시스템을 위한 응력 기반 굽힘 모드 한계 결정

문제 정의
Vera C. Rubin 천문대의 주거울-삼차거울(M1M3)은 156개의 공압 액추에이터로 지지되는 8.4미터 크기의 단일 구조 보로실리케이트 광학 부품입니다. 운영 중 능동 광학 시스템(AOS)은 중력, 열 구배 및 기계적 교란으로 인한 광학 수차를 보정하기 위해 액추에이터 힘을 지속적으로 조정합니다. 이러한 보정은 파면 측정으로부터 도출되지만, 결과로 발생하는 힘 구성이 유리 기판에 과도한 기계적 응력을 유발하지 않도록 보장하는 것이 매우 중요합니다. 거울은 17톤에 달하는 대체 불가능한 부품이며, 유리의 안전성은 최우선 제약 조건입니다.

표준 유한 요소 분석(FEA)은 정확한 응력 예측을 제공하지만, 상당한 계산 시간이 소요되어 텔레스코프 운영에 필요한 실시간 제어 결정에는 적합하지 않습니다. 따라서 모든 반복 계산마다 전체 FEA 시뮬레이션을 수행하지 않고도 후보 보정 벡터에 대한 거울 응력을 추정할 수 있는 신속한 예측 방법이 필요합니다.

방법론
본 논문은 사전 계산된 단위 굽힘 모드 응력의 제곱합 제곱근(RSS) 조합을 기반으로 한 신속한 응력 추정 방법론을 제안합니다. 이 접근 방식은 유리 기판의 선형 탄성(응력은 변형률에 비례하고, 변형률은 변위에 비례함)에 의존합니다. 굽힘 모드 변위는 변위 공간에서 직교하므로, 저자들은 이러한 모드들의 피크 응력 또한 RSS를 통해 결합될 것이라고 이론을 세웠습니다.

구현 워크플로우는 다음과 같습니다:

1. FEA 모델링: M1M3 거울 및 능동 지지 시스템에 대한 상세한 NASTRAN 모델을 개발했습니다. 처음 20개의 탄성 굽힘 모드(첫 6개의 강체 운동 제외)를 추출했습니다.
2. 단위 응력 계산: 각 20개 굽힘 모드에 대해, 1 µm RMS 광학 표면 변위를 달성하기 위해 정규화된 액추에이터 힘을 가하는 정적 분석을 수행했습니다. 각 모드의 양(+)과 음(-)의 방향에 대해 피크 주 주응력(인장) 및 부 주응력(압축)을 추출하여, 40개의 단위 응력 값으로 구성된 사전 계산된 룩업 테이블(LUT)을 생성했습니다.
3. RSS 예측 공식: 임의의 굽힘 모드 보정 계수 $\{c_i\}$에 대하여, 예측된 피크 유리 인장 응력($\sigma_{RSS}$)은 다음과 같이 계산됩니다:
   $$\sigma_{RSS} = \sqrt{\sum_{i=1}^{20} (c_i \cdot \sigma_{unit, i})^2}$$
   여기서 $\sigma_{unit, i}$는 계수 $c_i$의 부호에 따라 선택됩니다.
4. 안전 계수: 공간적으로 상관관계가 있는 모드들이 RSS(제곱합 방식)가 과소평가할 수 있는 선형 응력 합산을 일으킬 수 있는 경우를 다루기 위해, 1.25의 상수 곱셈 안전 계수를 제안했습니다.

주요 결과 및 검증
RSS 기반 예측은 다양한 굽힘 모드 조합에 대해 직접적인 NASTRAN FEA 시뮬레이션과 비교 검증되었습니다:

• 전체 모드 시험: 20개 모드가 균일한 진폭으로 동시에 활성화되었을 때, RSS 방법은 FEA 결과와 +3% 이내의 오차로 피크 응력을 예측했습니다. 반면, 직접적인 선형 응력 합산은 FEA 결과를 239~241% 과대평가했습니다.
• 실제 다중 모드 시험: 10가지의 현실적인 보정 시나리오(단일 지배 모드부터 밀집된 다중 모드 조합까지)를 테스트했습니다.
  • 10개 시험 중 7개에서 RSS 예측이 FEA와 ±5% 이내로 일치했습니다.
  • 단일 지배 보정(예: Trial 1)의 경우, 예측이 단순 스케일링으로 환원되어 거의 완벽한 일치(-1%)를 보였습니다.
  • 과소평가는 상당한 계수와 중첩된 공간 응력 패턴을 가진 중간 차수 모드가 포함된 경우 발생했습니다(예: Trial 5: -23%; Trial 10: -19%).
  • 과대평가는 공간적 상쇄가 FEA에는 반영되지만 RSS에는 반영되지 않는 희소하고 낮은 진폭의 시험에서 발생했습니다(예: Trial 3: +45%).
• 안전 계수의 효능: RSS 예측에 1.25배 안전 계수를 적용함으로써, 모든 테스트 케이스에 대해 보수적인 경계를 확보하였고, 이를 통해 과소평가 위험을 효과적으로 완화하면서도 직접 합산보다 훨씬 더 높은 정확도를 유지했습니다.

의의 및 운영 통합
본 논문은 이 RSS 기반 방법론이 M1M3 능동 광학 시스템에 필수적인 즉각적인 응력 마진 평가를 가능하게 한다고 주장합니다. 이 연구의 의의는 다음과 같은 세 가지 구체적인 운영 능력에 있습니다:

1. 실시간 안전 점검: 새로운 보정 벡터가 생성된 직후, 액추에이터 명령이 전달되기 전에 응력 추정치를 즉시 계산할 수 있습니다. 만약 응력이 운영 한계를 초과하면, 해당 보정을 표시하거나 거부하여 유리 파손을 방지할 수 있습니다.
2. 응력 제약 최적화: 이 방법은 최적화 알고리즘이 액추에이터 힘 최적화 루프에 응력 제약을 직접 포함할 수 있게 합니다. RSS 평가는 FEA에 비해 계산량이 무시할 수 있는 수준이므로, AOS는 반복적인 FEA 호출 없이도 광학 형상 요구 사항과 유리 응력 제한을 모두 만족하는 힘 솔루션을 반복적으로 찾을 수 있습니다.
3. LUT 설계 및 트레이드오프 분석: 사전 계산된 단위 응력 LUT는 응력 대비 성능 트레이드 곡선 구축과 다양한 텔레스코프 지향 구성에 따른 굽힘 모드의 진폭 한계 설정을 지원합니다.

저자들은 직접적인 모달 합산(101~302%의 오차로 인해 운영상 불가능)은 작업 불가능한 반면, 적절한 안전 계수로 보정된 RSS 방법은 물리적으로 의미 있고 계산 효율적인 실시간 응력 모니터링 접근 방식을 제공한다고 결론지었습니다. 향후 과제는 이 프레임워크를 M2 보조 거울에 적용하고, 실제 Rubin 천문대 운영에서 얻은 굽힘 모드 데이터를 활용하는 것입니다.