風力機程序設(shè)計報告

1、風力機程序設(shè)計報告翼型NACA4412指導老師： X X X 專 業(yè)：XXXXXXXXXXXXXXXX姓 名： X X X 學 號： XXXXXXXXXXXX 一 目的：1.根據(jù)課堂所學知識計算風輪設(shè)計直徑、轉(zhuǎn)速以及軸向和徑向誘導因子2.繪制風力機葉片弦長隨葉片展向長度的變化曲線3.繪制風力機扭角隨葉片展向長度的變化曲線二 1.給定原始參數(shù)：計功率P：7.6KW（由6KW+學號16*0.1KW）葉片數(shù)Z：3設(shè)計風速V：6m/s功率系數(shù)Cp：0.4給定最佳葉尖速比：5 空氣的密度p：1.225Kg/m3由給定翼型升阻比參數(shù)，得出最佳攻角為8度2.需要計算的參數(shù)：風輪的直徑D風輪葉片的轉(zhuǎn)速n各半徑

2、處葉素處的葉尖速比各半徑處葉素處的切向速度Vw各半徑處葉素處的合成速度Vr升力系數(shù)Cl，阻力系數(shù)Cd，最佳攻角，各半徑處葉素處的安裝角，各半徑處葉素處的入流角，各半徑處葉素處的扭角，各半徑處葉素處的軸向誘導因子a，周向誘導因子b，各半徑處葉素的弦長t四設(shè)計的原理與步驟：1，確定風輪的直徑D：計算風輪葉片的公式為：D=8PCpV3，其中P-為風輪的額定功率，Cp-為風能利用系數(shù)，-為空氣密度，V-為設(shè)計風速2.確定風輪的轉(zhuǎn)速n：記算風輪葉片的計算公式為：n=60*v*D，其中V-為設(shè)計風速，-為葉尖速比，D-為葉輪的設(shè)計直徑3.確定風輪各半徑處尖速比 10%=×10%V1 20%=&#

3、215;20%V1 30%=×30%V1 40%=×40%V1 50%=×50%V1 60%=×60%V1 70%=×70%V1 80%=×80%V1 90%=×90%V14.計算截面的弦長葉片半徑10%截面處弦長；葉片半徑20%截面處弦長；葉片半徑30%截面處弦長；葉片半徑40%截面處弦長；葉片半徑50%截面處弦長；葉片半徑60%截面處弦長；葉片半徑70%截面處弦長；葉片半徑80%截面處弦長；葉片半徑90%截面處弦長；5，翼型弦長的計算各葉素弦長公式其中：葉片數(shù)；：最大升阻比是的升力系數(shù)；：半徑為r處的葉素旋轉(zhuǎn)速度；：葉尖

4、速比。6.計算安裝角及扭角根據(jù)相關(guān)關(guān)系式可通過迭代方法求得軸向誘導因子a和周向誘導因子b，從而計算出安裝角的展向分布，具體迭代步驟如下：假設(shè)a和b的初值，一般可取0；（1）計算入流角；（2）計算安裝角 ；（3）根據(jù)翼型空氣動力特性曲線得到葉素的升力系數(shù)和阻力系數(shù)；（4）計算葉素的法向力系數(shù)和切向力系數(shù)；（5）計算a和b的新值其中，為葉素實度（6）比較新計算的a和b值與上一次的a和b值，如果誤差小于設(shè)定的誤差值（可取0.001），則迭代終止；否則，再回到（1）繼續(xù)迭代，實際過程中限于時間考慮會同時施加較大的步數(shù)限制。五設(shè)計結(jié)果1.程序界面2扭角隨葉片展向的變化3.弦長隨葉片展向的變化六，總結(jié) 經(jīng)

5、過這次的課程設(shè)計，我有很多的心得體會。由于VB程序設(shè)計是大一時學習的，關(guān)于編程的一些步驟以及程序中的一些字符、方法也忘記了不少，只能找出VB書重新學習一遍，好在先前的基礎(chǔ)，除了一些需要新學的功能（如畫圖），其他基本知識能很快記起來。正好趁著本門課程考試，這個程序中所有關(guān)于書上的知識點都連帶著一起復(fù)習了，這就為課程設(shè)計打好了知識基礎(chǔ)。由于初次進行課程設(shè)計，過程中存在一定的不足，一些公式由于太過繁瑣，報告就省略了部分輸入，另外在代碼中部分變量命名或者排列存在一定的缺陷，望老師見諒。一份努力，一分收獲，相信以后自己能夠做的更好。七.代碼'通用變量Const pi = 3.1415926Dim

6、 t(1 To 10) As Single '弦長Dim bt(1 To 10) As Single '扭角 Private Sub Command1_Click() List1.clear:list2.clear:list3.clear:list4.clearDim r1 As Single, i As Integer, vr As Single Dim cn As Single, ct As Single Dim difa As Single, difb As Single Dim a1 As Single, b1 As Single Dim z As Integer &#

7、39;葉片數(shù) Dim jsd As Single '角速度 Dim sd(1 To 10) As Single '實度 Dim r As Single, jsb As Single '半徑和尖速比 Dim fai As Single, afa As Single '入流角、攻角 Dim a As Double, b As Double 'a軸向誘導因子,b周向誘導因子 Dim As Single, p As Single, v As String, cp As Single, power As Single '葉尖速比;p 密度;v 來流風速;c

