Tổng quan luận án

Luận án tiến sĩ với đề tài "Động lực học và điều khiển tay máy có khâu đàn hồi chuyển động tuần hoàn" do nghiên cứu sinh Đinh Công Đạt thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Văn Khang và PGS. Nguyễn Quang Hoàng, bảo vệ tại Đại học Bách khoa Hà Nội năm 2023 thuộc ngành Cơ học (mã số: 9440109).

Tính cấp thiết và khoảng trống nghiên cứu

Trong kỹ thuật robot hiện đại, xu hướng thiết kế chuyển dịch từ các cơ cấu vật rắn nặng nề, cồng kềnh sang các tay máy có kết cấu thanh mảnh, trọng lượng nhẹ nhằm giảm năng lượng tiêu thụ, nâng cao vận tốc làm việc và tăng tính an toàn khi tương tác. Tuy nhiên, tính đàn hồi của các khâu thanh mảnh dẫn đến rung động kết cấu trong quá trình làm việc với vận tốc cao. Hai hệ quả trực tiếp xuất hiện khi bỏ qua yếu tố đàn hồi trong mô hình toán học là: sai số trong việc xác định mô-men/lực phát động của nguồn dẫn động và sai số định vị tại khâu thao tác cuối. Khi các cơ cấu này thực hiện chuyển động tuần hoàn ở tốc độ cao, phương trình vi phân chuyển động chứa các hệ số biến thiên theo thời gian có chu kỳ, dẫn đến hiện tượng dao động tham số và nguy cơ mất ổn định động lực (cộng hưởng tham số). Vấn đề đặt ra là phải xây dựng mô hình toán học chính xác, xác định vùng làm việc ổn định và điều khiển bù trừ dao động nhằm đáp ứng yêu cầu bám quỹ đạo mong muốn.

Mục tiêu nghiên cứu

  1. Xây dựng mô hình cơ học và mô hình toán học giải tích cho các lớp tay máy có khâu đàn hồi (tập tọa độ suy rộng bao gồm tọa độ khâu dẫn và tọa độ đàn hồi).
  2. Tuyến tính hóa hệ phương trình vi phân chuyển động phi tuyến xung quanh chuyển động cơ bản (chuyển động của mô hình rắn tương ứng) bằng khai triển Taylor của hàm ma trận và hàm véc tơ theo biến véc tơ.
  3. Xác định các tham số điều khiển tối ưu nhằm đảm bảo cơ hệ hoạt động trong vùng ổn định động lực, tránh hiện tượng cộng hưởng tham số.
  4. Tính toán dao động tuần hoàn và đề xuất phương pháp giải gần đúng bài toán động lực học ngược để xác định lực/mô-men phát động sao cho khâu thao tác bám sát quỹ đạo định trước.

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

  • Đối tượng nghiên cứu: Động lực học và điều khiển tay máy phẳng dạng chuỗi có chứa khâu đàn hồi, cụ thể gồm 3 mô hình:
    1. Tay máy một khâu đàn hồi chuyển động quay có mang tải trọng tại khâu thao tác cuối.
    2. Tay máy hai khâu T-R (khâu 1 tịnh tiến coi là rắn, khâu 2 quay coi là đàn hồi).
    3. Tay máy hai khâu R-R (khâu 1 quay coi là rắn, khâu 2 quay coi là đàn hồi mang tải trọng cuối).
  • Phạm vi nghiên cứu:
    • Khảo sát chuyển động tuần hoàn trong mặt phẳng.
    • Khâu đàn hồi được mô hình hóa theo lý thuyết dầm Euler-Bernoulli chịu uốn ngang (bỏ qua biến dạng dọc, xoắn và biến dạng cắt).
    • Khảo sát tính ổn định động lực học theo lý thuyết Floquet kết hợp phương pháp quy hoạch thực nghiệm Taguchi và tính toán nghiệm số trên phần mềm Maple, Matlab.

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

Các hướng nghiên cứu trước đây

