HOME > Browse Articles > Archive

Archive

Journal of the Korean Society for Precision Engineering - Vol. 39 , No. 2

[ REGULAR ]
Journal of the Korean Society for Precision Engineering - Vol. 39, No. 2, pp. 123-129
Abbreviation: J. Korean Soc. Precis. Eng.
ISSN: 1225-9071 (Print) 2287-8769 (Online)
Print publication date 01 Feb 2022
Received 17 Aug 2021 Revised 17 Nov 2021 Accepted 22 Nov 2021
DOI: https://doi.org/10.7736/JKSPE.021.087

무급유 터보 냉동기의 압축기를 지지하는 자기베어링 제어기의 효율적 최적화에 관한 연구
권은상1 ; 노명규1, # ; 이남수2 ; 백성기2 ; 박영우1
1충남대학교 메카트로닉스공학과
2LG전자 에어솔루션연구소

A Study on the Efficient Optimization of Controller for Magnetic Bearings Supporting Oil-Free Turbo-Chiller Compressor
Eunsang Kwon1 ; Myounggyu Noh1, # ; Namsoo Lee2 ; Seongki Baek2 ; Young-Woo Park1
1Department of Mechatronics Engineering, Chungnam National University
2H&A, LG Electronics Inc.
Correspondence to : #E-mail: mnoh@cnu.ac.kr, TEL: +82-42-821-6877


Copyright © The Korean Society for Precision Engineering
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Funding Information ▼

Abstract

Magnetic bearings are being actively adopted by the turbo-chiller industry because of their higher efficiency during partial load, quieter operation, and smaller footprint than that which machines with ball bearings provide. Since magnetic bearings are open-loop unstable, feedback control is necessary. In the industry, traditional PID-based control is preferred to model-based control, because of its simplicity. When traditional control algorithms are used, significant resources are required to obtain and tune control parameters, which is an impediment to the widespread use of magnetic bearing technology in the industry. In this paper, we propose a mixed optimization method by combining genetic algorithm and sequential quadratic programming. To obtain the initial guess to be used for the mixed optimization, a phase-margin maximization algorithm is also proposed, based on the rigid-body model of the system. Mixed optimization results in suitable control parameters in less than 2.8% of the time it takes a genetic algorithm only to find similar solutions. The proposed optimization also ensures the robustness of the control parameters. The output sensitivity measured from a prototype compressor with magnetic bearings confirms the validity of the control parameters.


Keywords: Active magnetic bearing, Turbo-Chiller, Control design, Optimization, Genetic algorithm, Sequential quadratic programming
키워드: 능동 자기베어링, 터보 냉동기, 제어기 설계, 최적화, 유전자 알고리즘, 순차적 2차 프로그래밍

1. 서론

자기베어링은 자기력을 이용해 회전축을 부상함으로써 비접촉 회전이 가능하도록 하며, 시스템의 신뢰도 및 효율 향상, 모니터링 기능 제공, 고장허용 작동 등 기존 베어링에 비해 여러 장점을 가지고 있다. 특히, 터보 냉동기 분야에서 부분 부하 효율 향상, 운전 소음 저감, 설치 면적의 감소 등으로 인해 자기베어링을 장착한 무급유 시스템에 대한 관심이 증대되고 있으며, 많은 터보 냉동기 업체가 자기부상 무급유 제품을 출시하고 있다.1

능동 자기베어링은 내재적 불안정성으로 궤환 제어가 필수적이다. 센서가 감지한 회전축의 위치를 입력으로 제어기가 제어 명령을 생성하면 증폭기가 그에 상응하는 전류를 베어링에 공급하여 비접촉 부상이 유지되도록 한다.2 자기베어링의 제어기를 대략적으로 분류하면 모델 기반과 경험 기반 제어기로 나눌 수 있다. 모델 기반은 시스템의 모델에 근거해 제어기 형태와 제어 파라미터를 결정하는 방법이다.3,4 경험 기반 제어기 설계는 실험을 통해 얻은 데이터와 경험을 기반으로 시행착오(Trial and Error)를 거쳐 제어 파라미터를 결정하는 방법이다.5,6 시스템을 정확하게 표현할 수 있는 수학적 모델 부재, 제어기 차수 증가로 인한 구현의 어려움 등으로 인해 현장에서는 경험 기반 제어기가 주로 사용된다.

