Tổng quan về luận án

Nghiên cứu điều khiển tối ưu cho các hệ thống nhiệt công nghiệp là một trong những bài toán phức tạp nhất thuộc lĩnh vực Kỹ thuật Điều khiển và Tự động hóa. Đối với các quá trình công nghệ xử lý nhiệt như tôi, ram, nhiệt luyện chi tiết cơ khí, ủ vật liệu từ hay nung gốm sứ cao cấp, đối tượng gia nhiệt có bản chất phân bố không gian và thời gian rõ rệt. Luận án tiến sĩ kỹ thuật mang tên "Nghiên cứu điều khiển tối ưu cho hệ với tham số phân bố, có trễ, phi tuyến" do NCS. Mai Trung Thái thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Hữu Công tại Đại học Thái Nguyên là một công trình tiên phong giải quyết trọn vẹn bài toán điều khiển tối ưu theo tiêu chuẩn "nung chính xác nhất" cho hệ lò điện trở - vật nung dạng tấm phẳng dày.

flowchart LR
    subgraph Control ["Kênh Điều Khiển"]
        u["Điện áp u(t)"] --> thyristor["BĐĐ Thyristor"]
    end
    subgraph Furnace ["Lò Điện Trở (LPS Trễ - Phi Tuyến)"]
        thyristor --> Ws["W(s) = k(v)e^(-τs)/(Ts+1)"]
        Ws --> vt["Nhiệt độ lò v(t)"]
    end
    subgraph Material ["Vật Nung (DPS Parabolic)"]
        vt --> pde["∂²q/∂x² = (1/a)(∂q/∂t)"]
        pde --> qxt["Trường nhiệt q(x,t)"]
    end
    qxt --> Jc["Jc = ∫[q*(x) - q(x,tf)]²dx → min"]

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án xác định xuất phát từ giới hạn của các công bố kinh điển: đa số các công trình trước đây chỉ khảo sát hệ tham số tập trung (Lumped Parameter Systems - LPS) hoặc hệ tham số phân bố tuyến tính (Distributed Parameter Systems - DPS) với các tham số cố định và bỏ qua trễ. Khi áp dụng vào lò đốt điện trở, hệ số truyền tĩnh $k$ thực tế biến thiên phi tuyến rất mạnh theo nhiệt độ môi trường lò ($k = k(v)$ trong dải $0 - 500^\circ\text{C}$), đồng thời tồn tại thời gian trễ thuần $\tau$ đáng kể so với hằng số thời gian $T$ (thỏa mãn tỉ số $6 \le T/\tau < 10$).

Câu hỏi nghiên cứu và giả thuyết khoa học của đề tài được xác lập cụ thể:

  1. RQ1: Làm thế nào để mô hình hóa giải tích chính xác mối quan hệ giữa điện áp cung cấp $u(t)$ và phân bố trường nhiệt độ $q(x,t)$ bên trong vật nung khi hệ số truyền tĩnh $k(v)$ mang tính phi tuyến và lò có trễ thuần lớn?
  2. RQ2: Cấu trúc thuật toán số nào cho phép giải tối ưu phiếm hàm mục tiêu nung chính xác nhất dưới các ràng buộc cơ lý nghiêm ngặt mà không rơi vào bế tắc nghiệm suy biến kiểu "bang-bang" của nhiệt độ lò?
  3. Hypothesis H1: Việc thay thế khâu trễ $e^{-\tau s}$ bằng phép xấp xỉ Padé bậc một trong miền $6 \le T/\tau < 10$ kết hợp với tuyến tính hóa từng đoạn hệ số $k(v)$ và biến đổi Laplace ngược sẽ triệt tiêu sai số mô hình hóa tốt hơn đáng kể so với khai triển Taylor truyền thống.
  4. Hypothesis H2: Chuyển đổi phiếm hàm mục tiêu tích phân liên tục về bài toán tối ưu hàm nhiều biến thông qua công thức tích phân số Simpson/hình thang sẽ xác định được quỹ đạo điện áp tối ưu $u^*(t)$ khả thi tuyệt đối trên phần cứng biến đổi điện trở thực tế.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa của Lý thuyết hệ tham số phân bố (Butkovskii, 1969; Wang & Tung, 1964), Phương trình truyền nhiệt Fourier dạng parabolic với điều kiện biên loại 3 (Robin boundary condition), và Lý thuyết xấp xỉ hàm hữu tỉ Padé.

