Tổng quan về luận án

Hệ thống bồn chứa chất lỏng ghép nối (Coupled-Tank System - CTS) là bài toán kiểm soát điểm chuẩn (benchmark control problem) kinh điển được Ủy ban Liên đoàn Quốc tế về Điều khiển Tự động (IFAC, 1990) xác lập nhằm thẩm định hiệu năng của các công cụ và giải thuật điều khiển tiên tiến. Động lực học của CTS phản ánh trực tiếp các hệ thống công nghiệp cốt lõi như tháp chưng cất phân đoạn, hệ thống bao hơi nhà máy nhiệt điện, dây chuyền lọc hóa dầu và các trạm xử lý nước thải quy mô lớn. Đặc tính vật lý của CTS tiềm ẩn những thách thức điều khiển nghiêm trọng bao gồm tính phi tuyến cao bắt nguồn từ định luật xả Torricelli, đặc tính động học pha không cực tiểu (non-minimum phase) do sự xuất hiện của các điểm không nằm bên phải mặt phẳng phức (Right Half Plane zeros), cùng sự tương tác đa biến (MIMO) phức tạp khi dòng chất lỏng lưu thông hai chiều giữa các bồn chứa.

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự bất đối xứng giữa các công trình nghiên cứu lý thuyết mô phỏng thuần túy và tính khả thi trong thực thi thời gian thực trên phần cứng thực nghiệm. Trong khi phần lớn các nghiên cứu trước đây (chẳng hạn như mô hình tuyến tính hóa chuỗi Taylor trong các công trình kinh điển) thường bỏ qua các thành phần phi tuyến bậc cao và không lượng hóa được độ bất định tham số, việc triển khai các giải thuật hiện đại trên giàn thử nghiệm vật lý chịu ràng buộc cơ cấu chấp hành nghiêm ngặt (điện áp điều khiển bơm 0–5V DC biến đổi thành tín hiệu PWM 24V) và nhiễu tải đột ngột (load disturbances) do mở van xả vẫn chưa được giải quyết triệt để.

Nghiên cứu tập trung giải quyết các câu hỏi và giả thuyết khoa học cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để trích xuất mô hình toán học thực nghiệm có độ chính xác cao cho CTS thông qua kỹ thuật nhận dạng hệ thống hộp đen thay vì các mô hình giải tích giải phóng từ nguyên lý đầu tiên vốn chứa nhiều sai số tham số?
  2. Câu hỏi nghiên cứu 2 (RQ2): Phương pháp xây dựng bộ điều khiển tích phân tỉ lệ (PI) dựa trên tối ưu hóa lồi Bất đẳng thức Ma trận Tuyến tính (LMI) có thể khắc phục hoàn toàn nhược điểm đáp ứng chậm và độ vọt lố lớn của phương pháp Ziegler-Nichols truyền thống hay không?
  3. Câu hỏi nghiên cứu 3 (RQ3): Cơ chế thích nghi mờ nào có thể triệt tiêu hiện tượng rung giật tần số cao (chattering) của Điều khiển Chế độ Trượt (SMC) mà vẫn đảm bảo tính bền vững (robustness) tuyệt đối khi xảy ra nhiễu tải thời gian thực?

Khung lý thuyết của công trình tích hợp bốn trụ cột học thuật: Định luật Bảo toàn Khối lượng và Thủy động lực học Bernoulli, Lý thuyết Điều khiển Tối ưu Tuyến tính phản hồi trạng thái giải bằng LMI, Lý thuyết Ổn định Lyapunov cho Hệ thống Cấu trúc Biến (VSC/SMC), và Hệ suy diễn mờ (FLC) thích nghi cấu trúc kép. Nghiên cứu thực hiện trên giàn phần cứng Advantech PCI-1711 DAQ, bộ khuếch đại PSUPA và môi trường MATLAB/Simulink Real-Time Workshop với bộ dữ liệu kích thích thực nghiệm gồm $N = 6.000$ mẫu thời gian thực (chu kỳ trích mẫu $T_s = 0.1\text{s}$, thời gian vận hành $600\text{s}$), tạo nên bước đột phá về tính chuẩn xác và khả năng ứng dụng thực tiễn.

Literature Review và Positioning

Lịch sử phát triển các giải thuật điều khiển hệ thống bồn chứa ghép nối chứng kiến bốn trào lưu nghiên cứu chính với những đóng góp và hạn chế rõ rệt:

Trào lưu thứ nhất khởi nguồn từ việc xây dựng mô hình giải tích dựa trên phương trình cân bằng khối lượng và nguyên lý Bernoulli (Johansson, 2000; Astrom & Bell, 2000). Các mô hình này phụ thuộc vào việc khai triển chuỗi Taylor quanh điểm làm việc tĩnh. Tuy nhiên, việc triệt tiêu các số hạng bậc cao dẫn đến sai số động học lớn khi hệ thống dịch chuyển khỏi điểm cân bằng, đồng thời không phản ánh được tính chất biến thiên thời gian của hệ số lưu lượng qua van.

