Tổng quan luận án

Vận tải đường sắt đóng vai trò quan trọng trong việc đáp ứng nhu cầu gia tăng cả về số lượng lẫn chất lượng của hệ thống giao thông vận tải phục vụ phát triển kinh tế - xã hội. Trong cơ cấu chi phí sản xuất của các xí nghiệp đầu máy đường sắt, chi phí năng lượng chiếm tỷ trọng chi phối, dao động từ 60% đến 65%. Sự biến động theo chiều hướng gia tăng của giá năng lượng trực tiếp làm tăng chi phí vận hành, gây áp lực lên hiệu quả khai thác toàn ngành. Việc nâng cao hiệu suất sử dụng năng lượng không chỉ phụ thuộc vào giải pháp công nghệ chế tạo phương tiện mà còn chịu tác động từ phương pháp vận dụng và điều khiển đoàn tàu trên tuyến. Việc vận hành tối ưu thông qua thao tác của tài xế và kế hoạch chạy tàu có thể tiết giảm tới 10% chi phí nhiên liệu mà không đòi hỏi kinh phí đầu tư lớn cho hạ tầng hay thay mới đầu máy.

Khoảng trống nghiên cứu được chỉ ra xuất phát từ thực tiễn và lý thuyết: Các công trình nghiên cứu kinh điển về điều khiển tối ưu chuyển động đoàn tàu thường xem xét một đoàn tàu đơn lẻ vận hành tự do trong khu gian với các mô hình động lực học giản lược (xem đoàn tàu như một chất điểm tập trung tại trọng tâm, lực kéo và lực cản là các biến liên tục hoặc không đổi). Khi áp dụng vào thực tế mạng lưới đường sắt—đặc biệt là đường sắt đơn có mật độ vận hành nhiều đoàn tàu, yêu cầu tác nghiệp tránh vượt tại các ga trung gian, hoặc hệ thống đường sắt đôi sử dụng phân khu đóng đường cố định/di động—các mô hình đơn lẻ không giải quyết được bài toán phối hợp tối ưu tổng thể. Bên cạnh đó, điều kiện hình học tuyến đường tại Việt Nam có nhiều đoạn đổi dốc ngắn với chiều dài nhỏ hơn chiều dài đoàn tàu, khiến giả thiết chất điểm gây sai số trong tính toán lực cản phụ. Đồng thời, các hệ số trong công thức tính lực cản chính thực tế có sự sai khác so với thông số do nhà chế tạo cung cấp.

Mục đích nghiên cứu của luận án bao gồm các mục tiêu cụ thể:

  1. Xây dựng mô hình toán học của bài toán điều khiển tối ưu nhiều đoàn tàu trên tuyến theo chỉ tiêu tiết kiệm năng lượng.
  2. Xây dựng phương pháp xác định các giá trị lực cản phụ có tính đến chiều dài và vị trí hiện tại của đoàn tàu trên tuyến.
  3. Xây dựng các thuật toán xác định vị trí của điểm chuyển chế độ tối ưu khi đoàn tàu đi qua các phân đoạn có độ dốc lên và dốc xuống.
  4. Áp dụng nguyên lý cực đại của Pontryagin để giải bài toán điều khiển tối ưu cục bộ cho từng đoàn tàu.
  5. Xây dựng đồ thị đặc tính điều khiển tối ưu đoàn tàu theo chỉ tiêu tiết kiệm năng lượng dựa trên các thông số đo ghi thực tế (tốc độ, vị trí, công suất) và các thông số cho trước (trắc dọc, trắc ngang của tuyến, mác tàu, công lệnh tốc độ).

Đối tượng nghiên cứu của luận án là đoàn tàu sử dụng đầu máy diesel truyền động điện cùng đoàn toa xe, kết hợp với đoạn đường sắt mà đoàn tàu vận hành qua trên tuyến Hà Nội – Lào Cai.

Phạm vi nghiên cứu bao gồm: Thiết lập mô hình toán học giải tích, phát triển thuật toán điều khiển tối ưu, đánh giá chất lượng thuật toán thông qua mô phỏng số và đối chứng với dữ liệu chạy tàu thực tế do tài xế điều khiển trên tuyến đường sắt Yên Viên – Lào Cai, đồng thời khảo sát mở rộng về trắc dọc trên tuyến Thống Nhất đoạn Vinh – Nha Trang.

Tổng quan tài liệu và vị trí của luận án

