FE 코드 NOSA ITACA를 위한 구조 해석 향상 모델 업데이트

초록

본 논문은 역사적 석조 구조물의 동적 특성을 정확히 재현하기 위해 NOSA‑ITACA 유한요소(FE) 코드에 구현된 모델 업데이트 절차를 제시한다. 자유 파라미터를 최적화하여 실험으로 얻은 고유진동수와 모드 형태를 일치시키는 신뢰구역(trust‑region) 기반 알고리즘과 Lanczos 기반 국부 축소 모델을 결합한다. 간단한 예제와 토레 델레 오레(루카) 및 마달레나 다리(보르고) 두 실제 사례에 적용해 계산 효율성, 민감도 분석, 잡음에 대한 강인성을 검증한다.

상세 분석

논문은 구조역학에서 모델 업데이트를 역문제(inverse problem)로 정의하고, 목표함수 Φ(x)를 고유진동수와 모드 형태의 차이를 동시에 최소화하도록 설계하였다. 파라미터 벡터 x는 재료 탄성계수, 질량, 경계조건 등 p 차원의 박스 구간 Ω 안에 제한되며, 각 파라미터에 대한 하한·상한을 명시한다. 목표함수는 주파수 차이의 제곱합에 가중치 w를 곱하고, 모드 형태 차이는 MAC(Modal Assurance Criterion) 기반 정규화된 스칼라 γ_i 로 표현한다. 이를 통해 주파수와 모드 형태를 균형 있게 고려할 수 있다.

수치해법은 신뢰구역 최적화 프레임워크 안에서 수행된다. 현재 파라미터 추정점 x_k에서 Lanczos 방법을 이용해 가장 작은 q개의 고유값·고유벡터를 효율적으로 계산하고, 이를 기반으로 국부적인 파라미터 의존 축소 모델을 만든다. 축소 모델은 선형 근사(1차 정확도)를 제공하므로, 저비용으로 목적함수의 감소 방향을 탐색할 수 있다. 신뢰구역 반경은 모델의 정확도에 따라 자동으로 조정되며, 모델이 충분히 정확하면 전역 최소점으로 수렴한다는 수학적 보장이 제시된다.

민감도 분석에서는 목표함수의 잔차 r(x)와 그 제이콥 행렬 J(x)의 특이값 분해(SVD)를 수행한다. 가장 큰 특이값 σ_i에 대응하는 특이벡터는 파라미터 변화에 가장 민감한 방향을 나타내며, 작은 특이값은 거의 영향을 주지 않는 방향을 의미한다. 조건수 κ = σ_max/σ_min을 통해 입력 데이터(실험 주파수·모드)의 작은 변동이 최적 파라미터에 미치는 영향을 정량화한다. 논문은 잡음이 있는 경우에도 해가 급격히 변하지 않도록 신뢰구역과 조건수를 활용한 안정성 검증 절차를 제시한다.

실제 적용 사례에서는 간단한 2차원 프레임 모델을 먼저 테스트해 알고리즘의 수렴 속도와 계산 시간 절감을 확인하였다. 이어서 토레 델레 오레와 마달레나 다리의 복합 구조에 대해 파라미터(콘크리트 탄성계수, 석조 압축강도, 경계 스프링 등)를 선택하고, 실험으로 얻은 5~7개의 고유진동수와 대응 모드 형태를 목표함수에 포함시켰다. 결과적으로 기존 상용 최적화기와 비교했을 때, 제안된 방법은 동일한 정확도를 유지하면서 평균 60 % 이상의 계산 시간을 단축하였다. 또한, 잡음 수준을 인위적으로 증가시킨 민감도 실험에서 조건수가 크게 상승하는 경우에만 파라미터 추정이 불안정해짐을 확인하고, 이를 사용자에게 경고하는 메커니즘을 구현하였다.

전반적으로 논문은 대규모 석조 구조물의 동적 해석에 필요한 모델 업데이트를 자동화하고, 효율적인 축소 모델과 신뢰구역 기반 최적화를 결합함으로써 실무 엔지니어가 제한된 실험 데이터만으로도 신뢰성 높은 FE 모델을 구축할 수 있음을 증명한다.