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

  1. (Dept. of Electrical Engineering, Chonnam National University, Republic of Korea.)



In-Wheel Motor, Power Factor, Permanent Magnet Vernier Motor, Torque Density, Winding Factor

1. 서 론

직접 구동 인휠 모터(Direct-Drive In-Wheel Motor)는 모터가 바퀴 내부에 직접 탑재되어 별도의 동력 전달 장치 없이 독립 구동이 가능한 구동 방식이다[1]. 직접 구동 인휠 모터는 감속기 등 동력 전달 장치의 부재로 추가적인 배터리 탑재 공간 및 차량 내부 공간을 확보할 수 있으며, 동력 전달 과정에서의 손실을 최소화할 수 있는 장점을 가진다. 하지만 바퀴 내부에 탑재되는 모터, 인버터, 브레이크 등의 구동 요소들로 인해 바퀴 자체의 무게가 증가하고, 제한된 공간에 구동 요소들을 배치해야 하므로 이를 보완하기 위한 높은 토크 및 출력 밀도가 요구되고, 동시에 낮은 토크 리플과 넓은 구동 속도 범위 또한 요구된다[2] [3].

따라서 영구자석 동기전동기(Permanent Magnet Synchronous Motor, PMSM)보다 높은 토크 밀도를 가지는 영구자석 버니어 전동기(Permanent Magnet Vernier Motor, PMVM)를 직접 구동 인휠 모터로 적용하기 위한 연구가 진행되고 있다[4] [5]. PMVM은 기어링 효과를 이용하여 높은 토크 밀도를 가지는 것에 반해, 낮은 역률 특성을 가지므로 인버터 용량 및 비용이 증가하는 단점이 존재한다[6].

기존 PMVM 연구에서는 다양한 극-슬롯 조합을 대상으로 토크 및 역률 특성에 대한 분석이 수행되어 왔다. 그러나 기존 연구의 상당수는 고정자 슬롯수와 회전자 극수가 동시에 변화하는 비교 방식을 채택하고 있어, 기어비 변화 자체가 전동기 성능에 미치는 영향을 분리하여 평가하는 데 한계가 있다.

따라서 본 논문에서는 고정자 슬롯수를 12로 고정한 상태에서 회전자 극수를 8극, 16극, 20극으로 변화시켜 각각 기어비 1, 2, 5에 해당하는 모델을 설계하고, 기어비 변화가 토크 밀도 및 역률 특성에 미치는 영향을 분석한다. 또한 기어비 5 모델의 경우 코일 피치 1의 집중권과 코일 피치 2, 3의 분포권을 각각 적용할 수 있으므로, 권선법 및 코일 피치 변화에 따른 버니어 모터의 특성 변화를 추가적으로 비교 분석한다. 이때 권선법에 따른 엔드부 높이 차이를 반영하여 토크 밀도를 산정하였으며, 인휠 구동계의 비현가질량을 고려한 질량 기준 토크 밀도까지 함께 검토함으로써 설계 유효성을 평가한다. 이후 기어링 효과, 권선 계수, 축 방향 길이 등 주요 설계 파라미터 간의 상관관계를 2D 유한 요소 해석(Finite Element Method, FEM) 결과를 바탕으로 분석하고, 수식 기반의 파라미터 도출 과정을 통해 PMVM의 낮은 역률 특성의 원인을 고찰한다. 이를 통해 최종적으로 비교 대상 모델 중 토크 밀도와 역률 측면에서 전반적으로 우수한 특성을 갖는 모델을 도출한다.

2. 버니어 모터의 이론 및 특성

2.1 토크 발생 원리

PMSM은 고정자 자계와 회전자 자계가 동일한 극쌍수를 가지며, 고정자에서 형성되는 회전자계와 회전자 자계가 동기화되어 토크를 발생시킨다.

반면 PMVM은 마그네틱 기어와 유사한 자계 변조 효과를 이용하기 위해 고정자 자계와 회전자 자계가 서로 다른 극쌍수를 갖도록 설계된다[7]. 이때 회전자 자계의 극쌍수를 고정자 자계의 극쌍수로 나눈 값인 기어비는 식(1)로 정의된다.

(1)
$G = \frac{p_r}{p_s}$

회전자 영구자석 극쌍수, 고정자 권선 자계 극쌍수, 자계 변조 극(Flux Modulation Pole, FMP) 개수가 식(2)의 관계를 만족하면 고정자 권선 자계는 자계 변조 극과의 상호작용을 통해 변조되어 회전자 영구자석 자계의 극쌍수와 동일한 차수의 고조파를 포함하는 다양한 공간 고조파를 생성한다[8]. 이때 생성된 고조파 중 회전자 영구자석 자계와 동일한 극쌍수를 갖는 유효 고조파 성분은, 원래의 고정자 권선 자계와는 반대 방향의 회전 특성을 가지며, 회전자 영구자석 자계와 동기화되어 토크를 발생시킨다.

(2)
$p_r = N_{FMP} \pm p_s$

$p_r$은 회전자 영구자석 극쌍수, $p_s$는 고정자 권선자계 극쌍수, $N_{FMP}$는 자계 변조 극수이다.

2.2 공극 퍼미언스법

버니어 모터의 자계 변조 현상은 공극 퍼미언스법을 이용하여 해석할 수 있다. 공극 퍼미언스법은 회전자 영구자석에 의해 형성되는 기자력과 고정자 슬롯 및 자계 변조 극 구조에 의해 형성되는 공극 퍼미언스를 각각 함수 형태로 표현하고, 이들의 곱으로부터 공극 자속밀도를 구하는 방법이다.

먼저, 회전자 영구자석에 의해 형성되는 기자력의 기본파 성분은 식(3)과 같이 나타낼 수 있다.

(3)
$F_m(\theta, \theta_m) \approx F_1 \cos p_r(\theta - \theta_m)$

$\theta_m$은 회전자 초기 각, $\theta$는 회전자 회전 각, $F_1$은 영구자석 기자력 기본파 크기이다.

한편, 고정자 슬롯과 자계 변조 극 구조에 의해 형성되는 공극 퍼미언스는 식(4)와 같이 표현할 수 있다. $P_0, P_1$은 각각 공극 퍼미언스의 평균값, 고조파이다.

