풍력터빈 PM형 동기발전기의 와전류손실과 열 해석에 관한 연구

A Study on Rotor Eddy Current Loss and Thermal Analysis of PM Synchronous Generator for Wind Turbine

Abstract

In this paper, eddy current loss, iron loss and heat transfer of PMSG with 2,000kW capacities were analyzed for wind turbine. The PMSG with 3 split magnet was analyzed using ansoft maxwell commercial program and, generator was tested by Back to Back converter with no load condition at laboratory. Rotor surface temperature was measured by Pt100 sensors for investigating heat transfer from rotor to atmosphere. The simulation results shows 27.4kW eddy current loss in no load condition and 50.2kW eddy current loss in rated load condition with 3 split magnet, and also shows 4.3kW iron loss in no load condition and 7.3kW iron loss rated load condition. The heat transfer coefficient of convection between rotor surface and atmosphere was investigated by $9.6W/m^2{\cdot}K$. Therefore the heat transfer from rotor to atmosphere was about 17kW(54%) and from rotor to air-gap was about 14.6kW(46%) in no load condition. It is identified that the cooling system for stator have to include the 46% of iron loss, and heat dissipation structure of rotor surface have to be suggested and designed for efficiency improvement of generator.

1. 서 론

최근 풍력발전기 산업에서는 고장이 적고 간단한 Drive Train의 요구가 점점 커지고 있으며, 이에 직결 방식의 풍력발전기가 현재 활발히 연구되고 있다. 직결방식 중에서도 크기와 무게 측면에서 장점이 많은 영구자석형 동기발전기의 개발이 네덜란드Delft 공대, Harkosan, Vensys사 등을 중심으로 더욱 활발히 이루어져 왔으며, 최근 그 기술들이 완성 단계에 이르러 제품들을 출시하고 있다[1].

최근 영구자석 발전기 설계에서 회전자 부분에서 발생하는 손실에 관한 연구가 많이 이루어지고 있으며, 회전자 부분에서 발생하는 손실은 주로 회전자의 도전체에서 발생하는 와전류 손실을 의미한다.[2] 와전류 손실은 회전자의 도전체, 즉 영구자석이나 영구자석의 비산을 방지하기 위해 삽입된 리테이너 등에서 발생하는 손실을 의미한다. 이러한 와전류 손실은 모두 열의 소스원이 되며 이는 영구자석의 감자 현상에 영향을 줄 뿐만 아니라 권선의 저항 등에도 영향을 미치게 되어 발전기의 안정적인 운전을 방해하는 요소가 되고 있다[3∼5].

자속이 도전체를 통과할 때 자속의 변화에 상응하는 와전류가 생성 되는데, 그 도전체의 면적이 클수록 큰 전압 불균형으로 인해 큰 와전류가 생성된다. 그러므로 자속이 통과하는 도전체는 자속 통과면적을 작게 할 필요가 있다. 영구자석 또한 동일한 영향으로 와전류가 생성되며 이 손실을 줄여야 발전기 효율을 높일 수 있다. 그러나, 현실적으로 영구자석을 무한하게 분할할 수 없고, 설계되는 발전기의 제작 가능성과 제작비용 등을 고려하여 분할 수량을 결정하여야 한다.

본 연구에서는 설계된 발전기는 와전류손을 줄이기 위하여 영구자석을 최대 3개로 분할하였다. 한편 영구자석의 와전류 손실계산에서는 영구자석의 길이와 넓이의 비례관계에서 2차원 해석결과는 보정을 통하여 3차원 해석과 유사한 결과를 확인할 수 있다.

따라서, Ansoft maxwell 2D 프로그램을 이용해 영구자석의 분할에 따른 와전류 손실을 2차원 유한요소 해석을 통하여 해석결과를 확인하고 로터 요크의 철손을 확인하고자 한다[6]. 한편 발전기에서 생성되는 손실은 모두 열로 발산이 되며, 열의 발산방식은 대류, 전도, 복사 3가지의 형태를 가진다. 열의 발산이 잘되지 않고 열에너지가 축적되는 곳은 온도가 상승하여 재질이 열화 혹은 소손이 발생하므로 열 방출을 위한 구조 및 냉각시스템의 설계가 매우 중요하다[7].

본 연구에서는 자석의 와전류분석, Rotor Yoke의 철손 분석 및 연계된 열해석을 통해 적절한 발전기 냉각시스템 설계의 방향을 제시하고자 한다.

 

