A Study on Air-gap Control for Transverse Flux Permanent Magnet Type Magnetic Levitation Electromagnet System

횡자속 영구자석형 자기부상전자석 시스템의 공극제어에 관한 연구

Abstract

In this paper, we proposes a study on air gap control for magnetic levitation of transverse flux permanent magnet electromagnets. In general, mechanical systems have a high failure rate of bearings. Bearings in particular are problematic because they have high surface wear rate and degradations. To solve this problem, replacing the bearing with a magnetic levitation electromagnet system can provide lightweight and efficiency improvements. However, precise air gap control is essential to control the magnetic levitation electromagnet system. Therefore, in this paper, we identify the instable cause of gap control through a mathematical modeling and verify through experiment a control algorithm that can use compensation.

1. 서 론

자기부상 전자석 시스템(magnetic levitation electromagnet system)은 전자기력을 이용하여 물체를 부상시키는 장치로서 1970년 이후 활발하 게 연구가 진행되어 왔다.

기술의 발전으로 회전기기에서 정지기기로 전기 접점의 무접점화, 수동제어에서 능동 제어 시스템 으로 발전되어 가고 있다. 이러한 현대산업의 요 구조건에 일치하는 자기부상 전자석 시스템으로 물체를 부상시키기 때문에 회전기기의 기계적 마 찰에 의한 마모를 줄여 시스템의 수명 증가 및 유 지보수 비용의 절감과 마찰저항, 마찰열 등을 무 시함으로써 고속에서도 높은 가속 성능을 낼 수 있다는 이점을 갖는다[1]-[3]. 이러한 이점으로 전 자기력을 통한 부상 및 추진하는 자기부상열차의 연구가 활발히 진행되고 있으며, 이송 및 분류 시 스템에서도 소규모 물류창고의 효율적인 공간 활 용과 장비 고장이나 노후화와 같은 유지 관리 비 용 문제를 자기부상 기술을 채택하여 개발하고 있 는 등 다양한 산업 분야에 응용되고 있다[4][5].

하지만 자기부상 전자석 시스템은 본질적으로 비선형성이 강한 불안정한 시스템으로 시스템 안 정화 및 신뢰성 확보를 위한 정밀한 공극 제어 알 고리즘이 필수적이며, 전자석과 부상체 사이에 작용하는 전자기력에 대한 정보를 명확히 알고 있어 야 부상체의 위치와 무게가 변해도 안정된 제어 성능을 가질 수 있다. 또한 자기부상 전자석 시스 템의 선형화 과정에서 필연적으로 모델링 오차가 발생하며 이러한 모델링 오차나 외란 등에 대하여 양호한 제어성능을 만족할 수 있어야 한다[6]-[10]. 따라서 자기부상 시스템의 실제 적용 가능성을 보 장하면서 동시에 장기적으로 발생할 수 있는 문제 를 해결할 수 있는 기술적인 해법이 필요하다.

Fig. 1 Overall system block diagram of prototype magnetic bearing

본 논문에서는 Fig. 1과 같이 자기부상 전자석 시스템을 구성하여 공극 제어 시 불안정성 원인의 분석 및 시스템의 자기력의 특성상 존재하는 비선 형성 해석을 위해 Taylor 급수 전개를 이용하여 선형화하여 수학적 모델을 수립하고 제어 알고리 즘 및 디지털 제어기 구현을 통해 공극 제어 알고 리즘 신뢰성을 검증하였다.

2. 본 론

2.1 자기부상 전자석의 수학적 모델링

2.1.1 자기회로의 동작 특성

자기부상 전자석 시스템의 자기 등가회로를 Fig. 2에 나타내었으며, 전자석 E-1에 의해 발생하는 자속은 식(1)과 같이 나타내었다.

Fig. 2 Equivalent circuit of magnetic levitation electromagnet system

\(\psi _ { 1 } = \frac { i N } { 2 R _ { g } } = \frac { i N } { 2 } \frac { \mu _ { 0 } A } { g - x }\)

여기서 은 전자석 E-1의 자속, 는 전자석 E-1 의 코일에 흐르는 전류, 은 코일의 턴 수, 은 진공 상태에서의 투자율, 는 자기저항, 는 전 자석의 면적, 는 자로의 길이, 그리고 는 자기 부상 전자석 시스템의 이동자 변위이다.

전자석 E-1 코일의 쇄교자속은 코일의 턴 수 및 코일에 흐르는 전류에 의해 식(2)로 표현된다.

\(\lambda _ { 1 } = \frac { N ^ { 2 } i } { 2 } \frac { \mu _ { 0 } A } { g - x }\)

