KR102465767B1 - 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템 - Google Patents

정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템 Download PDF

Info

Publication number
KR102465767B1
KR102465767B1 KR1020210070030A KR20210070030A KR102465767B1 KR 102465767 B1 KR102465767 B1 KR 102465767B1 KR 1020210070030 A KR1020210070030 A KR 1020210070030A KR 20210070030 A KR20210070030 A KR 20210070030A KR 102465767 B1 KR102465767 B1 KR 102465767B1
Authority
KR
South Korea
Prior art keywords
thrust
tower shadow
tower
equation
shadow
Prior art date
Application number
KR1020210070030A
Other languages
English (en)
Inventor
임채욱
Original Assignee
한밭대학교 산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 한밭대학교 산학협력단 filed Critical 한밭대학교 산학협력단
Priority to KR1020210070030A priority Critical patent/KR102465767B1/ko
Application granted granted Critical
Publication of KR102465767B1 publication Critical patent/KR102465767B1/ko

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2260/00—Function
    • F05B2260/84—Modelling or simulation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Wind Motors (AREA)

Abstract

본 발명은 본 발명은 타워 섀도우에 의한 풍력터빈의 추력 모델링 방법 및 시스템에 관한 것으로, 보다 상세하게는 타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값을 연산하는 추력연산부; 타워 섀도우에 의해 변동하는 추력값을 획득하는 타워 섀도우 추력변동획득부; 상기 추력연산부에서 연산된 추력값과, 상기 타워 섀도우 추력변동획득부에서 획득된 추력값을 기반으로 타워 섀도우 추력변동계수를 연산하는 타워 섀도우 추력변동계수 연산부; 상기 타워 섀도우 추력변동계수를 데이터베이스화하는 데이터베이스; 및 상기 추력연산부에서 연산된 추력값과, 상기 데이터베이스에 저장된 타워 섀도우 추력변동계수를 이용하여 타워 섀도우에 의해 각 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 블레이드 추력 모델링부;를 포함하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템에 관한 것이다.

Description

