Tổng quan luận án

Nghiên cứu về cơ cấu robot công nghiệp đóng vai trò quan trọng trong việc tự động hóa quá trình sản xuất, nâng cao năng suất, hạ giá thành và cải thiện điều kiện làm việc trong các môi trường độc hại. Robot song song là hệ thống cơ học có cấu trúc mạch vòng kín, trong đó khâu tác động cuối (bàn máy động) được liên kết với bệ máy cố định thông qua nhiều chuỗi động học (các chân). Nhờ cấu trúc liên kết song song, tải trọng được chia sẻ đồng đều giữa các chân, giúp robot đạt độ cứng vững cao, tỷ lệ tải trọng trên khối lượng bản thân lớn, khả năng phát động gia tốc cao và triệt tiêu sai số tích lũy so với robot dạng chuỗi nối tiếp. Tuy nhiên, trở ngại lớn nhất trong việc ứng dụng robot song song là sự xuất hiện của các điểm kỳ dị động học và động lực học trong không gian làm việc. Tại các cấu hình kỳ dị, ma trận Jacobi suy biến làm cơ hệ mất hoặc tăng bất thường số bậc tự do, dẫn đến việc mất khả năng điều khiển hoặc làm phát sinh lực và mô men cục bộ vô cùng lớn, gây nguy cơ phá hủy kết cấu cơ khí.

Bên cạnh đó, hầu hết các công trình nghiên cứu hiện hữu chỉ tập trung khảo sát mô hình cơ học thuần túy của các khâu khớp mà bỏ qua động lực học của hệ dẫn động (bao gồm động cơ điện một chiều và hộp giảm tốc). Trong thực tế hoạt động tốc độ cao, các đặc tính điện - cơ, quán tính của rô-to và mô men phi tuyến từ hệ dẫn động gây ảnh hưởng đáng kể đến độ chính xác bám quỹ đạo. Ngoài ra, việc thiếu hụt cảm biến đo biến khớp bị động cùng với các yếu tố bất định tham số và nhiễu ngoài đặt ra yêu cầu phải có giải pháp ước lượng động học và bộ điều khiển phi tuyến bền vững.

Luận án tiến sĩ cơ học của tác giả Vũ Đức Vương, với đề tài "Phân tích kỳ dị và điều khiển trượt mờ robot song song phẳng có kể đến hệ dẫn động", được thực hiện tại Trường Đại học Bách khoa Hà Nội dưới sự hướng dẫn khoa học của PGS. Nguyễn Quang Hoàng, tập trung giải quyết đồng bộ các vấn đề mô hình hóa, nhận diện kỳ dị, khắc phục kỳ dị và thiết kế điều khiển bền vững cho robot song song phẳng.

Mục tiêu tổng quát của luận án là góp phần hoàn thiện phương pháp khảo sát động học, xây dựng mô hình động lực học có kể đến hệ thống dẫn động và phát triển các giải pháp điều khiển tiên tiến nhằm nâng cao chất lượng hoạt động của robot song song phẳng. Luận án đặt ra 05 nhiệm vụ nghiên cứu cụ thể được đánh số như sau:

  1. Đưa ra cơ sở lý thuyết xây dựng mô hình động lực học cho robot song song có kể đến ảnh hưởng của động cơ dẫn động.
  2. Xác định kỳ dị và phát triển giải pháp vượt kỳ dị trong giải quyết bài toán động học, động lực học robot song song.
  3. Xây dựng thuật toán điều khiển trực tiếp robot song song trong không gian thao tác.
  4. Xây dựng bộ ước lượng động học xác định tham số thay thế cho thông tin thu được từ thiết bị đo thực trên các khớp bị động của robot song song.
  5. Xử lý các yếu tố bất định của môi trường và sai lệch tham số ảnh hưởng đến hoạt động của robot song song bằng bộ điều khiển trượt mờ.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Các cơ cấu robot song song được dẫn động bởi động cơ điện một chiều kết hợp hộp giảm tốc, tập trung vào các bài toán phân tích động học, mô hình hóa động lực học và các thuật toán điều khiển phi tuyến.
  • Phạm vi nghiên cứu: Khảo sát lý thuyết và tiến hành thực nghiệm số trên đối tượng đại diện cụ thể là robot song song phẳng 3 bậc tự do dạng 3RRR (mỗi chân gồm 3 khớp quay Revolute – Revolute – Revolute). Nghiên cứu giới hạn trong không gian chuyển động phẳng 2D đối với cơ cấu 3RRR đối xứng được dẫn động bởi 3 động cơ điện một chiều qua hộp giảm tốc.

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

