KR20160007838A - 풍력터빈 설계 및 해석방법 - Google Patents

풍력터빈 설계 및 해석방법 Download PDF

Info

Publication number
KR20160007838A
KR20160007838A KR1020140082993A KR20140082993A KR20160007838A KR 20160007838 A KR20160007838 A KR 20160007838A KR 1020140082993 A KR1020140082993 A KR 1020140082993A KR 20140082993 A KR20140082993 A KR 20140082993A KR 20160007838 A KR20160007838 A KR 20160007838A
Authority
KR
South Korea
Prior art keywords
blade
design
analysis
shape
data
Prior art date
Application number
KR1020140082993A
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 KR1020140082993A priority Critical patent/KR20160007838A/ko
Publication of KR20160007838A publication Critical patent/KR20160007838A/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
    • F03D1/00—Wind motors with rotation axis substantially parallel to the air flow entering the rotor 
    • F03D1/06—Rotors
    • F03D1/0608—Rotors characterised by their aerodynamic shape
    • F03D1/0633—Rotors characterised by their aerodynamic shape of the blades
    • 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
    • F03D3/00—Wind motors with rotation axis substantially perpendicular to the air flow entering the rotor 
    • F03D3/06—Rotors
    • F03D3/061—Rotors characterised by their aerodynamic shape, e.g. aerofoil profiles
    • 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
    • F05B2240/00—Components
    • F05B2240/20—Rotors
    • F05B2240/21—Rotors for wind turbines
    • 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
    • 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
    • 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/74—Wind turbines with rotation axis perpendicular to the wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (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

본 발명은 풍력터빈의 설계 및 해석을 간편하게 수행하기 위한 풍력터빈 설계 및 해석방법에 관한 것이다. 이를 위해 풍력터빈 핵심부품인 블레이드의 개념설계, 모델링(3차원 형상 모델링의 수행), 격자(유한요소 모델링의 수행), 하중조건, 응력평가 등에 대하여 후처리 결과로부터 자동으로 보고서를 생성해 주는 기능을 제공한다.

Description

풍력터빈 설계 및 해석방법{Method for designing and analyzing wind turbine blade}
본 발명은 풍력터빈의 주요부에 대하여 설계 및 해석을 간편하게 수행하기 위한 풍력터빈 설계 및 해석방법에 관한 것이다.
이를 위해 풍력터빈 핵심부품인 블레이드의 개념설계, 모델링(3차원 형상 모델링의 수행), 격자(유한요소 모델링의 수행), 하중조건, 응력평가 등에 대하여 후처리 결과로부터 자동으로 보고서를 생성해 주는 기능을 제공한다.
풍력산업은 바람을 이용한 발전단가가 저렴한 친환경 에너지로 기존의 화석연료 발전방식과 경쟁이 가능한 유망한 산업분야이다.
우리나라 에너지원 별 발전량을 비교하면 풍력은 수력 다음으로 높은 비중을 차지하며 향후에도 꾸준한 비율 상승이 예상된다.
하지만 풍력 부품 및 풍력 시스템의 시발점인 설계 SW, 해석 SW 및 엔지니어링 서비스 등은 실직적으로 해외 기술에 의존하고 있고 하드웨어와 소프트웨어의 조화로운 성장이 절실하다.
또한, 풍력터빈에서 가장 중요한 부품인 블레이드의 경우 전문적인 지식이 없으면, 설계 및 해석이 매우 어렵기 때문에 여기에 투자되는 시간과 인적자원이 매우 크다.
블레이드의 경우 일정한 섹션(Section)마다 에어포일(Airfoil)이라고 부르는 단면 형상을 지니게 된다. 이 에어포일의 경우 각 섹션마다 모두 동일한 것이 아니며, 블레이드의 설계 시 사용되는 단면 형상 데이타(Position, Chord length, Twist angle) 역시 모두 다르기 때문에, 모델링 시 각 섹션에 따른 반복적인 작업이 필수로 진행되게 된다.
블레이드 해석의 경우 격자를 생성하는 작업을 통해 유한요소 모델을 만들어주어야 하며, 하중과 구속 조건을 적용하는데 있어 전문적인 지식이 필요하다. 하지만, 블레이드의 설계 및 해석을 원하는 곳에서 항상 이러한 인적 자원이 마련되어 있지 않다는 문제점을 가진다.
블레이드의 설계 및 해석에 있어 섹션에 따른 형상의 변경이나 모델의 변경이 있을 때, 처음부터 그 모델링을 다시 해야만 하며, 이 때 유사한 형태의 모델링을 함에도 불구하고 시간의 투자가 많다는 문제점을 지닌다.
없음.
본 발명은 풍력터빈 주요부품에 대한 설계, 해석 시에 가지는 문제점을 해결하고자 안출된 것으로, 특히 블레이드 주요부품에 대하여 단순한 설계와 해석에 대한 정보의 입력 만으로도 간편하고 전문적인 지식이 없이도 설계와 해석을 수행하며 작업에 소모되는 시간을 대폭 절감하기 위한 엔지니어링 서비스 융합 소프트웨어 이다.
블레이드 형상설계용 공개프로그램과 연동됨과 동시에 블레이드형상을 결정하는 디자인모듈(Design Module), 블레이드 해석용 상용 프로그램과의 데이터 교환을 위한 연결모듈(Connection Module), 디자인모듈로부터 설계된 블레이드 형상에 대한 공력해석을 수행하는 블레이드 해석용 상용 프로그램을 준비해 두는 단계(s1);
상기 블레이드 형상설계용 공개프로그램의 경로를 미리 설정하여 상기 연결모듈에 의하여 데이터 교환이 가능하도록 하는 단계(s2);
상기 디자인 모듈에서,
입력된 에어포일(Airfoil)의 형상 이름과 레이놀즈수(Reynolds Number), 마하수(Mach Number), 유체밀도(Fluid Density), 동점성계수(Kinetic Viscosity) 데이터가 블레이드 형상설계용 공개프로그램으로 전송되는 단계(s3);
블레이드 형상설계용 공개프로그램에서, 상승력계수(CL), 저항계수(CD), 코드길이(Chord length), 포지션(Position), 트위스트앵글(Twist angle)값이 계산된 후, 다시 디자인모듈로 전송되는 단계(s4);
상기 디자인 모듈로 전송된, 상승력계수와 저항계수와 코드길이와 포지션데이터(에어포일의 섹션별 위치를 표현하는 데이터)과 트위스트앵글값을 이용하여, 블레이드 직경, 주속비 및 회전수를 계산하고, 날개 끝 손실개수와 축, 접선 방향의 유도계수를 계산하며 코드 길이, 유입각 및 트위스트앵글(비틀림각)을 계산하여 블레이드 형상이 모델링 되는 단계(s5);
상기 모델링된 블레이드 형상데이터가 연결모듈에 의하여, 블레이드 해석용 상용 프로그램으로 전송되는 단계(s6);
전송된 블레이드 형상에 대하여, 공력해석 모델을 생성한 뒤, 블레이드 공력해석을 수행(CFD)되어, 공력해석 결과와 최초 블레이드 설계치를 비교 검증하여, 수정된 동력계수를 피드백함으로써 블레이드 형상을 결정하는 단계(s7);
상기 s7 단계의 블레이드 공력해석결과 및 결정된 블레이드 형상에 대한 보고서가 출력되는 단계(s8); 을 포함하는 풍력터빈 블레이드 설계 및 해석방법이 제시된다.
본 발명의 풍력터빈 블레이드의 설계 및 해석을 위한 소프트웨어는 블레이드의 각 섹션에 따라 다르게 지정되는 Airfoil의 형상 및 설계 데이타(Position, Chord length, Twist angle)를 바탕으로 자동적으로 모델링이 이루어 지며, 이를 위해 사용자에 입력해야 할 정보는 기존의 작업에 비해 매우 간소하다.
기존의 블레이드 설계 시 사용자가 직접 블레이드의 섹션 마다 에어포일의 형상과 설계 데이타를 적용하여 모델링을 하고, 격자 작업 및 하중과 구속 조건의 적용을 수동적으로 해주어야 했던 것에 비하여 초기에 입력된 정보를 바탕으로 섹션 별 구성 및 형상의 구현과 격자 구현, 하중 및 구속의 적용까지 자동으로 이루어지기 때문에 기존에 발생하던 휴먼 에러(human error)의 발생을 최소한으로 줄일 수 있으며, 발생시에도 발생 원인을 쉽게 찾아 낼 수가 있다. 또한, 자동화로 인해 기존의 작업에 비하여 투자 시간이 비약적으로 줄어들게 된다.
기존의 블레이드 설계 및 해석 방식에 비하여 사용자에 의한 휴먼 에러가 줄고, 항상 일정한 규격에 따라 작업이 수행되기 때문에 해석상의 오차를 줄일 수 있는 효과를 가진다.
블레이드의 다양한 설계 데이타의 변경이나 에어포일 형상 변경에 대해서 손쉽고 빠르게 변경하여 해석 결과를 확인 할 수 있으므로, 다양한 설계 안의 검토가 가능하며 풍력터빈의 기타 주요부품(메인샤프트, 타워)에 대한 설계 등에도 적용가능하다.
도 1은 본 발명의 소프트웨어 전체 구성에 대한 개략적인 설명을 나타내는 흐름도
도 2a, 도 2b는 블레이드 설계 알고리즘
도 3은 풍력터빈 블레이드 설계식의 참고도
도 4는 본 발명의 해석 알고리즘 모듈 프로그래밍에 필요한 ANSYS WORKBENCH의 내부 구성언어 순서도
도 5는 본 발명의 플랫폼 모듈 설계도
도 6은 본 발명의 GUI 참고도
도 7은 본 발명의 해석 절차를 나타낸 순서도
도 8은 본 발명을 통해 구현된 블레이드의 모델링 결과
도 9는 구체적인 형상 설정 기능을 통하여 사용자가 원하는 코드길이, 포지션(Position), 트위스트앵글(Twist angle)을 가진 형태의 블레이드를 생성하는 장면
도 10은 본 발명을 통해 구현된 블레이드의 구조 및 유동 해석의 결과
도 11은 ANSYS를 통한 매뉴얼 방식의 블레이드 모델링 결과
도 12 및 도 13은 본 발명에 의한 오토매틱 방식의 블레이드 모델링 결과
도 14는 블레이드 출력성능 검증을 위한 자료이다.
도 15는 타워의 형상 결정후 CFD 해석과 구조해석이 수행된 결과 화면에 해당한다.
이하, 본 발명의 풍력터빈 블레이드의 설계 및 해석 소프트웨어에 대한 바람직한 실시예를 첨부된 도면을 참조하여 설명하도록 한다.
도 1은 본 발명의 소프트웨어 전체 구성에 대한 개략적인 설명을 나타내는 흐름도, 도 2a, 도 2b는 블레이드 설계 알고리즘, 도 3은 풍력터빈 블레이드 설계식의 참고도, 도 4는 본 발명의 해석 알고리즘 모듈 프로그래밍에 필요한 ANSYS WORKBENCH의 내부 구성언어 순서도, 도 5는 본 발명의 플랫폼 모듈 설계도, 도 6은 본 발명의 GUI 참고도, 도 7은 본 발명의 해석 절차를 나타낸 순서도, 도 8은 본 발명을 통해 구현된 블레이드의 모델링 결과, 도 9는 구체적인 형상 설정 기능을 통하여 사용자가 원하는 코드길이, 포지션(Position), 트위스트앵글(Twist angle)을 가진 형태의 블레이드를 생성하는 장면, 도 10은 본 발명을 통해 구현된 블레이드의 구조 및 유동 해석의 결과, 도 11은 ANSYS를 통한 매뉴얼 방식의 블레이드 모델링 결과, 도 12 및 도 13은 본 발명에 의한 오토매틱 방식의 블레이드 모델링 결과, 도 14는 블레이드 출력성능 검증을 위한 자료이다.
1. 본 발명의 구성
도 1은 본 발명의 풍력터빈 블레이드의 설계 및 해석 소프트웨어의 전체 구성에 대한 흐름도를 나타내고 있다
편의상, 본 발명의 풍력터빈 블레이드 설계 및 해석 소프트웨어를 "윈블레이드(WINBLADE)"로 명명하기로 하며, 도면에서도 이 용어를 사용한다.
본 발명에서의 소프트웨어는 크게 두 개의 외부 프로그램과 연동이 된다.
즉, 본 발명의 윈블레이드는, 블레이드 설계시 블레이드 형상설계용 공개 프로그램인 엑스포일(xfoil), 해석은 블레이드 해석용 상용 프로그램인 앤시스(ANSYS)와 연동할 수 있는 모듈을 구비하여야 한다. 기타 유사기능을 수행하는 상용 프로그램이 사용될 수도 있으며, 상기 프로그램들이 통상 많이 사용되므로 이를 사용하는 것은 일례이다.
본 발명의 윈블레이드(WINBLADE)와 엑스포일(xfoil) 사이에서는, 풍력터빈 블레이드의 설계 데이터를 주고 받게 되며, 여기서 도출된 블레이드 설계값을 ANSYS와 연동하여 계산을 수행 한 후 성능 해석 결과값을 확인할 수 있는 형태로 구성되어 있다.
따라서, 본 발명의 풍력터빈 블레이드의 설계 및 해석 소프트웨어의 구동에 있어서 공개소프트웨어인 xfoil과 상용 유한요소 해석 프로그램인 ANSYS 등이 필요하다. 기타 유사기능의 설계, 해석용 프로그램이 적용될 수도 있다.
상기 블레이드 형상설계용 공개프로그램(xfoil)의 경로를 미리 설정하여 상기 연결모듈을 통한 데이터 교환이 가능하도록 하는 것이 필요하다.
2. 플랫폼 모듈( Platform Module )
첨부된 도 5는 본 발명의 본 발명에서 구축된, 풍력터빈 블레이드 설계 및 해석 소프트웨어인 윈블레이드(WINBLADE)의 전체적인 프로그램 플랫폼 디자인 블록도로 나타낸 것이다
전체 본 발명에서 제시하는 윈블레이드 소프트웨어의 플랫폼은 디자인모듈(Design Module), ANSYS 연결모듈(Connection Module), 리포트 모듈(Report Module), 제어시스템 모듈(Conrol system Module) 및 수명주기모듈(Life cycle Module) 로 구성되어 있다.
디자인모듈(Design Module)의 여러 모듈 중에서, 블레이드 Module의 경우, 에어포일디자인(Airfoil Design)과 로터디자인(Rotor Design)으로 나뉜다.
사용자가 원하는 Airfoil의 형상 이름과 Reynolds Number, Mach Number, Fluid Density, Kinetic Viscosity, NCrit 와 같은 데이터를 대화창을 통하여 입력하게 되면, 공개 프로그램인 엑스포일(xfoil)을 통하여 결과 값으로 CL, CD, 코드길이(Chord length), 포지션(Position), 트위스트앵글(Twist angle) 값을 받을 수 있다.
본 발명의 윈블레이드(WINBLADE) 소프트웨어의 구동 절차를 설명한다.
도 6은 본 발명의 풍력터빈 블레이드 설계 및 해석 소프트웨어의 GUI 구성을 나타내고 있다.
먼저, 사용자가 원하는 블레이드의 설계 데이터를 윈블레이드 그래픽유저인터페이스(GUI) 창을 통해서 입력하게 되면, 자동적으로 xfoil과 연동하여 그래프 형태의 결과와 블레이드에 대한 3차원 형태의 모델링 결과를 표시해주게 된다.
여기서 사용자는 좀 더 구체적으로 원하는 형상의 블레이드 모양을 설정할 수 있으며, 설정된 데이터를 ANSYS와 연동하여 구조 해석 및 유동 해석을 수행하게 되며, 그 결과도 확인할 수 있게 된다.
도 6 상단에 프로그램 설정 및 구동에 필요한 메인메뉴(Main Menu)와 자주 사용하는 명령에 대한 툴 아이콘(Tool Icon)이 배치되어 있으며, 왼쪽에는 해석의 진행과 설정 상황을 보여주는 트리메뉴(Tree Menu)가 있다. 중앙에는 결과 그래프와 3D 모델을 보여주는 Graph & 3D View 메뉴가 구성되어 있다. 왼쪽 중앙부와 오른쪽 측면에는 선택한 기능에 대해 각각의 상세설정이 가능한 프라퍼티(Properties) 메뉴로 구성이 되어 있으며, 하단에는 작업의 진행 과정과 에러메시지(Error Message) 등의 출력이 가능한 Message & Output 창으로 구성되어 있다.
도 7에서는 본 발명의 풍력터빈 블레이드 설계 및 해석 소프트웨어의 해석 절차에 대한 순서를 표시하고 있다.
본 발명의 경우 필수적으로 구동에 있어서 엑스포일(xfoil)과 ANSYS를 필요로 하기 때문에 사용에 앞서 이들 프로그램들을 설치해주어야 하며, 프로그램 구동시 우선적으로 콘피규어(configure) 기능을 통해 설치된 xfoil의 경로를 정의해 주어야 한다. 도 7의 순서로 블레이드 설계 및 해석이 진행되며, 최종적인 블레이드 모델은 사용자의 필요에 따라서 CAD 파일 형태로의 엑스포트(Export) 또한 가능하다.
첨부된 도 8은 본 발명을 통해 구현된 블레이드의 모델링 결과이고, 도 9를 통한 구체적인 형상 설정 기능을 통하여 사용자가 원하는 코드길이, 포지션(Position), 트위스트앵글(Twist angle)을 가진 형태의 블레이드를 간편하게 생성할 수 있다. 여기서 구현된 블레이드의 설계 데이터를 ANSYS와 연동하여 3D 모델을 수행하게 되며, 다양한 정보에 대하여 보다 쉽고 빠르게 반복적인 모델링이 가능하다는 장점을 가질 수 있다.
도 10은 본 발명을 통해 구현된 블레이드의 구조 및 유동 해석의 결과이다. 사용자는 해석에 대한 지식이 전무하다 하더라도 자동화된 본 발명의 풍력터빈 블레이드 설계 및 해석 소프트웨어를 통하여 도 10과 같은 결과를 얻어낼 수 있다. 특히, 여러 모델에 대해서 해석의 수행 시 자동화된 모듈을 통하여 그 해석 수행 시간을 비약적으로 단축 시켜 줄 수 있다는 장점을 지닌다.
3. 블레이드 설계의 알고리즘
도 2a는 블레이드의 설계 알고리즘을 나타내고 있다. 도 2a를 참조하면 먼저 블레이드에 대한 설계식을 정립한다. 초기 가정한 동력계수로 모델링 및 해석을 진행하게 되고, 해석 결과의 동력계수로 반복역산을 수행하며, 계산된 동력계수를 본 발명의 풍력터빈 블레이드의 설계 및 해석 소프트웨어에 반영하게 된다. 그 후 설계 알고리즘을 통해 블레이드 설계 알고리즘을 수립하게 된다.
구체적인 블레이드 설계 순서는 도 2b에서 도시된 바와 같다.
블레이드의 정격출력, 블레이드 직경, 주속비 및 회전수를 계산하고 날개 끝 손실개수와 축, 접선 방향의 유도계수를 계산한다. 그리고 코드 길이, 유입각 및 트위스트앵글(비틀림각)을 계산한다. 이어서 블레이드 형상을 모델링을 하고 공력해석 모델을 생성한 뒤, 블레이드 공력해석을 수행(CFD)함으로써, 공력해석 결과와 최초 블레이드 설계치를 비교 검증하여, 수정된 동력계수를 반복 피드백함으로써 블레이드 형상을 결정한다.
본 발명의 풍력터빈 블레이드의 설계 및 해석 소프트웨어에 사용된 블레이드 설계식은 도 3에서 보여진다. 도 3과 같은 코드길이(Chord length), 트위스트 앵글(Twist angle) 등의 설계식을 통하여 자동화된 모델링을 수행하게 된다.
4. 블레이드 해석 알고리즘 모듈 프로그래밍
블레이드의 해석에는 ANSYS WorkBench 라는 상용 유한요소해석 프로그램을 사용할 수 있다. 도 4를 참조하면 ANSYS WorkBench의 경우 전체 프로그램은 Python, 3D 모델링과 구조 해석의 경우 Jscript, 유동 해석의 경우 Perl로 이루어져 있다.
본 발명의 풍력터빈 블레이드 설계 및 해석 소프트웨어인 윈블레이드(WINBLADE)와 ANSYS의 연결 모듈을 구성하기 위하여, ANSYS WorkBench의 내부 구성 언어를 파악하여 블레이드의 모델링과 구조 해석, 유동 해석에 해당하는 부분을 각 스크립트 언어를 통하여 자동화 하였다.
이러한 자동화의 기반을 만들기 위하여 블레이드 선행연구 논문의 레퍼런스 데이터를 이용하였다.
블레이드의 중요 설계 데이터인, 코드길이, 트위스트앵글 및 포지션 데이터(에어포일의 섹션별 위치를 표현하는 데이터임)를 변수처리하여 사용자가 원하는 형상의 블레이드를 자동적으로 구성할 수 있도록 하였다.
또한, 설계된 블레이드에 대하여 구조해석을 수행할 수 있는 FEM 모델을 만드는 작업과 함께 사용자가 원하는 구조하중과 구속조건을 적용할 수 있도록 Python과 Jscript를 통하여 해석 모듈을 구성하였다. 이와 함께, Ptython과 Perl스크립트를 통하여 유동해석을 위한 해석 모듈 또한 구성하였다. 두 모듈의 구동에 필요한 에어포일 데이터(Airfoil Data)는 Naca Series 중 임의의 형상을 통해 구성하였다.
5. 블레이드 모델 구현성능 검증
본 발명의 방법을 블레이드 모델 구현에 대하여 성능 검증을 하기 위하여 상용 프로그램인 ANSYS 와의 비교를 수행하였다.
각각 동일한 에어포일 형상(NACA63415)과 Position, Chord length, Twist angle)에 대하여 블레이드 모델링을 수행하였다.
먼저, ANSYS를 통해 매뉴얼 방식으로 블레이드를 모델링 하였고, 그 결과를 도 11에 나타내었다.
다음으로 본 발명에 따른 방식(WinBlade, 윈블레이드)을 통해 설계에 필요한 데이터를 오토매틱(Automatic) 방식으로 블레이드 모델을 구현하였고 그 결과를 도 12 및 도 13에 나타내었다.
결과를 살펴보면 두 방식 모두 올바른 형태의 블레이드 모델이 구현되었음을 확인할 수 있으며, 이는 본 발명 WinBlade를 통한 블레이드 모델링이 사용성에 있어 타당성이 있음을 확인할 수 있다.
6. 블레이드 출력성능 검증
본 발명에 의한 방법에 대한 블레이드 출력 성능 검증을 위하여 100kW 급 블레이드를 구성하였고, 해석에 사용된 레퍼런스 데이터는 다음 표와 같다.
Rated Power(kW) 100
Rated wind velocity(m/s) 10
Design TSR 6
Rotor diameter(m) 22.65
Rotational speed(rpm) 40.46
Air density(kg/m^3) 1.225
도 14에서 블레이드의 출력이 80.821 kW 됨을 확인할 수 있다. 이는 풍력 터빈의 유동해석의 경우 전이과정(Transient) 해석을 해야 정확한 결과 값을 알 수 있는 것에 비해 본 해석에서는 출력에 대한 검증을 위해 Steady 상태의 해석을 한 것과 더불어 해석 시간에 대한 효율을 높이고자 격자의 숫자를 제한한 것, 그리고 완전한 형상의 블레이드가 아닌 1/3 모델로 구성이 되어 난류의 영향이 제대로 적용되지 않았다는 점을 감안 할 때 80.821kW의 해석 결과는 CFD 자동화 모듈의 신빙성을 충분히 검증 해줄 수 있는 결과로 볼 수 있다.
7. 타워( Tower ), 메인샤프트( Main shaft ) 설계 등
상기에서는 블레이드 설계 해석에 대하여 설명하였지만, 타워 혹은 메인샤프트와 같이 풍력터빈의 주요부품에 대한 해석 작업이 필요할 경우에도 상기와 유사한 방법을 사용하여 상용프로그램과의 호환을 통한 동적 특성 해석 등이 가능하다.
가령 타워의 높이를 디자인 모듈에서 계산이 가능한데, 풍속, 공기밀도, 풍력터빈 효율 등의 데이터를 입력하면 타워의 높이가 계산되며, 타워 높이 계산후 타워 지름과 타워 구성파트 개수등을 입력하여 타워 형상디자인이 가능하다. 이어서 타워에 대한 공력 특성 등을 해석하기 위하여 ANSYS 프로그램을 이용할 수 있다. 도 15는 타워의 형상 결정후 CFD 해석과 구조해석이 수행된 결과 화면에 해당한다.
기타 본 발명에 의한 풍력터빈 해석후 각종 해석 데이터 자료 즉, 개념설계, 모델링, 격자 정보, 하중조건, 응력평가 결과 등의 후처리 결과를 자동으로 보고서 형태로 생성해 주는 리포트 제너레이터 모듈 등도 추가로 구현가능하다. 게다가 풍력터빈의 경제적 평가 예를 들어, 풍속에 따른 발전량 등을 계산하여 경제성을 고려한 평가를 할 수도 있을 것이다.
리포트 모듈에서는 블레이드 등의 공력해석 결과 및 블레이드 형상에 대한 결과 보고서가 출력되도록 구성될 수 있다.

