BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT VIỆN CƠ ĐIỆN NÔNG NGHIỆP VÀ CÔNG NGHỆ SAU TH HOẠCH ************************************** NGUYỄN ĐỨC THẬT NGHIÊN CỨU MỘT SỐ THÔNG SỐ CHÍNH CỦA BỘ PHẬN TÁC[.]
BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT VIỆN CƠ ĐIỆN NÔNG NGHIỆP VÀ CÔNG NGHỆ SAU TH HOẠCH ************************************** NGUYỄN ĐỨC THẬT NGHIÊN CỨU MỘT SỐ THƠNG SỐ CHÍNH CỦA BỘ PHẬN TÁCH NGỌN VÀ LÁ ỨNG DỤNG TRÊN MÁY LIÊN HỢP THU HOẠCH MÍA NGUYÊN CÂY LUẬN ÁN TIẾN SĨ CHUYÊN NGÀNH: KỸ THUẬT CƠ KHÍ MÃ SỐ: 952 01 03 NGƯỜI HƯỚNG DẪN KHOA HỌC: TSKH BẠCH QUỐC KHANG PGS.TS LƯƠNG VĂN VƯỢT Hà Nội – 2023 i LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi, kết nghiên cứu trình bày luận án trung thực, khách quan chưa dùng bảo vệ để lấy học vị Tôi xin cam đoan giúp đỡ trình thực luận án cám ơn, thơng tin trích dẫn luận án rõ nguồn gốc Hà Nội, ngày tháng năm 2023 Tác giả luận án Nguyễn Đức Thật ii LỜI CẢM ƠN Tôi xin chân thành cảm ơn tập thể thầy: TSKH Bạch Quốc Khang – Viện Cơ điện nông nghiệp Công nghệ sau thu hoạch; PGS.TS Lương Văn Vượt – Trường Đại học Chu Văn An; TS Đậu Thế Nhu – Viện Cơ điện nông nghiệp Công nghệ sau thu hoạch, hướng dẫn, giúp đỡ tơi hồn thành luận án Đồng thời xin chân thành cảm ơn lãnh đạo Viện Cơ điện nông nghiệp Công nghệ sau thu hoạch, tập thể cán chuyên gia Viện đơn vị, cá nhân quan tạo điều kiện cho tơi q trình thực luận án Tôi xin trân trọng cảm ơn nhà khoa học, bạn đồng nghiệp quan đặc biệt thành viên Gia đình giúp đỡ, ủng hộ động viên để tơi hồn thành luận án Xin trân trọng cảm ơn! Tác giả luận án Nguyễn Đức Thật iii MỤC LỤC Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục ký hiệu chữ viết tắt viii Danh mục bảng xii Danh mục hình xiv MỞ ĐẦU 1 Tính cấp thiết đề tài nghiên cứu Mục tiêu nghiên cứu 3 Đối tượng phạm vi nghiên cứu 3.1 Đối tượng nghiên cứu 3.2 Phạm vi nghiên cứu Nội dung nghiên cứu Ý nghĩa khoa học thực tiễn luận án 5.1 Ý nghĩa khoa học 5.2 Ý nghĩa thực tiễn Những đóng góp luận án Cấu trúc luận án Chương TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU 1.1 Đặc điểm sinh học mía 1.1.1 Thân mía 1.1.2 Ngọn mía 1.1.3 Lá mía 1.2 Kết nghiên cứu đặc tính cơ, lý mía 10 1.2.1 Thế giới 10 1.2.2 Việt Nam 14 iv 1.3 Tổng quan cơng trình nghiên cứu tách mía 15 1.3.1 Tách mía trước mía vào phận làm 16 1.3.2 Tách mía phận làm máy liên hợp THM 17 1.4 Kết nghiên cứu tách máy liên hợp thu hoạch mía 21 1.4.1 Cấu tạo phận tách máy liên hợp THM 21 1.4.2 Kết nghiên cứu tách mía giới 23 1.4.2.1 Một số loại kết cấu cánh lơ bóc 23 1.4.2.2 Ảnh hưởng dạng vật liệu làm bóc 24 1.4.2.3 Các nghiên cứu chế độ làm việc phận tách mía 26 1.4.3 Kết nghiên cứu tách mía Việt Nam 26 1.4.3.1 Kết nghiên cứu kết cấu cánh lơ bóc 26 1.4.3.2 Kết nghiên cứu chế độ làm việc lô 27 1.5 Đề xuất nguyên lý phận tách mía 30 1.5.1 Đề xuất nguyên lý phận bẻ 30 1.5.2 Sơ đồ nguyên lý phận tách mía 32 1.6 Các thông số ảnh hưởng tiêu đánh giá phận tách 32 mía 1.6.1 Các thông số ảnh hưởng 32 1.6.2 Các tiêu đánh giá phận tách mía 33 Kết luận chương 34 Chương PHƯƠNG PHÁP NGHIÊN CỨU 35 2.1 Phương pháp nghiên cứu lý thuyết 35 2.1.1 Phương pháp giải tích 36 2.1.2 Phương pháp số 36 2.2 Phương pháp nghiên cứu thực nghiệm 36 2.2.1 Cấu tạo mơ hình thí nghiệm 36 2.2.1.1 Cấu tạo, nguyên lý hoạt động mô đun tách 37 v 2.2.1.2 Cấu tạo, nguyên lý hoạt động mô đun bẻ 38 2.2.1.3 Sơ đồ truyền động phận tách ngọn, mía 39 2.2.2 Phân tích, lựa chọn thơng số cách điều chỉnh 40 2.2.2.1 Phân tích, lựa chọn thơng số 40 2.2.2.2 Cách điều chỉnh thơng số 42 2.2.3 Bố trí thí nghiệm 46 2.2.4 Phương pháp xác định tiêu 46 2.2.4.1 Tỷ lệ tạp chất 47 2.2.4.2 Tỷ lệ tổn thất 47 2.2.4.3 Chi phí lượng riêng 48 2.2.5 Thiết bị thí nghiệm dụng cụ đo đạc 48 2.2.5.1 Thiết bị dụng cụ thí nghiệm 48 2.2.5.2 Các bước tiến hành thí nghiệm 49 2.2.6 Phương pháp nghiên cứu thực nghiệm đơn yếu tố 50 2.2.7 Phương pháp quy hoạch thực nghiệm đa yếu tố 52 2.2.8 Phương pháp phân tích hồi quy 54 2.2.8.1 Đánh giá phương sai 56 2.2.8.2 Xác định dạng tổng quát mô hình hồi quy 56 2.2.8.3 Xác định hệ số mơ hình hồi quy 56 2.2.8.4 Phương pháp kiểm tra tính tương thích mơ hình tốn 56 2.2.8.5 Kiểm tra mức ý nghĩa hệ số hồi quy 57 2.2.8.6 Xác minh điều kiện tiên phân tích hồi quy 58 2.2.9 Phương pháp tối ưu hóa hàm tiêu 58 2.2.9.1 Tìm tọa độ điểm đặc biệt bề mặt hàm tiêu 58 2.2.9.2 Chuyển phương trình hồi quy sang dạng tắc 59 2.2.9.3 Phương pháp giải tốn thương lượng tìm cực trị có điều kiện 2.3 Phương pháp thu thập thông tin 60 61 vi 2.3.1 Phương pháp thu thập thông tin thứ cấp 61 2.3.2 Phương pháp thu thập thông tin qua nguồn sơ cấp 62 Kết luận chương 63 Chương NGHIÊN CỨU CƠ SỞ LÝ THUYẾT XÁC ĐỊNH MỘT 64 SỐ THÔNG SỐ CỦA BỘ PHẬN TÁCH NGỌN VÀ LÁ MÍA 3.1 Q trình tách mía bóc 64 3.2 Ảnh hưởng yếu tố động học đến q trình tách mía 66 bóc 3.2.1 Nghiên cứu chiều dài quét bóc 66 3.2.1.1 Xác định chiều dài quét bóc 66 3.2.1.2 Ảnh hưởng tỷ số chiều cao đặt lơ bóc đến chiều dài qt 68 3.2.1.3 Khảo sát ảnh hưởng hệ số động học đến chiều dài quét 69 3.2.2 Nghiên cứu hệ số quét lặp 69 3.2.2.1 Xác định hệ số quét lặp 69 3.2.2.2 Ảnh hưởng tỷ số chiều cao đặt lơ bóc đến hệ số qt lặp 70 3.2.2.3 Khảo sát ảnh hưởng số động học đến hệ số quét lặp 71 3.3 Nghiên cứu động lực học trình bẻ 72 3.3.1 Một số nhận xét q trình bẻ dàn thí nghiệm 72 3.3.2 Mơ hình nghiên cứu q trình bẻ mía 73 3.3.2.1 Phân tích q trình bẻ mía mô đun bẻ 73 3.3.2.2 Các giả thiết xây dựng mơ hình 74 3.3.2.3 Mơ hình biến dạng mía 75 3.3.2.4 Điều kiện bẻ gãy mía 84 3.4 Khảo sát mơ q trình bẻ mía 85 3.4.1 Thuật tốn mơ 85 3.4.2 Khảo sát trình động 88 3.4.3 Khảo sát trình tĩnh 89 vii 3.4.3.1 Phương trình mơ tả biến dạng tĩnh 89 3.4.3.2 Mơ men uốn mía theo chiều dài góc quay khác 90 bẻ 3.4.3.3 Ứng suất uốn mía theo chiều dài góc quay khác 91 bẻ 3.4.3.4 Ứng suất uốn cực đại mía theo góc quay bẻ 92 3.4.3.5 Khảo sát ảnh hưởng thơng số hình học 3.5 Xây dựng mơ hình khảo sát q trình bay mía sau khỏi 93 95 buồng làm 3.5.1 Xây dựng mô hình chuyển động (bay) thân mía 96 3.5.2 Mơ quỹ đạo chuyển động thân mía 99 3.5.3 Khảo sát ảnh hưởng vận tốc ban đầu đến quỹ đạo chuyển động 100 thân mía Kết luận chương 101 Chương KẾT QUẢ NGHIÊN CỨU THỰC NGHIỆM 103 4.1 Cơ sở chọn thông số nghiên cứu xác định miền biến thiên 103 4.1.1 Xác định thông số đầu vào tiêu chất lượng 103 4.1.2 Xác định miền biến thiên thông số đầu vào 105 4.2 Kết nghiên cứu thực nghiệm đơn yếu tố 107 4.2.1 Ảnh hưởng vận tốc cấp (Vk) tới tiêu 107 4.2.2 Ảnh hưởng số động học (λ) tới tiêu 112 4.2.3 Ảnh hưởng mật độ mía (q) tới tiêu 116 4.2.4 Ảnh hưởng khoảng cách tâm lô bẻ (Lm) đến tiêu 121 4.3 Xác định giá trị thông số tối ưu phương pháp quy hoạch hóa 126 thực nghiệm (QHTN) đa yếu tố 4.3.1 Mã hóa thơng số nghiên cứu lập ma trận thí nghiệm 127 viii 4.3.2 Kết xử lý số liệu thí nghiệm hàm tỷ lệ tạp chất 128 4.3.3 Kết xử lý số liệu thí nghiệm hàm tỷ lệ tổn thất 132 4.3.4 Kết xử lý số liệu thí nghiệm hàm chi phí lượng riêng 135 4.4 Tối ưu q trình tách ngọn, mía 138 4.5 Kết ứng dụng máy liên hợp THM nguyên SHC-0,2A 142 Kết luận chương 144 KẾT LUẬN VÀ ĐỀ NGHỊ 146 Danh mục cơng trình cơng bố liên quan đến Luận án 148 Tài liệu tham khảo 149 Phụ lục 160 ix DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT Ký hiệu, Chữ viết tắt Diễn giải Đơn vị Đạo hàm riêng góc theo thời gian J Đạo hàm riêng mơ men qn tính theo thời gian Đạo hàm riêng mô men uốn theo thời gian Đạo hàm riêng lực dọc theo thời gian Q Đạo hàm riêng lực cắt theo thời gian S Đạo hàm riêng phân tố S theo thời gian A Diện tích mặt cắt phân tố ax Gia tốc theo phương x m/s2 ay Gia tốc theo phương y m/s2 m2 b c CGH 1/2 chiều dài quét 1/2 khoảng cách hai trục lơ bóc Cơ giới hóa m m Chiều dài đường quét bóc m D2 dc d Đường kính lị xo Đường kính Góc ơm dây cung ds m m Độ dM Số gia mô men uốn mặt cắt phân tố Nm dm Khối lượng phân tố ds g dN Số gia lực dọc trục mặt cắt phân tố Áp lực đỉnh bẻ tác dụng lên chiều dài phân tố N dQ Số gia lực cắt mặt cắt phân tố N ds Chiều dài phân tố mía xét m E Es Mơ đun đàn hồi thép lị xo Mơ đun đàn hồi dóng mía Fgη Lực cản khí động tác dụng lên khối mía theo chiều η dP N MPa MPa N No 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Y1 17.72 14.42 9.01 3.49 29.14 24.26 17.89 11.56 14.9 11.57 7.12 2.69 24.05 19.99 14.02 9.29 15.56 11.35 19.33 8.86 9.19 18.6 13.07 10.16 Y2 17.86 14.01 9.1 3.49 29.22 23.95 17.45 12.04 15.17 12.11 6.84 2.62 24.16 20.37 14.26 8.97 15.61 11.59 19.67 9.75 9.3 18.17 13.63 10.59 Y3 17.35 14.29 8.59 3.66 29.58 23.36 17.85 11.12 14.73 11.73 6.58 2.72 24.11 20.04 14.79 9.79 15.61 11.31 19.29 9.71 9.73 17.89 13.56 10.47 Ytb 17.643 14.240 8.900 3.547 29.313 23.857 17.730 11.573 14.933 11.803 6.847 2.677 24.107 20.133 14.357 9.350 15.593 11.417 19.430 9.440 9.407 18.220 13.420 10.407 Ytt 17.954 14.064 8.703 3.632 29.035 24.011 17.847 11.642 14.710 11.842 6.848 2.800 24.177 20.175 14.377 9.195 15.773 11.237 19.493 9.377 9.444 18.182 13.336 10.491 Ytt - Ytb 0.311 -0.176 -0.197 0.085 -0.278 0.154 0.117 0.068 -0.224 0.039 0.001 0.123 0.070 0.042 0.021 -0.155 0.180 -0.180 0.063 -0.063 0.038 -0.038 -0.084 0.084 Tam cua mat quy hoach X = ( -2.586 , -1.787 , 1.810 , -0.253 ) Ytam = 1.49066 cac he so chinh tac W 0.3239 0.7713 1.4687 -1.2019 vec to rieng U ( A=U.W.UT) -0.8662 -0.4875 -0.0799 -0.0756 -0.2213 0.2488 0.9408 -0.0631 -0.4463 0.8314 -0.3188 0.0897 -0.0397 -0.0966 0.0826 0.9911 182 Si 0.06943 0.04390 0.07410 0.00963 0.05493 0.20903 0.05920 0.21173 0.04923 0.07693 0.07293 0.00263 0.00303 0.04263 0.15523 0.17080 0.00083 0.02293 0.04360 0.25270 0.08143 0.12790 0.09310 0.04923 Kết xử lý số liệu hàm chi phí lượng riêng DANH GIA DONG NHAT PHUONG SAI tieu chuan kohren G = 0.2491 he so tu m = 24 he so tu n-1= tieu chuan tra bangk ( 5%) G = 0.2871 KET QUA XU LY SO LIEU HOI QUY he so b0,0 b1,0 b1,1 b2,0 b2,1 b2,2 b3,0 b3,1 b3,2 b3,3 b4,0 b4,1 b4,2 b4,3 b4,4 Gia tri he so 2102.0028 225.9500 74.5472 308.7648 19.7958 113.5472 -289.3389 29.1000 34.7958 85.1806 -21.6694 30.4542 29.7000 24.8958 71.3289 Tieu chuan T student 126.5458 27.6273 3.4148 37.7532 2.2820 5.2013 -35.3780 3.3546 4.0112 3.9019 -2.6496 3.5107 3.4238 2.8700 3.2674 Phuong sai luong (lap) Sb = 1203.98125 So bac tu kb = 48 Phuong sai tuong thich Sa = 1645.61498 So bac tu ka = Tieu chuan FISHER F = 1.3668 Chuan Ro Ro = 0.9959 183 No 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Y1 2812.6 2261.5 2052.8 1539.9 3314.8 2697.6 2545.5 2211.9 2664.1 2209.4 2067 1658 3147.9 2821.8 2638.8 2367.9 2397.1 1959.9 2553.7 1907.1 1860.2 2501.7 2092.6 2198.1 Y2 2846.9 2239.9 2061 1506.6 3247.8 2736.4 2590.7 2196 2679.2 2253.2 2046.3 1567.4 3195.5 2882.9 2745 2348 2444.2 1914.3 2528.9 1872.8 1918.7 2479.4 2240.2 2212.5 Y3 2843 2210.1 2053.2 1524.4 3217.8 2729 2572.8 2213.5 2742.8 2242.4 2060.3 1656.8 3253 2844.9 2680.6 2437.5 2427.8 1916 2524.4 1906.4 1889.9 2473.2 2093.9 2202.6 Ytb 2834.2 2237.2 2055.7 1523.7 3260.2 2721 2569.7 2207.2 2695.4 2235 2057.9 1627.4 3198.8 2849.9 2688.1 2384.5 2423 1930.1 2535.7 1895.4 1889.6 2484.8 2142.3 2204.4 Ytt 2839.055 2228.455 2052.942 1521.525 3240.149 2745.949 2593.220 2178.203 2712.294 2223.510 2044.981 1635.381 3212.971 2840.588 2684.842 2391.642 2402.500 1950.600 2524.315 1906.785 1897.844 2476.522 2151.662 2195.001 Ytt - Ytb 4.848 -8.752 -2.765 -2.148 -20.024 24.909 23.513 -28.970 16.937 -11.480 -12.876 7.991 14.181 -9.269 -3.281 7.185 -20.533 20.533 -11.352 11.352 8.244 -8.244 9.389 -9.389 Tam cua mat quy hoach X = ( -1.844 , -1.610 , 2.271 , 0.484) Ytam = 1311.22172 cac he so chinh tac 134.6720 85.3044 57.1276 67.4998 vec to rieng U ( A=U.W.UT) -0.3206 0.5068 0.6536 -0.4616 -0.7606 -0.6402 0.1041 -0.0273 -0.4480 0.4887 -0.0725 0.7451 -0.3434 0.3072 -0.7461 -0.4806 184 Si 352.64333 666.09333 21.37333 277.66333 2466.33333 424.36000 518.12333 93.60333 1744.44333 520.68000 111.56333 2700.36000 2769.67000 951.80333 2862.17333 2208.40333 571.64333 668.24333 248.96333 384.32333 855.63000 224.66333 7198.52333 54.27000 Phụ lục Chương trình khảo sát, mơ Bài tốn khảo sát biến dạng uốn động mía clc; clear; global R dth Ia As ro Es Ln Nl hs smesh ps P EI miu Ep Ev ; global Rd hd xd nq omeg alfa0 R=0.1; dth=0.03; Ia=pi*dth^4/64; As=pi*dth^2/4; ro=1.2e3; miu=As*ro; Es=1e9; EI=Es*Ia; Ep=2e8; Ev=1e5; %giam chan Nl=200; Ln=0.3; hs=Ln/Nl; ps=100000; %================= thong so lo danh la Rd=0.2; hd=-0.05+Rd; xd=0.35; nq=800; omeg=pi*nq/30; % pt canh dap (xs-xd)^2+(ys-hd)^2=Rd^2 % xs=xd-Rd*cos(alfa); ys=hd-Rd*sin(alfa) alfa0 ys=0 alfa0=asin(hd/Rd); %====================================== smesh=linspace(0,Ln,Nl); %dung cho ve thi theo s %dien kien t=0 y1=zeros(1,2*Nl); dy_dt=zeros(1,2*Nl); P=zeros(1,Nl); %for i=15:20 % P(i)=ps; %end %for i=Nl-10:Nl-5 % P(i)=ps*(i-Nl+10)/5; %end Nt=301; tend=0.01; t_interval=linspace(0,tend,Nt); 185 %time interval and initial conditions y0=zeros(1,2*Nl); opts = odeset('RelTol',1e-4,'AbsTol',1e-6); %solution [t,y] = ode45(@p2dydt,t_interval, y0,opts); ys=y(:,Nl+1:2*Nl); for k=1:Nt xs(k,1)=(hs^2-ys(k,1)^2)^0.5; for i=2:Nl xs(k,i)=xs(k,i-1)+(hs^2-(ys(k,i)-ys(k,i-1))^2)^0.5; end end for k=1:Nt dy_ds2(k,1)=(ys(k,2)-2*ys(k,1))/hs^2; for i=2:Nl-1 dy_ds2(k,i)=(ys(k,i+1)-2*ys(k,i)+ys(k,i-1))/hs^2; end dy_ds2(k,Nl)=0; end Muon=dy_ds2*Es*Ia; sigma=Muon/Ia*dth/2; figure(1) yb=y(:,2*Nl); plot(t_interval*1000,yb); grid on figure(2) for i=1:10:Nt plot(smesh,y(i,Nl+1:2*Nl)); hold on end grid on % figure(3) % for i=1:10:Nt % plot(xs(i,:),ys(i,1:Nl)); % hold on %end % grid on figure(4) for i=1:20:Nt plot(smesh,Muon(i,1:Nl)); hold on end grid on function dydt = p1dydt(t,y) %fi dfidt dung o cac node 1->Nl global R dth Ia As ro Es Ln Nl hs smesh ps P EI miu Ep Ev; global Rd hd xd nq omeg alfa0 dydt= zeros(2*Nl,1); y1=y(Nl+1:2*Nl)'; yv1=y(1:Nl); 186 dy=y(Nl+1:2*Nl); ys=[0 y1 0]; % dieu kien bien ================= ys(1)=ys(3); ys(Nl+3)=2*ys(Nl+2)-ys(Nl+1); ys(Nl+4)=4*ys(Nl+2)-4*ys(Nl+1)+ys(Nl); %================================ % tinh toa be========================== alfa=alfa0-omeg*t; xt=xd-Rd*cos(alfa); yt=hd-Rd*sin(alfa); ytv=Rd*omeg*cos(alfa); %van toc P= zeros(1,Nl); kd=round(xt/hs); P(kd)=Ep*(yt-y1(kd))+Ev*(ytv-yv1(kd)); P(kd+1)=P(kd); P(kd-1)=P(kd); P(kd+2)=P(kd); P(kd-2)=P(kd); %lec=yt-y1(kd) %lecv=ytv-yv1(kd) %y(1->Nl)=y1 ; y(Nl+1->2*n)=dy; % Tinh dfi/ds, dfi^2/ds^2 o cac node for i=1:Nl %dyds(i)= (ys(i+3)-fi(i+1))/2/hs; %dyds2(i)=(ys(i+1)-2*ys(i+2)+ys(i+3))/hs^2; %dyds3(i)=(ys(i+4)-2*ys(i+3)+2*ys(i+1)-ys(i))/2/hs^3 dyds4(i)=(ys(i+4)-4*ys(i+3)+6*ys(i+2)-4*ys(i+1)+ys(i))/hs^4; end %Tinh dao ham theo thoi gian o cac node for i=1:Nl dydt(i)=(-EI*dyds4(i)+P(i))/miu; dydt(i+Nl)=y(i); end Bài tốn khảo sát biến dạng uốn tĩnh mía clc; clear; global R dth Ia As ro Es Ns hs smesh ps P EI miu Ep Ev Xc Yc ; global Rd hd xd nq omeg alfa0 R=0.1; dth=0.03; Ia=pi*dth^4/64; As=pi*dth^2/4; ro=1.2e3; miu=As*ro; Es=2.5e8; EI=Es*Ia; 187 %================= thong so lo danh la Rd=0.2; hd=0.025+Rd; xd=0.2+Rd; nq=800; omeg=pi*nq/30; % pt canh dap (xs-xd)^2+(ys-hd)^2=Rd^2 % xs=xd-Rd*cos(alfa); ys=hd-Rd*sin(alfa) alfa0 ys=0 alfa0=asin(hd/Rd); Ns=30; Smax=xd; %chieu dai ngon mia xtl=[]; ytl=[]; Mtl=[]; Ntl=[]; alfa=pi; sgraph=[]; alfaplot=[]; x=[Smax,50]; %dieu kien dau for alfa=0:pi/20:pi/2+pi/10 Xc=xd-Rd*sin(alfa); Yc=hd-Rd*cos(alfa); Smax=xd-0.05*alfa/(pi/2); %chieu dai giam theo thoi gian x0=x; [x, fval]=fsolve(@fmia4,x0) Ln=x(1); %chieu dai mia den diem dat luc Fb=x(2); %luc tac dung s_interval=linspace(Ln,0,Ns); %ma tran theo chieu dai s y0=[0 0,-Fb/EI,0,0,0]; %time interval and initial conditions %y(1)=fi; y(2)=dfi/ds; y(3)=d2fi/ds2; y(4)=N;y(5)=x; y(6)=y %solution [s,y] = ode45(@tinh4_dydt,s_interval, y0); fi0=-y(Ns,1); yfi=y(:,6); xfi=y(:,5); xl=xfi*cos(fi0)-yfi*sin(fi0); yl=xfi*sin(fi0)+yfi*cos(fi0); xl=xl-xl(Ns); yl=yl-yl(Ns); xmax=xl(1)+(Smax-Ln)*cos(fi0); ymax=yl(1)+(Smax-Ln)*sin(fi0); if Smax>Ln xl=[xmax;xl]; yl=[ymax;yl]; xtl=[xtl,xl]; ytl=[ytl,yl]; Mtl=[Mtl,y(:,2)*EI]; Ntl=[Ntl,y(:,4)]; 188 alfaplot=[alfaplot,alfa*180/pi]; figure(1) plot(xl,yl,'LineWidth',2); hold on grid on sgraph=[sgraph,(linspace(Ln,0,Ns))']; end end circle=rsmak('circle',Rd,[xd,hd]); fnplt(circle), axis equal hold off xlabel('X,m'); ylabel('Y,m'); sigma=(Mtl/Ia*dth/2-4*Ntl/pi/dth^2)/1e6; II=size(Mtl); str =num2str(alfaplot',3); figure(2); for i=1:II(2) plot(sgraph(:,i)*1000,Mtl(:,i),'LineWidth',2,'color',rand(1,3)); hold on end xlabel('s,mm'); ylabel('Mu,Nm'); grid on legend(str); figure(3); for i=1:II(2) plot(sgraph(:,i)*1000,sigma(:,i),'LineWidth',2,'color',rand(1,3)) ; hold on end xlabel('s,mm'); ylabel('\sigma,Mpa'); grid on legend(str); sigmamax=max(sigma); figure(4); plot(alfaplot,sigmamax,'LineWidth',2); grid on xlabel('\alpha^{O}'); ylabel('\sigma_{max},Mpa'); a=max(sigmamax) 189 Bài toán khảo sát quĩ đạo bay thân mía clc; clear; global Lm g v0 kf R fm fm fl kf vq beta hl rol rok kd t1 tend to at; g=0;%9.81; kf=1; v0=4.5; R=0.1; fm=0.4; Lm=2.0; rol=30; rok=1.1; hl=10; kd=1; vq=0; beta=pi*30/180; % goc cuar gio quat fl=0.6; fi0=-25*pi/180; %X=y1; Y=y2,fi=y3 %dx/dt=y4; dY/dt=y5; dfi/dt=y6 %smesh=linspace(0,Ln,Nl); %dung cho ve thi theo s %dien kien t=0 for v0=5:0.5:8 for vq=0:0.5:0 xm0=R*sin(fi0); ym0=-R*cos(fi0); dxm0=v0*cos(fi0); dym0=v0*sin(fi0); dfi0=0; to=0; dto=0; tend=0.8; Nt=10; dt=tend/(Nt-1); t_interval=linspace(0,tend,Nt); %time interval and initial conditions y0=[xm0 ym0 fi0 dxm0 dym0 dfi0 to dto]; opts = odeset('RelTol',1e-4,'AbsTol',1e-6); %solution [t,y] = ode45(@miabay_dydt,t_interval, y0,opts); yend=y(:,2)+Lm*sin(y(:,3)/2); yend1=y(:,2)-Lm*sin(y(:,3)/2); xend=y(:,1)+Lm*cos(y(:,3))/2; xend1=y(:,1)-Lm*cos(y(:,3))/2; 190 dxm=y(4); dym=y(5); xla=y(:,1)+y(:,7).*cos(y(:,3)); yla=y(:,2)+y(:,7).*sin(y(:,3)); fi=y(:,3); dxm=y(:,4); dym=y(:,5); dfi=y(:,6); to=y(:,7); dto=y(:,8); dxla=dxm+dto.*cos(fi)-to.*dfi.*sin(fi); dyla=dym+dto.*sin(fi)+to.*dfi.*cos(fi); %fi*180/pi %b1=dto.*cos(fi); %b2=dto.*sin(fi); %gocfi1=atan(b2./b1)*180/pi Nt1=Nt; t1=tend; for i=Nt:-1:1 if y(i,7)v0 kf=1; else kf=-1; end lx=((xm-R*sin(fi))^2+(ym+R*cos(fi))^2)^0.5; %khoang cach tu tam toi diem ty % Tinh phan luc Nm if lx