Tổng quan lịch sử phát triển cho thấy cấu trúc robot song song bắt đầu từ nguyên mẫu bệ đỡ rạp chiếu phim do Gwinnet phát minh năm 1928 [1], tiếp đó là cơ cấu kiểm tra lốp xe vòng kín 6 bậc tự do của Gough năm 1947 (hoàn thiện năm 1955) [2], và hệ thống mô phỏng buồng lái máy bay do Stewart đề xuất năm 1965 [3]. Về sau, các biến thể nổi bật như robot Delta không gian của Clavel [4][5][6] và các dạng robot song song phẳng 3RPR, 3RRR được ứng dụng rộng rãi trong công nghiệp đóng gói, định vị gia công (robot CMW380), đào tạo phi hành gia (hệ thống Rendezvous của NASA [12]) và phẫu thuật y tế chính xác 6 bậc tự do của M. Wapler [11].

So sánh đặc tính kỹ thuật cơ bản giữa robot chuỗi và robot song song được tổng hợp trong công trình nghiên cứu:

Tiêu chí so sánh Robot chuỗi (nối tiếp) Robot song song
Độ chính xác khâu chấp hành cuối Thấp hơn (do sai số tích lũy) Cao hơn (sai số được bù trung bình giữa các chân)
Độ chính xác khi lặp lại Trung bình Rất cao
Độ cứng vững cơ học Thấp (chịu uốn lớn) Cao (tải trọng chuyển thành lực kéo/nén)
Tỷ lệ khối lượng tải / robot Thấp Cao (tiết kiệm năng lượng khi mang tải)
Gia tốc khâu thao tác cuối Trung bình Lớn (quán tính chuyển động nhỏ)
Không gian làm việc Rộng, linh hoạt tránh vật cản Hẹp, giới hạn bởi các mạch vòng kín
Tính linh hoạt thao tác Cao Kém hơn
Chi phí thời gian tính toán Thấp Cao và phức tạp

Về phân tích kỳ dị, đối với robot chuỗi, các phương pháp bình phương tối thiểu có cản DLS (Damped Least Squares) và phân rã giá trị kỳ dị SVD (Singular Value Decomposition) thường được áp dụng [20]–[22]. Tuy nhiên, các kỹ thuật này không áp dụng trực tiếp được cho robot song song do đặc tính vòng kín. Merlet [23] đã dùng hình học Grassmann để phân tích kỳ dị bệ Stewart; Park và Kim [24] dùng hình học Riemann phân loại khả năng thao tác; Ma [16] phân tích cấu hình kỳ dị dựa trên kiến trúc hệ thống; Gosselin và Angeles [15], [19] phát triển phương pháp giải tích ma trận Jacobi phân loại kỳ dị thành ba nhóm chính: kỳ dị động học thuận, kỳ dị động học ngược và kỳ dị hỗn hợp. Về giải pháp vượt kỳ dị động lực học, các nghiên cứu của Oetomo và Ang [74], Ider [75], Jazar [76] và Schuler et al. [77] đã đề xuất việc hiệu chỉnh phương trình chuyển động bằng đạo hàm bậc cao.

Về mô hình hóa động lực học, các phương pháp truyền thống như phương trình Lagrange dạng nhân tử, phương trình Newton – Euler, nguyên lý công ảo Jourdain và phương trình Kane [27]–[31] đã được nhiều tác giả trong nước và quốc tế sử dụng. Tại Việt Nam, các nhóm nghiên cứu tại Viện Cơ học, Trường Đại học Bách khoa Hà Nội và Trường Đại học Bách khoa TP.HCM đã công bố nhiều nghiên cứu về động lực học hệ nhiều vật mạch vòng [22], [24], [27], [34]–[37], phương pháp số giải động lực học ngược [31], [49], biến đổi phương trình vi phân đại số (DAE) sang phương trình vi phân thường (ODE) [50], robot song song dùng trong máy CNC 5 trục [57] và phương pháp xoắn vít Plücker xác định kỳ dị [56], [58].

Về chiến lược điều khiển, các bộ điều khiển kinh điển như PD bù trọng lực, PID mô men tính toán [5], [48]–[50] dần bộc lộ hạn chế khi hệ thống có tính phi tuyến mạnh và chịu bất định. Các hướng tiếp cận mới kết hợp điều khiển mờ (Fuzzy Logic), mạng nơ-ron và giải thuật di truyền (GA) được triển khai trên các mô hình như robot Isoglide 3 bậc tự do [62], [63], [64]. Các nghiên cứu điều khiển trong không gian thao tác [65]–[71] mở ra khả năng tác động trực tiếp lên quỹ đạo bàn thao tác nhưng thường đơn giản hóa mô hình.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết:

  • Đa phần các công bố trong và ngoài nước chỉ thiết lập mô hình cơ học của các khâu khớp mà bỏ qua phương trình động lực học của hệ dẫn động (động cơ điện một chiều và hộp giảm tốc), dẫn đến sự sai lệch khi thiết kế điều khiển bám tốc độ cao.
  • Thiếu các giải pháp số hoàn chỉnh cho bài toán vượt qua các vùng kỳ dị động học và động lực học nhằm duy trì tính liên tục của quỹ đạo.
  • Việc thu nhận biến khớp bị động phục vụ hồi tiếp trạng thái gặp khó khăn về không gian lắp đặt và chi phí cảm biến, đòi hỏi một thuật toán ước lượng động học tin cậy không bị trôi dạt do sai số tích phân.
  • Chưa có cấu trúc điều khiển trượt mờ thích nghi tối ưu bằng GA vừa xử lý triệt để hiện tượng chattering vừa kháng được độ bất định tham số lớn trong không gian thao tác.

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