2. 동기발전기 설계 및 시뮬레이션

2.1 해석 모델

2,000kW 풍력터빈 표면부착형 영구자석 동기발전기의 설계 요구사항은 표 1과 같다. 설계출력과 회전속도는 각각 2,080kW, 17rpm 이며 부하 선전압은 392Vrms의 요건이다. 그림 1은 설계된 발전기로 영구자석은 1극당 3개로 분할한 형태이며, 동기발전기의 해석모델의 설계사양은 표 2와 같다.2.1 해석 모델

표 1풍력발전기 시스템 설계 사양 Table 1 Design specifications of WTG

그림 1영구자석 설계 Fig. 1 Design model of PMSG

표 2PMSG 발전기 설계사양 Table 2 Design specifications of PMSG model

본 논문에서 선정된 영구자석의 경우 다른 자석 보다 높은 에너지 밀도를 갖고 에너지 변환기기에 적용 시 높은 출력밀도를 낼 수 있는 NdFeB를 선정하였으며, NdFeB계의 영구자석 특징은 표 3과 같다.

표 3영구자석 재료 물성치 Table 3 Material Properties of PM

2.2 와전류 해석이론

일반적으로, 자계 분포 예측을 위한 해석은 자계 , 자속밀도 , 그리고 전류밀도 를 갖는 Maxwell 방정식으로부터 얻어진다. 자기 벡터 포텐셜 는 식 (1), (2)로 정의된다.

패러데이 법칙과 자기 벡터 포텐셜의 정의로부터 식(3)으로 표현되는 전계특성식을 얻을 수 있다.

도체에 유도되는 와전류는 로 표현 가능하므로 식 (3)은 식(5)로 다시 쓸 수 있다. E는 전계이고, σ 는 전기전도도이다.

2.3 와전류 시뮬레이션 및 결과

본 논문에서 Ansoft Maxwell 2D 프로그램을 활용하여 유한요소 해석 시뮬레이션을 수행하였다. 해석 조건은 무부하 운전일 때와 정격 운전 일 때 발생되는 와전류 손실을 분석하였다. 정격속도는 17rpm, 정격 주파수 13.88Hz, 정격 전류 2,300Arms이고 그림 2는 발전기 설계도면을 나타낸다.

그림 2발전기 설계된 모델 Fig. 2 Designed PMSG Model

와전류 발생 억제를 위해 3분할한 자석 모델의 해석 결과는 무부하 조건에서 와전류손실이 27.42kW이고, 정격부하에서는 50.21kW로 나타났다.

 

3. 발전기 Rotor Yoke의 철손

3.1 철손 해석이론

철손 Pv는 시간적으로 변하는 자계에 의해 철심내부에서의 자구가 변화할 때 발생되는 에너지 손실에 의한 히스테리시스 손실 Ph, 와전류 손실 Pc, 이상 와전류 손실분 Pe로 구분된다.

여기서 Ph는 히스테리시스손실, Pc는 와전 류 손실, Pe는 이상와전류 손실이며 Kh는 히스테리시스 손실상수, Kc는 와전류손실 상수, Ke는 이상와전류 손실상수를 말한다.

그림 3은 구상흑연주철로 설계된 로터 요크의 자기적 특성을 측정하기 위한 시편과 시험 장비를 나타내며, 측정한 결과 BH 곡선은 그림 4와 같고, 주파수에 따른 철손은 그림 5와 같이 나타난다.

그림 3B-H 실험 Fig. 3 B-H Test

그림 4로터 요크 B-H 곡선 Fig. 4 B-H Curve for Rotor Yoke

그림 5로터 요크의 주파수에 따른 철손 특성 Fig. 5 Iron loss characteristics according to frequency for Rotor Yoke

그림 4에서 로터 요크는 자계(H)에 따라 자속밀도(B)가 1.2T에서 급격히 포화되기 시작하므로 1.2T 이하가 되도록 설계하는 것이 필요하며, 그림 5는 사용 주파수에 따른 구상 흑연 주철의 손실 값을 측정한 것으로 주파수가 증가하면 손실이 증가한다.

3.2 철손 시뮬레이션 및 결과

Ansoft Maxwell 2D 프로그램을 활용하여 유한요소 해석 시뮬레이션을 수행하였으며, 해석 조건은 무부하 운전일 때와 정격 운전 일 때의 철손을 분석하였다. 이때의 정격속도는 17rpm, 정격 주파수 13.88Hz, 정격 전류 2,300Arms이다. 시뮬레이션 결과는 그림 6과 같다.

