Projekt turbiny akcyjnej

Dane:

ṁ= 637 $\frac{\text{kg}}{s}$


P0 = 9, 2 MPa


T0 = 422 C


dp = 0, 943 m


$$c_{0} = 48\ \frac{m}{s}$$

ρd = 0, 25 m


$$\frac{u}{c_{1}} = 0,526$$

n= 3000 min−1

ϕ= 0,96

ψ= 0,89


μ1 = 0, 98


μ2 = 0, 97

Κ= 1,3

ε= 1


τ1 = τ2 = 0, 9


α1 = 14

1. Obliczam prędkość obwodową:


$$u = \frac{\pi*d_{p}*n}{60} = \frac{\pi*0,943*3000}{60} = 148,1\frac{m}{s}$$

2. Obliczam energię kinetyczną początkową Δh0:


$$h_{0} = \frac{c_{0}^{2}}{2} = \frac{48^{2}}{2} = 1,15\frac{\text{kJ}}{\text{kg}}$$

3. Obliczam rzeczywistą prędkość bezwzględną c1:


$$c_{1} = \frac{u}{0,526} = \frac{148,1}{0,526} = 281,5\frac{m}{s}$$

4. Obliczam teoretyczną prędkość bezwzględną c1s:


c1 = c1s * φ


$$c_{1s} = \frac{c_{1}}{\varphi} = \frac{281,5}{0,96} = 293,2\frac{m}{s}$$

5. Obliczam spadek entalpii na kierownicy:


$$c_{1s} = \sqrt{2H^{k} + c_{0}^{2}}$$


$$H_{k} = \frac{c_{1}^{2}}{2*\varphi} - \frac{c_{0}^{2}}{2} = \frac{{281,5}^{2}}{2*0,96} - \frac{48^{2}}{2} = 40,1\frac{\text{kJ}}{\text{kg}}$$

6. Obliczam spadek entalpii w całym stopniu:


$$H = \frac{H^{k}}{(1 - \rho_{d})} = \frac{40,1}{(1 - 0,25)} = 53,5\frac{\text{kJ}}{\text{kg}}$$

7. Obliczam spadek entalpii na wirniku:


$$H^{w} = H - H^{k} = 53,5 - 40,1 = 13,4\frac{\text{kJ}}{\text{kg}}$$

8. Obliczam prędkość względną w1dla α1 = 14:

$w_{1} = \sqrt{u^{2} + c_{1}^{2} - 2*u*c_{1}*cos\alpha_{1}}$=$\sqrt{{148,1}^{2} + {281,5}^{2} - 2*148,1*285,1*cos14}$=

=142,4 $\frac{m}{s}$

9. Obliczam teoretyczną prędkość względną w2s:

$w_{2s} = \sqrt{2H^{w} + w_{1}^{2}\ }$=$\sqrt{2*13,4*10^{3} + {142,4}^{2}}$= 216,9 $\frac{m}{s}$

10. Obliczam prędkość względną w2:


$$w_{2} = \psi*w_{2s} = 0,89*216,9 = 193,04\frac{m}{s}$$

11. Obliczam kąt β1 wpływu do kanału wirnikowego:


$$\beta_{1} = arctg\frac{\sin\alpha_{1}}{\cos\alpha_{1} - \frac{u}{c_{1}}} = arctg\frac{sin14}{cos14 - 0,526} = 29$$

12. Odczytuje z wykresu i-s entalpie i objętość właściwą:


$$i_{0} = 3180\frac{\text{kJ}}{\text{kg}}$$


$$\nu_{0} = 0,031\ \frac{m^{3}}{\text{kg}}$$

13. Obliczam entalpię punktu 1s:


$$i_{1s} = i_{0} - H^{k} = 3180 - 40,1 = 3199,9\frac{\text{kJ}}{\text{kg}}$$

14. Odczytuję z wykresu i-s parametry dla punktu 1s:


$$\nu_{1s} = 0,035\ \frac{m^{3}}{s}$$


P1 = 7, 93 MPa


T1 = 399C

15. Obliczam entalpie dla punktu 2s:


$$i_{2s} = i_{0} - H = 3180 - 53,5 = 3126,5\frac{\text{kJ}}{\text{kg}}$$

16. Odczytuję z wykresu i-s parametry dla punktu 2s:


$$\nu_{2s} = 0,036\ \frac{m^{3}}{s}$$


P2 = 7, 57 MPa


T2 = 391C

17. Obliczam pole przepływu pary przez łopatkę kierowniczą A1:


$$A_{1} = \frac{m*\nu_{1s}}{c_{1}} = \frac{637*0,035}{281,5} = 0,0792\ m^{2}$$

18. Obliczam wysokość łopatki kierowniczej l1:


$$l_{1} = \frac{A_{1}}{\pi*d_{p}*\varepsilon*\tau_{1}*\mu_{1}*sin\alpha_{1}} = \frac{0,0792}{\pi*0,943*1*0,9*0,98*sin14} = 0,125\ m$$

19. Obliczam pole przepływu pary przez łopatkę wirnika A2:


$$A_{2} = \frac{m*\nu_{2s}}{w_{2}} = \frac{637*0,036}{193,04} = 0,119\ m^{2}$$

20. Obliczam wysokość łopatki wirnikowej l2:

lz – długość przykrycia zewnętrznego w łopatce wirnikowej (przyjmuję lz = 0, 002 m);

lw – długość przykrycia wewnętrznego w łopatce wirnikowej (przyjmuję lw = 0, 002 m)


lz = l1 + (lz+lw) = 0, 125 + 0, 004 = 0, 129 m

21. Obliczam kąt wypływu z kanału wirnikowego β2:


$$\beta_{2} = arcsin\frac{A_{2}}{\pi*d_{p}*l_{2}*\mu_{2}*\tau_{2}*\varepsilon} = arcsin\frac{0,119}{\pi*0,943*0,129*0,97*0,9*1} = 21$$

22. Obliczam rzeczywistą prędkość bezwzględną w wirniku c2:

$c_{2} = \sqrt{u^{2} + w_{2}^{2} - 2*u*w_{2}*cos\beta_{2}}$=$\sqrt{{148,1}^{2} + {193,04}^{2} - 2*148,1*193,04*cos21}$=

=76,3 $\frac{m}{s}$

23. Obliczam kąt wylotowy α2:


$$\alpha_{2} = arctg\frac{\sin\beta_{2}}{\cos\beta_{2} - \frac{u}{w_{2}}} = arctg\frac{sin21}{cos21 - \frac{148,1}{198,04}} = 65$$

24. Obliczam stratę wylotową hwyl:


$$h_{\text{wyl}} = \frac{c_{2}^{2}}{2} = \frac{{76,3}^{2}}{2} = 2,91\frac{\text{kJ}}{\text{kg}}$$

25. Obliczam straty w dyszy i na łopatce wirnika:


$$h_{d} = \frac{1}{2}c_{1}^{2}*\left( \frac{1}{\varphi^{2}} - 1 \right) = 3,37\frac{\text{kJ}}{\text{kg}}$$


$$h_{\tau} = \frac{1}{2}w_{2}^{2}*\left( \frac{1}{\psi^{2}} - 1 \right) = 4,89\frac{\text{kJ}}{\text{kg}}$$

26. Obliczam pracę obwodową hu:


hu = (H+h0) − (hd+hτ+hwyl) = (53,5+1,15) − (3,37+4,89+2,96)=

$= 43,43\frac{\text{kJ}}{\text{kg}}$

27. Obliczam sprawność obwodową ηu:


$$\eta_{d} = \frac{h_{u}}{H + h_{0}} = \frac{43,43}{53,5 + 1,15} = 0,79$$

28. Obliczam sprawność kierowniczą ηd:


ηd = φ2 = 0, 962 = 0, 92

29. Obliczam sprawność wirnikową ηw:


ηw = ψ2 = 0, 892 = 0, 79

30. Obliczam Lu:

c1u = cosα1 * c1 = cos14 * 281, 5=273,1 $\frac{m}{s}$

c2u = cosα2 * c2 = cos65 * 76, 3=32,2 $\frac{m}{s}$


$$L_{u} = u*\left( c_{1u} + c_{2u} \right) = 148,1*\left( 273,1 + 32,2 \right) = 45,21\frac{\text{kJ}}{\text{kg}}$$

31. Obliczam prędkość dźwięku na wylocie z kierownicy a1s:

$a_{1s} = \sqrt{\kappa*\nu_{1s}*P_{1}}$=$\sqrt{1,3*0,035*7,93*10^{6}}$=601 $\frac{m}{s}$

32. Obliczam liczbę Macha na wylocie z kierownicy Ma1:


$$M_{a1} = \frac{c_{1s}}{a_{1s}} = \frac{293,2}{601} = 0,48$$

33. Obliczam prędkość dźwięku na wylocie z wirnika w2s:

$a_{2s} = \sqrt{\kappa*\nu_{2s}*P_{2}}$=$\sqrt{1,3*0,036*7,57*10^{6}}$= 595 $\frac{m}{s}$

34. Obliczam liczbę Macha na wylocie z wirnika Ma2:


$$M_{a2} = \frac{w_{2s}}{a_{2s}} = \frac{216,9}{595} = 0,36$$

35. Obliczam składową obwodową prędkości względnej wlotowej w1u:


$$w_{1u} = w_{1}*cos\beta_{1} = 142,4*cos29 = 124,5\frac{m}{s}$$

36. Obliczam składową obwodową prędkości względnej wylotowej w2u:


$$w_{2u} = w_{2}*cos\beta_{2} = 193,04*cos21 = 180,2\frac{m}{s}$$

37. Obliczam składową promieniową prędkości względnej wlotowej w1a:


$$w_{1a} = w_{1}*sin\beta_{1} = 142,4*sin29 = 69,04\frac{m}{s}$$

38. Obliczam składową promieniową prędkości względnej wylotowej w2a:


$$w_{2u} = w_{2}*sin\beta_{2} = 193,04*sin21 = 69,18\frac{m}{s}$$

39. Obliczam siłę F działającą na łopatkę wirnika:


F = m * (w1+w2) = 637 * (142,4+193,04) = 213675 N

40. Obliczam siłę obwodową:


Fu = m * (w1u+w2u) = 637 * (124,5+180,2) = 194094 N

41. Obliczam siłę odśrodkową:


Fa = m * (w1a+w2a) + π * dp * l2 * (P1P2)=


=637 * (142,4+193,04) + π * 0, 943 * 0, 129 * (7,93*103−7,57*106) = 989956 N

42. Dobieram optymalną cięciwę profilu kierownicy:


$$\left( \frac{s}{l} \right)_{\text{opt}} = 0,8 \div 1,3$$

przyjmuję 0,8


$$s_{1} = l_{1}*\left( \frac{s}{l} \right)_{\text{opt}} = 0,125*0,8 = 0,1\ m$$

43. Obliczam smukłość łopatki kierowniczej:


$$\lambda_{1} = \frac{l_{1}}{s_{1}} = \frac{0,125}{0,1} = 1,250\ m$$

44. Obliczam współczynnik strat na kierownicy st1:


st1 = 1 − φ2 = 1 − 0, 962 = 0, 078

45. Obliczam współczynnik strat na wirniku st2:


st2 = 1 − ψ2 = 1 − 0, 892 = 0, 208

46. Dobieram optymalna cięciwę profilu łopatki wirnika:


$$\left( \frac{s}{l} \right)_{\text{opt}} = 0,3 \div 0,4$$

przyjmuję 0,4


$$s_{2} = l_{2}*\left( \frac{s}{l} \right)_{\text{opt}} = 0,129*0,4 = 0,0561\ m$$

47. Obliczam smukłość łopatek wirnika:


$$\lambda_{2} = \frac{l_{2}}{s_{2}} = \frac{0,129}{0,0516} = 2,5\ m$$

48. Obliczam podziałkę t1 dla kierownicy:

Dobieram podziałkę z Atlasu łopatki kierowniczej:


topt = 0, 7 ÷ 0, 85

Przyjmuję 0,7


t1 = s1 * topt = 0, 1 * 0, 7 = 0, 07 m

49. Obliczam liczbę łopatek kierowniczych:


$$z_{k} = \frac{\pi*d_{p}}{t_{1}} = \frac{\pi*0,943}{0,07} = 42$$

50. Obliczam podziałkę t1 dla wirnika:

Dobieram podziałkę z Atlasu łopatki wirnika:


topt = 0, 58 ÷ 0, 68

Przyjmuję 0,68


t2 = s2 * topt = 0, 0516 * 0, 68 = 0, 03588 m

51. Obliczam liczbę łopatek wirnika:


$$z_{w} = \frac{\pi*d_{p}}{t_{2}} = \frac{\pi*0,943}{0,03588} = 84$$

52. Obliczam wskaźnik minimalny Wmin1na zginanie profilu kierownicy:

Z Atlasu profili :


Wx1,   x1kr = 0, 5100 cm3


s = 51, 46 mm


$$W_{min1} = Wx_{1}x_{1}^{\text{kr}}*\left( \frac{s_{1}}{s^{'}} \right)^{3} = 0,5100*\left( \frac{100}{51,46} \right)^{3} = 3,74\ \text{cm}^{3}$$

53. Obliczam naprężenia zginające na łopatkę kierownicy σzg1:


$$\sigma_{zg1} = \frac{F_{u}*l_{1}}{W_{min1}*2*z_{k}} = \frac{194093*0,125}{3,74*2*42} = 77,2\ MPa$$

54. Obliczam wskaźnik minimalny Wminna zginanie profilu wirnika:

Z Atlasu profili :


Wx1,   x1kr = 0, 3000 cm3


s = 25, 63 mm


$$W_{min2} = Wx_{1}x_{1}^{\text{kr}}*\left( \frac{s_{2}}{s^{'}} \right)^{3} = 0,3000*\left( \frac{51,6}{25,63} \right)^{3} = 2,45\ \text{cm}^{3}$$

55. Obliczam naprężenia zginające na łopatkę wirnika σzg2:


$$\sigma_{zg2} = \frac{F_{u}*l_{2}}{W_{min2}*2*z_{w}} = \frac{194093*0,129}{2,45*2*84} = 60,8\ MPa$$


Wyszukiwarka