코일에 흐르는 전류와 쇄교자속의 비례상수로 인 덕턴스는 식(2)로부터 식(3)과 같이 다시 정리된다.

\(L = \frac { N ^ { 2 } } { 2 } \frac { \mu _ { 0 } A } { g - x }\)

이때, 자기부상 전자석 시스템의 공극에 비해 이동자 변위 가 작을 경우, 테일러급수 전개를 적용하여 식(4)와 같이 표현할 수 있으며, 이를 통 해 첫 번째와 두 번째 항만 고려한 인덕턴스는 식 (5)으로 근사화하여 표현이 가능하다.

\(\frac { 1 } { g - x } = \frac { 1 } { g } ( 1 + \frac { x } { g } + ( \frac { x } { g } ) ^ { 2 } + ( \frac { x } { g } ) ^ { 3 } + \cdots )\)

\(L = L _ { 0 } ( 1 + \frac { x } { g } ) , L _ { 0 } : \frac { N ^ { 2 } \mu _ { 0 } A } { 2 g }\)

여기서, 에 의해 공극에 저장된 자계 에너지  ′는 식(6)과 같이 표현할 수 있다.

\(W ^ { \prime } = \frac { 1 } { 2 } L i ^ { 2 }\)

식(6)의 자계 에너지  ′로부터 발생된 전자기 력 와 그에 의한 이동자 변위 에 대한 인덕턴 스 의 변화는 식(7) 및 식(8)과 같다.

\(F = \frac { \partial W ^ { \prime } } { \partial x } = \frac { \partial L } { \partial x } \frac { i ^ { 2 } } { 2 }\)

\(\frac { \partial L } { \partial x } = \frac { L _ { 0 } } { g }\)

식(7) 및 식(8)을 통하여 자계 에너지에서 발생 하는 전자기력의 크기를 식(8)로 나타낼 수 있다.

\(F = \frac { L _ { 0 } } { g } \frac { i ^ { 2 } } { 2 }\)

식(9)에서 볼 수 있듯이 자기부상 전자석 시스 템의 공극과 전류 변화에 따른 전자기력의 크기는 비선형 관계를 가짐을 알 수 있다. 공극에서 발생 한 전자기력의 경우 전자석 입력 전류의 제곱에 비례하며, 전류의 입력이 일정한 경우 공극 변위 에 따라 반비례함을 보여준다.

2.1.2 이동자 변위 에 대한 전자기력

앞 절에서 유도한 공극과 전자기력의 비선형 관계식은 자기부상 전자석 시스템의 이동자 변위 의 인덕턴스 변화를 무시하여 출력된 결과이며, 이동자 변위 는 전자기력에 직접적인 영향을 미 치므로 이에 대한 추가적인 분석이 필요하다.

식(4)에서 자기부상 전자석 시스템의 안정성에 영향을 미치는 세 번째 항 이상의 성분에 대한 인 덕턴스는 식(10)과 같이 정의할 수 있다.

\(L = L _ { 0 } ( 1 + \frac { x } { g } + ( \frac { x } { g } ) ^ { 2 } )\)

이에 따라 자계 에너지 및 자계 에너지의 변화 에 따른 전자기력은 식(11) 및 식(12)와 같다.

\(W = \frac { 1 } { 2 } i ^ { 2 } L _ { 0 } [ 1 + \frac { x } { g } + ( \frac { x } { g } ) ^ { 2 } ]\)

\(F = \frac { L _ { 0 } } { 2 g } i ^ { 2 } + \frac { L _ { 0 } } { g ^ { 2 } } i ^ { 2 } x\)

식(12)에서 볼 수 있듯이 이동자 변위 에 의 한 전자기력이 자기부상 전자석 시스템에 불안정 한 특성을 야기하며, 식(12)를 특정 동작점에서 선형화된 전자기력으로 표현하면 식(13)으로 나타 낼 수 있다.

\(F = k _ { i } i + k _ { x } x\)

여기서, 는 전류 강성이며 는 변위 강성이다.

2.1.3 자기부상 전자석의 운동 방정식

흡인식자기부상 전자석 시스템의 운동방정식과 블록도를 식(14) 및 Fig. 3에 나타내었다.

Fig. 3 Block diagram of magnetic levitation electromagnet system

\(F = m \frac { d x ^ { 2 } } { d t } + f _ { d }\)

여기서, 은 이동자 무게이며, 는 외란을 뜻한다.

3. 자기부상 전자석의 디지털 제어기 설계

3.1.1 자기부상 전자석의 전압 방정식

