Mobile QR Code QR CODE : Korean Journal of Air-Conditioning and Refrigeration Engineering
Korean Journal of Air-Conditioning and Refrigeration Engineering

Korean Journal of Air-Conditioning and Refrigeration Engineering

ISO Journal TitleKorean J. Air-Cond. Refrig. Eng.
  • Open Access, Monthly
Open Access Monthly
  • ISSN : 1229-6422 (Print)
  • ISSN : 2465-7611 (Online)

  1. ( Student, Graduate School of Refrigeration and Air-conditioning Engineering, Pukyong National University, Busan, 48513, Korea )
  2. ( Professor, Department of Refrigeration and Air-conditioning Engineering, Pukyong National University, Busan, 48513, Korea )



Variable speed refrigeration system(VSRS; 가변속 냉동시스템), Robust control(강인 제어), Sliding mode control(슬라이딩 모드 제어), Kalman filter(칼만 필터), Feedforward control(피드포워드 제어)

기호설명

$A$, $B$, $C$: 계수행렬 [-]
$d$: 부동작 시간 [sec]
$e$: 제어오차 [℃]
$F_{f}$: 피드포워드 제어기 게인 [-]
$F_{i}$: 인버터 주파수 [Hz]
$K$: 제어기 게인 [-]
$k$: DC gain [-]
$L$: 칼만 필터 게인 [-]
$p$: 위상평면의 상태변수 [-]
$T$: 제어량(온도 현재값) [℃]
$T^{*}$: $T$의 설정값 [℃]
$u_{c}$: 연속 제어입력 [V]
$u_{d}$: 불연속 제어입력 [V]
$u_{f}$: 피드포워드 제어입력 [V]
$V_{o}$: EEV 개도 지령 [step]
$\Delta$: 변화량 [-]
$\lambda$: 슬라이딩 라인의 기울기 [-]
$\sigma$: 슬라이딩 함수 [-]
$\tau$: 시정수 [sec]
$\Phi$: 경계층 두께 [-]

하첨자

$cp$: Compressor [-]
$ev$: EEV [-]
$o$: Oil outlet [-]
$s$: Superheat [-]

상첨자

$'$: 행렬의 전치 [-]
$\hat{}$: 추정값 [-]

1. 서 론

가변속 냉동시스템(VSRS)은 에너지 절약 성능이 뛰어나며 고정밀 온도 제어가 가능해 산업계 전반에 널리 사용되고 있다. VSRS는 압축기 회전수를 변화시켜 냉매의 질량유량을 제어함으로써 목표 온도를 제어하며, 동시에 COP를 최대로 유지하기 위해 과열도도 제어하므로 다입․출력계(MIMO) 시스템으로 취급된다.

VSRS의 온도 제어 방법은 수학적 모델 필요성 유무에 따라 모델 기반 제어와 AI 수법 제어로 구분된다. 모델 기반 제어는 전달함수와 상태 공간 모델 기반으로 나뉘며, 전자의 대표 제어기로는 PID, 후자는 최적제어 (optimal control)를 들 수 있다. PID는 단일 입․출력 시스템(SISO)일 경우, 설계가 쉽고 구현이 간편하면서도 비교적 양호한 제어성능을 가진다. 하지만 공칭모델 기반이므로 모델 불확실성 등에 대한 제어의 강인성 (robustness)을 확보하기 어렵다. LQR(Linear Quadratic Regulator), LQG(Linear Quadratic Gaussian)로 대표되는 최적제어는 MIMO 시스템의 제어기 설계에 적합하며 제어 정도(control accuracy)와 입력 에너지의 절충을 통해 최적성은 담보하지만 이 또한 모델 불확실성과 외란에 대한 제어의 강인성을 충분히 보장하지는 않는다. 한편, AI 수법의 Fuzzy 제어는 제어계 설계 시 제어대상의 동특성 모델이 불필요하므로 비선형성이 커 선형 모델을 얻기 어려운 VSRS의 제어법으로 적합하다. 그러나 이 제어법은 고정밀 제어일수록 제어 로직도 방대해져 고가의 마이크로프로세서를 필요로 한다. 특히, 핵심 구성 요소인 멤버십함수(membership function)와 규칙베이스 (rule base)의 설계가 전문가의 지식과 경험에 지나치게 의존함으로써 체계적인 설계가 어렵다.

모델 기반 제어를 위한 VSRS의 동특성 모델로는 열교환기에 대해 Navier-Stokes의 지배방정식을 적용하여 도출한 해석학적 모델이 제안되었다. 이 모델은 고차의 비선형 편미분방정식으로 선형화 및 저차원화 과정에서 모델의 불확실성(uncertainty)이 발생하며, 파라미터의 동정(identification) 과정에서 불확실성이 증대된다. 또한, 제어기 설계 시 상태관측기의 타당성 검증도 쉽지 않다. 또 다른 모델링 방법으로는 특정 동작점에서 동특성 실험을 통해 전달함수를 먼저 구하고, 이로부터 저차원의 실용적인 상태 공간 모델을 얻는 방법이 있다. 그러나 이 모델도 파라미터 동정 과정 및 모델링 시와 다른 운전 환경 조건에서는 모델의 불확실성이 커지므로 제어의 강인성을 보장할 수 없게 된다. 결국, PID나 최적제어와 같은 모델 기반의 제어는 모델의 불확실성과 같은 공칭모델이 갖는 한계성 및 잡음을 포함한 외란 유입으로 인해 설계한 제어기의 성능이 크게 열화(deterioration) 하므로 제어의 강인성(robustness)을 확보하기 어려운 문제점이 있다. 따라서 모델의 비선형성과 불확실성이 큰 VSRS의 정밀한 온도 제어를 위해서는 제어기의 게인이 공칭모델의 파라미터에 의존하지 않도록 설계하는 것이 바람직하며, 이를 만족하는 모델 기반형의 강인한(robust) 제어법의 적용이 요구된다. 이를 위해 본 논문에서는 칼만 필터와 피드포워드 제어기를 갖는 슬라이딩 모드 제어(Sliding Mode Control; SMC)를 제안한다. SMC란 제어 시스템 구조를 바꾸는 가변구조제어(Variable Structure Control; VSC)의 한 방법이며, 설계한 슬라이딩 평면 또는 선(sliding surface or line) 상에 슬라이딩 모드를 실현하는 제어를 말한다. 이 제어는 불연속 제어입력의 게인이 공칭모델의 파라미터에 독립적으로 설계되므로, 모델 불확실성과 외란에 강인한 제어성능은 물론 속응성과 정밀성도 동시에 보장한다.