Tác giả luận án đã hệ thống hóa các công trình nghiên cứu trong và ngoài nước theo các nhóm vấn đề chuyên sâu:

  1. Mô hình hóa phần tử đàn hồi và tay máy một khâu đàn hồi:

    • Sử dụng phương pháp khai triển theo các dạng riêng (Assumed Modes Method - AMM) kết hợp phương trình Lagrange loại 2: Canon và Schmitz (1984), Bayo (1987), Hastings và Book (1987), Tso và cộng sự (1999).
    • Sử dụng nguyên lý Hamilton mở rộng và phương trình Newton-Euler: Rakhsha và Goldenberg (1985), Barbieri và Ozguner (1988).
    • Khảo sát khâu đàn hồi kết hợp khớp tịnh tiến: Tabarrok và cộng sự (1988), Love và cộng sự (2004).
    • Khảo sát đồng thời biến dạng dài và uốn: Baruh và Tadikonda (1989).
    • Ứng dụng phương pháp phần tử hữu hạn (FEM): Nagarajan và Turcic (1990), Bricout và cộng sự (1990), Alberts và cộng sự (1992), Moulin và Bayo (1991), Mohamed và Tokhi (2004), Lee và Wang (2001), Theodore và Ghosal (1995).
    • Mô hình tham số tập trung (LPM): Zhu và cộng sự (1999). Ảnh hưởng của lực cắt và quán tính quay: Morris và Madani (1996), Bayo (1987). Tối ưu hóa hình dạng khâu: Asada và cộng sự (1991), Wang và Russel (1994).

  2. Tay máy hai khâu và nhiều khâu đàn hồi:

    • Ứng dụng phương pháp Lagrange: Ower và Vegte (1987), Sunada và Dubowsky (1981), Dado và Soni (1986), Fukuda (1985), Low và Vidyasagar (1988). Lee (1994) đã chỉ ra khiếm khuyết của mô hình Lagrange thông thường khi cho phép biến dạng dài không mong muốn và đề xuất mô hình cải tiến.
    • Ứng dụng nguyên lý Hamilton và phương pháp Kane: Singh (1993), Zhang và cộng sự (1996), Dogan và Iftar (2002), Meghdari (1991), Baruh (1989). Morris và Madani (1996) chỉ ra trong robot công nghiệp 6 bậc tự do, thường có hai khâu chính gây quán tính lớn cần mô hình hóa đàn hồi.

  3. Động lực học ngược và điều khiển tay máy đàn hồi:

    • Động lực học ngược trong miền tần số và thời gian: Bayo (1987, 1989), Kanaoka và Yoshikawa (1998), Green và Sasiadek (2001, 2003).
    • Điều khiển tay máy đàn hồi: Cannon và Schmitz (1984) khởi xướng thí nghiệm điều khiển phản hồi vị trí đầu mút; Benosman và Vey (2004) tổng quan các mục tiêu điều khiển (triệt tiêu dao động dư, chuyển trạng thái nghỉ, bám quỹ đạo).

  4. Tình hình nghiên cứu trong nước:

    • Nguyễn Văn Khang và cộng sự (1997, 2001, 2005) nghiên cứu dao động cơ cấu cam có cần đàn hồi và cơ cấu phẳng có khâu đàn hồi.
    • Vũ Văn Khiêm (2004) nghiên cứu dao động tuần hoàn cơ cấu phẳng bằng phương pháp tách cấu trúc.
    • Nguyễn Văn Khang, Nguyễn Sỹ Nam, Nguyễn Phong Điền (2004, 2006, 2008) phân tích động lực học và tuyến tính hóa hệ phương trình vi phân chuyển động của cơ cấu 4 khâu, 6 khâu bản lề phẳng đàn hồi.
    • Dương Xuân Biên (2018), Chu Anh Mỳ và cộng sự (2017, 2020) thiết lập phương trình vi phân cho tay máy đàn hồi bằng FEM với giả thiết biến dạng nhỏ.
    • Nguyễn Quang Hoàng và cộng sự (2018, 2021) ứng dụng FEM kết hợp điều khiển PD+ cho tay máy 1 và 2 khâu đàn hồi.

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

  • Áp dụng phương pháp hệ quy chiếu đồng hành (Floating Frame of Reference Formulation - FFRF) để thiết lập hệ phương trình vi phân chuyển động phi tuyến tường minh cho hệ tay máy phẳng có khâu đàn hồi dạng chuỗi.
  • Xây dựng công thức toán học tổng quát cho khai triển Taylor của hàm ma trận và hàm véc tơ theo biến véc tơ để tuyến tính hóa phương trình chuyển động quanh chuyển động cơ bản.
  • Giải quyết bài toán dao động tham số bằng cách kết hợp lý thuyết Floquet với phương pháp Taguchi để xác định bộ tham số điều khiển tối ưu ($K_P, K_D$) giúp hệ tránh được vùng cộng hưởng tham số.
  • Đề xuất quy trình giải gần đúng bài toán động lực học ngược cho tay máy chuyển động tuần hoàn nhằm xác định mô-men phát động bù trừ sai lệch đàn hồi.

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

