Tổng quan về luận án

Sự bùng nổ của cuộc Cách mạng Công nghiệp lần thứ tư (Công nghiệp 4.0) đã và đang tái cấu trúc toàn diện các nền tảng sản xuất và logistics toàn cầu. Trong hệ sinh thái nhà kho thông minh (Smart Warehouse) – trụ cột của chuỗi cung ứng hiện đại, các hệ thống tự động hóa như Unitload AS/RS (Automated Storage and Retrieval Systems) và mô hình GTP (Goods-to-People) đòi hỏi sự tham gia ngày càng sâu rộng của Robot tự hành (Autonomous Mobile Robot - AMR, Automated Guided Vehicle - AGV). Theo báo cáo "World Robotics" của Liên đoàn Robot Quốc tế (IFR), quy mô thị trường công nghệ robot toàn cầu năm 2021 đạt 79,5 tỷ USD, tăng lên 89,27 tỷ USD vào năm 2022 và được dự báo vượt ngưỡng 225,6 tỷ USD vào năm 2030 với tốc độ tăng trưởng kép hàng năm (CAGR) đạt 12,29%. Song song đó, phân tích từ Allied Market Research chỉ ra rằng riêng thị trường robot tự hành thương mại đã đạt giá trị 9.340 triệu USD vào năm 2018 và dự kiến cán mốc 39.585 triệu USD vào năm 2026 với CAGR ấn tượng 21,5%.

Theo tổng kết của McFarlane trong công trình "Intelligent logistics: Involving the customer", kho hàng thông minh không đơn thuần là sự số hóa mà là quá trình tích hợp sâu giữa lập kế hoạch, giám sát và vận hành thông minh theo thời gian thực. Wen Jinming trong nghiên cứu "Swarm robotics control and communications: Imminent challenges for next generation smart logistics" nhấn mạnh rằng thách thức lớn nhất của các hệ sinh thái logistics thông minh thế hệ mới chính là khả năng điều khiển và phối hợp chuyển động an toàn, chính xác của phi đội robot hoạt động trong không gian hẹp với mật độ vật cản cao.

Khoảng trống nghiên cứu (Research Gap) cốt lõi hiện nay nằm ở sự mất cân bằng giữa độ chính xác hình học và chi phí tính toán thời gian thực của các giải thuật tránh va chạm. Các tiếp cận lập kế hoạch toàn cục truyền thống (Global Path Planning) hoàn toàn bất lực trước các vật cản động xuất hiện ngẫu nhiên. Ngược lại, các phương pháp phản ứng cục bộ (Reactive Collision Avoidance) như Trường thế năng nhân tạo (Artificial Potential Field - APF) hay Dải đàn hồi truyền thống (Elastic Strips) của O. Khatib thường rơi vào điểm kỳ dị (singularities), cực tiểu cục bộ (local minima), hoặc gây ra hiện tượng bế tắc, dao động cơ học khi hướng chuyển động của robot trùng hoàn toàn với hướng của véc-tơ lực cản trực diện. Đồng thời, hầu hết các thuật toán hiện hành đều xử lý va chạm theo dạng nhị phân hình học thuần túy (va chạm hoặc không va chạm), thiếu hụt mô hình không gian xác suất rủi ro va chạm liên tục để tối ưu hóa quỹ đạo an toàn trong môi trường bất định.

Để giải quyết triệt để bài toán trên, luận án tiến sĩ công nghệ thông tin chuyên ngành Quản lý Hệ thống Thông tin (Mã số: 9480205.01QTD) của Nghiên cứu sinh Nông Minh Ngọc tại Viện Công nghệ Thông tin – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của PGS.TS Đỗ Năng Toàn và PGS.TS Vũ Việt Vũ, đã xác lập hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để phát hiện chính xác thời điểm và tọa độ tiếp xúc va chạm đa chiều giữa robot và vật cản phức tạp trong thời gian thực mà không làm bùng nổ độ phức tạp tính toán?
  • Câu hỏi nghiên cứu 2 (RQ2): Cần cải tiến cơ chế tổng hợp ngoại lực của giải thuật Elastic Strips như thế nào để triệt tiêu hoàn toàn điểm kỳ dị khi véc-tơ vận tốc đối đầu trực diện với điểm va chạm?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để thiết lập mô hình tính toán xác suất va chạm thể tích $P$ và trường phân vùng đồng mức phần bù xác suất $(1-P)$ nhằm tạo ra hành lang di chuyển mượt mà, thích nghi với độ bất định của cảm biến?
  • Giả thuyết nghiên cứu 1 (H1): Việc tích hợp cấu trúc cây phân lớp hệ bao BVH với hộp bao theo hướng OBB dựa trên Định lý Trục Tách biệt (SAT) 15 trục sẽ giảm bậc phức tạp tính toán kiểm tra va chạm xuống mức tuyến tính từng phần nhưng vẫn đảm bảo độ chính xác tiếp xúc tại thời điểm $T$.
  • Giả thuyết nghiên cứu 2 (H2): Cải tiến hàm ngoại lực thế năng bằng cách phân rã lực cản thành thành phần pháp tuyến và thành phần tiếp tuyến cưỡng bức sẽ giải phóng robot khỏi bế tắc chuyển động trong dải đàn hồi.
  • Giả thuyết nghiên cứu 3 (H3): Phân vùng đồng mức xác suất va chạm dựa trên tích phân hàm mật độ Gauss 2D cho phép robot tự hành thiết lập khoảng cách an toàn động (Dynamic Safety Margin), giảm thiểu đồng thời tỷ lệ cảnh báo sai (False Alarm) và tỷ lệ bỏ sót va chạm (Miss Detection).

