Tổng quan luận án

Quá trình nung vật nung dày trong công nghệ sản xuất kim loại và thép công nghiệp là công đoạn tiêu tốn lượng năng lượng lớn. Tại Việt Nam, theo báo cáo của Bộ Công Thương năm 2012 và Tổng công ty Thép Việt Nam năm 2014, các nhà máy luyện kim cán thép tiêu thụ khoảng 3,5 tỷ kWh điện hàng năm dù chỉ hoạt động dưới 50% công suất thiết kế; thời gian nung kéo dài từ 90 đến 180 phút và tiêu hao năng lượng từ 550 đến 690 kWh/tấn (so với mức trung bình thế giới từ 360 đến 430 kWh/tấn). Trở ngại kỹ thuật lớn nhất trong việc tối ưu hóa quá trình này nằm ở việc đo lường trực tiếp nhiệt độ bên trong lòng vật nung là không khả thi trong điều kiện công nghiệp (nhiệt độ cao, bụi bẩn, rung lắc). Đồng thời, các thuộc tính nhiệt vật lý của vật liệu biến thiên phi tuyến theo nhiệt độ và chuyển pha, khiến hệ thống mang bản chất phi tuyến, tham số phân bố rải và quán tính lớn.

Khoảng trống nghiên cứu được luận án chỉ ra là: phần lớn các nghiên cứu đã công bố tiếp cận hệ lò - vật nung theo các trường hợp riêng lẻ, sử dụng mô hình hàm truyền đạt quán tính bậc nhất có trễ hoặc mô hình thực nghiệm chia lớp đẳng nhiệt dạng phương trình vi phân thường (ODE) bậc thấp. Các nghiên cứu sử dụng phương trình đạo hàm riêng (PDE) thường giả thiết các tham số nhiệt vật lý là hằng số. Về mặt điều khiển, hầu hết các hệ thống chỉ dừng lại ở việc điều khiển nhiệt độ lò nung hoặc nhiệt độ bề mặt vật nung, chưa giải quyết được bài toán điều khiển trực tiếp trường nhiệt độ bên trong lòng (tâm) vật nung theo quỹ đạo tối ưu.

Mục tiêu và nhiệm vụ nghiên cứu của luận án được xác định gồm:

  1. Nghiên cứu khả năng sử dụng mô hình truyền nhiệt dạng phương trình đạo hàm riêng (PDE) có tham số thay đổi theo nhiệt độ bên trong vật nung để chuyển đổi về mô hình trạng thái phục vụ thiết kế bộ điều khiển.
  2. Xây dựng cấu trúc và thuật toán điều khiển có khả năng quan sát và điều khiển trường nhiệt độ trong lòng vật nung bám theo quỹ đạo định trước, đảm bảo độ chênh lệch nhiệt độ giữa bề mặt và tâm nằm trong giới hạn cho phép và đáp ứng các ràng buộc công nghệ.
  3. Kiểm chứng tính khả thi và độ tin cậy của mô hình cùng thuật toán điều khiển đề xuất thông qua mô phỏng trên phần mềm Matlab và công cụ mô phỏng vật lý Simscape.

Đối tượng nghiên cứu của luận án là lò nung với vật nung dày đặt trong lò (vật nung có tiêu chuẩn Biot $Bi \ge 0{,}5$). Phạm vi nghiên cứu giới hạn ở bài toán vật nung dạng tấm phẳng được cấp nhiệt đồng đều theo hai trục không gian $x, y$ và truyền nhiệt theo chiều dày trục $z$; tín hiệu điều khiển là lượng nhiệt cấp trực tiếp vào một hoặc cả hai bề mặt trên và dưới của vật nung; nhiệt độ trong lòng vật nung được giả thiết không đo được trực tiếp mà phải ước lượng thông qua mô hình quan sát trạng thái.

Tổng quan tài liệu và vị trí của luận án

Luận án đã tổng hợp và phân loại các công trình nghiên cứu trong và ngoài nước theo hai nhóm chính: xây dựng mô hình quá trình truyền nhiệt và thiết kế hệ thống điều khiển nung.