Cơ sở lý thuyết và khung phân tích

  • Phương pháp hệ quy chiếu đồng hành: Vị trí tuyệt đối của một điểm $P$ bất kỳ trên khâu đàn hồi thứ $i$ được xác định qua biểu thức: $$\mathbf{r}i = \mathbf{R}i + \mathbf{A}i (\mathbf{u}{i0} + \mathbf{u}{if})$$ Trong đó: $\mathbf{R}i$ là véc tơ vị trí gốc hệ tọa độ đồng hành gắn với khâu $i$; $\mathbf{A}i$ là ma trận cosin chỉ hướng; $\mathbf{u}{i0}$ là véc tơ vị trí khi chưa biến dạng; $\mathbf{u}{if}$ là véc tơ biến dạng đàn hồi. Khi triệt tiêu biến dạng ($\mathbf{u}{if} = \mathbf{0}$), mô hình trở về dạng vật rắn tuyệt đối.
  • Biểu diễn biến dạng bằng phương pháp Ritz - Galerkin: Khâu đàn hồi được mô hình hóa theo lý thuyết dầm Euler-Bernoulli. Chuyển vị uốn ngang tương đối $w(x, t)$ được biểu diễn dưới dạng chuỗi tách biến: $$w(x, t) = \sum_{i=1}^N X_i(x) q_{ei}(t)$$ Với $X_i(x)$ là các hàm dạng thỏa mãn điều kiện biên của dầm công-xôn (ngàm - tự do), suy ra từ nghiệm của phương trình đặc trưng: $1 + \cos(\beta_i l) \cosh(\beta_i l) = 0$. Tọa độ $q_{ei}(t)$ là các tọa độ suy rộng đàn hồi phụ thuộc thời gian.
  • Thiết lập phương trình chuyển động: Sử dụng phương trình vi phân Lagrange loại 2 cho hệ cơ học chịu liên kết holonom: $$\frac{d}{dt}\left(\frac{\partial T}{\partial \dot{\mathbf{q}}}\right) - \frac{\partial T}{\partial \mathbf{q}} + \frac{\partial \Pi}{\partial \mathbf{q}} = \mathbf{Q}^$$ Trong đó: $T$ là tổng động năng của hệ (khâu dẫn, đĩa khớp, dầm đàn hồi, tải trọng cuối); $\Pi$ là tổng thế năng (thế năng biến dạng đàn hồi dầm $\Pi_{dh} = \frac{1}{2}\int_0^l EI \left(\frac{\partial^2 w}{\partial x^2}\right)^2 dx$ và thế năng trọng trường $\Pi_g$); $\mathbf{Q}^$ là véc tơ lực suy rộng không có thế (bao gồm lực/mô-men phát động $\boldsymbol{\tau}$ và lực cản nhớt $\mathbf{F}_d, \mathbf{M}_d$).
Ký hiệu toán học Ý nghĩa vật lý / cơ học
$\mathbf{q} = [\mathbf{q}_a^T, \mathbf{q}_e^T]^T$ Véc tơ tọa độ suy rộng gồm tọa độ khâu dẫn $\mathbf{q}_a$ và tọa độ đàn hồi $\mathbf{q}_e$
$E, I, \rho, A$ Mô đun đàn hồi, mô-men quán tính tiết diện, khối lượng riêng, diện tích tiết diện
$C_i, D_i, n_{ij}, h_{ij}$ Các tích phân hình học dạng: $C_i = \int_0^l X_i dx$, $D_i = \int_0^l x X_i dx$, $n_{ij} = \int_0^l X_i X_j dx$, $h_{ij} = \int_0^l X''_i X''_j dx$
$\rho_j, \lambda_i$ Các nhân tử Floquet và số mũ Floquet đặc trưng cho tính ổn định của hệ vi phân tuần hoàn
$K_P, K_D$ Ma trận hệ số khuếch đại vị trí và vận tốc của bộ điều khiển phản hồi
$\text{SNR}$ Tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio) trong phương pháp Taguchi
  • Lý thuyết ổn định hệ phương trình vi phân tuyến tính hệ số tuần hoàn (Lý thuyết Floquet): Hệ phương trình vi phân tuyến tính biểu diễn sai lệch $\mathbf{y}(t)$ quanh quỹ đạo cơ bản có dạng: $$\dot{\mathbf{y}}(t) = \mathbf{A}(t)\mathbf{y}(t), \quad \text{với } \mathbf{A}(t + T) = \mathbf{A}(t)$$ Ma trận chuyển tiếp trạng thái qua một chu kỳ $\mathbf{\Phi}(T)$ có các giá trị riêng là các nhân tử Floquet $\rho_j$. Điều kiện ổn định tiệm cận của hệ chuyển động tuần hoàn là: $$|\rho_j| < 1 \quad \forall j$$
  • Phương pháp Taguchi: Sử dụng mảng trực giao (bảng $L_9$) và hàm mục tiêu tỷ số tín hiệu trên nhiễu $\text{SNR}$ theo đặc tính "càng nhỏ càng tốt" đối với chỉ số $|\rho|{\max}$ nhằm xác định cấu hình thông số điều khiển $K_P, K_D$ tối ưu: $$\text{SNR} = -10 \log{10}\left(\frac{1}{n}\sum_{k=1}^n y_k^2\right)$$

Phương pháp nghiên cứu và công cụ thực hiện

  • Phương pháp giải tích: Thiết lập biểu thức giải tích động năng, thế năng và các ma trận hệ số thông qua phương trình Lagrange loại 2; xây dựng công thức khai triển vi phân Taylor ma trận.
  • Phương pháp số và mô phỏng:
    • Tích phân số phương trình vi phân bằng phương pháp Runge-Kutta.
    • Sử dụng phần mềm Maple cho các biến đổi giải tích hình thức phức tạp và đạo hàm véc tơ/ma trận.
    • Sử dụng phần mềm Matlab / Simulink để tính toán nhân tử Floquet, tối ưu hóa Taguchi và mô phỏng số đáp ứng động lực học của các mô hình tay máy.

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

Chương 1: Tổng quan và cơ sở lý thuyết

Chương này tổng quan về các phương pháp mô hình hóa yếu tố đàn hồi trong cơ cấu robot (AMM, FEM, RFEM, LPM) và phân tích sâu sắc các công trình nghiên cứu trên thế giới và tại Việt Nam. Về mặt cơ sở lý thuyết, chương thiết lập 3 nền tảng toán học - cơ học chính:

  1. Phương pháp hệ quy chiếu đồng hành: Xác định hệ tọa độ đồng hành chuyển động cùng khâu đàn hồi, cho phép tách chuyển động vật rắn lớn và biến dạng đàn hồi nhỏ.
  2. Khai triển Taylor theo biến véc tơ: Đưa ra công thức khai triển tổng quát cho hàm véc tơ và hàm ma trận theo biến véc tơ làm công cụ giải tích phục vụ tuyến tính hóa phương trình chuyển động.
  3. Lý thuyết ổn định của hệ phương trình vi phân tuyến tính hệ số tuần hoàn: Trình bày định lý Floquet, ma trận đơn đạo (monodromy matrix), định nghĩa nhân tử Floquet $\rho_j$ và số mũ Floquet $\lambda_i$ làm chuẩn mực đánh giá ổn định động lực học.

Chương 2: Phương trình vi phân chuyển động của tay máy đàn hồi

