Advanced
Stabilization of Underwater Glider by Buoyancy and Moment Control: Feedback Linearization Approach
Stabilization of Underwater Glider by Buoyancy and Moment Control: Feedback Linearization Approach
Journal of Ocean Engineering and Technology. 2014. Dec, 28(6): 546-551
Copyright © 2014, Korean Society of Ocean Engineers
  • Received : July 02, 2014
  • Accepted : December 15, 2014
  • Published : December 31, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
성철 지
호재 이
문환 김
지현 문

Abstract
This paper addresses a feedback linearization control problem for the nonlinear dynamics of an underwater glider system. We consider the buoyancy and moment as control inputs, which come from the mass variation and elevator control, respectively. Moment-to-force coupling increases the nonlinearities, which make the controller design difficult. By using a feedback linearization technique, we convert the nonlinear underwater glider to an equivalent linear model and design a linear controller. The controller for the equivalent converted linear system is designed using sufficient conditions in terms of linear matrix inequalities. Then, the control input of the nonlinear model of an underwater glider is formulated from the linear control input. An experimental examination is implemented to verify the effectiveness of the proposed technique.
Keywords
1. 서 론
바다는 무한한 광물, 에너지, 공간 등의 자원을 보유하고 있지만 인간은 그것의 극히 일부분만을 활용하고 있으며 대부분은 오랜 세월동안 사용하지 못하고 있다. 해양은 인간이 접근하기에는 위험한 극한의 요소가 많기 때문에 직접적인 해양 탐사를 위해서는 무인 수중운동체의 운용이 필수적이다.
대표적인 해양탐사 시스템으로 무인잠수정(Autonomous underwater vehicle)이 잘 알려져 있다( Nakamura and Savant, 1992 ). 무인 잠수정은 운용자의 개입이 최소화되어 스스로 자율운항이 가능하여 군사적, 학술적 연구에 사용되는 제어시스템이다. 무인잠수정은 해양환경의 다양한 정보를 수집하거나 감시, 정찰활동 등의 중요한 임무를 수행하며 특히 사람이 접근하기 어려운 심해 탐사에 활용이 가능하다는 장점이 있다. 하지만 프로펠러를 사용하여 추진력을 생성하므로 선체의 무게가 무겁고 부피가 크다. 또한 추진동력을 생성하기 위한 에너지의 저장량이 한정되어 있으므로 장시간 탐사에 불리하다는 단점이 있다.
한편 수중글라이더는 바다의 심층과 표층을 오가면서 원하는 지점으로 이동할 수 있도록 고안된 무인 해양탐사 로봇으로서 역할은 무인잠수정과 크게 다르지 않지만 프로펠러 추진방식을 사용하는 무인잠수정과는 달리 별도의 추진체를 사용하지 않기때문에 에너지 측면에 있어서 무인잠수정보다 효율적이다( Park, 2013 ). 수중글라이더의 추진력은 유체의 유입과 배출에 의한 부력 조절과 내부의 질량 이동과 날개의 움직임을 이용한 모멘트 변화에 의해 발생된다. 따라서 동력 에너지의 사용이 비교적 적기 때문에 장시간, 장거리의 해양 탐사에 적합하다( Park, 2013 ; Park et al., 2012 ).
최근 수중글라이더의 동역학 해석 및 제어기 설계에 관한 연구가 크게 주목받기 시작했다. 그 중에서 특이섭동(Singular perturbation) 기법은 수중글라이더 동역학 해석에 효과적인 방법 중 하나로 알려져있다. 특이섭동 기법은 하나의 시스템을 응답이 느린 동역학과 빠른 동역학으로 분리하여 해석하는 방법으로 장주기와 단주기 모델로 나누어지는 수중글라이더의 동역학에 적합하다. 하지만 장주기의 축소모델(Reduced model)과 단주기의 경계층모델(Boundary-layer model)은 근사화된 비선형 모델로써 이를 기반으로 설계된 제어기는 수중글라이더를 국소적으로(Locally) 안정화할 수 있지만 전역적으로(Globally) 안정화하지는 못한다( Bhatta and Leonard, 2004 ; Bhatta and Leonard, 2008 ; Zhang et al., 2012 ).
논문 ( Bhatta and Leonard, 2002 )에서는 내부 질량이동에 의한 수중글라이더의 안정화 문제를 다루며 제어기 설계를 위해 피드백 선형화(Feedback linearization) 기법을 사용한다. 여기서, 선형화된 시스템은 상대차수(Relative degree)가 시스템의 차수보다 낮은 내부동역학으로 표현되며 전체 동역학은 등가의 삼각시스템(Triangular system)으로 표현된다. 이것은 다시 비선형의 영점동역학(Zero dynamics)과 부분적으로 선형화된 시스템(Partially linearized system)으로 분리되어 해석된다. 이 방법을 이용하여 분리된 각각의 시스템을 안정화하는 제어기를 설계할 수는 있지만 근사화 오차가 존재하기 때문에 전역 안정도를 보장하지 못한다.
본 논문은 비선형 수중글라이더의 안정화 문제를 다룬다. 수중글라이더의 제어입력으로서 유체의 유입과 배출에 의해 조절되는 부력과 승강타(Elevator)의 움직임에 의해 발생하는 모멘트를 고려한다. 상태변수들의 복잡한 비선형성과 모멘트 제어입력에 의한 추가적인 상호연결(Coupling) 힘은 제어기 설계를 어렵게 만든다. 수중글라이더의 제어기를 설계하기 위해 피드백 선형화 기법을 사용하여 등가의 선형시스템 제어기 설계문제로 변환한다. 차의 시스템 동역학을 부분적으로 선형화하여 차의 선형시스템으로 변환한 후 모멘트 제어기를 설계한다. 이어 모멘트 제어기를 포함한 나머지 차 시스템에 대한 추가적인 피드백 선형화를 수행하여 부력 제어기를 설계한다. 변환된 각 부분 선형시스템의 제어기 설계조건은 선형행렬부등식의 형태로 제시되며 설계된 선형 제어기로부터 비선형 수중글라이더의 제어입력을 만든다. 이 방법은 기존의 연구와는 달리 근사모델을 사용하지 않기 때문에 폐루프 시스템의 전역적 안정성을 보장한다.
2. 수중글라이더 모델링
Fig. 1 은 짧은 구간에서 이동하는 수중글라이더의 모형을 도식화하여 나타낸 것이다. 여기서, 아래첨자 ‘1’과 ‘2’는 이동하는 수중글라이더의 순간 위치를 구분하기 위해 쓰인다. 수중글라이더의 모델링을 위해 다음과 같은 가정을 도입한다.
PPT Slide
Lager Image
Hydrodynamnic force and moment of underwater glider
가정 1: 수중글라이더는 타원형의 강체이며 차원의 수직 평면에서 운항함을 가정한다.
가정 2: 수중글라이더는 그림 1 과 같이 중심각이 ∅ 1 ―∅ 2 인호(Arc)의 궤적을 따라 움직임을 가정한다
수중글라이더의 동역학 모델링을 위해 동체기준의 좌표계를 사용한다. 벡터들 ei ∈ℝ 3 , i j 3 ≔{1, 2, 3}는 동체의 기하학적 중심(Geometric center)으로부터 각각 동체의 머리 방향, 지면(paper)의 윗 방향, 동체의 수직 윗 방향의 단위벡터이며 오른손 법칙을 따른다.
M2 ∈ℝ ≥0 e 2 방향의 모멘트, L , D ∈ℝ >0 는 양력(Lift)과 항력(Drag)을 나타내며 Vk ∈ℝ, k j 2 는 경로방향의 속도, ∅ k , αk , θk , k j 2 는 각각 경로각(Path angle), 받음각 (Attack angle), 피치각(Pitch angle)을 나타낸다.
e 2 방향의 모멘트 제어입력 u 2 와 결합계수(Coupling factor) δ ∈ ℝ에 의한 e 3 방향의 상호작용 힘 f = δV 2 u 2 을 고려하자.
가정 3 ( Bhatta and Leonard, 2008 ; Zhang et al., 2012 ): 매우 짧은 구간 내에서 ∅ 1 ≈∅ 2 ≈∅, V 1 V 2 V 이며 주변 유체의 관성이 수중글라이더 동역학에 미치는 영향이 거의 없다고 가정한다.
가정 3에 의해 동체 기준축(Body axis) e 2 , e 3 의 부가질량(Added mass)과 e 2 방향의 부가모멘트는 고려하지 않는다. 뉴턴 제 2 법칙과 운동량방정식을 사용하면 수중글라이더의 동역학은 다음과 같이 나타낼 수 있다[자세한 유도는 ( Bhatta and Leonard, 2008 )를 참조].
PPT Slide
Lager Image
여기서, m 은 수중글라이더 동체의 질량, m 0 은 수중글라이더의 질량에서 수중글라이더 부피와 동일한 물의 질량을 뺀 값이며 부력조절기에 의해 조절 가능한 값이다( Bhatta and Leonard, 2008 ).
가정 4: 부력 조절에의한 수중글라이더의 내부 질량 이동과 모멘트는 발생하지 않음을 가정한다.
다음의 변수들을 정의하자.
PPT Slide
Lager Image
여기서, K D0 , K D 는 항력계수(Drag coefficient), K L0 , K L 는 양력계 수(Lift coefficient), M M0 , K M 은 피치모멘트(Pitch moment)계수, K q 는 피칭감쇠(Pitching damping)계수이다. 2 J 2 e 2 방향 의 각속도와 관성모멘트를 의미한다.
다음의 새로운 상태변수들과 제어입력을 고려하자.
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
이며 아래첨자 ′ e ′는 각각의 평형점을 나타내기 위해 사용된다. 새로운 상태변수들과 부력 제어입력 u 1 을 추가로 도입하면 (1)은 다음의 0-평형점을 갖는 5차 방정식 형태로 표현된다.
PPT Slide
Lager Image
3. 주요 결과: 피드백 선형화
정의 1: 벡터공간 f ( x ), g ( x )를 따라 h ( x )의 방향성미분(Lie derivative)을 각각 다음과 같이 정의한다.
PPT Slide
Lager Image
또한 고차의 방향성미분을 다음과 같이 정의한다.
PPT Slide
Lager Image
PPT Slide
Lager Image
상태변수
PPT Slide
Lager Image
와 제어입력
PPT Slide
Lager Image
를 정의하면 수중글라이더 동역학 (2)는 다음과 같은 형태로 표현된다.
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
이며
PPT Slide
Lager Image
이다. 다음의 표준행렬
PPT Slide
Lager Image
을 가지는 선형 변환 TL : ℝ 5 → ℝ 5 을 고려하자. 식 (3)의 상태변수 x ∈ℝ 5 에 대하여 변환 T L ( x )을 계산하면 변환된 상태변수의 동역학을 다음과 같이 기술할 수 있다.
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
PPT Slide
Lager Image
이며 비선형 수중글라이더 모델 (3)과 등가이다. (4)에서 ζ-동역학의 입력 v 2 의 비선형 요소를 상쇄시키고 상태궤환 폐루프 시스템을 구성하기 위해 다음 형태의 모멘트 제어입력을 사용하면 ζ-동역학은 다음과 같은 폐루프 시스템의 형태를 가지게 된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서, F 2 는 선형시스템의 제어이득행렬이다. A 2 + B 2 F 2 가 Hurwitz 행렬이 되도록 F 2 를 설계하면 ζ-동역학은 점근적으로 안정하나 시스템 (3) (동등하게 (2))의 안정성을 보장하지는 못한다. 만약 상태변수
PPT Slide
Lager Image
를 점근적으로 안정하게 만드는 제어입력 u 1 이 존재한다면 (4) (동등하게 (3) 혹은 (2)) 또한 점근적으로 안정하다.
이제
PPT Slide
Lager Image
-동역학을 점근적으로 안정화하는 u 1 을 설계하자. 이를 위해
PPT Slide
Lager Image
-동역학에 대한 추가적인 피드백 선형화를 수행한다. 다음의 비선형 변환 TN : ℝ 3 → ℝ 3 을 고려하자.
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
이다. 부력 제어입력
PPT Slide
Lager Image
을 정의하면 변환된 상태변수 η 에 대하여 다음의 폐루프 선형 시스템을 얻을 수 있다.
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
이며 식 (4)의 ζ-동역학과 등가이다. 이제 폐루프 시스템 (8)과 (6)으로 부터 확장된 시스템을 다음과 같이 쓸 수 있다.
PPT Slide
Lager Image
여기서,
PPT Slide
Lager Image
이며 (3)과 등가이다. 따라서 선형시스템 (9)가 점근적으로 안정하다면 시스템 (3)을 점근적으로 안정화하는 제어입력을 계산할 수 있다. 비선형 시스템 (3)의 제어입력 u 는 다음의 형태를 가진다.
PPT Slide
Lager Image
선형시스템 (9)의 제어이득행렬
PPT Slide
Lager Image
는 다음 정리로부터 구할 수 있다.
정리1 : 다음 선형행렬부등식
PPT Slide
Lager Image
을 만족하는 양한정 행렬
PPT Slide
Lager Image
과 비정방행렬
PPT Slide
Lager Image
이 존재한다면 선형시스템 (9)는 점근적으로 안정하다. 여기서, M = FX 이다.
증명: 양한정 행렬
PPT Slide
Lager Image
에 대하여 양한정 함수 V ( z )≔ z T P z 를 정의한다. 리아푸노프 안정도 이론에 의해 다음 부등식이 만족되면 (9)는 점근적으로 안정하다.
PPT Slide
Lager Image
여기서, P -1 에 의한 합동변환(Congruence transform)과 P -1 = X , FX = M 의 변수치환이 사용되었다.
참고1 : 피드백 선형화를 위한 상태변수들의 선형 변환 T L 과 비선형 변환 T N 은 미분동형사상(Diffeomorphism)이어야 한다. 즉, 새로운 상태변수의 동역학은 변환 이전의 비선형 동역학의 위상적 성질을 보존해야 하며 변환 함수는 연속함수로서 미분 가능하고 그것들의 역함수도 역시 미분가능해야한다.
참고 2: 정리 의 선형행렬부등식 (11)의 해공간은 컨벡스 집합(Convex set)의 형태를 가진다. 식 (11)의 최적 해를 구하기 위해 본 논문에서는 MATLAB robust control tool box의 LMI solver를 사용하였다.
4. 모의실험
다음 모델 파라미터를 갖는 (2) 형태의 시스템을 고려하자 ( Bhatta and Leonard, 2008 ).
PPT Slide
Lager Image
PPT Slide
Lager Image
정리 1을 만족하는 다음의 행렬들
PPT Slide
Lager Image
PPT Slide
Lager Image
로부터
PPT Slide
Lager Image
,
PPT Slide
Lager Image
를 계산하면 다음을 얻는다.
F 1 =[−2.1212 −5.2727 −2.7879]
F 2 =[−1.3125 −0.9375]
이제 제어이득행렬 F 1 , F 2 로부터 (10) 형태의 제어 입력 u 를 계산할 수 있다. 시스템 (2) 또는 (3)에 대하여 초기 값 V =0.70 m / s , ∅=0 rad, m 0 =1.46 kg, α =0.0175 rad, 2 =0 rad/s을 선정하여 t ∈[0.20]에서 모의실험을 수행한다.
Fig. 2 는 수중글라이더의 속도 V , 경로각 ∅, 질량 m 0 , 받음각 α , 피치 각속도 2 의 시간 응답을 보인다. 그림에서 알 수 있듯이 각 변수들의 궤적(파란 실선)은 초 이내에 평형점(빨간괘선)에 수렴하여 시스템이 안정함을 보이고 있다.
PPT Slide
Lager Image
Time responses of glider variables V[m/5], ∅[rad], m0[kg], α[rad], [rad/s] (solid blue) and their equilibrium points (dashed red)
Fig. 3 은 제어입력의 시간 응답이다. 상태변수들이 평형점에 도달한 후의 부력제어입력 u 1 과 모멘트 제어입력
PPT Slide
Lager Image
는 0에 수렴함을 보여준다
PPT Slide
Lager Image
Time responses of buoyancy control input u1 and moment control input
Fig. 4 는 수중글라이더의 이동궤적을 보여준다. 여기서, 빨간괘선은 평형점 상태에서 수중글라이더의 기준경로, 파란 괘선은 실제 이동경로이다. 받음각을 나타내기 위해 수중글라이더의 모형을 타원 형태로 나타내며 수중글라이더의 방향을 표시하기 위해 꼬리 부분을 십자(Cross)로 표시한다. 그림에서 알수 있듯이 제안한 기법을 적용한 수중글라이더는 기준 경로에 빠르게 수렴하여 안정한 움직임을 보이고 있다.
PPT Slide
Lager Image
Glider trajectory (dashed blue) and desired path (dashed red)
5. 결 론
본 논문은 부력과 모멘트 제어에 의한 수중글라이더의 안정화 문제를 논하였다. 안정한 시스템을 설계하기 위해 5차의 시스템을 부분 선형화하여 2차의 선형시스템으로 모델링 한 후 모멘트 제어기를 설계하였다. 이어 모멘트 제어기를 포함한 나머지 2차 시스템에 대한 선형화를 수행하여 부력 제어입력을 만들어냈다. 각 선형시스템의 제어 이득행렬은 선형행렬부등식의 설계조건으로부터 구하였으며 이로부터 수중글라이더의 제어입력을 만들어냈다. 모의실험을 통해 제안한 피드백 선형화 기법의 효용성을 입증하였다.
Acknowledgements
이 논문은 2014년도 정부(미래창조과학부)의 재원으로 한국연구 재단의 지원을 받아 수행된 연구임(No. 2014R1A2A2A01005664).
References
Bhatta P. , Leonard N.E. 2002 Stabilization and Coordination of Underwater Gliders Decision and Control, 2002, Proceedings of the 41st IEEE Conference on 2081 - 2086
Bhatta P. , Leonard N.E. 2004 A Lyapunov Function for Vehicles with Lift and Drag: Stability of Gliding Decision and Control, 2004. CDC. 43rd IEEE Conference on 4101 - 4106
Bhatta P. , Leonard N.E. 2008 Nonlinear Gliding Stability and Control for Vehicles with Hydrodynamic Forcing Automatica 44 (5) 1240 - 1250    DOI : 10.1016/j.automatica.2007.10.006
Zhang F. , Tan X. , Khalil H.K. 2012 Passivity-based Controller Design for Stabilzation of Underwater Gliders American Control Conference (ACC) 5408 - 5413
Nakamura Y. , Savant S. 1992 Nonlinear Tracking Control of Autonomous Underwater Vehicles Robotics and Automation, 1992. Proceedings., 1992 IEEE International Conference on A4 - A9
Park J. 2013 Underwater Glider: Its Applicability in the East/Japan Sea Ocean and Polar Research 35 (2) 107 - 121    DOI : 10.4217/OPR.2013.35.2.107
Park Y. , Lee S. , Lee Y. , Jung S. , Jang N. , Lee H. 2012 Report of East Sea Crossing by Underwater Glider Journal of the Korean Society of Oceanography 17 (2) 130 - 137