Bối cảnh và vấn đề nghiên cứu

Nguồn năng lượng hóa thạch phục vụ cho các ngành công nghiệp nói chung và ngành giao thông vận tải đường bộ nói riêng đang ngày càng cạn kiệt. Trên các phương tiện giao thông truyền thống, năng lượng từ nhiên liệu thực sự được chuyển hóa để dẫn động phương tiện chỉ chiếm dưới 20% (thường dao động trong khoảng 10% đến 16%), phần lớn còn lại bị tổn hao dưới dạng nhiệt lượng và khắc phục các lực cản. Một lượng lớn động năng liên tục bị tiêu hao qua hệ thống treo khi xe di chuyển trên bề mặt đường gồ ghề; các giảm chấn thủy lực truyền thống triệt tiêu dao động này bằng cách hấp thụ và chuyển hóa chúng thành nhiệt năng phát tán ra môi trường.

Theo số liệu thống kê từ Viện Chiến lược và Phát triển Giao thông Vận tải (TDSI) giai đoạn 2011–2016, hoạt động giao thông vận tải tại Việt Nam tiêu thụ khoảng 30% tổng nhu cầu năng lượng quốc gia với tốc độ tăng trưởng 10%/năm. Trong đó, vận tải đường bộ là lĩnh vực tiêu thụ nhiên liệu lớn nhất (chiếm 68% tổng tiêu thụ của ngành) và chịu trách nhiệm cho khoảng 85% tổng lượng phát thải khí nhà kính (tương đương khoảng 30 triệu tấn $\text{CO}_2$ mỗi năm). Tính đến năm 2023, theo thống kê của Bộ Giao thông Vận tải, Việt Nam có khoảng 45,5 triệu phương tiện đang lưu hành. Việc thu hồi năng lượng cơ học lãng phí từ hệ thống giảm xóc để chuyển đổi thành điện năng là hướng nghiên cứu có giá trị thực tiễn nhằm nâng cao hiệu quả sử dụng năng lượng và cắt giảm phát thải.

Mục tiêu và nhiệm vụ nghiên cứu của đề tài gồm:

  1. Nghiên cứu tổng quan nguyên lý hoạt động và các kết cấu thu hồi năng lượng từ hệ thống treo trên phương tiện giao thông.
  2. Thiết lập mô hình toán học và xây dựng mô hình mô phỏng hệ thống thu hồi năng lượng điện từ giảm xóc bằng phần mềm MATLAB/Simulink.
  3. Khảo sát, phân tích đáp ứng động lực học và hiệu suất phát điện của hệ thống dưới các kịch bản kích thích mặt đường khác nhau (tín hiệu bước và tín hiệu điều hòa ở các dải vận tốc).
  4. Đánh giá khả năng tích trữ điện năng vào mạch tụ điện và đề xuất cơ sở lý thuyết cho việc phát triển mẫu thử nghiệm thực tế.

Đối tượng nghiên cứu của khóa luận là hệ thống thu hồi năng lượng điện từ tích hợp trên giảm xóc sau của xe máy (sử dụng thông số kỹ thuật thực tế của dòng xe Sirius 2007 110cc). Phạm vi nghiên cứu giới hạn ở việc mô hình hóa dao động một bậc tự do (1-DOF), ứng dụng hiện tượng cảm ứng điện từ trên máy phát điện tuyến tính trực tiếp và khảo sát quá trình nạp/xả năng lượng qua mạch RC.

Cơ sở lý thuyết và phương pháp

Đề tài ứng dụng định luật cảm ứng điện từ Faraday để xác định điện áp sinh ra khi cuộn dây và hệ thống nam châm có chuyển động tương đối:

$$U = B \cdot l \cdot v$$

Trong đó:

  • $U$: Điện áp cảm ứng sinh ra (V).
  • $B$: Mật độ từ thông qua cuộn dây (Tesla - T).
  • $l$: Tổng chiều dài dây dẫn nằm trong vùng từ trường (m).
  • $v$: Vận tốc chuyển động tương đối giữa cuộn dây và từ trường (m/s).

