Giới thiệu dự án

Sự bùng nổ của thị trường xe điện (Electric Vehicle - EV) toàn cầu với tốc độ tăng trưởng kép hàng năm (CAGR) đạt hơn 21.7% là minh chứng rõ nét cho quá trình chuyển dịch tất yếu từ phương tiện sử dụng động cơ đốt trong (Internal Combustion Engine - ICE) sang năng lượng sạch nhằm giảm phát thải khí nhà kính và đáp ứng các tiêu chuẩn khí thải nghiêm ngặt như Euro 5, Euro 6 và QCVN 109:2021/BGTVT tại Việt Nam. Trong cấu trúc của một chiếc xe thuần điện (Battery Electric Vehicle - BEV), bộ pin cao áp (High-Voltage Traction Battery Pack) đóng vai trò là "trái tim", chiếm từ 30% đến 40% tổng chi phí sản xuất và trực tiếp quyết định đến tầm hoạt động, hiệu suất động lực học cũng như độ an toàn vận hành.

+-----------------------------------------------------------------------------------+
|                        TỔNG QUAN HỆ THỐNG PIN XE ĐIỆN                             |
+-----------------------------------------------------------------------------------+
|  +--------------------+    +----------------------+    +-----------------------+  |
|  | Samsung INR18650   | -> |  Mô hình 2-RC        | -> |  Bộ lọc Kalman mở     |  |
|  | Cell 3.6V - 2000mAh|    |  Thevenin Hộp xám    |    |  rộng (EKF) ước lượng |  |
|  +--------------------+    +----------------------+    +-----------------------+  |
|                                        |                                          |
|                                        v                                          |
|  +--------------------+    +----------------------+    +-----------------------+  |
|  | VinFast VF8 Eco    | <- |  Chu trình WLTP      | <- |  Tính toán Động lực   |  |
|  | Battery Pack ~82kWh|    |  Class 3 Dynamic Load|    |  học Xe & Năng lượng  |  |
|  +--------------------+    +----------------------+    +-----------------------+  |
+-----------------------------------------------------------------------------------+

Vấn đề cốt lõi đặt ra trong thiết kế hệ thống pin ô tô điện nằm ở các điểm nghẽn kỹ thuật:

  • Tính chất phi tuyến phức tạp của các phản ứng điện hóa bên trong cell pin dưới các điều kiện tải biến thiên liên tục.
  • Sai số tích lũy trong việc ước lượng trạng thái sạc (State of Charge - SOC) khi chỉ sử dụng phương pháp đếm Coulomb truyền thống.
  • Thách thức trong việc tối ưu hóa cách ghép nối hàng nghìn cell đơn lẻ để vừa đáp ứng điện áp danh định (~350V - 400V), vừa duy trì độ tin cậy và hạn chế tổn hao nhiệt do nội trở ($R_{int}$).

Đồ án tốt nghiệp chuyên ngành Công nghệ Kỹ thuật Ô tô tại Trường Đại học Sư phạm Kỹ thuật TP.HCM với đề tài "Thiết kế và mô phỏng hệ thống pin trên ô tô điện bằng phần mềm MATLAB/Simulink" do sinh viên Trần Nhật Duy và Phan Nguyễn Minh Nguyên thực hiện dưới sự hướng dẫn của TS. Nguyễn Mạnh Cường đã giải quyết trọn vẹn bài toán trên.

Mục tiêu cụ thể của dự án:

  1. Xây dựng mô hình toán học động lực học ô tô dựa trên thông số thực tế của dòng xe SUV điện VinFast VF8 Eco 2023 trên chu trình thử nghiệm tiêu chuẩn WLTP Class 3.
  2. Thiết kế và định cỡ chi tiết hệ thống pin từ cell đơn lẻ (Samsung INR18650-20R) đến cấu trúc pack pin hoàn chỉnh.
  3. Ứng dụng mô hình mạch tương đương Thevenin hai nhánh RC ($2$-RC Dual Polarization) kết hợp thuật toán lọc Kalman mở rộng (Extended Kalman Filter - EKF) để mô phỏng chính xác điện áp đầu cực ($V_t$) và ước lượng SOC.
  4. Kiểm chứng công suất tiêu thụ, tổn hao năng lượng và hiệu suất vận hành toàn hệ thống trên môi trường MATLAB/Simulink.