Trào lưu thứ hai tập trung vào việc hiệu chỉnh bộ điều khiển PID kinh điển. Ziegler và Nichols (1942) đã thiết lập nền tảng hiệu chỉnh thực nghiệm, nhưng bộc lộ độ vọt lố nghiêm trọng và thời gian xác lập kéo dài khi áp dụng cho hệ thống pha không cực tiểu. Skogestad (2003) đề xuất phương pháp hiệu chỉnh SIMC dựa trên mô hình nội (IMC), cải thiện khả năng đáp ứng nhưng vẫn nhạy cảm trước nhiễu ngẫu nhiên. Astrom và Wittenmark (1995) đưa ra giải pháp điều khiển thích ứng tham chiếu mô hình (MRAC) dựa trên quy tắc MIT, song cấu trúc này dễ mất ổn định khi xuất hiện các điểm kỳ dị trong ước lượng tham số trực tuyến.

Trào lưu thứ ba ứng dụng lý thuyết điều khiển bền vững và tối ưu hóa không gian trạng thái, điển hình là định lý cạnh (Edge Theorem) và phân vùng D của Neimark (1992), cùng các cấu trúc $H_\infty$ và điều khiển dự báo mô hình (MPC) (Camacho & Bordons, 2004). Mặc dù giải quyết được bài toán ràng buộc, chi phí tính toán trực tuyến của MPC gây khó khăn cho việc nhúng trên các vi điều khiển công nghiệp tốc độ cao.

Trào lưu thứ tư phát triển các hệ thống cấu trúc biến (Variable Structure Systems - VSC) và điều khiển chế độ trượt (Sliding Mode Control - SMC) (Utkin, 1977; Slotine & Li, 1991). Mặc dù SMC có ưu thế vượt trội về tính bất biến trước các bất định phù hợp (matched uncertainties), vấn đề rung giật (chattering) do hàm đóng cắt không liên tục gây tổn hại nghiêm trọng đến cơ cấu truyền động bơm và van điều tiết.

                  ┌─────────────────────────────────────────┐
                  │   Mô hình giải tích & Tuyến tính hóa   │
                  │        (Johansson, 2000)                │
                  └────────────────────┬────────────────────┘
                                       │ (Bỏ qua số hạng bậc cao)
                                       ▼
                  ┌─────────────────────────────────────────┐
                  │    Điều khiển PID / IMC / MRAC          │
                  │  (Ziegler-Nichols, 1942; Skogestad, 2003)│
                  └────────────────────┬────────────────────┘
                                       │ (Vọt lố lớn, nhạy cảm nhiễu)
                                       ▼
                  ┌─────────────────────────────────────────┐
                  │   Điều khiển Tối ưu & Chế độ Trượt SMC  │
                  │     (Utkin, 1977; Slotine & Li, 1991)    │
                  └────────────────────┬────────────────────┘
                                       │ (Hiện tượng Chattering)
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────────┐
│              ĐỊNH VỊ CỦA NGHIÊN CỨU HIỆN TẠI (Mahapatro & Subudhi)              │
│  - Nhận dạng Output Error (OE) với $N = 6.000$ mẫu trắng                       │
│  - Tối ưu hóa lồi LQR-LMI giải qua Bổ đề Schur (Solver MINCX)                  │
│  - Điều khiển Thích nghi Mờ - Chế độ Trượt (AFSMC) triệt tiêu Chattering       │
└─────────────────────────────────────────────────────────────────────────────────┘

Nghiên cứu định vị tại điểm giao thoa giữa điều khiển cấu trúc biến và kỹ thuật trí tuệ tính toán thích nghi. So với công trình của Johansson (2000) tại Đại học Lund về hệ Quadruple Tank Process vốn chỉ khảo sát các chế độ đa biến tuyến tính hóa và công trình của Skogestad (2003) tại NTNU Na Uy tập trung vào quy tắc điều chỉnh PI giải tích, nghiên cứu này vượt lên bằng việc thực thi đồng thời chu trình nhận dạng hệ thống thực nghiệm Output Error (OE) bậc hai và triển khai bộ điều khiển AFSMC thời gian thực có khả năng thích nghi tham số hạt nhân trực tuyến, loại bỏ hoàn toàn hiện tượng chattering và duy trì tính ổn định toàn cục.

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

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

Nghiên cứu mang lại những đóng góp nền tảng cho lý thuyết điều khiển tự động thông qua ba khía cạnh đột phá:

Thứ nhất, nghiên cứu mở rộng Khung Điều khiển Tối ưu Tuyến tính LQR bằng cách chuyển đổi bài toán thiết kế bộ điều khiển PI tăng cường trạng thái thành bài toán Tối ưu hóa Bất đẳng thức Ma trận Tuyến tính (LMI). Bằng việc bổ sung biến trạng thái tích phân sai số $x_I = \int_0^t (r - y)d\tau$, hệ thống mở rộng được mô tả qua phương trình vi phân ma trận: $$\dot{\tilde{x}} = \tilde{A}\tilde{x} + \tilde{B}u$$ Thông qua việc áp dụng Bổ đề Schur (Schur Complement Lemma), phương trình Đại số Riccati (ARE): $$A^T P + P A - P B R^{-1} B^T P + Q = 0$$ được tái cấu trúc thành bài toán tối ưu hóa lồi tìm kiếm ma trận đối xứng xác định dương $\hat{P} = P^{-1}$ và vector phản hồi $\hat{Y} = K P^{-1}$: $$\begin{bmatrix} \hat{A}\hat{P} + \hat{P}\hat{A}^T - \tilde{B}\hat{Y} - \hat{Y}^T\tilde{B}^T & \hat{P}Q^{1/2} & \hat{Y}^T R^{1/2} \ Q^{1/2}\hat{P} & -I & 0 \ R^{1/2}\hat{Y} & 0 & -I \end{bmatrix} < 0$$ Cách tiếp cận này loại bỏ hoàn toàn sự phụ thuộc vào phương pháp dò tìm thử - sai của các kỹ thuật kinh điển.

Thứ hai, nghiên cứu phát triển cấu trúc Bộ điều khiển PI Mờ Thích nghi (Adaptive Fuzzy PI - AFPI) song song độc đáo, phân tách hành vi điều khiển thành hai nhánh: nhánh sơ cấp (primary) giữ vai trò xác lập hệ số hạt nhân cố định $(k_{p0}, k_{i0})$ dựa trên đáp ứng Ziegler-Nichols, và nhánh thứ cấp (auxiliary) tạo ra các tham số biến thiên liên tục $(\Delta k_p(t), \Delta k_i(t))$ thông qua Hệ suy diễn mờ Mamdani: $$u(t) = \left[k_{p0} + \Delta k_p(t)\right]e(t) + \int_0^t \left[k_{i0} + \Delta k_i(t)\right]e(\tau)d\tau$$ Cấu trúc này tạo bước chuyển dịch mô hình (paradigm shift) từ điều chỉnh tham số tĩnh sang thích ứng động lực học phi tuyến theo thời gian thực.

Thứ ba, nghiên cứu hoàn thiện chứng minh Ổn định Lyapunov cho Bộ điều khiển Chế độ Trượt Mờ Thích ứng (AFSMC). Bề mặt trượt phi tuyến được xác định bởi: $$S(t) = \left(\frac{d}{dt} + \lambda\right)^{n-1} e(t)$$ Hàm Lyapunov $V(S) = \frac{1}{2}S^2$ được bảo đảm thỏa mãn điều kiện tiếp cận bền vững $\dot{V}(S) = S\dot{S} \le -\eta |S|$ ($\eta > 0$) thông qua việc sử dụng FLC để xấp xỉ liên tục thành phần đóng cắt gián đoạn $u_{sw} = K_{sw}\text{sgn}(S)$, triệt tiêu năng lượng tần số cao gây rung chấn cơ học.

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

Khung phân tích của luận án tích hợp bốn cấu trúc lý thuyết thành một hệ thống điều khiển đa tầng hoàn chỉnh:

                                  KHUNG PHÂN TÍCH ĐA TẦNG
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ TẦNG 1: NHẬN DẠNG HỆ THỐNG THỰC NGHIỆM (System Identification)                          │
│ - Tín hiệu kích thích: White Noise ($600\text{s}$, $T_s = 0.1\text{s}$, 6000 mẫu)     │
│ - Cấu trúc mô hình: Output Error (OE) bậc hai                                         │
│ - Thuật toán ước lượng: Least Squares (LS) tối thiểu hóa MSE                           │
│ - Thẩm định: Auto-correlation Analysis trong khoảng tin cậy 95%                        │
└───────────────────────────────────────────┬────────────────────────────────────────────┘
                                            ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ TẦNG 2: TỐI ƯU HÓA LỒI KHÔNG GIAN TRẠNG THÁI (LMI-LQR Formulation)                     │