(4)
$P(\theta) \approx P_0 - P_1 \cos(N_{FMP}\theta)$

공극 자속밀도는 영구자석 기자력과 공극 퍼미언스의 곱으로부터 구할 수 있으며 식(5)와 같이 표현된다.

(5)
$B_{ver} \approx \sum \left[ F_1 P_0 \cos p_r(\theta - \theta_m) - \left(\frac{F_1 P_1}{2}\right) \cos((p_r - N_{FMP})\theta - p_r \theta_m) \right]$

(5)의 첫 번째 항은 회전자 영구자석에 의한 기본 자속밀도 성분을 의미하고, 두 번째 항은 자계 변조 극과의 상호작용에 의해 생성되는 변조 자속밀도 성분을 의미한다.

이때 변조 자속밀도 성분은 고정자 권선과 쇄교하여 역기전력 및 토크를 발생시킨다. 상 권선에 쇄교되는 자속으로부터 역기전력은 시간에 대한 미분으로 구할 수 있으며, 이를 정리하면 버니어 모터의 상 역기전력은 식(6)과 같이 표현된다.

(6)
$E_{ph} = \frac{K_w N_{ph} D_g l_{stk} \omega}{\sqrt{2}} \left\{ F_{pm} \left[ P_0 + \frac{6q-1}{2} P_1 \right] \right\}$

$K_w$는 권선 계수, $l_{stk}$는 축 방향 길이, $N_{ph}$는 턴수, $\omega$는 전기 각속도, $D_g$는 공극 길이, $q$는 매극매상 슬롯수를 나타낸다.

(6)으로부터 버니어 모터의 역기전력은 공극 퍼미언스의 평균 성분 $P_0$뿐만 아니라 자계 변조에 의해 생성되는 고조파 성분 $P_1$의 영향을 함께 받음을 알 수 있다. 특히 (6q−1)항은 버니어 구조에서 자계 변조 효과가 강화될수록 역기전력이 증가할 수 있음을 나타내고, 토크는 상 역기전력과 전류의 관계로부터 식(7)과 같이 나타낼 수 있다[7] [9].

(7)
$T = \frac{3 E_{ph} I_{max}}{2 \omega}$

(7)로부터 역기전력이 증가할수록 동일 전류 조건에서 출력 토크 또한 증가함을 알 수 있다. 따라서 버니어 모터는 자계 변조 극에 의해 형성되는 변조 자속밀도 성분을 이용하여 일반적인 PMSM보다 높은 역기전력과 토크를 얻을 수 있으며, 이것이 버니어 모터가 높은 토크 밀도를 가지는 주요 원인이다. 본 논문에서 단위 축 방향 길이 당 토크를 식(8)과 같이 정의하여 모델 간 토크 성능을 비교하였다.

(8)
$T_{unit} = \frac{T_{average}}{l_{stk}}$

2.3 페이저도 및 역률

PMVM의 역률 특성을 분석하기 위해 먼저 영구자석 전동기의 페이저도를 살펴본다. 페이저도는 적용되는 제어 방식에 따라 두 가지 형태로 구분할 수 있다. 그림 1의 (a)는 릴럭턴스 토크를 고려하지 않는 경우의 페이저도이며, 그림 1의 (b)는 고려한 경우이다[10].

이론적으로 표면부착형 PMVM 또한 id=0 제어를 적용하지만, 자속 포화에 의한 d-q축 인덕턴스 차이에 의해 최대 토크를 발생시키는 전류 위상각 $\beta$가 0이 아니므로, PMVM의 역률 분석을 위한 페이저도는 그림 1의 (b)가 적합하다.

그림 1. 영구자석 전동기 페이저도

Fig. 1. Permanent Magnet Motor Phasor Diagram

../../Resources/kiee/KIEE.2026.75.8.1773/fig1.png

영구자석 전동기의 d-q축 전압 방정식은 식(9)와 식(10)으로 표현되고, 전류는 식(11), (12)와 같이 분해되며, 역률각 $\theta$와 역률은 식(13),(14)와 같이 정의된다. 따라서 역률은 전압과 전류 사이의 위상각에 의해 결정되며, 일반적으로 q축 리액턴스에 의한 전압 강하가 증가할수록 감소한다.

(9)
$V \cos\delta = E - I_d X_d + I_q R_a$
(10)
$V \sin\delta = I_q X_q + I_d R_a$
(11)
$I_d = I_s \cos\beta$
(12)
$I_q = I_s \sin\beta$
(13)
$\theta = \delta - \beta$
(14)
$pf = \cos\theta$

이러한 관점에서 PMVM은 PMSM 대비 회전자 영구자석 극수가 많고 기어비가 크므로 동일 기계 각속도 조건에서 높은 전기 각속도를 갖는다. 따라서 권선 리액턴스 성분이 증가해 전압과 전류 사이의 위상차가 커져 역률이 감소한다. 또한 높은 회전자 극쌍수로 인접 자석 간의 누설 자속이 증가하며, 고정자 코일 피치에 비해 회전자 극 피치가 작은 구조에서는 누설 자속이 추가로 증가해 역률 저하의 원인이 된다[11].

또한 PMVM에서는 자계 변조 극에 의해 고정자 권선 기자력이 여러 차수의 고조파로 변조되며, 이들 중 회전자 영구자석 기자력과 동일한 차수를 갖는 고조파만이 평균 토크 생성에 기여한다. 따라서 나머지 고조파 성분은 평균 토크 생성에 기여하지 않으므로 누설 리액턴스와 누설 자속을 증가시켜 역률 저하의 원인으로 작용한다[12].

즉, PMVM의 역률 저하는 단순히 극수 증가에 따른 전기 각속도 증가뿐 아니라, 자계 변조에 의해 발생하는 비유효 고조파 성분에 따른 누설 리액턴스 증가, 그리고 회전자 극수 증가에 따른 자석 간 누설 자속 증가가 복합적으로 작용한 결과이다. 따라서 PMVM의 역률 특성을 정확히 분석하기 위해서는 기어비, 극-슬롯 조합, 권선 구조에 따른 리액턴스 변화와 고조파 영향을 함께 고려할 필요가 있다.

(15)
$X_q = \omega L_q$
(16)
$\omega = 2\pi f = 2\pi \frac{2p_r \times RPM}{120}$

