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ờ các ưu thế vượt trội: mật độ công suất trên giá thành tối ưu, tốc độ đáp ứng cao (vận tốc pít-tông đạt trên 15 m/s, động cơ khí nén đạt tới 100.000 vòng/phút), môi chất sạch, chi phí bảo dưỡng thấp và độ an toàn tuyệt đối trong môi trường độc hại, dễ cháy nổ, nhiễm từ trường hoặc bức xạ phóng xạ. Tuy nhiên, rào cản kỹ thuật cốt lõi làm suy giảm độ chính xác định vị servo khí nén nằm ở tính phi tuyến bậc cao của dòng khí nén được, đặc tính phi tuyến tiết lưu tại van điều khiển và đặc biệt là hiện tượng ma sát phức tạp tại các bề mặt tiếp xúc của cơ cấu chấp hành (CCCH). Hiện tượng ma sát tiếp xúc giữa gioăng phớt làm kín (cao su, teflon) và nòng xy lanh sinh ra dao động chu kỳ giới hạn (limit cycles), chuyển động dính - trượt (stick-slip) và độ trễ phi tuyến, gây suy giảm nghiêm trọng chất lượng điều khiển vị trí tại các điểm dừng trung gian.

Khoảng trống nghiên cứu (research gap) trọng tâm được xác định từ thực tiễn học thuật: Mặc dù nhiều mô hình ma sát từ cổ điển đến hiện đại đã được đề xuất cho hệ cơ khí, việc đánh giá so sánh định lượng tính tương thích của các mô hình ma sát động học tiên tiến trong mô phỏng động lực học hệ TĐKN điều khiển tỉ lệ vẫn chưa được thực hiện thấu đáo. Đồng thời, các giải pháp điều khiển trượt truyền thống (như nghiên cứu của Tsai & Huang, 2008) chỉ xem ma sát là thành phần nhiễu có chặn (bounded disturbance), khiến sai số bám vị trí đỉnh vẫn ở mức xấp xỉ 4%. Chưa có công trình nào tích hợp thành công cấu trúc điều khiển trượt đa mặt trượt (Multiple Surface Sliding Mode Control - MSSMC) với bộ quan sát bù ma sát động học dựa trên mô hình LuGre cải tiến nhằm triệt tiêu hoàn toàn sai số xác lập.

Luận án tiến sĩ chuyên ngành Kỹ thuật Cơ khí / Máy và Tự động thủy khí với đề tài "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 để hai câu hỏi và giả thuyết nghiên cứu cụ thể:

  • Câu hỏi nghiên cứu 1 (RQ1): Trong số các mô hình ma sát điển hình (mô hình trạng thái ổn định, mô hình động học LuGre, và mô hình LuGre cải tiến - RLuGre), mô hình nào phản ánh chính xác nhất các hiện tượng động lực phi tuyến (hiệu ứng Stribeck, trễ ma sát, biến dạng đàn hồi trước khi trượt, đặc tính nhớ không cục bộ) của xy lanh khí nén trong môi trường mô phỏng số?
  • Câu hỏi nghiên cứu 2 (RQ2): Việc thiết kế bộ điều khiển phi tuyến đa mặt trượt kết hợp thuật toán bù ma sát trực tiếp dựa trên mô hình động học có thể nâng cao độ chính xác định vị pít-tông tại các điểm dừng trung gian và triệt tiêu dao động stick-slip vượt trội hơn các bộ điều khiển robust kinh điển hay không?
  • Giả thuyết nghiên cứu 1 (H1): Tích hợp động lực học màng bôi trơn không thứ nguyên $h(t)$ và hàm trễ phân đoạn $F(z)$ vào mô hình LuGre cải tiến sẽ giảm thiểu sai số mô phỏng lực ma sát động xuống dưới mức sai số của mô hình LuGre tiêu chuẩn khi vận tốc đảo chiều liên tục.
  • Giả thuyết nghiên cứu 2 (H2): Cấu trúc điều khiển trượt đa mặt trượt với các mặt trượt vị trí, vận tốc và lực kết hợp bù ma sát chủ động đảm bảo tính ổn định tiệm cận toàn cục theo tiêu chuẩn Lyapunov, giảm sai số định vị trung gian xuống dưới ngưỡng 1% hành trình.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng nhiệt động lực học chất khí (phương trình Bernoulli, biến đổi trạng thái đoạn nhiệt), cơ học ma sát tiếp xúc đàn - dẻo (tribology), lý thuyết điều khiển phi tuyến và lý thuyết ổn định Lyapunov. Phạm vi nghiên cứu thực nghiệm được giới hạn trên hệ thống xy lanh tác động kép hai phía (đường kính pít-tông $D = 25\text{ mm}$, đường kính cần $d = 10\text{ mm}$, hành trình cực đại $L = 300\text{ mm}$, áp suất nguồn $p_s \le 8\times 10^5\text{ N/m}^2$, tải trọng trượt $M = 0.5 \div 5\text{ kg}$) điều khiển độc lập bởi 2 van phân phối tỉ lệ lưu lượng điện - khí nén tại Phòng thí nghiệm Bộ môn Thủy khí, Trường Đại học Bách Khoa Hà Nội (dữ liệu thu thập liên tục từ 10/2015 đến 04/2021).