Đóng góp đột phá của công trình được định lượng thông qua việc giảm thiểu sai số trường nhiệt mục tiêu $J_c$ xuống mức tối thiểu (sai số nhiệt độ bề mặt và tâm vật ở cuối chu kỳ nung đạt dưới $1 - 2^\circ\text{C}$ trên dải nhiệt độ khảo sát $300^\circ\text{C} - 700^\circ\text{C}$), rút ngắn thời gian tính toán mô phỏng và đảm bảo loại bỏ hoàn toàn hiện tượng nứt vỡ ứng suất nhiệt. Phạm vi nghiên cứu tập trung vào vật nung dạng tấm phẳng dày tiêu chuẩn ($\text{Bi} \ge 0.5$) gồm hai loại vật liệu chịu lửa và cách nhiệt đặc trưng là Samot và Diatomite, với chu kỳ nung kéo dài từ $3300\text{s}$ đến $4500\text{s}$.

Literature Review và Positioning

Lịch sử phát triển lý thuyết điều khiển hệ tham số phân bố khởi nguồn từ công trình nền tảng của Butkovskii & Lerner (1960) và Butkovskii (1969) khi thiết lập Nguyên lý cực đại cho hệ mô tả bởi phương trình tích phân phi tuyến. Tiếp đó, Wang & Tung (1964) đã hoàn thiện mô tả không gian trạng thái của hệ DPS qua phương trình đạo hàm riêng (PDEs), phân tích tính điều khiển được (controllability) và quan sát được (observability). Trong cùng giai đoạn, Sakawa (1964) áp dụng phép tính biến phân đưa bài toán điều khiển nhiệt độ về phương trình tích phân Fredholm loại một.

Các tranh luận khoa học cốt lõi trong y văn quốc tế xoay quanh hai luồng quan điểm đối nghịch về phương pháp luận tổng hợp bộ điều khiển cho hệ DPS:

  • Trường phái "Thiết kế sau đó xấp xỉ" (Design-then-approximate): Đại diện bởi Curtain & Pritchard (1974, 1995), Fattorini (1993), Gibson (1979, 1981) với việc xây dựng toán tử Riccati trong không gian Hilbert vô hạn chiều và nửa nhóm (semigroup). Ưu điểm là tính chặt chẽ toán học tuyệt đối nhưng nhược điểm lớn là cực kỳ khó áp dụng vào các đối tượng phi tuyến có ràng buộc trạng thái thực tế do hiện tượng tràn chiều sai số (spillover effect) như Yoshida (1989) cảnh báo.
  • Trường phái "Xấp xỉ sau đó thiết kế" (Approximate-then-design): Phát triển bởi Banks & Kunisch (1984), Axelband (1966), Sage & Chaudhuri (1967), Padhi & Balakrishnan (2005) sử dụng các phép chiếu Galerkin, phân tích trực giao thích hợp (Proper Orthogonal Decomposition - POD) hoặc sai phân hữu hạn (FDM) để rút gọn hệ PDE về hệ phương trình vi phân thường (ODEs).
graph TD
    A["Hệ Tham Số Phân Bố (DPS)"] --> B["Design-then-Approximate<br/>(Curtain, Fattorini, Gibson)"]
    A --> C["Approximate-then-Design<br/>(Banks, Padhi, Sirovich)"]
    B --> D["Toán tử vô hạn chiều / Riccati<br/>Nhược điểm: Tràn chiều sai số"]
    C --> E["Rút gọn Galerkin / POD / FDM<br/>Nhược điểm: Mất thông tin biên liên tục"]
    A --> F["Phương Pháp Luận Án (Mai Trung Thái)<br/>Laplace + Padé 1 + Tuyến tính hóa từng đoạn k(v) + Tối ưu số"]

Tại Việt Nam, các nghiên cứu của Nguyễn Hữu Công (2006, 2010) đã bước đầu giải quyết bài toán nung tối ưu cho phôi cán trong lò nung đốt bằng dầu nặng FO. Tuy nhiên, công trình của Nguyễn Hữu Công (2010) chỉ dừng lại ở mô hình trễ nhỏ ($T/\tau \ge 10$) sử dụng xấp xỉ Taylor cấp 1 và coi hệ số truyền tĩnh $k$, hệ số dẫn nhiệt $\lambda$, hệ số khuếch tán nhiệt độ $a$ là các hằng số không đổi.

Vị trí của luận án được định vị chính xác tại điểm nghẽn của các nghiên cứu tiền nhiệm:

  1. So sánh với nghiên cứu của Subasi (2006) và Qian Li et al. (2012): Các tác giả quốc tế nghiên cứu điều khiển điểm tối ưu cho hệ parabolic nhưng coi cơ cấu chấp hành là tuyến tính tức thời. Luận án tích hợp trễ vận tải lớn ($\tau$) và động học phi tuyến của lò điện vào phương trình trạng thái biên.
  2. So sánh với Behzad Talaei et al. (2014): Nghiên cứu của Talaei dùng mạng nơ-ron sai phân xấp xỉ quy hoạch động (ADP) đòi hỏi tài nguyên tính toán rất lớn và khó cài đặt trên vi điều khiển công nghiệp; trong khi đó, giải pháp giải tích số kết hợp biến đổi Laplace của luận án cho phép tính toán trước quỹ đạo điện áp điều khiển ngoại tuyến (offline trajectory) với độ tin cậy và tính ổn định tuyệt đối.

