Tổng quan luận án
Luận án tiến sĩ kỹ thuật với đề tài "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 nghiên cứu sinh 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, chuyên ngành Kỹ thuật điều khiển và Tự động hóa tại Trường Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên (bảo vệ năm 2018). Luận án giải quyết bài toán điều khiển tối ưu cho đối tượng nhiệt động lực học trong công nghiệp luyện kim và xử lý nhiệt vật liệu.
Tính cấp thiết và khoảng trống nghiên cứu
Lý thuyết điều khiển tối ưu cho hệ có tham số tập trung (Lumped Parameter Systems – LPS) đã được phát triển tương đối hoàn chỉnh, nhưng các nghiên cứu dành cho hệ với tham số phân bố (Distributed Parameter Systems – DPS) còn gặp nhiều rào cản khi triển khai vào thực tế. Hệ DPS được ứng dụng rộng rãi trong các quá trình gia nhiệt công nghiệp như tôi, ram, nhiệt luyện chi tiết cơ khí, ủ vật liệu từ, nung gạch men, cán thép hay sản xuất nhôm kính.
Trong các công nghệ trước đây sử dụng lò đốt nhiên liệu dầu nặng FO, hàm truyền của lò được mô tả bằng khâu quán tính có trễ với các hệ số truyền nhiệt được coi là hằng số tuyến tính và thời gian trễ tương đối nhỏ so với hằng số thời gian ($T/\tau \ge 10$). Tuy nhiên, đối với công nghệ nhiệt luyện chi tiết máy, tôi ram hoặc ủ vật liệu từ, lò nung chủ yếu là lò điện trở đốt nóng bằng dây điện trở. Khảo sát thực nghiệm nhận dạng lò điện trở trong dải nhiệt độ từ $0^\circ\text{C}$ đến $500^\circ\text{C}$ cho thấy hệ số truyền tĩnh $k$ thay đổi rất lớn theo nhiệt độ của lò ($k = k(v)$), mang bản chất phi tuyến rõ rệt. Đồng thời, khi thời gian trễ $\tau$ của lò lớn đáng kể so với hằng số thời gian $T$ ($6 \le T/\tau < 10$), việc sử dụng phép xấp xỉ chuỗi Taylor bậc nhất thông thường sẽ gây ra sai số lớn, làm giảm chất lượng điều khiển.
Khoảng trống nghiên cứu mà tác giả xác định là: chưa có công trình trong và ngoài nước đưa ra lời giải số hoàn chỉnh cho bài toán điều khiển tối ưu hệ tham số phân bố có trễ khi đồng thời xét đến tính phi tuyến của hệ số truyền tĩnh $k(v)$ và trường hợp trễ lớn ($6 \le T/\tau < 10$) trên đối tượng lò điện trở nung vật dày theo tiêu chuẩn nung chính xác nhất.
Mục tiêu nghiên cứu
Luận án xác định 4 mục tiêu nghiên cứu cụ thể được đánh số như sau:
- Xây dựng mô hình toán của đối tượng với tham số phân bố, có trễ, phi tuyến; xét cả trường hợp có hệ số trễ lớn (tỷ số $T/\tau$ nằm trong khoảng $6 \le T/\tau < 10$).
- Nghiên cứu thay thế khâu trễ $e^{-\tau s}$ trong khâu quán tính bậc nhất có trễ bằng phép xấp xỉ Pade bậc một khi đối tượng có thời gian trễ ($\tau$) là lớn đáng kể so với hằng số thời gian ($T$) của nó, tức là khi tỷ số $T/\tau$ thỏa mãn điều kiện $6 \le T/\tau < 10$.
- Tìm ra lời giải cho bài toán điều khiển tối ưu cho hệ với tham số phân bố, có trễ, phi tuyến bằng phương pháp số. Hệ này được đặc trưng bằng quá trình gia nhiệt một phía trong lò điện trở đối với vật dày. Trong đó quan tâm nhất tới tính phi tuyến (thay đổi) của hệ số truyền tĩnh $k$ của lò điện trở và trường hợp thời gian trễ $\tau$ là lớn đáng kể so với hằng số thời gian $T$ của lò.
- Mô phỏng và thực nghiệm để chứng minh tính chính xác và tính ổn định của nghiệm tối ưu.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Hệ thống điều khiển nhiệt độ lò điện trở và vật nung, cấu thành một hệ thống với tham số phân bố, có trễ và phi tuyến.
- Phạm vi nghiên cứu:
- Hình học và cơ chế truyền nhiệt: Vật nung có dạng tấm phẳng với bề dày $L$ nhỏ hơn nhiều so với chiều rộng và chiều dài; quá trình truyền nhiệt được khảo sát là truyền nhiệt một phía theo phương toạ độ không gian một chiều $x$.
- Đặc tính nhiệt động: Không gian lò đủ nhỏ so với kích thước vật nung; vật nung là vật dày có tiêu chuẩn Biot $\text{Bi} \ge 0{,}5$; vật nung được chia thành các lớp mỏng để coi hệ số dẫn nhiệt $\lambda$ và hệ số dẫn nhiệt độ $a$ là tuyến tính trên từng lớp.
- Đặc tính trễ và vùng công nghệ: Tỷ số giữa hằng số thời gian và thời gian trễ thỏa mãn $6 \le T/\tau < 10$; vùng nhiệt độ khảo sát không cao, giới hạn ở mức $\le 700^\circ\text{C}$.
- Công cụ và môi trường triển khai: Lập trình tính toán và mô phỏng trên phần mềm MATLAB/Simulink; thực nghiệm kiểm chứng trên mô hình vật lý tại phòng thí nghiệm 308 – Nhà Thí nghiệm, Khoa Điện tử, Trường Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên.
Tổng quan tài liệu và vị trí của luận án
Luận án đã hệ thống hóa lịch sử phát triển của lý thuyết điều khiển tối ưu cho hệ tham số phân bố (DPS) từ thập niên 1960 đến nay, phân tích các trường phái tiếp cận lý thuyết và thực nghiệm:
| Giai đoạn / Năm |
Tác giả tiêu biểu |
Hướng nghiên cứu và phương pháp tiếp cận |
| 1960 – 1969 |
Buttkovskii & Lerner (1960, 1969), Wang & Tung (1964), Sakawa (1964), Wang (1966), Axelband (1966), Sage & Chaudhuri (1967), Kim & Erzberger (1965), Alvarado & Mukundan (1969), Goldwyn, Sriram & Graham (1969) |
Đặt nền móng nguyên lý cực đại cho DPS; mô tả toán học bằng phương trình vi phân đạo hàm riêng (PDEs); phương pháp biến phân dẫn đến phương trình tích phân Fredholm loại 1; quy hoạch lồi; ma trận Riccati và phương trình Hamilton - Jacobi cho truyền sóng và truyền nhiệt tấm kim loại; ứng dụng biến đổi Laplace. |
| 1970 – 1989 |
Hassan & Solberg (1972), Julio (1972), Zone & Chang (1972), Davis & Perkins (1972), Yu & Seinfeld (1973), Curtain & Pritchard (1974, 1975), Balas (1978), Gibson (1979, 1981), Banks & Kunisch (1984), Pourki (1988), Taasan (1989), Shimemura (1989), Yoshida (1989), Matsumoto & Yoshida (1989) |
Điều khiển tối ưu toàn phương trong miền thời gian rời rạc; điều khiển tách kênh; vị trí cảm biến tối ưu; toán tử Riccati trong không gian Hilbert vô hạn chiều; phương pháp đa lưới (multi-grid); thiết kế bộ điều khiển toàn phương tuyến tính (LQR); xử lý hiện tượng tràn tín hiệu do xấp xỉ hữu hạn chiều. |
| 1990 – 2005 |
Bernstein & Rosen (1990), Li & Fadali (1991), Fattorini (1993), Burns & King (1994, 1998), Alli & Singh (1995, 2000), Joshi, Speyer & Kim (1995), Choe & Chang (1995, 1998), Farahi (1996), Banaszuk (1997), Drakunov & Barbieri (1997), Christofides & Daoutidis (1997), D’Andrea (1998), Godasi, Karakas & Palazoglu (1999, 2000), Padhi & Balakrishnan (2005) |
Lý thuyết chiếu tối ưu Bernstein/Hyland; phương pháp chiếu Galerkin và thiết kế bộ điều khiển bền vững LQG; so sánh phương pháp PDEs và phương pháp phương trình tích phân (IE) cho lò phản ứng hình ống; điều khiển cuốn chiếu (backstepping) trong không gian Sobolev; điều khiển trượt; phương pháp phân tích trực giao thích hợp (Proper Orthogonal Decomposition – POD / Karhunen - Loeve); thiết kế bộ điều khiển động học ngược tối ưu. |
| 2006 – 2014 |
Yadav (2006), Qian Li, Ning Li & Shaoyuan Li (2011), Natori (2012), Wu Qun Wei (2013), Behzad Talaei, Hao Xu & Jagannathan (2014), Pleshivtseva (2014) |
Điều khiển tối ưu thích nghi dựa trên mạng nơron; quy hoạch động xấp xỉ (ADP); xấp xỉ sai phân hữu hạn (FDM); ứng dụng biến đổi Fourier cho DPS bậc hai đa biến; áp dụng xấp xỉ Pade cho hệ có trễ với bộ quan sát nhiễu; tối ưu hóa tham số bài toán giá trị biên. |
| Nghiên cứu trong nước |
Nguyễn Hữu Công và cộng sự (tài liệu [7], [8], [9], [10], [11]) |
Giải bài toán điều khiển tối ưu DPS có trễ cho lò gia nhiệt đốt dầu FO bằng phương pháp số. Xử lý tác động điều khiển công suất $p(t)$, khâu trễ xấp xỉ bằng chuỗi Taylor với giả thiết trễ nhỏ ($T/\tau \ge 10$) và các tham số $a, \lambda, k$ là hằng số. |
Khoảng trống luận án lựa chọn giải quyết
Các công trình trước đây áp dụng nguyên lý cực đại Pontryagin hoặc phương pháp biến phân thường dẫn đến nghiệm điều khiển tối ưu nhiệt độ lò $v^*(t)$ có dạng đóng cắt nhảy cấp (bang-bang). Tín hiệu điều khiển này không thể thực thi trực tiếp trên lò nhiệt do quán tính nhiệt của môi trường không khí trong lò rất lớn. Nhằm khắc phục nhược điểm này, luận án không chọn $v(t)$ làm biến điều khiển trực tiếp mà chuyển sang điều khiển điện áp cung cấp cho lò $u(t)$ thông qua bộ biến đổi xoay chiều/xoay chiều dùng thyristor có quán tính điện rất nhỏ. Luận án tập trung giải quyết đồng thời hai thách thức kỹ thuật chưa được xử lý trong các công trình tiền nhiệm:
- Xử lý tính phi tuyến của hệ số truyền tĩnh $k(v)$ thay đổi theo nhiệt độ lò trong quá trình gia nhiệt.
- Ứng dụng xấp xỉ Pade bậc một để thay thế khâu trễ $e^{-\tau s}$ cho các đối tượng có độ trễ lớn trong dải $6 \le T/\tau < 10$, khắc phục sai số lớn của xấp xỉ Taylor.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết và mô hình toán học
Quá trình truyền nhiệt dẫn một chiều trong vật nung dạng tấm phẳng được mô tả bởi phương trình đạo hàm riêng parabolic bậc hai (phương trình dẫn nhiệt Fourier):
$$\frac{\partial^2 q(x,t)}{\partial x^2} = \frac{1}{a}\frac{\partial q(x,t)}{\partial t} \quad (0 \le x \le L, ; 0 \le t \le t_f)$$
Trong đó:
- $q(x,t)$ là nhiệt độ tại toạ độ $x$ ở thời điểm $t$ ($^\circ\text{C}$);
- $L$ là bề dày của vật nung ($\text{m}$);
- $t_f$ là thời gian nung cho phép ($\text{s}$);
- $a = \frac{\lambda}{C \rho}$ là hệ số dẫn nhiệt độ ($\text{m}^2/\text{s}$), với $\lambda$ là hệ số dẫn nhiệt ($\text{W}/\text{m}\cdot^\circ\text{C}$), $C$ là nhiệt dung riêng ($\text{J}/\text{kg}\cdot^\circ\text{C}$), và $\rho$ là khối lượng riêng ($\text{kg}/\text{m}^3$).
Các điều kiện biên và điều kiện ban đầu:
- Phân bố nhiệt độ ban đầu đồng nhất: $q(x,0) = q_0(x) = \text{const}$.
- Điều kiện biên loại 3 tại bề mặt tiếp xúc nhiệt ($x = 0$): $\lambda \frac{\partial q(x,t)}{\partial x}\Big|_{x=0} = -\alpha [v(t) - q(0,t)]$, với $\alpha$ là hệ số trao đổi nhiệt giữa không khí trong lò và bề mặt vật ($\text{W}/\text{m}^2\cdot^\circ\text{C}$), $v(t)$ là nhiệt độ môi trường không khí trong lò.
- Điều kiện biên cách nhiệt tại mặt đáy ($x = L$): $\frac{\partial q(x,t)}{\partial x}\Big|_{x=L} = 0$.
Mối quan hệ giữa điện áp cấp $u(t)$ và nhiệt độ môi trường không khí trong lò $v(t)$ là một khâu quán tính bậc nhất có trễ với hệ số truyền tĩnh $k$ phi tuyến phụ thuộc vào nhiệt độ $v$:
$$T \frac{dv(t)}{dt} + v(t) = k(v) \cdot u(t-\tau)$$
Để giải quyết tính phi tuyến, luận án tuyến tính hóa từng đoạn hệ số $k(v)$ thành $N$ đoạn với các giá trị hằng số $k_1, k_2, \dots, k_N$ tương ứng với các khoảng thời gian $t_1, t_2, \dots, t_N$ và dải nhiệt độ $v_1, v_2, \dots, v_N$.
Phiếm hàm mục tiêu và các điều kiện ràng buộc
Luận án giải bài toán nung chính xác nhất, hướng tới cực tiểu hóa sai lệch nhiệt độ giữa phân bố thực tế ở cuối quá trình gia nhiệt $q(x,t_f)$ và phân bố mong muốn $q^(x)$:
$$J_c = \int_0^L \left[ q^(x) - q(x,t_f) \right]^2 dx \to \min$$
Hệ thống chịu 5 điều kiện ràng buộc công nghệ trong toàn bộ thời gian $t \in [0, t_f]$:
- Giới hạn điện áp cung cấp cho lò: $U_1 \le u(t) \le U_2$ (hoặc $0 \le u(t) \le U_2$).
- Giới hạn nhiệt độ tối đa của môi trường lò (bảo vệ vật liệu xây lò): $v(t) \le U_3$.
- Giới hạn nhiệt độ bề mặt vật nung (chống quá nhiệt, nứt vỡ, nóng chảy): $q(0,t) \le U_4$.
- Giới hạn độ chênh lệch nhiệt độ giữa các lớp vật nung (giảm thiểu nội ứng suất nhiệt đối với vật dày $\text{Bi} \ge 0{,}5$ trong dải nhiệt độ nguy hiểm dưới $500^\circ\text{C}$): $\left| \frac{\partial q}{\partial x} \right| \le U_5$.
Phương pháp xấp xỉ Pade và tích phân số
Khâu trễ $e^{-\tau s}$ được xấp xỉ Pade bậc một (với $u=r=1$) cho tỷ số $6 \le T/\tau < 10$:
$$e^{-\tau s} \approx \frac{1 - \frac{\tau}{2}s}{1 + \frac{\tau}{2}s}$$
Phương pháp số chuyển đổi tích phân xác định trong hàm mục tiêu $J_c$ bằng công thức hình thang và công thức Simpson. Phép biến đổi Laplace thuận và ngược được sử dụng để xác định hàm truyền đạt toán tử $Q(x,s)$, từ đó tìm ra quan hệ giải tích giữa trường nhiệt độ $q(x,t)$ và các giá trị điện áp rời rạc $u_i$. Phiếm hàm $J_c$ liên tục được đưa về dạng hàm số của nhiều biến số thực, giải bằng các thuật toán tối ưu hóa quy hoạch toán học trên máy tính.
Nội dung chính theo từng chương
Chương 1: Tổng quan về điều khiển tối ưu cho hệ với tham số phân bố, có trễ, phi tuyến
Chương 1 phân loại bản chất toán học của hệ DPS (phương trình elliptic, hyperbolic và parabolic), làm rõ tính chất phức tạp của hệ DPS có trễ và phi tuyến khi các biến trạng thái phụ thuộc đồng thời vào thời gian $t$ và tọa độ không gian $x$.
Tác giả phân tích các phương pháp truyền thống (biến phân, nguyên lý cực đại Pontryagin) và chỉ ra hạn chế cốt lõi: việc giải trực tiếp hệ phương trình vi tích phân phi tuyến dẫn đến phương trình tích phân Fredholm loại 1 rất khó xác định nghiệm giải tích; tín hiệu điều khiển nhiệt độ lò $v^*(t)$ tối ưu nhận được có dạng xung vuông nhảy cấp (bang-bang), không tương thích với quán tính nhiệt lớn của lò. Từ đó, chương 1 xác lập hướng đi mới: chuyển đổi tín hiệu điều khiển sang điện áp $u(t)$ qua thyristor, tuyến tính hóa từng đoạn hệ số truyền tĩnh $k(v)$ và ứng dụng xấp xỉ Pade bậc một cho đối tượng có độ trễ lớn.
Chương 2: Đề xuất và giải bài toán điều khiển tối ưu cho hệ với tham số phân bố, có trễ, phi tuyến bằng phương pháp số sử dụng phép biến đổi Laplace
Chương 2 là nội dung trọng tâm về mặt lý luận của luận án, thực hiện các bước giải toán học và thực nghiệm nhận dạng đối tượng:
- Nhận dạng mô hình lò điện trở: Thực nghiệm cấp điện áp $u = 220\text{V}$ vào lò, thu thập dữ liệu đáp ứng nhiệt độ thời gian thực qua hệ thống cảm biến và xác định được hằng số thời gian $T$, thời gian trễ $\tau$, cùng bảng giá trị thực nghiệm của hệ số truyền tĩnh $k$ theo các dải nhiệt độ $v(t)$.
- Thiết lập công thức xấp xỉ Pade bậc một: Chứng minh toán học việc chuyển đổi hàm truyền đạt của lò $W(s) = \frac{k e^{-\tau s}}{Ts + 1}$ sang dạng phân thức hữu tỷ bậc hai khi thay $e^{-\tau s} \approx \frac{1 - \frac{\tau}{2}s}{1 + \frac{\tau}{2}s}$.
- Xác lập lời giải trường nhiệt độ: Sử dụng phép biến đổi Laplace giải phương trình đạo hàm riêng truyền nhiệt kết hợp điều kiện biên loại 3, thu được biểu thức phân bố nhiệt độ $q(x,t)$ và nhiệt độ bề mặt $q(0,t)$ theo điện áp $u(t)$.
- Rời rạc hóa và thuật toán tối ưu: Chia không gian thành $n$ lớp và thời gian $t_f$ thành $m$ khoảng rời rạc; áp dụng công thức hình thang và Simpson để tính xấp xỉ tích phân mục tiêu $J_c$; xây dựng thuật toán kiểm tra và ép các điều kiện ràng buộc kỹ thuật ($U_1, U_2, U_3, U_4, U_5$).
Chương 3: Các chương trình tính và các kết quả mô phỏng
Chương 3 trình bày cấu trúc các chương trình tính toán số và kết quả mô phỏng trên phần mềm MATLAB/Simulink:
- Xây dựng chương trình tính nghiệm riêng $\mu_i$, chương trình tính các hàm cơ sở $g_\mu(x,t)$, $g_\mu(x, t-\tau)$, và chương trình tối ưu hóa giải phiếm hàm $J_c$.
- Tiến hành khảo sát mô phỏng trên hai loại mẫu vật liệu nung điển hình có thông số nhiệt vật lý khác nhau: mẫu Samot (gạch chịu lửa samot) và mẫu Diatomite (vật liệu diatomite).
- Đánh giá sự phụ thuộc của sai số hàm mục tiêu $e = J_c$ vào số bước chia không gian $n$, số bước chia thời gian $m$, và tổng thời gian nung $t_f$:
| Mẫu vật liệu |
Tham số khảo sát ($n, m, t_f$) |
Kết quả sai số $e$ (Xấp xỉ Taylor) |
Kết quả sai số $e$ (Xấp xỉ Pade 1) |
Nhận xét từ dữ liệu mô phỏng |
| Mẫu Samot |
$n=6, m=40, t_f = 4200\text{s}$ |
$e = 5{,}4$ |
$e = 3{,}4$ |
Pade 1 giảm đáng kể sai số nhiệt độ so với Taylor |
| Mẫu Samot |
$n=8, m=64, t_f = 4200\text{s}$ |
$e = 3{,}6$ |
$e = 2{,}5$ |
Khi tăng số bước chia $n, m$, sai số hàm mục tiêu giảm dần |
| Mẫu Samot |
$n=10, m=100, t_f = 4200\text{s}$ |
$e = 2{,}8$ |
$e = 2{,}1$ |
Lưới tính $n=10, m=100$ cho độ chính xác cao nhất |
| Mẫu Samot |
$t_f = 3900\text{s}$ |
$e = 6{,}9$ |
$e = 5{,}10$ |
Giảm thời gian nung $t_f$ làm tăng sai số do quán tính nhiệt |
| Mẫu Samot |
$t_f = 3600\text{s}$ |
$e = 4{,}12$ |
$e = 1{,}11$ |
Pade 1 duy trì độ chính xác tốt hơn khi thời gian nung rút ngắn |
| Mẫu Samot |
$t_f = 3300\text{s}$ |
$e = 0{,}91$ |
$e = 0{,}13$ |
Phù hợp với chế độ nung tối ưu đặt ra |
| Mẫu Diatomite |
$t_f = 4200\text{s}$ |
$e = 2{,}16$ |
$e = 1{,}17$ |
Kiểm chứng thành công trên vật liệu diatomite có tính dẫn nhiệt kém |
| Mẫu Diatomite |
$t_f = 3900\text{s}$ |
$e = 0{,}19$ |
$e = 0{,}18$ |
Chất lượng phân bố nhiệt đạt độ đồng đều cao |
| Mẫu Diatomite |
$t_f = 3600\text{s}$ |
$e = 1{,}21$ |
$e = 1{,}20$ |
Nghiệm tối ưu ổn định ở các điều kiện nung khác nhau |
Kết quả mô phỏng khẳng định việc áp dụng phép xấp xỉ Pade bậc một luôn cho sai số hàm mục tiêu nhỏ hơn và đáp ứng nhiệt bám sát nhiệt độ đặt hơn so với phép xấp xỉ Taylor khi đối tượng có trễ lớn ($6 \le T/\tau < 10$).
Chương 4: Thực nghiệm kiểm chứng chất lượng phương pháp đã đề xuất trên mô hình hệ thống thực
Chương 4 thực hiện thí nghiệm vật lý trên hệ thống lò điện trở thực tế để kiểm tra tính đúng đắn của thuật toán:
- Cấu hình hệ thống thực nghiệm: Bao gồm lò nhiệt điện trở, mẫu nung Samot và Diatomite, các cặp cảm biến nhiệt độ (can nhiệt), mạch khuếch đại tín hiệu can nhiệt, Card thu thập dữ liệu NI USB-6008 của hãng National Instruments, bộ biến đổi điện áp xoay chiều/xoay chiều 1 pha sử dụng 2 thyristor đấu song song ngược và mạch điều khiển pha xung.
- Kịch bản thí nghiệm thực:
- Mẫu Samot: Thực hiện với nhiệt độ đặt $q^* = 300^\circ\text{C}$ ($t_f = 4200\text{s}$) và $q^* = 400^\circ\text{C}$ ($t_f = 4200\text{s}$).
- Mẫu Diatomite: Thực hiện với nhiệt độ đặt $q^* = 300^\circ\text{C}$ ($t_f = 4500\text{s}$) và $q^* = 400^\circ\text{C}$ ($t_f = 4500\text{s}$).
- Kết quả đo đạc thực nghiệm cho thấy quy luật biến thiên của điện áp tối ưu $u^*(t)$, nhiệt độ lò $v(t)$ và phân bố nhiệt độ tại các lớp trong mẫu nung hoàn toàn trùng khớp với quy luật thu được từ tính toán lý thuyết và mô phỏng trong Chương 2 và Chương 3, đảm bảo các điều kiện giới hạn công nghệ $U_1 \dots U_5$.
Kết quả và những đóng góp mới
Đóng góp mới về mặt lý luận khoa học
- Thiết lập thành công mô hình toán học và tìm ra lời giải tường minh xác định trường nhiệt độ $q(x,t)$ trong vật nung tấm phẳng truyền nhiệt một phía khi biết tín hiệu điều khiển điện áp $u(t)$ cấp cho lò điện trở, xét trong điều kiện hệ số truyền tĩnh $k$ của lò là phi tuyến ($k=k(v)$) và đối tượng có thời gian trễ lớn ($6 \le T/\tau < 10$).
- Đề xuất phương pháp giải bài toán điều khiển tối ưu cho hệ tham số phân bố, có trễ, phi tuyến theo tiêu chuẩn nung chính xác nhất bằng phương pháp số sử dụng phép biến đổi Laplace kết hợp với phép xấp xỉ Pade bậc một.
Đóng góp mới về mặt thực tiễn
- Luận án là một trong những công trình khoa học đầu tiên tại Việt Nam giải quyết trọn vẹn bài toán điều khiển tối ưu hệ tham số phân bố, có trễ, phi tuyến trên đối tượng lò điện trở nung vật dày dạng tấm phẳng.
- Xây dựng bộ chương trình tính toán số và mô phỏng trên môi trường MATLAB/Simulink, cho phép tính toán quỹ đạo điện áp điều khiển tối ưu $u^*(t)$ đáp ứng các ràng buộc công nghệ phức tạp.
- Chế tạo và tích hợp thành công hệ thống phần cứng thực nghiệm với card NI USB-6008 và mạch điều khiển pha xung thyristor; kiểm chứng thành công trên hai mẫu vật liệu thực tế (Samot và Diatomite), cung cấp cơ sở kỹ thuật tin cậy để chuyển giao ứng dụng vào các hệ thống lò nhiệt luyện công nghiệp.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của luận án
- Mô hình truyền nhiệt chỉ khảo sát một chiều theo phương bề dày $x$ trên vật nung dạng tấm phẳng, chưa xét truyền nhiệt đa chiều (2 chiều hoặc 3 chiều).
- Phạm vi nhiệt độ khảo sát bị giới hạn ở mức $\le 700^\circ\text{C}$, chưa nghiên cứu vùng nhiệt độ cao hơn (nơi bức xạ nhiệt chiếm ưu thế tuyệt đối và các hệ số $a, \lambda$ biến thiên mạnh).
- Các hệ số nhiệt vật lý của vật nung ($a, \lambda$) được giả thiết tuyến tính (hằng số) trên từng lớp mỏng trong suốt quá trình nung; luận án mới chỉ tập trung xử lý tính phi tuyến của hệ số truyền tĩnh $k$ và độ trễ của lò điện trở.
- Luận án mới áp dụng phép xấp xỉ Pade bậc một cho khoảng tỷ số $6 \le T/\tau < 10$, chưa giải chi tiết cho các bậc Pade cao hơn (bậc hai, bậc ba) khi tỷ số $T/\tau < 6$.
Hướng nghiên cứu tiếp theo
- Mở rộng bài toán điều khiển tối ưu cho các vật nung có cấu trúc không gian 2 chiều, 3 chiều và hình dạng phức tạp (hình trụ, khối hộp).
- Nghiên cứu bài toán tối ưu khi đồng thời tính đến sự biến thiên phi tuyến của cả hệ số dẫn nhiệt $\lambda(t)$ và hệ số dẫn nhiệt độ $a(t)$ của vật nung theo nhiệt độ.
- Nghiên cứu thuật toán điều khiển tối ưu khi sử dụng phép khai triển Pade bậc hai (với $2 \le T/\tau < 6$) và Pade bậc ba (với $1 \le T/\tau < 2$) cho các hệ thống có độ trễ cực lớn.
- Phát triển hệ thống điều khiển tối ưu vòng kín có phản hồi trạng thái nhiệt độ thực để nâng cao khả năng triệt tiêu nhiễu trong môi trường công nghiệp.
Giá trị tham khảo
Luận án mang lại giá trị tham khảo thiết thực cho nhiều nhóm đối tượng:
- Nghiên cứu sinh và học viên cao học ngành Kỹ thuật điều khiển, Tự động hóa và Toán ứng dụng: Có thể tham khảo phương pháp giải số cho phương trình đạo hàm riêng parabolic kết hợp phương trình vi phân có trễ phi tuyến, phương pháp chuyển đổi phiếm hàm mục tiêu về hàm nhiều biến, và kỹ thuật xấp xỉ Pade trong miền toán tử Laplace.
- Giảng viên và nhà nghiên cứu: Sử dụng tài liệu Chương 1 và Chương 2 làm học liệu tham khảo chuyên sâu cho các học phần Điều khiển hệ tham số phân bố, Điều khiển tối ưu nâng cao và Nhận dạng hệ thống.
- Kỹ sư thiết kế và chuyên gia công nghệ nhiệt: Chương 2, Chương 3 và Chương 4 cung cấp sơ đồ mạch lực thyristor, cấu hình thu thập dữ liệu với NI USB-6008, cùng thuật toán điều khiển điện áp nung chính xác áp dụng trực tiếp cho các dây chuyền tôi ram, ủ vật liệu từ, sấy gốm sứ và sản xuất vật liệu chịu lửa.
Câu hỏi thường gặp
1. Tại sao tác giả chọn điều khiển điện áp cung cấp $u(t)$ qua bộ biến đổi thyristor thay vì điều khiển trực tiếp nhiệt độ lò $v(t)$?
Trong các nghiên cứu giải bằng phương pháp biến phân hoặc nguyên lý cực đại Pontryagin, nghiệm tối ưu nhiệt độ lò $v^*(t)$ thường có dạng nhảy cấp (bang-bang). Tuy nhiên, môi trường không khí trong lò có quán tính nhiệt rất lớn nên không thể biến thiên tức thời theo dạng xung vuông. Khi chuyển sang điều khiển điện áp $u(t)$ cấp cho dây đốt thông qua bộ biến đổi thyristor, tín hiệu điều khiển có thể thay đổi nhanh chóng do quán tính của các phần tử điện tử công suất là không đáng kể, đảm bảo tính khả thi khi triển khai thực tế.
2. Sự khác biệt cốt lõi giữa luận án này và công trình của Nguyễn Hữu Công (tài liệu [10]) là gì?
Công trình [10] giải bài toán điều khiển tối ưu cho lò đốt dầu nặng FO với biến điều khiển là công suất cấp $p(t)$, giả thiết trễ nhỏ ($T/\tau \ge 10$) sử dụng xấp xỉ Taylor, và coi các hệ số $a, \lambda, k$ là hằng số tuyến tính. Luận án của Mai Trung Thái nghiên cứu lò điện trở điều khiển bằng điện áp $u(t)$, giải quyết bài toán với hệ số truyền tĩnh $k(v)$ phi tuyến thay đổi theo nhiệt độ lò và đối tượng có độ trễ lớn ($6 \le T/\tau < 10$) sử dụng phép xấp xỉ Pade bậc một.
3. Tiêu chuẩn chọn bậc xấp xỉ Pade cho khâu trễ $e^{-\tau s}$ được quy định như thế nào trong luận án?
Căn cứ theo tài liệu [7], việc lựa chọn phương pháp xấp xỉ khâu trễ $e^{-\tau s}$ phụ thuộc vào tỷ số giữa hằng số thời gian $T$ và thời gian trễ $\tau$:
- $T/\tau \ge 10$: Sử dụng khai triển chuỗi lũy thừa Taylor;
- $6 \le T/\tau < 10$: Sử dụng khai triển Pade bậc một (Pade 1, $r=1$);
- $2 \le T/\tau < 6$: Sử dụng khai triển Pade bậc hai (Pade 2, $r=2$);
- $1 \le T/\tau < 2$: Sử dụng khai triển Pade bậc ba (Pade 3, $r=3$).
Luận án tập trung giải quyết trường hợp $6 \le T/\tau < 10$ với xấp xỉ Pade bậc một.
4. Luận án đã thực hiện mô phỏng và thực nghiệm trên những mẫu vật liệu nào và dùng thiết bị phần cứng gì?
Luận án thực hiện mô phỏng và thực nghiệm trên 2 mẫu vật liệu nung điển hình là mẫu Samot (gạch chịu lửa) và mẫu Diatomite (vật liệu cách nhiệt). Hệ thống phần cứng thực nghiệm sử dụng card chuyển đổi NI USB-6008 kết nối máy tính chạy MATLAB/Simulink, mạch điều khiển pha xung, bộ biến đổi xoay chiều/xoay chiều 1 pha dùng 2 thyristor song song ngược, cảm biến can nhiệt và mạch khuếch đại tín hiệu can nhiệt tại phòng thí nghiệm 308, Khoa Điện tử – Trường Đại học Kỹ thuật Công nghiệp Thái Nguyên.
5. Tại sao cần có ràng buộc $|\partial q / \partial x| \le U_5$ trong quá trình gia nhiệt vật dày?
Đối với vật nung dày ($\text{Bi} \ge 0{,}5$), sự chênh lệch nhiệt độ giữa bề mặt và tâm vật là rất lớn. Trong giai đoạn nhiệt độ biến thiên từ nhiệt độ môi trường đến nhiệt độ biến dạng dẻo ($500^\circ\text{C}$), gradient nhiệt độ quá lớn sẽ sinh ra nội ứng suất nhiệt, gây nứt vỡ hoặc biến dạng vật liệu. Điều kiện ràng buộc $|\partial q / \partial x| \le U_5$ bắt buộc thuật toán điều khiển phải khống chế độ chênh lệch nhiệt độ giữa các lớp không vượt quá giới hạn an toàn cho phép.
Kết luận
Luận án tiến sĩ của tác giả Mai Trung Thái đã giải quyết thành công 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ệ tham số phân bố, có trễ và phi tuyến trên đối tượng lò điện trở gia nhiệt cho vật dày dạng tấm phẳng. Bằng cách kết hợp phép biến đổi Laplace, kỹ thuật tuyến tính hóa từng đoạn hệ số truyền tĩnh phi tuyến $k(v)$, phép xấp xỉ Pade bậc một cho khâu trễ lớn ($6 \le T/\tau < 10$) và phương pháp số, tác giả đã xây dựng thuật toán xác định chính xác quỹ đạo điện áp điều khiển tối ưu $u^*(t)$. Các kết quả mô phỏng số và thực nghiệm trên mô hình vật lý với mẫu Samot và Diatomite qua card NI USB-6008 đã chứng minh tính đúng đắn, độ tin cậy và khả năng ứng dụng thực tiễn của phương pháp đã đề xuất.