Claims (1)

  1. 블레이드 형상설계용 공개프로그램과 연동됨과 동시에 블레이드형상을 결정하는 디자인모듈(Design Module), 블레이드 해석용 상용 프로그램과의 데이터 교환을 위한 연결모듈(Connection Module), 디자인모듈로부터 설계된 블레이드 형상에 대한 공력해석을 수행하는 블레이드 해석용 상용 프로그램을 준비해 두는 단계(s1);
    상기 블레이드 형상설계용 공개프로그램의 경로를 미리 설정하여 상기 연결모듈에 의하여 데이터 교환이 가능하도록 하는 단계(s2);

    상기 디자인 모듈에서,
    입력된 에어포일(Airfoil)의 형상 이름과 레이놀즈수(Reynolds Number), 마하수(Mach Number), 유체밀도(Fluid Density), 동점성계수(Kinetic Viscosity) 데이터가 블레이드 형상설계용 공개프로그램으로 전송되는 단계(s3);

    블레이드 형상설계용 공개프로그램에서, 상승력계수(CL), 저항계수(CD), 코드길이(Chord length), 포지션(Position), 트위스트앵글(Twist angle)값이 계산된 후, 다시 디자인모듈로 전송되는 단계(s4);

    상기 디자인 모듈로 전송된, 상승력계수와 저항계수와 코드길이와 포지션데이터(에어포일의 섹션별 위치를 표현하는 데이터)과 트위스트앵글값을 이용하여, 블레이드 직경, 주속비 및 회전수를 계산하고, 날개 끝 손실개수와 축, 접선 방향의 유도계수를 계산하며 코드 길이, 유입각 및 트위스트앵글(비틀림각)을 계산하여 블레이드 형상이 모델링 되는 단계(s5);

    상기 모델링된 블레이드 형상데이터가 연결모듈에 의하여, 블레이드 해석용 상용 프로그램으로 전송되는 단계(s6);
    전송된 블레이드 형상에 대하여, 공력해석 모델을 생성한 뒤, 블레이드 공력해석을 수행(CFD)되어, 공력해석 결과와 최초 블레이드 설계치를 비교 검증하여, 수정된 동력계수를 피드백함으로써 블레이드 형상을 결정하는 단계(s7);

    상기 s7 단계의 블레이드 공력해석결과 및 결정된 블레이드 형상에 대한 보고서가 출력되는 단계(s8); 을 포함하는 풍력터빈 설계 및 해석방법
