ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN VĂN THƠNG THÍ NGHIỆM HẦM GIĨ ẢO CHO CƠNG TRÌNH NHÀ BẰNG MƠ PHỎNG CFD Chun ngành: Kỹ thuật xây dựng cơng trình dân dụng công nghiệp Mã số ngành: 60 58 02 08 LUẬN VĂN THẠC SĨ Tp Hồ Chí Minh, 02 - 2018 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: Cán hướng dẫn 1: TS Nguyễn Hữu Thành Cán hướng dẫn 2: PGS.TS Nguyễn Quốc Ý Cán chấm nhận xét 1: TS Lê Văn Phước Nhân Cán chấm nhận xét 2: TS Lê Trung Kiên Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG Tp HCM, ngày 07 tháng 02 năm 2018 Thành phần Hội đồng đánh giá Luận văn thạc sĩ gồm: Chủ tịch: PGS TS Nguyễn Văn Hiệp Thư kí: PGS TS Nguyễn Minh Long Phản biện (thành viên): TS Lê Văn Phước Nhân Phản biện (thành viên): TS Lê Trung Kiên Thành viên: TS Hồ Hữu Chỉnh CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc Lập - Tự Do - Hạnh Phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: NGUYỄN VĂN THÔNG MSHV: 1570658 Ngày, tháng, năm sinh: 20/08/1992 Nơi sinh: Bình Định Chuyên ngành: KTXD cơng trình dân dụng cơng nghiệp MN: 60 58 02 08 I TÊN ĐỀ TÀI: Thí nghiệm hầm gió ảo cho cơng trình nhà mơ CFD II NHIỆM VỤ VÀ NỘI DUNG Tìm hiểu sở lý thuyết CFD, phần mềm Ansys Fluent chất thơng số phân tích gió Sử dụng phần mềm Ansys Fluent mơ phân tích nghiên cứu công bố Sử dụng phương pháp CFD để xác định tải trọng gió tĩnh tác dụng lên cơng trình III NGÀY GIAO NHIỆM VỤ : 10/07/2017 IV NGÀY HOÀN THÀNH NHIỆM VỤ : 07/02/2018 V HỌ VÀ TÊN CBHD: TS Nguyễn Hữu Thành, PGS.TS Nguyễn Quốc Ý Tp HCM, ngày… tháng… Năm … CHỦ TỊCH HỘI ĐỒNG NGÀNH CÁN BỘ HƯỚNG DẪN TS Nguyễn Hữu Thành CÁN BỘ HƯỚNG DẪN TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG PGS TS Nguyễn Quốc Ý LỜI CẢM ƠN Luận văn thạc sĩ Xây dựng cơng trình dân dụng công nghiệp nằm hệ thống luận cuối khóa nhằm trang bị cho Học viên cao học khả tự nghiên cứu, biết cách giải vấn đề cụ thể đặt thực tế xây dựng… Đó trách nhiệm niềm tự hào học viên cao học Để hoàn thành luận văn này, cố gắng nỗ lực thân, học viên nhận giúp đỡ từ tập thể cá nhân Học viên xin ghi nhận tỏ lòng biết ơn đến tập thể cá nhân dành cho học viên giúp đỡ quý báu Đầu tiên Học viên xin bày tỏ lòng biết ơn sâu sắc đến thầy TS Nguyễn Hữu Thành thầy PGS.TS Nguyễn Quốc Ý Thầy đưa gợi ý để hình thành nên ý tưởng đề tài Thầy góp ý nhiều cách nhận định đắn vấn đề nghiên cứu, cách tiếp cận nghiên cứu hiệu giúp đỡ lúc khó khăn trình thực đề tài Học viên xin chân thành cảm ơn quý Thầy Cô Khoa Kỹ Thuật Xây dựng, trường Đại học Bách Khoa Tp.HCM truyền dạy kiến thức quý giá, kiến thức thiếu đường nghiên cứu khoa học nghiệp sau Luận văn thạc sĩ hoàn thành thời gian quy định với nỗ lực thân, nhiên khơng có thiếu sót Kính mong q Thầy Cơ dẫn thêm để học viên bổ sung kiến thức hồn thiện Luận văn Xin trân trọng cảm ơn Tp HCM, ngày tháng năm 2018 Nguyễn Văn Thơng TĨM TẮT LUẬN VĂN THẠC SĨ Hiện với phát triển ngành xây dựng hàng loạt cơng trình có chiều cao lớn, hình dáng phức tạp xây chen đời dẫn đến việc xác định tải trọng gió lên cơng trình tương đối khó khăn, tốn Vì vậy, luận văn học viên nghiên cứu đề xuất phương pháp xác định tải trọng gió, đặc biệt cho cơng trình khơng có bảng tra tiêu chuẩn Để kiểm chứng khả ứng dụng phần mềm, khả mô đánh giá gió xung quanh cơng trình nên học viên tiến hành mơ cơng trình có dạng hình trụ (CAARC) nghiên cứu trước có kết mơ thí nghiệm hầm gió Cả hai mơ áp dụng mơ hình k có điều chỉnh (mơ hình MMK) so với mơ hình standard k , mơ hình điều chỉnh đề cập nghiên cứu Murakami cộng (1997) Sau khảo sát mơ hình, mật độ lưới xác định thông số dùng mô phỏng, học viên tiếp tục áp dụng phương pháp CFD phần mềm Ansys Fluent để xác định phân tích tải trọng gió tĩnh cho cơng trình thực tế có hình dạng mặt mặt đứng phức tạp, khơng có bảng tra tiêu chuẩn thiết kế gió Từ có nhận định quan trọng phương pháp CFD tính tốn theo tiêu chuẩn xác định lực gió tĩnh tác dụng lên cơng trình LỜI CAM ĐOAN Ngoại trừ kết tham khảo từ cơng trình nghiên cứu khác trích dẫn luận văn, tơi xin cam đoan cơng việc tơi thực hướng dẫn PGS TS Nguyễn Quốc Ý TS Nguyễn Hữu Thành Kết Luận văn chưa công bố nghiên cứu khác Học viên xin chịu trách nhiệm công việc thực Tp HCM, ngày tháng năm 2018 Nguyễn Văn Thơng MỤC LỤC TÓM TẮT LUẬN VĂN THẠC SĨ LỜI CẢM ƠN LỜI CAM ĐOAN MỤC LỤC .5 DANH MỤC HÌNH ẢNH .8 DANH MỤC BẢNG BIỂU 10 DANH MỤC CHỮ VIẾT TẮT .11 DANH MỤC KÍ HIỆU 12 CHƯƠNG TỔNG QUAN 13 1.1 Đặt vấn đề 13 1.2 Tình hình nghiên cứu 19 1.3 Sự cần thiết, mục tiêu, phạm vi nghiên cứu ý nghĩa đề tài 21 1.4 Nội dung nghiên cứu 21 CHƯƠNG CƠ SỞ LÍ THUYẾT .23 2.1 Cơ sở lý thuyết CFD khí động lực học 23 Định luật bảo toàn khối lượng 23 Định luật II Newton 24 2.2 Phương pháp giải phương trình chủ đạo 25 2.3 Mơ hình Reynolds-Averaged Navier-Stokes (RANS) mơ hình MMK 25 Mơ hình Reynolds-Averaged Navier-Stokes (RANS) .27 Mơ hình MMK 28 CHƯƠNG KIẾM CHỨNG PHƯƠNG PHÁP CFD 30 3.1 Mô nghiên cứu Huang cộng (2007) 30 Miền tính tốn mật độ lưới 30 Mơ hình rối 34 Điều kiện biên .35 Kết nhận xét 37 Kết luận .41 3.2 Mô báo Huang cộng (2007) có thay đổi miền tính tốn .41 Miền tính tốn áp dụng mô 41 Kết phân tích 42 Kết luận .46 3.3 Mô nghiên cứu Kim cộng (2003) 46 Thông số ban đầu 47 Kết nhận xét 48 Kết luận .53 CHƯƠNG TÍNH TỐN TẢI TRỌNG GIĨ TĨNH CHO CƠNG TRÌNH NHÀ BẰNG MƠ PHỎNG CFD 54 4.1 Tổng quan 54 4.2 Thông số ban đầu .56 Hàm vận tốc theo phương đứng 56 Cường độ rối kích thước xốy .58 Mật độ lưới mơ hình 60 4.3 Kết nhận xét 61 Gió phương X .61 Gió phương Y .69 CHƯƠNG KẾT LUẬN VÀ KIẾN NGHỊ 76 5.1 Kết luận 76 5.2 Kiến nghị hướng phát triển 77 TÀI LIỆU THAM KHẢO .78 PHỤ LỤC 82 Phụ lục A Xác định miền tính tốn (computational domain) .82 Phụ lục B Khảo sát hệ số lực 95 Phụ lục C Kết hệ số lực tính tốn theo tiêu chuẩn phân tích CFD 101 Phụ lục D Mơ hình MMK 110 LÝ LỊCH TRÍCH NGANG 112 DANH MỤC HÌNH ẢNH Hình 1.1 The Everrich (Quận – TP.HCM) .15 Hình 1.2 Chung cư CT3 (Hoàng Mai – Hà Nội) 15 Hình 1.3 Tổ hợp cơng trình Keangnam 16 Hình 1.4 Cơng trình Platinum Park (Kuala Lumpur, Malaysia) 16 Hình 1.5 Cơng trình Shanghai Tower (Thượng Hải, Trung Quốc) 17 Hình 1.6 Hầm gió CPP 18 Hình 1.7 Mắt cắt ngang mơ hình hầm gió (Jaeyong Chung cộng sự, 2012) .18 Hình 2.1 Các phương pháp giải sử dụng CFD phần mềm .25 Hình 3.1 Miền tính tốn điều kiện biên (Huang cộng sự, 2007) 30 Hình 3.2 Miền tính tốn khảo sát 31 Hình 3.3 Hệ số lực tính tốn cho trường hợp 32 Hình 3.4 Cách chia lưới 32 Hình 3.5 Miền tính tốn dùng khảo sát chiều dài phía sau cơng trình zone (L) .33 Hình 3.6 Ảnh hưởng chiều dài L đến hệ số lực 34 Hình 3.7 Đồ thị phân bố vận tốc 35 Hình 3.8 Phân bố cường độ rối 36 Hình 3.9 Hệ số áp suất Cp chiều cao Z = 2/3 H CAARC building 38 Hình 3.10 Phân bố vận tốc gió mặt phẳng X-Y Z = 2/3H 39 Hình 3.11 Phân bố vận tốc gió mặt phẳng X-Z y = 40 Hình 3.12 Miền tính toán điều kiện biên (tác giả) 42 Hình 3.13 Hệ số áp suất Cp chiều cao Z = 2/3 H CAARC building 43 Hình 3.14 Phân bố vận tốc gió mặt phẳng X-Y Z = 2/3H 44 Hình 3.15 Phân bố vận tốc gió mặt phẳng X-Z y = 45 Hình 3.16 Hệ thống tạo lớp biên 46 Hình 3.17 Phân bố vận tốc cường độ rối (Kim cộng sự, 2003) [19] 47 Hình 3.18 Các mặt phẳng khảo sát .48 Hình 3.19 Dịng khơng khí mặt phẳng trung tâm 49 Hình 3.20 Dịng khơng khí mặt phẳng 0.5H 49 Hình 3.21 Dịng khơng khí mặt phẳng 1.0H 50 +1.88 +1.80 +1.70 +1.56 -0.07 +1.60 - - +1.88 +0.2 +2.05 +2.10 +1.66 -0.04 +1.60 - - +1.77 +2.0 +2.0 +1.89 -0.25 +2.0 +2.0 +1.54 +0.20 +1.55 +1.60 Kết Bảng B.1 thể qua Hình B.33 Hình B4 98 2.5 CFD IS 875 GS Manabu Ito Hệ số lực CD 1.5 0.5 Hình dáng tiết diện Hình B.3 Biểu đồ hệ số lực theo phương tác dụng gió 0.25 CFD 0.2 IS 875 GS Manabu Ito 0.15 Hệ số lực CL 0.1 0.05 -0.05 -0.1 -0.15 -0.2 -0.25 -0.3 Hình dáng tiết diện Hình B.4 Biểu đồ hệ số lực theo phương vng góc với phương tác dụng gió 99 B.3 Kết luận Dựa kết mơ tiết diện thể Hình B.33 Hình B, sau so sánh với kết tài liệu đưa số nhận xét sau: - Hệ số lực dọc theo phương tác động gió (IS 875 (Part 3), 2013) GS Manabu Ito cho kết chênh lệch với vài tiết diện kết chênh lệch không đáng kể - Kết hệ số CD thu từ phân tích CFD có chênh lệch với tiêu chuẩn tư liệu GS Manabu Ito chênh lệch nhỏ, chênh lệch lớn so với (IS 875 (Part 3), 2013) 11% cho tiết diện chênh lệch nhỏ 2.5% ứng với tiết diện - Hệ số lực theo phương vng góc với phương tác dụng gió xác định mô CFD cho kết khác với (IS 875 (Part 3), 2013) tư liệu GS Manabu Ito, đồng thời ứng với tiết diện đối xứng kết từ CFD cho khác khơng, cịn tiêu chuẩn tài liệu cho khơng Từ nhận định kết luận mơ CFD thu hệ số CL ≠ - Vùng ảnh hưởng xác định khả thi áp dụng vào phân tích khảo sát Luận Văn - Mô CFD có khả áp dụng vào xác định tải trọng gió tĩnh tác dụng vào cơng trình, cho kết gần tương đương với thí nghiệm, nhiên có chênh lệch kết Sự chênh lệch mơ chưa kể đến độ nhám bề mặt cơng trình so với cơng trình thực tế, thay đổi mơi trường xung quanh ứng với thời tiết thực tế, mô xét đến độ nhớt khối lượng riêng khơng khí 68 F atm 100 Phụ lục C Kết hệ số lực tính tốn theo tiêu chuẩn phân tích CFD Bảng C.1 Hệ số lực tính tốn theo TCVN 2737-1995 Tầng số (STT) ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET Chiều cao tầng (m) 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 3.2 4.5 4.5 3.5 Cao độ tầng (m) 127.2 124 120.8 117.6 114.4 111.2 108 104.8 101.6 98.4 95.2 92 88.8 85.6 82.4 79.2 76 72.8 69.6 66.4 63.2 60 56.8 53.6 50.4 47.2 44 40.8 37.6 34.4 31.2 28 24.8 21.6 18.4 15.2 12 4.5 Bề rộng đón gió theo phương Tải trọng gió tĩnh quy lực tập trung Phương X Phương Y Phương X Phương Y (m) 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 (m) 26.4 33 40.05 41.6 46.65 49.95 53.25 56.83 56.83 60.4 60.4 63.98 63.98 67.55 67.55 71.13 71.13 74.7 74.7 78.28 78.28 81.85 81.85 85.43 85.43 89 89 89 89 89 89 89 89 89 89 89 89 89 89 89 (daN) 16436.667 20467.5 24745 25603.333 28600.833 30505 32393.333 34433.333 34297.5 36294.167 36114.167 38061.667 37870.833 39786.667 39585.833 41458.333 41205 43009.167 42742.5 44508.333 44229.167 45956.667 45567.5 47150.833 46744.167 48277.5 47853.333 47430 46847.5 46212.5 45576.667 44809.167 43961.667 43114.167 42161.667 41102.5 37354.167 53890.833 48981.667 36196.667 (daN) 13075 13025 12975 12925 12875 12825 12775 12725 12675 12619.167 12556.667 12494.167 12431.667 12369.167 12306.667 12240.833 12165.833 12090.833 12015.833 11940.833 11865.833 11790.833 11690.833 11590.833 11490.833 11390.833 11291.667 11191.667 11054.167 10904.167 10754.167 10572.5 10373.333 10173.333 9948.3333 9698.3333 8814.1667 12715.833 11557.5 8540.8333 101 Hệ số áp lực phân bố theo chiều cao Phương X Phương Y 2.34 2.34 2.33 2.32 2.31 2.30 2.29 2.28 2.27 2.26 2.25 2.24 2.23 2.22 2.21 2.19 2.18 2.17 2.15 2.14 2.13 2.11 2.10 2.08 2.06 2.04 2.02 2.01 1.98 1.95 1.93 1.90 1.86 1.82 1.78 1.74 1.69 1.62 1.47 0.00 2.34 2.34 2.33 2.32 2.31 2.30 2.29 2.28 2.27 2.26 2.25 2.24 2.23 2.22 2.21 2.19 2.18 2.17 2.15 2.14 2.13 2.11 2.10 2.08 2.06 2.04 2.02 2.01 1.98 1.95 1.93 1.90 1.86 1.82 1.78 1.74 1.69 1.62 1.47 0.00 Bảng C.2 Hệ số lực từ phân tích CFD theo phương X block D Wind X_Trường hợp (CFD) LEVEL ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET Fd (N) Fl (N) h(m) z (m) B(m) v (z) (m/s) W0 (daN/m2) Cd Cl 0.02239 0.04069 0.05297 0.06219 0.06901 0.07654 0.08306 0.09015 0.09195 0.09742 0.09989 0.10401 0.10516 0.10960 0.11083 0.11435 0.11331 0.11572 0.11574 0.11809 0.11787 0.12054 0.12029 0.12385 0.12331 0.12157 0.12013 0.11706 0.11469 0.11205 0.10901 0.10565 0.10259 0.09938 0.09635 0.09271 0.08685 0.11609 0.13308 0.00000 -0.00827 -0.01053 -0.00904 -0.00930 -0.00977 -0.00893 -0.00966 -0.00867 -0.00944 -0.00824 -0.00826 -0.00658 -0.00774 -0.00616 -0.00754 -0.00553 -0.00636 -0.00543 -0.00638 -0.00455 -0.00581 -0.00520 -0.00610 -0.00158 -0.00184 -0.00220 -0.00315 -0.00341 -0.00369 -0.00307 -0.00318 -0.00295 -0.00326 -0.00309 -0.00304 -0.00311 -0.00302 -0.00038 -0.00245 0.00000 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.0075 0.01125 0.01125 0.00875 0.318 0.31 0.302 0.294 0.286 0.278 0.27 0.262 0.254 0.246 0.238 0.23 0.222 0.214 0.206 0.198 0.19 0.182 0.174 0.166 0.158 0.15 0.142 0.134 0.126 0.118 0.11 0.102 0.094 0.086 0.078 0.07 0.062 0.054 0.046 0.038 0.03 0.0225 0.01125 0.066 0.0825 0.100125 0.104 0.116625 0.124875 0.133125 0.142075 0.142075 0.151 0.151 0.15995 0.15995 0.168875 0.168875 0.177825 0.177825 0.18675 0.18675 0.1957 0.1957 0.204625 0.204625 0.213575 0.213575 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 10 0.76 1.11 1.20 1.36 1.35 1.40 1.44 1.47 1.50 1.50 1.55 1.53 1.55 1.54 1.57 1.55 1.54 1.51 1.52 1.49 1.50 1.47 1.48 1.47 1.48 1.42 1.41 1.39 1.38 1.37 1.35 1.33 1.31 1.30 1.29 1.28 1.32 1.26 1.68 0.00 11 -0.28 -0.29 -0.20 -0.20 -0.19 -0.16 -0.17 -0.14 -0.15 -0.13 -0.13 -0.10 -0.11 -0.09 -0.11 -0.07 -0.09 -0.07 -0.08 -0.06 -0.07 -0.06 -0.08 -0.02 -0.02 -0.03 -0.04 -0.04 -0.04 -0.04 -0.04 -0.04 -0.04 -0.04 -0.04 -0.04 -0.05 0.00 -0.03 0.00 102 Bảng C.3 Hệ số lực từ phân tích CFD theo phương X block D Wind X_Trường hợp (CFD) ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET Fd (N) Fl (N) h(m) z (m) B(m) v (z) (m/s) W0 (daN/m2) Cd Cl 0.023639 0.040788 0.053897 0.061675 0.068928 0.075269 0.08068 0.087138 0.089848 0.093906 0.095471 0.100072 0.101359 0.104837 0.106363 0.109682 0.110805 0.114504 0.115771 0.119525 0.119889 0.124469 0.126221 0.130426 0.131702 0.130696 0.130149 0.128637 0.12763 0.125675 0.123262 0.120191 0.118006 0.115484 0.114608 0.112652 0.10921 0.149858 0.16846 -0.00275 -0.00492 -0.00297 -0.00339 -0.00236 -0.00353 -0.00447 -0.00321 -0.00539 -0.00368 -0.00513 -0.00412 -0.00397 -0.00379 -0.00438 -0.00179 -0.00304 -0.00159 -0.0038 -0.00271 -0.00201 -0.00274 -0.00173 0.001095 0.000597 0.001482 0.000647 0.001167 0.001553 0.001353 0.001162 0.000899 0.00025 0.000324 0.00027 -0.00134 -0.00268 0.000207 -0.00439 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.0075 0.01125 0.01125 0.00875 0.318 0.31 0.302 0.294 0.286 0.278 0.27 0.262 0.254 0.246 0.238 0.23 0.222 0.214 0.206 0.198 0.19 0.182 0.174 0.166 0.158 0.15 0.142 0.134 0.126 0.118 0.11 0.102 0.094 0.086 0.078 0.07 0.062 0.054 0.046 0.038 0.03 0.0225 0.01125 0.066 0.0825 0.100125 0.104 0.116625 0.124875 0.133125 0.142075 0.142075 0.151 0.151 0.15995 0.15995 0.168875 0.168875 0.177825 0.177825 0.18675 0.18675 0.1957 0.1957 0.204625 0.204625 0.213575 0.213575 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 0.80 1.12 1.22 1.35 1.35 1.38 1.39 1.42 1.47 1.45 1.48 1.47 1.50 1.48 1.51 1.48 1.51 1.49 1.52 1.51 1.52 1.52 1.56 1.55 1.58 1.52 1.53 1.53 1.54 1.53 1.52 1.51 1.51 1.51 1.53 1.55 1.67 1.62 2.12 0.00 10 -0.09 -0.13 -0.07 -0.07 -0.05 -0.06 -0.08 -0.05 -0.09 -0.06 -0.08 -0.06 -0.06 -0.05 -0.06 -0.02 -0.04 -0.02 -0.05 -0.03 -0.03 -0.03 -0.02 0.01 0.01 0.02 0.01 0.01 0.02 0.02 0.01 0.01 0.00 0.00 0.00 -0.02 -0.04 0.00 -0.06 0.00 103 Wind X_Trường hợp (CFD) ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET FD (N) FL (N) h(m) z (m) B(m) v (z) (m/s) W0 (daN/m2) CD CL 0.0240 0.0379 0.0510 0.0611 0.0670 0.0735 0.0801 0.0866 0.0888 0.0949 0.0970 0.1019 0.1037 0.1091 0.1107 0.1148 0.1168 0.1201 0.1221 0.1259 0.1273 0.1319 0.1328 0.1399 0.1411 0.1408 0.1401 0.1394 0.1377 0.1364 0.1341 0.1330 0.1311 0.1281 0.1267 0.1250 0.1220 0.1610 0.1814 -0.0024 -0.0022 -0.0013 -0.0040 -0.0027 -0.0032 -0.0034 -0.0022 -0.0029 -0.0026 -0.0042 -0.0039 -0.0033 -0.0026 -0.0041 -0.0025 -0.0032 -0.0033 -0.0036 -0.0039 -0.0031 -0.0025 -0.0021 -0.0022 -0.0030 -0.0026 -0.0024 -0.0024 -0.0025 -0.0033 -0.0024 -0.0026 -0.0034 -0.0047 -0.0047 -0.0052 -0.0053 -0.0070 -0.0073 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.0075 0.01125 0.01125 0.00875 0.318 0.31 0.302 0.294 0.286 0.278 0.27 0.262 0.254 0.246 0.238 0.23 0.222 0.214 0.206 0.198 0.19 0.182 0.174 0.166 0.158 0.15 0.142 0.134 0.126 0.118 0.11 0.102 0.094 0.086 0.078 0.07 0.062 0.054 0.046 0.038 0.03 0.0225 0.01125 0.066 0.0825 0.100125 0.104 0.116625 0.124875 0.133125 0.142075 0.142075 0.151 0.151 0.15995 0.15995 0.168875 0.168875 0.177825 0.177825 0.18675 0.18675 0.1957 0.1957 0.204625 0.204625 0.213575 0.213575 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 10 0.82 1.04 1.15 1.34 1.31 1.35 1.38 1.41 1.45 1.46 1.50 1.50 1.53 1.54 1.57 1.55 1.59 1.57 1.60 1.59 1.62 1.61 1.64 1.66 1.70 1.64 1.65 1.66 1.66 1.66 1.66 1.67 1.68 1.67 1.70 1.72 1.86 1.74 2.29 0.00 11 -0.08 -0.06 -0.03 -0.09 -0.05 -0.06 -0.06 -0.04 -0.05 -0.04 -0.07 -0.06 -0.05 -0.04 -0.06 -0.03 -0.04 -0.04 -0.05 -0.05 -0.04 -0.03 -0.03 -0.03 -0.04 -0.03 -0.03 -0.03 -0.03 -0.04 -0.03 -0.03 -0.04 -0.06 -0.06 -0.07 -0.08 -0.08 -0.09 0.00 104 Bảng C.3 Hệ số lực từ phân tích CFD theo phương X block D Wind X_Trường hợp (CFD) LEVEL ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET Fd (N) Fl (N) h(m) z (m) B(m) v (z) (m/s) W0 (daN/m2) Cd Cl 0.12278 0.16468 0.20840 0.22183 0.23103 0.24171 0.25297 0.27153 0.26645 0.27570 0.26402 0.28042 0.26556 0.27183 0.25581 0.25965 0.24079 0.24079 0.22725 0.23339 0.21660 0.21849 0.20905 0.22418 0.21758 0.21223 0.20829 0.20046 0.19756 0.19429 0.19213 0.18943 0.19069 0.19178 0.19293 0.19379 0.18867 0.24649 0.26448 0.00000 -0.04596 -0.05599 -0.04938 -0.04634 -0.05392 -0.05316 -0.06464 -0.02646 -0.03534 -0.02736 -0.04250 -0.01906 -0.03744 -0.01701 -0.04625 -0.01998 -0.02892 -0.01868 -0.02843 -0.01691 -0.02627 -0.03150 -0.03472 -0.02085 -0.01965 -0.02008 -0.02014 -0.02041 -0.02005 -0.02025 -0.02009 -0.02014 -0.02006 -0.01963 -0.01919 -0.01879 -0.01766 -0.02526 -0.02278 0.00000 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.0075 0.01125 0.01125 0.00875 0.318 0.31 0.302 0.294 0.286 0.278 0.27 0.262 0.254 0.246 0.238 0.23 0.222 0.214 0.206 0.198 0.19 0.182 0.174 0.166 0.158 0.15 0.142 0.134 0.126 0.118 0.11 0.102 0.094 0.086 0.078 0.07 0.062 0.054 0.046 0.038 0.03 0.0225 0.01125 0.066 0.0825 0.100125 0.104 0.116625 0.124875 0.133125 0.142075 0.142075 0.151 0.151 0.15995 0.15995 0.168875 0.168875 0.177825 0.177825 0.18675 0.18675 0.1957 0.1957 0.204625 0.204625 0.213575 0.213575 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 0.2225 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 10 4.18 4.50 4.71 4.85 4.52 4.43 4.37 4.41 4.35 4.26 4.09 4.13 3.93 3.83 3.62 3.51 3.27 3.14 2.98 2.94 2.75 2.67 2.58 2.67 2.61 2.47 2.45 2.38 2.38 2.37 2.38 2.38 2.44 2.51 2.58 2.67 2.88 2.67 3.34 0.00 11 -1.57 -1.53 -1.12 -1.01 -1.05 -0.98 -1.12 -0.43 -0.58 -0.42 -0.66 -0.28 -0.55 -0.24 -0.65 -0.27 -0.39 -0.24 -0.37 -0.21 -0.33 -0.39 -0.43 -0.25 -0.24 -0.23 -0.24 -0.24 -0.24 -0.25 -0.25 -0.25 -0.26 -0.26 -0.26 -0.26 -0.27 -0.27 -0.29 0.00 105 Bảng C.4 Hệ số lực từ phân tích CFD theo phương Y block D Wind Y_Trường hợp (CFD) Fd (N) ROOF 0.012206 T37 0.009252 T36 0.013736 T35 0.020163 T34 0.016849 T33 0.022788 T32 0.015958 T31 0.029083 T30 0.024679 T29 0.029909 T28 0.023911 T27 0.029104 T26 0.026116 T25 0.029286 T24 0.027301 T23 0.029057 T22 0.027345 T21 0.029111 T20 0.026003 T19 0.027573 T18 0.027479 T17 0.024664 T16 0.020494 T15 0.027888 T14 0.030057 T13 0.029896 T12 0.02852 T11 0.027492 T10 0.026717 T9 0.026441 T8 0.025982 T7 0.02608 T6 0.024639 T5 0.02393 T4 0.023057 T3 0.021338 T2 0.019965 T.KT 0.018701 T.LUNG 0.019512 T.TRET Fl (N) h(m) z (m) B(m) -0.004 -0.0073 -0.01008 -0.00094 -0.01243 -0.01184 -0.01599 -0.01451 -0.01356 -0.01522 -0.01457 -0.017 -0.0143 -0.01893 -0.01711 -0.0197 -0.01766 -0.02022 -0.01879 -0.02033 -0.01795 -0.02139 -0.01974 -0.02218 -0.02143 -0.02036 -0.01998 -0.01959 -0.02041 -0.01917 -0.01891 -0.01876 -0.02007 -0.01997 -0.01847 -0.01854 -0.01706 -0.02106 -0.02574 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.01125 0.01125 0.00875 0.318 0.310 0.302 0.294 0.286 0.278 0.270 0.262 0.254 0.246 0.238 0.230 0.222 0.214 0.206 0.198 0.190 0.182 0.174 0.166 0.158 0.150 0.142 0.134 0.126 0.118 0.110 0.102 0.094 0.086 0.078 0.070 0.062 0.054 0.046 0.038 0.030 0.02250 0.01125 0.00000 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 106 v (z)(m/s) 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 W0 (daN/m2) CD CL 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 0.52 0.40 0.59 0.87 0.73 0.99 0.70 1.28 1.09 1.33 1.07 1.30 1.18 1.33 1.24 1.33 1.26 1.35 1.21 1.29 1.30 1.18 0.98 1.35 1.47 1.47 1.42 1.39 1.36 1.37 1.36 1.39 1.34 1.32 1.31 1.25 1.29 0.86 1.04 0.00 10 -0.17 -0.31 -0.43 -0.04 -0.54 -0.52 -0.70 -0.64 -0.60 -0.68 -0.65 -0.76 -0.64 -0.86 -0.78 -0.90 -0.81 -0.94 -0.88 -0.95 -0.85 -1.02 -0.95 -1.07 -1.05 -1.00 -1.00 -0.99 -1.04 -0.99 -0.99 -1.00 -1.09 -1.11 -1.05 -1.08 -1.10 -0.97 -1.38 0.00 Bảng C.5 Hệ số lực từ phân tích CFD theo phương Y block D Wind Y_Trường hợp (CFD) ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET FD (N) FL (N) h(m) z (m) B(m) 0.009404 0.005887 0.017597 0.018179 0.021227 0.023362 0.018834 0.027247 0.027859 0.030179 0.026002 0.030415 0.025816 0.03073 0.025422 0.028153 0.026541 0.026671 0.021978 0.024647 0.023839 0.019144 0.016234 0.023148 0.024525 0.024715 0.024462 0.025472 0.025817 0.025666 0.025605 0.025469 0.02566 0.025356 0.02435 0.023836 0.022367 0.024604 0.023399 -0.02382 -0.03466 -0.04459 -0.05091 -0.05731 -0.06369 -0.07038 -0.07774 -0.07572 -0.08434 -0.084 -0.09427 -0.09258 -0.10179 -0.1006 -0.11078 -0.10842 -0.1185 -0.1163 -0.12567 -0.12103 -0.13005 -0.12607 -0.14068 -0.13687 -0.13658 -0.13441 -0.13579 -0.13441 -0.13492 -0.13408 -0.13482 -0.13448 -0.135 -0.13467 -0.13407 -0.12889 -0.16525 -0.18248 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.0075 0.01125 0.01125 0.00875 0.318 0.31 0.302 0.294 0.286 0.278 0.27 0.262 0.254 0.246 0.238 0.23 0.222 0.214 0.206 0.198 0.19 0.182 0.174 0.166 0.158 0.15 0.142 0.134 0.126 0.118 0.11 0.102 0.094 0.086 0.078 0.07 0.062 0.054 0.046 0.038 0.03 0.0225 0.01125 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 107 v (z)(m/s) 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 W0 (daN/m2) CD CL 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 0.40253 0.252909 0.758781 0.786847 0.9224 1.019261 0.825126 1.198827 1.231181 1.3398 1.159817 1.36327 1.162959 1.391621 1.15752 1.289141 1.222508 1.236103 1.025154 1.157434 1.127455 0.912178 0.779674 1.121042 1.198333 1.219143 1.219012 1.283382 1.31637 1.325864 1.341833 1.356283 1.391371 1.403782 1.381348 1.392823 1.44554 1.127722 1.250814 10 -1.01983 -1.48911 -1.92274 -2.20382 -2.49041 -2.77896 -3.08322 -3.42054 -3.34614 -3.74423 -3.74661 -4.22521 -4.17075 -4.60945 -4.58063 -5.07286 -4.99373 -5.49191 -5.42486 -5.90161 -5.7242 -6.19668 -6.0548 -6.8132 -6.68779 -6.73753 -6.69823 -6.8419 -6.85319 -6.96943 -7.02624 -7.17942 -7.2921 -7.47382 -7.63955 -7.83433 -8.3296 -7.57413 -9.75452 Wind Y_Trường hợp (CFD) ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET FD (N) FL (N) h(m) z (m) B(m) 0.0061 0.0070 0.0092 0.0110 0.0108 0.0105 0.0097 0.0110 0.0102 0.0103 0.0089 0.0103 0.0090 0.0110 0.0101 0.0113 0.0105 0.0115 0.0117 0.0114 0.0120 0.0116 0.0107 0.0086 0.0070 0.0055 0.0038 0.0020 0.0016 0.0013 0.0006 0.0001 0.0000 0.0002 0.0008 0.0011 0.0016 0.0030 0.0033 -0.0225 -0.0293 -0.0379 -0.0430 -0.0454 -0.0498 -0.0539 -0.0581 -0.0548 -0.0600 -0.0574 -0.0635 -0.0615 -0.0662 -0.0642 -0.0704 -0.0667 -0.0724 -0.0693 -0.0758 -0.0725 -0.0783 -0.0738 -0.0833 -0.0798 -0.0759 -0.0744 -0.0730 -0.0708 -0.0700 -0.0675 -0.0645 -0.0619 -0.0593 -0.0567 -0.0540 -0.0499 -0.0613 -0.0654 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.0075 0.01125 0.01125 0.00875 0.318 0.31 0.302 0.294 0.286 0.278 0.27 0.262 0.254 0.246 0.238 0.23 0.222 0.214 0.206 0.198 0.19 0.182 0.174 0.166 0.158 0.15 0.142 0.134 0.126 0.118 0.11 0.102 0.094 0.086 0.078 0.07 0.062 0.054 0.046 0.038 0.03 0.0225 0.01125 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 108 v (z)(m/s) 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 W0 (daN/m2) CD CL 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 0.259 0.303 0.395 0.474 0.471 0.459 0.425 0.485 0.450 0.458 0.395 0.463 0.403 0.500 0.460 0.518 0.484 0.533 0.546 0.535 0.567 0.551 0.516 0.415 0.340 0.273 0.188 0.101 0.081 0.068 0.030 0.008 0.000 0.010 0.045 0.064 0.103 0.136 0.178 10 -0.965 -1.258 -1.633 -1.860 -1.974 -2.172 -2.362 -2.558 -2.421 -2.665 -2.559 -2.846 -2.772 -2.998 -2.924 -3.223 -3.073 -3.353 -3.233 -3.559 -3.427 -3.729 -3.545 -4.033 -3.899 -3.745 -3.710 -3.677 -3.612 -3.617 -3.536 -3.435 -3.355 -3.280 -3.219 -3.155 -3.225 -2.808 -3.496 Bảng C.6 Hệ số lực từ phân tích CFD theo phương Y block D Wind Y_Trường hợp (CFD) ROOF T37 T36 T35 T34 T33 T32 T31 T30 T29 T28 T27 T26 T25 T24 T23 T22 T21 T20 T19 T18 T17 T16 T15 T14 T13 T12 T11 T10 T9 T8 T7 T6 T5 T4 T3 T2 T.KT T.LUNG T.TRET FD (N) FL (N) h(m) z (m) B(m) 0.010215 0.009444 0.01123 0.012346 0.011719 0.010523 0.009461 0.011169 0.010615 0.011527 0.009976 0.012105 0.010917 0.012105 0.011493 0.014067 0.013271 0.016387 0.014996 0.016749 0.016276 0.016436 0.016257 0.011973 0.002268 -9.70E-05 0.001749 0.002126 0.002239 0.002421 0.00257 0.002761 0.002876 0.003024 0.003292 0.003492 0.003735 0.005128 0.00648 -0.01303 -0.01703 -0.02014 -0.02159 -0.02345 -0.0249 -0.02591 -0.02674 -0.02548 -0.02742 -0.0258 -0.0283 -0.02662 -0.02849 -0.02708 -0.02862 -0.02785 -0.02901 -0.02811 -0.02981 -0.02903 -0.03034 -0.02827 -0.02892 -0.02556 -0.02347 -0.02191 -0.0214 -0.02055 -0.02016 -0.01975 -0.01916 -0.01882 -0.01817 -0.01755 -0.01642 -0.01419 -0.01772 -0.01799 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.01125 0.01125 0.00875 0.318 0.310 0.302 0.294 0.286 0.278 0.270 0.262 0.254 0.246 0.238 0.230 0.222 0.214 0.206 0.198 0.190 0.182 0.174 0.166 0.158 0.150 0.142 0.134 0.126 0.118 0.110 0.102 0.094 0.086 0.078 0.070 0.062 0.054 0.046 0.038 0.030 0.02250 0.01125 0.00000 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 0.0525 109 v (z)(m/s) 9.53 9.51 9.49 9.48 9.46 9.44 9.42 9.40 9.38 9.36 9.34 9.31 9.29 9.26 9.24 9.21 9.19 9.16 9.13 9.10 9.07 9.03 9.00 8.96 8.92 8.88 8.83 8.78 8.73 8.67 8.61 8.54 8.47 8.38 8.28 8.16 8.01 7.77 7.19 7.07 W0 (daN/m2) CD CL 55.62 55.42 55.22 55.01 54.79 54.57 54.35 54.11 53.88 53.63 53.38 53.12 52.85 52.58 52.29 52.00 51.69 51.37 51.04 50.70 50.34 49.97 49.58 49.16 48.73 48.27 47.78 47.26 46.70 46.09 45.43 44.71 43.91 43.01 41.97 40.75 39.30 36.94 31.67 30.58 0.44 0.41 0.48 0.53 0.51 0.46 0.41 0.49 0.47 0.51 0.44 0.54 0.49 0.55 0.52 0.64 0.61 0.76 0.70 0.79 0.77 0.78 0.78 0.58 0.11 0.00 0.09 0.11 0.11 0.13 0.13 0.15 0.16 0.17 0.19 0.20 0.24 0.24 0.35 0.00 10 -0.56 -0.73 -0.87 -0.93 -1.02 -1.09 -1.14 -1.18 -1.13 -1.22 -1.15 -1.27 -1.20 -1.29 -1.23 -1.31 -1.28 -1.34 -1.31 -1.40 -1.37 -1.45 -1.36 -1.40 -1.25 -1.16 -1.09 -1.08 -1.05 -1.04 -1.03 -1.02 -1.02 -1.01 -1.00 -0.96 -0.92 -0.81 -0.96 0.00 Phụ lục D Mơ hình MMK Đoạn code tham khảo nghiên cứu Huang cộng (2007) #include "udf.h" DEFINE_TURBULENT_VISCOSITY(NMK_mu_t, c, t) { real mu_t,S,O; real Sxx,Sxy,Sxz; /*Sij*/ real Syx,Syy,Syz; real Szx,Szy,Szz; real Oxx,Oxy,Oxz; /*Vorticity tensor*/ real Oyx,Oyy,Oyz; real Ozx,Ozy,Ozz; real rho = C_R(c,t); /*density*/ real k = C_K(c,t); /*k*/ real d = C_D(c,t); /*e*/ Sxx=0.5*(C_DUDX(c,t)+C_DUDX(c,t)); Sxy=0.5*(C_DVDX(c,t)+C_DUDY(c,t)); Sxz=0.5*(C_DWDX(c,t)+C_DUDZ(c,t)); Syx=0.5*(C_DUDY(c,t)+C_DVDX(c,t)); Syy=0.5*(C_DVDY(c,t)+C_DVDY(c,t)); Syz=0.5*(C_DWDY(c,t)+C_DVDZ(c,t)); Szx=0.5*(C_DUDZ(c,t)+C_DWDX(c,t)); Szy=0.5*(C_DVDZ(c,t)+C_DWDY(c,t)); Szz=0.5*(C_DWDZ(c,t)+C_DWDZ(c,t)); S=sqrt(2*(Sxx*Sxx+Sxy*Sxy+Sxz*Sxz+Syx*Syx+Syy *Syy+Syz*Syz+Szx*Szx+Szy*Szy+Szz*Szz)); Oxx=0.5*(C_DUDX(c,t)-C_DUDX(c,t)); Oxy=0.5*(C_DVDX(c,t)-C_DUDY(c,t)); Oxz=0.5*(C_DWDX(c,t)-C_DUDZ(c,t)); Oyx=0.5*(C_DUDY(c,t)-C_DVDX(c,t)); Oyy=0.5*(C_DVDY(c,t)-C_DVDY(c,t)); 110 Oyz=0.5*(C_DWDY(c,t)-C_DVDZ(c,t)); Ozx=0.5*(C_DUDZ(c,t)-C_DWDX(c,t)); Ozy=0.5*(C_DVDZ(c,t)-C_DWDY(c,t)); Ozz=0.5*(C_DWDZ(c,t)-C_DWDZ(c,t)); O=sqrt(2*(Oxx*Oxx+Oxy*Oxy+Oxz*Oxz+Oyx*Oyx+Oyy *Oyy+Oyz*Oyz+Ozx*Ozx+Ozy*Ozy+Ozz*Ozz)); if((O/S-1)