Luận Văn Thạc Sĩ Kỹ Thuật Xây Dựng.pdf

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ CÔNG NGHIỆP LONG AN MAI LƯU QUỐC ĐẠT ỨNG DỤNG CHƯ NG TR NH PH N T CH SỐ TRONG THIẾT KẾ M THANH CHO HỘI TRƯỜNG LUẬN VĂN THẠC SĨ Chuyên ngành Kỹ Thuật Xây D[.]

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ CÔNG NGHIỆP LONG AN MAI LƯU QUỐC ĐẠT ỨNG DỤNG CHƯ NG TR NH PH N T CH SỐ TRONG THIẾT KẾ M THANH CHO HỘI TRƯỜNG LUẬN VĂN THẠC SĨ Chuyên ngành Kỹ Thuật Xây Dựng Long An, năm 2019 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ CÔNG NGHIỆP LONG AN - MAI LƯU QUỐC ĐẠT ỨNG DỤNG CHƯ NG TR NH PH N T CH SỐ TRONG THIẾT KẾ M THANH CHO HỘI TRƯỜNG DIGITAL ANALYSIS PROGRAM APPLICATION IN THE SOUND DESIGN FOR THE HALL LUẬN VĂN THẠC SĨ Chuyên ngành Kỹ Thuật Xây Dựng Mã ngành: 8.58.02.01 Người hướng dẫn khoa học: PGS TS TRƯ NG T CH THIỆN Long An, năm 2019 i LỜI CAM ĐOAN Ngoài kết tham khảo từ cơng trình khác nhƣ đƣợc ghi luận văn, xin cam đoan luận văn tơi thực luận văn đƣợc nộp Trƣờng Đại học Kinh tế Công nghiệp Long An Số liệu kết nghiên cứu luận văn hoàn toàn trung thực chƣa đƣợc sử dụng công bố cơng trình khác Mọi giúp đỡ cho việc thực luận văn đƣợc cảm ơn thơng tin trích dẫn luận văn đƣợc ghi rõ nguồn góc HỌC VIÊN THỰC HIỆN Mai Lưu Quốc Đạt ii LỜI CẢM N Luận văn Cao học hồn thành kết q trình học tập nghiên cứu học viên trƣờng Đại học kinh tế Công nghiệp Long An Bên cạnh nỗ lực thân học viên, việc hồn thành chƣơng trình luận văn khơng thể thiếu giảng dạy, quan tâm, giúp đỡ tập thể quý thầy, cô Khoa Kiến trúc Xây dựng ( Trƣờng Đại học Kinh tế Công nghiệp Long An) trình học tập nhƣ hồn thành luận văn cao học Tôi xin chân thành cảm ơn thầy giáo hƣớng dẫn: PGS.TS Trương Tích Thiện tập thể q thầy, cơ, đồng nghiệp tận tình quan tâm hƣớng dẫn truyền đạt kiến thức kinh nghiệp, tạo điếu kiện thuận lợi để tơi hồn thành luận văn cao học HỌC VIÊN THỰC HIỆN Mai Lưu Quốc Đạt iii MỤC LỤC MỤC LỤC iii DANH MỤC TỪ VIẾT TẮT TIẾNG VIỆT v DANH MỤC VIẾT TẮT TIẾNG ANH v DANH MỤC BẢNG BIỂU vi DANH MỤC ĐỒ THỊ VÀ HÌNH VẼ vi CHƢƠNG TỔNG QUAN VỀ ĐỀ TÀI 1.1 Tính cần thiết đề tài 1.2 Mục tiêu nghiên cứu 1.2.1 Mục tiêu chung 1.2.2 Mục tiêu cụ thể 1.3 Đối tƣợng nghiên cứu 1.4 Phạm vi nghiên cứu 1.5 Câu hỏi nghiên cứu 1.6 Phƣơng pháp nghiên cứu CHƢƠNG CƠ SỞ LÝ THUYẾT 2.1 Lý thuyết âm 2.1.1 Các phƣơng trình chủ đạo 2.1.2 Cơ sở phƣơng trình sóng phần tử hữu hạn 2.1.3 Phƣơng trình chủ đạo với ảnh hƣởng lƣu lƣợng trung bình 2.1.4 Cơ sở phần tử hữu hạn cho phƣơng trình sóng đối lƣu 2.1.5 Ma trận truyền âm 10 2.2 Bài toán truyền âm .12 2.2.1 Điều kiện biên âm 12 2.2.2 Điều kiện biên hấp thu 14 2.3 Thơng số đầu tốn truyền âm 15 2.3.1 Áp lực âm 15 2.3.2 Công suất âm 17 2.5.5 Mức áp suất âm 18 2.4 Âm học ngẫu nhiên - Khuếch tán trƣờng âm (Acoustic Diffuse Sound Field) 18 2.5 Âm phòng 20 2.5.1 Mơ hình khuếch tán cho âm phòng 20 2.5.2 Điều kiện biên 21 2.5.3 Tiêu chuẩn âm 21 iv 2.5.4 Sự kết hợp phòng đƣợc ngăn tƣờng 22 2.6 Âm hội trƣờng .23 2.6.1 Đánh giá chất lƣợng âm hội trƣờng theo chủ quan 23 2.6.2 Đánh giá chất lƣợng âm hội trƣờng theo khác quan 23 2.6.3 Áp dụng nguyên lý âm hình học để thiết kế hình dạng phịng 24 2.6.4 Tránh tƣợng xấu âm học 24 2.6.5 Thiết kế tạo trƣờng khuếch tán âm 25 2.6.6 Thiết kế hội trƣờng theo thời gian âm vang 26 2.6.7 Công thức xác định thời gian âm vang 26 2.6.8 Thời gian âm vang tối ƣu 27 CHƢƠNG CÁC MƠ HÌNH TÍNH TỐN 29 3.1 Phân tích tốn truyền âm chƣơng trình ANSYS 29 3.1.1 Bài toán hấp thụ âm truyền qua 29 3.1.2 Bài toán truyền âm hai phòng theo tiêu chuẩn ASTM [6] 34 3.2 Phân tích tốn truyền âm chƣơng trình Odeon 36 3.2.1 Mơ hình toán 37 3.2.2 Thông số vật liệu 37 3.2.3 Điều kiện biên 42 3.2.4 Kết tính tốn phần mềm Odeon 43 3.2.5 Kết tính tốn theo lý thuyết 48 3.3 Phân tích tốn truyền âm chƣơng trình COMSOL .49 3.3.1 Mơ hình toán 49 3.3.2 Điều kiện biên 51 3.3.3 Kết mô 54 3.4 Kết luận 55 CHƢƠNG KẾT LUẬN – KIẾN NGHỊ .57 4.1 Kết luận 57 4.1.1 Những ƣu điểm luận văn 57 4.1.2 Những thiếu sót luận văn 57 4.2 Kiến nghị - hƣớng phát triển 58 DANH MỤC TÀI LIỆU THAM KHẢO 58 v DANH MỤC TỪ VIẾT TẮT TIẾNG VIỆT STT TỪ VIẾT TẮT VIẾT ĐẦY ĐỦ VN Việt Nam PTHH Phần tử hữu hạn TCNV Tiêu chuẩn Việt Nam DANH MỤC VIẾT TẮT TIẾNG ANH STT TỪ VIẾT TẮT VIẾT ĐẦY ĐỦ vi DANH MỤC BẢNG BIỂU Bảng 3.1 Thông tin vật liệu thƣ viện Odeon 37 Bảng 3.2 Thông số chi tiết vật liệu sử dụng mơ hình 38 Bảng 3.3 Tọa độ nguồn âm 43 Bảng 3.4 Tọa độ nơi thu âm 43 Bảng 3.5 Kích thƣớc phịng 50 DANH MỤC ĐỒ THỊ VÀ H NH VẼ Hình 1.1 Nhà hát Opera Paris – Garnier Hình 2.1 Điều kiện hấp thu 15 Hình 3.1 Mơ tả ống truyền âm ngăn mỏng có đục lỗ với điều kiện biên [3] 29 Hình 3.2 Mơ tả bề dày panel đƣờng kính lỗ phần lỗ đƣợc khoét 30 Hình 3.3 Hệ số hấp thu âm theo tần số số mẫu vật liệu phƣơng pháp đo lƣờng thực nghiệm 31 Hình 3.4 Số lƣợng lỗ panel 31 Hình 3.5 Chia lƣới mơ hình ansys 32 Hình 3.6 Khai báo toàn điều kiện biên điều kiện đầu vào ANSYS 32 Hình 3.7 Kết mô mức áp suất âm hay mức cƣờng độ âm ống ANSYS 32 Hình 3.8 So sánh kết thực nghiệm kết Ansys hệ số hấp thụ âm 33 Hình 3.9 Đánh giá sai số hệ số hấp thụ âm thực nghiệm tính tốn 33 Hình 3.10 Căn phòng thực nghiệm đo mát truyền âm [4] 34 Hình 3.11 Mơ hình phịng đƣợc xây dựng ANSYS 34 Hình 3.12 Kết mesh lƣới mơ hình size 0.05m 35 Hình 3.13 Mơ tả điều kiện đầu vào cho toán 35 Hình 3.14 Sự phân bố mức cƣờng độ âm phòng 36 vii Hình 3.15 4.1 Các bƣớc tiến hành Odeon 36 Hình 3.16 Mơ hình khán phịng dùng tốn 37 Hình 3.17 Hệ số hấp thụ vật liệu 42 Hình 3.18 Vị trí nguồn âm nơi thu âm 43 Hình 3.19 Kết tính tốn Odeon 44 Hình 3.20 Thời gian hồi âm hệ thống 44 Hình 3.21 Thời gian hồi âm T30 T20 45 Hình 3.22 Mức áp suất âm theo dải tần số 45 Hình 3.23 Các giá trị thu đƣợc R01 45 Hình 3.24 Các giá trị thu đƣợc R02 46 Hình 3.25 Thời gian hồi âm trung bình T30 46 Hình 3.26 Các giá trị T30 theo dải tần số 46 Hình 3.27 Mức áp suất âm theo mạng cân A, tuyến tính mạng cân C47 Hình 3.28 Kết tính toán theo lý thuyết 48 Hình 3.29 So sánh kết tính tốn theo lý thuyết kết theo Odeon 49 Hình 3.30 Mơ hình hình học toán 50 Hình 3.31 Mơ hình hình học phần tử 3D 51 Hình 3.32 Mơ hình chia lƣới mơ hình 51 Hình 3.33 Điều kiện biên cho ghế sofa, bàn, tủ tivi giƣờng 52 Hình 3.34 Điều kiện biên cho vách – trần phịng 52 Hình 3.35 Điều kiện biên cửa phòng ngủ 53 Hình 3.36 Điều kiện biên nguồn âm 53 Hình 3.37 Sự phân bố mức cƣờng độ âm khơng gian 02 phịng 54 Hình 3.38 Dịng lƣơng âm phòng 55 CHƯ NG TỔNG QUAN VỀ ĐỀ TÀI 1.1 Tính cần thiết đề tài Mơi trƣờng âm dễ chịu rõ ràng có đƣợc biết cách kết hợp hài hòa cấu trúc âm hình học vật liệu nội thất sử dụng Ngoài vật liệu nội thất hút âm lý tƣởng cách bố trí, lắp đặt hút âm nhƣ quan trọng Âm rõ liên quan đến phù hợp âm với hoạt động diễn không gian định, cách âm tƣơng tác cách ngƣời tiếp nhận tƣợng âm Từ đó, cần phải lƣu ý đến yếu tố then chốt nhƣ truyền âm, hút âm, phản âm khuếch tán âm Điều phụ thuộc nhiều vào hình dạng, kích thƣớc, bề mặt nội thất đồ dùng trang trí gian phịng Chính vậy, thiết kế mơi trƣờng âm học phòng, cần cân yếu tố để hỗ trợ tối đa cho mục đích sử dụng Mỗi yếu tố âm tƣơng ứng với phân bổ âm khác liên quan đến tƣợng dội âm, độ rõ lời nói, khả quan thính giác độ suy yếu không gian Tùy chức khu vực, có yêu cầu khác môi trƣờng âm cần đƣợc cân nhắc trƣớc kỹ trƣớc vào thiết kế, đáp ứng tiêu chuẩn Việt Nam Do đó, thiết kế mơi trƣờng âm học phịng q trình phức tạp, dễ xảy sai sót Để giảm bớt sai sót nhầm lẫn q trình thiết kế, ta cần có giải pháp cụ thể rõ ràng để thử nghiệm thiết kế trƣớc thi công Một giải pháp thử nghiệm hiệu áp dụng chƣơng trình tính tốn số, vừa giảm đƣợc thời gian thử nghiệm, vừa giảm đƣợc chi phí thử nghiệm Tại Việt Nam, vấn đề tính toán tối ƣu cho toán truyền âm chƣa đƣợc quan tâm nhiều, có nghiên cứu liên quan đến vấn đề - Trong tiêu chuẩn Việt Nam liên quan đến lĩnh vực truyền âm khơng gian kín có tiêu chuẩn TCVN 9369: 2012 – Tiêu chuẩn thiết kế nhà hát [5] Tuy nhiên tiêu chuẩn tập trung vào nội dung thiết kế kết cấu bố cục sân khấu, chƣa có quy định cụ thể lan truyền âm qua kết cấu, khả hút âm vật liệu, nơi bố trí vật liệu phản âm tán âm … - Nguyễn Minh Nhật [4] (2011) có thực luận văn đại học đề tài tính tốn truyền âm khơng gian kín với hƣớng dẫn PGS.TS Trƣơng Tích Thiện 2 Luận văn sử dụng chƣơng trình tốn số để tính tốn q trình lan truyền âm khán phòng nhỏ Tuy chƣa có so sánh với thực nghiệm nhƣng kết luận văn cho thấy tiềm ứng dụng chƣơng trình tính tốn số vào tốn thiết kế âm kiến trúc - Nguyễn Thanh Nhã, Trƣơng Tích Thiện [6] công bố kết nghiên cứu ứng dụng phƣơng pháp khơng lƣới giải tốn truyền âm chiều Nghiên cứu đƣợc thực dựa phƣơng pháp số khơng lƣới, kết tính tốn đƣợc so sánh với thực nghiệm kết tính tốn từ ANSYS (chƣơng trình tính tốn đƣợc xây dựng dựa phƣơng pháp số: phƣơng pháp phần tử hữu hạn, phƣơng pháp thể tích hữu hạn,…) Các kết nghiên cứu chứng minh hiệu phƣơng pháp số toán truyền âm Đối với toán truyền âm, theo tác giả giới có số phƣơng pháp số đƣợc phát triển thƣơng mại hóa thành chƣơng trình mơ Đó phƣơng pháp phần tử hữu hạn (chƣơng trình ANSYS) phƣơng pháp phần tử biên (chƣơng trình SYSNOISE), phƣơng pháp Hybrid reflection method (chƣơng trình ODEON ROOM ACOUSTICS) Hiện nay, toán truyền âm chƣa đƣợc quan tâm nhiều Việt Nam Do việc nghiên cứu áp dụng chƣơng trình tính tốn số để giải toán liên quan đến truyền âm cần đƣợc quan tâm nghiên cứu Đó lý tác giả chọn thực đề tài “Ứng dụng chương trình phân tích số thiết kế âm cho hội trường” với hƣớng dẫn PGS TS Trƣơng Tích Thiện 1.2 Mục tiêu nghiên cứu 1.2.1 Mục tiêu chung Tính tốn mơ lan truyền âm mơi trƣờng kín mà cụ thể khơng gian phịng kín, hội trƣờng 1.2.2 Mục ti u cụ thể Mục tiêu 1: Tìm hiểu phân tích lý thuyết tính tốn tốn truyền âm mơi trƣờng kín, từ chọn lựa phƣơng pháp phù hợp để giải toán luận văn 3 Mục tiêu 2: Mô lan truyền âm hội trƣờng kín với cấu trúc vị trí nguồn âm đƣợc bố trí tùy ý Mục tiêu 3: Từ kết tính tốn mơ phỏng, luận văn tập trung đánh giá, phân tích để đƣa quy trình phân tích q trình lan truyền âm chƣơng trình tính tốn số 1.3 Đối tượng nghiên cứu Âm vừa mang đến thoải mái vừa gây căng thẳng cho ngƣời Do đó, kỹ thuật có nhiều dạng tốn truyền âm cần đƣợc tính tốn, mơ phỏng, kể nhƣ số tốn sau: - Bài tốn cách âm cho cơng trình kiến trúc - Bài tốn phân tích lan truyền âm cho cơng trình nhà hát, rạp hát - Bài tốn khử rung động, làm giảm tiếng ồn động xe ô tô tác động lên hành khách - Bài toán giảm tiếng ồn phát từ nhà xƣởng tác động lên công nhân dân cƣ xung quanh… Trong luận văn tác giả chọn đối tƣợng nghiên cứu đề tài tốn phân tích q trình truyền âm hội trƣờng kín 1.4 Phạm vi nghiên cứu Trong phạm vi giới hạn đề tài, tác giả tập trung ứng dụng chƣơng trình phân tích số phù hợp để thực mơ q trình truyền âm khơng gian hội trƣờng kín với điều kiện khác vị trí đặt nguồn âm, số lƣợng nguồn âm, vị trí bố trí vật liệu cách âm 4 Hình 1.1 Nhà hát Opera Paris – Garnier 1.5 Câu hỏi nghiên cứu - Phƣơng pháp tính tốn phù hợp với tốn truyền âm mơi trƣờng kín - Cơ sở đánh giá tính xác kết mô - Các tiêu chuẩn thiết kế âm trƣờng phù hợp với Việt Nam 1.6 Phương pháp nghi n cứu Nghiên cứu phƣơng pháp nhƣ chƣơng trình tính tốn cho tốn truyền âm Chọn lựa chƣơng trình phù hợp để tính tốn cho mơ hình luận văn Tìm kiếm tài liệu, báo, mơ hình thí nghiệm có liên quan đến lĩnh vực nghiên cứu luận văn để làm sở so sánh đánh giá độ tin cậy phƣơng pháp nghiên cứu Sử dụng chƣơng trình tính tốn đƣợc lựa chọn để mơ tốn có kết đƣợc cơng bố từ nhà nghiên cứu trƣớc để kiểm chứng tính xác chƣơng trình Sử dụng chƣơng trình đƣợc kiểm chứng để tính tốn cho mơ hình luận văn Phân tích kết tính tốn để đƣa kết luận kiến nghị luận văn 5 CHƯ NG C SỞ LÝ THUYẾT 2.1 Lý thuyết âm 2.1.1 Các phương trình chủ đạo Lý thuyết âm lĩnh vực khoa học liên quan đến mơ tả sóng âm Nó xuất phát từ động lực học lƣu chất Sự lan truyền sóng âm lƣu chất đƣợc mơ tả bởi: - Phƣơng trình liên tục     (  v )  Q t (2.1) Trong đó: v vector vận tốc;  khối lƣợng riêng; Q hệ số tổ hợp khối lƣợng t thời gian - Phƣơng trình bảo tồn động lƣợng  dv dt     p   S   b (2.2)  Trong đó: - Tensor ứng suất nhớt; a   p a  0c t  Q (2.3) 0 Phƣơng trình Navier-Stokes dạng tuyến tính:  v a t   0 pa    30   c  p a 4 Trong đó: v a vận tốc âm - áp lực; b Lực thân p Phƣơng trình liên tục tuyến tính  v - S  Q   (2.4)   t pa áp suất âm Phƣơng trình sóng âm đƣợc viết dƣới dạng:    p         0c  pa t  4           c p      t  t   4  Q   Q                  0  (2.5) Trong đó: c vận tốc truyền âm;  khối lƣợng riêng lƣu chất môi trƣờng truyền âm; module khối lƣu chất; K  hệ số nhớt động học Phƣơng trình (2.5) cịn đƣợc gọi phƣơng trình tiêu hao sóng truyền âm Trong áp suất hàm điều hòa đƣợc định nghĩa bởi: p ( r , t )  R e  p ( r )e j t    (2.6) Phƣơng trình (2.5) đƣợc viết gọn lại thành phƣơng trình không đồng Helmholtz:    p     4     p  j      30   c   0c   2.1.2 Cơ sở phương trình s ng   4  Q   Q                 0  p    j      (2.7) ng phần t h u hạn Phƣơng pháp số Galerkin sử dụng để xây dựng lại cho phƣơng trình (2.5), nhân (2.5) với hàm dạng w tích phân miền thể tích với số thao tác biến đổi ta đƣợc:   F   F  c w  p t  4     2c   w F 1 Q  t dv  2  4  w     2c F     nˆ  p d S  t   dV   w F  F  4 30  4 Q  0  30  p  d V  t     w p    F  nˆ  Q d S  d V  (2.8)  d V   Trong đó: - dV : Vi phân thể tích miền âm - dS : Vi phân bề mặt miền điều kiện biên - nˆ F F : Vector pháp tuyến hƣớng ngồi điều kiện biên F Từ phƣơng trình cân động lƣợng, vận tốc pháp tuyến điều kiện biên miền âm thanh, đƣợc cho bởi: v n ,F  nˆ t  4      t   0c   4 nˆ  Q  nˆ  p  t    v (2.9) Phƣơng trình (2.9) vào phƣơng trình (2.8) tạo dạng “yếu phƣơng trình (2.5), đƣợc cho bởi:   F c v   w  p w t n ,F t F dV  dS   4  w     2c F   w F Q  t dV    F p  d V  t     w p    F   4 w Q  0  30   d V  (2.10)  d V   Gia tốc pháp tuyến vi phân mơi trƣờng lƣu chất tính vi phân trực giao chuyển vị dòng lƣu chất, đƣợc cho bởi: v n ,F  nˆ t  uF t (2.11) Trong đó, u F chuyển vị phần tử lƣu chất Sau kết hợp phƣơng trình (2.11) phƣơng trình (2.10) cho ra:   F   w F c  p w t  u n ,F t 2 dV  dS   4  w     2c F   w F Q  t dV    F p  d V  t     w p   F  0  4 w Q  0  30   d V  (2.12)  d V   2.1.3 Phương trình chủ đạo với ảnh hưởng ưu ượng trung ình Phƣơng trình động lƣợng (Navier-Stokes) phƣơng trình liên tục đƣợc đơn giản hóa để đến phƣơng trình sóng đối lƣu sử dụng giả định sau: - Lƣu chất nén đƣợc - Lƣu chất khơng xốy - Nhiễu loạn áp suất nhỏ lƣu chất - Khơng có ứng suất nhớt - Khí lý tƣởng, trình đoạn nhiệt (adiabatic) thuận nghịch (reversible) - Ảnh hƣởng lƣu lƣợng trung bình Khi phƣơng trình liên tục tuyến tính (2.3) là:   p a  v  p a     v  c  t  a Q (2.13) Phƣơng trình Navier - Stokes tuyến tính là:  v a  v  v   t  a   p a  b      (2.14) Trong đó: v a vận tốc trung bình lƣu chất b lực khối, hai đại lƣợng đƣợc xác định thông qua vận tốc lƣu lực thế: v a   0 (2.15) b   Với: (2.16)  ( r ,t ) : Vận  ( r ,t ) : Lực tốc lƣu Phƣơng trình sóng đối lƣu thứ hai đƣợc cho bởi:   D 0  0c t d t  v0         0c  D 0  q  d t   (2.17) Với: D0 - dt   q Q  -  v  t c  j     v  (   )  0   p suất âm đƣợc xác định theo phƣơng trình (2.18): p a   ( j    v   ) (2.18) 2.1.4 Cơ sở phần t h u hạn cho phương trình s ng đối ưu Sử dụng phép biến đổi Galerkin cho phƣơng trình (2.17): nhân với hàm dạng tích phân miền thể tích với số thao tác biến đổi ta đƣợc: w  v0  w v 0         c F     0c F   F   w t dV    F w c  v0 w    v 0     c   d V    v 0  t   w F v  0c t dV dV   nˆ d S    (2.19)   F w qdV Hàm điều hòa vận tốc tiềm đƣợc cho bởi:  ( r ,t )  R e [  ( r )e Với  j t (2.20) ] : Biên độ vận tốc lƣu Lƣu lƣợng trung bình tĩnh thu đƣợc với giả định dịng chảy khơng xốy khơng nén đƣợc Phƣơng trình Laplace, liên quan đến vận tốc đƣợc cho bởi:   0 (2.21) Giá trị lƣu lƣợng trung bình bề măt miền đƣợc sử dụng điều kiên Neumann khác khơng cho phƣơng trình Laplace Điều kiện Neumann khơng áp dụng cho phần lại miền bề mặt mà lƣu lƣợng trung bình tiếp tuyến với bề mặt đƣợc áp dụng cho điều kiệm Neumann túy  0 Với giải định dòng chảy với lƣu lƣợng trung bình tiếp tuyến với tƣờng phƣơng trình sóng đối lƣu đƣợc giải miền tần số, điều kiện biên điều kiện chuyển động đƣợc cho bởi: nˆ v a  j  u n  v  u n  u n nˆ ( nˆ  )v (2.22) Với: nˆ : vector đơn vị điện biên un : điều kiện biên chuyển vị Miền tích phân miền chuyển động (2.19) đƣợc truyền bởi:  M 0 0 w   nˆ d S   M 0 0 u n ( j  w  v  w )d S   cM 0 0 w u n nˆ v d l (2.23) 10 Nếu vận tốc đơn vị v n điều kiện biên đƣợc biết, là: j u n  v n (2.24) Miền tích phân (2.19) đƣợc biến đổi:  M 0 0 w   nˆ d S   M 0 0 v nw dS  j 0 v  n 0  M  ( v 0 w ) d S  j C M 0 0 w v n  nˆ v d l (2.25) Nếu điều kiện trở kháng đƣợc tồn phần vùng chuyển động, đƣợc viết: j u n  Y p a (2.26) Miền tích phân (2.19) đƣợc viết lại:  M 0 0  j w   nˆ d S   Y w  d S Z  j   Y (v  w ) ( w   )d S  Z   Y  w ( v   )  ( v  w )  d s    Z  w Y (v   ) nˆ v d l (2.27) cM  C Y w  nˆ v d l Z 2.1.5 Ma trận truyền âm Phƣơng trình (2.12) bao gồm áp suất p lƣu chất thành phần chuyển vị kết cấu u x ,F ,u y ,F ,u z ,F biến phụ thuộc để đến giải toán Phần tử hữu hạn xấp xỉ hàm dạng cho thay đổi vùng chứa thành phần áp suất chuyển vị đƣợc viết: T P  { N } { Pe } T u  {N '} {u e } Với - {N } : phần tử hàm dạng áp suất - {N '} : phần tử hàm dạng chuyển vị - { Pe } : node vector áp suất (2.28) (2.29) 11 { u e } = { u xe } , { u ye } , { u ze } : - thành phần vector node chuyển vị Từ phƣơng trình (2.28) (2.29):  P t  T  { N } { Pe } (2.30) T {u }  {N } {u e } t (2.31) T  P  { N } {  Pe } (2.32) Thay phƣơng trình (2.28) (2.29), phƣơng trình sóng (2.5) trở thành:   F + 4 T c { N } { N } d V { p e } +  F T  { N  0c T } { n } { N ' } d S { u e ,F } =  F T F 0 [ N ] [ N ]d V { p e }    F T { N } { N } d V { q } +  F 4 30 T { N }{ N } d V { p e } (2.33) T [  N ] [  N ] d V {q } Với {n } : Vector pháp tuyến hƣớng lƣu chất {q } : nodal mass source vector {q } : đạo hàm bậc nodal mass source vector Các điều kiện khác đƣợc định nghĩa sở truyền âm Phƣơng trình (2.33) đƣợc viết gọn lại dƣới dạng ma trận để diễn dãi phƣơng trình sóng: M  F   p   C   p    K e e   F   M  F   F   p    R  e    T u    f  e ,F F (2.34) Với =   F C  F   =   F K  F   =   F  0c T 4  0c 0 {N }{N } d V T Ma trận khối lƣợng lƣu chất âm [ N ] [ N ]d V T [ N ] [ N ]d V Ma trận giảm chấn lƣu chất âm : Matrix độ cứng lƣu chất âm ... từ cơng trình khác nhƣ đƣợc ghi luận văn, xin cam đoan luận văn tơi thực luận văn đƣợc nộp Trƣờng Đại học Kinh tế Công nghiệp Long An Số liệu kết nghiên cứu luận văn hoàn toàn trung thực chƣa đƣợc... TRƯỜNG DIGITAL ANALYSIS PROGRAM APPLICATION IN THE SOUND DESIGN FOR THE HALL LUẬN VĂN THẠC SĨ Chuyên ngành Kỹ Thuật Xây Dựng Mã ngành: 8.58.02.01 Người hướng dẫn khoa học: PGS TS TRƯ NG T CH THIỆN... 54 3.4 Kết luận 55 CHƢƠNG KẾT LUẬN – KIẾN NGHỊ .57 4.1 Kết luận 57 4.1.1 Những ƣu điểm luận văn 57 4.1.2 Những thiếu sót luận văn