Khung lý thuyết của nghiên cứu được xây dựng trên sự giao thoa của Lý thuyết Hình học Tính toán Phân cấp (Bounding Volume Hierarchies - BVH, Gottschalk et al.), Lý thuyết Trường Thế năng Nhân tạo (Khatib, 1985), Lý thuyết Dải Đàn hồi Động (Elastic Strips Theory, Khatib, 1999), và Lý thuyết Đánh giá Rủi ro Va chạm Xác suất (Probabilistic Collision Risk Assessment). Phạm vi nghiên cứu tập trung vào hệ thống robot di động tự hành bánh xe hoạt động trong môi trường 2D/3D của kho hàng thông minh, được kiểm chứng giải thuật thông qua mô phỏng toán học chuyên sâu trên MATLAB và lập trình C/C++.

Literature Review và Positioning

Lịch sử nghiên cứu điều hướng và tránh va chạm cho robot tự hành chứng kiến sự phân hóa rõ nét thành hai trường phái chính:

                                 HỆ THỐNG TRÁNH VA CHẠM ROBOT TỰ HÀNH
                                                   |
                 +---------------------------------+---------------------------------+
                 |                                                                   |
     Trường phái Dựa trên Kế hoạch                                      Trường phái Dựa trên Phản ứng
         (Planning-Based)                                                   (Reactive-Based)
                 |                                                                   |
     +-----------+-----------+                                       +---------------+---------------+
     |           |           |                                       |               |               |
Dijkstra (1959)  A* (1968)  RRT (2001)                         Certainty Grid    APF / VFF     Elastic Strips
   Offline/Tĩnh  Heuristic   Cây ngẫu nhiên                     Moravec (1988)  Khatib / Koren   Khatib (1999)
                                                                Nhiễu cảm biến  Local Minima    Kỳ dị trực diện
                                                                                                     |
                                                                                    +----------------+----------------+
                                                                                    |   ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN   |
                                                                                    +----------------+----------------+
                                                                                    | • BVH-OBB + 15 Trục cô lập SAT  |
                                                                                    | • Elastic Strips Cải tiến       |
                                                                                    | • Phân vùng đồng mức (1-P)      |
                                                                                    +---------------------------------+

Trường phái tiếp cận dựa trên lập kế hoạch (Planning-based Approach) khởi đầu từ các giải thuật tìm đường kinh điển như thuật toán Dijkstra (1959) và giải thuật A* (Hart, Nilsson & Raphael, 1968). Các phương pháp này tìm kiếm quỹ đạo tối ưu toàn cục dựa trên bản đồ topo cố định. Mở rộng hướng tiếp cận này, LaValle (2001) đề xuất thuật toán Cây ngẫu nhiên khám phá nhanh (Rapidly-exploring Random Tree - RRT) nhằm xử lý không gian cấu hình nhiều chiều, trong khi Chen và cộng sự phát triển hệ thống Định vị và Lập bản đồ đồng thời (SLAM) tự động cập nhật nút tọa độ. Gần đây, Shini và cộng sự (2018) ứng dụng Điều khiển Dự báo Mô hình (Model Predictive Control - MPC) để dẫn đường cho hệ thống đa robot. Tuy nhiên, điểm hạn chế chí tử của trường phái này là phụ thuộc tuyệt đối vào tri thức bản đồ tĩnh tiên nghiệm, tiêu tốn tài nguyên tính toán khổng lồ khi tái cấu hình bản đồ và hoàn toàn không thích ứng kịp với các tương tác động bất ngờ trong nhà kho.

Trường phái tiếp cận dựa trên phản ứng (Reactive-based Approach) vận hành theo thời gian thực dựa trên dữ liệu cảm biến cục bộ. Moravec và Elfes (1988) giới thiệu mô hình Lưới xác định (Certainty Grid) biểu diễn không gian ô lưới (Cell Value - CV) để xử lý nhiễu siêu âm. Borenstein và Koren (1989, 1991) phát triển phương pháp Trường lực ảo (Virtual Force Field - VFF) tính toán lực đẩy tỷ lệ nghịch với khoảng cách vật cản. O. Khatib (1985) đặt nền móng lý thuyết với Trường thế năng nhân tạo (APF), sau đó được Krogh và Thorpe (1986) mở rộng thành Trường thế năng tổng quát (GPF), Newman và Hogan (1987) bổ sung đặc tính hình học vật cản, và Brooks (1986), Arkin (1987) hiện thực hóa trên phần cứng cảm biến. Năm 1999, O. Khatib cùng Brock nâng tầm lý thuyết với phương pháp Dải đàn hồi (Elastic Strips), xem quỹ đạo chuyển động như một sợi dây đàn hồi chịu tác động của nội lực co giãn và ngoại lực đẩy từ vật cản.

Cuộc tranh luận học thuật sâu sắc nảy sinh giữa hai quan điểm đối nghịch:

  1. Quan điểm tất định hình học: Cho rằng kiểm tra va chạm phải dựa trên giao cắt khối đa diện chính xác để đảm bảo an toàn tuyệt đối, bất chấp chi phí tính toán cao.
  2. Quan điểm trường lực thế năng: Cho rằng cần đơn giản hóa hình học thành các điểm khối chất điểm để tối ưu hóa thời gian tính toán phản ứng mili-giây, chấp nhận rủi ro va chạm ở biên do mất thông tin góc cạnh.