│ - Biến đổi trạng thái tăng cường (Augmented State-Space)                               │
│ - Tối ưu hóa hàm mục tiêu Riccati qua Bổ đề Schur                                     │
│ - Bộ giải điểm trong (Interior-Point MINCX Solver) xác định ma trận độ lợi $K$         │
└───────────────────────────────────────────┬────────────────────────────────────────────┘
                                            ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ TẦNG 3: ĐIỀU KHIỂN BỀN VỮNG & CẤU TRÚC THÍCH NGHI PHI TUYẾN (AFSMC Architecture)       │
│ - Bề mặt trượt phi tuyến $S(t) = \dot{e}(t) + c e(t)$ bảo đảm ổn định tiệm cận       │
│ - Hệ suy diễn mờ Mamdani 25 luật làm trơn hàm chuyển mạch $\text{sgn}(S)$             │
│ - Thẩm định Lyapunov đa biến $\dot{V}(S) \le -\eta |S|$                                │
└───────────────────────────────────────────┬────────────────────────────────────────────┘
                                            ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ TẦNG 4: THỰC THI THỜI GIAN THỰC (Hardware-In-The-Loop Implementation)                   │
│ - Real-Time Workshop (RTW) sinh mã C tự động qua Target Language Compiler (TLC)        │
│ - Bo ghép nối PCI-1711 DAQ và khối công suất PSUPA (Biến đổi 0-5V DC sang PWM 24V)    │
│ - Kiểm soát bù trừ nhiễu tải đột ngột khi mở van xả thực nghiệm                        │
└────────────────────────────────────────────────────────────────────────────────────────┘

Điều kiện biên hoạt động (boundary conditions) được xác lập nghiêm ngặt: dải điện áp bơm danh định $u \in [0, 5]\text{V}$, mức nước giới hạn $h_1, h_2 \in [0, 30]\text{cm}$, tiết diện bồn $A = 0.01389\text{ m}^2$, diện tích lỗ xả $a_1 = a_2 = 0.1245\text{ cm}^2$, gia tốc trọng trường $g = 9.8\text{ m/s}^2$, và hệ số dòng bơm $\eta = 0.1194$.

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ẽ lập trường nhận thức luận Thực chứng Kỹ thuật Thực nghiệm (Empirical-Deductive Engineering Positivism). Toàn bộ các quy luật động học và giải thuật điều khiển đều bắt buộc phải vượt qua quy trình kiểm chứng kép: mô phỏng toán học giải tích trên MATLAB và thực thi vật lý thời gian thực trên giàn thiết bị thực nghiệm.

Thiết kế nghiên cứu đa cấp (multi-level design) bao gồm:

  • Cấp độ 1 (Mô hình hóa động học): Thiết lập hệ phương trình vi phân phi tuyến liên tục mô tả động lực học cân bằng chất lỏng giữa hai bồn: $$\frac{dh_1(t)}{dt} = \frac{\eta u(t)}{A} - \frac{a_1}{A}\sqrt{2gh_1(t)}$$ $$\frac{dh_2(t)}{dt} = \frac{a_1}{A}\sqrt{2gh_1(t)} - \frac{a_2}{A}\sqrt{2gh_2(t)}$$
  • Cấp độ 2 (Nhận dạng tham số thực nghiệm): Kích thích hệ thống bằng tín hiệu ngẫu nhiên trắng để thu thập quan hệ vào - ra không phụ thuộc giả định lý thuyết.
  • Cấp độ 3 (Tổng hợp và Nhúng giải thuật): Lập trình hóa các cấu trúc LMI, AFPI, SMC, và AFSMC thông qua mã nguồn C được biên dịch trực tiếp vào nhân thời gian thực của hệ điều hành.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu thực nghiệm tuân thủ tiêu chuẩn sáu bước nghiêm ngặt:

[Thiết kế kích thích White Noise] ──► [Thu thập dữ liệu DAQ PCI-1711] ──► [Tiền xử lý & Lọc nhiễu]
                                                                                   │
[Thẩm định tương quan 95%] ◄── [Ước lượng Bình phương nhỏ nhất] ◄── [Lựa chọn mô hình OE]

Dữ liệu được thu thập từ cảm biến mức áp suất đặt tại đáy bồn, tín hiệu đo chuyển đổi về dạng $0–5\text{V}$ DC qua mạch điều hòa tín hiệu. Để đảm bảo tính giá trị cấu trúc (construct validity), $6.000$ mẫu dữ liệu thực nghiệm được chia tách độc lập cho pha ước lượng (estimation) và pha thẩm định (validation).

Phương pháp kiểm định phi tham số Tự tương quan dư sai (Auto-correlation analysis of residuals) được áp dụng. Kết quả chỉ ra toàn bộ các hệ số tương quan của sai số ước lượng đều nằm hoàn toàn bên trong dải giới hạn tin cậy 95% ($\pm 1.96/\sqrt{N}$), chứng minh mô hình không bị bỏ sót bất kỳ thành phần động học có quy luật nào.

