Luận án tiến sĩ kỹ thuật cơ khí động lực nghiên cứu thiết bị biến đổi năng lượng sóng điện ứng dụng cho khai thác năng lượng sóng tại vùng biển việt nam

LỜI CẢM ƠN Trong thời gian học tập nghiên cứu hoàn thành luận án “Nghiên cứu thiết bị biến đổi lượng sóng - điện ứng dụng cho khai thác lượng sóng vùng biển Việt Nam” trường Đại học Bách khoa Hà nội tác giả nhận nhiều giúp đỡ từ tổ chức cá nhân: Tác giả xin bày tỏ lòng biết ơn sâu sắc đến tập thể nhà khoa học là: GS TS Nguyễn Thế Mịch PGS.TS Đặng Thế Ba tận tình hướng dẫn bảo vạch định hướng khoa học để tác giả hoàn thành luận án Tác giả xin gửi lời cảm ơn chân thành đến thầy, cô Bộ môn Máy tự động thủy khí, Viện khí động lực, trường đại học Bách khoa Hà Nội Tập thể giảng viên Bộ môn Cơ kỹ thuật – Đại học Công nghệ - ĐHQGHN; Hội học Thủy khí; Hội máy thủy khí tạp chí khoa học động viên giúp đỡ hỗ trợ tinh thần vật chất thời gian tác giả thực luận án Tác giả xin cảm ơn đến anh em đồng nghiệp nói chung anh em đồng nghiệp Viện khoa học Thủy lợi Miền Trung Tây Nguyên động viên, giúp đỡ, tạo điều kiện thuận lợi cho tác giả q trình nghiên cứu cơng tác Cuối tác giả xin bày tỏ lòng biết ơn sâu sắc đến người thân gia đình động viên, tạo điều kiện khuyến khích cho tác giả vượt qua khó khăn để hồn thành luận án tiến sĩ Nghiên cứu sinh Phùng Văn Ngọc i LỜI CAM ĐOAN Tên đề tài luận án: “Nghiên cứu thiết bị biến đổi lượng sóng - điện ứng dụng cho khai thác lượng sóng vùng biển Việt Nam” Tôi xin cam đoan đề tài luận án tơi hồn tồn tơi thực hướng dẫn tận tình tập thể nhà khoa học gồm GS.TS Nguyễn Thế Mịch, PGS.TS Đặng Thế Ba Những kết nghiên cứu luận án số liệu đo đạc thí nghiệm Phịng thí nghiệm trường trường Đại học Công Nghệ ĐHQGHN hồn tồn trung thực chưa có cơng bố Hà nội, ngày 22 tháng 08 năm 2019 Thay mặt tập thể hƣớng dẫn khoa học Ngƣời cam đoan GS TS Nguyễn Thế Mịch Phùng Văn Ngọc ii MỤC LỤC LỜI CẢM ƠN i LỜI CAM ĐOAN ii MỤC LỤC iii DANH MỤC CÁC KÝ HIỆU vi DANH MỤC CHỮ CÁI LATINH viii DANH MỤC CÁC CHỮ VIẾT TẮT ix DANH MỤC CÁC BẢNG BIỂU x DANH MỤC CÁC HÌNH VÀ ĐỒ THỊ xi MỞ ĐẦU Chƣơng TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Nguồn lƣợng từ biển 1.2 Năng lƣợng sóng biển 1.2.1 Các đặc trưng sóng biển 1.2.2 Năng lượng vùng hình thành sóng 1.2.3 Đặc trưng lượng sóng biển .11 1.3 Tiềm năng lƣợng sóng biển Việt Nam .15 1.4 Tổng quan tình hình nghiên cứu phát triển thiết bị chuyển đổi lƣợng sóng biển 21 1.4.1 Sơ lược số nguyên lý thiết bị chuyển đổi lượng sóng 21 1.4.2 Tổng quan thiết bị chuyển đổi lượng sóng Thế giới 25 1.4.3 Các thiết bị chuyển đổi lượng sóng xa bờ .28 1.4.4 Các thiết bị chuyển đổi lượng sóng gần bờ 34 1.5 Tổng quan tình nghiên cứu thiết bị chuyển đổi lƣợng nƣớc 37 1.6 Kết luận chƣơng 39 Chƣơng CƠ SỞ NGHIÊN CỨU LỰA CHỌN THIẾT BỊ CHUYỂN ĐỔI NĂNG LƢỢNG SÓNG BIỂN 41 2.1 Cơ sở tính tốn dịng lƣợng sóng biển 41 2.2 Các đặc trƣng vùng nghiên cứu phát triển thiết bị 41 2.3 Cơ sở tính tốn tƣơng tác sóng biển với cơng trình 43 2.2.1 Các giả thuyết sóng cơng trình .44 iii 2.3.2 Xác định tải trọng sóng tác dụng lên cơng trình 45 2.4 Nghiên cứu lựa chọn thiết bị chuyển đổi lƣợng sóng biển .46 2.4.1 Lựa chọn nguyên lý chuyển đổi lượng sóng phù hợp mục tiêu để nghiên cứu phát triển 47 2.4.2 Lựa chọn dạng phao tiếp nhận cho thiết bị .49 2.5 Phƣơng pháp tính tốn mơ hình phao dạng hộp chữ nhật 51 2.6 Nghiên cứu nguyên lý làm việc thiết bị sử dụng dạng phao trụ chóp nón .54 2.6.1 Nguyên lý chuyển đổi lượng sóng phao dạng chóp nón 54 2.6.2 Mơ hình cấu máy phát 56 2.6.3 Nguyên lý hoạt động mơ hình 57 2.7 Kết luận chƣơng 57 Chƣơng NGHIÊN CỨU ĐỘNG LỰC HỌC THIẾT BỊ CHUYỂN ĐỔI NĂNG LƢỢNG SÓNG BIỂN 58 3.1 Xây dựng mơ hình tính tốn thiết bị 58 3.1.1 Các phương trình học thiết bị 58 3.1.2 Phương trình chuyển động phao 1: 59 3.1.3 Phương trình chuyển động phao 60 3.2 Phƣơng trình cho hệ thống máy phát 63 3.2.1 Mô hình điện từ trường máy phát 63 3.2.2 Áp dụng cấu hình máy phát nam châm chuyển động thẳng 64 3.3 Phƣơng trình lực điện từ máy phát 68 3.4 Tính tốn mơ hoạt động hệ thống chuyển đổi lƣợng sóng sang điện 69 3.4.1 Giới thiệu phầm mềm mô 69 3.4.2 Tính phân bố từ trường nam châm vĩnh cửu FlexPDE 70 3.4.3 Tính tốn mơ hoạt động máy phát Matlab 77 3.5 Kết luận chƣơng 81 Chƣơng THIẾT KẾ CHẾ TẠO MƠ HÌNH THỬ NGHIỆM THIẾT BỊ CHUYỂN ĐỔI NĂNG LƢỢNG SÓNG BIỂN 83 4.1 Tính tốn khảo sát đặc tính, tính tốn thiết kế cho thiết bị có P= 50W 83 4.2 Thiết kế, chế tạo hệ thống chuyển đổi lƣợng sóng thử nghiệm 87 4.2.1 Thiết kế, chế tạo máy phát điện 87 4.2.2 Thiết kế, chế tạo hệ thống phao thu nhận lượng 89 iv 4.2.3 Thử nghiệm, hoàn thiện thiết bị 91 4.3 Đề xuất mơ hình chuyển đổi lƣợng sóng 5kW để phát triển thiết bị quy mơ thực tế 101 4.4 Khảo sát số đặc tính 104 4.5 Kết luận chƣơng .110 a) Kết luận chung 111 b) Khuyến nghị 112 TÀI LIỆU THAM KHẢO .114 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA LUẬN ÁN 118 v DANH MỤC CÁC KÝ HIỆU Awp: Phần diện tích tiếp xúc với nước biển với Awp = BL ( B chiều rộng L chiều dài phao) B: Bề rộng phao (m) BT: Cảm ứng từ (T) Dc: Điện cảm D : Đường kính phao thứ cấp (m) D1: Đường kính phao sơ cấp (m) E : Năng lượng sóng tổng ( lb.ft N.m ) Ekz : Động toàn phần dao động phao (KJ) Epz : Thế toàn phần dao động phao (KJ) ET: Cường độ điện trường (V/m) HT: Cường độ từ trường (A/m) Hs : Chiều cao sóng (m) g: gia tốc trọng trường g = 9,81 (m/s2) m: khối lượng kèm nước chuyển động theo m = ρdAwp Z0: Biên độ dao động nhấp nhô cực đại phao (m) F0: Biên độ lực dao động (N) ω: Tần số góc (rad/s) J: Mật độ dòng điện (A/m2) ρc : mật độ điện tích (C/m3) M : độ từ hóa (A/m) n : Số vòng dây T : Chu kỳ (s) f : Tần số (Hz) c : Vận tốc lan truyền (m/s) ρ : Khối lượng riêng (kg/m3) Plt : Công suất lý thuyết thiết bị (kW) Ptt : Công suất thực tế thiết bị (kW) V : Vận tốc truyền sóng (m/s) P : Cơng suất sóng ( ft.lb/sec W, kW) f : Tần số ( 1/T ) ( Hz ) c : Vận tốc sóng ( ft/s m/s) h : Độ sâu mực nước biển ( ft m ) vi L: Chiều dài phao dn: Phần ngập nước phao (m) (m) Fz: Lực kích thích sóng lên phao (KN) F0: Là biên độ lực dao động (KN) Pz: Công suất thu từ phao (KW) Tz: Chu kỳ dao động nhấp nhô tự nhiên (s) mw: Khối lượng phần nước biển tác động vào phao m: khối lượng kèm nước chuyển động theo m = ρdAwp Sb(t) dịch chuyển đứng phao khỏi vị mặt nước tĩnh ban đầu mb1 khối lượng phao sb (t ) gia tốc chuyển dịch đứng phao Fe,b(t) lực kích thích theo chiều đứng sóng tới Fr,b(t) lực theo chiều đứng sóng phát xạ Fb,b(t) tổng lực đẩy thủy tĩnh lực trọng trường Fb,v(t) lực cản Fb,f(t) lực ma sát nhớt nước Fb,u(t) lực điện từ trường máy phát Z: Biên độ dao động phao hình hộp chữ nhật (m) vii DANH MỤC CHỮ CÁI LATINH λ: bước sóng ( ft m) γ : Trọng lượng riêng (N/m3) η : Hiệu suất thiết bị (%) εo : Độ điện thẩm chân không (εo = 8,854.10-12 F/m) ε : Độ điện thẩm môi trường μo : Độ từ thẩm chân không (μo = 4π 10-7 T.m/A) μ : Độ từ thẩm mơi trường ωz: Tần số góc dao động (rad/s) viii DANH MỤC CÁC CHỮ VIẾT TẮT Ký hiệu Tiếng anh Tiếng việt CERT Commission Energy Research and Technology Ủy ban công nghệ lượng EVN Vietnam Electricity Tập đoàn điện lực viện nam HEB Hydro Electric Barrel Hydro Electric Barrel IEA International Energy Agency Ủy ban lượng giới LPM Linear power management Thiết bị chuyển đổi lượng sóng biển chuyển động tịnh tiến NL Energy Năng lượng NLTT Renewable energy Năng lượng tái tạo NLSB Wave ocean energy Năng lượng sóng biển NTB Nam Trung Bộ OWAP Ocean Wave Air Piston Chuyển đổi lượng sóng biển piston nén khí OWC Oscillating Water Column Thiết bị chuyển đổi lượng sóng nhờ động cột nước OES Ocean energy systems Hệ thống lượng biển R&D Research & Development Research and Development VKHTLVN Viện khoa học thủy lợi Việt nam VN Việt nam ix DANH MỤC CÁC BẢNG BIỂU Bảng 1.1 Mức lượng sóng biển số vùng biển giới Bảng 1.2 Năng lượng TB tháng số vùng có tiềm (kW/m) 19 Bảng 1.3 Các thiết bị chuyển đổi lượng sóng xa bờ thực 33 Bảng 1.4 Các thiết bị gần bờ phát triển giới 36 Bảng Các thơng số kích thước khối lượng phụ thuộc vào dạng phao 50 Bảng Các thông số thiết kế ban đầu thiết bị 84 Bảng Thống kê chi tiết phận cấu thành phao1 phao 90 Bảng 4.3 Năng lượng thu theo kích thước phao (cả ma sát) 100 Bảng 4.4 Thông số đặc trưng chủa máy phát điện LPM 103 Bảng 4.5 Cơng suất trung bình tối đa ống dây ứng với điều kiện sóng vùng Nam Trung Bộ 108 x Tháng I II III IV V VI VII VIII IX X XI XII Cả năm TRAM SO: 66 DO SAU: 10,0m 17,7 14,9 8,4 7,0 3,1 7,1 15,5 6,8 6,1 3,4 6,1 24,4 10,0 48,6 45,0 21,6 8,5 1,7 5,2 10,6 4,5 2,9 4,0 12,1 38,5 16,9 56,5 53,9 22,7 8,6 1,7 5,3 10,6 4,6 2,9 4,2 14,7 49,2 19,6 37,4 34,8 19,4 7,0 2,1 5,3 16,3 17,6 13,9 7,0 3,5 9,5 12,2 7,5 5,7 3,3 5,7 19,6 10,2 28,9 27,1 19,4 6,8 2,4 5,7 TOA DO: 9o15E 106o00N TRAM SO: 67 DO SAU: 8,0m TOA DO: 9o00E 105o30N TRAM SO: 68 DO SAU: 15,0m TOA DO: 9o00E 105o45N TRAM SO: 69 DO SAU: 15,0m 7,9 5,5 2,8 3,7 10,1 26,9 13,6 TOA DO: 8o45E 105o30N TRAM SO: 70 DO SAU: 15,0m TOA DO: 8o30E 105o15N TRAM SO: 71 DO SAU: 15,0m TOA DO: 8o30E 105o00N PL - 12 9,2 5,7 4,0 3,7 8,3 25,4 12,2 Tháng I II III IV V VI VII VIII IX X XI XII Cả năm TRAM SO: 72 DO SAU: 26,0m 0,0 0,2 0,4 1,2 6,4 18,0 17,4 22,5 15,2 3,9 0,2 0,1 7,1 0,1 0,0 0,1 0,2 0,2 0,6 0,6 0,7 0,5 0,2 0,2 0,2 0,3 0,0 0,0 0,1 0,2 0,9 1,8 1,9 2,4 1,6 0,4 0,0 0,0 0,8 0,0 0,1 0,2 0,6 3,5 9,2 9,1 11,3 7,7 2,0 0,1 0,0 3,6 0,3 0,1 0,2 0,2 0,3 0,7 0,7 0,9 0,6 0,3 0,6 0,8 0,5 0,4 0,2 0,5 1,0 5,2 12,2 12,8 16,8 11,0 2,7 1,0 1,3 5,4 TOA DO: 10o15E 103o45N TRAM SO: 73 DO SAU: 5,0m TOA DO: 10o15E 104o15N TRAM SO: 74 DO SAU: 5,0m TOA DO: 10o15E 104o30N TRAM SO: 75 DO SAU: 13,0m TOA DO: 10o00E 104o15N TRAM SO: 76 DO SAU: 13,0m TOA DO: 10o00E 104o30N TRAM SO: 77 DO SAU: 14,0m TOA DO: 9o45E 104o30N PL - 13 Tháng I II III IV V VI VII VIII IX X XI XII Cả năm TRAM SO: 78 DO SAU: 15,0m 0,7 0,3 0,5 1,0 5,5 13,4 13,5 17,8 11,5 3,0 1,7 2,0 5,9 0,4 0,2 0,5 1,1 5,7 14,9 14,8 18,5 12,7 3,4 1,0 1,2 6,2 0,3 0,2 0,5 1,2 6,4 18,0 17,5 22,5 15,4 4,0 0,8 0,8 7,3 0,4 0,2 0,5 1,2 6,4 17,6 17,4 22,4 15,1 4,0 1,0 1,1 7,3 0,4 0,2 0,5 1,2 6,5 18,0 17,6 22,6 15,4 4,1 1,0 1,0 7,4 0,0 0,2 0,4 1,2 5,7 14,9 14,9 18,6 12,7 3,3 0,2 0,1 6,0 TOA DO: 9o30E 104o30N TRAM SO: 79 DO SAU: 20,0m TOA DO: 9o15E 104o30N TRAM SO: 80 DO SAU: 30,0m TOA DO: 9o00E 104o30N TRAM SO: 81 DO SAU: 18,0m TOA DO: 8o45E 104o30N TRAM SO: 82 DO SAU: 30,0m TOA DO: 8o30E 104o30N TRAM SO: 83 DO SAU: 10,0m TOA DO: 8o30E 104o45N PL - 14 Phụ lục 2: Phụ lục công suất phát điện khả phát điện theo tháng Các phương án sử dụng phao khác Đơn vị 10 11 12 tính Bình qn Phao cao 13m Công suất phát điện MW Khả Triệu phát điện KWh 438,6 295,2 220,2 94,7 155,0 315,3 318,7 381,6 272,8 162,1 157,6 339,6 262,6 326,3 198,4 163,8 68,2 115,3 227,0 237,1 283,9 196,4 120,6 113,5 252,6 191,9 2, Phao cao 3,5m Mức giảm công suất % 15,0 4,4 4,5 0,5 0,2 3,7 5,9 PL - 15 7,6 2,7 0,8 2,4 7,6 4,6 Công suất phát điện MW Khả Triệu phát điện KWh 372,4 282,1 210,2 94,2 154,7 303,8 299,8 352,6 265,5 160,8 153,8 313,9 247,0 277,3 189,6 156,4 67,8 115,1 218,7 223,0 262,3 191,1 119,7 110,7 233,6 180,4 Phao cao 3,5m Mức giảm công suất Công suất phát điện % MW Khả Triệu phát điện KWh 21,1 8,4 8,0 1,8 1,4 8,4 10,5 13,8 6,6 1,8 4,0 12,7 8,2 316,7 269,6 200,7 93,7 154,4 292,7 282,0 325,7 258,3 159,6 150,1 290,2 232,8 235,6 181,2 149,3 67,5 114,9 210,7 209,8 242,3 186,0 118,7 108,0 215,9 170,0 PL - 16 Phụ lục 2: Sơ đồ mật độ lượng vùng biển khu vực tỉnh Khánh Hòa PL - 17 Phụ lục 3: Sơ đồ mật độ lượng vùng biển khu vực tỉnh Bình Thuận PL - 18 Phụ lục 4: Sơ đồ mật độ lượng vùng biển khu vực tỉnh Ninh Thuận PL - 19 Phụ lục 5: Chương trình tính tốn thơng số thiết bị phần mềm Mathlab % Giai he phuong trinh vi phan mo ta hoat dong cua thiet bi phao % PHUNG VAN NGOC, DANG THE BA AND TA QUANG HUNG - UET %% clear all clc global omega global mb1 Feb Rb Rb2 Sb global md1 Fed Rd Rd2 Rdb Sd global Rngoai Rtrong L nturn lpm dpm phih ncoil nslice lamda %% two_bouy_lpm_dat ; %% cac dieu kien ban dau zb0=0; ub0=0; zd0=0; ud0=0; i0=0; n=0; %% Thoi gian mo phong tmax=scks*tm; dt=0.01; %% Dac trung cua song bien omega=2*pi/tm; k=omega*omega/g ; buocsong=2*pi/k; % tinh tan so goc omega % So song % buoc song %% phuong trinh dao dong song bien % figure(1) PL - 20 % hold on % t=[0:dt:tmax]; % y=ksia*cos(omega.*t+pi); %% van toc dich chuyen cua song % v=-omega*ksia*sin(omega*t+pi); % plot(t,y,'-r') % plot(t,v,'g') % title('ptrinh dao dong va van toc dich chuyen song bien' ) % grid on %% Lap tim dai long toi uu ltubeii=[]; Ptbii=[]; Ayb=[]; Febii=[]; Rbii=[]; mb1ii=[]; Ptbkk=[]; ii=0; ltube0=ltube; mtube0=mtube; lb0=lb; mb0=mb; kk=0; for Rngoai=50:2:50 kk=kk+1; for ltubei=9:0.5:12.5 ii=ii+1 ltube=ltubei; ltubeii(ii)=ltube; PL - 21 %% Tinh he so cho phuong trinh phao % He so luc day Archimedes Swb = pi*rb*rb ; %tiet dien phao Sb=row*g*Swb; %% So hang luc quan tinh phao % Khoi luong kem cua phao mb33=2/3*pi*rb*rb*rb*muyb33*row; % Khoi luong cua long mtube= pi*ltube*dtube*(2*rtube+dtube)*rotube; dmtube=mtube-mtube0; % Khoi luong tang them c long % Khoi luong nuoc kem theo long mtube33=pi*ltube*dtube*(2*rtube+dtube)*row ; % khoi luong kem long % mtube33 = 0; % Tong khoi luong quan tinh cua phao va long mb1= mb+mtube+mb33+mtube33; %% Cac he so luc can len phao % Luc can kem dlb=dmtube/(pi*rb*rb*row); %Do chim them tang khoi luong long Rb33 = exp(-2*k*(lb+dlb))*omega*row * 2/3 *pi*rb*rb*rb*eb33; % - luc can nuoc len thiet dien phao Rbv = 0.5*row*cvb*pi*rb*rb ; % - Luc can nuoc len thiet dien long (Tube) Rtv = 0.5*row*cvb*pi*dtube*(2*rtube+dtube); % Rbt = ; rrfm = cfm; %he_so_luc_ma_sat_co_khi ??? Rbs = 1200; % he so ma sat khac cua nuoc len phao % ???? PL - 22 Rb = Rb33+Rbs; %tong luc can tuyen tinh len phao Rb2=Rbv+Rtv; % tong can phi tuyen len phao %% - Luc tac dong cua song len phao Feb=exp(-k*(lb+dlb))*kb33*row*g*Swb*ksia; %% PHAO 2-d % So hang dan hoi - luc Archimet Swd = pi*rphao2*rphao2 ; %tiet_dien_cot_ong Sd=row*g*Swd; %% -So hang quan tinh-tong KL phao 2= Khoi luong tinh + kem cua ong noi + kem cua may phat muyd33 = muyb33; md33=pi*rphao2*rphao2*lphao2*muyd33*row; %khoi kem cua doan ong noi may phat md33t=pi*rtube*rtube*ltube*row; %khoi luong kem cua nuoc tube md1=md+md33+md33t; %tong khoi luong c/d cua dia %% So hang can - tong he so luc can len phao ed33 = eb33 ; Rd33 = exp(-2*k*lphao2)*omega*row*pi*rphao2*rphao2*lphao2*ed33; %he so luc can kem cvd = cvb ; Rdv = 0.5*row*cvd*pi*rgene*rgene ; %can nhot len vo may phat %he so can nhot len long phao Rds =10; % He so ma sat khac len phao Rd = Rd33+Rds; Rd2=Rdv; %tong luc can tuyen tinh len phao %Tong can phi tuyen len phao 2; %% - Luc song tac dung len phao kd33 = kb33 ; Fed=exp(-k*lphao2)*kd33*row*g*Swd*ksia; PL - 23 % Fed=0; % Kiem tra o co luc ngoai %% tong luc can sinh co cau lien ket phao voi dia qua may phat Rfbd = 100; % Ma sat long va dia Rfg = 100; % Ma sat cua piston may phat Rabs=0; % he so hap thu chuyen doi dien Rdb=Rfbd+Rfg+Rabs; %% - Ve va kiem tra luc song % lucsong ; %% MAY PHAT DIEN l2pm = 2*dpm +2*dw ; % di chu ky tu truong gom nam cham l1cuon=(l2pm-3*wcachpha)/3 ; % dai cuon day - quan c?a pha nturn = floor(l1cuon/(2*rd))*floor(hcoil/(2*rd)); %So vong day cuon lamda=2*dw+2*dpm; % buoc song tu truong % tiet dien day quan Sd=pi*rd^2; % chieu dai day quan cua mot cuon day cua pha lday=nturn*2*pi*drotor*drotor ; %_tinh gia tri dien tro cua cuon day Rtrong=dientrosuat*lday/Sd; %_tinh he so tu cam L ccua cuon ?ay scuon=pi*drotor*drotor; L=nuy0*(nturn^2)*scuon*l1cuon; %% mac noi tiep ncoil vong day Rtrong=ncoil*Rtrong; L=ncoil*L; %% toi uu cong suat % Rngoai=Rtrong; %% GIAI HE PHUONG TRINH VI PHAN % Kiem tra chuyen dong cua phao PL - 24 Febii(ii)=Feb; Rbii(ii)=Rb; mb1ii(ii)=mb1; %[t,z]=ode23('twobouy_ode1',[0,tmax],[0,0]); %% Phuong trinh day du voi he PT [t,z]=ode23('twobouy_odebdgFU',[0,tmax],[0,0,0,0,0]); zz=z(:,1); Ayb(ii)=max(abs(zz)); IIa=z(:,5)/Rngoai; Pout=IIa.*z(:,5); Ptb=mean(Pout) Ptbii(ii)=Ptb; end Ptbkk(kk)=Ptb; end %% Ve va ghi du lieu file figure(1); plot(ltubeii,Ayb); title('Ayb'); % plot(ltubeii,Ptbii); figure(2); plot(Febii); title('Feb'); %plot(t,z(:,1),t,z(:,3)); %title('Dich chuyen '); %legend('phao 1', 'phao 2'); %grid; figure(3); plot(Rbii); PL - 25 title('Rb'); %plot(t,z(:,2),t,z(:,4)); %title('Van toc'); %legend('phao 1', 'phao 2'); %grid; figure(4); plot(mb1ii); title('mb1'); %plot(t,z(:,5)); %title('Dien the'); %legend('tren tai'); %grid; %% Kiem tra luc dien tu %lucdientu ; %% Tính cơng suat %congsuat ; PL - 26