Nghiên cứu xác định một số thông số làm việc chính của máy gieo hạt đậu tương kết hợp bón phân

HỌC VIỆN NÔNG NGHIỆP VIỆT NAM NGUYỄN CHUNG THÔNG NGHIÊN CỨU XÁC ĐỊNH MỘT SỐ THƠNG SỐ LÀM VIỆC CHÍNH CỦA MÁY GIEO HẠT ĐẬU TƯƠNG KẾT HỢP BÓN PHÂN LUẬN ÁN TIẾN SĨ NHÀ XUẤT BẢN HỌC VIỆN NÔNG NGHIỆP - 2021 HỌC VIỆN NÔNG NGHIỆP VIỆT NAM NGUYỄN CHUNG THÔNG NGHIÊN CỨU XÁC ĐỊNH MỘT SỐ THÔNG SỐ LÀM VIỆC CHÍNH CỦA MÁY GIEO HẠT ĐẬU TƯƠNG KẾT HỢP BÓN PHÂN N : Mã so: N ười ướ Kỹ t uật k í 52 01 03 dẫ k oa ọc: PGS.TS Lê Mi Lư TS N uyễ Xuâ T iết HÀ NỘI - 2021 LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng tơi, kết nghiên cứu trình bày luận án trung thực, khách quan chưa dùng để bảo vệ lấy học vị 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 Tôi xin cam đoan việc sử dụng kết nghiên cứu công bố vào phần nội dung luận án đồng tác giả chủ nhiệm đề tài đồng ý, cho phép sử dụng Hà Nội, ngày tháng năm 2021 Tác giả luận án Nguyễn Chung Thông iii LỜI CẢM ƠN Trong suốt thời gian học tập, nghiên cứu hồn thành luận án, tơi hướng dẫn, bảo tận tình thầy cô giáo, giúp đỡ của, động viên bạn bè, đồng nghiệp gia đình Nhân dịp hồn thành luận án, cho phép tơi bày tỏ lịng kính trọng biết ơn sâu sắc tới hai thầy hướng dẫn: PGS.TS Lê Minh Lư TS Nguyễn Xuân Thiết khoa Cơ – Điện, Học viện Nông nghiệp Việt Nam, người tận tình hướng dẫn, dành nhiều công sức, thời gian điều kiện cho suốt trình học tập thực đề tài Tôi xin chân thành cảm ơn Bộ môn Cơ học kỹ thuật, Khoa Cơ Điện, Ban Quản lý đào tạo, Ban Tổ chức cán bộ, Ban Giám đốc Học viện Nông nghiệp Việt Nam tạo điều kiện thuận lợi cho tơi học nâng cao trình độ đồng thời nơi đào tạo thời gian qua Tôi xin chân thành cảm ơn Khoa, Ban, Viện Học viện giúp đỡ chuyên môn tạo điều kiện thuận lợi cho q trình nghiên cứu, thực luận án Tơi xin chân thành cảm ơn động viên, ủng hộ, đóng góp ý kiến đồng nghiệp, thầy giáo, nhà khoa học nước q trình thực luận án Tơi xin cảm ơn thành viên gia đình, bạn bè ủng hộ, động viên, tạo điều kiện thuận lợi để hoàn thành luận án Xin trân trọng cảm ơn! Hà Nội, ngày tháng năm 2021 Nghiên cứu sinh Nguyễn Chung Thông MỤC LỤC Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục chữ viết tắt vi Danh mục ký hiệu vii Danh mục bảng x Danh mục hình xii Trích yếu luận án xv Thesis abstract xvii Phần Mở đầu 1.1 Tính cấp thiết đề tài 1.2 Mục tiêu nghiên cứu 1.2.1 Mục tiêu tổng quát 1.2.2 Mục tiêu cụ thể 1.3 Phạm vi nghiên cứu 1.4 Những đóng góp luận án .4 1.5 Ý nghĩa khoa học thực tiễn đề tài 1.5.1 Ý nghĩa khoa học 1.5.2 Ý nghĩa thực tiễn Phần Tổng quan vấn đề nghiên cứu .6 2.1 Tình hình sản xuất giới hóa canh tác đậu tương 2.1.1 Tình hình sản xuất giới hóa canh tác đậu tương giới 2.1.2 Tình hình sản xuất giới hóa canh tác đậu tương nước 2.1.3 Kỹ thuật canh tác đậu tương truyền thống 2.2 Quy trình hệ thống máy canh tác đậu tương theo hướng giới hóa đồng 11 2.2.1 Quy trình canh tác đậu tương 11 2.2.2 Nguồn động lực hệ thống máy canh tác .13 2.3 Phương án canh tác sơ đồ kết cấu tổng thể máy GBĐ 15 2.3.1 Phương án gieo hạt, bón phân 15 2.3.2 Sơ đồ kết cấu nguyên lý làm việc máy GBĐ 17 2.4 Tổng quan nguyên lý gieo hạt bón phân, số mẫu máy gieo hạt bón phân 19 2.4.1 Một số nguyên lý phận gieo hạt 19 2.4.2 Một số nguyên lý phận bón phân 23 2.4.3 Một số mẫu máy gieo 26 2.5 Lựa chọn phận gieo phận bón phân 34 2.5.1 Lựa chọn phận gieo 34 2.5.2 Lựa chọn phận bón phân 37 Kết luận phần 40 Phần Phương pháp nghiên cứu 42 3.1 Nội dung nghiên cứu 42 3.2 Phương pháp nghiên cứu 42 3.2.1 Phương pháp giải tích 42 3.2.2 Phương pháp số 43 3.2.3 Phương pháp mơ hình hóa mơ hệ thống kỹ thuật 43 3.2.4 Phương pháp phần tử rời rạc số (Discrete element method) .44 3.2.5 Phương pháp nghiên cứu xác định số thơng số hình học lý tính .44 3.2.6 Phương pháp nghiên cứu thực nghiệm đơn yếu tố đa yếu tố .50 3.2.7 Phương pháp thử nghiệm máy 58 Phần Kết nghiên cứu 59 4.1 Mơ hình tốn xác định thông số phận gieo 59 4.1.1 Mơ hình xác định quy luật chuyển động hạt 60 4.1.2 Mơ hình xác định khả lấy hạt 66 4.1.3 Mơ hình xác định khả phân ly giữ hạt 71 4.1.4 Mô hình xác định tác động lưỡi gạt phân ly .73 4.1.5 Mơ hình tốn xác định khả nhả hạt đĩa gieo 75 4.1.6 Đĩa cố định ống dẫn hạt 77 4.2 Mơ hình tốn xác định thơng số phận bón phân 79 4.2.1 Xác định lượng phân qua cửa xả .79 4.2.2 Năng suất hệ số nạp trục .81 4.3 Nghiên cứu ổn định liên hợp máy gbđ 86 4.3.1 Mơ hình ổn định liên hợp máy GBĐ 86 4.3.2 Kết mô .92 4.4 Kết nghiên cứu thực nghiệm số thông số hình học lý tính 99 4.4.1 Kết nghiên cứu thực nghiệm đậu tương 99 4.4.2 Kết nghiên cứu thực nghiệm phân bón 109 4.4.3 Kết nghiên cứu số tham số nghiên cứu ổn định máy 110 4.5 Kết nghiên cứu thực nghiệm đơn đa yếu tố 111 4.5.1 Kết nghiên cứu thực nghiệm đơn yếu tố phận gieo 111 4.5.2 Kết nghiên cứu thực nghiệm đơn yếu tố phận bón phân 119 4.5.3 Kết nghiên cứu thực nghiệm đa yếu tố phận gieo 123 4.6 Kết sơ thử nghiệm máy 133 4.6.1 Thời gian, địa điểm điều kiện thử nghiệm 133 4.6.2 Kết thử nghiệm 135 Kết luận phần 135 Phần Kết luận kiến nghị .138 5.1 Kết luận 138 5.2 Kiến nghị 140 Danh mục cơng trình cơng bố có liên quan đến luận ấn 141 Tài liệu tham khảo 142 Phụ lục 146 DANH MỤC CHỮ VIẾT TẮT Từ viết tắt Nghĩa tiếng Việt DEM Discrete element method (Phương pháp phần tử rời rạc) ĐTĐL Đề tài độc lập FAOSTAT The Food and Agriculture Organization Corporate Statistical Database (Tổ chức Lương thực Nông nghiệp Liên Hiệp Quốc) GBĐ Tên máy gieo hạt đậu tương kết hợp bón phân GMO Genetically Modified (Cải biến di truyền) HP Mã lực LHM Liên hợp máy NN-PTNT Nông nghiệp Phát triển nông thôn USD United States Dollar (Đô la Mỹ) TTXVN Thông xã Việt Nam VNĐ Việt Nam đồng DANH MỤC KÝ HIỆU Ký hiệu Đơn vị Ý nghĩa ρvl kg/m3 Khối lượng thể tích hạt đậu tương phân bón mvl kg Khối lượng hạt chiếm chỗ cốc thí nghiệm Vvl m3 độ λ Thể tích chiếm chỗ cốc thí nghiệm Góc ma sát vật liệu với hạt đậu tương phân bón (tgλ = f) vm m/s Tốc độ làm việc máy MSF Trung bình bình phương nhân tố MSE Trung bình bình phương phương sai α Mức ý nghĩa f1 Bậc tự nhân tố f2 Bậc tự sai số N Tổng số thí nghiệm ̅ Giá trị trung bình chung ̅ Giá trị trung bình mẫu i (mức i) Xi Giá trị thí nghiệm thơng số vào thứ i X0i Mức sở thông số vào thứ i xi Giá trị mã chuẩn hóa thông số vào thứ i ɛi Khoảng biến thiên thông số thứ i a0, ai, aij Hệ số phương trình hồi quy (i ≠ j) b0, bi, bij Các ước lượng hệ số hồi quy Phương sai S2y S S2b0, S2bi, S,2bij B bii Sai số hệ số hồi quy theo phương sai m Bề rộng mặt luống H m Chiều cao luống R m Chiều rộng rãnh luống Zh hàng Số hàng gieo Zp hàng Số hàng phân bk m nhh hạt/hốc ah m Qp kg/ha Lượng phân ωđ rad/s Vận tốc góc đĩa gieo Dc m Khoảng cách hai hàng kép Số hạt hốc Khoảng cách hốc Đường kính cốc định lượng fc Hệ số ma sát hạt với cánh gạt vlđ m/s Vận tốc lỗ đĩa δđ m Chiều dầy đĩa gieo Δvm m Chiều dài cạnh vát mép dlđ m Đường kính lỗ đĩa dh m Đường kính hạt đậu tương Zlđ lỗ Số lỗ đĩa gieo e m Khe hở hai đĩa Phh N Lực phá hủy hạt Ppl N Lực phân ly hạt Kpl N/m fx Độ cứng phân ly Hệ số chảy cửa xả phân dp kg/m Khối lượng thể tích phân ex m Chiều rộng cửa xả phân bx m Chiều dài cửa xả phân Ψn % Hệ số nạp đầy Zt rãnh dr m Đường kính rãnh trục Lt m Chiều dài rãnh trục nt vòng/phút Số vòng quay trục Dt m Đường kính ngồi trục Qt kg/h Số rãnh trục suất làm việc trục 15 Lê Minh Lư (2002) Nghiên cứu dao động máy kéo bánh có tính đến đặc trưng phi tuyến phần tử đàn hồi Luận án tiến sĩ, Trường Đại học Nông nghiệp I 16 Lê Minh Lư (2019) Nghiên cứu, thiết kế, chế tạo thử nghiệm hệ thống máy, thiết bị phục vụ giới hóa đồng khâu canh tác đậu tương Báo cáo đề tài nghiên cứu khoa học độc lập cấp Nhà nước ĐTĐL.CN – 21/15 Bộ Khoa học Công nghệ 17 Lê Quyết Tiến, Nguyễn Xuân Biên, Hoàng Nghĩa Đạt & Nguyễn Xuân Thủy (2018) Kết nghiên cứu thiết kế chế tạo máy gieo ngơ, bón phân viên nén đất lúa hiệu thích ứng với biến đổi khí hậu đồng Sơng cứu long Tạp chí Cơng nghiệp nơng thơn 30: 35-42 18 Lê Sĩ Hùng (2003) Máy gieo ngô MGBN - 04, thuộc đề tài cấp bộ: Nghiên cứu số máy móc phục vụ cho canh tác trồng cạn đất dốc Bộ Nông nghiệp Phát triển Nông thôn 19 Mai Loan (2009) Khoa học/Máy canh tác đa - sản phẩm Truy cập từ: https://vietnamnet.vn/khoahoc/2009/07/860326/ ngày 27/7/2009 20 Ngơ Thế Dân, Trần Đình Long, Trần Văn Lài, Đỗ Thị Dung & Phạm Thị Đào (1999) Cây đậu tương NXB Nông nghiệp, Hà Nội 21 Nguyễn Bảng & Đồn Văn Điện (1978) Cấu tạo máy nơng nghiệp NXB Đại học Trung học chuyên nghiệp 22 Nguyễn Bảng (1995) Lý thuyết máy canh tác NXB trường Đại học Nông nghiệp I 23 Nguyễn Chung Thông (2019) Nghiên cứu giới hóa khâu gieo hạt đậu tương kết hợp với bón phân Báo cáo đề tài nghiên cứu khoa học, đề tài nhánh cấp Nhà nước Học viện Nông nghiệp Việt Nam 24 Nguyễn Huy Bích (2016) Cơng tác đào tạo nguồn nhân lực khí nơng nghiệp Kỷ yếu hội thảo CLB Khoa học – Công nghệ trường Đại học Kỹ thuật lần thứ 48 NXB Đại học Cần thơ 25 Nguyễn Khắc Thường & Đặng Thế Huy (1983) Nguyên lý máy NXB Nông thôn 26 Nguyễn Văn Đồn (1997) Nghiên cứu sở tính tốn thiết kế máy bóc phân ly vỏ hạt trẩu, lai Luận án tiến sĩ, Viện Cơ điện Nông nghiệp 27 Nguyễn Văn Mai, Nguyễn Đức Thật, Ngơ Vi Bình & Bùi Thanh Kỳ (2013) Tiến kỹ thuật, công nghệ đề nghị cơng nhận:Máy gieo ngơ kết hợp bón phân hai hàng (MGH-02D), ứng dụng vùng trung du, miền núi Viện Cơ điện Nông Công nghệ sau thu hoạch 28 Nguyễn Văn Muốn, Nguyễn Viết Lầu, Trần Văn Nghiễn & Hà Đức Thái (1999) 29 30 31 32 Máy canh tác nông nghiệp NXB Giáo dục Phạm Thượng Hàn, Nguyễn Trọng Quế & Nguyễn Văn Hòa (2002) Kỹ thuật đo lường đại lượng vật lý NXB Giáo dục Phạm Văn Lang & Bạch Quốc Khang (1998) Cơ sở lý thuyết quy hoạch thực nghiệm ứng dụng kỹ thuật nông nghiệp NXB Nông nghiệp, Hà Nội Phạm Văn Tờ & Lương Văn Vượt (2004) Cơ học lý thuyết NXB Đại học Nông nghiệp, Hà Nội Sở Khoa học Công nghệ tỉnh Thừa Thiên Huế (2009) Hà Tây: sáng chế thành công máy gieo hạt đậu tương Truy cập từ https://skhcn.thuathienhue.gov.vn/? gd=6&cn=293&tc=4832 ngày 28/12/2009 33 Tạp chí Cơng thương – Bộ Cơng thương (2019) Thị trường đậu tương giới tháng 2/2019 Truy cập từ http://www.tapchicongthuong.vn/bai-viet/thi-truong-dautuong-the-gioi-thang-22019-60122.htm ngày 10/02/2019 34 Tổng cục Thống kê (2017) Niên giám thống kê năm 2016 NXB Thống kê 35 Tổng cục Thống kê (2020) Niên giám thống kê năm 2019 NXB Thống kê 36 Trần Thị Nhị Hường & Đặng Thế Huy (1987) Một số phương pháp toán học học nông nghiệp NXB Nông nghiệp, Hà Nội 37 Trần Văn Điền (2007) Cây đậu tương NXB Nông nghiệp Hà Nội 38 Trung tâm Thông tin Công nghiệp Thương mại (2020) Sản lượng đậu tương giới niên vụ 2019-2020 Truy cập từ http://www.asemconnectvietnam.gov.vn/default.aspx?ID1=1&ZID1=21&ID8=936 39 ngày 25/01/2020 39 Trung tâm Thông tin Công nghiệp Thương mại Bộ Công Thương (2019) Nhập khô đậu tương năm 2019 dự kiến tăng nhẹ so với năm 2018 Truy cập từ http://thongtincongthuong.vn/nhap-khau-kho-dau-tuong-trong-nam-2019du-kien-tang-nhe-so-voi-nam-2018/ngày 11/11/2019 40 Trung tâm Thông tin Công nghiệp Thương mại Bộ Công Thương (2020) Nhập đậu tương Việt Nam tăng mạnh Truy cập từ http://thongtincongthuong.vn/nhap-khau-dau-tuong-ve-viet-nam-tang-manh/ ngày 6/5/2020 41 Vũ Bá Minh & Hồng Minh Nam (1998) Q trình thiết bị cơng nghệ hóa học (Tập 2: Cơ học vật liệu rời) NXB Khoa học Kỹ thuật 42 Vũ Đình Chính & Bùi Thị Cúc (2010) Nghiên cứu thời vụ gieo trồng cho đậu tương Xuân đất đồi gị huyện Chương Mỹ, Hà Nội Tạp chí Kinh tế Sinh thái 36: 33-47 43 Vũ Đình Chính & Lê Thị Lý (2011) Ảnh hưởng liều lượng phân bón đến sinh trưởng, phát triển suất số giống đậu tương đất phù sa đê huyễn Vĩnh Tường, tỉnh Vĩnh Phúc Tạp chí Khoa học Phát triển 1(9): 42-58 Tiếng Anh: 44 Adam B Flaugh (2011) Studies on Power Hop Using a onlinear Model of an Agricultural Tractor Degree of Master of Science University of Nebraska, US 45 Agritechno Yazaki (2016) Product/Technology/Services Retrieved at https://jgoodtech.jp/en/web/page/corp/-/info/JC0000000000751/appeal/eng on 16 January, 2019 46 Claas Harvest Centre (2017) New Machine/Seeding Retrieved at https://www.claasharvestcentre.com/new-machinery/seeding/rapid-a-400-800s on 12 January, 2019 47 Deere & Company (2017) Equipment/Agriculture/Seeding Equipment Retrieved at https://www.deere.com/en/seeding-equipment/730-air-disk-drill/# on 22 January, 2019 48 Dekenag (2017) Planter & Harvester/Seeding machine Retrieved at http://www.dekenag.com/productMechanies/2BYF-series-Corn-Planter-265.html 49 FAOSTAT (2019) Retrieved at http://FAOSTAT.fao.org on 10 January, 2020 50 Hoffman J.D (1994) Numerical methods for engineers and scientists Mc GrawHill, Inc, New York St Louis San Fancico 51 Jang Automation Co.,Ltd (2017) Product/Air Seeding Machine Retrieved at http://www.jangauto.com/vie/sub.php?code=02_prod0401&mode=view&no=36&c ate=20&page=1 on 10 January, 2019 52 John B (2005) Runge - Kutta methods for ordinary differential equations COE Workshop on Numerical Analysis Kyushu University The University of Auckland New Zealand 53 Kuhn (2017) Farm machinery, Seeding, Pneumatic seed drills Retrieved at https://www.kuhn.com/com_en/range/seeding.html on 10 January, 2019 54 Cundall P A & Strack O D L (1979) A discrete numerical model for granular assemblies Geotechnique 29(1): 47–65 https://doi.org/10.1680/geot.1979.29.1.47 55 Wang Z., Xi-wen L., Xiang T., Gou-zhong Z & Shan Z (2010) Precision Rice Hill-Direct – Seeding Technology and Machine Based on the Combination of Agricultural Machinery and Agronomic Technology Joural of South China Agriculture University 31(1) Tiếng Nga: 56 А.Б.ЛУРЬЕ И А.А.ГРОМБЧЕВСКИЙ (1977) РАСЧЕТ И КОНСТРУИРОВАНИЕ СЕЛЬСКОХОЗЯЙСТВЕННЫХ МАШИН ЛЕНИНГРАД «МАШИНОСТРОЕНИЕ» ЛЕНИНГРАДСКОЕ ОТДЕЛЕНИЕ PHỤ LỤC Các bước tính phân tích phương sai (một nhân tố) Tính tỷ số: F= MSF 1P MSE Tra bảng tra theo tiêu chuẩn thống kê Fisher (C) mức ý nghĩa α với bậc tự (r - 1, N - r), F < C giả thiết H (có sai khác) chấp nhận ngược lại bác bỏ giả thiết H0 Trong công thức 1P, MSF MSE trung bình bình phương nhân tố trung bình bình phương phương sai, xác định sau: SSF MSF r − 2P = 3P SS T2 MSE E = T2 T2 T2 N− r r SSF = ∑ N − y̅ ) = + + + r − (y̅ i i N1 N2 Nr N i=1 T2 T2 SSE = ∑ ∑ y2 − ( + ij N1 N2 T2 + + Nr r ) 4P 5P đó: (f1= r - 1) bậc tự nhân tố, (f2= N – k) bậc tự sai số, SSF tổng bình phương nhân tố, SSE tổng bình phương sai số, T tổng giá trị, N tổng số thí nghiệm (N = m.r), y̅ giá trị trung bình chung, y̅ i giá trị trung bình mẫu i (mức i) Các bước tính phân tích tương quan Phương pháp phân tích tương quan cho phép xác định mối quan hệ phụ thuộc hàm tiêu Y vào thông số vào X phương trình hồi quy Từ xây vẽ đồ thị quan hệ để phân tích, đánh giá lựa chọn vùng làm việc thông số vào để hàm tiêu có giá trị gần tối ưu Nếu Y có hồi quy tuyến tính theo X, đường thẳng có dạng: y = ax+ b 6P đó: a b gọi hệ số hồi quy, xác định theo công thức: N ∑ xy − (∑ x) (∑ y) a= N (∑ x2) − (∑ x)2 7P ∑ y − a ∑ x b= N 8P Theo phương pháp đoán giá trị hàm mục tiêu Y theo giá trị thông số vào X Các bước tính nghiên cứu thực nghiệm đa yếu tố Bước Xây dựng mơ hình hồi quy xác định hệ số Hàm hồi quy mô tả phụ thuộc hàm tiêu Y vào thông số ảnh hưởng X biểu diễn phương trình: Y 2= a0 + a1 X1 + ⋯ + Xi + aij Xi Xj + ⋯ + a11X2 + ⋯ + aii X 9P i đó: Y giá trị hàm tiêu, Xi giá trị thí nghiệm thơng số vào thứ i; ai, aij hệ số phương trình hồi quy (i ≠ j) Để thuận lợi cho việc xử lý số liệu việc xác định hệ số mơ hình, thay giá trị thực Xi giá trị mã hóa xi theo cơng thức: x Xi − X0i i = ε i , (i = 1, 2, k) đó: X0i mức sở thơng số vào thứ i, xi giá trị mã chuẩn hóa thông số vào thứ i đại lượng không thứ nguyên, ɛi khoảng biến thiên thông số thứ i Giá trị thí nghiệm thơng số X i chia làm mức: mức trên, sở mức Tương ứng với giá mã hóa xi nhận trị số là: xit = 1, x0i = xid = -1 Phương trình hồi quy thực nghiệm dạng mã viết theo công thức: k k-1 Y = b + ∑b x + ∑ x2 k ∑b i i i=1 i=1 j=i+1 k x x + + ∑ b ij i j i=1 ii 10P i đó: b0, bi, bij ước lượng hệ số hồi quy (gọi tắt hệ số hồi quy), k số thông số vào ảnh hưởng xi Số hệ số phương trình hồi quy xác định theo cơng thức: (k + 1)(k + 2) ktn = Tổng số thí nghiệm N xác định phần: + Phần sở (hạt nhân) phần thực nghiệm mức với số thí nghiệm 2k 2k-p (p hệ số rút gọn) + Phần thực nghiệm điểm có số thí nghiệm 2k + Phần thực nghiệm tâm hay mức với số thí nghiệm N0 ≥ Các giá trị p, α N0 xác định theo kế hoạch thực nghiệm cụ thể Số lần lặp lại điểm thí nghiệm thường m = Như vây, tổng số thí nghiệm kế hoạch thí nghiệm: N = 2k−p + k + N0 11P Bước Mức ý nghĩa hệ số hồi quy Mức ý nghĩa hệ số hồi quy kiểm tra theo tiêu chuẩn thống kê Student với mức ý nghĩa α = 0,05, bậc tự f = N(m – 1) tra bảng tìm tb Các hệ số hồi quy phương trình có nghĩa hệ số tb0, tbi, tbij tbii lớn tb, ngược lại hệ số ảnh hưởng không rõ rệt lên hàm mục tiêu ta loại bỏ hệ số Các bước tính bao gồm: tính phương sai S 2y sai số S2b0, S2bi, S2 bij, S2 bii hệ số hồi quy theo phương sai sai số thí nghiệm Phương sai tính theo công thức: N m S = ∑ ∑(y − y̅ ) y uj N(m − u u=1 j=1 1) Sai số hệ số hồi quy tính theo cơng thức: Sb0 = Sy k1 ; Sbi = Sy k5 Sbij = Sy k6 ; Sbii = Sy k3  k4 Các hệ số tb0, tbi, tbij tbii theo công thức: tb0 = b0 ; b0 S t = bij ; bij S bi j t bi = bi Sbi = bii t bii S bii 12 P Bước Kiểm tra tính thích ứng mơ hình tốn Sự thích ứng mơ hình tốn kiểm tra theo tiêu chuẩn thống kê Fisher theo cơng thức 13P Tra bảng tìm F b với α = 0,05 bậc tự do: f1 = N – k* f2 = N(m -1) S2 tu F = Sy2 13P đó: S phương sai thích ứng xác định theo cơng thức 14P, tu k* số hệ số hồi quy có nghĩa Nếu F < Fb mơ hình bậc thích ứng Nếu F > Fb mơ hình bậc khơng thích ứng, ta phải chuyển sang phương án bậc cao N S2 = tu ∑ (y u=1 u − yu )2 14P N−k* đó: yu giá trị tính tốn từ phương trình hồi quy Bước Chuyển dạng mã dạng thực Để xử lý số liệu nghiên cứu quy hoạch thực nghiệm cần sử dụng giá trị thông số ảnh hưởng dạng mã x i Tuy nhiên, để xác định giá trị hàm tiêu phụ thuộc vào thông số vào dạng tự nhiên ta thường chuyển hàm tiêu từ dạng mã dạng thực k k−1 Y = a + ∑a X + ∑ o i i i=1 k k ∑a X X + ∑a X2 ij i i=1 j=i+1 j ii i i=1 Các hệ số xác định theo công thức: ao = bo bi − ∑ εi X + 0i b ∑∑ k i ij j=1 i=1 X0 bi b ∑ j a k =i εi − ε j ε i j= b aij = ij εiεj ; i ε2 j = i εiεj X − 2b ii i bii εi2 X0i X0 j 0i (với i ≠ j) 15P Bảng 1P Thơng tin máy thiết bị thí nghiệm TT Tên thiết bị Động điện AC2 32-21 Ảnh thiết bị Thông số kỹ thuật - Số vịng quay: 2880 vịng/phút - Cơng suất: kW - Hiệu suất: 83 % - Khối lượng: 43 kg Biến tần VicRuns VD120-4T-4.0GB Máy đo tốc độ quay cầm tay EXTECH 461920 Thiết bị đo góc ma sát Gunt TM225 - Công suất: kW - Điện áp: ÷ 380 V - Tần số: ÷ 650 Hz - Phạm vi đo: ÷ 9999,9 vịng/phút - Độ phân giải: 0,1 vòng/phút - Khoảng cách mục tiêu: 500 mm - Phạm vi góc đo: ÷ 900 - Ma sát: nghỉ ma sát lăn - Kích thước DxRxC: 100x20x70 cm - Phạm vi đo: ÷ 20 kN Máy kéo nén Shimadzu AGS20kNG - Điện áp: ÷ 380 V - Kích thước DxRxC: 80x50x220 cm - Động điêzen gián tiếp - Chiều chiều dài sở: 1845 mm Máy kéo Kubota L4508VN - Công suất thực: 43,6 Hp - Khối lượng tổng: 1370 kg - Vận tốc: ÷ 24,2 km/h - Thiết bị thủy lực: điểm treo, sức nâng điểm nâng 1300 kg Bảng 2P Bảng giá trị tần số biến tần số vòng quay trục TT Trục đĩa gieo Tần số Số vòng quay (Hz) nđ (vòng/phút) Trục trục Tần số Số vòng quay (Hz) nt (vòng/phút) 6,0 10 4,5 6,6 11 5,4 7,2 12 6,3 7,8 13 7,2 8,4 14 8,1 9,0 15 9,0 10 9,6 16 9,9 11 10,2 17 10,8 12 10,8 18 11,7 13 10 11,4 19 12,6 14 11 12,0 20 13,5 15 12 12,6 21 14,4 16 13 13,2 22 15,3 17 14 13,8 23 16,2 18 15 14,4 24 17,1 19 16 15,0 25 18,0 20 17 15,6 26 18,9 21 18 16.2 27 19,8 22 19 16,8 28 20,7 23 20 17,4 29 21,6 24 21 18,0 30 22,5 25 22 18,6 31 23,4 26 23 19,2 32 24,3 27 24 19,8 33 25,2 28 25 20,4 34 26,1 29 26 21,0 35 27,0 30 Bảng 3P Kết thí nghiệm sơ lựa chọn phận bón phân Kiểu trục đặt đáy thùng Thông số Kiểu trục đặt cửa xả Rãnh Rãnh Rãnh cắt lệch Rãnh Rãnh Rãnh cắt lệch 79,4 68,7 81,5 76,5 60,4 72,2 Hệ số nạp đầy (%) 136 59 70 131 52 62 Tỷ lệ tăng khối lượng thể tích (%) 36 22 16 31 14 Năng suất (kg/h) trục - Đường kính ngồi Dt: 110 mm - Chiều dài Lt: 50 mm - Đường kính rãnh dr: 18 mm Thông số cấu tạo - Bề rộng cửa nạp cn: mm trục - Số rãnh Zt: 10 rãnh cửa xả thùng - Số vòng quay trục nt: 12 vịng/phút - Kích thước cửa xả bx.ex: 15x12 mm - Trường hợp: dùng cánh khuấy Bảng 4P Thông số kỹ thuật làm việc LHM TT I 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 1.10 1.11 1.12 1.13 II 2.1 2.2 2.3 III 3.1 3.2 Tên thông số Máy BPĐ Khối lượng: m2 Chiều rộng máy: Bm Chiều dài máy: Lm Chiều cao máy: Hm Số hàng gieo: Zh Số hàng phân: Zp Khoảng cách hàng đầu với mép luống: m Khoảng cách hai hàng kép: bk Đường kính bánh xe: Dbx Tỷ số truyền trung trục đĩa gieo: ig Tỷ số truyền trung trục trục cuốn: it Bộ phận gieo: 02 phận Đường kính ngồi đĩa gieo: Dđ Đường kính vịng trịn chia lỗ đĩa gieo: Dl Chiều dầy đĩa gieo: δđ Đường kính lỗ đĩa: dđ Khe hở đĩa: e Bộ phận bón phân: 02 phận Chiều dài cửa xả phân: bx Chiều rộng cửa xả phân: ex Đường kính ngồi trục cuốn: Dt Chiều dài trục cuốn: Lt Đường kính rãnh trục cuốn: dr Bề rộng cửa nạp phân rãnh trục cuốn: cn Số rãnh trục cuốn: Zr Kích thước luống Bề rộng trung bình mặt luống: B Chiều cao trung bình luống: H Chiều rộng rãnh luống: R Máy kéo Kubota L4508VN: bảng 1P sách hướng dẫn sử dụng Vận tốc làm việc thử nghiệm Đơn vị Giá trị kg mm mm mm hàng hàng mm mm mm 285 1300 1200 1100 75 150 800 0,81 0,79 mm mm mm mm mm 255 240 9,5 mm mm mm mm mm mm rãnh 15 14 110 50 18 10 mm mm mm 800 250 400 m/s 0,6; 0,9 1,3 Bảng 5P Kết thử nghiệm LHM mật độ khoảng cách Vị trí Vị trí Vị trí Mật độ gieo (cây/m2) 142 143 141 47,6 8,5 135 145 144 47,3 8,5 Luống 138 135 146 46,8 8,6 TB 138 141 144 47,2 8,5 Luống 129 133 124 43,1 9,3 137 129 133 44,6 9,0 Luống 135 131 127 43,9 9,2 TB 134 131 128 43,8 9,0 Luống 113 117 121 39,2 10,3 108 106 115 36,7 10,9 Luống 120 119 114 39,4 10,2 TB 114 114 117 38,4 10,5 TT Vận tốc máy (m/s) Tổng số (cây) Luống Luống Luống Luống 0,6 0,9 1,3 Khoảng cách ah (cm) Ghi chú: số liệu lấy sau 10 ngày gieo Bảng 6P Kết thử nghiệm LHM suất máy, sản lượng hệ số nạp phân Thùng (kg) Thùng (kg) Năng suất máy (kg/h) Sản lượng (kg/ha) Hệ số nạp (%) Luống 0,92 0,93 199,8 513,9 58,3 Luống 0,91 0,94 199,8 513,9 58,3 Luống 0,93 0,89 196,6 505,6 57,3 Luống 0,95 0,92 202,0 519,4 58,9 Luống 0,93 0,95 203,0 522,2 59,2 Luống 0,96 0,94 205,2 527,8 59,9 Luống 0,96 0,93 204,1 525,0 59,5 Luống 0,91 0,89 194,4 500,0 56,7 Luống 0,92 0,95 202,0 519,4 58,9 Luống 10 0,89 0,92 195,5 502,8 57,0 200,2 515,0 58,4 TT luống TB Ghi chú: Trường hợp vận tốc máy kéo 0,9 m/s Bảng 7P Một số thông số máy kéo sử dụng nước ta Thông số Đơn vị T-125 Khối lượng kg 6450 Mô men quán tính khối kg.m2 lượng với trục ngang T-25 T-150K K-700 MTZ-5 T-40 1640 3200 2660 10105 160000 260000 43127 26000 5692 11640 Zetop MTSZ MTZ 4011 -50 -80 1851 2800 3200 2360 2370 Chiều dài sở mm 2720 1778 2860 3150 Khoảng cách trọng tâm tới cầu trước mm 950 658 935 1100 Độ cứng cầu trước N/m Độ cứng bánh lốp trước N/m 251000 172000 834000 580000 460000 525000 260000 320000 358000 Độ cứng bánh lốp sau N/m 254000 210000 686000 494000 565000 610000 260000 309000 432000 560000 1500 1525 1612 600000 266000 255000 Nguồn: Lê Minh Lư (2002) Bảng 8P Kết tính vị trí tâm quay tức thời I Tọa độ điểm I Trường hợp Đơn vị I0 I1 I2 I3 I4 I5 I6 I7 I8 xI m 0,656 zI m 0,655 0,654 0,652 0,648 0,655 0,654 0,652 0,650 -0,024 -0,047 -0,070 -0,094 0,024 0,047 0,070 0,094 Sơ đồ vị trí Ghi chú: Tọa độ điểm I tọa độ so với điểm C; I0 treo vị trí nằm ngang; I1, I2, I3, I4, I5, I6, I7 I8: trường hợp tọa độ điểm D theo phương thẳng đứng so với trường hợp treo vị trí nằm ngang là: 0,025, 0,05, 0,075, 0,1, -0,025, -0,05, -0,075 -0,1 m ... gieo hạt đậu tương kết hợp với bón phân làm sở cho việc thiết kế, chế tạo máy gieo đậu tương kết hợp với bón phân Lựa chọn kết cấu chế độ làm việc phù hợp cho liên hợp máy gieo đậu tương kết hợp. .. phận gieo phận bón phân) máy gieo hạt đậu tương kết hợp với bón phân làm sở cho việc thiết kế, chế tạo máy gieo đậu tương kết hợp với bón phân Lựa chọn kết cấu chế độ làm việc phù hợp cho liên hợp. .. hợp máy gieo đậu tương kết hợp với bón phân liên hợp máy làm việc số điều kiện mặt đồng khơng phẳng điển hình 1.3 PHẠM VI NGHIÊN CỨU 1.3.1 Phạm vi đối tượng nghiên cứu Máy gieo hạt đậu tương kết