경험 기반 제어기 설계는 설계자의 전문성과 경험이 요구되며, 시행착오에 많은 시간과 자원이 필요하다. 이러한 문제를 해결하기 위해 최적화 방법을 이용하여 경험 기반 제어기 설계를 자동화하고자 하는 연구가 진행되었다.7,8 제어기 설계를 위한 최적화는 최적점이 다수 존재하는 전역 최적(Global Optimization) 문제로서 기존의 연구는 모두 유전자 알고리즘(Genetic Algorithm)을 사용하였다. 전역 최적화의 가장 큰 문제점은 전역 최적점으로의 수렴 여부를 확신하기 어렵다는 것이다. 설령 전역 최적점으로 수렴하더라도 계산에 소요되는 시간이 매우 길어 현장 적용이 어려울 수 있다. 예를 들어, 유전자 알고리즘의 경우 인구수(Population)와 세대수(Generation)를 증가시키면 전역 최적점으로 수렴하도록 할 수 있지만 계산 비용의 증가를 감수해야만 한다.

최적화에서의 첫 단계는 초기값 설정이다. 자기베어링 시스템은 개루프 불안정(Open-Loop Unstable)하므로, 최적화가 제대로 시작되기 위해 안정성을 제공하는 초기 제어기 선정이 필요하다. 본 논문의 첫 번째 목표는 회전축의 강체 운동 모델에 기반하여 안정성이 확보된 초기 부상 제어기를 적절히 설계할 수 있는 과정을 제시하는 것이다. 그 다음으로 유전자 알고리즘과 기울기 기반 최적화 방법을 혼합하여 강인성이 보장된 제어 파라미터를 효율적으로 탐색하는 방법을 제시하고자 한다. 혼합 최적화 방법으로 최적 수렴 시간이 2.8% 이하로 감소함을 확인하였고, 실험을 통해 최적 제어기 적합성을 검증하였다.


2. 이론적 배경
2.1 Plant 전달함수

자기베어링 시스템에서 제어기의 제어대상인 Plant는 회전축, 베어링 액츄에이터, 전류증폭기, 센서 등으로 구성된다. 자기베어링에 의해 지지된 회전축의 동역학을 전달함수 형태로 표현하면 식(1)과 같다.

Grotor=Kact1ms2+Ks+k=1NflexAr,ks2+2ζr,kωr,ks+ωr,k2(1) 

회전축의 유연모드를 고려하는 경우 Nflex개의 유연모드에 대해 모드강도 Ar, 모드댐핑 ζr, 모드주파수 ωr 등으로 유연모드를 모델링할 수 있다. 유연모드에 대한 파라미터는 유한요소 모델로부터 구하거나 실험을 통한 시스템 식별과정(System Identification)으로부터 구할 수 있다. 회전축의 강체 운동은 음의 부호를 가지는 개루프(Open-Loop) 스프링 상수 Ks 와 질량 m 으로 모델링된다. 코일에 전류를 공급하는 증폭기는 대부분 2차 시스템의 특성을 가지며 전달함수는 다음 식(2)와 같이 표현할 수 있다.

Gamp=Kampωamp2s2+2ζampωamps+ωamp2(2) 

여기서 ωamp. 는 증폭기의 주파수 대역을 의미하며, ζamp. 는 댐핑이다. 센서의 주파수 대역은 제어범위보다 매우 크므로 고정된 이득 Ksensor로 모델링할 수 있다. 따라서, Plant 전달함수는 식(3)으로 같이 구할 수 있다.

