[논문 리뷰] A numerical scheme for impact problems

Lætitia Paoli, Michelle Schatzman|ArXiv.org|1999. 10. 04.

Dynamics and Control of Mechanical SystemsEngineering참고 문헌 14인용 수 20

한 줄 요약

이 논문은 n 자유도를 가진 기계적 충격 문제를 해결하기 위한 새로운 수치적 스킴을 제안한다. 여기서 운동은 C³-정규성 조건을 만족하는 닫힌 집합 K 내에 제약된다. 이 스킴은 충격 시간을 명시적으로 탐지하지 않고, 복원계수 e ∈ [0,1]을 사용한 국소 유한차분 접근 방식을 통해 이를 해결하며, 사전 추정치를 바탕으로 수렴성을 증명하고, 국소적 및 전역적 해의 존재성을 확립한다.

ABSTRACT

We consider a mechanical system with impact and n degrees of freedom, written in generalized coordinates. The system is not necessarily Lagrangian. The representative point of the system must remain inside a set of constraints K; the boundary of K is three times differentiable. At impact, the tangential component of the impulsion is conserved, while its normal coordinate is reflected and multiplied by a given coefficient of restitution e between 0 and 1. The orthognality is taken with respect to the natural metric in the space of impulsions. We define a numerical scheme which enables us to approximate the solutions of the Cauchy problem: this is an ad hoc scheme which does not require a systematic search for the times of impact. We prove the convergence of this numerical scheme to a solution, which yields also an existence result. Without any a priori estimates, the convergence and the existence are local; with some a priori estimates, the convergence and the existence are proved on intervals depending exclusively on these estimates. This scheme has been implemented with a trivial and a non trivial mass matrix.

연구 동기 및 목표

• 충격 시간을 명시적으로 탐지하지 않는 기계적 시스템의 충격 문제를 위한 수치적 방법을 개발한다.
• 일반적인 소산력과 C³-정규 집합 K로 정의된 제약 조건을 가진 라그랑주가 아닌 시스템을 다룰 수 있다.
• 연속 시간 충격 문제의 해에 수렴하는 수치적 스킴의 수렴성을 확립한다.
• 사전 추정치와 스킴 수렴성을 이용해 국소적 및 전역적 해의 존재성을 증명한다.
• 경계 직선화에 의해 발생하는 2차 오차를 포함한 유한차분 스킴에서의 경계 국소화 문제에 철저히 대응한다.

제안 방법

• 스킴은 일반화된 좌표와 질량 행렬 M(u)을 사용하여 역학을 모델링하며, 경계 ∂K에서의 충격 반응 μ는 λdϕ(u)로 표현된다.
• 시간에 대한 유한차분 이산화를 사용하며, 경계 ∂K 근처에서 국소 재구성 기법을 적용하여 비연속적인 충격 역학을 처리한다.
• 사영된 탄성 반동 계수 e를 활용한 코탄젠트 메트릭 기반 수직성 조건을 통해 접선 운동량의 보존과 법선 운동량의 반사 조건을 강제한다.
• 핵심적 혁신은 ∂K 근처에서 동역학을 근사하면서 경계 직선화에 의해 발생하는 2차 오차를 제어하는 국소 스킴을 사용하는 것이다.
• 수렴성 증명은 사전 추정치와 컴actness 추론에 기반하며, 이로써 이산 해 수열이 연속 해로 수렴함을 보여준다.
• 이산 속도 노름 |V^m|_{U^m}^2 에 대한 추정치를 이용하여, 이산 속도 노름이 연속 속도 노름의 본질적 상한으로 수렴함을 증명한다.

실험 결과

연구 질문

• RQ1 충격 시간을 명시적으로 탐지하지 않아도 되는 충격 문제를 위한 수치적 스킴을 구성할 수 있는가?
• RQ2 경계가 곡선일 경우, 경계 ∂K에서 비연속적인 동역학을 유한차분으로 어떻게 근사할 수 있는가?
• RQ3 수치적 스킴이 진정한 충격 문제의 해로 수렴하기 위한 조건은 무엇인가?
• RQ4 운동량과 위치에 대한 사전 추정치가 해의 전역 존재성에 어떻게 영향을 주는가?
• RQ5 복원계수 e는 수치적 스킴의 안정성과 수렴성에 어떤 역할을 하는가?

주요 결과

• 수치적 스킴은 연속 충격 문제의 해로 수렴함을 증명하며, 허용 가능한 집합 내의 초기 조건에 대해 해의 존재성을 입증한다.
• 사전 추정치 없이 국소적으로 수렴성과 존재성을 확립하였고, 이러한 추정치가 존재할 경우 전역적으로 성립한다.
• 스킴이 이산 속도 노름 |V^m|_{U^m}^2 가 연속 속도 노름의 본질적 상한으로 수렴함을 증명함으로써 수렴성이 입증된다.
• 초기 운동량과 K의 기하학적 구조에 의해 결정되는 시간 τ(R) 이후로 스킴이 수렴하지 않으면 모순이 유도된다.
• 증명은 ∂K 근처에서 경계 직선화에 의해 발생하는 2차 오차를 추정치를 통해 철저히 제어하는 국소화 기법에 기반한다.
• 운동량과 위치에 대한 사전 경계를 조합함으로써 국소 수렴성과 함께 전역 존재 결과를 확보하였으며, 이는 스킴이 오직 이러한 경계에 의해 결정되는 시간 τ(R)까지 연장될 수 있음을 보장한다.