Data và phân tích

Phân tích nhận dạng hệ thống đã xác định cấu trúc mô hình sai số đầu ra (Output Error - OE) bậc hai cho hai bồn chứa: $$y(k) = \frac{B(q)}{F(q)}u(k) + e(k) = \frac{b_1 q^{-1} + b_2 q^{-2}}{1 + f_1 q^{-1} + f_2 q^{-2}}u(k) + e(k)$$

Hàm truyền liên tục tương đương thu được qua phép biến đổi song tuyến tính Tustin:

  • Bồn 1: $$T_1(s) = \frac{0.004711 s + 0.004877}{s^2 + 0.035 s + 0.00021}$$
  • Bồn 2: $$T_2(s) = \frac{0.002145 s + 0.003124}{s^2 + 0.042 s + 0.00035}$$

Trong thiết kế LMI-LQR, các ma trận trọng số tối ưu được lựa chọn thông qua ma trận bán xác định dương $Q$ và xác định dương $R$:

  • Bồn 1: $Q_1 = \begin{bmatrix} 2.4628 & 0 \ 0 & 2.9373 \end{bmatrix}$, $R_1 = 16.371$
  • Bồn 2: $Q_2 = \begin{bmatrix} 4.1938 & 0 \ 0 & 1.3268 \end{bmatrix}$, $R_2 = 0.0115$

Sử dụng bộ giải MINCX trong LMI Control Toolbox, vector độ lợi phản hồi trạng thái tối ưu được xác định chuẩn xác:

  • $K_{\text{Tank 1}} = \begin{bmatrix} 0.2761 & 0.4542 \end{bmatrix}$
  • $K_{\text{Tank 2}} = \begin{bmatrix} 0.8265 & 0.6640 \end{bmatrix}$

Hệ mờ trong AFSMC sử dụng 2 biến đầu vào là sai số mức $e(t)$ và đạo hàm sai số $\dot{e}(t)$, cùng 1 biến đầu ra là tín hiệu điều khiển bù $u_{fuzzy}$. Mỗi biến được định nghĩa qua 5 hàm liên thuộc tam giác chuẩn hóa trên vũ trụ đàm luận: Negative Large (NL), Negative Small (NS), Zero (ZE), Positive Small (PS), Positive Large (PL). Cơ sở quy tắc mờ gồm 25 luật IF-THEN được xây dựng dựa trên nguyên lý điều khiển chuyên gia và giải mờ bằng phương pháp trọng tâm (Centroid method).

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

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

Các phát hiện thực nghiệm lượng hóa sự vượt trội mang tính bước ngoặt của các giải thuật đề xuất so với phương pháp kinh điển:

Thuật toán điều khiển Đối tượng Thời gian tăng $t_r$ (s) Thời gian xác lập $t_s$ (s) Thời gian đỉnh $t_p$ (s) Độ vọt lố $M_p$ (%) ISE IAE
PI (Ziegler-Nichols) Bồn 1 8.65 29.17 10.50 4.00 43.55 64.20
PI dựa trên LMI Bồn 1 2.12 21.60 8.65 2.00 27.81 62.82
PI (Ziegler-Nichols) Bồn 2 16.46 67.74 27.48 25.00 27.33 39.54
PI dựa trên LMI Bồn 2 1.78 26.78 20.93 9.00 17.81 43.22
SMC Truyền thống Bồn 1 & 2 1.45 18.20 5.20 1.50 14.20 31.40
AFSMC Đề xuất Bồn 1 & 2 1.10 12.40 4.10 0.00 8.65 19.80 
So sánh Thời gian xác lập (Settling Time - ts) tại Bồn 2 (Pha không cực tiểu):
Ziegler-Nichols PI : [==================================================] 67.74s
LMI-based PI       : [====================] 26.78s  (-60.5%)
AFSMC Đề xuất      : [=========] 12.40s             (-81.7%)

Năm phát hiện cốt lõi từ dữ liệu thực nghiệm bao gồm:

Thứ nhất, đối với Bồn 1, bộ điều khiển PI dựa trên LMI rút ngắn thời gian tăng $75.5%$ (từ $8.65\text{s}$ xuống $2.12\text{s}$), giảm thời gian xác lập $25.9%$ (từ $29.17\text{s}$ xuống $21.60\text{s}$), triệt tiêu $50%$ độ vọt lố (từ $4%$ xuống $2%$), và giảm chỉ số tích phân bình phương sai số (ISE) từ $43.55$ xuống $27.81$ (giảm $36.14%$).

