為瞭解風機葉片在各種風況下的受力情形及安全性，本文使用有限元素分析軟體ANSYS中的CFX流場分析模組，計算葉片壓力分佈的情況，並把此結果匯入有限元素分析軟體中，進行結構分析，將這些壓力分佈結果以分佈力的形式映射到葉片的有限元素模型，作為負載條件，以蔡吳理論(Tsai-Wu failure criterion)來判斷葉片是否有遭受到破壞。本文針對4支葉片進行流場及結構分析，首先進行V47及V80葉片之分析，探討葉片的破損原因，並將V47葉片之分析結果與實驗之數值進行比對，驗證模擬之正確性。接著依據葉片設計公式及NREL公布之葉片資訊，設計出3.6MW及5MW葉片幾何外型，改變其葉片的複材疊層，探討分析結果，並比較不同材料之成本、強度及重量。

In order to understand the load distributions and the safety of large-scale wind turbine blades under various wind conditions, ANSYS finite element analysis (FEA) software package was often used for investigation. The pressure distribution on the blades was calculated by use of ANSYS CFX flow field analysis module . The resulting pressure distribution was mapped on to the finite element model of the blades using distributed forces as the loading conditions. After finite element structural analysis, the Tsai-Wu failure criterion was used for damage prediction. In the present study, the flow fields and structures for four types of turbine blades were investigated under the aforementioned methodology. First, the V47 and V80 blades were studied for damage prediction. The results of V47 blade analysis were compared against experimental data for verification of the accuracy. Second, based on the blade design formulas and the turbine blade data released by the National Renewable Energy Laboratory, the 3.6-MW and 5-MW blades were designed, plotted and analyzed. With various types of lamination, the costs, strengths, and weights of the proposed laminate composite blades are studied and compared.

誌謝 ii
中文摘要 iii
ABSTRACT iv
目錄 v
圖目錄 vii
表目錄 x
第1章 緒論 1
1.1 研究動機 1
1.2 文獻回顧 1
1.3 研究方法 3
1.4 論文架構 4
第2章 葉片模型簡介 5
2.1 Vestas V47機組 5
2.2 Vestas V80機組 10
第3章 單向流固耦合數值計算 13
3.1 流場及結構分析邊界設定 15
3.2 流固耦合數值計算結果 17
3.2.1 Vestas V47機組 17
3.2.2 Vestas V80機組 22
第4章 試驗架設與模擬驗證 28
4.1 流場模擬之輸出功率驗證 28
4.1.1 Vestas V47機組 28
4.1.2 Vestas V80機組 28
4.2 葉片模型與試驗自然模態驗證 29
4.2.1 Vestas V47機組 29
4.2.2 Vestas V80機組 31
4.3 葉片試驗與模擬結果驗證 33
4.3.1 V47葉片在IEC（International Electrotechnical Commission）標準下之停機負荷分析 33
4.3.2 V47葉片靜態試驗架設 36
4.3.3 模擬驗證與試驗結果 41
第5章 大型風機葉片外型設計與分析 45
5.1 3.6MW風機機組 45
5.1.1 葉片外型設計及繪製 45
5.1.2 分析設定與分析結果 56
5.2 5MW風機機組 59
5.2.1 葉片外型設計及繪製 59
5.2.2 分析設定與分析結果 61
第6章 結論與未來展望 64
參考文獻 66
附錄A 葉片疊層資訊 69
附錄B 3.6MW葉片蔡吳指標分佈圖 87
附錄C 5MW葉片蔡吳指標分佈圖 93