KR1020140082993A 2014-07-03 2014-07-03 풍력터빈 설계 및 해석방법 KR20160007838A (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020140082993A KR20160007838A (ko) 2014-07-03 2014-07-03 풍력터빈 설계 및 해석방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020140082993A KR20160007838A (ko) 2014-07-03 2014-07-03 풍력터빈 설계 및 해석방법

Publications (1)

Publication Number Publication Date
KR20160007838A true KR20160007838A (ko) 2016-01-21

Family

ID=55308393

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020140082993A KR20160007838A (ko) 2014-07-03 2014-07-03 풍력터빈 설계 및 해석방법

Country Status (1)

Country Link
KR (1) KR20160007838A (ko)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109066772A (zh) * 2018-08-01 2018-12-21 广州供电局有限公司 风电爬坡率的平抑方法、装置和系统
CN110375947A (zh) * 2019-07-05 2019-10-25 浙江大学 一种输电塔体型系数的现场实测装置和测试方法
KR20200021768A (ko) 2018-08-21 2020-03-02 두산중공업 주식회사 유동 해석을 최적화하기 위한 장치 및 이를 위한 방법
KR20200021766A (ko) 2018-08-21 2020-03-02 두산중공업 주식회사 유동 해석 장치 및 이를 위한 방법
KR20200021767A (ko) 2018-08-21 2020-03-02 두산중공업 주식회사 학습된 모델을 이용한 해석 장치 및 이를 위한 방법
CN113945357A (zh) * 2021-10-28 2022-01-18 南京航空航天大学 一种风力机叶片颤振临界风速预测方法及装置
KR20230103843A (ko) 2021-12-30 2023-07-07 영남대학교 산학협력단 유동 해석 방법 및 이를 수행하기 위한 컴퓨팅 장치
KR20230145871A (ko) 2022-04-11 2023-10-18 영남대학교 산학협력단 유동 해석 방법 및 이를 수행하기 위한 컴퓨팅 장치
CN117540613A (zh) * 2024-01-10 2024-02-09 西安羚控电子科技有限公司 一种气动载荷二次插值处理方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
없음.

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109066772A (zh) * 2018-08-01 2018-12-21 广州供电局有限公司 风电爬坡率的平抑方法、装置和系统
US11416654B2 (en) 2018-08-21 2022-08-16 Dosan Enerbility Co., Ltd. Analysis apparatus using learned model and method therefor
KR20200021768A (ko) 2018-08-21 2020-03-02 두산중공업 주식회사 유동 해석을 최적화하기 위한 장치 및 이를 위한 방법
KR20200021766A (ko) 2018-08-21 2020-03-02 두산중공업 주식회사 유동 해석 장치 및 이를 위한 방법
KR20200021767A (ko) 2018-08-21 2020-03-02 두산중공업 주식회사 학습된 모델을 이용한 해석 장치 및 이를 위한 방법
US11531795B2 (en) 2018-08-21 2022-12-20 Dosan Enerbility Co., Ltd. Apparatus for optimizing flow analysis and method therefor
US11544428B2 (en) 2018-08-21 2023-01-03 Dosan Enerbility Co., Ltd. Flow analysis apparatus and method therefor
CN110375947A (zh) * 2019-07-05 2019-10-25 浙江大学 一种输电塔体型系数的现场实测装置和测试方法
CN113945357A (zh) * 2021-10-28 2022-01-18 南京航空航天大学 一种风力机叶片颤振临界风速预测方法及装置
KR20230103843A (ko) 2021-12-30 2023-07-07 영남대학교 산학협력단 유동 해석 방법 및 이를 수행하기 위한 컴퓨팅 장치
KR20230145871A (ko) 2022-04-11 2023-10-18 영남대학교 산학협력단 유동 해석 방법 및 이를 수행하기 위한 컴퓨팅 장치
CN117540613A (zh) * 2024-01-10 2024-02-09 西安羚控电子科技有限公司 一种气动载荷二次插值处理方法及装置
CN117540613B (zh) * 2024-01-10 2024-04-05 西安羚控电子科技有限公司 一种气动载荷二次插值处理方法及装置

Similar Documents

Publication Publication Date Title
KR101418445B1 (ko) 풍력터빈 블레이드 설계 및 해석방법
KR20160007838A (ko) 풍력터빈 설계 및 해석방법
Bonaiuti et al. On the coupling of inverse design and optimization techniques for the multiobjective, multipoint design of turbomachinery blades
Marten et al. Qblade guidelines
Kim et al. Developement and verification of a performance based optimal design software for wind turbine blades
CN101923467A (zh) Cad软件、ug软件与cae软件之间的模型传输方法
Rubechini et al. A redesign strategy to improve the efficiency of a 17-stage steam turbine
Chaudhary et al. Modeling and optimal design of small HAWT blades for analyzing the starting torque behavior
CN106354967A (zh) 一种涡轮盘结构设计分析一体化方法
He et al. Three-dimensional aerodynamic optimization for axial-flow compressors based on the inverse design and the aerodynamic parameters
CN113742861A (zh) 一种适用于风力发电机风洞试验的叶片模型优化设计方法
CN117634091A (zh) 一种基于OpenFOAM软件的水轮机数值模拟平台搭建方法及系统
Safari et al. A novel combination of adaptive tools for turbomachinery airfoil shape optimization using a real-coded genetic algorithm
Zangeneh et al. 3D multi-disciplinary inverse design based optimization of a centrifugal compressor impeller
Sairam The influence of radial area variation on wind turbines to the axial induction factor
Krumme Performance prediction and early design code for axial turbines and its application in research and predesign
Bonnet et al. Parametric modelling of large wind turbine blades
Ouellet et al. A preliminary design system for turbine discs
Obaida et al. Loss reduction in a 1.5 stage axial turbine by computer-driven stator hub contouring
Reising Steady and unsteady performance of a transonic compressor stage with non-axisymmetric end walls
CN117454721B (zh) 基于数字仿真实验的风电场尾流叠加效应评估方法和介质
Kadhim Effect of non-axisymmetric casing on flow and performance of an axial turbine
CN109766584B (zh) 一种自定义式发电机控制器电磁暂态自动建模方法及系统
Strofylas et al. Wind turbine blade structure parameterization using T4T
Abbasi et al. Effects of different turbulence models in simulation of unsteady tip leakage flow in axial compressor rotor blades row

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E601 Decision to refuse application