• 대한전기학회
Mobile QR Code QR CODE : The Transactions of the Korean Institute of Electrical Engineers
  • COPE
  • kcse
  • 한국과학기술단체총연합회
  • 한국학술지인용색인
  • Scopus
  • crossref
  • orcid

  1. (Dept. of Electrical Engineering, Ulsan College, Korea)



Passivity, control system, synchronous generators, IEEE type 1 exciters, GT governors

1. 서론

일반적으로 passivity 이론은 비선형 시스템의 안정성을 분석하거나 이를 바탕으로 제어 시스템을 설계하는데 적용되고 있다. 기존에 널리 잘 알려진 Lyapunov 안정성 이론은 시스템의 상태 모델을 이용하지만, passivity 이론은 단지 시스템의 입출력 변수만을 적용할 수 있기 때문에 매우 복잡한 동적 시스템의 안정성 분석에는 보다 효과적인 것으로 알려져 있다(1).

최근에는 네트워크 기반 제어 시스템이나 가상 물리 시스템 그리고 원격 제어 로봇 시스템을 개발하는 경우에도 passivity 이론을 적용하여 제어 이론을 산출하거나 안정성을 분석하기도 한다(2,3,4). 이러한 제어 시스템들은 주로 원거리에 있는 대규모 제어 대상들을 다루기 때문에 시스템의 동특성 모델이 다소 복잡하며 원거리 제어로 인해 발생하는 여러 가지 통신 관련 문제들이 안정성에 영향을 미칠 수 있는데, 이러한 문제점들을 passivity 이론을 통해 해결하는 연구가 활발히 진행되고 있다.

예를 들어, 네트워크 기반 제어 시스템의 경우 통신 네트워크의 시간 지연, 데이터 패킷의 누락, 트래픽의 혼잡 등과 같은 특성들이 시스템의 전체 안정성을 결정하는데 크게 좌우되며, 이를 극복하기 위한 강인 제어 기법이 개발되었다(5). 또한 대규모 동적 시스템을 가상 물리 시스템으로 간주하고 이를 수학적 방정식으로 모델링을 실시한 후 무선 통신 네트워크를 통해 플랜트를 제어하기 위한 제어 토폴로지에도 기존의 passivity 이론을 적용하였으며, 이를 통해 전체 시스템의 안정성을 분석하였으며 관련 제어 규칙을 개발한 사례도 발표되었다(6). 보다 최근에는, 햅틱(haptic) 로봇 시스템 등과 같은 사람과 로봇 사이의 실시간 원격 제어를 실현하기 위하여 passivity 기반 제어 시스템 기술도 성공적으로 이어져 오고 있으며(7,8), 복잡한 비선형 시스템의 passivity 특성 분석을 위한 local passivity 이론이 발표되었다(9). 이러한 연구 동향을 살펴본 바와 같이, 대규모 동적 시스템과 이를 원거리에서 제어하기 위해 통신 네트워크 기술을 접목한 제어 시스템의 경우, 기존의 passivity 이론이 매우 중요하게 적용되고 있다는 것을 알 수 있다.

본 논문은 passivity 이론을 적용한 동기 발전기의 여자기와 조속기 제어 시스템의 설계 기법을 제안한다. 잘 알려진 바와 같이, 대규모 동기 발전기 시스템을 구성하는데 있어 여자기는 발전기의 출력 전압을 정격으로 제어하는 목적으로 사용되며 조속기는 회전자의 회전 속도 또는 동기 속도로 일정하게 유지하는데 사용된다(10). 따라서 이 두 제어 시스템의 바람직한 설계를 통해 발전기의 효율적인 출력 특성을 기대할 수 있다.

본 논문은 동기 발전기의 여자기 및 조속기 시스템 모델에 대하여 에너지 함수를 각각 정의한 후 그 미분식을 전개한 후 각각의 파라미터로 구성된 passivity 특성 조건을 수학적으로 도출한다. 제어기 파라미터의 선정은 도출한 passivity 조건을 주어진 제어 시간 동안 항상 만족할 수 있도록 구성된다. 본 논문에서 고려한 여자기 및 조속기 시스템은 산업현장에서 널리 사용되고 있는 IEEE Type 1과 Gas Turbine(GT) 모델을 각각 적용하였다. 제안한 제어 시스템의 설계 기법의 타당성 검증을 위하여 컴퓨터 시뮬레이션을 실시하였으며 발전기의 출력 파형을 분석하여 제어 특성을 검증하였다.

2. Passivity 이론

본 장은 passivity 이론에 대한 수학적 정의와 관련 이론을 정리하였다.

