Tổng quan về luận án

Hệ thống truyền động khí nén (TĐKN) giữ vai trò then chốt trong nền sản xuất công nghiệp hiện đại nhờ ưu thế vượt trội về mật độ công suất, chi phí đầu tư thấp, khả năng vận hành ở vận tốc cao (pít-tông đạt trên 15 m/s, động cơ quay đạt 100.000 vòng/phút), độ sạch môi trường và tính an toàn cháy nổ cao trong các môi trường khắc nghiệt chịu phóng xạ, hóa chất hay điện từ trường. Tuy nhiên, rào cản kỹ thuật lớn nhất trong việc mở rộng phạm vi ứng dụng của hệ thống TĐKN điều khiển tỉ lệ vào các hệ thống đòi hỏi độ chính xác cao (như máy công cụ CNC, robot phẫu thuật y tế, thiết bị đo lường chuẩn) nằm ở "đặc tính động lực học phức tạp và phi tuyến bậc cao do tính nén được của không khí, đặc tính phi tuyến của van khí và lực ma sát trong các cơ cấu chấp hành (CCCH) khí nén". Hiện tượng ma sát phi tuyến tiếp xúc giữa bề mặt phớt làm kín (gioăng cao su, teflon) và thân xi lanh dẫn đến chu kỳ giới hạn (limit cycles), hiện tượng dính - trượt (stick-slip), và độ trễ phi tuyến, gây suy giảm nghiêm trọng độ chính xác định vị tại các điểm dừng trung gian.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ việc: dù nhiều mô hình ma sát đã được đề xuất cho các cơ cấu cơ khí và thủy lực, việc đánh giá định lượng và so sánh có hệ thống hiệu quả mô phỏng của các mô hình ma sát động học hiện đại đối với CCCH khí nén vẫn chưa được thực hiện toàn diện. Đồng thời, các giải pháp điều khiển trượt (sliding mode control) truyền thống thường xem ma sát như nhiễu bất định có chặn trên mà chưa tích hợp cấu trúc bù ma sát động học thích ứng theo mô hình.

Luận án tiến sĩ chuyên ngành Kỹ thuật Cơ học / Tự động hóa và Cơ điện tử giải quyết hai câu hỏi và giả thuyết khoa học chính:

  1. Câu hỏi nghiên cứu 1 (RQ1): Mô hình ma sát nào trong số các mô hình kinh điển và hiện đại (mô hình trạng thái ổn định, mô hình LuGre, mô hình LuGre cải tiến - RLuGre) phản ánh chính xác nhất các hành vi động học ma sát phi tuyến trong cả giai đoạn trước khi trượt (pre-sliding) và giai đoạn trượt (sliding) của xy lanh khí nén?
    • Giả thuyết 1 (H1): Mô hình LuGre cải tiến (RLuGre) tích hợp động lực học độ dày màng bôi trơn $h$ và hàm trễ phi tuyến $F(z)$ với bộ nhớ không cục bộ sẽ tái hiện chính xác nhất động lực học hệ thống thực tế so với mô hình LuGre chuẩn và mô hình ma sát tĩnh.
  2. Câu hỏi nghiên cứu 2 (RQ2): Cấu trúc điều khiển phi tuyến nào triệt tiêu được hiện tượng trượt dính và đảm bảo độ chính xác bám quỹ đạo vị trí tại các điểm dừng trung gian dưới sự biến thiên tải trọng và áp suất?
    • Giả thuyết 2 (H2): Bộ điều khiển trượt đa mặt trượt (Multi-Surface Sliding Mode Controller - MSSMC) kết hợp khâu bù ma sát thích nghi dựa trên quan sát trạng thái LuGre sẽ triệt tiêu sai số xác lập, hạ sai số vị trí đỉnh xuống dưới 4% và triệt tiêu dao động chattering.

Phạm vi thực nghiệm được giới hạn nghiêm ngặt trên hệ thống TĐKN tỉ lệ sử dụng xy lanh tác động hai phía SMC CM2L25-300 (đường kính nòng $D = 25\text{ mm}$, đường kính cần $d = 10\text{ mm}$, hành trình $L = 300\text{ mm}$), áp suất nguồn khí nén làm việc biến thiên đến $8 \times 10^5\text{ N/m}^2$ ($8\text{ bar}$), và tải trọng ngoài biến thiên $M = 0{,}5 \div 5{,}0\text{ kg}$, được thu thập dữ liệu tại Phòng thí nghiệm 207-C8, Trường Đại học Bách Khoa Hà Nội trong giai đoạn 10/2015 đến 04/2021.