Luận án đã lấp đầy khoảng trống bằng cách kết hợp xấp xỉ Padé bậc một cho khâu trễ lớn ($6 \le T/\tau < 10$), phân đoạn hóa phi tuyến $k(v)$ theo nhiệt độ lò và giải tích hóa trường nhiệt $q(x,t)$ bằng biến đổi Laplace ngược, chuyển phiếm hàm tối ưu liên tục thành bài toán quy hoạch toán học nhiều biến giải bằng phương pháp số.

Đó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à bổ sung sâu sắc cho Lý thuyết Điều khiển Tối ưu Hệ tham số phân bố (Butkovskii, 1969) và Lý thuyết Truyền nhiệt phi tuyến thông qua các đóng góp lý thuyết cốt lõi:

Thứ nhất, luận án thiết lập thành công mô hình toán học giải tích tường minh liên kết giữa tín hiệu đầu vào điện áp $u(t)$ của bộ biến đổi nguồn và trường nhiệt độ phân bố $q(x,t)$ của tấm vật liệu dày. Mô hình vượt qua giới hạn của phương pháp cổ điển bằng cách mô tả chính xác sự phụ thuộc phi tuyến của hệ số truyền tĩnh lò $k(v)$ theo từng dải nhiệt độ: $$T \frac{dv(t)}{dt} + v(t) = k(v) u(t-\tau)$$ với $k(v)$ được tuyến tính hóa từng đoạn thành tập hợp các hằng số ${k_1, k_2, k_3, \dots, k_N}$ tương ứng với các khoảng thời gian ${t_1, t_2, t_3, \dots, t_N}$.

Thứ hai, luận án giải quyết triệt để sự mâu thuẫn vật lý trong các công trình kinh điển của Pontryagin khi áp dụng cho quá trình nhiệt. Trong các nghiên cứu trước (như [29], [46], [72]), tín hiệu điều khiển tối ưu tối ưu tìm được trực tiếp là nhiệt độ lò $v^(t)$ mang dạng đóng mở "bang-bang" (xung vuông nhảy cấp). Điều này bất khả thi trong thực tế vì quán tính nhiệt của lò rất lớn. Bằng cách bổ sung khâu động học lò điện trở vào mạch điều khiển, luận án chuyển hàm điều khiển tối ưu về điện áp cung cấp $u^(t)$ hoặc công suất $p^*(t)$. Tín hiệu điện áp nhảy cấp hoàn toàn khả thi do các linh kiện bán dẫn công suất (Thyristor) có quán tính điện vô cùng nhỏ (tính bằng micro giây).

graph LR
    subgraph Classical ["Giải Pháp Kinh Điển"]
        Jc1["Hàm Mục Tiêu"] --> Pontryagin["Nguyên Lý Cực Đại"]
        Pontryagin --> vt_bang["v*(t) Dạng Bang-Bang (Xung Vuông)"]
        vt_bang --> Fail["BẤT KHẢ THI<br/>Quán tính nhiệt lò quá lớn"]
    end
    subgraph Proposed ["Giải Pháp Luận Án"]
        Jc2["Hàm Mục Tiêu"] --> NumLaplace["Số Hóa + Biến Đổi Laplace"]
        NumLaplace --> ut_opt["u*(t) Điện Áp Tối Ưu"]
        ut_opt --> Thyristor["Điều Khiển Góc Mở Thyristor α"]
        Thyristor --> Success["KHẢ THI VẬT LÝ<br/>Quán tính bán dẫn cực nhỏ"]
    end

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết toán - kỹ thuật:

  1. Lý thuyết biến đổi tích phân Laplace cho PDE Parabolic: Đưa phương trình vi phân đạo hàm riêng bậc hai $\frac{\partial^2 q(x,t)}{\partial x^2} = \frac{1}{a} \frac{\partial q(x,t)}{\partial t}$ về phương trình vi phân thường theo toán tử Laplace $s$ trong miền không gian $x$, cho phép biểu diễn nghiệm trường nhiệt $Q(x,s)$ dưới dạng giải tích giải được qua các hàm hypebolic ($\cosh, \sinh$).
  2. Lý thuyết xấp xỉ hữu tỉ Padé bậc một: Thay thế khâu trễ thuần $e^{-\tau s} \approx \frac{2 - \tau s}{2 + \tau s}$ cho các đối tượng có độ trễ tương đối lớn ($6 \le T/\tau < 10$), khắc phục hoàn toàn sự mất ổn định pha và sai số biên độ lớn của phép xấp xỉ Taylor ($e^{-\tau s} \approx \frac{1}{1 + \tau s}$).
  3. Quy hoạch phi tuyến và Tối ưu hóa số: Chuyển đổi phiếm hàm mục tiêu bình phương tích phân: $$J_c = \int_0^L \left[ q^*(x) - q(x, t_f) \right]^2 dx \to \min$$ thành hàm mục tiêu đại số nhiều biến $J_c(u_1, u_2, \dots, u_m)$ thông qua rời rạc hóa không gian $n$ lớp và thời gian $m$ bước bằng công thức Simpson và hình thang.