Gplant=KsensorGrotorGamp(3) 

Fig. 1은 Plant 전달함수의 한 예로서 유한요소해석을 통해 구한 결과이다. 회전축은 총 3개의 유연모드를 가지는 유연축으로 모델링되었으며, 전류증폭기의 주파수 대역은 680 Hz, 댐핑은 0.7이다.


Fig. 1 
An example of a plant transfer function which includes the dynamics of rotor, amplifier, and sensor

2.2 제어기

자기베어링의 내재적 불안정성을 극복하기 위해 다양한 형태의 제어기가 제시되었으나, 상업용 제품에서 가장 널리 사용되는 제어기는 비례-적분-미분(Proportional-Integral-Derivative, PID) 형태이다. 미분기의 고주파 이득을 제한하고 구현이 가능하도록 비례-미분기를 변형하면 식(4)와 같이 제어기를 표현할 수 있다.

GPID=KPτzs+1τps+1+KIs(4) 

식(4)의 첫 항은 위상선행(Phase Lead) 필터라고도 하며, α = τz/τp 가 위상선행과 고주파 이득을 결정한다. 위상선행은 진동 관점에서 댐핑과 관련이 있으며, 위상선행을 증가하기 위해 2차 위상선행 필터를 사용하기도 한다.9,10

GPID=KPτz2s2+2ζzτzs+1τp2s2+2ζpτps+1+KIs(5) 

일반적으로 PID 제어기와 함께, 센서의 고주파 잡음을 저감하는 저역 통과(Lowpass) 필터, 유연모드를 제어하기 위한 노치 필터가 제어기에 사용된다.

Gcon=GPIDGLPFGNF(6) 

Fig. 2는 자기베어링 제어 시스템을 블록 선도로 보여준다.


Fig. 2 
Block diagram of magnetic bearing control. The plant consists of bearing actuator (Including amplifier), rotor and sensor


3. 강체모드 제어기 설계
3.1 강체모드 제어기 설계 과정

회전축의 운전속도가 유연모드에 비해 현저하게 낮거나 초기 부상을 위한 제어기 설계의 경우 회전축의 강체 운동만을 고려할 수 있다. 또한 유연모드 제어기 설계의 초기값으로도 강체 모드 제어기가 활용될 수 있다. 제어기 설계의 목표는 ISO 14839-3에서 제시하는 출력민감도의 최소화인데,11 단일입력-단일출력(SISO)인 경우 위상여유도(Phase Margin)가 출력민감도와 직접적 관계가 있으므로, 위상여유도를 최대한 확보하도록 설계한다. 위상여유도는 Gain Cross-Over 주파수, 즉 시스템 개루프 전달함수의 크기가 1이 되는 주파수에서 개루프 전달함수의 위상과 180o 간의 차이를 의미한다.

Plant 전달함수의 전체 이득을 Kf = KactKampKsensor 로 정의하고, 시상수 형태로 표현하면 식(7)과 같다.

Gplant=Kf/Ksτm2s2-11τamp2s2+2ζampτamps+1(7) 

강체 운동의 시상수와 증폭기의 시상수는 다음 식(8)과 같이 정의된다.

τm=mKs,τamp=1ωamp(8) 

Gain Cross-Over 주파수는 식(9)를 만족하는 주파수이다.

GplantjωgcGconjωgc=1(9) 

수치적 방법으로 식(9)의 해를 구할 수 있겠으나, 분수함수 및 제곱근 형태이므로 보다 수월하게 해를 구하기 위해 식(9) 좌변의 분모와 분자를 각각 D(gc), N(gc) 라고 하면 식(10)이 되어, 다항식 형태로 정리할 수 있다.

Njωgc2=Djωgc2(10) 