VSRS를 대상으로 한 기존의 SMC 설계법으로는 MBM(Moving Boundary Model)과 열교환기 지배방정식에 의한 해석학적 모델로부터 Hankel의 특이치 분해를 통해 얻은 5차원 상태 공간 모델에 기반한 사례가 있다.이 모델 또한 여전히 고차원이어서 파라미터의 동정과 상태관측기의 타당성을 검증하기가 어려웠다. 이를 해결하고자 동특성 실험으로 구한 전달함수 모델로부터 Pade 근사를 거쳐 실용적인 2차원 상태 공간 모델을 얻고, 최적 절환 초평면을 갖는 SMC가 최근 제안되었다. 최적 절환 초평면을 갖는 SMC는 모델링과 제어기 설계가 쉬우면서도 우수한 제어성능을 보였다. 하지만, 이 설계법은 최종 서보계가 3차원의 위상평면(phase plane)을 갖게 됨으로써 상태변수의 거동을 직관적으로 파악하기 어려웠다. 또한, 기존의 방법은 최적 절환 초평면의 설계 파라미터인 하중함수를 대각행렬로 표현하므로 제어입력과 설계 파라미터 간의 인과관계를 파악하기 어려울 뿐만 아니라 하중함수 설정 시 번거로운 반복 시행을 요하는 문제가 있다.

본 논문에서는 VSRS에 대한 기존 SMC의 문제점들을 해결하기 위해 저차원의 실용적인 상태 공간 모델을 토대로 2차원 위상평면을 갖는 슬라이딩 평면을 설계한다. 이를 통해 제안한 제어기는 상태변수의 거동을 직관적으로 이해하기 쉬우며, 제어성능 개선이 용이한 장점을 가짐과 동시에 하중함수의 번거로운 반복 시행 문제를 줄일 수 있는 장점이 있다. 뿐만 아니라 SMC의 주요 설계 파라미터 및 제어기 게인이 스칼라(scalar)로 표현되므로 제어입력과 설계 파라미터들 간의 인과관계를 쉽게 파악할 수 있다. 특히, 엄밀한 상태변수 추정과 제어기 설계의 일반성이 담보되도록 칼만필터로 상태관측기를 구성함으로써 프로세스 잡음 및 관측 잡음에 강인한 설계가 되도록 하였다. 불연속 제어입력으로 인한 채터링(chattering) 발생을 억제하기 위해 부호함수를 포화함수로 대체하고, 이로 인해 목표 제어량의 과도특성이 열화하는 부작용을 줄이기 위해 피드포워드 (feedforward) 제어기를 추가로 설계한다. 포화함수의 사용은 sigmoid 함수에 비해 마이크로프로세서의 처리 부담을 줄여 실용성 향상에 기여한다. 제안된 제어기는 VSRS 기반의 오일쿨러 시스템을 제어대상으로 한 시뮬레이션과 실험을 통해 그 타당성이 검증되며, PI 제어기와의 제어성능을 비교함으로써 그 유효성을 확인한다.

2. 칼만 상태관측기를 갖는 슬라이딩 모드 제어기 설계

2.1 VSRS의 2차원 상태 공간 모델링

Fig. 1은 VSRS로 구성된 오일쿨러 시스템(Oil-Cooler System; OCS)의 모식도이다. OCS는 정밀 공작기계의 공작물 가공 과정에서 발생하는 열을 신속히 제거, 공작물의 열 변형과 가공 정밀도 저하를 방지한다. 이 장치는 공작물 가공 부위에 공급되는 오일 온도를 일정하게 제어하는 역할을 한다. OCS의 오일출구온도는 압축기의 회전수를 가변시켜 냉매의 질량유량을 조절함으로써 제어된다. 이때 냉매 질량유량의 급격한 과다 또는 과소로 인한 액백현상이나 과열 증기 압축 등의 부작용을 최소화 하고, COP가 최대인 운전점에서 시스템을 작동시키기 위해 과열도(superheat)도 보조적으로 제어한다. 과열도는 증발기 입․출구의 압력차가 작을 경우, 증발기 출구와 입구 측 온도 차로 구해지며 EEV의 개도를 조절함으로써 제어된다.

Fig. 2는 VSRS의 압축기와 EEV의 전달함수 모델을 구하기 위한 입․출력 관계를 나타낸 블록선도이다.

Fig. 1 Schematic diagram of the oil cooler system based on VSRS.
../../Resources/sarek/KJACR.2021.33.10.506/fig1.png

Fig. 2 Input and output variables of transfer functions for a compressor and EEV.
../../Resources/sarek/KJACR.2021.33.10.506/fig2.png

OCS의 주된 제어량은 오일출구온도 $T_{o}$이며, 보조 제어량은 과열도 $T_{s}$이다. 두 제어량에 대한 제어입력 (조작량)은 인버터 주파수 $F_{i}$와 EEV 개도(opening angle) 지령 $V_{o}$이다. Fig. 2에서 실선으로 표시된 전달함수는 압축기와 EEV의 동특성 모델로서 제어기를 설계하기 위한 용도이고, 점선으로 표시된 전달함수는 시뮬레이션 시에 간섭계인 제어대상의 실제 거동을 엄밀히 모사하기 위한 용도이다. 전달함수 $G_{d}(s)$는 외란(disturbance)인 열부하 변동에 따른 $T_{o}$의 거동, $G_{icp}(s)$와 $G_{iev}(s)$는 압축기 회전수 및 EEV 개도 변화가 $T_{o}$와 $T_{s}$에 미치는 간섭 영향을 각각 나타낸다.

식(1)Fig. 1의 OCS를 대상으로 Fig. 2의 전달함수 $P(s)$를 얻기 위해 동작점 근방에서 압축기와 EEV의 조작량인 $F_{i}$와 $V_{o}$를 미소 변동시켜 얻은 전달함수이다. 동특성 실험 시, 간섭 영향을 주는 조작량은 일정(constant; c) 값으로 고정하였다. 식(1)로부터 전달함수는 부동작 시간을 갖는 1차계임을 알 수 있다.