Điều kiện biên (Boundary conditions) được xác lập nghiêm ngặt:

  • Biên truyền nhiệt mặt trước ($x = 0$): $-\lambda \left.\frac{\partial q(x,t)}{\partial x}\right|_{x=0} = \alpha [v(t) - q(0,t)]$ (Điều kiện biên loại 3).
  • Biên cách nhiệt mặt sau ($x = L$): $\left.\frac{\partial q(x,t)}{\partial x}\right|_{x=L} = 0$ (Điều kiện biên loại 2, thông lượng nhiệt triệt tiêu).
  • Ràng buộc điện áp: $U_1 \le u(t) \le U_2$.
  • Ràng buộc nhiệt độ môi trường lò: $v(t) \le U_3$ (bảo vệ vật liệu xây lò).
  • Ràng buộc nhiệt độ bề mặt vật nung: $q(0,t) \le U_4$ (chống cháy cục bộ, biến tính hóa học).
  • Ràng buộc gradient nhiệt độ: $\left| \frac{\partial q}{\partial x} \right| \le U_5$ (triệt tiêu ứng suất nhiệt gây nứt vỡ vật liệu).

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

Thiết kế nghiên cứu

Luận án lựa chọn lập trường triết học khoa học thực chứng (Positivism) kết hợp với chủ nghĩa duy thực phản biện (Critical Realism). Thiết kế nghiên cứu áp dụng phương pháp hỗn hợp định lượng (Quantitative Mixed-Design) bao gồm: Mô hình hóa toán học giải tích $\to$ Thuật toán tối ưu hóa số $\to$ Mô phỏng số trên môi trường chuyên dụng $\to$ Thực nghiệm kiểm chứng trên mô hình vật lý thời gian thực.

Thiết kế đa tầng (Multi-level design) được cấu trúc chặt chẽ qua 3 cấp độ:

  • Cấp độ 1 (Mô hình vật lý lò nhiệt): Nhận dạng thực nghiệm tham số động học lò điện trở, xác định hằng số thời gian $T$, thời gian trễ $\tau$ và hàm truyền đạt phi tuyến $k(v)$.
  • Cấp độ 2 (Mô hình dẫn nhiệt 1D): Rời rạc hóa vật thể nung theo chiều dày $L$ thành $n$ phân lớp hữu hạn ($n = 6, 8, 10$).
  • Cấp độ 3 (Bộ điều khiển tối ưu): Tối ưu hóa chuỗi điện áp điều khiển rời rạc ${u_j}$ trên $m$ khoảng thời gian ($m = 40, 64, 100$).

Quy trình nghiên cứu rigorous

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

  1. Giao thức nhận dạng tham số lò: Cung cấp điện áp bước nhảy $u = 220\text{V}$, ghi nhận đường cong đáp ứng quá độ nhiệt độ $v(t)$ thông qua cảm biến can nhiệt K kết nối card thu thập dữ liệu NI USB-6008 của hãng National Instruments. Áp dụng phương pháp Ziegler-Nichols và đồ thị tiếp tuyến xác định điểm uốn để trích xuất $T$ và $\tau$.
  2. Kỹ thuật tam giác đạc (Triangulation): Đối chiếu chéo kết quả trường nhiệt giữa (i) Công thức giải tích biến đổi Laplace, (ii) Kết quả mô phỏng số trên MATLAB & Simulink, và (iii) Dữ liệu đo đạc thực tế từ các đầu đo nhiệt gắn tại bề mặt ($x=0$), tâm ($x=L/2$) và đáy ($x=L$) của mẫu vật nung.
  3. Kiểm soát độ tin cậy và độ giá trị (Validity & Reliability): Hiệu chuẩn các kênh đo nhiệt độ bằng thiết bị chuẩn, bảo đảm sai số đo của can nhiệt và mạch khuếch đại không vượt quá $\pm 0.5^\circ\text{C}$. Độ ổn định của nghiệm tối ưu được kiểm tra tính nhạy (sensitivity analysis) khi thay đổi bước chia lưới không gian và thời gian.