Nhóm hướng nghiên cứu Các công trình tiêu biểu được trích dẫn Đặc điểm và cách tiếp cận
Mô hình thực nghiệm chia lớp Pedersen & Wittenmark [53], [10, 23, 25, 27, 29, 34, 35, 44, 48] Chia vật nung thành $n$ lớp đẳng nhiệt, mô tả bằng hệ ODE phi tuyến; tham số được nhận dạng qua dữ liệu thực nghiệm đo đạc từ các cặp nhiệt điện gắn trên mẫu thử.
Phương pháp giải tích và số giải PDE [14, 15, 16, 17, 18, 19, 20], [30, 36, 45, 63, 74] Sử dụng 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) để tính toán trường nhiệt; một số nghiên cứu xấp xỉ PDE thành ODE nhưng giả thiết tham số dẫn nhiệt là hằng số.
Nhận dạng và tối ưu hóa tham số [3, 13, 32, 38, 41, 49, 51, 70, 76, 77] Ứng dụng quy hoạch phi tuyến, mạng nơ-ron hồi quy với thuật toán Levenberg-Marquardt, giải thuật di truyền (GA) và tối ưu bầy đàn (PSO) để tìm tham số mô hình.
Hệ thống và thuật toán điều khiển [9, 15, 18, 23, 25, 27, 28, 34, 44, 48, 53, 65, 77, 78] Phân chia thành điều khiển hở và điều khiển vòng kín; sử dụng các thuật toán kinh điển như PI, PID, điều khiển dự báo mô hình (MPC) và điều khiển tối ưu để ổn định nhiệt độ vùng lò hoặc bề mặt phôi.

Khoảng trống nghiên cứu được luận án lựa chọn giải quyết là: khắc phục nhược điểm tăng bậc mô hình khi chia nhiều lớp đẳng nhiệt của phương pháp thực nghiệm, đồng thời xử lý tính phi tuyến của các tham số nhiệt dung riêng $c(T)$ và hệ số dẫn nhiệt $\lambda(T)$ trong mô hình lý thuyết PDE. Luận án đề xuất phương pháp chuyển đổi mô hình PDE sang mô hình không gian trạng thái bậc cố định thông qua công thức tách biến Galerkin kết hợp suy luận mờ, từ đó thiết kế bộ điều khiển phản hồi trạng thái trượt dọc trục thời gian có xét đến ràng buộc và nhiễu.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án xây dựng trên các nền tảng lý thuyết và công cụ phân tích sau:

  • Định luật nhiệt động lực học và phương trình vi phân truyền nhiệt: Sử dụng định luật thứ nhất về nhiệt động lực học, định luật bảo toàn năng lượng và định lý phân kỳ Gauss (Ostrogradsky) để thiết lập phương trình vi phân đạo hàm riêng parabolic mô tả sự biến thiên nội năng và phân bố nhiệt độ theo không gian và thời gian: $$\rho c \frac{\partial T}{\partial t} = \nabla \cdot (\lambda \nabla T)$$
  • Lý thuyết không gian hàm Sobolev và tính đặt chỉnh (well-posedness): Sử dụng không gian Sobolev $W^{k,p}$ và không gian Hilbert $H^1$ để phân tích tính tồn tại, duy nhất và liên tục theo điều kiện biên của nghiệm phương trình truyền nhiệt.
  • Phương pháp xấp xỉ biến phân Galerkin: Sử dụng hệ hàm cơ sở trực giao không gian ba chiều ${h_1(z), h_2(z), h_3(z)}$ để chiếu phương trình vi phân đạo hàm riêng vào không gian trạng thái hữu hạn chiều.
  • Lý thuyết tập mờ và hệ suy luận mờ Mamdani: Sử dụng kỹ thuật mờ hóa để biểu diễn hàm phi tuyến thực nghiệm của hệ số dẫn nhiệt $\lambda(T)$ và nhiệt dung riêng $c(T)$ phụ thuộc vào nhiệt độ $T$, sử dụng hàm liên thuộc đầu vào dạng tam giác/hình thang, đầu ra singleton và phương pháp giải mờ trọng tâm (centroid).
  • Lý thuyết điều khiển tối ưu dạng toàn phương (LQR): Xây dựng thuật toán giải phương trình đại số Riccati (ARE) cho mô hình tuyến tính dừng cục bộ tại từng bước thời gian trượt, kết hợp khâu đạo hàm - quán tính bậc nhất để ước lượng bất định tham số và nhiễu.