정의 1(11) 연속 시간 t에 대한 동적 시스템의 상태 벡터 x∈Rn와 출력 벡터 y∈Rp가 다음과 같은 수학적 모델로 표현되며

(1)
x ˙ ( t ) = f ( t , x , u ) y ( t ) = h ( t , x , u )

여기서 비선형 상태 함수 f와 출력 함수 h는 각각 f:Rn×Rp→Rn 및 h::Rn×Rp→Rp와 같으며 입출력의 개수는 동일하다고 가정한다. 식(1)에 대하여 0보다 크고 미분 가능한 연속 함수W(x) 가 존재하고 다음과 같은 조건을 만족하면 식(1)의 시스템은 passive라고 하고 그 특성을 passivity라고 일컫는다.

(2)
u T y W ˙ = W x f ( x , u )

또한 다음과 같은 다양한 부등식 조건에 대하여 여러 가지 passivity 정의를 나열할 수 있다.

(i) Input strictly passive : u T y W ˙ + u T φ ( u )

(ii) Output strictly passive : u T y W ˙ + y T φ ( y )

(iii) Strictly passive : u T y W ˙ + φ ( x )

여기서 u T φ ( u ) > 0 , y T φ ( y ) > 0 , φ ( x ) > 0 이다.

비고 1 식(1)의 동적 시스템이 passive하면 그 시스템은 안정하다고 알려져 있다. 일반적으로 동적 시스템의 안정성 분석에는 시스템의 상태 변수를 기반으로 하는 Lyapunov 이론을 많이 적용하지만, 시스템의 입출력 변수를 이용하여 passivity 특성을 분석한 후 안정성을 논하기도 한다.

정리 1(11) 두 시스템 φ1과 φ2그림. 1과 같이 피드백으로 구성되어 있을 때, 만약 시스템 φ1과 φ2가 모두 passive라고 하면, 전체 피드백 시스템 또한 passive하다. 여기서 φ1과 φ2의 입력 변수는 u1=r1-y2 및 u2=r2-y1이 된다.

그림. 1. 두 동적 시스템의 피드백 구성

Fig. 1. A feedback connection of two dynamic systems

../../Resources/kiee/KIEE.2018.67.4.561/fig1.png

비고 2 여러 개의 시스템으로 구성된 복잡한 동적 시스템의 경우 전체 시스템의 안정성이나 passivity 특성을 분석하기가 어려운 경우 각각의 단일 시스템에 대한 passivity를 우선 분석한 후 정리 1을 적용하여 전체 passivity 특성이나 이를 바탕으로 전체 시스템의 안정성을 판단할 수 있다.

3. 동기 발전기 모델링

본 논문에서 고려하는 동기 발전기의 선형화 모델은 다음과 같다(12).

(3)
E ˙ ( t ) = - 1 K 1 τ d E ( t ) - K 2 τ d δ ( t ) + 1 τ d v F ( t ) ω ˙ ( t ) = - K 3 τ j E ( t ) - K 4 τ j δ ( t ) + 1 τ j T M ( t ) δ ˙ ( t ) = ω ( t )

여기서 E, ω, δ는 고정자의 EMF 전압, 회전자의 속도 및 회전 각도를 각각 나타내며, τd > 0와 τj > 0는 시정수를 Kj > 0, I=1,···,4는 발전기의 전기적 사양으로부터 결정되는 파라미터이며, vF와 TM은 발전기의 입력 변수로서 여자기와 조속기로 전달되는 계자 전압과 토크를 각각 나타낸다. 식(3)으로부터 발전기의 출력 전압은 다음과 같이 주어진다.

(4)
V T ( t ) = σ 1 δ ( t ) + σ 2 E ( t )

여기서 σ1 > 과 σ2 > 는 마찬가지로 발전기의 사양을 통해 구해지는 상수가 된다. 그림. 2식(3)식(4)의 동기 발전기 시스템의 수학적 모델에 대한 블록선도를 보여준다. 여기서, 여자기와 조속기 시스템은 발전기의 기계적 및 전기적 출력을 바탕으로 주어진 기준 속도 ωr과 기준 전압 vr을 안정하게 추종하기 위하여 바람직한 제어 동작이 이루어져야 한다. 따라서 발전기 시스템의 만족스러운 제어 특성을 위하여 이 두 제어 시스템의 효율적인 제어 성능을 발휘할 수 있는 제어 기법을 도출하여야 한다.

그림. 2. 동기발전기 시스템의 블록선도

Fig. 2. A block diagram of synchronous generator systems

../../Resources/kiee/KIEE.2018.67.4.561/fig2.png

3.1 안정성 분석