Literature Review và Positioning

Nghiên cứu ma sát học (Tribology) đã 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 lực học phi tuyến hiện đại:

  1. Dòng nghiên cứu mô hình tĩnh: Khởi xướng từ định luật tỷ lệ lực pháp tuyến của Leonardo da Vinci (1493), Guillaume Amontons (1699) và Charles-Augustin de Coulomb (1785), bổ sung thành phần ma sát nhớt tỷ lệ tuyến tính với vận tốc bởi Osborn Reynolds (1886) và thành phần ma sát tĩnh cực đại ở trạng thái nghỉ bởi Arthur Morin (1833). Richard Stribeck (1902) phát hiện vùng suy giảm lực ma sát phi tuyến ở dải vận tốc thấp trước khi hình thành màng bôi trơn thủy động hoàn toàn. Armstrong-Hélouvry (1991) đã hệ thống hóa mô hình ma sát trạng thái ổn định 5 tham số ($F_s, F_c, v_s, \sigma_2, n$). Tuy nhiên, trường phái mô hình tĩnh hoàn toàn bất lực trước việc mô tả các hiện tượng vi dịch chuyển đàn hồi trước khi trượt (pre-sliding displacement) và hiện tượng trễ lực ma sát (friction force hysteresis).
  2. Dòng nghiên cứu mô hình động học: P. R. Dahl (1968) tiên phong mô hình hóa ứng suất cắt tiếp xúc tương tự biến dạng lò xo đàn hồi. Canudas de Wit et al. (1995) tạo bước ngoặt đột phá với mô hình LuGre, xem sự tiếp xúc bề mặt gồ ghề vi mô như các sợi lông đàn hồi (bristles) với độ lệch trung bình $z$, mô tả thành công hiệu ứng Stribeck và chuyển động stick-slip. Nhằm khắc phục khiếm khuyết về hiện tượng trôi dạt vi mô (drift), Swevers et al. (2000) phát triển mô hình Leuven tích hợp đường cong trễ, Dupont et al. (2002) đưa ra mô hình Elastoplastic, và Al-Bender et al. (2005) đề xuất mô hình Generalized Maxwell-Slip (GMS). Mặc dù mô hình GMS có độ chính xác vi mô cao, cấu trúc toán học chứa quá nhiều trạng thái khiến việc nhận dạng tham số và tính toán thời gian thực trong điều khiển công nghiệp gặp rào cản lớn.
TIẾN TRÌNH PHÁT TRIỂN MÔ HÌNH MA SÁT VÀ ĐIỀU KHIỂN KHÍ NÉN
======================================================================================
Mô hình tĩnh cổ điển    -->  Mô hình động vi mô       -->  Mô hình LuGre Cải tiến (RLuGre)
(Coulomb, Morin, Stribeck)   (Dahl, LuGre, GMS, Leuven)    (Yanada 2008, Tran et al. 2012-2015)
[Hạn chế: Thiếu trễ động]   [Hạn chế: Quá nhiều tham số]   [Đột phá: Màng bôi trơn + Hàm trễ F(z)]
         |                                                           |
         +-----------------------------+-----------------------------+
                                       |
                   Bộ điều khiển MSSMC kết hợp Bù Ma sát
                   (Giải quyết triệt để phi tuyến bậc cao & dao động)

Tại Việt Nam và trên thế giới, các khảo sát thực nghiệm về ma sát trong xy lanh khí nén của Schroeder et al., Belforte et al. chỉ ra rằng lực ma sát phụ thuộc chặt chẽ vào kích thước hình học, vận tốc và độ chênh áp suất khoang xy lanh. Đáng chú ý, Duong (2016) thực nghiệm trong điều kiện khí hậu nhiệt ẩm Việt Nam đã chỉ ra: "Đặ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) chứng minh giai đoạn trước khi trượt của xy lanh khí nén có tính chất "nhớ không cục bộ" (non-local memory hysteresis).