+-----------------------------------------------------------------------------------+
|                           HỆ THỐNG TRUYỀN ĐỘNG KHÍ NÉN TỈ LỆ                      |
+-----------------------------------------+-----------------------------------------+
                                          |
        +---------------------------------+---------------------------------+
        |                                                                   |
        v                                                                   v
+---------------------------------------+   +---------------------------------------+
|  MỤC TIÊU 1: MÔ HÌNH HÓA ĐỘNG LỰC HỌC |   |     MỤC TIÊU 2: ĐIỀU KHIỂN PHI TUYẾN   |
+---------------------------------------+   +---------------------------------------+
| - Mô hình Trạng thái ổn định          |   | - Bộ điều khiển Đa mặt trượt (MSSMC)  |
| - Mô hình LuGre động học chuẩn        |   | - Tích hợp bộ bù ma sát LuGre động học|
| - Mô hình LuGre cải tiến (RLuGre)     |   | - Chứng minh ổn định theo Lyapunov    |
+-------------------+-------------------+   +-------------------+-------------------+
                    |                                           |
                    +---------------------+---------------------+
                                          v
+-----------------------------------------------------------------------------------+
|   ĐÁNH GIÁ THỰC NGHIỆM ĐỐI CHỨNG (HUST LAB 207-C8: 1.16 ms, Filter 32 Hz)         |
|   -> Nâng cao độ chính xác mô phỏng và hạ sai số định vị trung gian đỉnh < 4%     |
+-----------------------------------------------------------------------------------+

Literature Review và Positioning

Nghiên cứu về cơ chế ma sát trong các hệ cơ điện tử trải qua quá trình tiến hóa từ các mô hình tĩnh cổ điển đến các mô hình động học vi cấu trúc phức tạp. Nền tảng ma sát học khởi nguồn từ Leonardo da Vinci (1493), Guillaume Amontons (1699) và Charles Augustin de Coulomb (1785) với định luật ma sát khô tỉ lệ tuyến tính với lực pháp tuyến $F_C = \mu N$. Tiếp đó, Osborn Reynolds (1886) đưa vào thành phần ma sát nhớt $\sigma_2 v$, và Arthur Morin (1833) bổ sung khái niệm lực ma sát tĩnh cực đại $F_s > F_C$. Richard Stribeck (1902) phát hiện hiệu ứng phi tuyến vận tốc thấp: lực ma sát suy giảm phi tuyến khi vận tốc trượt tăng nhẹ trước khi bước vào chế độ bôi trơn thủy động hoàn toàn. Armstrong-Hélouvry (1991) đã hệ thống hóa mô hình trạng thái ổn định 5 thông số kết hợp Stribeck, Coulomb và nhớt.

Tuy nhiên, các mô hình tĩnh thất bại hoàn toàn khi mô tả các hiệu ứng động học: biến dạng đàn hồi trước khi trượt (pre-sliding displacement), hiện tượng trễ lực ma sát (friction force hysteresis), sự thay đổi của lực đứt-gãy (break-away force), và dao động trượt-dính (stick-slip). P.R. Dahl (1968) tiên phong mô hình hóa chuyển vị đàn hồi vi mô tựa lò xo giữa các điểm tiếp xúc gồ ghề (asperities) cho ứng dụng hàng không vũ trụ. Canudas de Wit et al. (1995) tạo bước ngoặt với mô hình LuGre, tích hợp hàm Stribeck $g(v)$ vào động lực học độ võng vi mô trung bình $z$, cho phép mô tả cả giai đoạn trước khi trượt và trượt liên tục. Nhằm khắc phục hiện tượng trôi dạt (drift) và mô tả chính xác vòng trễ, Swevers et al. (2000) phát triển mô hình Leuven, Dupont et al. (2002) đề xuất mô hình Elastoplastic, và Al-Bender et al. (2005) xây dựng mô hình Maxwell-Slip tổng quát (GMS). Mặc dù có độ chính xác hình thức cao, các mô hình Leuven và GMS chứa quá nhiều thông số phức tạp, gây khó khăn lớn cho việc nhận dạng thời gian thực trên vi xử lý công nghiệp.

Tiến trình phát triển mô hình ma sát cơ cấu chấp hành:
[Coulomb (1785)] -> [Stribeck (1902)] -> [Dahl (1968)] -> [LuGre (1995)] -> [RLuGre (Tran et al., Yanada 2008-2013)]
                                                                                |
                                                  +-----------------------------+-----------------------------+
                                                  v                                                           v
                                    [Động lực học màng bôi trơn h]                             [Hàm trễ F(z) Non-local memory]