Phạm vi nghiên cứu tập trung vào mô hình hóa 1D/hộp xám của pin Lithium-ion NMC, kiểm nghiệm trên chu trình vận hành tiêu chuẩn, không bao gồm phân tích va chạm cơ học 3D hay mô phỏng nhiệt chuyên sâu CFD.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi lựa chọn công nghệ Lithium-ion NMC dạng trụ (Cylindrical Cell), các giải pháp lưu trữ năng lượng phổ biến trên phương tiện giao thông đã được tổng hợp và đánh giá qua ma trận kỹ thuật:

Tiêu chí kỹ thuật Pin Axit - Chì (Lead-Acid) Pin Niken-Kim loại (NiMH) Pin Lithium-ion (NMC 18650)
Mật độ năng lượng ($Wh/kg$) $30 - 50$ $60 - 120$ $150 - 250$
Hiệu suất sạc/xả ($%$) $70 - 85%$ $65 - 75%$ $95 - 99%$
Tuổi thọ chu kỳ (Cycles @ 80% DOD) $300 - 500$ $500 - 1,000$ $1,000 - 2,500$
Điện áp danh định cell ($V$) $2.0 - 2.1$ $1.2$ $3.6 - 3.7$
Tốc độ tự xả ($%/tháng$) $5 - 15%$ $20 - 30%$ $1 - 3%$
Hiệu ứng nhớ (Memory Effect) Không Rõ rệt Không có

Yêu cầu kỹ thuật hệ thống pin được phân loại theo mô hình MoSCoW:

  • Must have (Bắt buộc): Đáp ứng công suất kéo cực đại ($>150,kW$), dung lượng khả dụng trên $80,kWh$, sai số ước lượng SOC $<2%$, mạch ngắt bảo vệ quá áp ($4.2V$) và dưới áp ($2.5V$).
  • Should have (Nên có): Tích hợp thuật toán lọc Kalman mở rộng EKF tự hiệu chỉnh sai số trôi cảm biến dòng, mô phỏng tổn hao nhiệt Joule ($I^2R_0$).
  • Could have (Có thể mở rộng): Tích hợp ước lượng trạng thái sức khỏe (State of Health - SOH) dựa trên điện trở trong suy giảm theo chu kỳ.
  • Won't have (Chưa thực hiện): Quản lý nhiệt chủ động đa kênh chất lỏng 3D trong giai đoạn này.

Thiết kế hệ thống

Kiến trúc mô phỏng được xây dựng trên nguyên lý Model-Based Design (MBD), phân tách thành 4 phân hệ chính liên kết chặt chẽ qua bus dữ liệu Simulink:

graph TD
    A["Chu trình lái WLTP Class 3<br/>(Vận tốc vs Thời gian)"] --> B["Khối Động lực học Ô tô<br/>(VinFast VF8 Specs)"]
    B --> C["Công suất yêu cầu tại bánh xe<br/>(Ptong & Fcan)"]
    C --> D["Khối Quản lý Năng lượng & Phụ tải<br/>(Auxiliary Load: 1200W)"]
    D --> E["Mô hình Pack Pin 2-RC Thevenin<br/>(Cell Samsung 18650-20R)"]
    E --> F["Thuật toán EKF Ước lượng SOC<br/>(Look-up Table OCV-SOC)"]
    F -->|Feedback SOC & Terminal Voltage| E
    E --> G["Ghi nhận Dữ liệu & Đánh giá Hiệu năng<br/>(Tổn hao, Sụt áp, Wh/km)"]