(1)
$\begin{cases} \left . P_{cp}(s)=\dfrac{\Delta T_{o}}{\Delta F_{i}}\right |_{V_{o}=c}=\dfrac{-0.43}{1680s+1}e^{-51s}\\ P_{ev}(s)=\dfrac{\Delta T_{s}}{\Delta V_{o}}\left |_{F_{i}=c}=\dfrac{-0.045}{67s+1}e^{-5s}\right . \end{cases}$

식(2)는 열부하 변동에 대한 $T_{o}$의 거동, 식(3)은 압축기의 회전수($F_{i}$) 변동이 $T_{s}$에 미치는 간섭 영향을 각각 나타낸다. Fig. 2의 $G_{iev}(s)$는 그 영향이 미미하여 시뮬레이션 시에는 제외되었다. 이는 본 논문이 비간섭제어 보다 외란을 포함한 모델의 불확실성에 강인한 제어기 설계에 초점을 맞추고 있기 때문이다.

(2)
$G_{d}(s)=\dfrac{\Delta T_{o}}{\Delta T_{d}}=\dfrac{19.9}{1790s+1}$

(3)
$G_{icp}(s)=\dfrac{\Delta T_{s}}{\Delta F_{i}}=\dfrac{348.1s-0.467}{885s+1}$

SMC의 설계를 위한 OCS의 2차원 상태 공간 모델은 식(1)을 Pade 근사를 통해 선형 2차계로 표현한 후, 역 라플라스 변환을 통해 제어계의 해석과 설계가 쉬운 가제어 표준형인 식(4), 식(5)와 같이 유도된다.

(4)
$\begin{cases} \begin{bmatrix}\dot x_{1cp}(t)\\\dot x_{2cp}(t)\end{bmatrix}=\begin{bmatrix}0& 1\\-2.33\times 10^{-5}& -0.0398\end{bmatrix}\begin{bmatrix}x_{1cp}(t)\\x_{2cp}(t)\end{bmatrix}+\begin{bmatrix}0\\1\end{bmatrix}u_{cp}(t)\\ y_{cp}(t)=[-1.0059\times 10^{-5}2.5594\times 10^{-4}]\begin{bmatrix}x_{1cp}(t)\\x_{2cp}(t)\end{bmatrix} \end{cases}$

(5)
$\begin{cases} \begin{bmatrix}\dot x_{1ev}(t)\\\dot x_{2ev}(t)\end{bmatrix}=\begin{bmatrix}0& 1\\-0.006& -0.4150\end{bmatrix}\begin{bmatrix}x_{1ev}(t)\\x_{2ev}(t)\end{bmatrix}+\begin{bmatrix}0\\1\end{bmatrix}u_{ev}(t)\\ y_{ev}(t)=[-2.682\times 10^{-4}6.705\times 10^{-4}]\begin{bmatrix}x_{1ev}(t)\\x_{2ev}(t)\end{bmatrix} \end{cases}$

여기서 상태변수 $x(t)$는 가제어 표준형 변환 과정에서 얻어진 가상의 물리량임에 주목할 필요가 있다. 또한, 시간의 함수인 상태변수 $x(t)$, 출력변수 $y(t)$, 제어입력 $u(t)$의 ‘$t$’는 표기의 편의상 이후에는 생략한다.

2.2 칼만 상태관측기를 갖는 슬라이딩 모드 제어기 설계

Fig. 3은 SMC의 위상평면도이다. SMC는 Fig. 3(a)에서 상태변수 $p_{1}$과 $p_{2}$를 설계자가 설정한 슬라이딩 평면을 따라 원점 $O$에 수렴시킨다. Fig. 3(b)는 서보계를 위한 SMC의 위상평면도이다. 서보계는 $p_{1}$과 $p_{2}$를 제어오차 $e(=T-T^{*})$와 그 미분값 $\dot e$으로 각각 치환함으로써 구성된다. 여기서 $T^{*}$는 $T$의 설정값이다.

Fig. 3 Phase plane of sliding mode control.
../../Resources/sarek/KJACR.2021.33.10.506/fig3.png

2차원 슬라이딩 함수 $\sigma$는 일반적으로 $\sigma =\lambda p_{1}+p_{2}$로 설정되며, 이때 $\lambda(\lambda >0)$는 슬라이딩 라인의 기울기로서 $p_{1}$이 0에 수렴하는 속도를 의미한다. 본 논문에서는 SMC에 의한 서보계 구축을 목표로 하므로 상태변수를 제어오차 $e$, 즉 $e_{o}(=T_{o}-T_{o}^{*})$와 $e_{s}(=T_{s}-T_{s}^{*})$로 하여 $\sigma$를 식(6)과 같이 설계한다.

(6)
$\sigma =\lambda e+\dot e$

제어량 $T_{o}$와 $T_{s}$를 제어하기 위한 SMC의 제어입력 $u$는 슬라이딩 모드 입력인 연속 제어입력 $u_{c}$와 도달 모드(reaching mode) 입력인 불연속 제어입력 $u_{d}$의 합으로 식(7)과 같이 구해진다.

(7)
$u=u_{c}+u_{d}$

여기서 $u_{c}$는 슬라이딩 모드 입력이므로 식(6)에서 $\dot\sigma =0$일 때의 제어입력 $u$로부터 구한다. 우선 $T_{o}$를 제어하기 위한 $u_{c}$를 설계하기 위해 식(4)의 상태 공간 모델을 일반화 하면, 식(8)로 표현된다.

(8)
$\begin{cases} \begin{bmatrix}\dot x_{1}\\\dot x_{2}\end{bmatrix}=\begin{bmatrix}0& 1\\-A_{21}& -A_{22}\end{bmatrix}\begin{bmatrix}x_{1}\\x_{2}\end{bmatrix}+\begin{bmatrix}0\\1\end{bmatrix}u\\ y =[C_{11}C_{12}]\begin{bmatrix}x_{1}\\x_{2}\end{bmatrix} \end{cases}$

식(8)의 출력방정식에서 상태변수 $x_{i}(i=1,\:2)$가 출력 $y(=T_{o})$에 미치는 기여도를 Hankel의 특이치 분해로 분석한 결과를 반영하면 근사적으로 식(9)가 얻어진다.

