Mục đích nghiên cứu của Luận văn nhằm xây dựng và kiểm tra mô hình tính toán tấm bằng ngôn ngữ lập trình matlab. Dựa vào các kết quả phân tích đưa ra một số kết luận quan trọng và hướng phát triển tiếp tục cho Luận văn. Mời các bạn cùng tham khảo!
BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC TP HỒ CHÍ MINH TRẦN HUỲNH BẢO TRUNG PHÂN TÍCH ỨNG XỬ CỦA KẾT CẤU TẤM CHỊU TẢI TRỌNG NỔ DÙNG PHẦN TỬ TỨ GIÁC TRƠN MISQ20 TÓM TẮT LUẬN VĂN THẠC SĨ KỸ THUẬT XÂY DỰNG Thành Phố Hồ Chí Minh - 2018 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC TP HỒ CHÍ MINH TRẦN HUỲNH BẢO TRUNG PHÂN TÍCH ỨNG XỬ CỦA KẾT CẤU TẤM CHỊU TẢI TRỌNG NỔ DÙNG PHẦN TỬ TỨ GIÁC TRƠN MISQ20 Chuyên ngành: Kỹ thuật xây dựng Mã số: 8.58.02.01 TÓM TẮT LUẬN VĂN THẠC SĨ KỸ THUẬT XÂY DỰNG NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS NGUYỄN VĂN HIẾU Thành phố Hồ Chí Minh - 2018 i DANH MỤC CHƯƠNG 1.TỔNG QUAN 1.1.GIỚI THIỆU CHUNG: 1.2.TÌNH HÌNH NGHIÊN CỨU ĐÃ ĐƯỢC CƠNG BỐ VÀ SỰ ĐĨNG GĨP CỦA ĐỀ TÀI 1.2.1.Tình hình nghiên cứu ngồi nước: 1.2.2.Các cơng trình nghiên cứu nước: 1.2.3.Sự đóng góp đề tài: 1.2.4.Mục tiêu phạm vi nghiên cứu đề tài: CHƯƠNG 2.CƠ SỞ LÝ THUYẾT 2.1.NỔ TRONG MÔI TRƯỜNG KHÔNG KHÍ 2.2.LÝ THUYẾT TẢI TRỌNG NỔ 2.3.MÔ PHỎNG TẢI TRỌNG NỔ TÁC DỤNG LÊN KẾT CẤU TẤM BÊ TÔNG CỐT THÉP 2.3.1.Các giả thuyết mô tải trọng nổ - Tác nhân gây nổ “đủ xa” cơng trình để đảm bảo nghiệm vật lý toán 2.3.2.Mô tải trọng nổ tác dụng lên bê tông cốt thép 2.4.CÁC GIẢ THUYẾT PHÂN TÍCH ĐỘNG LỰC HỌC TẤM BÊ TƠNG CỐT THÉP CHỊU UỐN 2.4.1.Lý thuyết dày 2.4.1.Phương trình chuyển động: 2.5.PHƯƠNG PHÁP PHẦN TỬ TỨ GIÁC TRƠN MISQ20 2.5.1.Biến dạng trơn uốn 2.5.2.Biến dạng trượt ii 2.5.3.Ma trận cản: 2.6.PHƯƠNG PHÁP TÍCH PHÂN NEWMARK 10 2.6.1.Phương pháp giải: 10 2.6.2.Bài tốn dộng lực học tuyến tính 10 CHƯƠNG 3.KẾT QUẢ PHÂN TÍCH SỐ 11 3.1.Bài toán 1: Khảo sát hội tụ tốn việc mơ hình hố tải trọng nổ dạng tải trọng nổ 11 3.2.Bài toán 2: Khảo sát thời gian phân tích ứng xử bê tông cốt thép chịu tải trọng nổ phần tử tứ giác trơn MISQ20 so với phần thử hữu hạn Q4 thông thường 12 3.3.Bài toán 3: Khảo sát theo kết thực nghiệm từ “bài báo cáo phân tích động lực học tải trọng nổ” A C Jacinto, R D Ambrosini, R F Danesi (2002) [38] 13 3.4.Bài toán 4: Khảo sát ảnh hưởng kích thước đến ứng xử bê tơng cốt thép 15 3.5.Bài tốn 5: Phân tích ảnh hưởng vật liệu đến ứng xử bê tông cốt thép 16 3.6.Bài tập 6: Khảo sát ảnh hưởng khoảng cách nổ tới ứng xử bê tông cốt thép 17 3.7.Bài tập 7: Khảo sát ảnh hưởng khối lượng thuốc nổ tới ứng xử bê tông cốt thép 18 CHƯƠNG 4.KẾT LUẬN VÀ KIẾN NGHỊ 19 4.1 Kết luận: 19 4.2 Kiến nghị 19 iii DANH MỤC HÌNH ẢNH, ĐỒ THỊ Hình 2.8 Lịch sử áp suất – thời gian sóng nổ Hình 2.11 Mơ hình phân tích tốn bê tơng cốt thép Hình 2.12 Mơ hình tính r(j,i) Hình 3.1 Chuyển vị nút trọng tâm bước thời gian t=0.04s chia nhỏ phần tử 11 Hình 3.2 Biểu đồ chuyển so sánh chuyển vị nút theo thời gian phần tử MISQ20 Q4 12 Hình 3.3 Biểu đồ so sánh thời gian phân tích số phần tử tứ giác trơn MISQ20 phần tử thông thường Q4 13 Hình 3.4 Kết đo áp lực thực nghiệm theo A C Jacinto, R D Ambrosini, R F Danesi (2002) [36] 13 Hình 3.6 Biểu đồ áp lực sóng nổ phần tử tâm thép 14 Hình 3.7 Biểu đồ so sánh gia tốc chuyển động thép sử dụng phương pháp số phần tử MISQ20, phần mềm ABAQUS phương pháp thực nghiệm 14 Hình 3.8 Chuyển vị lớn vị trí tâm vị trí tọa độ tâm nổ thay đổi bề dày 15 Hình 3.9 Chuyển vị theo thời gian vị trí tâm vị trí tọa độ tâm nổ ứng với cấp độ bền bê tông khác 16 Hình 3.10 Chuyển vị theo thời gian vị trí tâm vị trí tọa độ tâm nổ thay đổi khoảng cách nổ R0 17 Hình 3.11 Chuyển vị theo thời gian vị trí tâm vị trí tọa độ tâm nổ thay đổi khối lượng thuốc nổ 18 iv DANH MỤC BẢNG Bảng 3.3.So sánh chuyển vị bước thời gian lặp t = 0.04s tăng số lượng phần tử bê tông cốt thép 11 Bảng 3.4.Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ thay đổi chiều dày bê tông cốt thép 15 Bảng 3.5 Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ ứng với giá trị cấp độ bền khác 16 Bảng 3.6.Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ ứng với giá trị khoảng cách thuốc nổ R0 khác 17 Bảng 3.7 Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ ứng với giá trị khối lượng thuốc nổ thay đổi 18 v TÓM TẮT Trong năm gần đây, vấn đề ngăn ngừa phòng chống thiệt hại vũ khí đại gây quan tâm nghiên cứu nhà khoa học nước Những vụ nổ thường gây thiệt hại nghiêm trọng cho cơng trình, ảnh hưởng tính mạng tài sản người Việc nghiên cứu ảnh hưởng kết cấu chịu tải trọng nổ cần thiết quan trọng đảm bảo an tồ cho cơng trình, tính mạng tài sản người Trong luận văn này, việc phân tích ứng xử động lực học kết cấu bê tông cốt thép chịu tải trọng nổ nghiên cứu Sử dụng phương pháp số để mơ hình hố tải trọng nổ dạng tải trọng động, mơ hình hố bê tơng cốt thép ngàm cạnh phần tử tứ giác trơn với bốn nút bậc tự MISQ20 thiết lập phương trình chuyển động giải phương pháp Newmark Phân tích ứng xử tuyến tính theo bước thời gian lần lược khảo sát thay đổi bề dày bê tông cốt thép, khoảng cách thuốc, cao trình đặt chất nổ khối lượng thuốc nồ TNT tương đương So sánh đồ thị chuyển vị lớn điểm khảo sát bê tông cốt thép Xây dựng kiểm tra mơ hình tính tốn ngơn ngữ lập trình matlab Dựa vào kết phân tích đưa số kết luận quan trọng hướng phát triển tiếp tục cho đề tài CHƯƠNG TỔNG QUAN 1.1 GIỚI THIỆU CHUNG: Trong bối cảnh giới hoạt động khủng bố mối đe dọa xung đột tranh chấp chủ quyền lãnh thổ quốc gia, xung đột sắc tộc trở thành vấn nạn gia tăng khắp giới Để phòng tránh loại bom đạn, nhiều quốc gia nghiên cứu, xây dựng cơng trình quốc phịng, mục đích che giấu, bảo vệ lực lượng quân vũ khí trước phá hoại bom đạn Trong luận văn này, nghiên cứu ứng xử kết cấu bê tông cốt thép chịu tải trọng nổ dùng phần tử tứ giác trơn bốn nút với bậc tự MISQ20 Nguyễn Văn Hiếu, từ tốn phân tích tuyến tính, phương pháp đề xuất đánh giá tính hiệu quả, độ xác độ tin cậy Luận văn nghiên cứu ảnh hưởng thuốc nổ TNT để đảm bảo toán cho kết tin cậy, với khoảng cách hợp lý để tránh phá hoại đột ngột bê tơng cốt thép 1.2 TÌNH HÌNH NGHIÊN CỨU ĐÃ ĐƯỢC CƠNG BỐ VÀ SỰ ĐĨNG GĨP CỦA ĐỀ TÀI 1.2.1 Tình hình nghiên cứu ngồi nước: Đã có nhiều cơng trình nghiên cứu mơ hình hóa áp lực nổ tác dụng lên vật thể kết cấu thực hiện: - Smith (1994) [3], Mơ hình hoá áp suất phản xạ đỉnh tương tác sóng nổ với bề mặt mục tiêu tác động tĩnh - Các tác giả nghiên cứu khác: M.Remennikov (2003) [4] mô 3D AUTODYN [5], Nelson Lam đồng nghiệp (2004) [6] M.Remennikov Rose (2005) [7], Ngo tác giả (2007) [8], Tavakoli Kiakojouri (2012) [9] - Các vấn đề dùng mơ số để khảo sát phi tuyến hình học cho vỏ phương pháp phần tử hữu hạn trơn nghiên cứu: Giáo sư G R Liu Nguyễn Thời Trung (2010) [13] methodSFEM) Nguyễn Văn Hiếu (2009) [14], Nguyễn Xuân Hùng cộng (2008) [15, 16] , Nguyễn Văn Hiếu (2009) [14] nghiên cứu phi tuyến hình học dựa phần tử MISQ20, X.Y Cui cộng (2008) [17] 1.2.2 Các cơng trình nghiên cứu nước: - Ở việt nam, Đỗ Kiến Quốc Lê Đức Tuấn (2007) [18] nghiên cứu ứng xử động lực học kết cấu bê tông cốt thép tác dụng tải trọng nổ Nguyễn Đức Duyến (2010) [19] nghiên cứu phương pháp tính tốn cơng trình phịng thủ dân chịu tác động tải trọng nổ - Nguyễn Trọng Phước Trần Minh Thi (2011) [20] phân tích động lực học kết cấu tường dọc tác động tải trọng nổ - Đỗ Ngọc Thuận (2014) [21], phân tích ứng xử dày đàn hồi chịu tải trọng nổ Các nghiên cứu phân tích phi tuyến hình học kết cấu dạng tấm/vỏ phần tử hữu hạn trơn phát triển như: Nguyễn Hoài Nam, Nguyễn Văn Hiếu, Lương Văn Hải, Châu Đình Thành (2013) [22], Phạm Hồng Cơng, Nguyễn Đình Đức (2013) [23], Nguyễn Văn Hiếu, Nguyễn Hoài Nam, Trần Đồng Kiếm Lam, Lê Văn Thơng (2013) [24] 1.2.3 Sự đóng góp đề tài: Tính mới: Lần áp dụng phần tử tứ giác trơn MISQ20 cho phân tích kết cấu chịu tải trọng nổ Tính thời sự: Phương pháp phần tử hữu hạn trơn MISQ20 tác giả Nguyễn Văn Hiếu [14] phát triển sử dụng rộng rãi mơ phỏng, phân tích ứng xử với độ phức tạp khác Góp phần nâng cao hiểu biết ứng xử kết cấu xảy tượng nổ Ý nghĩa khoa học: Mở rộng phạm vi ứng dụng phần tử hữu hạn trơn MISQ20 cho tính tốn kết cấu chịu tải trọng biến đổi theo thời gian 1.2.4 Mục tiêu phạm vi nghiên cứu đề tài: Mô tải trọng nổ dạng tải trọng động dựa số giả thiết để đảm bảo nghiệm vật lí tốn động lực học Xây dựng ma trận tính kết cấu: ma trận độ cứng, ma trận khối lượng, ma trận cản véc tơ tải trọng nút, không xét phi tuyến Thiết lập phương trình chuyển động giải phương pháp tích phân trực tiếp Newmark, sử dụng trình tự phân tích bước thời gian để xác định ứng xử bê tông cốt thép Mô tải trọng nổ tác dụng lên mơ hình kiểm tra so sánh tải trọng nổ, ứng xử với nghiên cứu công bố Khảo sát hội tụ, sai số lưới chia phần tử Phân tích ứng xử bê tơng cốt thép chịu tải trọng nổ phần tử MISQ20 so với phần tử hữu hạn Q4 thông thường Khảo sát yếu tố ảnh hưởng tới ứng xử bê tơng cốt thép: kích thước, vật liệu khoảng cách nổ, khối lượng thuốc nổ 12 Nhận xét: Qua kết bảng 3.3 hình 3.1 ta thấy việc phân tích ứng xử tác dụng tải trọng nổ phần tử tứ giác trơn MISQ20 ngơn ngữ lập trình matlab hội tụ lưới chia phần tử 1414 Điều chứng tỏ chia nhỏ với lưới phần tử 1414 đủ nghiệm toán bắt đầu hội tụ 3.2 Bài toán 2: Khảo sát thời gian phân tích ứng xử bê tơng cốt thép chịu tải trọng nổ phần tử tứ giác trơn MISQ20 so với phần thử hữu hạn Q4 thông thường Kết phân tích chuyển vị nút bê tông cốt thép lưới phần tử 16 16 hình 3.2 Hình 3.2 Biểu đồ chuyển so sánh chuyển vị nút theo thời gian phần tử MISQ20 Q4 Nhận xét: Qua hình 3.2 đồ thị chuyển vị nút bê tông cốt thép tải trọng nổ theo thời gian cách sử dụng phần tử MISQ20 với phần tử hữu hạn Q4 thông thường, cho ta thấy đường đồ thị gần trùng nhau, nên sai số kết nhỏ độ tin cậy hai phần tử Kết so sánh thời gian xử lý phân tích ứng xử bê tông cốt thép tải trọng nổ phần tử tứ giác trơn MISQ20 phần tử hữu hạn thông thường Q4 hình 3.3 13 Hình 3.3 Biểu đồ so sánh thời gian phân tích số phần tử tứ giác trơn MISQ20 phần tử thông thường Q4 Nhận xét: Từ kết biểu đồ thời gian phân tích liệu cho ta thấy phần tử tứ giác trơn MISQ20 có tốc độ xử lý nhanh phần tử hữu hạn Q4 thông thường 3.3 Bài toán 3: Khảo sát theo kết thực nghiệm từ “bài báo cáo phân tích động lực học tải trọng nổ” A C Jacinto, R D Ambrosini, R F Danesi (2002) [38] Hình Kết đo áp lực thực nghiệm theo A C Jacinto, R D 14 Ambrosini, R F Danesi (2002) [38] Kết phân tích ứng xử thép tải trọng nổ phần tử hữu hạng MISQ20 xét nút tâm hình 3.6 Hình Biểu đồ áp lực sóng nổ phần tử tâm thép Hình Biểu đồ so sánh gia tốc chuyển động thép sử dụng phương pháp số phần tử MISQ20, phần mềm ABAQUS phương pháp thực nghiệm Nhận xét: Qua kết biểu đồ áp lực nổ hình 3.6 dùng phần tử tứ giác 15 trơn MISQ20 giống với biểu đồ áp nổ hình 3.5 trích từ báo nghiên cứu thực nghiệm A C Jacinto, R D Ambrosini, R F Danesi (2002) [38] hình 3.7 biểu đồ so sánh gia tốc chuyển động thép phần tử MISQ20 [14], phần mềm ABAQUS thực nghiệm [38] gần trùng giai đoạn dao động cưỡng áp lực nổ Các kết khảo sát có sai số khơng đáng kể 3.4 Bài tốn 4: Khảo sát ảnh hưởng kích thước đến ứng xử bê tông cốt thép Bảng 3.4.Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ thay đổi chiều dày bê tơng cốt thép Hình 3.8 Chuyển vị lớn vị trí tâm vị trí tọa độ tâm nổ thay đổi bề dày 16 Nhận xét: Qua kết khảo sát bảng 3.4 hình 3.8 ta tăng chiều dày tâm vị trí toạ độ thuốc nổ có chuyển vị giảm dần Với thông số trên, để đảm bảo bê tông phù hợp quy định độ võng giới hạn đàn hồi, không bị phá hoại nên xây bê tơng có bề dày từ 0.56m trở lên 3.5 Bài tốn 5: Phân tích ảnh hưởng vật liệu đến ứng xử bê tông cốt thép Bảng 3.5 Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ ứng với giá trị cấp độ bền khác Hình 3.9 Chuyển vị theo thời gian vị trí tâm vị trí tọa độ tâm nổ ứng với cấp độ bền bê tông khác 17 Nhận xét: Qua kết khảo sát bảng 3.5 hình 3.9 cho ta thấy chuyển vị bê tơng cốt thép khơng có chênh lệch đáng kể thay đổi cấp độ bền bê tông Đường biểu diễn chuyển vị có độ dốc không lớn ta tăng cấp độ bền bê tông Để giảm chuyển vị khơng nên dùng biện pháp tăng cấp độ bền bê tông 3.6 Bài tập 6: Khảo sát ảnh hưởng khoảng cách nổ tới ứng xử bê tông cốt thép Bảng 3.6 Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ ứng với giá trị khoảng cách thuốc nổ R0 khác Hình 3.10 Chuyển vị theo thời gian vị trí tâm vị trí tọa độ tâm nổ thay đổi khoảng cách nổ R0 18 Nhận xét: Qua kết khảo sát bảng 3.4 hình 3.9 cho ta thấy khoảng cách trường hợp bê tơng cốt thép xa tác nhân gây nổ chuyển vị giảm dần Do để đảm bảo chuyển vị bê tơng cốt phải có giải pháp giảm khoảng cách nổ từ 3.0m trở lên, có phương án làm tăng bề dày bê tơng cốt thép bề dày 3.7 Bài tập 7: Khảo sát ảnh hưởng khối lượng thuốc nổ tới ứng xử bê tông cốt thép Bảng 3.7 Kết chuyển vị lớn tâm vị trí tọa độ tâm nổ ứng với giá trị khối lượng thuốc nổ thay đổi Hình 3.11 Chuyển vị theo thời gian vị trí tâm vị trí tọa độ tâm nổ thay đổi khối lượng thuốc nổ 19 Nhận xét: Qua kết khảo sát bảng 3.4 hình 3.9 cho ta thấy tăng thuốc nổ chuyển vị tăng Hình 3.10 cho ta thấy, đường biểu đồ chuyển vị độ võng tâm có độ dốc lớn thẳng, chứng tỏ khối lượng thuốc nổ tăng chuyển vị tăng tuyến tính CHƯƠNG KẾT LUẬN VÀ KIẾN NGHỊ 4.1 Kết luận: (1) Phương pháp phần tử hữu hạn tứ giác trơn MISQ20 có xác cao, khối lượng tính tốn lớn, thời gian phân tích nhanh nhiều so với phần tử hữu hạn Q4 thơng thường Bài tốn phần tử tứ diện bốn nút với bậc tự MISQ20 hình thành để mơ hình phân tích sâu kết cấu (2) Dưới tác dụng tải trọng nổ với khối lượng thuốc nổ phù hợp khoảng cách định, Bài toán kết cấu bê tơng cốt thép chịu tải trọng nổ đưa dạng tốn phân tích ứng xử động lực học kết cấu bê tông cốt thép (3) Các yếu tố kích thước chiều dày bê tơng cốt thép, tính chất vật liệu bê tơng cốt thép, khoảng cách nổ, khối lượng thuốc nổ có ảnh hưởng nhiều đến ứng xử bê tơng cốt thép Tuy nhiên, yếu tố tính chất vật liệu bê tông cụ thể module đàn hồi E có ảnh hưởng khơng rõ rệch đến chuyển vị bê tông cốt thép chịu tải trọng nổ Vì để giảm chuyển vị bê tơng cốt thép phương pháp tăng cấp độ bền bê tông không mang lại hiệu quà cao Phương án tăng chiều dày bê tông cốt thép chịu tải trọng nổ hiệu tốt 4.2 Kiến nghị 20 Những vấn đề định có khả phát triển nghiên cứu tốt tương lai sau: (1) Phương pháp phần tử hữu hạn tứ giác trơn MISQ20 có khả phân tích biên nên đảm bảo độ xác với lưới thơ phần tử tứ giác có dạng lồi lõm với toán phức tạp, điều kiện biên xác so với phần tử dựa kỹ thuật tích phân miền phần tử thơng thường Do vậy, cần phát triển nghiên cứu toán có hình dạng lưới thơ hay phần tử tứ giác có dạng lồi lõm để khảo sát tính ưu việt, hiệu độ xác cao phần tử tứ giác trơn MISQ20 (2) Khi phân tích ứng xử tải trọng nổ dùng phần tử tứ giác trơn MISQ20, để đảm bảo việc mô kết cấu tính chất vật liệu làm việc phù hợp với thực tế yếu tố phi tuyến cần nghiên cứu phát triển (3) Dưới tác dụng thuốc nổ lên kết cấu bê tông cốt thép, nhiệt lượng lớn toả tác dụng lên kết cấu bê tông cốt thép ảnh hưởng đến tính chất vật liệu làm thay đổi khả chịu lực bê tông cốt thép Do để đảm bảo tính xác cao hơn, cần nghiên cứu phân tích yếu tố nhiệt độ sau nổ Bài tốn bê tơng cốt thép chịu tải trọng nổ xét trường hợp điều kiện biên bê tông cốt thép ngàm 01 cạnh, chưa xét trường hợp bê tông cốt thép gàm 04 cạnh Trong trường hợp tác động sóng nổ lên bê tông cốt thép liên kết cứng bị cản trở dao động dẫn tới ứng xử có nhiều thay đổi khác Trang TÀI LIỆU THAM KHẢO [1] Alexander M Remennikov, “A Review of Methods for Predicting Bomb Blast Effects on Buildings”, Faculty of Engineering, University of Wollongong, NSW.2522, Australia [2] Kingery CN and Blulmash G, (April 1984), “Airblast parameters from TNT Sherical Air Blast and Hemispherical Surface Blast”,Technical Report ARBRL – TR – 02555, US Armament Reseach and Development Cennical, Ballistic Reseach Laboratory, Aberdeen proving Ground, MD [3] Smith PD and Hetherington JG, (1994), “Blast and Ballistic Loading of structures”, Butterowrth-Heinemenn, Oxford [4] Alexander M Remennikov, (2003), “A Review of Methods for Predicting Bomb Blast Effects on Buildings”, Faculty of Engineering, University of Wollongong, NSW 2522, Australia [5] Century Dynamics, (2002), “AUTODYN-2D & 3D User Documentation”, [6] Nelson Lam, Priyan Mendis and Tuan Ngo (2004), Response Spectrum Solution for Blast Loading, The University of Melbourne, Parkville 3010, Victoria, Australia [7] Remenennikow, AM and Rose, AT (2005), “Modelling blast loads on buildins in complex city geometries”, Faculty of Engineering, Trang University in complex city geometries”, Faculty of Engineering, University of Wollongong, NSW.2522, Australia [8] Ngo, T, Mendis, P, Gupta, A and Ramsay, J (2007), “Blast loading and blast effects on structure – An Overview”, EJSE Special Inssure: Loading on Structrures Ejse Internation, 16p [9] Tavakoli, HR, Kiakojouri, F (2012), “Numerical dynamic anasis of stiffened plate under blast loading’’, Latin American Journal oSolids and Structures, 15p [10] Gilmour J R (1998), “Numerical Modelling of the Progressive Collapse of Structructures as a result of Impact or Explosion”, PhD Thesis, City University, London, UK [11] S W Park, Q Xia and M Zhou, “Dynamic behavior of concrete at high strain rates and pressures”, Internationak Journal of Impact Engineering 25 (2001) 887-910 [12] Anupam Sharma, Lyle N Long and Ted Krauthammer, “Using the direct simulation Monto Carlo approach for the Blast – Impact problem”, The Pennsylvania State University, University Park, PA 16802, USA [13] GR, L and N TT (2010), “Smoothed finite element methods”, CRC Press: Taylor and Francis Group, New York [14] Nguyen-Van, H (2009), “Development and application of assumed strain smoothing finite element technique for composite Trang plate/shell structures”, Ph.D Thesis, University of southern queensland’s dissertation, Australia [15] H Nguyen-Xuan, T Rabczuk, Ste´phane Bordas, J.F Debongnie, "A smoothed finite element method for plate analysis," Computer Methods in Applied Mechanics and Engineering, 197, pp 1184–1203, 2008 [16] N.Nguyen-Thanh, Timon Rabczuk, H.Nguyen-Xuan, Ste´phane Bordas (2008), "A smoothed finite element method for shell analysis", Computer Methods in Applied Mechanics and Engineering, 198, pp 165-177 [17] X.Y Cui, G R Liu, G Y Li, X Zhao, T.T Nguyen, G.Y Sun (2008), “A Smoothed Finite Element Method (SFEM) for Linear and Geometrically Nonlinear Analysis of Plates and Shells’’, Computer Modeling in Engineering & Sciences, 28, pp 109-126 [18] Đỗ Kiến Quốc, Lê Đức Tuấn, Ứng xử động lực học kết cấu bê tông cốt thép tác dụng tải trọng nổ, T/C Khảo sát thiết kế, số 2, 2007 [19] Nguyễn Đức Duyến, (2009), “Tính cơng trình phịng thủ dân chịu tác dụng sóng nổ có kể đến cấp bảo vệ hợp lý cơng trình”, Luận văn Tiến sĩ kỹ thuật, Học viện Kĩ thuật Quân [20] Nguyễn Trọng Phước Trần Minh Thi, 2011 “Response of vertical wall structure under under blast loading by dynamic Trang analysis” ,Procedia Engineering, Vol 14, pp 3308-3316 (doi:10.1016/j.proeng.2011.07.418) [21] Đỗ Ngọc Thuận, (2014),“Phân ứng xử đàn hồi tác dụng tải trọng nổ”, Luận văn thạc sĩ, Trường đại học mở TP.HCM [22] Nguyễn Hoài Nam, Nguyễn Văn Hiếu, Lương Văn Hải, Châu Đình Thành (2013) “Phân tích ứng xử phi tuyến hình học kết cấu vỏ chịu tải trọng tĩnh phương pháp phần tử hữu hạn trơn”, Tạp chí Xây Dựng, Số 11, trang 105-108 [23] Phạm Hồng Công, Nguyễn Đình Đức, “Phân tích phi tuyến đối xứng S-FGM đàn hồi sử dụng lý thuyết biến dạng trượt bậc phương pháp đặt hàm chuyển vị”, Hội nghị Khoa học toàn quốc Cơ học Vật rắn biến dạng lần thứ XI Thành phố Hồ Chí Minh, 7-9/11/2013 272-281 [24] Nguyễn Văn Hiếu, Nguyễn Hồi Nam, Trần Đồng Kiếm Lam, Lê Văn Thơng "Mơ hình phân tích phi tuyến hình học kết cấu tấm/vỏ composite sử dụng phần tử tứ giác trơn" Hội nghị Khoa học toàn quốc Cơ học Vật rắn biến dạng lần thứ XI Thành phố Hồ Chí Minh, 7-9/11/2013 [25] Wilson Fickett (1985), “Detonation in miniature”, In John D Buckmaster, Editor, The Mathematics of Combustion, pp 133-181, SIAM, Philadenphia, 1985 Trang [26] Bernar Lewis and Guenther von Elbe (1961), “Combustion Flames and Explosis of Gasses”, Academic Press, New York, Second Edition [27] Krauthammer T., (2000).“Modern Protective Structures, Design, Analysic and Evaluation, Course notes”, The Pennsylvania State University, 358pp [28] Henrych J (1979), “The Dynamics of Explosion and Its Use”, Elsevier, Amsterdam [29] Rankine W.J.H (1870), “Philosophical Transactions of the Royal Society”, Vol 160, 277 [30] Brode HL (1955), “A Numerical Solutions of Spherical Blast Waves”, Journal of Applied Physics, American Institute of Physics No.6 [31] W E Baker, P A Cox, P S Westine, J J Kulesz and R A Stehlow, (1983), “Explosion Hazards and Evaluation”, Elsevier Scientific Publishing Company, New York [32] Bulson P.S, (1997), “Explosive loading of engineering structures”, E & FN Spon, London, 233 pp [33] Nelson Lam, Priyan Mendis and Tuan Ngo, “Response Spectrum Solution for Blast Loading”, The University of Mebourne, Parkville 3010, Victoria, Australia Trang [34] Gilbert Ford Kinney (1962), “Explosive Shocks in Air”, The MacMinllan Company, New York [35] Chu Quốc Thắng (1997), “Phương Pháp Phần Tử Hữu Hạn”, Nhà Xuất Bản Khoa Học Kỹ Thuật, [36] Đỗ Kiến Quốc (2003), “Bài Giảng Động Lực Học Kết Cấu”, Phòng Đào Tạo Sau Đại Học, Trường Đại Học Bách Khoa TP Hồ Chí Minh [37] Đỗ Kiến Quốc Nguyễn Trọng Phước (2010), “Các phương pháp số trọng động lực học kết cấu” ,TP Hồ Chí Minh: NXB ĐHQG Tp.HCM [38] A C Jacinto, R D Ambrosini, R F Danesi (2002), “Dynamic response of plates subjected to blast loading”, Proceedings of the Institution of Civil Engineers, Structures & Buildings 152, August 2002 Issue 3, Pages 269-276 ... TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC TP HỒ CHÍ MINH TRẦN HUỲNH BẢO TRUNG PHÂN TÍCH ỨNG XỬ CỦA KẾT CẤU TẤM CHỊU TẢI TRỌNG NỔ DÙNG PHẦN TỬ TỨ GIÁC TRƠN MISQ20 Chuyên ngành: Kỹ thuật xây dựng... gian phân tích số phần tử tứ giác trơn MISQ20 phần tử thông thường Q4 Nhận xét: Từ kết biểu đồ thời gian phân tích liệu cho ta thấy phần tử tứ giác trơn MISQ20 có tốc độ xử lý nhanh phần tử hữu... gần trùng nhau, nên sai số kết nhỏ độ tin cậy hai phần tử Kết so sánh thời gian xử lý phân tích ứng xử bê tông cốt thép tải trọng nổ phần tử tứ giác trơn MISQ20 phần tử hữu hạn thơng thường Q4