(15), (16)에 따르면 Xq는 q축 인덕턴스 Lq뿐 아니라 전기 각속도 $\omega$에도 비례한다. 따라서 PMVM은 SPMSM에 비해 높은 회전자 극수를 가지므로 동일한 기계 각속도 조건에서도 더 큰 전기 각속도를 갖게 되며, 그 결과 Xq가 증가하는 경향을 나타낸다. 또한 PMVM은 높은 회전자 극수에 의해 인접 자석 간 누설 자속이 증가하고, 자계 변조로 인해 발생하는 비유효 고조파 성분 역시 리액턴스 성분을 증가시키므로 전반적으로 역률이 낮아지는 특성을 갖는다.

2.4 권선 계수

권선 계수는 공극에서 기자력을 결정하는 요소로 역기전력 기본파에 비례한다. 일반적인 영구자석 동기기의 경우의 권선 계수 Kw는 분포 계수 Kd와 단절 계수 Kp의 곱을 이용하여 구할 수 있으며 아래의 식은 각 계수의 관계 식이다[13]. $m$은 상의 개수, $\beta$는 극 피치와 코일 피치의 비율이다.

(17)
$K_d = \frac{\sin(\pi/2m)}{q \sin(\pi/2mq)}$
(18)
$K_p = \sin\frac{\beta\pi}{2}$
(19)
$K_w = K_d \times K_p$

일반적인 영구자석 동기기의 경우 회전자 영구자석 기자력이 역기전력에 영향을 미치지만, 버니어 모터의 경우 회전자 영구자석에 의한 기자력과 자계 변조에 의한 공간 고조파 두 성분에 의해 역기전력이 결정된다. 따라서 PMVM의 역기전력 수식을 이용하여 권선 계수를 유도할 수 있다[14].

(20)
$E_{ph}(t) = \omega D_g N_{ph} l_{stk} \left( B_{PM0} K_{w0} \mp \frac{p_s}{p_r} B_{PM1} K_w \right) \times \cos(p_r \theta_m \mp p_s \theta_{ph}) \\ = \omega D_g N_{ph} l_{stk} F_{pm} P_{0max} K_w \cos(p_r \theta_m \mp p_s \theta_{ph}) \times \cos(p_r \theta_m \mp p_s \theta_{ph})$

여기서 P0max는 형상에 따른 퍼미언스 최댓값이고, 최종적으로 버니어 전동기의 권선 계수는 식(21)과 같이 도출된다.

(21)
$K_w = \frac{P_0}{P_{0max}} K_{w0} + \frac{G_m P_1}{2 P_{0max}} K_{w1}$
(22)
$P_0 = \frac{\mu_0}{g_{m+a}} (1 - 1.6\beta c_0)$
(23)
$P_1 = \frac{\mu_0}{g_{m+a}} \frac{2\beta}{\pi} \left( \frac{0.78125}{0.78125 - 2c_0^2} \right) \sin(1.6\pi c_0)$
(24)
$\beta = \frac{1}{2} - \left( \sqrt{4 + \left( \frac{o}{g_{m+a}} \right)^2} \right)^{-1}$
(25)
$c_0 = \frac{o}{t_d}$

표 1. 모델별 권선 계수

Table 1. Winding Factors for Each Model

Model (a) (b) (c) (d) (e)
Kw 0.87 0.90 0.79 1.37 1.58

Kw0는 회전자 극수가 pr이고 권선은 ps를 기준으로 설계하였을 때 권선 계수, Kw1은 회전자 극수가 ps일 때 계산한 권선 계수이다. 이때 퍼미언스값은 식(22)-(25)를 통해 구할 수 있으며, 모델별 권선 계수를 표 1에 정리하였다. 이때 $t_d$는 슬롯 피치, $o$은 슬롯 개구폭, $g_{m+a}$는 유효 공극 길이이다.

3. 모델링 및 전자계 특성 분석

3.1 분석 모델 선정

PMVM의 전자계 특성을 분석하기 위해 표 2의 모델 설계 제원을 만족하는 분석 모델을 선정하였다. 2.2절에서 언급한 바와 같이 기어링 효과에 의해 기자력이 기존보다 높거나 낮은 차수의 고조파로 변조되어, 목적에 따라 PMVM을 저속 고토크 또는 고속 저토크로 구동할 수 있게 한다. 직접 구동 인휠 모터는 감속기 없이 작동하므로 저속 고토크 구동이 요구되며, 이를 위해 G>1을 만족하는 극-슬롯 조합을 선정해야 한다.

본 연구에서는 모델 간 비교 분석을 위해 회전자 단면 자석 사용량을 동일하게 설정하였다. 또한 설계 가능한 극-슬롯 조합 중 회전자 극수가 과도하게 커지지 않도록 고정자 슬롯수를 12로 고정하여 모델별 극-슬롯 조합을 선정하였다.

8극 12슬롯 SPMSM (a)모델을 기준으로, 일반적으로 PMVM에 적용되는 4:3 극-슬롯 조합을 적용하여 고정자 권선 자계의 극쌍수와 회전자 영구자석의 극쌍수의 비가 1:2가 되도록 16극 12슬롯 (b)모델을 설계하였다. 또한, 고정자 극쌍수와 회전자 극쌍수의 비가 1:5가 되도록 20극 12슬롯 PMVM을 설계하였다. 20극 12슬롯(G=5) 모델의 경우, 코일 피치를 1, 2, 3으로 다르게 설정할 수 있어 코일 피치 변화에 따른 특성 분석을 위해 코일 피치에 따라 (c), (d), (e)모델을 설계하였고 전체 모델의 형상을 그림 2에 정리하였다.

표 2. 모델 설계 제원

Table 2. Model Design Parameters

Motor Type SPMSM PMVM PMVM
Pole/Slot 8P12S 16P12S 20P12S
Gear Ratio - 2 5
Rotor O.D. 330 [mm]
Air-gap Length 1 [mm]
Fill Factor 45 [%]
Current Density 10 [Arms/mm2]
Base Speed 500 [rpm]
Average Torque 955 [N·m]
Rated Voltage 630 [Vll_peak]
PM Usage Area 6443.45 [mm2]

그림 2. SPMSM, PMVM 모델 형상

Fig. 2. Configurations of the SPMSM, PMVM Models