Để tối ưu hóa mật độ từ trường $B$, tác giả khảo sát lý thuyết mảng nam châm Halbach (Halbach Array) – cấu hình đổi hướng cực từ 90 độ liên tiếp nhằm tập trung mật độ từ thông về một phía khe hở không khí. Các công thức xác định thông số hình học cuộn dây gồm:

$$\text{Chiều dài một vòng dây: } l = N \cdot D_a \cdot \pi$$

$$\text{Số vòng dây: } N = \frac{\pi \cdot A_c}{2\sqrt{3} \cdot A_w}$$

$$\text{Bề rộng cuộn dây: } H = d \cdot N$$

Trong đó $D_a$ là đường kính trung bình của cuộn dây, $A_c$ là diện tích mặt cắt ngang của cuộn dây, $A_w$ là diện tích mặt cắt của dây quấn, $H$ là bề rộng cuộn dây và $d$ là đường kính sợi dây đồng.

Phương pháp nghiên cứu được triển khai bao gồm:

  • Phương pháp phân tích - tổng hợp lý thuyết: Thu thập tài liệu chuyên khảo, bài báo khoa học quốc tế về giảm chấn tái sinh năng lượng, phân loại các cơ cấu dẫn động trực tiếp (Direct Drive) và gián tiếp (Indirect Drive), các giải pháp chỉnh lưu chuyển động (Motion Rectifier) cơ học và thủy lực.
  • Phương pháp mô hình hóa và mô phỏng số: Thiết lập hệ phương trình vi phân mô tả dao động hệ thống khối lượng - lò xo - giảm chấn (Mass-Spring-Damper - MSD) 1 bậc tự do trong miền thời gian; sử dụng công cụ MATLAB/Simulink để xây dựng sơ đồ khối mô phỏng hệ thống thu hồi điện từ tuyến tính, phân tích đáp ứng tần số và đáp ứng quá độ.
  • Nguồn dữ liệu: Dữ liệu cấu hình vật liệu từ tính; thông số kỹ thuật thực nghiệm của xe máy Sirius 2007 110cc (khối lượng toàn tải $m = 120\text{ kg}$, hệ số cản giảm chấn $c = 200\text{ N}\cdot\text{s/m}$); dữ liệu tham chiếu từ các công trình nghiên cứu quốc tế của P. Zhang, Z. Wang, T. Ai-Yafeai, F. Silva.

Thiết kế và triển khai

Hệ thống được thiết kế theo nguyên lý máy phát điện tuyến tính trực tiếp (Linear Electromagnetic Regenerative Shock Absorber - LERA) gắn đồng trục với bộ giảm chấn.

Khối chức năng trong Simulink Chi tiết thiết kế và thông số kỹ thuật
Khối kích thích mặt đường (Road Profile) Sử dụng Signal Editor, Pulse Generator và Transport Delay để tạo 2 dạng tín hiệu: tín hiệu xung bước (mô phỏng gờ giảm tốc cao 5 cm) và tín hiệu điều hòa (sóng sin tuần hoàn).
Khối dao động cơ học (MSD Block) Mô hình hóa hệ thống treo 1 bậc tự do với khối lượng $m = 120\text{ kg}$, hệ số cản $c = 200\text{ N}\cdot\text{s/m}$, nhận kích thích mấp mô để xuất tín hiệu vị trí $x(t)$ và vận tốc $\dot{x}(t)$.
Khối máy phát tuyến tính (Linear Harvester) Chuyển đổi vận tốc tương đối $\dot{x}(t)$ thành sức điện động (EMF) dựa trên mật độ từ thông mảng Halbach 2 lớp ($B \approx 0{,}8\text{ T}$).
Khối lưu trữ điện năng (RC Circuit) Mạch chỉnh lưu kết hợp tụ điện tích trữ dung lượng cao $C = 10\text{ F}$, mô phỏng hai pha đóng/ngắt công tắc để khảo sát quá trình sạc và xả điện.
Kích thích mặt đường (Road Profile)
Khối dao động cơ học (MSD Block: m = 120 kg, c = 200 N.s/m)
Máy phát tuyến tính LERA (Mảng Halbach 2 lớp, B = 0.8 T)
Mạch nạp / xả và lưu trữ điện năng (Tụ điện C = 10 F)

