Nghiên cứu xác định chế độ làm việc hợp lý của liên hợp máy thu hoạch khoai tây

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - ĐỖ ĐÌNH THI NGHIÊN CỨU XÁC ĐỊNH CHẾ ĐỘ LÀM VIỆC HỢP LÝ CỦA LIÊN HỢP MÁY THU HOẠCH KHOAI TÂY LUẬN VĂN THẠC SĨ KỸ THUẬT Hà Nội - Năm 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - ĐỖ ĐÌNH THI NGHIÊN CỨU XÁC ĐỊNH CHẾ ĐỘ LÀM VIỆC HỢP LÝ CỦA LIÊN HỢP MÁY THU HOẠCH KHOAI TÂY CHUYÊN NGÀNH KỸ THUẬT CƠ KHÍ ĐỘNG LỰC LUẬN VĂN THẠC SĨ KỸ THUẬT NGƢỜI HƢỚNG DẪN 1: TS ĐÀM HOÀNG PHÚC NGƢỜI HƢỚNG DẪN 2: TS BÙI VIỆT ĐỨC Hà Nội – Năm 2016 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu tơi Những nội dung đƣợc trình bày luận văn thực với hƣớng dẫn khoa học thầy giáo TS Đàm Hoàng Phúc – giảng viên Bộ mơn tơ – Viện Cơ khí động lực – Trƣờng Đại học Bách Khoa Hà Nội thầy giáo TS Bùi Việt Đức – giảng viên Bộ môn ô tô – Khoa Cơ điện- Học Viện Nông nghiệp Việt Nam Toàn nội dung luận văn hoàn toàn phù hợp với nội dung đƣợc đăng ký phê duyệt Hiệu trƣởng Trƣờng Đại học Bách Khoa Hà Nội Các số liệu, kết luận văn trung thực Hà Nội, ngày 29 tháng 03 năm 2016 Tác giả Đỗ Đình Thi i LỜI CẢM ƠN Có đƣợc kết nghiên cứu hồn thành luận văn Thạc sĩ mình, ngồi nỗ lực thân, nhận đƣợc giúp đỡ nhiệt tình nhiều cá nhân, tổ chức nhà khoa học Trƣớc tiên xin bày tỏ lòng biết ơn chân thành sâu sắc đến thầy giáo TS Đàm Hoàng Phúc – giảng viên Bộ mơn tơ –Viện Cơ khí động lực – Trƣờng Đại học Bách Khoa Hà Nội thầy giáo TS Bùi Việt Đức – giảng viên Bộ môn ô tô – Khoa Cơ điện- Học Viện Nông nghiệp Việt Nam, ngƣời tận tình bảo, truyền đạt cho tơi kiến thức kinh nghiệm suốt trình thực nghiên cứu khoa học Nhân dịp này, xin chân thành cảm ơn thầy cô giáo Viện Cơ khí động lực – Trường Đại học Bách Khoa Hà Nội thầy cô Khoa Cơ điện- Học Viện Nơng nghiệp Việt Nam,, người có ý kiến quý báu tạo điều kiện giúp đỡ tơi hồn thành luận văn Tơi xin chân thành cảm ơn Ban Giám hiệu trường Đại học Bách Khoa Hà Nội, Viện Đào tạo sau Đại học, lãnh đạo cán Viện Cơ khí động lực tạo điều kiện giúp đỡ trình nghiên cứu khoa học Tơi xin chân thành cảm ơn người thân gia đình, bạn bè, đồng nghiệp động viên, khích lệ tơi suốt q trình học tập, nghiên cứu hồn thành luận văn Xin trân trọng cảm ơn! Hà Nội, tháng năm 2016 Tác giả luận văn Đỗ Đình Thi ii 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 vi DANH MỤC HÌNH vii ĐẶT VẤN ĐỀ 1 Sự cấp thiết vấn đề nghiên cứu Mục tiêu đề tài CHUƠNG 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU .2 1.1 Đối tƣợng nghiên cứu 1.1.1 Tình hình sản xuất khoai tây nƣớc ta 1.1.2 Đặc tính sinh trƣởng phát triển yêu cầu nông học khoai tây 1.1.3 Các phƣơng pháp thu hoạch khoai tây nƣớc ta 1.1.4 Tình hình giới hóa thu hoạch khoai tây ngồi nƣớc 1.1.5 Tình hình nghiên cứu ứng dụng giới hóa thu hoạch khoai tây CHƢƠNG 2: CƠ SỞ LÝ THUYẾT TÍNH TỐN THÀNH LẬP LIÊN HỢP MÁY 2.1 Tính chất lý đất 2.1.1 Khả chống nén đất .8 2.1.2 Khả chống cắt đất 10 2.2 Máy kéo bánh lốp .17 2.3 Các yếu tố ảnh hƣởng đến hiệu sử dụng LHM thu hoạch khoai tây 19 2.3.1 Động học máy kéo bánh lốp 19 2.3.2 Bánh xe lăn không trƣợt (lăn túy) 19 2.3.3 Các tƣợng trƣợt bánh xe 21 2.4 Các tính chất bánh (bánh lốp có săm) 24 2.5 Bán kính bánh xe 26 2.6 Cân công suất hiệu suất 28 2.7 Ảnh hƣởng lực cản kéo đến độ trƣợt 30 iii 2.8 Xác định phản lực pháp tuyến mặt đƣờng tác dụng lên bánh xe máy kéo .33 2.9 Tính kéo máy kéo .37 2.10 Đồ thị cân công suất đƣờng đặc tính kéo 39 2.10.1 Đồ thị cân công suất 39 2.10.2 Đƣờng đặc tính kéo .41 2.11 Đƣờng đặc tính kéo máy kéo dùng hộp số học 42 2.11.1 Khái niệm chung đƣờng đặc tính kéo dùng hốp số học .42 2.12 Máy thu hoạch khoai tây 46 2.12.1 Tính tốn thơng số máy đào khoai tây 47 2.12.2 Lực cản máy thu hoạch khoai tây 50 2.13 Các thông số chế độ làm việc tối ƣu Liên hợp máy .53 2.13.1 Xác định bề rộng làm việc tối ƣu Btƣ 54 2.13.2 Xác định vận tốc làm việc tối ƣu Vtu .55 CHƢƠNG III XÂY DỰNG MƠ HÌNH TÍNH TỐN CHO LHM THU HOẠCH 58 3.1 Các giả thiết để xây dựng mơ hình 58 3.2 Mơ hình hóa liên hợp máy 58 3.2.1 Phần tử động .59 3.2.2 Phần tử truyền lực (TL) .61 3.2.3 Phần tử hệ thống di động (DĐ) 61 3.2.4 Phần tử máy thu hoạch khoai tây 62 3.2.5 Phần tử điều kiện sử dụng .62 3.3 Các hàm mục tiêu .62 3.4 Mơ hình tốn LHM thu hoạch khoai tây 63 3.5 Thuật giải toán tối ƣu 65 CHƢƠNG KHẢO SÁT XÁC ĐỊNH BỀ RỘNG VÀ VẬN TỐC LÀM VIỆC TỐI ƢU CỦA LHM THU HOẠCH KHOAI TÂY CHỊU ẢNH HƢỞNG CỦA MỘT SỐ YẾU TỐ SỬ DỤNG 67 4.1 Ảnh hƣởng chiều dài đƣờng làm việc 67 iv 4.2 Ảnh hƣởng lực cản riêng máy thu hoạch khoai tây .68 4.3 Ảnh hƣởng độ sâu đào 69 4.4 Ảnh hƣởng tốc độ làm việc 70 4.5 Quan hệ bề rộng vận tốc tối ƣu .71 KẾT LUẬN VÀ KIẾN NGHỊ 73 TÀI LIỆU THAM KHẢO 74 PHỤ LỤC 76 v DANH MỤC BẢNG Bảng Tình hình sản xuất khoai tây Việt Nam 2000 – 2010 Bảng Hệ thống kiểu cỡ máy kéo nông nghiệp 18 Bảng Phân loại máy kéo theo lớp lực kéo .42 vi DANH MỤC HÌNH Hình 1.1: Máy thu hoạch khoai tây giới Hình 1.2: Máy thu hoạch khoai tây nhóm nghiên cứu thuộc Học Viện Nơng nghiệp Việt Nam thiết kế chế tạo Hình 2.1 Đặc tính nén đất Hình 2.2 Đặc tính cắt đất 10 Hình 2.3 Ảnh hƣởng tải trọng pháp tuyến đến khả chống cắt đất .12 Hình 2.4 Sự phụ thuộc lực cắt đất T vào tải trọng pháp tuyến N 13 Hình 2.5 Ứng suất sinh đất tác dụng mẫu bám bánh xe 13 Hình 2.6 Sự phụ thuộc hệ số ma sát nghỉ fn hệ số ma sát trƣợt fδ vào áp suất p .14 Hình 2.7 Sự phụ thuộc ứng suất cắt vào biến dạng 15 Hình 2.8 Sự phụ thuộc ứng suất cắt giới hạn ηδ ứng suất pháp ζ 15 Hình 2.9 Sơ đồ vận tốc điểm vành ngồi bánh xe lăn khơng trƣợt 20 Hình 2.10 Sơ đồ xác định toạ độ điểm bánh xe lăn không trƣợt 21 Hình 2.11 Sơ đồ vận tốc quĩ đạo chuyển động bánh xe .23 Hình 2.12 Đặc tính biến dạng pháp tuyến biến dạng tiếp tuyến lốp .24 Hình 2.13 Đặc tính biến dạng góc lốp 25 Hình 2.14 Sơ đồ xác định loại bán kính bánh xe 27 Hình 2.15 Ảnh hƣởng độ trƣợt đến hệ số cản lăn 30 Hình 2.16 Sơ đồ nghiên cứu độ trƣợt bánh chủ động 31 Hình 2.17 Đặc tính trƣợt máy kéo 33 Hình 2.18 Sơ đồ máy kéo với máy nơng nghiệp treo vận chuyển 34 Hình 2.19 Sơ đồ lực tác dụng lên phận làm việc cày 35 Hình 2.20 Sơ đồ liên hợp máy cày treo làm việc 36 Hình 2.21 Quan hệ hiệu suất kéo lực kéo móc 39 Hình 2.22 Đồ thị cân công suất máy kéo 40 vii Hình 2.23 Đƣờng đặc tính kéo máy kéo .40 Hình 2.24 Đƣờng đặc tính kéo máy kéo dùng hộp số học 43 Hình 2.25 Máy đào khoai tây 47 Hình 2.26 Sơ đồ lƣỡi đào phẳng .48 Hình 2.27 Lực tác động lên lƣỡi đào 49 Hình 2.28 Bộ phận sàng tách 50 Hình 2.29 Đƣờng đặc tính hiệu suất kéo máy kéo 54 Hình 2.30 Sự phụ thuộc hiệu suất kéo hk, chi phí nhiên liệu riệng ge suất W vào vận tốc làm việc V 56 Hình 3.1 Sơ đồ khối mơ hình liên hợp máy (5 phần tử) 59 Hình Đƣờng đặc tính tự điều chỉnh động .59 Hình 3.3 Đồ thị xác định bề rộng vận tốc làm việc tối ƣu liên hợp máy 65 Hình 4.1 Ảnh hƣởng chiều dài đƣờng làm việc tới suất chi phí nhiên liệu riêng 67 Hình 4.2 Ảnh hƣởng lực cản riêng đến suất chi phí nhiên liệu riêng 68 Hình 4.3 Ảnh hƣởng độ sâu làm việc tới suất chi phí nhiên liệu riêng .69 Hình 4.4 Ảnh hƣởng vận tốc làm việc đến chi phí nhiên liệu riêng 70 Hình 4.5 Ảnh hƣởng vận tốc làm việc đến suất LHM 71 Hình 4.6 Quan hệ bề rộng vận tốc tối ƣu LHM 72 viii Từ đặc tính hình 4.2 ta xác định đƣợc bề rộng vận tốc làm việc tối ƣu liên hợp máy ứng với lực cản riêng khác nhận đƣợc suất chi phí nhiên liệu riêng tƣơng ứng 4.3 Ảnh hƣởng độ sâu đào Các quan hệ vận tốc, suất chi phí nhiên liệu riêng với bề rộng làm việc thay đổi giá trị h khác Nếu chọn trƣớc lực cản riêng K chiều dài LR bề rộng vận tốc làm việc tối ƣu phụ thuộc vào độ sâu đào h nhƣ hình 4.3 Qua ta thấy ảnh hƣởng lực cản riêng h đến thông số tối ƣu Btu Vtu đáng kể thay đổi theo quy luật gần nhƣ tuyến tính Từ đồ thị hình 4.3 ta xác định đƣợc bề rộng vận tốc làm việc tối ƣu liên hợp máy ứng với độ sâu đào khác nhận đƣợc suất chi phí nhiên liệu riêng tƣơng ứng ANH HUONG DO SAU DAO DEN NANG SUAT va CHI PHI NHIEN LIEU RIENG K0 = 40kN/m2; Lr = 50m; B = 0.85 m 16 0.8 Wmax [ha/h] Wmax gcmin 0.7 gcmin [kg/ha] 14 0.6 12 12 0.5 10 10 0.4 8 0.3 6 17 0.2 10 11 12 13 Do sau dao h[cm] 14 15 16 Hình 4.3 Ảnh hƣởng độ sâu làm việc tới suất chi phí nhiên liệu riêng 69 4.4 Ảnh hƣởng tốc độ làm việc Các quan hệ suất chi phí nhiên liệu riêng với vận tốc làm việc thay đổi Nếu chọn trƣớc lực cản riêng K0 chiều dài LR, độ sâu đào h, suất chi phí nhiên liệu phụ thuộc vào tốc độ làm việc nhƣ hình 4.5, 4.6 Qua ta thấy ảnh hƣởng tốc độ làm việc đến thông số suất tối đa chi phí nhiên liệu tối thiểu đáng kể có cực trị tƣơng ứng với bề rộng làm việc khác Từ đồ thị khảo sát, ta xác định đƣợc suất chi phí nhiên liệu liên hợp máy ứng với vận tốc làm việc khác ANH HUONG VAN TOC DEN CHI PHI NHIEN LIEU RIENG gc K0 = 40kN/m2; Lr = 50m; h = 12 cm 13 Chi phi nhien lieu rieng gc [kg/ha] 12 11 B=0.4 10 B=0.6 B=0.7 B=0.8 B=1.0 0.5 1.5 Van toc V [m/s] 2.5 3.5 Hình 4.4 Ảnh hƣởng vận tốc làm việc đến chi phí nhiên liệu riêng 70 ANH HUONG VAN TOC DEN NANG SUAT K0 = 40kN/m2; Lr = 50m; h = 12 cm 0.6 Nang suat W [ha/h] 0.5 B=1.0 B=0.8 B=0.7 B=0.6 0.4 B=0.4 0.3 0.2 0.1 0 0.5 1.5 Van toc V [m/s] 2.5 3.5 Hình 4.5 Ảnh hƣởng vận tốc làm việc đến suất LHM 4.5 Quan hệ bề rộng vận tốc tối ƣu Quan hệ bề rộng vận tốc tối ƣu LHM thu hoạch khoai tây đƣợc xác định sở kết tính tốn ảnh hƣởng điều kiện, chế độ làm việc nhƣ chiều dài đƣờng làm việc, hệ số cản, độ sâu đào vv, tới suất tối đa chi phí nhiên liệu riêng tối thiểu Đồ thị 4.6 biểu diễn quan hệ bề rộng vận tốc tối ƣu LHM thu hoạch khoai tây điều kiện làm việc cụ thể Dựa vào quan hệ xác định chế độ làm việc hợp lý sở để tính tốn thành lập LHM 71 QUAN HE GIUA BE RONG VA VAN TOC TOI UU K0 = 40kN/m2; Lr = 50m; h = 12 cm 2.5 Vtu [m/s] 1.5 0.5 0.4 0.5 0.6 0.7 Be rong B[m] 0.8 0.9 Hình 4.6 Quan hệ bề rộng vận tốc tối ƣu LHM 72 KẾT LUẬN VÀ KIẾN NGHỊ - Các kết luận từ kết nghiên cứu Quá trình làm việc liên hợp máy thu hoạch khoai tây diễn phức tạp phụ thuộc vào nhiều yếu tố Trong phạm vi nghiên cứu đề tài có giới hạn lựa chọn số yếu tố mơ hình nghiên cứu, đáp ứng u cầu đơn giản hóa việc tính tốn mà đảm bảo độ xác cần thiết Mơ hình tốn đề tài mơ tả đƣợc q trình làm việc LHM thu hoạch khoai tây, với liên hệ chặt chẽ phần tử: Máy kéo – Máy công tác – Điều kiện sử dụng Mơ hình sử dụng để nghiên cứu chế độ làm việc, khảo sát thông số kết cấu điều kiện sử dụng LHM thu hoạch khoai tây Phƣơng pháp xác định tối ƣu bề rộng làm việc vận tốc LHM thu hoạch khoai tây mở rộng phạm vi nghiên cứu ứng dụng để phục vụ trình thiết kế máy thu hoạch khoai tây thành lập LHM thu hoạch khoai tây Dựa kết nghiên cứu quan hệ bề rộng vận tốc làm việc (hình 4.6), điều kiện làm việc cụ thể (độ sâu đào 12 cm, lực cản riêng đất 40 kN/m2, chiều dài ruộng 50 m), xác định vận tốc làm việc tối ƣu Liên hợp máy Kubota B2420 máy đào khoai bề rộng 0.85 m 1.75 m/s, tƣơng đƣơng với số truyền chế độ làm việc định mức động - Các kiến nghị nghiên cứu Tiếp tục hồn thiện mơ hình nghiên cứu theo hƣớng bổ sung yếu tố ảnh hƣởng thuộc phạm vi tổ chức sử dụng cho công việc xác định tối ƣu bề rộng, vận tốc làm việc LHM thu hoạch khoai tây Tiếp tục tổ chức nghiên cứu thực nghiệm để xác định xác số hệ số công thức thực nghiệm, đồng thời kiểm chứng lại mơ hình nghiên cứu 73 TÀI LIỆU THAM KHẢO Tiếng Việt 1.Võ Tấn Thặng (1994), Sử dụng máy nông nghiệp, NXB Nông nghiệp Hà Nội Nguyễn Hữu Cẩn, Dƣ Quốc Thịnh, Phạm Minh Thái, Nguyễn Văn Tài, Lê Thị Vàng (2000), Lý thuyết ô tô – máy kéo, NXB Khoa học kỹ thuật, Hà Nội Vũ Liêm Chính, Phan Nguyên Di (2001), Giáo trình Động lực học máy, NXB Khoa học & Kỹ thuật, Hà Nội Nguyễn Phùng Quang (2003), Matlab & Simulink, NXB Khoa học kỹ thuật Nguyễn Ngọc Quế(2007), Giáo trình tơ máy kéo xe chun dùng, Đại học Nông nghiệp Hà Nôi, Hà Nội Nguyễn Ngọc Quế, Ảnh hưởng sơ đồ truyền động đến tính chất kéo bám máy kéo làm việc dốc ngang, Tạp chí Nơng nghiệp Phát triển nông thôn, số 09/2006, trang 42-44 Nguyễn Anh Tuấn (2006) Xác định thời điểm sang số tối ưu hệ thống truyền lực thủy xe xích quân sự, Luận án tiến sỹ kĩ thuật, HVKTQS, Hà Nội Nguyễn Khắc Trai (2007), Tính điều khiển quỹ đạo chuyển động ô tô, NXB Giao thông vận tải Đỗ Sanh (2004), Cơ học, tập hai, động lực học, NXB Giáo dục Hà Nội 10 Nông Văn Vìn (2000), Động lực học chuyển động máy kéo – ô tô, giáo trình, Đại học Nông nghiệp I, Hà Nội 11 Nơng Văn Vìn, Hàn Trung Dũng, Phương pháp xây dựng đường đặc tính kéo lý thuyết – thực nghiệm máy nơng nghiệp, Tạp chí khoa học kỹ thuật Nông nghiệp số & 5-2006, trang 244 - 251 12 Đào Hữu Đồn, Nguyễn Ngọc Quế, Nơng Văn Vìn, Phạm Thị Thu Hằng, Phương pháp xác định tiêu kéo bám phận di động xích Tạp chí Khoa học Phát triển 2013, Tập 11, Số 8: Tr 1142-1149 13 Nguyễn Ngọc Oanh (2012), Nghiên cưu thành lập Liên hợp mayskeos phù hợp với điều kiện làm đất vùng đồng bắc bộ, Luận văn Thạc sỹ kỹ thuật, 74 Trƣờng Đại học Giao thông vận tải Tiếng nƣớc 14 M G BEKKER (1968), Introduction to Terrain-Vehicle Systems, The University of Michigan Press 15 J.Y Wong, Ph.D D.Sc (2001), Theory of ground vehicles, Department of Mechanical and A erospace Engineering Carleton University, Ottawa Canada 16 Gheorghe Brasov (1974), Dynamische Belastungen des Antriebs von Kettenschleppern bei konstanten Betriebsbedingungen Grundl Landtechnik Bd 24, Nr.6 17 Kheiralla A.F Yousif G.A (2012) Conceptual Design of a Rubber Tracked Minivehicle for Small Holders Using Off-Road Vehicle Engineering Techniques International Conference on Trends in Industrial and Mechanical Engineering, March 24-25, Dubai 75 PHỤ LỤC Liên hợp máy thu hoạch khoai tây Máy kéo Kubota B2420 - Động Diesel xi lanh công suất 24 mã lực - Hộp số khí số tiến, số lùi - Bộ phận di động bánh lốp - Trọng lƣợng máy kéo 595 kg - Khoảng tốc độ tiến từ 1.4 đến 18.8 km/h, lùi từ 1.4 đến 9.1 km/h 76 Máy thu hoạch khoai tây Thông số kỹ thuật máy đào khoai tây Cấu trúc liên kết với máy kéo Treo Kích thƣớc Dài x Rộng x cao 1770 × 1150 × 900mm Trọng lƣợng máy 230Kg Máy kéo hỗ trợ làm việc 22.1-29.4 Hp Lƣỡi đào đất củ Lƣỡi phẳng Chiều rộng làm việc 850mm Năng suất làm việc 0.10-0.2ha / h Độ sâu đào 150-250 mm Kiểu tách đất củ Lƣới chuỗi (Thanh lƣới lắp đĩa xích quay kết hợp rung) Khoảng cách rãnh luống 800-1200 mm Tốc độ hoạt động máy kéo 1.2-5.0 Km / h 77 Chƣơng trình tính tốn Matlab-Simulink % ANH HUONG CUA KET CAU VA DIEU KIEN SU DUNG DEN CHI PHI NHIEN LIEU VA NANG SUAT LHM THU HOACH KHOAI TAY % LV 2016 % Hoc vien: Do Dinh Thi %===================================================== clc clear all % CAC SO LIEU BAN DAU G = 5.95 ; % KN luong mk etam=0.85; f = 0.1; % He so can lan eps1= 0; eps2= 4.07; % Cac he so de tinh luc can may dao eps1= 0; eps2= 2.07; % Cac he so de tinh luc can may dao Ad= 16.6953; Bd= 1813.7/100 ;%kN Tinh truot Pfi=Bd; Vqv= 2.0 ; Lr = 50 ; % m/s - Van toc quay vong % m - chieu dai ruong a1= 10.96; a2= 45.40 ; % Cac he so de tinh He so sd thoi gian lam viec GeH= 5.5 ; %kg/h - chi phi nhien lieu gio cuc dai NeH= 17.65 ; % kW - Cong suat danh nghia cua dong co 78 geH= 0.213 ; % kg/kWh chi phi nhien lieu rieng cua dong co gama1= 1.2386; gama2 = 2.3753; % ====== TINH TOAN KHAO SAT ================ Lr = 50 ; %m hc=12 ; % cm Do sau dao h=hc/100 ; % m K0= 40 ; % kN/m^2 Luc can rieng cua may dao ch1=['K0 = ' num2str(K0) 'kN/m^2']; ch2=['; Lr = ' num2str(Lr) 'm']; ch3=['; h = ' num2str(hc) chtit= [ch1 ' cm'] ch2 ch3]; % Ky hieu dieu kien khao sat no=312 % toc quay PTO v/ph n1=173 % toc quay sang v/ph B0=0.4: 0.15: 1.05; % Cho B thay doi m=length(B0); %================================== for j= 1:m B= B0(j); Pk=0; Ne=0; i=1; dV=0.01; V=-dV; W=[]; gcay=[]; V0=[]; while (Pk