(Luận văn) 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 lu an NGHIÊN CỨU ẢNH HƯỞNG CỦA MỘT SỐ THÔNG va n SỐ ĐẾN KHẢ NĂNG DI ĐỘNG CỦA LIÊN HỢP MÁY to p ie gh tn THU HOẠCH NGHÊU nl w CHUYÊN NGÀNH: KỸ THUẬT CƠ KHÍ d oa MÃ SỐ: 60520103 va an lu ll u nf LUẬN VĂN THẠC SỸ KỸ THUẬT CƠ KHÍ oi m z at nh NGƯỜI HƯỚNG DẪN KHOA HỌC z m co l gm @ TS BÙI VIỆT ĐỨC an Lu Hà Nội, 2014 n va ac th si 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ố lồi thủy hải sản, ni 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 lu Tiềm định hướng phát triển nghề nuôi nghêu thể rõ qua an n va định phê duyệt quy hoạch phát triển ni nhuyễn thể hàng hóa tập trung đến vùng ni nghêu, sị đến năm 2015 năm 2020 cho tỉnh đồng sông Cửu gh tn to năm 2020, Bộ Nông nghiệp Phát triển nông thơn; có Quy hoạch ie Long (ĐBSCL) như: Cà Mau, Bạc Liêu, Tiền Giang, Sóc Trăng, Kiên Giang, Bến p Tre Trà Vinh, với tổng diện tích đến năm 2015 28.110 năm 2020 phát nl w triển lên 35.690 d oa Quyết định Quy hoạch nêu rõ: đến năm 2015, diện tích nuôi nghêu 15.950 ha, va 206.300 an lu sản lượng 142.700 đến năm 2020, diện tích ni nghêu 20.590 ha, sản lượng m nuôi nghêu khâu thu hoạch ll u nf Một công việc tiêu tốn nhiều sức lao động trình oi Hiện phương pháp thu hoạch nghêu chủ yếu vùng chăn ni nghêu hồn z at nh 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 z hoạch l gm @ 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 m co 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 an Lu yêu cầu kỹ thuật - công nghệ, thay cho công việc thu hoạch thủ công n va ac th si 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ệ hồn chỉnh, đáp ứng tốt 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 lu 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 an n va vụ sản xuất cho vùng ni nghêu ven biển phía Nam Việt Nam Trường Đồng Nai Với lý trình bầy chúng tơi chọn thực đề tài: ie gh tn to Đại học Nông nghiệp Hà Nội Trường Trung cấp nghề Cơ Điện Đông Nam Bộ p “Nghiên cứu ảnh hưởng số thông số đến khả di động LHM nl w thu hoạch nghêu ” d oa Mục tiêu nghiên cứu va nghiên cứu là: an lu Xuất phát từ lý thực đề tài nêu đặt mục tiêu u nf Xác định ảnh hưởng số thông số đến khả di động LHM ll 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 - m oi 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 z at nh 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 l gm @ Nô ̣i dung nghiên cứu của đề tài z suất hiệu sử dụng máy m co Do thời gian nghiên cứu có hạn, đề tài giới hạn nội dung sau: - Tính toán thông số bô ̣ phâ ̣n di đô ̣ng an Lu - Xây dựng mô hình tổ ng quát LHM thu hoa ̣ch nghêu n va ac th si - 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ủ lu công hiê ̣n vẫn đảm bảo suấ t và chấ t lươ ̣ng sản phẩ m an n va Tính khoa học đề tài: tn to Đề tài phát triển vấn đề khoa học lý luận thực tiễn gh việc giới hóa khâu thu hoạch nghêu p ie Ý nghĩa thực tiễn đề tài: w Đề 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 d oa nl kinh tế, kỹ thuật, đảm bảo tính thời vụ khâu thu hoạch nghêu ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si 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 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 lu Do đó, nghêu coi đối tượng chủ lực cho phát triển nuôi biển an n va kỷ XXI Theo thống kê FAO (2006) tổng sản lượng động vật thân mềm 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 gh tn to 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 ie tỉnh miền Tây Nam Bộ như: Gị Cơng Đơng (Tiền Giang), Cầu Ngang, Duyên Hải p (Trà Vinh), Bình Đại, Ba Tri, Thạnh Phú (Bến Tre), Vĩnh Châu (Sóc Trăng); nl w nghêu Bến Tre người dân đưa vào nuôi thử nghiệm số vùng cửa sông d oa ven biển cho kết tốt an lu Nghêu đối tượng người dân quan tâm giá trị dinh dưỡng va Thịt thơm ngon có nhiều chất dinh dưỡng Prơtêin chiếm 15,66%, Lipit u nf chiếm 3,43%, khống chiếm 3-13% (Nguyễn Chính CTV) ll mặt hàng hải sản có giá trị xuất cao Đặc biệt tháng 10/2008 hội đồng m oi 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 z at nh (Marine Sterwarship Council) cho nghêu Bến Tre trở thành đặc sản biển z khu vực Đông Nam Á Sự công nhận giúp nghêu Bến Tre và các giới l gm @ tỉnh khu vực có nhiều hội xuất sang nhiều thị trường khu vực m co 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 an Lu n va ac th si 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 đồ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ê ̣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 lu hàng xuấ t khẩ u lớn ( Lê Xuân Sinh, 2007 ) an n va 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 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 gh tn to 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 ie hàng năm đưa biển lượng phù sa lớn, bồi thành bãi triề u rộng p tương đối phẳng, cấu tạo trầm tích baĩ bồi chủ yếu dạng cát, nl w sa, sét thích hợp cho ni nghêu ni trồng loại thủy sản khác d oa Với vị trí địa lý điều kiện tự nhiên thuận lợi, tiềm để phát triển kinh an lu tế thủy sản, nghề nuôi trồng thủy sản đó có nghề nuôi nghêu Những va năm gần với việc phát triển chung ngành kinh tế, nghề ni trồng u nf thủy sản có bước phát triển nhanh chóng diện tích, sản lượng giá trị ll 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 m oi tế khu vực góp phần cải thiện đời sống cho nhân dân xã ven biển z at nh Mặc dù vùng bãi triều ven biển phiá nam có nhiều tiềm phát triển ni trồng thủy sản, cịn phần lớn diện tích chưa khai z gm @ thác có hiệu quả; Đặc biệt nghề ni 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ã l m co 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 an Lu chế sách nên việc ni nghêu chưa đạt kết n va ac th si 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 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ượ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; lu đả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 an n va Kỹ thuâ ̣t nuôi nghêu: gh tn to 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 p ie 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 w 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 oa nl trường ô nhiễm d 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 an lu đáy cứng, nhiều bùn, dẫn đến chết nghêu u nf va Sau cho ̣n đươ ̣c baĩ nuôi tiế n hành vê ̣ sinh mă ̣t baĩ , ̣n ta ̣p chấ t, cối ở ll 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 oi m 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 z at nh - 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 z chế nghiêng theo hướng nước chảy và gió lớn @ l gm Thả giố ng: Thả làm ̣t vào tháng 2, tháng và tháng 7, tháng thưa để nuôi lớn m co Đợt 1: Thả giố ng kić h cỡ 3.000-5000 con/1kg Sau tới tháng ương thì san an Lu Đợt 2: Thả giố ng kích cỡ 500 - 800 con/kg Nên thả nuôi luôn, không cầ n san n va ac th si 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 lu Thu hoa ̣ch nghêu an Nghêu đươ ̣c thu hoạch vào lúc triều rút, lúc chúng ăn no thải va n 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 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 ie gh tn to lên, nghêu thường ngậm cát, bùn) p để thu hoa ̣ch vì lúc này nghêu béo nhấ t nl w 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 d oa và kỹ thuâ ̣t cao có thể đa ̣t 50 - 55 tấ n/ha an lu Hiê ̣n khâu thu hoa ̣ch ngêu chỉ là thủ công dùng cào tới hoă ̣c va 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 u nf đồ ng thời không đảm bảo đươ ̣c tin ́ h thời vu ̣ ll Xuất phát từ vấn đề trên: “Nghiên cứu ảnh hưởng số thông số m oi đến khả di động LHM thu hoạch nghêu ” góp phầ n hoàn thiê ̣n LHM thu z at nh 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 z gm @ 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 l m co 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ó an Lu đáy cát hay cát bùn đó tỷ lệ cát 68 - 75%, kích thước cỡ ̣t từ 0,006 - n va ac th si 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 ) lu Theo nghiên cứu của Võ Si ̃ Tấ n ( 1999 ) ở Gò Công Đông cho thấ y nghêu tâ ̣p an Nghêu lồi rộng nhiệt: Chúng sống điều kiện nhiệt độ từ - n va trung ở đô ̣ cao từ khoảng 0,8 - 1,5 m so với hải đồ gh tn to 35oC, khoảng nhiệt độ thích hợp cho nghêu phát triển 28-31oC, tốc độ dòng chảy ie 0,1 - 0,25 m/s, hàm lượng oxy hoà tan khoảng - mg/l, pH 6-9, độ mặn dao động p từ 19-26% độ mặn phù hợp cho nghêu phát triển 22-25% nl w Tập tính sống: Nghêu lồi sống đáy gặp điều kiện không thuận d oa lợi nhiệt độ giảm, độ mặn thay đổi đột ngột chúng thường lên khỏi đáy an lu di chuyển nơi khác có điều kiện thích hợp va Theo Trương Quốc Phú (1999) độ mặn giảm xuống 5% hầu hết nghêu u nf 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 ll 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 m oi 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 z at nh động vật thân mềm hai mảnh vỏ nói chung, trình sống trải qua giai đoạn z gm @ - 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ố l m co 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 an Lu khép vỏ Do giai đoạn cần đáy phẳng cần có vật bám cát to, sỏi n va ac th si - 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 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, lu an vùng trung triều hay hạ triều, khác cá thể mà di n va truyền tạo to tn Tốc độ tăng trưởng nghêu phụ thuộc vào vùng phân bố nhiều hay thức ie gh ăn Nghêu phân bố vùng cửa sông phong phú thành phần thực vật phù du p 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 w vùng triều cao Ngồi chúng cịn phụ thuộc vào nhiều nhân tố độ mặn, oa nl nhiệt độ, sóng gió Sự sinh trưởng thay đổi từ năm đến năm khác d khu vực mà có nhiệt độ biến đổi theo mùa Vào mùa xuân mùa hè nhiệt độ lu va an 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 u nf ll Trương Quốc Phú (1999) nghiên cứu đặc điểm sinh trưởng nghêu cho m oi thấy, nghêu sinh trưởng nhanh từ tháng - sinh trưởng chậm từ tháng 10-4 z at nh năm sau (tháng 10-11) nước từ sơng Cửu Long đổ mạnh làm độ mặn z xuống thấp, ngao vùi để tránh nước làm giảm thời gian bắt mồi gm @ Tốc độ sinh trưởng trung bình 1,72 mm/tháng 789,47 mg/tháng Tốc độ l sinh trưởng khối lượng tương đối nhanh so với tốc độ sinh trưởng chiều dài Ở m co 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, an Lu Tiền Giang ngao cỡ 25 mm nặng 2,8 g/con n va ac th si 81 Chương trin ̀ h mô phỏng 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: % % lu an W=10507; % Trong luong xe n va kw=W/W00; % He so tang luong de tinh lai Do lun hg=0.25; % Chieu cao tam D= 140; % chieu dai tiep dat 140 p ie gh tn to chW=['(Khi G=', num2str(W), 'N)']; B=25; Rf=6.5; % ban kinh banh truoc 25 Rk=13; d oa nl w % % be rong xich B=35 ; Rr=9 ; % Ban kinh banh sau u nf va an % Chieu cao mau bam cm ll oi m H=3 ; lu Rp=D/5; %% Ban kinh banh sa0 chu dong 25 z at nh H0=6000; % [N]= Luc cang xich Ban dau; H0=6000 Gp=9; % Buoc vau xich z % 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 m co l gm @ e=0.02; an Lu Vphay= 1.00; % m/s Ttoc vong cua banh chu dong n va ac th si 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; lu an %mf=0.623; n va mf=0.3630; % n1=0.8664 k3=30.3; k4=30.3; n4=0.632; p ie gh tn to n1=0.8664*0.5; % nl w teta0=0.000; d oa tam1=1000; va i=0; an lu tam2=0; dT=0.001; % Buoc bien thien ll m i=i+1; u nf while (abs(tam1-tam2)>50) oi % - z at nh % sf0=4.5556; ds=0.001; sm=5.5; %G=18000 p0= k1*H^n1; z an Lu %pf0=k1*(sf0)^n1; m co e0=e+(hg+H)*tan(teta0)/D; l teta0= atan((sr0-sf0)/D); gm sr0=sm+i*ds; @ %for i=1: 1540 n va ac th si 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); lu an tam2=2*B*D*p0+(2*k2*B*D/(n2+1))*((sr0-H)^(n2+1)-(sf0- n va H)^(n2+1))/(sr0-sf0); to test = tam1-tam2; thu1=W*cos(teta0)*(0.5+e0)*D; p ie gh tn end thu2=((2*B*D^2)/(sr0-sf0)^2)*((k2/(n2+2)*((sr0-H)^(n2+2)-(sf0-H)^(n2+2)))- nl w (k2/(n2+1)*(sf0-H)*((sr0-H)^(n2+1)-(sf0-H)^(n2+1)))+(0.5*p0*(sr0^2-sf0^2)- d oa p0*sf0*(sr0-sf0))); e0 sf0 M=[ e e0 sr0 va Mc=' e an lu KT=thu1-thu2; pf0 pr0 Teta'; pr0=%6.4f \n',e0*D,sr0,pr0) z at nh fprintf(' tam1= %6.4f sr0=%6.4f oi fprintf(' test= %6.4f pf0=%6.4f \n', e*D,sf0,pf0) m fprintf(' e0D= %6.4f sf0=%6.4f ll fprintf(' eD= %6.4f u nf ;sf0 sr0; pf0 pr0 ]; Teta= %6.4f \n',test, Teta0) tam2= %6.4f p0= %6.4f z \n',KT) x=0:0.01:D; an Lu s0x= sf0 + (sr0-sf0)*x./D; m co fprintf(' thu1= %6.4f thu2=%6.4f \n',thu1,thu2) l gm fprintf(' KT= %6.4f @ pm= %6.4f\n',tam1,tam2,p0,pm) n va ac th si 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; lu an id= id0/100; n va tetai=teta0; jB=id*D/(1-id); gh tn to idphay=1-(1-id)/cos(tetati); %sO LIEU CU (dX=0.1) p ie dX=0.05; N= round(D/dX); TPtmX=0; s0iX = sf0+(sr0-sf0)*(X/D); oi m X= i2*dX; XV(i2)=X; ll u nf for i2=1:N va an MangtmX=[0]; lu MangTP=[0]; d oa nl w TPpiX=0; z at nh piX = k1*H.^n1 + k2*(s0iX - H).^n2; PIX1(i2)=piX; jx = jB + idphay*X; z TmX(i2)=2*B*tmX; l gm @ tmX = (mc + mf*piX)*(1-exp(-a*(jB + idphay*X))); %Luu vao mang an Lu TPpiX=TPpiX+piX*dX; m co TmS(i2)= 4*H*(mc + mf*piX/pi*acot(H/B))*(1-exp(-a*(jB+idphay*X))); n va ac th si 85 TPtmX=TPtmX+tmX*dX; MangTP=[MangTP TPpiX]; MangtmX=[MangtmX TPtmX]; end MangHieuTP=[0]; MangTX=[0]; MangKQ=[0]; MangKQ2=[0]; MangS0ix=[3]; lu an MphayiX=[0]; n va MT0X=[0]; XM1=round(0.492*Rp/dX); gh tn to MtphaymX=[0]; p ie Xm1=XM1*dX; Mangpphay=[0]; nl w tpphayix=0; ll for ih=2:N z at nh MangTpphaymS=[0]; oi TpphaymS=0; m MangtphaymX=[0]; u nf TpphaymX=0; va an Mangtam=[0]; lu tpiX2=0; d oa MangX2=[0]; z HieuTP=B*(MangTP(XM1)-MangTP(ih)); an Lu MangTX=[MangTX TX]; m co TX=H0+B*MangtmX(ih); l MangHieuTP=[MangHieuTP HieuTP]; gm @ Xh=ih*dX; XVV(ih)=X; n va ac th si 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 lu an KQ=(Xm1-Xh)*(MangHieuTP(ih)/MangTX(ih)); n va end to MangKQ=[MangKQ KQ]; p ie gh tn %Luu KQ phep chia vao mang KQ if (ih round(Rp/dX) ll m end u nf tam=ih; va if mod((ih-1),round(Rp/dX))==0 oi s0iXm1 = (((sr0-sf0)*Xm1)/D+sf0)-MangKQ(tam); z at nh end z Thay doi gia tri XM sau moi mot vong+0.0124999999*Rp/dX if mod(ih,round(Rp/dX))==0 an Lu end m co Xm1=XM1*dX; l XM1=round(XM1+Rp/dX); gm @ % n va ac th si 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)); lu an TpphaymX =real(TpphaymX + tphaymX*dX); n va MangtphaymX=[MangtphaymX TpphaymX]; to MangtphaymX(ih)=2*B*MangtphaymX(ih); TphaymS= 4*H*(mc + mf*pphayiX/pi*acot(H/B))*(1-exp(-a*(jf))); p ie gh tn MtphaymX=[MtphaymX TpphaymX]; TpphaymS = real(TpphaymS+TphaymS*dX); nl w MangTpphaymS=[MangTpphaymS TpphaymS]; an lu end d oa T0X(ih)=H0+B*MtphaymX(ih); va MT0X=[MT0X T0X]; %N0=1400 u nf TphaymXX(i) = MangtphaymX(N); ll TphaymSS(i) = MangTpphaymS(N); z at nh - oi % m T0Xx(i)=MT0X(N) Tmbb=trapz(XV,TmX);% -1 z an Lu Tmb(i)=real(Tmbb); m co id0V(i)=id0; l XX(i2)=X; gm TPtmXX(i2)=real(Tmbb); @ Tmss=trapz(XV,TmS);% -2 n va ac th si 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); lu an for i3= 1:M n va teta= i3*dT; tetamV(i3)= teta; %luu lai to sf0iteta = (Rf+H)*(cos(teta+tetat0i)-cos(tetaf0-tetat0i))/cos(teta+tetat0i); pfiteta=k1*(sf0iteta)^n1; p ie gh tn %pfiteta = k1*H^n1 + k2*(sf0iteta -H)^n2; pfitetaV(i3)= pfiteta; d oa nl w jfteta = (Rf+H)*((tetaf0 - teta)-(1-id)*(sin(tetaf0+teta)- sin(teta + tetat0i))); an lu tfiteta = (mc + mf*pfiteta)*(1-exp(-a*jfteta)); va tfitetatphan(i3) = tfiteta*cos(teta) - pfiteta*sin(teta); u nf tfsteta = 4*H*(mc+mf*pfiteta/pi*acot(H/B))*(1-exp(-a*jfteta)); ll tfstetatphan(i3)= tfsteta*cos(teta); oi z at nh end m tetatphan(i3)= teta; Tfbb = 2*B*(Rf+H)*trapz(tetatphan,tfitetatphan); % -3 z Tfb(i)=real(Tfbb); an Lu tetar0i= teta0; %arctan((sr0-sf0)/D); m co % - l Tfs(i)=real(Tfss); % -4 gm @ Tfss=(Rf+H)*trapz(tetatphan,tfstetatphan); n va ac th si 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)); lu an pridelta = k1*H^n1 + k2*(sr0idelta-H)^n2; n va trdelta = (mc+mf*pridelta)*(1-exp(-a*jrdelta)); to trdeltatphan(i4)= trdelta*cos(tetat0i-delta) + pridelta*sin(tetat0i-delta); trsdeltatphan(i4)=trsdelta*cos(tetar0i - delta); p ie gh tn trsdelta= 4*H*(mc+mf*pridelta*acot(H/B)/pi)*(1-exp(-a*jrdelta)); end % -5 nl w Trbb= 2*B*(Rr+H)*trapz(deltatphan,trdeltatphan); % -6 d oa Trss= (Rr+H)*trapz(deltatphan,trsdeltatphan); va Trs(i)=real(Trss); an lu Trb(i)=real(Trbb); ll u nf T3(i)= Tmb(i)+Tms(i)+Tfb(i)+ Tfs(i)+Trb(i)+Trs(i); oi m z at nh %============================ z gm @ tetat0i= teta0; %arctan((sr0-sf0)/D); tetaf0 =acos(cos(tetat0i)-sf0*cos(tetat0i)/(Rf+H))-tetat0i; an Lu M= round(tetaf0/dT); m co Tongf=0; dT=0.01; l jfs= (Rr +H)*sin(tetaf0)*id/(1-id); n va ac th si 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)); lu an end n va sfs= Tongf; idphay= 1- (1-id)/cos(teta0); gh tn to % p ie js= idphay*D/(1-idphay); Tongr=0; %dX=0.1; % CHUA DUNG nl w for n=1:N N= round(D/dX); d oa X= n*D/N; an lu s0iX = sf0+(sr0-sf0)*(X/D); va piX = k1*H.^n1 + k2*(s0iX - H).^n2; PIX(n)=piX; oi m srs = sfs + Tongr; ll end u nf Tongr = Tongr + c0*(piX*n*D/N)^c1 *((n*js/N)^c2 - ((n-1)*js/N)^c2); z at nh sri= (sr0 + srs)*cos(teta0); T2(i)= 2*k1*B*H^(n1+1)/(n1+1) + 2*k1*B*H^n1*(sri - H) + z m co l gm %T2(i)=real(T2); @ 2*k2*B/(n2+2)*(sri-H)^(n2+1); %========================================== an Lu n va ac th si 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))- lu an W*(hg*sin(tetati)-D*(0.5-e)*cos(tetati)))/(T2(i)+W*cos(tetati)/ n va tan(tetati)); p ie gh tn to % ========================================== E1(i)= T3(i)*Vt/1000; % kW - Cong suat chu dong nl w E2(i)= T2(i)*V/1000 ; % kW - Cong suat chi phi cho luc can d oa E3(i)= (T3(i)*(1/(1-id) - 1/cos(teta0))*V + W*V*tan(teta0))/1000; an lu if T4