Nội dung chính theo từng chương

Chương 1: Giới thiệu tổng quan (Introduction)

Chương 1 làm rõ tính cấp thiết của việc tái tạo năng lượng từ hệ thống treo trên phương tiện giao thông. Tác giả trình bày ý nghĩa khoa học trong việc mô hình hóa máy phát điện tuyến tính và ý nghĩa thực tiễn đối với mục tiêu giảm tiêu thụ nhiên liệu. Điểm mới của đề tài là ứng dụng mảng nam châm tuyến tính để khai thác trực tiếp dao động thẳng đứng do mặt đường gây ra mà không qua cụm biến đổi chuyển động phức tạp.

Phần tổng quan tài liệu (Literature Review) phân tích các nghiên cứu tiêu biểu:

  • Nghiên cứu của Z. Wang (2022) về bộ thu hồi năng lượng áp điện không tiếp xúc (PEH) kết hợp trục vít me - đai ốc, tạo ra công suất cực đại $24{,}28\text{ W}$ ở tốc độ $60\text{ km/h}$ trên đường ngẫu nhiên và đạt $3346\text{ W}$ ở tốc độ $30\text{ km/h}$ trên đường có xung kích thích khi xe không tải.
  • Nghiên cứu của T. Ai-Yafeai (2021) mô hình hóa nửa xe tích hợp lớp vật liệu áp điện trên lò xo giảm xóc trước và sau.
  • Nghiên cứu của F. Silva (2013) về bộ truyền động lực tái sinh (RFA) hoạt động trong dải tần $0{,}5 - 20\text{ Hz}$, cho công suất thu hồi tối đa đạt $984{,}4\text{ W}$ ở tần số kích thích cao nhất.

Chương 2: Kết cấu và nguyên lý hoạt động của các hệ thống thu hồi năng lượng (Construction & Operating Principles of Energy Harvesting Systems)

Chương này đi sâu vào phân tích cấu tạo hệ thống treo truyền thống và phân loại các công nghệ thu hồi năng lượng:

  • Hệ thống treo xe máy: Phân tích đặc tính các loại lò xo (lò xo trụ, lò xo côn, lò xo xoắn) và cấu tạo giảm chấn thủy lực dạng vỏ (shell damper) gồm ống giảm xóc (fork tube), thanh trượt (slider), van tiết lưu (damper valve), dầu giảm chấn (fork oil) và phớt chắn dầu.
  • Chế độ dẫn động trực tiếp (Direct Drive Mode - LERA): Khảo sát độ từ cảm của các nhóm vật liệu (Diamagnetic, Paramagnetic, Ferromagnetic, Ferrimagnetic). Đánh giá 4 cấu hình sắp xếp nam châm: nam châm dọc trục có đệm, nam châm hướng tâm có đệm, nam châm hướng tâm thuần túy và mảng Halbach. Dữ liệu từ nghiên cứu của P. Zhang chỉ ra rằng mảng Halbach 2 lớp cho mật độ từ thông trung bình đạt $0{,}8\text{ Tesla}$, vượt trội so với cấu hình dọc trục có đệm ($0{,}4375\text{ Tesla}$).
  • Chế độ dẫn động gián tiếp (Indirect Drive Mode): Phân tích cơ chế hoạt động của trục vít me bi (Ball Screw), cơ cấu thanh răng - bánh răng (Rack and Pinion), tấm áp điện (Piezoelectric Plate) và hệ thống thủy lực. Điểm nhấn của chương là nguyên lý bộ chỉnh lưu chuyển động cơ học (Mechanical Motion Rectifier) sử dụng hai ly hợp con lăn (roller clutches) và bộ chỉnh lưu thủy lực (Hydraulic Motion Rectifier) với việc tối ưu hóa van một chiều (giúp tăng hiệu suất chuyển đổi công suất từ 12% lên đến 40%).