Tổng quan tài liệu trong luận án phân loại các giải pháp tiết kiệm năng lượng trong giao thông đường sắt theo các nhóm giải pháp công nghệ và quản lý vận hành. Theo các tài liệu tổng hợp chuyên ngành, năng lượng tiêu thụ của đoàn tàu gồm bốn thành phần chính: năng lượng để gia tốc và vượt dốc; năng lượng thắng lực cản chuyển động (ma sát và khí động học); năng lượng phụ trợ phục vụ tiện nghi (điều hòa, quạt gió, chiếu sáng trên tàu khách); và năng lượng tổn hao trên hệ thống truyền động (động cơ điện kéo, bộ biến đổi, cơ cấu cơ khí).

Nhóm hướng nghiên cứu Tác giả, công trình và mốc thời gian tiêu biểu Nội dung và kết quả chính
Nghiên cứu trong nước Tổng cục Đường sắt Việt Nam & Bộ môn Đầu máy ĐH Giao thông vận tải (thập niên 1980) Xây dựng định mức tiêu hao nhiên liệu cho đầu máy diesel sau giai đoạn chuyển đổi sức kéo diesel (bắt đầu từ 1963 ở miền Nam, 1968 ở miền Bắc).
Công ty Tư vấn Đường sắt (Dự án năm 2007) Đề xuất giải pháp khoa học công nghệ: hợp lý hóa vận chuyển, lập biểu đồ chạy tàu, đào tạo tài xế, thay thế bộ điều tốc cơ khí - thủy lực bằng điều tốc điện tử.
Mô hình động lực học và các chế độ điều khiển cơ bản Kokotovic và Singh (1970-x) Xây dựng mô hình điều khiển tối ưu phi tuyến bậc hai để giảm tiêu hao năng lượng thông qua kiểm soát dòng điện kéo.
Milroy (1980-x) Xác lập đặc tuyến chuyển động tối ưu cho hành trình ngắn gồm 3 giai đoạn: gia tốc cực đại, chạy đà (thả trôi), và phanh hãm cực đại.
Asnis và cộng sự; Lee và cộng sự; P. Howlett Bổ sung giai đoạn chạy ổn tốc cho hành trình dài, hoàn thiện mô hình đặc tuyến 4 giai đoạn: gia tốc cực đại, ổn tốc, chạy đà, và phanh với lực cực đại.
Điều khiển tối ưu theo nguyên lý Pontryagin cho 1 đoàn tàu P. Howlett và cộng sự; Golovitcher; Cheng và P. Howlett Sử dụng nguyên lý cực đại Pontryagin để tìm điều kiện cần của đặc tuyến liên tục; chia tuyến thành các đoạn có độ dốc đơn để tìm điểm chuyển đổi cục bộ; chứng minh giai đoạn chạy đều bị gián đoạn khi gặp dốc.
Khmelnitsky (2000) Nghiên cứu bài toán điều khiển tối ưu với độ dốc biến đổi, hiệu suất lực kéo biến đổi và giới hạn tốc độ tùy ý.
Albrecht và cộng sự Phát triển phương pháp số dựa trên nguyên lý cực đại Pontryagin, ứng dụng vào các hệ thống trợ giúp lái tàu thương mại như Freightmister, Energymister2.
Điều khiển nhiều đoàn tàu và tối ưu hóa biểu đồ chạy tàu Wang và Goverde Điều khiển tối ưu hai đoàn tàu chạy cùng hướng trên một khu đoạn theo nguyên tắc "làn sóng xanh", cân đối giữa chi phí nhiên liệu và thời gian hành trình.
Yang và cộng sự Thiết kế cấu hình tốc độ cho mạng đường sắt đô thị nhiều đoàn tàu với đặc tuyến 3 giai đoạn; đồng bộ hóa lịch trình để tận dụng năng lượng hãm tái sinh.
Liu và Golovitcher Phát triển phương pháp giải tích xác định điểm chuyển trạng thái điều khiển với hiệu suất lực kéo không đổi và giới hạn tốc độ tùy ý.
Chen; HU, Yin-Ping, Chao HU; Li và cộng sự Ứng dụng thuật toán di truyền (GA), thuật toán tối ưu bầy đàn (PSO), thuật toán chia đôi để tối ưu hóa biểu đồ chạy tàu và tìm điểm chuyển chế độ tối ưu.
Albrecht và Oettich; Wang và cộng sự Thiết kế thuật toán giảm thời gian chờ tại các ga trung gian, áp dụng phương pháp phổ giả và quy hoạch tuyến tính hỗn hợp cho nhóm đoàn tàu.
Li và Lo; Gupta và cộng sự; Su và cộng sự Tối ưu hóa phân phối thời gian bổ sung và đồng bộ hóa hành động tăng tốc - hãm giữa các đoàn tàu để tận dụng năng lượng hãm tái sinh trên hệ thống tàu điện ngầm.