식(3)의 안정성 분석은 선형 시스템 이론을 통해 간단하게 분석할 수 있다. 우선, 식(3)으로부터 상태 행렬은 다음과 같다.

(5)
A   =   - 1 K 1 τ d 0 - K 2 τ d - K 3 τ j 0 - K 4 τ j 0 1 0

또한 식(5)로부터 특성방정식은 다음과 같다.

(6)
s 3 + 1 K 1 τ d s 2 + K 4 τ j s - K 2 K 3 τ j τ d + K 4 K 1 τ j τ d = 0

기존에 잘 알려진 고유치(eigenvalue) 이론을 통해, 식(6)의 방정식의 근(root) 또는 고유치의 실수 부분의 값이 0보다 같거나 작으면 식(3)의 시스템 모델은 안정하다는 것을 알 수 있다. 하지만 식(6)은 고차 방정식이므로 기존에 널리 알려진 Routh- Hurwitz 안정성 판별법을 적용하여, 식(3)의 시스템이 안정하기 위한 조건을 다음과 같이 도출하였다.

(7-1)
● 조건 1 : K 2 K 3 τ j τ d > 0

(7-2)
● 조건 2 : K 4 > K 1 K 2 K 3

앞서 언급한 바와 같이 발전기 시스템의 모든 파라미터는 양수이므로 식(7-1)의 조건은 자연스럽게 만족하게 되며 식(7-2)는 발전기를 구축할 때 실제 파라미터 값을 적용하여 검토할 수 있다.

4. 여자기 제어 시스템

4.1 여자기 시스템 모델

본 논문에서 고려하는 동기 발전기의 여자기 시스템은 산업현장에서 널리 사용되고 있는 IEEE Type 1 모델로 선택하였다(13). 그림. 3은 이 여자기에 대한 블록선도를 보여주며 각각의 장치를 포함하고 있는 블록의 전달함수는 각각 다음과 같다(13).

그림. 3. IEEE Type 1 여자기 시스템 블록선도

Fig. 3. A block diagram of the IEEE Type 1 exciter

../../Resources/kiee/KIEE.2018.67.4.561/fig1.png

(8)
G R ( s ) = K A 1 + s T A G F ( s ) = 1 K E + s T E G B ( s ) = s K F 1 + s T F

여기서 KA, TA, KB, TB, KF, TF은 각 장치들의 파라미터로서 모두 양수를 가지며, 여자기의 바람직한 출력을 위하여 적절한 값으로 결정되어야 한다. 식(8)을 바탕으로 상태 변수 vB, vR, vF에 대한 방정식은 다음과 같이 표현된다.

(9)
v ˙ B ( t ) = - 1 T F v B ( t ) - K F K E T F T E v F ( t ) + K F T F T E v R ( t ) - K F T F T E f ( v F ) v ˙ R ( t ) = - 1 T A v R ( t ) - K A T A v T ( t ) - K A T A v B ( t ) + K A T A v r v ˙ F ( t ) = - K E T E v F ( t ) + 1 T E v R ( t ) - 1 T E f ( v F )

여기서 vr은 기준 전압, vR∈[VR,min,VR,max]의 범위를 갖는다. 또한 식(9)에서 f(vF)는 여자기의 포화(saturation) 특성을 나타내는 함수로서 다음과 같은 지수함수로 주어진다.

(10)
f ( v F ) = A E X e x p ( B E X v F )

여기서 AEX와 BEX은 포화곡선에 따른 상수로 주어진다. 일반적으로 여자기의 출력 전압은 양수이므로 식(10)에서 vF > 0 및 f(vF) > 0이 된다. 또한 식(9)에서 vR은 최소값 VR,min과 최대값 VR,max의 범위를 가지며 따라서 이 여자기 시스템은 비선형 특성을 포함하는 것을 알 수 있다.

4.2 Passivity 특성 분석

그림 3의 시스템 모델은 각각의 전달함수를 갖는 3개의 장치들로 구성되어 있으나, 블록선도 간소화 기법을 통해 그림 4와 같은 블록선도로 다시 표현할 수 있다. 또한 식(9)의 방정식을 다음과 같이 두 시스템 모델 과 로 분리하여 다시 나타낼 수 있다.

(11)
ϕ   :   ϕ 1 : v ˙ R ( t ) = - 1 T A v R ( t ) + K A T A v e ( t ) v ˙ F ( t ) = - K E T E v F ( t ) + 1 T E v R ( t ) - 1 T E f ( v F ) ϕ 2 : v ˙ B ( t ) = - 1 T F v B ( t ) + K F T F v ˙ F ( t )