Đối với CCCH khí nén, Schroeder et al. và Belforte et al. xác định thực nghiệm lực ma sát phụ thuộc vào kích thước nòng, vận tốc và áp suất buồng khí. Duong (2013) nghiên cứu ảnh hưởng của môi trường nhiệt ẩm Việt Nam, chứng minh: "đặc tính ma sát thay đổi theo tốc độ dịch chuyển có dạng đường cong Stribeck với giá trị nhỏ nhất của vận tốc trong khoảng 25 ÷ 35 mm/s. Đồng thời lực ma sát giảm khoảng 10% ÷ 17% khi nhiệt độ tăng từ 15 ÷ 50 °C, lực ma sát giảm khoảng 14% ÷ 24% khi độ ẩm tương đối tăng từ 51% ÷ 99%". Nouri (2000) phát hiện: "đặc tính trễ của ma sát trong giai đoạn trước khi trượt được trình bày bằng đặc tính 'nhớ không cục bộ'... giá trị tương lai của hàm ma sát ở một thời điểm bất kỳ t (t > t0) phụ thuộc không chỉ vào giá trị hiện tại ở thời điểm t0 và đối số của nó mà còn giá trị cực đại trước đó của hàm ma sát". Tiếp đó, Yanada & Sekikawa (2008) cùng Tran et al. (2012, 2013) phát triển mô hình LuGre cải tiến (RLuGre) bổ sung độ dày màng bôi trơn $h$ và khâu vi phân dẫn pha bậc một $\sigma_2(v + T \dot{v})$.

Về điều khiển, các công trình của Wang et al., Richer & Hurmuzlu, và Takosoglu et al. đã mô hình hóa phi tuyến dòng khí nhưng bỏ qua ma sát động. Tsai & Huang (2008) đề xuất bộ điều khiển trượt đa mặt trượt (MSSMC) trên các mặt trượt vị trí, vận tốc và lực, đạt sai số bám đỉnh 4% nhưng chỉ coi ma sát là nhiễu có chặn. Các công trình bù ma sát của Situm et al. (Fuzzy-PID), Gao & Feng (2010), hay Drakunov et al. (Sliding Mode bù nhớt) chỉ dựa trên mô hình Coulomb hoặc tĩnh.

Luận án này định vị chính xác vào khoảng trống tri thức: thiết lập cầu nối tích hợp giữa mô hình ma sát động học RLuGre hoàn chỉnh vào hệ thống mô phỏng tổng thể 7 phương trình vi phân phi tuyến bậc cao và hiện thực hóa bộ điều khiển MSSMC tích hợp bộ quan sát bù ma sát LuGre trên hệ thực nghiệm thời gian thực.

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

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

Luận án đóng góp vào lý thuyết điều khiển tự động và cơ học chất khí - ma sát học thông qua ba luận điểm đột phá:

  • Mở rộng lý thuyết ma sát học tiếp xúc trong môi trường bôi trơn mỡ: Luận án chứng minh sự tương thích của lý thuyết bôi trơn màng mỡ (Li et al., Sugimura et al.) với động lực học biến thiên của màng bôi trơn không thứ nguyên $h$, giải thích cơ chế suy giảm lực đứt-gãy (break-away force) sau các khoảng thời gian dừng ngắn thông qua số hạng $F_s(1-h)$.
  • Khung toán học động học RLuGre hoàn chỉnh: Thiết lập hệ phương trình vi sai tích hợp khâu trễ cục bộ phi tuyến $F(z)$, mở rộng phương trình vi phân trạng thái võng sợi tiếp xúc:

$$\frac{dz}{dt} = \frac{F(z)}{g(v,h)} v$$

$$g(v,h) = F_c + (1-h)(F_s - F_c)e^{-(v/v_s)^n}$$

$$F_r = F(z) + \sigma_1 \frac{dz}{dt} + \sigma_2\left(v + T \frac{dv}{dt}\right)$$

trong đó động lực học màng bôi trơn được điều khiển bởi hằng số thời gian $\tau_h$ phân tách theo giai đoạn: tăng tốc ($h \le h_{ss}, \tau_h = \tau_{hp}$), giảm tốc ($h > h_{ss}, \tau_h = \tau_{hn}$), và dừng ($\tau_h = \tau_{h0}$).

  • Mô hình toán học kết nối tầng 5 biến trạng thái: Xây dựng mô hình không gian trạng thái tổng quát biểu diễn quan hệ phi tuyến giữa tín hiệu điện áp điều khiển $u_1, u_2$ của hai van tỉ lệ lưu lượng độc lập với vị trí pít-tông $x$, vận tốc $v$, áp suất hai buồng $p_1, p_2$, và trạng thái võng ma sát $z$.

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

Khung phân tích của luận án tích hợp liên ngành giữa Cơ học thủy khí nén (Định luật bảo toàn khối lượng, Định luật bảo toàn năng lượng Bernoulli, Định luật I Nhiệt động lực học cho quá trình đoạn nhiệt khí lý tưởng), Ma sát học vi mô (Asperity contact micro-mechanics), và Lý thuyết điều khiển phi tuyến hiện đại (Lyapunov stability theory).

+---------------------------------------------------------------------------------------------------+
|                                KHUNG PHÂN TÍCH TOÁN HỌC LIÊN NGÀNH                                |
+---------------------------------------------------------------------------------------------------+
|                                                                                                   |
|  [ĐIỀU KHIỂN ĐA MẶT TRƯỢT (MSSMC)]           [ĐỘNG LỰC HỌC CHẤT KHÍ]        [MA SÁT HỌC VI MÔ]    |
|  - Mặt trượt 1: S1 = e = x - x_d             - Dòng nén đoạn nhiệt qua van: - Asperities võng z   |
|  - Mặt trượt 2: S2 = v - alpha1                qm1(u1, p1), qm2(u2, p2)     - Màng mỡ bôi trơn h  |
|  - Mặt trượt 3: S3 = (p1*A1-p2*A2) - alpha2  - Biến thiên áp suất khoang:   - Trễ non-local F(z)  |
|                                                dp1/dt, dp2/dt               - Vi phân pha dẫn T   |
|                                                                                                   |
+-------------------------------------------------+-------------------------------------------------+
                                                  |
                                                  v
+---------------------------------------------------------------------------------------------------+
|              LUẬT ĐIỀU KHIỂN ĐIỆN ÁP TỔNG HỢP u1, u2 TÍCH HỢP BÙ MA SÁT DỰA TRÊN LuGre           |
|                               V = 1/2*(S1^2 + S2^2 + S3^2) <= 0 (Lyapunov)                        |
+---------------------------------------------------------------------------------------------------+

Điều kiện biên xác định: Giới hạn nhiệt độ phòng ổn định $T_0 = 293\text{ K}$, bỏ qua rò rỉ khí qua khe hở pít-tông, cấu trúc xy lanh giả định rắn tuyệt đối không biến dạng đàn hồi thể tích, và tần số dao động kích thích không vượt quá ngưỡng trễ khí nén của van ($< 2\text{ Hz}$).

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ triết lý thực chứng thực nghiệm (Positivist / Empirical Paradigm) kết hợp mô hình hóa giải tích diễn dịch (deductive analytical modeling). Thiết kế đa mức (multi-level design) bao gồm:

  1. Mức 1 (Component Level): Nhận dạng thông số tĩnh và động của từng phần tử riêng biệt (van tỉ lệ SMC VEF3121, cảm biến áp suất PSE540, cảm biến vị trí Novotechnik LWH0300).
  2. Mức 2 (Subsystem Friction Level): Khảo sát độc lập hành vi ma sát của xy lanh SMC CM2L25-300 dưới các điều kiện kích thích vận tốc và áp suất khác nhau.
  3. Mức 3 (System Control Level): Tích hợp toàn hệ thống kín, thực thi thuật toán điều khiển và đánh giá chất lượng bám quỹ đạo thời gian thực.
        +-------------------------------------------------------------+
        |                 HỆ THỐNG THỰC NGHIỆM ĐỒNG BỘ                |
        +------------------------------+------------------------------+
                                       |
       +-------------------------------+-------------------------------+
       |                               |                               |
       v                               v                               v
+---------------+             +------------------+             +---------------+
|  KHÍ ĐỘNG LỰC |             |  CƠ KHÍ CHẤP HÀNH|             | ĐO LƯỜNG - XỬ LÝ|
| - Máy nén khí |             | - Xy lanh khí nén|             | - Cảm biến LWH|
|   PUMA PK1090A|             |   CM2L25-300     |             |   (0-300 mm)  |
| - 02 Van tỉ lệ|             | - Tải M = 0.5-5kg|             | - 02 Cảm biến |
|   VEF3121-1-02|             | - Ray dẫn hướng  |             |   PSE540      |
| - 02 Driver   |             |   chính xác      |             | - DAQ USB4711 |
|   VEA250      |             |                  |             | - Visual C++  |
+---------------+             +------------------+             +---------------+

Quy trình nghiên cứu rigorous

Hệ thống thực nghiệm được thiết kế đối xứng kép sử dụng hai van tỉ lệ lưu lượng điện - khí nén 3 cửa (SMC VEF3121-1-02, diện tích cửa thông $12\text{ mm}^2$, dòng điều khiển $0 \div 1\text{ A}$ tương ứng điện áp $0 \div 5\text{ VDC}$) điều khiển độc lập hai khoang xy lanh, khắc phục nhược điểm mất đối xứng lưu lượng của van 5/3 truyền thống.

Giao thức đo lường và nhận dạng thông số ma sát:

  • Cảm biến vị trí tuyến tính: Novotechnik LWH0300 (dải đo $0 \div 300\text{ mm}$, độ tuyến tính độc lập $\pm 0{,}05%\text{ F.S.}$, độ lặp lại $0{,}01\text{ mm}$).
  • Cảm biến áp suất: Hai bộ SMC PSE540 (dải đo $0 \div 1\text{ MPa}$, độ chính xác $\pm 2%\text{ F.S.}$).
  • Thu thập dữ liệu: Card DAQ Advantech USB-4711A (16-bit A/D, 12-bit D/A, tốc độ lấy mẫu $150\text{ kS/s}$), chu kỳ trích mẫu thời gian thực $\Delta t = 1{,}16\text{ ms}$ trên nền tảng Microsoft Visual C++.
  • Lọc tín hiệu vi phân: Vận tốc $v$ và gia tốc $a$ được tính toán qua đạo hàm số bậc một và bậc hai từ vị trí $x$, kết hợp bộ lọc số phi nhân quả thông thấp (acausal low-pass filter) tần số cắt $32\text{ Hz}$ nhằm triệt tiêu hoàn toàn nhiễu lượng tử hóa mà không gây trễ pha.

Lực ma sát thực nghiệm $F_r$ được trích xuất động lực học theo phương trình cân bằng D'Alembert:

$$F_r = p_1 A_1 - p_2 A_2 - M a$$

trong đó $A_1 = \frac{\pi D^2}{4} = 4{,}909 \times 10^{-4}\text{ m}^2$, $A_2 = \frac{\pi (D^2 - d^2)}{4} = 4{,}123 \times 10^{-4}\text{ m}^2$.

Data và phân tích

Quy trình nhận dạng thông số ma sát RLuGre được tiến hành qua hai bước phân lập:

  1. Thông số tĩnh ($F_s, F_c, v_s, n, \sigma_2, v_b$): Thực hiện chạy pít-tông ở các vận tốc không đổi từ $0{,}001\text{ m/s}$ đến $0{,}5\text{ m/s}$. Áp dụng phương pháp tối ưu hóa phi tuyến bình phương tối thiểu (Nonlinear Least Squares) để khớp đường cong thực nghiệm với phương trình Stribeck. Vận tốc chuyển tiếp $v_b$ được xác định chính xác bằng $1{,}5$ lần vận tốc tại điểm cực tiểu của lực ma sát Stribeck.
  2. Thông số động ($\sigma_0, \sigma_1, \tau_{hp}, \tau_{hn}, \tau_{h0}, T, F(z)$):
    • Hằng số thời gian màng bôi trơn $\tau_{hp} = 0{,}035\text{ s}$ và $\tau_{hn} = 0{,}058\text{ s}$ được xác định từ phản ứng suy giảm hàm mũ của lực ma sát khi thay đổi vận tốc dạng bước nhảy.
    • Hằng số thời gian pha dẫn $T$ được nhận dạng từ quan hệ tuyến tính giữa sai số dư $\Delta F$ và gia tốc $dv/dt$: đạt $T = 0{,}33\text{ s}$ cho hành trình thuận và $T = 0{,}07\text{ s}$ cho hành trình nghịch.
    • Độ cứng tiếp xúc $\sigma_0$ và độ cản $\sigma_1$ được tinh chỉnh và kiểm chứng thông qua so sánh độ vọt lố của lực đứt-gãy thực nghiệm và mô phỏng trên MATLAB/Simulink.

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 khoa học mang tính bước ngoặt:

  1. Sự vượt trội của mô hình RLuGre trong việc tái hiện trễ động học: So sánh định lượng giữa mô hình tĩnh, LuGre chuẩn và RLuGre chỉ ra rằng mô hình RLuGre là mô hình duy nhất mô phỏng chính xác vòng trễ ma sát - vận tốc ở dải vận tốc cao ($> 0{,}1\text{ m/s}$) và hiện tượng suy giảm lực đứt-gãy theo chu kỳ vận tốc đảo chiều hình sin ($0{,}5 \div 2\text{ Hz}$). Mô hình LuGre chuẩn đánh giá thấp độ trễ ở vận tốc lớn do thiếu số hạng vi phân bậc một $\sigma_2 T \dot{v}$.
  2. Khả năng mô phỏng vòng trễ con trong chế độ trước khi trượt: Trong giai đoạn chuẩn bị trượt ($|x| < 15\text{ }\mu\text{m}$), khi vận tốc đảo chiều liên tiếp nhiều lần, lực ma sát thể hiện đặc tính lò xo phi tuyến có bộ nhớ không cục bộ. Hàm phân đoạn $f_i(z) = f_i(z_i) + c_i(1 - e^{-k_i(z - z_i)})$ trong RLuGre đã mô phỏng khớp hoàn toàn các vòng lặp trễ bên trong (internal hysteresis loops), điều mà mô hình LuGre truyền thống hoàn toàn bất lực.
  3. Triệt tiêu hiện tượng dính - trượt bằng bộ điều khiển MSSMC bù ma sát LuGre: Khi thử nghiệm đáp ứng vị trí với tín hiệu đặt dạng bước và hình sin tại các tần số $0{,}25\text{ Hz}, 0{,}5\text{ Hz}$ và $1{,}0\text{ Hz}$, bộ điều khiển MSSMC kết hợp bù ma sát LuGre triệt tiêu hoàn toàn dao động dính - trượt ở vùng vận tốc gần 0, không xuất hiện hiện tượng dao động tham số (chattering) trên tín hiệu điện áp van $u_1, u_2$.
  4. Độ chính xác bám vị trí vượt bậc: Sai số vị trí đỉnh lớn nhất của hệ thống được kiểm soát nghiêm ngặt dưới mức $4%$ đối với quỹ đạo hình sin tần số cao và đạt độ chính xác định vị dừng trung gian với sai số xác lập tiến về $0\text{ mm}$ (dưới ngưỡng sai số đo của cảm biến $0{,}05\text{ mm}$), vượt trội so với các thuật toán PID tuyến tính và SMC thông thường không có khâu bù ma sát động.