Chương 2 tập trung thiết lập phương trình vi phân chuyển động phi tuyến và thực hiện tuyến tính hóa giải tích quanh quỹ đạo chuyển động cơ bản cho 3 cơ hệ cụ thể:

  1. Tay máy một khâu đàn hồi chuyển động quay:

    • Cấu tạo gồm động cơ quay khâu $OE$ (chiều dài $l$, tiết diện $A$), tải trọng tập trung $m_E$ tại đầu $E$.
    • Tọa độ suy rộng: $\mathbf{q} = [q_a, q_{e1}, \dots, q_{eN}]^T$. Phương trình vi phân chuyển động khi lấy $N=1$ dạng: $$\begin{cases} \left[J_1 + m_E l^2 + \frac{1}{3}\rho A l^3 + (m_E X_1^2(l) + \rho A n_{11})q_{e1}^2\right]\ddot{q}a + 2(m_E X_1^2(l) + \rho A n{11})\dot{q}a q{e1}\dot{q}{e1} \ \quad + (\rho A D_1 + m_E l X_1(l))\ddot{q}{e1} + \left(\frac{1}{2}m_{OE} + m_E\right)gl\cos q_a - m_E g X_1(l)q_{e1}\sin q_a - \rho A g C_1 q_{e1}\sin q_a + M_d = \tau \ (m_E X_1^2(l) + \rho A n_{11})\ddot{q}{e1} + (m_E l X_1(l) + \rho A D_1)\ddot{q}a - (m_E X_1^2(l) + \rho A n{11})\dot{q}a^2 q{e1} \ \quad + (m_E X_1(l) + \rho A C_1)g\cos q_a + EI h{11} q_{e1} = 0 \end{cases}$$

  2. Tay máy hai khâu đàn hồi T-R:

    • Gồm khâu 1 tịnh tiến theo phương thẳng đứng (tọa độ $q_{a1}$, khối lượng $m_1$), đầu mang đĩa $B$ (khối lượng $m_B$, bán kính $r$). Khâu 2 quay quanh trục qua đĩa $B$ (tọa độ góc $q_{a2}$, chiều dài $l_2$, khối lượng $m_2$) là thanh đàn hồi.
    • Tọa độ suy rộng: $\mathbf{q} = [q_{a1}, q_{a2}, q_{e1}, \dots, q_{eN}]^T$.
    • Thiết lập tường minh biểu thức động năng, thế năng và hệ phương trình vi phân phi tuyến liên kết giữa chuyển động tịnh tiến - quay - uốn đàn hồi.

  3. Tay máy hai khâu đàn hồi R-R:

    • Gồm khâu 1 quay quanh gốc cố định $O$ (rắn, chiều dài $l_1$, tọa độ $q_{a1}$), đầu mang đĩa $B$ (khối lượng $m_B$, mô men quán tính $J_B$), khâu 2 quay quanh khớp $B$ (đàn hồi, chiều dài $l_2$, tọa độ góc $q_{a2}$, tải mút $m_E$).
    • Thiết lập hệ phương trình chuyển động đầy đủ mô tả tương tác phi tuyến phức tạp giữa hai chuyển động quay và các mode dao động uốn.

  4. Tuyến tính hóa quanh chuyển động cơ bản:

    • Chuyển động cơ bản $\mathbf{q}^0(t) = [\mathbf{q}_a^0(t)^T, \mathbf{0}^T]^T$ được định nghĩa là chuyển động của tay máy khi coi các khâu hoàn toàn rắn tuyệt đối.
    • Đặt $\mathbf{q}(t) = \mathbf{q}^0(t) + \mathbf{y}(t)$ với $\mathbf{y}(t) = [\Delta \mathbf{q}_a^T, \mathbf{q}_e^T]^T$ là véc tơ sai lệch nhỏ.
    • Sử dụng khai triển Taylor ma trận giải tích, hệ phương trình được đưa về dạng vi phân ma trận cấp hai tuyến tính biến thiên theo thời gian: $$\mathbf{M}_0(t)\ddot{\mathbf{y}} + \mathbf{C}_0(t)\dot{\mathbf{y}} + \mathbf{K}_0(t)\mathbf{y} = \mathbf{h}(t)$$ Khi chuyển động cơ bản có tính chu kỳ, các ma trận $\mathbf{M}_0, \mathbf{C}_0, \mathbf{K}_0$ là các hàm tuần hoàn chu kỳ $T$.

Chương 3: Điều khiển ổn định động lực và tính toán động lực học ngược tay máy có khâu đàn hồi

Chương 3 giải quyết trọn vẹn 3 bài toán ứng dụng cốt lõi:

  1. Điều khiển ổn định động lực học dựa trên phương pháp Taguchi:

    • Đưa tín hiệu điều khiển phản hồi $\Delta \boldsymbol{\tau} = -\mathbf{K}_P \Delta \mathbf{q}_a - \mathbf{K}_D \Delta \dot{\mathbf{q}}_a$ vào hệ thống tuyến tính hóa.
    • Chuyển hệ về dạng phương trình trạng thái $\dot{\mathbf{x}} = \mathbf{A}(t)\mathbf{x}$. Tính toán ma trận chuyển tiếp chu kỳ $\mathbf{\Phi}(T)$ để trích xuất các nhân tử Floquet $\rho_j$.
    • Xây dựng bài toán tối ưu thông số điều khiển $K_P, K_D$ sử dụng bảng trực giao $L_9$ của Taguchi. Lựa chọn mức tham số tối ưu để cực tiểu hóa $|\rho|_{\max} < 1$, đảm bảo mọi nghiệm nhiễu tắt dần nhanh theo thời gian.
    • Khảo sát mô phỏng và bảng số liệu nhân tử Floquet trên các dải vận tốc góc làm việc:
      • Tay máy 1 khâu: Khảo sát tại $\omega = 2, 4, 6, 8, 10\text{ rad/s}$.
      • Tay máy 2 khâu T-R: Khảo sát tại $\omega = 1, 2, 4, 8, 9\text{ rad/s}$.
      • Tay máy 2 khâu R-R: Khảo sát tại $\omega = 4, 6, 8, 10, 12\text{ rad/s}$.

  2. Tính toán dao động tuần hoàn:

    • Thiết lập thuật toán số tìm nghiệm tuần hoàn cưỡng bức $\mathbf{y}_p(t) = \mathbf{y}_p(t+T)$ của hệ phương trình vi phân tuyến tính có hệ số tuần hoàn.
    • Xác định chính xác điều kiện đầu $\mathbf{y}(0)$ để hệ nhanh chóng xác lập trạng thái dao động ổn định chu kỳ mà không bị bùng nổ biên độ do kích động tham số.

  3. Tính toán gần đúng động lực học ngược:

    • Xác định quỹ đạo thực tế của khâu thao tác cuối $E(x_E, y_E)$ chịu ảnh hưởng của biến dạng đàn hồi: $$\begin{cases} x_E(t) = x_E^R(t) + \Delta x_E(w) \ y_E(t) = y_E^R(t) + \Delta y_E(w) \end{cases}$$
    • Đề xuất thuật toán xác định gần đúng mô-men/lực dẫn động bổ sung $\Delta \boldsymbol{\tau}(t)$ kết hợp với mô-men chuyển động rắn $\boldsymbol{\tau}_R(t)$ để tạo thành mô-men phát động $\boldsymbol{\tau}_e(t) = \boldsymbol{\tau}_R(t) + \Delta \boldsymbol{\tau}(t)$, giúp điểm thao tác cuối bám sát quỹ đạo thiết kế với sai số triệt tiêu.

Chương 4: Kết luận chung, kiến nghị và phương hướng phát triển

Chương này tổng hợp các kết quả lý thuyết và mô phỏng đã đạt được trong toàn bộ luận án, đồng thời chỉ rõ các hạn chế và đề xuất các phương hướng nghiên cứu tiếp theo của đề tài.

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

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

  1. Thiết lập mô hình động lực học tường minh bằng phương pháp hệ quy chiếu đồng hành: Luận án đã ứng dụng thành công phương pháp hệ quy chiếu đồng hành kết hợp phương trình Lagrange loại 2 để xây dựng hệ phương trình vi phân chuyển động phi tuyến hoàn chỉnh cho các dạng tay máy phẳng có khâu đàn hồi (1 khâu quay, 2 khâu T-R, 2 khâu R-R). Mô hình này đảm bảo tính nhất quán cao, khi triệt tiêu biến dạng đàn hồi sẽ suy biến chính xác về mô hình tay máy rắn cổ điển.
  2. Phương pháp giải tích tuyến tính hóa Taylor tổng quát: Xây dựng thành công công thức khai triển Taylor của hàm ma trận và hàm véc tơ theo biến véc tơ. Đây là đóng góp có ý nghĩa phương pháp luận, giúp chuyển đổi hệ phương trình vi phân chuyển động phi tuyến phức tạp sang hệ phương trình vi phân tuyến tính hệ số tuần hoàn quanh quỹ đạo cơ bản một cách chặt chẽ và có hệ thống.
  3. Quy trình tối ưu hóa ổn định động lực học kết hợp Floquet - Taguchi: Đề xuất một giải pháp kết hợp lý thuyết ổn định Floquet và phương pháp quy hoạch thực nghiệm Taguchi để tối ưu hóa ma trận hệ số điều khiển ($K_P, K_D$). Thuật toán cho phép tìm được bộ tham số điều khiển tối ưu giúp hệ thống luôn có môđun nhân tử Floquet cực đại nhỏ hơn 1 ($|\rho|_{\max} < 1$), đảm bảo tay máy vận hành ổn định ngoài vùng cộng hưởng tham số.

Đóng góp mới về mặt thực tiễn và ứng dụng

  1. Phương pháp giải gần đúng bài toán động lực học ngược: Đề xuất giải pháp giải gần đúng bài toán động lực học ngược cho tay máy có khâu đàn hồi chuyển động tuần hoàn. Phương pháp cho phép tính toán chính xác độ võng đàn hồi và quy luật mô-men phát động bù trừ, giúp cơ cấu bám sát quỹ đạo làm việc mong muốn với độ chính xác cao.