Luận án của NCS. Nông Minh Ngọc định vị tại điểm giao thoa đột phá: Không hy sinh độ chính xác hình học cũng không đánh đổi tốc độ phản ứng thời gian thực. Bằng cách dung hợp Cây phân cấp hệ bao (BVH) chứa các Hộp bao định hướng (Oriented Bounding Box - OBB) với giải thuật Dải đàn hồi cải tiến, kết hợp cùng Lý thuyết Trường xác suất rủi ro va chạm, luận án đã định vị vượt trội so với các công trình quốc tế đương đại. Cụ thể:

  • So với nghiên cứu của Rostami và cộng sự (2019) về cải tiến giải thuật APF giải quyết cực tiểu cục bộ nhưng vẫn dựa trên chất điểm đơn giản, luận án giải quyết trọn vẹn mô hình vật lý 3D của cả robot lẫn vật cản bằng OBB 15 trục cô lập.
  • So với công trình của Sang và cộng sự (2021) kết hợp lai A* và APF cho phương tiện không người lái trên mặt nước (USV) vốn nặng về tính toán tĩnh, giải pháp của luận án phản ứng trơn tru với vật thể động và triệt tiêu hoàn toàn sự cố dừng đột ngột của dải đàn hồi Khatib khi tiếp cận trực diện.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mang lại ba đột phá lý thuyết mang tính nền tảng cho lĩnh vực Điều khiển học và Robot học:

Luận án mở rộng trực tiếp Lý thuyết Dải đàn hồi (Elastic Strips Theory) của O. Khatib (1999). Trong mô hình gốc của Khatib, khi robot tiếp cận vật cản có véc-tơ vị trí tương đối ngược hướng hoàn toàn với véc-tơ vận tốc mong muốn, tổng hợp ngoại lực $F_{ext}$ bị triệt tiêu thành phần vuông góc, dẫn đến việc robot bị bế tắc động học (Deadlock/Kinematic Singularity) hoặc dao động cơ học tần số cao không thể dập tắt. Luận án đã tái thiết lập công thức trường ngoại lực cản thông qua việc tích hợp hệ số trượt tiếp tuyến cưỡng bức, phá vỡ thế đối xứng lực và đảm bảo quỹ đạo đàn hồi luôn tự động uốn cong về phía có gradien thế năng thấp nhất.

Luận án hoàn thiện khung giải tích của Định lý Trục Tách biệt (Separating Axis Theorem - SAT) trong kiểm tra va chạm OBB đa diện lồi do Gottschalk, Lin và Manocha khởi xướng. Tác giả không chỉ dừng lại ở bài toán nhị phân "có/không va chạm" mà thiết lập hệ phương trình đóng để suy biến chính xác nhãn thời gian tiếp xúc đầu tiên $T$ và tọa độ tiếp điểm duy nhất $P(x, y, z)$ dựa trên 15 trục cô lập ứng viên:

$$\mathbf{V} \in \left{ \mathbf{A}_i, \mathbf{B}_j, \mathbf{A}_i \times \mathbf{B}_j \right} \quad (i, j \in {0, 1, 2})$$

Luận án kiến tạo Lý thuyết Trường Phân vùng Đồng mức Phần bù Xác suất $(1-P)$. Thay vì xem rủi ro là một hàm bước nhảy rời rạc, nghiên cứu mô hình hóa rủi ro va chạm thành một trường vô hướng liên tục trên mặt phẳng $E^2$, chứng minh một cách toán học rằng việc điều hướng robot dọc theo đường biên đồng mức xác suất tối ưu sẽ tối thiểu hóa năng lượng điều khiển cơ cấu chấp hành và duy trì cự ly an toàn động hoàn hảo.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột lý thuyết: Hình học phân cấp tính toán (BVH-OBB), Động lực học dải đàn hồi (Modified Elastic Strips), và Lý thuyết Đo lường Xác suất Rủi ro (Collision Probability Field).

                                 KHUNG PHÂN TÍCH TÍCH HỢP ĐỘC ĐÁO
+---------------------------------------------------------------------------------------------------+
| 1. TẦNG HÌNH HỌC PHÂN CẤP (BVH-OBB Filtering Layer)                                               |
|    • Lớp thô: Hộp bao trục AABB loại trừ 80%+ vật thể không có nguy cơ.                           |
|    • Lớp tinh: Hộp bao hướng OBB + SAT 15 trục chiếu xác định khoảng cách R, R0, R1.              |
|    • Giải tích tiếp điểm: Khi R = R0 + R1, truy xuất nghiệm thời gian T và tọa độ va chạm P.      |
+---------------------------------------------------------------------------------------------------+
+---------------------------------------------------------------------------------------------------+
| 2. TẦNG ĐỘNG LỰC DẢI ĐÀN HỒI CẢI TIẾN (Modified Elastic Strips Dynamics Layer)                    |
|    • Nội lực: F_int duy trì tính liên tục, độ trơn và độ co ngắn của quỹ đạo.                     |
|    • Ngoại lực cải tiến: F_ext phân rã lực đẩy + lực dịch góc khi vector hướng trùng điểm va chạm.|
|    • Cập nhật quỹ đạo thời gian thực: uốn cong dải đàn hồi, giải phóng bế tắc động học.           |
+---------------------------------------------------------------------------------------------------+
+---------------------------------------------------------------------------------------------------+
| 3. TẦNG ĐIỀU HƯỚNG ĐỒNG MỨC XÁC SUẤT (ISO-Probability Navigation Layer)                           |
|    • Mô hình sai số Gauss thể tích: Phân tích hiệp phương sai định vị và kích thước vật cản.      |
|    • Trường xác suất va chạm P và Trường phần bù an toàn (1-P).                                   |
|    • Lập quỹ đạo bám biên đồng mức (1-P): Hành lang di chuyển tối ưu, loại bỏ giật cục cơ học.   |
+---------------------------------------------------------------------------------------------------+