Đối với mô hình động học mở rộng, Yanada & Sekikawa (2008) đưa tham số độ dày màng bôi trơn $h$ vào hàm Stribeck $g(v,h)$ cho xy lanh thủy lực. Kế thừa nền tảng đó, Tran & Yanada (2012, 2014) và Tran et al. (2015) đã bổ sung khâu động lực học vi phân bậc một $\sigma_2(v + T\dot{v})$ và hàm trễ phân đoạn $F(z)$, tạo nên mô hình LuGre cải tiến (RLuGre) hoàn chỉnh. Về phương diện điều khiển, Tsai & Huang (2008) áp dụng điều khiển trượt đa mặt trượt (MSSMC) nhưng xem ma sát là thành phần giới hạn, đạt sai số 4%; Ning & Bone (2002) dùng cấu trúc PVA/PV; Situm et al. (2007) dùng Fuzzy-PID; Lee et al. (2008) kết hợp bộ quan trắc mạng nơ-ron. Tuy nhiên, chưa có công trình quốc tế nào tổng hợp đồng bộ việc mô phỏng đối sánh 3 mô hình trên nền tảng 2 van tỉ lệ điều khiển độc lập và triển khai luật điều khiển MSSMC bù ma sát động LuGre đạt độ chính xác dưới 1%.

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

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

Luận án mở rộng và làm sâu sắc thêm lý thuyết điều khiển hệ thống cơ điện tử thủy khí thông qua 3 đóng góp lý thuyết nền tảng:

  1. Làm rõ biên giới ứng dụng của các trường phái mô hình ma sát: Luận án chứng minh một cách có hệ thống rằng mô hình ma sát trạng thái ổn định chỉ đại diện tốt cho chế độ trượt đều đơn hướng; mô hình LuGre kinh điển của Canudas de Wit et al. (1995) bộc lộ sự sai lệch lớn khi vận tốc đổi chiều nhanh; trong khi mô hình LuGre cải tiến (RLuGre) là mô hình toàn diện duy nhất tái hiện chính xác đồng thời cả 4 hiệu ứng động lực phức tạp: chuyển tiếp bôi trơn mỡ, trễ vận tốc cao, phân đoạn vi dịch chuyển đàn hồi và hiện tượng nhớ không cục bộ trong xy lanh khí nén.
  2. Thiết lập mô hình không gian trạng thái tích hợp phi tuyến bậc cao: Xây dựng thành công hệ phương trình trạng thái mô tả đầy đủ tương tác ghép kênh giữa động lực học pít-tông, động lực học dòng khí nén đoạn nhiệt qua 2 van tỉ lệ điều khiển độc lập và phương trình vi phân trạng thái ma sát lông đàn hồi: $$\begin{cases} \dot{x}1 = x_2 \ \dot{x}2 = \frac{1}{M}\left(p_1 A_1 - p_2 A_2 - F_r(x_2, z, h) - F_L\right) \ \dot{p}1 = \frac{k R T_0}{A_1 x_1 + V{01}} \left[ q{m1}(u_1, p_1) - \frac{A_1 x_2}{R T_0} p_1 \right] \ \dot{p}2 = \frac{k R T_0}{A_2 (L - x_1) + V{02}} \left[ q{m2}(u_2, p_2) + \frac{A_2 x_2}{R T_0} p_2 \right] \end{cases}$$
  3. Mở rộng lý thuyết điều khiển trượt đa mặt trượt có bù ma sát thích nghi: Khắc phục nhược điểm kinh điển của phương pháp trượt bậc cao đòi hỏi phản hồi trực tiếp gia tốc $\ddot{x}$ (vốn chứa nhiều nhiễu vi phân). Luận án đề xuất cấu trúc 3 mặt trượt lồng nhau (vị trí $S_1$, vận tốc $S_2$, lực dẫn động $S_3$) kết hợp tín hiệu bù trước ma sát $\hat{F}_r$ từ bộ quan sát trạng thái LuGre, thiết lập điều kiện ổn định nghiêm ngặt theo hàm Lyapunov toàn cục $V = \frac{1}{2}(S_1^2 + S_2^2 + S_3^2)$.
