Tổng quan về luận án

Nghiên cứu về tự động hóa quá trình nhiệt luyện công nghiệp đóng vai trò then chốt trong việc tối ưu hóa hiệu suất năng lượng và nâng cao chất lượng cơ tính sản phẩm thép. Luận án tiến sĩ kỹ thuật với đề tài "Phương pháp mô hình hóa và điều khiển trường nhiệt độ vật nung dày" do nghiên cứu sinh Nguyễn Việt Dũng thực hiện dưới sự hướng dẫn khoa học của PGS. Hoàng Minh Sơn và TS. Nguyễn Văn Hòa tại Đại học Bách khoa Hà Nội (2018), thuộc chuyên ngành Kỹ thuật Điều khiển và Tự động hóa, là một công trình tiên phong giải quyết bài toán điều khiển phi tuyến phân bố rải cho hệ thống gia nhiệt phôi kim loại.

Trong bối cảnh thực tiễn của ngành luyện kim Việt Nam, theo báo cáo của Bộ Công Thương và Tổng công ty Thép Việt Nam, cả nước có gần 70 dự án luyện cán thép công suất trên 100.000 tấn/năm nhưng chỉ vận hành dưới 50% công suất thiết kế, tiêu thụ tới 3,5 tỷ kWh điện mỗi năm. Do công nghệ điều khiển nhiệt truyền thống còn nhiều hạn chế: "thời gian nung, chi phí luyện một mẻ phôi thép cao gần gấp đôi so với trung bình trên thế giới (từ 90-180 phút, tiêu thụ từ 550-690 kWh/ tấn so với trung bình thế giới vào khoảng 360-430 kWh)". Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các hệ thống hiện hữu chủ yếu điều khiển vòng hở hoặc điều khiển bám nhiệt độ bề mặt và nhiệt độ khí lò, trong khi nhiệt độ vùng lõi (tâm phôi) – yếu tố quyết định trực tiếp tới ứng suất nhiệt, sự hình thành lớp xỉ oxy hóa và độ biến dạng dẻo – lại không thể đo trực tiếp bằng cảm biến công nghiệp trong môi trường khắc nghiệt (>1000°C, rung lắc, bụi bẩn).

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  1. RQ1: Làm thế nào để quy đổi chính xác phương trình đạo hàm riêng (PDE) phi tuyến mô tả quá trình truyền nhiệt phân bố không gian với các tham số nhiệt vật lý phụ thuộc nhiệt độ về mô hình không gian trạng thái hữu hạn chiều mà không làm bùng nổ số chiều của hệ thống?
  2. RQ2: Cấu trúc điều khiển tối ưu nào có khả năng ước lượng và điều khiển nhiệt độ tâm phôi bám chính xác theo quỹ đạo công nghệ đặt trước trong điều kiện có ràng buộc về công suất gia nhiệt, tồn tại nhiễu đo và sai lệch mô hình?
  • Giả thuyết H1: Mô hình xấp xỉ không gian Galerkin kết hợp hệ suy luận mờ Mamdani (Galerkin-Mờ) có khả năng chuyển hóa phương trình PDE truyền nhiệt 1D/2D thành hệ phương trình vi phân thường (ODE) song tuyến bậc thấp với độ chính xác tương đương các mô hình phần tử hữu hạn đa lớp nhưng giảm thiểu đáng kể chi phí tính toán thời gian thực.
  • Giả thuyết H2: Thuật toán điều khiển tối ưu toàn phương tuyến tính (LQR) có tham số biến đổi trượt dọc theo trục thời gian (moving-horizon LQR) kết hợp cấu trúc điều khiển mô hình nội (IMC) và khâu bù bất định sẽ đảm bảo tính ổn định tiệm cận của hệ thống kín, triệt tiêu sai lệch bám nhiệt độ tâm vật nung trong phạm vi sai số dưới 10-15%.