Khoảng trống mà luận án lựa chọn giải quyết tập trung vào ba điểm mấu chốt:

  • Phần lớn các nghiên cứu áp dụng nguyên lý Pontryagin chỉ khảo sát một đoàn tàu chạy tự do trong khu gian. Luận án mở rộng việc thiết lập hàm mục tiêu và giải pháp điều khiển cho nhiều đoàn tàu trên tuyến, bao gồm tình huống hai đoàn tàu tránh nhau tại ga trung gian trên tuyến đường đơn và các đoàn tàu kế tiếp nhau trong hệ thống đóng đường cố định hoặc di động.
  • Khắc phục hạn chế của mô hình chất điểm truyền thống bằng cách đề xuất phương pháp tính toán lực cản phụ có xét đến sự phân bố chiều dài đoàn tàu khi vận hành qua các đoạn đổi dốc liên tiếp có chiều dài ngắn hơn chiều dài đoàn tàu.
  • Kết hợp lý thuyết điều khiển tối ưu với việc hiệu chuẩn thực nghiệm các hệ số lực cản chính ($a, b, c$) thông qua thuật toán di truyền dựa trên số liệu đo đạc thực tế tại Việt Nam, thay vì sử dụng hoàn toàn thông số lý thuyết định mức.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án sử dụng nền tảng lý thuyết điều khiển tối ưu kết hợp với phương pháp giải tích và thuật toán số để giải quyết bài toán chuyển động của đoàn tàu:

Khung lý thuyết và mô hình toán học

Cơ sở lý thuyết trung tâm là Nguyên lý cực đại Pontryagin áp dụng cho hệ thống phi tuyến, không Autonom (do phương trình vi phân chuyển động phụ thuộc trực tiếp vào biến không gian quãng đường $s$).

Phương trình vi phân chuyển động của đoàn tàu lấy quãng đường $s$ làm biến số độc lập: $$\frac{dv}{ds} = \frac{1}{\gamma v} \left[ \frac{F - B_T}{P + Q} - \omega_0(v) - \omega_d(s) \right]$$ $$t(s) = \int_{s_d}^{s} \frac{1}{v} ds$$

Trong đó:

  • $v(s)$ và $t(s)$ là vận tốc và thời gian chuyển động tại vị trí $s$.
  • $P$ là trọng lượng đầu máy, $Q$ là tổng trọng lượng các toa xe.
  • $F$ là lực kéo của đầu máy, $B_T$ là lực hãm đoàn tàu.
  • $\omega_0(v)$ là lực cản chính đơn vị của đoàn tàu, biểu diễn dưới dạng hàm bậc hai theo vận tốc: $\omega_0(v) = a + bv + cv^2$.
  • $\omega_d(s)$ là lực cản phụ đơn vị do độ dốc và bán kính đường cong.
  • $\gamma$ là hệ số quán tính của đoàn tàu.

Phương pháp hệ số nhân Lagrange

Tiêu chuẩn tối ưu là cực tiểu hóa tổng chi phí năng lượng kéo $A_e = \int_{s_d}^{s_c} \frac{F}{\eta} ds$ đồng thời thỏa mãn ràng buộc đẳng thức về thời gian chạy tàu quy định trên khu gian $T_c$. Luận án sử dụng phương pháp nhân tử Lagrange để đưa bài toán tối ưu có ràng buộc thời gian về bài toán tối ưu không ràng buộc với hàm mục tiêu mở rộng: $$A' = \int_{s_d}^{s_c} \left[ u_f f_m(v) + \frac{\lambda}{v} \right] ds \to \min$$ Trong đó $u_f = \frac{F(v)}{F_{max}(v)}$ và $u_b = \frac{B_T(v)}{B_{Tmax}(v)}$ là các biến điều khiển tương đối ($0 \le u_f, u_b \le 1$), $\lambda$ là hệ số nhân Lagrange điều chỉnh tương ứng với thời gian chạy tàu, $f_m(v)$ và $b_m(v)$ là gia tốc kéo và gia tốc hãm đơn vị lớn nhất.

Thiết lập hàm Hamilton và biến phụ trợ

Hàm Hamilton mở rộng của hệ thống được xác lập cùng với hệ phương trình đồng trạng thái. Bằng cách định nghĩa biến phụ trợ $p(s) = \frac{\psi_1(s)}{\psi_0 v(s)}$, điều kiện cực đại của hàm Hamilton được chuyển đổi thành bài toán phân tích giá trị của hàm $p(s)$ dọc theo hành trình, từ đó xác lập các điều kiện chuyển đổi trạng thái giữa các chế độ vận hành.

Nguồn tư liệu và dữ liệu thực tế

  • Dữ liệu tuyến đường: Trắc dọc, trắc ngang, vị trí ga, độ dốc, bán kính đường cong và công lệnh tốc độ trên tuyến đường sắt Hà Nội – Lào Cai (khảo sát thực tế đoạn Yên Viên – Lào Cai) và tuyến Thống Nhất (đoạn Vinh – Nha Trang).
  • Dữ liệu phương tiện: Thông số kỹ thuật của đầu máy diesel truyền động điện, thành phần đoàn xe, chiều dài đoàn tàu tiêu chuẩn (khảo sát trường hợp $L = 500\text{ m}$).
  • Dữ liệu đo ghi thực nghiệm: Tốc độ, vị trí, công suất động cơ diesel được thu thập từ thiết bị đo ghi chuyên dụng lắp đặt trên cabin đầu máy, phục vụ việc nhận dạng và hiệu chỉnh các hệ số lực cản $a, b, c$ thông qua thuật toán di truyền.

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

Chương 1: Tổng quan về vấn đề điều khiển tối ưu năng lượng trong giao thông đường sắt

Chương 1 hệ thống hóa các nhân tố ảnh hưởng đến mức tiêu hao nhiên liệu của đầu máy và tổng quan các phương pháp giải quyết bài toán tiết kiệm năng lượng.

  • Các nhân tố ảnh hưởng được phân tích thành 4 nhóm chính: trạng thái kỹ thuật và phương pháp vận dụng động cơ diesel - máy phát - động cơ điện kéo; hoạt động của hệ thống phụ trợ (làm mát, bôi trơn, cấp nhiên liệu, nén khí); phương pháp vận dụng và thao tác lái tàu của tài xế; công tác kiểm tra, sửa chữa phương tiện.
  • Phân tích kỹ thuật vận dụng chỉ ra rằng thao tác của tài xế (lựa chọn vị trí tay máy phù hợp với công suất động cơ, lợi dụng động năng khi lên dốc, xuống dốc và chạy đà) quyết định trực tiếp tới chi phí nhiên liệu đơn vị ($kg/10^4\text{ t.km}$). Các hoạt động hợp lý hóa vận dụng có tiềm năng giảm tới 10% chi phí nhiên liệu. Nghiên cứu thực nghiệm được trích dẫn cho thấy cùng một thành phần đoàn tàu và khối lượng kéo, chi phí nhiên liệu đơn vị giữa các điều kiện vận hành có thể chênh lệch nhau tới 11%.
  • Tổng kết tình hình nghiên cứu quốc tế từ thập niên 1970 đến nay: phân tích các mô hình đặc tuyến chuyển động 3 giai đoạn (Milroy) và 4 giai đoạn (Asnis, Howlett); các ứng dụng nguyên lý cực đại Pontryagin giải bài toán điều khiển liên tục và gián đoạn; các phương pháp giải tích và thuật toán lai (quy hoạch động, quy hoạch tuyến tính, thuật toán di truyền GA, thuật toán bầy đàn PSO, tối ưu đàn kiến) trong xây dựng biểu đồ chạy tàu và tận dụng năng lượng hãm tái sinh.
  • Xác định nhiệm vụ nghiên cứu: xây dựng hàm mục tiêu tổng quát cho nhiều đoàn tàu trên tuyến; ứng dụng nguyên lý cực đại Pontryagin cho bài toán cục bộ của từng tàu; xử lý ảnh hưởng của chiều dài đoàn tàu qua các đoạn dốc ngắn; và hiệu chỉnh các hệ số lực cản đơn vị bằng thực nghiệm.

Chương 2: Phân tích, xây dựng bài toán điều khiển tối ưu chuyển động đoàn tàu