../../Resources/kiee/KIEE.2026.75.8.1773/fig2.png

표 3. 모델별 q축 파라미터

Table 3. q-Axis Parameters of Each Model

Model (a) (b) (c) (d) (e)
I[A] 5927
β[°] 26 12 6 4 4
Lq[μH] 5.04 2.49 2.13 3.98 4.94
Xq[mΩ] 1.0 1.0 1.1 2.1 2.6
XqIq[V] 5.3 5.9 6.6 12.3 15.3

설계 모델의 2D FEM 해석을 통해 도출한 파라미터를 표 3에 정리하였다. 기어비가 높은 모델일수록 Xq가 증가하고, Xq에 Iq를 곱한 q축 리액턴스 전압 강하 성분인 XqIq 또한 증가하는 것을 확인할 수 있다. 이는 Lq를 통해 비교할 경우 그 경향이 정확히 확인되지 않으나, 회전자 극수 증가에 따른 전기 각속도의 증가가 Xq 증가의 주요 원인임을 알 수 있다.

추가로 G=5 모델 내부에서 코일 피치가 1에서 3으로 증가할수록 Xq 및 XqIq가 증가하는 것을 통해 권선법에 따라 리액턴스 성분에 영향을 줌을 알 수 있고 이는 2.3절에서 언급한 고정자 코일 피치에 비해 회전자 극 피치가 작은 경우 누설 자속이 증가하여 리액턴스 성분이 증가한 것으로 볼 수 있다.

3.2 토크 특성 분석

선정 모델에 대한 자속 밀도 Contour Plot을 그림 3에 정리하였다. 앞서 2.2절에서 언급했듯이 PMVM은 기어링 효과에 의해 형성된 자계 변조 성분이 토크 생성에 기여하므로, 동일한 슬롯수, 단면 자석 사용량, 전류 밀도 조건에서도 SPMSM 대비 더 높은 토크를 발생시킬 수 있고, 이를 그림 4표4에서 확인할 수 있다.

특히 (a)모델인 SPMSM과 (b)모델인 G=2 PMVM은 권선 계수가 각각 0.87과 0.90으로 유사한 수준임에도 (b)모델의 Tunit이 더 크게 나타났다. 이는 단순한 권선 계수 차이에 의한 결과라기보다, PMVM의 기어링 효과에 의해 고정자 권선 기자력이 자계 변조 극과의 상호작용을 통해 회전자 영구자석 자계 차수에 대응하는 유효 고조파 성분을 형성하고, 이 성분이 평균 토크 생성에 기여하기 때문으로 해석할 수 있다. 이는 권선 계수가 거의 동일한 조건에서 기어링 효과가 평균 토크 증가에 미치는 영향을 직접적으로 보여주는 결과이다.

또한 (c)모델은 (a)모델보다 권선 계수가 더 낮음에도 Tunit이 더 크게 나타났다. 이는 권선 계수 감소에 따른 기자력 감소 효과보다 PMVM의 기어링 효과에 의한 토크 증가 효과가 더 크게 작용하였기 때문으로 해석할 수 있다. 즉, (c)모델의 결과는 PMVM에서 평균 토크가 단순히 권선 계수의 크기에 의해서만 결정되는 것이 아니라, 자계 변조에 의해 형성되는 유효 고조파 성분의 영향도 함께 반영된다는 점을 보여준다.

G=5 모델의 경우에는 코일 피치를 다르게 설정할 수 있으므로, 권선법 변화에 따른 전자계 특성 차이를 추가적으로 비교하였다. 코일 피치가 증가하면 권선 계수가 증가하게 되며, 이에 따라 공극에서 형성되는 기자력의 크기 또한 증가한다. 본 연구에서는 식(17)-(25)를 이용하여 각 권선법에 대한 권선 계수를 계산하였으며, 이를 바탕으로 코일 피치 변화가 평균 토크에 미치는 영향을 비교하였다. 동일한 G=5 조건에서 모델의 코일 피치 증가에 따라 권선 계수가 증가하고, 이에 따라 $T_{unit}$ 또한 증가하는 경향을 나타냈다.

그림 3. 모델별 자속 밀도 분포도

Fig. 3. Flux Density Contour Plot of Each Model

../../Resources/kiee/KIEE.2026.75.8.1773/fig3.png

그림 4. 모델별 Tunit 파형

Fig. 4. Tunit Waveforms of Each Model

../../Resources/kiee/KIEE.2026.75.8.1773/fig4.png

표 4. 2D FEM 해석 결과

Table 4. 2D FEM Analysis Results

Model (a) (b) (c) (d) (e)
Tunit[N·m] 6.3 10.0 8.1 10.7 11.0
Tripple[%] 19.3 9.7 7.1 5.7 6.1
THD[%] 15.9 9.0 7.0 7.2 7.5
lstk[mm] 151.8 95.2 117.9 89.1 86.8

부피 대비 토크 밀도는 권선 엔드부를 반영한 전체 축 방향 길이를 기준으로 산정하였다. (a)-(c)모델은 집중권을 적용해 코일 엔드부 높이를 40mm로 (d), (e)모델은 분포권을 적용해 70mm로 반영하였다.

표 5에 따르면 (b)모델은 (a)모델 대비 평균 토크 증가와 축 방향 길이 감소가 동시에 나타나 토크 밀도가 크게 향상되었다. (c)모델도 (a)모델 대비 높은 토크 밀도를 보였다. (d)모델은 분포권에 따른 엔드부 길이 증가로 (c)모델보다 낮은 토크 밀도를, (e)모델은 권선 계수 증가에 따른 평균 토크 향상이 엔드부 증가에 따른 부피 증가를 일부 상쇄해 (c)모델보다 소폭 높은 토크 밀도를 나타냈다. 이는 G=5 내에서 코일 피치 증가에 따른 권선 계수 향상이 토크 밀도 증가로 직결되는 것은 아니며, 코일 엔드부 높이를 함께 고려해야 함을 의미한다.

