Nghiên cứu xác định các tham số nhiệt động lực học của động cơ tua bin phản lực trên cơ sở các kích thước hình học

i LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng Các số liệu, kết trình bày luận án trung thực chưa công bố công trình tác giả khác Tác giả luận án Bùi Văn Thƣởng ii LỜI CẢM ƠN Lời cảm ơn chân thành xin gửi tới PGS.TS Phạm Vũ Uy, PGS.TS Phạm Thế Phiệt tận tình hướng dẫn, giúp đỡ suốt trình thực hoàn thành luận án Tôi xin bày tỏ lòng biết ơn đến PGS.TS Vũ Quốc Trụ, PGS.TS Đặng Ngọc Thanh giúp đỡ, động viên nhiều trình nghiên cứu Xin chân thành cảm ơn nhà khoa học Quân đội cho ý kiến đóng góp quý báu Cảm ơn Đảng ủy, Ban giám đốc, Phòng Sau đại học Học viện Kỹ thuật Quân sự; Đảng ủy Thủ trưởng khoa Hàng không Vũ trụ, Bộ môn Động phản lực khoa Hàng không Vũ trụ; Viện Tên lửa viện KH&CNQS tạo điều kiện cho hoàn thành luận án Lời cuối cùng, xin cảm ơn bạn bè, đồng nghiệp gia đình quan tâm, chia xẻ, giúp đỡ suốt thời gian qua Hà Nội, Ngày 20 tháng năm 2016 Bùi Văn Thƣởng iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC CHỮ VIẾT TẮT vii DANH MỤC CÁC KÝ HIỆU viii DANH MỤC CÁC BẢNG x DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ xi MỞ ĐẦU Chƣơng TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU VỀ ĐỘNG CƠ TUA BIN PHẢN LỰC 1.1 Động hành trình TBPL tên lửa đối hải 1.2 Tình hình nghiên cứu ĐCTBPL giới nước liên quan đến nghiên cứu xác định tham số nhiệt động lực học 1.2.1 Tình hình nghiên cứu giới 1.2.2 Tình hình nghiên cứu nước 17 1.3 Luận giải việc đặt mục tiêu nội dung cần nghiên cứu đề tài 19 1.3.1 Phân tích khả thiết kế chế tạo ĐCTBPLN 19 1.3.2 Phân tích việc lựa chọn phương pháp xác định tham số nhiệt động lực học 23 1.3.3 Bố cục luận án 27 1.4 Kết luận chương 28 iv Chƣơng CÁC KÍCH THƢỚC HÌNH HỌC CƠ BẢN VÀ CƠ SỞ LÝ THUYẾT ĐỂ TÍNH TOÁN CÁC THAM SỐ NHIỆT ĐỘNG LỰC HỌC CỦA ĐỘNG CƠ 2.1 Xây dựng mô hình hình học động 29 2.1.1 Mô tả khảo sát nguyên lý hoạt động động 29 2.1.2 Xây dựng mô hình máy nén 31 2.1.3 Xây dựng mô hình buồng đốt 34 2.1.4 Xây dựng mô hình tua bin 35 2.1.5 Xây dựng mô hình thiết bị (miệng phun) 36 2.1.6 Các kích thước hình học động mẫu nghiên cứu 36 2.2 Các liệu xác định tham số nhiệt động lực học 39 2.2.1 Các tam giác tốc độ cấp máy nén tua bin chế độ tính toán 39 2.2.2 Các liệu để tính toán tham số nhiệt động lực học 41 2.3 Xây dựng mô hình toán cho sơ đồ động TBPLHD hai máy nén rotor có buồng trộn 47 2.3.1 Mô hình toán chế độ tính toán động 47 2.3.2 Kết tính toán sơ tham số nhiệt động lực học 52 2.4 Kết luận chương 55 Chƣơng MÔ PHỎNG XÁC ĐỊNH CÁC ĐẶC TÍNH HOẠT ĐỘNG CỦA CÁC BỘ PHẬN CHÍNH CỦA ĐỘNG CƠ 3.1 Các phần mềm ứng dụng tính toán mô 56 3.1.1 Phần mềm ANSYS 56 3.1.2 Phần mềm ASTRA 58 3.2 Ứng dụng tính toán mô ANSYS-CFX 59 3.2.1 Xây dựng mô hình mô 60 v 3.2.2 Chia lưới mô hình mô 63 3.2.3 Định nghĩa phần tử, mặt, xác định điều kiện ban đầu, điều kiện biên ―Setup‖ 65 3.2.4 Thực tính toán mô ―Solution‖ 69 3.2.5 Đưa kết tính toán mô ―Results‖ 71 3.3 Ứng dụng phần mềm ASTRA tính toán nhiệt động buồng đốt 80 3.3.1 Thiết lập tính toán nhiệt động buồng đốt 80 3.3.2 Phân tích xử lý kết 80 3.4 Xây dựng đặc tính phận động 83 3.4.1 Xây dựng đặc tính máy nén 83 3.4.2 Xây dựng đặc tính tua bin 87 3.4.3 Xây dựng đặc tính miệng phun 92 3.5 Kết luận chương 93 Chƣơng ĐỒNG BỘ CÁC ĐẶC TÍNH CỦA CÁC BỘ PHẬN ĐỘNG CƠ ĐỂ XÁC ĐỊNH CÁC THAM SỐ NHIỆT ĐỘNG LỰC HỌC 4.1 Các chế độ hoạt động cân chế độ tính toán 94 4.1.1 Các chế độ hoạt động cân đồ thị đặc tính 94 4.1.2 Các điều kiện chế độ hoạt động cân 96 4.2 Xác định đường hoạt động cân đặc tính máy nén động mẫu nghiên cứu 99 4.2.1 Các điều kiện đồng hoạt động phận động mẫu 100 4.2.2 Các bước xác định đường làm việc cân bằng cách đồng đặc tính theo vòng quay 100 4.3 Xác định tham số nhiệt động lực học 112 4.3.1 Xác định tham số nhiệt động lực học ĐTT (Hb =0, Mb =0) 112 4.3.2 Xác định tham số nhiệt động lực học theo Mb (Hb =0) 114 vi 4.4 Kết luận chương 117 KẾT LUẬN VÀ KIẾN NGHỊ 118 DANH MỤC CÔNG TRÌNH CỦA TÁC GIẢ 121 TÀI LIỆU THAM KHẢO 122 PHỤ LỤC 132 vii DANH MỤC CÁC CHỮ VIẾT TẮT Chữ viết tắt Ý nghĩa ĐCTBPL Động tua bin phản lực ĐCTBPLHD Động tua bin phản lực hai dòng ĐCTBPLN Động tua bin phản lực nhỏ ĐCDT Động dòng thẳng ĐTT Điểm tính toán KCB Khí cụ bay KHKT Khoa học kỹ thuật KH&CN Khoa học công nghệ MBKNL (UAV) Máy bay không người lái MNCA Máy nén cao áp MNTA Máy nén thấp áp TBB Thiết bị bay KCB Khí cụ bay TBĐL Thiết bị động lực TBPL Tua bin phản lực TBPLHD Tua bin phản lực hai dòng Chỉ số dƣới Ý nghĩa knd, kta Máy nén thấp áp kvd, kca Máy nén cao áp vent Quạt nén vca, rca Vào cao áp, cao áp tb Tua bin lp… Loa phụt… viii DANH MỤC CÁC KÝ HIỆU Ký hiệu Ý nghĩa p Áp suất T Nhiệt độ P Lực đẩy  Mật độ n Số vòng quay q Mật độ dòng G Lưu lượng g oxl Hệ số xả khí điều khiển làm mát ổ bi gt Hệ số cấp nhiên liệu m Mức phân dòng i Entanpi S Entropi  Tỷ số nén, tỷ số giãn nở  Hệ số bảo toàn áp suất toàn phần cp Nhiệt dung riêng đẳng áp N Công suất  Hiệu suất F Diện tích u Tốc độ vòng c, ca Tốc độ tốc độ dọc trục mv Hệ số lưu lượng chất khí  Hệ số tốc độ dòng q    ,     ,     Các hàm khí động Q Nhiệt lượng ix Hu Nhiệt trị thấp nhiên liệu Hb Độ cao bay Mb Số M bay (tốc độ bay) L Công X Tham số tua bin  Góc tốc độ vòng tốc độ tuyệt đối  Góc tốc độ vòng tốc độ tương đối  Khối lượng riêng  Trọng lượng riêng x DANH MỤC CÁC BẢNG Bảng 1.1 So sánh số lượng chi tiết rotor hai động 12 Bảng 1.2 Các tham số đặc tính ĐCTBPLN hãng 16 Bảng 1.3 Kết hoàn thiện đặc tính động 17 Bảng 2.1 Các đường kính máng chảy mặt cắt 37 Bảng 2.2 Các góc profil đường kính máy nén thấp áp 38 Bảng 2.3 Các góc β α vòng làm việc vòng chỉnh dòng 38 Bảng 2.4 Số làm việc tầng máy nén (rotor) 38 Bảng 2.5 Số chỉnh dòng dẫn dòng tầng máy nén (stator) 39 Bảng 2.6 Các góc profil hướng chiều làm việc tua bin 39 Bảng 2.7 Số lượng vòng làm việc hướng chiều tua bin 39 Bảng 2.8 Các tham số nhiệt động lực học (tính toán sơ bộ) 53 Bảng 2.9 Các tham số động theo số Mb (Hb = 0) 54 Bảng 3.1 Ký hiệu số tham số ―Expressions‖ 72 Bảng 3.2 Kết tính tham số máy nén ANSYS-CFX 75 Bảng 3.3 Kết tính tham số tua bin ANSYS-CFX 77 Bảng 3.4 Kết tính tham số miệng phun (loa phụt) 79 Bảng 3.5 Kết tính toán ASTRA 80 Bảng 3.6 Các tham số đơn vị đo ASTRA 82 Bảng 4.1 Kết tính đường hoạt động cân 110 Bảng 4.2 Các tham số NĐLH động mặt cắt 113 Bảng 4.3 Sự thay đổi tham số theo Mb (phụ lục 4) 114 Bảng 4.4 So sánh số tham số động (tài liệu tính toán) 116 Bảng 4.5 Các tham số riêng động (tính toán sơ thử nghiệm số) 116 154 Phụ lục 3.2 Một số kết tính ANSYS-CFX cho tua bin D:\Thuong\HOSO-NCS-THUONG\TUABIN100_1.8_files\dp0\CFX\CFX\Fluid Flow CFX_001.res Biểu diễn tốc độ dòng khí qua tua bin Biểu diễn thay đổi entanpi toàn phần tua bin 155 Biểu diễn thay đổi nhiệt độ toàn phần tua bin Biểu diễn thay đổi áp suất toàn phần tua bin Nhận xét: Trên hình biểu diễn tham số dòng khí cháy qua tua bin thấy rõ tham số áp suất, nhiệt độ giảm dần, tốc độ qua cấp tăng dần tua bin thành phần biến lượng nhiệt hàm sản phẩm cháy thành động sinh công quay rotor động 156 Phụ lục 3.3 Kết tính ANSYS-CFX cho loa Kết mô loa đưa hình ảnh số liệu tham số dòng sản phẩm khí thoát khỏi loa Biến đổi áp suất tổng thiết bị Biến đổi tốc độ dòng thiết bị Nhận xét: Tại loa dòng sản phẩm cháy tăng tốc độ thoát môi trường tạo lực đẩy phản lực áp suất, nhiệt độ mật độ dòng giảm xuống 157 Phụ lục 3.4 Kết tính toán ASTRA cho buồng đốt ┌───────────────────────────── Исходные данные ──────────────────────────────┐ I=0, p=0.7693,0.793113, Fuel=(C1H1.996[-1958]), Ox=(N21.9095 O6.699[261.5] Ar0.377 C1O20.0145), Alpha=5.989,5.311; └─────────────────────────────────────────── ─────────────────────────────────┘ Брутто-формула раб.тела: C 762689 H 1.52233 N 52.3479 O 16.0058 Теоретическое соотношение ок./гор = 13.2145 └─────────────────────────────────────────── ───────────────────────────── Характеристики равновесия - СИ P=0.7693 T=1028.67 V=0.384231 S=7.68192 I=233.805 U=-61.7808 M=34.5592 Cp=1.17654 k=1.32316 Cp"=1.17787 k"=1.32267 A=625.258 Mu=0.0000417 Lt=0.0702818 Lt"=0.070358 MM=28.9358 Cp.г=1.17654 k.г=1.32316 MM.г=28.9358 R.г=287.348 Z=0 Пл=0 Bm=0.101272 Содержание компонентов - мoль/кг O2 6.65376 OH 0.2042e-5 H2O 0.88805 N2 26.126 NO 0.0013796 NO2 0.0001634 CO2 0.889832 Характеристики равновесия - СИ P=0.793113 T=1028.67 V=0.372695 S=7.67316 I=233.805 U=-61.7808 M=34.5592 Cp=1.17654 k=1.32316 Cp"=1.17787 k"=1.32266 A=625.258 Mu=0.0000417 Lt=0.0702818 Lt"=0.070358 MM=28.9358 Cp.г=1.17654 k.г=1.32316 MM.г=28.9358 R.г=287.348 Z=0 Пл=0 Bm=0.101272 Содержание компонентов - мoль/кг O2 6.65376 OH 0.2026e-5 H2O 0.88805 N2 26.126 NO 0.0013796 NO2 0.0001659 CO2 0.889832 Характеристики равновесия - СИ P=0.7693 T=1084.18 V=0.404984 S=7.75116 I=230.319 U=-81.2318 M=34.5608 Cp=1.19011 k=1.31832 Cp"=1.19215 k"=1.31759 A=640.689 Mu=0.0000431 Lt=0.0734762 Lt"=0.0735973 158 MM=28.9345 Cp.г=1.19011 k.г=1.31832 R.г=287.361 Z=0 Пл=0 Bm=0.101785 Содержание компонентов - мoль/кг O2 6.47265 OH 0.5693e-5 H2O 0.999824 NO 0.00235 NO2 0.0001942 HNO2 0.1185e-5 Характеристики равновесия - СИ P=0.793113 T=1084.18 V=0.392825 S=7.7424 U=-81.2317 M=34.5608 Cp=1.19011 k=1.31832 k"=1.31759 A=640.689 Mu=0.0000431 Lt"=0.0735972 MM=28.9345 Cp.г=1.19011 k.г=1.31832 R.г=287.361 Z=0 Пл=0 Bm=0.101785 Содержание компонентов - мoль/кг O2 6.47265 OH 0.5650e-5 H2O 0.999824 NO 0.00235 NO2 0.0001971 HNO2 0.1212e-5 MM.г=28.9345 N2 26.084 CO2 1.00183 I=230.319 Cp"=1.19215 Lt=0.0734762 MM.г=28.9345 N2 26.084 CO2 1.00183 Qua tính toán ASTRA cho buồng đốt chế độ n max nhận kết tham số nhiệt động thành phần sản phẩm cháy dải áp suất 7,7- 7,93 at dải nhiệt độ 1028- 1084 K 159 Phụ lục ĐỒNG BỘ CÁC ĐẶC TÍNH Ở CHẾ ĐỘ VÒNG QUAY LỚN NHẤT (n = 100%) 1) Cho trước M.b = 0, H.b = 0: Th  288 K p h  101325Pa  k  1.4 kg  1.33 , R  287.3 J J Rg  288 kg K kg K m Va  340.352 Va  k R Th s kg h  1.225 m Vb  M b  Va  kc  0.97 J J Cpcp  1212 Cpv  1004 kg K kg K l0  14.9 II  0.9 zm  0.993 J Cpg  1160 kg K  g oxl  0.027 2) Thiết bị vào âm: v  0.98 Tzh  Th  1  ( k  1)  p zv  p h  v pzv  0.98 ph n k100  36000  Mb   k 3) Tính cho TBV: Tzv  Tzh Tzv 1 Th n k90  32400 n k80  28800 4) Chọn nhánh: nk=100% , điểm 2: Xác định giá trị đồ thị GvII kg kg m'  n k  2   GvI  4.1392 GvII  3.57207 GvI 60s s s  zvent  2.51667 Tzk  617.787K m'  0.863 Tzvent  396.225K p zvent  251443Pa p zk  793113Pa n k100 Tính giá trị:  zk  7.93818 Tzv 4 0.5 GvIqc  GvI  7.074 10 m K  s pzv Tzk  Tzv T kqc   1.1451 Tzv Tzv 0.5 Gvtaqc  GvI (  m')   0.0013179m K  s pzv Gnl Gt  t Q  Gt  Hu  g    Hu  42.910   J kg 160 T  Q GvI Cpcp  g t1  Cpcp  ( T ) Hu  g Tính hệ số khí dư cho buồng đốt:   1 1  g t1 l01 g t  0.01627 g t  l0 g t1  0.01264 2  g t  l0 l0  14.9 l01  13.2145 Lượng khí sau máy nén lấy lượng nhỏ điều khiển khoảng 0.025-0.03%, ta lấy giá trị 0.027 5) Xác định:  Gvbd  GvI  g oxl  Tzk 0.5 Gvbdqc  GvI  goxl   0.000126216m K  s pzk   Có thể tính theo công thức:   Tzv Tzk  Tzv 4 0.5 Gvbdqc1  GvI  goxl    1  1.2699453 10 m K  s pzv  zk Tzv   Chọn nhiệt độ sau buồng đốt: thay đổi cho đảm bảo cân lưu lượng: Tzg  Tzk  438.3415K p zg  p zk  kc  7.693 10 Pa Tính T.zg gần điểm cân nhất: Gvtbqc_max  0.00017013m K 0.5 s Chọn T.3 gần nhất: p zg   T3   Gvtbqc_max   1056.12885K Gvbd   DeltaTbd  1056.12885K  Tzk  438.342K 6) Chọn:   Tzg Tzv  3.66711 Có thể tính gần đúng: Tzg GvI     3.622    Gvbd    kc  zk  p zg Gvbd Tzv   GvI   p zv   7) Buồng đốt:  bd  Tzg Tzk  1.70954 Có thể tính Θ.bd theo công thức:  Tzg      1.71  Tzv  T kqc  8) Dòng sản phẩm cháy qua buồng đốt không trích khí làm mát mà toàn vào tua bin: Gvtb  Gvbd 161 Tzg 0.5 Gvtbqc  Gvtb  0.00017012997 m K  s pzg Có thể tính gần đúng: Tzv Tzg  1  Gvtbqc1  GvI  goxl      pzv  zk kc Tzv    0.5 Gvtbqc1  0.00017118m K  s n t  n k nt nk Tzv 1  116.004   116.004 0.5 0.5 Tzg Tzv Tzg K s K s Có thể tính: 9) Đọc đặc tính tua bin: Gvtbqcmax  0.00017013 m K 0.5 s  11 Gvtbqc  Gvtbqcmax  2.819 10 m K 0.5 s Điều kiện lưu lượng cân chọn tỷ số T.zg/T.zv: 6 0.5 G vtbqc  10 m K s Tzvcm  Tzvent p zvent  2.514 10 Pa p zvcm  p zvent  II Tzvent  396.225K Từ đồ thị đặc tính tua bin xác định:  zt  4.0 Cpcm  Rcm  T tbqc  0.22 m' Cpv  Cpg  m' m' R  Rg  m' Tzt  Tzg  T tbqc  Tzg Tzt  823.78K J J J Cpv  1004 Cpg  1160 Cpcm  1.088 10  kg K kg K kg K J Rcm  287.676 kg K Cpcm kcm  Cpcm  Rcm p zt  p zg  zt kcm  kcm  1  ckp      kcm 1 p'zvcm   ckp  1.869 p zvcm p zt p'zt  Pa Pa Xác định theo công thức thực nghiệm: p zcm  cm Pa  p'zt  p'zvcm  p zcm  202195.550841  Pa m'  1 m' Tzcm  m' Cpv  Tzvent  Cpg  Tzt Tzcm  640.972K (  m')  Cpcm mv  R      k  1 2 K s mv  0.0404 m k 1 Tzg  Tzt Tzg T tbqc  0.22 kg k k T tbqc  kg mg    Rg kg   K s mg  0.03964 m  1  kg kcm kcm mcm    Rcm kcm   K s mcm  0.03996 m  1  kcm  162 10) Xác định cho miệng phun 5 p zvcm  2.263 10 Pa Tzvent  396.225K p zcm  2.022 10 Pa p zt  1.923 10 Pa  zt  Tzvcm  396.225K Tzcm 0.5 Gvmpqc  Gvtb  GvI m'    0.00097595m K  s cm pzcm   Giá trị cần so sánh đồ thị đặc tính miệng phun  cm  p zt pzvcm  vcm  p zcm  cm  0.951 T cmqc  p zcm  vcm  1.119 Tzg  Tzcm T cmqc  0.393 Tzg T cmqcII  Tzvcm  Tzcm Tzvcm T cmqcII  0.618 Hoặc tính theo biểu thức từ tổng hai dòng: Tzg Tzvcm   1 G'mpqc   Gvtb   zt  cm   T cmqc  GvI  m'  vcm   T cmqcII   pzg cm pzvcm cm G'mpqc  0.00097595m K 0.5 s  c  1.996 11) Tính mức giãn loa phụt: Tính tỷ số giảm áp miệng phun theo hai dòng khác nhau:  c2   c2'  p zv ph p zv ph   zvent  II cm   zk  kc  zt  vcm  cm  cm 12) Trên đặc tính tua bin ta so sánh mức giãn nở đọc theo thông số lưu lượng dòng tìm π.c1 giá trị giãn nở tính toán theo mục 11) π.c2 Thể hình, theo trường hợp 10a), b), c) 13) Nếu tường hợp tua bin chưa bị nghẽn dòng, tức là: Gvtbqc  Gvtb Tzg pzg   Gvtbqc  max  Gvtb Tzg   pzg     max Sẽ xét đặc tính miệng phun:  ckp  1.869  c   ckp Trường hợp miệng phun chế độ tới hạn Bằng tính toán đảm bảo điều kiện cân lưu lượng vòng quay 14) Xét điều kiện cân công suất:  mp   c T ktaqc   Tzvent  Tzvcm Tzvent  Tzv  Tzv Tzg Tzg GvI Tzv Tzv  N qc  zm Cpg GvI  goxl     zk  kc  T tbqc   Cpv  T kqc  Cpv GvI m'  T k   p  Tzv pzv pzv zg   m N qc  0.407895 0.5 K s Tại điểm ΔW.qc ta thu điểm hoạt động cân nhánh 15) Dựng đặc tính máy nén đường HĐCB 163 16) Tính tham số động gian no hoan toan: kcm  1    kcm   kcm   ph   Cc1   c 2  Rcm Tzcm      kcm  p zcm     Cc1  472.9153m s 1 Nếu giãn nở không hoаn toаn:  c   ckp kcm  1    kcm    kcm  Cc2   c 2 R T  1     kcm  cm zcm    ckp   Cc2  451.7481m s 1 Phần phía trước loa làm việc chế độ giãn nở chưa hoàn toàn, sau giãn nở tiếp lõi trung tâm (vỏ phận sinh khí) nên tốc độ khí thoát tính theo tốc độ tương đương: Cc2 Tc2  Tzcm  2 kcm kcm  Tc2  547.164K Cc2  c  2  Rcm kcm kcm   Rcm Tc2 p c2  p zcm  ckp p c2  1.082 10 Pa  c  1.06 Mật độ dòng    kcm  q c     cs  ph p zcm 1 kcm 1    kcm   kcm    c    cs  c 2  kcm 1      q c  0.996 Cc2  cs  0.501 kcm Rcm Tc2  cs  0.977 kcm  k  1  2  cm c  c    kcm  1   kcm 1 p cs  p c2  c  c  0.491 p cs  2.203 10 Pa Van toc tuong duong Cc  Cc2  pc2  ph  Tc2  mcm qc  pcs m Cc  470.124 s Tc  539.38K Cc Tc  Tzcm  2 kcm kcm   Rcm p c  1.013 10 Pa Vì giãn nở loa giãn nở hoàn toàn Xác định tham số riêng động cơ, lưu lương không khí vа hiệu suất: kg Gv  GvI  GvII  7.71127 s Lực đẩy riêng lực đẩy: Pud   g t   g oxl   1    Cc  Vb  m'     P  Pud  Gv pc  pzcm c 164 Hiệu suất hiệu dụng động cơ: e  Pud  (  m')   2 g t Hu   g oxl e  0.3083 Lưu lượng nhiên liệu riêng   g t   g oxl Cud  Pud  (  m') kg Cud  0.06452 N hr kg Cud  0.645 10 N hr  kg   kG hr    Cong hieu dung cua chu trinh 1  goxl 2  1  g t   C  Vb   m'  c    Le  J Le  1.11447 10  kg Hiệu suất toàn phần hiệu suất bay: Vb 0  Pud   g t  g oxl        Hu  m'  Vb v  Pud  Le Xác định khối lượng riêng trọng lượng riêng xác định H=0 M=0 chế độ tính toán M dc  95kg  dc  M dc Gv  12.32 kg kg  dc  M dc P  0.026 kg N s  p zh     p zh  v  p   zventII   pzventII     p zk     p zg   p   zt   p zcm     pc   p  c   p x  Pa  Tzh     Tzh  T   zventII   TzventII     Tzk     Tzg   T   zt   Tzcm     Tc   T  c   Tx  K  c1a     c1a  c   aventII   cakvd     cak     C1a   C   2a   ctvdr     Cc   C  c   Cx  m s 165 288 196 396.23 2.51·10 158 396.23 2.51·10 169 617.79 7.93·10 134 1.06·10 7.69·10 209 823.78 1.92·10 153 640.97 2.02·10 153 539.38 1.01·10 470.12 539.38 1.01·10 470.12  7 8 196 9.73·10 3 augment xi1Txp x Cx   9.93·10 2 288 10 T x 10 px Cx 10 xi1 M  ( 0.2 0.4 0.6 0.8 0.85 1.0) V  M  Va V  ( 68.07 136.14 204.21 272.28 289.3 340.35)   g t  g oxl    Pud  1    Cc  Vb  m'       m s     ( k  1)  M2  Tzh  Th  1     Tzh  ( 288 290.304 297.216 308.736 324.864 329.616 345.6) K Tzv  Tzh   M  b_M  ( k  1) 2 21  M    b_M  ( 0.2182 0.4313 0.6348 0.8251 0.8704 ) k1  V  b  k 2 R  Tzh   b  ( 0.2182 0.4313 0.6348 0.8251 0.8704 ) k1 166   k  k 1  ( k  1)   2 b  b  1   k1     b   b  1  ( k  1)   k  1   b  ( 0.992 0.969 0.933 0.887 0.874 0.833)  b  ( 0.972 0.896 0.784 0.656 0.624 0.528)   1      ( k  1) k 1  ( k  1) 2 k 1 qb    b 1      k b     q b  ( 0.337 0.629 0.842 0.963 0.98 )   b  1  ( k  1)   b  k 1  b  ( 0.98 0.924 0.84 0.74 0.714 0.634) k1   k    k  1   Tzh   p zh  p h       Th    5 p zh  1.013 10 1.042 10 1.131 10 1.292 10 5 1.545 10 1.625 10 1.918 10 p zv  p zh  v  p zv  9.93  10 1.021 10 1.109 10 Dvvent  258mm d vvent  134mm Fvvent d rvent  178mm 2 D vvent  d vvent      Frvent FrventII 2 D vvent  Dgiua      1.267 10 FrventI 1.514 10 Frca  Pa Dgiua  224.8mm 2 D  d giua rvent      Fvvent  0.03818m Frvent  0.0274m FrventII  0.0126m d vca  176mm d rca  204mm Dca  220mm Fvca 1.88  10 2 D vvent  d rvent       D  d 2 vca   ca    1.593 10  Pa  D  d 2 rca   ca    FrventI  0.0148m Fvca  0.013685m 3 Frca  5.3281 10 m Nhiệt độ không khí cửa vào MN:   k 1 Tzh  Th     Mb     Áp suất không khí  k    k  1   Tzh   p zh  p h       Th   Tzh  ( 288 290.304 297.216 308.736 324.864 329.616 345.6) K 167 p zh  ( 101325 104190.5846 113134.6279 129240.4201 154453.7515 162506.936 191801.047)  Pa Tzv  ( 288 290.304 297.216 308.736 324.864 329.616 345.6) K p zv  ( 99298.5 102106.77291 110871.93538 126655.61173 151364.67644 159256.79728 187965.02607)  Pa Tốc độ tương đối  1a   c1a 2 k k1  R Tzv  1a  ( 0.632 0.63 0.623 0.611 0.595 0.591 0.577) Mật độ dòng tương đối  1    k 1  k  1 k 1  k 1 q1a     1a     1a    k1      q 1a  ( 0.84 0.837 0.831 0.821 0.807 0.803 0.789)   q1a    Gv1  mv 0.0389m  pzv   Tzv   kg Gv1  ( 7.7197 7.8862 8.3985 9.2957 10.6448 11.0633 12.5429) s   P  Pud  Gv1 PM0  3645N    3 3 3  P  3.645 10 3.225 10 2.983 10 2.883 10 2.903 10 2.925 10 3.031 10 N   Pud  (  m') e        e  ( 0.306 0.229 0.173 0.132 0.102 0.096 0.08 ) 2 g t Hu   g oxl   3600g t   g oxl s Cud  Cud  ( 0.065 0.075 0.086 0.099 0.112 0.116 0.127) Pud  (  m') m    g oxl    1  g t   Cc2  V2  m'    Le      V     g t  Cc V  V  0  Pud     b   g t  g oxl    2    Hu  g t   Cc  V   m'         0  ( 0.076 0.133 0.174 0.204 0.21 0.226)    b  ( 0.127 0.226 0.307 0.375 0.391 0.438) 168  stack e b 0  T  0.306   0.229  0.173   0.132   0.102  0.096  0.08  0   0.127 0.076 0.226 0.133 0.307 0.174  0.375 0.204 0.391 0.21   0.438 0.226