Chương 2 là phần nội dung lý thuyết cốt lõi của luận án, tập trung mô hình hóa toán học và phát triển các phương pháp giải.

  • Tiêu chuẩn tối ưu và mô hình hóa đối tượng: Xác lập hàm chi phí vận chuyển tổng quát $C_{vc} = C_n A_n + C_w A_w + C_t A_t + C_g T_c$. Khi tích hợp chi phí hao mòn vào chi phí nhiên liệu, hàm mục tiêu rút gọn về dạng $C_{vc} = C_n A_n + C_g T_c$. Thiết lập hệ phương trình vi phân chuyển động không Autonom với biến độc lập $s$ cùng các ràng buộc biên: tốc độ $0 < v \le V_{max}(s)$, lực kéo $0 \le F \le F_{max}$, dòng điện động cơ $0 \le I_d \le I_{dmax}(v)$, lực hãm $0 \le B_T \le B_{Tmax}(v)$.
  • Phân tích hàm Hamilton và 5 chế độ điều khiển tối ưu: Dựa vào biến phụ trợ $p(s)$, luận án chứng minh điều kiện cực đại của hàm Hamilton dẫn đến 5 chế độ vận hành cơ bản:
    1. $p < 0$: $u_f = 0, u_b = 1$ — Chế độ hãm lớn nhất (ký hiệu: TM).
    2. $p = 0$: $u_f = 0, u_b \in [0; 1]$ — Chế độ hãm duy trì tốc độ (ký hiệu: CT), giữ vận tốc bằng vận tốc giới hạn $V_{max}$.
    3. $p \in (0; 1)$: $u_f = 0, u_b = 0$ — Chế độ chạy đà / quán tính (ký hiệu: QT).
    4. $p = 1$: $u_f \in [0; 1], u_b = 0$ — Chế độ ổn tốc (ký hiệu: C), vận tốc đoàn tàu được duy trì tại giá trị không đổi $V_{on}$ phụ thuộc vào hệ số Lagrange $\lambda$ và các thông số $a, b, c$ theo biểu thức $v^3 \omega_0'(v) = \lambda$.
    5. $p > 1$: $u_f = 1, u_b = 0$ — Chế độ kéo lớn nhất (ký hiệu: LK).
  • Cấu trúc đặc tính tối ưu trên các dạng địa hình:
    • Đối với khu gian phẳng ($v_d < v_c$): Đặc tuyến gồm 4 chế độ tuần tự LK $\to$ C $\to$ QT $\to$ TM. Trên các khu đoạn ngắn, chế độ C bị triệt tiêu, đặc tuyến trở thành LK $\to$ QT $\to$ TM. Trường hợp $v_c > 0$, chế độ TM được loại bỏ (LK $\to$ C $\to$ QT). Khi $v_d > v_c$, cấu trúc là QT $\to$ C $\to$ QT $\to$ TM.
    • Phân tích bước nhảy của hàm $p(s)$ và tính toán chuyển đổi chế độ thông qua việc bảo toàn giá trị hàm Hamilton $H_{n-1} = H_n$ qua các phân đoạn có độ dốc thay đổi.
  • Bài toán điều khiển tối ưu nhiều đoàn tàu trên tuyến:
    • Xét trường hợp đường sắt đơn điển hình ở Việt Nam: hai đoàn tàu chạy ngược chiều phải tránh nhau tại ga trung gian C tại thời điểm quy định $T_3$. Hành trình từ ga A đến ga D được phân chia thành hai khoảng thời gian $T_1 = T_3 - T_1 = T_a$ và $T_2 = T_8 - T_3 = T_b$. Phiếm hàm mục tiêu được gán hai hệ số nhân Lagrange $\lambda_1$ (cho đoạn từ ga A đến ga C) và $\lambda_2$ (cho đoạn từ ga C đến ga D).
    • Khi vượt qua ranh giới ga trung gian $S_1$, hàm $p(s)$ xuất hiện bước nhảy bậc từ $p^-$ sang $p^+$, dẫn đến sự thay đổi vận tốc ổn định từ $V_{on1}$ sang $V_{on2}$.
    • Luận án phân tích chi tiết điều kiện kiểm tra tính tối ưu khi có độ dốc lên hoặc dốc xuống tại lân cận điểm $S_1$ cho cả hai trường hợp $V_{on2} > V_{on1}$ và $V_{on2} < V_{on1}$ thông qua biểu thức quan hệ giữa các giá trị vận tốc và lực cản.
    • Mở rộng bài toán cho hành trình qua $k$ ga trung gian với vector hệ số nhân $\lambda_1, \lambda_2, \dots, \lambda_{k+1}$.
    • Phân tích mở rộng cho tuyến đường đôi chạy cùng hướng: trong hệ thống đóng đường cố định, các hệ số $\lambda_k$ của tàu đi sau phụ thuộc vào thời điểm vượt phân khu của tàu đi trước; trong hệ thống đóng đường di động, vị trí phân ranh $S_k(t)$ và hệ số $\lambda_k(t)$ là các biến số phụ thuộc động vào tọa độ tức thời của tàu đi trước.
  • Xét ảnh hưởng của chiều dài đoàn tàu:
    • Khảo sát thực tế trên tuyến Hà Nội – Lào Cai và Vinh – Nha Trang cho thấy tồn tại nhiều đoạn dốc ngắn liên tiếp có chiều dài nhỏ hơn chiều dài đoàn tàu ($L = 500\text{ m}$).
    • Đạo hàm hàm Hamilton theo quãng đường: $$\frac{dH}{ds} = \psi_1 \left( \frac{F - B_T}{P + Q} \right)' - \psi_1 \omega_0'(v) - \psi_1 g'(s)$$
    • Khi xem đoàn tàu là chất điểm, trong mỗi đoạn dốc $g'(s) = 0$ nên $H = \text{const}$. Tuy nhiên, khi xét đoàn tàu có chiều dài phân bố khối lượng đều dọc theo thân tàu, trọng tâm đoàn tàu đi qua nhiều độ dốc cục bộ cùng lúc khiến giá trị độ dốc tương đương biến thiên liên tục ($g'(s) \ne 0$), dẫn đến $H$ không còn là hằng số. Luận án đề xuất phương pháp tích phân lực cản phụ dọc theo chiều dài $L$ để hiệu chỉnh lực cản dốc tác dụng lên đoàn tàu.