Đặc tính so sánh Dẫn động trực tiếp (Direct Drive) Dẫn động gián tiếp: Thanh răng - Bánh răng Dẫn động gián tiếp: Vít me bi Dẫn động gián tiếp: Tấm áp điện
Nguyên lý hoạt động Máy phát điện từ tuyến tính chuyển đổi chuyển động nén/nhả trực tiếp thành điện năng. Thanh răng biến chuyển động tịnh tiến thành chuyển động quay để kéo máy phát điện. Trục vít me bi biến đổi dao động tịnh tiến thành chuyển động quay ma sát thấp. Vật liệu áp điện tạo ra điện áp khi bị biến dạng/nén cơ học.
Công suất tạo ra Trung bình (có thể nâng cao theo kích thước) Khá cao ($50 - 100\text{ W}$) Rất cao ($100\text{ W} - 3\text{ kW}$) Thấp ($1 - 5\text{ W}$)
Hiệu suất chuyển đổi Cao ($60 - 80%$) Trung bình ($50 - 60%$) Cao (lên đến $70%$) Trung bình
Độ phức tạp cơ khí Thấp (ít chi tiết chuyển động) Trung bình Trung bình đến Cao Thấp
Độ êm dịu chuyển động Tốt, dễ điều chỉnh Có thể giảm do ma sát ăn khớp Tốt do ma sát lăn thấp Phù hợp với xe chú trọng độ êm dịu
Chi phí chế tạo Trung bình Thấp Khá đắt Thấp

Dựa trên bảng so sánh, nhóm tác giả lựa chọn phương án dẫn động trực tiếp (LERA) để triển khai mô hình hóa vì cấu trúc nhỏ gọn, độ tin cậy cơ học cao và phù hợp không gian lắp ráp trên giảm xóc xe máy.

Chương 3: Mô hình hóa, mô phỏng và phân tích kết quả (Modeling, Simulation and Results Analysis)

Chương 3 trình bày các bước xây dựng mô hình tính toán trên MATLAB/Simulink:

  • Giới thiệu tổng quan hệ thống phần mềm MATLAB, cú pháp dòng lệnh, cấu trúc file .m, .mat, và các công cụ chuyên ngành như Simulink, Stateflow, Fast Fourier Transform (FFT).
  • Thiết lập kịch bản mô phỏng xung bước: Xe vượt qua gờ giảm tốc cao $5\text{ cm}$, khảo sát ở hai dải tốc độ $6\text{ m/s}$ ($21{,}6\text{ km/h}$) và $24\text{ m/s}$ ($86{,}4\text{ km/h}$). Thời gian chạy qua gờ được đặt là 20 giây kèm thời gian nghỉ 8 giây giữa hai xung mấp mô.
  • Kết quả mô phỏng động lực học khối lượng - lò xo - giảm chấn (MSD): Vị trí dịch chuyển và vận tốc tương đối của giảm chấn dao động mạnh khi gặp xung kích thích, sau đó nhanh chóng ổn định nhờ lực cản giảm chấn ($c = 200\text{ N}\cdot\text{s/m}$).
  • Kết quả phát điện: Điện áp ngõ ra tỷ lệ thuận với vận tốc dao động của giảm xóc. Ở vận tốc di chuyển $24\text{ m/s}$, biên độ điện áp đỉnh và công suất tức thời sinh ra cao hơn rõ rệt so với ở vận tốc $6\text{ m/s}$.
  • Mô phỏng quá trình nạp/xả tụ điện: Tích hợp mạch sạc với tụ điện dung lượng $C = 10\text{ F}$. Đồ thị quá trình nạp cho thấy điện áp trên tụ tăng dần theo thời gian khi hệ thống treo dao động và duy trì mức năng lượng dự trữ để cung cấp ngược lại cho tải tiêu thụ trong chu kỳ xả.