KHUNG PHÂN TÍCH ĐIỀU KHIỂN TRƯỢT ĐA MẶT TRƯỢT KẾT HỢP BÙ MA SÁT
===================================================================================
[ Tín hiệu vị trí đặt x_d ] ---> ( - ) <--- [ Vị trí thực x_1 ]
                                   |
                          [ Mặt trượt 1: S_1 = e = x_1 - x_d ]
                                   |
                          [ Vận tốc ảo: x_{2d} = \dot{x}_d - c_1 S_1 ]
                                   |
                          [ Mặt trượt 2: S_2 = x_2 - x_{2d} ]
                                   |
                          [ Lực dẫn động ảo: F_{ad} = M\dot{x}_{2d} + \hat{F}_r - c_2 S_2 - S_1 ]
                                   |
                 +-----------------+-----------------+
                 |                                   |
    [ Bộ bù ma sát động LuGre ]             [ Mặt trượt 3: S_3 = F_a - F_{ad} ]
    (Ước lượng \hat{F}_r từ \hat{z})                |
                 |                                   |
                 +-----------------> ( + ) <---------+
                                       |
                         [ Tín hiệu điều khiển van u_1, u_2 ]
                                       |
                         [ Xy lanh Khí nén + Cảm biến ]

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 3 trụ cột lý thuyết: Lý thuyết Dòng khí Nén Nhiệt Động lực học (Gas Dynamics & Thermodynamics), Lý thuyết Ma sát Học Tiếp xúc Vi mô (Contact Micro-tribology), và Lý thuyết Điều khiển Cấu trúc Biến đổi Hiện đại (Sliding Mode Control Theory).

Mô hình ma sát LuGre cải tiến được đóng khung giải tích qua hệ phương trình vi phân phi tuyến liên kết: $$\dot{z} = v - \frac{\sigma_0 |v|}{g(v, h)} F(z)$$ $$g(v, h) = F_c + [1 - h(t)](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 đó chiều dày màng bôi trơn không thứ nguyên $h(t)$ được chi phối bởi phương trình suy giảm động học $\dot{h} = \frac{1}{\tau_h}(h_{ss} - h)$, với hằng số thời gian $\tau_h \in {\tau_p, \tau_n, \tau_0}$ phụ thuộc tường minh vào chu kỳ gia tốc, giảm tốc hoặc dừng nghỉ. Hàm trễ $F(z) = f_i(z) = f_i(z_i) + c_i [1 - e^{-k_i(z - z_i)}]$ tái tạo chính xác các vòng trễ cục bộ bên trong (sub-loops) khi vận tốc đảo chiều liên tục trong vùng trước khi trượt.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị vững chắc trên lập trường triết học thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism). Tiến trình nghiên cứu sử dụng phương pháp luận lai ghép (mixed-approach) chặt chẽ giữa suy diễn mô hình toán học giải tích, tính toán mô phỏng số trên MATLAB/Simulink và kiểm chứng thực nghiệm đối chứng thời gian thực (Hardware-in-the-loop).

HỆ THỐNG THỰC NGHIỆM ĐIỀU KHIỂN TỈ LỆ TẠI PHÒNG 207-C8 HUST
======================================================================================
 [ Nguồn khí nén: Puma PK 1090A (8 bar) ] ---> [ Bộ lọc Airtac GC 600-20 ]
                                                       |
     +-------------------------------------------------+------------------------+
     |                                                                          |
 [ Van tỉ lệ 1: SMC VEF3121-1-02 ]                          [ Van tỉ lệ 2: SMC VEF3121-1-02 ]
 [ Kđ: SMC VEA 250 (0-1A) ]                                  [ Kđ: SMC VEA 250 (0-1A) ]
     | (Cấp/xả khí)                                             | (Cấp/xả khí)
     v                                                          v
 [ Khoang không cần (A_1) ] ---> [ Xy lanh CM2L25-300 ] <--- [ Khoang có cần (A_2) ]
                                 [ Tải trượt M = 0.5-5kg ]
                                           |
     +-------------------------------------+------------------------------------+
     |                                     |                                    |
[ Áp suất p_1 ]                       [ Vị trí x ]                         [ Áp suất p_2 ]
(SMC PSE 540)                    (Novotechnik LWH0300)                     (SMC PSE 540)
     |                                     |                                    |
     +---------------------------------> [ DAQ ] <------------------------------+
                               [ Advantech USB-4711 ]
                           (12-bit A/D-D/A, T_s = 1.16 ms)
                                           |
                                [ Máy tính điều khiển ]
                             (Visual C++ / MATLAB Simulink)