Chương 3: Xây dựng các thuật toán điều khiển; mô phỏng thuật toán và kiểm nghiệm thuật toán trên một tuyến cụ thể

Theo tóm tắt tổng thể trong cấu trúc luận án, nội dung Chương 3 giải quyết các vấn đề thực nghiệm và ứng dụng:

  • Xây dựng các thuật toán điều khiển tối ưu hóa điểm chuyển chế độ dựa trên hệ thống công thức giải tích và nguyên lý cực đại đã thiết lập ở Chương 2.
  • Xây dựng chương trình tính toán hiệu chuẩn các hệ số $a, b, c$ trong hàm lực cản chính $\omega_0(v) = a + bv + cv^2$ bằng thuật toán di truyền (Genetic Algorithm - GA) kết hợp với dữ liệu đo thực nghiệm từ thiết bị ghi thông số vận hành trên đoàn tàu.
  • Mô phỏng thuật toán trên máy tính với điều kiện biên và trắc dọc cụ thể của tuyến đường sắt Yên Viên – Lào Cai.
  • Kiểm nghiệm thực tế: So sánh mức tiêu hao nhiên liệu và quỹ đạo vận tốc của phương pháp điều khiển đề xuất với kết quả chạy tàu thực tế do tài xế điều khiển trong cùng điều kiện tải trọng và thời gian quy định.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và khoa học

  1. Thiết lập mô hình toán học và phương pháp giải bài toán điều khiển tối ưu cho nhiều đoàn tàu vận hành trên tuyến theo chỉ tiêu tiết kiệm năng lượng, giải quyết sự tương tác giữa các đoàn tàu tránh nhau tại ga trung gian trên tuyến đường đơn thông qua việc phân đoạn hệ số nhân Lagrange.
  2. Xây dựng phương pháp phân tích bài toán điều khiển tối ưu hai đoàn tàu chạy cùng hướng trong cả hai cấu hình điều khiển: hệ thống đóng đường tự động cố định và hệ thống đóng đường di động.
  3. Đề xuất phương pháp tính toán lực cản phụ có xét đến sự phân bố chiều dài của đoàn tàu khi vận hành qua các đoạn đổi dốc liên tiếp có chiều dài ngắn hơn chiều dài đoàn tàu, chỉ ra sự biến thiên của hàm Hamilton khi loại bỏ giả thiết chất điểm.
  4. Xây dựng thuật toán giải tích và điều kiện biên để xác định chính xác tọa độ điểm chuyển đổi giữa 5 chế độ vận hành tối ưu (kéo lớn nhất, ổn tốc, chạy đà, hãm duy trì tốc độ, hãm lớn nhất) trên các biên dạng trắc dọc phức tạp (dốc lên, dốc xuống).

Đóng góp mới về mặt thực tiễn và giải pháp kỹ thuật

  1. Đề xuất phương pháp hiệu chuẩn các hệ số $a, b, c$ của hàm lực cản chính đơn vị $\omega_0(v)$ cho đoàn tàu cụ thể trên tuyến đường sắt Việt Nam bằng cách kết hợp thuật toán di truyền với dữ liệu đo ghi thực tế từ đầu máy, thay thế các hệ số định mức cũ.
  2. Cung cấp cơ sở khoa học để thiết kế và chế tạo hệ thống thiết bị hỗ trợ lái tàu lắp đặt trực tiếp trên cabin đầu máy diesel, hiển thị khuyến nghị về vị trí tay máy và thời điểm chuyển đổi chế độ vận hành nhằm tối ưu hóa chi phí nhiên liệu mà vẫn đảm bảo an toàn và biểu đồ chạy tàu.
  3. Cung cấp tài liệu tham khảo khoa học phục vụ công tác giảng dạy, nghiên cứu chuyên ngành Đầu máy - Toa xe, Kỹ thuật điều khiển và Vận tải đường sắt.

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