Luận án kết hợp chặt chẽ giữa cơ sở lý thuyết giải tích cơ học, giải thuật tối ưu hóa tiến hóa và phương pháp mô phỏng số trên máy tính.

Khung lý thuyết và khái niệm sử dụng:

  • Cơ học giải tích: Sử dụng phương trình Lagrange loại hai dạng nhân tử kết hợp biến đổi tọa độ suy rộng tối thiểu để thiết lập mô hình toán học cho hệ nhiều vật có liên kết mạch vòng kín. Ổn định hóa nghiệm phương trình vi phân đại số bằng phương pháp Baumgarte và phương pháp chiếu hiệu chỉnh nhằm triệt tiêu sai số ràng buộc.
  • Mô hình hệ dẫn động điện cơ: Khảo sát hệ thống động cơ điện một chiều kích từ vĩnh cửu kết nối qua hộp giảm tốc tỉ số truyền $r_g$. Phương trình mạch phần ứng và phương trình cân bằng mô men rô-to: $$u_i = R_{a,i} i_{a,i} + L_{a,i} \frac{d i_{a,i}}{dt} + K_{e,i} \omega_{m,i}$$ $$M_{m,i} = K_{m,i} i_{a,i} = I_{m,i} \dot{\omega}{m,i} + b{m,i} \omega_{m,i} + \frac{\tau_i}{r_g}$$ trong đó $u_i$ là điện áp điều khiển, $R_a, L_a$ là điện trở và điện cảm phần ứng, $K_e, K_m$ là hệ số sức điện động và hệ số mô men, $I_m$ là mô men quán tính rô-to, $\tau_i$ là mô men phát động tại khớp chủ động.
  • Lý thuyết kỳ dị ma trận Jacobi: Dựa trên phương trình liên kết tổng quát cấp vận tốc $J_x \dot{x} + J_q \dot{q} = 0$:
    • Kỳ dị loại 1 (động học ngược): $\det(J_q) = 0$, khâu chấp hành mất bậc tự do di chuyển theo một số hướng, các chân đạt biên vùng làm việc (chân duỗi thẳng hoặc gấp hoàn toàn).
    • Kỳ dị loại 2 (động học thuận): $\det(J_x) = 0$, bàn máy động tăng thêm bậc tự do cục bộ, mất tính cứng vững (ba khâu dẫn sau đồng quy hoặc song song).
    • Kỳ dị loại 3 (kỳ dị hỗn hợp): Đồng thời thỏa mãn $\det(J_x) = 0$ và $\det(J_q) = 0$.
  • Kỹ thuật ước lượng động học và ổn định vi phân: Áp dụng phương trình hiệu chỉnh đạo hàm liên kết $\dot{f}(q) = -K f(q)$ (với ma trận đường chéo $K > 0$) để thiết lập bộ ước lượng biến khớp bị động $\hat{y}$ với nghiệm sai số dạng $e_i(t) = e_i(0)\exp(-k_{ii}t)$ hội tụ tiệm cận về 0, loại bỏ hiện tượng trôi dạt do tích phân số.
  • Thuật toán tiến hóa và điều khiển phi tuyến:
    • Giải thuật di truyền (GA): Sử dụng quần thể $4n$ cá thể, toán tử chọn lọc, lai ghép và đột biến để tìm điều kiện ban đầu $q(0)$ tối ưu và tối ưu hóa tham số mặt trượt.
    • Điều khiển trượt (SMC): Thiết lập mặt trượt $s(e, \dot{e}) = \dot{e} + \lambda e$, sử dụng hàm Lyapunov $V = \frac{1}{2}s^T s$ bảo đảm điều kiện trượt $\dot{V} \le -\eta |s|$.
    • Điều khiển logic mờ (Fuzzy Logic): Sử dụng luật hợp thành Max-Min, cấu trúc hàm thuộc tam giác và giải mờ bằng phương pháp trọng tâm (Center of Gravity) nhằm mờ hóa hàm dấu $\text{sign}(s)$, triệt tiêu hiện tượng rung chấn (chattering) và thích nghi với độ bất định tham số khối lượng/mô men quán tính.

Dữ liệu mô hình mô phỏng: Nghiên cứu sử dụng mô hình hình học và thông số vật lý của robot song song phẳng 3RRR với các giá trị chuẩn:

Thành phần kết cấu Ký hiệu tham số Giá trị danh định Đơn vị
Khoảng cách giữa các bệ cố định $L_0$ 1,072 $\text{m}$
Chiều dài khâu chủ động 1 $l_{1,i}$ 0,250 $\text{m}$
Chiều dài khâu bị động 2 $l_{2,i}$ 0,250 $\text{m}$
Kích thước cạnh bàn máy động $b$ 0,100 $\text{m}$
Khối lượng khâu 1 $m_{1,i}$ 0,850 $\text{kg}$
Khối lượng khâu 2 $m_{2,i}$ 0,185 $\text{kg}$
Khối lượng bàn máy động $m_3$ 0,500 $\text{kg}$
Mô men quán tính khâu 1 $I_{C1,i}$ $0,0044$ $\text{kg}\cdot\text{m}^2$
Mô men quán tính khâu 2 $I_{C2,i}$ $0,00096$ $\text{kg}\cdot\text{m}^2$
Mô men quán tính bàn máy động $J_{C3}$ $0,00672$ $\text{kg}\cdot\text{m}^2$
Tỉ số truyền hộp giảm tốc $r_g$ 10,01 -
Hằng số mô men động cơ DC $K_m$ 3,00 $\text{N}\cdot\text{m}/\text{A}$
Hằng số sức điện động động cơ $K_e$ 0,02 $\text{V}\cdot\text{s}/\text{rad}$
Mô men quán tính rô-to động cơ $I_m$ 0,0001 $\text{kg}\cdot\text{m}^2$
Điện trở mạch phần ứng $R_a$ 1,50 $\Omega$

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

Chương 1: Tổng quan

Chương 1 giới thiệu tổng quát về nguồn gốc, cấu trúc, phân loại và ứng dụng của các hệ thống robot song song trong công nghiệp và kỹ thuật y sinh. Tác giả phân tích chi tiết đặc thù động học mạch vòng kín, so sánh tính năng kỹ thuật với robot chuỗi, đồng thời nêu rõ các thách thức cốt lõi liên quan đến sự suy biến không gian làm việc tại các vị trí kỳ dị.

Chương này phân loại hệ thống kỳ dị thành 3 nhóm cơ bản dựa trên cấu trúc hình học và ma trận đạo hàm liên kết: kỳ dị loại 1 do cấu hình chân duỗi thẳng hoặc gập gọn làm $\det(J_q) = 0$; kỳ dị loại 2 do các khâu sau đồng quy tại một điểm hoặc song song làm $\det(J_x) = 0$; kỳ dị loại 3 là trạng thái hỗn hợp thỏa mãn đồng thời $\det(J_q) = 0$ và $\det(J_x) = 0$.

Tổng quan các nghiên cứu trong và ngoài nước được tổng hợp cụ thể, qua đó chỉ ra các thiếu sót về mô hình hóa tích hợp hệ thống dẫn động, bài toán vượt kỳ dị và phương pháp điều khiển không gian thao tác. Trên cơ sở đó, chương 1 xác định rõ mục tiêu, đối tượng khảo sát là robot song song phẳng 3RRR dẫn động bằng động cơ điện một chiều và xác lập hệ thống 05 nhiệm vụ nghiên cứu cho các chương tiếp theo.

Chương 2: Phân tích kỳ dị động học robot song song

Chương 2 tập trung giải quyết bài toán động học và khảo sát các cấu hình kỳ dị của robot song song:

  • Thiết lập hệ phương trình liên kết động học ở mức vị trí: $$f_i(q_{ai}, b_i) = l_i(x), \quad i=1,2,3$$ Đối với robot 3RRR, hệ phương trình liên kết tọa độ giữa các gốc bệ $O_i(x_{Oi}, y_{Oi})$, khớp nối $A_i, B_i$ và tâm bàn máy $P(x_P, y_P)$ với góc nghiêng $\alpha$ được biểu diễn: $$x_{Oi} + l_{1,i}\cos\theta_i + l_{2,i}\cos(\theta_i + \phi_i) = x_P + b\cos(\alpha + \gamma_i)$$ $$y_{Oi} + l_{1,i}\sin\theta_i + l_{2,i}\sin(\theta_i + \phi_i) = y_P + b\sin(\alpha + \gamma_i)$$
  • Phân tích động học ngược: Xây dựng phương pháp giải tích dựa trên cấu trúc hình học tam giác, áp dụng định lý hàm số cosin để tính toán góc quay khớp chủ động $\theta_i$ và góc phụ $\phi_i$. Nghiên cứu chỉ ra rằng ứng với một vị trí bàn thao tác trong vùng làm việc, hệ thống tồn tại $2^3 = 8$ cấu hình hình học phân biệt, ký hiệu $\sigma_j = [\pm 1, \pm 1, \pm 1]$ ($j=1,\dots,8$).
  • Thuật toán số cho bài toán động học ngược tổng quát: Sử dụng ma trận Jacobi kết hợp hàm ổn định hóa sai số vi phân $e = -Ke$ để triệt tiêu sai số tích phân số theo thời gian: $$\dot{q} = -J_q^{-1} \left[ J_x \dot{x} + K f(q, x) \right]$$ Tìm nghiệm điều kiện ban đầu $q(0)$ bằng giải thuật di truyền (GA) kết hợp thuật toán lặp Newton - Raphson nhằm tránh nghiệm rơi vào cực tiểu địa phương ngoài vùng giới hạn khớp.
  • Xây dựng bộ ước lượng động học: Đề xuất thuật toán ước lượng trực tuyến các biến khớp bị động $\hat{y}$ và vận tốc $\dot{\hat{y}}$ mà không cần lắp đặt thêm cảm biến đo vật lý: $$\dot{\hat{y}} = -J_y^{-1}(q_a, \hat{y}) \left[ J_{qa}(q_a, \hat{y})\dot{q}_a + K f(q_a, \hat{y}) \right]$$
  • Bản đồ kỳ dị và giải pháp vượt kỳ dị động học: Khử các biến phụ để xây dựng phương trình ràng buộc bậc hai $(r_{Bi} - r_{Ai})^2 - l_2^2 = 0$. Thiết lập ma trận Jacobi phân tích định thức $\det(J_x)$ và $\det(J_q)$. Thuật toán vượt kỳ dị động học cho phép quỹ đạo chuyển động đi cắt qua điểm kỳ dị một cách liên tục mà không làm gián đoạn chương trình tính toán số.

Chương 3: Động lực học robot song song có kể đến hệ dẫn động

Chương 3 thiết lập mô hình động lực học hoàn chỉnh tích hợp giữa hệ cơ cấu nhiều vật và hệ dẫn động cơ điện:

  • Mô hình hóa hệ thống cơ điện: Ghép nối phương trình vi phân chuyển động của cơ cấu cơ khí phẳng với hệ phương trình vi phân mô tả mạch điện phần ứng và chuyển động quay của rô-to động cơ DC.
  • Phương trình vi phân chuyển động tổng quát với tọa độ suy rộng tối thiểu: $$M(q)\ddot{q} + C(q, \dot{q})\dot{q} + G(q) = \tau$$ Khi tích hợp hệ dẫn động qua hộp giảm tốc tỉ số truyền $r_g$, ma trận khối lượng mở rộng $M_e(q)$ và véc-tơ lực suy rộng $u$ được tính đến quán tính $I_m$ và các hệ số điện cơ $K_m, K_e, R_a$: $$M_e(q)\ddot{q} + C_e(q, \dot{q})\dot{q} + G(q) = B u$$
  • Giải quyết bài toán động lực học thuận và ngược: Thực hiện các mô phỏng kiểm chứng so sánh mô hình cơ học thuần túy với mô hình tích hợp hệ dẫn động. Kết quả chỉ ra rằng việc bỏ qua động cơ DC dẫn đến sai số mô men phát động đáng kể ở các giai đoạn tăng tốc nhanh.
  • Ổn định hóa phương trình liên kết: Ứng dụng phương pháp Baumgarte và phương pháp chiếu hiệu chỉnh để loại bỏ sự trôi dạt nghiệm liên kết hình học trong quá trình giải tích phân số dài hạn.
  • Thuật toán vượt kỳ dị động lực học: Sử dụng kỹ thuật mở rộng không gian trạng thái và hiệu chỉnh phương trình vi phân chuyển động bằng các đạo hàm bậc cao lân cận điểm kỳ dị, bảo đảm mô men điều khiển và điện áp không bị tăng vọt vô hạn khi robot chuyển động qua cấu hình có $\det(J_x) \to 0$.

Chương 4: Điều khiển trượt robot song song có kể đến hệ dẫn động

Chương 4 nghiên cứu và thiết kế các thuật toán điều khiển chuyển động cho robot song song tích hợp hệ dẫn động:

  • Khảo sát các bộ điều khiển truyền thống: Thực hiện mô phỏng đối chứng bộ điều khiển vị trí PD + bù trọng lực và bộ điều khiển PID dựa trên mô men tính toán, đánh giá chất lượng bám quỹ đạo và độ nhạy trước nhiễu.
  • Thiết kế bộ điều khiển trượt (SMC) trong không gian khớp:
    • Thiết lập véc-tơ mặt trượt trong không gian biến khớp: $s = \dot{e} + \Lambda e$ với $e = q - q_d$.
    • Tổng hợp luật điện áp điều khiển $u$ dựa trên hàm tiêu chuẩn Lyapunov nhằm bảo đảm tính ổn định tiệm cận toàn cục.
    • Tích hợp bộ ước lượng động học đã xây dựng ở Chương 2 vào vòng hồi tiếp điều khiển, kiểm chứng tính khả thi khi không có cảm biến đo góc khớp bị động.
  • Thiết kế bộ điều khiển trượt trực tiếp trong không gian thao tác:
    • Biến đổi hệ phương trình động lực học về hệ tọa độ tối thiểu đại diện cho vị trí và góc nghiêng của bàn máy động $x = [x_P, y_P, \alpha]^T$.
    • Xây dựng luật điều khiển trượt trực tiếp trên sai lệch không gian làm việc: $e_x = x - x_d$, $s_x = \dot{e}_x + \Lambda_x e_x$.
    • So sánh mô phỏng giữa điều khiển SMC và điều khiển PD trong không gian thao tác: Kết quả khẳng định điều khiển trượt trong không gian thao tác triệt tiêu sai lệch bám nhanh hơn, góc nghiêng bàn máy động đạt độ ổn định cao hơn, tuy nhiên vẫn còn tồn tại hiện tượng rung chấn tần số cao (chattering) do thành phần hàm dấu $\text{sign}(s)$.