Khung lý thuyết của công trình tích hợp lý thuyết phương trình đạo hàm riêng parabolic, giải tích hàm trên không gian Sobolev $W^{k,p}$, xấp xỉ không gian biến phân Ritz-Galerkin, lý thuyết tập mờ và điều khiển tối ưu tối thiểu hóa phiếm hàm mục tiêu toàn phương với phương trình vi phân đại số Riccati. Phạm vi nghiên cứu tập trung vào lớp vật nung dày có chuẩn số Biot $Bi \ge 0,5$, bề mặt gia nhiệt phân bố đều trên trục không gian $x, y$ và truyền nhiệt chủ đạo theo phương chiều dày $z$.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về mô hình hóa quá trình nung phôi kim loại phản ánh sự phân nhánh rõ rệt giữa hai trường phái tiếp cận chính:

  1. Trường phái mô hình hóa thực nghiệm chia lớp đẳng nhiệt (Lumped Parameter Isothermal Models): Điển hình là nghiên cứu của Pedersen và Wittenmark (1998, 2003) cùng các công trình của Depree et al. Mô hình này chia vật nung thành $N$ lớp lát cắt đẳng nhiệt dọc theo chiều dày và thiết lập hệ phương trình vi phân thường (ODE) bậc $N$. Mặc dù thuận tiện cho việc áp dụng các công cụ nhận dạng tham số offline thông qua thuật toán Levenberg-Marquardt hoặc giải thuật di truyền (GA), phương pháp này vấp phải mâu thuẫn cơ bản: để tăng độ mịn và độ chính xác của phân bố trường nhiệt độ, số lớp $N$ bắt buộc phải tăng lên, dẫn đến hiện tượng bùng nổ số chiều trạng thái (curse of dimensionality), gây quá tải cho bộ vi xử lý điều khiển thời gian thực và làm mất đi tính liên tục vật lý của phương trình truyền nhiệt.
  2. Trường phái mô hình hóa giải tích và số trị phân tán (Distributed Parameter PDE Systems): Đại diện bởi các nghiên cứu của W.H. Ray (1981), A.G. Butkovskiy (1982), và sau đó là Jaluria (2003). Các tác giả trực tiếp giải phương trình nhiệt Fourier phi tuyến thông qua phương pháp sai phân hữu hạn (FDM), phần tử hữu hạn (FEM) hoặc thể tích hữu hạn (FVM). Tuy nhiên, các phần mềm thương mại tính toán số như ANSYS, COMSOL hay các bộ giải PDE bậc cao chỉ đóng vai trò mô phỏng tĩnh hoặc giám sát (offline/open-loop prediction), không thể nhúng trực tiếp vào vòng lặp điều khiển phản hồi kín có chu kỳ quét (sampling time) tính bằng giây.

Vị trí học thuật của luận án được xác lập tại điểm giao thoa đột phá: chuyển đổi phương trình PDE parabolic phi tuyến về dạng mô hình không gian trạng thái song tuyến (bilinear state-space model) bậc 3 cố định nhờ phương pháp chiếu trực giao Galerkin, đồng thời khắc phục triệt để giả định phi thực tế của các nghiên cứu quốc tế trước đây (vốn coi hệ số dẫn nhiệt $\lambda$ và nhiệt dung riêng $c$ là hằng số) bằng cách tích hợp suy luận mờ Mamdani theo thời gian thực.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  • So với mô hình chia lớp của Pedersen & Wittenmark: Luận án duy trì bậc mô hình cố định ($n=3$) bất kể chiều dày vật nung, giảm hơn 60% số lượng biến trạng thái cần ước lượng nhưng vẫn đảm bảo tính toán liên tục trường nhiệt độ tại mọi tọa độ $z$.
  • So với mô hình xấp xỉ Galerkin thuần túy của các công trình luyện kim quốc tế (như các nghiên cứu được trích dẫn trong tài liệu [14, 16]): Luận án đã giải quyết được điểm nghẽn toán học khi thay thế phép tích phân không gian phi tuyến phức tạp $\bar{\lambda}_i(t), \bar{c}_i(t)$ bằng một khâu mờ hóa tĩnh Mamdani 1 ngõ vào - 2 ngõ ra cực kỳ gọn nhẹ, mở đường cho việc tổng hợp bộ điều khiển tối ưu trực tuyến.