BẢNG SO SÁNH NĂNG LỰC CỦA CÁC MÔ HÌNH MA SÁT TRONG MÔ PHỎNG CCCH KHÍ NÉN:
+-------------------------------+-------------------+-------------------+-------------------+
| Tiêu chí khảo sát             | Mô hình Tĩnh      | Mô hình LuGre     | Mô hình RLuGre    |
+-------------------------------+-------------------+-------------------+-------------------+
| Hiệu ứng Stribeck tĩnh        | Xuất sắc          | Xuất sắc          | Xuất sắc          |
| Hiện tượng trượt - dính       | Không thể         | Khá               | Xuất sắc          |
| Độ trễ ở vận tốc cao          | Không thể         | Kém (bị hẹp)      | Xuất sắc (với T)  |
| Trễ bộ nhớ không cục bộ       | Không thể         | Không thể         | Xuất sắc (với F(z)|
| Biến thiên màng bôi trơn h    | Không thể         | Không thể         | Xuất sắc          |
+-------------------------------+-------------------+-------------------+-------------------+

Implications đa chiều

  • Về lý thuyết điều khiển: Khẳng định tính tất yếu của việc tích hợp mô hình vi cấu trúc ma sát vào các bộ điều khiển phi tuyến hiện đại, đặt nền móng cho các kiến trúc điều khiển thích nghi đa tầng cho hệ thống truyền động mềm (soft actuators).
  • Về phương pháp luận: Chuẩn hóa quy trình thực nghiệm nhận dạng phân tách hai pha (tĩnh - động) cho CCCH khí nén bằng cách sử dụng cấu hình van tỉ lệ kép, loại bỏ hoàn toàn ảnh hưởng gây nhiễu do mất đối xứng áp suất.
  • Về kỹ thuật công nghiệp: Cung cấp giải pháp phần mềm cho phép biến các hệ truyền động khí nén chi phí thấp thành các cơ cấu định vị chính xác cao tương đương hệ thống servo thủy lực hoặc động cơ bước, giảm thiểu điện năng tiêu thụ do triệt tiêu được hiện tượng xả khí thừa khi dao động quanh điểm dừng.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Giới hạn về mô hình nhiệt ẩm: Mặc dù tham chiếu kết quả của Duong (2013), mô hình toán RLuGre trong luận án chưa tích hợp trực tiếp phương trình biến thiên của nhiệt độ và độ ẩm môi trường theo thời gian thực ($T_a, RH$) vào hàm $g(v,h)$.
  2. Giới hạn về độ mòn phớt làm kín: Các thông số ma sát được nhận dạng trong trạng thái phớt còn mới; ảnh hưởng của quá trình thoái hóa cơ học và mài mòn vật liệu sau hàng triệu chu kỳ làm việc chưa được khảo sát dài hạn.
  3. Đơn giản hóa nhiệt động học buồng khí: Quá trình nhiệt động trong xy lanh được coi là đoạn nhiệt lý tưởng ($k = 1{,}4$), bỏ qua quá trình truyền nhiệt qua thành vỏ xy lanh bằng nhôm/thép ra môi trường khi vận hành ở tần số cao liên tục.
  4. Khối lượng tính toán của RLuGre trong điều khiển thời gian thực: Trong Chương 4, do giới hạn năng lực xử lý tức thời của vi xử lý điều khiển tại chu kỳ $1{,}16\text{ ms}$, bộ bù ma sát thời gian thực mới tích hợp mô hình LuGre động học chuẩn thay vì toàn bộ mô hình RLuGre đầy đủ (vốn được dùng trong mô phỏng giải tích Chương 3).

Chương trình nghiên cứu 5-10 năm tiếp theo định hướng:

  • Mở rộng thuật toán tự chỉnh định thông số ma sát online (Adaptive Online Parameter Identification) sử dụng mạng nơ-ron học sâu (Deep Neural Networks) hoặc bộ lọc Kalman mở rộng (EKF).
  • Tích hợp mô hình RLuGre hoàn chỉnh trực tiếp vào chip FPGA/DSP để thực thi bộ bù ma sát vi mô ở tần số trích mẫu $> 10\text{ kHz}$.
  • Khảo sát ứng dụng trên các kết cấu xy lanh không cần (rodless cylinders) và động cơ khí nén quay chuyển vị góc.

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

  • Ảnh hưởng học thuật: Thiết lập chuẩn đối sánh (benchmark) mới trong mô hình hóa ma sát cho lĩnh vực Tự động thủy khí tại Việt Nam và khu vực, mở ra các công bố quốc tế chất lượng cao trên các tạp chí chuyên ngành uy tín (IEEE/ASME Transactions on Mechatronics, Mechanism and Machine Theory).
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp cho các doanh nghiệp chế tạo máy tự động hóa, dây chuyền đóng gói thực phẩm, dược phẩm, và thiết bị y tế; cho phép thay thế các cụm servo thủy lực đắt tiền bằng hệ khí nén tỉ lệ tiết kiệm chi phí $40 \div 60%$.
  • Lợi ích kinh tế - xã hội: Thúc đẩy sử dụng nguồn năng lượng sạch (khí nén), giảm thiểu ô nhiễm rò rỉ dầu khoáng của hệ thủy lực, tối ưu hóa hiệu suất năng lượng của hệ thống khí nén công nghiệp quốc gia.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ điện tử / Tự động hóa: Tiếp cận bộ dữ liệu thực nghiệm chuẩn, hệ phương trình toán học giải tích chi tiết và cấu trúc mô phỏng MATLAB/Simulink hoàn chỉnh của hệ TĐKN tỉ lệ.
  • Kỹ sư R&D trong các tập đoàn chế tạo robot và tự động hóa công nghiệp: Bản thiết kế hoàn chỉnh của bộ điều khiển MSSMC kết hợp bộ bù ma sát, sẵn sàng nhúng vào các bộ điều khiển khả trình PLC/PAC thương mại.
  • Các nhà hoạch định chính sách công nghệ và quản lý năng lượng: Cơ sở khoa học để xây dựng tiêu chuẩn định mức tiêu hao năng lượng cho hệ thống tự động hóa khí nén trong sản xuất xanh.

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 là việc hoàn thiện và chứng minh thực nghiệm tính vượt trội của mô hình ma sát LuGre cải tiến (RLuGre) trên cơ cấu chấp hành khí nén. Luận án mở rộng trực tiếp lý thuyết mô hình hóa LuGre cổ điển của Canudas de Wit et al. (1995) và phát triển các kết quả của Yanada & Sekikawa (2008) bằng cách tích hợp đồng thời: i) Động lực học độ dày màng bôi trơn $h$ phụ thuộc vận tốc và gia tốc; ii) Khâu vi phân pha dẫn bậc một $\sigma_2(v + T \dot{v})$ cho dải vận tốc cao; và iii) Hàm trễ phi tuyến $F(z)$ với bộ nhớ không cục bộ cho chế độ chuẩn bị trượt.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu quốc tế trước đây?