Nguồn dữ liệu và công cụ phục vụ nghiên cứu bao gồm:

  • Các bộ thông số nhiệt vật lý của vật liệu thép (hệ số dẫn nhiệt $\lambda$, nhiệt dung riêng $c$, khối lượng riêng $\rho$) trích xuất từ tài liệu chuyên ngành luyện kim [14, 16, 37, 59].
  • Công cụ phần mềm tính toán và mô phỏng Matlab/Simulink và Toolbox Simscape (mô hình hóa vật lý hệ lò nung và vật nung truyền nhiệt đa lớp).

Nội dung chính theo từng chương

Chương 1: Bài toán điều khiển quá trình nung

Chương 1 phân tích các đặc trưng công nghệ và yêu cầu điều khiển trong quá trình nung kim loại công nghiệp, tập trung vào ba chỉ tiêu kỹ thuật:

  1. Nhiệt độ bề mặt vật nung: Phải được nâng theo giản đồ công nghệ thích hợp; nếu vượt ngưỡng cho phép sẽ làm tăng tốc độ tạo xỉ, nóng chảy lớp xỉ làm dính phôi vào đáy lò, gây mất mát kim loại do oxy hóa và khử cacbon.
  2. Độ đồng đều nhiệt độ theo tiết diện: Đối với vật nung dày ($Bi \ge 0{,}5$), gradient nhiệt độ giữa bề mặt và tâm xuất hiện rõ rệt. Trong giai đoạn nhiệt độ dưới $500^\circ\text{C}$ (trước khi kim loại đạt trạng thái biến dạng dẻo), tốc độ gia nhiệt phải được kiểm soát chặt chẽ để khống chế hiệu nhiệt độ giữa bề mặt và tâm, ngăn ngừa ứng suất nhiệt gây nứt vỡ phôi.
  3. Thời gian nung: Cần được tối thiểu hóa để giảm tiêu hao năng lượng nhiệt và hạn chế quá trình oxy hóa bề mặt, đảm bảo nhịp độ cấp phôi cho dây chuyền cán tiếp theo.

Tác giả chỉ ra thực trạng các hệ thống điều khiển lò nung truyền thống chủ yếu là điều khiển hở hoặc điều khiển vòng kín nhiệt độ khí lò/bề mặt vật nung. Từ đó, chương 1 đề xuất cấu trúc điều khiển dạng mô hình nội (Internal Model Control - IMC): kết hợp vòng điều khiển nhiệt độ lò nung bên ngoài với vòng điều khiển trường nhiệt độ bên trong lòng vật nung, lấy tín hiệu đặt là nhiệt độ tâm $T(0,t)$ và tín hiệu điều khiển $u(t)$ là cường độ dòng nhiệt cấp vào bề mặt.

Chương 2: Xây dựng mô hình quá trình truyền nhiệt trong vật nung

Chương 2 trình bày tổng quan các phương pháp mô hình hóa truyền nhiệt và phát triển mô hình trạng thái đề xuất:

1. Tổng quan các phương pháp mô hình hóa

  • Mô hình lý thuyết PDE: Thiết lập phương trình truyền nhiệt 1D, 2D và 3D dựa trên định luật cân bằng năng lượng. Phương trình truyền nhiệt 1D theo phương chiều dày $z \in [-Z/2, Z/2]$ có dạng: $$\frac{\partial T(z,t)}{\partial t} = \frac{1}{\rho c(z)} \frac{\partial}{\partial z} \left( \lambda(z) \frac{\partial T(z,t)}{\partial z} \right)$$
  • Các phương pháp số giải gần đúng PDE:
    • Phương pháp sai phân hữu hạn (FDM): Rời rạc hóa không gian và thời gian với tỷ lệ bước $M = \frac{a \Delta t}{(\Delta y)^2}$. Với bài toán 2D chọn $M = 0{,}25$ và bài toán 1D chọn $M = 0{,}5$, nhiệt độ nút tại thời điểm mới bằng trung bình cộng các nút lân cận.
    • Phương pháp phần tử hữu hạn (FEM-Galerkin): Sử dụng hàm hình dạng toàn cục $\phi_i(z)$, đưa phương trình truyền nhiệt về hệ phương trình vi phân đại số dạng ma trận: $C \frac{dT}{dt} + K T = f$.
    • Phương pháp thể tích hữu hạn (FVM): Rời rạc hóa tấm vật nung thành $N$ thể tích hữu hạn và tích phân phương trình bảo toàn trên từng phân miền không gian.
  • Mô hình thực nghiệm chia lớp đẳng nhiệt: Phân tích mô hình của Pedersen và Wittenmark [53] chia độ dày vật nung thành $n$ lớp, mỗi lớp mô tả bằng phương trình vi phân bậc nhất phi tuyến: $$\frac{dT_i}{dt} = g(T_i, h_i) \left[ a_{i,1}(T_{i-1}^4 - T_i^4) + a_{i,2}(T_{i-1} - T_i) + a_{i,3} T_i + a_{i,4} \right]$$ Nhược điểm của phương pháp này là số biến trạng thái tăng tuyến tính theo số lớp chia, gây khó khăn cho việc thiết kế điều khiển phản hồi liên tục.