Điều kiện biên (Boundary Conditions) được xác lập rõ ràng: Không gian hoạt động là mặt phẳng 2D có giới hạn của nhà kho thông minh; các vật cản động và tĩnh là các khối đa diện lồi hoặc tập hợp các khối lồi; cảm biến đo đạc (LiDAR, siêu âm, camera) có sai số chuẩn Gauss trong dải đo hữu hạn; robot tự hành có giới hạn về vận tốc tuyến tính tối đa $v_{max}$, gia tốc góc $\omega_{max}$ và bán kính quay vòng tối thiểu.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt thế giới quan khoa học thực chứng (Positivism) kết hợp phương pháp luận diễn dịch - thực nghiệm định lượng. Nghiên cứu triển khai thiết kế đa tầng (Multi-level Analytical Design) bao gồm: Tầng lọc hình học thô (Coarse Geometric Filter), Tầng phân tích tiếp xúc chính xác (Exact Contact Analysis), Tầng biến dạng quỹ đạo đàn hồi (Elastic Trajectory Deformation), và Tầng kiểm soát rủi ro ngẫu nhiên (Stochastic Risk Regulation).

Quy mô thực nghiệm mô phỏng bao gồm 4 kịch bản động lực học tiêu chuẩn với sự tham gia của robot tự hành đối đầu với 1 vật cản, 2 vật cản song song, 3 vật cản hỗn hợp và hệ thống đa vật thể phân tán ngẫu nhiên, tái hiện các cấu trúc luồng di chuyển điển hình trong kho hàng tự động.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu kiểm tra và tránh va chạm được thiết kế hai pha khép kín:

                            QUY TRÌNH NGHIÊN CỨU & XỬ LÝ VA CHẠM HAI PHA
                                                 |
         +---------------------------------------+---------------------------------------+
         |                                                                               |
   PHA 1: BROAD-PHASE (Lọc nhanh)                                           PHA 2: NARROW-PHASE (Tính toán tinh)
         |                                                                               |
• Thiết lập hộp bao AABB theo trục toạ độ                               • Thiết lập hộp bao OBB theo hướng (C, A_i, a_i)
• Xác định toạ độ cực trị (xmin, ymin, zmin)                            • Chiếu 8 đỉnh lên trục cô lập d = C0 + t*V
• Kiểm tra điều kiện không giao nhau qua 4 bất đẳng thức                • Khảo sát 15 trục cô lập: 3 A_i, 3 B_j, 9 A_i x B_j
• Loại bỏ tức thì các cặp vật thể không nguy cơ                          • Tính bán kính chiếu R0, R1 và khoảng cách tâm R
                                                                        • Kiểm tra điều kiện phân tách: R > R0 + R1
                                                                                         |
                                                                        +----------------+----------------+
                                                                        |                                 |
                                                                  R > R0 + R1                       R <= R0 + R1
                                                                  (An toàn tuyệt đối)              (Xác nhận va chạm)
                                                                                                          |
                                                                                        • Lưu nhãn thời gian tiếp xúc T
                                                                                        • Giải tích tìm tiếp điểm P
                                                                                        • Truyền toạ độ P vào Elastic Strips
  1. Pha 1 (Broad-phase Collision Detection): Sử dụng hộp bao dạng trục (AABB) bao quanh robot và vật cản xác định bởi tâm $C_1, C_2$ và độ dài các cạnh $a_0, a_1, a_2$ và $b_0, b_1, b_2$. Tọa độ cực trị được tính toán:

$$x_{min1} = C_0 - \frac{a_0}{2}, \quad x_{max1} = C_0 + \frac{a_0}{2}$$

Kiểm tra nhanh điều kiện không giao nhau qua 4 hệ bất đẳng thức tọa độ. Pha này loại bỏ hơn 80% các cặp đối tượng không có nguy cơ va chạm trong không gian làm việc mà chỉ tiêu tốn chi phí tính toán cực tiểu $O(1)$.

  1. Pha 2 (Narrow-phase Exact Computation): Khi hộp bao AABB giao nhau, thuật toán kích hoạt cấu trúc Cây phân lớp hệ bao BVH chứa các hộp bao theo hướng (OBB). OBB được định nghĩa bởi tâm $C$, 3 véc-tơ trực chuẩn chỉ hướng $A_0, A_1, A_2$ và 3 bán trục độ dài $a_0, a_1, a_2 > 0$. Tám đỉnh của OBB được biểu diễn giải tích:

$$C + \sum_{i=0}^{2} s_i a_i A_i \quad \text{với } |s_i| = 1$$

Áp dụng Định lý Trục Tách biệt (SAT), tập hợp 15 trục cô lập ứng viên được thiết lập gồm: 3 trục hướng của hình thứ nhất ($A_i$), 3 trục hướng của hình thứ hai ($B_j$), và 9 trục tích có hướng ($A_i \times B_j$). Hình chiếu của 8 đỉnh lên trục cô lập $d = C_0 + tV$ tạo ra bán kính chiếu cực đại:

$$R_0 = a_0 |A_0 \cdot V| + a_1 |A_1 \cdot V| + a_2 |A_2 \cdot V|$$

$$R_1 = b_0 |B_0 \cdot V| + b_1 |B_1 \cdot V| + b_2 |B_2 \cdot V|$$

Khoảng cách hình chiếu giữa hai tâm:

$$R = |(C_1 - C_0) \cdot V|$$

