비선형 시스템 안정성 분석을 위한 기술함수의 정성적 그래픽 표현 방법

초록

본 논문은 조각식 비선형 특성을 갖는 시스템의 기술함수(Describing Function)를 복잡한 수식 없이 손쉽게 그래픽으로 그릴 수 있는 정성적 방법을 제시한다. 기존의 정확한 기술함수 계산이 요구하는 적분과 복잡한 식을 대신해, 비선형 특성의 구간별 기울기와 불연속점만을 이용해 두 가지 기본 형태(데드존·릴레이)의 합으로 표현하고, 이를 기반으로 두 개의 규칙과 알고리즘을 제안한다. 제안된 방법은 교육 현장에서 빠른 한눈에 파악 가능한 한계 진동 예측을 가능하게 하며, 사례 연구를 통해 기존 정확한 방법과 동일한 정성적 결과를 얻음이 검증된다.

상세 분석

논문은 기술함수 방법이 비선형 시스템의 한계 진동을 예측하는 데 유용하지만, 정확한 기술함수 (F(X))를 구하기 위해서는 비선형 특성 (y(x))에 대한 적분식(식 3‑4)을 풀어야 하며, 이는 특히 조각식·불연속성을 가진 경우 계산량이 크게 증가한다는 점을 지적한다. 이를 해결하기 위해 저자는 두 가지 기본 비선형 모델, 즉 ‘데드존(dead zone)’과 ‘임계값을 가진 릴레이(relay with threshold)’의 기술함수를 각각 식 5‑8 로 정리하고, 이 두 함수를 선형 결합함으로써 임의의 조각식 비선형을 재구성할 수 있음을 보인다(속성 1). 핵심 아이디어는 각 구간의 기울기 차이 ((m_{i+1}-m_i))와 불연속 높이 (Y_j)를 해당 기본 함수의 파라미터로 매핑하는 것이다.

Algorithm 1은 입력으로 구간 경계점 ((x_i,y_i))와 진폭 벡터 (X)를 받아, 각 구간에 대해 (\Phi)와 (\Psi) 함수를 적용해 (F(X))를 수치적으로 계산한다. 이 과정은 MATLAB 코드로 공개되어 재현성을 높인다. 그러나 저자는 교육적 목적에서 ‘정성적 그래픽’이 더 유용하다고 주장한다. 속성 2와 Algorithm 2는 (F(X))의 형태를 ‘연속적이면서 구간별 지수‑유사 상승((\tilde\Phi))’과 ‘임펄시브‑유사 하강((\tilde\Psi))’ 함수로 근사한다. 구체적으로, 초기 구간에서는 첫 번째 기울기 (m_0)가 그대로 유지되고, 각 구간 전이점 (X_j)에서는 (\tilde\Phi)가 0에서 시작해 점차 (m_j)에 수렴한다. 동시에 불연속점이 존재하면 (\tilde\Psi)가 추가되어 급격한 상승을 모사한다. 이러한 규칙은 손으로도 쉽게 그릴 수 있도록 설계되었으며, 실제 사례(그림 4·5)에서 정량적 (F(X))와 정성적 (\tilde F(X))가 거의 일치함을 보인다.

핵심 기여는 다음과 같다. 첫째, 복잡한 적분 대신 기울기와 불연속만으로 기술함수를 구성하는 체계적 방법을 제시했다. 둘째, 정성적 그래픽 규칙을 수학적으로 근거(속성 2)와 알고리즘(Algorithm 2)으로 명시해 교육용 툴킷을 제공했다. 셋째, MATLAB 구현과 사례 분석을 통해 기존 정확한 방법과 정성적 방법이 한계 진동 존재 여부와 대략적인 진폭·주파수를 동일하게 예측한다는 점을 실증했다. 한계점으로는 비대칭·동적 비선형(예: 히스테리시스)이나 고차 고조파 영향을 무시한다는 가정이 남아 있어, 정밀 설계 단계에서는 여전히 정확한 기술함수 계산이 필요하다는 점을 인정한다.