Nghiên cứu chế tạo và thử nghiệm mô hình sấy buồng theo nguyên lý bơm nhiệt

Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN THỊ ÚT HIỀN NGHIÊN CỨU CHẾ TẠO VÀ THỬ NGHIỆM MƠ HÌNH SẤY BUỒNG THEO NGUN LÍ BƠM NHIỆT Chun ngành : CƠNG NGHỆ NHIỆT LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 12 năm 2011 TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SAU ĐẠI HỌC Độc lập - Tự - Hạnh phúc -oOo - TP HCM, ngày tháng năm 2010 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: NGUYỄN THỊ ÚT HIỀN Phái: Nữ Ngày, tháng, năm sinh: 10/08/1983 Nơi sinh: Bình Chun ngành: CƠNG NGHỆ NHIỆT MSHV: Định 00609414 1- TÊN ĐỀ TÀI: NGHIÊN CỨU CHẾ TẠO VÀ THỬ NGHIỆM MƠ HÌNH SẤY BUỒNG THEO NGUN LÍ BƠM NHIỆT 2- NHIỆM VỤ LUẬN VĂN: - Nghiên cứu lý thuyết cơng nghệ sấy bơm nhiệt - Tính tốn thiết kế, chế tạo mơ hình sấy buồng theo nguyên lí bơm nhiệt - Thử nghiệm sấy hành mơ hình với nhiều chế độ - Phân tích, đánh giá kết khảo nghiệm - Lập trình tính tốn q trình sấy lí thuyết thực tế 3- NGÀY GIAO NHIỆM VỤ : 14-02-20110 4- NGÀY HOÀN THÀNH NHIỆM VỤ : 02-12-2011 5- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: GS TS LÊ CHÍ HIỆP Nội dung đề cương LV thạc sĩ Hội đồng Chuyên ngành thông qua CÁN BỘ HƯỚNG DẪN GS TS LÊ CHÍ HIỆP CHỦ NHIỆM BỘ MƠN KHOA QL CHUN QL CHUN NGÀNH NGÀNH GS TS LÊ CHÍ HIỆP CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học : GS TS LÊ CHÍ HIỆP ……………… (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét : (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét: (Ghi rõ họ, tên, học hàm, học vị chữ ký) Luận văn thạc sĩ bảo vệ Trường Đại Học Bách Khoa, ĐHQG Tp.HCM ngày…….tháng … năm …… Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: (Ghi rõ họ, tên, học hàm, học vị Hội đồng chấm bảo vệ luận văn thạc sĩ) …………………………………………………………………… …………………………………………………………………… …………………………………………………………………… …………………………………………………………………… …………………………………………………………………… Xác nhận Chủ tịch Hội đồng đánh giá LV Bộ môn quản lý chuyên ngành sau luận văn sửa chữa (nếu có) Chủ tịch Hội đồng đánh giá LV Bộ môn quản lý chuyên ngành LỜI CẢM ƠN Tác giả xin bày tỏ lòng biết ơn sâu sắc chân thành đến GS.TS Lê Chí Hiệp tận tình hướng dẫn đóng góp nhiều ý kiến quan trọng quý giá cho nội dung luận văn Tác giả xin gửi lời cảm ơn đến thầy cô anh em Bộ môn Công nghệ Nhiệt trường Đại học Công nghiệp Thực Phẩm Thành phố Hồ Chí Minh, thầy anh em Bộ môn Điện Lạnh trường Kinh tế- Kỹ thuật Phú Lâm Thành phố Hồ Chí Minh giúp đỡ tác giả góp phần vào thành cơng luận văn Thành cơng luận văn cịn thể biết ơn sâu sắc đến người thân, thương bạn bè đồng nghiệp hỗ trợ động viên tác giả suốt thời gian thực luận văn Tp Hồ Chí Minh, ngày tháng 12 năm 2011 Tác giả Nguyễn Thị Út Hiền TÓM TẮT LUẬN VĂN THẠC SĨ Trong nghiên cứu này, tiến hành tính tốn, thiết kế sơ mơ hình sấy buồng theo nguyên lí bơm nhiệt Mục đích trình tính tốn, thiết kế phục vụ cho cơng việc thử nghiệm, đánh giá kỹ thuật sấy bơm nhiệt hành Đồng thời xây dựng chương trình tính tốn q trình sấy lí thuyết thực tế cho mơ hình sấy buồng theo nguyên lí bơm nhiệt phần mềm visua basic Chúng tiến hành thử nghiệm sấy hành nhiều chế độ sấy khác thông qua việc điều chỉnh thông số nhiệt độ tác nhân sấy vận tốc sấy nhằm tìm chế độ hoạt động hợp lí mơ hình sấy buồng theo ngun lí bơm nhiệt chúng tơi so sánh kiểm tra tính khả thi mơ hình, hiệu lượng chất lượng sản phẩm sấy Thơng qua việc xử lí phân tích kết thí nghiệm để đánh giá yếu tố nhiệt độ tác nhân sấy, tốc độ sấy ảnh hưởng đến chất lượng sản phẩm sấy (màu sắc, mùi vị, thành phần vitamin C) hiệu hệ thống sấy (thời gian sấy, điện tiêu thụ) GVHD: GS.TS LÊ CHÍ HIỆP - i - MỤC LỤC Chương 1: Mở đầu 1.1 Lý chọn đề tài 1.1.1 Sơ lược tình hình sấy 1.1.2 Tầm quan trọng việc sấy nông sản thực phẩm 1.2 Mục đích, đối tượng phạm vi nghiên cứu 1.2.1 Mục đích nghiên cứu 1.2.2 Đối tượng nghiên cứu 1.2.3 Phạm vi nghiên cứu 1.3 Nội dung phương pháp nghiên cứu 1.4 Ý nghĩa khoa học thực tiễn đề tài nghiên cứu Chương 2: Tổng quan 2.1 Tình hình sản xuất bn bán sử dụng hành 2.2 Giới thiệu hành 2.3 Thực trạng kỹ thuật công nghệ sấy hành 10 2.4 Tìm hiểu chung trình sấy 12 2.4.1 Khái niệm trình sấy 12 2.4.2 Ẩm vật liệu sấy 12 2.4.3 Phân loại vật liệu ẩm đặc tính xốp 13 2.4.3.1 Phân loại vật liệu ẩm 13 2.4.3.2 Đặc tính xốp vật liệu ẩm 13 2.4.4 Các dạng liên kết lượng liên kết ẩm 13 2.4.5 Truyền nhiệt truyền chất động học trình sấy 14 2.4.5.1 Truyền nhiệt truyền chất trình sấy 14 2.4.5.2 Động học trình sấy 15 2.4.6 Các yếu tố ảnh hưởng đến trình sấy 17 2.4.6.1 Ảnh hưởng nhiệt độ trình sấy 17 2.4.6.2 Ảnh hưởng độ ẩm trình sấy 17 2.4.6.3 Ảnh hưởng lượng tác nhân sấy 17 2.4.6.4 Ảnh hưởng độ dày lớp vật liệu sấy 18 HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - ii - 2.5 Các phương pháp sấy 18 2.5.1 Sấy tự nhiên 18 2.5.2 Sấy thiết bị 18 2.5.2.1 Phương pháp sấy nóng 18 2.5.2.2 Phương pháp sấy lạnh 20 2.6 Lý thuyết sấy nhiệt độ thấp 22 Chương 3: Tính tốn thiết kế mơ hình thực nghiệm 26 3.1 Lựa chọn mơ hình thiết kế 26 3.2 Các thơng số tính tốn 27 3.3 Tính tốn kích thước buồng sấy 28 3.4 Xây dựng trình sấy lý thuyết đồ thị t-d 28 3.4.1 Đồ thị t-d 28 3.4.2 Tính tốn q trình sấy 29 3.5 Xây dựng trình sấy thực tế đồ thị t-d 32 3.5.1 Cân nhiệt cho trình sấy thực tế 32 3.5.2 Đồ thị t-d 39 3.5.3 Tính tốn q trình sấy thực tế 39 3.6 Tính tốn thiết kế máy sấy bơm nhiệt 41 3.6.1 Chọn môi chất nạp thông số môi chất 41 3.6.2 Chọn tính tốn chu trình bơm nhiệt 42 3.6.3 Xây dựng đồ thị lập bảng xác định giá trị điểm nút 43 3.6.4 Tính tốn chu trình 44 3.6.5 Tính chọn máy nén 45 3.6.6 Tính tốn dàn ngưng 47 3.6.7 Tính tốn dàn lạnh 54 3.6.8 Chọn đường ống dẫn môi chất 60 3.6.9 Tính tốn dàn ngưng phụ 61 3.7 Chương trình tính tốn 61 3.7.1 Tính thể tích buồng sấy 62 3.7.2 Tính tốn q trình sấy lí thuyết 63 HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - iii - 3.7.3 Tính tốn q trình sấy thực 64 Chương : Thí nghiệm bàn luận 69 4.1 Giới thiệu mơ hình thí nghiệm 69 4.1.1 Mơ hình thực nghiệm 69 4.1.2 Thơng số mơ hình 69 4.1.3 Vận hành điều chỉnh mơ hình 70 4.2 Phương pháp tiến hành thí nghiệm 73 4.2.1 Phương pháp đo đạc 73 4.2.2 Phương pháp phân tích cảm quan 76 4.3 Kết thơng số máy chạy không tải 77 4.3.1 Các thông số thực nghiệm điều chỉnh tỷ số nén bơm nhiệt 77 4.3.2 Các thông số thực nghiệm điều chỉnh tốc độ gió sau buồng sấy 78 4.3.3 Các thông số thực nghiệm để điều chỉnh nhiệt độ TNS trước buồng 79 4.4 Lựa chọn thông số khảo sát trình sấy 80 4.4.1 Nhiệt độ TNS (TNS ) 80 4.4.2 Vận tốc TNS ( lưu lượng TNS) 80 4.4.3 Khối lượng vật liệu sấy mẻ 81 4.5 Phương pháp xử lý số liệu 82 4.5.1 Quy hoạch thực nghiệm 82 4.5.2 Phương pháp xử lý kết thực nghiệm 83 4.6 Thực nghiệm sấy hành máy sấy bơm nhiệt 83 4.6.1 Mục đích yêu cầu 83 4.6.2 Vật liệu khảo nghiệm 83 4.6.3 Xác định tốc độ giảm ẩm sấy hành 85 4.6.3.1 Quy trình sấy 85 4.6.3.2 Số lượng thí nghiệm 85 4.6.3.3 Kết khảo nghiệm sấy hành bàn luận 86 4.7 Nhận xét đánh giá 97 4.7.1 Tốc độ giảm ẩm mẫu thay đổi tốc độ tác nhân sấy nhiệt độ 350C HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ 97 LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - iv - 4.7.2 Tốc độ giảm ẩm mẫu thay đổi tốc độ tác nhân sấy nhiệt độ 400C 98 4.7.3 Tốc độ giảm ẩm mẫu sấy tốc độ tác nhân sau buồng sấy 99 4.8 Đánh giá cảm quan 100 4.9 Lượng ẩm tách đơn vị lượng 101 4.10 So sánh chất lượng sản phẩm sấy phương pháp sấy khác 103 Chương 5: Kết luận kiến nghị 105 5.1 Kết luận 105 5.2 Kiến nghị 106 TÀI LIỆU THAM KHẢO 107 Phụ lục1: Bảng số liệu thí nghiệm sấy hành nhiệt độ 400 C,tốc độ TNS sau buồng sấy  = 0,5 m/s, khối lượng VLS G=1220g 108 Phụ lục 2: Bảng số liệu thí nghiệm sấy hành nhiệt độ 400 C, tốc độ 109 TNS sau buồng sấy  = 1,0 m/s, khối lượng VLS G=1250g Phụ lục 3: Bảng số liệu thí nghiệm sấy hành nhiệt độ 400 C,tốc độ 110 TNS sau buồng sấy  = 1,5 m/s, khối lượng VLS G=1250g Phụ lục 4: Bảng số liệu thí nghiệm sấy hành nhiệt độ 350 C, tốc độ TNS sau buồng sấy  = 0,5 m/s, khối lượng vật liệu sấy G=1250 g 111 Phụ lục 5: Bảng số liệu thí nghiệm sấy hành nhiệt độ 350 C, tốc độ TNS sau buồng sấy  = 1,0 m/s, khối lượng VLS G=1150g 112 Phụ lục 6: Bảng số liệu thí nghiệm sấy hành nhiệt độ 350 C, tốc độ TNS sau buồng sấy  = 1,5 m/s, khối lượng VLS G= 1100g 113 Phụ lục 7: Một số hình ảnh mẫu sấy bơm nhiệt 114 Phụ lục 8: Một số hình ảnh mẫu sấy phương pháp khác kết phân tích thành phần vitaminC 115 Phụ lục 9: Một số hình ảnh mơ hình sấy buồng theo ngun lí bơm nhiệt 116 Phụ lục 10: Một số hình ảnh dụng cụ đo HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ 118 LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - v - Phụ lục 11: Thông số kỹ thuật máy nén 119 Phụ lục 12: Chương trình lập trình tính tốn nhiệt mơ hình sấy buồng theo ngun lí bơm nhiệt 120 DANH MỤC CÁC BẢNG BIỂU Tr ang HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 119 - Phụ lục 10: Một số hình ảnh dụng cụ đo Hình dụng cụ đo độ ẩm Hình dụng cụ đo nhiệt độ Hình dụng cụ đo tốc độ gió HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 120 - Phụ lục 12: Chương trình lập trình tính tốn nhiệt mơ hình sấy buồng theo ngun lí bơm nhiệt GIAO DIEN CHINH Private Sub MENUSAYLYTHUYET_Click() QUATRINHSAYLYTHUYET.Show End Sub Private Sub MENUSAYTHUC_Click() QUATRINHSAYTHUC.Show End Sub Private Sub menutinh_Click() B1.Show End Sub Private Sub MENUTONTHATNHIET_Click() CHUONG.Show End Sub TÍNH TOÁN THIẾT KẾ BUỒNG SẤY 'Nut Tinh' Private Sub Command1_Click() Dim a, b, c, d, TTHD, TTTB If ((IsNumeric(Text1.Text)) And (IsNumeric(Text2.Text)) And (IsNumeric(Text3.Text)) And (IsNumeric(Text4.Text))) Then a = CSng(Text1.Text) b = CSng(Text2.Text) c = CSng(Text3.Text) d = CSng(Text4.Text) TTHD = Round(a / (b * c), 2) Text5.Text = CStr(TTHD) TTTB = TTHD + d Text6.Text = CStr(TTTB) End If End Sub 'Chuyen tu nut Tab qua nut Enter' Private Sub Form_KeyPress(KeyAscii As Integer) If KeyAscii = 13 Then SendKeys "{TAB}" KeyAscii = End If End Sub TÍNH TỐN Q TRÌNH SẤY LÝ THUYẾT 'Nut Tinh' Private Sub Command1_Click() trangthai0_Click trangthai1_Click trangthai2_Click trangthai3_Click tinhG2_Click End Sub 'Chuyen tu nut Tad qua nut Enter' Private Sub Form_KeyPress(KeyAscii As Integer) If KeyAscii = 13 Then HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 121 - SendKeys "{TAB}" KeyAscii = End If End Sub 'Thong so trang thai 0' Private Sub trangthai0_Click() Dim p, t0, phi0, d0, i0, pb0 p = CSng(Txtp.Text) t0 = CSng(Txtt0.Text) phi0 = CSng(Txtphi0.Text) pb0 = Exp((12 - (4026.42 / (235.5 + t0)))) d0 = Round(0.621 * (phi0 / 100) * (pb0) / (p - (phi0 / 100) * pb0), 4) i0 = Round(1.004 * t0 + d0 * (2500 + 1.842 * t0), 4) Txtd0.Text = CStr(d0) Txti0.Text = CStr(i0) End Sub 'Thong so trang thai 1' Private Sub trangthai1_Click() Dim p, t1, phi1, d1, i1, pb1 p = CSng(Txtp.Text) t1 = CSng(Txtt1.Text) phi1 = CSng(Txtphi1.Text) pb1 = Exp((12 - (4026.42 / (235.5 + t1)))) d1 = Round(0.621 * (phi1 / 100) * (pb1) / (p - (phi1 / 100) * pb1), 4) i1 = Round(1.004 * t1 + d1 * (2500 + 1.842 * t1), 4) Txtd1.Text = CStr(d1) Txti1.Text = CStr(i1) End Sub 'Thong so trang thai 2' Private Sub trangthai2_Click() Dim p, t2, phi2, d2, i2, pb2 p = CSng(Txtp.Text) t2 = CSng(Txtt2.Text) d2 = CSng(Txtd1.Text) pb2 = Exp((12 - (4026.42 / (235.5 + t2)))) phi2 = Round(d2 * p / (pb2 * (0.621 + d2)), 4) * 100 i2 = Round(1.004 * t2 + d2 * (2500 + 1.842 * t2), 4) Txtphi2.Text = CStr(phi2) Txti2.Text = CStr(i2) Txtd2.Text = CStr(d2) End Sub 'Thong so trang thai 3' Private Sub trangthai3_Click() Dim p, t3, phi3, d3, i3, pb3 p = CSng(Txtp.Text) t3 = CSng(Txtt3.Text) i3 = CSng(Txti2.Text) d3 = Round((i3 - 1.004 * t3) / (2500 + 1.842 * t3), 4) pb3 = Exp((12 - (4026.42 / (235.5 + t3)))) phi3 = Round(d3 * p / (pb3 * (0.621 + d3)), 4) * 100 Txtphi3.Text = CStr(phi3) Txtd3.Text = CStr(d3) Txti3.Text = CStr(i3) HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 122 - End Sub 'Tinh G2, w' Private Sub tinhG2_Click() Dim W, G1, w1, w2, G2, t, d2, d3, Ltt, G, i1, i2, qtt, Q, Qk, i3, qttl, Qll, Q0 w1 = CSng(Txtdoam1.Text) w2 = CSng(Txtdoam2.Text) G1 = CSng(Txtg1.Text) W = Round(G1 * (w1 - w2) / (100 - w2), 3) G2 = Round(G1 - W, 4) Txtg2 = CStr(G2) Txtla = CStr(W) 'tinh Luong KK tuan hoan Qt say' d3 = CSng(Txtd3.Text) d2 = CSng(Txtd2.Text) t = CSng(Txth.Text) Llt = / (d3 - d2) L = Round(W * Llt, 2) TxtL.Text = CStr(L) G = Round(L / (t * 3600), 4) TxtLs.Text = CStr(G) 'Tinh nang suat nhiet dan nong' i1 = CSng(Txti1.Text) i2 = CSng(Txti2.Text) qtt = (i2 - i1) / (d3 - d2) Q = W * qtt Qk = Round(Q / (t * 3600), 4) Txtdn.Text = CStr(Qk) 'Tinh nag suat nhiet dan lanh' i3 = CSng(Txti3.Text) qttl = Llt * (i3 - i1) Qll = W * qttl Q0 = Round(Qll / (t * 3600), 4) Txtdl.Text = CStr(Q0) End Sub TÍNH TỔN THẤT NHIỆT 'Chuyen nut tad qua nut Enter' Private Sub Form_KeyPress(KeyAscii As Integer) If KeyAscii = 13 Then SendKeys "{TAB}" KeyAscii = End If End Sub 'Tinh nhiet trung binh buong say' Private Sub Txttss_Change() Dim tf1, t2, t3 If IsNumeric(Txttts.Text) And IsNumeric(Txttss.Text) Then t2 = CSng(Txttts.Text) t3 = CSng(Txttss.Text) tf1 = (t2 + t3) / Txttbtbs.Text = CStr(tf1) End If End Sub HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 123 - 'Tinh nhiet trung binh buong say' Private Sub Txttts_Change() Dim tf1, t2, t3 If IsNumeric(Txttts.Text) And IsNumeric(Txttss.Text) Then t2 = CSng(Txttts.Text) t3 = CSng(Txttss.Text) tf1 = (t2 + t3) / Txttbtbs.Text = CStr(tf1) End If End Sub 'Tinh anpha 1' Private Sub Txtvt_Change() Dim v, d, nuy, Re, Nu, lamda, anpha1 If IsNumeric(Txtvt.Text) And IsNumeric(Txtd1.Text) And IsNumeric(Txtnuy1.Text) And IsNumeric(Txtlamda1.Text) Then v = CSng(Txtvt.Text) d = CSng(Txtd1.Text) nuy = CSng(Txtnuy1.Text) lamda = CSng(Txtlamda1.Text) Re = v * (d / 1000) * ((nuy) ^ (-1)) If (Re > 100000) Then Nu = 0.032 * ((Re) ^ 0.8) End If If (Re 100000) Then Nu = 0.032 * ((Re) ^ 0.8) End If If (Re 100000) Then Nu = 0.032 * ((Re) ^ 0.8) End If If (Re 100000) Then Nu = 0.032 * ((Re) ^ 0.8) End If If (Re q1) Then dentaq = ((Q - q1) / Q) * 100 End If If (q1 > Q) Then dentaq = ((q1 - Q) / q1) * 100 End If Txtdentaq = CStr(dentaq) HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 125 - End Sub 'Tinh mat dong nhiet q' Private Sub Txttbtbs_Change() Dim anpha1, tf1, tw2, Q If IsNumeric(Txtanpha1.Text) And IsNumeric(Txttbtbs.Text) And IsNumeric(Txttw2.Text) Then anpha1 = CSng(Txtanpha1.Text) tf1 = CSng(Txttbtbs.Text) tw2 = CSng(Txttw2.Text) Q = Round(anpha1 * (tf1 - tw2), 4) Txtq = CStr(Q) End If End Sub 'Tinh mat dong nhiet q' Private Sub Txtanpha1_Change() Dim anpha1, tf1, tw2, Q If IsNumeric(Txtanpha1.Text) And IsNumeric(Txttbtbs.Text) And IsNumeric(Txttw2.Text) Then anpha1 = CSng(Txtanpha1.Text) tf1 = CSng(Txttbtbs.Text) tw2 = CSng(Txttw2.Text) Q = Round(anpha1 * (tf1 - tw2), 4) Txtq = CStr(Q) End If End Sub 'Tinh nhiet vach ngoai va nhiet trung binh vach ngoai' Private Sub Txthsdn3_Change() Dim cd1, hsdn1, cd2, hsdn2, cd3, hsdn3, tw1, Q, tf2, ttb2, tw2, t0 If IsNumeric(Txtchieuday1.Text) And IsNumeric(Txthsdn1.Text) And IsNumeric(Txtchieuday2.Text) And IsNumeric(Txthsdn2.Text) And IsNumeric(Txtchieuday3.Text) And IsNumeric(Txthsdn3.Text) Then cd1 = CSng(Txtchieuday1.Text) hsdn1 = CSng(Txthsdn1.Text) cd2 = CSng(Txtchieuday2.Text) hsdn2 = CSng(Txthsdn2.Text) cd3 = CSng(Txtchieuday3.Text) hsdn3 = CSng(Txthsdn3.Text) tw2 = CSng(Txttw2.Text) Q = CSng(Txtq.Text) t0 = CSng(Txttn.Text) tw1 = Round(tw2 - Q * ((cd1 / (1000 * hsdn1)) + (cd2 / (1000 * hsdn2)) + (cd3 / (1000 * hsdn3))), 1) ttb2 = (tw1 + t0) / Txttvn = CStr(tw1) Txttbbn = CStr(ttb2) End If End Sub 'Tinh nhiet vach ngoai va nhiet trung binh vach ngoai' Private Sub Txtchieuday3_Change() Dim cd1, hsdn1, cd2, hsdn2, cd3, hsdn3, tw1, Q, tf2, ttb2, tw2, t0 If IsNumeric(Txtchieuday1.Text) And IsNumeric(Txthsdn1.Text) And IsNumeric(Txtchieuday2.Text) And IsNumeric(Txthsdn2.Text) And IsNumeric(Txtchieuday3.Text) And IsNumeric(Txthsdn3.Text) Then cd1 = CSng(Txtchieuday1.Text) hsdn1 = CSng(Txthsdn1.Text) cd2 = CSng(Txtchieuday2.Text) hsdn2 = CSng(Txthsdn2.Text) HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 126 - cd3 = CSng(Txtchieuday3.Text) hsdn3 = CSng(Txthsdn3.Text) tw2 = CSng(Txttw2.Text) Q = CSng(Txtq.Text) t0 = CSng(Txttn.Text) tw1 = Round(tw2 - Q * ((cd1 / (1000 * hsdn1)) + (cd2 / (1000 * hsdn2)) + (cd3 / (1000 * hsdn3))), 1) ttb2 = (tw1 + t0) / Txttvn = CStr(tw1) Txttbbn = CStr(ttb2) End If End Sub 'Tinh nhiet vach ngoai va nhiet trung binh vach ngoai' Private Sub Txthsdn2_Change() Dim cd1, hsdn1, cd2, hsdn2, cd3, hsdn3, tw1, Q, tf2, ttb2, tw2, t0 If IsNumeric(Txtchieuday1.Text) And IsNumeric(Txthsdn1.Text) And IsNumeric(Txtchieuday2.Text) And IsNumeric(Txthsdn2.Text) And IsNumeric(Txtchieuday3.Text) And IsNumeric(Txthsdn3.Text) Then cd1 = CSng(Txtchieuday1.Text) hsdn1 = CSng(Txthsdn1.Text) cd2 = CSng(Txtchieuday2.Text) hsdn2 = CSng(Txthsdn2.Text) cd3 = CSng(Txtchieuday3.Text) hsdn3 = CSng(Txthsdn3.Text) tw2 = CSng(Txttw2.Text) Q = CSng(Txtq.Text) t0 = CSng(Txttn.Text) tw1 = Round(tw2 - Q * ((cd1 / (1000 * hsdn1)) + (cd2 / (1000 * hsdn2)) + (cd3 / (1000 * hsdn3))), 1) ttb2 = (tw1 + t0) / Txttvn = CStr(tw1) Txttbbn = CStr(ttb2) End If End Sub 'Tinh nhiet vach ngoai va nhiet trung binh vach ngoai' Private Sub Txtchieuday2_Change() Dim cd1, hsdn1, cd2, hsdn2, cd3, hsdn3, tw1, Q, tf2, ttb2, tw2, t0 If IsNumeric(Txtchieuday1.Text) And IsNumeric(Txthsdn1.Text) And IsNumeric(Txtchieuday2.Text) And IsNumeric(Txthsdn2.Text) And IsNumeric(Txtchieuday3.Text) And IsNumeric(Txthsdn3.Text) Then cd1 = CSng(Txtchieuday1.Text) hsdn1 = CSng(Txthsdn1.Text) cd2 = CSng(Txtchieuday2.Text) hsdn2 = CSng(Txthsdn2.Text) cd3 = CSng(Txtchieuday3.Text) hsdn3 = CSng(Txthsdn3.Text) tw2 = CSng(Txttw2.Text) Q = CSng(Txtq.Text) t0 = CSng(Txttn.Text) tw1 = Round(tw2 - Q * ((cd1 / (1000 * hsdn1)) + (cd2 / (1000 * hsdn2)) + (cd3 / (1000 * hsdn3))), 1) ttb2 = (tw1 + t0) / Txttvn = CStr(tw1) Txttbbn = CStr(ttb2) End If End Sub 'Tinh nhiet vach ngoai va nhiet trung binh vach ngoai' HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 127 - Private Sub Txthsdn1_Change() Dim cd1, hsdn1, cd2, hsdn2, cd3, hsdn3, tw1, Q, tf2, ttb2, tw2, t0 If IsNumeric(Txtchieuday1.Text) And IsNumeric(Txthsdn1.Text) And IsNumeric(Txtchieuday2.Text) And IsNumeric(Txthsdn2.Text) And IsNumeric(Txtchieuday3.Text) And IsNumeric(Txthsdn3.Text) Then cd1 = CSng(Txtchieuday1.Text) hsdn1 = CSng(Txthsdn1.Text) cd2 = CSng(Txtchieuday2.Text) hsdn2 = CSng(Txthsdn2.Text) cd3 = CSng(Txtchieuday3.Text) hsdn3 = CSng(Txthsdn3.Text) tw2 = CSng(Txttw2.Text) Q = CSng(Txtq.Text) t0 = CSng(Txttn.Text) tw1 = Round(tw2 - Q * ((cd1 / (1000 * hsdn1)) + (cd2 / (1000 * hsdn2)) + (cd3 / (1000 * hsdn3))), 1) ttb2 = (tw1 + t0) / Txttvn = CStr(tw1) Txttbbn = CStr(ttb2) End If End Sub 'Tinh nhiet vach ngoai va nhiet trung binh vach ngoai' Private Sub Txtchieuday1_Change() Dim cd1, hsdn1, cd2, hsdn2, cd3, hsdn3, tw1, Q, tf2, ttb2, tw2, t0 If IsNumeric(Txtchieuday1.Text) And IsNumeric(Txthsdn1.Text) And IsNumeric(Txtchieuday2.Text) And IsNumeric(Txthsdn2.Text) And IsNumeric(Txtchieuday3.Text) And IsNumeric(Txthsdn3.Text) Then cd1 = CSng(Txtchieuday1.Text) hsdn1 = CSng(Txthsdn1.Text) cd2 = CSng(Txtchieuday2.Text) hsdn2 = CSng(Txthsdn2.Text) cd3 = CSng(Txtchieuday3.Text) hsdn3 = CSng(Txthsdn3.Text) tw2 = CSng(Txttw2.Text) Q = CSng(Txtq.Text) t0 = CSng(Txttn.Text) tw1 = Round(tw2 - Q * ((cd1 / (1000 * hsdn1)) + (cd2 / (1000 * hsdn2)) + (cd3 / (1000 * hsdn3))), 1) ttb2 = (tw1 + t0) / Txttvn = CStr(tw1) Txttbbn = CStr(ttb2) End If End Sub 'Tinh anpha2' Private Sub Txtlamda2_Change() Dim g, beta, dentat, d, nuy2, pr, nu2, c, n, a, t0, lamda2, anpha2 If IsNumeric(Txtnuy2.Text) And IsNumeric(Txtlamda2.Text) And IsNumeric(Txtpr2.Text) Then ttb2 = CSng(Txttbbn.Text) t0 = CSng(Txttn.Text) nuy2 = CSng(Txtnuy2.Text) lamda2 = CSng(Txtlamda2.Text) pr = CSng(Txtpr2.Text) d = CSng(Txtd1.Text) g = 9.81 beta = / (ttb2 + 273) dentat = (ttb2 - t0) a = g * beta * dentat * ((d / 1000) ^ 3) * pr / (nuy2 ^ 2) HVTH: NGUYỄN THỊ ÚT HIỀN THẠC SĨ LUẬN VĂN GVHD: GS.TS LÊ CHÍ HIỆP - 128 - If (0.0001