자기부상 전자석 시스템은 자체적으로 부상 상태의 유지가 불가능하기 때문에 공극 센서를 이용 하여 파악한 이동자의 위치를 통한 피드백 제어기 의 사용이 필수적이다. 이러한 피드백 제어기를 포함한 자기부상 전자석 시스템의 전체 제어 블록 도는 Fig. 4과 같다. Fig. 4에서 알 수 있듯이 자 기부상 전자석 시스템의 전압 방정식은 식(15)와 같이 정리할 수 있다.

Fig. 4 Control block diagram of magnetic levitatioin electromagnet system

\(V _ { s } = R _ { s } i _ { s } + \frac { d \lambda _ { s } } { d t }\)

여기서, 는 인가 전압, 는 쇄교자속, 는 코 일에 흐르는 상전류, 는 고정자 권선 저항이다.

이때, 식(15)의 자기회로가 선형적이고 누설자 속이 없다고 가정하면 식(16)과 같이 표현할 수 있다.

\(V _ { s } = R _ { s } i _ { s } + L _ { s } ( x ) \frac { d i _ { s } } { d t } + i _ { s } \frac { d L _ { s } ( x ) } { d x } \frac { d x } { d t }\)

자기부상 전자석 시스템의 전류에 따른 자속 포화가 없을 경우 전류 변화에 따른 인덕턴스를 일정한 계수로 표현이 가능하며 이동자 변위량이 작을 때 전압 방정식은 식(17)과 같이 간략화를 할 수 있으며, 이를 라플라스 변환하면 식(18)을 얻을 수 있다.

\(V _ { s } = R _ { s } i _ { s } + L _ { s } ( x ) \frac { d i _ { s } } { d t }\)

\(\frac { I _ { s } ( s ) } { V _ { s } ( s ) } = \frac { 1 } { R _ { s } + L _ { s } s } = \frac { 1 / R _ { s } } { 1 + \tau _ { s } s }\)

여기서,   로 전자석의 전기적 시정 수이다.

3.1.2 자기부상 전자석의 공극 제어기

식(17)을 보면 자기부상 전자석 시스템의 공극 제어를 위해서는 미분 요소가 필수적이다. 이를 반영한 공극 제어기의 블록도는 Fig. 5와 같다. PD 제어기 적용 시 공극 제어기를 위한 전달함수 는 다음과 같다.

Fig. 5 PID Control block diagram of magnetic levitation electromagnet system

\(G _ { d s ) } = K _ { p q } + T _ { d } s\)

공극 제어기를 적용한 자기부상 전자석 시스템 안정성 계산을 위한 폐루프 전달함수는 식(20)과 같다.

\(\left. \begin{array} { l }{ G _ { c ( s ) } ^ { c } = \frac { x ( s ) } { x ^ { * } ( s ) } = \frac { ( K _ { p q } + T _ { d } s ) \times k _ { i } } { ( m s ^ { 2 } - k _ { x } ) + ( K _ { p q } + T _ { d } s ) \times k _ { i } } } \\ { = \frac { ( K _ { p q } + T _ { d } s ) \times k _ { i } } { m s ^ { 2 } + s T _ { d } k _ { i } + ( K _ { p q } k _ { i } - k _ { x } ) } } \end{array} \right. \)

식(20)으로부터 특성 방정식은 식(21)과 같다.

\(m s ^ { 2 } + s T _ { d } k _ { i } + ( K _ { p q } k _ { i } - k _ { x } ) = 0\)

식(21)을 통해 자기부상 전자석 시스템의 안정 성을 알아보기 위한 조건은 아래와 같이 표현할 수 있다.

1.  및 가 0인 경우

특성 방정식의 근은 식(22)로 나타나며, 식(22) 및 Fig. 6에서 알 수 있듯이 조건1의 경우 우반 평면에 극점이 존재하므로 시스템이 불안정하게 된다.

Fig. 6 Root locus analysis. (a) condition, (b) condition2

\(s = \pm \sqrt { 4 m k _ { x } }\)

2. 가 0이고     인 경우

특성 방정식의 근은 식(23)과 같이 나타나며, 식(23) 및 Fig. 6에서 볼 수 있듯이 조건2의 경우

두 극점이 허수축에 위치하기 때문에 안정된 영역 에서의 사용이 불가능하다.

\(s = \pm j \sqrt { 4 m ( K _ { p q } k _ { i } - k _ { x } ) }\)

3. 가 양수이고     인 경우

특성 방정식의 근은 식(24)와 같고, 식(24) 및 Fig. 7에서 볼 수 있듯이 시스템의 두 극점이 좌 반 평면에 위치하기 때문에 안정적인 시스템의 동 작이 가능하다.

Fig. 7 Root locus analysis of condition3

