Nghiên cứu thuật toán ước lượng soc cho modul pin lithium

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP ĐỖ THỊ THÙY DƯƠNG NGHIÊN CỨU THUẬT TOÁN ƯỚC LƯỢNG SoC CHO MODUL PIN LITHIUM LUẬN VĂN THẠC SỸ KHOA HỌC KỸ THUẬT ĐIỀU KHIỂN VÀ TỰ ĐỘNG HÓA Thái Nguyên – Năm 2019 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP ĐỖ THỊ THÙY DƯƠNG NGHIÊN CỨU THUẬT TOÁN ƯỚC LƯỢNG SoC CHO MODUL PIN LITHIUM Chuyên ngành: Kỹ thuật điều khiển tự động hóa LUẬN VĂN THẠC SỸ KHOA HỌC KỸ THUẬT ĐIỀU KHIỂN VÀ TỰ ĐỘNG HÓA NGƯỜI HƯỚNG DẪN KHOA HỌC PGS.TS NGUYỄN VĂN CHÍ Thái Nguyên - Năm 2019 i Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc BẢN XÁC NHẬN CHỈNH SỬA LUẬN VĂN THẠC SĨ Họ tên tác giả luận văn: Đỗ Thị Thùy Dương Đề tài luận văn: Nghiên cứu thuật toán ước lượng SoC cho modul pin lithium Chuyên ngành: Kỹ thuật điều khiển tự động hóa Mã số: Tác giả, Cán hướng dẫn khoa học Hội đồng chấm luận văn xác nhận tác giả sửa chữa, bổ sung luận văn theo biên họp Hội đồng ngày 16/6/2019 với nội dung sau: - Sửa lỗi tả trang 14, 33, 35, 42 Chỉnh sửa số thứ tự hình vẽ luận văn - Chương bổ sung tên thành: Xây dựng mơ hình tốn pin Lithium - Sửa lại mã số - Đánh lại thứ tự công thức luận văn - Thống ký hiệu SoC (State of Charge) - Viết lại phần kết luận chung luận văn, đổi thành “ Kết luận đề xuất” Thái Nguyên, ngày… tháng ….năm 2019 Cán hướng dẫn Tác giả luận văn PGS.TS Nguyễn Văn Chí Đỗ Thị Thùy Dương CHỦ TỊCH HỘI ĐỒNG PGS.TS Võ Quang Lạp ii Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn LỜI CAM ĐOAN Tên là: Đỗ Thị Thùy Dương Sinh ngày 31 tháng 10 năm 1986 Học viên lớp cao học khoá 20 chuyên ngành Kỹ thuật điều khiển tự động hóa Trường đại học kỹ thuật Cơng nghiệp Thái Nguyên Hiện công tác : Trường Cao đẳng nghề số – Bộ quốc phịng Tơi xin cam đoan luận văn “Nghiên cứu thuật toán ước lượng SoC cho modul Pin Lithium” thầy giáo PGS.TS Nguyễn Văn Chí hướng dẫn nghiên cứu tơi với tất tài liệu tham khảo có nguồn gốc, xuất xứ rõ ràng Thái Nguyên, ngày…….tháng … năm 2019 Học viên Đỗ Thị Thùy Dương iii Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn LỜI CẢM ƠN Sau thời gian nghiên cứu, làm việc khẩn trương hướng dẫn tận tình giúp đỡ thầy giáo PGS.TS Nguyễn Văn Chí, luận văn với đề tài “Nghiên cứu thuật tốn ước lượng SoC cho modul Pin Lithium” hoàn thành Tác giả xin bày tỏ lòng biết ơn sâu sắc tới: Thầy giáo hướng dẫn PGS.TS Nguyễn Văn Chí tận tình dẫn, giúp đỡ tác giả hồn thành luận văn Các thầy cô giáo Trường Đại học kỹ thuật công nghiệp Thái Nguyên, số đồng nghiệp, quan tâm động viên, giúp đỡ tác giả suốt q trình học tập để hồn thành luận văn Mặc dù cố gắng hết sức, nhiên điều kiện thời gian kinh nghiệm thực tế thân cịn ít, đề tài khơng thể tránh khỏi thiếu sót Vì vậy, tác giả mong nhận đóng góp ý kiến thầy giáo, cô giáo bạn bè đồng nghiệp cho luận văn tơi hồn thiện Tơi xin chân thành cảm ơn! Thái Nguyên, ngày……tháng……năm 2019 Tác giả luận văn Đỗ Thị Thùy Dương iv Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn MỤC LỤC TRANG PHỤ BÌA…………………………………………………………… ………… ……… i BẢN XÁC NHẬN CHỈNH SỬA LUẬN VĂN THẠC SĨ ii LỜI CAM ĐOAN iii LỜI CẢM ƠN iv MỤC LỤC v DANH MỤC HÌNH ẢNH vi DANH MỤC CÁC CHỮ VIẾT TẮT viii LỜI NÓI ĐẦU 1 Tính cấp thiết đề tài Mục tiêu nghiên cứu Nội dung luận văn CHƯƠNG 1: GIỚI THIỆU VỀ PIN LITHIUM VÀ THAM SỐ SOC 1.1 Giới thiệu pin lithium – Ion 1.2 Ưu điểm Pin Lithium ứng dụng .6 1.3 Vai trò tham số SoC số phương pháp xác định SoC 1.4 Kết luận chương 10 CHƯƠNG 2: XÂY DỰNG MƠ HÌNH TỐN CHO MODUL PIN LITHIUM11 2.1 Mơ hình tốn Pin Lithium 11 2.2 Rời rạc hóa mơ hình pin Lithium Ion 16 2.3 Mơ hình ESC pin Lithium Ion 17 2.4 Xác định tham số mơ hình ESC 19 2.5 Kết xác định tham số mơ hình ESC cho loại Pin 31 2.6 Kết luận chương 39 CHƯƠNG 3: QUAN SÁT SOC CỦA PIN LITHIUM SỬ DỤNG BỘ LỌC KALMAN MỞ RỘNG 40 3.1 Quan sát SoC sử dụng quan sát Kalman mở rộng .40 3.2 Kết quan sát SoC cho pin Lithum Ion Samsung INR18650-25R 20/35A 2500mAh 18650 44 3.2 Kết luận chương 53 KẾT LUẬN VÀ ĐỀ XUẤT 54 TÀI LIỆU THAM KHẢO 55 v Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn DANH MỤC HÌNH ẢNH Hình 1-1: Một loại Pin Lithium 32650 hãng Panasonic: Hình 1-2: Quá trình sạc xả Pin Lithium- Ion: Hình 1-3: Minh họa SoC đồ thị minh họa thay đổi SoC điện áp hở mạch trình nạp xả Hình 1-4: Cấu trúc hệ thống BMS .9 Hình 2-1: Mơ hình điện áp hở mạch Pin – Lithium Ion 11 Hình 2-2: Quan hệ OCV SoC nhiệt độ 250C cho 03 loại Pin Lithium 12 Hình 2-3: Mơ hình pin Lithium kể đến nội trở 13 Hình 2-4: Hiện tượng điện áp khuếch tán Pin 14 Hình 2-5: Mơ hình Pin có kể đến tượng phân cực tuyến tính 14 Hình 2-6: Hiện tượng điện áp trễ 15 Hình 2-7: Mơ hình ESC pin Lithium - Ion 17 Hình 2-8: Sự thay đổi điện áp hai cực pin theo kịch .20 Hình 2-9: Sự thay đổi điện áp hai cực pin theo kịch .22 Hình 2-10: Hiệu suất Coulomb cho loại cell pin Lithium Ion khác 24 Hình 2-11: Quan hệ OCV SoC tương ứng với trình nạp xả bước bước cho loại pin ứng với nhiệt độ cố định .25 Hình 2-12: Quan hệ OCV SoC tương ứng nhiệt độ 00C (bên trái) nhiệt độ thay đổi (bên phải) 26 Hình 2-13: Điện áp OCV kịch 28 Hình 2-14: Điện áp OCV kịch 28 Hình 2-15: Điện áp OCV kịch 29 Hình 2-16: Pin Lithium Ion SAMSUNG INR18650-25R 20/35A 2500mAh 18650 32 Hình 2-17: Quan hệ SoC OCV Pin SAMSUNG nhiệt độ -250C 34 Hình 2-18: Quan hệ SoC OCV Pin SAMSUNG nhiệt độ -150C 34 Hình 2-19: Quan hệ SoC OCV Pin SAMSUNG nhiệt độ -50C .35 Hình 2-20: Quan hệ SoC OCV Pin SAMSUNG nhiệt độ 50C .35 vi Số hóa Trung tâm Học liệu Cơng nghệ thơng tin – ĐHTN http://lrc.tnu.edu.vn Hình 2-21: Quan hệ SoC OCV Pin SAMSUNG nhiệt độ 150C .35 Hình 2-22: Quan hệ SoC OCV Pin SAMSUNG nhiệt độ 350C .36 Hình 2-23: Quan hệ SoC OCV Pin SAMSUNG nhiệt độ 450C .36 Hình 2-24: Quan hệ tham số mơ hình Pin theo nhiệt độ 37 Hình 2-25: So sánh điện áp đầu mơ hình hai trường hợp test xác định OCV test xác định tham số với điện áp thực nghiệm 450C .38 Hình 2-26: Sai lệch điện áp OCV mơ hình ESC thực nghiệm 450C 39 Hình 3-1: Dữ liệu động học cho kịch 48 Hình 3-2: Dữ liệu động học cho kịch 48 Hình 3-3: Dữ liệu động học cho kịch 49 Hình 3-4: Kết quan sát SoC sai lệch quan sát SoC cho kịch 50 Hình 3-5: Kết quan sát SoC sai lệch quan sát SoC cho kịch .51 Hình 3-6: Kết quan sát SoC sai lệch quan sát SoC cho kịch 52 vii Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn DANH MỤC CÁC CHỮ VIẾT TẮT Từ viết tắt Tên tiếng Anh Tên tiếng Việt LTHI Lithium Ion Battery Pin Ion Battery SoC State of Charge Trạng thái nạp OCV Open Circuit Voltage Điện áp hở mạch ESC Enhanced Self Correct Circuit Mơ hình mạch điện tương Model đương pin có kể đến ảnh hưởng điện áp trễ, nhiệt độ, điện trở trong, phân cực điện áp v.v BMS Battery Management Systems Hệ thống quản lý pin viii Số hóa Trung tâm Học liệu Cơng nghệ thơng tin – ĐHTN http://lrc.tnu.edu.vn LỜI NĨI ĐẦU Tính cấp thiết đề tài Ngày Pin Lithium Ion ứng dụng rộng rãi nhiều thiết bị điện quan trọng thiết bị viễn thông, điện tử, điện thoai v.v Và đặc biệt mở xu hướng sử dụng cho ngành cơng nghiệp tơ điện Pin Lithium Ion có có mật độ lượng lớn hơn, chi phí thấp hơn, tốc độ tự xả lâu vòng đời sử dụng dài Phần lớn nghiên cứu gần BMS (Hệ thống quản lý pin) tập trung vào vào loại Pin Ước lượng trạng thái nạp hay cịn gọi ước lượng trạng thái tích điện (SoC) đóng vai trị cực quan trọng nghiên cứu liên quan đến BMS, cho phép cung cấp thơng tin liên quan đến dung lượng lại Pin Ước lượng SoC xác ngăn chặn tượng nạp, tượng xả mức gây tổn hại đến pin ngăn chặn tượng nhiệt, cháy nổ, cân cell pin modul pin Nhiễu yếu tố không mong muốn ảnh hướng đến trình xác định SoC Pin, thay đổi nhiễu ảnh hưởng nhiều đến q trình Có số phương pháp xác định SoC phương pháp đếm Coulomb, phương pháp điện áp hở mạch Hai phương pháp thực theo ngun tắc vịng hở khơng thể tự hiệu chỉnh sai lệch xác định SoC trở lên lớn Chính năm gần tác giả tập trung sử dụng lọc Kalman để ước lượng SoC lọc Kalman có khả ước lượng trạng thái điều kiện có nhiễu Vậy toán ước lượng trạng thái SoC cho Pin Lithium Ion sử dụng nguyên tắc lọc Kalman mang tính cấp thiết Mục tiêu nghiên cứu - Mục tiêu chung Ước lượng trạng thái SoC Pin Lithium Ion sử dụng lọc Kalman mở rộng dựa mơ hình mạch điện tương đương Pin Lithium Ion có xét đến tượng động học Pin nhiệt độ làm việc Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn x k » f (xˆ k+- 1, u k - 1, wk - ) + df (x k - 1, u k - 1, wk - ) dx k - x k - = xˆ k+- 144444444442 44444444443 (x k- - xˆ k+- Aˆ k - + df (x k - 1, u k - 1, wk - ) dwk - wk - = wk - 144444444442 44444444443 (w k- ) (3.4) - wk - ) Bˆ k - Vậy ta có x%k- » Aˆ k - 1xˆ k+- + Bˆ k - 1wˆ k - (3.5) Ma trận hiệp phương sai sai lệch ước lượng T ù é S -x%,k = E ê(x%k )(x%k ) ú» Aˆ k - 1S x+%,k Aˆ kT- + Bˆ k - 1S w%Bˆ kT- ú ëê û (3.6)  Bước 1c: Ước lượng đầu Đầu hệ ước lượng yˆ = E éêh(x k , uk , vk ) Y k - ù » h xˆ k- , uk , vk ú ë û ( ) (3.7) vk = E éêvk ù ë ú û  Bước 2a: Xác định ma trận khuếch đại trạng thái ( y%k = yk - yˆ = h(x k , uk , vk ) - h xˆ k- , uk , vk ) (3.8) Sử dụng khai triển Taylor cho thành phần thứ phương trình ta có ( ) yk » h xˆ k- , u k , vk + dh (x k , u k , vk ) x k = xˆ k- dx k 144444442 44444443 Cˆ k + dh (x k , u k , vk ) vk = vk dvk 144444442 44444443 (v k (x k - xˆ k- ) (3.9) - vk ) Dˆ k Ma trận hiệp phương sai sai lệch đầu S y%,k » Cˆ k - 1S x+%,kCˆ kT- + Dˆ k - 1S v%Dˆ kT- (3.10) Ma trận hiệp phương sai sai lệch ước lượng trạng thái sai lệch đầu 41 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn S -xy%%,k » S x+%,kCˆ kT (3.11) Ma trận khuếch đại trạng thái xác định - Lk = S -x%,kCˆ kT éêCˆ k S -x%,kCˆ kT + Dˆ k S v%,k Dˆ kT ù ú ë û (3.12)  Bước 2b: Hiệu chỉnh trạng thái quan sát xˆ k+ = xˆ k- + Lk (yk - yˆk ) (3.13)  Bước 2c: Cập nhật ma trận hiệp phương sai sai lệch quan sát S x+%,k = S -x%,k - Lk S y%,k LTk (3.14) Tổng kết lại ta có thuật tốn tính tốn lọc Kalman mở rộng sau: Bảng 3-1: Thuật tốn tính tốn lọc Kalman mở rộng Mơ hình khơng gian trạng thái ìï x = f (x , u , w ) k- k- k- ï k í ïï yk = h (x k , u k , vk ) ỵ Định nghĩa ma trận df (x k - 1, u k - 1, wk - ) Aˆ k = dx k - x Cˆ k = dh (x k , u k , vk ) dx k Khởi tạo: x k = xˆ k- k- = xˆ k+- , Dˆ k = df (x k - 1, u k - 1, wk - ) , Bˆ k = dwk - w k- = wk - dh (x k , u k , vk ) vk = vk dvk é + k=0 đặt xˆ 0+ = E éêx ù , S = E ê x - xˆ 0+ ú % x ,0 ë û êë ( )(x - xˆ 0+ T ù ) úûú Tính toán: Từ k =1,2, Ước lượng trạng thái: xˆ k- = f (xˆ k+- 1, u k - 1, wk - ) Ước lượng ma trận hiệp phương sai sai lệch quan sát: T ù é S -x%,k = E ê(x%k )(x%k ) ú» Aˆ k - 1S x+%,k Aˆ kT- + Bˆ k - 1S w%Bˆ kT- ú ëê û ( Ước lượng đầu ra: yˆ = h xˆ k- , uk , vk ) 42 Số hóa Trung tâm Học liệu Cơng nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn Xác định ma trận khuếch đại quan sát: - Lk = S -x%,kCˆ kT éêCˆ k S -x%,kCˆ kT + Dˆ k S v%,k Dˆ kT ù ú ë û Cập nhật trạng thái quan sát: xˆ k+ = xˆ k- + Lk (yk - yˆk ) Cập nhật ma trận hiệp phương sai sai lệch quan sát: S x+%,k = S -x%,k - Lk S y%,k LTk 3.1.2 Quan sát SoC sử dụng quan sát Kalman mở rộng Từ phương trình mơ hình ESC pin Lirthium Ion xác định chương II ìï é- h(k )D t ù ïï ê ú é ùé ù úêz (k ) ú ê Q ïï úé ê1 ê úê i (k ) ù ïï ê úê ú ú úê úêiR (k )ú+ ê B R C ïï X (k + 1) = ê0 AR C ú sgn( i ( k )) ê ú ïí ê úê ú ê êë ú û ú 0 A ( k ) h ( k ) A ( k ) êë úê ú H H ïï ûë û ê ú êë ú ïï û ïï ïï v(k ) = OCV (z (k ),T (k )) + M 0s (k ) + Mh (k ) - å R j i R (k ) - R 0i (k ) j ïïỵ j (3.15) cần xác định ma trận Aˆ k , Bˆ k , Cˆ k , Dˆ k Trước hết ta xác định thành phần để tính tốn ma trận Aˆ k , Bˆ k Giả thiết dòng điện qua pin đo bị nhiễu ik + wk , đo ik Để đơn giản giả thiết hk = Phương trình SoC viết zk + = zk - Dt (i - wk ) Q k (3.16) ¶ zk + ¶ zk z k = zˆk+ = 1, ¶ zk + ¶ wk wk = w = - Dt , Q (3.17) Phương trình dịng điện qua điện trở pin 43 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn iR ,k + ét ù é1 - t ù 1ú ê1 Lú ê ê1 - t ú i + w = êê0 t + ú R ,k + ê ú( k k) ê ú ê ú Oú êëM êë M ú 144442 44443û 1442 443û AR C t j = e - Dt R jC j (3.18) BR C , đạo hàm riêng theo dòng điện nhiễu dòng điện ¶ iR ,k + ¶ iR ,k iR ,k = iˆk+ ¶ i R ,k + = AR C , ¶ wk wk = w (3.19) = BRC Phương trình trễ điện áp hk + = AH ,k hk + (AH ,k - 1)sgn (ik + wk ) - AH ,k = e ¶ hk + ¶ hk hk = hˆk+ wk = w =e (3.20) (ik + wk )gD t Q , đạo hàm riêng theo trễ điện áp nhiễu (ik + wk )gD t Q = AH ,k , ¶ hk + ¶ wk hk = hˆk+ wk = w =- gD t AH ,k + sgn (ik + wk )hˆk+ Q ( ) (3.21) Tiếp theo ta xác định ma trận Cˆ k , Dˆ k Phương trình điện áp đầu pin y (k ) = OCV (z (k )) + M 0s(k ) + Mh(k ) - å R j iR (k ) - R 0i (k ) (3.22) j j Ta có hệ số ¶ yk ¶ sk = M 0, ¶ yk ¶ zk z k = zˆk- = ¶ yk ¶ hk = M, ¶ yk ¶ iR j ,k = - Rj, ¶ yk ¶ vk = (3.23) ¶ OCV (z k ) ¶ zk z k = zˆk- 3.2 Kết quan sát SoC cho pin Lithum Ion Samsung INR18650-25R 20/35A 2500mAh 18650  Dữ liệu mơ hình pin xác định chương 2, sau: 44 Số hóa Trung tâm Học liệu Cơng nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn Bảng 3-2 Dữ liệu mô hình Pin Tên biến OCV0 (1x201) OCVrel (1x201) SoC (1x201) Dữ liệu 2.8685 3.1006 3.2236 3.5534 3.5783 3.6562 3.6573 3.6683 3.6708 3.6964 3.7005 3.7273 3.7306 3.7535 3.7564 3.7647 3.7654 3.7706 3.7715 3.7791 3.7802 3.7897 3.7910 3.8026 3.8045 3.8181 3.8202 3.8336 3.8355 3.8505 3.8530 3.8720 3.8753 3.9029 3.9070 3.9344 3.9378 3.9615 3.9648 3.9886 3.9923 4.0196 4.0237 4.0540 4.0582 4.0866 4.0908 4.1220 4.1266 4.1595 4.1644 -0.00163 -0.00093 -0.00065 -0.00042 -0.00041 0.00012 0.00013 0.00003 0.00000 -0.00015 -0.00015 -0.00015 -0.00015 -0.00018 -0.00018 -0.00013 -0.00010 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00001 -0.00002 -0.00003 -0.00003 -0.00002 -0.00002 0.00001 0.00001 0.00016 0.00020 0.00027 0.00026 0.00012 0.00012 0.00008 0.00008 0.00008 0.00008 0.00005 0.00005 -0.00001 -0.00001 0.00002 0.00002 0.00001 0.00001 0.00001 0.00001 0.00000 0.00500 0.01000 0.04500 0.05000 0.08500 0.09000 0.12500 0.13000 0.16500 0.17000 0.20500 0.21000 0.24500 0.25000 0.28500 0.29000 0.32500 0.33000 0.36500 0.37000 0.40500 0.41000 3.3034 3.6009 3.6587 3.6735 3.7046 3.7335 3.7587 3.7661 3.7726 3.7816 3.7927 3.8062 3.8224 3.8376 3.8552 3.8787 3.9112 3.9414 3.9683 3.9959 4.0280 4.0623 4.0952 4.1311 4.1693 -0.00058 -0.00040 0.00011 -0.00002 -0.00015 -0.00016 -0.00019 -0.00008 0.00000 0.00000 0.00000 -0.00002 -0.00004 -0.00002 0.00003 0.00024 0.00023 0.00011 0.00008 0.00008 0.00004 -0.00001 0.00002 0.00002 0.00001 0.01500 0.05500 0.09500 0.13500 0.17500 0.21500 0.25500 0.29500 0.33500 0.37500 0.41500 3.3649 3.6211 3.6599 3.6765 3.7085 3.7364 3.7602 3.7664 3.7735 3.7828 3.7943 3.8083 3.8246 3.8396 3.8579 3.8823 3.9153 3.9450 3.9715 3.9998 4.0322 4.0664 4.0995 4.1358 4.1742 -0.00053 -0.00039 0.00011 -0.00005 -0.00015 -0.00017 -0.00019 -0.00005 0.00000 0.00000 -0.00001 -0.00002 -0.00004 -0.00002 0.00003 0.00028 0.00021 0.00009 0.00008 0.00007 0.00003 -0.00001 0.00002 0.00001 0.00000 0.02000 0.06000 0.10000 0.14000 0.18000 0.22000 0.26000 0.30000 0.34000 0.38000 0.42000 3.4149 3.6383 3.6614 3.6797 3.7123 3.7393 3.7615 3.7672 3.7746 3.7841 3.7960 3.8102 3.8264 3.8416 3.8605 3.8861 3.9194 3.9482 3.9748 4.0036 4.0366 4.0705 4.1040 4.1403 4.1797 -0.00049 -0.00035 0.00010 -0.00007 -0.00016 -0.00017 -0.00019 -0.00003 0.00000 0.00000 -0.00001 -0.00002 -0.00004 -0.00001 0.00005 0.00029 0.00019 0.00009 0.00008 0.00007 0.00003 -0.00001 0.00002 0.00002 -0.00001 0.02500 0.06500 0.10500 0.14500 0.18500 0.22500 0.26500 0.30500 0.34500 0.38500 0.42500 45 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN 3.4572 3.6506 3.6628 3.6840 3.7159 3.7425 3.7625 3.7680 3.7755 3.7855 3.7976 3.8122 3.8281 3.8438 3.8632 3.8900 3.9231 3.9516 3.9784 4.0074 4.0411 4.0744 4.1083 4.1452 4.1849 -0.00047 -0.00025 0.00009 -0.00011 -0.00015 -0.00016 -0.00018 -0.00002 0.00001 0.00000 -0.00001 -0.00003 -0.00004 -0.00001 0.00007 0.00031 0.00017 0.00009 0.00008 0.00007 0.00001 0.00000 0.00002 0.00001 0.00000 0.03000 0.07000 0.11000 0.15000 0.19000 0.23000 0.27000 0.31000 0.35000 0.39000 0.43000 3.4938 3.6543 3.6644 3.6881 3.7198 3.7461 3.7633 3.7686 3.7769 3.7868 3.7992 3.8141 3.8299 3.8459 3.8661 3.8942 3.9270 3.9549 3.9816 4.0112 4.0454 4.0783 4.1129 4.1498 4.1907 -0.00046 -0.00003 0.00007 -0.00012 -0.00015 -0.00017 -0.00017 0.00000 0.00000 0.00000 -0.00001 -0.00003 -0.00003 -0.00001 0.00009 0.00031 0.00016 0.00008 0.00008 0.00006 0.00000 0.00001 0.00002 0.00002 -0.00001 0.03500 0.07500 0.11500 0.15500 0.19500 0.23500 0.27500 0.31500 0.35500 0.39500 0.43500 3.5255 3.6551 3.6664 3.6921 3.7237 3.7498 3.7642 3.7698 3.7778 3.7884 3.8008 3.8160 3.8318 3.8481 3.8690 3.8986 3.9306 3.9584 3.9853 4.0154 4.0497 4.0825 4.1174 4.1546 4.2025 -0.00044 0.00010 0.00005 -0.00014 -0.00015 -0.00017 -0.00015 -0.00001 0.00000 -0.00001 -0.00001 -0.00003 -0.00003 0.00001 0.00012 0.00029 0.00014 0.00008 0.00008 0.00006 0.00000 0.00001 0.00002 0.00001 -0.00012 0.04000 0.08000 0.12000 0.16000 0.20000 0.24000 0.28000 0.32000 0.36000 0.40000 0.44000 http://lrc.tnu.edu.vn OCV (1x178) SoC0 (1x178) OCVrel (1x178) 0.44500 0.45000 0.45500 0.46000 0.46500 0.47000 0.47500 0.48000 0.48500 0.49000 0.49500 0.50000 0.50500 0.51000 0.51500 0.52000 0.52500 0.53000 0.53500 0.54000 0.54500 0.55000 0.55500 0.56000 0.56500 0.57000 0.57500 0.58000 0.58500 0.59000 0.59500 0.60000 0.60500 0.61000 0.61500 0.62000 0.62500 0.63000 0.63500 0.64000 0.64500 0.65000 0.65500 0.66000 0.66500 0.67000 0.67500 0.68000 0.68500 0.69000 0.69500 0.70000 0.70500 0.71000 0.71500 0.72000 0.72500 0.73000 0.73500 0.74000 0.74500 0.75000 0.75500 0.76000 0.76500 0.77000 0.77500 0.78000 0.78500 0.79000 0.79500 0.80000 0.80500 0.81000 0.81500 0.82000 0.82500 0.83000 0.83500 0.84000 0.84500 0.85000 0.85500 0.86000 0.86500 0.87000 0.87500 0.88000 0.88500 0.89000 0.89500 0.90000 0.90500 0.91000 0.91500 0.92000 0.92500 0.93000 0.93500 0.94000 0.94500 0.95000 0.95500 0.96000 0.96500 0.97000 0.97500 0.98000 0.98500 0.99000 0.99500 1.00000 2.49 2.50 2.51 2.52 2.53 2.54 2.55 2.56 2.57 2.58 2.59 2.60 2.61 2.62 2.63 2.64 2.65 2.66 2.67 2.68 2.69 2.70 2.71 2.72 2.73 2.74 2.75 2.76 2.77 2.78 2.79 2.80 2.81 2.82 2.83 2.84 2.85 2.86 2.87 2.88 2.89 2.90 2.91 2.92 2.93 2.94 2.95 2.96 2.97 2.98 2.99 3.00 3.01 3.02 3.03 3.04 3.05 3.06 3.07 3.08 3.09 3.10 3.11 3.12 3.13 3.14 3.15 3.16 3.17 3.18 3.19 3.20 3.21 3.22 3.23 3.24 3.25 3.26 3.27 3.28 3.29 3.30 3.31 3.32 3.33 3.34 3.35 3.36 3.37 3.38 3.39 3.40 3.41 3.42 3.43 3.44 3.45 3.46 3.47 3.48 3.49 3.50 3.51 3.52 3.53 3.54 3.55 3.56 3.57 3.58 3.59 3.60 3.61 3.62 3.63 3.64 3.65 3.66 3.67 3.68 3.69 3.70 3.71 3.72 3.73 3.74 3.75 3.76 3.77 3.78 3.79 3.80 3.81 3.82 3.83 3.84 3.85 3.86 3.87 3.88 3.89 3.90 3.91 3.92 3.93 3.94 3.95 3.96 3.97 3.98 3.99 4.00 4.01 4.02 4.03 4.04 4.05 4.06 4.07 4.08 4.09 4.10 4.11 4.12 4.13 4.14 4.15 4.16 4.17 4.18 4.19 4.20 4.21 4.22 4.23 4.24 4.25 4.26 -0.0082 -0.0080 -0.0078 -0.0076 -0.0074 -0.0072 -0.0069 -0.0067 -0.0065 -0.0063 -0.0061 -0.0059 -0.0056 -0.0054 -0.0052 -0.0050 -0.0048 -0.0046 -0.0043 -0.0041 -0.0039 -0.0037 -0.0035 -0.0033 -0.0030 -0.0028 -0.0026 -0.0024 -0.0022 -0.0020 -0.0017 -0.0015 -0.0013 -0.0011 -0.0009 -0.0006 -0.0004 -0.0002 0.0001 0.0003 0.0006 0.0008 0.0011 0.0014 0.0017 0.0020 0.0022 0.0025 0.0028 0.0031 0.0034 0.0037 0.0039 0.0042 0.0045 0.0048 0.0051 0.0054 0.0056 0.0059 0.0062 0.0065 0.0068 0.0071 0.0073 0.0076 0.0079 0.0082 0.0085 0.0087 0.0090 0.0093 0.0097 0.0101 0.0105 0.0111 0.0119 0.0126 0.0133 0.0140 0.0147 0.0154 0.0161 0.0168 0.0175 0.0182 0.0190 0.0197 0.0206 0.0216 0.0227 0.0238 0.0249 0.0260 0.0270 0.0281 0.0293 0.0305 0.0319 0.0333 0.0348 0.0363 0.0377 0.0393 0.0409 0.0427 0.0446 0.0465 0.0484 0.0505 0.0527 0.0551 0.0574 0.0599 0.0630 0.0682 0.0742 0.0992 0.1282 0.1450 0.1576 0.1695 0.1821 0.1953 0.2096 0.2258 0.2407 0.2618 0.3212 0.3692 0.4062 0.4375 0.4645 0.4895 0.5150 0.5409 0.5639 0.5842 0.6025 0.6177 0.6299 0.6413 0.6530 0.6655 0.6789 0.6928 0.7075 0.7226 0.7377 0.7523 0.7668 0.7802 0.7932 0.8055 0.8174 0.8288 0.8403 0.8522 0.8645 0.8770 0.8891 0.9005 0.9118 0.9228 0.9337 0.9445 0.9551 0.9655 0.9757 0.9854 0.9930 0.9987 1.0035 1.0079 1.0123 1.0167 1.0211 1.0255 0.000057 0.000056 0.000056 0.000055 0.000054 0.000054 0.000053 0.000053 0.000052 0.000051 0.000051 0.000050 0.000049 0.000049 0.000048 0.000048 0.000047 0.000046 0.000046 0.000045 0.000044 0.000044 0.000043 0.000043 0.000042 0.000041 0.000041 0.000040 0.000040 0.000039 0.000038 0.000038 46 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn OCVQ (1 x 8) OCVeta (1 x 8) Các nhiệt độ xét T (1 x 8) h(T ) (1 x 8) Q(T ) (1 x 8) g(T ) (1 x 8) 0.000038 0.000038 0.000038 0.000039 0.000039 0.000040 0.000040 0.000041 0.000041 0.000041 0.000041 0.000040 0.000039 0.000038 0.000038 0.000037 0.000036 0.000036 0.000035 0.000034 0.000033 0.000033 0.000032 0.000031 0.000030 0.000030 0.000029 0.000028 0.000027 0.000027 0.000026 0.000025 0.000024 0.000024 0.000023 0.000022 0.000021 0.000021 0.000020 0.000021 0.000027 0.000034 0.000040 0.000044 0.000043 0.000043 0.000042 0.000042 0.000041 0.000040 0.000040 0.000039 0.000039 0.000038 0.000041 0.000047 0.000053 0.000056 0.000055 0.000054 0.000054 0.000053 0.000052 0.000053 0.000058 0.000065 0.000068 0.000067 0.000067 0.000066 0.000067 0.000073 0.000079 0.000081 0.000080 0.000078 0.000083 0.000090 0.000093 0.000092 0.000095 0.000104 0.000101 0.000025 -0.000012 -0.000397 -0.000029 0.000121 0.000162 0.000180 0.000200 0.000203 0.000258 0.000256 0.000280 0.000604 0.000028 -0.000005 0.000010 0.000035 0.000066 0.000080 0.000088 0.000038 -0.000016 -0.000082 -0.000203 -0.000302 0.000338 -0.000341 -0.000318 -0.000272 -0.000217 0.000165 -0.000133 -0.000119 -0.000123 -0.0001130.000114 -0.000097 -0.000083 -0.000063 -0.000043 0.000017 0.000002 0.000009 0.000008 -0.000008 -0.000023 -0.000022 -0.000020 -0.000018 -0.000015-0.000013 0.000014 -0.000010 -0.000006 -0.000001 0.000025 0.000068 0.000110 0.000163 0.000217 0.000270 0.000324 0.000377 2.16059783836010 2.16352610129109 2.17686662338975 2.17058839729510 2.17510662791514 2.17529962872337 2.18858227281753 2.17285900776249 0.987641340252708 0.989120659632985 0.991699131770744 0.987604769535184 0.988428318746643 0.985439396412454 0.991240379867523 0.988478517031051 -25 -15 -5 15 25 35 45 2.1497 2.1378 2.1597 2.1877 2.1944 2.1507 2.1515 2.1524 2.150 2.138 2.160 2.188 2.194 2.151 2.152 2.152 250 17.2157 250 250 250 250 250 250 M (T ) (1 x 8) 0.0043 0.0619 0.0086 0.0060 0.0057 0.0042 0.0036 0.0023 M (T ) (1 x 8) 0.1142 0.1301 0.0380 0.0305 0.0227 0.0244 0.0213 0.0196 R (T ) (1 x 8) 0.5967 0.4166 0.2645 0.1622 0.1149 0.0898 0.0781 0.0718 RC (T ) RC (1 x 8) R (T ) (1 x 8) 1.3781 0.8885 0.8188 0.7706 0.9576 2.0578 2.1560 2.6972 0.4655 0.2373 0.0652 0.0219 0.0137 0.0105 0.0093 0.0087 47 Số hóa Trung tâm Học liệu Cơng nghệ thơng tin – ĐHTN http://lrc.tnu.edu.vn Các liệu động học cho kịch 1,2,3 Kịch 1: Hình 3-6: Dữ liệu động học cho kịch Kịch 2: Hình 3-7: Dữ liệu động học cho kịch 48 Số hóa Trung tâm Học liệu Cơng nghệ thơng tin – ĐHTN http://lrc.tnu.edu.vn Kịch 3: Hình 3-8: Dữ liệu động học cho kịch  Các ma trận hiệp phương sai chọn T xˆ 0+ = éê0 0ù ú, ë û é1e - 0 ê S x+%,0 = êê 1e - ê 0 1e êë S v%,k = 2e - S w% = ù ú ú ú 3ú ú û  Kết quan sát cho kịch sau: 49 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn Kịch Hình 3-9: Kết quan sát SoC sai lệch quan sát SoC cho kịch 50 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn Kịch Hình 3-10: Kết quan sát SoC sai lệch quan sát SoC cho kịch 51 Số hóa Trung tâm Học liệu Cơng nghệ thơng tin – ĐHTN http://lrc.tnu.edu.vn Kịch Hình 3-11: Kết quan sát SoC sai lệch quan sát SoC cho kịch Nhận xét: Đối với kịch 1: Lượng xả nhiều lượng nạp (thể hình 3-1), điện áp giảm dần từ 4.2V xuống 3.5V, tổng lượng nạp vào pin 0.8Ah, tổng lượng xả 2.8Ah Kết quan sát SoC hình 3-4 cho thấy SoC giảm từ 100% xuống 5% Đường quan sát mơ hình ESC gần trùng với đường SoC thực tế xác định liệu thực nghiệm Sai lệch quan sát mơ tả hình 3-4 có giá trị khoảng ±1% Đối với kịch 2: Đây kịch cho nạp đầy sau cho xả với tốc độ C/30 đến điện áp mức trì mức này, nghĩa SoC nhỏ 52 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn biểu thị Hình 3-2 Kết quan sát SoC Hình 3-5 cho thấy SoC giảm nhanh từ 100% xuống 2.5% giai đoạn trì điện áp SoC giảm xuống 0% Đường quan sát mơ hình ESC gần trùng với đường SoC thực tế xác định liệu thực nghiệm Sai lệch quan sát mô tả hình 34 có giá trị khoảng ±0.1% Đối với kịch 3: Đây kịch cho nạp đầy đến điện áp mức max trì mức này, nghĩa SoC gần 100% biểu thị Hình 3-3 Kết quan sát SoC Hình 3-6 cho thấy SoC tăng đến 100% Đường quan sát mơ hình ESC gần trùng với đường SoC thực tế xác định liệu thực nghiệm Sai lệch quan sát mô tả hình 3-4 có giá trị khoảng ±1% 3.2 Kết luận chương Trong chương luận văn, tác giả tìm hiểu lọc Kalman mở rộng áp dụng cho hệ phi tuyến Sau áp dụng ngun tắc cho mơ hình ESC xây dựng chương nhằm quan sát SoC cho pin Lithium Ion dựa 03 kịch dòng áp nhiệt độ Pin Kết cho thấy SoC pin xác định xác với sai số đảm bảo yêu cầu nhỏ 2% 53 Số hóa Trung tâm Học liệu Cơng nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn KẾT LUẬN VÀ ĐỀ XUẤT a Kết luận Luận văn với tên đề tài “Nghiên cứu thuật toán ước lượng trạng thái SoC cho modul pin Lithium” tìm hiểu pin Lithium Ion, nguyên lý hoạt động ứng dụng Pin Lirthium Ion Một tham số quan trọng cần xác định q trình làm việc Pin SoC Một phương pháp thực nghiên cứu năm gần phương pháp dựa quan sát SoC từ biến dòng điện, điện áp hở mạch nhiệt độ làm việc pin mang lại độ xác cao Xây dựng mơ hình cho cell pin Lithium Ion dựa phương pháp mạch điện tương đương, mơ hình ESC Mơ hình ESC kể đến tượng phân cực tuyến tính, trễ điện áp q trình nạp xả Qua xác định quan hệ điện áp hở mạch OCV Các tham số mô hình pin xác định thơng qua kịch liệu thực nghiệm Pin Tìm hiểu lọc Kalmam mở rộng áp dụng cho hệ phi tuyến Sau áp dụng ngun tắc cho mơ hình ESC Kết cho thấy SoC pin xác định xác với sai số đảm bảo yêu cầu nhỏ 2% b Đề xuất hướng nghiên cứu tiếp theo Trong luận văn xác định tham số pin SAM (SAMSUNG), INR18650-25R 20/35A 2500mAh 18650 Lithium ion, pin sạc cho thiết bị điện tử cầm tay Hiện nước phát triển, LIB trọng phát triển quân đội, ứng dụng cho phương tiện di chuyển chạy điện kĩ thuật hàng khơng Nó kì vọng thay cho ắc qui chì tơ, xe máy loại xe điện Hơn nữa, việc thay cho ắc qui chì cịn hứa hẹn việc đảm bảo mơi trường sạch, nâng cao an tồn sử dụng tránh việc sử dụng dung dịch điện ly chứa axit, hạn chế phát thải kim loại nặng môi trường, pin Li-ion đảm bảo điện ngang với ắc qui 54 Số hóa Trung tâm Học liệu Cơng nghệ thơng tin – ĐHTN http://lrc.tnu.edu.vn TÀI LIỆU THAM KHẢO [1] https://vi.wikipedia.org/wiki/Pin_Li-ion [2] https://www.sparkfun.com/products/13813 [3] Ruifeng Zhang 1,2,3,* ID , Bizhong Xia , Baohua Li , Libo Cao , Yongzhi Lai , Weiwei Zheng , Huawen Wang and Wei Wang 3, State of the Art of Lithium-Ion Battery SOC Estimation for Electrical Vehicles, Energies 2018, 11, 1820; doi:10.3390/en11071820 [4] Gregory L.Plett, “Battery Modeling”, Vol 1, Artech House, London, 2015 [5] M Sami Fadali, Antonio Visioli, “Digital Control Engineering, Analysis and Design” Elsevier Press, 2009 [6] Gregory L.Plett, “Applied Kalman Filtering”, Artech House, London, 2017 [7] Nguyễn Vĩnh Thụy, Nguyễn Văn Chí, “Ước lượng SoC Pin Lithium Ion sử dụng quan sát bình phương cực tiểu hồi quy”, Tạp chí Khoa học Cơng nghệ, Đại học Thái Nguyên, 2018 55 Số hóa Trung tâm Học liệu Công nghệ thông tin – ĐHTN http://lrc.tnu.edu.vn ... ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP ĐỖ THỊ THÙY DƯƠNG NGHIÊN CỨU THUẬT TOÁN ƯỚC LƯỢNG SoC CHO MODUL PIN LITHIUM Chuyên ngành: Kỹ thuật điều khiển tự động hóa LUẬN VĂN THẠC SỸ KHOA HỌC KỸ THUẬT ĐIỀU... định SoC trở lên lớn Chính năm gần tác giả tập trung sử dụng lọc Kalman để ước lượng SoC lọc Kalman có khả ước lượng trạng thái điều kiện có nhiễu Vậy toán ước lượng trạng thái SoC cho Pin Lithium. .. Sau thời gian nghiên cứu, làm việc khẩn trương hướng dẫn tận tình giúp đỡ thầy giáo PGS.TS Nguyễn Văn Chí, luận văn với đề tài ? ?Nghiên cứu thuật tốn ước lượng SoC cho modul Pin Lithium? ?? hoàn