2. Xây dựng mô hình Galerkin-Mờ đề xuất

Để khắc phục các hạn chế trên, tác giả áp dụng công thức xấp xỉ tách biến Galerkin, biểu diễn nghiệm nhiệt độ dưới dạng tổ hợp tuyến tính của ba hàm cơ sở không gian: $$T(z,t) = \sum_{i=1}^3 h_i(z) x_i(t) = h_1(z)x_1(t) + h_2(z)x_2(t) + h_3(z)x_3(t)$$ Trong đó hệ trục tọa độ không gian được định nghĩa: $$h_1(z) = 1, \quad h_2(z) = \frac{2z}{Z}, \quad h_3(z) = \left(\frac{2z}{Z}\right)^2 - \frac{1}{3}$$ Ý nghĩa vật lý của 3 biến trạng thái ảo $x(t) = [x_1(t), x_2(t), x_3(t)]^T$:

  • $x_1(t)$: Nhiệt độ trung bình của vật nung.
  • $x_2(t)$: Độ không đối xứng của biểu đồ phân bố nhiệt độ giữa hai bề mặt.
  • $x_3(t)$: Độ không đồng đều nhiệt độ tức thời bên trong vật nung.

Thay biểu thức khai triển vào phương trình truyền nhiệt, thu được hệ phương trình không gian trạng thái song tuyến bậc 3: $$\dot{x}(t) = \begin{bmatrix} 0 & 0 & 0 \ 0 & 0 & 0 \ 0 & 0 & -\frac{12\lambda(T)}{c(T)\rho Z^2} \end{bmatrix} x(t) + \begin{bmatrix} \frac{1}{c(T)\rho Z} & \frac{1}{c(T)\rho Z} \ -\frac{3}{c(T)\rho Z} & \frac{3}{c(T)\rho Z} \ \frac{5}{c(T)\rho Z} & \frac{5}{c(T)\rho Z} \end{bmatrix} u(t)$$ $$T(z,t) = \begin{bmatrix} h_1(z) & h_2(z) & h_3(z) \end{bmatrix} x(t)$$ Với $u(t) = [u_1(t), u_2(t)]^T$ là cường độ dòng nhiệt cấp vào mặt dưới và mặt trên ($W/m^2$).

Để xác định các tham số phụ thuộc nhiệt độ $\lambda(T)$ và $c(T)$, tác giả tích hợp khâu quan hệ mờ Mamdani:

  • Tín hiệu đầu vào nhiệt độ $T$ được mờ hóa thành 9 tập mờ $Tl_j$ cho $\lambda(T)$ và 12 tập mờ $Tc_i$ cho $c(T)$.
  • Tập mờ đầu ra là các singleton $ld_j$ và $c_i$.
  • Giải mờ theo phương pháp trọng tâm (centroid).
Thành phần thông số mờ Số lượng tập mờ đầu vào Dạng tập mờ đầu ra Phương pháp giải mờ
Hệ số dẫn nhiệt $\lambda(T)$ 9 tập mờ ($Tl_1 \dots Tl_9$) Singleton ($ld_j$) Điểm trọng tâm (Centroid)
Nhiệt dung riêng $c(T)$ 12 tập mờ ($Tc_1 \dots Tc_{12}$) Singleton ($c_i$) Điểm trọng tâm (Centroid)

Kết quả mô phỏng so sánh mô hình Galerkin-Mờ với mô hình lưới sai phân và mô hình vật lý Simscape cho thấy quỹ đạo nhiệt độ tâm $T(0,t)$ và phân bố trường nhiệt độ theo chiều dày $z$ có độ trùng khớp cao.

Chương 3: Điều khiển quá trình nung với bộ điều khiển LQR có ràng buộc trượt dọc trên trục thời gian