제어기에서 위상선행 필터만을 고려한다면 이 식은 Ωgc=ωgc2 에 대해 6차 방정식이며, Newton-Raphson 방법과 같은 과정을 통해 해를 구할 수 있다. 만일 증폭기의 주파수 대역이 매우 높고 (ωampωgc), 위상선행 필터에서 극점이 영점에 비해 매우 크다면 (τpτz), 식(10)은 다음 식(11)과 같이 2차 방정식으로 근사화할 수 있다.

aΩgc+bΩgc+c=0a=κf2Kp2Tz2-Tm2b=2κf2Kp22Zz-1Tz-2Tmc=κf2Kp2-1(11) 

식(11)에서 κf = Kf/|Ks| 이며, 대문자변수는 소문자변수의 제곱을 의미한다. 즉, Tm=τm2, Tz=τz2, Zz=ζz2 이다.

Gain Cross-Over 주파수를 구하면, 다음 식(12)로 위상여유도를 구할 수 있다.

ϕPM=Gampjωgc+GLFjωgc(12) 

만일 위상여유도가 0보다 작으면 시스템의 불안정성을 의미하며, 최적화의 구속 조건으로 활용할 수 있다.

3.2 강체모드 제어기 설계 결과

3.1절의 과정을 통해 강체모드 제어기의 설계를 수행하였다. 대상은 50마력 압축기의 반경 방향 베어링이며 회전축의 질량은 7.3 kg, 개루프 스프링 상수 -871 N/m, 액츄에이터 상수 126 N/A이다. 전류증폭기의 주파수 대역은 800 Hz로 가정하였다. 비례이득은 12이며, 적분기는 사용하지 않았다. 초기 제어기는 시행착오를 거쳐 결정하였다.

Table 1Fig. 3은 제어기 설계의 결과이다. 제어기 설계의 목표는 위상여유도를 최대화하여 안정성을 높이는 것이지만, 이와 함께 외란에 대한 강인성을 향상하는 것이다. 안정성과 강인성을 동시에 평가하는 도구로 출력민감도를 사용한다. 출력민감도함수는 식(13)으로 정의되며, SO14839-3은 출력민감도 함수의 최댓값인 ∥S()∥가 3을 넘지 않도록 권고하고 있다.

S=11+GplantGcpm(13) 
Table 1 
Design result for rigid-body controller
Lead filter parameter Initial Optimized
Zero time constant (τz) [ms] 1.768 1.276
Pole time constant (τp) [ms] 0.589 0.284
Zero damping (ζz) 0.7 1.037
Pole damping (ζp) 0.7 1.999


Fig. 3 
Output sensitivity function of phase-margin optimized control compared with initial control

Fig. 3은 초기 제어기와 최적화된 제어기의 출력민감도를 비교하고 있으며, 최댓값은 3.67에서 2.32로 감소하여 기준을 만족한다.


4. 유연모드 제어기 최적화
4.1 유연모드 제어기 최적화 과정

많은 경우 회전축의 유연모드는 자기베어링 제어기 설계에 지대한 영향을 미친다. 운전 속도가 유연모드보다 낮더라도 제어 주파수 범위 내에 있는 유연모드는 센서 잡음, 운전속도의 배수 성분 등에 의해 진동을 유발하며, 제어 안정성을 저하한다. 유연모드에 대응하는 가장 일반적인 방법은 노치 필터의 사용이다. 노치 필터는 식(14)와 같은 전달함수로 표현된다.

GNF=s2+2ζnzωnfs+ωnf2s2+2ζnpωnfs+ωnf2(14) 

노치 주파수는 유연모드 주파수에 따라 결정되며, 노치 필터의 설계는 노치의 폭을 결정하는 ζnp 와 노치의 깊이를 결정하는 ζnz 를 정하는 과정이라고 할 수 있다.