Hai hộp bao OBB hoàn toàn không giao nhau khi và chỉ khi tồn tại ít nhất một trục cô lập thỏa mãn $R > R_0 + R_1$. Ngược lại, nếu trên toàn bộ 15 trục đều thỏa mãn $R \le R_0 + R_1$, va chạm chắc chắn xảy ra.

Điểm tiếp xúc va chạm $P$ tại thời điểm tới hạn $T$ (khi $R = R_0 + R_1$) được giải mã chính xác thông qua hệ phương trình:

$$P = C_0 + T v_0 + \sum_{i=0}^{2} x_i A_i = C_1 + T v_1 + \sum_{j=0}^{2} y_j B_j$$

Với $V = A_i$, nghiệm tọa độ tiếp xúc giải tích trên trục $A_i$ đạt giá trị:

$$x_i = \varepsilon \cdot \left[ a_i + \sum_{j=0}^{2} b_j |c_{ij}| \right] + \sum_{j=0}^{2} y_j c_{ij} \quad (\varepsilon = \pm 1, c_{ij} = A_i \cdot B_j)$$

Data và phân tích

Dữ liệu hình học và động học được chuẩn hóa với các thông số vật lý chính xác:

  • Bán kính an toàn cơ sở của robot: $R_{robot} = 0.35\text{ m}$ đến $0.5\text{ m}$.
  • Kích thước chướng ngại vật mô phỏng: từ $0.2\text{ m} \times 0.2\text{ m}$ đến $1.2\text{ m} \times 0.8\text{ m}$.
  • Vận tốc di chuyển của robot: $v_0 = 0.5\text{ m/s}$ đến $1.2\text{ m/s}$.
  • Vận tốc vật cản động: $v_1 = 0.2\text{ m/s}$ đến $0.6\text{ m/s}$.
  • Chu kỳ trích xuất tín hiệu cảm biến thời gian thực: $\Delta t = 50\text{ ms}$.

Môi trường tính toán và kiểm chứng giải thuật được triển khai trên nền tảng MATLAB R2022b kết hợp module tối ưu hóa C++ MEX functions để kiểm tra độ trễ thuật toán.

Phát hiện đột phá và implications

Những phát hiện then chốt

Nghiên cứu mang lại 4 phát hiện đột phá có giá trị khoa học và ứng dụng thực tiễn cao:

+---------------------------------------------------------------------------------------------------+
|                                BỐN PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                                  |
+---------------------------------------------------------------------------------------------------+
| 1. Triệt tiêu hoàn toàn điểm kỳ dị bế tắc chuyển động của Elastic Strips                          |
|    • Phương pháp Khatib (1999): Triệt tiêu lực tiếp tuyến khi đối đầu trực diện -> Robot dừng kẹt.|
|    • Cải tiến của Luận án: Bổ sung lực trượt góc tiếp tuyến -> Quỹ đạo uốn cong mượt 100%.       |
+---------------------------------------------------------------------------------------------------+
| 2. Cắt giảm đột biến độ trễ tính toán phát hiện va chạm qua Bảng tra 15 trục OBB                 |
|    • Giải tích hóa phép chiếu SAT giúp thời gian xử lý < 2.5 ms (chu kỳ điều khiển 50 ms).        |
|    • Tích hợp BVH giúp giảm hơn 80% phép thử OBB không cần thiết từ tầng Broad-phase.             |
+---------------------------------------------------------------------------------------------------+
| 3. Xây dựng Hành lang an toàn mượt mà thông qua Phân vùng đồng mức (1-P)                          |
|    • Thay thế biên giới nhị phân cứng bằng trường vô hướng xác suất liên tục.                     |
|    • Loại bỏ hoàn toàn hiện tượng rung lắc/giật cục (Chattering) cơ học tại cơ cấu chấp hành bánh.|
+---------------------------------------------------------------------------------------------------+
| 4. Cân bằng tối ưu giữa Cảnh báo sai (False Alarm) và Bỏ sót va chạm (Miss Detection)             |
|    • Phân tích độ nhạy xác suất thể tích giúp tối ưu hóa ngưỡng kích hoạt né tránh.               |
|    • Đảm bảo xác suất an toàn vận hành P_safe > 99.85% trong mọi kịch bản mô phỏng đa vật thể.   |
+---------------------------------------------------------------------------------------------------+
  • Phát hiện 1: Khắc phục triệt để hiện tượng dừng kẹt do kỳ dị đối xứng trong giải thuật Elastic Strips. Trong kịch bản vật cản di chuyển trực diện ngược chiều với robot trên cùng trục tọa độ, giải thuật Elastic Strips nguyên bản của Khatib làm triệt tiêu thành phần ngoại lực tiếp tuyến, khiến robot đứng yên hoặc rung lắc mất kiểm soát. Thuật toán Elastic Strips cải tiến của luận án với thành phần lực đẩy phân rã góc đã tạo ra mô-men né tránh tự nhiên, bẻ cong dải đàn hồi sang vùng không gian trống với độ lệch quỹ đạo tối thiểu và đưa robot về đích an toàn.
  • Phát hiện 2: Tối ưu hóa thời gian tính toán kiểm tra va chạm OBB nhờ bảng tra trước 15 trục. Việc giải tích hóa toàn bộ các biểu thức tính $R, R_0, R_1$ và tọa độ tiếp điểm $P$ giúp thời gian tính toán cho mỗi chu kỳ kiểm tra va chạm OBB giảm xuống dưới $2,5\text{ ms}$, đáp ứng hoàn hảo yêu cầu khắt khe của hệ thống thời gian thực (Real-time Systems) với tần số lấy mẫu $20\text{ Hz}$ ($\Delta t = 50\text{ ms}$).
  • Phát hiện 3: Trường phân vùng đồng mức phần bù xác suất $(1-P)$ triệt tiêu hiện tượng giật cục cơ học (Chattering). Khi áp dụng ranh giới an toàn hình học cứng, robot liên tục đổi hướng gấp khúc khi chạm biên an toàn. Ngược lại, quỹ đạo bám theo đường đẳng mức xác suất $(1-P)$ tạo ra độ cong chuyển động cấp $C^2$ trơn tru, giúp bảo vệ động cơ chấp hành và duy trì vận tốc ổn định.
  • Phát hiện 4: Sự cân bằng tối ưu giữa Báo động sai và Bỏ sót va chạm trong không gian xác suất thể tích. Phân tích độ nhạy xác suất va chạm chỉ ra rằng ngưỡng xác suất tới hạn $P_{threshold} = 0.05$ đến $0.08$ là điểm cân bằng Pareto tối ưu, loại bỏ hoàn toàn hiện tượng phanh gấp không cần thiết khi vật cản chỉ lướt qua vùng biên ngoài.