Kiến nghị và đề xuất từ luận án

  • Ứng dụng quy trình tính toán ổn định và bù trừ mô-men dẫn động vào quá trình thiết kế, chế tạo các cơ cấu robot công nghiệp tốc độ cao và robot nhẹ trong thực tế.
  • Tích hợp công cụ giải tích tuyến tính hóa vào các phần mềm mô phỏng động lực học để hỗ trợ kỹ sư tự động hóa việc đánh giá vùng mất ổn định tham số.

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

Hạn chế của luận án

  • Phạm vi không gian chuyển động: Luận án chỉ giới hạn khảo sát các mô hình tay máy chuyển động trong mặt phẳng hai chiều (2D).
  • Mô hình biến dạng đàn hồi: Chỉ xét biến dạng uốn phẳng theo lý thuyết dầm Euler-Bernoulli; chưa tính đến biến dạng xoắn, kéo - nén dọc trục và các hiệu ứng bậc cao (quán tính quay, biến dạng trượt cắt theo mô hình dầm Timoshenko).
  • Phương pháp kiểm chứng: Toàn bộ các kết quả đánh giá ổn định động lực học và điều khiển bám quỹ đạo mới dừng lại ở mức độ tính toán giải tích và mô phỏng số thông qua phần mềm Maple và Matlab, chưa được kiểm chứng thực nghiệm trên mô hình vật lý chế tạo thực tế.

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

  1. Mở rộng bài toán động lực học và điều khiển sang không gian ba chiều (3D) cho các robot nhiều bậc tự do có khâu đàn hồi.
  2. Xây dựng mô hình tính toán kết hợp đồng thời các dạng biến dạng: uốn trong không gian, xoắn và dãn dọc trục.
  3. Ứng dụng các phương pháp phần tử hữu hạn nâng cao kết hợp các thuật toán điều khiển phi tuyến hiện đại (điều khiển thích nghi, điều khiển trượt, điều khiển bền vững).
  4. Thiết kế và chế tạo mô hình thực nghiệm tay máy có khâu đàn hồi để đo đạc dao động thực tế và kiểm chứng tính chính xác của thuật toán điều khiển đề xuất.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị cao trong lĩnh vực Cơ học ứng dụng, Động lực học hệ nhiều vật và Kỹ thuật điều khiển robot:

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên:
    • Là tài liệu tham khảo chi tiết về kỹ thuật thiết lập phương trình Lagrange loại 2 cho hệ nhiều vật đàn hồi bằng phương pháp hệ quy chiếu đồng hành.
    • Cung cấp công cụ toán học giải tích rõ ràng về khai triển vi phân Taylor cho hàm ma trận/véc tơ theo biến véc tơ.
    • Minh họa phương pháp ứng dụng lý thuyết Floquet trong khảo sát ổn định dao động tham số của các hệ kỹ thuật thực tế.
  • Đối với kỹ sư thiết kế máy và chuyên gia robot:
    • Cung cấp quy trình tính toán, lựa chọn thông số điều khiển ($K_P, K_D$) tối ưu thông qua phương pháp Taguchi nhằm triệt tiêu rung động và tránh vùng phá hủy do cộng hưởng tham số.
    • Cung cấp thuật toán tính toán mô-men bù trừ cho động cơ khâu dẫn để nâng cao độ chính xác vị trí khâu công tác cuối khi robot vận hành ở tốc độ cao.

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

1. Luận án nghiên cứu những mô hình tay máy cụ thể nào và sử dụng lý thuyết dầm nào để mô hình hóa khâu đàn hồi?

Luận án nghiên cứu 3 mô hình tay máy phẳng dạng chuỗi:

  1. Tay máy 1 khâu đàn hồi chuyển động quay có mang tải trọng tại đầu mút.
  2. Tay máy 2 khâu T-R (khâu 1 chuyển động tịnh tiến coi là rắn, khâu 2 chuyển động quay coi là đàn hồi).
  3. Tay máy 2 khâu R-R (khâu 1 quay coi là rắn, khâu 2 quay coi là đàn hồi có mang tải trọng đầu mút).

Biến dạng của khâu đàn hồi được mô hình hóa theo lý thuyết dầm Euler-Bernoulli chịu uốn ngang (bỏ qua biến dạng cắt và quán tính quay), rời rạc hóa chuyển vị bằng chuỗi hàm riêng Ritz-Galerkin với điều kiện biên ngàm - tự do.

2. Khái niệm "chuyển động cơ bản" dùng trong tuyến tính hóa phương trình vi phân chuyển động được định nghĩa như thế nào?