(9)
$y\approx C_{11}x_{1}$

제어오차 $e$는 $e=T-T^{*}$이므로 $e=C_{11}x_{1}-T^{*}$로 된다. 이 $e$와 $\dot e$을 구해 식(6)에 대입하여 $\dot\sigma$을 구하면 식(10)이 얻어진다. $\dot\sigma =0$일 때의 제어입력 $u$가 $u_{c}$이므로, 식(10)으로부터 식(11)이 구해진다. $T_{s}$를 제어하기 위한 $u_{c}$도 위에서 기술한 $T_{o}$의 경우와 동일한 방법으로 구해진다. 최종적으로 $u_{c}$는 시스템의 파라미터 값들을 대입하여 식(12)와 같이 유도된다.

(10)
$\dot\sigma =C_{11}(\lambda x_{2}-A_{21}x_{1}-A_{22}x_{2}+u)$

(11)
$u_{c}=A_{21}x_{1}+(A_{22}-\lambda)x_{2}$

(12)
$\begin{cases} u_{c₋cp}=2.33\times 10^{-5}x_{1cp}+(0.0398-\lambda_{cp})x_{2cp}\\ u_{c₋ev}=0.006x_{1ev}+(0.4150-\lambda_{ev})x_{2ev} \end{cases}$

한편, 불연속 제어입력 $u_{d}$는 Lyapunov의 제2법칙인 도달법칙으로부터 식(13)과 같이 설계된다. 여기서 계수 $K$는 스위칭 게인으로서 반복 시행으로 구해지며, 부호(signum)함수 $sgn(\sigma)$는 식(14)로 정의된다.

(13)
$\begin{cases} u_{d₋cp}=-K_{cp}sgn(\sigma_{cp})\\ u_{d₋ev}=-K_{ev}sgn(\sigma_{ev}) \end{cases}$

(14)
$sgn(\sigma)= 1;\sigma >0,\:0;\sigma =0,\:-1;\sigma <0$

식(12)의 연속 제어입력 $u_{c}$는 시스템의 파라미터와 연동되어 있지만, 식(13)의 불연속 제어입력 $u_{d}$는 시스템의 파라미터와 무관하다는 사실에 주목할 필요가 있다. 시스템의 상태변수가 외란 또는 모델의 불확실성에 의해 슬라이딩 라인을 이탈하면, 이 불연속 제어입력 $u_{d}$가 상태변수를 강제적으로 슬라이딩 라인에 구속시킴으로써 SMC는 외란과 모델의 불확실성에도 강인(robust)한 제어성능을 보이게 된다.

설계된 SMC 서보 제어계는 Lyapunov 안정 조건을 만족한다. 슬라이딩 함수 $\sigma$에 관한 Lyapunov 함수를 $V=\sigma^{2}/2$로 정의하면, $V>0$이고 $\dot V <0$이므로, 설계된 제어계는 점근 안정임을 알 수 있다.

식(12)의 $u_{c}$는 상태변수 $x$의 정보를 필요로 하지만, 이 $x$는 제어대상 모델을 식(4), 식(5)에서와 같이 가제어 표준형으로 변환된 가상의 물리량이므로 제어 시에는 상태관측기를 설계하여 이 값들을 추정한다. OCS의 압축기와 EEV의 출력방정식에는 잡음을 포함한 열부하 등의 외란이 인가되므로, 상태관측기는 이에 강인한 칼만 필터(Kalman filter)로 구성한다. 칼만 필터 설계는 백색잡음(white noise)으로 가정된 시스템 잡음 $\omega$와 관측 잡음 $\nu$를 갖는 식(15)의 선형 시스템을 대상으로 출력 $y$로부터 상태 $x$의 최소 평균 제곱 오차(minimum mean square error)의 추정값을 찾는 문제이다. 이때, 식(15)의 $A$, $B$, $C$는 상태 공간 모델인 식(4)식(5)에서 정의된 행렬이다.

(15)
$\begin{cases} \dot x =Ax+Bu+\omega \\ y=Cx+\nu \end{cases}$

식(16)은 칼만 필터를 적용한 상태관측기이다. 상첨자 기호 ‘$\hat{}$’은 추정값, $L$은 칼만 필터의 게인으로서 $L=SC'V^{-1}$로부터 구해진다. 여기서 행렬 $S$는 식(17)의 리카티 방정식(Riccati equation)으로부터 얻어진다. 식(17)에서 $V$와 $W$는 식(15)의 시스템 및 관측 잡음인 $\omega$와 $\nu$의 파워 스펙트럼 밀도이다.

