맥놀이는 매우 근접한 두 주파수의 진동이 중첩될 때, 진동이 주기적으로 커졌다 작아지는 현상이다. 회전축이나 회전쉘, 풍력 타워 같은 축 대칭형 구조가 질량이나 강성의 분포 또는 경계조건에서 미세한 비대칭성을 갖는 경우, 하나의 진동모드가 비슷한 주파수를 갖는 진동모드 쌍으로 분리되고 그 간섭으로 맥놀이가 발생한다. 맥놀이가 중요한 의미를 갖는 진동계로 한국 종이 있다. 적절한 주기로 우~웅 우~웅 반복되는 맥놀이는 종소리의 웅장함과 생동감을 높인다⁽¹⁾. 한국 종에서 맥놀이 진동의 감쇠비는 종소리가 얼마나 오래 지속되는가를 결정하는 중요한 음향 성능 인자이며, 내부 균열 등의 구조적 변화를 감지하는 데에도 사용되는 물리적 인자다. 종의 진동과 같이 감쇠가 작은 진동계의 감쇠비를 구하는 데에는 대수감소법(logarithmic decrement method)이 자주 사용된다. 대수감소법은 단일 주파수의 진동 파형으로부터 점성감쇠비를 추출하는 오래된 방법이나 근래에 들어서도 감쇠비 추출 오차를 줄이려는 많은 연구가 있었다. Huang 등⁽²⁾은 부분적으로 존재하는 잡음이 대수감소에 주는 큰 오차를 제거하기 위해 면적감소(area decrement)의 개념을 도입했다. Tweten 등⁽³⁾은 잡음 오차 감소를 위해 대수감소의 평균화 과정에서 오차를 최소화시키는 최적 싸이클 수를 제시했다. Little 등⁽⁴⁾은 대수감소의 불확도를 최소화시키는 최적 싸이클 수에 대한 이론식을 제시했다. 그러나 이러한 연구들은 모두 단일 주파수 파형에 대한 연구이다. Fig. 1과 같이 매우 근접한 두 개의 주파수 운동이 간섭하여 만드는 맥놀이 파동에는 기존의 대수감소법을 바로 적용할 수가 없다. 주파수 필터링을 통해 두 주파수 성분을 완전히 분리해야 기존의 대수감소법을 적용할 수 있다. 문제는 두 주파수 성분이 매우 근접한 경우, 주파수 필터링이 개별 파형을 상당히 왜곡시킨다는 점이다. Liao and Wells⁽⁵⁾에 의하면, 필터의 대역폭은 목표 주파수 성분을 포함하도록 충분히 커야 함과 동시에 인접한 주파수 성분을 완전히 배제할 정도로 작아야 한다. 따라서 두 주파수 성분이 가까이 붙어있는 종의 맥놀이에서 필터링은 상당한 왜곡을 일으킨다. 대책으로 Nakutis 등⁽⁶⁾은 맥놀이의 포락선 데이터를 사용하는 확장된 대수감소법을 제시하였으나, 이 방법은 Fig. 1(a)와 같이 두 주파수 성분의 진폭이 같을 때 나오는 완벽하게 선명한 맥놀이에만 한정되었다. 인접한 두 주파수 성분의 진폭이 다른 일반적인 경우, (b)와 같이 흐린 맥놀이 파동을 얻게 된다. 선행연구에서 웨이블릿 변환⁽⁷⁾을 사용해서 다자유도 진동계 모드의 감쇠비를 동시에 추출하는 방법을 제시하였다. 현장에서 보다 쉽게 사용 가능한 방법으로 최근에는 일반적인 맥놀이 파형의 감쇠비 추출을 위한 포락선 대수감소법(envelope logarithmic decrement method)⁽⁸⁾를 제시하고 평창 동계올림픽대종에 적용하였다. 이 방법은 필터링 과정 없이 다양한 패턴의 맥놀이 포락선 데이터로부터 감쇠비를 추출할 수 있다. 그러나 포락선 대수감소법을 사용할 때, 측정 조건에 따라 오차가 발생할 수가 있다. 파형 측정시 사용하는 샘플링 주파수, 측정 파형의 선명도, 맥놀이를 만드는 주파수 쌍의 차이가 중요한 오차 요인이다.