위상선행 필터와 노치 필터를 동시에 사용하면 강체모드 제어기에서와 같은 위상여유도 최적화 방법을 사용할 수 없다. 유연모드에 의해 다수의 Gain Cross-Over 주파수가 있을 수 있으며, 노치 필터로 인한 차수의 증가로 수치적 불안정성이 증가한다. 또한 출력민감도 최소화는 여러 개의 지역 최소점(Local Minima)을 가질 수 있어, 전역 최소화(Global Minimization)가 필요하다.

흔히 사용되는 전역 최소화 방법으로 유전자 알고리즘이 있다. 일정 크기의 인구에 대해 비용함수에 해당하는 적합성(Fitness)을 계산하고, 진화와 돌연변이를 거쳐 세대를 반복함으로써 최적을 찾는 방법이다. 유전자 알고리즘의 가장 큰 문제점은 전역 최적치를 찾기 위해 인구가 크거나 세대 수가 많아야 하기 때문에 수렴하기 위해 많은 계산 시간의 필요하다는 것이다. 다음의 사례 연구를 통해 최적화 과정의 문제점을 파악한다.

대상은 340 kW 냉매압축기용 자기베어링으로서 회전축 질량 72 kg, 개루프 스프링 상수 -1,958 N/mm, 액츄에이터 상수103 N/A, 전류증폭기 주파수 대역 680 Hz의 특성을 가진다. Plant 전달함수는 Fig. 1과 같으며 3개의 유연모드를 확인할 수 있다. 유연모드의 주파수는 각각 393, 766, 1,550 Hz이다. 자기베어링 제어기는 식(6)과 같이 PID 기반 위상선행 필터, 노치 필터, 저역 통과 필터로 구성되어 있으며, 저역 통과 필터의 주파수는 1,500 Hz로 설정되었다. 제어기 변수는 위상선행 필터의 4개, 각 유연모드별 노치 필터당 2개로 총 10개이다.

이를 유전자 알고리즘을 이용해 최적화할 때 걸리는 계산 시간을 인구수에 대비해 구한 결과가 Fig. 4에 나와 있다. 계산 시간은 거의 인구수에 비례한다. 인구수가 50일 때 46.5초인거에 비해 인구수가 25,600일 때 거의 10,000초의 시간이 걸려 최적값을 찾아낸다. 계산은 i9-7920X (2.96 GHz) CPU, 64 GB RAM, 64 bit Win-10 컴퓨터에서 MATLAB Global Optimization Toolbox의 유전자 알고리즘을 사용하였다.


Fig. 4 
Convergence time with respect to the size of population (Red). Normalized standard deviation of solutions with respect to the population size

유전자 알고리즘의 특성상 인구수가 작을 때는 지역 최솟값으로 수렴할 가능성이 존재한다. 유전자 알고리즘을 여러 번 수행하여 각 파라미터의 평균과 표준편차를 구하면 정규화된 표준편차의 최댓값을 다음 식(15)로 구할 수 있다.

σ¯max=maxkstdipk,imeanipk,i(15) 

여기서, pk,ik 번째의 제어 파라미터의 i 번째 최적해이다. 즉, σ¯max는 유전자 알고리즘의 해가 얼마나 퍼져 있는지를 나타낸다. Fig. 4의 파란색 그래프는 유전자 알고리즘을 10회 수행하여 얻은 σ¯max로 인구수가 작을 때는 최적화 결과값이 40% 이상 다를 수 있다. 수렴 시간의 증가에도 불구하고 인구수가 10,000을 넘어야 최적값의 표준편차가 20% 이하로 줄어들어 전역 최솟값에 수렴하는 경향을 보인다.