Implications đa chiều

  • Ý nghĩa lý thuyết: Cung cấp phương pháp giải tích tường minh cho bài toán tìm tiếp điểm va chạm trong không gian 3 chiều của các vật thể OBB lồi; đồng thời mở rộng lý thuyết điều khiển dải đàn hồi sang miền xác suất ngẫu nhiên.
  • Đổi mới phương pháp luận: Thiết lập quy trình chuẩn 3 tầng (Lọc AABB $\rightarrow$ Tính tiếp điểm OBB/SAT $\rightarrow$ Uốn dải đàn hồi bám đường đồng mức xác suất), có khả năng chuyển giao áp dụng cho nhiều lớp phương tiện tự hành khác nhau.
  • Ứng dụng thực tiễn: Triển khai trực tiếp vào phần mềm điều khiển trung tâm của robot vận chuyển hàng hóa tự động AGV/AMR trong các nhà máy thông minh, kho thương mại điện tử quy mô lớn, giúp nâng cao mật độ hoạt động của robot trên một đơn vị diện tích sàn mà không lo ngại xung đột luồng giao thông.
  • Khuyến nghị chính sách: Cung cấp cơ sở khoa học để các cơ quan quản lý tiêu chuẩn hóa quy chuẩn an toàn kỹ thuật (Safety Standards) cho phương tiện tự hành công nghiệp hoạt động trong môi trường có sự tương tác gần với con người (Human-Robot Interaction - HRI).

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 hạn chế khoa học cần tiếp tục hoàn thiện:

  1. Xấp xỉ hình học cho vật thể phi lồi phức tạp: Kỹ thuật OBB đơn lẻ vẫn tạo ra khoảng không gian trống dư thừa đối với các vật cản có hình dạng lõm sâu (ví dụ: kệ hàng dạng chữ U, xe nâng có càng vươn). Cần sử dụng phân rã đa diện lồi (Convex Decomposition) hoặc k-DOPs (k-Discrete Oriented Polytopes) để khớp sát biên dạng hơn.
  2. Giả định phân phối sai số Gauss: Mô hình xác suất va chạm thể tích trong luận án hiện giả định sai số vị trí tuân theo phân phối chuẩn Gauss. Trong môi trường công nghiệp thực tế với hiện tượng phản xạ đa đường của sóng siêu âm hoặc bề mặt kim loại làm lóa cảm biến LiDAR, nhiễu có thể mang đặc tính phi Gauss phức tạp.
  3. Ràng buộc động học phi toàn chỉnh cấp cao: Nghiên cứu tập trung giải quyết bài toán quỹ đạo hình học và lực phản ứng đàn hồi, chưa tích hợp sâu các ràng buộc động lực học phi toàn chỉnh (Non-holonomic constraints) cấp cao như giới hạn gia tốc giật (Jerk limitation) và hiện tượng trượt bánh xe trên bề mặt sàn trơn trượt.
  4. Quy mô kiểm chứng phần cứng: Các kết quả thực nghiệm mới dừng lại ở mức mô phỏng toán học chi tiết và kiểm chứng thuật toán trên dữ liệu số hóa; việc thử nghiệm thực địa trên phi đội hàng chục robot vật lý hoạt động liên tục 24/7 cần nguồn lực triển khai lớn hơn.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) bao gồm:

  • Tích hợp kỹ thuật biểu diễn hình học hàm khoảng cách có dấu (Signed Distance Fields - SDF) kết hợp k-DOPs để xử lý vật thể phi lồi thời gian thực.
  • Ứng dụng Bộ lọc Kalman không mùi (Unscented Kalman Filter - UKF) hoặc Lọc hạt (Particle Filter) để ước lượng ma trận hiệp phương sai sai số phi tuyến trong mô hình xác suất.
  • Kết hợp học tăng cường sâu (Deep Reinforcement Learning) để tự động tinh chỉnh hệ số đàn hồi nội lực $k_{int}$ và ngoại lực $k_{ext}$ theo mật độ kho bãi.
  • Mở rộng thuật toán cho không gian 3 chiều hoàn chỉnh ứng dụng cho thiết bị bay không người lái (UAV) và phương tiện lặn ngầm tự hành (AUV).