질량 대비 토크 밀도는 엔드부 길이에 따른 구리 질량 변화를 반영해 산정하고, 재료별 무게와 함께 표 6에 정리하였다. (b)모델은 (a)모델 대비 자석 및 강판 사용량이 감소하면서도 높은 평균 토크로 토크 밀도가 크게 향상되었다. (c)-(e)모델도 SPMSM 대비 높은 토크 밀도를 보였으나, 최대 평균 토크 모델인 (e)모델이 최대 토크 밀도를 갖지 않았다. 질량 증가까지 함께 고려한 경우 (d)모델이 가장 높은 토크 밀도를 가진다. 따라서 인휠 모터 설계 시 단순 평균 토크뿐 아니라 엔드부에 따른 체적 및 권선 질량 증가를 함께 고려해야 한다.

표 5. 모델별 부피 및 토크 밀도

Table 5. Volume and Torque Density of Each Model

Model (a) (b) (c) (d) (e)
Volume[m3] 0.0164 0.0116 0.0135 0.0136 0.0134
TDensity[kN·m/m3] 58.2 82.6 70.7 70.2 71.2

표 6. 모델별 재료 무게 및 토크 밀도

Table 6. Material Mass and Torque Density of Each Model

Model (a) (b) (c) (d) (e)
Copper[Kg] 8.1 6.6 7.3 8.4 10.3
Magnet[Kg] 51.9 32.6 40.3 30.5 29.7
Steel[Kg] 7.3 4.6 5.7 4.3 4.2
Total[Kg] 67.3 43.8 53.3 43.2 44.2
TDensity[N·m/Kg] 14.2 21.8 17.9 22.1 21.6

3.3 역률 특성 분석

앞서 2.3절에서 설명한 바와 같이 영구자석 전동기의 역률은 무부하 역기전력 E와 q축 리액턴스 전압 강하 성분인 XqIq의 상대적인 크기에 결정된다. 그림 1의 페이저도에서 확인할 수 있듯 XqIq 성분이 작을수록, 역기전력이 클수록 전압과 전류 사이의 위상차가 감소하여 역률이 개선된다. 따라서 PMVM의 역률 특성을 분석하기 위해서는 리액턴스 성분과 역기전력의 변화를 함께 검토할 필요가 있다.

모터 설계 시 목표 토크를 만족하는 적층 길이를 결정한 후, 전압 제한을 고려하여 턴수를 결정한다. 이때 전압은 무부하 역기전력과 저항 및 인덕턴스에 의한 전압 강하 성분의 합으로 해석할 수 있다. 식(6)에 따르면 무부하 역기전력은 회전자 극쌍수, 권선 계수, 축 방향 길이에 비례한다. 전압 강하 성분은 인덕턴스 및 저항에 비례하며, 저항은 축 방향 길이와 코일 피치에 의해 결정되고, 인덕턴스에 의한 전압 강하는 n2과 각속도에 비례한다. 따라서 동일 목표 토크와 전압 제한 조건에서 회전자 극쌍수가 클수록 무부하 역기전력이 커지는 만큼 턴수가 제한되며, 동시에 인덕턴스에 의한 전압 강하가 증가하여 서로 상반된 방향으로 작용한다. 이러한 상충관계로 인해 PMVM에서는 두 성분을 함께 고려한 설계가 요구된다.

그림 5. 무부하 역기전력 파형

Fig. 5. No-Load Back-EMF Waveforms of Each Model

../../Resources/kiee/KIEE.2026.75.8.1773/fig5.png

표 7. 선간 역기전력 및 전압 관련 파라미터

Table 7. Line-to-Line Back-EMF and Voltage Parameters

Model (a) (b) (c) (d) (e)
lstk[mm] 151.8 95.2 117.9 89.1 86.8
Ell_1[V] 11.1 11.5 11.6 14.4 15.8
Vll_1[V] 12.4 13.4 14.7 23.3 28.0
n[Turn] 51 47 42.75 27 22.5
Ell[V] 565.6 542.4 496.3 389.1 354.6
Vll[V] 630

표 8. 설계 모델 전자계 파라미터

Table 8. Model Parameters

Model (a) (b) (c) (d) (e)
Eph[V] 327.7 314.5 287.8 225.6 205.6
Is[A] 116.2 126.0 138.6 219.4 263.3
β[°] 26 12 6 4 4
R[mΩ] 209.4 159.3 145.5 67.5 70.6
Ld[mH] 9.3 5.1 3.8 2.9 2.5
Lq[mH] 13.1 5.5 3.9 2.9 2.5
Xq[Ω] 2.7 2.3 2.0 1.5 1.3
XqIq[V] 268.8 275.2 281.5 332.3 343.8

표 7에서 확인할 수 있듯이, 모델 (a)에서 (e)로 갈수록 1턴 기준 선간 역기전력 Ell_1이 증가하고, 표 3에서 확인할 수 있듯 XqIq가 증가하여 Vll_1또한 증가하게 되어 전압 제한 조건에 의해 결정되는 턴수는 감소한다. 그 결과 최종 선간 무부하 역기전력 Ell은 오히려 감소하는 경향을 보인다. 이는 PMVM 모델의 1턴 역기전력이 증가하더라도, q축 인덕턴스에 의한 전압 증가로 인해 전압 제한을 만족하는 턴수가 감소하여 최종 역기전력은 작아질 수 있음을 의미한다.

표 8은 정격속도 500RPM에서 설계 모델의 전자계 파라미터를 정리한 것이다. Iq는 턴수 n에 반비례하고, Lq는 n2에 비례하므로 Xq 역시 n2에 비례한다. 따라서 XqIq는 n에 비례하기 때문에 PMVM은 높은 회전자 극수에 의해 표 3와 같이 동일 턴수 기준에서는 Xq가 증가하더라도, 실제 전압 제한 조건에서 턴수가 감소하면 Xq는 오히려 감소할 수 있다.

그러나 표 8에서 확인할 수 있듯, Iq의 증가와 Xq의 감소가 동시에 작용한 결과, 최종적으로 XqIq는 모델 (a)에서 (e)로 갈수록 증가하는 경향을 나타낸다.

표 9그림 6은 앞서 구한 파라미터를 기반으로 모델 (a)-(e)의 역률, 피상전력, 페이저도를 나타낸 것이다. (a)모델이 가장 높은 역률을 가지며, 기어비 증가 및 코일 피치 증가에 따라 PMVM 모델의 역률이 점차 감소하는 경향을 확인할 수 있다. 특히 (d), (e)모델에서는 권선 계수 증가로 평균 토크 특성은 향상되었지만, 동시에 XqIq 성분의 상대적 영향이 커져 역률은 더욱 저하되는 결과를 보였다.