8、p 功率系數(shù);power 功率 afa = 0.14 z = 3 cl = 1.2424 cd = 0.0222 power = Val(textpower.Text) = Val(text.Text) p = Val(textp.Text) v = Val(textv.Text) cp = Val(textcp.Text) '賦值 r = 0.5 * Sqr(8 * power / (3.14 * p * cp * v 3) '葉片半徑 n = (60 * * v) / (pi * 2 * r) '轉(zhuǎn)速 Textd.Text = Format(2 * r, "

9、;0.0000") Textn.Text = Format(n, "0.0000") For i = 1 To 10 r1 = 0.1 * i * r jsb = n * 3.14 * 2 * r1 / (60 * v) jsd = jsb * v / r1 vr = Sqr(v * v + jsd * r1 * jsd * r1) t(i) = (8 * 3.14 * r1 / (3 * z * cl) * (v / (jsb * vr) Form1.List1.AddItem m & "0%半徑處" & Space(5) &

10、amp; Format(t(i), "0.0000") sd(i) = t(i) * z / (2 * 3.14 * r) a = b = 0 a1 = b1 = 0 Do fai = Atn(v * (1 - a) / (jsd * r1 * (1 + b) cn = cl * Cos(fai) + cd * Sin(fai) ct = cl * Sin(fai) - cd * Cos(fai) a1 = 1 / (4 * Sin(fai) * Sin(fai) / (sd(i) * cn) + 1) b1 = 1 / (4 * Sin(fai) * Cos(fai) /

11、 (sd(i) * ct) - 1) difa = Abs(a - a1): difb = Abs(b - b1) a = a1: b = b1 Loop While Abs(a - a1) > 0.00001 And Abs(b - b1) > 0.00001 Form1.List3.AddItem i & "0%半徑處" & Space(7) & Format(a, "0.0000") Form1.List4.AddItem m & "0%半徑處" & Space(7) &

12、; Format(b, "0.0000") fai = Atn(v * (1 - a) / (jsd * r1 * (1 + b) bt(i) = (fai - afa) * 180 / 3.14 Form1.List2.AddItem i & "0%半徑處" & Space(7) & Format(bt(i), "0.000") Next Command2.Enabled = True Command5.Enabled = True End Sub Private Sub Command5_Click()Lo

13、ad Form2Form2.ShowForm2.Pic1.Scale (-20, 55)-(120, -10)For i = 0 To 50 Step 5 Form2.Pic1.Line (0, i)-(100, i), vbBlack Form2.Pic1.CurrentX = -5 Form2.Pic1.CurrentY = i If i > 0 Then Form2.Pic1.Print i End IfNext iFor i = 0 To 100 Step 10Form2.Pic1.CurrentX = i - 5Form2.Pic1.CurrentY = -1Form2.Pic

14、1.Print i; "%"Next iForm2.Pic1.CurrentX = 21: Form2.Pic1.CurrentY = -5Form2.Pic1.Print "葉弦隨展向長度變化"Form2.Pic1.CurrentX = -10: Form2.Pic1.CurrentY = 50Form2.Pic1.Print "扭角/度"For i = 0 To 100 Step 10Form2.Pic1.Line (i, 0)-(i, 50), vbBlackNext iForm2.Pic1.DrawWidth = 2Form2

15、.Pic1.Line (0, 0)-(110, 0)Form2.Pic1.Line (0, 0)-(0, 55)Form2.Pic1.CurrentX = 10: Form2.Pic1.CurrentY = bt(1)For i = 1 To 10Form2.Pic1.Line -(10 * i, bt(i), vbYellowNext iFor i = 1 To 8Form2.Pic1.CurrentX = 10 * i: Form2.Pic1.CurrentY = bt(i) + 2Form2.Pic1.Print Format(bt(i), "0.0000")Next

16、 iEnd SubPrivate Sub Command2_Click()Dim m As IntegerLoad Form3Form3.ShowForm3.Pic1.Scale (-20, 200)-(120, -20)For i = 0 To 200 Step 10Form3.Pic1.Line (0, i)-(100, i), vbBlackNext iFor i = 0 To 100 Step 10Form3.Pic1.CurrentX = i - 5Form3.Pic1.CurrentY = -1Form3.Pic1.Print i; "%"Next iForm3

17、.Pic1.CurrentX = 21: Form3.Pic1.CurrentY = -10Form3.Pic1.Print "葉弦隨展向長度變化"Form3.Pic1.CurrentX = -10: Form3.Pic1.CurrentY = 200Form3.Pic1.Print "弦長m"For i = 0 To 100 Step 10 Form3.Pic1.Line (i, 0)-(i, 200), vbBlackNext iFor i = 0 To 200 Step 10Form3.Pic1.CurrentX = -7Form3.Pic1.CurrentY = iIf i > 0 ThenForm3.Pic1.Print (i / 50)End IfNext iForm3.Pic1.DrawWidth = 2Form3.Pic1.Line (0, 0)-(110, 0)Form3.Pic1.Line (0, 0)-(0, 200)For