Đó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 quan trọng vào lý thuyết điều khiển hệ tham số phân bố rải (Distributed Parameter Systems - DPS) thông qua 4 nội dung trụ cột:

                     [Phép chiếu tách biến Galerkin]
                      [Tích hợp Khâu Mờ Hóa Mamdani]
                     [Tối ưu hóa LQR trượt trục thời gian]
  1. Chứng minh tính đặt chỉnh (Well-posedness) của phương trình truyền nhiệt phi tuyến: Sử dụng công cụ giải tích hàm trong không gian Sobolev $W^{k,p}(\Omega)$, luận án chứng minh sự tồn tại, tính duy nhất và tính phụ thuộc liên tục vào điều kiện biên của nghiệm phương trình vi phân đạo hàm riêng truyền nhiệt phi tuyến có tham số biến thiên theo nhiệt độ.
  2. Xây dựng không gian trạng thái ảo 3 chiều giàu ý nghĩa vật lý: Thông qua phép chiếu trực giao với hệ hàm cơ sở $h_1(z) = 1$, $h_2(z) = \frac{2z}{Z}$, $h_3(z) = \frac{3}{2}\left(\frac{2z}{Z}\right)^2 - \frac{1}{2}$, nghiệm nhiệt độ được xấp xỉ dưới dạng tổng tuyến tính: $$T(z, t) = \sum_{i=1}^{3} h_i(z) x_i(t)$$ Các biến trạng thái được sáng tỏ về mặt bản chất vật lý:
    • $x_1(t)$: Đại diện cho nhiệt độ trung bình toàn thể tích vật nung.
    • $x_2(t)$: Thể hiện mức độ bất đối xứng của biểu đồ nhiệt độ giữa mặt trên và mặt dưới phôi.
    • $x_3(t)$: Thể hiện độ không đồng đều nhiệt độ tức thời giữa bề mặt ngoài và tâm vật nung.
  3. Mở rộng lý thuyết điều khiển tối ưu LQR cho hệ song tuyến phi dừng: Thiết lập cấu trúc xấp xỉ tuyến tính dừng cục bộ $H_k$ trượt dọc theo trục thời gian (moving-horizon linear time-invariant approximations), cho phép giải phương trình đại số Riccati tại từng chu kỳ lấy mẫu để xác định ma trận độ lợi phản hồi $K_k$, khắc phục hạn chế của lý thuyết LQR kinh điển vốn chỉ áp dụng cho hệ tuyến tính dừng (LTI).
  4. Đồng hóa cấu trúc điều khiển mô hình nội (IMC) với khâu lọc nhân quả: Tích hợp bộ quan sát trạng thái Galerkin-Mờ bên trong mạch vòng điều khiển kín, kết hợp khâu đạo hàm - quán tính bậc nhất có hằng số thời gian $\tau$ để ước lượng và triệt tiêu toàn bộ sai lệch mô hình và nhiễu tải bất định.

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

Khung phân tích của luận án kết hợp 3 lý thuyết độc lập: lý thuyết truyền nhiệt nhiệt động lực học, lý thuyết suy luận mờ mờ hóa tham số thực nghiệm và lý thuyết tối ưu hóa toàn phương Riccati.

Cấu trúc toán học của mô hình không gian trạng thái song tuyến thu được: $$\dot{x}(t) = \begin{bmatrix} 0 & 0 & 0 \ 0 & 0 & 0 \ 0 & 0 & -\frac{12\lambda(T)}{\rho c(T) Z^2} \end{bmatrix} x(t) + \begin{bmatrix} \frac{1}{\rho c(T) Z} & \frac{1}{\rho c(T) Z} \ -\frac{2}{\rho c(T) Z} & \frac{2}{\rho c(T) Z} \ \frac{1}{\rho c(T) Z} & \frac{1}{\rho c(T) Z} \end{bmatrix} \begin{bmatrix} u_1(t) \ u_2(t) \end{bmatrix}$$ $$y(t) = T(z, t) = [h_1(z) \quad h_2(z) \quad h_3(z)] x(t)$$

Điều kiện biên và phạm vi ứng dụng (boundary conditions): Áp dụng cho vật nung hình học dạng tấm phẳng, đối xứng nhiệt hai mặt hoặc gia nhiệt đơn mặt, thông lượng nhiệt cấp $u_1(t), u_2(t) \in [u_{\min}, u_{\max}]$ bị chặn trên và chặn dưới bởi công suất mỏ đốt và giới hạn chịu nhiệt của bề mặt kim loại nhằm chống chảy xỉ.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng định lượng kết hợp diễn dịch toán học (deductive-positivist paradigm). Cấu trúc nghiên cứu phân tầng gồm 3 giai đoạn:

Giai đoạn Mục tiêu Phương pháp / Công cụ Đầu ra chính
Giai đoạn 1 Mô hình hóa lý thuyết Tách biến Galerkin + Không gian Sobolev Hệ trạng thái song tuyến bậc 3
Giai đoạn 2 Nhận dạng & Mờ hóa Hệ mờ Mamdani + Dữ liệu nhiệt vật lý thép Khâu ước lượng $\lambda(T)$ và $c(T)$
Giai đoạn 3 Tổng hợp điều khiển & Kiểm chứng LQR trượt trục thời gian + Matlab/Simscape Bộ điều khiển tối ưu chống sai lệch mô hình

Quy trình nghiên cứu rigorous

  • Quy trình mờ hóa Mamdani: Biến đầu vào nhiệt độ $T$ được chia thành 9 tập mờ $Tl_j$ ($j=1..9$) cho hệ số dẫn nhiệt $\lambda(T)$ và 12 tập mờ $Tc_i$ ($i=1..12$) cho nhiệt dung riêng $c(T)$. Giá trị đầu ra sử dụng các tập mờ dạng singleton: $ld_j \in [60, 57.8]$ W/(m·K) và $c_i \in [480, 680]$ J/(kg·K). Cơ chế suy luận max-min và giải mờ theo phương pháp trọng tâm (centroid method) đảm bảo tính trơn và khả vi của các tham số khi nhiệt độ thay đổi từ 20°C đến 1250°C.
  • Xác thực chéo (Triangulation): Kiểm chứng mô hình Galerkin-Mờ bằng cách so sánh độc lập với 3 đối chứng:
    1. Mô hình sai phân hữu hạn (FDM) đa lớp lưới mịn.
    2. Mô hình thực nghiệm chia lớp đẳng nhiệt cân bằng năng lượng dạng ODE.
    3. Mô hình vật lý đa miền xây dựng trên Toolbox Simscape của Matlab (mô phỏng chính xác hiện tượng bức xạ khí lò, truyền nhiệt đối lưu và dẫn nhiệt phi tuyến 3D).
  • Tính vững và triệt tiêu nhiễu: Kiểm tra độ nhạy của bộ điều khiển dưới tác động đồng thời của nhiễu tải đầu vào $d(t)$ (biến động nhiệt trị nhiên liệu) và nhiễu đo ngẫu nhiên $n(t)$ (nhiễu trắng phân bố chuẩn trên tín hiệu nhiệt độ bề mặt).