그림 6. 모델별 페이저도

Fig. 6. Phasor Diagrams for Each Model

../../Resources/kiee/KIEE.2026.75.8.1773/fig6.png

표 9. 모델별 역률 및 피상전력

Table 9. Power Factor and Apparent Power of Each Model

Model (a) (b) (c) (d) (e)
PF[%] 93.8 85.4 77.7 60.4 56.3
S[kVA] 59.8 64.3 71.4 93.2 102.4

표 10. 비유효 고조파 비율 및 자속 활용도 지표

Table 10. Non-working Harmonic Ratio and Flux Utilization Index

Model (a) (b) (c) (d) (e)
$K_v$ 1.68 1.47 1.64 1.39 1.29
$K_\phi$ 2.10 1.68 1.25 1.18 1.15

추가로 비유효 고조파 및 인접 자석 간 누설에 의한 역률 저감 분석을 위해 부하 공극 자속 밀도 고조파를 유효 고조파와 비유효 고조파로 구분하고, 비유효 고조파 대비 유효 고조파의 상대 크기 $K_v$를 식(26)과 같이 정의하였다. $B_v$의 차수는 모델별로 (a) 4차, (b) 8차, (c)∼(e) 10차로 설정하였다.

(26)
$K_v = \frac{B_v}{\sqrt{\sum B_{v,non}^2}}$

$B_v$ : 유효 고조파 공극 자속 밀도

$B_{v,non}$ : 비유효 고조파 공극 자속 밀도

3.4 손실 특성 분석

또한 동일 자석 사용 조건에서 인접 자석 간 누설을 간접적으로 비교하기 위해 1턴 기준 선간 역기전력 Ell_1을 이용하여 PM 자속 활용 계수 $K_\phi$를 식(27)과 같이 정의하였다.

(27)
$K_\phi = \frac{E_{ll\_1}}{p_r K_w l_{stk}} \times 100$

표 10과 같이, $K_v$와 $K_\phi$ 모두 (a)에서 (e)로 갈수록 감소한다. 이는 기어링 효과 및 회전자 극수 증가에 따라 공극 자속 밀도 고조파 중 유효 고조파 비율이 줄고, 인접 자석 간 누설로 유효 PM 자속 활용도가 저하되기 때문이다.

역률 저하는 시스템 차원의 인버터 용량 증가로 이어진다. 모델 (a)에서 (e)로 갈수록 역률이 감소함에 따라 요구 인버터 용량은 증가한다. 특히 G=5 모델의 경우 평균 토크 측면에서 유리한 결과를 보였으나, 표 9에 따르면 낮은 역률로 피상전력이 증가하여 인버터 용량 증가로 인해 시스템 측면에서 불리할 수 있다. 또한 동일 출력, 동일 전압 제한 조건 내에서 낮은 역률 특성은 출력 보상을 위한 스택 또는 입력 전류의 증가로 이어져 부피 증가와 동손 증가를 초래한다. 따라서 직접 구동 인휠 모터 설계에서는 평균 토크 향상뿐 아니라 역률 특성까지 함께 고려하여 최적 모델을 선정할 필요가 있다.

표 11은 정격속도 500RPM에서 모델별 손실 성분과 효율을 나타낸다. Wc는 동손, Wi는 철손, Weddy는 자석 와류손, Weddy/lstk는 적층 길이당 자석 와류손, η는 효율을 의미한다.

자석 와류손은 자석을 축 방향으로 5층 분할하고 각 층을 절연한 조건에서 3D FEM 해석을 진행하였다. 표 11에 따르면 모든 모델에서 동손이 가장 큰 비중을 차지하며, 철손과 자석 와류손은 상대적으로 작은 값을 나타낸다.

(b)모델(G=2 PMVM)은 (a)모델(SPMSM) 대비 평균 토크가 증가하면서도 동손이 감소하였다. 이는 기어링 효과로 목표 토크에 필요한 축 방향 길이가 줄어 전체 권선 길이가 감소했기 때문이다. 또한 자석 및 철심 사용량 감소에 따라 자석 와류손과 철손도 낮은 값을 가져 (b)모델이 정격속도에서 효율이 가장 높은 모델로 확인되었다.

표 12는 정격속도에서 부하 시 주요 공극 자속 밀도 Bg를 크기순으로 정리한 것으로, (b)모델은 (a)모델보다 높은 차수의 고조파를 포함한다. 이 고조파는 식(2)에 따라 고정자 권선 자계와 회전자 영구자석 자계가 자계 변조 극과 상호작용하며 형성된 성분이다. 이때 (b)모델의 손실이 크게 증가하지 않은 것은, 축 방향 길이 감소와 재료 사용량 감소 효과가 고조파 증가에 따른 손실 증가 요인을 상쇄하였기 때문이다.

(c)-(e)모델은 기어비 증가와 권선법 변화에 따라 손실 특성이 다르게 나타났다. (c)모델은 (a)모델 대비 높은 평균 토크를 확보하면서 집중권을 유지해 축 방향 길이 감소의 이점을 가지나, 회전자 극수 증가와 자계 변조로 고조파 자속 성분이 늘어 철손과 자석 와류손은 (a), (b)모델보다 증가한다.

즉 높은 기어비는 평균 토크 향상에는 유리하지만, 주요 공극 자속밀도 고조파 차수를 높여 자속 변화 주파수를 증가시키므로 철손과 자석 와류손에는 불리하다. 표 12에서 G=5 모델들이 SPMSM 및 G=2모델보다 높은 차수의 주요 고조파를 포함하는 것이 손실 증가의 원인으로 해석할 수 있다.

(d), (e)모델은 슬롯피치 2, 3의 분포권을 적용해 집중권 대비 엔드부 길이가 증가한다. 코일 피치 증가에 따른 권선 계수 향상으로 평균 토크가 증가하고 축 방향 길이는 감소하지만, 분포권에 따른 권선 길이 증가로 동손이 증가한다. 또한 높은 회전자 극수로 인해 공극 자속 밀도 주요 고조파 차수가 증가해 철손 및 자석 와류손 또한 증가하였다. 그 결과 (e)모델은 가장 높은 평균 토크에도 동손과 자석 와류손이 가장 크게 나타난다. 따라서 권선 계수 향상으로 평균 토크는 증가해도 손실 측면에서 반드시 유리하지는 않다.