Dựa trên các giả thiết và phạm vi nghiên cứu được văn bản công bố, các hạn chế kỹ thuật của luận án gồm:

  • Mô hình tính toán coi hệ số hiệu dụng $\eta$ của cơ cấu truyền động là hằng số không đổi trong toàn bộ dải làm việc, chưa xét đến sự biến thiên của hiệu suất theo nhiệt độ và các chế độ phụ tải động cơ.
  • Giả thiết lực kéo và lực hãm là các biến liên tục trong phạm vi giới hạn ($0 \le u_f, u_b \le 1$), trong khi trên đầu máy diesel truyền động điện thực tế, lực kéo bị gián đoạn theo từng vị trí tay máy (nấc tay máy) tương ứng với các mức công suất động cơ diesel.
  • Do khối lượng tính toán lớn, bài toán tối ưu chưa thể giải trực tiếp theo thời gian thực trên cabin mà phải phân tách thành hai giai đoạn: tính toán ngoại tuyến (offline) các quỹ đạo đặc tính mẫu, sau đó hệ thống tự động/hỗ trợ lái tàu thực hiện việc bù sai lệch trong quá trình vận hành thực tế.
  • Khảo sát thực nghiệm mới chỉ thực hiện trên đối tượng đầu máy diesel truyền động điện trên tuyến Hà Nội – Lào Cai; việc mở rộng cho các loại sức kéo khác (như đầu máy điện, đoàn tàu điện động lực phân tán EMU) chưa được thực hiện chi tiết trong khuôn khổ nghiên cứu này.

Giá trị tham khảo

Kết quả của luận án có giá trị tham khảo chuyên môn đối với các nhóm đối tượng:

  • Nghiên cứu sinh và học viên sau đại học chuyên ngành Kỹ thuật Giao thông, Điều khiển Tự động: Tham khảo phương pháp vận dụng Nguyên lý cực đại Pontryagin cho hệ thống không Autonom và kỹ thuật sử dụng hệ số nhân Lagrange gián đoạn để giải bài toán tối ưu có ràng buộc thời gian và vận tốc.
  • Giảng viên và cán bộ nghiên cứu tại các trường đại học khối kỹ thuật giao thông: Sử dụng nội dung Chương 2 làm tài liệu giảng dạy chuyên sâu về sức kéo đường sắt, động lực học đoàn tàu và tối ưu hóa vận tải.
  • Kỹ sư vận dụng đầu máy và chuyên viên kỹ thuật Tổng công ty Đường sắt Việt Nam: Tham khảo thuật toán và cấu trúc đặc tuyến 5 chế độ vận hành để xây dựng sổ tay hướng dẫn lái tàu tiết kiệm năng lượng, lập kế hoạch phân bổ thời gian biểu đồ chạy tàu và thiết kế hệ thống màn hình hỗ trợ tài xế trên cabin.

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

1. Tại sao hàm Hamilton trong bài toán điều khiển đoàn tàu lại không phải là một hằng số theo quãng đường khi xét chiều dài đoàn tàu?

Theo chứng minh ở Chương 2, đạo hàm của hàm Hamilton theo quãng đường $s$ có dạng $\frac{dH}{ds} = \psi_1 \left( \frac{F - B_T}{P + Q} \right)' - \psi_1 \omega_0'(v) - \psi_1 g'(s)$. Khi coi đoàn tàu là chất điểm chuyển động trong một đoạn dốc cố định, độ dốc là hằng số nên $g'(s) = 0$, dẫn đến $H = \text{const}$. Tuy nhiên, khi xét đoàn tàu có chiều dài phân bố thực tế (ví dụ $L = 500\text{ m}$) đi qua khu vực có các đoạn dốc ngắn liên tiếp nhỏ hơn chiều dài tàu, trọng tâm và các phần của đoàn tàu đồng thời chịu nhiều độ dốc khác nhau. Khi đó $g'(s) \ne 0$, khiến giá trị $H$ biến thiên liên tục theo quãng đường và đòi hỏi phải xét lại điều kiện cực đại của hàm Hamilton.

2. Biến phụ trợ $p(s)$ được định nghĩa như thế nào và xác định những chế độ vận hành nào?