Chương 3 xây dựng và kiểm chứng thuật toán điều khiển tối ưu cho đối tượng truyền nhiệt:

1. Nguyên lý thiết kế bộ điều khiển

  • Cấu trúc điều khiển mô hình nội (IMC) được sử dụng để quan sát vector trạng thái ảo không đo được $x(t)$ từ mô hình Galerkin-Mờ.
  • Tại mỗi thời điểm lấy mẫu $t_k$, mô hình trạng thái song tuyến phi tuyến được tuyến tính hóa xung quanh điểm làm việc hiện tại thành mô hình tuyến tính dừng cục bộ $H_k(A_k, B_k, C_k)$ trên khoảng thời gian $[t_k, t_{k+1})$.
  • Tiêu chuẩn tối ưu tích phân bình phương sai lệch trạng thái và tín hiệu điều khiển được thiết lập: $$J_k = \int_{t_k}^{t_{k+1}} \left( \Delta x_k^T(t) Q_k \Delta x_k(t) + \Delta u_k^T(t) R_k \Delta u_k(t) \right) dt$$ Với $Q_k, R_k$ là các ma trận xác định dương được hiệu chỉnh bởi các hệ số $\alpha, \beta$.
  • Ma trận phản hồi trạng thái tối ưu $K_k = R_k^{-1} B_k^T L_k$ được tính toán thông qua việc giải nghiệm đối xứng xác định dương $L_k$ của phương trình đại số Riccati (ARE): $$A_k^T L_k + L_k A_k - L_k B_k R_k^{-1} B_k^T L_k + Q_k = 0$$

2. Xử lý bất định, nhiễu và ràng buộc

  • Thành phần bất định tham số $\Delta A, \Delta B$ và nhiễu tác động $d(t), n(t)$ được ước lượng thông qua khâu causal đạo hàm - quán tính bậc nhất với hằng số thời gian $\tau$.
  • Bộ điều khiển bổ sung thuật toán kiểm tra và xử lý điều kiện ràng buộc công nghệ về biên độ dòng nhiệt $u(t) \in [u_{\min}, u_{\max}]$.
  • Luận án phân tích và chứng minh tính ổn định tiệm cận của hệ thống dữ liệu lấy mẫu (sample-data system) theo lý thuyết.

3. Kết quả mô phỏng kiểm chứng

Hệ thống điều khiển được cài đặt và thử nghiệm mô phỏng trên đối tượng lò nung - vật nung xây dựng bằng Toolbox Simscape trong các kịch bản:

  • Bài toán bám quỹ đạo nhiệt độ tâm $T(0,t)$ và bám nhiệt độ bề mặt vật nung.
  • Thử nghiệm trong điều kiện có đồng thời nhiễu đầu vào $d(t)$, nhiễu đo $n(t)$ và sai lệch mô hình - đối tượng.
  • Kết quả mô phỏng xác nhận hệ thống điều khiển bám sát quỹ đạo đặt trước, triệt tiêu sai lệch tĩnh, giữ được chênh lệch nhiệt độ bề mặt và tâm trong giới hạn cho phép mà không vi phạm ràng buộc công suất cấp nhiệt.

Kết quả và những đóng góp mới

Luận án tự nêu các đóng góp mới về mặt khoa học và thực tiễn như sau:

Đóng góp về mặt lý luận

  1. Đề xuất phương pháp mô hình hóa quá trình truyền nhiệt trong vật nung dày bằng cách kết hợp công thức tách biến Galerkin bậc thấp với hệ suy luận mờ Mamdani (mô hình Galerkin-Mờ). Phương pháp này chuyển đổi thành công phương trình vi phân đạo hàm riêng (PDE) phi tuyến tham số phân bố rải thành hệ phương trình trạng thái song tuyến bậc 3 cố định, giữ nguyên được bản chất phụ thuộc nhiệt độ của các tham số $\lambda(T)$ và $c(T)$ mà không làm tăng bậc mô hình như phương pháp chia lớp đẳng nhiệt.
  2. Thiết lập cơ sở phân tích tính điều khiển được, tính quan sát được và tính đặt chỉnh (well-posed) của mô hình truyền nhiệt trong không gian Sobolev, làm tiền đề vững chắc cho việc thiết kế bộ quan sát và điều khiển trạng thái.
  3. Xây dựng thuật toán điều khiển LQR có tham số biến đổi trượt dọc trục thời gian cho hệ song tuyến, tích hợp bộ ước lượng sai lệch/nhiễu qua khâu lọc quán tính bậc nhất và thuật toán xử lý ràng buộc tín hiệu điều khiển, đồng thời chứng minh tính ổn định tiệm cận của hệ sample-data.