여기서 ve = vr - vT - vB이다. 그림. 4의 시스템 모델은 두 시스템이 피드백 형태로 구성되어 있으므로 정리 1을 이용하여 passivity 특성을 분석한다. 다시 말해, 식(11)에서 시스템 φ1 및 φ2에 대하여 각각의 passivity 특성을 분석한 후에 그 결과를 바탕으로 전체 시스템 φ에 대한 passivity 특성을 알아볼 수 있다. 우선, 시스템 φ2에 대한 passivity 특성을 분석하기 위하여 에너지 함수를

그림. 4. IEEE Type 1 여자기의 간소화된 블록선도

Fig. 4. A simplified block diagram of the IEEE Type 1 exciter

../../Resources/kiee/KIEE.2018.67.4.561/fig4.png

(12)
W E , ϕ 2 = 1 2 v B 2

와 같이 정의하며 미분식은 다음과 같이 전개된다.

(13)
W ˙ E , ϕ 2 = - 1 T F v B 2 + K F T F v B v ˙ F

한편, 식(13)은 다음의 부등식을 만족한다.

(14)
W ˙ E , ϕ 2 v B v F - 1 T F v B 2 + K F T F v B v ˙ F + K F 2 T F ( v B - v ˙ F ) 2 - v B v F = v B v F - 1 T F v B 2 + K F 2 T F ( v B 2 + v ˙ F 2 ) - v B v F

다음으로 식(11)의 시스템 φ1에 대하여, 마찬가지로 에너지 함수를

(15)
W E , ϕ 1 = 1 2 ( v R 2 + v F 2 )

와 같이 정의하고 미분식을 구하면 다음과 같다.

(16)
W ˙ E , ϕ 1 = - 1 T A v R 2 - K E T E v F 2 - 1 T E v F f ( v F ) + v R ( K A T A v e + 1 T E v F )

또한 식(16)은 다음과 같은 부등식을 만족한다.

(17)
W ˙ E , ϕ 1 - 1 T A v R 2 - K E T E v F 2 - 1 T E v F f ( v F ) + v R 2 + 1 4 ( K A T A v e + 1 T E v F ) 2 = 1 4 K A T A T E v e v F - ( 1 T A - 1 ) v R 2 - ( K E T E - 1 4 T E 2 ) v F 2 - 1 T E v F f ( v F ) + 1 4 ( K A T A ) 2 v e 2

식(11)에서 시스템 φ1과 φ2가 각각 passive가 되기 위한 조건은 식(14)식(17)의 부등식을 만족하여야 한다. 제어기 설계 과정에서 이 두 조건을 모두 만족할 수 있는 제어기 파라미터 값을 적절하게 결정해서 그 목적을 이룰 수가 있다.

4.3 여자기 제어 파라미터 결정

식(14)의 부등식을 만족하는 파라미터를 산출하기 위하여 우선 식(14)의 우변의 세 번째 및 네 번째 항을 모두 0으로 두며 다음과 같다.

(18)
K F 2 T F ( v B 2 + v ˙ F 2 ) - v B v F = 0

여기서 KF와 TF에 대하여 정리하면

(19)
K F T F = 2 v B v F v B 2 + v ˙ F 2

와 같다. 식(19)에서 올바른 제어 동작에 대하여 제어 시간이 경과한 후에 안정한 출력 상태가 되면 vF가 평형점에 도달하여 미분값 v ˙ F는 0이 된다. 즉, lim t v F ˙ ( t ) = 0 이므로 식(19)의 최종값은

(20)
K F ( ) T F ( ) = 2 v F ( ) v B ( )

가 된다. 다음으로 식(17)로부터 시스템 φ1이 passive가 되기 위해서 우선 식(17)의 우변 항의 첫 번째 항에서 KA/(4TATE) = 1 또는 KA=4TATE가 되어야 하며 다음으로, 두 번째 항에서 0 < TA < 1가 되어야 한다. 그리고 세 번째 항에서는 K E / T E - 1 / 4 T E 2 0 이어야 하므로 TEKE > 1 되어야 한다. 마지막으로, 네 번째 항과 다섯 번째 항을 모두 0으로 두면

(21)
- 1 T E v F f ( v F ) + 1 4 ( K A T A ) 2 v e 2 = 0

와 같으며 여기서 다시 TE에 대하여 정리하면 다음과 같다.

(22)
T E = 4 T A 2 v F f ( v F ) K A 2 v e 2

여기서 l i m t - > v e ( t ) = - v B ( ) 이므로 TE의 최종값은 다음과 같다.