Thứ hai, tại Bồn 2 (nơi chịu tác động nặng nề của đặc tính pha không cực tiểu và ghép kênh thủy lực), giải thuật Ziegler-Nichols bộc lộ sự suy giảm hiệu năng nghiêm trọng với thời gian xác lập lên tới $67.74\text{s}$ và độ vọt lố $25%$. Ngược lại, bộ điều khiển LMI-PI kiểm soát thời gian xác lập ở mức $26.78\text{s}$ (giảm $60.47%$), khống chế độ vọt lố dưới $9%$, và hạ thấp chỉ số ISE từ $27.33$ xuống $17.81$ (giảm $34.83%$).

Thứ ba, bộ điều khiển PI Mờ Thích nghi (AFPI) chứng minh khả năng tự chỉnh định hệ số trực tuyến vượt trội, duy trì đáp ứng bám quỹ đạo chính xác trên toàn dải điểm đặt (setpoint tracking) mà không bị bão hòa cơ cấu chấp hành.

Thứ tư, thử nghiệm thực tế đối với bộ điều khiển SMC truyền thống ghi nhận hiện tượng chattering dao động điện áp tần số cao quanh bề mặt trượt $S=0$, gây rung giật cơ học liên tục ở động cơ bơm DC.

Thứ năm, giải thuật Điều khiển Chế độ Trượt Mờ Thích ứng (AFSMC) đề xuất đã triệt tiêu $100%$ hiện tượng chattering nhờ cơ chế làm trơn bề mặt trượt của FLC, đạt độ vọt lố $0%$, đưa thời gian xác lập về mức kỷ lục $12.40\text{s}$, đồng thời triệt tiêu hoàn toàn sai lệch mức nước khi có sự cố mở van xả cưỡng bức (load disturbance rejection).

Implications đa chiều

Về mặt học thuật, công trình thiết lập chuẩn mực phương pháp luận trong việc kết hợp tối ưu hóa lồi LMI và hệ mờ thích ứng cho lớp đối tượng phi tuyến pha không cực tiểu.

Về phương diện kỹ thuật công nghệ, quy trình tích hợp giữa MATLAB Real-Time Workshop, bộ chuyển đổi Target Language Compiler (TLC), và bo mạch DAQ cung cấp một khung thiết kế tạo mẫu nhanh (Rapid Control Prototyping) có khả năng chuyển giao trực tiếp cho các hệ thống nhúng công nghiệp như PLC, PAC hay vi xử lý DSP.

Về ứng dụng thực tiễn, giải thuật AFSMC giải quyết triệt để bài toán kiểm soát mức trong các hệ thống công nghiệp đòi hỏi độ an toàn cao: hệ thống cấp nước bao hơi nhà máy nhiệt điện (ngăn ngừa sự cố cạn nước gây nổ ống sinh hơi hoặc tràn nước sang tuabin), tháp tách phân đoạn dầu mỏ, và các bồn phản ứng hóa chất liên tục.

Về chính sách và tiêu chuẩn, nghiên cứu cung cấp cơ sở định lượng để xây dựng các quy chuẩn tiết kiệm năng lượng và bảo vệ thiết bị chấp hành, giúp kéo dài tuổi thọ của bơm và van điều tiết công nghiệp nhờ việc giảm thiểu tần suất đóng cắt dòng điều khiển.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận bốn hạn chế kỹ thuật:

  1. Mô hình thực nghiệm được thẩm định trên hệ thống hai bồn chứa ghép nối đơn dòng; các hiệu ứng động học phức tạp hơn của cấu hình bốn bồn đa biến (Quadruple Tank) với nhiều điểm không RHP dịch chuyển chưa được tích hợp đầy đủ.
  2. Giả định dòng chảy lý tưởng không xét đến độ trễ thời gian lớn (large time-delay) do chiều dài đường ống truyền dẫn trong các nhà máy quy mô siêu lớn.
  3. Ràng buộc cơ cấu chấp hành mới chỉ xét đến giới hạn biên độ điện áp bơm $(0–5\text{V})$, chưa mô hình hóa phi tuyến tính vùng chết (dead-zone) và giới hạn tốc độ biến thiên (rate limits) của động cơ.
  4. Cảm biến mức dạng áp suất thủy tĩnh nhạy cảm với bọt khí và biến thiên tỷ trọng chất lỏng khi nhiệt độ môi trường thay đổi.

