최근 공학 분야에서는 현대 인간의 삶의 수준을 고려한 다양한 기술이 개발되고 있다. 그러나 기계 시스템과 장비의 사용에는 불가피한 요소가 작용하는데, 그 중 하나는 진동과 소음이다. 이 연구에서는 다양한 시스템의 진동과 소음을 억제하는데 효과적인 입자 충돌 댐퍼(particle impact damper; PID)를 고려하였다.

입자 충돌 댐퍼는 구조물이 진동할 때 진폭이 가장 크게 발생하는 부분에 입자를 담은 용기를 부착하여 구조물의 진동을 억제하는 댐퍼이다. 구조물이 운동할 때 입자가 운동 에너지를 흡수하여 자유롭게 이동하며 용기내 충돌이 발생된다. 이때 충돌로 인하여 마찰력이 발생하고 진동이 감쇠된다.

입자 충돌 댐퍼는 점탄성(visco elastic) 댐퍼와 달리 광범위한 온도에서 작동할 수 있기 때문에 사람의 손길이 닿지 못하는 가혹한 환경에서 사용할 수 있고, 시스템의 특성과 용도에 따라 설계를 할 수 있기 때문에 넓은 범위의 주파수에 적용이 가능하며, 반영구적인 사용이 가능하기 때문에 경제적이면서도 효과적인 댐퍼라고 할 수 있다.

입자 충돌 댐퍼는 1940년대 Lieber and Jensen⁽¹⁾이 항공기 플러터 현상 제어를 연구하면서 단일 입자를 담은 충돌 댐퍼로 처음 소개되었다. 이후 Masri^(2,3)는 주기당 양쪽으로 입자가 두 번 충돌할 때 충돌 댐퍼가 안정적으로 동작하는 것을 주장하였고 정현파에 의해 감쇠되는 1자유도 시스템의 엄밀해를 정의하고 그 안정성을 분석하였다. Masri and Ibrahim⁽⁴⁾은 Runge Kutta 방법을 사용하여 입자 충돌 댐퍼의 동작을 수치적으로 해석하고 임의의 가진에 의한 1자유도계 응답을 분석했다. Bapat and Sankar⁽⁵⁾는 자유 및 강제 진동 하에서 단일 입자 충격 댐퍼의 응답을 연구하였다. 시뮬레이션을 위하여 엄밀해를 이용한 수학적 모델을 정의하였고 실험을 진행하였다. 입자와 시스템의 질량비가 댐퍼의 감쇠력과 비례하는 것과 반발 계수 및 용기 내 간격이 일정 수치보다 커지면 감쇠 효과가 줄어드는 것을 입증했다. Papelou and Masri^(6~8)는 다중 입자를 담은 입자 충돌 댐퍼의 운동 방성식을 풀기 위해 다중의 입자를 동일한 질량을 가진 동등한 단일 입자로 대체하여 분석적 해를 계산하여 시스템의 수치 해석을 진행하였다. 마찬가지로 Friend and Kinra⁽⁹⁾는 입자를 덩어리 질량으로 표현하여 다중 입자 댐퍼의 운동에 대한 분석적 접근 방식을 개발하고 중력, 진동 진폭 및 입자와 용기 사이의 간격 거리의 영향을 제시하였다. Yang 등⁽¹⁰⁾은 단일 충돌 댐퍼를 활용하여 레일 진동 저감 연구를 진행하였고 시스템의 수치해석을 위하여 Runge-Kutta 방법을 사용하였다. Lu 등⁽¹¹⁾은 이산요소법(DEM)을 이용한 입자 충돌 댐퍼의 수치해석법을 제안하였다. 이산요소법은 각각의 입자의 모션 및 시스템과 충돌로 인한 상호작용을 모두 반영하여 거동을 수치적으로 분석하는 해석 방법이다. 그러나 이산요소법으로 시스템 해석을 구현하려면 많은 리소스와 시간이 소요된다.

충돌로 인한 비선형 운동으로 인해 입자 충돌 댐퍼의 운동을 정확하게 해석하기에는 어려움이 있다. 이산요소법을 이용한 해석법은 오랜 시간이 걸리기 때문에 다각적인 방면에서 시뮬레이션을 거쳐야하는 입자 충돌 댐퍼 설계에 응용하기에는 아직 적합하지 않다. 그렇기 때문에 다중입자를 단일 입자 모델로 변환하고 시스템의 거동을 시뮬레이션하는 방법으로 연구가 진행되어왔다.

이 연구에서는 여러 개의 입자가 일렬로 정렬된 입자 충돌 댐퍼를 장착한 시스템의 거동을 비교적 간단한 방법으로 예측할 수 있는 알고리즘을 개발하였다. 알고리즘을 입증하기 위하여 볼베어링 구슬을 담은 입자 충돌 댐퍼를 장착한 1자유도계 실험 모델을 제작하여 실험을 진행하고 시뮬레이션 결과와 실험 데이터를 비교하였다.