Tác động và ảnh hưởng

  • Tác động học thuật: Luận án mở ra hướng nghiên cứu mới trong việc kết hợp hình học phân cấp tính toán với trường thế năng dải đàn hồi và trường xác suất rủi ro. Các công bố khoa học từ luận án (công trình CT1, CT2) đóng góp tài liệu tham khảo giá trị cao cho cộng đồng nghiên cứu Robot học và Tự động hóa trong và ngoài nước.
  • Chuyển đổi công nghiệp logistics: Việc ứng dụng thuật toán giúp các doanh nghiệp logistics tối ưu hóa hiệu suất vận hành của hệ thống kho thông minh, nâng cao tốc độ luân chuyển hàng hóa thêm 15–20% nhờ giảm thời gian chờ và triệt tiêu tình trạng tắc nghẽn cục bộ của robot tại các lối đi hẹp.
  • Nâng cao an toàn lao động và xã hội: Trong bối cảnh mô hình cộng tác Người - Robot (Cobot/HRI) ngày càng phổ biến, thuật toán tránh va chạm dựa trên phân vùng đồng mức xác suất cung cấp "lá chắn vô hình" bảo vệ tuyệt đối cho công nhân vận hành trong kho, giảm thiểu rủi ro tai nạn lao động công nghiệp.
  • Tầm vóc quốc tế: Nghiên cứu giải quyết trực tiếp bài toán cốt lõi của công nghệ Robot tự hành theo các tiêu chuẩn quốc tế của IEEE Robotics and Automation Society, tạo tiền đề để công nghệ tự hành nội địa của Việt Nam hội nhập và cạnh tranh sòng phẳng với các giải pháp toàn cầu.

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

+---------------------------------------------------------------------------------------------------+
|                               HỆ THỐNG ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                              |
+---------------------------------------------------------------------------------------------------+
| 1. Nghiên cứu sinh & Nhà khoa học chuyên ngành:                                                   |
|    • Tiếp cận khung toán học tường minh về SAT 15 trục và giải tích tọa độ tiếp điểm P.           |
|    • Khai thác hướng nghiên cứu giao thoa giữa BVH, Dải đàn hồi và Trường xác suất (1-P).         |
+---------------------------------------------------------------------------------------------------+
| 2. Kỹ sư R&D & Doanh nghiệp chế tạo Robot tự hành (AMR/AGV):                                      |
|    • Nhận chuyển giao thuật toán Elastic Strips cải tiến không điểm kỳ dị, dễ dàng nạp vào chip.  |
|    • Cắt giảm chi phí phần cứng xử lý nhờ tối ưu hóa thời gian tính toán kiểm tra va chạm < 2.5ms.|
+---------------------------------------------------------------------------------------------------+
| 3. Quản lý kho bãi & Chuyên gia Tự động hóa Chuỗi cung ứng:                                       |
|    • Gia tăng mật độ luồng xe tự hành trên cùng diện tích mặt bằng kho thông minh.                |
|    • Triệt tiêu hoàn toàn sự cố dừng kẹt dây chuyền vận chuyển hàng hóa tự động AS/RS và GTP.    |
+---------------------------------------------------------------------------------------------------+
| 4. Cơ quan hoạch định & Ban Quản lý Tiêu chuẩn Công nghệ:                                         |
|    • Sở hữu bộ công cụ toán học định lượng rủi ro va chạm để xây dựng khung an toàn kỹ thuật số.  |
+---------------------------------------------------------------------------------------------------+

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất của luận án là việc mở rộng Lý thuyết Dải đàn hồi (Elastic Strips Theory) của O. Khatib (1999) kết hợp với việc hình thức hóa nghiệm giải tích của Định lý Trục Tách biệt (SAT) trên cấu trúc OBB 15 trục cô lập. Luận án đã giải quyết trọn vẹn điểm kỳ dị đối xứng – một khiếm khuyết nội tại của mô hình Khatib khiến dải đàn hồi bị triệt tiêu lực bẻ góc khi tiếp cận trực diện vật cản. Bằng cách thiết lập cơ chế phân rã lực cản và áp dụng trường phân vùng đồng mức phần bù xác suất $(1-P)$, luận án đã biến dải đàn hồi từ một mô hình thuần cơ học tất định thành một hệ thống điều khiển thích nghi ngẫu nhiên có độ mượt động học vượt trội.

2. Đổi mới phương pháp luận của luận án so với các công trình quốc tế tiêu biểu?

So với nghiên cứu của Rostami và cộng sự (2019) (cải tiến giải thuật Trường thế năng nhân tạo APF nhưng xem robot và vật cản là các chất điểm hình cầu đơn giản), luận án duy trì mô hình hình học 3D đầy đủ thông qua hộp bao định hướng OBB, bảo toàn nguyên vẹn thông tin góc cạnh và hướng của vật thể mà không làm tăng bậc phức tạp tính toán.

So với công trình của Sang và cộng sự (2021) (kết hợp lai giữa thuật toán tĩnh A* và APF cho phương tiện hàng hải USV), phương pháp của luận án vượt trội ở khả năng thích ứng thời gian thực với vật cản động ngẫu nhiên: Sự kết hợp giữa Cây phân lớp hệ bao BVH ở tầng lọc thô và giải tích tiếp điểm va chạm ở tầng lọc tinh cho phép trích xuất chính xác thời điểm va chạm $T$ và tọa độ tiếp xúc $P$ với độ trễ tính toán dưới $2,5\text{ ms}$, trong khi giải thuật của Sang et al. tiêu tốn tài nguyên đáng kể để tái lập kế hoạch toàn cục khi có biến động môi trường.

3. Phát hiện gây bất ngờ nhất trong quá trình tính toán và mô phỏng là gì?

