기존 유한요소 코드로 주기재료 분산관계 계산하는 새로운 요소 기반 방법

초록

본 논문은 상용 유한요소(FE) 프로그램에서 전역 행렬에 접근하지 못하는 제약을 극복하고, 블록 경계조건을 요소 수준에서 직접 적용함으로써 주기재료의 분산관계를 효율적으로 계산하는 방법을 제시한다. 사용자 정의 요소 서브루틴만으로 구현 가능하며, 실수·복소수 연산 모두 지원한다. 다양한 차원·물리 현상에 적용 가능함을 확인하기 위해 여러 재료 모델과 셀 형상에 대해 해석·검증을 수행하였다.

상세 분석

이 연구는 블록 이론에 기반한 주기재료의 파동 전파 해석을 기존 상용 유한요소 코드에 통합하는 데 초점을 맞춘다. 전통적인 접근법은 전역 강성·질량 행렬에 변환 행렬을 곱하거나 행·열 연산을 수행해 블록 경계조건(BBC)을 적용하는데, 이는 상용 코드가 전역 행렬에 직접 접근할 수 없기 때문에 실현이 어렵다. 저자들은 이러한 한계를 극복하기 위해 BBC를 요소 수준에서 구현하는 전략을 제안한다. 구체적으로, 각 요소의 연결 정보(노드 집합)를 두 종류의 연결 연산자 C_c 와 C_a 로 구분한다. C_c 는 기존과 동일하게 요소 행렬을 계산하는 데 사용되고, C_a 는 어셈블 과정에서 BBC가 반영된 형태로 전역 행렬에 기여하도록 설계된다. 이렇게 하면 전역 행렬을 직접 수정할 필요 없이, 어셈블 단계에서 이미 위상 변이를 포함한 축소된 행렬 K_R, M_R 이 생성된다.

이 방법의 핵심은 위상 변이 e^{i k·a} 를 요소 행렬에 직접 삽입함으로써, 동일한 물리적 셀을 복제하는 인접 셀 간의 자유도 관계를 자동으로 만족시키는 것이다. 따라서 복소수 고유값 문제를 실수 연산 기반 2‑mesh 기법이나 복소수 연산을 지원하는 솔버에 그대로 전달할 수 있다. 또한, 요소 수준에서 BBC를 적용하므로 스칼라·벡터 문제, 2‑D·3‑D 차원, 고전 탄성부터 코시어·마이크로폴라와 같은 고차원 이론까지 확장 가능하다.

검증 사례에서는 (1) 평면 파동에 대한 전통적인 탄성 매질, (2) 코시어 연성 매질, (3) 다양한 셀 기하학(직사각형, 삼각형, 3‑D 격자) 등을 대상으로 해석을 수행하였다. 각 사례에서 분석된 분산곡선은 문헌에 보고된 해석적 해와 수치 결과와 일치함을 보였으며, 셀 크기 변화에 대한 민감도 테스트에서도 안정적인 결과를 얻었다. 이러한 검증을 통해 제안된 요소 기반 BBC 적용 방법이 일반적인 상용 FE 코드(예: ABAQUS, ANSYS, FEAP‑pv)에서 구현 가능하고, 추가적인 전처리·후처리 없이도 정확한 밴드 구조를 도출할 수 있음을 입증하였다.

결과적으로, 이 접근법은 전통적인 전역 행렬 조작이 불가능한 환경에서도 주기재료의 파동 특성을 정밀히 분석할 수 있는 실용적인 도구를 제공한다. 사용자 정의 요소 서브루틴만으로 구현 가능하므로, 기존 코드의 라이선스·보안 제약을 유지하면서도 연구·산업 현장에서 광범위하게 활용될 수 있다.