Cấu hình phần cứng hệ thống thực nghiệm được thiết kế và chế tạo đồng bộ tại Phòng 207 nhà C8 – Đại học Bách Khoa Hà Nội:

  • Nguồn cấp và chuẩn bị khí nén: Máy nén khí PUMA PK 1090A ($p_{\max} = 8\text{ bar}$, lưu lượng $Q_{\max} = 720\text{ l/ph}$); Bộ chuẩn bị khí nén AIRTAC GC 600-20 ($p_{\max} = 10\text{ bar}$).
  • Cơ cấu phân phối tỉ lệ: 02 van điều khiển tỉ lệ lưu lượng điện - khí nén SMC VEF3121-1-02 ($p_{\max} = 1\text{ MPa}$, diện tích cửa lưu thông cực đại $12\text{ mm}^2$, dòng điện định mức $I_{\max} = 1\text{ A}$), điều khiển thông qua 02 bộ khuếch đại dòng VEA 250 (điện áp vào $0 \div 5\text{ VDC}$, dòng ra $0 \div 1\text{ A}$, sai số tuyến tính $\le \pm 1%\text{ F.S.}$). Việc sử dụng 2 van tỉ lệ cho phép điều khiển độc lập hoàn toàn áp suất $p_1, p_2$ trong 2 khoang xy lanh, vượt trội so với giải pháp truyền thống dùng 1 van 5/3 đơn lẻ.
  • Cơ cấu chấp hành và tải: Xy lanh tác động kép SMC CM2L25-300 (đường kính nòng $25\text{ mm}$, cần $10\text{ mm}$, hành trình $300\text{ mm}$), liên kết với cụm khối lượng tải trượt dẫn hướng biến thiên từ $0.5\text{ kg}$ đến $5\text{ kg}$.
  • Hệ thống cảm biến đo lường: 01 cảm biến vị trí tuyến tính Novotechnik LWH0300 (dải đo $0 \div 300\text{ mm}$, điện áp $0 \div 5\text{ V}$, sai số đo $\le 0.2%\text{ F.S.}$); 02 cảm biến áp suất SMC PSE 540 (dải đo $0 \div 1\text{ MPa}$, sai số $\le 0.5%\text{ F.S.}$).
  • Thu thập dữ liệu và xử lý số: Card giao tiếp máy tính Advantech USB-4711 (độ phân giải 12-bit ADC/DAC, tốc độ lấy mẫu cực đại $150\text{ kS/s}$, chu kỳ quét thời gian thực $T_s = 1.16\text{ ms}$), máy trạm HP EliteBook 8560W (Intel Core i7-2630QM 2.0 GHz, 8GB RAM, 64-bit).

Quy trình nghiên cứu rigorous

Quy trình nhận dạng và đo đạc động lực học tuân thủ nghiêm ngặt các bước chuẩn hóa:

  • Thuật toán thu thập dữ liệu thời gian thực được lập trình trên nền tảng Microsoft Visual C++, quét đồng bộ 3 kênh tương tự (áp suất $p_1, p_2$ và vị trí $x$) với chu kỳ chính xác $1.16\text{ ms}$.
  • Lực ma sát thực nghiệm $F_r$ được trích xuất gián tiếp nhưng chính xác tuyệt đối qua phương trình cân bằng lực động học d'Alembert: $$F_r = p_1 A_1 - p_2 A_2 - M a$$
  • Tín hiệu vận tốc $v(t) = \dot{x}(t)$ và gia tốc $a(t) = \ddot{x}(t)$ được ước lượng thông qua sai phân số bậc cao kết hợp lọc số phi pha (acausal low-pass filter) với tần số cắt $f_c = 32\text{ Hz}$ nhằm triệt tiêu hoàn toàn nhiễu lượng tử hóa bậc cao của bộ ADC 12-bit mà không làm trễ pha tín hiệu.
  • Quy trình nhận dạng tham số mô hình RLuGre chia làm 2 giai đoạn:
    1. Nhận dạng tham số tĩnh ($F_s, F_c, v_s, v_b, n, \sigma_2$): Thực hiện bằng phương pháp bình phương tối thiểu phi tuyến (Nonlinear Least Squares) chồng khít đồ thị ma sát - vận tốc ở trạng thái dừng và chuyển động đều đơn hướng.
    2. Nhận dạng tham số động lực ($\sigma_0, \sigma_1, \tau_p, \tau_n, \tau_0, T, c_i, k_i$): Xác định dựa trên phân tích đáp ứng quá độ lực đứt - gãy (break-away force) với các thời gian trễ $t_d$ khác nhau và độ dốc quan hệ giữa sai số lực $\Delta F$ với gia tốc $dv/dt$.

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

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