Technology Stack triển khai:

  • Nền tảng mô phỏng: MATLAB R2023b / Simulink Version 10.8.
  • Toolboxes chuyên dụng: Simscape Electrical, Powertrain Blockset, Control System Toolbox, Signal Processing Toolbox.
  • Cơ sở dữ liệu đặc tính thực nghiệm: Bảng tra cứu đa chiều (2D Look-up Tables) chứa các thông số điện áp mạch hở $U_{OCV}(SOC)$, điện trở thuần $R_0(SOC)$, các nhánh phân cực $R_1, C_1$ (phản ứng truyền điện tích nhanh) và $R_2, C_2$ (khuếch tán nồng độ chậm) trích xuất từ dữ liệu thử nghiệm xung dòng HPPC (Hybrid Pulse Power Characterization).

Methodology

Quy trình nghiên cứu áp dụng chu trình chữ V (Automotive V-Model):

  1. Phân tích yêu cầu & Chuẩn hóa đầu vào: Trích xuất thông số khí động học và tự trọng của VinFast VF8 Eco ($m = 2600,kg$, $C_x = 0.28$, $A = 2.8,m^2$, $f = 0.015$).
  2. Thiết kế toán học & Xây dựng thuật toán: Lập trình phương trình trạng thái không gian rời rạc cho mạch $2$-RC và thuật toán ma trận hiệp phương sai EKF.
  3. Mô hình hóa khối trên Simulink: Tích hợp khối động lực học kéo, khối cell pin đơn lẻ và cấu hình module tổ hợp.
  4. Mô phỏng MIL (Model-in-the-Loop): Chạy thử nghiệm chu trình WLTP Class 3 với bước thời gian $T_s = 0.1,s$.

Implementation và kết quả

Development process

Quá trình tính toán động lực học xác định lực kéo tổng cộng tác dụng lên ô tô tại mỗi thời điểm:

$$F_{tong} = F_f + F_w + F_i + F_j$$

Trong đó:

  • $F_f = m \cdot g \cdot f \cdot \cos(\alpha)$: Lực cản lăn của lốp xe trên mặt đường.
  • $F_w = \frac{1}{2} \cdot \rho \cdot C_x \cdot A \cdot v^2$: Lực cản khí động học không khí.
  • $F_i = m \cdot g \cdot \sin(\alpha)$: Lực cản dốc mặt đường.
  • $F_j = \delta \cdot m \cdot \frac{dv}{dt}$: Lực cản quán tính khi xe tăng/giảm tốc ($\delta \approx 1.05$).

Công suất tức thời yêu cầu từ hệ thống pin:

$$P_{pack}(t) = \frac{F_{tong}(t) \cdot v(t)}{\eta_{trans} \cdot \eta_{motor}} + P_{aux}$$

function [P_batt, I_batt] = calc_battery_power(v, a, params)
    % Thông số đầu vào
    m = params.mass;          % 2600 kg
    f = params.rolling_coef;   % 0.015
    g = 9.81;
    rho = 1.225;              % Mật độ không khí kg/m3
    Cd = params.drag_coef;     % 0.28
    A = params.frontal_area;   % 2.8 m2
    eta = params.efficiency;   % 0.92
    P_aux = 1200;              % Phụ tải phụ trợ 1.2 kW

    % Lực thành phần
    F_roll = m * g * f;
    F_aero = 0.5 * rho * Cd * A * (v^2);
    F_accel = 1.05 * m * a;
    F_total = F_roll + F_aero + F_accel;

    % Công suất yêu cầu
    if F_total >= 0
        P_wheel = F_total * v;
        P_batt = (P_wheel / eta) + P_aux;
    else
        % Phanh tái sinh (Regenerative Braking) với hiệu suất thu hồi 65%
        P_wheel = F_total * v;
        P_batt = (P_wheel * 0.65 * eta) + P_aux;
    end
end

Mô hình mạch Thevenin 2-RC mô tả điện áp đầu cực tế bào pin:

$$V_t(t) = U_{OCV}(SOC) - I(t) \cdot R_0 - U_1(t) - U_2(t)$$

$$\frac{dU_1(t)}{dt} = -\frac{U_1(t)}{R_1 C_1} + \frac{I(t)}{C_1}, \quad \frac{dU_2(t)}{dt} = -\frac{U_2(t)}{R_2 C_2} + \frac{I(t)}{C_2}$$