sequenceDiagram
    participant ID as "Nhận Dạng Hệ Thống"
    participant Math as "Mô Hình Toán Laplace + Padé"
    participant Sim as "Mô Phỏng MATLAB/Simulink"
    participant Hardware as "Thực Nghiệm Lò Thật & NI USB-6008"

    ID->>Math: Tham số T, τ, k(v)
    Math->>Sim: Quỹ đạo tối ưu u*(t), q(x,t)
    Sim->>Hardware: Áp dụng u*(t) qua bộ biến đổi Thyristor
    Hardware->>Sim: Phản hồi q_thực(x,t) & v_thực(t)
    Hardware-->>Math: Đối chiếu kiểm chuẩn tam giác đạc (Sai số < 1.8%)

Data và phân tích

Đặc tính mẫu vật liệu nghiên cứu:

  • Mẫu Samot (Gạch chịu lửa): Kích thước tấm phẳng dày $L = 0.065\text{m}$, khối lượng riêng $\rho = 2100\text{kg/m}^3$, nhiệt dung riêng $C = 1050\text{J/(kg}\cdot^\circ\text{C)}$, hệ số dẫn nhiệt $\lambda = 1.05\text{W/(m}\cdot^\circ\text{C)}$, hệ số khuếch tán nhiệt $a = 4.76 \times 10^{-7}\text{m}^2/\text{s}$, hệ số trao đổi nhiệt $\alpha = 45\text{W/(m}^2\cdot^\circ\text{C)}$, hệ số $\text{Bi} = \alpha L / \lambda = 2.785$ (vật dày, $\text{Bi} > 0.5$).
  • Mẫu Diatomite (Gạch cách nhiệt): Kích thước tấm phẳng dày $L = 0.05\text{m}$, khối lượng riêng $\rho = 600\text{kg/m}^3$, nhiệt dung riêng $C = 920\text{J/(kg}\cdot^\circ\text{C)}$, hệ số dẫn nhiệt $\lambda = 0.18\text{W/(m}\cdot^\circ\text{C)}$, hệ số khuếch tán nhiệt $a = 3.26 \times 10^{-7}\text{m}^2/\text{s}$, hệ số trao đổi nhiệt $\alpha = 35\text{W/(m}^2\cdot^\circ\text{C)}$, hệ số $\text{Bi} = 9.722$ (vật rất dày).

Kỹ thuật phân tích tính toán nâng cao:

  • Xây dựng chương trình chuyên dụng trên môi trường MATLAB giải hệ phương trình xác định các hàm ảnh $g_\mu(x,t)$ và $g_\mu(x, t-\tau)$ ($\mu=1,2,3$).
  • Thuật toán giải bài toán tối ưu hóa ràng buộc phi tuyến (Nonlinear Programming Solver) dựa trên hàm mục tiêu đại số cực tiểu hóa sai số $e = J_c$.
  • Kiểm tra tính bền vững (Robustness checks) bằng cách so sánh hiệu năng của nghiệm tối ưu khi xấp xỉ hàm trễ theo chuỗi Taylor so với khai triển Padé bậc một ở các độ phân giải lưới khác nhau: $n=6, m=40$; $n=8, m=64$; $n=10, m=100$.

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

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

  1. Hiệu quả vượt trội của phép xấp xỉ Padé bậc một so với Taylor: Kết quả mô phỏng trên mẫu Samot với cùng điều kiện thời gian nung $t_f = 4200\text{s}$ cho thấy, khi số bước chia lưới tăng từ $(n=6, m=40)$ lên $(n=10, m=100)$, sai số hàm mục tiêu $J_c$ theo xấp xỉ Taylor giảm từ $3.45$ xuống $2.82$, trong khi theo xấp xỉ Padé 1, sai số $J_c$ giảm mạnh từ $2.15$ xuống chỉ còn $0.84$. Khai triển Padé bậc một giúp nâng cao độ chính xác lên hơn $65%$ nhờ bảo toàn đáp ứng pha của khâu trễ lớn ($6 \le T/\tau < 10$).
  2. Quy luật biến thiên phi tuyến của hệ số truyền tĩnh $k(v)$: Nhận dạng thực nghiệm chứng minh rằng hệ số $k$ không phải là hằng số mà suy giảm phi tuyến theo nhiệt độ lò: trong dải nhiệt độ từ môi trường đến $200^\circ\text{C}$, $k_1 \approx 1.85$; dải $200^\circ\text{C} - 350^\circ\text{C}$, $k_2 \approx 1.62$; và dải $350^\circ\text{C} - 500^\circ\text{C}$, $k_3 \approx 1.38$. Việc phân đoạn $k$ thành 3 khoảng thời gian $t_1, t_2, t_3$ là điều kiện tiên quyết để triệt tiêu sai số mô hình.
  3. Tính chất phản trực giác của quỹ đạo điện áp tối ưu $u^*(t)$: Khác với trực giác vận hành thông thường là cấp điện áp tối đa ở giai đoạn đầu để nung nhanh, quỹ đạo tối ưu $u^(t)$ bắt buộc phải khởi đầu ở mức điện áp trung bình để khống chế gradient nhiệt độ $\partial q/\partial x \le U_5$, sau đó tăng vọt lên mức trần $U_2$ ở pha giữa và chủ động giảm áp sớm trước thời điểm kết thúc $t_f$ nhằm tận dụng quán tính nhiệt và trễ truyền nhiệt, đưa nhiệt độ toàn bộ các lớp vật nung hội tụ đồng đều về giá trị đặt $q^$.
  4. Kiểm chứng thực nghiệm chính xác cao: Trên hệ thống thí nghiệm thực tế với mẫu Samot nung đến $q^* = 300^\circ\text{C}$ ($t_f = 4200\text{s}$) và $q^* = 400^\circ\text{C}$ ($t_f = 4200\text{s}$), nhiệt độ bề mặt đo được đạt $301.2^\circ\text{C}$ và tâm vật đạt $298.6^\circ\text{C}$ (độ chênh nhiệt giữa các lớp chỉ $2.6^\circ\text{C}$). Với mẫu Diatomite ($q^* = 400^\circ\text{C}, t_f = 4500\text{s}$), nhiệt độ thực nghiệm bám sát đường cong lý thuyết với sai số cực đại không quá $1.8%$.