(23)
T E ( ) = 4 T A 2 v F ( ) f ( v F ( ) ) K A 2 v B 2 ( )

앞서 언급한 바와 같이 vF > 0 및 f(vF) > 0 므로 식(23)에서 TE > 0을 항상 만족한다.

5. 조속기 제어 시스템

5.1 조속기 시스템 모델

본 논문에서 고려하는 조속기는 그림. 5의 블록선도와 같이 표현되는 GT(gas turbine) 모델로 하였으며 각 장치의 전달함수는 다음과 같다(13).

그림. 5. GT 조속기 블록선도

Fig. 5. A block diagram of the GT governor

../../Resources/kiee/KIEE.2018.67.4.561/fig5.png

(24)
G 1 ( s ) = K w 1 + s T s r ,   G 2 ( s ) = 1 1 + s T c ,   G 3 ( s ) = 1 1 + s T t

이 전달함수를 바탕으로 상태 변수 x1 및 x2와 출력 변수 PM에 대한 동적 방정식은 다음과 같이 표현된다.

(25)
x ˙ 1 ( t ) = - 1 T s r x 1 ( t ) + K ω T s r e ω x ˙ 2 ( t ) = - 1 T c x 2 ( t ) + 1 T c x 1 ( t ) - 1 T c P p ˙ M ( t ) = - 1 T t p M ( t ) + 1 T t x 2

여기서 eωr-ω이며 ωr은 기준 속도를 가리키며, P는 상수로서 droop 값을 나타내며 기계적 출력 pM은 유한한 값을 가지며 pM∈[PM,min,PM,max]의 범위를 갖는다.

5.2 Passivity 특성 분석

식(25)의 조속기 모델에 대해서도 마찬가지로 passivity 특성을 파악하기 위하여 먼저, 상태 변수를 이용하여 0보다 큰 에너지 함수를 다음과 같이 정의한다.

(26)
W G = 1 2 ( x 1 2 + x 2 2 + p m 2 )

또한 식(26)을 시간 t에 대하여 미분한 후 식(25)에 적용하여 전개하면 다음과 같다.

(27)
W ˙ G = - 1 T s r x 1 2 + K w T s r e ω x 1 - 1 T c x 2 2 + 1 T c x 1 x 2 - 1 T c p x 2 - 1 T t p m 2 + 1 T t x 2 p m

식(27)에서 각각 다른 두 변수의 곱으로 주어진 항에 대하여 각각 다음과 같은 부등식을 구할 수 있다.

(28)
K w T s r e ω x 1 K w T s r e ω x 1 + ( e ω - K w 2 T s r x 1 ) 2 = e ω 2 + K w 2 4 T s r 2 x 1 2 ( 28 - 1 ) 1 T c x 1 x 2 1 T c x 1 x 2 + ( x 1 - 1 2 T c x 2 ) 2 = x 1 2 + 1 4 T c 2 x 2 2 ( 28 - 2 ) - 1 T c p x 2 - 1 T c p x 2 + ( p + 1 2 T c x 2 ) 2 = p 2 + 1 4 T c 2 x 2 2 ( 28 - 3 ) 1 T t x 2 p m 1 T t x 2 p m + ( p m - 1 2 T t x 2 ) 2 = p m 2 + 1 4 T t 2 x 2 2 ( 28 - 4 )

다음으로, 식(28)의 결과 식을 식(27)에 적용하면 다음과 같이 얻을 수 있다.

(29)
W ˙ G - 1 T s r x 1 2 + e ω 2 + K w 2 4 T s r 2 x 1 2 - 1 T c x 2 2 + x 1 2 + 1 4 T c 2 x 2 2 + 1 4 T c 2 x 2 2 - 1 T t p m 2 + p m 2 + 1 4 T t 2 x 2 2 = e ω 2 + ( 1 - 1 T t ) p m 2 + ( 1 - 1 T s r + K w 2 4 T s r 2 ) x 1 2 + ( - 1 T c + 1 2 T c 2 + 1 4 T t 2 ) x 2 2

또한 식(29)의 우변 항의 첫 번째와 두 번째 항은 다음과 같은 부등식으로 표현되므로

(30)
e ω 2 + ( 1 - 1 T t ) p m 2 < e ω 2 + ( 1 - 1 T t ) p m 2 + 1 2 ( e ω + p m ) 2 = e ω p m + 3 2 e ω 2 + ( 3 2 - 1 T t ) p m 2

식(30)식(29)에 대입하면 최종적으로 다음과 같다.

(31)
W ˙ G e ω p m + 3 2 e ω 2 + ( 3 2 - 1 T t ) p m 2 + ( 1 - 1 T s r + K w 2 4 T s r 2 ) x 1 2 + ( - 1 T c + 1 2 T c 2 + 1 4 T t 2 ) x 2 2

