Nghiên cứu xác định một số thông số tối ưu của máy chữa cháy rừng bằng đất cát

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC LÂM NGHIỆP NGHIÊN CỨU XÁC ĐỊNH MỘT SỐ THÔNG SỐ TỐI ƯU CỦA MÁY CHỮA CHÁY RỪNG BẰNG ĐẤT CÁT LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜ NG ĐẠI HỌC LÂM NGHIỆ P - NGHIÊN CỨU XÁC ĐỊNH MỘT SỐ THÔNG SỐ TỐI ƯU CỦA MÁY CHỮA CHÁY RỪNG BẰNG ĐẤT CÁT LUẬN ÁN TIẾN SĨ KỸ THUẬT Ngành : Kỹ thuật khí Mã số : 52 01 03 Người hướng dẫn khoa học HÀ NỘI, NĂM 2022 i LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng tơi, hướng dẫn khoa học PGS TS Dương Văn Tài TS Trần Văn Tưởng Các kết nghiên cứu trình bày luận án trung thực, khách quan chưa cơng bố cơng trình khác Tôi xin cam đoan giúp đỡ cho việc thực luận án cảm ơn, thơng tin trích dẫn luận án rõ nguồn gốc Hà Nội, ngày… tháng……năm 2022 Hướng dẫn khoa học Tác giả luận án ii LỜI CẢM ƠN Xin chân thành cảm ơn nhà khoa học, quan nhiệt tình giúp đỡ tơi hồn thành luận án khoa học Trước hết xin bầy tỏ lòng biết ơn sâu sắc tới PGS TS với ý kiến đóng góp quan trọng dẫn khoa học quý giá trình thực cơng trình nghiên cứu Trân trọng cảm ơn Ban Giám hiệu Trường Đại học Lâm nghiệp, Phòng Đào tạo sau đại học giúp đỡ tạo điều kiện để tơi hồn thành luận án Trân trọng cảm ơn Lãnh đạo Trường Đại học Đồng Nai tạo điều kiện thuận lợi để tơi hồn thành luận án Trân trọng cảm ơn Khoa Cơ điện Công trình, Bộ mơn Cơng nghệ máy chun dùng Trường Đại học Lâm nghiệp giúp đỡ, tạo điều kiện thuận lợi cho tơi hồn thành nhiệm vụ học tập nghiên cứu Trân trọng cảm ơn nhà khoa học Trường Đại học Lâm nghiệp, Học viện Nông nghiệp Việt Nam, Học viện Kỹ thuật quân đóng góp ý kiến q báu để tơi hồn thành luận án Hà Nội, ngày… tháng……năm 2022 Tác giả luận án iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC BẢNG vii DANH MỤC HÌNH viii DANH MỤC CÁC KÝ HIỆU VÀ VIẾT TẮT xi MỞ ĐẦU Chương TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU 1 Tình hình cháy rừng giới Việt Nam 1 Tình hình cháy rừng giới 1 Tình hình cháy rừng Việt Nam Khái quát công nghệ chữa cháy rừng Các cơng trình nghiên cứu thiết bị chữa cháy rừng giới 10 Các cơng trình nghiên cứu thiết bị chữa cháy rừng Việt Nam 14 Mục tiêu nghiên cứu luận án 19 Nội dung nghiên cứu luận án 19 Nghiên cứu lý thuyết 19 Nghiên cứu thực nghiệm 20 Đối tượng nghiên cứu 20 Cấu tạo máy phun đất cát chữa cháy rừng Nguyên lý hoạt động 20 21 Thông số kỹ máy chữa cháy rừng đất cát 22 Một số tồn máy chữa cháy rừng đất cát 23 Khái quát đất rừng phục vụ cho chữa cháy 23 Đặc điểm đất rừng 23 Phân loại đất rừng: 24 Thành phần đất rừng 25 Một số tính chất lý đất rừng 26 Nguyên lý trình chữa cháy [5]; [22] Bản chất trình cháy 29 29 iv Các điều kiện trình cháy Các chế lý, hoá dập tắt đám cháy rừng 29 31 10 Phương pháp nghiên cứu 35 10 Phương pháp nghiên cứu lý thuyết 35 10 Phương pháp nghiên cứu thực nghiệm Kết luận chương 35 36 Chương CƠ SỞ XÁC ĐỊNH CÁC THÔNG SỐ TỐI ƯU CỦA MÁY CHỮA CHÁY RỪNG BẰNG ĐẤT CÁT 37 Đặc điểm yêu cầu hệ thống cắt đất 37 1 Đặc điểm trình cắt đất 37 2 Yêu cần kỹ thuật đất sau cắt phục vụ cho chữa cháy 38 2 Xây dựng mơ hình tính tốn hệ thống cắt đất, hút đất phun đất vào đám cháy 38 2 Mơ hình tính tốn hệ thống cắt đất, hút phun đất vào đám cháy 38 2 Cấu tạo nguyên lý hoạt động hệ thống cắt đất 40 Cơ sở lý thuyết tính tốn hệ thống cắt đất Nguyên lý cắt đất dạng búa 41 42 Quan hệ động học qua trình cắt đất 44 3 Quan hệ động lực học trình cắt đất Lực cắt đất 45 46 Thiết lập hệ phương trình chuyển động hệ: Đĩa – dao cắt kéo văng đất 48 Tính tốn rung động hệ thống cắt đất 60 Xây dựng mơ hình tính dao động hệ máy cắt đất dạng búa 60 Thiết lập hệ phương trình dao động 61 Phân tích lực kích động tác dụng lên trục đĩa thép 63 Cơ sở lý thuyết tính tốn hệ thống hút phun đất vào đám cháy 68 Xác định vận tốc dịng khơng khí đường ống thẳng đứng 68 Xác định vận tốc dịng khơng khí ống dẫn nằm ngang 70 Tính tốn quạt hút phun đất 71 Tính tốn hệ thống hút đất 78 v 5 Tính tốn hệ thống phun đất Kết luận chương 82 84 Chương KHẢO SÁT CÁC YẾU TỐ ẢNH HƯỞNG ĐẾN ĐỘNG LỰC HỌC CỦA MÁY CHỮA CHÁY RỪNG BẰNG ĐẤT CÁT 85 Khảo sát yếu tố ảnh hưởng đến công suất cắt- kéo văng đất 85 1 Ảnh hưởng bán kính động học đĩa thép đến công suất cắt -kéo văng đất 86 Ảnh hưởng chiều dài dao cắt đất đến công suất cắt- kéo văng đất 87 3 Ảnh hưởng mật độ khối lượng dao cắt đến công suất cắt - kéo văng đất 88 Khảo sát ảnh hưởng vận tốc góc đĩa thép lắp dao cắt đến đến công suất cắt - kéo văng đất 89 Khảo sát rung động hệ thống cắt đất 90 Xác định thông số đầu vào để khảo sát rung động hệ thống 90 2 Giải mơ phương trình vi phân dao động hệ thống cắt đất 91 3 Giải pháp giảm rung cho hệ thống cắt đất 97 3 Khảo sát thông số ảnh hưởng đến vận tốc áp lực quạt hút phun đất 98 3 Khảo sát ảnh hưởng góc lắp ráp đầu đến vận tốc quạt hút phun đất 98 3 Khảo sát ảnh hưởng góc lắp ráp đầu đến áp lực phun quạt hút 99 Xác định số thông số hợp lý máy chữa cháy rừng đất cát 99 Chương NGHIÊN CỨU THỰC NGHIỆM KIỂM CHỨNG MÔ HÌNH TÍNH TỐN LÝ THUYẾT VÀ XÁC ĐỊNH MỘT SỐ THÔNG SỐ TỐI ƯU CỦA MÁY CHỮA CHÁY RỪNG BẰNG ĐẤT CÁT 102 Mục tiêu, nhiệm vụ nghiên cứu thực nghiệm 102 1 Mục tiêu nghiên cứu thực nghiệm Nhiệm vụ nghiên cứu 102 103 Đối tượng nghiên cứu thực nghiệm 103 4 Phương pháp xác định khối lượng đất phun 108 Thiết bị, dụng cụ đo áp lực khối lượng đất phun vào đám cháy 109 vi Thiết bị khuếch đại 109 4 Kết nghiên cứu thực nghiệm kiểm chứng mơ hình tính tốn lý thuyết 110 4 Chuẩn bị thí nghiệm 110 4 Tổ chức tiến hành thí nghiệm 110 Xác định số thông số tối ưu máy chữa cháy rừng đất cát Chọn phương pháp nghiên cứu 115 115 Phương pháp xử lý kết thí nghiệm Kết thí nghiệm đơn yếu tố 118 123 Kết thực nghiệm đa yếu tố 133 Xác định giá trị tối ưu tham số ảnh hưởng 140 Phương pháp tìm giá trị tối ưu thơng số đầu vào 140 Kết giải toán tối ưu theo phương pháp hàm tỷ lệ tối ưu tổng quát 142 Xác định công suất động 143 Xác định công suất động hệ thống cắt đất 143 Tính tốn cơng suất động hệ thống hút phun đất 143 Thông số kỹ thuật máy chữa cháy rừng đất cát tính tốn tối ưu hồn thiện 144 Kết luận chương 146 KẾT LUẬN VÀ KIẾN NGHỊ 147 Kết luận 147 Kiến nghị 148 TÀI LIỆU THAM KHẢO 149 vii DANH MỤC BẢNG Bảng 1: Tổng hợp tình hình cháy rừng Việt Nam từ năm 2010 - 2020 Bảng 2: Thông số kỹ thuật máy chữa cháy rừng đất cát 22 Bảng 3: Trạng thái độ chặt đất, ứng với độ ẩm 27 Bảng Khả chịu tải loại đất khác 52 Bảng So sánh lực cắt đất tính tốn theo lý thuyết với kết thực nghiệm 114 Bảng Bảng so sánh gia tốc rung động cực đại hệ thống cắt đất mơ hình tính toán lý thuyết kết thực nghiệm 115 Bảng 3: Kết thí nghiệm ảnh hưởng đường kính đĩa thép D đến khối lượng đất phun 124 Bảng 4: Kết thí nghiệm ảnh hưởng đường kính đĩa thép D 125 Bảng 5: Kết thí nghiệm ảnh hưởng chiều dài dao cắt L đến hối lượng đất phun 126 Bảng 6: Kết thí nghiệm Ảnh hưởng chiều dài dao cắt L đến áp lực phun 127 Bảng 7: Kết thí nghiệm ảnh hưởng góc lắp ráp đầu β2 đến khối lượng đất phun 128 Bảng 8: Kết thí nghiệm ảnh hưởng góc lắp ráp đầu β2 130 Bảng 9: Kết thí nghiệm ảnh hưởng số lượng cánh quạt gió đến khối lượng đất phun 131 Bảng 10: Kết thí nghiệm ảnh hưởng số lượng cánh Z quạt gió đến áp lực phun 132 Bảng 11 Mức thí nghiệm thơng số đầu vào 134 Bảng 12 Ma trận thí nghiệm theo kế hoạch Boks - Benken 135 Bảng 13: Kết ảnh hưởng tham số đến hàm 136 khối lượng đất phun 136 Bảng 14: Kết ảnh hưởng tham số đến hàm áp lực phun 137 137 Bảng 15: Thông số kỹ thuật máy phun đất cát chữa cháy rừng chế tạo theo kết tính tốn thơng số tối ưu 145 CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA TÁC GIẢ Trần Văn Tưởng, Dương Văn Tài, , “Khảo sát động lực học hệ thống cắt đất máy chữa cháy rừng đất cát”, Tạp chí Cơng nghiệp nông thôn số 31/2018 trang 17-23 , Dương Văn Tài, “Dao động máy chữa cháy rừng đất cát”, Tạp chí Cơng nghiệp nơng thơn số 42/2021 trang 48-55 , Dương Văn Tài, “Xác định số thông số tối ưu máy chữa cháy rừng đất cát”, Tạp chí Cơng nghiệp nơng thơn số 42/2021 trang 56-61 PHỤ LỤC PHỤ LỤC 01 KẾT QUẢ THÍ NGHIỆM THĂM DỊ ẢNH HƯỞNG CỦA CÁC THAM SỐ ĐẾN KHỐI LƯỢNG VÀ ÁP LỰC ĐẤT PHUN VÀO ĐÁM CHÁY Chiều dài dao cắt L(cm) Góc lắp ráp đầu β2 (độ) 100 Số cánh quạt Z Khối lượng đất phun Kg/phút Áp lực đất phun (N) Đường kính đĩa thép D (cm) 15 18 1,78 15 100 18 1,69 15 100 18 1,82 15 100 18 1,73 15 100 18 1,81 15 100 18 1,71 15 100 18 1,74 15 100 18 1,86 15 100 18 1,72 10 15 100 18 1,81 11 15 100 18 1,71 12 15 100 18 1,83 13 15 100 18 1,75 14 15 100 18 1,78 15 15 100 18 1,79 16 15 100 18 1,88 17 15 100 18 1,73 18 15 100 18 1,81 19 15 100 18 1,76 20 15 100 18 1,89 21 15 100 18 1,71 22 15 100 18 1,81 23 15 100 18 1,71 24 15 100 18 1,88 25 15 100 18 1,69 26 15 100 18 1,72 17,8 18,4 17,5 16,1 17,2 16,7 18,1 17,9 17,7 17,7 18,6 18,1 17,8 17,3 18,1 17,8 17 17,2 18,4 18,5 17,7 17,1 18,1 17,8 18,4 17,5 Số TN 27 15 100 18 1,83 28 15 100 18 1,71 29 15 100 18 1,81 30 15 100 18 1,74 31 15 100 18 1,76 32 15 100 18 1,82 33 15 100 18 1,71 34 15 100 18 1,81 35 15 100 18 1,73 36 15 100 18 1,85 37 15 100 18 1,78 38 15 100 18 1,79 39 15 100 18 1,88 40 15 100 18 1,73 41 15 100 18 1,81 42 15 100 18 1,76 43 15 100 18 1,79 44 15 100 18 1,81 45 15 100 18 1,71 46 15 100 18 1,81 47 15 100 18 1,76 48 15 100 18 1,89 49 15 100 18 1,71 50 15 100 18 1,81 16,1 17,2 16,7 18,1 17,9 17,7 17,7 18,6 18,1 17,8 17,3 18,1 17,8 17 17,2 18,4 18,5 17,7 17,1 18,1 18,4 18,5 17,7 17,1 PHỤ LỤC 02 Q TRÌNH THÍ NGHIỆM ĐO MƠ MEN VÀ RUNG ĐỘNG CỦA MÁY ĐÀO ĐẤT PHỤ LỤC 03 Q TRÌNH THÍ NGHIỆM ĐO KHỐI LƯỢNG VÀ ÁP LỰC ĐẤT PHUN VÀO ĐÁM CHÁY PHỤ LỤC 04 DỤNG CỤ VÀ THIẾT BỊ THÍ NGHIỆM PHỤ LỤC 05 COD CỦA CHƯƠNG TRÌNH TÍNH CÔNG SUẤT CẮT- KÉO VĂNG ĐẤT function Dao_dong_Bua_va_Cong_suat_cat_dat clc; Ten_ghi = 'So_lieu_khao_sat_Hoa xlsx'; g=9 81 ; % m/s^2 gia toc Trong truong R =0 ; % m L0= 03; %m ro_b = 78 ; % kg/m khoi luong rieng bua m_b = 55 g r0_d = 40 ; % kg/m2 khoi luong rieng dia h = 01 ; % m sau nhat cat mG = 003 ; % kg khoi luong dat cat - keo - vang 2,7 10^3kg/m^3 Xic_dat = 150000; % N/m2 ung suat gioi han Dat Beta = pi; % rad goc cat Omega = 262 ; % rad/s J1 = r0_d*pi*R^4/2 ; % mo men quan tinh dia Pc = 2*pi*h^2*Xic_dat/2 % luc cat da % L = 08 ; % m %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ cla; T_gian =10 ; t =0:0 001:T_gian; Bang=zeros(); Bang1(1,1)= "Do dai bua - cm" ; Bang1(1,2)= " Cong suat cat_KVang_W"; % Bang2(1,1)= "Khoi luong rieng bua - kg/m" ; Bang2(1,2)= " Cong suat cat_KVang_W"; %========================================== for j = 1: % ro_b =ro_b0 + (j-1)*0 2*ro_b0; L = L0+ (j-1)*0 01 b=R+L ; %m [t,z]=ode45(@rhs1,t,[0 0 ]); % BD_Max= max(z(:,3)*10^3*R2) bay = size(z(:,1)); N_c = zeros(); tg1 = Pc*b*Omega ; tg2 = mG*b^3*Omega^2/(Beta*b-h) for k =1:bay g1 = z(k,2); g2 =z(k,4); tgt1 =2*(R^2*L +L^3/3)*Omega; tgt2 = (L^3/3)*(g1+g2); tgt3 = (R*L^2/2) *((2*Omega + g1)*cos(z(k,1)) +(2*Omega + g2)*cos(z(k,3))); TS = ro_b*(tgt1 + tgt2 + tgt3) + J1*Omega - Pc*h/Omega; tgb1 = 2*(L*R^2 + L^3/3) + R*L^2*(cos(z(k,1)) + cos(z(k,3))); MS = J1 + mG*b^2 + ro_b*tgb1; tg = tg1+ tg2*(TS/MS); N_c(k) = tg; end Cong_suat_cat_W_Max= max(N_c); Bang(j,1) = L*100 ; %Bang(j,1) = ro_b ; Bang(j,2)=Cong_suat_cat_W_Max; end % xlswrite(['So_lieu_khao_sat_Hoa xlsx'], Bang1,2,['A' num2str(2)]); % xlswrite(['So_lieu_khao_sat_Hoa xlsx'],Bang,2,['A' num2str(3)]); %========================================== % axis off plot(Bang(:,1),Bang(:,2),'R-'); %,'linewidth',1); % axis([9 10 -0 5]); axis on hold on xlabel("Do dai bua (cm)"); % xlabel("Khoi luong rieng bua (kg/m)" ylabel("Cong suat cat keo vang (W) "); grid on %********************************************* function xdot=rhs1(t,z) G1= z(2); G2 = z(4); dxdt_1 = z(2); dxdt_3 = z(4); tgT1 =2*(R^2*L +L^3/3)*Omega; tgT2 = (L^3/3)*(G1+G2); tgT3 = (R*L^2/2) *((2*Omega + G1)*cos(z(1)) +(2*Omega + G2)*cos(z(3))); A_TS = ro_b*(tgT1 + tgT2 + tgT3) + J1*Omega - Pc*h/Omega; tgM1 = 2*(L*R^2 + L^3/3) + R*L^2*(cos(z(1)) + cos(z(3))); A_MS = J1 + mG*b^2 + ro_b*tgM1; tgVP = (R/L)*3*Omega^2*b^3*mG/(Beta*b-h); tgVP = tgVP*t; tgVP2t= 5*(R/L)*Omega^2*sin(z(1))+ 5*(g/L)*sin(Omega*t + z(1)); tgVP4t= 5*(R/L)*Omega^2*sin(z(3))+ 5*(g/L)*sin(Omega*t + pi+ z(3)); tgVP2p = sin(z(1))*( (Omega +0 5*G1)*A_MS -A_TS)/A_MS^2; tgVP4p = sin(z(3))*( (Omega +0 5*G2)*A_MS -A_TS)/A_MS^2; VP2 = tgVP *tgVP2p -tgVP2t ; VP4 = tgVP *tgVP4p -tgVP4t ; %+++++++++++++++++++++++++++++++++++ dxdt_2 = VP2; dxdt_4 = VP4 xdot=[dxdt_1 ; dxdt_2 ; dxdt_3 ; dxdt_4 ]; enD %************************************** end PHỤ LỤC 06 COD CỦA CHƯƠNG TRÌNH TÍNH DAO ĐỘNG CỦA MÁY CẮT ĐẤT function Dao_dong_may_cat_dat clc; Ten_ghi = 'So_lieu_khao_sat_Hoa xlsx'; g=9 81 ; % m/s^2 gia toc Trong truong R =0 ; % m L= 07; %m ro_b = 78 ; % kg/m khoi luong rieng bua m_b = 55 g m_b = 055 ; % kg khoi luong bua r0_d = 40 ; % kg/m2 khoi luong rieng dia %r0_d = 80; h = 01 ; % m sau nhat cat mG = 003 ; % kg khoi luong dat cat - keo - vang 2,7 10^3kg/m^3 Xic_dat = 150000; % N/m2 ung suat gioi han Dat Beta = pi/2; % rad goc cat Omega = 262 ; % rad/s J1 = r0_d*pi*R^4/2 ; % mo men quan tinh dia Pc = 2*pi*h^2*Xic_dat/2 % luc cat dat b=R+L; %m %&&&&&&&&&&&&&&&&&&&&&&&&&&&&& L1= 13 ; % m a=0 10 ; % m L3=0 17; %m L2= 27 ; %m C10=200000 ; % N/m C2= 50000 ;% N/m m =9 ; % kg khoi luong máy Te_ta = pi/4; % Goc nghieng T0 = Beta/Omega ; % thoi gian xung luc T = pi/Omega ; % Chu ky xung luc s0 = T0/2 ; S = T/2; N = 5; a0 = T0 /T; ak=zeros(); for k=1:N ak(k) = 2*T/(T0*k^2*pi^2)*(1-cos(k*pi*T0/T)); end %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ cla; T_gian =10 ; t =0:0 001:T_gian; [t,z]=ode45(@Dao_dong_bua,t,[0 0 ]); bay = size(z(:,1)) CS_c = zeros(); N_TB=zeros(); tg1 = Pc ; tg2 = mG*b^2*Omega/(pi*b-h); for k =1:bay g1 = z(k,2); g2 =z(k,4); tgt1 =2*(R^2*L +L^3/3)*Omega; tgt2 = (L^3/3)*(g1+g2); tgt3 = (R*L^2/2) *((2*Omega + g1)*cos(z(k,1)) +(2*Omega + g2)*cos(z(k,3))); TS = ro_b*(tgt1 + tgt2 + tgt3) + J1*Omega - Pc*h/Omega; tgb1 = 2*(L*R^2 + L^3/3) + R*L^2*(cos(z(k,1)) + cos(z(k,3))); MS = J1 + mG*b^2 + ro_b*tgb1; tgD = tg1+ tg2*(TS/MS); tg = tgD*b*Omega; CS_c(k) = tg; XL_TB(k) = tgD; end Cong_suat_cat= max(CS_c); Xung_luc_TB = max(XL_TB) P_q = m_b *Omega^2*(R + 5*L); P_T = sqrt(P_q^2 + Xung_luc_TB^2); Bang=zeros(); T_gian1= 10; Bang1(1,1) = " C1- kN/m "; Bang1(1,2) = " Bien_do_Dia - (mm) " ; Bang1(1,3)= " Luc cam tay - (kG)"; for j = 1: 14 C1= C10+(j-1)*50000; t =0:0 001:T_gian1; [t,x]=ode45(@Dao_dong_May,t,[0 L2 ]); Bien_do_Max = max(x(:,1)*L2); % tinh theo m Luc_lo_xo = Bien_do_Max *C1; Bien_do_dia = Bien_do_Max*(L1+L2)/L2; Luc_Tay_cam = Luc_lo_xo*L3/L2; Bang(j,1) = C1/1000; Bang(j,2)= Bien_do_dia* 1000; Bang(j,3) = Luc_Tay_cam /10; end xlswrite(['So_lieu_khao_sat_Hoa xlsx'],Bang1,1,['A' num2str(2)]); xlswrite(['So_lieu_khao_sat_Hoa xlsx'],Bang,1,['A' num2str(3)]); plot(Bang(:,1),Bang(:,3),'R-'); %,'linewidth',1); xlabel("Do cung C1 (kN/m)"); ylabel("Luc tay cam tren (kG) "); % plot(t(:),x(:,1)*(L2+L1)*1000,'R-'); %,'linewidth',1); % axis([9 10 -1 15]); % axis on % hold on % xlabel("Thoi gian (s)") % ylabel("Bien cua truc dia theo phuong dung - (mm) ") grid on %********************************************* function Sidot=Dao_dong_bua(t,z) G1= z(2); G2 = z(4); dSidt_1 = z(2); dSidt_3 = z(4); tgT1 =2*(R^2*L +L^3/3)*Omega; tgT2 = (L^3/3)*(G1+G2); tgT3 = (R*L^2/2) *((2*Omega + G1)*cos(z(1)) +(2*Omega + G2)*cos(z(3))); A_TS = ro_b*(tgT1 + tgT2 + tgT3) + J1*Omega - Pc*h/Omega; tgM1 = 2*(L*R^2 + L^3/3) + R*L^2*(cos(z(1)) + cos(z(3))); A_MS = J1 + mG*b^2 + ro_b*tgM1; tgVP = (R/L)*3*Omega^2*b^3*mG/(pi*b-h); tgVP = tgVP*t; tgVP2t= 5*(R/L)*Omega^2*sin(z(1))+ 5*(g/L)*sin(Omega*t + z(1)); tgVP4t= 5*(R/L)*Omega^2*sin(z(3))+ 5*(g/L)*sin(Omega*t + pi+ z(3)); tgVP2p = sin(z(1))*( (Omega +0 5*G1)*A_MS -A_TS)/A_MS^2; tgVP4p = sin(z(3))*( (Omega +0 5*G2)*A_MS -A_TS)/A_MS^2; VP2 = tgVP *tgVP2p -tgVP2t ; VP4 = tgVP *tgVP4p -tgVP4t ; %+++++++++++++++++++++++++++++++++++ dSidt_2 = VP2; dSidt_4 = VP4 ; Sidot=[dSidt_1 ; dSidt_2 ; dSidt_3 ; dSidt_4 ]; end %************************************** function xdot=Dao_dong_May(t,x) phi = x(1); phi1 =x(2) ; y = x(3) ; y1 =x(4); tg11 = m*y^2 ; tg12 =m*y*phi; tg22 = m*(phi^2+1); Del = m^2*y^2; tgP1 = a0/2; for k=1:N tgP1 = tgP1 + ak(k)*cos(k*pi*t/S); end tgP1 = 2*P_T*tgP1; tgP22 = tgP1*(phi*cos(Te_ta)+sin(Te_ta)) - 2*m*phi*phi1*y1 m*y*phi1^2-C1*(y-a)*phi^2 -C2*(y-L2)+m*g*phi; tgP11 = tgP1*(y+L1)*cos(Te_ta)-2*m*y*y1*phi1 -C1*(y-a)^2*phi+m*g*y; Del_1= tgP11*tg22-tgP22*tg12; Del_2 = tg11*tgP22-tg12*tgP11; %+++++++++++++++++++++++++++++++++++ if Del~=0 dxdt_1 = x(2); dxdt_3 = x(4); dxdt_2 = Del_1/Del; dxdt_4 = Del_2/Del ; xdot=[dxdt_1 ; dxdt_2 ; dxdt_3 ; dxdt_4 ]; end end end ... 99 Xác định số thông số hợp lý máy chữa cháy rừng đất cát 99 Chương NGHIÊN CỨU THỰC NGHIỆM KIỂM CHỨNG MƠ HÌNH TÍNH TỐN LÝ THUYẾT VÀ XÁC ĐỊNH MỘT SỐ THÔNG SỐ TỐI ƯU CỦA MÁY CHỮA CHÁY RỪNG BẰNG ĐẤT... Cấu tạo máy phun đất cát chữa cháy rừng Nguyên lý hoạt động 20 21 Thông số kỹ máy chữa cháy rừng đất cát 22 Một số tồn máy chữa cháy rừng đất cát 23 Khái quát đất rừng phục vụ cho chữa cháy 23... cứu tính tốn tối ưu số thông số máy Xuất phát từ lý luận án tiến hành nghiên cứu với tên đề tài: "Nghiên cứu xác định số thông số tối ưu máy chữa cháy rừng đất cát? ?? Mục đích nghiên cứu luận án