Fig. 3에 보이는 것과 같이 수평 운동을 하는 1자유도 알루미늄 구조물을 제작하고 윗부분에 볼베어링의 개수를 조절하여 담을 수 있는 플라스틱 용기를 부착하여 실험을 진행하였다. 입자를 제외하고 자유 진동 실험을 하여 시스템의 강성값과 감쇠 계수를 산출하였다. 이후 입자의 충돌실험을 진행하여 입자와 용기의 충돌 반발계수, e_s와 입자간의 반발계수, e_p를 측정하여 Table 1과 같은 값을 시뮬레이션 알고리즘에 반영하였다.

Fig. 3 
PID system with three particles

한 개부터 세 개까지 볼베어링 개수를 바꿔가며 시스템의 자유진동실험을 진행했다. 볼베어링의 개수가 많을수록 용기 내 입자 지름의 합을 제외한 빈공간(gap)의 거리는 짧아지며 그 길이는 Table 2와 같다.

실험을 진행하며 시스템과 입자의 초기 변위와 진동시 거동을 영상으로 기록하였고 영상내 물체의 변위를 측정할 수 있는 소프트웨어인 tracker⁽¹¹⁾를 이용하여 변위를 측정하였다. 입자가 움직일 때 미끄러짐 조건을 반영하지 않고 시스템과 입자의 변위를 시뮬레이션하여 실험 값과 비교하였다.

입자의 움직임을 제한하여 충돌 댐퍼가 작동하지 않을 때의 시스템의 변위를 측정하였고, 실험값과 시뮬레이션 값이 Fig. 4에서 나타나듯이 매우 유사한 것을 확인하였다.

Fig. 4 
Displacement when particles are immovable

Fig. 5에 나타난 것과 같이 볼베어링 한 개를 이용한 입자 충돌 댐퍼 시스템의 자유 진동과 볼베어링의 변위를 측정하고 시뮬레이션한 것과 비교하였다. 시뮬레이션으로 나타난 시스템의 진동 감쇠와 입자의 움직임이 실험 결과와 유사한 것을 확인하였다.

Fig. 5 
Displacements of the main system (a) and single particle (b)

Fig. 6에 나타난 것과 같이 볼베어링 두 개를 이용한 입자 충돌 댐퍼 시스템의 자유 진동 실험에서도 시스템의 감쇠와 두 볼베어링이 충돌하기까지 변위가 시뮬레이션 결과와 유사한 것을 확인하였다.

Fig. 6 
Displacements of the main system (a) and two particles (b)

Fig. 7에 나타난 것과 같이 볼베어링 세 개를 이용한 입자 충돌 댐퍼 시스템의 자유 진동 실험에서도 시스템과 볼베어링의 거동이 시뮬레이션한 것과 유사하게 진동이 억제되는 것을 알 수 있다.

Fig. 7 
Displacements of the main system (a) and three particles (b)

이 연구에서 소개한 알고리즘은 다수의 입자가 일렬로 배열된 유형의 입자 충돌 댐퍼 시스템을 예측하는 알고리즘이다. 입자의 회전을 반영한 입자 충돌 댐퍼 시스템의 운동 방정식을 유도하고 입자와 시스템의 변위를 이용하여 충돌 발생 시점을 정의하였다. 충돌을 반영한 동적 해석을 수행하기 위하여 이산 시간 상태 공간 방정식을 활용하였다. 세 가지 유형의 입자 충돌 댐퍼 시스템에 대한 자유 진동 실험을 수행하고 시스템과 입자의 변위 값을 측정하여 시뮬레이션 데이터와 비교하였다. Figs. 5 ~ 7에서 알 수 있듯이 실험과 시뮬레이션 결과에서 전반적인 유사성을 확인하였다. 그러나 시뮬레이션으로 구현한 시스템의 진동이 다소 빠른것을 Fig. 5, Fig. 6에서 볼 수 있다. 이는 용기 바닥의 거칠기에 의한 미끄러짐 특성 조건을 반영하지 않았기 때문인 것으로 사료된다. Fig. 7에서 비교적 높은 유사도가 나타난것은 입자를 추가할수록 입자와 용기 사이의 빈 공간이 짧아지며 미끄러짐에 의한 영향이 상대적으로 적기 때문인것으로 이해할 수 있다.

여러 입자의 충돌에 의한 시스템의 비선형적인 거동을 시뮬레이션에 모두 반영하기에는 어려움이 있는것이 사실이다. 그렇기 때문에 이전의 입자 충돌 댐퍼 연구에서 대부분 동적 해석을 수행할 때 여러 입자를 하나의 입자로 변환하여 엄밀해를 정의하여 분석하거나 Runge-Kutta 방법을 적용하였다. 이 연구에서 소개한 수치 해석법은 이산 시간 상태 공간 방정식을 활용하여 주 시스템과 각 입자의 운동과 충돌을 상대적으로 간단하고 정확하게 시뮬레이션하고 시스템의 특성을 반영하여 직관적으로 해석할 수 있는 장점이 있다. 앞으로 입자 충돌 댐퍼 설계를 위한 활용이 가능할 것으로 판단된다.