Vật liệu mẫu Thời gian nung $t_f$ (s) Nhiệt độ đặt $q^*$ ($^\circ\text{C}$) Sai số $J_c$ (Taylor) Sai số $J_c$ (Padé 1) Chênh lệch nhiệt độ bề mặt - đáy ($^\circ\text{C}$)
Samot ($L=0.065\text{m}$) $3300$ $300$ $4.12$ $1.45$ $4.2$
Samot ($L=0.065\text{m}$) $4200$ $300$ $2.82$ $0.84$ $2.6$
Diatomite ($L=0.05\text{m}$) $3600$ $400$ $3.95$ $1.28$ $5.1$
Diatomite ($L=0.05\text{m}$) $4500$ $400$ $2.41$ $0.72$ $3.0$

Implications đa chiều

  • Về mặt học thuật: Cung cấp phương pháp luận tổng quát hóa giải tích - số kết hợp, giải quyết thành công bài toán điều khiển tối ưu hệ DPS parabolic có đồng thời cả trễ vận tải và tham số phi tuyến mà không làm phức tạp hóa cấu trúc tính toán.
  • Về mặt phương pháp: Quy trình xấp xỉ Padé kết hợp biến đổi Laplace ngược và tích phân Simpson có thể chuyển giao trực tiếp để giải quyết các hệ phân bố khác như quá trình khuếch tán hóa học, truyền sóng trong môi trường đàn hồi hay phân bố nồng độ phản ứng sinh học.
  • Về mặt công nghệ và thực tiễn: Cho phép các nhà máy luyện kim, chế tạo gốm sứ và sản xuất vật liệu từ thiết kế các bộ điều khiển quá trình nhiệt hiệu suất cao, tiết kiệm từ $12% - 18%$ năng lượng tiêu thụ điện so với phương pháp điều khiển PID truyền thống, đồng thời triệt tiêu tỷ lệ phế phẩm nứt vỡ do sốc nhiệt.
  • Về chính sách và tiêu chuẩn: Định hình tiêu chuẩn tiết kiệm năng lượng và kiểm soát phát thải nhiệt trong công nghiệp nặng, thúc đẩy lộ trình tự động hóa các lò nung công nghiệp tại Việt Nam theo chuẩn công nghệ hiện đại.

Limitations và Future Research

Mặc dù đạt được những kết quả lý thuyết và thực nghiệm đột phá, công trình vẫn tồn tại một số giới hạn nghiên cứu:

  1. Giới hạn một chiều không gian: Luận án mới chỉ xét quá trình truyền nhiệt 1D theo phương $x$ cho tấm phẳng (chiều dày $L$ rất nhỏ so với chiều dài và chiều rộng). Trong thực tế, các chi tiết cơ khí phức tạp có truyền nhiệt 2D hoặc 3D với hình học phi đối xứng.
  2. Tuyến tính hóa tham số vật liệu: Luận án giả thiết hệ số dẫn nhiệt $\lambda$ và hệ số khuếch tán nhiệt độ $a$ của vật nung là hằng số dọc theo chiều dày các lớp mỏng; thực tế ở các dải nhiệt độ rất cao ($> 1000^\circ\text{C}$), tính phi tuyến của $\lambda(q)$ và nhiệt dung riêng $C(q)$ trở nên đáng kể.
  3. Cơ cấu nung tĩnh một phía: Đối tượng nghiên cứu là lò tĩnh truyền nhiệt một phía; chưa mở rộng cho hệ lò liên tục (continuous reheating furnace) có phôi chuyển động hoặc gia nhiệt bức xạ hai phía đối xứng.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng chiến lược:

  • Hướng 1: Mở rộng thuật toán cho hệ tham số phân bố đa chiều (2D/3D Heat Conduction PDEs) áp dụng cho phôi trụ tròn hoặc hình khối chữ nhật.
  • Hướng 2: Tích hợp tính phi tuyến đồng thời của cả lò và vật liệu ($\lambda(q), a(q)$) kết hợp các bậc xấp xỉ Padé cao hơn (Padé bậc 2, Padé bậc 3 khi tỉ số $1 \le T/\tau < 6$).
  • Hướng 3: Phát triển cấu trúc điều khiển vòng kín thích nghi (Closed-loop Adaptive Control) với cơ chế bù nhiễu thời gian thực khi điện áp lưới hoặc nhiệt dung lò thay đổi đột ngột.
  • Hướng 4: Cài đặt thuật toán tối ưu hóa trực tiếp lên các chip điều khiển nhúng DSP/FPGA công nghiệp để tính toán tối ưu thích nghi trực tuyến (Online Real-time Optimization).

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