Bốn hướng nghiên cứu mở rộng trong giai đoạn tới:

  • Hướng 1: Mở rộng giải thuật AFSMC sang cấu hình hệ thống phân tán đa biến (Decentralized Multi-variable CTS) có xét đến trễ mạng truyền thông công nghiệp (Networked Control Systems).
  • Hướng 2: Tích hợp giải thuật Điều khiển Chế độ Trượt Bậc Phân số (Fractional-Order SMC) kết hợp mạng Nơ-ron Mờ Tự thích nghi (ANFIS) để tối ưu hóa năng lượng điều khiển.
  • Hướng 3: Phát triển giải thuật Điều khiển Dự báo Mô hình Thần kinh - Mờ Bền vững (Robust Neuro-Fuzzy MPC) có khả năng xử lý đồng thời các ràng buộc trạng thái và cơ cấu chấp hành trong thời gian thực.
  • Hướng 4: Xây dựng giải pháp bản sao số (Digital Twin) trên nền tảng điện toán biên IIoT và giao thức OPC-UA phục vụ chẩn đoán lỗi thời gian thực (Fault Tolerant Control).

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

Về phương diện học thuật, nghiên cứu đóng góp một cơ sở dữ liệu thực nghiệm chuẩn mực cho cộng đồng điều khiển tự động quốc tế, với tiềm năng trích dẫn cao trong các tạp chí chuyên ngành hàng đầu như IEEE Transactions on Control Systems Technology, Control Engineering Practice, và Journal of Process Control.

Đối với ngành công nghiệp chế biến, việc triển khai giải thuật AFSMC giúp tối ưu hóa hiệu suất truyền nhiệt trong các lò hơi công nghiệp, giảm thiểu dao động mức nước, qua đó tiết kiệm ước tính $3–5%$ nhiên liệu tiêu thụ và giảm thiểu nguy cơ dừng máy ngoài kế hoạch.

Về mặt xã hội và môi trường, kiểm soát mức dòng chảy chính xác tại các nhà máy xử lý nước thải công nghiệp đảm bảo lưu lượng xả thải đạt chuẩn môi trường, ngăn ngừa nguy cơ tràn hóa chất độc hại ra nguồn nước tự nhiên.

Trên bình diện quốc tế, công trình đóng góp giải pháp thực thi khả thi cho các quốc gia đang phát triển trong việc hiện đại hóa các nhà máy công nghiệp cũ thông qua giải pháp điều khiển tiên tiến mà không cần thay thế toàn bộ phần cứng đắt tiền.

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

Nghiên cứu mang lại giá trị gia tăng trực tiếp cho bốn nhóm đối tượng chuyên môn:

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận hoàn chỉnh từ nhận dạng hệ thống hộp đen OE, biến đổi LMI qua Bổ đề Schur, đến kỹ thuật chứng minh ổn định Lyapunov cho hệ thích ứng mờ.
  • Kỹ sư R&D và Tự động hóa: Sở hữu bộ thông số hiệu chỉnh tối ưu $(K_{\text{LMI}}, k_{p0}, k_{i0}, \lambda, \eta)$, sơ đồ khối Simulink, và quy trình sinh mã C tự động qua Real-Time Workshop để ứng dụng ngay vào sản xuất.
  • Doanh nghiệp Vận hành Nhà máy Hóa chất & Nhiệt điện: Nâng cao độ tin cậy của hệ thống, giảm thiểu chi phí bảo trì cơ cấu bơm/van nhờ triệt tiêu hiện tượng chattering, và đảm bảo an toàn vận hành tuyệt đối.
  • Các Cơ quan Quản lý Tiêu chuẩn Công nghiệp: Cung cấp số liệu thực nghiệm chuẩn xác để xây dựng các khung định mức kỹ thuật về độ ổn định và chỉ số sai số tích phân (ISE, IAE) cho hệ thống điều khiển mức chất lỏng.

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ì? Đó là việc thiết lập cấu trúc tích hợp giữa Phương pháp Tối ưu hóa Lồi Bất đẳng thức Ma trận Tuyến tính (LMI) và Hệ Thích nghi Mờ - Chế độ Trượt (AFSMC). Nghiên cứu đã mở rộng Lý thuyết Điều khiển Tối ưu Kalman-Yakubovich-Popov thông qua việc chuyển đổi cấu trúc PI tăng cường sang bài toán phản hồi trạng thái giải bằng Bổ đề Schur, đồng thời chứng minh chặt chẽ tính ổn định tiệm cận Lyapunov cho mặt trượt được làm trơn bằng FLC.

  2. Đổi mới phương pháp luận thể hiện như thế nào khi so sánh với hai nghiên cứu quốc tế tiêu biểu? So với nghiên cứu của Johansson (2000) (vốn dựa trên mô hình tuyến tính hóa giải tích cục bộ) và nghiên cứu của Skogestad (2003) (dựa trên quy tắc điều chỉnh PI giải tích SIMC cho mô hình FOPDT), nghiên cứu này tạo bước đột phá phương pháp luận bằng việc: (1) Nhận dạng thực nghiệm thông qua mô hình Output Error (OE) bậc hai kích thích bởi tiếng ồn trắng $6.000$ mẫu; (2) Triển khai thực thi phần cứng thời gian thực (Hardware-in-the-loop) với bo DAQ PCI-1711 thay vì mô phỏng số thuần túy.

  3. Phát hiện thực nghiệm gây bất ngờ nhất được hỗ trợ bởi dữ liệu là gì? Đó là sự phân kỳ hiệu năng nghiêm trọng của bộ điều khiển PI Ziegler-Nichols tại Bồn 2 so với Bồn 1 do ảnh hưởng của đặc tính pha không cực tiểu. Trong khi tại Bồn 1, PI truyền thống đạt thời gian xác lập $29.17\text{s}$ và vọt lố $4%$, thì tại Bồn 2, thời gian xác lập bị kéo dài đột biến lên tới $67.74\text{s}$ cùng độ vọt lố $25%$. Bộ điều khiển LMI-PI đã giải quyết triệt để sự suy biến này khi duy trì thời gian xác lập ổn định ở mức $26.78\text{s}$ và khống chế độ vọt lố ở mức $9%$.

  4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không? Có. Luận án cung cấp đầy đủ thông số kích thước vật lý của hệ thống ($A = 0.01389\text{ m}^2, a_1 = a_2 = 0.1245\text{ cm}^2$), sơ đồ phần cứng kết nối qua Advantech PCI-1711 DAQ, thông số khối PSUPA (biến đổi $0–5\text{V}$ sang PWM $24\text{V}$), tập 25 luật mờ IF-THEN chi tiết, và các ma trận trọng số $Q, R$ dùng trong thuật toán tối ưu MINCX.

  5. Chương trình nghiên cứu 10 năm được phác thảo ra sao? Chương trình 10 năm định hướng tiến trình: (Năm 1–3) Mở rộng giải thuật AFSMC cho hệ Quadruple Tank đa biến phi tập trung; (Năm 4–6) Tích hợp giải thuật SMC bậc phân số và nơ-ron mờ thích nghi ANFIS xử lý trễ mạng công nghiệp; (Năm 7–10) Xây dựng nền tảng Bản sao số (Digital Twin) công nghiệp tích hợp AI chẩn đoán lỗi thời gian thực trên nền tảng đám mây và giao thức IIoT OPC-UA.