Data và phân tích

  • Tham số đối tượng vật nung tiêu chuẩn: Chiều dày tấm thép $Z = 0,2$ m; khối lượng riêng $\rho = 7850$ kg/m³; dải nhiệt độ khảo sát từ $T_0 = 20^\circ\text{C}$ (293 K) đến nhiệt độ ra lò mục tiêu $T_{\text{target}} = 1150^\circ\text{C}-1200^\circ\text{C}$ (1423 K - 1473 K).
  • Thiết lập thuật toán LQR trượt: Tại mỗi chu kỳ lấy mẫu $\Delta t$, hệ thống tuyến tính hóa cục bộ xung quanh điểm làm việc danh định: $$\dot{\xi}k(t) = A_k \xi_k(t) + B_k \mu_k(t)$$ Phiếm hàm mục tiêu toàn phương cần cực tiểu hóa: $$J_k = \int{0}^{\infty} \left( \xi_k^T(t) Q_k \xi_k(t) + \mu_k^T(t) R_k \mu_k(t) \right) dt$$ với các ma trận trọng số xác định dương $Q_k = \text{diag}(q_1, q_2, q_3) > 0$ và $R_k = \text{diag}(r_1, r_2) > 0$.
  • Xử lý ràng buộc tín hiệu điều khiển: Tín hiệu điều khiển $u(t) = u_s[k] + \mu_k(t)$ được giới hạn trong tập lồi $U = [u_{\min}, u_{\max}]$ thông qua thuật toán chiếu bão hòa có hiệu chỉnh ma trận Riccati trực tiếp theo sơ đồ lưu đồ thuật toán tối ưu.

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

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

  1. Sự tương thích hoàn hảo giữa mô hình Galerkin-Mờ và mô hình phần tử hữu hạn: Kết quả mô phỏng số chỉ ra rằng sai lệch nhiệt độ tại tâm vật nung $T(0, t)$ giữa mô hình Galerkin-Mờ bậc 3 và mô hình sai phân hữu hạn (FDM) lưới 20 nút là dưới 1,2% trong suốt quá trình nâng nhiệt từ 20°C lên 1200°C. Đặc biệt, thời gian thực thi thuật toán giảm hơn 8 lần so với phương pháp FDM.
  2. Khắc phục nghịch lý điều khiển vùng biến dạng dẻo (Critical Plastic Region at 500°C): Luận án chứng minh bằng thực nghiệm mô phỏng rằng giai đoạn nhiệt độ đạt xấp xỉ 500°C là vùng nhạy cảm nhất đối với ứng suất nhiệt. "Đối với vật nung dày là vật nung có chỉ số Bi >= 0,5... việc chọn chế độ nung bao gồm nhiệt độ ban đầu của lò và nhiệt độ ban đầu của vật nung có ý nghĩa quan trọng, đặc biệt trong giai đoạn nhiệt độ vật nung từ giá trị ban đầu đến nhiệt độ bắt đầu biến dạng dẻo của kim loại (khoảng 500°C) là vùng nhiệt độ đặc biệt quan trọng đối với thép". Thuật toán LQR trượt đề xuất đã tự động ghìm tốc độ nâng nhiệt bề mặt, giữ hiệu nhiệt độ mặt - tâm $\Delta T = |T_{\text{surface}} - T_{\text{core}}| \le 85^\circ\text{C}$, triệt tiêu hoàn toàn nguy cơ nứt nẻ phôi do ứng suất nhiệt cục bộ.
  3. Triệt tiêu sai lệch mô hình nhờ cấu trúc mô hình nội mở rộng: Khi áp đặt sai số cố tình giữa mô hình đối tượng Simscape và mô hình toán học Galerkin-Mờ lên tới 20% về tham số dẫn nhiệt $\lambda$, bộ điều khiển LQR kết hợp khâu bù đạo hàm - quán tính bậc nhất vẫn điều khiển nhiệt độ tâm vật nung bám sát quỹ đạo đặt với sai số xác lập $e_{ss} < 1,5^\circ\text{C}$.
  4. Khả năng kháng nhiễu kép vượt trội: Dưới tác động đồng thời của nhiễu xung dòng nhiệt đầu vào (biến thiên $\pm 15%$) và nhiễu đo cảm biến nhiệt độ bề mặt (nhiễu trắng biên độ $\pm 10^\circ\text{C}$), hệ thống duy trì độ ổn định tiệm cận vững, không xảy ra hiện tượng dao động cưỡng bức tín hiệu điều khiển van cấp nhiên liệu.

