Ⅰ서 론
최근 무인기(드론), 자율주행차, 무인선 등 무인 이동체 산업이 미래 신 성장 분야로 급부상하고 있 다. 그중에서도 드론은 개발과 보급 속도가 기대를 뛰어 넘으면서 ICT 뿐만 아니라 산업전반에 걸쳐 주목을 받으면서 빠르게 용도 확대가 진행되고 있 는데, 특히, 재해예방, 농업용, 산업현장관리, 스포 츠 활동과 방송 등의 영상제작에 효용성이 입증되 고 있다.
드론의 짐발 시스템은 외란과 장착부위 진동의 영향을 받음에도 불구하고 카메라 영상의 흔들림을 제거하고 원하는 표적 방향으로 영상장비를 유지하 는 장치로서 카메라 모듈과 안정화 구동장치로 나 눌 수 있다. 카메라 모듈은 직접적으로 영상을 획득 하는 부분으로 광학계 및 이를 제어하는 전자보드 로 구성되고, 안정화 구동장치는 rolling/ pitching/ yawing 운동을 하여 주변의 동적인 외란에도 불구 하고 카메라의 시선을 지향하거나 추적하는 목적으 로 사용된다.
짐발 시스템의 안정화 구동장치에 대한 제어는 PID(Propotional-Integral-Derivative Controller, 비례-적 분-미분 제어기)제어 기법이 가장 많이 사용되고 있 다. PID제어 기법은 조정을 위한 파라미터의 수가 적고, 그 구조가 단순하며, 각 요소의 작용에 대한 이해가 쉬우며, 비례이득과 적분시간 및 미분시간 등의 변화에 대한 응답특성의 변화가 직관적이고 적용방법이 용이하다는 특징을 가지고 있다[1]. 그 러나 제어대상에 시간지연이 존재하거나 혹은 시스 템 등의 특성이 변화하면 PID 파라미터의 재조정이 필요성이 대두되어, 이러한 문제의 해결책으로 비 선형제어의 간단화 등을 위해 Intelligent-PID제어라 는 새로운 제어기법이 M. Fliess, C. Join 등에 의해 제안되었다[2-4]. 이 기법은
-
제어대상의 모델을 거의 필요로 하지 않는다
-
구조가 PID제어와 동일할 정도로 간단하다.
-
고차의 시스템에 대해서도 파라미터의 조정이 간단하다.
-
제어대상의 특성이 변화하여도 제어기의 파라 미터를 재조정할 필요가 없다.
등의 특징을 가지는 것으로 기존 PID를 대체할 수 있는 유망한 제어기법의 하나이며 현재 다양한 종 류의 플랜트에 적용한 실험결과가 보고되고 있다 [5].
본 연구에서는 짐발 시스템의 3 축에 대한 모델 링 과정이 동일하므로, 중복성을 피하기 위하여 한 축에 대한 동적모델을 유도한다. 모델링된 플랜트 에 대하여 기존의 PID 및 I-PD제어기(Integral Propotioanl-Derivative Controller)를 설계하여 외란에 대한 각속도 응답특성을 MATLAB/Simulink를 통하 여 시뮬레이션 한다. 또한 본 연구에서 제안하는 Intelligent-PID제어기를 시스템 모델링을 기반으로 하지 않고 오차방정식을 안정화하는 조건으로 설계 한 후, PID, I-PD 및 I-PID 제어기법과 비교하여 응 답특성 및 파리미터의 변화에 대하여도 강인성을 확인한다.
Ⅱ짐발시스템의 수학적 모델
드론으로부터 깨끗하고 안정적인 영상 획득을 위해 카메라는 짐발시스템에 설치되는데, 짐발시스 템은 <Fig. 1>과 같이 3축(롤/피치/요) 운동을 할 수 있는 링크 구조물, 각 축의 회전운동을 위한 직접구 동 모터, 오일러 자세각 측정센서, 안정화 제어 알 고리즘 연산을 위한 DSP 제어기로 구성된다. 일반 적으로 3축 운동은 대칭적이므로 본 연구에서는 1 축에 대한 동적모델을 세우고 이에 대한 안정화 제 어기를 설계한다.
드론의 좌표계를 D(i, j, k)라 하고 짐발의 좌표계 를 G(r, p, y)라 하고, 각 좌표계의 축에 대한 각속도 벡터들을
라 하자.
<Fig. 2>와 같이 짐발 좌표계의 p축에 있는 구동 모터에 의한 회전운동 방정식을 유도해 본다. p축 이 각도 є만큼 회전할 때 짐발 축의 각속도를 드론 축의 각속도 함수로 나타내면[6-9],
이다.
한편, 짐발의 관성행렬을 다음과 같이 두자.
여기서, Jr, Jp, Jy는 각각 r, p, y축에 대한 관성모 멘트이고, Jrp, Jry, Jyp를 관성 곱모멘트 라고 하자. Newton의 법칙에 의하여 토크방정식은
이 되는데, 여기서 를 각 운동량이라 하고 다음과 같이 정의한다.
따라서 짐발의 각 축에 대한 토크는 식 (4)에서 각속도벡터와 운동량벡터의 벡터곱을 통하여
이 된다. 특히, 짐발의 p축에 대한 토크방정식은
와 같이 두는데, 여기서 τ 는 모터로부터 짐발 링크 에 인가되는 토크이며, 식 (5)의 항들을 대입하여 다음과 같이 정리한다.
만일, 짐발의 축에 대하여 무게분포가 대칭적이 라면 관성 곱모멘트가 0, 즉 Jrp = Jry = Jyp = 0 이므로, 식 (8)은 다음과 같이 된다.
여기에 식 (2)의 ωr , ωy를 대입하여, 드론의 기본 축 외란에 대한 짐발 p축의 회전운동 방정식을 다 음과 같이 얻는다.
한편 식 (10)에서 축을 회전시키는 토크 τ 가 기 어 없는 직접구동 DC모터에 의해 발생된다고 하자. <Fig. 3>은 직접구동 DC모터의 등가회로를 나타내 고 있으며 이 회로에 대한 동적방정식은 다음과 같 은 미분방정식으로 표현된다[10].
여기서 Jm = Ja+ JL 로서 모터의 관성과 링크의 관성을 합한 것이며, bm = B+ BL 은 모터의 마찰계 수와 짐발링크의 마찰계수를 합한 값인데, 본 연구 에서는 무시한다.
이상과 같이 유도된 짐발 링크에 대한 회전운동 방정식 (10)과 링크의 회전토크를 발생시키는 구동 방정식 (11)-(14)로 구성된 짐발시스템에 대한 블록 선도를 <Fig. 4>에 나타낸다.
ⅢIntelligent-PID 제어기 설계
시스템을 해석하고 설계하기 위하여 우선 수행 할 일은 제어하고자 하는 동적 시스템을 수학적으 로 표현하는 모델링 작업이다. 일반적으로 시간역 모델링 방법인 상태공간모델식과 에너지 접근법인 Lagrange 방정식을 이용한 모델식을 기반으로 제어 기를 설계하고 있지만, 본 연구에서는 비선형 요소 와 불확실한 요소가 많은 복잡한 시스템에 대하여 모델링 없이 제어기를 설계하는 방법에 대하여 기 술한다.
일반적으로 시스템 동적방정식은 입력 u(t)와 시 스템 출력 y(t)의 함수로
와 같이 표현되는데, 이 식은 아주 짧은 시간동안 유효한
로 나타낼 수 있다. 여기서 n은 일반적으로 1 또는 2가 되며, α는 설계자가 결정하는 값이고, F(t)는 y(t), u(t) 및 이들의 고계미분의 함수로서 제어대 상의 모든 정보가 들어있는 미지의 항으로 임의 시 간에 입력 u(t)와 측정값 y(n) (t)을 이용하여 실시간 으로 측정할 수 있다.
만일, F(t)를 추정할 수 있다면 폐루프 제어시스 템에서 제어대상의 바람직한 출력을 위하여 제어입 력은 다음과 같이 결정한다. 특히, n=2일 경우 식 (17)로부터[2-4],
여기서
y* (t)는 제어대상의 출력이 추종해야 하는 목표 신호
e(t)= y* (t)- y(t) 는 추종오차
KP , KI, KD는 일반적인 PID 제어기 이득 이다. 한편 F(t)는 식 (17)로부터 y(n) (t) - αu(t) 이 므로 식(18)에 대입하여 폐루프 응답의 오차방정식
을 얻으며, 양변에 라플라스변환을 적용하면 다음 과 같은 폐루프 특성방정식을 얻는다.
그러므로 식(20)이 안정으로 되도록 KP , KI, KD 를 설정하면 e(t) →0 (r →∞)로 된다.
실제로, 식 (18)의 F(t)는 미지이기 때문에 식 (17)를 이용하여 실시간으로 계산해야 할 필요가 있 다. 이를 위해 식 (17)의 u(t) 대신에 충분히 작은 값을 가지는 시간 h(>0) 만큼 이전의 값 u(t-h)를 사용하여 다음과 같이 F(t)를 추정한다.
이 값을 식 (18)에 대입하여 폐루프 제어시스템 에서 제어대상의 바람직한 출력을 위하여 실제의 제어입력을 다음과 같이 결정한다.
이를 Intelligent PID 제어기라고 하며, 이 제어기 에 의해 구현된 제어시스템 구조를 <Fig. 5>에 나타 낸다.
Ⅳ시뮬레이션 및 결과
2장에서 유도한 짐발 동역학 플랜트에 대한 안정 화 제어기로 많이 사용되는 고전적인 PID 제어기들 과 본 연구에서 제안하는 Intelligent PID제어기와의 성능을 비교하여 본다. 특히, 플랜트의 환경 변화에 따른 안정화 성능을 비교하여 본다. 제어기들의 성 능을 시뮬레이션하기 위해서 DC모터의 파라미터를 <Table 1>과 같이 둔다.
또한 카메라가 부착된 짐발시스템이 무게가 1[Kg]이고 반경이 14[cm]인 회전체라 가정하면 관 성모멘트는 라 둘 수 있다. 따라서 모터의 관성모멘트와 합하여 시스템 의 관성모멘트는 이 된다. 본 시물레이션에서 마찰계수는 무시한다.
먼저, 고전적인 PID제어기는 오차신호(e(t))나 플 랜트의 출력신호(y(t))에 비례(KP), 적분(KI), 미분 (KD) 상수를 곱하여 다음과 같은 조합의 알고리즘 으로 제어신호 u(t)를 형성할 수 있다.
여기서 식(23)으로 표현된 제어알고리즘은 <Fig. 6>의 (a)로, 식(24)로 표현된 제어알고리즘은 <Fig. 6>의 (b)와 같은 제어기 구조를 형성한다. 여기서 고전적인 PID 제어기의 게인값들을 Zigler-Nichols 기법에 의해 KP = 0.96 , KI = 10.7 , KD = 0.0001 라 둔다. 목표 입력 각속도(y*(t) = ωc)를 0.35[rad/sec] 로 할 때 드론 본체의 i축의 각속도가 ωi =0.09[rad/sec] 로 변하고, k축의 각속도가 ωk =0.17[rad/sec] 로 변동하며, j축의 각속도가 ωj = 0~0.61[rad/sec]로 변하는 변동 외란에 대해 짐발 p축의 각속도 출력(y(t) = ωp)를 시뮬레이션한 결과를 <Fig. 7>에 나타내었다.
다음은 본 논문에서 제시한 Intelligent-PID제어기 의 설계에 따라 오차방정식 (20)의 안정화 조건으로 부터 KP = 12 , KI = 8 , KD = 6 을 얻고, 여기서 시 간지연 h를 0.01로 두고 α=1로 하여 제어기를 설계 한다. 이를 <Fig. 5>와 같은 구조의 제어시스템으로 구성한 후, 앞서 제시한 실험 조건을 대입하여 시뮬 레이션한 결과 <Fig. 8>과 같은 응답특성을 얻었다. 이 결과로부터 PID제어기나 I-PD제어기의 경우 변 동외란에 대한 드론의 각속도 응답특성이 악화됨을 볼 수 있다. 그러나 Intelligent PID제어기의 경우 환 경변화에 매우 강인한 응답특성을 보여주고 있다. 또한 Intelligent PID 제어기의 설계인 경우, 동적 모 델을 사용하지 않고도 제어기 파라미터를 선택할 수 있었고 제어성능도 우수함을 확인할 수 있다.
Ⅴ결 론
최근 드론은 재해예방, 농업용, 산업현장관리, 스 포츠 활동, 방송 등의 영상제작에 효용성이 입증되 고 있다. 드론의 짐발 시스템은 외란과 장착부위 진 동의 영향을 받음에도 불구하고 카메라 영상의 흔 들림을 제거하고 원하는 표적 방향으로 영상장비를 유지하는 장치이다. 짐발 시스템의 안정화 구동장 치에 대한 제어는 PID제어 기법이 가장 많이 사용 되고 있는데, 제어대상에 시간지연이 존재하거나 혹은 외란 등의 특성이 변화하면 파라미터의 재조 정이 필요해 진다. 이러한 문제점을 극복하기 위하 여 본 연구에서는 Intelligent-PID제어 기법을 짐발 시스템의 안정화 구동제어에 적용하여 보았다. 기 존 PID제어 기법과의 응답특성을 비교하기 위하여 드론의 3축 운동 중에서 한 축에 대한 동적 모델을 유도하였고, 이 모델을 기반으로 PID 및 I-PD제어 의 각속도 응답특성을 시뮬레이션 및 비교 검토 하 였다. 이 결과, 본 논문에서 제안한 Intelligent-PID 제어에 의한 응답특성은 외부환경의 변화에도 불구 하고 강인한 Servo 특성을 보이는 우수성을 확인하 였다. 앞으로 연구결과를 제작된 드론의 짐발시스 템 안정화 제어기로 구현하여 안정되고 깨끗한 영 상이미지를 얻는 실험을 진행할 예정이다.