Nghiên cứu ảnh hưởng của một số thông số đến khả năng di động của liên hợp máy thu hoạch nghêu

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 ̃ NGUYÊN VĂN VỮNG ́ NGHIÊN CỨU ẢNH HƯỞNG CỦA MÔṬ SÔ THÔNG ́ ́ SÔ ĐÊN KHẢ NĂNG DI ĐÔNGG̣ CỦA LIÊN HỢP MÁY THU HOẠCH NGHÊU CHUYÊN NGÀNH: KỸ THUẬT CƠ KHÍ MÃ SỐ: 60520103 LUẬN VĂN THẠC SỸ KỸ THUẬT CƠ KHÍ NGƯỜI HƯỚNG DẪN KHOA HỌC TS BÙI VIỆT ĐỨC Hà Nội, 2014 ̉̉ ̉̀ MƠ ĐÂU Tính cấp thiết vấn đề nghiên cứu Nước ta có bờ biển dài, có tiềm lớn phát triển ni trồng thủy hải sản Trong số lồi thủy hải sản, nuôi ngao (miền Bắc) hay nghêu (miền Nam) nghề đem lại hiệu kinh tế lớn, làm hồi sinh nhiều vùng đất ven biển trước bị bỏ hoang, tạo việc làm thu nhập ổn định cho phận cư dân nghèo sinh sống ven biển, góp phần xóa đói giảm nghèo đảm bảo an sinh xã hội Tiềm định hướng phát triển nghề nuôi nghêu thể rõ qua định phê duyệt quy hoạch phát triển ni nhuyễn thể hàng hóa tập trung đến năm 2020, Bộ Nông nghiệp Phát triển nơng thơn; có Quy hoạch vùng ni nghêu, sò đến năm 2015 năm 2020 cho tỉnh đồng sông Cửu Long (ĐBSCL) như: Cà Mau, Bạc Liêu, Tiền Giang, Sóc Trăng, Kiên Giang, Bến Tre Trà Vinh, với tổng diện tích đến năm 2015 28.110 năm 2020 phát triển lên 35.690 Quyết định Quy hoạch nêu rõ: đến năm 2015, diện tích ni nghêu 15.950 ha, sản lượng 142.700 đến năm 2020, diện tích ni nghêu 20.590 ha, sản lượng 206.300 Một công việc tiêu tốn nhiều sức lao động q trình ni nghêu khâu thu hoạch Hiện phương pháp thu hoạch nghêu chủ yếu vùng chăn ni nghêu hồn tồn sử dụng phương pháp thu hoạch thủ cơng, suất thấp chi phí lao động cao Ngồi phương phương pháp người dân dùng phương pháp khác dùng bơm tạo dịng nước có áp lực để phân loại làm nghêu để thu hoạch Với diện tích sản lượng lớn, tăng hàng năm vậy, việc giới hóa khâu thu hoạch nghêu, mà cụ thể thiết kế, chế tạo máy thu hoạch nghêu đáp ứng yêu cầu kỹ thuật - công nghệ, thay cho công việc thu hoạch thủ công việc làm cấp thiết, nhằm nâng cao suất chất lượng sản phẩm, giảm chi phí sản xuất, thỏa mãn tối đa tính thời vụ yêu cầu ngày cao cho tiêu dùng nước xuất Để tạo LHM thu hoạch nghêu hồn chỉnh với đầy đủ tính kinh tế - kỹ thuật – cơng nghệ hồn chỉnh, đáp ứng tốt yêu cầu công việc, cần phải thực đầy đủ nhiều công đoạn phức tạp khác liên quan đến đối tượng tác động (con nghêu), điều kiện làm việc (nền đất), nguồn động lực (máy kéo) phận công tác (máy thu hoạch) Trên sở thực phần công việc nằm đề án thiết kế chế tạo Liên hợp máy thu hoạch nghêu, thuộc chương trình hợp tác nghiên cứu khoa học phục vụ sản xuất cho vùng nuôi nghêu ven biển phía Nam Việt Nam Trường Đại học Nơng nghiệp Hà Nội Trường Trung cấp nghề Cơ Điện Đông Nam Bộ Đồng Nai Với lý trình bầy chúng tơi chọn thực đề tài: “Nghiên cứu ảnh hưởng số thông số đến khả di động LHM thu hoạch nghêu ” Mục tiêu nghiên cứu Xuất phát từ lý thực đề tài nêu đặt mục tiêu nghiên cứu là: Xác định ảnh hưởng số thông số đến khả di động LHM thu hoạch nghêu làm việc nhằm tính tốn thơng số, xây dựng quan hệ đất bánh để lựa chọn hệ thống di động hợp lý cho LHM thu hoạch nghêu nhằm nâng cao khả di động cho LHM, góp phần hồn chỉnh LHM với đầy đủ tính kinh tế - kỹ thuật – công nghệ, đáp ứng tốt yêu cầu công việc nhằm tăng suất hiệu sử dụng máy Nôịdung nghiên cứu đềtài Do thời gian nghiên cứu có hạn, đề tài giới hạn nội dung sau: - Xây dưngG̣ mô hinhh̀ tổng quát LHM thu hoacḥ nghêu - Tinh́ tốn thơng sốbơ G̣phâṇ di đơngG̣ - Xây dưngG̣ mô hinhh̀ bô pG̣ hâṇ di đôngG̣ LHM + Mô hinhh̀ xich́ cao su không biến dangG̣ + Mô hinhh̀ xich́ cao su biến dangG̣ - Xây dưngG̣ mơ hinhh̀ tính tốn - Khảo sát ảnh hưởng yếu tốkết cấu tG̣ hống di đơngG̣ Tính đề tài: Liên hợp máy thu hoạch nghêu nghiên cứu với tính đưa giới hóa khâu thu hoacḥ vào vùng chăn ni nghêu thay thếcho viêcG̣ thu hoacḥ thủ công hiêṇ vẫn đảm bảo suất vàchất lươngG̣ sản phẩm Tính khoa học đề tài: Đề tài phát triển vấn đề khoa học lý luận thực tiễn việc giới hóa khâu thu hoạch nghêu Ý nghĩa thực tiễn đề tài: Đề tài góp phần giới hóa khâu thu hoạch nghêu theo hướng nâng cao hiệu kinh tế, kỹ thuật, đảm bảo tính thời vụ khâu thu hoạch nghêu ̉̉ Chương ́ ̀ TÔNG QUAN CÁC VÂN ĐÊ NGHIÊN CỨU 1.1 Tởng quan vềtình hình sản xuất, chăn ni và thu hoạch nghêu khu vực phía Nam Trong chiến lược phát triển Nuôi trồng thuỷ sản nay, nghêu xem đối tượng ưu đầy triển vọng Với vai trò quan trọng làm thực phẩm, góp phần làm hồi sinh nhiều vùng đất ven biển trước bị bỏ hoang, tạo việc làm thu nhập ổn định cho phận cư dân nghèo sinh sống ven biển, góp phần xóa đói giảm nghèo đảm bảo an sinh xã hội Do đó, nghêu coi đối tượng chủ lực cho phát triển nuôi biển kỷ XXI Theo thống kê FAO (2006) tổng sản lượng động vật thân mềm ni tính đến năm 2004 13,25 triệu tấn, chiếm 22,3% tổng sản lượng thuỷ sản nuôi Ở nước ta nghêu phân bố tự nhiên khu vực vùng triều cửa sông ven biển tỉnh miền Tây Nam Bộ như: Gị Cơng Đơng (Tiền Giang), Cầu Ngang, Dun Hải (Trà Vinh), Bình Đại, Ba Tri, Thạnh Phú (Bến Tre), Vĩnh Châu (Sóc Trăng); nghêu Bến Tre người dân đưa vào nuôi thử nghiệm số vùng cửa sông ven biển cho kết tốt Nghêu đối tượng người dân quan tâm giá trị dinh dưỡng Thịt thơm ngon có nhiều chất dinh dưỡng Prơtêin chiếm 15,66%, Lipit chiếm 3,43%, khống chiếm 3-13% (Nguyễn Chính CTV) mặt hàng hải sản có giá trị xuất cao Đặc biệt tháng 10/2008 hội đồng bảo tồn biển quốc tế vừa cấp giấy chứng nhận tiêu chuẩn thương hiệu MSC (Marine Sterwarship Council) cho nghêu Bến Tre trở thành đặc sản biển khu vực Đông Nam Á Sự công nhận giúp nghêu Bến Tre vàcác tinhh̉ khu vưcG̣ có nhiều hội xuất sang nhiều thị trường khu vực giới Tính đến hết tháng 2-2011, tỉnh Bến Tre có 10 hợp tác xã thủy sản với gần 9000 hộ xã viên 35 tập đồn ni nghêu Các HTXTS tập đoàn biến gần 16.000 bãi bồi thành khu ni nghêu, có khoảng 9650 khai thác Riêng huyêṇ Thạnh Phú, Bình Đại Ba Tri tinhh̉ Bến Tre tập trung hầu hết HTX tập đồn ni nghêu với trữ lượng từ 32.000 – 40.000 tấn/năm, sản lượng thu hoạch 20.000 tấn/năm Ở TràVinh códiêṇ tich́ ni khoảng 1000 cho sản lươngG̣ 5.000 với suất binhh̀ quân tấn/ha Đến năm 2006 đầu năm 2007, diêṇ tich́ nuôi nghêu 20.000 với tổng sản lươngG̣ ởtồn bơ G̣khu vưcG̣ Nam Bơ đG̣ aṭ70 - 80 ngàn tấn/năm vàđa ̃trởthành măṭ hàng xuất lớn ( Lê Xuân Sinh, 2007 ) Khu vưcG̣ phiá nam vùng có bờ biển dài gần 6000 Km với diêṇ tich́ bãi bồi rôngG̣ lớn Ven biển thuộc địa phận tinhh̉ Bến tre, Tràvinh, CàMau, Bạc Liêu, Tiền Giang, Kiên Giang có nhiều cửa sơng lớn đổ biển Các cửa sông hàng năm đưa biển lượng phù sa lớn, bồi thành bãi triều rộng tương đối phẳng, cấu tạo trầm tích baĩ bồi chủ yếu dạng cát, sa, sét thích hợp cho nuôi nghêu nuôi trồng loại thủy sản khác Với vị trí địa lý điều kiện tự nhiên thuận lợi, tiềm để phát triển kinh tế thủy sản, nghề nuôi trồng thủy sản đócónghềni nghêu Những năm gần với việc phát triển chung ngành kinh tế, nghề nuôi trồng thủy sản có bước phát triển nhanh chóng diện tích, sản lượng giá trị Kết phát triển ni trồng thủy sản đóng góp đáng kể vào phát triển kinh tế khu vưcG̣ góp phần cải thiện đời sống cho nhân dân xã ven biển Mặc dù vùng bãi triều ven biển phia ́ nam có nhiều tiềm phát triển nuôi trồng thủy sản, vẫn cịn phần lớn diện tích chưa khai thác có hiệu quả; Đặc biệt nghề nuôi nghêu - Nghề mà vốn đầu tư sở hạ tầng không nhiều, hiệu kinh tế lại cao Một số năm gần đây, nhân dân xã ven biển tinhh̉ Tiền Giang, TràVinh, Bến Tre, Sóc Trăng tự phát đầu tư nuôi nghêu Song, chưa tổ chức quy hoạch, thiếu kỹ thuật chuyên ngành chế sách nên việc ni nghêu chưa đạt kết Theo quyết đinḥ phê duyêṭcủa Bộ Nơng nghiệp Phát triển nơng thơn; có Quy hoạch vùng ni nghêu, sị đến năm 2015 năm 2020 cho tỉnh đồng sông Cửu Long với tổng diện tích đến năm 2015 28.110 năm 2020 phát triển lên 35.690 với sản lươngg̣ tới năm 2020 357.300 Để khai thác có hiệu tiềm diện tích đất vùng bãi bồi ven biển phia ́ nam vấn đề đặt phải có nghiên cứu tổng thể từ khâu quy hoacḥ baĩ nuôi, kỹthuâṭnuôi tới khâu thu hoacḥ chi tiết để đưa vào nuôi nghêu, đáp ứng nhu cầu thị trường nguyện vọng nhân dân xã ven biển việc làm cần thiết, phù hợp với chủ trương phát triển nuôi trồng thủy sản năm tới; đảm bảo cân sinh thái phát triển nghề nuôi trồng thủy sản bền vững Ky ̃thuâṭnuôi nghêu: Kỹthuâṭ nuôi: Choṇ baĩ nuôi: Choṇ baĩ triều gần cửa sơng, it́ sóng gió, phẳng Khi triều xuống thấp cócảnăng phơi baĩ vài giờ/ngày Chất đáy nên choṇ làcát bùn Khơng choṇ nơi cónguy bi G̣nước thải nhiễm cơng vìnghêu nhaỵ cảm với môi trường ô nhiễm Bãi nuôi vùng triều tương đối cao, thời gian nước rút tương đối dài, chất đáy cứng, nhiều bùn, dẫn đến chết nghêu Sau choṇ đươcG̣ baĩ nuôi tiến hành vê G̣ sinh măṭbai,̃ doṇ tapG̣ chất, cối đáy Nến đất cứng phải xới lên cho tơi xốp, mềm giúp nghêu vùi minhh̀ tốt Dùng đăng tre hoăcG̣ lưới chắn xung quanh, cao 0,6 tới 0,7m, chân đăng vùi sâu 0,2 - 0,3m Dùng cocG̣ cắm cho đăng lưới đứng nghiêng vềphiá trong, cocG̣ cách 1,2m tới 1,5m Trong bai,̃ căng nhiều dây ngang docG̣ theo chiều gióđểhaṇ chếnghiêng theo hướng nước chảy vàgiólớn Thảgiống: Thảlàm đơṭ vào tháng 2, tháng vàtháng 7, tháng Đơṭ 1: Thảgiống kich́ cỡ3.000-5000 con/1kg Sau tới tháng ương thìsan thưa đểni lớn Đơṭ 2: Thảgiống kích cỡ500 - 800 con/kg Nên thảni ln, khơng cần san thưa, mâṭđô 3G̣ 00-450 con/kg Chúýkhi thảlúc triều lên đểnghêu dễvùi lớp bùn đáy Chăm sóc: Khi rải nghêu baĩ ni, vỏnghêu cịn yếu, chưa thểvùi sâu xuống đáy nên haṇ chếđi laịtrên baĩ nghêu Bắt hết loaịốc mỡtrơn baĩ triều thảnghêu giống Trong qtrình ni thường xun kiểm tra bắt ốc haịnghêu con, thường xuyên kiểm tra rào chắn, chân rào đểngăn chăṇ nghêu ngồi Nếu nghêu tâpG̣ trung ởmơṭchỡphải san giañ vìmâṭđơ G̣cao se ̃làm cho nghêu châm lớn Thu hoacḥ nghêu Nghêu đươcG̣ thu hoạch vào lúc triều rút, lúc chúng ăn no thải thức ăn thừa thãi, giữ lại nước nên thịt nghêu (ngược lại, lúc triều lên, nghêu thường ngậm cát, bùn) Sau ni - 10 tháng cóthểthu hoacḥ nghêu Nên lưạ choṇ cuối mùa mưa đểthu hoacḥ vìlúc nghêu béo Năng suất nghêu binhh̀ quân đaṭ20 - 25 tấn/ha Nếu người ni cókinh nghiêm vàkỹthṭcao cóthểđaṭ50 - 55 tấn/ha Hiêṇ khâu thu hoacḥ ngêu chỉlàthủcông dùng cào tới hoăcG̣ dùng tia nước áp lưcG̣ đểthu hoacḥ ngêu phícho khâu thu hoacḥ làrất lớn đồng thời không đảm bảo đươcG̣ tinh́ thời vu.G̣ Xuất phát từ vấn đề trên: “Nghiên cứu ảnh hưởng của số thông số đến khả di động của LHM thu hoạch nghêu ” góp phần hồn thiêṇ LHM thu hoacḥ nghêu áp dungG̣ cho vùng chăn nuôi nghêu thời gian tới vấn đề cần thiết cóýnghiã thưcG̣ tiễn 1.2 Đặc điểm tự nhiên vùng chăn nuôi nghêu Nghêu loài sống vùi phân bố bãi biển, eo vịnh có đáy cát pha bùn Theo nghiên cứu tiêu môi trường nghêu đồng sông Cửu Long Nguyễn Hữu phungG̣ ( 1996 ) cho thấy nghêu thường phân bố vùng có đáy cát hay cát bùn đótỷ lệ cát 68 - 75%, kich́ thước cỡhaṭtừ 0,006 - 0,25 mm, tỷ lệ sét 21 - 31%, đất thịt có tỷ trọng thấp < 7% Theo Trương Quốc Phú( 1999 ) Bãi nghêu phân bố thường gần cửa sơng có sóng gió nhẹ, có nguồn nước chảy vào, chủyếu ởvùng trung triều, hạ triều, nơi cóđơ dG̣ ốc tương đối phẳng Theo Nguyễn Tác An vàNguyễn Văn LucG̣ ( 1994 ) cho nghêu phân bốở vùng thời gian phơi baĩ từ - giờ/ngày Đô G̣ sâu cưcG̣ đaịtimh̀ thấy nghêu lúc nước ròng là2,5m Nghêu phân bốởvùng cónền đáy cát miṇ đến cát trung cópha lẫn hàm lươngG̣ bùn lỏng vàxác hữu ( 10 - 18 % ), vào mùa mưa bùn lỏng bao phủ đáy baĩ nghêu ( 1,5 - 2,5 cm ) Theo nghiên cứu VõSi ̃ Tấn ( 1999 ) ởGịCơng Đơng cho thấy nghêu tâpG̣ trung ởđơ cG̣ ao từ khoảng 0,8 - 1,5 m so với hải đồ Nghêu lồi rộng nhiệt: Chúng sống điều kiện nhiệt độ từ o o 35 C, khoảng nhiệt độ thích hợp cho nghêu phát triển 28-31 C, tốc độ dòng chảy 0,1 - 0,25 m/s, hàm lượng oxy hoà tan khoảng - mg/l, pH 6-9, độ mặn dao động từ 19-26% độ mặn phù hợp cho nghêu phát triển 22-25% Tập tính sống: Nghêu lồi sống đáy gặp điều kiện không thuận lợi nhiệt độ giảm, độ mặn thay đổi đột ngột chúng thường lên khỏi đáy di chuyển nơi khác có điều kiện thích hợp Theo Trương Quốc Phú (1999) độ mặn giảm xuống 5% hầu hết nghêu di chuyển tới vùng nơi có độ muối cao hơn, bãi nghêu có xu hướng dịch chuyển từ bờ xa vào mùa mưa mùa khơ có xu hướng tiến lại gần bờ nên diện phân bố nghêu vào mùa khô rộng so với mùa mưa Nghêu nói riêng động vật thân mềm hai mảnh vỏ nói chung, q trình sống trải qua giai đoạn - Giai đoạn ấu trùng sống trơi sống đáy: Sau trứng nở chuyển sang giai đoạn ấu trùng phù du Giai đoạn này, ấu trùng trôi nước Sự phân bố chúng phụ thuộc lớn vào dòng chảy thuỷ triều Kết thúc giai đoạn ấu trùng phù du giai đoạn sống đáy Lúc hình thành chân, màng áo khép vỏ Do giai đoạn cần đáy phẳng cần có vật bám cát to, sỏi - Giai đoạn trưởng thành: nghêu sống vùi đáy Để hơ hấp lấy thức ăn nước, nghêu thị ống hút nước lên mặt bãi ống thoát hút nước nghêu ngắn nên nghêu chui sâu, thường cách đáy vài cm Vào mùa lạnh nghêu vùi xuống sâu, khơng q 10 cm Nghêu lồi ăn lọc, chúng bắt mồi theo hình thức thụ động Khi triều dâng lên nghêu thị ống hút nước lên cát để lọc mồi ăn Nghêu thường chọn mảnh vụn hữu cơ, loài thực vật phù du Đặc điểm sinh trưởng: Sự sinh trưởng động vật thân mềm hai mảnh vỏ nói chung nghêu nói riêng có thay đổi theo lồi, vị trí địa lý phân bố, thời tiết, vùng trung triều hay hạ triều, khác mỗi cá thể mà di truyền tạo Tốc độ tăng trưởng nghêu phụ thuộc vào vùng phân bố nhiều hay thức ăn Nghêu phân bố vùng cửa sông phong phú thành phần thực vật phù du mùn bã hữu cơ, nghêu lớn nhanh Nghêu sống vùng triều thấp thuờng lớn nhanh vùng triều cao Ngồi chúng cịn phụ thuộc vào nhiều nhân tố độ mặn, nhiệt độ, sóng gió Sự sinh trưởng thay đổi từ năm đến năm khác khu vực mà có nhiệt độ biến đổi theo mùa Vào mùa xuân mùa hè nhiệt độ nước ấm lên, thức ăn dư thừa sinh trưởng tăng lên nhanh chóng Tổng thời gian từ sinh đến lúc thu hoạch trải qua 18 - 20 tháng Trương Quốc Phú (1999) nghiên cứu đặc điểm sinh trưởng nghêu cho thấy, nghêu sinh trưởng nhanh từ tháng - sinh trưởng chậm từ tháng 10-4 năm sau (tháng 10-11) nước từ sông Cửu Long đổ mạnh làm độ mặn xuống thấp, ngao vùi để tránh nước làm giảm thời gian bắt mồi Tốc độ sinh trưởng trung bình 1,72 mm/tháng 789,47 mg/tháng Tốc độ sinh trưởng khối lượng tương đối nhanh so với tốc độ sinh trưởng chiều dài Ở vùng biển khác nhau, tốc độ sinh trưởng nghêu khác Cùng năm tuổi, nghêu Trà Vinh cỡ 20 mm nặng trung bình 2,7 g/con, Duyên Hải 3,7 g/con, Tiền Giang ngao cỡ 25 mm nặng 2,8 g/con 81 Chương trinh̀ mô Chương trình tính tốn Matlab % KHAO SAT TINH CHAT KEO BAM CUA MAY KEO XICH clc clear all % SO LIEU DAU VAO: % % W=10507; % Trong luong xe kw=W/W00; % He so tang luong de tinh lai Do lun chW=['(Khi G=', num2str(W), 'N)']; hg=0.25; D= 140; B=35 ; % % Chieu cao tam % chieu dai tiep dat 140 % be rong xich B=25; Rf=6.5; % ban kinh banh truoc 25 Rk=13; %% Ban kinh banh sa0 chu dong 25 Rr=9 ; % Ban kinh banh sau Rp=D/5; H=3 ; % Chieu cao mau bam cm H0=6000; % [N]= Luc cang xich Ban dau; H0=6000 Gp=9; % Buoc vau xich e=0.02; % e/D Do lech tam ld=88 ; % khoang cach tu diem moc keo den duong thang di qua diem giua cua may keo hd=38; % m Chieu cao diem moc keo Vphay= 1.00; % m/s Ttoc vong cua banh chu dong 82 %k1= 01.167/10; %k2= 24.17/10; % N/cm^(n2+2) %n1=0.188; %n2=0.798181; %mc= 4.5; %c0=0.100; pm=W/(2*B*D); pm0=pm; k1=0.49147; k2=0.56107; mc=0.46387; %mf=0.623; mf=0.3630; n1=0.8664*0.5; % n1=0.8664 k3=30.3; k4=30.3; n4=0.632; % -teta0=0.000; tam1=1000; tam2=0; dT=0.001; % Buoc bien thien i=0; while (abs(tam1-tam2)>50) i=i+1; % % sf0=4.5556; ds=0.001; sm=5.5; %G=18000 p0= k1*H^n1; %for i=1: 1540 sr0=sm+i*ds; teta0= atan((sr0-sf0)/D); e0=e+(hg+H)*tan(teta0)/D; %pf0=k1*(sf0)^n1; 83 %pr0=p0 + k2*(sr0-H)^n2; %X1=(H-sf0)/tan(teta0); %tam1=W*cos(teta0); %tam2=(2*B*D/(sr0-sf0))*((k1/(n1+1))*((H^(n1+1))-sf0^(n1+1))+(sr0H)*(p0+(k2/(n2+1))*(sr0-H)^n2)); Teta0=teta0*180/pi; pf0=p0 + k2*(sf0-H)^n2; pr0=p0 + k2*(sr0-H)^n2; tam1=W*cos(teta0); tam2=2*B*D*p0+(2*k2*B*D/(n2+1))*((sr0-H)^(n2+1)-(sf0H)^(n2+1))/(sr0-sf0); test = tam1-tam2; end thu1=W*cos(teta0)*(0.5+e0)*D; thu2=((2*B*D^2)/(sr0-sf0)^2)*((k2/(n2+2)*((sr0-H)^(n2+2)-(sf0-H)^(n2+2)))- (k2/(n2+1)*(sf0-H)*((sr0-H)^(n2+1)-(sf0-H)^(n2+1)))+(0.5*p0*(sr0^2-sf0^2)p0*sf0*(sr0-sf0))); KT=thu1-thu2; Mc=' e M=[ e fprintf(' fprintf(' fprintf(' test= %6.4f Teta= %6.4f \n',test, Teta0) fprintf(' tam1= %6.4f pm= %6.4f\n',tam1,tam2,p0,pm) fprintf(' KT= %6.4f fprintf(' thu1= %6.4f thu2=%6.4f x=0:0.01:D; s0x= sf0 + (sr0-sf0)*x./D; 84 p0x=p0 +k2*(s0x - H).^n2; % tt %====================================================== tetati=teta0; T3i=0.001; i=0; dR=0.1 ; %=1; % truot dR=1 [%] phan tram %while (abs(T3-T3i)>epsilon) for i=1: 600 %=60 diem tinh id0= i*dR; id= id0/100; tetai=teta0; idphay=1-(1-id)/cos(tetati); jB=id*D/(1-id); dX=0.05; %sO LIEU CU (dX=0.1) N= round(D/dX); TPpiX=0; TPtmX=0; MangTP=[0]; MangtmX=[0]; for i2=1:N X= i2*dX; XV(i2)=X; s0iX = sf0+(sr0-sf0)*(X/D); piX = k1*H.^n1 + k2*(s0iX - H).^n2; PIX1(i2)=piX; jx = jB + idphay*X; tmX = (mc + mf*piX)*(1-exp(-a*(jB + idphay*X))); TmX(i2)=2*B*tmX; TmS(i2)= 4*H*(mc + mf*piX/pi*acot(H/B))*(1-exp(-a*(jB+idphay*X))); %Luu vao mang TPpiX=TPpiX+piX*dX; 85 TPtmX=TPtmX+tmX*dX; MangTP=[MangTP TPpiX]; MangtmX=[MangtmX TPtmX]; end MangHieuTP=[0]; MangTX=[0]; MangKQ=[0]; MangKQ2=[0]; MangS0ix=[3]; MphayiX=[0]; MT0X=[0]; MtphaymX=[0]; XM1=round(0.492*Rp/dX); Xm1=XM1*dX; Mangpphay=[0]; tpphayix=0; MangX2=[0]; tpiX2=0; Mangtam=[0]; TpphaymX=0; MangtphaymX=[0]; TpphaymS=0; MangTpphaymS=[0]; for ih=2:N Xh=ih*dX; XVV(ih)=X; HieuTP=B*(MangTP(XM1)-MangTP(ih)); MangHieuTP=[MangHieuTP HieuTP]; TX=H0+B*MangtmX(ih); MangTX=[MangTX TX]; 86 %Doi dau neu ih >XM if ih>Xm1 KQ=(Xm1-Xh)*(MangHieuTP(ih)/((-1)*MangTX(ih))); end if MangTX(ih)==0 KQ=0; else KQ=(Xm1-Xh)*(MangHieuTP(ih)/MangTX(ih)); end %Luu KQ phep chia vao mang KQ MangKQ=[MangKQ KQ]; if (ih round(Rp/dX) if mod((ih-1),round(Rp/dX))==0 tam=ih; end s0iXm1 = (((sr0-sf0)*Xm1)/D+sf0)-MangKQ(tam); end % Thay doi gia tri XM sau moi mot vong+0.0124999999*Rp/dX if mod(ih,round(Rp/dX))==0 XM1=round(XM1+Rp/dX); Xm1=XM1*dX; end 87 sphay0iX=(MangKQ(ih)+s0iXm1-((sr0-sf0)/D)*(Xm1-Xh)); MangS0ix=[MangS0ix sphay0iX]; % Tinh tich phan pphayiX va luu vao mang Mangpphay pphayiX=k1*H^n1+k2*(sphay0iX-H)^n2; PIX11(ih)=pphayiX; ppphayXX(ih)=real(pphayiX); MphayiX=[MphayiX ppphayXX(ih)]; jf= jB + idphay*X; tphaymX = (mc + mf*pphayiX)*(1-exp(-a*jf)); TpphaymX =real(TpphaymX + tphaymX*dX); MangtphaymX=[MangtphaymX TpphaymX]; MangtphaymX(ih)=2*B*MangtphaymX(ih); MtphaymX=[MtphaymX TpphaymX]; TphaymS= 4*H*(mc + mf*pphayiX/pi*acot(H/B))*(1-exp(-a*(jf))); TpphaymS = real(TpphaymS+TphaymS*dX); MangTpphaymS=[MangTpphaymS TpphaymS]; T0X(ih)=H0+B*MtphaymX(ih); end MT0X=[MT0X T0X]; %N0=1400 TphaymSS(i) = MangTpphaymS(N); T0Xx(i)=MT0X(N) % Tmbb=trapz(XV,TmX);% -1 Tmss=trapz(XV,TmS);% -2 TPtmXX(i2)=real(Tmbb); XX(i2)=X; id0V(i)=id0; Tmb(i)=real(Tmbb); 88 Tms(i)=real(Tmss); % pphayiXX(ih)= real(pphayiX); XXX(ih)=X; tetat0i= teta0; %arctan((sr0-sf0)/D); tetaf0 =acos(cos(tetat0i)-sf0*cos(tetat0i)/(Rf+H))-tetat0i; dT=0.01; kk=tetaf0/dT; M=round(kk); for i3= 1:M teta= i3*dT; tetamV(i3)= teta; %luu lai sf0iteta = (Rf+H)*(cos(teta+tetat0i)-cos(tetaf0-tetat0i))/cos(teta+tetat0i); %pfiteta = k1*H^n1 + k2*(sf0iteta -H)^n2; pfiteta=k1*(sf0iteta)^n1; pfitetaV(i3)= pfiteta; jfteta = (Rf+H)*((tetaf0 - teta)-(1-id)*(sin(tetaf0+teta)- sin(teta + tetat0i))); tfiteta = (mc + mf*pfiteta)*(1-exp(-a*jfteta)); tfitetatphan(i3) = tfiteta*cos(teta) - pfiteta*sin(teta); tfsteta = 4*H*(mc+mf*pfiteta/pi*acot(H/B))*(1-exp(a*jfteta)); tfstetatphan(i3)= tfsteta*cos(teta); tetatphan(i3)= teta; end Tfbb = 2*B*(Rf+H)*trapz(tetatphan,tfitetatphan); % -3 Tfss=(Rf+H)*trapz(tetatphan,tfstetatphan); Tfb(i)=real(Tfbb); Tfs(i)=real(Tfss); % tetar0i= teta0; %arctan((sr0-sf0)/D); % -4 89 tetar0 = tetar0i; ddelta =0.001; n= round(tetar0/ddelta); for i4= 1:n deltatphan(i4) = i4*ddelta; delta= i4*ddelta; sr0idelta=((Rr+H)*(cos(tetat0i)-cos(tetaf0+tetat0i))+ D*sin(tetat0i)+(Rr+H)*(cos(delta)-cos(tetat0i)))/cos(delta); jrdelta =( Rr+H)*((tetat0i-delta)-(1-id)*sin(tetat0i)- sin(delta)); pridelta = k1*H^n1 + k2*(sr0idelta-H)^n2; trdelta = (mc+mf*pridelta)*(1-exp(-a*jrdelta)); trdeltatphan(i4)= trdelta*cos(tetat0i-delta) + pridelta*sin(tetat0i-delta); trsdelta= 4*H*(mc+mf*pridelta*acot(H/B)/pi)*(1-exp(-a*jrdelta)); trsdeltatphan(i4)=trsdelta*cos(tetar0i - delta); end Trbb= 2*B*(Rr+H)*trapz(deltatphan,trdeltatphan); Trss= (Rr+H)*trapz(deltatphan,trsdeltatphan); Trb(i)=real(Trbb); Trs(i)=real(Trss); T3(i)= Tmb(i)+Tms(i)+Tfb(i)+ Tfs(i)+Trb(i)+Trs(i); %============================ tetat0i= teta0; %arctan((sr0-sf0)/D); tetaf0 =acos(cos(tetat0i)-sf0*cos(tetat0i)/(Rf+H))-tetat0i; jfs= (Rr +H)*sin(tetaf0)*id/(1-id); Tongf=0; dT=0.01; M= round(tetaf0/dT); 90 for m=1:M % CHUA DUNG tetam = tetaf0*(1-m/M); sf0itetam= (Rf + H)*(cos(tetam + teta0) - cos(tetaf0 + teta0))/cos(tetam + teta0); sfi= (sf0 + sf0itetam)*cos(teta0); pfitetam= k1*H^n1 + k2*(sfi - H)^n2; PFteta(m)=pfitetam; Tongf = Tongf + c0*((pfitetam*cos(tetam))^c1 * ((m/M*jfs)^c2 - ((m-1)/M*jfs)^c2)); end sfs= Tongf; % -idphay= 1- (1-id)/cos(teta0); js= idphay*D/(1-idphay); Tongr=0; %dX=0.1; for n=1:N N= round(D/dX); % CHUA DUNG X= n*D/N; s0iX = sf0+(sr0-sf0)*(X/D); piX = k1*H.^n1 + k2*(s0iX - H).^n2; PIX(n)=piX; Tongr = Tongr + c0*(piX*n*D/N)^c1 *((n*js/N)^c2 - ((n-1)*js/N)^c2); end srs = sfs + Tongr; sri= (sr0 + srs)*cos(teta0); T2(i)= 2*k1*B*H^(n1+1)/(n1+1) + 2*k1*B*H^n1*(sri - H) + 2*k2*B/(n2+2)*(sri-H)^(n2+1); %T2(i)=real(T2); %========================================== 91 Vt= Vphay ; % m/s - Van toc ly thuyet V= Vt*(1- id) ; % m/s - Van toc thuc te T1 (i) = T3(i); T4(i)= T3(i)/cos(teta0) - W*tan(teta0) - T2(i); P(i)= W/cos(teta0) - T3(i)*tan(teta0); % ================================= tetati=asin((sri-sfi)/D); z(i)=sri-(T2(i)*(Rr+(hd-Rr)*cos(tetati)-(ld0.5*D)*sin(tetati))- W*(hg*sin(tetati)-D*(0.5-e)*cos(tetati)))/(T2(i)+W*cos(tetati)/ tan(tetati)); % ========================================== E1(i)= T3(i)*Vt/1000; % kW - Cong suat chu dong E2(i)= T2(i)*V/1000 ; % kW - Cong suat chi phi cho luc can E3(i)= (T3(i)*(1/(1-id) - 1/cos(teta0))*V + W*V*tan(teta0))/1000; if T4