NGHIÊN CỨU GIẢI PHÁP THIẾT BỊ ĐIỀU KHIỂN TỰ ĐỘNG HỆ THỐNG SAN PHẲNG MẶT ĐỒNG ỨNG DỤNG CÔNG NGHỆ LASER VÀ MÔ PHỎNG BỀ MẶT SAN PHẲNGNGHIÊN CỨU GIẢI PHÁP THIẾT BỊ ĐIỀU KHIỂN TỰ ĐỘNG HỆ THỐNG SAN PHẲNG MẶT ĐỒNG ỨNG DỤNG CÔNG NGHỆ LASER VÀ MÔ PHỎNG BỀ MẶT SAN PHẲNG
BỘ GIÁO DỤC VÀ ĐÀO TẠO ĐẠI HỌC NÔNG LÂM THÀNH PHỐ HỒ CHÍ MINH **************** TRẦN MINH LỘC NGHIÊN CỨU GIẢI PHÁP THIẾT BỊ ĐIỀU KHIỂN TỰ ĐỘNG HỆ THỐNG SAN PHẲNG MẶT ĐỒNG ỨNG DỤNG CÔNG NGHỆ LASER VÀ MÔ PHỎNG BỀ MẶT SAN PHẲNG LUẬN VĂN THẠC SĨ KHOA HỌC KỸ THUẬT Thành Phố Hồ Chí Minh Tháng 02/2011 BỘ GIÁO DỤC VÀ ĐÀO TẠO ĐẠI HỌC NƠNG LÂM THÀNH PHỐ HỒ CHÍ MINH **************** TRẦN MINH LỘC NGHIÊN CỨU GIẢI PHÁP THIẾT BỊ ĐIỀU KHIỂN TỰ ĐỘNG HỆ THỐNG SAN PHẲNG MẶT ĐỒNG ỨNG DỤNG CÔNG NGHỆ LASER VÀ MÔ PHỎNG BỀ MẶT SAN PHẲNG Chuyên ngành : Kỹ thuật Cơ khí Mã số : 60 52 14 LUẬN VĂN THẠC SĨ KHOA HỌC KỸ THUẬT Hướng dẫn khoa học TS NGUYỄN VĂN HÙNG Thành Phố Hồ Chí Minh Tháng 02/2011 NGHIÊN CỨU GIẢI PHÁP THIẾT BỊ ĐIỀU KHIỂN TỰ ĐỘNG HỆ THỐNG SAN PHẲNG MẶT ĐỒNG ỨNG DỤNG CÔNG NGHỆ LASER VÀ MÔ PHỎNG BỀ MẶT SAN PHẲNG TRẦN MINH LỘC Hội đồng chấm luận văn: TS NGUYỄN NHƯ NAM Chủ tịch: Đại học Nông Lâm TP HCM TS LÊ ANH ĐỨC Thư ký: Đại học Nông Lâm TP HCM Phản biện 1: TS PHAN HIẾU HIỀN Đại học Nông Lâm TP HCM Phản biện 2: PGS TS TRẦN THỊ THANH Đại học Nông Lâm TP HCM Uỷ viên: TS NGUYỄN VĂN HÙNG Đại học Nông Lâm TP HCM TRƯỜNG ĐẠI HỌC NƠNG LÂM TP HỒ CHÍ MINH HIỆU TRƯỞNG i LÝ LỊCH CÁ NHÂN Tôi tên Trần Minh Lộc, sinh ngày 20 tháng 09 năm 1978, huyện Đức Phổ, tỉnh Quảng Ngãi Con Ông Trần Minh Bảo Bà Nguyễn Thị Bờ Tốt nghiệp tú tài trường Trung học phổ thông Đức Phổ I, tỉnh Quảng Ngãi năm 1998 Tốt nghiệp Đại học năm 2003, ngành Cơ khí chế tạo máy, hệ dài hạn tập trung, trường Đại học Sư Phạm Kỹ Thuật, thành phố Hồ Chí Minh Sau làm việc liên tục đơn vị sau: Tháng năm 2003 đến 8/2005 làm việc công ty TNHHTM&DV VPS, chức vụ kỹ sư thiết kế Tháng 9/2005 đến nay, Giảng viên khoa Cơ Khí Trường Cao Đẳng Cơng Nghệ Thủ Đức, thành phố Hồ Chí Minh Tháng năm 2008 theo học Cao học ngành Cơ khí Đại học Nơng Lâm, Thủ Đức, thành phố Hồ Chí Minh Tình trạng gia đình: Độc thân Địa liên lạc: 103/9A đường Dương Đình Hội, khu phố 5, phường Phước Long B quận 9, thành phố Hồ Chí Minh Điện thoại DĐ: 0983.383.337 Cơ quan: Số 53 Võ Văn Ngân, phường Linh Chiểu, quận Thủ Đức, thành phố Hồ Chí Minh Điện thoại: 08.22158645 Email: minhloc09@gmail.com loc.tranminh@gmail.com ii LỜI CAM ĐOAN Tôi cam đoan cơng trình nghiên cứu tơi Các số liệu kết qủa luận văn trung thực chưa công bố cơng trình khác Tác giả luận văn Trần Minh Lộc iii CẢM TẠ Để hoàn thành luận văn thạc sĩ tơi xin chân thành bày tỏ lòng kính trọng biết ơn sâu sắc đến: Tiến sĩ Nguyễn Văn Hùng, Trưởng môn Cơ điện tử Trường Đại học Nơng lâm thành phố Hồ Chí Minh – Giáo viên hướng dẫn đề tài Thầy tận tình bảo, giúp đỡ tạo điều kiện thuận lợi cho tơi suốt q trình làm đề tài Qua thời gian làm việc thầy, tơi có kiến thức nghiên cứu khoa học, cách nhận định đánh giá vấn đề,… Đó tảng cho tiếp bước vững công tác giảng dạy sau Ban giám hiệu, Phòng sau đại học, Cô PGS.TS Trần Thị Thanh, Chủ nhiệm Khoa tập thể Giảng viên Khoa Cơ khí - Cơng nghệ, Trường Đại học Nơng Lâm thành phố Hồ Chí Minh tận tình giảng dạy, truyền đạt kiến thức quí báu tạo điều kiện tốt trình học Cao học thực luận văn Ban lãnh đạo quý thầy cô Trung tâm Năng Lượng Máy Nông Nghiệp (CAEM), Trường Đại Học Nơng Lâm Tp.Hồ Chí Minh cung cấp thông tin liên quan tới đề tài, hướng dẫn, tạo điều kiện để đề tài triển khai Ban giám đốc công ty Cổ Phần Nông Nghiệp Lý Tưởng hỗ trợ kinh phí cho đề tài Xin cảm ơn quý thầy, cô phản biện đề tài cho lời nhận xét qúi báu Qua phản hồi tơi để hồn thiện q trình nghiên cứu Các anh chị học viên lớp cao học khóa 2007 2008 hỗ trợ, giúp đỡ, động viên chia sẻ kinh nghiệm làm việc Cảm ơn cha mẹ, anh em, bạn bè đồng nghiệp tạo điều kiện tốt suốt trình học tập làm đề tài iv TÓM TẮT Đề tài “Nghiên cứu giải pháp thiết bị điều khiển tự động hệ thống san phẳng mặt đồng ứng dụng công nghệ laser mô bề mặt san phẳng ” tiến hành Khoa Cơ khí – Cơng nghệ Trung tâm Năng lượng Máy Nông nghiệp Trường Đại học Nơng lâm Thành phố Hồ Chí Minh, thời gian từ tháng 03 năm 2010 đến tháng 12 năm 2010 với kết tóm tắt sau: Nghiên cứu, cải tiến cụm van thủy lực từ hệ thống ngoại nhập, khắc phụ tượng nóng dầu trình hoạt động hệ thống Nghiên cứu chế tạo thành công điều khiển PLC thay cho điều khiển ngoại nhập, góp phần nội địa hóa thiết bị điều khiển Nghiên cứu, thiết kế, chế tạo giảm chấn cấu nâng hạ gàu san nhằm khắc phục tượng rung động nâng bánh xe cụm gàu san lên khỏi mặt đồng gàu san vùng đất cao Một chương trình mơ bề mặt đồng viết phần mềm Matlab 7.5 để mơ tính tốn cho q trình san phẳng Kết tính tốn mơ thể giao diện trực quan Chương trình giúp xác định độ nhấp nhơ ban đầu khu đất, dễ dàng tìm đường chạy hợp lý cho máy kéo đồng thời tính tốn thời gian ước lượng chi phí cho q trình san phẳng v ABSTRACT The thesis "Study on automatic controlling solutions of a laser land leveling system and simulation the field surface applied leveling" was done at the Faculty of Engineering and the Center for Energy and Agriculture Machinery, Nong Lam University of Ho Chi Minh City, from March to December 2010 with the summary results as following: Improved the hydraulic system, resolved phenomenon the hot oil happening during the operational process Designed and manufactured the PLC controller in order to replace the imported one and reduce the equipment cost Designed and manufactured the damping and lifting device for the scraper to reduce the “shock” phenomenon A program was written based on the Matlab 7.5 software to simulate and calculate for the operational process Results of calculation and simulation are represented by visual interface This program is very helpful for determining the first undulating of land, finding a logical route of tractors and calculating leveling cost vi MỤC LỤC TRANG Trang chuẩn Y i Lý lịch cá nhân ii Lời cam đoan iii Cảm tạ iv Tóm tắt v Tóm tắt tiếng Anh vi Mục lục vii Bảng ký hiệu chữ viết tắt xi Danh sách liệt kê ký hiệu .xii Danh sách hình xiv Danh sách bảng xvii MỞ ĐẦU 1.1 Đặt vấn đề 1.2 Mục tiêu đề tài 1.3 Phạm vi nghiên cứu TỔNG QUAN 2.1 Đối tượng nghiên cứu 2.1.1 Cấu tạo nguyên lý hoạt động hệ thống san phẳng laser 2.1.2 Bộ phát laser 2.1.3 Bộ nhận laser 2.1.4 Hộp điều khiển 2.1.5 Hệ thống thủy lực 2.1.6 Gàu san… .11 2.2 Cơ sở tính tốn thiết kế hệ thống thủy lực 12 2.2.1 Cơ sở tính tốn bơm 12 2.2.1.1 Các đại lượng đặc trưng bơm 12 vii 2.2.1.2 Tiêu chuẩn chọn bơm 15 2.2.2 Các loại van thủy lực ứng dụng mạch điều khiển tự động 15 2.2.2.1 Van đảo chiều 15 2.2.2.2 Van solenoid ứng dụng điều khiển tự động 17 2.2.2.3 Van servo 18 2.2.3 Cơ sở tính tốn thiết kế xi lanh thủy lực .20 2.3 Một số nghiên cứu hệ thống san phẳng laser 23 2.3.1 Tại Việt Nam 23 2.3.2 Một số nghiên cứu hệ thống san phẳng laser giới 24 NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU .29 3.1 Nội dung nghiên cứu 29 3.2 Phương pháp nghiên cứu .29 3.2.1 Phương pháp nghiên cứu lý thuyết 29 3.2.1.1 Phương pháp tiếp cận hệ thống 29 3.2.1.2 Phương pháp thiết kế hệ thống thủy lực 30 3.2.1.3 Phương pháp thiết kế phận điều khiển tự động 31 3.2.2 Phương pháp nghiên cứu thực nghiệm 31 3.2.2.1 Phạm vi khảo nghiệm thiết bị thí nghiệm 31 3.2.2.2 Thiết bị dụng cụ đo 32 3.2.2.3 Phương pháp đo 32 3.2.3 Phương pháp đo mô bề mặt đồng 33 3.2.3.1 Phương pháp đo bề mặt đồng 33 3.2.3.2 Phương pháp mô bề mặt đồng 34 3.2.4 Phương pháp nội suy Spline bậc hai tìm hàm mơ bề mặt đồng .35 3.2.5 Phương pháp tính tốn thể tích đất vùng 37 3.2.6 Phương pháp khảo nghiệm đánh giá phận chế tạo 39 KẾT QUẢ VÀ THẢO LUẬN 40 4.1 Giải pháp mơ bề mặt đồng tính tốn thơng số q trình san phẳng 40 viii Bảng PL3.4 Chiều cao đo thực tế lô thứ hai 2033 2020 2018 1998 1978 1997 2016 2017 2018 2033 2012 1949 1966 1975 1971 1995 1965 1982 1999 2014 2007 1981 1965 1975 1979 1995 1968 1983 1999 2001 2001 2013 1963 1975 1987 1995 1971 1984 1998 1988 2001 2004 2005 1983 1962 1989 1985 1986 1987 1914 1986 1981 1993 1981 1970 1980 1996 1987 1978 1939 1972 1959 1981 1979 1978 1972 2007 1988 1969 1965 1976 1969 1983 1983 1983 1958 1982 1972 1963 1979 1953 1978 1974 1974 1974 1959 1981 1974 1966 1979 1929 1987 1964 1964 1965 1960 1980 1975 1970 1980 1932 1982 1980 1970 1964 1990 1985 1985 1984 1982 1943 1958 1985 1965 1952 1975 1969 1965 1962 1984 1954 1934 1990 1960 1940 1960 1954 1945 1940 1987 1954 1934 1990 1960 1940 1960 1954 1945 1940 1987 Bảng PL3.5 Chiều cao nhấp nhô so với mặt phẳng trung bình lơ thứ hai 56.3 43.3 41.3 21.3 1.3 20.3 39.3 40.3 41.3 56.3 35.3 -27.7 -10.7 -1.7 -5.7 18.3 -11.7 5.3 22.3 37.3 30.3 4.3 -11.7 -1.7 2.3 18.3 -8.7 6.3 22.3 24.3 24.3 36.3 -13.7 -1.7 10.3 18.3 -5.7 7.3 21.3 11.3 24.3 27.3 28.3 6.3 -14.7 12.3 8.3 9.3 10.3 -62.7 9.3 4.3 16.3 4.3 -6.7 3.3 19.3 10.3 1.3 -37.7 -4.7 -17.7 4.3 2.3 1.3 -4.7 30.3 11.3 -7.7 -11.7 -0.7 -7.7 6.3 6.3 6.3 -18.7 5.3 -4.7 -13.7 2.3 -23.7 1.3 -2.7 -2.7 -2.7 -17.7 4.3 -2.7 -10.7 2.3 -47.7 10.3 -12.7 -12.7 -11.7 -16.7 3.3 -1.7 -6.7 3.3 -44.7 5.3 3.3 -6.7 -12.7 13.3 8.3 8.3 7.3 5.3 -33.7 -18.7 8.3 -11.7 -24.7 -1.7 -7.7 -11.7 -14.7 7.3 -22.7 -42.7 13.3 -16.7 -36.7 -16.7 -22.7 -31.7 -36.7 10.3 -22.7 -42.7 13.3 -16.7 -36.7 -16.7 -22.7 -31.7 -36.7 10.3 100 Bảng PL3.6 Chiều cao đo thực tế lô thứ ba 2047 2047 2040 2057 2012 2012 2023 2028 2024 2033 2054 2054 2050 2051 2015 2015 2031 2034 2027 2012 2054 2054 2053 2051 2034 2034 2029 2033 2023 2007 2053 2053 2055 2051 2052 2052 2027 2031 2018 2001 2044 2044 2072 2052 2076 2076 1998 1991 1999 2001 2030 2030 2061 2062 2059 2059 2002 1979 1988 1986 2016 2016 2050 2073 2042 2042 2006 1968 1977 1972 2024 2024 2034 2016 2028 2028 2042 1974 1963 1976 2019 2019 2028 2017 2023 2023 2028 1952 1969 1953 2014 2014 2022 2018 2018 2018 2014 1930 1974 1929 2038 2038 2052 2036 2082 2082 2008 1982 1948 1932 2025 2025 2031 2018 2024 2024 1984 1967 1939 1943 2012 2012 2010 2000 1966 1966 1960 1953 1930 1954 2012 2012 2010 2000 1966 1966 1960 1953 1930 1954 Bảng PL3.7 Chiều cao nhấp nhô so với mặt phẳng trung bình lơ thứ ba 32.7 32.7 25.7 42.7 -2.3 -2.3 8.7 13.7 9.7 18.7 39.7 39.7 35.7 36.7 0.7 0.7 16.7 19.7 12.7 -2.3 39.7 39.7 38.7 36.7 19.7 19.7 14.7 18.7 8.7 -7.3 38.7 38.7 40.7 36.7 37.7 37.7 12.7 16.7 3.7 -13.3 29.7 29.7 57.7 37.7 61.7 61.7 -16.3 -23.3 -15.3 -13.3 15.7 15.7 46.7 47.7 44.7 44.7 -12.3 -35.3 -26.3 -28.3 1.7 1.7 35.7 58.7 27.7 27.7 -8.3 -46.3 -37.3 -42.3 9.7 9.7 19.7 1.7 13.7 13.7 27.7 -40.3 -51.3 -38.3 4.7 4.7 13.7 2.7 8.7 8.7 13.7 -62.3 -45.3 -61.3 -0.3 -0.3 7.7 3.7 3.7 3.7 -0.3 -84.3 -40.3 -85.3 23.7 23.7 37.7 21.7 67.7 67.7 -6.3 -32.3 -66.3 -82.3 10.7 10.7 16.7 3.7 9.7 9.7 -30.3 -47.3 -75.3 -71.3 -2.3 -2.3 -4.3 -14.3 -48.3 -48.3 -54.3 -61.3 -84.3 -60.3 -2.3 -2.3 -4.3 -14.3 -48.3 -48.3 -54.3 -61.3 -84.3 -60.3 101 Bảng PL3.8 Chiều cao đo thực tế lô thứ tư 2135 2140 2145 2140 2115 2116 2118 2119 2115 2088 2137 2137 2146 2140 2114 2114 2111 2118 2111 2087 2129 2129 2146 2127 2145 2145 2119 2117 2107 2098 2121 2121 2147 2115 2176 2176 2128 2116 2104 2109 2088 2088 2032 2110 2110 2110 2110 2006 2018 2108 2102 2102 2027 2115 2101 2101 2091 2043 2052 2097 2116 2116 2022 2120 2092 2092 2073 2081 2087 2087 2118 2118 2098 2170 2188 2188 2089 2078 2081 2080 2115 2115 2107 2150 2184 2184 2088 2063 2072 2061 2111 2111 2116 2131 2180 2180 2086 2048 2063 2042 2140 2140 2128 2111 2115 2115 2025 2038 2046 2048 2100 2100 2097 2095 2086 2086 2035 2041 2037 2045 2061 2061 2066 2080 2058 2058 2046 2044 2028 2041 2047 2047 2040 2057 2012 2012 2023 2028 2024 2033 Bảng PL3.9 Chiều cao nhấp nhơ so với mặt phẳng trung bình lô thứ tư 39.1 44.1 49.1 44.1 19.1 20.1 22.1 23.1 19.1 -7.9 41.1 41.1 50.1 44.1 18.1 18.1 15.1 22.1 15.1 -8.9 33.1 33.1 50.1 31.1 49.1 49.1 23.1 21.1 11.1 2.1 25.1 25.1 51.1 19.1 80.1 80.1 32.1 20.1 8.1 13.1 -7.9 -7.9 -63.9 14.1 14.1 14.1 14.1 -89.9 -77.9 12.1 6.1 6.1 -68.9 19.1 5.1 5.1 -4.9 -52.9 -43.9 1.1 20.1 20.1 -73.9 24.1 -3.9 -3.9 -22.9 -14.9 -8.9 -8.9 22.1 22.1 2.1 74.1 92.1 92.1 -6.9 -17.9 -14.9 -15.9 19.1 19.1 11.1 54.1 88.1 88.1 -7.9 -32.9 -23.9 -34.9 15.1 15.1 20.1 35.1 84.1 84.1 -9.9 -47.9 -32.9 -53.9 44.1 44.1 32.1 15.1 19.1 19.1 -70.9 -57.9 -49.9 -47.9 4.1 4.1 1.1 -0.9 -9.9 -9.9 -60.9 -54.9 -58.9 -50.9 -34.9 -34.9 -29.9 -15.9 -37.9 -37.9 -49.9 -51.9 -67.9 -54.9 -48.9 -48.9 -55.9 -38.9 -83.9 -83.9 -72.9 -67.9 -71.9 -62.9 102 Phụ lục Kết khảo nghiệm đánh giá phận thiết kế Ngày khảo nghiệm: 15/11/2010 1/12/2010 Địa điểm: Lô đất sau khoa Cơ Khí, Trường ĐH Nơng Lâm TPHCM, phường Linh Trung, quận Thủ Đức 4.1 Khảo nghiệm xác định nhiệt độ dầu hệ thống thủy lực trình hoạt động Ngày khảo nghiệm: 15/11/2010 Người khảo nghiệm: Trần Minh Lộc, Nguyễn Đức Cảnh Bảng PL4.1 Nhiệt độ ống dầu hệ thống thủy lực Thời gian đo Giờ qui đổi ngày (phút) cụm van 3/2 (0C) cụm van 4/3 (0C) 14 45 35 35 14 55 10 46 38 15 05 20 54 42 15 15 30 65 44 15 25 40 67 45 15 35 50 68 46 16 45 60 70 46 15 55 70 71 46 TT Nhiệt độ ống dầu Nhiệt độ ống dầu 4.2 Khảo nghiệm hệ thống điều khiển PLC Ngày khảo nghiệm: 15/11/2010 Người khảo nghiệm: Trần Minh Lộc, Nguyễn Đức Cảnh 103 Bảng PL4.2 Kết đo thời gian nâng hạ hệ thống điều khiển Bộ điều khiển ngoại nhập Bộ điều khiển PLC Số Độ Thời gian Thời gian hạ từ Thời gian nâng Thời gian hạ từ lần cao nâng từ mức mức trung bình từ mức trung mức trung bình đo (mm) trung bình lên xuống mức bình lên mức xuống mức thấp mức cao(giây) thấp (giây) cao (giây) (giây) 1,08 Lần 50 0,68 Lần 50 0,58 Lần 50 0,67 Lần 80 0,98 Lần 80 0,99 Lần 80 1,01 TB 1,01 0,64 0.96 1,36 TB 1,24 0,99 TB 1,01 TB 1,29 1,27 0,67 0,71 0,69 TB 0,69 0,99 1,05 TB 1,02 1,04 1,15 1,12 1,08 1,38 1,33 TB 1,11 TB 1,35 1,35 4.3 Xác định áp suất dầu trình nâng hạ gàu san Ngày khảo nghiệm: 1/12/2010 Người khảo nghiệm: Trần Minh Lộc, Nguyễn Đức Cảnh Bảng PL4.3 Áp suất nâng hạ gàu san chế độ khơng tải khơng có giảm chấn Thời gian giây giây giây giây Lần 10 12,5 15 25 Lần 12,5 15 17,5 20 Lần 12,5 15 20 22,5 Lần 10 12,5 15 17,5 Lần 10 12,5 15 22,5 Áp suất trung bình 11 13,5 16,5 21,5 Áp suất(bar) 104 Bảng PL4.4 Áp suất nâng hạ gàu san chế độ khơng tải có giảm chấn Thời gian giây giây giây giây Lần 12,5 15 17,5 22,5 Lần 10 15 17,5 22,5 Lần 7,5 12,5 15 20 Lần 10 10 15,5 20 Lần 12,5 12,5 15 22,5 Áp suất trung bình 10,5 13 16,1 21,5 Áp suất (bar) Bảng PL4.5 Áp suất nâng hạ gàu san chế độ có tải khơng có giảm chấn Thời gian giây giây giây giây 10 giây Lần 22,5 25 27,5 30 35 Lần 25 25 30 40 40 Lần 25 27,5 35 37,5 40 Lần 22,5 25 27,5 30 37,5 Lần 22,5 25 30 35 35 Áp suất trung bình 23,5 25,5 30 34,5 37,5 Áp suất(bar) Bảng PL4.6 Áp suất nâng hạ gàu san chế độ có tải có giảm chấn Thời gian Áp suất(bar) Lần giây giây giây giây 10 giây 20 22,5 25 25 30 Lần 20 25 25 27,5 32,5 Lần 22,5 27,5 30 30 35 Lần 20 25 27,5 37,5 40 Lần 17,5 20 22,5 24,5 30 27,5 35 31 27,5 33 Áp suất trung bình 105 Phụ lục Chương trình Matlab mơ mặt đồng PLC LOGO 5.1 Chương trình Matlab function varargout = mpmatdong(varargin) gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, 'gui_Singleton', gui_Singleton, 'gui_OpeningFcn', @mpmatdong_OpeningFcn, 'gui_OutputFcn', @mpmatdong_OutputFcn, 'gui_LayoutFcn', [] , 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end function mpmatdong_OpeningFcn(hObject, eventdata, handles, varargin) handles.output = hObject; guidata(hObject, handles); function varargout = mpmatdong_OutputFcn(hObject, eventdata, handles) varargout{1} = handles.output; function edit1_Callback(hObject, eventdata, handles) user_entry = str2double(get(hObject,'string')); if isnan(user_entry) errordlg('You must enter a numeric value','Bad Input','modal') end function edit1_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end Tương tự cho edit lại …………………………… function edit100_Callback(hObject, eventdata, handles) 106 user_entry = str2double(get(hObject,'string')); if isnan(user_entry) errordlg('You must enter a numeric value','Bad Input','modal') end set(handles.checkbox3,'Value',0); function pushbutton1_KetQuaDoThi_Callback(hObject, eventdata, handles) rotate3d off; cla reset; val1=get(handles.edit1,'String');val1=str2num(val1);h(1,1)=val1; …………………………………… val100=get(handles.edit100,'String');val100=str2num(val100);h(10,10)=val100; a=get(handles.radiobutton1,'Value'); b=get(handles.radiobutton2,'Value'); c=get(handles.radiobutton3,'Value'); d=get(handles.radiobutton4,'Value'); e=get(handles.radiobutton5,'Value'); f=get(handles.checkbox1,'Value'); g=get(handles.checkbox2,'Value'); i=get(handles.checkbox3,'Value'); htb=mean(mean(h)); Ht=(h-htb)*(-1); [x,y]=meshgrid(1:10,1:10); [x1,y1]=meshgrid(1:0.18:10); % tao mang xap xi cho matran vuong set(handles.axes1,'Visible','off'); set(handles.axes2,'Visible','off'); xlabel('X-axis, m'); ylabel('Y-axis, m'); zlabel('do cao mat dong, m'); title('Mo phong mat dong'); if a==1 %Ve thi luoi set(handles.axes2,'Visible','off'); set(handles.axes1,'Visible','on'); axes(handles.axes1); mesh(Ht); az = -45; el = 60; view(az, el); end if b==1 %Ve thi cot set(handles.axes2,'Visible','off'); set(handles.axes1,'Visible','on'); axes(handles.axes1); bar3(Ht); end if c==1 %Ve thi noi suy bac ba set(handles.axes2,'Visible','off'); set(handles.axes1,'Visible','on'); axes(handles.axes1); 107 h1=interp2(x,y,Ht,x1,y1,'spline'); meshz(h1); az = -45; el = 60; view(az, el); colorbar('ZTickLabel', {'Freezing','Cold','Cool','Neutral', 'Warm','Hot','Burning','Nuclear'}) end if d==1 %Ve thi noi suy tuyen tinh set(handles.axes2,'Visible','off'); set(handles.axes1,'Visible','on'); axes(handles.axes1); h1=interp2(x,y,Ht,x1,y1,'linear'); mesh(h1); az = -45; el = 60; view(az, el); end if e==1 %Ve thi noi suy diem gan nhat set(handles.axes2,'Visible','off'); set(handles.axes1,'Visible','on'); axes(handles.axes1); h1=interp2(x,y,Ht,x1,y1,'nearest'); mesh(h1); az = -45; el = 60; view(az, el); end if g==1 % Ve thi noi suy bac ba phong to set(handles.uipanel1,'Visible','off'); set(handles.uipanel2,'Visible','off'); set(handles.uipanel3,'Visible','off'); set(handles.axes1,'Visible','off'); set(handles.axes2,'Visible','on'); axes(handles.axes2); H1=interp2(x,y,Ht,x1,y1,'spline'); meshz(H1); colorbar('ZTickLabel', {'Freezing','Cold','Cool','Neutral', 'Warm','Hot','Burning','Nuclear'}) az = -45; el = 60; view(az, el); end if i==1 % Hien trang mat bang khu dat set(handles.uipanel1,'Visible','off'); set(handles.uipanel2,'Visible','off'); set(handles.uipanel3,'Visible','off'); set(handles.axes1,'Visible','off'); set(handles.axes2,'Visible','on'); axes(handles.axes2); 108 [x2,y2]=meshgrid(1:0.18:10); H1=interp2(x,y,Ht,x1,y1,'spline');%ham tinh tich phan H2=interp2(x,y,Ht,x2,y2,'spline');% Ham ve noi suy duong dong muc %meshz(H1); syms x y real% khai bao dang toan hoc cho bien x va y tp=(double(int(int(H1,x,0,1),y,0,1))/1000)+0.0005; lg=(tp>0); tp1=lg.*tp; V=sum(sum(tp1)) b1=[ones(25) zeros(25,26);zeros(26,25) zeros(26)].*tp; b2=[zeros(25) zeros(25,26);ones(26,25) zeros(26)].*tp; b3=[zeros(25) zeros(25,26);zeros(26,25) ones(26)].*tp; b4=[zeros(25) ones(25,26);zeros(26,25) zeros(26)].*tp; V1=sum(sum(b1)); V2=sum(sum(b2)) V3=sum(sum(b3)); V4=sum(sum(b4)); kt=sum(sum(tp)) [C,h] = contourf(x2,y2,H2,9); clabel(C,h,'FontSize',13,'Color','k','Rotation',0); colorbar('ZTickLabel', {'Freezing','Cold','Cool','Neutral', 'Warm','Hot','Burning','Nuclear'}) hold on x=[1 10]; y=[5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 ]; plot(x,y,'color','m','LineWidth',3) hold on x=[5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 ]; y=[1 10]; plot(x,y,'color','m','LineWidth',3) xlabel('truc x (met)'); ylabel('truc y (met)'); text(2.5,3,'Vung 1','FontSize',18,'Color','r'); text(3,2.5, ['V1 = ',num2str(V1)],'FontSize',18, 'HorizontalAlignment','center', 'BackgroundColor',[.7 7]); ………………………… text(7.5,3,'Vung 4','FontSize',18,'Color','r'); text(8,2.5, ['V4 = ',num2str(V4)],'FontSize',18, 'HorizontalAlignment','center', 'BackgroundColor',[.7 7]); 109 end if f==1 rotate3d on; end function pushbutton2_Close_Callback(hObject, eventdata, handles) close all; function pushbutton4_Callback(hObject, eventdata, handles) set(handles.uipanel1,'Visible','on'); set(handles.uipanel2,'Visible','on'); set(handles.uipanel3,'Visible','on'); set(handles.axes2,'Visible','off'); colorbar off; function checkbox3_Callback(hObject, eventdata, handles) set(handles.radiobutton1,'Value',0); set(handles.radiobutton2,'Value',0); set(handles.radiobutton3,'Value',0); set(handles.radiobutton4,'Value',0); set(handles.radiobutton5,'Value',0); set(handles.checkbox2,'Value',0); set(handles.checkbox3,'Value',1); function pushbutton3_Ketquatinhtoan_Callback(hObject, eventdata, handles) clc val1=get(handles.edit1,'String');val1=str2num(val1);h(1,1)=val1; val2=get(handles.edit2,'String');val2=str2num(val2);h(1,2)=val2; Tương tự cho edit lại …………………………………… val100=get(handles.edit100,'String');val100=str2num(val100);h(10,10)=val100; val116=get(handles.edit116,'String');val116=str2num(val116);a=val116 val117=get(handles.edit117,'String');val117=str2num(val117);b=val117 Vt=0.9; v=4; bgs=2; htb=mean(mean(h)) % Tinh gia tri trung binh Ht=h-htb [x,y]=meshgrid(1:10,1:10); [x1,y1]=meshgrid(1:0.18:10); H1=interp2(x,y,Ht,x1,y1,'spline'); syms x y real % khai bao dang toan hoc cho bien x va y tp=(double(int(int(H1,x,0,1),y,0,1))/1000)-0.0005; lg=(tp>0); tp1=lg.*tp; V=sum(sum(tp1)) % Tim the tich dat can phai san phang b1=[ones(25) zeros(25,26);zeros(26,25) zeros(26)].*tp; b2=[zeros(25) zeros(25,26);ones(26,25) zeros(26)].*tp; 110 b3=[zeros(25) zeros(25,26);zeros(26,25) ones(26)].*tp; b4=[zeros(25) ones(25,26);zeros(26,25) zeros(26)].*tp; V1=sum(sum(b1)); V2=sum(sum(b2)) V3=sum(sum(b3)); V4=sum(sum(b4)); kt=sum(sum(tp)); % Tinh thoi gian san tho if ((V10)&&(V3>0)&&(V4>0))||((V1>0)&&(V20))||((V2>0)&&(V10))||((V3>0)&&(V10))||((V4>0)&&(V10)&&(V30)&&(V4>0)&&(V20)&&(V3>0)&&(V2