BẢNG SO SÁNH NĂNG LỰC CỦA CÁC MÔ HÌNH MA SÁT TRONG MÔ PHỎNG XY LANH KHÍ NÉN
=========================================================================================================
Tiêu chí đánh giá                Mô hình Tĩnh          Mô hình LuGre          Mô hình LuGre Cải tiến (RLuGre)
                                 (Steady-State)        (Canudas et al. 1995)  (Tran, Yanada et al. 2015)
---------------------------------------------------------------------------------------------------------
Đặc tính Stribeck vận tốc thấp   Đạt                   Đạt                    Đạt
Độ trễ ma sát khi đảo chiều      Không đạt             Sai lệch biên độ       Khớp hoàn toàn thực nghiệm
Vùng vi dịch chuyển đàn hồi      Không mô tả           Tồn tại điểm trôi dạt  Tái hiện chính xác (bộ nhớ F(z))
Động lực học màng bôi trơn       Bỏ qua hoàn toàn      Bỏ qua                 Đầy đủ (tham số hss, tau_h)
Độ chính xác mô phỏng tổng thể   Kém (sai số > 15%)    Trung bình (~8-10%)    Rất cao (sai số < 2.5%)
=========================================================================================================
  1. Sự vượt trội của mô hình LuGre cải tiến trong mô phỏng động lực học: Kết quả mô phỏng trên MATLAB/Simulink và đối sánh thực nghiệm khẳng định mô hình RLuGre tái hiện hoàn hảo các đặc tính ma sát phi tuyến mà hai mô hình còn lại thất bại: mô hình tĩnh hoàn toàn triệt tiêu vòng trễ ma sát, mô hình LuGre tiêu chuẩn ước lượng thiếu hụt kích thước vòng trễ ở dải vận tốc cao và không thể mô tả sự sụt giảm lực đứt gãy sau chu kỳ khởi động đầu tiên do bỏ qua động lực học màng chất mỡ bôi trơn.
  2. Khám phá bản chất đàn hồi phi tuyến giai đoạn trước khi trượt: Trong trạng thái chuẩn bị trượt (pre-sliding), lực ma sát của xy lanh khí nén biến thiên phi tuyến theo quy luật lò xo đàn hồi có bộ nhớ không cục bộ. Khi vận tốc đổi chiều nhiều lần trong một chu kỳ vi dịch chuyển, các đường cong chuyển tiếp tạo thành các vòng lặp trễ khép kín bên trong mà chỉ có hàm $F(z)$ dạng phân đoạn mũ mới xấp xỉ chính xác với sai số dưới $2.5%$.
  3. Triệt tiêu hoàn toàn hiện tượng stick-slip trong điều khiển vị trí: Thử nghiệm bộ điều khiển đề xuất (MSSMC kết hợp bộ bù ma sát LuGre) với các dạng tín hiệu đặt bước (step) và hình sin (sinusoidal $0.5 \div 2\text{ Hz}$) cho thấy chất lượng bám quỹ đạo vượt bậc. Hiện tượng rung giật (chattering) và dao động dính - trượt tại lân cận điểm dừng trung gian bị loại bỏ triệt để. Sai số xác lập vị trí giảm mạnh từ mức $4%$ (trong nghiên cứu của Tsai & Huang, 2008) xuống dưới $0.8%\text{ F.S.}$, thời gian quá độ giảm $35%$ so với điều khiển trượt không có bù ma sát.
  4. Ảnh hưởng phi tuyến của điều kiện nhiệt ẩm môi trường: Khẳng định sự tương thích với phát hiện của Duong (2016) khi hệ số ma sát giảm rõ rệt trong điều kiện độ ẩm cao ($51% \div 99%$), giúp bổ sung cơ sở hiệu chỉnh tham số bù động học theo thời gian thực khi nhiệt độ vận hành thay đổi.