Biến phụ trợ được định nghĩa qua tỷ số giữa biến liên hợp $\psi_1(s)$ và vận tốc: $p(s) = \frac{\psi_1(s)}{\psi_0 v(s)}$. Căn cứ vào giá trị của $p(s)$, điều kiện cực đại của hàm Hamilton xác lập 5 chế độ điều khiển tối ưu:

  • $p < 0$: Chế độ hãm lớn nhất (TM) với $u_f = 0, u_b = 1$.
  • $p = 0$: Chế độ hãm duy trì tốc độ (CT) với $u_f = 0, u_b \in [0; 1]$.
  • $p \in (0; 1)$: Chế độ chạy đà / quán tính (QT) với $u_f = 0, u_b = 0$.
  • $p = 1$: Chế độ ổn tốc (C) với $u_f \in [0; 1], u_b = 0$, duy trì vận tốc ổn định $V_{on}$.
  • $p > 1$: Chế độ kéo lớn nhất (LK) với $u_f = 1, u_b = 0$.

3. Vận tốc ổn tốc $V_{on}$ phụ thuộc vào các yếu tố nào và được xác định ra sao?

Tại chế độ ổn tốc ($p = 1$), đạo hàm $\frac{dp}{ds} = 0$, dẫn đến phương trình cân bằng $\frac{1}{v^3} \left[ \lambda - v^2 \omega_0'(v) \right] = 0$, hay tương đương $\lambda = v^3 \omega_0'(v)$. Do lực cản chính có dạng $\omega_0(v) = a + bv + cv^2$, phương trình trở thành một hàm đơn điệu tăng theo vận tốc. Nghiệm của phương trình này xác định giá trị vận tốc ổn định $V_{on} = f(\lambda, a, b, c)$. Vận tốc này phụ thuộc trực tiếp vào hệ số nhân Lagrange $\lambda$ (yếu tố đại diện cho ràng buộc thời gian chạy tàu) và các hệ số lực cản đơn vị của đoàn tàu.

4. Luận án giải quyết bài toán hai đoàn tàu tránh nhau tại ga trung gian trên đường đơn như thế nào?

Đối với tuyến đường đơn có hai đoàn tàu tránh nhau tại ga C vào thời điểm $T_3$, hành trình từ ga A đến ga D được chia thành hai khoảng thời gian độc lập $T_1 = T_3 - T_1$ và $T_2 = T_8 - T_3$. Luận án đưa vào phiếm hàm tiêu chuẩn hai hệ số nhân Lagrange riêng biệt: $\lambda_1$ cho đoạn trước ga C ($s \le S_1$) và $\lambda_2$ cho đoạn sau ga C ($s > S_1$). Tại vị trí ranh giới $S_1$, hàm $p(s)$ có bước nhảy bậc từ $p^-$ sang $p^+$, kéo theo sự thay đổi mức vận tốc ổn định từ $V_{on1}$ sang $V_{on2}$, từ đó cho phép tối ưu hóa cục bộ từng phân đoạn mà vẫn đáp ứng thời điểm tránh tàu bắt buộc.

5. Phương pháp hiệu chuẩn các hệ số trong hàm lực cản chính của đoàn tàu được thực hiện như thế nào?

Thay vì sử dụng các thông số mặc định do nhà sản xuất cung cấp vốn có sai lệch so với điều kiện khai thác thực tế tại Việt Nam, luận án xây dựng chương trình tính toán nhận dạng lại các hệ số $a, b, c$ trong hàm $\omega_0(v) = a + bv + cv^2$. Quá trình hiệu chuẩn sử dụng thuật toán di truyền (GA) để tìm kiếm bộ tham số tối ưu, dựa trên sự khớp nối với dữ liệu đo thực tế về tốc độ, vị trí và công suất được ghi nhận từ thiết bị đo chuyên dụng gắn trên đầu máy chạy trên tuyến Yên Viên – Lào Cai.

Kết luận

Luận án đã giải quyết bài toán điều khiển tối ưu năng lượng chuyển động đoàn tàu bằng cách áp dụng Nguyên lý cực đại Pontryagin cho hệ phi tuyến không Autonom, mở rộng từ bài toán một đoàn tàu sang bài toán nhiều đoàn tàu có xét đến giao cắt tránh nhau tại ga trung gian và tương tác giãn cách đóng đường. Nghiên cứu đã làm rõ ảnh hưởng của sự phân bố chiều dài đoàn tàu khi vượt qua các đoạn dốc ngắn liên tiếp đối với tính bất biến của hàm Hamilton, đồng thời đề xuất giải pháp xác định điểm chuyển chế độ tối ưu trong 5 trạng thái vận hành cơ bản. Kết quả nghiên cứu kết hợp giữa mô hình hóa toán học, thuật toán di truyền hiệu chuẩn lực cản thực nghiệm và mô phỏng trên tuyến Yên Viên – Lào Cai tạo cơ sở khoa học để phát triển các hệ thống hỗ trợ lái tàu tiết kiệm nhiên liệu trong ngành đường sắt.