Chương 5: Điều khiển trượt mờ robot song song có kể đến hệ dẫn động

Chương 5 phát triển bộ điều khiển trượt mờ (Fuzzy Sliding Mode Control - FSMC) nhằm loại bỏ chattering và thích nghi với các yếu tố bất định tham số:

  • Cơ sở cấu trúc mờ hóa: Xây dựng hệ thống mờ Mamdani với hai đầu vào là biến mặt trượt $s$ và đạo hàm mặt trượt $\dot{s}$, đầu ra là biến điều khiển hiệu chỉnh $v$. Sử dụng 7 tập mờ chuẩn cho mỗi biến: Âm Lớn (NB), Âm Vừa (NM), Âm Nhỏ (NS), Không (ZE), Dương Nhỏ (PS), Dương Vừa (PM), Dương Lớn (PB) cùng luật hợp thành Max-Min và giải mờ trọng tâm.
  • Bộ điều khiển FSMC bù nhiễu và bất định tham số: Thiết lập mô phỏng kiểm chứng độ bền vững khi thông số khối lượng và mô men quán tính của robot bị sai lệch so với giá trị thiết kế thực tế ở 2 kịch bản:
    • Trường hợp 1: Nhận dạng tham số đạt 80% giá trị chính xác.
    • Trường hợp 2: Nhận dạng tham số chỉ đạt 50% giá trị chính xác (bất định 50%). Kết quả mô phỏng cho thấy FSMC bù trừ hoàn toàn sai lệch mô hình, duy trì quỹ đạo bám chính xác của bàn máy động theo các trục $x, y$ và góc $\alpha$.
  • Tối ưu hóa bộ điều khiển FSMC bằng giải thuật di truyền (GA):
    • Ứng dụng GA để tìm kiếm bộ thông số tối ưu cho ma trận độ dốc mặt trượt $\Lambda$ và các trọng số hàm thuộc của bộ điều khiển mờ.
    • Mờ hóa thành phần gián đoạn $\text{sign}(s)$ thành biến trơn sinh bởi bộ suy diễn mờ.
    • Đồ thị mô phỏng điện áp điều khiển và mô men phát động chứng minh hiện tượng rung giật chattering bị triệt tiêu hoàn toàn, tín hiệu điều khiển chuyển tiếp êm thuận, sai số bám quỹ đạo của các khớp chủ động giảm rõ rệt so với bộ điều khiển SMC kinh điển.

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

Các kết quả và đóng góp mới của luận án được tổng hợp cụ thể:

  • Về mặt lý luận:
    1. Xây dựng hoàn chỉnh mô hình toán học động lực học cho robot song song phẳng có kể đến đầy đủ các đặc tính phi tuyến của hệ dẫn động cơ điện (động cơ điện một chiều và hộp giảm tốc), khắc phục sự thiếu hụt của các mô hình cơ học thuần túy trước đây.
    2. Đề xuất phương pháp phân tích, nhận diện và phân vùng kỳ dị dựa trên ma trận Jacobi đại số khử biến phụ, đồng thời phát triển thuật toán số cho phép quỹ đạo động học và động lực học vượt qua các điểm kỳ dị mà không gây mất ổn định hay phân kỳ toán học.
    3. Xây dựng thành công thuật toán ước lượng động học cho robot song song có khả năng tự động triệt tiêu sai số trôi dạt bằng kỹ thuật phản hồi vi phân ổn định, cho phép xác định chính xác trạng thái vị trí và vận tốc của các khớp bị động.
    4. Phát triển cấu trúc điều khiển trượt mờ (FSMC) thích nghi kết hợp giải thuật di truyền (GA), giải quyết đồng thời hai nhược điểm cốt lõi của điều khiển trượt truyền thống: loại bỏ hiện tượng chattering và bảo đảm tính bền vững trước sự bất định tham số lên tới 50%.
  • Về mặt thực tiễn:
    1. Cung cấp bộ công cụ thuật toán và chương trình tính toán mô phỏng hoàn chỉnh trên máy tính phục vụ trực tiếp cho công tác tính toán, kiểm nghiệm và thiết kế các dạng robot song song phẳng tương tự.
    2. Giải pháp ước lượng động học giúp giảm số lượng cảm biến vật lý cần lắp ráp trên các khớp bị động, giúp tối ưu hóa kết cấu cơ khí, giảm khối lượng chuyển động và hạ giá thành chế tạo thiết bị.
    3. Thuật toán điều khiển trực tiếp trong không gian thao tác tạo tiền đề cho việc tích hợp các hệ thống đo lường hiện đại (như cảm biến thị giác máy tính, cảm biến laser) vào điều khiển thời gian thực cho robot công nghiệp.

Kiến nghị và giải pháp đề xuất:

  • Khuyến nghị ứng dụng mô hình động lực học tích hợp hệ dẫn động khi thiết kế hệ thống điều khiển bám quỹ đạo tốc độ cao nhằm tránh quá tải động cơ.
  • Đề xuất áp dụng bộ điều khiển trượt mờ tối ưu bằng GA cho các cơ cấu robot song song làm việc trong môi trường chịu tải trọng biến thiên liên tục hoặc có thông số động học khó xác định chính xác.

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

Hạn chế của nghiên cứu:

  • Toàn bộ các kết quả khảo sát, đánh giá tính đúng đắn của mô hình động lực học, thuật toán ước lượng động học và các bộ điều khiển trượt mờ mới chỉ dừng lại ở mức độ thực nghiệm số thông qua mô phỏng trên máy tính, chưa được kiểm chứng trên mô hình thực nghiệm vật lý thực tế.
  • Đối tượng nghiên cứu giới hạn trên hệ thống robot song song phẳng 3 bậc tự do (3RRR), chưa mở rộng khảo sát cho các dạng robot song song không gian có cấu trúc phức tạp hơn (như robot Stewart 6 bậc tự do hay robot Delta 3D).
  • Hệ dẫn động khảo sát chuyên biệt cho loại động cơ điện một chiều thông thường kết hợp hộp giảm tốc cơ khí, chưa tính đến các dòng động cơ xoay chiều không chổi than (BLDC, PMSM) hoặc hệ dẫn động trực tiếp (Direct-drive).

Hướng nghiên cứu tiếp theo:

  • Chế tạo mô hình thực nghiệm vật lý robot song song phẳng 3RRR để đo đạc, kiểm chứng và hiệu chỉnh các kết quả mô phỏng số.
  • Mở rộng thuật toán phân tích kỳ dị, ước lượng động học và bộ điều khiển FSMC tối ưu sang các hệ thống robot song song không gian nhiều bậc tự do.
  • Ứng dụng xử lý ảnh và công nghệ đo laser trực tiếp trong không gian công tác để hoàn thiện hệ thống điều khiển bám thời gian thực.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị chuyên môn cao đối với các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và sinh viên các ngành Kỹ thuật Cơ học, Cơ kỹ thuật, Cơ điện tử và Tự động hóa: Tham khảo phương pháp xây dựng phương trình động lực học hệ nhiều vật vòng kín, kỹ thuật ổn định hóa ràng buộc liên kết (Baumgarte, phép chiếu hiệu chỉnh) và quy trình thiết kế các bộ điều khiển phi tuyến bền vững (SMC, FSMC, GA).
  • Giảng viên và nhà nghiên cứu: Sử dụng nội dung Chương 2 và Chương 3 để làm tài liệu giảng dạy chuyên đề về Động lực học Robot nâng cao, phân tích cấu trúc ma trận Jacobi và xử lý các trạng thái kỳ dị động học/động lực học.
  • Kỹ sư thiết kế và phát triển hệ thống tự động: Ứng dụng thuật toán ước lượng động học để giảm thiểu cảm biến phần cứng trong thiết kế cơ cấu máy song song; khai thác các bảng thông số động cơ và quy trình thiết kế bộ điều khiển trượt mờ để ứng dụng trên các máy công cụ CNC, cơ cấu định vị vi mô hoặc thiết bị mô phỏng chuyển động.

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

1. Luận án phân loại các dạng kỳ dị của robot song song phẳng 3RRR như thế nào và điều kiện toán học tương ứng là gì? Luận án phân loại kỳ dị thành 3 nhóm cơ bản dựa trên ma trận Jacobi của phương trình liên kết $J_x \dot{x} + J_q \dot{q} = 0$:

  • Kỳ dị loại 1 (kỳ dị động học ngược): Xảy ra khi $\det(J_q) = 0$, xuất hiện khi một chân robot duỗi thẳng hoặc gập lại hoàn toàn, tâm bàn máy nằm trên biên không gian làm việc, robot mất bậc tự do di chuyển.
  • Kỳ dị loại 2 (kỳ dị động học thuận): Xảy ra khi $\det(J_x) = 0$, xuất hiện khi ba khâu sau của các chân hoặc đồng quy tại một điểm hoặc song song với nhau, làm bàn máy động tăng thêm bậc tự do cục bộ và mất khả năng chịu lực.
  • Kỳ dị loại 3 (kỳ dị hỗn hợp): Xảy ra khi thỏa mãn đồng thời $\det(J_x) = 0$ và $\det(J_q) = 0$, kết hợp giữa hiện tượng khâu sau song song/đồng quy và chân robot bị duỗi thẳng.