\(s = - T _ { d } k _ { i } \pm \sqrt { ( T _ { d } k _ { i } ) ^ { 2 } - 4 m ( K _ { p q } k _ { i } - k _ { x } ) }\)

이에 따라 자기부상 전자석 시스템의 안정적인 동작을 위해서는 미분 성분이 필요함을 알 수 있다.

4. 실 험

4.1.1 실험장치 구성

횡자속 영구자석형 자기부상 전자석의 공극 제 어 성능 평가를 위해 자기부상 전자석 시스템을 구성하였다. 실험 장치는 Fig. 8에서 볼 수 있듯 이 이동자의 공극 제어를 위한 링 형태의 권선을 갖는 2개의 전자석, 단상 AD/DC/AC 컨버터, 공 극 위치 측정을 위한 갭 센서 및 제어 보드로 구 성하였다. 본 실험에 사용된 제어 보드의 프로세 서는 TMS320VC33을 이용하였으며, A/D 컨버터 를 통해 공극 위치 및 전자석의 입력 전류를 디지 털 값으로 변환하였다. 또한 갭 센서는 와전류 형 센서로 자세한 사양은 Table 1과 같다.

Fig. 8 Overall experiment setup

Table 1. System of gap sensor

4.1.2 실험 및 결과

본 논문에서는 횡자속 영구자석형 자기부상 전 자석의 PID 공극 제어기 적용 시의 공극 제어 성 능 평가를 위하여 일정 공극 제어, 공극 지령 변동 및 바이어스 전류 변동 시의 실험을 진행하였 으며, 실험은 이동자 무게 29 [kg], 고정자 저항 3.39[Ω], 전류 강성() 76.878 [N/A], 변위 강성 () 38,439.02 [N/m], 정격 전류 8 [A], 전류 제 어기 주기 100 [us], 스위칭 주파수 10 [kHz]의 조 건에서 초기 위치 7.6 [mm], 기준 공극 4 [mm]로 설정하여 진행했다.

Fig. 9 Constant air gap control test of magnetic levitation electromagnet system

Fig. 10 Air gap reference change test of magnetic levitation electromagnet system

초기 위치 7.6 [mm]에서 기준 공극 4 [mm]로 지령을 주었을 때 공극과 전류의 파형을 Fig. 9에 나타내었다. Fig. 9에서 알 수 있듯이, 실제 공극이 지령 공극을 잘 추종하고 있으며 실제 전류 역 시 지령 전류에 수렴하고 있음을 알 수 있다. Fig. 10는 운전 중 공극 지령 변동에 따른 실제 공극 의 응답 특성을 나타내었으며, Fig. 11은 일정 공 극 제어 중 바이어스 전류 변동 시 공극과 전류 파형을 보여주고 있다. Fig. 10과 Fig. 11에서 알 수 있듯이 공극 지령이 스텝 형태로 변동이 생겨 도 실제 공극이 값을 잘 추종하고 있으며 일정 공 극 제어 시 바이어스 전류의 변동에도 불구하고 기준 공극 값인 4 [mm]로 수렴하고 있음을 확인 할 수 있다.

Fig. 11 Bias current variation test of magnetic levitaion electromagnet system

5. 결 론

본 논문에서는 자기부상 전자석 시스템을 통해 횡자속 영구자석형 자기부상 전자석의 공극 제어 를 위한 기본 개념을 설계하였으며, 특히 디지털 제어기 설계 및 시스템의 안정성과 관련된 내용의 수학적 모델링을 수립하고 실험을 진행하였다.

그 결과, 자기부상 전자석 시스템 공극 제어를 위한 수학적 모델 구성 시 시스템의 불안정성의 원인을 분석하여 미분 요소를 포함한 제어 알고리 즘 수립 및 디지털 제어기를 구현하고자 하였고, 또한 자기부상 전자석 시스템의 성능 평가를 통해 일정 공극 제어 및 공극 지령 변동, 바이어스 전 류 변동에 따른 공극 제어 성능 평가 등의 실험을 통해 공극 제어 알고리즘의 신뢰성을 검증하였다.

사 사

본 연구는 연구개발특구진흥재단의 ‘지역혁신 메 가프로젝트’ 사업(과제명: 수소연료전지 기반 하이 브리드 분산 전기추진시스템을 활용한 커뮤터기 기 술, 과제번호: 2023-DD-UP-0026 )과 2023년도 정부(산업통상자원부)의 재원으로 한국에너지기술평 가원의 지원을 받아 수행된 연구임(141518667, 최 대이륙중량 200kg급 비행체용 순정격출력 30kw 급 연료전지 파워팩 시스템 개발).