[논문 리뷰] Simple model of bouncing ball dynamics. Displacement of the limiter assumed as a cubic function of time

Andrzej Okniński, B. Radziszewski|arXiv (Cornell University)|2012. 03. 19.

Quantum chaos and dynamical systems참고 문헌 14인용 수 4

한 줄 요약

이 논문은 시간에 따라 주기적인 삼차함수로 변위하는 리미터와 충돌하는 뛰는 공의 해석적 접근이 가능한 모델을 제시한다. 이는 영향 시간과 역학을 정확히 계산할 수 있도록 하며, 주요 기여는 고정점의 주기 두배 이후에 발생하는 모서리 유형의 분기를 통해 혼돈대가 나타나는 것을 해석적으로 유도한 것이다.

ABSTRACT

Nonlinear dynamics of a bouncing ball moving vertically in a gravitational field and colliding with a moving limiter is considered and the Poincare map, describing evolution from an impact to the next impact, is described. Displacement of the limiter is assumed as periodic, cubic function of time. Due to simplicity of this function analytical computations are possible. Several dynamical modes, such as fixed points, 2 - cycles and chaotic bands are studied analytically and numerically. It is shown that chaotic bands are created from fixed points after first period doubling in a corner-type bifurcation. Equation for the time of the next impact is solved exactly for the case of two subsequent impacts occurring in the same period of limiter's motion making analysis of chattering possible.

연구 동기 및 목표

• 주기적으로 움직이는 리미터를 갖는 뛰는 공의 단순하면서도 해석적으로 풀 수 있는 모델을 개발하기 위해.
• 기본적인 사인 함수 리미터 모델에서 영향 시간 방정식을 풀기 어려운 문제를 해결하기 위해 삼차 변위 함수를 사용함으로써 해석적 타당성을 확보하기 위해.
• 비연속 동역학계에서 채터링, 그레이징, 혼돈과 같은 복잡한 역학적 행동의 해석적 조사를 가능하게 하기 위해.
• 한 주기 동안 다중 영향이 발생하는 정확한 조건을 유도하고 안정성 분석을 수행하기 위해.
• 특히 주기 두배 이후의 모서리 유형 분기로 인한 혼돈의 발생을 분기 분석을 통해 특성화하기 위해.

제안 방법

• 리미터의 변위를 주기적인 삼차함수 $ Y_c(T) = 12\tilde{T}(\tilde{T} - 1/2)(\tilde{T} - 1) $ 로 모델링하여 해석적 타당성을 확보한다. 여기서 $ \tilde{T} = T - \lfloor T \rfloor $ 이다.
• 비차원 중력 $ \gamma $ 와 복원계수 $ R $ 를 사용하여 영향 후 속도 $ V_i $ 와 영향 시간 $ T_i $ 를 연결하는 Poincaré 맵 (2.1) 을 구성한다. $ \Delta_{i+1} = T_{i+1} - T_i $ 로 정의한다.
• 일정한 가속도 하에서 운동 방정식을 풀어 다음 영향 시간 $ \Delta_{i+1} $ 의 정확한 해석적 표현을 도출한다. 이는 폐쇄형 비선형 맵 (3.12c,d) 를 이끌어내며 해석적 처리를 가능하게 한다.
• 역학을 단순화하고 그레이징 및 채터링 분석을 가능하게 하기 위해 상대 속도 변수 $ v_i = V_i - G(3T_i^2 - 3T_i + 1/2) $ 를 도입한다.
• 수렴 조건 $ \lambda_i < 1 $ 을 가정하여, $ v_i \to 0 $ 일 때 채터링의 최종 단계를 모델링하기 위해 점근적 전개 (3.14) 와 근사 맵 (3.16) 을 사용한다.
• 영향 시간이 무한히 반복되는 그레이징 조건을 유도하기 위해 $ T_{(\infty)} < 1 $ 을 풀며, 이로부터 $ T_B < T_i < T_C $ 의 범위를 도출한다 (식 3.21–3.22).

실험 결과

연구 질문

• RQ1 삼차 리미터 변위 함수의 사용이 뛰는 공 시스템에서 영향 시간을 정확히 계산하는 데 어떻게 기여하는가?
• RQ2 한 리미터 주기 내에서 다중 영향(채터링)이 발생하는 조건은 무엇이며, 이를 어떻게 해석적으로 특성화할 수 있는가?
• RQ3 혼돈대는 시스템에서 어떻게 발생하며, 그 형성 배경이 되는 분기 메커니즘은 무엇인가?
• RQ4 그레이징 만능의 역할은 무엇이며, 채터링의 안정성과 수렴성에 어떤 영향을 미치는가?
• RQ5 채터링 중에 N번째 영향 시간에 대해 정확하거나 근사적인 해석적 표현을 도출할 수 있는가?

주요 결과

• 혼돈대는 첫 번째 주기 두배 이후에 발생하는 모서리 유형의 분기로 인해 고정점에서 기인하며, $ R = 0.85 $ 와 $ \gamma \in [0, 0.06] $ 를 사용한 분기 다이어그램으로 확인되었다.
• 다음 영향 시간 $ \Delta_{i+1} $ 는 이차 방정식의 해를 정확히 구함으로써 계산되며, 이는 폐쇄형 비선형 맵 (3.12c,d) 를 이끌어내어 전면적인 해석적 처리를 가능하게 한다.
• 채터링은 근사식 (3.16) 을 통해 해석적으로 모델링되며, $ \Delta_{i+1} \approx \frac{2v_i}{-3G + 2 + 6GT_i} $ 를 도출한다. 수렴은 $ \lambda_i < 1 $ 일 때 보장된다.
• 그레이징 만능은 $ T_* \geq \max(T_{cr}, 0) $ 일 때만 존재하며, 여기서 $ T_{cr} = \frac{1}{2} - \frac{1}{3G} $ 이다. 이는 한 방향에서는 안정적이며 고유값 $ \Lambda_1 = 1 $ 을 가지며, 다른 방향에서는 중립적이다.
• 무한히 많은 영향(그레이징) 조건은 $ \max(T_{A_1}, T_{A_2}, T_B) < T_i < T_C $ 로 유도되며, $ T_B $ 와 $ T_C $ 는 각각 식 (3.21) 과 (3.22) 에 의해 주어지며, 이는 $ T_{(\infty)} < 1 $ 을 보장한다.
• 채터링 중 N번째 영향 시간은 정확히 $ T_{(N)} = T_i + \frac{2v_i}{-3G + 2 + 6GT_i} \cdot \frac{1 - R^N}{1 - R} $ 로 계산되며, 이는 영향 간격의 지수 감쇠를 보여준다.