Implications đa chiều

  • Về mặt lý thuyết: Luận án mở ra hướng tiếp cận mới cho việc thiết kế bộ điều khiển phản hồi trạng thái cho lớp hệ thống tham số phân bố rải phi tuyến parabolic tổng quát mà không cần qua các phép biến đổi tọa độ không gian phức tạp.
  • Về mặt phương pháp luận: Khẳng định tính khả thi của việc tích hợp lý thuyết mờ (Fuzzy Logic) vào các phương pháp giải tích cổ điển (Galerkin/Ritz) để biến đổi các hệ phi tuyến tham số biến thiên phức tạp thành hệ song tuyến có cấu trúc rõ ràng.
  • Về mặt ứng dụng công nghiệp: Cung cấp giải pháp phần mềm điều khiển tối ưu có thể cài đặt trực tiếp vào hệ thống PLC/DCS hiện có của các nhà máy luyện kim mà không đòi hỏi đầu tư thay mới toàn bộ phần cứng cảm biến đắt tiền.
  • Về mặt chính sách và môi trường: Đóng góp giải pháp kỹ thuật trực tiếp nhằm hiện thực hóa mục tiêu giảm tiêu hao năng lượng từ 550-690 kWh/tấn xuống mức tiệm cận thế giới (360-430 kWh/tấn), giảm phát thải khí nhà kính và xỉ thải kim loại cho ngành công nghiệp nặng quốc gia.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions):

  1. Giả thiết phân bố nhiệt đồng nhất trên mặt phẳng $x, y$: Luận án giới hạn bài toán truyền nhiệt một chiều theo phương chiều dày $z$. Trong thực tế các lò nung phôi tấm siêu trường, sự mất mát nhiệt qua các cạnh bên ($x, y$) có thể tạo ra gradien nhiệt độ 2D/3D đáng kể.
  2. Bỏ qua động học chuyển pha vi mô của tinh thể thép: Mô hình mờ hóa Mamdani dựa trên các đường đặc tính nhiệt vật lý vĩ mô $\lambda(T)$ và $c(T)$, chưa tính đến hiệu ứng tỏa nhiệt ẩn (latent heat) trong các điểm chuyển pha tế vi austenit-ferit diễn ra tức thời.
  3. Giả định lượng nhiệt cấp bề mặt là tín hiệu điều khiển lý tưởng: Luận án giả định vòng điều khiển nhiệt độ khí lò và độ mở van nhiên liệu bên ngoài đã được tối ưu hóa hoàn hảo và đáp ứng tức thời thông lượng nhiệt $u(t)$.

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

  • Mở rộng mô hình Galerkin-Mờ sang không gian 2 chiều và 3 chiều cho các loại phôi đúc có tiết diện phức tạp (phôi tròn, phôi dầm chữ I).
  • Tích hợp kỹ thuật điều khiển dự báo mô hình phi tuyến (NMPC) trên nền tảng mô hình trạng thái song tuyến để xử lý tối ưu hơn nữa các ràng buộc đa biến đồng thời.
  • Tiến hành thực nghiệm trực tiếp trên lò nung công nghiệp quy mô lớn có trang bị hệ thống thu thập dữ liệu chuyên dụng PTM1-010 hoặc hộp ghi dữ liệu PhoenixTM chịu nhiệt độ cao.

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

  • Ảnh hưởng học thuật: Các bài báo khoa học xuất bản từ luận án trên các tạp chí chuyên ngành uy tín (Tạp chí Khoa học & Công nghệ các trường đại học kỹ thuật, Hội nghị toàn quốc về Điều khiển và Tự động hóa - VCCA) đã xác lập một hướng nghiên cứu điển hình về điều khiển hệ DPS tại Việt Nam, mở ra tiềm năng trích dẫn cao trong lĩnh vực điều khiển quá trình nhiệt luyện.
  • Chuyển dịch công nghệ công nghiệp: Luận án cung cấp thuật toán lõi giúp ngành thép và gia công áp lực nâng cao chất lượng cán phôi, giảm tỷ lệ nứt phế phẩm, kéo dài tuổi thọ trục cán và tiết kiệm ước tính 10-15% tổng năng lượng đốt lò.
  • Đóng góp phát triển bền vững: Nâng cao hiệu suất năng lượng của lò nung công nghiệp là đóng góp thiết thực cho cam kết giảm phát thải carbon quốc gia trong khuôn khổ chuyển đổi xanh ngành công nghiệp sản xuất vật liệu xây dựng.

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

  • Nghiên cứu sinh & Nhà khoa học điều khiển: Tiếp cận một phương pháp luận chuẩn xác về chuyển đổi mô hình PDE thành ODE song tuyến và kỹ thuật giải Riccati trượt trục thời gian.
  • Kỹ sư R&D & Thiết kế hệ thống tự động hóa nhà máy thép: Sở hữu thuật toán quan sát nhiệt độ tâm phôi không đo được và bộ điều khiển LQR tối ưu dễ lập trình trên các dòng PLC công nghiệp thế hệ mới (Siemens S7-1500, ABB AC800M).
  • Chủ doanh nghiệp luyện kim: Nâng cao năng lực cạnh tranh nhờ cắt giảm tiêu hao nhiên liệu/điện năng nung phôi và giảm chi phí mài mòn cơ khí trục cán.
  • Cơ quan quản lý năng lượng (Bộ Công Thương): Có thêm căn cứ khoa học và giải pháp kỹ thuật để ban hành các định mức tiêu hao năng lượng bắt buộc trong ngành sản xuất thép.

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 độc đáo nhất là việc thiết lập mô hình không gian trạng thái song tuyến bậc 3 Galerkin-Mờ từ phương trình đạo hàm riêng truyền nhiệt parabolic phi tuyến. Luận án đã mở rộng lý thuyết tách biến biến phân Ritz-Galerkin cổ điển bằng cách ghép nối với hệ suy luận mờ Mamdani, giải quyết triệt để bài toán phụ thuộc phi tuyến của hệ số dẫn nhiệt $\lambda(T)$ và nhiệt dung riêng $c(T)$ mà không làm tăng bậc của hệ phương trình vi phân trạng thái.

2. Tính đổi mới về mặt phương pháp luận khi so sánh với các công trình đi trước?