So với nghiên cứu của Tsai & Huang (2008) (chỉ sử dụng MSSMC thuần túy và coi ma sát là nhiễu chặn) và Drakunov et al. (chỉ bù ma sát nhớt), luận án đã thiết kế cấu hình phần cứng van tỉ lệ lưu lượng kép điều khiển độc lập hai khoang, kết hợp quy trình nhận dạng tách biệt thông số ma sát tĩnh/động với bộ lọc số phi nhân quả $32\text{ Hz}$, đồng thời chứng minh ổn định tiệm cận toàn cục hệ kín theo tiêu chuẩn Lyapunov.

So sánh cấu trúc điều khiển:
Tsai & Huang (2008): [Tín hiệu đặt] -> [Bộ điều khiển MSSMC] ---------------> [Hệ Khí Nén] (Ma sát bị bỏ qua)
Luận án này:         [Tín hiệu đặt] -> [Bộ điều khiển MSSMC] -> (+) ---------> [Hệ Khí Nén]
                                                 ^               ^
                                                 |               | (Tín hiệu bù u_comp)
                                                 +-- [BỘ BÙ MA SÁT LUGRE]

3. Phát hiện bất ngờ nhất từ thực nghiệm được lý giải như thế nào?

Phát hiện bất ngờ nhất là hằng số thời gian pha dẫn $T$ có sự bất đối xứng rất lớn giữa hai chiều chuyển động: $T = 0{,}33\text{ s}$ ở hành trình thuận nhưng giảm xuống $T = 0{,}07\text{ s}$ ở hành trình nghịch. Hiện tượng này được giải thích về mặt cơ học do sự bất đối xứng tiết diện hữu ích của pít-tông ($A_1 \ne A_2$) kết hợp với sự mất cân bằng lực căng của cơ cấu phớt làm kín cần pít-tông khi chịu tải nén so với tải kéo.