제어의 강인성 측면에서 제어 파라미터의 변화에 대한 최적점의 변화가 크지 않아야 한다. 파라미터의 섭동(Perturbation)에 따라 비용함수가 크게 바뀐다면, 제한 조건 근처에 있거나 지역 최소점을 둘러싸는 영역이 매우 작음을 의미하며 파라미터 변화에 따라 제어 결과가 크게 바뀔 수 있다. 이를 확인하기 위해 유전자 알고리즘을 수행하여 얻은 최적점에서 무작위로 10% 이내 만큼 파라미터를 변화시켜 초기값을 설정하고, 기울기 기반 최소화 방법인 Sequential Quadratic Programming (SQP)를 사용하여 최적화를 수행하였다. 만일 SQP에 의해 도달한 최소점이 유전자 알고리즘 결과와 유사하면, 최소점을 둘러싸는 파라미터 영역이 10% 정도의 파라미터 변화를 감당할 수 있음을 의미한다.

Table 2는 인구수 50으로 유전자 알고리즘을 10회 수행하여, 각 최적점의 강인성을 확인한 결과이다. 유전자 알고리즘으로 구한 출력민감도의 최솟값이 두 번째 열에 나와 있다. 유전자 알고리즘을 통해 얻은 파라미터값을 10% 이내로 무작위 변화시켜 초기값으로 설정한 후, SQP를 통해 최소화를 진행하였다. Table 2의 세 번째 열은 SQP에 의해 구한 파라미터와 유전자 알고리즘에서 얻은 파라미터 간의 차이를 비로 나타낸 것이며 다음 식(16)으로 정의된다.

Table 2 
Robustness of local optima found by GA
Sol. No. GA Optim. Param. Var. SQP Optim.
1 2.01 0.467 18.5
2 2.06 0.093 2.58
3 1.82 0.442 6.00
4 1.95 0.081 2.06
5 1.98 0.333 3.63
6 2.04 0.276 6.76
7 1.91 0.238 3.63
8 2.00 0.164 3.05
9 1.93 0.511 16.1
10 1.94 0.084 2.28

εmax=maxkpk,GA-pk,SQPpk,GA(16) 

여기서 pk,GA 는 유전자 알고리즘에 의해 구한 k 번째 제어 파라미터이며, pk,SQP 는 SQP로 구한 k 번째 제어 파라미터이다. 1번 해의 경우 파라미터가 46% 이상 변하였으며 출력민감도는 거의 불안정성에 가까운 18.5의 결과를 얻었다. 유전자 알고리즘 해가 강인성을 보인 경우는 2, 4, 10번으로 모두 파라미터의 변화가 10% 이내이며, SQP로 구한 출력민감도가 유전자 알고리즘의 출력민감도와 크게 다르지 않았다.

이상의 결과에 근거하여, 본 논문에서는 유전자 알고리즘과 기울기 기반 최적화 방법을 혼합한 효율적 최적화 과정을 제안한다. Fig. 5는 혼합 최적화 방법의 개략적인 순서도로서, 시작 단계에서 최소한의 인구수로 유전자 알고리즘을 수행하여 최적 파라미터를 구한다. 최적 파라미터로부터 일정범위(예를 들어 10%) 내에 무작위로 파라미터를 변화시켜 다수의 초기값을 구하고, SQP를 통해 해를 구한다. 식(15)를 이용해 정규화된 표준편차의 최댓값을 구하여 0.2 이하가 되는지 판단하여, 만족하지 않는 경우 다시 유전자 알고리즘으로 새로운 최적해를 찾아 반복한다. Fig. 4에서 확인할 수 있는 바와 같이 정규화된 표준편차의 최댓값이 0.2 이하로 되려면 인구수가 최소 10,000을 넘어야 하고, 수렴에 필요한 계산 시간은 4,495초이다. 본 논문에서 제안하는 혼합 최적화 과정을 통해 동일한 수렴 결과를 얻는데 걸리는 계산 시간은 평균 126초로 유전 알고리즘만을 사용하였때에 비해 35.7배 빠르게 결과를 얻을 수 있었다.


Fig. 5 
Flow chart of mixed optimization algorithm

4.2 유연모드 제어기 최적화 결과