정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템{Evaluation of Dynamic Thrust under Tower Shadow in Wind Turbine below Rated Wind Speed}
본 발명은 타워 섀도우에 의한 풍력터빈의 추력 모델링 방법 및 시스템에 대한 것이다. 보다 상세하게는, 타워 섀도우 추력변동계수를 정의하여 이러한 타워 섀도우 추력변동계수를 이용하여 타워 섀도우에 의해 블레이드에 발생하는 추력변동 동하중을 모델링할 수 있는 방법 및 시스템에 관한 것이다.
바람은 유용하게 이용할 수 있는 재생에너지 중 하나이다. 풍력터빈은 발전기의 기술적인 발전과 블레이드를 이용한 양력기술과 역사를 같이하고 있다. 풍력터빈은 바람이 가진 유용한 에너지를 로터 블레이드의 회전을 통하여 기계적 에너지로 전환하고, 발전기 제어를 통하여 전기를 생산하는 기계이다. 최근 풍력터빈은 대형화되는 추세로 수 MW급이 시장의 주류를 형성하고 있다
도 1은 통상의 풍력터빈시스템(1)의 부분 사시도를 도시한 것이다. 도 1에 도시된 바와 같이, MW급 대형 풍력터빈시스템(1)은 블레이드(3), 허브(6), 로터(2), 기어박스(5), 타워(7), 베어링, 발전기(4), 피치시스템, 요잉시스템 등 다양한 부품들로 구성되어 있다.
풍력터빈에서는 발전기의 제어를 통하여 생산되는 전기 파워의 질을 높이는 것도 중요하지만, 각 부품에 작용하는 정적 및 동적 하중에 대하여 20년의 내구성을 갖도록 설계 및 제작하는 것 또한 중요하다.
풍력터빈이 대형화되면서 블레이드의 길이가 40~50미터 이상으로 길어지게 되어, 로터 블레이드가 회전할 때 블레이드에 발생하는 비대칭 하중이 증가하게 되었다. 윈드쉬어(Wind Shear), 타워새도우(Tower Shadow), 난류풍속 등의 영향으로 블레이드에는 비대칭 하중이 발생한다. 이 비대칭 하중이 커지면 동적 하중이 증가하여 풍력터빈의 설계 시에 이를 반영해야 하며, 최근에는 이 비대칭 하중을 줄이기 위하여 세 개의 블레이드의 피치각을 독립적으로 제어하는 개별피치제어 방법을 적용하고 있다.
타워 섀도우에 의해 블레이드에 작용하는 추력의 변동을 구하기 위해서는 양력계수와 항력계수와 같은 공력계수를 이용하여 블레이드의 길이 방향에 대하여 적분하는 BEM(Blade Element Momentum) 방법을 적용하면 된다. 블레이드에 작용하는 비대칭 동하중을 계산하기 위하여 지금까지의 연구들에서는 BEM 방법을 주로 적용하였다. Bladed와 같은 상용 소프트웨어에서도 BEM 방법을 적용하여 추력변동에 의한 비대칭 동하중을 계산하고 있다.
풍력터빈의 추력은 일반적으로 공력계수의 정보를 이용하여 BEM 방법을 적용하여 계산할 수 있지만, 공력계수 정보 대신에 추력계수를 이용하여도 구할 수 있다. 풍력터빈의 추력 계산 시에 추력계수를 이용하는 것이 실용적이고 용이하지만, 추력계수를 이용한 추력 표현식을 적용하면 허브에서의 풍속 변화에 따른 추력의 변화만 표현할 뿐, 타워 섀도우 등에 의해 블레이드에 작용하는 비대칭 동하중에 의한 추력 변동을 표현할 수 없다. 타워 섀도우에 의한 추력의 변동은 일반적으로 BEM 방법을 적용하여 계산할 수 있는데, 기존의 추력계수를 이용한 추력 표현식을 적용하여도 동일하게 표현할 수 있다면 추력 변동의 모델이 유용할 것이다. 더 나아가 추력계수를 이용하여 타워 섀도우의 영향에 의한 추력변동을 표현할 수 있다면, 향후 이를 이용하여 타워 섀도우를 추정하거나 타워제어와 개별피치제어에 적용하기 위한 기초 기술이 될 수 있을 것이다.
한국 공개특허 제2015-0126127호 일본 공개특허 제2010-281297호 한국 공개특허 제2014-0041223호 일본 공개특허 제2009-114977호
따라서 본 발명은 상기와 같은 종래의 문제점을 해결하기 위하여 안출된 것으로서, 본 발명의 실시예에 따르면, 타워 섀도우 추력변동계수를 정의하여 이러한 타워 섀도우 추력변동계수를 이용하여 정격풍속 이하에서, 타워 섀도우에 의해 블레이드에 발생하는 추력변동 동하중을 간편하고 신속하게 모델링할 수 있는, 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈 블레이드의 추력 동하중 모델링 방법 및 시스템을 제공하는데 그 목적이 있다.
또한, 본 발명의 일실시예에 따르면, 타워 섀도우에 대한 추력변동에 의한 동하중을 모델링하는데 있어, 공력계수를 이용한 BEM 방법이 아닌, 정의되는 타워 섀도우 추력변동계수와 추력계수를 이용하여 신속하고 효율적으로 타워 섀도우에 대한 추력변동을 모델링할 수 있는, 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈 블레이드의 추력 동하중 모델링 방법 및 시스템을 제공하는데 그 목적이 있다.
한편, 본 발명에서 이루고자 하는 기술적 과제들은 이상에서 언급한 기술적 과제들로 제한되지 않으며, 언급하지 않은 또 다른 기술적 과제들은 아래의 기재로부터 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.
본 발명의 제1목적은, 정격풍속 이하에서 풍력터빈의 추력 동하중 모델링 시스템에 있어서, 타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값을 연산하는 추력연산부; 타워 섀도우에 의해 변동하는 추력값을 획득하는 타워 섀도우 추력변동획득부; 상기 추력연산부에서 연산된 추력값과, 상기 타워 섀도우 추력변동획득부에서 획득된 추력값을 기반으로 타워 섀도우 추력변동계수를 연산하는 타워 섀도우 추력변동계수 연산부; 상기 타워 섀도우 추력변동계수를 데이터베이스화하는 데이터베이스; 및 상기 추력연산부에서 연산된 추력값과, 상기 데이터베이스에 저장된 타워 섀도우 추력변동계수를 이용하여 타워 섀도우에 의해 각 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 블레이드 추력 모델링부;를 포함하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템으로서 달성될 수 있다.
그리고 상기 블레이드 추력 모델링부에 의해 모델링된 각 블레이드의 동하중 모델링을 기반으로 로터 추력 동하중을 모델링하는 로터추력 모델링부;를 더 포함하는 것을 특징으로 할 수 있다.
또한 상기 타워 섀도우에 의한 풍속 변동은 이하의 수학식 1 및 수학식 2에 의해 정의되는 것을 특징으로 할 수 있다.
[수학식 1]
Figure 112021062585564-pat00001
[수학식 2]
Figure 112021062585564-pat00002
상기 수학식 1, 2에서, Vh는 허브 높이에서의 평균 풍속, vts는 타워 섀도우에 의해 변동하는 풍속, Dt는 타워의 지름, x는 타워의 중심부터 타워 앞에 있는 블레이드 중심까지의 거리, y는 풍속이 통과하는 로터의 회전면에서의 타워 중심부터 블레이드까지의 횡 방향 거리이다.
그리고 상기 추력연산부에 의해 연산되는 추력값은 이하의 수학식 5에 의해 연산되는 것을 특징으로 할 수 있다.
[수학식 5]
Figure 112021062585564-pat00003
[수학식 6]
Figure 112021062585564-pat00004
수학식 5에서, Vm(t)는 로터 허브에서의 풍속, β(t)는 블레이드의 피치각, λ(t)는 주속비로서 수학식 6으로 정의되고, CT(λ(t),β(t))는 추력계수이고, 수학식 6에서 Ωr(t)는 로터의 회전 각속도이다.
또한 상기 타워 섀도우 추력변동계수 연산부는 이하의 수학식 9에 의해 상기 타워 섀도우 추력변동계수를 연산하며, 상기 추력 모델링부는 이하의 수학식 10에 의해 타워 섀도우에 의해 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 것을 특징으로 할 수 있다.
[수학식 9]
Figure 112021062585564-pat00005
[수학식 10]
Figure 112021062585564-pat00006
수학식 9에서, FTi,ts(θ,λ)는 타워 섀도우 추력변동획득부에서 획득된 추력값이고, FTi(λ)는 상기 수학식 5에 의해 연산된 추력값이고, 수학식 10에서,
Figure 112021062585564-pat00007
는 타워 섀도우 추력변동계수이고, R은 블레이드의 길이이며, ρ는 공기밀도,
Figure 112021062585564-pat00008
는 추력계수이고, Vm(t)는 로터 허브에서의 풍속이다.
또한 상기 타워 섀도우 추력변동계수의 최대값은 최소값이 일어나는 로터 방위각에서 전후 30도 범위 내에서 각각에서 발생되는 것을 특징으로 할 수 있다.
그리고 상기 타워 섀도우 추력변동계수의 최대값과 최소값은 최적 주속비 구간에서 일정한 것을 특징으로 할 수 있다.
또한 최적주속비 구간보다 작은 풍속에서 주속비는 커지고, 최적 주속비 구간보다 큰 풍속에서는 주속비가 작아지며, 상기 타워 섀도우 추력변동계수는, 주속비가 최적 주속비보다 작아짐에 따라 추력변동비율이 작아지고, 주속비가 최적 주속비보다 커짐에 따라 추력변동 비율이 커지는 것을 특징으로 할 수 있다.
본 발명의 제2목적은 정격풍속 이하에서 풍력터빈의 추력 동하중 모델링 방법에 있어서, 추력연산부가 타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값을 연산하는 단계; 타워 섀도우 추력변동획득부가 타워 섀도우에 의해 변동하는 추력값을 획득하는 단계; 타워 섀도우 추력변동계수 연산부가 상기 추력연산부에서 연산된 추력값과, 상기 타워 섀도우 추력변동획득부에서 획득된 추력값을 기반으로 타워 섀도우 추력변동계수를 연산하는 단계; 데이터베이스가 상기 타워 섀도우 추력변동계수를 데이터베이스화하는 단계; 및 블레이드 추력 모델링부가 추력연산부에서 연산된 추력값과, 상기 데이터베이스에 저장된 타워 섀도우 추력변동계수를 이용하여 타워 섀도우에 의해 각 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 단계;를 포함하고, 상기 타워 섀도우에 의한 풍속 변동은 이하의 수학식 1 및 수학식 2에 의해 정의되고, 상기 추력연산부에 의해 연산되는 추력값은 이하의 수학식 5에 의해 연산되며, 상기 타워 섀도우 추력변동계수 연산부는 이하의 수학식 9에 의해 상기 타워 섀도우 추력변동계수를 연산하고, 상기 추력 모델링부는 이하의 수학식 10에 의해 타워 섀도우에 의해 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법로서 달성될 수 있다.
[수학식 1]
Figure 112021062585564-pat00009
[수학식 2]
Figure 112021062585564-pat00010
[수학식 5]
Figure 112021062585564-pat00011
[수학식 6]
Figure 112021062585564-pat00012
[수학식 9]
Figure 112021062585564-pat00013
[수학식 10]
Figure 112021062585564-pat00014
상기 수학식 1, 2에서, Vh는 허브 높이에서의 평균 풍속, vts는 타워 섀도우에 의해 변동하는 풍속, Dt는 타워의 지름, x는 타워의 중심부터 타워 앞에 있는 블레이드 중심까지의 거리, y는 풍속이 통과하는 로터의 회전면에서의 타워 중심부터 블레이드까지의 횡 방향 거리이고, 수학식 5에서, Vm(t)는 로터 허브에서의 풍속, β(t)는 블레이드의 피치각, λ(t)는 주속비로서 수학식 6으로 정의되고, CT(λ(t),β(t))는 추력계수이고, 수학식 6에서 Ωr(t)는 로터의 회전 각속도이며, 수학식 9에서, FTi,ts(θ,λ)는 타워 섀도우 추력변동획득부에서 획득된 추력값이고, FTi(λ)는 상기 수학식 5에 의해 연산된 추력값이고, 수학식 10에서,
Figure 112021062585564-pat00015
는 타워 섀도우 추력변동계수이고, R은 블레이드의 길이이며, ρ는 공기밀도,
Figure 112021062585564-pat00016
는 추력계수이고, Vm(t)는 로터 허브에서의 풍속이다.
그리고 터추력 모델링부가 상기 블레이드 추력 모델링부에 의해 모델링된 각 블레이드의 동하중 모델링을 기반으로 로터 추력 동하중을 모델링하는 단계;를 더 포함하는 것을 특징으로 할 수 있다.
또한 상기 타워 섀도우 추력변동계수의 최대값은 최소값이 일어나는 로터 방위각에서 전후 30도 범위 내에서 각각에서 발생되고, 상기 타워 섀도우 추력변동계수의 최대값과 최소값은 최적 주속비 구간에서 일정하며, 최적주속비 구간보다 작은 풍속에서 주속비는 커지고, 최적 주속비 구간보다 큰 풍속에서는 주속비가 작아지고, 상기 타워 섀도우 추력변동계수는, 주속비가 최적 주속비보다 작아짐에 따라 추력변동비율이 작아지고, 주속비가 최적 주속비보다 커짐에 따라 추력변동 비율이 커지는 것을 특징으로 할 수 있다.
본 발명의 실시예에 따르면, 타워 섀도우 추력변동계수를 정의하여 이러한 타워 섀도우 추력변동계수를 이용하여 정격풍속 이하에서, 타워 섀도우에 의해 블레이드에 발생하는 추력변동 동하중을 간편하고 신속하게 모델링할 수 있 는 효과를 갖는다.
또한, 본 발명의 일실시예에 따르면, 타워 섀도우에 대한 추력변동에 의한 동하중을 모델링하는데 있어, 공력계수를 이용한 BEM 방법이 아닌, 정의되는 타워 섀도우 추력변동계수와 추력계수를 이용하여 신속하고 효율적으로 타워 섀도우에 대한 추력변동을 모델링할 수 있는 효과를 갖는다.
한편, 본 발명에서 얻을 수 있는 효과는 이상에서 언급한 효과들로 제한되지 않으며, 언급하지 않은 또 다른 효과들은 아래의 기재로부터 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.
본 명세서에 첨부되는 다음의 도면들은 본 발명의 바람직한 실시예를 예시하는 것이며, 발명의 상세한 설명과 함께 본 발명의 기술적 사상을 더욱 이해시키는 역할을 하는 것이므로, 본 발명은 그러한 도면에 기재된 사항에만 한정되어 해석 되어서는 아니 된다.
도 1은 통상의 풍력터빈시스템의 부분 사시도,
도 2는 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 모델링 시스템의 구성을 나타낸 블록도,
도 3a 및 도 3b는 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 모델링 방법의 흐름도,
도 4는 포텐셜 흐름 모델의 개략도,
도 5는 2MW 풍력터빈 데모 모델의 물리적인 파라미터값 표,
도 6은 타워 섀도우가 없는 경우에 대하여 정격 이하의 풍속(4 ~ 11m/s)에서 풍속에 따른 블레이드 하나에서의 추력 그래프,
도 7은 허브에서의 풍속이 7m/s인 경우 방위각에 따른 타워 섀도우에 의한 추력의 변동 그래프,
도 8은 풍속이 7m/s인 경우 방위각에 따른 타워 섀도우 추력 변동계수 그래프,
도 9는 타워 섀도우가 있는 경우 정격 이하의 풍속(4 ~ 11m/s)에서 풍속에 따른 타워 섀도우 추력 변동계수의 최대값, 최소값 그래프,
도 10은 타워 섀도우가 있는 경우 정격 이하의 풍속(4 ~ 11m/s)에서 주속비에 따른 타워 섀도우 추력 변동계수의 최대값, 최소값 그래프,
도 11은 타워 섀도우가 있는 경우, 주속비와 방위각에 따른 섀도우 추력 변동계수,
도 12는 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 모델링 시스템의 구성을 나타낸 블록도,
도 13 및 도 14는 풍속 7m/s인 경우 타워 섀도우 영향으로 인해 세 블레이드에서 발생하는 추력변동(도 13)과, 로터 허브에서 발생하는 추력변동(도 14)에 대한 수치실험 결과를 도시한 것이다.
이상의 본 발명의 목적들, 다른 목적들, 특징들 및 이점들은 첨부된 도면과 관련된 이하의 바람직한 실시예들을 통해서 쉽게 이해될 것이다. 그러나 본 발명은 여기서 설명되는 실시예들에 한정되지 않고 다른 형태로 구체화될 수도 있다. 오히려, 여기서 소개되는 실시예들은 개시된 내용이 철저하고 완전해질 수 있도록 그리고 통상의 기술자에게 본 발명의 사상이 충분히 전달될 수 있도록 하기 위해 제공되는 것이다.
본 명세서에서, 어떤 구성요소가 다른 구성요소 상에 있다고 언급되는 경우에 그것은 다른 구성요소 상에 직접 형성될 수 있거나 또는 그들 사이에 제 3의 구성요소가 개재될 수도 있다는 것을 의미한다. 또한 도면들에 있어서, 구성요소들의 두께는 기술적 내용의 효과적인 설명을 위해 과장된 것이다.
본 명세서에서 기술하는 실시예들은 본 발명의 이상적인 예시도인 단면도 및/또는 평면도들을 참고하여 설명될 것이다. 도면들에 있어서, 막 및 영역들의 두께는 기술적 내용의 효과적인 설명을 위해 과장된 것이다. 따라서 제조 기술 및/또는 허용 오차 등에 의해 예시도의 형태가 변형될 수 있다. 따라서 본 발명의 실시예들은 도시된 특정 형태로 제한되는 것이 아니라 제조 공정에 따라 생성되는 형태의 변화도 포함하는 것이다. 예를 들면, 직각으로 도시된 영역은 라운드지거나 소정 곡률을 가지는 형태일 수 있다. 따라서 도면에서 예시된 영역들은 속성을 가지며, 도면에서 예시된 영역들의 모양은 소자의 영역의 특정 형태를 예시하기 위한 것이며 발명의 범주를 제한하기 위한 것이 아니다. 본 명세서의 다양한 실시예들에서 제1, 제2 등의 용어가 다양한 구성요소들을 기술하기 위해서 사용되었지만, 이들 구성요소들이 이 같은 용어들에 의해서 한정되어서는 안 된다. 이들 용어들은 단지 어느 구성요소를 다른 구성요소와 구별시키기 위해서 사용되었을 뿐이다. 여기에 설명되고 예시되는 실시예들은 그것의 상보적인 실시예들도 포함한다.
본 명세서에서 사용된 용어는 실시예들을 설명하기 위한 것이며 본 발명을 제한하고자 하는 것은 아니다. 본 명세서에서, 단수형은 문구에서 특별히 언급하지 않는 한 복수형도 포함한다. 명세서에서 사용되는 '포함한다(comprises)' 및/또는 '포함하는(comprising)'은 언급된 구성요소는 하나 이상의 다른 구성요소의 존재 또는 추가를 배제하지 않는다.
아래의 특정 실시예들을 기술하는데 있어서, 여러 가지의 특정적인 내용들은 발명을 더 구체적으로 설명하고 이해를 돕기 위해 작성되었다. 하지만 본 발명을 이해할 수 있을 정도로 이 분야의 지식을 갖고 있는 독자는 이러한 여러 가지의 특정적인 내용들이 없어도 사용될 수 있다는 것을 인지할 수 있다. 어떤 경우에는, 발명을 기술하는 데 있어서 흔히 알려졌으면서 발명과 크게 관련 없는 부분들은 본 발명을 설명하는데 있어 별 이유 없이 혼돈이 오는 것을 막기 위해 기술하지 않음을 미리 언급해 둔다.
이하에서는 본 발명에서는 추력계수를 이용한 추력 표현식을 적용할 때 타워 섀도우에 의해 블레이드에 발생하는 추력변동 동하중을 표현하는 방법에 대하여 제시한다. 이를 위하여 “타워 섀도우 추력변동 계수”를 정의 및 도입하고, 2MW 풍력터빈을 대상으로 정격이하의 풍속에서 “타워 섀도우 추력변동 계수”를 구하고 분석한 후 Matlab/Simulink에서 타워 섀도우 동하중 모델을 구현하고 검증한다.
이하에서는 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템(100)의 구성, 기능 그리고, 모델링 방법에 대해 설명하도록 한다. 먼저, 도 2는 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템(100)의 구성을 나타낸 블록도를 도시한 것이다. 그리고, 도 3a 및 도 3b는 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법의 흐름도를 도시한 것이다.
본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템(100)은 정의되는 타워 섀도우 추력변동계수와 추력계수를 이용하여 기존 공력계수를 이용한 BEM 방법을 적용하는 것보다 신속하고 효율적으로 타워 섀도우가 고려된 각 블레이드와 로터의 추력을 모델링할 수 있게 된다.
이러한 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템(100)은, 도 2에 도시된 바와 같이, 추력연산부(10)와, 타워 섀도우 추력변동획득부(20)와, 타워 섀도우 추력변동계수 연산부(30)와, 블레이드 추력 모델링부(40)와, 데이터베이스(50)와, 로터 추력 모델링부(60) 등을 포함하여 구성될 수 있음을 알 수 있다.
도 3a는 본 발명의 실시예에 따른 타워 섀도우 추력변동계수를 데이터베이스화하는 방법에 대한 흐름도를 도시한 것이고, 도 3b는 본 발명의 실시예에 따라 데이터베이스에 저장된 타워 섀도우 추력변동계수와 추력연산부(10)에서 연산되는 추력값을 기반으로 타워 섀도우에 의해 블레이드에 발생되는 추력변동에 의한 동하중과 로터에서 발생되는 추력변동에 의한 동하중을 모델링하는 방법의 흐름도를 도시한 것이다.
추력연산부(10)는 후에 상세하게 설명되는 바와 같이, 추력계수를 이용하여, 타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값을 연산하고(S1), 타워 섀도우 추력변동획득부(20)는 타워 섀도우에 의해 변동하는 추력값을 획득하게 된다(S2).
그리고, 타워 섀도우 추력변동계수 연산부(30)는 후에 상세히 설명되는 바와 같이, 추력연산부(10)에서 연산된 추력값과, 타워 섀도우 추력변동획득부(20)에서 획득된 추력값을 기반으로 타워 섀도우 추력변동계수를 연산하게 된다(S3). 그리고, 데이터베이스(50)는 추력계수와 연산된 타워 섀도우 추력변동계수를 저장하여 데이터베이스화하게 된다(S4).
그리고, 모델링방법은 도 3b에 도시된 바와 같이, 추력연산부(10)는 추력계수를 이용하여, 타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값을 연산하고(S11), 데이터베이스(50)로부터 타워 섀도우 추력변동계수를 결정하고(S12), 추력 모델링부(40)는 추력연산부(10)에서 연산되는 타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값(S11)과, 데이터베이스(50)에서 결정된 타워 섀도우 추력변동계수를 이용하여 타워 섀도우에 의해 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하게 된다(S13). 그리고 로터추력 모델링부(60)는 블레이드 추력 모델링부(40)에 의해 모델링된 각 블레이드의 동하중 모델링을 기반으로 로터 추력 동하중을 모델링하게 된다(S14).
이하에서는 공력계수를 이용한 기존 BEM 방법과, 종래 추력계수를 이용한 추력 연산방법에 대해 설명한 후, 본 발명의 실시예에 따른 타워 섀도우에 의한 풍력터빈 블레이드의 추력 모델링을 보다 상세하게 설명하도록 한다.
먼저, 풍력터빈에서 발생하는 타워 섀도우는 바람과 타워의 간섭에 의한 것이다. 바람이 풍력터빈으로 불어 들어오면 블레이드를 거친 후에 타워를 거쳐 지나가게 된다. 이때 타워는 바람의 흐름을 방해하게 되고, 타워를 지나는 지점의 풍속에 변화를 주게 된다. 이에 따라 타워 섀도우에 의해 추력의 변동이 발생하고 블레이드에는 동하중이 발생하게 된다. 로터가 회전할 때 타워 섀도우에 의해 블레이드에 발생하는 하중은 로터 회전각에 따라 달라진다. 세 블레이드는 로터에 120도 간격으로 결합되어 있으므로 세 블레이드에는 비대칭 하중이 발생하게 된다.
타워 섀도우에 의한 풍속의 변동을 계산하기 위한 방법으로는 대표적으로 포텐셜 흐름 모델(potential flow model), 타워 섀도우 결손 모델(tower-shadow-deficit model), 실증 모델(empirical model)이 있다. 특히 전방향형(upwind) 풍력터빈에 대하여 풍속의 변동을 계산할 때에는 포텐셜 흐름 모델을 적용할 경우 더 정확한 결과를 얻을 수 있는 것으로 알려져 있다. 따라서 본 연구에서는 전방향형 2MW 풍력터빈 대상으로 포텐셜 흐름 모델을 적용하여 수치실험을 진행하였다
도 1은 포텐셜 흐름 모델의 개략도이고, 포텐셜 흐름 모델에서 타워 섀도우에 의한 풍속의 변동은 수학식 1과 수학식 2로 표현된다.
[수학식 1]
Figure 112021062585564-pat00017
[수학식 2]
Figure 112021062585564-pat00018
수학식 1, 2에서, Vh는 허브 높이에서의 평균 풍속을 뜻하며, vts는 타워 섀도우에 의해 변동하는 풍속을 뜻하고, Dt는 타워의 지름, x는 타워의 중심부터 타워 앞에 있는 블레이드 중심까지의 거리를 뜻한다. y는 풍속이 통과하는 로터의 회전면에서의 타워 중심부터 블레이드까지의 횡 방향 거리를 뜻한다.
타워 섀도우가 고려되지 않은 경우에는 블레이드의 길이에 따라서 풍속은 변동이 없이 일정하며 허브 높이에서의 풍속과 같다. 하지만 타워 섀도우가 고려된 경우에는 블레이드의 길이와 로터 회전 시 변화하는 로터 방위각에 따라서 거리 y가 달라지기 때문에 풍속의 변동이 발생한다. 이로 인해서 타워 섀도우가 고려된 경우에는 로터 회전 시 로터 방위각이 변화함에 따라 각각 세 블레이드에는 풍속 변동에 의한 추력의 변동이 발생하게 된다. 또한 각각 세 블레이드에는 120도의 위상차가 존재하기 때문에 로터 회전 시 로터 방위각이 변화함에 따라 블레이드마다 각기 다른 비대칭 하중이 유발하게 된다.
풍력터빈이 바람에 의해 후방으로 밀리는 힘인 추력(FT)은 BEM(Blade Element Momentum) 방법을 이용하여 계산이 가능하다. 바람이 불어 풍력터빈의 로터-블레이드가 회전할 때 블레이드 단면에 작용하는 공력은 블레이드 익형의 양력계수(Cl)와 항력계수(Cd)의 함수이다. 블레이드 익형의 각 단면에서 발생하는 양력과 항력의 성분은 공력 토크와 추력을 발생시킨다. 블레이드 하나에서 발생하는 전체 추력은 블레이드 각 단면에서 발생하는 추력을 블레이드 길이에 따라 적분함으로 얻을 수 있는데 수학식 3과 같이 표현된다.
[수학식 3]
Figure 112021062585564-pat00019
여기서
Figure 112021062585564-pat00020
는 공기밀도, c는 익형 시위선의 길이, Vrel는 익형 앞전에서의 상대속도,
Figure 112021062585564-pat00021
는 익형 받음각과 블레이드 피치각의 합이다.
블레이드가 세 개인 풍력터빈의 경우 로터 허브에 작용하는 전체 추력은 각 블레이드의 추력(FTi)의 합으로 수학식 4와 같이 표현된다.
[수학식 4]
Figure 112021062585564-pat00022
타워 섀도우가 없는 경우에는 블레이드가 세 개인 풍력터빈의 로터 허브 작용하는 전체 추력은 블레이드 하나에 작용하는 추력인 수학식 3을 세 배하면 된다. 타워 섀도우가 없는 경우에는 풍속이 블레이드 길이에 따라서 일정하다. 그러나 타워 섀도우를 고려하면 풍속은 블레이드 길이에 따라 다르고 또한 로터의 회전각에 따라서 다르기 때문에, 각 블레이드에는 서로 다른 비대칭 하중이 작용하게 된다.
회전날개의 형상과 블레이드 익형에 관한 공력 특성이 알려져 있다면 수학식 3으로부터 추력계수(CT)를 계산할 수 있다. 바람에 의해 블레이드 하나에서 발생하는 풍력터빈의 추력은 수학식 5와 같이 표현되는데, 풍속뿐만 아니라 주속비와 피치각에 대하여 비선형적으로 표현된다. 본 발명의 실시예에 따른 추력연산부(에 의해 연산되는 추력값은 이하의 수학식 5에 의해 연산되게 된다.
[수학식 5]
Figure 112021062585564-pat00023
여기서 Vm은 로터 허브에서의 풍속, β는 블레이드의 피치각, λ는 주속비로 수학식 6과 같이 정의된다. 그리고 CT는 추력계수로 수학식 3과 4의 BEM 방법에서의 추력(FT) 표현을
Figure 112021062585564-pat00024
로 나눔으로 수학식 7과 같이 정의되며, 블레이드의 각 단면에서의 익형의 양력계수, 항력계수, 형상 등이 주어져야 계산이 가능하다.
[수학식 6]
Figure 112021062585564-pat00025
수학식 6에서 Ωr(t)는 로터의 회전 각속도이다.
[수학식 7]
Figure 112021062585564-pat00026
블레이드가 세 개인 풍력터빈의 경우 로터 허브에 작용하는 전체 추력은 각 블레이드의 추력(FTi)의 합으로 수학식 4와 같이 동일하게 표현된다. 타워 섀도우가 없는 경우 블레이드가 세 개인 풍력터빈에 작용하는 전체 추력은 블레이드 하나에 작용하는 추력인 수학식 5를 세 배하면 된다. 그러나 타워 섀도우가 있는 경우의 풍속은 블레이드 길이에 따라 다르고 또한 로터의 회전각에 따라서 다르기 때문에, 수학식 5를 이용하여 각 블레이드에 작용하는 서로 다른 비대칭 하중을 계산할 수 없다.
이하에서는 현장에 적용되는 2MW 풍력터빈 시스템에서 타워 섀도우에 의한 추력 특정을 살펴보도록 한다. GL Garrad Hassan의 Bladed S/W에서는 2MW 풍력터빈 데모 모델을 제공하고 있다. Bladed S/W에서는 블레이드, 로터, 타워, 파워 트레인 등 풍력터빈의 각 구성요소들을 표현할 수 있는 창들이 있는데, 거기에서 제시된 2MW 풍력터빈 데모 모델의 물리적인 파라미터값들을 정리하면 도 5와 같다. 최적 주속비 구간의 풍속은 5.3~9.3m/s이고 정격풍속은 11.7m/s이다. 풍력터빈의 출력계수(CP)는 최적 주속비(λopt)가 8.1일 때 최대값(CP,max)이 0.4662이다.
풍력터빈은 풍속과 블레이드의 상호작용이 비선형적인 특성을 가져 풍속에 따른 작동점이 제어기의 동특성에 영향을 받으므로, 토크 제어기와 피치 제어기의 설계가 우선적으로 고려되어야 한다. 정격풍속 이하에서는 블레이드의 피치각이 일정한 값으로 유지되므로 피치 제어기의 고려는 필요 없으나 토크 제어기의 고려는 반드시 필요하다. 정격풍속 이하에서는 풍력터빈의 응답(로터 속도, 추력 등) 특성은 토크 제어기의 형태에 따라 영향을 받는다. 토크 제어기는 풍력터빈에서 실용적으로 많이 사용되는 최적모드게인을 이용하는 토크제어 방법을 고려하여 적용하였다(. 발전기의 토크 크기는 최적모드게인에 발전기의 속도를 제곱하여 얻게 된다. 토크 제어기의 작동으로 인해 최적 주속비 구간에서는 주속비가 최적 주속비를 유지하게 되고 최대의 출력파워를 얻게 된다.
기어박스가 포함된 풍력터빈인 경우 이하의 수학식 8에서와 같이, 발전기 회전 각속도(Ωg)는 기어박스의 증속비(Ngb)와 상관관계가 있다.
[수학식 8]
Figure 112021062585564-pat00027
타워 섀도우가 없는 경우 풍력터빈의 세 개의 블레이드에 작용하는 추력은 모두 같고 각 블레이드의 추력이 합해져서 로터 허브에 추력을 발생시킨다. 앞에서 설명한 토크 제어기를 적용하여 로터 허브에서의 풍속에 따른 각 블레이드에서의 추력의 크기를 BEM(Blade Element Momentum)방법을 이용하여 계산하였다. 도 6은 타워 섀도우가 없는 경우에 대하여 정격 이하의 풍속(4 ~ 11m/s)에서 풍속에 따른 블레이드 하나에서의 추력의 정상상태 값을 보여주고 있는데, 풍속이 증가함에 따라 추력이 증가하는 것을 확인할 수 있다.
타워 섀도우가 있는 경우 풍력터빈의 세 개의 블레이드에 작용하는 추력은 시간에 따라 모두 다르고 각 블레이드에 작용하는 추력에 변동이 발생한다. 로터가 회전할 때 각 블레이드의 추력은 로터의 회전각에 따라서 변동하게 된다. 앞에서 설명한 토크 제어기를 적용하여 로터 허브에서의 풍속에 따른 각 블레이드에서의 추력의 크기를 BEM방법을 이용하여 계산하였다.
도 7은 허브에서의 풍속이 7m/s인 경우 타워 섀도우에 의한 추력의 변동을 보여주고 있다. 타워 섀도우에 의한 추력 변동의 최소값은 로터의 회전 방위각(θ)이 180도일 때 나타남을 확인할 수 있고, 추력 변동의 최대값은 최소값이 일어나는 로터 방위각 180도에서 전후 약 20도(약 160도와 200도)부근에서 발생함을 확인할 수 있다.
추력계수를 이용하는 수학식 5는 허브에서의 풍속의 변화에 따른 추력의 변화를 표현할 뿐이다. 수학식 5를 사용할 경우 허브에서의 풍속이 일정할 경우 평균 풍속에 대한 블레이드 하나에 작용하는 추력의 정상상태 값을 구할 수 있다.
그러나 추력계수를 이용하는 수학식 5는 타워 섀도우가 있는 경우 타워 섀도우에 의한 추력의 변동을 표현할 수 없다. 본 발명의 실시예에서는 추력계수를 이용하는 수학식 5를 기반으로 타워 섀도우에 대한 추력의 변동을 표현하고자 한다.
수학식 10과 같이 무차원 파라미터인 “타워 섀도우 추력변동 계수(Ctvts)”를 도입하여 타워 섀도우 동하중 모델을 표현하는 방법을 제안하고자 한다.
“타워 섀도우 추력변동 계수(Ctvts)”는 수학식 9와 같이 타워 섀도우에 의해 변동하는 추력값(FTi,ts)을 타워 섀도우가 없을 때의 추력값(FTi)으로 나눈 비로 정의하고자 한다. 피치각이 일정한 정격이하의 풍속인 경우 “타워 섀도우 추력변동 계수(Ctvts)”는 로터의 회전 방위각(θ)과 주속비(λ)의 함수로 표현할 수 있다.
[수학식 9]
Figure 112021062585564-pat00028
[수학식 10]
Figure 112021062585564-pat00029
정의된 “타워 섀도우 추력변동 계수(Ctvts)”는 타워 섀도우로 인한 추력의 변동이 타워 섀도우가 없는 경우의 추력보다 어느 정도의 비율로 증가 또는 감소하는지를 알려주는 것이라 볼 수 있다.
수학식 9에서 FTi는 수학식 5를 사용하여 구할 수 있고, FTi,ts는 수학식 3으로 구한 값을 이용하면 된다. 예를 들어 평균 풍속이 7m/s일 때, 타워 섀도우가 없는 경우의 추력(FTi)은 40.22kNm이고, 타워 섀도우가 있는 경우의 추력(FTi,ts)은 도 7과 같이 로터의 회전 방위각(θ)에 따라 변동한다. 수학식 9로부터 “타워 섀도우 추력변동 계수(Ctvts)”를 구하면 도 8과 같이 구할 수 있다.
도 9는 타워 섀도우가 있는 경우 정격 이하의 풍속(4 ~ 11m/s)에서 풍속에 따른 “타워 섀도우 추력변동 계수(Ctvts)”의 최대값과 최소값을 보여주고 있다.
도 9에 도시된 바와 같이, 풍속이 증가함에 따라 최대값과 최소값 모두 조금씩 커지는 경향임을 알 수 있다. 특히 풍속 6 m/s, 7 m/s, 8 m/s, 9 m/s에서는 최대값들과 최소값들은 변하지 않고 동일함을 주목할 필요가 있다.
이는 최적 주속비 구간(풍속 5.3 ~ 9.3 m/s)에서는 주속비의 값이 최적 주속비(λ0 = 8.1)로 일정하므로, “타워 섀도우 추력변동 계수(Ctvts)”가 동일함을 의미한다.
최적 주속비 구간에서 타워 섀도우가 있는 경우의 추력은 타워 섀도우가 없는 경우보다 최대 1.4% 정도 증가하고 최소 21.7% 정도 감소함을 확인할 수 있다.
도 10은 타워 섀도우가 있는 경우 정격 이하의 풍속(4 ~ 11m/s)에서 주속비 따른 “타워 섀도우 추력변동 계수(Ctvts)”의 최대값과 최소값을 보여준다. 최적 주속비 구간보다 작은 풍속에는 주속비가 커지고, 최적 주속비 구간보다 큰 풍속에는 주속비가 작아진다. 주속비가 최적 주속비보다 작아짐에 따라 추력변동 비율이 조금 작아지고(최대값 약 2.4% 증가, 최소값 약 18.8% 감소), 주속비가 최적 주속비보다 커짐에 따라 추력변동 비율이 약간 커짐(최대값 약 1.4% 증가, 최소값 약 22.2% 감소)을 확인할 수 있다.
정격 풍속이하에서는 제안된 수학식 9의 “타워 섀도우 추력변동 계수(Ctvts)”와 추력계수를 이용하는 수학식 10을 이용하여, 로터의 회전 방위각(θ)과 주속비(λ)에 따라 블레이드에 타워 섀도우로 인해 발생하는 추력 변동을 구현할 수 있다.
구해진 “타워 섀도우 추력변동 계수(Ctvts)”는 도 11과 같은데, 로터의 회전 방위각(θ)과 주속비(λ)의 함수로 표현된다. 세 블레이드는 120도 간격으로 설치되므로, 각 블레이드의 회전 방위각(θ1,θ2,θ3)은 수학식 11과 같다. 각 블레이드에 발생하는 추력 변동은 “타워 섀도우 추력변동 계수(Ctvts)”의 사용 시에 로터의 방위각에 0도와 120도와 240도를 각각 더하여 사용하면 된다.
[수학식 11]
Figure 112021062585564-pat00030
그리고 타워 섀도우로 인해 세 블레이드에 발생하는 비대칭 추력 변동으로 인해 로터 허브에 발생하는 추력 변동(FT,ts)은 수학식 12와 같이 세 블레이드의 각 추력을 더하면 된다.
[수학식 12]
Figure 112021062585564-pat00031
제안된 수학식 9의 “타워 섀도우 추력변동 계수(Ctvts)”와 추력계수를 이용하는 수학식 10을 이용하여 Matlab/Simulink를 이용하여 타워 섀도우의 동하중 모델을 구현하였다(도 12).
구현된 Matlab /Simulink를 이용하여 풍속에 따른 추력변동을 확인하였다. 도 13 및 도 14는 풍속 7m/s인 경우 타워 섀도우의 영향으로 인해 세 블레이드에서 발생하는 추력 변동(FT1,ts, FT2,ts, FT3,ts)과 로터 허브에 발생하는 추력 변동(FT,ts)에 대한 수치실험 결과를 각각 보여주고 있다. 추력계수를 이용하지 않고 공력계수를 직접 이용하는 BEM 방법의 결과와 일치하는 동일한 응답을 확인할 수 있었다.
본 발명의 실시예에서는 추력계수를 이용하여 타워 섀도우에 의한 추력변동 동하중을 표현하는 방법을 제안하였다. 이를 위하여 “타워 섀도우 추력변동 계수”를 도입하였고, 이를 타워 섀도우에 의해 변동하는 추력값을 타워 섀도우가 없을 때의 추력값으로 나눈 비로 정의하였다. 피치각이 일정한 정격이하의 풍속에서 “타워 섀도우 추력변동 계수”를 로터의 회전 방위각과 주속비의 함수로 표현할 수 있음을 확인하였다.
2MW 육상용 풍력터빈을 대상으로 “타워 섀도우 추력변동 계수”를 구하였는데, 최적 주속비 구간에서 이의 최대값은 약 1.014이고 최소값은 약 0.783임을 확인하였다. 구해진 “타워 섀도우 추력변동 계수”와 추력계수를 이용하여 Matlab/Simulink에서 타워 섀도우 동하중 모델을 구현하였고, 블레이드에 작용하는 추력변동을 추력계수와 “타워 섀도우 추력변동 계수”를 동시에 이용하여 표현할 수 있음을 보였다.
또한, 상기와 같이 설명된 장치 및 방법은 상기 설명된 실시예들의 구성과 방법이 한정되게 적용될 수 있는 것이 아니라, 상기 실시예들은 다양한 변형이 이루어질 수 있도록 각 실시예들의 전부 또는 일부가 선택적으로 조합되어 구성될 수도 있다.
1:풍력터빈 시스템
2:로터
3:블레이드
4:발전기
5:기어박스
6:허브
7:타워
10:추력연산부
20:타워 섀도우 추력변동획득부
30:타워 섀도우 추력변동계수 연산부
40:블레이드 추력 모델링부
50:데이터베이스
60:로터추력 모델링부
100:타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템

Claims (11)

  1. 정격풍속 이하에서 풍력터빈의 추력 동하중 모델링 시스템에 있어서,
    타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값을 연산하는 추력연산부;
    타워 섀도우에 의해 변동하는 추력값을 획득하는 타워 섀도우 추력변동획득부;
    상기 추력연산부에서 연산된 추력값과, 상기 타워 섀도우 추력변동획득부에서 획득된 추력값을 기반으로 타워 섀도우 추력변동계수를 연산하는 타워 섀도우 추력변동계수 연산부;
    상기 타워 섀도우 추력변동계수를 데이터베이스화하는 데이터베이스; 및
    상기 추력연산부에서 연산된 추력값과, 상기 데이터베이스에 저장된 타워 섀도우 추력변동계수를 이용하여 타워 섀도우에 의해 각 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 블레이드 추력 모델링부;를 포함하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템.
  2. 제 1항에 있어서,
    상기 블레이드 추력 모델링부에 의해 모델링된 각 블레이드의 동하중 모델링을 기반으로 로터 추력 동하중을 모델링하는 로터추력 모델링부;를 더 포함하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템.
  3. 제 1항에 있어서,
    상기 타워 섀도우에 의한 풍속 변동은 이하의 수학식 1 및 수학식 2에 의해 정의되는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템:
    [수학식 1]
    Figure 112021062585564-pat00032

    [수학식 2]
    Figure 112021062585564-pat00033

    상기 수학식 1, 2에서, Vh는 허브 높이에서의 평균 풍속, vts는 타워 섀도우에 의해 변동하는 풍속, Dt는 타워의 지름, x는 타워의 중심부터 타워 앞에 있는 블레이드 중심까지의 거리, y는 풍속이 통과하는 로터의 회전면에서의 타워 중심부터 블레이드까지의 횡 방향 거리이다.
  4. 제 3항에 있어서,
    상기 추력연산부에 의해 연산되는 추력값은 이하의 수학식 5에 의해 연산되는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템:
    [수학식 5]
    Figure 112022104290600-pat00034

    [수학식 6]
    Figure 112022104290600-pat00035

    수학식 5에서, Vm(t)는 로터 허브에서의 풍속, β(t)는 블레이드의 피치각, λ(t)는 주속비로서 수학식 6으로 정의되고, CT(λ(t),β(t))는 추력계수이고, 수학식 6에서 Ωr(t)는 로터의 회전 각속도이고, R은 블레이드의 길이이며, ρ는 공기밀도이다.
  5. 제 4항에 있어서,
    상기 타워 섀도우 추력변동계수 연산부는 이하의 수학식 9에 의해 상기 타워 섀도우 추력변동계수를 연산하며,
    상기 추력 모델링부는 이하의 수학식 10에 의해 타워 섀도우에 의해 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템:
    [수학식 9]
    Figure 112021062585564-pat00036

    [수학식 10]
    Figure 112021062585564-pat00037

    수학식 9에서, FTi,ts(θ,λ)는 타워 섀도우 추력변동획득부에서 획득된 추력값이고, FTi(λ)는 상기 수학식 5에 의해 연산된 추력값이고,
    수학식 10에서,
    Figure 112021062585564-pat00038
    는 타워 섀도우 추력변동계수이고, R은 블레이드의 길이이며, ρ는 공기밀도,
    Figure 112021062585564-pat00039
    는 추력계수이고, Vm(t)는 로터 허브에서의 풍속이다.
  6. 제 5항에 있어서,
    상기 타워 섀도우 추력변동계수의 최대값은 최소값이 일어나는 로터 방위각에서 전후 30도 범위 내에서 각각에서 발생되는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템.
  7. 제 6항에 있어서,
    상기 타워 섀도우 추력변동계수의 최대값과 최소값은 최적 주속비 구간에서 일정한 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템.
  8. 제 7항에 있어서,
    최적주속비 구간보다 작은 풍속에서 주속비는 커지고, 최적 주속비 구간보다 큰 풍속에서는 주속비가 작아지며,
    상기 타워 섀도우 추력변동계수는, 주속비가 최적 주속비보다 작아짐에 따라 추력변동비율이 작아지고, 주속비가 최적 주속비보다 커짐에 따라 추력변동 비율이 커지는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 시스템.
  9. 정격풍속 이하에서 풍력터빈의 추력 동하중 모델링 방법에 있어서,
    추력연산부가 타워 섀도우가 없을 때의 풍력터빈 블레이드의 추력값을 연산하는 단계;
    타워 섀도우 추력변동획득부가 타워 섀도우에 의해 변동하는 추력값을 획득하는 단계;
    타워 섀도우 추력변동계수 연산부가 상기 추력연산부에서 연산된 추력값과, 상기 타워 섀도우 추력변동획득부에서 획득된 추력값을 기반으로 타워 섀도우 추력변동계수를 연산하는 단계;
    데이터베이스가 상기 타워 섀도우 추력변동계수를 데이터베이스화하는 단계; 및
    블레이드 추력 모델링부가 추력연산부에서 연산된 추력값과, 상기 데이터베이스에 저장된 타워 섀도우 추력변동계수를 이용하여 타워 섀도우에 의해 각 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 단계;를 포함하고,
    상기 타워 섀도우에 의한 풍속 변동은 이하의 수학식 1 및 수학식 2에 의해 정의되고,
    상기 추력연산부에 의해 연산되는 추력값은 이하의 수학식 5에 의해 연산되며,
    상기 타워 섀도우 추력변동계수 연산부는 이하의 수학식 9에 의해 상기 타워 섀도우 추력변동계수를 연산하고,
    상기 추력 모델링부는 이하의 수학식 10에 의해 타워 섀도우에 의해 블레이드에 발생되는 추력변동에 의한 동하중을 모델링하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법:
    [수학식 1]
    Figure 112021062585564-pat00040

    [수학식 2]
    Figure 112021062585564-pat00041

    [수학식 5]
    Figure 112021062585564-pat00042

    [수학식 6]
    Figure 112021062585564-pat00043

    [수학식 9]
    Figure 112021062585564-pat00044

    [수학식 10]
    Figure 112021062585564-pat00045

    상기 수학식 1, 2에서, Vh는 허브 높이에서의 평균 풍속, vts는 타워 섀도우에 의해 변동하는 풍속, Dt는 타워의 지름, x는 타워의 중심부터 타워 앞에 있는 블레이드 중심까지의 거리, y는 풍속이 통과하는 로터의 회전면에서의 타워 중심부터 블레이드까지의 횡 방향 거리이고,
    수학식 5에서, Vm(t)는 로터 허브에서의 풍속, β(t)는 블레이드의 피치각, λ(t)는 주속비로서 수학식 6으로 정의되고, CT(λ(t),β(t))는 추력계수이고, 수학식 6에서 Ωr(t)는 로터의 회전 각속도이며,
    수학식 9에서, FTi,ts(θ,λ)는 타워 섀도우 추력변동획득부에서 획득된 추력값이고, FTi(λ)는 상기 수학식 5에 의해 연산된 추력값이고,
    수학식 10에서,
    Figure 112021062585564-pat00046
    는 타워 섀도우 추력변동계수이고, R은 블레이드의 길이이며, ρ는 공기밀도,
    Figure 112021062585564-pat00047
    는 추력계수이고, Vm(t)는 로터 허브에서의 풍속이다.
  10. 제 9항에 있어서,
    로터추력 모델링부가 상기 블레이드 추력 모델링부에 의해 모델링된 각 블레이드의 동하중 모델링을 기반으로 로터 추력 동하중을 모델링하는 단계;를 더 포함하는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법.
  11. 제 9항에 있어서,
    상기 타워 섀도우 추력변동계수의 로터 방위각에서 전후 30도 범위 내에서 각각에서 발생되고,
    상기 타워 섀도우 추력변동계수의 최대값과 최소값은 최적 주속비 구간에서 일정하며,
    최적주속비 구간보다 작은 풍속에서 주속비는 커지고, 최적 주속비 구간보다 큰 풍속에서는 주속비가 작아지고,
    상기 타워 섀도우 추력변동계수는, 주속비가 최적 주속비보다 작아짐에 따라 추력변동비율이 작아지고, 주속비가 최적 주속비보다 커짐에 따라 추력변동 비율이 커지는 것을 특징으로 하는 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법.
KR1020210070030A 2021-05-31 2021-05-31 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템 KR102465767B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020210070030A KR102465767B1 (ko) 2021-05-31 2021-05-31 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020210070030A KR102465767B1 (ko) 2021-05-31 2021-05-31 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템

Publications (1)

Publication Number Publication Date
KR102465767B1 true KR102465767B1 (ko) 2022-11-10

Family

ID=84101867

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020210070030A KR102465767B1 (ko) 2021-05-31 2021-05-31 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템

Country Status (1)

Country Link
KR (1) KR102465767B1 (ko)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009114977A (ja) 2007-11-07 2009-05-28 Shinko Electric Co Ltd 風力発電装置
JP2010281297A (ja) 2009-06-08 2010-12-16 Chia Yuan Li 風力発電機
KR20140041223A (ko) 2012-09-27 2014-04-04 한국전력공사 최소 하중에 의한 풍력터빈의 출력 제어방법
KR20150126127A (ko) 2014-05-02 2015-11-11 (주)제이피엠엔지니어링 해상풍력타워의 모노파일 기초에 대한 다층 지반해석방법
KR101734757B1 (ko) * 2016-03-29 2017-05-11 한밭대학교 산학협력단 윈드쉬어에 의한 풍력터빈 블레이드의 추력 모델링 방법 및 시스템
JP2020139427A (ja) * 2019-02-27 2020-09-03 株式会社日立製作所 風力発電システム、風量発電システムの制御方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009114977A (ja) 2007-11-07 2009-05-28 Shinko Electric Co Ltd 風力発電装置
JP2010281297A (ja) 2009-06-08 2010-12-16 Chia Yuan Li 風力発電機
KR20140041223A (ko) 2012-09-27 2014-04-04 한국전력공사 최소 하중에 의한 풍력터빈의 출력 제어방법
KR20150126127A (ko) 2014-05-02 2015-11-11 (주)제이피엠엔지니어링 해상풍력타워의 모노파일 기초에 대한 다층 지반해석방법
KR101734757B1 (ko) * 2016-03-29 2017-05-11 한밭대학교 산학협력단 윈드쉬어에 의한 풍력터빈 블레이드의 추력 모델링 방법 및 시스템
JP2020139427A (ja) * 2019-02-27 2020-09-03 株式会社日立製作所 風力発電システム、風量発電システムの制御方法

Similar Documents

Publication Publication Date Title
Lindenburg Investigation into rotor blade aerodynamics
KR101734757B1 (ko) 윈드쉬어에 의한 풍력터빈 블레이드의 추력 모델링 방법 및 시스템
Hansen Aerodynamics of wind turbines
do Rio et al. An extension of the blade element momentum method applied to diffuser augmented wind turbines
Lanzafame et al. Horizontal axis wind turbine working at maximum power coefficient continuously
US7393177B2 (en) Vertical axis wind turbine with optimized blade profile
Liu et al. Development and application of an improved blade element momentum method model on horizontal axis wind turbines
Martinez et al. An improved BEM model for the power curve prediction of stall‐regulated wind turbines
Kim et al. Aerodynamic design and performance analysis of multi-MW class wind turbine blade
Boorsma et al. Rotor experiments in controlled conditions continued: New Mexico
Hauptmann et al. Comparison of the lifting-line free vortex wake method and the blade-element-momentum theory regarding the simulated loads of multi-MW wind turbines
Madsen Low frequency noise from wind turbines mechanisms of generation and its modelling
Maniaci An investigation of WT_Perf convergence issues
CN108138747B (zh) 确定和控制定速风力涡轮机叶片的攻角的方法
Lee et al. Design optimization of wind turbine blades for reduction of airfoil self-noise
Hamlaoui et al. Improved BEM method for HAWT performance predictions
Schaffarczyk et al. Comparison of 3D transitional CFD simulations for rotating wind turbine wings with measurements
KR102465767B1 (ko) 정격풍속 이하에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템
Buck et al. Redefining the design objectives of large offshore wind turbine rotors
KR102494317B1 (ko) 정격풍속 이상에서 타워 섀도우에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템
Mesquita et al. An improved approach for performance prediction of HAWT using the strip theory
KR102519428B1 (ko) 정격풍속 이상에서 윈드쉬어에 의한 풍력터빈의 추력 동하중 모델링 방법 및 시스템
JP6865672B2 (ja) 風力発電装置の評価方法および設計方法
Xu et al. Development and application of a dynamic stall model for rotating wind turbine blades
Debbache et al. Analytical Performance Study of Fixed Speed Wind Turbine

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant