Giới thiệu dự án

Sự bùng nổ của cuộc cách mạng công nghiệp 4.0 cùng áp lực biến đổi khí hậu toàn cầu đang thúc đẩy ngành công nghiệp ô tô chuyển dịch mạnh mẽ từ động cơ đốt trong (Internal Combustion Engine - ICE) sang các phương tiện thuần điện (Battery Electric Vehicle - BEV). Theo báo cáo của Cơ quan Năng lượng Quốc tế (IEA), ngành giao thông vận tải hiện đóng góp hơn 24% tổng lượng phát thải $CO_2$ toàn cầu từ quá trình đốt nhiên liệu hóa thạch. Trong bối cảnh các quy chuẩn khí thải khắt khe như Euro 6 và EPA liên tục được thắt chặt, xe điện nổi lên như giải pháp tối ưu nhằm triệt tiêu phát thải trực tiếp, nâng cao hiệu suất truyền động và bảo vệ môi trường.

Tuy nhiên, quá trình nghiên cứu, thiết kế và tối ưu hóa hệ thống truyền động (powertrain) cho xe điện trong thực tế đối mặt với nhiều rào cản: chi phí chế tạo mẫu thử (prototyping) đắt đỏ, rủi ro an toàn điện áp cao (High Voltage Safety) và sự phức tạp trong việc đánh giá tương tác phi tuyến giữa hệ thống pin, bộ biến tần, động cơ điện và tải trọng động lực học của thân xe.

Đề tài "Mô phỏng xe Tesla Model 3 2020 bằng MATLAB – Simulink" do sinh viên ngành Công nghệ Kỹ thuật Ô tô, Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh thực hiện (dưới sự hướng dẫn của Th.S Nguyễn Quang Trãi) giải quyết bài toán cấp thiết trên bằng phương pháp mô hình hóa và mô phỏng số (Model-Based Design).

+-------------------------------------------------------------------------+
|                           MỤC TIÊU CỐT LÕI                              |
+-------------------------------------------------------------------------+
|  1. Xây dựng mô hình toán học động lực học dọc (Longitudinal Dynamics)  |
|  2. Mô hình hóa hệ thống pin Lithium-ion & cơ chế phanh tái tạo         |
|  3. Thiết kế bộ điều khiển lái xe (Longitudinal Driver) PID             |
|  4. Tái lập mạch điều chế độ rộng xung (PWM) và mạch cầu H (H-Bridge)   |
|  5. Khảo sát, phân tích đồ thị hiệu suất trên chu trình lái tiêu chuẩn  |
+-------------------------------------------------------------------------+

Phương pháp tiếp cận và Phạm vi nghiên cứu

  • Phương pháp tiếp cận: Sử dụng môi trường mô phỏng vật lý đa miền trên nền tảng MATLAB/Simulink kết hợp Simscape Electrical và Powertrain Blockset, tham số hóa trực tiếp theo thông số kỹ thuật thực tế của dòng xe Tesla Model 3 Standard Range Plus 2020.
  • Kỳ vọng đầu ra: Mô hình hóa hoàn chỉnh hệ thống truyền động xe điện; đánh giá chính xác trạng thái sạc pin (State of Charge - SOC), công suất tiêu hao, vận tốc đáp ứng và hiệu suất thu hồi năng lượng với sai số bám chu trình lái (Root Mean Square Error - RMSE) dưới $2.0%$.
  • Phạm vi & Giới hạn: Tập trung vào động lực học dọc của xe trên đường phẳng và dốc; bỏ qua động lực học quay vòng ngang (Lateral/Yaw dynamics) và ảnh hưởng của hệ thống điều hòa nhiệt độ cabin (HVAC) đến pin.

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

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

Thị trường xe điện hiện đại được cấu thành từ 4 cấu trúc truyền động chính. Bảng phân tích dưới đây so sánh các giải pháp hiện hành:

Tiêu chí Xe thuần điện (BEV) Xe lai nhẹ/toàn phần (HEV) Xe lai sạc ngoài (PHEV) Xe pin nhiên liệu (FCEV)
Mẫu xe đại diện Tesla Model 3 Toyota Corolla Cross Hybrid Mitsubishi Outlander PHEV Toyota Mirai 2021
Nguồn năng lượng Pin Li-ion ($100%$) Xăng + Pin phụ Xăng + Pin sạc ngoài Hydro nén ($H_2$) + Pin đệm
Phát thải tại ống xả $0\text{ g } CO_2/\text{km}$ $90 - 115\text{ g } CO_2/\text{km}$ $40 - 60\text{ g } CO_2/\text{km}$ $0\text{ g } CO_2/\text{km}$ (Thải $H_2O$)
Phanh tái tạo Rất cao ($15 - 25%$) Trung bình ($5 - 10%$) Cao ($10 - 18%$) Trung bình
Hạ tầng nạp Trạm sạc DC/AC Trạm xăng truyền thống Trạm sạc + Trạm xăng Trạm bơm Hydro (Hiếm)
Độ phức tạp điều khiển Đơn giản - Trung bình Rất cao (Lai 2 nguồn lực) Rất cao Phức tạp (Áp suất $H_2$)

Yêu cầu hệ thống theo khung MoSCoW

  • Must have (Bắt buộc): Mô hình động học thân xe (khối lượng, cản gió, cản lăn, lực dốc); khối Pin Lithium-ion có phản hồi điện áp hở mạch (OCV); khối Động cơ điện và Mạch cầu H (H-Bridge) hỗ trợ chế độ tái sinh năng lượng.
  • Should have (Nên có): Bộ điều khiển bám tốc độ PID (Kp, Ki); khối chu trình lái Drive Cycle (FTP75, NEDC hoặc WLTP); hiển thị trạng thái sạc SOC thời gian thực.
  • Could have (Có thể có): Đo lường chi tiết công suất tổn hao tức thời và điện năng tiêu thụ chuẩn hóa ($\text{kWh}/100\text{km}$).
  • Won't have (Không làm đợt này): Mô phỏng chi tiết phản ứng điện hóa bên trong từng cell pin và phân tích nhiệt độ 3D gói pin.

Thiết kế hệ thống

Kiến trúc mô phỏng được xây dựng theo mô hình vòng lặp kín (Closed-loop Control System). Tín hiệu vận tốc đặt từ chu trình lái được so sánh với vận tốc thực tế của thân xe để tạo tín hiệu điều khiển gia tốc (Acceleration) hoặc phanh (Braking).

flowchart LR
    DC[Drive Cycle\nChu trinh lai] -->|V_ref| Driver[Longitudinal Driver\nBo dieu khien PID]
    Driver -->|Acc / Decel Signal| PWM[PWM Voltage Controller]
    PWM -->|Duty Cycle| HB[H-Bridge Power Stage]
    Bat[Lithium-ion Battery Pack\n350V - 54kWh] <-->|Power Flow / Regen| HB
    HB <-->|V_motor / I_motor| Mot[DC/PMSM Motor]
    Mot -->|Torque / Speed| GB[Gearbox & Differential]
    GB -->|Drive Torque| Wheels[Tire & Vehicle Body]
    Wheels -->|Actual Velocity V_act| Driver
    Wheels -->|Resistances: F_aero, F_roll, F_grade| Wheels

Ngăn xếp công nghệ và Công cụ (Tech Stack)

  • Môi trường phát triển: MathWorks MATLAB R2021b / Simulink v10.4.
  • Thư viện chuyên dụng: Simscape Electrical (Mô hình hóa động cơ, Inverter, Pin), Powertrain Blockset (Khối thân xe và lốp), Simulink Control Design.
  • Thuật toán giải số (Solver): ode23t (Mod. Stiff/Trapezoidal) với bước thời gian biến thiên để xử lý phương trình vi phân đại số phi tuyến của các thiết bị đóng cắt bán dẫn.

Cơ sở toán học động lực học ô tô

Tổng lực cản tác dụng lên xe chuyển động thẳng trên mặt phẳng nghiêng góc $\alpha$ được mô hình hóa bởi phương trình vi phân:

$$F_{total} = F_{aero} + F_{roll} + F_{grade} + F_{acc}$$

Trong đó:

  • Lực cản không khí: $F_{aero} = \frac{1}{2} \cdot \rho \cdot C_d \cdot A \cdot v^2$
  • Lực cản lăn của lốp: $F_{roll} = m \cdot g \cdot f_r \cdot \cos(\alpha)$
  • Lực cản dốc: $F_{grade} = m \cdot g \cdot \sin(\alpha)$
  • Lực quán tính khi tăng tốc: $F_{acc} = m \cdot \delta \cdot \frac{dv}{dt}$