앞서 제시한 혼합 최적화 과정을 냉매압축기 자기베어링 제어기에 적용하였다. 자기베어링이 장착된 압축기는 Fig. 6에 나와 있다. 시스템 식별을 위해 Sine Sweep 방법이 사용되었다. 제어기에서 특정한 주파수의 입력을 인가하고, 출력에서 해당 주파수 성분을 DFT를 이용해 추출하여 주파수 응답을 구하는 방식이다.


Fig. 6 
Prototype compressor for 340 kW chiller. The cooling capacity is 500RT

제어기는 위상선행 필터, 노치 필터, 저역 통과 필터로 구성되어 있다. Table 3은 시행착오 과정을 거쳐 수동으로 설정한 제어 파라미터와 최적화 과정을 통해 얻은 제어 파라미터를 비교한다. Fig. 7은 시제품에서 Sine Sweep 방법을 이용하여 측정한 출력민감도이다. 경험에 근거하여 수동으로 설정한 제어기의 경우 출력민감도의 최댓값이 3.43로 ISO 14839-3의 기준을 초과하는데 반해 혼합 최적화 과정으로 얻은 제어기를 적용한 결과 최댓값이 2.13로 기준을 만족하였다.

Table 3 
Optimization result for flexible-mode control
Parameter Initial Optimized
Lead
filter
Zero T.C. (τz) [ms] 1.768 0.798
Pole T.C. (τp) [ms] 0.589 0.080
Zero damp. (ζz) 0.7 1.414
Pole damp. (ζp) 0.7 0.623
Notch 1 Zero damp. (ζnz1) 0.005 0.05
Pole damp. (ζnp1) 0.5 0.148
Notch 2 Zero damp. (ζnz2) 0.005 0.004
Pole damp. (ζnp2) 0.5 0.391
Notch 3 Zero damp. (ζnz3) 0.005 0.001
Pole damp. (ζnp3) 0.5 0.5


Fig. 7 
Output sensitivity measured from a prototype compressor test rig using sinesweep tests. Manually tuned control (Blue) vs. Optimized control (red)


5. 결론

자기베어링 제어기로 산업 현장에서 널리 사용되고 있는 경험 기반 제어기의 설계는 전문성과 경험이 요구된다. 경험 기반 제어기를 최적화하기 위한 기존의 연구는 수렴에 걸리는 시간으로 인해 현장 적용에 어려움이 있다. 전역 최적화와 지역 최적화를 혼합한 효율적 최적화를 통해 수렴 시간을 2.8% 이내로 단축하였다. 혼합 최적화의 초기값으로 활용할 수 있는 강체모드 제어기의 위상여유도 최적화 방안도 함께 제시하였다. 본 논문에서 제시한 제어기 설계 자동화 과정이 보다 유용성을 갖기 위해서는 압축기 이외에도 펌프나 송풍기와 같은 다양한 시스템에서의 검증이 필요하다.


NOMENCLATURE
G : Transfer Function
S : Output Sensitivity Function
s : Laplace Variable
m : Rotor Mass
Kact : Actuator Constant
Ks : Open-Loop Stiffness
Nflex : Number of Flexible Modes
Ar : Modal Strength
ζr : Modal Damping
ωr : Mode Frequency
Kamp : Amplifier Gain
KP : Proportional Gain
KI : Integral Gain
τz : Zero Time Constant of Lead Filter
τp : Pole Time Constant of Lead Filter
ζz : Zero Damping of Lead Filter
ζp : Pole Damping of Lead Filter
ζnz : Zero Damping of Notch Filter
ζnp : Pole Damping of Notch Filter
ϕPM : Phase Margin
σ¯max : Maximum Normalized Standard Deviation
εmax : Maximum Normalized Parameter Difference
p : Optimizing Parameters

Acknowledgments

이 논문은 2021년도 산업통상자원부의 재원으로 한국에너지기술평가원의 지원을 받아 수행된 연구임(No. 20208901010010).