Fig. 1 
Envelopes of beat waves

이 연구에서는 주어진 맥놀이 조건과 감쇠비 참값을 사용하여 맥놀이 파형을 만든 후, 위 세 가지 오차 요인에 따른 감쇠비 추출 오차를 시뮬레이션을 통해 분석한다. 그 결과로부터 포락선 대수감소법의 오차 범위를 파악하고 측정 오차를 감소시키기 위한 정보를 제공하는 것이 이 연구의 목적이다.

미세한 주파수 차이를 갖는 두 개의 점성감쇠 진동이 중첩될 때 발생하는 맥놀이 파동은 식 (1)로 표현된다.

ω_L은 미세하게 낮은 주파수 성분이고, ω_H는 미세하게 높은 주파수 성분이다. ζ_L과 ζ_H는 ω_L, ω_H 두 주파수 진동의 감쇠비다. 동일한 진동모드가 ω_L, ω_H의 매우 근접한 주파수로 분리된 경우, 두 감쇠비는 같다고 보고, 식 (1)을 식 (2)로 단순화시킬 수 있다.

여기서 ζ_L과 ζ_H는 동일한 값 ζ로 표시하였고, 평균 주파수 ω=(ω_L+ω_H)/2를 사용하였다.

식 (2)의 괄호 안의 맥놀이 파동을 미세한 주파수 차이 Δω=ω_H-ω_L을 사용하여 표시하면 식 (3)과 같다.

식 (3)은 삼각함수 합성 공식에 의해 식 (4)로 표시 가능하다.

식 (4)는 식 (5) ~ 식 (7)로 정리된다.

식 (6)은 완벽하게 선명한 코싸인 포락선 맥놀이고 식 (7)은 싸인 포락선 맥놀이다. 즉, Fig. 2(a)의 임의의 맥놀이는 (b)의 코싸인 맥놀이와 (c)의 싸인 맥놀이의 합이다. 이때 맥놀이 포락선의 정점에서의 진폭은 코싸인 맥놀이 진폭이고, 골에서의 진폭은 싸인 맥놀이 진폭이 된다.

Fig. 2 
Cosine and sine envelopes

포락선 대수감소를 구하기 위해 Fig. 3의 확대된 맥놀이 파형을 사용한다. Fig. 3의 포락선 간격은 T_b = 2π/Δω이고, t_N = t₀ + NT_b 이다. 두 시점 t₀, t_N 에서 포락선 정점의 진폭비를 취하면, 식 (8)과 같다.

Fig. 3 
Beat wave and logarithmic decrement

따라서 포락선 대수감소를 식 (9)와 같이 구할 수 있다.

식 (9)로부터 감쇠비를 식 (10)과 같이 구할 수 있다.

포락선 대수감소법으로 감쇠비를 추출하는 과정에서 발생 가능한 오차를 분석하였다. 성덕대왕신종의 1차 고유주파수 쌍과 감쇠비를 사용해서 다양한 패턴의 맥놀이 파형을 만들어 포락선 대수감소법을 적용하였다. 샘플링 주파수, 선명도, 주파수 쌍의 차이에 의한 오차를 검토했다. 목표주파수의 10배 이상으로 샘플링하면 오차율은 0.2 % 이내로 수렴하였다. 포락선의 선명도가 높아지면서 오차율은 영으로 수렴하였고, 선명도가 떨어져도 오차율은 0.25 % 수준이었다. 주파수 쌍의 차이가 증가하면 오차율도 증가하나, 성덕대왕신종의 주파수와 감쇠 수준에서는 임의의 맥놀이 파형을 사용하여 감쇠비를 추출하더라도 0.25 % 이내의 작은 오차율을 보였다.