Ⅰ. 서론
영구자석 동기전동기는 높은 출력밀도와 효율 확보가 가능하나, 제작에 필요한 재료의 가격이 비싸고 설계가 유도전동기 대비 다소 어려운 문제점을 가지고 있다. 따라서 효율 및 유지보수 편의성 등이 모두 고려된 최적 의전 동기 개발 및 관련 제어연구가 필요한 실정이다. 뿐만 아니라, 실질적인 전동기에 의한 구동은 좁은 정격영역에서의 최고효율의 증대 요구와 전체 전기구동 영역에서의 평균 효율 증대, 평균 출력의 증대 요구로 이어지고 있다. 이러한 움직임에 의해 영구자석이 필요 없는 릴럭턴스 전동기 (Reluctance Motor)가 하나의 대안으로 검토되고 있다.
그림 1~2의 개발 방향 및 이슈(Issue)와 같이 현재의 희토류 영구자석 전동기를 대체할 수 있는 미래기술 개발을 준비하게 되었고, 희토류 저감형 전동기와 탈 희토류 전동기의 기술 선점을 요구하고 있다[1].
그림 1. 구동용 전동기의 개발 방향
Fig. 1. DC Direction of development of motors for drive
그림 2. 전동기 성능향상을 위한 개발 이슈
Fig. 2. Development issues for improving the performance of motors
본 논문에서는 희토류 영구자석 전동기를 대체할 수 있는 미래기술 개발과, 희토류 저감형 전동기와 탈 희토류 전동기의 기술 선점을 요구하는 시대적 이슈 (Issue) 에 맞춰 영구자석이 필요 없는 스위치드 릴럭턴스 전동기 (Switched Reluvtance Motor, SRM)에 대해 기계언어(C 언어)를 통한 모델링과 그에 따른 SRM의 특성을 연구하고자 한다.
II. PSIM 프로그램에서의 SRM Modeling
PSIM 시뮬레이션 프로그램에서는 3상 SRM을 모델링한 library를 그림 3과 같이 제공하여 SRM의 변수값들을 입력할 수 있고 구동용 컨버터와 회전자 위치검출 회로를 추가적으로 이용하여 그림 4와 같이 SRM 구동 및 제어회로를 모델링할 수 있다[2].
그림 3. PSIM 프로그램에서 제공하는 3상 SRM library
Fig. 3. SRM library provided by PSIM program
그림 4. PSIM 프로그램에서 제공된 SRM library를 이용한 전동기 구동 및 제어 모델링
Fig. 4. Motor drive and control modeling using the SRM library provided in the PSIM program
따라서 PSIM 및 특정 프로그램에서 제공된 library를 이용하여 전동기 구동 및 제어회로를 설계하여 그에 따른 결과를 확인할 수 있다. 그러나 제공된 특정 library 만 이용할 수 있기 때문에 여러 변수값을 가지고 있는 전동기를 적용하는데 한계가 있으며 프로그램 구입 및 업데이트에 해당되는 경제적 비용도 고려해야 되는 문제점이 있다.
III. 기계언어를 이용한 SRM Modeling
본 논문에서는 II절에서 언급한 기존 시뮬레이션용 프로그램의 문제점을 고찰하여 기계언어(C 언어)를 통해 비선형 인덕턴스 특성을 갖는 SRM을 모델링 하였으며 다음의 내용들을 가정하였다. 또한 모델링 과정에서 미분 방정식의 해를 구하기 위해 Runge-Kutta 방법을 사용하였으며 제어부의 샘플링 시간은 50μs로 하였다.
①SRM의 입력측 DC링크 전압의 맥동은 무시한다.
②SRM은 자기포화가 없다고 가정한다.
③비대칭 브리지 컨버터에서 스위치 및 다이오드의 손실은 무시한다.
1. Modeling을 위한 수식화[3-5]
기계언어(C 언어)를 통해 모델링에 사용된 SRM의 한 상의 전압 방정식은 다음과 같다.
\(\begin{aligned} V_{s}=& R_{e} i+L(\theta, i) \frac{d i}{d t}+i \frac{d L(\theta, i)}{d t} \\ & \pm M\left(\theta, i^{\prime}\right) \frac{d i^{\prime}}{d t} \end{aligned}\) (1)
\(\frac{d L(\theta, i)}{d t}=\frac{d L(\theta, i)}{d \theta} \omega\) (2)
\(\frac{d i}{d t}=\frac{d i}{d \theta} \omega\) (3)
식 (2), (3)을 식 (1)에 대입하면,
\(\begin{aligned} V_{s}=& R_{e} i+L(\theta, i) \frac{d i}{d \theta} \omega \\ &+i \frac{d L(\theta, i)}{d \theta} \omega \pm M\left(\theta, i^{\prime}\right) \frac{d i^{\prime}}{d \theta} \omega \end{aligned}\) (4)
여기서,
Vs : DC 인가 전압.
Re : 고정자의 권선 저항 및 컨버터 소자 저항의 합.
L(θ,i) : 인덕턴스.
i : 여자상의 순시전류.
i' : 앞상의 여자 순시전류.
θ : 고정자와 회전자의 상대위치.
M(θ, i) : 인접한 상의 상호 인덕턴스.
ω : 회전자의 회전 각속도.
이고, 이상적인 스위치를 사용하며, 자기포화가 없는 인덕턴스를 가정하면 식 (4)는 다음과 같이 간략화 할 수 있다.
\(V_{s}=R i+L(\theta) \frac{d i}{d t}+i \frac{d L(\theta)}{d \theta} \omega\) (5)
(R : 고정자 권선 저항.)
회전자 극수 Nr이 4개일 때 한 상의 인덕턴스의 변환 주기는 90°가 된다. 자기포화를 무시하여, 전류의 크기와 관계없이 인덕턴스가 항상 일정한 최대값 Lmax와 최소값 Lmin을 갖는다고 가정하면, 한 상의 인덕턴스는 고정자와 회전자의 극호(pole-arc) βs, βr에 따라 고정자극과 회전자극이 겹치기 시작하는 시점을 θ = 0으로 하였을 때 다음과 같다.
\(\begin{aligned} &L_{1}(\theta)= \\ & \begin{cases}L_{\min }+K \theta & \left(\theta \leqq \beta_{s}\right) \\ L_{\max } & \left(\beta_{s} \leqq \theta \leqq \beta_{r}\right) \\ L_{\min }-K\left(\theta-\beta_{r}-\beta_{s}\right) & \left(\beta_{r} \leqq \theta \leqq \beta_{r}+\beta_{s}\right) \\ L_{\min } & \left(\beta_{r}+\beta_{s} \leqq \theta \leqq \tau\right)\end{cases} \end{aligned}\) (6)
또한 j번째 상의 인덕턴스는 기준 인덕턴스로부터 식 (7)과 같이 얻을 수 있다.
\(L_{j}(\theta)=L_{1}\left[\theta-\frac{\tau}{N_{r}}(j-1)\right]\) (7)
여기서, K는 인덕턴스 파형 중 증가하는 부분의 인덕턴스의 기울기이며, 다음과 같다.
\(K=\frac{L_{\max }-L_{\min }}{\beta_{s}}\) (8)
\(\begin{aligned} T_{j} &=\frac{1}{2} \frac{d L(\theta)}{d \theta} i_{j}^{2} \\ &=\left\{\begin{array}{cl} \frac{1}{2} K i_{j}^{2} & \left(0 \leqq \theta \leqq \beta_{s}\right) \\ 0 & \left(\beta_{s} \leqq \theta \leqq \beta_{r}\right) \\ -\frac{1}{2} K i_{j}^{2} & \left(\beta_{r} \leqq \theta \leqq \beta_{r}+\beta_{s}\right) \\ 0 & \left(\beta_{r}+\beta_{s} \leqq \theta \leqq \tau\right) \end{array}\right. \end{aligned}\) (9)
또한, 자기포화가 없다고 가정하면 j번째 상에서 발생되는 한 상의 토크는 식 (9)와 같이 표현되며, 3상 SRM 의 합성 토크는 식 (10)과 같이 구해진다.
\(T=\sum_{j=1}^{3} T_{j}=\frac{1}{2} \sum_{j=1}^{3} \frac{d L_{j}(\theta)}{d \theta} i_{j}^{2}\) (10)
선형적인 부하를 갖는 전동기의 기계 방정식은
\(T=T_{L}+B \omega+J \frac{d \omega}{d t}\) (11)
로 표현되며, 식 (5), (10), (11)로부터 SRM의 동특성 방정식을 구하면 식 (12)와 같다.
\(\begin{aligned} \frac{d i_{j}}{d t} &=\frac{V_{j}}{L_{j}(\theta)}-\frac{R i_{j}}{L_{j}(\theta)}-\frac{d L_{j}(\theta)}{d \theta} \frac{i_{j}}{L_{j}(\theta)} \\ \frac{d \omega}{d t} &=\frac{1}{J}\left[-B \omega+\frac{1}{2} \sum_{j=1}^{3} \frac{d L_{j}(\theta)}{d \theta} i_{j}^{2}-T_{L}\right] \\ \frac{d \theta}{d t} &=\omega \end{aligned}\) (12)
2. 기계언어를 이용한 SRM Modeling 결과[6]
기계언어(C 언어)를 통해 모델링한 SRM의 특성을 살펴보기 위해 OriginLab 프로그램을 이용하여 데이터 들을 그래프화 하였다. 그림 5는 본 논문에서 3상 6/4 pole의 특정 변수를 가지고 있는 SRM의 모델링 결과로인덕턴스 Profile을 보여주고 있다[4]. 식 (6)과 (7)을 통해 회전자 위치에 따른 인덕턴스의 최대값 Lmax와 최소값 Lmin을 가지고 있는 SRM의 인덕턴스 특성을 확인할 수 있다. 또한 3상 SRM의 합성 토크 식 (10)에 대한 결과로 그림 6과 같이 무부하시 구동에 따른 토크 리플이 포함된 SRM의 전형적인 토크 그래프를 보여주고 있다.
그림 5. 3상 6/4 pole SRM의 인덕턴스 profile
Fig. 5. Inductance profile of 3 phase 6/4 pole SRM
그림 6. SRM의 토크 파형 (무부하)
Fig. 6. Torque waveform of SRM (no load)
그림 7과 8은 기준속도(Wref) 증감 및 감속에 대한속도 변화시 속도응답(Wreal)과 기준전류(Iref) 파형의 결과를 보여주고 있다. 1, 000→1, 500rpm, 1, 500→ 1, 000rpm으로 기준속도 변화시 설계된 제어기에 의해 기준전류의 변화량과 기준속도를 추종하는 응답의 결과를 보여주고 있다.
그림 7. 속도 변화시 속도 및 I_ref(기준전류) 파형 (1, 000rpm → 1, 500rpm)
Fig. 7. DC motor closed loop control system
그림 8. 속도 변화시 속도 및 I_ref(기준전류) 파형 (1, 500rpm → 1, 000rpm)
Fig. 8. DC motor closed loop control system
모델링된 SRM을 이용하여 구동 및 제어특성을 확인하기 위해 제어부 및 구동부를 포함한 전체 시스템에서 히스테리시스 전류제어 방식과 전압 PWM제어 방식을 수행하여 그에 따른 특성을 살펴보았다. 그림 9, 10은 기준속도 2, 000rpm (중첩각 0°)일 때 히스테리시스 전류제어 방식과 전압 PWM제어 방식에 대한 각각의 전류 파형 모델링 결과를 보여주고 있다.
그림 9. SRM의 상전류 파형 (히스테리시스 전류제어)
Fig. 9. Phase current waveform of SRM (hysteresis current control)
그림 10. SRM의 상전류 파형 (전압 PWM 제어)
Fig. 10. Phase current waveform of SRM (voltage PWM control)
그림 9에서는 히스테리시스 전류제어기를 사용했을 때 한 상의 전류가 Flat-Top하게 흐르는 것을 보여주므로써 설계된 히스테리시스 전류제어기의 Upper-Lower band 안에서 전류제어가 이루어지고 있음을 알 수 있다. 전류가 히스테리시스 하위 밴드(Lower band)에 도달하게 되면 전류를 평활하게 하기위해 제어기는 스위치를 On하게 되고 다시 전류가 히스테리시스 상위 밴드 (Upper band)에 도달하게 되면 스위치를 Off하는 방식으로 전류가 Flat-Top 하게 흐르고 있다. 반면 그림 10 에서는 전압 PWM시 그림 9의 상전류 결과와 달리 전류가 Flat-Top 하지 못하고 감소하는 부분이 발생되는 결과를 볼 수 있다.
그림 11, 12는 각 제어방식의 구동에 따른 SRM의 상전압 파형으로 인덕턴스 Profile이 양(+)의 기울기가 존재하는 영역에서 스위칭 소자의 On-Off에 따른 양의 값을 가지는 상전압을 볼 수가 있다. 또한 인덕턴스 Profile이 음의 기울기를 가지는 영역에서 빠른 전류의 소호를 위해 스위칭 소자를 Off 시킴으로써 음(-)의 상전압 특성을 확인할 수 있다.
그림 11 SRM의 상전압 파형 (히스테리시스 전류제어)
Fig. 11. Phase voltage waveform of SRM (hysteresis current control)
그림 12. SRM의 상전압 파형 (전압 PWM 제어)
Fig. 12. Phase voltage waveform of SRM (voltage PWM control)
SRM의 토크는 식 (9)와 같이 나타낼 수 있으므로 인덕턴스의 기울기가 일정하다면 전류를 Flat-Top하게 흐르게 함으로써 일정한 토크를 얻을 수가 있게 된다. 따라서 그림 9~12의 리모델링의 결과로 전압 PWM 제어 방식에 비해 히스테리시스 전류제어를 통해 토크 리플의 저감의 효과를 얻을 수 있다.
V. 결론
본 논문에서는 기존 시뮬레이션용 프로그램의 문제점을 고찰하여 기계언어(C 언어)를 통해 비선형 인덕턴스특성을 갖는 SRM을 모델링을 하였고 제어부와 구동부를 포함한 전체 시스템에 대한 특성에 대해 고찰하였다. 따라서 기존 특정 프로그램에서 제공하는 전동기 library 만 이용하는 제약과 여러 변수들을 가지고 있는 전동기들의 특성을 확인할 수 없는 문제들을 본 연구를 통해 연구된 모델링 프로그램을 통해 다양한 특성의 연구 결과를 얻을 수 있을 것으로 사료된다.
※ 이 연구는 2021년도 광주대학교 대학 연구비의 지원을 받아 수행되었음.
참고문헌
- Y. S. Chae, Technology Development Trend of Reluctance Motors for Automobiles, Global Auto News 2013.
- PSIM User's Manual, Powersim Inc.
- T. J. E Miller, Switched Reluctance Motors and Their Control, Oxford University press, 1993.
- J. M. Noh and J. H. Kim, "Automated Inductance Measurement of a Switched Reluctance Motor Using Voltage Integration Method," Trans. KIEE, Vol. 64, No. 8, pp. 1180~1185, 2015. DOI: http://dx.doi.org/10.5370/KIEE.2015.64.8.1180
- S. I. Jeong, "Modeling of Switched Reluctance Motor (SRM) Drive and Control System using Rotor Position Information Sensor", The Journal of The Institute of Internet, Broadcasting and Communication (JIIBC), Vol. 21, No. 3, pp.137-142, 2021. DOI: https://doi.org/10.7236/JIIBC.2021.21.3.137
- Y. H. Yoon, "Control of SRM with Modified C-dump Converter in Cooling System of Automobiles," Trans. KIEE, Vol. 66, No. 8, pp. 1297~1302, 2017. DOI: http://doi.org/10.5370/KIEE.2017.66.8.1297