So với phương pháp chia lớp đẳng nhiệt kinh điển của Pedersen & Wittenmark (phải tăng số phương trình vi phân khi tăng độ dày phôi), phương pháp của luận án duy trì số chiều trạng thái bất biến ($n=3$). So với các phương pháp giải số FDM/FEM/FVM (chỉ dùng mô phỏng tĩnh do khối lượng tính toán quá lớn), phương pháp của luận án cho phép tính toán trực tuyến thời gian thực với chu kỳ trích mẫu micro-giây, tương thích hoàn toàn với cấu trúc điều khiển tối ưu LQR.

3. Phát hiện gây bất ngờ nhất có số liệu thực nghiệm/mô phỏng chứng minh?

Phát hiện về khả năng ổn định tiệm cận của bộ điều khiển LQR trượt ngay cả khi mô hình toán học bị sai lệch tham số tới 20% so với đối tượng vật lý Simscape. Biểu đồ mô phỏng cho thấy sai lệch nhiệt độ tâm phôi chỉ dao động trong biên độ cực hẹp $\pm 1,5^\circ\text{C}$ nhờ tác động bù sai lệch tức thời của cấu trúc điều khiển mô hình nội (IMC) kết hợp bộ lọc quán tính đạo hàm bậc nhất.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) không?

Có. Luận án trình bày chi tiết toàn bộ hệ tọa độ hàm cơ sở $h_i(z)$, bảng quy tắc mờ Mamdani với 9 tập mờ cho $\lambda(T)$ và 12 tập mờ cho $c(T)$, các ma trận trạng thái giải tích $A(T), B(T)$, lưu đồ thuật toán điều khiển LQR trượt giải phương trình Riccati từng bước tại Hình 3.7 và mã nguồn chương trình mô phỏng trên Matlab/Simulink trong phần Phụ lục.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Định hướng phát triển tập trung vào: (1) Mở rộng hệ điều khiển sang bài toán phân bố rải 3 chiều cho các hình học phức tạp; (2) Tích hợp học máy/mạng nơron hồi quy sâu để tự chỉnh định ma trận trọng số $Q, R$ của bộ điều khiển LQR theo thời gian thực; (3) Triển khai thử nghiệm pilot trên lò nung công nghiệp quy mô lớn và thương mại hóa gói phần mềm điều khiển tối ưu nung thép.

Kết luận

Luận án "Phương pháp mô hình hóa và điều khiển trường nhiệt độ vật nung dày" của tác giả Nguyễn Việt Dũng là một công trình nghiên cứu khoa học xuất sắc, kết hợp nhuần nhuyễn giữa toán học giải tích hàn lâm và kỹ thuật điều khiển ứng dụng:

  1. Chuẩn hóa phương pháp mô hình hóa Galerkin-Mờ: Đưa ra công thức chuyển đổi chuẩn xác từ phương trình PDE truyền nhiệt sang mô hình trạng thái song tuyến bậc 3, giữ nguyên ý nghĩa vật lý của các biến trạng thái.
  2. Giải quyết bài toán ước lượng trạng thái không đo được: Xây dựng thành công bộ quan sát nhiệt độ tâm và phân bố trường nhiệt độ theo chiều dày vật nung trong điều kiện môi trường công nghiệp khắc nghiệt.
  3. Phát triển thuật toán LQR tối ưu trượt dọc trục thời gian: Tối ưu hóa quỹ đạo nhiệt độ bám sát yêu cầu công nghệ, triệt tiêu ứng suất nhiệt phá hủy tại vùng biến dạng dẻo 500°C.
  4. Chứng minh tính vững và khả năng kháng sai lệch mô hình: Tích hợp thành công nguyên lý điều khiển mô hình nội (IMC) giúp hệ thống triệt tiêu sai số mô hình lên tới 20% và kháng nhiễu đo hiệu quả.
  5. Kiểm chứng đa tầng bằng công cụ mô phỏng vật lý Simscape: Đảm bảo độ tin cậy tuyệt đối của kết quả nghiên cứu trước khi triển khai ứng dụng thực tế.
  6. Mở ra 3 nhánh nghiên cứu chuyên sâu mới: Bao gồm điều khiển hệ DPS đa chiều 3D, điều khiển tối ưu hóa đa mục tiêu ràng buộc phi tuyến cho lò nung công nghiệp, và ứng dụng trí tuệ nhân tạo trong nhận dạng tham số nhiệt luyện thời gian thực.