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ỘT SỐ THÔNG SỐ ĐẾN KHẢ NĂNG DI ĐỘNG 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 nuôi trồng thủy hải sản Trong số loà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 nuôi 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 nuôi 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 nuôi nghêu 15.950 ha, sản lượng 142.700 đến năm 2020, diện tích nuôi 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 trình nuôi 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 nuôi nghêu hoàn toàn sử dụng phương pháp thu hoạch thủ công, suất thấp chi phí lao động cao Ngoà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 hoàn chỉnh với đầy đủ tính kinh tế - kỹ thuật – công nghệ hoà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ọ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 toá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 hoà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ô ̣i dung nghiên cứu của đề 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ựng mô hình tổ ng quát LHM thu hoa ̣ch nghêu - Tính toán thông số bô ̣ phâ ̣n di đô ̣ng - Xây dựng mô hiǹ h bô ̣ phâ ̣n di đô ̣ng LHM + Mô hình xích cao su không biế n da ̣ng + Mô hiǹ h xić h cao su biế n da ̣ng - Xây dựng mô hiǹ h tính toán - Khảo sát ảnh hưởng các yế u tố kế t cấ u ̣ thố ng di đô ̣ng 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 hoa ̣ch vào vùng chăn nuôi nghêu thay thế cho viê ̣c thu hoa ̣ch thủ công hiê ̣n vẫn đảm bảo suấ t và chấ t lươ ̣ng 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 nuôi 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 nuôi 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, Duyên 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%, khoá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 tỉnh khu vực 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 đoàn nuôi nghêu Các HTXTS tập đoàn biến gần 16.000 bãi bồi thành khu nuôi nghêu, có khoảng 9650 khai thác Riêng các huyê ̣n Thạnh Phú, Bình Đại Ba Tri của tin̉ h Bế n Tre tập trung hầu hết HTX tập đoàn nuôi 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ê ̣n tích nuôi khoảng 1000 cho sản lươ ̣ng 5.000 tấ n với suấ t bình quân tấ n/ha Đế n năm 2006 đầ u năm 2007, diê ̣n tích nuôi nghêu 20.000 với tổ ng sản lươ ̣ng ở toàn bô ̣ khu vực Nam Bô ̣ đa ̣t 70 - 80 ngàn tấ n/năm và đã trở thành mă ̣t hàng xuấ t khẩ u lớn ( Lê Xuân Sinh, 2007 ) Khu vực phiá nam vùng có bờ biển dài gầ n 6000 Km với diê ̣n tić h bãi bồi rấ t rô ̣ng lớn Ven biển thuộc địa phận các tỉnh 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ề nuôi 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 nuôi trồng thủy sản đóng góp đáng kể vào phát triển kinh tế khu vực 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 phiá nam có nhiều tiềm phát triển nuôi trồng thủy sả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 các tin̉ h 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 nuôi nghêu chưa đạt kết Theo quyế t ̣nh phê duyê ̣t của Bộ Nông nghiệp Phát triển nông thôn; có Quy hoạch vùng nuôi 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ượng tới năm 2020 là 357.300 Để khai thác có hiệu tiềm diện tích đất vùng bãi bồi ven biển phiá nam vấn đề đặt phải có nghiên cứu tổng thể từ khâu quy hoa ̣ch baĩ nuôi, kỹ thuâ ̣t nuôi tới khâu thu hoa ̣ch 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 Kỹ thuâ ̣t nuôi nghêu: Kỹ thuật nuôi: Cho ̣n baĩ nuôi: Cho ̣n baĩ triề u gầ n cửa sông, it́ sóng gió, bằ ng phẳ ng Khi triề u xuố ng thấ p có cả phơi baĩ vài giờ/ngày Chấ t đáy nên cho ̣n là cát bùn Không cho ̣n nơi có nguy bi ̣ nước thải ô nhiễm tấ n công vì nghêu rấ t nha ̣y 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 ̣n đươ ̣c baĩ nuôi tiế n hành vê ̣ sinh mă ̣t baĩ , ̣n ta ̣p chấ t, cối ở nề n đáy Nế n nề n đấ t cứng phải xới lên cho tơi xố p, mề m giúp nghêu vùi miǹ h tố t Dùng đăng tre hoă ̣c 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 co ̣c cắ m cho đăng lưới đứng nghiêng về phiá trong, các co ̣c cách 1,2m tới 1,5m Trong baĩ , căng nhiề u dây ngang ̣c theo chiề u gió để ̣n chế nghiêng theo hướng nước chảy và gió lớn Thả giố ng: Thả làm ̣t vào tháng 2, tháng và tháng 7, tháng Đợt 1: Thả giố ng kić h cỡ 3.000-5000 con/1kg Sau tới tháng ương thì san thưa để nuôi lớn Đợt 2: Thả giố ng kích cỡ 500 - 800 con/kg Nên thả nuôi luôn, không cầ n san thưa, mâ ̣t đô ̣ 300-450 con/kg Chú ý thả lúc triề u lên để nghêu dễ vùi lớp bùn đáy Chăm sóc: Khi mới rải nghêu baĩ nuôi, vỏ nghêu còn yế u, chưa thể vùi sâu xuố ng đáy nên ̣n chế la ̣i baĩ nghêu Bắ t hế t các loa ̣i ố c mỡ trơn bãi triề u rồ i mới thả nghêu giố ng Trong quá trình nuôi thường xuyên kiể m tra bắ t ốc ̣i nghêu con, thường xuyên kiể m tra rào chắ n, chân rào để ngăn chă ̣n nghêu thoát ngoài Nế u nghêu tâ ̣p trung ở mô ̣t chỗ phải san gian ̃ vì mâ ̣t đô ̣ cao sẽ làm cho nghêu châ ̣m lớn Thu hoa ̣ch nghêu Nghêu đươ ̣c 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 nuôi - 10 tháng có thể thu hoa ̣ch nghêu Nên lựa cho ̣n cuố i mùa mưa để thu hoa ̣ch vì lúc này nghêu béo nhấ t Năng suấ t nghêu biǹ h quân đa ̣t 20 - 25 tấ n/ha Nế u người nuôi có kinh nghiê ̣m và kỹ thuâ ̣t cao có thể đa ̣t 50 - 55 tấ n/ha Hiê ̣n khâu thu hoa ̣ch ngêu chỉ là thủ công dùng cào tới hoă ̣c dùng tia nước áp lực để thu hoa ̣ch ngêu phí cho khâu thu hoa ̣ch là rấ t lớn đồ ng thời không đảm bảo đươ ̣c tin ́ h thời vu ̣ Xuất phát từ vấn đề trên: “Nghiên cứu ảnh hưởng số thông số đến khả di động LHM thu hoạch nghêu ” góp phầ n hoàn thiê ̣n LHM thu hoa ̣ch nghêu áp du ̣ng cho những vùng chăn nuôi nghêu thời gian tới vấn đề cần thiết có ý nghiã thực 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 của Nguyễn Hữu phu ̣ng ( 1996 ) cho thấy nghêu thường phân bố những vùng có đáy cát hay cát bùn đó tỷ lệ cát 68 - 75%, kích thước cỡ ̣t 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ó đô ̣ dố c tương đố i bằ ng phẳ ng Theo Nguyễn Tác An và Nguyễn Văn Lu ̣c ( 1994 ) cho rằ ng nghêu phân bố ở vùng thời gian phơi baĩ từ - giờ/ngày Đô ̣ sâu cực đa ̣i tim ̀ 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ươ ̣ng bùn lỏng và xác hữu ( 10 - 18 % ), vào mùa mưa bùn lỏng bao phủ nề n đáy baĩ nghêu ( 1,5 - 2,5 cm ) Theo nghiên cứu của Võ Si ̃ Tấ n ( 1999 ) ở Gò Công Đông cho thấ y nghêu tâ ̣p trung ở đô ̣ cao từ khoảng 0,8 - 1,5 m so với hải đồ Nghêu loài rộng nhiệt: Chúng sống điều kiện nhiệt độ từ 35oC, khoảng nhiệt độ thích hợp cho nghêu phát triển 28-31oC, 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 loà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, 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 thoát 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 10 cm Nghêu loà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 thoát 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 loài, vị trí địa lý phân bố, thời tiết, vùng trung triều hay hạ triều, khác 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 Ngoài chúng 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 trin ̀ h mô phỏng Chương trình tính toá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; % Chieu cao tam D= 140; % chieu dai tiep dat 140 B=35 ; % % 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; mf = 0.41 ; a = 0.1; %c0=0.100; c1=0.855; c2=0.4446; pm=W/(2*B*D); pm0=pm; k1=0.49147; k2=0.56107; n2=0.8001; mc=0.46387; a=0.2011; c0=0.1481; c1=0.8491; c2=0.3272; %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 e0 sf0 M=[ e e0 sr0 pf0 pr0 Teta'; ;sf0 sr0; pf0 pr0 ]; fprintf(' eD= %6.4f sf0=%6.4f pf0=%6.4f \n', e*D,sf0,pf0) fprintf(' e0D= %6.4f sr0=%6.4f pr0=%6.4f \n',e0*D,sr0,pr0) fprintf(' test= %6.4f fprintf(' tam1= %6.4f Teta= %6.4f \n',test, Teta0) tam2= %6.4f p0= %6.4f pm= %6.4f\n',tam1,tam2,p0,pm) fprintf(' KT= %6.4f \n',KT) fprintf(' thu1= %6.4f thu2=%6.4f \n',thu1,thu2) 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 TphaymXX(i) = MangtphaymX(N); 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); % -5 % -6 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)-(ld-0.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