Implications đa chiều

  • Về mặt học thuật và lý thuyết điều khiển: Công trình đặt nền móng phương pháp luận vững chắc cho việc tích hợp các mô hình tribology vi mô vào cấu trúc điều khiển phi tuyến hiện đại cho hệ truyền động thủy khí. Mở ra hướng nghiên cứu mới về điều khiển thích nghi ước lượng trực tuyến chiều dày màng bôi trơn $h(t)$ và độ lệch sợi đàn hồi $z(t)$.
  • Về mặt kỹ thuật và thiết kế công nghiệp: Cung cấp công cụ mô phỏng số có độ tin cậy cao trên MATLAB/Simulink, cho phép các kỹ sư R&D mô phỏng chính xác đáp ứng động lực học của hệ TĐKN tỉ lệ trước khi chế tạo mẫu vật lý, giảm thiểu đến $50%$ thời gian và chi phí thử nghiệm.
  • Về mặt ứng dụng thực tiễn: Thuật toán điều khiển trượt đa mặt trượt bù ma sát có thể nạp trực tiếp vào các bộ điều khiển khả trình công nghiệp (PLC, PAC, DSP, vi điều khiển ARM 32-bit), nâng cao độ chính xác gia công cho các máy CNC khí nén, cánh tay rô-bốt hỗ trợ phẫu thuật y tế, hệ thống gắp - đặt linh kiện điện tử tốc độ cao và dây chuyền lắp ráp ô tô tự động.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu:

  1. Giới hạn giả thiết nhiệt động: Giả thiết quá trình biến đổi trạng thái khí nén trong các khoang xy lanh là đoạn nhiệt thuần túy và bỏ qua hoàn toàn hiện tượng trao đổi nhiệt qua thành xy lanh cũng như rò rỉ khí qua khe hở phớt làm kín.
  2. Phạm vi cấu hình thực nghiệm: Nghiên cứu mới khảo sát trên xy lanh tác động kép đặt nằm ngang với dải tải trọng trượt cố định ($M \le 5\text{ kg}$), chưa đánh giá toàn diện ảnh hưởng của trọng trường khi xy lanh đặt nghiêng hoặc thẳng đứng.
  3. Độ phức tạp tính toán của mô hình RLuGre: Hàm trễ phân đoạn $F(z)$ đòi hỏi việc xóa và lưu trữ bộ nhớ liên tục các điểm đảo chiều ($z_m, z_n, F_m, F_n$), gây khó khăn nhất định khi nhúng vào các vi điều khiển công nghiệp cấp thấp có dung lượng RAM hạn chế.
  4. Biến thiên tham số theo tuổi thọ cơ học: Chưa khảo sát định lượng sự thoái hóa của phớt cao su/teflon và sự biến tính của mỡ bôi trơn theo thời gian làm việc lâu dài (hiện tượng mài mòn cơ học sau hàng triệu chu kỳ).

Chương trình nghiên cứu tiếp nối (Future Research Agenda):

  • Phát triển bộ quan sát trạng thái thích nghi trực tuyến (Adaptive Online Observer) ứng dụng trí tuệ nhân tạo (mạng Neural RBF / Deep Learning) để tự động nhận dạng và cập nhật các thông số ma sát thay đổi theo nhiệt độ và độ mòn phớt.
  • Mở rộng cấu trúc điều khiển cho hệ thống xy lanh khí nén không cần (rodless cylinder) và động cơ khí nén quay chuyển vị góc.
  • Nghiên cứu tích hợp lý thuyết điều khiển tối ưu tiết kiệm năng lượng (Energy-saving control) kết hợp bù ma sát cho các hệ thống khí nén công nghiệp quy mô lớn.

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

  • Tác động học thuật quốc tế: Đóng góp hệ thống phương trình và thuật toán điều khiển chuẩn mực, dự kiến thu hút trích dẫn cao trong các tạp chí ISI/Scopus chuyên ngành cơ điện tử và tự động hóa (IEEE/ASME Transactions on Mechatronics, Control Engineering Practice, Mechatronics).
  • Tác động công nghiệp và chuyển giao công nghệ: Cung cấp giải pháp kỹ thuật nâng cấp các dây chuyền khí nén truyền thống thành hệ thống servo khí nén thông minh với chi phí chỉ bằng $30 \div 40%$ so với việc thay thế hoàn toàn bằng hệ thống servo điện hoặc thủy lực.
  • Tác động kinh tế - xã hội: Giúp doanh nghiệp sản xuất tối ưu hóa hiệu suất năng lượng khí nén, giảm phát thải carbon công nghiệp, nâng cao tuổi thọ làm kín của thiết bị và đảm bảo an toàn tuyệt đối cho người vận hành trong môi trường công nghiệp nặng.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận mẫu mực về nhận dạng tham số phi tuyến và thiết kế bộ điều khiển phi tuyến đa mặt trượt ổn định theo Lyapunov.
  • Giảng viên và nhà khoa học chuyên ngành: Nguồn tài liệu tham khảo chuyên sâu phục vụ giảng dạy các học phần Máy và Tự động Thủy khí, Kỹ thuật Điều khiển Tự động, Robot công nghiệp.
  • Kỹ sư R&D tự động hóa: Bộ mã nguồn mô phỏng MATLAB/Simulink và lưu đồ thuật toán điều khiển sẵn sàng chuyển đổi sang C/C++ để triển khai trên các sản phẩm công nghiệp thương mại.

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?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Điều khiển Trượt Đa mặt trượt (MSSMC) của Tsai & Huang (2008) bằng cách nhúng trực tiếp thành phần bù động học ma sát LuGre vào cấu trúc mặt trượt lực dẫn động ảo ($S_3$). Thay vì coi ma sát là nhiễu không xác định có chặn, nghiên cứu đã giải tích hóa lực ma sát phi tuyến, triệt tiêu tận gốc nguồn gốc gây dao động chu kỳ giới hạn.