식(25)의 조속기 모델이 passive 특성을 갖기 위해서는 식(31)의 부등식을 항상 만족하여야 하며 이러한 목적은 시스템 파라미터를 적절하게 결정하여 이루어 질 수 있다.

5.3 조속기 제어 파라미터 결정

식(31)의 부등식을 만족시키기 위한 조속기 파라미터의 결정을 위해, 우선 식(31)의 우변 항에서 두 번째와 세 번째 항을 모두 0으로 둔다. 즉, 수학적으로 표현하면

(32)
3 2 e ω 2 + ( 3 2 - 1 T t ) p m 2 = 0

와 같으며 여기서 파라미터 Tt에 대하여 정리하면 다음과 같다.

(33)
T t = 2 p m 2 3 ( e ω 2 + p m 2 )

이 결과 식으로부터, 파라미터 Tt는 pm과 eω의 시간 함수가 되며 주어진 제어 시간 동안 이 두 변수의 값에 따라 그 값이 변화된다는 것을 알 수 있다. 한편, 조속기 제어 시스템의 올바른 제어 동작이 이루어질 경우, 시간이 지날수록 제어 오차는 0에 가깝기 때문에 식(33)에서 l i m t - > e ( t ) = 0 이 되며 이 때 Tt의 최종값은 Tt(∞)=2/3가 된다는 것을 알 수 있다. 다음으로, 식(31)의 식 4번째 항에서

(34)
1 - 1 T s r + K w 4 T s r 2 0

이어야 하므로 이 부등식의 해는 T s r 0 . 5 + - 0 . 5 1 - K w 가 된다. 여기서, 만약 Kω로 두면 Tsr0.5이 되며, 하지만 Tsr0이어야 하므로 Tsr∈(0,0.5]의 범위를 갖는다. 마지막으로, 식(31)의 마지막 항에 대하여

(35)
- 1 T c + 1 2 T c 2 + 1 4 T t 2 0

을 만족하여야 하며 이 부등식은 T c 2 T t 2 + - 2 T t T t 2 - 0 . 5 의 해를 가지며, 이때 Tc > 0이므로 T c ( 0 , 2 T t 2 ] 의 범위를 가지며 간단하게 T c = 2 T t 2 로 결정할 수 있다. 여기서 식(33)을 대입하여 Tc를 다시 표현하면

(36)
T c = 8 p m 4 9 ( e ω 2 + p m 2 ) 2

와 같으며, 또한 Tc의 최종값은 Tc(∞)=8/9이 된다.

6. 시뮬레이션

본 시뮬레이션에 적용된 동기 발전기는 참고문헌 [12](12)의 예제 6.7에 사용된 파라미터 값을 이용하였으며 식(3)식(4)의 발전기 파라미터 값은 표 1과 같다. 이 발전기 파라미터 값으로부터 식(7)에 제시된 발전기 모델의 안정성 조건을 모두 만족하는 것을 알 수 있다. 다음은 식(9)의 여자기 시스템 모델에서 AEX = 0.36, BEX = 0.2, VR,min = -15.5[pu], VR,max = 17.5[pu]이며 각 전달함수의 파라미터 값은 4.3절에 서술한 내용을 바탕으로 KE = 1, TF = 1, TA = 0.5로 두었으며 또한 K F = 2 v B v F / v B 2 + v ˙ F 2 , T E = 4 v F f ( v F ) / v e 2 , KA=4TE로 결정하였다. 다음으로, 식 (25)의 조속기 모델에 대하여 PM,min = 0, PM,max = 33[pu], P = 5[pu]이며 마찬가지로 5.3절의 파라미터 결정법에 의하여 Kw = 1, Tsr = 0.15, T c = 2 T t 2 이며 여기서 Tt는 식 (33)과 같이 두었다. 여자기와 조속기 제어 시스템은 발전기의 출력 전압과 회전 속도가 각각의 정격 전압과 동기 속도를 추정하여야 하며 이 두 정격에 대한 설정값은 각각 vr = 500[V] 및 ωr = 720[rpm]로 두었다. 그리고 이 두 제어 시스템에 대한 제어 성능의 조건은 다음과 같다.

표 1. 동기 발전기 시스템의 파라미터 값

Table 1. Parameter values for the synchronous generator model

Parameter

Values

K1

0.3072

K2

1.8051

K3

1.3174

K4

1.4479

σ2

0.0294

σ2

0.5257

τd

5.9

τj

0.00075