4. Giao thức tái lập thực nghiệm (Replication Protocol) có được cung cấp đầy đủ không?

Toàn bộ thông số phần cứng (mã van SMC VEF3121-1-02, xy lanh CM2L25-300, DAQ USB-4711A), thuật toán lọc vi phân $32\text{ Hz}$, chu kỳ lấy mẫu $\Delta t = 1{,}16\text{ ms}$, mã nguồn Visual C++ và khối mô phỏng MATLAB/Simulink đều được lập tài liệu chi tiết, cho phép các nhóm nghiên cứu độc lập tái lập $100%$ kết quả.

5. Kế hoạch nghiên cứu 10 năm tới được phác thảo ra sao?

Kế hoạch 10 năm tập trung vào: i) Nhúng trực tiếp toàn bộ mô hình RLuGre vào vi điều khiển biên (Edge AI Microcontrollers); ii) Tích hợp thuật toán bù thích nghi tự học sâu với điều kiện môi trường nhiệt ẩm công nghiệp; iii) Chuẩn hóa thành gói thư viện điều khiển mở cho các hệ chấp hành khí nén y tế và robot công nghiệp.

Kết luận

Luận án "Nghiên cứu ứng dụng mô hình ma sát trong mô phỏng và điều khiển xy lanh khí nén" đã giải quyết triệt để vấn đề khoa học cốt lõi tồn tại nhiều năm trong ngành Máy và Tự động Thủy khí thông qua 5 đóng góp cụ thể:

  1. Đánh giá định lượng toàn diện ảnh hưởng của ba mô hình ma sát (Trạng thái ổn định, LuGre, và RLuGre) đến độ chính xác mô phỏng động lực học hệ thống TĐKN tỉ lệ.
  2. Khẳng định mô hình LuGre cải tiến (RLuGre) là mô hình tối ưu nhất, có khả năng mô tả trọn vẹn các hiệu ứng phi tuyến tĩnh, động, trễ vận tốc cao và trễ bộ nhớ không cục bộ trước khi trượt.
  3. Phát triển thành công bộ điều khiển phi tuyến mới kết hợp phương pháp trượt đa mặt trượt (MSSMC) với bộ bù ma sát động học LuGre, chứng minh toán học tính ổn định theo Lyapunov.
  4. Kiểm chứng thực nghiệm nghiêm ngặt trên hệ thống phần cứng chuyên dụng, chứng minh khả năng triệt tiêu trượt - dính và kiểm soát sai số vị trí đỉnh dưới $4%$.
  5. Xây dựng bộ công cụ mô phỏng chuẩn hóa trên MATLAB/Simulink và phần mềm điều khiển thời gian thực trên C++, mở ra hướng phát triển các hệ thống servo khí nén thông minh thế hệ mới.