Phân tích động lực học tấm bê tông cốt thép chịu tác dụng tải trọng nổ

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA TRẦN MINH THI PHÂN TÍCH ĐỘNG LỰC HỌC TẤM BÊ TÔNG CỐT THÉP CHỊU TÁC DỤNG TẢI TRỌNG NỔ CHUYÊN NGÀNH: XÂY DỰNG DÂN DỤNG VÀ CÔNG NGHIỆP MÃ SỐ NGÀNH: 60.58.20 LUẬN VĂN THẠC SĨ KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH, THÁNG 12 NĂM 2009 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 : TS Nguyễn Trọng Phước Cán chấm nhận xét 1: TS Ngô Hữu Cường Cán chấm nhận xét 2: TS Nguyễn Minh Long Luận văn thạc sĩ bảo vệ HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày 26 tháng 01 năm 2010 NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN   ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  NHẬN XÉT CỦA GIÁO VIÊN PHẢN BIỆN   ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐  ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HOÀ XÃ HỘI CHỦ NGHIÃ VIỆT NAM Độc Lập - Tự Do - Hạnh Phúc -oOo Tp HCM, ngày 04 tháng 12 năm 2009 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: Trần Minh Thi Giới tính : Nam Ngày, tháng, năm sinh : 14-02-1984 Nơi sinh : Bình Định Chuyên ngành : … Xây dựng dân dụng công nghiệp… Khoá (Năm trúng tuyển) : 2007 1- TÊN ĐỀ TÀI: PHÂN TÍCH ĐỘNG LỰC HỌC TẤM BÊ TƠNG CỐT THÉP CHỊU TẢI TRỌNG NỔ 2- NHIỆM VỤ LUẬN VĂN: – Tổng quan tình hình nghiên cứu tải trọng nổ, tìm hướng phát triển cho đề tài – Mô tải trọng nổ tác dụng lên bê tơng cốt thép – Thuyết lập ma trận có tính chất kết cấu: ma trận khối lượng, ma trận độ cứng, ma trận cản có xét đến ảnh hưởng phi tuyến – Giải phương trình chuyển động phương pháp tích phân số Newmark – Phân tích thí dụ số cho kết ứng xử bê tông cốt thép chịu tải trọng nổ, kiểm chứng tính đắn chương trình lập biểu đồ so sánh kết – Xây dựng chương trình phân tích động lực học bê tơng cốt thép chịu tải trọng nổ có xét phi tuyến nhờ hổ trợ ngơn ngữ lập trình kỹ thuật Matlab – Từ kết phân tích, rút số kết luận hướng phát triển cho nghiên cứu sau 3- NGÀY GIAO NHIỆM VỤ : 4- NGÀY HOÀN THÀNH NHIỆM VỤ : 04-12-2009 5- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN (Ghi đầy đủ học hàm, học vị ): TS Nguyễn Trọng Phước Nội dung đề cương Luận văn thạc sĩ Hội Đồng Chuyên Ngành thông qua CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) TRƯỞNG BAN QUẢN LÝ CHUYÊN NGÀNH (Họ tên chữ ký) i MỤC LỤC DANH MỤC CÁC HÌNH VẼ DANH MỤC CÁC BẢNG BIỂU LỜI CẢM ƠN TÓM TẮT CHƯƠNG MỞ ĐẦU 1.1 GIỚI THIỆU 1.2 TỒNG QUAN 1.3 MỤC TIÊU VÀ PHẠM VI CỦA ĐỀ TÀI CHƯƠNG TẢI TRỌNG NỔ 2.1 GIỚI THIỆU 2.2 SĨNG NỔ TRONG KHƠNG KHÍ 2.3 THIẾT LẬP CÁC PHƯƠNG TRÌNH SĨNG NỔ 2.4 CÁC THƠNG SỐ MẶT ĐẦU SĨNG 2.5 LÝ THUYẾT TẢI TRỌNG NỔ 2.6 MƠ PHỎNG TẢI TRỌNG NỔ 2.7 KẾT LUẬN CHƯƠNG PHÂN TÍCH ĐỘNG LỰC HỌC TẤM B Ê TÔNG CỐT THÉP 3.1 GIỚI THIỆU 3.2 THIẾT LẬP TRƯỜNG CHUYỂN VỊ CỦA PHẦN TỬ TẤM CHỊU UỐN 3.3 CÁC GIẢ THIẾT PHÂN TÍCH ĐỘNG LỰC HỌC TẤM CHỊU UỐN 3.4 PHƯƠNG TRÌNH CHUYỂN ĐỘNG 3.5 XÁC ĐỊNH CÁC MA TRẬN TÍNH CHẤT KẾT CẤU 3.6 PHƯƠNG PHÁP TÍCH PHÂN NEWMARK 3.7 THUẬT TỐN PHÂN TÍCH ĐLH TẤM BTCT CHỊU TẢI TRỌNG NỔ 3.8 GIỚI THIỆU CHƯƠNG TRÌNH BLASTSHELL 3.9 HƯỚNG DẪN SỬ DỤNG CHƯƠNG TRÌNH 3.10 KẾT LUẬN CHƯƠNG CÁC THÍ DỤ SỐ 4.1 GIỚI THIỆU 4.2 KIỂM CHỨNG CHƯƠNG TRÌNH BLASTSHELL 4.3 PHÂN TÍCH ĐLH TẤM BTCT VÀ KIỂM TRA SỰ HỘI TỤ BÀI TOÁN 4.4 ẢNH HƯỞNG CỦA BỀ DÀY TẤM ĐẾN CHUYỂN VỊ NGANG 4.5 ẢNH HƯỞNG CỦA TỈ SỐ CẢN ĐẾN ỨNG XỬ CỦA TẤM 4.6 ẢNH HƯỞNG CÁC THÔNG SỐ ĐẶC TR ƯNG CỦA TẢI TRỌNG NỔ 4.7 KẾT LUẬN CHƯƠNG KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 5.1 KẾT LUẬN 5.2 HƯỚNG PHÁT TRIỂN TÀI LIỆU THAM KHẢO MỤC LỤC ii iv vi vii 1 10 10 10 12 17 21 26 33 34 34 34 36 37 38 43 46 48 48 51 52 52 52 56 68 70 72 83 84 84 85 86 ii DANH MỤC CÁC HÌNH VẼ Hình 2.1 Mơi trường nổ khơng khí theo Krauthammer (2000) 11 Hình 2.2 Mơi trường nổ bề mặt theo Krauthammer (2000) 11 Hình 2.3 Bên tác nhân theo Henrych (1979) 12 Hình 2.4 Thể tích kiểm sốt phần tử theo Henrych (1979) 13 Hình 2.5 Bên tác nhân theo Henry ch (1979) 16 Hình 2.6 Mối quan hệ áp suất động lực học v áp suất tĩnh định 18 Hình 2.7 Áp suất tĩnh định lớn nhỏ với vụ nổ dạng cầu 20 Hình 2.9 Lịch sử áp suất – thời gian sóng nổ 21 Hình 2.10 Mối quan hệ b Z theo Nelson Lam [30] 22 Hình 2.11 Mơ hình sóng n ổ tác dụng vào mục tiêu 23 Hình 2.12 Mơ hình phân tích tốn t ấm 27 Hình 2.13 Mơ hình tính r (i ) 28 Hình 2.14 Mơ hình tính c ận  28 Hình 2.15 Mơ hình tính c ận   31 Hình 3.1 Phần tử dạng chữ nhật chịu uốn c ùng bậc tự 34 Hình 3.2 Mơ hình vật liệu 46 Hình 4.1 Sơ đồ hình học chia lưới 16x16 56 Hình 4.2 Đồ thị khảo sát 1000x3000x100mm dao động tự 57 Hình 4.4 Đồ thị khảo sát 1000x3000x100mm dao động tự 59 Hình 4.5 Đồ thị khảo sát 1000x3000x150 (mm) dao động tự 60 Hình 4.6 Đồ thị chuyển vị ngang nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ theo bước thời gian t 61 Hình 4.7 Đồ thị chuyển vị ngang nút đỉnh 1000x3000x150 (mm) chịu tải trọng nổ theo bước thời gian t 62 Hình 4.8 Đồ thị chuyển vị ngang nút 1000x3000x100 (mm) chịu tải trọng nổ theo bước thời gian t 63 Hình 4.9 Đồ thị chuyển vị ngang nút 1000x3000x150 (mm) chịu tải trọng nổ theo bước thời gian t 64 Hình 4.10 Áp lực nổ gán vào nút đỉnh 1000x3000x100 (mm) theo t 65 Hình 4.11 Đồ thị chuyển vị ngang nút đỉnh 000x3000x100 (mm) chịu tải trọng nổ theo thời gian t 65 Hình 4.12 Đồ thị chuyển vị ngang nút 1000x3000x100 (mm) chịu tải trọng nổ theo theo thời t 66 Hình 4.13 Đồ thị chuyển vị ngang lớn nút đỉnh 1000x3000 (mm) chịu tải tr ày khác 68 Hình 4.14 Đồ thị ảnh hưởng bề dày đến chuyển vị ngang lớn nút lân cận cao trình đặt thuốc nổ xét khơng có phi tuyến 69 Hình 4.15 Đồ thị giá trị chuyển vị ngang nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị tỉ số cản khác 70 Hình 4.16 Đồ thị chuyển vị ngang nút lân cận cao tr ình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, ứng với giá trị tỉ số cản khác 71 Hình 4.17 Đồ thị giá trị chuyển vị ngang lớn nút đỉnh v 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị R0 khác 72 DANH MỤC CÁC HÌNH VẼ iii Hình 4.18 Đồ thị chuyển vị ngang nút lân cận cao tr ình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, ứng với giá trị R0 khác 73 Hình 4.19 Đồ thị giá trị chuyển vị ngang lớn nút đỉnh v 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá t rị hno khác 74 Hình 4.20 Đồ thị giá trị chuyển vị ngang nút gần 1000x3000x150 (mm) chịu tải trọng nổ, ứng với giá trị hno khác 75 Hình 4.21 Đồ thị giá trị chuyển vị ngang lớn nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị bno khác 76 Hình 4.22 Đồ thị giá trị chuyển vị ngang nút lân cận cao tr ình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, ứng với giá trị bno khác 77 Hình 4.23 Đồ thị giá trị chuyển vị ngang lớn nút đỉnh v nút 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị W khác 78 Hình 4.24 Đồ thị giá trị chuyển vị ngang nút lân cận cao tr ình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, ứng với giá trị W khác 79 Hình 4.25 Đồ thị giá trị chuyển vị ngang lớn nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị R0 W khác 81 Hình 4.26 Đồ thị giá trị chuyển vị ngang lớn nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị hno W khác 82 DANH MỤC CÁC HÌNH VẼ iv DANH MỤC CÁC BẢNG BIỂU Bảng 2.1 Kết biến đổi áp suất tĩnh đỉnh v áp suất động 18 Bảng 2.2 Hệ số chuyển đổi cho loại thuốc nổ theo Baker v cộng 19 Bảng 4.1 Bảng so sánh chuyển vị ngang nút đỉnh 53 Bảng 4.2 Bảng so sánh chuyển vị ngang nút 53 Bảng 4.3 Bảng so sánh chu kỳ dao động tự nhi ên (tấm 1000x3000x100 (mm)) 54 Bảng 4.4 Bảng so sánh chu kỳ dao động tự nhi ên (tấm 1000x3000x150 (mm)) 55 Bảng 4.5 Chuyển vị ngang U z lớn nút đỉnh 1000x3000x100 (mm) dao động tự với chuyển vị vận tốc ban đầu u0  , v0  m/s 57 Bảng 4.6 Chuyển vị ngang U z lớn nút đỉnh 1000x3000x150 (mm) dao động tự do, với chuyển vị vận tốc ban đầu u0  , v0  m/s 58 Bảng 4.7 Chuyển vị ngang U z lớn nút 1000x3000x100 (mm) dao động tự với chuyển vị vận tốc ban đầu u0  , v0  m/s 59 Bảng 4.8 Chuyển vị ngang U z lớn nút 1000x3000x150 (mm) dao động tự do, với chuyển vị vận tốc ban đầu u0  , v0  m/s 60 Bảng 4.9 Bảng giá trị chuyển vị ngang U z lớn nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ 61 Bảng 4.10 Bảng giá trị chuyển vị ngang U z lớn nút đỉnh 1000x3000x150 (mm) chịu tải trọng nổ 62 Bảng 4.11 Chuyển vị ngang U z lớn nút 1000x3000x100 (mm) chịu tải trọng nổ, chia theo lưới 16x16 62 Bảng 4.12 Chuyển vị ngang U z lớn nút 1000x3000x150 chịu tải trọng nổ, chia theo lưới 16x16 63 Bảng 4.13 Giá trị mô men giới hạn M u 67 Bảng 4.14 Chuyển vị ngang U z lớn nút lân cận cao tr ình đặt thuốc nổ 67 Bảng 4.15 Chuyển vị ngang U z lớn nút đỉnh 1000x3000 (m m) chịu tải trọng nổ, ứng với giá trị bề d ày khác 68 Bảng 4.16 Chuyển vị ngang U z lớn nút lân cận cao tr ình hno 69 Bảng 4.17 Chuyển vị ngang U z lớn nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị tỉ số cản khác 70 Bảng 4.18 Chuyển vị ngang lớn U z nút lân cận cao trình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, với giá trị tỉ số cản  khác 71 Bảng 4.19 Chuyển vị ngang U z lớn nút đỉnh nút 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị R0 khác 72 Bảng 4.20 Chuyển vị ngang lớn U z nút lân cận cao trình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, với giá trị R0 khác 73 Bảng 4.21 Chuyển vị ngang U z lớn nút đỉnh nút 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị hno khác DANH MỤC CÁC BẢNG BIỂU 74 v Bảng 4.22 Chuyển vị ngang lớn U z nút lân cận cao trình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, với giá trị hno khác Bảng 4.23 Chuyển vị ngang U z lớn nút đỉnh nút 1000x3000x100 (mm) chịu tải trọng nổ, ứng với g iá trị bno 75 76 Bảng 4.24 Chuyển vị ngang U z lớn nút lân cận cao trình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, với giá trị bno khác 77 Bảng 4.25 Chuyển vị ngang U z lớn nút đỉnh nút 1000x3000x100 (mm) chịu tải trọng nổ, ứng với giá trị W khác 78 Bảng 4.26 Chuyển vị ngang U z lớn nút lân cận cao trình đặt thuốc nổ 1000x3000x150 (mm) chịu tải trọng nổ, với giá trị W khác 79 Bảng 4.27 Chuyển vị ngang U z lớn nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ, với giá trị R0 W khác 80 Bảng 4.28 Chuyển vị ngang lớn nút đỉnh 1000x3000x100 (mm) chịu tải trọng nổ, với giá trị hno W khác 82 DANH MỤC CÁC BẢNG BIỂU PL - 19 Ux40(it)=q(3*40 -2); if abs(Ux40(it))>abs(Ux40max) Ux40max=Ux40(it); end Uy40(it)=q(3*40 -1); if abs(Uy40(it))>abs(Uy40max) Uy40max=Uy40(it); end Uz40(it)=q(3*40); if abs(Uz40(it))>abs(Uz40max) Uz40max=Uz40(it); end Ux41(it)=q(3*41 -2); if abs(Ux41(it))>abs(Ux41max) Ux41max=Ux41(it); end Uy41(it)=q(3*41 -1); if abs(Uy41(it))>abs(Uy41max) Uy41max=Uy41(it); end Uz41(it)=q(3*41); if abs(Uz41(it))>abs(Uz41max) Uz41max=Uz41(it); end Ux42(it)=q(3*42 -2); if abs(Ux42(it))>abs(Ux42max) Ux42max=Ux42(it); end Uy42(it)=q(3*42 -1); if abs(Uy42(it))>abs(Uy42max) Uy42max=Uy42(it); end Uz42(it)=q(3*42); if abs(Uz42(it))>abs(Uz42max) Uz42max=Uz42(it); end Ux49(it)=q(3*49 -2); if abs(Ux49(it))>abs(Ux49max) Ux49max=Ux49(it); end Uy49(it)=q(3*49 -1); if abs(Uy49(it))>abs(Uy49max) Uy49max=Uy49(it); end Uz49(it)=q(3*49); if abs(Uz49(it))>abs(Uz49max) Uz49max=Uz49(it); end Ux50(it)=q(3*50 -2); if abs(Ux50(it))>abs(Ux50max) Ux50max=Ux50(it); end Uy50(it)=q(3*50 -1); if abs(Uy50(it))>abs(Uy50max) Uy50max=Uy50(it); end Uz50(it)=q (3*50); if abs(Uz50(it))>abs(Uz50max) Uz50max=Uz50(it); end Ux51(it)=q(3*51 -2); if abs(Ux51(it))>abs(Ux51max) Ux51max=Ux51(it); Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 20 end Uy51(it)=q(3*51 -1); if abs(Uy51(it))>abs(Uy51max) Uy51max=Uy51(it); end Uz51(it)=q(3*51); if abs(Uz51(it))>abs(Uz51max) Uz51max=Uz51(it); end end Ux31max; Uy31max; Uz31max; Ux32max; Uy32max; Uz32max; Ux33max; Uy33max; Uz33max; Ux40max; Uy40max; Uz40max; Ux41max; Uy41max; Uz41max; Ux42max; Uy42max; Uz42max; Ux49max; Uy49max; Uz49max; Ux50max; Uy50max; Uz50max; Ux51max; Uy51max; Uz51max; disp('******************************************************************' ); disp('***********XAC DINH MOMENT CUA MOI PHAN TU TA M********************' ); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************'); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('************************ ******************************************' ); fprintf('\n'); for i=1:ny for j=1:nx ie=nx*(i-1)+j; if ie==28 Qe28=[ Ux31max; Uy31max; Uz31max; Ux32max; Uy32max; Uz32max; Ux40max; Uy40max; Uz40max; Ux41max; Uy41max; Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 21 Uz41max;]; QQe=double(Qe28); Se=matrantinhmomen(B,H,thk,muy,E); Mo28=double(abs(Se*QQe)); elseif ie==29 Qe29=[ Ux32max; Uy32max; Uz32max; Ux33max; Uy33max; Uz33max; Ux41max; Uy41max; Uz41max; Ux42max; Uy42max; Uz42max;]; QQe=double(Qe29); Se=matrantinhmome n(B,H,thk,muy,E); Mo29=double(abs(Se*QQe)); elseif ie==36 Qe36=[ Ux40max; Uy40max; Uz40max; Ux41max; Uy41max; Uz41 max; Ux49max; Uy49max; Uz49max; Ux50max; Uy50max; Uz50max;]; QQe=double(Qe36); Se=matrantinhmomen(B,H,thk,muy,E) ; Mo36=double(abs(Se*QQe)); elseif ie==37 Qe37=[ Ux41max; Uy41max; Uz41max; Ux42max; Uy42max; Uz42max; Ux50max; Uy50max; Uz50max; Ux51max; Uy51max; Uz51max;]; QQe=double(Qe37); Se=matrantinhmomen(B,H,thk,muy,E); Mo37=double(abs(Se*QQe)); end end end figure(1) plot(tii,Uz50,'LineWidth',2); xlabel('Thoi gian (s)'); ylabel('Chuyen vi ngang nut lan can cao trinh dat thuoc no cua tam (m)' ); title('Do thi chuyen vi ngang nut lan can cao trinh dat thuoc no cua tam theo t Bai toan tuyen tinh' ); grid on; fprintf('\n'); Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 22 disp('**********PHAN TICH TAM CO XET PHI TUYEN VAT LIEU **********' ); if thk==0.15 Muth=30;%kN.m Nho su ho tro cua phan mem GALA elseif thk==0.2 Muth=70; elseif thk==0.25 Muth=110; elseif thk==0.3 Muth=140; end if(Loaitaidong ==1) %Truong hop phan tich tai no Taino=subs(Pni,t,t0); for i=1:(ny+1) for j=1:(nx+1) nutgn=(i -1)*(nx+1)+j; if nutgn==1 % (i==1 & j==1) P(3*nut gn-2)=0; P(3*nutgn -1)=0; P(3*nutgn)=Taino(nutgn); elseif nutgn==(nx+1) % (i==1 & j==nx+1) P(3*nutgn -2)=0; P(3*nutgn -1)=0; P(3*nutgn)=Taino(nutgn); elseif nutgn==ny*(nx+1)+1 % (i==ny+1&j==1) P(3*nutgn -2)=0; P(3*nutgn -1)=0; P(3*nutgn)=Taino(nutgn); elseif nutgn==(ny+1)*(nx+1) % (i==(ny+1)& j==(nx+1)) P(3*nutgn -2)=0; P(3*nutgn-1)=0; P(3*nutgn)=Taino(nutgn); elseif nutgn==j && 1 Muth K=zeros(sobtd1tam,sobtd1tam); elseif ie==36 && Mo36(1) > Muth K=zeros(sobtd1tam,sobtd 1tam); elseif ie==37 && Mo37(1) > Muth K=zeros(sobtd1tam,sobtd1tam); else K=matrandocungphantu(B,H,thk,muy,E); M=matrankhoiluongphantu(ro,B,thk,H); end phtu(ie,1)=ie+(i -1); phtu(ie,2)=ie+i; phtu(ie,3)=ie+i+nx+1; phtu(ie,4)=ie+i+nx; ix=chisotam(phtu(ie,:),sonutphtu,sobtdnut); KOS=matrandocungtam(KOS,K,ix); MOS=matrankhoiluongtam(MOS,M,ix); end end disp('******************* ***********************************************' ); disp('************************** gandkxoa.m **************************' ); disp('******************** MA TRAN DIEU KIEN XOA *********************' ); disp('********** Hoc vien :TRAN MINH THI ****** *****************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************'); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 24 if Mo28(1) > Muth && Mo29(1) > Muth && Mo36(1) > Muth && Mo37(1) > Muth nutxoa=41; chisoxoa(nutxoa*3 -2)=nutxoa*3-2; giatribtdchisoxoa(nutxoa*3 -2)=0; chisoxoa(nutxoa*3 -1)=nutxoa*3-1; giatribtdchisoxoa(nutxoa*3 -1)=0; chisoxoa(nutxoa*3)=nutxoa*3; giatribtdchisoxoa(nutxoa*3)=0; [KOS,MOS,P]=gandkxoa(KOS,MOS,P,chisoxoa,giatribtdchisoxoa,nutxoa,sobtdnut); %gan dieu kien bien cho nut chinh giua end disp('******************************************************************' ); disp('************************** gandk bien.m **************************' ); disp('******************** MA TRAN DIEU KIEN BIEN *********************' ); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************'); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); for j=1:(nx+1) nutkhop1=j; chisorb1(nutkhop1*3)=nutkhop1*3; giatribtdchiso(nutkhop1*3)=0; nutkhop2=ny*(nx+1)+j; chisorb2(nutkhop2*3)=nutkhop2*3; giatribtdchiso(nutkhop2*3)=0; end for i=1:(ny+1) nutkhop3=9*i-8; chisorb3(nutkhop3*3)=nutkhop3*3; giatribtdchiso(nutkhop3*3)=0; nutkhop4=9*i; chisorb4(nutkhop4*3)=nutkhop4*3 ; giatribtdchiso(nutkhop4*3)=0; end [KOS,MOS,P]=gandkbienpt(KOS,MOS,P,chisorb1,chisorb2,chisorb3,chisorb4,giatr ibtdchiso,nx,ny);%gan dieu kien bien disp('******************************************************************' ); disp('********************** ***** tanso.m ******************************' ); disp('************************* TAN SO DAO DONG ***********************' ); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUO C ***************'); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); omega=[];omega1=[];omega2=[];omega3=[]; KK=KOS; MM=MOS/9.81; [V,D]=eig(KK,MM); [lamda,h]=sort(diag(D)); omega=sqrt(lamda); T=2*pi./omega; disp('******************************************************************' ); disp('************************* matrancan.m ****************************' ); Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 25 disp('*********************** MA TRAN CAN TONG THE ********** **********'); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************'); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('*********************************** *******************************' ); fprintf('\n'); if Mo28(1) > Muth && Mo29(1) > Muth && Mo36(1) > Muth && Mo37(1) > Muth omega1=omega((nx+1)*4); omega2=omega((nx+1)*4+1); omega3=omega((nx+1)*4+2); else omega1=omega((nx+1)*4 -3); omega2=omega((nx+1)*4-2); omega3=omega((nx+1)*4 -1); end T1=2*pi./omega1; T2=2*pi./omega2; T3=2*pi./omega3; aa0=2*xuy*omega1*omega2/(omega1+omega2); aa1=2*xuy/(omega1+omega2); COS=aa0*MM+aa1*KK; CC=COS; MM=MOS/9.81; KK=KOS; PP=P; % -% -GAN CAC DIEU KIEN BAN DAU -% -ut=u0*ones(3*sonut,1); %Vecto chuyen vi ban dau vt=v0*ones(3*sonut,1);%Vecto van toc ban dau at=inv(MM)*(PP'-KK*ut-CC*vt);%Vecto gia toc ban dau disp('******************************************************************' ); disp('*********************** tichphanNewmark.m ************************' ); disp('**************************TICH PHAN NEWMARK***********************' ); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************'); disp('********** Ten chuong trinh :BLASTSHEL L **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG *************************** '); disp('******************************************************************' ); fprintf('\n'); delta=0.25; theta=0.5; a0=1/(delta*Dt^2); a1=theta/(delta*Dt); a2=1/(delta*Dt); a3=1/(delta*2)-1; a4=theta/delta-1; a5=Dt*(theta/delta -2)/2; a6=Dt*(1-theta); a7=Dt*theta; Uz50pt=[]; Uz50pt(it)=ut(3*50); Uz50ptmax=Uz50pt(it); % -% - CAC VONG LAP THEO THOI GIAN -% Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 26 tii=[];tii(1)=0; it=1;tt=0; fprintf('\n Vong lap thu %d' ,it); while tt Muth nutxoa=41; chisoxoa(nutxoa*3 -2)=nutxoa*3-2;%xoa diem chinh giu a giatribtdchisoxoa(nutxoa*3 -2)=0; chisoxoa(nutxoa*3 -1)=nutxoa*3-1; giatribtdchisoxoa(nutxoa*3 -1)=0; chisoxoa(nutxoa*3)=nutxoa*3; giatribtdchisoxoa(nutxoa*3)=0; [KOS,MOS,P]=gandkxoa(KOS,MOS, P,chisoxoa,giatribtdchisoxoa,nutxoa,sobtdnut); %gan dieu kien xoa end [KOS,MOS,P]=gandkbienpt(KOS,MOS,P,chisorb1,chisorb2,chisorb3,chisorb4,giatr ibtdchiso,nx,ny);%gan dieu kien bien MM=MOS/9.81; KK=KOS; COS=aa0*MM+aa1 *KK; CC=COS; PP=P; PP=double(PP'); EE=KK+a0*MM+a1*CC; FF=PP+MM*(a0*ut+a2*vt+a3*at)+CC*(a1*ut+a4*vt+a5*at); uti=inv(EE)*FF; ati=a0*(uti -ut)-a2*vt-a3*at; vti=vt+a6*at+a7*ati; ut=uti; vt=vti; at=ati; ut=double(ut); vt=double(vt); at=double(at); dem=0; for i=1:(ny+1) for j=1:(nx+1) nutgn=(i -1)*(nx+1)+j; dem=dem+1; q(3*nutgn -2)=ut(3*dem-2); q(3*nutgn-1)=ut(3*dem-1); q(3*nutgn)=ut(3*dem); end Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 28 end for j=1:(nx+1) nutkhop1=j; %lien ket khop tren canh y=0 nutkhop2=ny*(nx+1)+j; %lien ket khop tren canh y=Ho dem=dem+1; q(3*nutkhop1)=0; q(3*nutkhop2)=0; end for i=1:(ny+1) nutkhop3=9*i -8;%lien ket khop tren canh x=0 nutkhop4=9*i; %lien ket khop tren canh x=Bo dem=dem+1; q(3*nutkhop3)=0; q(3*nutkhop4)=0; end Uz50pt(it)=q(3*50) ; if abs(Uz50pt(it))>abs(Uz50ptmax) Uz50ptmax=Uz50pt(it); end end figure(2) plot(tii,Uz50pt,'LineWidth',2); xlabel('Thoi gian (s)'); ylabel('Chuyen vi ngang nut lan can cao trinh dat thuoc no cua tam(m)' ); title('Do thi chuyen vi ng ang nut lan can cao trinh dat thuoc no cua tam theo t Bai toan phi tuyen vat lieu' ); grid on; end end disp('******************************************************************' ); fprintf('***********XUAT DU LIEU NHAP VA KET QUA PHAN TICH************ \n'); fprintf('***************BAI TOAN TAM CHIU TAI TRONG NO**************** \n'); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************'); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************'); disp('******************************************************************' ); fprintf('\n'); fprintf('*******************SO LIEU HINH HOC CUA TAM *******************' ); fprintf('Kich thuoc cua dinh theo phuong X: Bo =%g \n',Bo); fprintf('Kich thuoc cua dinh theo phuong Z: Ho =%g \n',Ho); fprintf('Be day cua dinh: thk =%g \n',thk); fprintf('************************************************************* \n'); fprintf('\n'); fprintf('*******************THONG SO VAT LIEU CUA TAM******************' ); fprintf('Khoi luong rieng cua be tong (kN/m3): ro =%g \n',ro); fprintf('Module dan hoi cua be tong (kN/m2): E =%g \n',E); fprintf('He so no hong cua be tong: muy =%g \n',muy); fprintf('Ti so can cua be tong: xuy =%g \n',xuy); fprintf('\n'); fprintf('*********************SO DOAN CHIA LUOI TAM********************' ); fprintf('So doan chia dinh theo phuong X: nx =%g \n',nx); fprintf('So doan chia dinh theo phuong Y: ny =%g \n',ny); fprintf('\n'); fprintf('*********************THONG SO TAI TRONG NO ********************' ); fprintf('Khoang cach tu dinh no den dinh (m): Ro =%g \n',Ro); fprintf('Trong luong thuoc no (kg): W =%g \n',W); fprintf('Chieu cao so voi mat dat cua vi tri thuoc no (m): hno =%g \n',hno); fprintf('Khoang cach so voi dinh O cua vi tri thuoc no (m): bno =%g\n',bno); Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 29 fprintf('Ap suat khong (kN/m2): po =%g \n',po); fprintf('\n'); if Baitoan==0 fprintf('********************************************************* \n'); fprintf('*******PHAN TICH TAM KHONG XET PHI TUYEN *************** *\n'); fprintf('********************************************************* \n'); elseif Baitoan==1 fprintf('********************************************************* \n'); fprintf('*********PHAN TICH TAM CO XET PHI TUYEN VAT LIEU********* \n'); fprintf('********************************************************* \n'); end if Baitoan==0 fprintf('**********************CHUYEN VI CUA NUT ************************' ); fprintf('Chuyen vi cua nut dinh giua bai toan tuyen tinh (m): Uzdinhmax =%g\n',Uzdinhmax); fprintf('Chuyen vi cua nut co cung cao trinh hno bai toan tuyen tinh (m): Uzhnomax =%g \n',Uzhnomax); fprintf('Chuyen vi cua nut chinh giua tam bai toan tuyen tinh (m): Uzgiuamax =%g \n',Uzgiuamax); elseif Baitoan==1 fprintf('**********************CHUYEN VI CUA NUT ************************' ); fprintf('Chuyen vi cua nut gan giua tam bai toan tuyen tinh (m): Uz50max =%g\n',Uz50max); fprintf('Chuyen vi cua nut gan giua tam bai toan phi tuyen vat lieu (m): Uz50ptmax =%g \n',Uz50ptmax); end if Baitoan==0 fprintf('**********************CHU KY DAO DONG TU DO ************************' ); fprintf('Chu ky dao dong : T1 =%g \n',T1); fprintf('Chu ky dao dong : T2 =%g \n',T2); fprintf('Chu ky dao dong : T3 =%g\n',T3); end Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 30 Các chương trình disp('******************************************************************' ); disp('*********************** matrandocungphantu.m *********************' ); disp('******************************************************************' ); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************' ); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************** **********'); fprintf('\n'); function [K]=matrandocungphantu(B,H,thk,muy,E) syms x y Qk=[1 x y x^2 x*y y^2 x^3 x^2*y x*y^2 y^3 x^3*y x*y^3]; Ak=inv([0 0 0 0 0 0; -1 0 0 0 0 0; 0 0 0 0 0 0; 0 B 0 B^2 0 B^3 0; -1 -2*B 0 -3*B^2 0 0 0; B B^2 0 B^3 0 0 0; 0 B 2*H B^2 2*B*H 3*H^2 B^3 3*B*H^2; -1 -2*B -H -3*B^2 -2*B*H -H^2 -3*B^2*H H^3; B H B^2 B*H H^2 B^3 B^2*H B*H^2 H^3 B^3*H B*H^3; 0 0 2*H 0 3*H^2 0; -1 0 -H 0 -H^2 0 H^3; H 0 H^2 0 H^3 0;]); Bk=[diff(Qk,x,2); diff(Qk,y,2); 2*diff(diff(Qk,x),y)]; D=[1 muy 0; muy 0; 0 (1 -muy)/2]; I=simplify(E*thk^3/(12*(1 -muy^2))*int(int(Bk.'*D*Bk,y,0,H),x,0,B)); K=simplify(Ak.'*I*Ak); fprintf('\n'); disp('******************************************************************' ); disp('******************** matrankhoiluongphantu.m *********************' ); disp('*************************************************************** ***'); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************' ); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :05/12/2009**************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************** **********************' ); fprintf('\n'); function M=matrankhoiluongphantu(ro,B,thk,H) syms x y z Qm=[1 x y x^2 x*y y^2 x^3 x^2*y x*y^2 y^3 x^3*y x*y^3]; Am=inv([ 0 0 0 0 0 0; -1 0 0 0 0 0; 0 0 0 0 0 0; 0 B 0 B^2 0 B^3 0; -1 -2*B 0 -3*B^2 0 0 0; B B^2 0 B^3 0 0 0; 0 B 2*H B^2 2*B*H 3*H^2 B^3 3*B*H^2; Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 31 -1 -2*B -H B H B^2 0 -1 0 0 -3*B^2 -2*B*H -H^2 -3*B^2*H - H^3; B*H H^2 B^3 B^2*H 0 0 B*H^2 H^3 B^3*H B*H^3; -H 2*H 0 -H^2 3*H^2 0 0; - H^3; H 0 H^2 0 H^3 0;]); Bm=[-z*diff(Qm,x); -z*diff(Qm,y); Qm]; M=simplify(ro*Am.'*int(int(int(Bm.'*Bm,y,0,H),x,0,B),z, thk/2,thk/2)*Am); fprintf('\n'); disp('******************************** **********************************' ); disp('********************** matrandocungtam.m **** *********************' ); disp('******************************************************************' ); disp('********** Hoc vien :TRAN MINH THI ***************** ******'); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************' ); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :0 5/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); function [KOS]=matrandocungtam( KOS,K,ix) for i=1:length(ix) is=ix(i); for j=1:length(ix) js=ix(j); KOS(is,js)=KOS(is,js)+K(i,j); end end fprintf('\n'); disp('***************** *************************************************' ); disp('******************* matrankhoiluongtam.m *************************' ); disp('******************************************************************' ); disp('********** Hoc vien :TRAN MINH THI ** *********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************' ); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :0 5/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); function [MOS]=matrankhoiluongtam(MOS,M,ix) for i=1:length(ix) is=ix(i); for j=1:length(ix) js=ix(j); MOS(is,js)=MOS(is,js)+M(i,j); end end fprintf('\n'); disp('******************************************************************' ); disp('************************* chisotam.m *****************************' ); disp('******************************************************************' ); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************' ); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :0 5/12/2009 ************** ************'); Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 32 disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); function [ix]=chisotam(phtu,sonutphtu,sobtdnut) k=0; for i=1:sonutphtu s=(phtu(i)-1)*sobtdnut; for j=1:sobtdnut k=k+1; ix(k)=s+j; end end fprintf('\n'); disp('****************************************************************** '); disp('************************ gandkbien.m *****************************' ); disp('******************************************************************' ); disp('********** Hoc vien :TRAN MINH THI ***********************' ); disp('********** GV Huon g dan :TS.NGUYEN TRONG PHUOC ***************' ); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :0 5/12/2009 **************************' ); disp('********************************************** ********************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); function [KOS,MOS,P]=gandkbien(KOS,MOS,P,chisorb,giatribtdchiso) m=length(chisorb); n=length(KOS); for i=1:m c=chisorb(i); for j=1:n KOS(c,j)=0; KOS(j,c)=0; MOS(c,j)=0; MOS(j,c)=0; end KOS(c,c)=1; MOS(c,c)=1; P(c)=giatribtdchiso(i); end disp('********************** ********************************************' ); disp('************************ gandkbienpt.m ***************************' ); disp('******************************************************************' ); disp('********** Hoc vien :TRAN MINH THI ******* ****************'); disp('********** GV Huong dan :TS NGUYEN TRONG PHUOC ***************' ); disp('********** Ten chuong trinh :BLASTSHELL **************************' ); disp('********** Ngay hoan :0 5/12/2009 **************************' ); disp('******************************************************************' ); disp('******************** He don vi kN -m-kG ***************************' ); disp('******************************************************************' ); fprintf('\n'); function [KOS,MOS,P]=gandkbienpt(KOS,MOS,P,chisorb1,chisorb2,chisorb3,chisorb4,giatr ibtdchiso,nx,ny) n=length(KOS); for i=1:(nx+1) c=chisorb1(3*i); for j=1:n KOS(c,j)=0; KOS(j,c)=0; MOS(c,j)=0; MOS(j,c)=0; end Phụ lục: Mã nguồn chương trình BLASTSHELL PL - 33 KOS(c,c)=1; MOS(c,c)=1; P(c)=giatribtdchiso(i); end for i=1:(nx+1) d=chisorb2((ny*(nx+1)+i)*3); for j=1:n KOS(d,j)=0; KOS(j,d)=0; MOS(d,j)=0; MOS(j,d)=0; end KOS(d,d)=1; MOS(d,d)=1; P(d)=giatribtdchiso(d); end for i=1:(ny+1) e=chisorb3((9*i -8)*3); for j=1:n KOS(e,j)=0; KOS(j,e)=0; MOS(e,j)=0; MOS(j,e)=0; end KOS(e,e)=1; MOS(e,e)=1; P(e)=giatribtdchiso(i); end for i=1:(ny+1) f=chisorb4(9*i*3); for j=1:n KOS(f,j)=0; KOS(j,f)=0; MOS(f,j)=0; MOS(j,f)=0; end KOS(f,f)=1; MOS(f,f)=1; P(f)=giatribtdchiso(f); end Phụ lục: Mã nguồn chương trình BLASTSHELL ... 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 [14] Bài báo nghiên cứu ứng xử động lực học khung bê tông cốt thép chịu tác động tải trọng nổ. .. [30] phân tích động lực học tuyến tính bậc tự bê tông cốt thép chịu tải trọng nổ Đỗ Kiến Quốc, Lê Đức Tuấn (2/2007) nghiên cứu ứng xử động lực học khung bê tông cốt thép tác dụng tải trọng nổ [14]... đề tải trọng nổ đ ược quan tâm nhiều nhà khoa học giới Trong luận văn này, ứng xử động lực học b ê tông cốt thép chịu tải trọng nổ đ ược nghiên cứu Đầu tiên, tải trọng nổ mô dạng tải trọng động