Công trình của NCS. Mai Trung Thái tạo ra tác động học thuật sâu sắc trong cộng đồng Tự động hóa và Khoa học Điều khiển tại Việt Nam, đóng vai trò tài liệu tham khảo nền tảng cho các chương trình đào tạo sau đại học về điều khiển hệ tham số phân bố. Phương pháp luận của luận án dự báo sẽ thu hút trích dẫn quan trọng từ các nhóm nghiên cứu về tối ưu hóa quá trình nhiệt công nghiệp.

graph TD
    Thesis["Luận Án Tiến Sĩ (Mai Trung Thái)"] --> Acad["Học Thuật & Đào Tạo<br/>Giáo trình Sau đại học, Bài báo quốc tế"]
    Thesis --> Ind["Chuyển Giao Công Nghiệp<br/>Ủ thép cán, Nhiệt luyện chi tiết máy, Gốm sứ"]
    Thesis --> Energy["Môi Trường & Năng Lượng<br/>Giảm 12-18% điện năng, Giảm phát thải CO2"]
    Thesis --> Policy["Tiêu Chuẩn Hóa<br/>Định mức tiêu hao năng lượng công nghiệp"]

Về mặt công nghiệp, giải pháp điều khiển tối ưu đã chứng minh khả năng ứng dụng thực tế vượt trội trên các hệ thống lò điện trở nhiệt luyện, nhà máy gạch ceramic và tôi ram thép. Khả năng kiểm soát chính xác gradient nhiệt độ giúp nâng cao cơ tính sản phẩm, tăng tuổi thọ chi tiết máy và giảm thiểu lượng tiêu hao điện năng quốc gia. Kết quả nghiên cứu có tính tương thích toàn cầu, dễ dàng áp dụng cho các hệ thống nhiệt công nghiệp tại bất kỳ quốc gia nào đang trong tiến trình hiện đại hóa sản xuất.

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

  • Nghiên cứu sinh và Học viên cao học: Tiếp cận khung phương pháp luận hoàn chỉnh từ mô hình hóa toán học vi tích phân, xấp xỉ Padé đối tượng trễ đến thuật toán số hóa giải tích biến đổi Laplace trên máy tính.
  • Các nhà khoa học và Giảng viên chuyên ngành Điều khiển: Sở hữu bộ công cụ lý thuyết chặt chẽ và kết quả đối chứng tin cậy để tiếp tục phát triển các bài toán điều khiển biên bền vững cho hệ PDEs.
  • Kỹ sư R&D và Chuyên gia Công nghệ Nhiệt: Ứng dụng ngay sơ đồ điều khiển góc mở Thyristor qua card thu thập dữ liệu (DAQ) và quy trình lập trình tối ưu hóa quỹ đạo nhiệt vào các dự án nâng cấp lò công nghiệp thực tế.
  • Các nhà quản lý năng lượng và Doanh nghiệp sản xuất: Có cơ sở khoa học để hoạch định dự án đầu tư tự động hóa quá trình gia nhiệt, tiết kiệm chi phí năng lượng và giảm phát thải carbon.

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 mở rộng Lý thuyết Điều khiển Tối ưu Hệ tham số phân bố (Butkovskii, 1969) bằng cách tích hợp đồng thời trễ vận tải lớn ($6 \le T/\tau < 10$) và tính phi tuyến của hệ số truyền tĩnh lò $k(v)$ vào bài toán truyền nhiệt biên loại 3. Luận án đã tìm ra lời giải giải tích - số tường minh cho trường nhiệt $q(x,t)$ theo điện áp cung cấp $u(t)$, khắc phục triệt để nghịch lý điều khiển dạng xung vuông bất khả thi của các công trình kinh điển.

2. Điểm mới về phương pháp luận so với ít nhất hai nghiên cứu tiền nhiệm?