Kết luận

Luận án khẳng định sáu đóng góp khoa học và kỹ thuật nổi bật:

  1. Xây dựng thành công mô hình thực nghiệm sai số đầu ra (OE) bậc hai cho hệ thống bồn chứa ghép nối thông qua phương pháp nhận dạng hộp đen với $6.000$ mẫu kích thích tiếng ồn trắng, được kiểm định chuẩn xác qua phân tích tự tương quan dư sai đạt khoảng tin cậy 95%.
  2. Đề xuất quy trình chuyển đổi bài toán thiết kế bộ điều khiển PI sang bài toán phản hồi trạng thái giải bằng tối ưu hóa lồi LMI qua Bổ đề Schur, giảm $75.5%$ thời gian tăng và giảm $36.14%$ chỉ số sai số tích phân ISE tại Bồn 1 so với phương pháp Ziegler-Nichols.
  3. Giải quyết triệt để vấn đề suy biến đáp ứng động học do đặc tính pha không cực tiểu tại Bồn 2, rút ngắn thời gian xác lập từ $67.74\text{s}$ xuống $26.78\text{s}$ ($60.47%$) và khống chế độ vọt lố dưới $9%$.
  4. Phát triển cấu trúc điều khiển PI Mờ Thích nghi (AFPI) song song cho phép tự động điều chỉnh hệ số khuếch đại theo thời gian thực trên dải điểm đặt rộng lớn.
  5. Thiết kế và nhúng thành công giải thuật Điều khiển Chế độ Trượt Mờ Thích ứng (AFSMC) thời gian thực, triệt tiêu $100%$ hiện tượng rung giật chattering cơ khí, đạt độ vọt lố $0%$ và loại bỏ hoàn toàn sai lệch khi có nhiễu tải mở van đột ngột.
  6. Hoàn thiện quy trình tạo mẫu nhanh (Rapid Control Prototyping) đồng bộ giữa Simulink Real-Time Workshop, bộ biên dịch TLC và phần cứng Advantech PCI-1711 DAQ.

Công trình đánh dấu bước chuyển dịch mô hình nghiên cứu từ mô phỏng lý thuyết sang thực thi phần cứng thời gian thực bền vững, mở ra ba nhánh nghiên cứu đột phá về điều khiển trượt thích ứng thông minh, điều khiển phi tuyến đa biến cho các hệ thống pha không cực tiểu, và hệ thống nhúng điều khiển quá trình thời gian thực trong kỷ nguyên Công nghiệp 4.0.