Chuyển động cơ bản $\mathbf{q}^0(t) = [\mathbf{q}_a^0(t)^T, \mathbf{0}^T]^T$ được định nghĩa là quy luật chuyển động của tay máy khi giả thiết tất cả các khâu của cơ hệ đều là vật rắn tuyệt đối (biến dạng đàn hồi $\mathbf{q}_e = \mathbf{0}$). Hệ phương trình vi phân phi tuyến sau đó được tuyến tính hóa giải tích xung quanh chuyển động cơ bản này bằng khai triển Taylor ma trận/véc tơ, đưa về hệ phương trình vi phân tuyến tính có hệ số biến thiên tuần hoàn theo thời gian.

3. Lý thuyết Floquet và phương pháp Taguchi được kết hợp như thế nào để điều khiển ổn định động lực cho cơ hệ?

Hệ phương trình vi phân tuyến tính hóa sau khi áp dụng luật điều khiển phản hồi vị trí và vận tốc ($K_P, K_D$) được đưa về hệ phương trình trạng thái chu kỳ $\dot{\mathbf{x}} = \mathbf{A}(t)\mathbf{x}$. Lý thuyết Floquet được sử dụng để xác định ma trận chuyển tiếp trạng thái chu kỳ $\mathbf{\Phi}(T)$ và các nhân tử Floquet $\rho_j$. Điều kiện ổn định đòi hỏi mọi môđun nhân tử Floquet phải thỏa mãn $|\rho_j| < 1$. Phương pháp Taguchi với bảng trực giao $L_9$ và tỷ số tín hiệu trên nhiễu $\text{SNR}$ (tiêu chí "càng nhỏ càng tốt") được áp dụng để tối ưu hóa việc chọn các mức tham số $K_P, K_D$ nhằm cực tiểu hóa $|\rho|_{\max}$, bảo đảm hệ thống luôn ổn định tiệm cận và tránh xa vùng cộng hưởng tham số.

4. Những dải vận tốc góc nào đã được khảo sát mô phỏng đối với ba mô hình tay máy trong luận án?

Trong các khảo sát số và mô phỏng nghiệm chuyển tiếp ở Chương 3, các dải vận tốc góc $\omega$ được thực hiện cụ thể:

  • Tay máy một khâu đàn hồi: Khảo sát tại các vận tốc góc $\omega = 2, 4, 6, 8, 10\text{ rad/s}$.
  • Tay máy hai khâu T-R đàn hồi: Khảo sát tại các vận tốc góc $\omega = 1, 2, 4, 8, 9\text{ rad/s}$.
  • Tay máy hai khâu R-R đàn hồi: Khảo sát tại các vận tốc góc $\omega = 4, 6, 8, 10, 12\text{ rad/s}$.

5. Phương pháp gần đúng giải bài toán động lực học ngược trong luận án giải quyết mục tiêu kỹ thuật gì?

Bài toán động lực học ngược giải quyết việc xác định quy luật lực/mô-men phát động $\boldsymbol{\tau}_e(t)$ tại các khớp dẫn động sao cho khâu thao tác cuối $E$ chuyển động chính xác theo quỹ đạo định trước. Thuật toán của luận án tính toán quy luật mô-men dựa trên sự kết hợp giữa mô-men chuyển động rắn $\boldsymbol{\tau}_R(t)$ và thành phần mô-men bù trừ $\Delta \boldsymbol{\tau}(t)$ suy ra từ nghiệm dao động đàn hồi tuần hoàn, từ đó triệt tiêu sai số vị trí do độ võng uốn của khâu đàn hồi gây ra.

Kết luận

Luận án tiến sĩ của tác giả Đinh Công Đạt đã giải quyết một cách có hệ thống bài toán động lực học và điều khiển cho lớp tay máy phẳng có khâu đàn hồi chuyển động tuần hoàn. Bằng việc kết hợp phương pháp hệ quy chiếu đồng hành, công cụ khai triển Taylor ma trận giải tích, lý thuyết ổn định Floquet và phương pháp tối ưu hóa Taguchi, tác giả đã xây dựng hoàn chỉnh mô hình toán học, thuật toán điều khiển ổn định ngoài vùng cộng hưởng tham số và phương pháp giải gần đúng bài toán động lực học ngược. Các kết quả nghiên cứu đóng góp nền tảng lý luận vững chắc và mang lại giải pháp tính toán hữu ích cho việc thiết kế, điều khiển các thế hệ robot thanh mảnh, tốc độ cao và độ chính xác cao trong kỹ thuật hiện đại.