Weddy/lstk를 제시한 이유는 모델별 적층 길이가 서로 다르기 때문이다. Weddy에는 축 방향 길이 차이와 고조파 영향이 동시에 반영되므로, 적층 길이당 자석 와류손 Weddy/lstk를 통해 자계 고조파와 회전자 극수 증가가 자석 손실에 미치는 영향을 직접 비교할 수 있다. 표 11에서 Weddy/lstk는 (a)모델에서 (e)모델로 갈수록 증가하며, 이는 기어비 증가 및 코일 피치 변화에 따라 자석 내부의 와전류 손실이 증가함을 의미한다.

표 11. 정격점 손실 특성

Table 11. Loss Components at Rated Speed

Model (a) (b) (c) (d) (e)
Wc[W] 5037.3 3793.8 4191.4 4870.4 5964.4
Wi[W] 89.1 197.1 273.2 300.7 341.6
Weddy[W] 964.2 914.0 1031.4 1106.5 1347.0
Weddy/lstk[W/m] 6.3 9.6 11.4 11.6 12.7
η[%] 93.0 94.8 93.7 94.4 94.0

표 12. 크기에 따른 주요 공극 자속 밀도 고조파 분석

Table 12. Major Air-Gap Flux Density Harmonics

Model (a) (b) (c) (d) (e)
1st Dominant Bg
(FFT order)
0.753
(4)
0.880
(8)
0.887
(10)
0.825
(10)
0.807
(10)
2nd Dominant Bg
(FFT order)
0.319
(8)
0.436
(4)
0.350
(2)
0.459
(2)
0.502
(2)
3rd Dominant Bg
(FFT order)
0.272
(12)
0.326
(24)
0.326
(30)
0.302
(30)
0.307
(30)
4th Dominant Bg
(FFT order)
0.162
(28)
0.249
(20)
0.253
(22)
0.223
(22)
0.218
(22)

그림 7. 모델별 효율맵 및 범례

Fig. 7. Efficiency Maps of Each Model and Legend

../../Resources/kiee/KIEE.2026.75.8.1773/fig7.png

그림 7의 효율맵을 보면, (b)모델은 비교적 넓은 운전 영역에서 높은 효율을 보인다. 이는 토크 특성 향상에 따른 적층 길이 감소와 상대적으로 낮은 손실 성분이 동시에 반영된 결과이다. (c)모델 역시 (a)모델 대비 우수한 효율 분포를 나타내지만, 기어비 증가에 따른 철손 및 자석 와류손 증가로 인해 G=2 모델만큼의 효율 향상은 나타나지 않는다. 한편 (d), (e)모델은 분포권 적용에 따라 저속·고토크 영역에서는 높은 효율을 유지하더라도, 권선 길이 증가와 고조파 성분 증가의 영향으로 일부 운전 영역에서 낮은 효율 특성을 보인다.

4. 결 론

본 논문에서는 직접 구동 인휠 모터의 최적 모델을 선정하기 위해 고정자 슬롯수 12, 회전자 자석 단면적, 전류밀도를 동일하게 고정하고, 8극 12슬롯 SPMSM과 16극 12슬롯, 20극 12슬롯 PMVM을 설계하여 전자계 특성을 비교하였다.

토크 특성에서 PMVM 모델이 SPMSM 모델 대비 높은 토크 특성을 보였다. 특히 (a)모델과 (b)모델은 권선 계수가 유사함에도 (b)모델의 Tunit이 더 크게 나타나, 기어링 효과가 평균 토크 향상에 직접 기여함을 확인하였다. G=5 모델에서는 코일 피치 증가에 따라 권선 계수가 증가하고 평균 토크 역시 증가하였으나, 엔드부 길이 증가를 고려하면 부피 및 질량 대비 토크 밀도 향상은 제한적이었다. 이는 PMVM이 동일 목표 토크 조건에서 더 짧은 축 방향 길이로 설계될 수 있어, 인휠 모터의 체적 저감과 토크 밀도 향상에 유리함을 보여준다.

역률 측면에서는 역기전력과 XqIq성분의 경향을 분석하였다. 1턴 및 n턴 기준 파라미터를 구분하여 분석한 결과, 기어비 및 코일 피치 증가는 동일 턴수에서 기준 역기전력을 증가시키나, 실제 전압 제한 조건에서는 역기전력 감소와 XqIq 증가가 동시에 나타나 PMVM의 역률을 저하시킨다. 특히 G=5 모델은 평균 토크 향상에도 낮은 역률로 요구되는 인버터 용량이 크게 증가하여 시스템 측면에서 불리할 수 있다.

최적 모델 선정을 위해 토크 밀도, 역률, 손실 및 인버터 용량을 종합적으로 검토하였다. (a)-(c)모델은 집중권, 코일 피치 2, 3인 (d), (e)모델은 분포권을 적용하였다. 분포권 모델은 권선 계수 향상에 따라 평균 토크가 증가하였으나, 엔드부 길이 증가에 따른 체적 및 권선 질량 증가, 그리고 손실 증가가 함께 나타났다. 또한 손실 특성 분석 결과, 모든 모델에서 동손이 지배적이었으며, G=5 모델은 주요 공극 자속밀도 고조파 차수 증가와 분포권 적용의 영향으로 철손과 자석 와류손이 증가하였다. 반면 G=2 PMVM은 평균 토크 향상과 축 방향 길이 감소가 동시에 나타나 정격점 효율이 가장 높고, 운전 영역 전반에서도 가장 균형적인 효율 특성을 보였다.

역률 측면에서 (b)모델은 PF = 0.854, (a)모델은 PF = 0.938로 역률이 약 0.084 감소하였으나 축 방향 길이를 56.6mm 줄이고 전체 부피를 약 29.3% 감소시킬 수 있었다. 이에 따라 (b)모델은 (a)모델 대비 더 높은 토크 밀도와 출력 밀도를 확보하면서 G=5 모델보다 우수한 역률과 낮은 인버터 용량 증가를 나타냈다. 따라서 본 연구에서 16극 12슬롯 PMVM이 직접 구동 인휠 모터에 가장 적합한 모델로 판단된다.

References