(16)
$$ \left\{\begin{array}{l} \hat{\hat{x}}=A \hat{x}+B u+L(y-\hat{y}) \\ \hat{y}=C \hat{x} \end{array}\right. $$

(17)
$AS+SA'-SC'V^{-1}CS+W=0$

한편, 식(13)의 $u_{d}$는 식(14)의 부호함수로 인해 제어입력과 제어량에 채터링을 발생시킨다. 본 논문에서는 이를 억제하기 위해 경계층 두께 $\Phi(\Phi >0)$를 갖는 식(18)의 포화(saturation; sat)함수를 사용하였다.

(18)
$sat(\sigma)= 1;\sigma SUCC\Phi ,\:k\sigma ; |\sigma |\le\Phi ,\:-1;\sigma PREC -\Phi$ (단, $k =1/\Phi$)

포화함수는 $\Phi $로 인해 $T_{o}$의 과도응답을 열화시킨다. 이를 개선하고자 피드포워드 제어기를 추가로 설계한다. 피드포워드 제어기는 $\Phi $의 영향을 포함한 외란을 $\hat d$으로 식(19)와 같이 추정하여 식(20)과 같이 설계한다.

(19)
$\hat d =y-\hat y$

(20)
$u_{f}=F_{f}\hat d$

결과적으로 주 제어량인 $T_{o}$를 정밀하게 제어하기 위한 최종 제어입력 $u$는 식(8)식(20)을 더한 식(21)로 설계된다. 반면, 보조 제어량인 $T_{s}$의 경우에는 피드포워드 제어기를 생략한 $u = u_{c}+u_{d}$로 설계하였다.

(21)
$u=u_{c}+u_{d}+u_{f}$

Table 1은 본 논문에서 설계한 SMC 제어기와 설계 파라미터들을 각각 나타낸다.

Table 1. Designed parameters and gains of SMC with feedforward controller

Plant

Designed parameters and gains

$\lambda$

$K$ $\Phi$ $L$ $F_{f}$

Compressor

0.023

-6.72×106

0.93×104

[-0.0542 3.21×10-5]´

0.6280

EEV

0.0595

-2.62×106

7.0×104

[-0.0054 8.08×10-5]´

-

Fig. 4 Block diagrams for MATLAB simulations and real experiments to control $T_{o}$.
../../Resources/sarek/KJACR.2021.33.10.506/fig4.png

Fig. 4는 본 논문에서 제안한 칼만 상태관측기와 피드포워드 제어기를 갖는 압축기 제어용 SMC의 MATLAB 기반 시뮬레이션 및 실험 블록선도이다. EEV의 경우에는 $u_{f}$를 제외하면 Fig. 4와 동일한 구조이다.

3. 시뮬레이션 및 실험 결과

3.1 MATLAB 시뮬레이션 및 실험 결과

Fig. 5는 설계한 제어기의 타당성 검증에 사용된 VSRS로 구성된 OCS 실험장치의 개략도이다. Table 2Table 3은 시뮬레이션과 실험에 사용된 OCS의 주요 구성 요소에 대한 사양을 나타낸다.

Fig. 5 Diagram of experimental system.
../../Resources/sarek/KJACR.2021.33.10.506/fig5.png

Table 2. Specifications of the test unit

Component

Note

Compressor

Rotary type, 30-90 Hz, 0.86 kW

EEV

0~2,000 step, 12 V

Condenser

Air-cooled fin and tube type, 5.24 kW

Evaporator

Bare tube coil type, 2.1 kW(max.)

Refrigerant

R-22, 0.9 kg(max.)

Table 3. Specifications of the attached devices and oil

Component

Note

Inverter

4.5 kVA, 3phase, PWM, V/f=$c$ type

EEV drive

4 W, 24 V, Bipolar type

Heater

4.5 kW(max.)

Oil tank

Immersion type, 400mm$\times$400mm$\times$385mm

Oil

ISO VG 10, Velocite oil No. 6, 40 ℓ

가변속 압축기는 3상 농형 유도전동기에 의해 구동되는 로터리식 압축기로서 ‘$V/f=c$’ 타입의 인버터 주파수 변화를 통해 회전수가 제어된다. EEV의 개도는 EEV 드라이브로부터 지령을 받아 내부의 스텝모터를 조절하여 제어된다. MATLAB 기반의 실시간 제어장치(Real Time Controller; RTC)의 CPU는 $T_{o}$와 $T_{s}$의 설정값과 상태관측기를 통해 추정된 상태 추정값으로부터 제어입력을 연산한다. 연산된 제어입력은 RTC 내부의 D/A 변환기를 통해 아날로그 전압 지령으로 변환되어 인버터와 EEV 드라이브에 각각 입력된다. 실험에서 공작기계의 공작물 가공 과정에서 발생하는 열부하는 전기히터 장치를 통해 인가하였다.

시뮬레이션 및 실험에서는 지령값 변경과 계단 입력의 열부하 증․감에 따른 $T_{o}$와 $T_{s}$의 응답을 통해 SMC의 제어 강인성을 분석하였다. $T_{o}$의 설정값은 1,000 초에 30℃에서 25℃로 변경하였으며, $T_{s}$의 설정값은 $T_{s}$와 COP의 관계를 분석한 정특성 실험 결과를 토대로 선정되었으며, COP가 최대로 유지됨과 동시에 제어가 용이한 7℃로 설정하였다. 열부하 변동의 경우, 4,000 초에 정격 열부하(1.68 kW)의 10%를 증가시켰으며, 6,000 초에 다시 정격 열부하의 10%를 감소시켰다.

Fig. 6Fig. 7은 지령값 변경 및 열부하 증․감 시, 본 논문에서 설계한 SMC의 시뮬레이션 및 실험 결과를 각각 나타낸다. (a)는 $T_{o}$와 $T_{s}$의 응답을 나타내며, (b)와 (c)는 (a)의 제어입력인 인버터 주파수 $F_{i}$와 EEV 개도 지령 $V_{o}$를 각각 나타낸다. 이 결과들로부터 우선 시뮬레이션과 실제 실험에서의 제어 응답이 매우 일치함을 볼 수 있다. 이는 본 논문에서 설계한 제어기가 타당함을 보여준다.

Fig. 6의 시뮬레이션에서 $T_{o}$와 $T_{s}$는 지령값 변경 및 열부하 증․감 시에도 정상상태오차 없이 목표값에 엄밀히 수렴하였다. 다만, $T_{s}$는 지령값 변경 시 압축기의 급격한 회전수 변화로 인한 간섭 영향으로 큰 오버슈트를 갖는 과도특성을 보였다.

Fig. 7의 실험 결과들은 열부하 변동 및 실험 장치의 주변 온도가 모델링 시와 다르게 변동하는 환경에서 얻어진 것이다. 따라서 제안한 SMC는 모델 불확실성에도 강인한 제어성능을 가지는 것을 알 수 있다.

Fig. 6 Results of simulations applied SMC.
../../Resources/sarek/KJACR.2021.33.10.506/fig6.png

Fig. 7 Results of experiments applied SMC.
../../Resources/sarek/KJACR.2021.33.10.506/fig7.png

특히, 피드포워드 제어기는 열부하 변동에 의한 외란 인가 시에도 $T_{o}$를 정상상태오차 없이 목표 온도 25℃에 신속히 수렴시킴을 확인하였다. 또한 주된 제어량인 $T_{o}$가 지령값 변경 및 열부하 증․감 시에도 허용 오차 범위 $+- 0.1^{\circ}\mathrm{C}$ 이내로 엄밀히 제어됨을 확인하였다. 이를 통해 제안한 SMC는 모델의 불확실성 및 열부하 변동에 강인함을 확인하였다.

3.2 위상평면에서의 상태변수 거동 해석 및 경계층 두께 $\Phi$에 따른 제어량 응답 특성

Fig. 8(a)는 압축기, Fig. 8(b)는 EEV를 제어하기 위한 SMC의 2차원 위상평면도이며, 내부에 원점 부근을 상세히 확대하였다. Fig. 9에서 제어오차 $e$와 $\dot e$은 지령값 변경 및 열부하 증․감 시에도 위상평면의 원점(0,0) 으로 수렴함을 확인할 수 있다. 이를 통해 설계된 2차원 위상평면은 상태변수의 거동을 직관적으로 이해하기 쉬웠고, 이 거동을 분석함으로써 SMC가 모델 불확실성 및 외란에 강인함을 재확인할 수 있다.

불연속 제어입력 $u_{d}$에 포화함수를 사용하면, 경계층 두께 $\Phi$로 인해 제어량의 과도특성이 열화한다. 채터링 방지와 더불어 양호한 과도응답을 얻기 위해서는 이 두께 $\Phi$를 적절히 선정할 필요가 있다. Fig. 9Fig. 10은 $u_{d}$에 부호함수와 포함함수를 각각 적용한 시뮬레이션 결과이다. 그림에서 (a)는 $T_{o}$ 응답, (b)는 $e$와 $\dot e$의 거동을 각각 나타낸다. Fig. 9(a)에서 $T_{o}$는 설정값으로 신속히 수렴해 가지만 채터링이 발생한다. 이는 Fig. 9(b)에서 부호함수로 인해 원점 근방에서 변동하는 $e$와 $\dot e$의 거동을 분석함으로써 알 수 있다.

Fig. 8 Behaviors of state variables on phase plane.
../../Resources/sarek/KJACR.2021.33.10.506/fig8.png

Fig. 9 Results for applying signum function.
../../Resources/sarek/KJACR.2021.33.10.506/fig9.png

반면에, Fig. 10(a)는 포화함수를 적용함으로써 $T_{o}$의 응답에 채터링이 억제되었음을 보여준다. 이는 Fig. 10(b)에서 포화함수로 인해 원점으로 수렴하는 $e$와 $\dot e$의 거동을 분석함으로써 명확히 알 수 있다. 특히, 포화함수는 $|\sigma |\le\Phi$에서 $sat(\sigma)=k\sigma(단,\: k=1/\Phi)$로 정의되는데, $\Phi$가 클수록 식(13)의 $u_{d}$가 작아져 정착시간이 길어지므로 과도응답이 열화한다. 본 논문에서는 이를 보완하기 위해 피드포워드 제어기를 부가함으로써 $T_{o}$를 허용 오차 범위 $+- 0.1^{\circ}\mathrm{C}$ 이내로 신속하게 제어하였다.

Fig. 10 Results for applying saturation function.
../../Resources/sarek/KJACR.2021.33.10.506/fig10.png

3.3 칼만 상태관측기에 의한 외란 추정의 타당성

Fig. 11은 $T_{o}$ 제어기에 부가된 외란의 인가 위치 및 외란 추정 방법을 나타낸 시뮬레이션용 블록선도이다. Fig. 12는 인가된 외란(적색) $\overline{d}$와 추정된 외란(청색) $\hat d$을 시뮬레이션으로 비교한 결과이다. $\overline{d}=\hat d$이 엄밀히 성립하므로, 식(19)에 의한 외란 추정이 타당함을 알 수 있다.

Fig. 11 Disturbance estimation by simulation.
../../Resources/sarek/KJACR.2021.33.10.506/fig11.png

Fig. 12 Comparison between $\overline{d}$ and $\hat d$ by simulation.
../../Resources/sarek/KJACR.2021.33.10.506/fig12.png

3.4 PI 제어기와의 제어성능 비교 시뮬레이션 및 실험

본 논문에서는 제안한 SMC의 제어성능을 확인하기 위해 PI 제어기의 시뮬레이션 및 실험을 진행하고 그 응답들과 비교 분석하였다. PI 제어기의 시뮬레이션 및 실험은 SMC와 동일한 조건 하에서 진행하였으며, 지령값 변경 및 계단 입력의 열부하 증․감에 따른 $T_{o}$와 $T_{s}$의 응답을 분석하였다.

PI 제어기는 식(1)의 전달함수 모델을 기반으로 MATLAB tuner에 의해 설계된다. 압축기 제어용 PI 제어기 게인은 지령값 변경 시 SMC와 유사한 과도특성을 갖도록 $T_{o}$의 정착시간 $t_{s}[+- 2%]$를 1,700 초, 최대 언더슈트는 0.25℃(5.00%) 이내가 되도록 설계하였다. EEV 제어용 PI 제어기 게인은 과열도 $T_{s}$를 7℃로 제어함과 동시에 SMC에 의한 $T_{s}$의 과도응답과 유사하도록 미세 조정을 통해 구하였다. PI 제어기 적용 시 $T_{o}$의 제어기에는 적분 누적으로 인한 조작량의 포화를 방지하기 위해 안티-와인드업(anti-windup) 제어기를 별도로 부가하였다. Table 4는 시뮬레이션과 실험을 위해 설계한 PI 제어기와 안티-와인드업 게인을 각각 나타낸다.

Table 4. Designed gains of PI controller and anti-windup

Plant

Designed gains

P gain $K_{p}$ I gain $K_{i}$

Anti-windup gain $K_{a}$

Compressor

-17.0

-0.080

-8.50

EEV

-10.0

-0.200

-5.00

Fig. 13Fig. 14는 지령값 변경 및 열부하 증․감 시, PI 제어기의 시뮬레이션 및 실험 결과를 각각 나타 낸다. 이들 그림의 (a)는 $T_{o}$와 $T_{s}$의 응답, (b)와 (c)는 (a)에 대응하는 조작량 $F_{i}$와 $V_{o}$를 각각 나타낸다. 시뮬레이션과 실험의 제어 응답이 매우 일치하므로 설계한 PI 제어기의 타당성이 확인되었다. Fig. 14(a)의 $T_{s}$ 및 Fig. 14(b)에 보이는 미소한 변동(fluctuation)은 SMC의 실험 결과에서와 마찬가지로 잡음에 의한 영향으로 파악되었다.

Fig. 13 Results of simulations applied PI controller.
../../Resources/sarek/KJACR.2021.33.10.506/fig13.png

Fig. 14 Results of experiments applied PI controller.
../../Resources/sarek/KJACR.2021.33.10.506/fig14.png

Fig. 6Fig. 13의 시뮬레이션 결과를 우선 비교해 보면, SMC와 PI 제어기에 의한 $T_{o}$와 $T_{s}$는 지령값 변경 및 열부하 증․감 시에도 정상상태오차 없이 목표값에 엄밀히 수렴하였다. Fig. 6(a)Fig. 13(a)의 지령값 변경 시의 $T_{o}$를 비교한 결과, 과도응답이 설계 의도대로 유사하게 나타남을 확인하였다. 또한, 이 그림들에서 지령값 변경 시의 $T_{s}$의 응답은 압축기의 급격한 회전수 변화로 인한 간섭 영향으로 큰 오버슈트를 갖는 과도특성을 보였다. 다음으로 외란에 대한 제어 강인성 평가를 위해 열부하 증․감 시의 응답을 정량적으로 분석해 보면, 열부하 증가 시 SMC에 의한 $T_{s}$의 정착시간($t_{s}$)은 342 초, PI의 경우는 1,209 초로 확인되었다. 열부하 감소 시에는 SMC에 의한 $T_{s}$의 정착시간($t_{s}$)은 714 초였으며, PI의 경우 1,301 초로 나타났다. 이들 비교를 통해, 외란인 열부하 증․감 시 SMC가 PI 제어기보다 $T_{s}$를 목표 온도에 더 신속히 수렴시킴을 알 수 있다.

Fig. 7Fig. 14의 실험 결과를 비교해 보면, SMC와 PI 제어기의 $T_{o}$와 $T_{s}$는 지령값 변경 및 열부하 증․감 시에도 정상상태오차 없이 목표값에 수렴하였다. 외란에 대한 두 제어기의 강인성 비교 평가를 위해 열부하 증․감 시의 응답을 정량적으로 분석하였다. 우선 열부하 증가 시 SMC에 의한 $T_{o}$의 정착시간은 882 초, PI의 경우는 1,196 초였다. 최대 오버슈트는 SMC가 0.14℃, PI 제어기는 0.30℃로 확인되었다. 열부하 감소 시 SMC에 의한 $T_{o}$의 정착시간은 1,341 초, PI의 경우는 1,484 초였다. 최대 언더슈트는 SMC가 0.15℃, PI는 0.78℃로 확인되었다. 이 결과들을 통해 열부하 증․감 시의 주 제어량 $T_{o}$의 과도특성은 SMC가 PI 제어기보다 더 양호하게 나타남을 알 수 있다. 다음으로 지령값 변경 시 SMC와 PI의 $T_{s}$ 응답을 비교 분석한다. 두 제어기의 응답 모두 압축기의 급격한 회전수 변화로 인한 간섭 영향으로 큰 오버슈트를 갖는 과도특성을 보인다. SMC의 최대 오버슈트는 14.4℃, PI 제어기의 최대 오버슈트는 16.3℃로 SMC보다 13.2% 크게 나타났다. 이를 통해 SMC가 지령값 변경 시 냉매 질량유량의 급격한 변화에도 더욱 강인한 제어성능을 가지는 것을 알 수 있다. 한편, 외란인 열부하 증가 시 SMC에 의한 $T_{s}$의 정착시간은 116 초, PI의 경우 512 초로 확인되었다. 열부하 감소 시 SMC에 의한 $T_{s}$의 정착시간은 395 초, PI의 경우 959 초로 나타났다. 이를 통해 열부하 증․감 시에도 SMC가 PI 제어기보다 $T_{s}$를 목표 온도에 더 신속히 수렴시킴을 알 수 있다. 결과적으로, SMC가 PI 제어기보다 더 강인한 제어성능을 보임을 알 수 있다. 본 논문에서의 SMC와 PI 제어기 시뮬레이션 및 실험 결과들에 대한 엄밀한 정량적 비교 평가는 상대적 우수성 평가보다는 제안한 SMC의 강인한 제어성능을 보다 명확히 나타내기 위한 것임에 주목할 필요가 있다.

특히, Fig. 7Fig. 14의 지령값 변경 및 열부하 증․감 시, SMC의 $F_{i}$와 $V_{o}$가 PI 제어기의 조작량에 비해 더욱 안정적임을 알 수 있다. SMC에 의한 조작량의 안정성은 불연속 제어입력 $u_{d}$에 의해 보장된다. 불연속 제어입력 $u_{d}$는 외란과 모델 불확실성에 의해 제어오차가 발생한 경우 상태변수를 신속히 슬라이딩 라인에 구속시킨다. 이러한 제어 원리로 SMC는 제어오차 발생 시 짧은 시간 내에 큰 조작량을 인가함으로써 PI 제어기보다 더 양호한 과도특성을 갖는다.

뿐만 아니라 실험장치 주변의 외기 온도는 실험이 진행되는 8,000 초 동안 압축기의 발열로 인해 지속적 으로 상승하였다. 이로 인해 제어대상의 공칭모델은 모델 불확실성을 포함함에도, 설계한 SMC는 제어량을 설계사양에 잘 추종시키고 있음을 보인다. 모델 불확실성에 대한 SMC와 PI 제어기의 제어성능을 비교한 결과, SMC가 PI 제어기보다 더 강인함을 보였다. 이를 통해 본 논문에서 설계한 SMC는 외란과 모델 불확실성에 대한 제어 강인성은 물론, 속응성도 동시에 확보할 수 있음을 확인하였다. 이러한 SMC의 강인성은 공칭모델의 파라미터에 독립적으로 설계되는 SMC의 불연속 제어입력 $u_{d}$에 의해 보장된다.

본 논문에서는 OCS를 대상으로 SMC의 외란과 모델 불확실성에 대한 제어 강인성을 면밀히 검토하였다. SMC의 제어 강인성 확인을 위해 산업계에서 널리 사용되는 PI 제어기를 설계한 후, SMC와 동일한 조건 하에서 시뮬레이션과 실험을 진행하고, 그 결과를 상호 비교하였다. PI 제어기는 지령값 변경 시의 주 제어량 $T_{o}$의 응답이 SMC와 유사하도록 설계된 후, 열부하 변동과 같은 외란 인가 시의 제어 강인성에 주목하였다. 두 제어기의 강인성을 정량적으로 비교 분석한 결과, 본 논문에서 제안한 SMC가 PI 제어기보다 외란 및 모델 불확실성에 더 강인함을 확인하였다. 특히, 제안한 SMC 설계법은 기존의 SMC 설계법과는 달리 2차원 위상평면을 갖도록 설계함으로써 상태변수 거동 분석이 한층 용이해, 채터링 방지 및 정밀한 응답 특성 확보가 가능함을 확인하였다.

4. 결 론

본 논문에서는 VSRS에 대한 기존 SMC의 문제점을 해결하여 강인한 제어성능을 확보하기 위해 실용적인 상태 공간 모델을 토대로 2차원 위상평면을 갖는 슬라이딩 라인을 설계하였다. 또한, 상태 추정을 위해 칼만 상태관측기를 구성하였다. 특히, 불연속 제어입력의 포화함수가 목표 제어량의 과도특성을 열화시키는 부작용을 줄이기 위해 피드포워드 제어기를 추가로 설계하였다. 제안된 제어기는 VSRS 기반의 OCS를 제어대상으로 한 시뮬레이션과 실험을 통해 그 타당성이 검증되었으며, PI 제어기와의 제어성능을 비교함으로써 그 유효성을 확인하였다. 이를 통해 얻은 주요 결론은 다음과 같다.

(1) 2차원 슬라이딩 라인을 갖는 SMC는 위상평면에서 상태변수의 거동 분석이 용이해 채터링 방지 및 정밀한 제어 응답을 확보하였다.

(2) 칼만 필터를 갖는 상태관측기는 제어입력 연산에 필요한 상태변수 및 외란을 엄밀하게 추정하였다.

(3) 피드포워드 제어기는 포화함수로 인한 정착시간 증가 문제를 해결함으로써 제어량의 과도특성을 현저히 개선하였다.

(4) SMC와 PI 제어기의 제어성능을 비교한 결과 SMC가 PI 제어기보다 외란 및 모델 불확실성에 더 강인함을 확인하였다.

References

1 
Kim J. G., Han C. H., Jeong S. K., 2020, Disturbance observer-based robust control against model uncertainty and disturbance for a variable speed refrigeration system, International Journal of Refrigeration, Vol. 116, pp. 49-58DOI
2 
Kim J. G., Han C. H., Jeong S. K., 2019, Robust Control for Variable Speed Refrigeration System Using a Disturbance Observer Based on Optimization, Korean Journal of Air-Conditioning and Refrigeration Engineering, Vol. 31, No. 5, pp. 193-205Google Search
3 
Jeong S. K., Kwon T. E., 2019, Robust Linear Quadratic Gaussian Controller Design for Oil Coolers Based on a State Space Model, Korean Journal of Air-Conditioning and Refrigeration Engineering, Vol. 31, No. 3, pp. 130-139Google Search
4 
Jeong S. K., 2019, Optimum robust control for an oil cooler system with variable speed drive of machine tools, Journal of the Korean Society for Power System Engineering, Vol. 23, No. 3, pp. 18-26DOI
5 
Oh H. R., Yang S. W., Jeong S. K., 2020, Optimal Robust Servo Control for a VSRS with Model Uncertainty and Disturbance, 2020 SAREK Summer Annual Conference, pp. 66-69Google Search
6 
Han C. H., Kim J. G., Jeong S. K., 2018, Model Reduction Using the Hankel Singular Value Decomposition for Variable Speed Refrigeration System, 2018 SAREK Winter Annual Conference, pp. 173-176Google Search
7 
Yang S. W. Jeong S. K., 2020, Robust Temperature Controller Design Based on a Sliding Mode with Optimal Switching Hyper-plane for a Variable Speed Refrigeration System, Korean Journal of Air-Conditioning and Refrigeration Engineering, Vol. 33, No. 3, pp. 101-112Google Search
8 
Lee J. E., Jeong S. K., 2020, Robust Temperature Control for Variable Speed Refrigeration System Based on Sliding Mode Control Using a State Observer with Kalman Filter, 2020 SAREK Winter Annual Conference, pp. 23-26DOI
9 
Li H., Jeong S. K., 2006, An Empirical Model for Decoupling Control of a Variable Speed Refrigeration System, The Korean Society for Power System Engineering, Vol. 10, No. 3, pp. 81-87Google Search
10 
Li H., Jeong S. K., Yoon J. I., You S. S., 2008, An empirical model for independent control of variable speed refrigeration system, Applied Thermal Engineering, Vol. 28, No. 14-15, pp. 1918-1924DOI
11 
Li H., Jeong S. K., You S. S., 2009, Feedforward control of capacity and superheat for a variable speed refrigeration system, Applied Thermal Engineering, Vol. 29, No. 5-6, pp. 1067-1074DOI
12 
Jeong S. K., Han C. H., Li H., Wahyu K. W., 2018, Systematic design of membership functions for fuzzy logic control of variable speed refrigeration system, Applied Thermal Engineering, Vol. 142, pp. 303-310DOI
13 
Cao J. P., Jeong S.K., Jung Y. M., 2014, Fuzzy logic controller design with unevenly distributed membership function for high performance chamber cooling system, Journal of Central South University, Vol. 21, No. 7, pp. 2684-2692DOI
14 
Aprea C., Mastrullo R., Renno C., 2004, Fuzzy control of the compressor speed in a refrigeration plant, International Journal of Refrigeration, Vol. 27, No. 6, pp. 639-648DOI
15 
Vadim U. I., 1977, Variable structure systems with sliding modes, IEEE Transactions on Automatic Control, Vol. 22, No. 2, pp. 212-222DOI
16 
Vadim U. I., 1993, Sliding Mode Control Design Principles and Applications to Electric Drives, IEEE Transactions on Industrial Electronics, Vol. 40, No. 1, pp. 23-36DOI
17 
Alfonso D., Gianluca L., Ignazio M., Alessandro P., 2004, Second-Order Sliding-Mode Control of DC Drives, IEEE Transactions on Industrial Electronics, Vol. 51, No. 2, pp. 364-373DOI
18 
Kim N. H., Park Y. K., Sun J. E., Shin S. J., Min B. H., Park H. J., Kang S. H., Hur H., Ha M. Y., Lee M. C., 2018, Robust sliding mode control of a vapor compression cycle, International Journal of Control, Automation and Systems, Vol. 16, No. 1, pp. 62-78DOI