% Thuật toán Lọc Kalman Mở Rộng (EKF) cập nhật trạng thái SOC
function [SOC_est, P_cov] = EKF_SOC_Update(I_k, V_meas, SOC_prior, P_prior, Q, R, dt, Cap_Ah, params_LUT)
    % 1. Dự đoán trạng thái (Time Update)
    SOC_pred = SOC_prior - (I_k * dt) / (Cap_Ah * 3600);
    P_pred = P_prior + Q;

    % 2. Tính toán đạo hàm riêng OCV theo SOC (Jacobian H)
    OCV_pred = interp1(params_LUT.SOC_grid, params_LUT.OCV_grid, SOC_pred, 'pchip');
    dOCV_dSOC = (interp1(params_LUT.SOC_grid, params_LUT.OCV_grid, SOC_pred + 0.001, 'pchip') - ...
                 interp1(params_LUT.SOC_grid, params_LUT.OCV_grid, SOC_pred - 0.001, 'pchip')) / 0.002;
    H = dOCV_dSOC;

    % 3. Cập nhật đo lường (Measurement Update)
    V_pred = OCV_pred - I_k * params_LUT.R0;
    y_tilde = V_meas - V_pred; % Residual
    S = H * P_pred * H' + R;
    K = (P_pred * H') / S;     % Kalman Gain

    % Trạng thái ước lượng tối ưu
    SOC_est = SOC_pred + K * y_tilde;
    P_cov = (1 - K * H) * P_pred;
end

Testing và validation

Hệ thống được kích hoạt chạy kiểm thử qua chu trình WLTP Class 3 với tổng thời gian $1800,s$, quãng đường di chuyển $23.26,km$, vận tốc cực đại đạt $131.3,km/h$:

Tham số kiểm chứng Giá trị thiết kế mục tiêu Kết quả mô phỏng đạt được Đánh giá sai số / Độ chính xác
Điện áp Pack Pin tối đa ($V$) $403.2,V$ ($96S$) $401.8,V$ Sai số $<0.35%$
Năng lượng tiêu hao ($E_{avg}$) $180 - 220,Wh/km$ $198.4,Wh/km$ Nằm trong dải danh định
Sai số điện áp đầu cực ($V_t$) RMSE $< 25,mV$ RMSE $= 11.8,mV$ Đạt chuẩn chính xác cao
Sai số ước lượng SOC (EKF) $\le 2.0%$ Tối đa $1.15%$ (TB $0.48%$) Khử hoàn toàn trôi tích phân
Nhiệt Joule tổn thất trung bình $< 3.5%$ công suất $2.14%$ tổng năng lượng Phù hợp tản nhiệt khí/lỏng

Đổi mới và đóng góp

  1. Mô hình hóa hộp xám Thevenin 2-RC đa tham số động: Khác biệt với các nghiên cứu sử dụng mô hình $R_{int}$ tĩnh chỉ tính đến nội trở thuần, đề tài đã tách biệt hiện tượng quá thế kích hoạt (Activation Polarization qua nhánh $R_1-C_1$) và quá thế nồng độ (Concentration Polarization qua nhánh $R_2-C_2$), tái hiện chính xác hiện tượng phục hồi điện áp (Voltage Relaxation) sau khi xe ngừng tăng tốc hoặc kết thúc quá trình phanh tái sinh.
  2. Khắc phục triệt để điểm yếu trôi sai số của phương pháp đếm Coulomb: Thuật toán EKF liên tục cập nhật ma trận hiệp phương sai nhiễu đo lường $R$ và nhiễu hệ thống $Q$, tự động nắn chỉnh giá trị SOC về tiệm cận thực tế dựa trên gradient điện áp mạch hở $dU_{OCV}/dSOC$, loại bỏ sai số tích lũy do độ nhạy cảm biến dòng điện.
  3. Định lượng năng lượng trực quan trên dòng xe thương mại tại Việt Nam: Cung cấp bộ thông số động học chuẩn xác của VinFast VF8 Eco kết hợp cùng cell pin 18650, đóng góp nguồn tài liệu tham khảo có giá trị cao cho các kỹ sư phát triển phần mềm nhúng BMS trong nước.

Ứng dụng thực tế và triển khai

Khả năng ứng dụng thực tiễn

Mô hình Simulink được đóng gói dạng module chuẩn hóa, sẵn sàng phục vụ các kịch bản triển khai:

  • Tích hợp kiểm thử phần mềm điều khiển BMS (Software-in-the-Loop - SIL): Nhúng trực tiếp mã nguồn C/C++ sinh tự động từ Simulink Coder vào các dòng vi điều khiển ô tô như STM32G474 hoặc TI TMS570 (chuẩn an toàn ISO 26262 ASIL-D).
  • Thiết kế định cỡ hệ thống lưu trữ tĩnh (BESS): Tái sử dụng mô hình cho các trạm sạc xe điện kết hợp pin lưu trữ năng lượng mặt trời (Solar-Storage EV Charging Stations).
  • Công cụ giảng dạy thực hành ảo: Phục vụ đào tạo chuyên sâu về kỹ thuật xe điện tại các trường đại học chuyên ngành Cơ khí Động lực và Kỹ thuật Điều khiển.

Phân tích hiệu quả kinh tế & Kỹ thuật

+-----------------------------------------------------------------------------------+
|                        LỘ TRÌNH TRIỂN KHAI VÀ HIỆU QUẢ KINH TẾ                   |
+-----------------------------------------------------------------------------------+
|  Giai đoạn 1 (Tháng 1-2): Mô phỏng MIL & Hiệu chỉnh tham số trên Simulink        |
|  Giai đoạn 2 (Tháng 3-4): Sinh mã nhúng tự động C/C++ & Kiểm thử SIL             |
|  Giai đoạn 3 (Tháng 5-6): Kiểm thử HIL trên Hardware Platform (TI/Infineon)       |
+-----------------------------------------------------------------------------------+
|  Lợi ích kinh tế:                                                                 |
|  - Rút ngắn 60% thời gian hiệu chỉnh thuật toán SOC trên xe thực tế.              |
|  - Tiết kiệm hàng trăm triệu đồng chi phí thử nghiệm phá hủy pack pin cao áp.     |
+-----------------------------------------------------------------------------------+

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Tham số mạch $R_0, R_1, C_1, R_2, C_2$ hiện được nội suy ở nhiệt độ danh định ($25^\circ C$), chưa phản ánh toàn diện hiện tượng suy giảm dung lượng ở dải nhiệt độ cực trị ($<0^\circ C$ hoặc $>45^\circ C$).
  • Mô hình coi tất cả các cell trong pack có mức độ lão hóa đồng đều, chưa mô phỏng hiện tượng lệch áp cell (Cell Imbalance) và thuật toán cân bằng pin chủ động/thụ động (Active/Passive Balancing).

Hướng phát triển tiếp theo

  • Ghép nối mô hình nhiệt khối tập trung (Lumped Thermal Model) vào phương trình $2$-RC để cập nhật điện trở $R(T)$ theo nhiệt độ tức thời.
  • Nghiên cứu thuật toán lọc Kalman thích nghi kép (Dual Extended Kalman Filter - DEKF) để ước lượng đồng thời cả SOC và điện trở suy thoái để tính toán SOH theo thời gian thực.
  • Ứng dụng mạng nơ-ron học sâu (Deep Neural Networks / LSTM) kết hợp vật lý (Physics-Informed Neural Networks - PINN) để tăng tốc độ xử lý trên chip nhúng biên.

Đối tượng hưởng lợi

  • Sinh viên ngành Kỹ thuật Ô tô / Điện tử: Nắm bắt phương pháp luận thiết kế động lực học phương tiện, hiểu sâu nguyên lý cấu tạo và đặc tính sạc/xả của công nghệ pin Lithium-ion hiện đại.
  • Kỹ sư phát triển phần mềm nhúng BMS: Tiếp cận trực tiếp cấu trúc thuật toán EKF hoàn chỉnh bằng mã nguồn MATLAB, sẵn sàng chuyển đổi mã sang C/C++ cho các dòng vi điều khiển chuyên dụng.
  • Doanh nghiệp sản xuất và tích hợp xe điện: Sở hữu bộ công cụ mô phỏng định cỡ pack pin với độ tin cậy cao, giảm thiểu chi phí R&D và số lượng nguyên mẫu vật lý cần chế tạo.
  • Nhà nghiên cứu học thuật: Cung cấp cơ sở dữ liệu mở về mô hình tương đương $2$-RC và chu trình WLTP Class 3 cho các công trình nghiên cứu nâng cao về quản lý năng lượng xe điện.

Câu hỏi thường gặp

1. Cấu hình phần cứng tối thiểu để triển khai mô phỏng này là gì?

Máy tính cần trang bị tối thiểu CPU 4 nhân (Intel Core i5 Gen 10 hoặc AMD Ryzen 5 tương đương), 8GB RAM (khuyến nghị 16GB) và cài đặt MATLAB/Simulink phiên bản từ R2021a trở lên có tích hợp Control System Toolbox và Simscape.

2. Tại sao mô hình Thevenin 2-RC lại được ưu tiên hơn mô hình 1-RC và mô hình điện hóa hộp trắng?

Mô hình 1-RC không phản ánh được hiện tượng hồi phục điện áp hai giai đoạn (ngắn hạn và dài hạn). Trong khi đó, mô hình điện hóa hộp trắng (P2D Model) đòi hỏi giải hệ phương trình đạo hàm riêng phi tuyến cực kỳ phức tạp, không thể chạy thời gian thực trên các bộ điều khiển nhúng BMS. Mô hình 2-RC là điểm cân bằng hoàn hảo giữa độ chính xác động học ($>98%$) và tải tính toán của vi xử lý.

3. Thuật toán EKF xử lý nhiễu cảm biến dòng điện thực tế như thế nào?

EKF sử dụng ma trận hiệp phương sai nhiễu đo lường $R$ để gán trọng số độ tin cậy cho cảm biến điện áp và ma trận $Q$ cho sai số tích phân dòng điện. Khi cảm biến dòng bị nhiễu hoặc trôi giá trị không (Zero-point drift), thuật toán sẽ tự động dựa vào độ chênh lệch giữa điện áp đo được $V_{meas}$ và $U_{OCV}(SOC)$ để hiệu chỉnh giá trị trạng thái SOC về điểm hội tụ thực.

4. Mô hình này có thể áp dụng cho các loại hóa học pin khác như LFP hay Natri-ion không?

Hoàn toàn có thể. Người dùng chỉ cần đo đạc đường đặc tính thực nghiệm $U_{OCV}(SOC)$ và bộ thông số trở kháng xung qua thử nghiệm HPPC của cell pin tương ứng, sau đó cập nhật lại các bảng Look-up Table trong phân hệ mô phỏng của Simulink.

5. Chi phí bản quyền và tính khả thi khi chuyển giao thương mại?

Giải pháp được thiết kế trên môi trường chuẩn công nghiệp MATLAB/Simulink, cho phép đóng gói thành dạng C-code độc lập thông qua Embedded Coder, loại bỏ hoàn toàn chi phí bản quyền runtime khi nạp vào phần cứng thương mại của xe điện.

Kết luận

Đồ án "Thiết kế và mô phỏng hệ thống pin trên ô tô điện bằng phần mềm MATLAB/Simulink" đã hoàn thành xuất sắc các mục tiêu nghiên cứu, chứng minh tính khả thi vượt trội của phương pháp thiết kế hướng mô hình (MBD) trong công nghiệp ô tô hiện đại. Bằng việc kết hợp chính xác giữa mô hình động lực học xe VinFast VF8, cấu trúc pack pin đa module từ cell Samsung 18650, mô hình tương đương Thevenin $2$-RC và thuật toán lọc Kalman mở rộng EKF, công trình đã giải quyết triệt để bài toán ước lượng trạng thái năng lượng với độ chính xác cao (sai số dưới $1.15%$).

Đây là tài liệu kỹ thuật có giá trị thực tiễn cao, đóng góp nền tảng vững chắc cho hoạt động R&D xe điện trong nước và mở ra hướng tiếp cận chuẩn mực cho các kỹ sư, nhà nghiên cứu và sinh viên đam mê công nghệ phương tiện xanh.