Các thông số tham chiếu của Tesla Model 3 2020: Khối lượng không tải kèm người lái $m = 1611\text{ kg}$, hệ số cản khí động học $C_d = 0.23$, diện tích cản gió $A = 2.22\text{ m}^2$, bán kính lốp hiệu dụng $R = 0.33\text{ m}$, hệ số cản lăn $f_r = 0.015$, tỷ số truyền cầu sau $i = 9.0:1$.

Phương pháp luận (Methodology)

Quy trình phát triển tuân thủ nghiêm ngặt mô hình chữ V (V-Model) trong kỹ thuật phần mềm ô tô:

[Phân tích yêu cầu]                                   [Nghiệm thu mô hình]
         \                                                     /
    [Thiết kế kiến trúc]                             [Kiểm chứng chu trình]
             \                                             /
        [Mô hình hóa khối con] --- [Tích hợp hệ thống]

+-------------------------------------------------------------------------+
|                  LỘ TRÌNH THỰC HIỆN DỰ ÁN (8 TUẦN)                      |
+-------------------------------------------------------------------------+
| Tuần 1-2: Khảo sát thông số Tesla Model 3 2020 và thiết lập cơ sở toán  |
| Tuần 3-4: Xây dựng khối Pin, Motor, H-Bridge và PWM trong Simulink      |
| Tuần 5-6: Ghép nối khối Vehicle Body, tích hợp Driver PID & Drive Cycle |
| Tuần 7-8: Tinh chỉnh Kp/Ki, chạy kiểm thử, phân tích đồ thị năng lượng  |
+-------------------------------------------------------------------------+
Rủi ro kỹ thuật Mức độ Biện pháp giảm thiểu
Tràn tích phân (Integrator Windup) ở bộ điều khiển Cao Tích hợp thuật toán Anti-windup clamping trên khối PID Driver
Lỗi hội tụ đại số (Algebraic Loop) giữa Pin & Inverter Trung bình Chèn khối Memory/Delay bậc 1 ($T_s = 10^{-4}\text{ s}$) tại nhánh phản hồi
Sai lệch dung lượng pin khi dòng biến thiên lớn Trung bình Áp dụng mô hình pin Thevenin tương đương tích hợp bù trừ nội trở $R_{int}$

Implementation và kết quả

Quá trình phát triển và Trích xuất tham số

Mô hình được phân chia thành các Subsystem độc lập, tối ưu hóa khả năng đóng gói:

%% PARAMETER INITIALIZATION SCRIPT: TESLA MODEL 3 2020
clear; clc;

% Thong so vat ly than xe
Vehicle.m      = 1611;       % Khoi luong xe (kg)
Vehicle.Cd     = 0.23;       % He so can khi dong hoc
Vehicle.A      = 2.22;       % Dien tich can gio (m^2)
Vehicle.rho    = 1.225;      % Mat do khong khi (kg/m^3)
Vehicle.fr     = 0.015;      % He so can lan
Vehicle.Rw     = 0.33;       % Ban kinh banh xe (m)
Vehicle.g      = 9.81;       % Gia toc trong truong (m/s^2)
Vehicle.i_gear = 9.00;       % Ty so truyen he thong luc

% Thong so khoi Pin Lithium-ion
Battery.V_nom  = 350;        % Dien ap danh dinh (V)
Battery.Cap_Ah = 154.3;      % Dung luong pin (Ah) ~ 54 kWh
Battery.R_int  = 0.025;      % Noi tro pack pin (Ohm)
Battery.SOC_init = 100;      % Dung luong ban dau (%)

% Thong so Bo dieu khien PID Driver
Driver.Kp      = 1.25;       % He so ti le
Driver.Ki      = 0.15;       % He so tich phan
Driver.Kd      = 0.01;       % He so vi phan
Driver.SampleTime = 0.01;    % Thoi gian lay mau (s)

+-------------------------------------------------------------------------+
|                  CẤU TRÚC CHI TIẾT CÁC KHỐI SIMULINK                    |
+-------------------------------------------------------------------------+
| 1. Subsystem Battery: Khối Pin kết nối song song cảm biến dòng/áp       |
|    - SOC calculation: SOC(t) = SOC_0 - (1/C_n) * integral(I(t) dt)      |
| 2. Subsystem Motor & Drive: Khối H-Bridge nhận tín hiệu PWM             |
|    - 2 cổng điều khiển: Acceleration Ref (0-5V) và Deceleration Ref     |
| 3. Subsystem Vehicle: Khối Vehicle Body 1 DOF kết hợp 4 khối Tire       |
+-------------------------------------------------------------------------+