Chương 4: Kết luận và khuyến nghị (Conclusion and Recommendation)

Tổng kết toàn bộ kết quả mô hình hóa, khẳng định tính đúng đắn của phương pháp tiếp cận và đề xuất các giải pháp kỹ thuật cần hoàn thiện ở giai đoạn chế tạo mẫu thực nghiệm.

Kết quả và đóng góp

Khóa luận đạt được các kết quả cụ thể sau:

  1. Thiết lập mô hình động lực học hoàn chỉnh: Xây dựng thành công sơ đồ khối trên MATLAB/Simulink mô phỏng tương tác giữa kích thích mặt đường, đáp ứng dao động của cụm giảm xóc xe máy và đặc tính phát điện tuyến tính.
  2. Chứng minh hiệu quả cấu hình mảng Halbach: Khẳng định việc bố trí nam châm mảng Halbach hai lớp đạt mật độ từ thông $0{,}8\text{ Tesla}$, giúp tăng cường điện áp cảm ứng mà không làm tăng kích thước cơ cấu.
  3. Đánh giá định lượng tiềm năng kinh tế - môi trường: Với giả định 45,5 triệu xe máy tại Việt Nam vận hành trung bình 4 giờ/ngày và mỗi xe thu hồi công suất dao động ở mức $0{,}2375\text{ W}$, hệ thống có khả năng sản xuất liên tục $45{,}225\text{ kWh}$ điện năng. Tổng sản lượng điện tái tạo ước tính đạt $43{,}125\text{ triệu kWh/năm}$, mang lại giá trị kinh tế hơn 55 tỷ đồng (tính theo đơn giá điện kinh doanh $3.500\text{ VNĐ/kWh}$) và góp phần giảm phát thải khí $\text{CO}_2$ từ các nhà máy nhiệt điện.

Hạn chế và hướng nghiên cứu tiếp

Khóa luận chỉ ra một số điểm hạn chế trong phạm vi thực hiện:

  • Nghiên cứu dừng lại ở mức độ mô hình hóa lý thuyết và mô phỏng số học (Numerical Simulation) trên phần mềm MATLAB/Simulink, chưa tiến hành gia công chế tạo mẫu thử nghiệm vật lý (hardware prototype) để kiểm chứng trên băng thử dao động thực tế.
  • Mô hình dao động mới chỉ khảo sát hệ dao động một bậc tự do (1-DOF) trên giảm xóc sau xe máy, chưa xét đến tương tác động lực học toàn xe (Full-car model) hoặc các yếu tố nhiệt độ môi trường ảnh hưởng đến độ nhớt của dầu giảm chấn và từ tính của nam châm vĩnh cửu.

Hướng nghiên cứu tiếp theo được tác giả đề xuất:

  • Thiết kế, gia công chi tiết cơ khí và chế tạo mẫu máy phát điện tuyến tính LERA thực tế dựa trên các thông số đã tính toán.
  • Tiến hành thực nghiệm trên máy thử giảm chấn chuyên dụng để đo đạc chính xác công suất phát, hiệu suất chuyển đổi và độ bền cơ học.
  • Nghiên cứu tích hợp mạch điều khiển quản lý năng lượng thông minh (Power Management Circuit) nhằm ổn định điện áp đầu ra trước khi nạp vào ắc-quy hoặc siêu tụ điện của phương tiện.

Giá trị tham khảo

Đồ án tốt nghiệp là tài liệu tham khảo hữu ích cho sinh viên, học viên cao học và kỹ sư thuộc các chuyên ngành Kỹ thuật Ô tô, Kỹ thuật Cơ khí Động lực, Cơ điện tử và Kỹ thuật Năng lượng tái tạo.