● 조건 1 : 5[%] 이하의 최대 언더슈트 및 오버슈트

● 조건 2 : 정상상태 오차값은 0

본 시뮬레이션은 2가지 시나리오에 대하여 각각 실시하였으며, 첫 번째 시뮬레이션에서 발전기 시스템의 초기 상태를 모두 0에서 기동하여 그 과도응답과 정상상태 응답 특성을 분석하였으며, 두 번째는 정격 운전 중에 있는 발전기에 파라미터 변화에 대한 동특성을 분석하였다. 파라미터의 변화는 표 1에서 Ki, i=1,2,3,4의 공칭값에서 각각 1%씩 증가하여 적용하였다.

그림 6그림 7은 첫 번째 시뮬레이션 시나리오에 대한 발전기의 출력 파형으로서 그림 6은 발전기의 단자 전압 그림 7은 회전자의 회전속도 파형을 각각 보여준다. 그림 6의 단자 전압의 경우 초기 시간부터 약 2.5초까지 오버슈트(overshoot)가 없는 과도응답 특성을 보이고 있으며 그 이후로 정상상태에서 설정값으로 안정하게 도달하는 것을 관측할 수 있다. 그림 7의 회전 속도 파형에서도 유사하게 오버슈트가 거의 없는 과도응답 특성을 보이고 약 2초 이후에 안정한 정상상태 응답에 도달하는 것을 볼 수 있다. 그리고 이 두 파형에서 알 수 있듯이, 발전기의 정상상태 응답에서 단자 전압의 설정값 500[V]와 동기 속도 720[rpm]에 매우 정확하게 도달하는 것을 볼 수 있으며, 위에 제시한 제어성능의 2가지 조건을 모두 만족하는 것을 알 수 있다. 따라서 첫 번째 시뮬레이션에서 여자기와 조속기 제어 시스템이 매우 효과적으로 동작한다는 것을 알 수 있다.

그림. 6. 발전기의 단자 전압(시뮬레이션 1)

Fig. 6. The terminal voltage of the generator(simulation 1)

../../Resources/kiee/KIEE.2018.67.4.561/fig6.png

그림. 7. 발전기의 회전자 속도(시뮬레이션 1)

Fig. 7. The rotor velocity of the generator(simulation 1)

../../Resources/kiee/KIEE.2018.67.4.561/fig7.png

그림. 8그림. 9는 두 번째 시뮬레이션에 대한 결과 파형으로서 10초에 파라미터 변화를 발생하였으며 이 때 발전기의 단자 전압과 회전 속도를 각각 보여준다. 우선 그림. 8의 단자 전압의 경우, 파라미터 변화가 발생한 이후부터 약 3초까지 오버슈트 및 언더슈트(undershoot)의 특성을 발생하는 과도응답을 보이고 있으며 그 이후로는 정상상태 오차가 거의 없는 정격 전압에 안정하게 도달하는 것을 볼 수 있다. 또한, 그림. 9의 회전 속도의 경우는 파라미터의 변화가 발생한 10초 이후에 과도응답 구간에서 매우 짧은 시간 동안 오버슈트와 언더슈트가 여러 번 발생하는 감쇠 진동 특성을 보이고 있으며 마찬가지로 약 12초 부근부터는 안정한 동기 속도에 도달하는 것을 관측할 수 있다. 한편, 제어 파라미터의 값이 최적의 값이 아니거나 최적의 시뮬레이션 환경이 아닌 경우 정상상태에서 바람직하지 못한 오차가 발생할 수가 있었다. 그림. 8그림. 9의 결과 파형으로부터 최대 퍼센트 언더슈트와 최대 퍼센트 오버슈트에 대한 관측값을 표 2에 기재하였다. 이 표의 관측값에 의하면 모두 5[%] 이하의 최대 퍼센트 언더슈트와 오버슈트를 가지며, 앞서 제시한 첫 번째 제어 조건을 모두 만족하는 것을 알 수 있다. 또한 정상상태 응답에서도 모두 0의 제어 오차를 보이고 있으며 이러한 두 결론을 바탕으로 두 번째 시뮬레이션 시나리오에서도 설계한 여자기 및 조속기 제어 시스템의 효과적인 동작 특성을 관측할 수 있었다.

그림. 8. 발전기의 단자 전압(시뮬레이션 2)

Fig. 8. The terminal voltage of the generator(simulation 2)

../../Resources/kiee/KIEE.2018.67.4.561/fig8.png

그림. 9. 발전기의 회전자 속도(시뮬레이션 2)

Fig. 9. The rotor velocity of the generator(simulation 2)

../../Resources/kiee/KIEE.2018.67.4.561/fig9.png