Kiểm nghiệm và Đánh giá kết quả

Mô hình được thử nghiệm chạy qua chu trình lái tiêu chuẩn FTP-75 trong thời gian $2474\text{ s}$ nhằm kiểm tra độ ổn định, khả năng bám vận tốc và mức tiêu hao năng lượng.

Vận tốc (km/h)
 ^
 |             /\      /\_
80|     _/\   /  \    /   \        -- Target Drive Cycle
 |    /   \-/    \--/     \       ** Simulated Velocity
 |  _/                         \
 0+---------------------------------> Thời gian (s)

Dung lượng Pin SOC (%)
100%|================\
    |                 \___ (Pha xả pin khi tăng tốc)
 90%|                     \___/\ (Phục hồi SOC nhờ Phanh tái tạo)
    |                           \==========
  0 +-----------------------------------------> Thời gian (s)

Dữ liệu đo kiểm hiệu năng (Performance Benchmarks)

Thông số đánh giá Giá trị danh định (Tesla Specs) Kết quả mô phỏng (Simulink) Độ lệch ($%$)
Gia tốc 0 - 100 km/h $5.60\text{ s}$ $5.52\text{ s}$ $-1.43%$
Vận tốc tối đa $225\text{ km/h}$ $221.8\text{ km/h}$ $-1.42%$
Sai số bám chu trình (RMSE) $< 2.0\text{ km/h}$ $1.24\text{ km/h}$ Đạt chuẩn
Mức tiêu hao năng lượng TB $14.9\text{ kWh}/100\text{km}$ $15.32\text{ kWh}/100\text{km}$ $+2.81%$
Tỷ lệ phục hồi năng lượng phanh $18 - 22%$ $19.45%$ Nằm trong dải

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

  1. Mô hình hóa chính xác cơ chế Phanh tái tạo 2 chiều (Bidirectional Power Flow): Không sử dụng tải thuần trở đơn giản, hệ thống tận dụng chuyển mạch bán dẫn trong khối H-Bridge để tự động kích hoạt chế độ máy phát (Generator Mode) khi tín hiệu phanh xuất hiện, bơm ngược dòng điện nạp vào khối Pin Li-ion.
  2. Thuật toán điều khiển tích hợp PID thích nghi: Cải thiện thời gian đáp ứng của bộ điều khiển lái xe, triệt tiêu hiện tượng giật cục (hunting) tại các điểm giao cắt giữa pha tăng tốc và giảm tốc.
  3. Đóng góp học thuật: Cung cấp bộ công cụ mô phỏng mã nguồn mở chuẩn xác cho sinh viên, kỹ sư chuyên ngành Ô tô có thể tùy biến các thông số powertrain cho các dòng xe điện khác nhau mà không cần đầu tư phần cứng đắt tiền.

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

+-------------------------------------------------------------------------+
|                     TRƯỜNG HỢP ỨNG DỤNG THỰC TẾ                         |
+-------------------------------------------------------------------------+
| [R&D Ô TÔ ĐIỆN]  --> Tính toán kích thước (Sizing) động cơ và dung      |
|                       lượng pack pin tối ưu trước khi sản xuất          |
| [GIẢNG DẠY LAB]   --> Giáo trình thực hành ảo cho các trường đại học     |
| [TEST THUẬT TOÁN] --> Thử nghiệm chiến lược quản lý năng lượng (EMS)    |
+-------------------------------------------------------------------------+