Đóng góp về mặt thực tiễn

  1. Đề xuất giải pháp điều khiển trực tiếp trường nhiệt độ bên trong lòng (tâm) vật nung bám theo quỹ đạo công nghệ đặt trước, thay thế cách tiếp cận truyền thống chỉ điều khiển nhiệt độ bề mặt hoặc nhiệt độ lò. Giải pháp này giúp kiểm soát chính xác độ chênh lệch nhiệt độ theo tiết diện, hạn chế tối đa ứng suất nhiệt và nguy cơ nứt phôi trong giai đoạn nung dưới $500^\circ\text{C}$.
  2. Xây dựng mô hình vật lý và cài đặt hoàn chỉnh thuật toán điều khiển trên công cụ mô phỏng Matlab/Simscape, cung cấp công cụ kiểm chứng mô phỏng tin cậy phục vụ nghiên cứu và chuyển giao công nghệ điều khiển tự động hóa quá trình nung công nghiệp.

Hạn chế và hướng nghiên cứu tiếp

Dựa trên các nội dung được trình bày trong văn bản luận án, các điểm hạn chế và hướng phát triển gồm:

  • Hạn chế:
    • Mô hình nghiên cứu mới chỉ xét sự truyền nhiệt 1D theo phương chiều dày $z$, dựa trên giả thiết nhiệt lượng cấp đều theo hai trục $x, y$.
    • Tín hiệu điều khiển được giả thiết trực tiếp là dòng nhiệt cấp vào bề mặt vật nung, chưa ghép nối đầy đủ với động học của hệ thống đốt nhiên liệu, cơ cấu chấp hành van mở khí và quá trình khí động học phức tạp trong không gian lò.
    • Các kết quả đánh giá mới dừng lại ở mức độ mô phỏng số và kiểm chứng trên mô hình vật lý Simscape, chưa tiến hành thử nghiệm trực tiếp trên lò nung công nghiệp thực tế.
  • Hướng nghiên cứu tiếp:
    • Mở rộng phương pháp mô hình hóa Galerkin-Mờ cho không gian hai chiều (2D) và ba chiều (3D).
    • Tích hợp vòng điều khiển nhiệt độ lò nung (điều khiển cơ cấu chấp hành van cấp nhiên liệu/không khí) với vòng điều khiển trường nhiệt độ vật nung thành hệ thống điều khiển tầng hoàn chỉnh.
    • Ứng dụng và thử nghiệm thực nghiệm thuật toán trên hệ thống thiết bị lò nung công nghiệp tại các nhà máy luyện kim.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu phục vụ các đối tượng sau:

  • Nghiên cứu sinh và học viên cao học ngành Kỹ thuật Điều khiển và Tự động hóa: Tham khảo phương pháp tiếp cận mô hình hóa hệ tham số phân bố rải (PDE) về hệ trạng thái phi tuyến/song tuyến (ODE) thông qua phương pháp biến phân Galerkin kết hợp logic mờ; kỹ thuật thiết kế bộ điều khiển LQR không dừng trượt dọc trục thời gian có ràng buộc và ước lượng nhiễu.
  • Giảng viên và nhà nghiên cứu kỹ thuật nhiệt, luyện kim: Tham khảo các nội dung tổng quan về phương pháp tính toán số truyền nhiệt (FDM, FEM, FVM) và mô hình chia lớp đẳng nhiệt trong vật nung dày.
  • Kỹ sư thiết kế hệ thống tự động hóa công nghiệp: Tham khảo quy trình xây dựng mô hình vật lý hệ "Lò nung - Vật nung" trên công cụ Matlab/Simscape phục vụ kiểm nghiệm thuật toán trước khi triển khai thực tế.

Câu hỏi thường gặp