2. Tại sao việc mô hình hóa động lực học cần phải tính đến hệ dẫn động (động cơ DC và hộp giảm tốc)? Hầu hết các nghiên cứu truyền thống chỉ mô tả chuyển động của các khâu khớp cơ học thuần túy. Tuy nhiên, hệ dẫn động bao gồm rô-to động cơ và hộp giảm tốc có quán tính quay riêng ($I_m$), tổn hao ma sát và các đặc tính điện - cơ phi tuyến (sức điện động cảm ứng $K_e$, điện trở $R_a$, mô men phát động $K_m$). Khi robot chuyển động ở tốc độ và gia tốc cao, các yếu tố này làm biến đổi ma trận khối lượng và mô men quán tính biểu kiến của toàn hệ thống, dẫn đến sai lệch lớn trong tính toán lực/mô men điều khiển nếu chỉ dùng mô hình cơ học thuần túy.

3. Thuật toán ước lượng động học đề xuất trong luận án giải quyết vấn đề gì và dựa trên nguyên lý nào? Trong thực tế, việc lắp đặt cảm biến (encoder) tại tất cả các khớp bị động làm tăng khối lượng, phức tạp kết cấu và tăng chi phí chế tạo. Thuật toán ước lượng động học cho phép xác định gián tiếp trạng thái của các khớp bị động $\hat{y}$ và vận tốc $\dot{\hat{y}}$ chỉ từ thông tin cảm biến của khớp chủ động $q_a, \dot{q}a$. Để chống lại sự tích lũy sai số gây trôi dạt nghiệm khi tích phân số, thuật toán sử dụng phương trình hiệu chỉnh sai số vi phân $\dot{f}(q) = -K f(q)$ (với $K$ là ma trận đường chéo dương), bảo đảm sai số liên kết hình học hội tụ tiệm cận về 0 theo hàm mũ $e_i(t) = e_i(0)\exp(-k{ii}t)$.

4. Giải pháp điều khiển trượt mờ (FSMC) tối ưu bằng giải thuật di truyền (GA) có ưu điểm gì so với bộ điều khiển trượt (SMC) thông thường? Bộ điều khiển trượt kinh điển (SMC) sử dụng hàm gián đoạn $\text{sign}(s)$ nên thường gây ra hiện tượng rung giật tần số cao (chattering), gây nóng động cơ và kích thích các dao động cơ khí không mong muốn. Bộ điều khiển FSMC mờ hóa thành phần gián đoạn này thông qua bảng suy diễn mờ Mamdani và luật hợp thành Max-Min, tạo ra tín hiệu điều khiển chuyển tiếp trơn tru. Đồng thời, giải thuật di truyền (GA) tối ưu hóa các tham số mặt trượt và hàm thuộc, giúp hệ thống duy trì khả năng bám quỹ đạo chính xác ngay cả khi thông số mô hình bị bất định tới 50%.

5. Phương pháp nào được sử dụng để giải bài toán động học ngược cho robot 3RRR và tìm điều kiện ban đầu $q(0)$? Đối với robot phẳng 3RRR, bài toán động học ngược được giải bằng hai cách:

  • Phương pháp giải tích: Dựa trên phân tích hình học tam giác và định lý hàm số cosin, xác định được 8 cấu hình nghiệm khả dĩ $\sigma_j = [\pm 1, \pm 1, \pm 1]$.
  • Phương pháp số: Đạo hàm hệ phương trình liên kết để tìm vận tốc qua ma trận Jacobi và hàm ổn định hóa vi phân. Giá trị đầu $q(0)$ được xác định thông qua việc giải bài toán tối ưu hóa hàm khoảng cách giới hạn khớp bằng giải thuật di truyền (GA) kết hợp với thuật toán lặp Newton - Raphson để đạt độ chính xác cao.

Kết luận

Luận án tiến sĩ của tác giả Vũ Đức Vương đã giải quyết một cách có hệ thống và hoàn chỉnh các vấn đề cốt lõi về phân tích động học, mô hình hóa động lực học có kể đến hệ dẫn động và thiết kế điều khiển phi tuyến cho robot song song phẳng 3RRR. Công trình đã đóng góp các giải pháp học thuật quan trọng bao gồm: phương pháp vượt kỳ dị động học và động lực học, bộ ước lượng động học chống trôi dạt liên kết và cấu trúc điều khiển trượt mờ thích nghi tối ưu bằng giải thuật di truyền. Toàn bộ các đề xuất lý thuyết đã được kiểm chứng tính đúng đắn và hiệu quả thông qua hệ thống mô phỏng số chi tiết. Kết quả của luận án là nguồn tài liệu tham khảo có giá trị cao cho công tác nghiên cứu, đào tạo và ứng dụng thiết kế robot song song trong thực tiễn.