그림 6철손 해석 결과 Fig. 6 The analysis result of iron loss

무부하 조건에서 Rotor yoke에 생기는 철손은 4.2kW이며, 정격운전에서는 7.3kW이다. 그림 7은 정격운전에서 철손 분포도를 나타낸 것이다.

그림 7부하시 로터요크의 철손 분포 Fig. 7 Iron loss distribution of Rotor Yoke in load operation

 

4. 발전기 열발생

발전기에서 생성된 와전류 및 Rotor 철손은 열로 변환된다. 본 연구의 발전기는 그림 8과 같이 외전형으로 Rotor가 외측에서 회전하며 직접 대기와 접촉하여 방열하는 구조이고 내부의 Stator에 수냉식 냉각구조를 가지고 있다.

그림 8발전기의 작동부 Fig. 8 Active Part of Generator

회전자(로터) 냉각 효율의 관점에서 볼 때 종래의 대형 발전기가 내부 회전 로터 방식이여서 로터 냉각이 어렵고, 따라서 로터에 부착된 마그넷이 쉽게 고온이 되어 탈자화되는 단점을 방지하는 장점을 가지고 있다. 본 연구에서 외부에서 회전하는 로터의 열발생 근원인 영구자석 와전류를 포함한 철손이 열로 변환된 후 Rotor를 통해 대기로 방출되는 열과 발전기 내부로 전달되는 비율을 분석하였다. 특히, Stator 내부로 전달되는 열은 발전기의 효율과 수명을 단축 시키는 큰 영향을 주기 때문에 발전기 설계시 반드시 고려해야 할 중요한 요인이다. 그러나 발전기 내부의 공극을 통한 대류현상과 복사에 의한 방열량을 실험으로 계측하거나 이론적으로 추정하기가 어렵다.

본 연구에서 발전기의 방열상태를 해석하기 위해서 Maxwell 2D로 해석한 무부하시 이 정격속도의 와전류손 및 로터 철손이 모두 열로 변환된다고 가정한다. 실험을 통하여 무부하 정격속도로 발전기를 운전하여 로터 외측의 온도를 측정하고 외부로 전달되는 열전달을 분석한다, 이로 인해 Stator 내부로 전달되는 열전달을 추정할 수 있다. 제작된 MG(Motor-Generator)세트는 그림 9와 같이 2개의 발전기를 제작하여 하나는 모터로 다른 하나는 발전기로 동작시켜 시험하였다. 모터의 속도제어는 2MW 용량의 인버터를 사용하여 정격속도로 운전하였다.

그림 9모터-발전기 테스트 세트 Fig. 9 MG Test Set

4.1 열전달 해석이론

자석에서 발생한 열은 열전도율이 높은 로터 요크로의 많은 열전달이 예상되지만 공극 5mm를 통한 Stator 방향으로의 복사와 대류 열전달도 예상할 수 있다. 각 부분에 대한 열전달 계수는 표 4와 같다.

표 4열전달 계수 Table 4 The coefficient of heat transfer

전도에 의한 열전달 효율은 식 (7)에 의해 지배되고, 복사에 의한 열전달 효율은 식 (8)에 의해 지배된다.

여기서 α는 열전달 효율, K는 열전도도, T는 온도, y는 거리이다. ∊는 방사율이며 σ는 Stefan-Boltzmann 상수이다. 대류에 의한 열전달 효율은 식(9)에 의해 지배된다.

그리고 로터 요크의 표면 온도 상승시 외부로 방출되는 열전달은 대류현상으로 Newton의 냉각법칙 식(10)을 따른다.

여기의 는 단위시간당 방출 열량이고, A는 열방출 면적이다.

17rpm으로 회전하는 로터 표면의 대류 현상을 난류유동으로 적용하면, 표면에서의 온도구배를 의미하는 Nusselt수는

이다. 여기서 는 Reynolds 수이며, 는 Prandtl 수를 말한다.

U, X, 𝜐는 각각 기체의 속도, 방열부의 길이, 기체의 동점성(kinematic viscosity)을 나타낸다. 그리고 μ, Cp, α는 각각 역학적점성 (dynamic viscosity), 비열, 열확산계수를 나타낸다. 한편, Nu수는 열경계층 이론에 의해 다음과 같이 식(12)로 정의한다.

4.2 Rotor 온도측정