Các nội dung có giá trị ứng dụng cao bao gồm:

  • Phương pháp luận xây dựng mô hình toán học và sơ đồ khối mô phỏng hệ dao động cơ - điện từ tích hợp trong môi trường MATLAB/Simulink.
  • Bảng tổng hợp so sánh chi tiết giữa các cơ cấu thu hồi năng lượng trực tiếp và gián tiếp (vít me bi, thanh răng, áp điện, thủy lực).
  • Dữ liệu tính toán kích thước cuộn dây và cấu hình từ trường mảng Halbach ứng dụng cho các máy phát điện tuyến tính kích thước nhỏ gọn.

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

1. Đề tài lựa chọn phương án dẫn động trực tiếp (Direct Drive) hay gián tiếp (Indirect Drive) để mô phỏng?
Tác giả lựa chọn phương án dẫn động trực tiếp (Direct Drive - LERA). Lý do là kết cấu này có độ phức tạp cơ khí thấp, không cần bộ truyền động trung gian biến đổi chuyển động tịnh tiến thành chuyển động quay, giảm thiểu tổn hao ma sát và dễ dàng tích hợp vào không gian giảm xóc xe máy.

2. Cấu hình nam châm mảng Halbach có ưu thế gì so với cấu hình bố trí nam châm thông thường?
Theo số liệu nghiên cứu của P. Zhang được trích dẫn trong đồ án, cấu hình mảng Halbach 2 lớp tập trung mật độ từ thông về phía khe hở không khí chứa cuộn dây, đạt mật độ từ thông trung bình $0{,}8\text{ Tesla}$, cao hơn gần gấp đôi so với cấu hình nam châm dọc trục có đệm ($0{,}4375\text{ Tesla}$).

3. Thông số phương tiện thực tế nào được sử dụng trong mô hình mô phỏng Simulink?
Nhóm tác giả sử dụng thông số kỹ thuật của xe máy Yamaha Sirius 2007 110cc, bao gồm khối lượng ướt toàn tải $m = 120\text{ kg}$ và hệ số cản của bộ giảm chấn $c = 200\text{ N}\cdot\text{s/m}$.

4. Kịch bản kích thích mặt đường dạng bước (Step Signal) được thiết lập với các thông số nào?
Kịch bản mô phỏng gờ giảm tốc có độ cao $5\text{ cm}$, khảo sát ở hai mức vận tốc xe chạy là $6\text{ m/s}$ ($21{,}6\text{ km/h}$) và $24\text{ m/s}$ ($86{,}4\text{ km/h}$), thời gian vượt qua gờ là 20 giây và thời gian nghỉ giữa các gờ là 8 giây.

5. Mạch lưu trữ năng lượng trong mô phỏng sử dụng loại linh kiện nào?
Mô hình sử dụng mạch RC tích hợp tụ điện có điện dung $C = 10\text{ F}$ để nạp năng lượng sinh ra từ bộ máy phát tuyến tính và thực hiện xả năng lượng khi đóng/ngắt khóa chuyển mạch.

Kết luận

Khóa luận tốt nghiệp của sinh viên Nguyễn Trần Bình và Đỗ Việt Lâm đã hoàn thành việc nghiên cứu, tính toán và mô phỏng hệ thống thu hồi năng lượng điện từ giảm chấn xe máy sử dụng máy phát điện tuyến tính mảng Halbach trên phần mềm MATLAB/Simulink. Kết quả mô phỏng khẳng định tính khả thi của việc chuyển đổi dao động cơ học do mấp mô mặt đường thành điện năng hữu ích, với điện áp đầu ra tỷ lệ thuận với tốc độ dao động của hệ thống treo. Đề tài cung cấp cơ sở lý thuyết và mô hình tính toán hoàn chỉnh, mở ra tiền đề thuận lợi cho các nghiên cứu chế tạo mẫu thử nghiệm thực tế và ứng dụng công nghệ tái sinh năng lượng trên các phương tiện giao thông tại Việt Nam.