Phát hiện bất ngờ nhất là: Việc gia tăng kích thước biên an toàn hình học thuần túy (Geometry Safety Margin) không đồng nghĩa với việc gia tăng độ an toàn thực tế của robot, mà ngược lại còn làm gia tăng nguy cơ bế tắc động học và mất ổn định cơ cấu chấp hành. Khi nới rộng hộp bao hình học một cách cơ học, không gian tự do khả dụng bị thu hẹp đột ngột, khiến giải thuật dải đàn hồi tạo ra các biến dạng quá mức (over-deformation) tại các khúc cua hẹp, dẫn đến hiện tượng dao động vận tốc cực mạnh. Ngược lại, khi chuyển sang điều hướng bám theo đường biên của trường phân vùng đồng mức xác suất $(1-P)$, robot tự động lựa chọn quỹ đạo uốn lượn mượt mà với năng lượng tiêu hao thấp nhất, loại bỏ hoàn toàn các đỉnh xung lực đột ngột trên động cơ chấp hành.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp đầy đủ hệ thống công thức giải tích đóng cho toàn bộ 15 trường hợp trục cô lập SAT (Bảng các giá trị $R, R_0, R_1$), hệ phương trình xác định tọa độ tiếp điểm $P$ theo tham số $T, v_0, v_1$, thuật toán phân vùng đồng mức xác suất dạng mã giả (pseudocode), cùng hệ thống bảng tham số hình học và động lực học chi tiết (Bảng tham số tính toán xác suất va chạm, dữ liệu tính toán xác suất phần bù va chạm). Các nhà nghiên cứu hoàn toàn có thể lập trình tái lập 100% kết quả trên bất kỳ môi trường mô phỏng nào (như MATLAB, ROS/Gazebo, C/C++).

5. Lộ trình nghiên cứu 10 năm (10-Year Research Agenda) được định hình như thế nào?

Lộ trình 10 năm hướng tới sự phát triển toàn diện của Hệ thống Tự hành Thông minh:

  • Giai đoạn 2024–2026: Tích hợp bộ giải thuật BVH-OBB và Elastic Strips cải tiến vào hệ điều hành Robot chuyên dụng (ROS 2), thử nghiệm trên các đội xe AGV/AMR vật lý trong môi trường kho công nghiệp thực tế; mở rộng sang mô hình phi lồi k-DOPs.
  • Giai đoạn 2027–2030: Kết hợp Trí tuệ nhân tạo tăng cường sâu (Deep Reinforcement Learning) và Mạng nơ-ron đồ thị (Graph Neural Networks - GNN) để tối ưu hóa phi đội đa robot phân tán (Swarm Robotics) trong không gian kho 3D đa tầng.
  • Giai đoạn 2031–2034: Mở rộng lý thuyết trường rủi ro xác suất đồng mức cho các phương tiện tự hành không gian 3 chiều phức tạp hoạt động ở tốc độ cao như Thiết bị bay giao hàng đô thị (UAM/UAV) và Robot cứu hộ cứu nạn trong môi trường khắc nghiệt.

Kết luận

Luận án tiến sĩ của NCS. Nông Minh Ngọc đại diện cho một công trình nghiên cứu hàn lâm mẫu mực, có cấu trúc chặt chẽ và tính ứng dụng thực tiễn vượt trội, tạo ra bước đột phá quan trọng trong kỹ thuật phát hiện và tránh va chạm cho robot tự hành.

Năm đóng góp cụ thể, mang tính cốt lõi của công trình bao gồm:

  1. Xây dựng thành công thuật toán phát hiện va chạm chính xác hai mức dựa trên Cây phân lớp hệ bao BVH và Hộp bao theo hướng OBB, giải tích hóa thành công hệ phương trình tìm nhãn thời gian $T$ và tọa độ tiếp điểm $P$ trên 15 trục cô lập SAT với chi phí tính toán thời gian thực tối ưu.
  2. Đề xuất giải thuật Dải đàn hồi (Elastic Strips) cải tiến, loại bỏ hoàn toàn điểm kỳ dị bế tắc chuyển động khi robot đối đầu trực diện với vật cản thông qua cơ chế phân rã ngoại lực thế năng có thành phần dịch chuyển góc cưỡng bức.
  3. Thiết lập mô hình tính toán xác suất va chạm thể tích dựa trên tích phân hàm mật độ Gauss, làm cơ sở khoa học định lượng hóa rủi ro va chạm trong môi trường tồn tại sai số cảm biến.
  4. Phát triển thuật toán phân vùng đồng mức phần bù xác suất $(1-P)$, tạo ra hành lang điều hướng tối ưu, triệt tiêu hiện tượng giật cục cơ học và đảm bảo quỹ đạo chuyển động đạt độ trơn hình học cao.
  5. Kiểm chứng hoàn chỉnh giải thuật qua hệ thống kịch bản mô phỏng đa dạng từ đơn vật thể đến đa vật thể chuyển động phức tạp, chứng minh tính vượt trội về độ ổn định, tính an toàn và khả năng đáp ứng thời gian thực so với các giải pháp kinh điển.

Sự chuyển dịch mô thức (Paradigm Shift) mà luận án mang lại chính là việc thay thế tư duy né tránh va chạm hình học nhị phân cứng nhắc bằng Tư duy Điều khiển Dải đàn hồi Tối ưu bám Trường Rủi ro Xác suất Liên tục. Công trình mở ra ít nhất ba nhánh nghiên cứu mới: Điều khiển phi đội robot bầy đàn thích nghi xác suất; Dung hợp cảm biến thời gian thực trên cấu trúc hình học đa diện lồi; và Thiết lập chuẩn an toàn động học cho các hệ thống tương tác Người - Robot trong kỷ nguyên Công nghiệp 4.0 và 5.0.