로터가 정격속도 17rpm으로 운전시 로터의 표면에는 온도 측정용 Pt100Ω 센서 2개를 부착하여 측정하였다. 그리고, 실험에 사용된 발전기의 Rotor 사양은 다음의 표 5과 같다.

표 5로터의 사양 Table 5 Specifications of Rotor

로터 표면은 방열 효율을 위해 요철 구조로 되어있으며, 마루부와 골부에 각각 한 개씩 총 2개의 Pt100 센서를 설치하여 약 3시간에 걸쳐 측정하였다. 외부 온도는 11.5℃이며, 로터 표면은 80℃까지 상승하였다. 온도측정 결과 2개의 센서 데이터는 유사한 경향을 보였고, 그 결과는 그림 10에 나타내고 있다.

그림 10로터 요크의 온도 변화 Fig. 10 Temperature of Rotor Yoke

공기의 물성치는 온도에 따라 다르므로 기준 온도는 막온도(film temperature)를 기준으로 다음과 같다.

막온도 이다. 이때의 공기 물성치는 표 6과 같다.

표 6공기의 열 특성 Table 6 Thermal Properties of Air

표 6에서 언급된 물성치를 적용하여 열전달에 적용된 Parameter 및 로터 표면을 통한 외부로의 대류 열전달은 표 7과 같다. 표 7에서 열전달 16.97 kW는 로터의 표면으로부터 외부 대기로 방출되는 열량을 나타내고 있다.

표 7열전달량 Table 7 Value of Heat Transfer

4.3 결과 분석

앞에서 언급한 Maxwell 2D로 해석한 무부하 및 정격부하에서의 영구자석의 와전류 손실과 Rotor Yoke의 철손은 표 8과 같이 해석되었다. 표 8에서 보면 정격부하시의 총손실은 57.5kW로 무부하시 대비 약 82%가 증가를 보여주고 있다. 그리고 로터의 손실중에는 자석의 와전류 손실이 약 87%로 많은 비중을 차지함을 알 수 있다.

표 8로터 손실 Table 8 Loss of Rotor

표 9와 같이 자석의 와전류 및 로터 철손에 의한 손실 중 약 54%가 로터의 표면을 통해 외부로 열방출되며, 약 46%가 공극을 통해 발전기 내부로 전달됨을 확인할 수 있었다. 또한, 로터 표면으로부터 방출되는 열전달 효율계수 αCV는 9.6 W/m2·K가 됨을 알 수 있었다.

표 9열전달 경로 Table 9 Path of Heat Transfer

 

5. 결 론

본 논문에서는 영구자석의 와전류 손실 및 로터의 철손을 2차원 유한 요소법을 이용하여 구하고, 실험적으로 로터의 온도를 측정함으로서 열전달의 경향을 분석하였다.

자계해석을 통한 로터의 총 손실이 무부하시 31.6kW, 정격부하시 57.5kW임을 확인할 수 있었으며, 정격부하시의 손실은 무부하시 대비 약 82% 증가하였고, 이는 발전기 정격출력 2,000kW의 약 2.9% 임을 확인 할 수 있었다. 그리고, 로터의 손실은 자석의 와전류 손실이 대부분으로 약 87% 차지함을 알 수 있었다. 따라서 로터의 손실을 줄이기 위해서는 자석의 와전류 손실을 줄이는 연구방향이 매우 효과적이라는 것을 알 수 있었다.

한편 로터에서 발생한 손실로 로터 표면을 통해 외부로 방출되는 열전달은 약 54% 이고, 46%의 열은 공극을 통해 발전기 내부 Stator로 전달됨을 확인하였다. 따라서 발전기 냉각시스템 설계시 Stator 냉각 용량 외에도 로터로부터 전달되는 46%의 열용량을 고려한 냉각시스템 설계가 되어야 함을 알 수 있다.

또한, 로터 표면으로부터 방출되는 열전달 효율계수 αCV가 9.6 W/m2·K가 됨을 확인 할 수 있었다. 따라서 로터 표면을 통한 열전달의 비율을 높이기 위해서, 방열 효율이 높은 구조로 로터 표면을 연구하여 열전달 효율계수 αCV를 향상 시킬 필요가 있다.

향후 연구는 자석의 와전류 손실을 낮추기 위해 추가적인 자석 분할에 대한 해석이 필요하며, 발전기 내부 온도 측정을 통한 열전달의 분석이 필요하다.