1. Tại sao luận án lại đặt mục tiêu điều khiển nhiệt độ tại tâm vật nung thay vì nhiệt độ bề mặt như các nghiên cứu trước đây?
Đối với vật nung dày ($Bi \ge 0{,}5$), gradient nhiệt độ giữa bề mặt và tâm là rất lớn. Trong giai đoạn nhiệt độ dưới $500^\circ\text{C}$, nếu gia nhiệt quá nhanh, chênh lệch nhiệt độ lớn sẽ sinh ra ứng suất nhiệt gây nứt vỡ kim loại; ngược lại nếu nâng nhiệt quá chậm sẽ gây tiêu tốn năng lượng và tăng oxy hóa tạo xỉ. Điều khiển trực tiếp nhiệt độ tâm bám theo quỹ đạo đặt trước cho phép kiểm soát chính xác hiệu nhiệt độ giữa bề mặt và tâm trong giới hạn an toàn công nghệ.

2. Mô hình Galerkin-Mờ giải quyết được nhược điểm gì của các phương pháp mô hình hóa truyền thống?
Mô hình thực nghiệm chia lớp đẳng nhiệt làm tăng số phương trình vi phân và bậc trạng thái khi cần tăng độ mịn phân bố nhiệt. Trong khi đó, mô hình giải tích PDE lại rất phức tạp và khó áp dụng thuật toán điều khiển trạng thái khi tham số $\lambda(T), c(T)$ phụ thuộc nhiệt độ. Mô hình Galerkin-Mờ cố định bậc trạng thái ở bậc 3 với 3 biến trạng thái ảo, đồng thời sử dụng khâu mờ Mamdani để cập nhật phi tuyến các tham số nhiệt vật lý theo nhiệt độ mà không làm tăng bậc mô hình.

3. Ba biến trạng thái ảo $x_1(t), x_2(t), x_3(t)$ trong mô hình Galerkin-Mờ biểu thị đại lượng vật lý nào?

  • $x_1(t)$ biểu thị nhiệt độ trung bình của toàn bộ thể tích vật nung.
  • $x_2(t)$ biểu thị độ không đối xứng của biểu đồ nhiệt độ giữa mặt trên và mặt dưới của vật nung.
  • $x_3(t)$ biểu thị độ không đồng đều của nhiệt độ tức thời trong lòng vật nung. Trong giai đoạn đồng nhiệt cuối quá trình nung, $x_2(t)$ và $x_3(t)$ sẽ tiệm cận dần về 0.

4. Thuật toán LQR trượt dọc trên trục thời gian xử lý tính phi tuyến của đối tượng như thế nào?
Tại mỗi bước thời gian lấy mẫu $t_k$, thuật toán tuyến tính hóa mô hình trạng thái song tuyến quanh điểm làm việc hiện tại thành một mô hình tuyến tính dừng cục bộ $H_k(A_k, B_k, C_k)$. Bộ điều khiển giải phương trình đại số Riccati tại mỗi bước để cập nhật ma trận phản hồi $K_k$, kết hợp khâu lọc đạo hàm - quán tính bậc nhất để bù sai lệch mô hình và triệt tiêu ảnh hưởng của nhiễu.

5. Luận án đã sử dụng công cụ nào để kiểm chứng độc lập chất lượng của mô hình và bộ điều khiển đề xuất?
Tác giả sử dụng Toolbox Simscape của Matlab để xây dựng mô hình vật lý đa lớp độc lập của lò nung và vật nung. Mô hình Simscape đóng vai trò như đối tượng thực tế để kiểm chứng khả năng bám quỹ đạo nhiệt độ tâm và nhiệt độ bề mặt của bộ điều khiển khi chịu tác động của nhiễu đầu vào và nhiễu đo.

Kết luận

Luận án của tác giả Nguyễn Việt Dũng đã giải quyết bài toán mô hình hóa và điều khiển trường nhiệt độ cho vật nung dày trong lò luyện kim thông qua việc đề xuất mô hình trạng thái song tuyến Galerkin-Mờ bậc 3 kết hợp khâu suy luận mờ Mamdani. Trên nền tảng cấu trúc điều khiển mô hình nội, tác giả đã phát triển thuật toán điều khiển tối ưu LQR có tham số biến đổi trượt dọc trục thời gian, tích hợp khả năng ước lượng nhiễu và xử lý điều kiện ràng buộc công nghệ. Các kết quả mô phỏng trên nền tảng Matlab/Simscape đã xác nhận tính đúng đắn của phương pháp lý thuyết và khả năng điều khiển chính xác trường nhiệt độ trong lòng vật nung, mở ra hướng nghiên cứu ứng dụng nâng cao hiệu quả năng lượng trong ngành công nghiệp luyện kim.