REFERENCES
1. Jin, S., Cha, J., Chang, J., Choi, S., and Ahn, H., “Controlling the Fluid Induced Instability of A Supercritical CO2 Compressor Supported by Magnetic Bearing,” Journal of the Korean Society for Precision Engineering, Vol. 37, No. 10, pp. 737-742, 2020.
2. Schweitzer, G. and Maslen, E. H., “Magnetic Bearings. Theory, Design, and Application to Rotating Machinery,” Springer, 2009.
3. Maslen, E. H. and Sawicki, J. T., “Mu-Synthesis for Magnetic Bearings: Why Use such a Complicated Tool” Proc. of the ASME International Mechanical Engineering Congress and Exposition, pp. 1103-1112, 2007.
4. Noshadi, A., Shi, J., Lee, W. S., Shi, P., and Kalam, A., “Robust Control of an Active Magnetic Bearing System Using H and Disturbance Observer-Based Control,” Journal of Vibration and Control, Vol. 23, No. 11, pp. 1857-1870, 2017.
5. Spirig, M., Schmied, J., Jenckel, P., and Kanne, U., “Three Practical Examples of Magnetic Bearing Control Design Using a Modern Tool,” Journal of Engineering for Gas Turbines and Power, Vol. 124, No. 4, pp. 1025-1031, 2002.
6. Sun, J., Zhou, H., Ma, X., and Ju, Z., “Study on PID Tuning Strategy based on Dynamic Stiffness for Radial Active Magnetic Bearing,” ISA Transactions, Vol. 80, pp. 458-474, 2018.
7. Wei, C. and Söffker, D., “Optimization Strategy for PID-Controller Design of AMB Rotor Systems,” IEEE Transactions on Control Systems Technology, Vol. 24, No. 3, pp. 788-803, 2015.
8. Liu, B., Sjöberg, J., and Laiho, A., “Optimization-Based Radial Active Magnetic Bearing Controller Design and Verification for Flexible Rotors,” Proceedings of the Institution of Mechanical Engineers, Part I: Journal of Systems and Control Engineering, Vol. 230, No. 4, pp. 339-351, 2016.
9. Messner, W. C., Bedillion, M. D., Xia, L., and Karns, D. C., “Lead and Lag Compensators with Complex Poles and Zeros Design Formulas for Modeling and Loop Shaping,” IEEE Control Systems Magazine, Vol. 27, No. 1, pp. 44-54, 2007.
10. Kurvinen, E., Fittro, R., and Maslen, E., “Improving Compressor Surge Performance with Advanced Control,” Proceedings of the Institution of Mechanical Engineers, Part I: Journal of Systems and Control Engineering, Vol. 230, No. 7, pp. 672-679, 2016.
11. ISO 14839-3:2006, “Mechanical Vibration-Vibration of Rotating Machinery Equipped with Active Magnetic Bearings-Part 3: Evaluation of Stability Margin,” 2006.

Eunsang Kwon

MS Student in the Department of Mechatronics Engineering, Chungnam National University. His research interest is Automatic control.

E-mail: eunsang284@gmail.com

Myounggyu Noh

Professor in the Department of Mechatronics Engineering, Chungnam National University. His research interest is magnetic bearings and control.

E-mail: mnoh@cnu.ac.kr

Namsoo Lee

Chief Research Engineer at Air Solution Research Center, H&A, LG Electronics Inc.. His research interest is turbo-chiller design.

E-mail: ns.lee@lge.com

Seongki Baek

Senior Research Engineer at Air Solution Research Center, H&A, LG Electronics Inc.. His research interest is turbo-chiller control.

E-mail: baek.seongki@lge.com

Young-Woo Park

Professor in the Department of Mechatronics Engineering, Chungnam National University. His research interest is smart sensing and actuator.

E-mail: ywpark@cnu.ac.kr