1 
A. Takahashi, T. Suto, M. Ito, R. Iwano, T. Hara, "Comparison of Gear-Drive and Direct-Drive Systems for In-Wheel Motors," Torino, Italy, pp. 1-7, 2024. Google Search
2 
S. P. Nikam, V. Rallabandi, B. G. Fernandes, "A High-Torque-Density Permanent-Magnet Free Motor for in-Wheel Electric Vehicle Application," IEEE Transactions on Industry Applications, vol. 48, no. 6, pp. 2287-2295, 2012. DOI
3 
S. Qiu, Y. Lu, J. Li, "Multi-Physics Design and Analysis of Axial Flux Permanent Magnet Motors for In-Wheel Direct-Drive," Torino, Italy, pp. 1-7, 2024. Google Search
4 
R. Wang, B. Wang, S. Dai, M. Cheng, W. Hua, "Design and Analysis of High Torque Density Dual Three-Phase Fault-Tolerant Permanent Magnet Vernier Machine for In-Wheel Drive," IEEE Transactions on Transportation Electrification, vol. 11, no. 6, pp. 13104-13114, 2025. DOI
5 
Y. Yu, S. Xie, F. Chai, Y. Pei, X. Zhang, C. H. T. Lee, "Analysis and Suppression of Core and PM Losses in Fractional-Slot Permanent Magnet Vernier Machine and PMSM for in-Wheel Application," IEEE Transactions on Industry Applications, vol. 61, no. 4, pp. 6305-6316, 2025. DOI
6 
Y. Zhao, X. Ren, X. Fan, D. Li, R. Qu, "A High Power Factor Permanent Magnet Vernier Machine With Modular Stator and Yokeless Rotor," IEEE Transactions on Industrial Electronics, vol. 70, no. 7, pp. 7141-7152, 2023. DOI
7 
N. Kim, C. Kim, T. Ji, H. Yoon, H. Oh, S. Jung, "Comparative Study of Electromagnetic Performances and Loss of Vernier Motor and Synchronous Motor According to Eddy Current Effect," The Transactions of the Korean Institute of Electrical Engineers, vol. 72, no. 1, pp. 57-64, 2023. DOI
8 
D. Jang, Ph.D. dissertation, "A Study on Analysis and Optimization of Magnetic Field Modulation in Permanent Magnet Vernier Motor," Dong-A University, Republic of Korea, 2019. Google Search
9 
B. Kim, T. A. Lipo, "Operation and Design Principles of a PM Vernier Motor," IEEE Transactions on Industry Applications, vol. 50, no. 6, pp. 3656-3663, 2014. DOI
10 
Y. Yu, F. Chai, Y. Pei, M. Doppelbauer, C. H. T. Lee, "Investigation of Power Factor Characteristic in Permanent Magnet Motors for In-Wheel Direct Drive Application," IEEE Transactions on Energy Conversion, vol. 39, no. 3, pp. 1699-1710, 2024. DOI
11 
Y. Zhao, D. Li, Z. Liang, R. Qu, "A High Power Factor PM Vernier Machine With Segmented Stator," IEEE Transactions on Transportation Electrification, vol. 10, no. 4, pp. 9294-9303, 2024. DOI
12 
D. K. Kana Padinharu, G. J. Li, Z. Q. Zhu, R. Clark, Z. Azar, A. Thomas, "Investigation of scaling effect on power factor of permanent magnet Vernier machines for wind power application," IET Electric Power Applications, vol. 14, no. 11, pp. 2136-2145, 2020. DOI
13 
H.-J. Kim, S.-W. Baek, "Comparison of Characteristics of Concentrated Winding Ferrite Permanent Magnet Synchronous Motors according to Pole/Slot Combinations," The Journal of the Korea Institute of Electronic Communication Sciences, vol. 20, no. 3, pp. 495-504, 2025. Google Search
14 
H.-S. Kwon, Ph.D. dissertation, "Optimal Design of Vernier Permanent Magnet Motor for Robot Joint," Seoul National University, Republic of Korea, 2020. Google Search

저자소개

명준 (Jun Myeong)
../../Resources/kiee/KIEE.2026.75.8.1773/au1.png

He received the B.S. degree in electrical engineering from Chonnam National University, Gwangju, South Korea, in 2025. He is currently pursuing the integrated M.S.- Ph.D. degree with the Department of Electrical Engineering, Chonnam National University, Gwangju, South Korea. His research interests include numerical analysis and design optimization of electric machines.

임재성 (Jae Seong Lim)
../../Resources/kiee/KIEE.2026.75.8.1773/au2.png

He received the B.S. degree in electrical engineering from Chonnam National University, Gwangju, South Korea, in 2025. He is currently pursuing the M.S. degree in the Department of Electrical Engineering, Chonnam National University, Gwangju, South Korea. His research interests include optimal design, electromagnetic analysis, and performance improvement of electric machines.

이현식 (Hyeon Sik Lee)
../../Resources/kiee/KIEE.2026.75.8.1773/au3.png

He received the B.S. degree in electrical engineering from Chonnam National University, Gwangju, South Korea, in 2025. He is currently pursuing the integrated M.S.- Ph.D. degree with the Department of Electrical Engineering, Chonnam National University, Gwangju, South Korea. His research interests include thermal analysis and design optimization of electric machines, with particular interest in axial-flux motors.

이환규 (Hwan Gyu Lee)
../../Resources/kiee/KIEE.2026.75.8.1773/au4.png

He received the B.S. degree in electrical engineering from Dongshin University, Gwangju, South Korea, in 2025. He is currently pursuing the M.S. degree with the Department of Electrical Engineering, Chonnam National University, Gwangju, South Korea. His research interests include the design, analysis, and optimization of various electric machine topologies for traction and high-torque applications.

이진환 (Jin Hwan Lee)
../../Resources/kiee/KIEE.2026.75.8.1773/au5.png

He received the B.S. and Ph.D. degrees in electronic and electrical engineering from Sungkyunkwan University, Suwon, South Korea, in 2013, and 2018, respectively. He was a Senior Research Engineer in the R&D Division at Hyundai-Transys. He is currently an assistant professor in department of electrical engineering at Chonnam National University. His research interests include numerical analysis and design optimization of electric machines, development of optimization algorithm customized for electric machines, and development on analysis precision improvement technique of electric machines.