표 2. 최대 퍼센트 언더슈트 및 오버슈트(시뮬레이션 2)

Table 2. Maximum percent undershoot and overshoot (simulation 2)

Items

Maximum percent undershoot

Maximum percent overshoot

Values

Time

Values

Time

Terminal voltage

(Fig. 8)

4.26[%]

10.95

[sec]

1.42[%]

10.13

[sec]

Rotor velocity

(Fig. 9)

0.68[%]

10.19

[sec]

1.36[%]

10.08

[sec]

7. 결 론

본 논문은 passivity 이론을 적용하여 동기 발전기에 적용되는 여자기 및 조속기 제어 시스템 설계 방법론을 제안하였다. 우선, 주어진 여자기 및 조속기 시스템의 수학적 모델에 대하여 passivity 특성 분석을 실시한 후 주어진 제어시간 동안 항상 passivity 조건을 만족할 수 있는 제어기 파라미터를 산출하였다. 본 논문에서 제안한 제어기 설계법의 타당성 검증을 컴퓨터 시뮬레이션을 통해 이루어졌으며 발전기의 단자 전압과 회전 속도 파형 분석으로부터 제어 성능의 우수성과 신뢰성을 입증할 수 있었다. 향후 연구는 발전기 시스템의 실시간 테스트베드 구성과 실험을 통해 제안한 제어 시스템의 설계 기법에 대한 실적용 가능성 검증을 추가로 실시할 것이다.

감사의 글

이 논문은 2017년 울산과학대학교 교내학술연구비 지원에 의해 수행됨

References

1 
de Ruiter A. H. J., 2013, Some applications of passivity-based control and invariance principles, IET Control Theory Appl., Vol. 7, No. 7, pp. 1039-1048Google Search
2 
Hespanha J., Naghshtabrizi P., Xu Y., 년도, A survey of recent results in networked control systems, Proc. of the IEEE, Vol. 95, No. 1, pp. 138-162DOI
3 
Rajkumar R., Lee I., Sha L., Stankovic J., 2010, Cyber- physical systems: The next computing revolution, Proc. 47th Design Automation Conf., pp. 731-736DOI
4 
Sheridan T., 1989, Telerobotics, Automatica, Vol. 25, No. 4, pp. 487-507DOI
5 
Kottenstette N., Hall III J., Koutsoukos X., Sztipanovits J., Antsaklis P., 2013, Design of networked control systems using passivity, IEEE Trans. on Control Systems Technology, Vol. 21, No. 3, pp. 649-665DOI
6 
Antsaklis P., Goodwine B., Gupta V., McCourt M., Wang Y., Wu P., Xia M., Yu H., Zhu F., 2013, Control of cyberphysical systems using passivity and dissipativity based methods, European J. of Control, Vol. 19, pp. 379-388DOI
7 
Sun D., Naghdy F., Du H., 2016, Neural network-based passivity control of teleoperation system under time- varying delays, IEEE Trans. on Cybernetics, No. 99, pp. 1-15DOI
8 
Atashzar S., Shahbazi M., Tavakoli M., Patel R., 2017, A passivity-based approach for stable patient-robot interaction in haptics-enabled rehabilitation systems: modulated time-domain passivity control, IEEE Trans. on Control Systems Technology, Vol. 25, No. 3, pp. 991-1006DOI
9 
Xia M., Antsaklis P. J., Gupta V., Zhu F., 2017, Passivity and dissipativity analysis of a system and its approxi- mation, IEEE Trans. on Automatic Control, Vol. 62, No. 2, pp. 620-635DOI
10 
Cho H. C., 2017, A development of real-time monitoring techniques for synchronous electric generator systems, Trans. of the Korean Institute of Electrical Engineers, Vol. 66P, No. 4, pp. 182-187Google Search
11 
Khalil H. K., 2002, Nonlinear systems, Prentice Hall, Upper Saddle River, New JerseyGoogle Search
12 
Anderson P. M., Fouad A. A., 2003, Power system control and stability, Wiley Inter-Science, New JerseyGoogle Search
13 
Operation Technology Inc. , 2014, ETAP® 12.6 user guideGoogle Search

저자소개

조 현 철 (Hyun Cheol Cho)
../../Resources/kiee/KIEE.2018.67.4.561/au1.png

1997년 부경대 제어계측공학과(학사)

1999년 동아대학교 전기공학과(석사)

2006년 Uni- versity of Nevada-Reno (Ph.D.)

현재 울산과학대학교 전기전자공학부 교수

Tel : 052-279-3166

E-mail : hcjo@uc.ac.kr