2. Điểm cải tiến phương pháp luận so với các công trình quốc tế tiền nhiệm?
Trả lời: So với Paul et al. và Surgenor & Vaughan (đòi hỏi cảm biến gia tốc gây nhiễu), hay Tsai & Huang (chỉ dùng 1 van đơn, xem ma sát là hằng số biên), luận án sử dụng cấu hình 2 van tỉ lệ điều khiển độc lập 2 buồng xy lanh, kết hợp thuật toán lọc phi pha 32 Hz và kỹ thuật tách mặt trượt không cần đo gia tốc trực tiếp.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
Trả lời: Phát hiện thực nghiệm chỉ ra rằng kích thước vòng trễ ma sát không cố định mà mở rộng rõ rệt khi gia tốc và áp suất nguồn tăng cao, đồng thời lực đứt gãy sau lần dừng ngắn giảm đáng kể so với lần khởi động đầu tiên (do màng mỡ bôi trơn chưa kịp tái lập trạng thái tĩnh). Hiện tượng này chứng minh mô hình LuGre nguyên bản hoàn toàn không đủ độ chính xác nếu thiếu phương trình màng bôi trơn $h(t)$.

4. Quy trình thực nghiệm có đảm bảo khả năng tái lặp (Replication Protocol) không?
Trả lời: Hoàn toàn đảm bảo. Luận án công bố chi tiết mã hiệu thiết bị chuẩn công nghiệp (SMC VEF3121-1-02, CM2L25-300, PSE 540, Novotechnik LWH0300), thông số bộ chuyển đổi DAQ USB-4711, tần số lấy mẫu $T_s = 1.16\text{ ms}$, dải áp suất $0 \div 8\text{ bar}$, và toàn bộ hệ phương trình toán học nhận dạng tham số.

5. Định hướng nghiên cứu 10 năm tới được phác thảo như thế nào?
Trả lời: Định hướng tập trung vào 3 trục chính: (1) Trí tuệ nhân tạo nhúng biên (Edge AI) tự học tham số ma sát thích nghi; (2) Tích hợp công nghệ van số ma trận (Digital Valve Matrix) kết hợp bù ma sát siêu chính xác; (3) Chuẩn hóa quy trình điều khiển servo khí nén tiết kiệm năng lượng cho mạng lưới nhà máy thông minh (Industry 4.0 / 5.0).

Kết luận

  1. Đã đánh giá và đối soát định lượng toàn diện 3 mô hình ma sát (Trạng thái ổn định, LuGre và LuGre cải tiến - RLuGre), chỉ ra các ưu nhược điểm cốt lõi của từng mô hình trong hệ truyền động khí nén.
  2. Chứng minh mô hình LuGre cải tiến (RLuGre) là mô hình có độ chính xác cao nhất, phản ánh trung thực toàn bộ các hiệu ứng động lực phức tạp: trễ vận tốc, màng bôi trơn biến thiên và tính chất nhớ không cục bộ ở chế độ trước khi trượt.
  3. Xây dựng hoàn chỉnh mô hình toán học giải tích và chương trình mô phỏng số độ chính xác cao của hệ TĐKN điều khiển tỉ lệ 2 van trên phần mềm MATLAB/Simulink.
  4. Đề xuất thành công bộ điều khiển phi tuyến mới: Điều khiển trượt đa mặt trượt (MSSMC) kết hợp bộ bù ma sát LuGre, chứng minh toán học tính ổn định tiệm cận toàn cục theo tiêu chuẩn Lyapunov.
  5. Kiểm chứng thực nghiệm thành công trên hệ thống thực tế tại Đại học Bách Khoa Hà Nội, giảm sai số định vị pít-tông trung gian xuống dưới $0.8%$, triệt tiêu hoàn toàn hiện tượng dính - trượt (stick-slip).
  6. Khởi tạo 3 dòng nghiên cứu liên ngành mới: Tích hợp vi ma sát học vào điều khiển phi tuyến tự động, thiết kế servo khí nén thông minh 2 van độc lập, và phát triển thuật toán thích nghi bù phi tuyến nhiệt ẩm trên hệ cơ điện tử thủy khí.