So với nghiên cứu của Nguyễn Hữu Công (2010) (chỉ dùng xấp xỉ Taylor cho trễ nhỏ và giả định $k = \text{const}$) và nghiên cứu của Subasi (2006) (bỏ qua trễ và động học phi tuyến của cơ cấu chấp hành), luận án đã ứng dụng thành công khai triển Padé bậc một kết hợp phân đoạn tuyến tính hóa tham số $k(v)$ và số hóa tích phân Simpson. Điều này giúp giảm sai số hàm mục tiêu $J_c$ hơn $65%$ so với phương pháp Taylor ở cùng cấp lưới rời rạc.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ và có ý nghĩa thực tiễn cao nhất là tính phi tuyến suy giảm mạnh của $k(v)$ theo nhiệt độ (từ $1.85$ ở nhiệt độ thường xuống còn $1.38$ ở $500^\circ\text{C}$) và quy luật phân bố điện áp tối ưu $u^*(t)$. Để đạt độ đồng đều nhiệt cao nhất ở thời điểm kết thúc $t_f$, điện áp không được cấp cực đại ở đầu chu kỳ mà phải cấp tăng dần và ngắt sớm trước $t_f$ hàng trăm giây để dòng nhiệt tự cân bằng từ bề mặt vào tâm phôi.

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

Luận án cung cấp đầy đủ và chi tiết toàn bộ giao thức tái lập thực nghiệm:

  • Cấu hình phần cứng: Lò điện trở công suất thí nghiệm, mẫu thử Samot ($L=0.065\text{m}$) và Diatomite ($L=0.05\text{m}$), can nhiệt K hiệu chuẩn, mạch khuếch đại tín hiệu đo, card thu thập dữ liệu NI USB-6008, bộ biến đổi xoay chiều - xoay chiều 1 pha dùng 2 Thyristor đấu song song ngược và mạch tạo xung điều khiển pha.
  • Cấu hình phần mềm: Mã nguồn chương trình tính toán các hàm ảnh $g_\mu(x,t)$, thuật toán giải bài toán tối ưu phi tuyến trên MATLAB/Simulink và giao diện đồ họa giám sát thực nghiệm.

5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Lộ trình phát triển 10 năm tiếp theo bao gồm: (1) Hoàn thiện lý thuyết điều khiển tối ưu cho trường nhiệt 2D/3D không gian phi đối xứng; (2) Ứng dụng xấp xỉ Padé bậc 2 và bậc 3 cho các hệ thống nhiệt siêu trễ ($1 \le T/\tau < 6$); (3) Xây dựng hệ thống điều khiển vòng kín bền vững thích nghi trên nền tảng AI/Machine Learning; (4) Chế tạo thiết bị điều khiển nhúng công nghiệp chuyên dụng (Dedicated Industrial Controller) thương mại hóa cho ngành luyện kim và nhiệt luyện.

Kết luận

  1. Luận án đã giải quyết trọn vẹn và xuất sắc bài toán điều khiển tối ưu theo tiêu chuẩn "nung chính xác nhất" cho hệ lò điện trở - vật nung dạng tấm phẳng, đại diện tiêu biểu cho lớp hệ thống với tham số phân bố, có trễ lớn và phi tuyến.
  2. Thiết lập thành công mô hình toán học giải tích tường minh xác định trường phân bố nhiệt độ $q(x,t)$ theo điện áp điều khiển $u(t)$, tích hợp khâu quán tính có trễ với hệ số truyền tĩnh phi tuyến $k(v)$ thay đổi theo nhiệt độ lò.
  3. Ứng dụng đột phá phép xấp xỉ Padé bậc một cho khâu trễ $e^{-\tau s}$ trong miền tham số $6 \le T/\tau < 10$, chứng minh tính vượt trội về độ chính xác và tính ổn định so với phương pháp xấp xỉ chuỗi Taylor truyền thống.
  4. Xây dựng thuật toán số hóa tối ưu chuyển đổi phiếm hàm mục tiêu tích phân liên tục về bài toán quy hoạch toán học nhiều biến, cho phép tính toán chính xác quỹ đạo điện áp tối ưu $u^*(t)$ đáp ứng đầy đủ các ràng buộc công nghệ nhiệt phức tạp.
  5. Kiểm chứng toàn diện và nhất quán giữa lý thuyết, mô phỏng số trên MATLAB và thực nghiệm vật lý trên hai loại mẫu vật liệu Samot và Diatomite qua hệ thống ghép nối NI USB-6008, với sai số bám trường nhiệt thực tế dưới $1.8%$.
  6. Công trình mở ra bước tiến mới trong kỹ thuật điều khiển các quá trình công nghệ nhiệt, đóng góp to lớn vào kho tàng lý thuyết điều khiển tự động và mang lại giá trị gia tăng cao cho nền sản xuất công nghiệp hiện đại.