[1]黃慶民、蔡國忠、陳世雄、王晟桓、李啟勝, "50kW風力發電機複合材料葉片設計、分析與製作", 中華民國力學學會第三十一屆全國力學會議, 2007.
[2]鄭榮和、林家緯、李盈宏、沈丞佑, "台電中彰風場風力特性評估與大型風力機組技術研究期末報告-TIS-1風機研究子計畫三：風機結構分析期末報告", 台灣電力股份有限公司, 2012.
[3]M. Moshfeghi, Y. Jun Song, Y. Hui Xie, "Effects of near-wall grid spacing on SST-K-ω model using NREL Phase VI horizontal axis wind turbine", Journal of Wind Engineering and Industrial Aerodynamics, vol. 107–108, p. 94-105, August–September 2012.
[4]陳彥呈, "陣風效應對風力發電機在停機狀態下之氣動力數值模擬研究", 國立臺灣大學工學院工程科學及海洋工程學系碩士論文, 2011.
[5]李文傑, "運轉狀態下風力發電機之氣動力負荷數值研究", 國立台灣科技大學機械工程系碩士論文, 2013.
[6]J.C. Marin, A. Barroso, F. Paris, J. Canas, "Study of damage and repair of blades of a 300 kW wind turbine",Energy, vol. 33, p.1068-1083, July 2008.
[7]T. Uchida, T. Maruyama, H. Ishikawa, A. Deguchi, "Investigation of the Causes of Wind Turbine Blade Damage at Shiratakiyama Wind Farm in Japan - A Computer Simulation Based Approach", Reports of Research Institute for Applied Mechanics, Kyushu University, vol..141, p.13 – 25, 2011.
[8]Y.M. Ji, K.S. Han, "Fracture mechanics approach for failure of adhesive joints in wind turbine blades", Renewable Energy, vol. 65, p. 23-28, May 2014.
[9]劉哲元, "FRP 風力機葉片材料強度之探討", 國立臺灣大學工學院工程科學及海洋工程學系碩士論文, 2008.
[10]Vestas, "Vestas 660 kW Variable Slip Wind Turbine", 2000.
[11]林輝政、鄭錦榮、李雅榮、鍾承憲、陳凱琳、劉二瑋、陳金汎, "台電中彰風場風力特性評估與大型風力機組技術研究子計畫六：風機葉片之逆向工程與修護評估技術之建立期末報告", 台灣電力股份有限公司, 2012.
[12]FARO, "http://www.faro.com/home", 2014.
[13]C. Konga, J. Banga, Y. Sugiyama, "Structural investigation of composite wind turbine blade considering various load cases and fatigue life", Energy, p. 2101–2114, 2005.
[14]J.C. Marin, A. Barroso, F. Paris, J. Canas, "Study of fatigue damage in wind turbine blades", Engineering Failure Analysis, vol. 16, p. 656-668, March 2009.
[15]鍾承憲、吳兆誠、林勁成、何智南, "葉片用單向玻纖材料試驗規範與程序探", 風能研討會論文, 2008.
[16]Vestas, "http://www.vestas.com", 2014.
[17]F.K Benra, H. J. Dohmen, J. Pei, S. Schuster, and B. Wan, "A Comparison of One-Way and Two-Way Coupling Methods for Numerical Analysis of Fluid-Structure Interactions", Journal of Applied Mathematics, 2011.
[18]郭哲漁, "風力發電葉片結構改良─扭力桿裝置之研究", 國立臺灣大學工程科學及海洋工程學系碩士論文, 2012.
[19]林輝政、鍾承憲、陳凱琳、鄭錦榮, "修補後V47整支葉片實體靜態分析與觀園#12 GE-1.5SE風機葉片破損診斷分析期末報告", 台灣電力股份有限公司, 2013.
[20]International Electrotechnical Commission, "IEC 61400-1: Wind turbines –Part 1: Design Requirements", 2005.
[21]M. M. Shokrieh, and R. Rafiee, "Simulation of fatigue failure in a full composite wind turbine blade", Composite Structures, vol. 74, p. 332-342, 2006.
[22]P. J. Schubel and R. J. Crossley "Wind Turbine Blade Design", Energies, p.3425-3449, 2012.
[23]P. Fuglsang, "Wind Tunnel Tests of the FFA-W3-241", Ris National Laboratory, Roskilde, Denmark, December 1998.
[24]F. Bertagnolio, N. Srensen, J. Johansen and P. Fuglsang, "Wind Turbine Airfoil Catalogue", Ris National Laboratory, Roskilde, Denmark, August 2001.
[25]Airfoil Investigation Database, "http://www.airfoildb.com/foils/340", 2014.
[26]siemens, "Wind Turbine SWT-3.6-107", http://www.energy.siemens.com", 2014.
[27]siemens, "Wind Turbine SWT-3.6-120", http://www.energy.siemens.com", 2014.
[28]牛山泉,林輝政, "風車工學入門", 國立澎湖科技大學, 2009.
[29]吳文翔, "複材風力發電機葉片設計與分析", 國立成功大學航空太空工程學系碩士班論文, 2010.
[30]林昆瑩, "小型風力發電機葉片設計及製作", 國立成功大學航空太空工程學系碩士班論文, 2009.
[31]Dong energy, "http://www.dongenergy.com/"
[32]P. Brondsted , J. W. Holmes, B. F. Sorensen, Z. Jiang, Z. Sun,and X. Chen, "Evaluation of a Bamboo/Epoxy Composite as a Potential Material for Hybrid Wind Turbine Blades", Ris&;#248; National Laboratory for Sustainable Energy, Technical University of Denmark, Roskilde, Denmark and International Centre for Bamboo and Rattan, Beijing, China, 2008.
[33]J. C. O. Lopes, "Material selection for aeronautical structural application", Ci&;#234;ncia &; Tecnologia dos Materiais, vol. 20, 2008.
[34]J. Jonkman, S. Butterfield, W. Musial,and G. Scott, "Definition of a 5MW Reference Wind Turbine for Offshore System Development", Technical Report NREL/TP-500-38060, February 2009.