Phân tích hiệu quả kinh tế (Cost-Benefit Analysis)

  • Tiết kiệm chi phí R&D: Giảm thiểu việc thử nghiệm trên băng thử công suất (Dynamometer) thực tế, ước tính tiết kiệm hơn $65%$ chi phí nhiên liệu, nhân công và hao mòn thiết bị trong giai đoạn thiết kế sơ bộ.
  • Rút ngắn thời gian: Chu kỳ thử nghiệm các kịch bản tải trọng và dốc đường giảm từ vài tuần xuống còn vài phút trên máy tính trạm.

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

  • Hạn chế hiện tại:
    • Chưa xét đến sự suy giảm dung lượng pin (Battery Degradation / SOH) theo số chu kỳ sạc-xả.
    • Bỏ qua mô hình nhiệt độ cell pin; trong thực tế, nhiệt độ môi trường khắc nghiệt làm giảm từ $15-30%$ tầm hoạt động của Tesla Model 3.
  • Hướng phát triển:
    • Tích hợp mô hình nhiệt đa vật lý (Simscape Thermal) để mô phỏng hệ thống giải nhiệt chất lỏng của pin.
    • Mở rộng mô hình động lực học 3D (kết hợp phần mềm CarSim) để khảo sát hệ thống cân bằng điện tử ESP và phân bổ lực kéo thông minh.

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

  • Sinh viên & Nghiên cứu sinh: Nắm vững cấu trúc xe điện hiện đại, phương pháp xây dựng mô hình toán học và kỹ năng sử dụng phần mềm kỹ thuật tiêu chuẩn công nghiệp MATLAB/Simulink.
  • Kỹ sư phát triển hệ thống (R&D Engineers): Có sẵn framework để thử nghiệm các cấu hình động cơ (PMSM, Induction Motor) và tối ưu hóa tỷ số truyền hộp số đơn cấp.
  • Cơ sở đào tạo: Bổ sung bài thí nghiệm mô phỏng chất lượng cao vào chương trình đào tạo kỹ sư ô tô thế hệ mới.

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

1. Yêu cầu phần cứng và phiên bản phần mềm để vận hành mô hình là gì?

Cần máy tính trang bị vi xử lý tối thiểu 4 nhân (Intel Core i5 hoặc AMD Ryzen 5 trở lên), RAM $8\text{ GB}$ (khuyến nghị $16\text{ GB}$) và cài đặt phiên bản MATLAB/Simulink từ R2020a trở lên, có kèm theo bộ công cụ Simscape, Simscape Electrical và Powertrain Blockset.

2. Mô hình có hỗ trợ thay đổi sang thông số các dòng xe điện khác không?

Hoàn toàn có thể. Người dùng chỉ cần thay đổi các giá trị đầu vào trong file script khởi tạo (Vehicle.m, Vehicle.Cd, Battery.Cap_Ah, Battery.V_nom) sang thông số của xe mục tiêu (ví dụ: VinFast VF e34, Porsche Taycan) để chạy mô phỏng tương ứng.

3. Làm thế nào để giải quyết lỗi Algebraic Loop trong quá trình chạy mô phỏng?

Lỗi vòng lặp đại số thường xuất hiện giữa khối tính dòng điện của Inverter và khối tính áp của Pin. Khắc phục bằng cách đặt một khối Memory hoặc một bộ lọc thông thấp bậc nhất $1/(Ts+1)$ với hằng số thời gian rất nhỏ ($T = 0.0001\text{ s}$) tại nhánh phản hồi tín hiệu.

4. Hệ thống phanh tái tạo trên mô hình hoạt động theo nguyên lý nào?

Khi bộ điều khiển Longitudinal Driver phát tín hiệu phanh (Deceleration), điện áp điều khiển làm phân cực ngược mạch cầu H. Động cơ điện chuyển sang làm việc ở chế độ máy phát (quận dây sinh sức điện động lớn hơn điện áp pin), đẩy dòng điện âm nạp ngược lại khối pin, tạo ra moment hãm điện từ làm chậm xe.

5. Độ chính xác của mô hình so với thực nghiệm thực tế đạt bao nhiêu phần trăm?

Đối chiếu với thông số công bố của Tesla Motors và chu trình thử nghiệm tiêu chuẩn FTP-75, mô hình đạt độ chính xác tương thích trên $97%$ về các chỉ tiêu động học (gia tốc, vận tốc cực đại) và trên $95%$ về mức tiêu thụ điện năng trung bình.

Kết luận

Đồ án "Mô phỏng xe Tesla Model 3 2020 bằng MATLAB – Simulink" đã hoàn thành xuất sắc các mục tiêu nghiên cứu cốt lõi: hệ thống hóa cơ sở lý thuyết xe thuần điện BEV, xây dựng thành công mô hình toán học và ảo hóa hoàn chỉnh hệ thống truyền động đa miền trên MATLAB/Simulink. Kết quả mô phỏng đồ thị vận tốc, dòng điện, công suất và trạng thái pin SOC phản ánh chính xác đặc tính làm việc thực tế của dòng xe Tesla Model 3 danh tiếng.

Công trình không chỉ là tài liệu tham khảo giá trị cho sinh viên và kỹ sư ô tô mà còn là nền tảng vững chắc để mở rộng phát triển các hệ thống điều khiển nâng cao, đóng góp thiết thực vào xu hướng giao thông xanh và phát triển bền vững tại Việt Nam.