Tóm tắt Luận án tiến sĩ Kỹ thuật: 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

BỘ GIÁO DỤC VÀ ĐÀO TẠO
ĐẠI HỌC THÁI NGUYÊN
-----------------***-----------------

Mai Trung Thá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

Chuyên ngành: Kỹ thuật điều khiển và tự động hóa
Mã số: 9 52 02 16

TĨM TẮT LUẬN ÁN TIẾN SĨ KỸ THUẬT

Thái Nguyên - 2018


Cơng trình được hồn thành tại:
Trường Đại học Kỹ thuật Công Nghiệp Thái Nguyên

Người hướng dẫn khoa học: PGS.TS Nguyễn Hữu Công

Phản biện 1: ………………………………
Phản biện 2: ………………………………
Phản biện 2: ………………………………

Luận án sẽ được bảo vệ trước Hội đồng đánh giá luận án tiến sĩ cấp Trường
họp tại Trường Đại học Kỹ thuật Cơng nghiệp Thái ngun
Vào hồi……., giờ…….ngày…….tháng…….năm…………

Có thể tìm hiểu luận án tại:

1. Thư viện Quốc gia Việt Nam
2. Trung tâm học liệu Đại học Thái Nguyên
3. Thư viện Trường Đại học Kỹ thuật Công nghiệp Thái Nguyên


DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC CĨ
LIÊN QUAN ĐẾN LUẬN ÁN

1. Mai Trung Thái, Nguyễn Thị Mai Hương (2013), "Điều khiển tối ưu cho một
hệ với tham số phân bố sử dụng phương pháp Gradient ", Tạp chí Khoa học
Công nghệ - Đại học Thái Nguyên, Số 10, tập 110, tr. 45 - 52.
2. Mai Trung Thai, Nguyen Huu Cong, Nguyen Van Chi, Vu Van Dam,
(2017), “Applying Pade approximation model in optimal control problem for a
distributed parameter system with time delay”, International Journal of
Computing and Optimization, HIKARI Ltd, Vol.4, no.1, 2017, pp. 19-30
3. Mai Trung Thái, Nguyễn Thị Mai Hương (2017), “Hai phương pháp thay
thế đối tượng có trễ trong bài toán điều khiển tối hệ với tham số phân bố”,
ISSN 1859-1531, Tạp chí Khoa học và Cơng nghệ Đại học Đà nẵng, số 5
(114). 2017 – quyển 1.
4. Cong Huu Nguyen, Mai Trung Thai (2018), “Optimal control for a distributed
parameter system with time-delay, non-linear using the numerical method.
Application to one-sided heat conduction system”, ISSN 2395-0250,
International Journal of Thermal Engineering (IJTE), Vol 4, Issue 1, Jan-Feb
2018


1

MỞ ĐẦU
1. Đặt vấn đề

Lý thuyết điều khiển tối ưu đã được nghiên cứu từ lâu song cho tới nay các tác giả chủ
yếu nghiên cứu bài toán điều khiển tối ưu cho hệ có tham số tập trung mà chưa quan tâm nhiều
tới bài toán điều khiển tối ưu cho hệ với tham số phân bố. Điều khiển tối ưu cho hệ với tham
số phân bố được ứng dụng trong nhiều lĩnh vực khác nhau như: tôi, ram, nhiệt luyện, ủ vật liệu
từ, nung gạch men, cán thép,….Trong một số cơng nghệ, q trình gia nhiệt được thực hiện
trong lị nung thường bằng dầu nặng FO, ví dụ như q trình nung trong cán thép hay nung
phơi khi sản xuất nhơm kính. Trong trường hợp này, hàm truyền của lị nung là khâu qn tính
có trễ, cịn mối quan hệ giữa nhiệt độ lị là các phương trình đạo hàm riêng dạng parabolic với
điều kiện biên loại 3. Nếu ta xét bài toán điều khiển tối ưu cho quá trình “nung chính xác
nhất”, lúc này đối tượng điều khiển là hệ với tham số phân bố, có trễ. Với bài toán này, đã
được một số tác giả quan tâm và tìm được lời giải bằng phương pháp biến phân, phương pháp
dùng nguyên lý cực đại của Pontryagin hay phương pháp số như trong [8,10,72]. Trong đó
phương pháp số tỏ ra ưu việt hơn cả. Tuy nhiên trong một số cơng nghệ khác, lị nung là lị
điện, tức là đốt bằng dây điện trở như q trình tơi, ram, nhiệt luyện các chi tiết cơ khí, ủ vật
liệu từ, v.v…Lúc này hàm truyền của lò điện trở cũng là khâu qn tính bậc nhất có trễ dạng:
W( s ) 

Y( s )
k .e  s

X( s ) (s  1)

(0.1)

Nhưng, lúc này k là hệ số phụ thuộc vào nhiệt độ trong lò. Thực tế qua việc nhận dạng
lò điện trở thì k thay đổi khá nhiều, ví dụ như trong lò điện trở với dải nhiệt độ thay đổi từ
0-5000C. (Việc này sẽ được chứng minh ở phần sau). Vậy nếu vẫn xét bài toán điều khiển tối
ưu cho q trình “nung chính xác nhất” thì đây là bài toán điều khiển tối ưu cho đối tượng với
tham số phân bố, có trễ, phi tuyến. Chính sự phi tuyến của k làm cho lời giải của bài toán trở
nên rất phức tạp. Do vậy để bài tốn có thể được ứng dụng trong thực tế, luận án này tìm cách

đưa ra lời giải cho bài toán với điểm khác biệt lớn nhất là tính phi tuyến của k . Bài toán điều
khiển tối ưu vẫn được thực hiện bằng phương pháp số. Lời giải cho trường hợp xét tới tính phi
tuyến của k chưa được các tác giả trong và ngoài nước nghiên cứu. Ngoài ra, để mở rộng bài
toán điều khiển tối ưu, luận án cũng xét thêm trường hợp hệ số trễ () của lò điện trở là lớn
đáng kể so với hằng số thời gian (T) của nó.
2. Tính cấp thiết của luận án
Điều khiển tối ưu theo tiêu chuẩn nung chính xác nhất cho hệ với tham số phân bố được
ứng dụng trong nhiều lĩnh vực khác nhau trong các lĩnh vực công nghiệp. Các nghiên cứu
trước đây [10,79] cũng đã giải quyết bài toán điều khiển tối ưu cho hệ với tham số phân bố, có
trễ. Nếu trong lĩnh vực lị nung thì bài tốn này đã được áp dụng cho các cơng nghệ lò đốt
bằng dầu nặng FO. Tuy nhiên, với một số công nghệ như ủ vật liệu từ, tôi ram nhiệt luyện chi
tiết máy thì lị nung được thực hiện bằng lị điện. Vì vậy đây là bài tốn điều khiển tối ưu cho
hệ với tham số phân bố, có trễ, phi tuyến. Với bài tốn này, hiện nay chưa có sự nghiên cứu
của các tác giả trong và ngoài nước, vì vậy đề tài này có tính cấp thiết và nếu được giải quyết
sẽ một mặt bổ sung vào lý thuyết điều khiển cho hệ có tham số phân bố, mặt khác cũng mở ra
khả năng ứng dụng vào thực tế.


2

3. Mục tiêu của luận án
Xây dựng mơ hình tố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ìm ra lời giải cho bài tố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 q 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ở. Ngồi ra cịn quan tâm tới trường hợp thời
gian trễ () 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.
4. Đối tượng, phạm vi và phương pháp nghiên cứu của luận án
- Đối tượng nghiên cứu: Hệ thống điều khiển nhiệt độ lò điện trở và vật nung, đó là một hệ

với tham số phân bố, có trễ, phi tuyến.
- Phạm vi nghiên cứu: Nghiên cứu đối tượng động học có trễ mà có thời gian trễ () là lớn
đáng kể so với hằng số thời gian (T) của nó, tức là khi đối tượng có tỷ số T/ thỏa mãn điều
kiện 6  T/ < 10 [7], sau đó áp dụng vào việc 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. Lập trình mơ phỏng trên máy tính và thí nghiệm trên mơ hình vật
lý cụ thể.
- Phương pháp nghiên cứu: Phân tích, đánh giá các cơng trình nghiên cứu đã được cơng bố
trên các bài báo, tạp chí, các luận án về điều khiển tối ưu cho hệ với tham số phân bố, có trễ,
phi tuyến nhằm xác định chắc chắn các mục tiêu đề ra. Nghiên cứu các phương pháp 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. Lập trình mơ phỏng trên
Matlab & Simulink để kiểm chứng lại lý thuyết. Tiến hành thí nghiệm kiểm chứng kết quả
nghiên cứu lý thuyết và kết quả mơ phỏng trên mơ hình vật lý hệ thống phi tuyến cụ thể (lò
điện trở và vật nung).
5. Những đóng góp mới về mặt khoa học và thực tiễn của luận án
Luận án đã có những đóng góp mới như sau:
- Ứng dụng phép biến đổi Laplace, phương pháp tuyến tính hóa từng đoạn và khai triển
Pade đã đưa ra được lời giải tường minh cho bài toán tìm được trường nhiệt độ của vật nung khi
biết điện áp cung cấp cho lò - xét cả trong trường hợp hệ số truyền tĩnh k của lò là phi tuyến và
thời gian trễ của lị lớn.
- Đã tìm ra lời giải cho bài toán điều khiển tối ưu hệ với tham số phân bố, có trễ, phi
tuyến bằng phương pháp số, có tính tới cả điều kiện giới hạn pha – Hệ này ứng dụng cho bài
tốn nung chính xác nhất trong q trình gia nhiệt.
- Đã mơ phỏng và tiến hành thực nghiệm để kiểm chứng kết quả nghiên cứu lý thuyết
làm cơ sở cho việc triển khai điều khiển thực tế nung chính xác ở lị điện trở khi gia nhiệt cho
vật nung dạng tấm phẳng.
Những đóng góp trên có ý nghĩa khoa học, có giá trị thực tiễn, có thể áp dụng cho bài
tốn nung chính xác nhất trong q trình gia nhiệt ở lị điện trở, ví dụ áp dụng trong một số lĩnh
vực như q trình tơi, ram, nhiệt luyện các chi tiết máy, ủ vật liệu từ, v.v…
6. Cấu trúc của luận án:
Luận án được trình bày trong 4 chương chính và phần kết luận như sau:

3

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 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 3: Các chương trình tính tốn và các kết quả mô phỏng
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.
- Kết luận và kiến nghị
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.
1.1. Tổng quan chung
1.2. Tổng quan các cơng trình nghiên cứu về điều khiển tối ưu cho hệ với tham số phân
bố, có trễ, phi tuyến trong và ngồi nước.
Lý thuyết về điều khiển tối ưu cho hệ với tham số phân bố (DPS) đã được nghiên cứu từ
thập niên 60 của thế kỷ trước. Buttkovskii và Lerner đã đưa ra bài báo đầu tiên trong lĩnh vực
này vào năm 1960 [36], bắt đầu từ nguyên lý cực đại cho một lớp các hệ thống tham số phân
bố. Điều này đã cho ra một loạt các bài báo từ Butkovskii [32,34,35]. Các nghiên cứu này đã
đề cập đến việc mơ tả bài tốn và ngun lý cực đại cho một hệ tham số phân bố được mô tả
bởi một tập các phương trình tích phân phi tuyến. Các kết quả nghiên cứu [29,46,72,86], các
tác giả đã dùng nguyên lí cực đại của Pôntriagin và phương pháp biến phân để đưa ra lời giải
cho bài toán điều khiển tối ưu một hệ với tham số phân bố (cụ thể là bài tốn truyền nhiệt một
phía trong lị gia nhiệt). Khi tìm nghiệm tối ưu thường dẫn đến phải giải các phương trình
Fredhom loại một nên rất khó giải, khó khẳng định được sự tồn tại nghiệm. Đặc biệt khi hàm
điều khiển v(t) có kèm theo điều kiện ràng buộc thì việc tìm nghiệm của các phương trình trên
là rất khó khăn. Việc ứng dụng các kết quả này vào thực tế gặp khó khăn vì nghiệm tối ưu v*(t)
là các hàm đổi dấu tức thời, tức là tín hiệu điều khiển thuộc dạng “bang-bang”- xung vng.

Song v*(t) ở đây chính là nhiệt độ lị nên khơng thể thay đổi đột ngột được vì có qn tính nhiệt
khá lớn, nhất là khi tần suất thay đổi nhiều.
Ở Việt Nam, vấn đề này đã được một số học giả tiếp cận, nghiên cứu trong khoảng hơn
một thập niên trở lại đây. Các kết quả nghiên cứu chủ yếu là các bài báo, luận văn thạc sỹ của
một số học viên cao học thuộc một số trường Đại học trong cả nước. Theo hiểu biết của tác
giả, hiện nay các cơng trình nghiên cứu về hệ với tham số phân bố, có trễ, phi tuyến, đặc biệt
khi hệ này được áp dụng cho bài tốn truyền nhiệt một phía trong lị nung là lò điện (tức là đốt
bằng dây điện trở áp dụng cho các q trình tơi, ram, nhiệt luyện các chi tiết cơ khí, ủ vật liệu
từ,…) để điều khiển nhiệt độ cho vật nung có dạng tấm phẳng theo tiêu chuẩn nung chính xác
nhất ở trong nước chưa có tác giả nào nghiên cứu, chủ yếu mới dừng lại ở nghiên cứu về hệ
với tham số phân bố, có trễ, điển hình là một số bài báo và luận án tiến sĩ như [7,8,9,10,11]
Đặc biệt, theo [10], luận án nghiên cứu giải bài toán điều khiển tối ưu cho hệ với tham
số phân bố có trễ, hệ này được áp dụng cho hệ thống truyền nhiệt một phía trong lị nung để
điều khiển nhiệt độ cho phôi tấm theo tiêu chuẩn nung chính xác nhất, q trình gia nhiệt cho


4

lị nung bằng dầu nặng FO, ví dụ như q trình nung trong cán thép hay nung phơi khi sản xuất
nhơm kính. [10] đã đưa ra hướng khắc phục nhược điểm của [29,46,72,86] bằng cách không
dùng tác động điều khiển là nhiệt độ lị v(t) mà là cơng suất p(t) cung cấp cho lị thơng qua một
khâu chuyển đổi. Cơng suất cung cấp cho lò nếu dùng năng lượng điện thì việc đóng mở các
hệ thống cung cấp điện (ví dụ như rơle, công tắc tơ, bộ biến đổi thyristor…) được thực hiện rất
nhanh vì các thiết bị đóng cắt có qn tính nhỏ. Điều này hồn tồn có thể thực hiện được
trong thực tế. Khâu chuyển đổi biểu thị mối quan hệ giữa cơng suất cung cấp cho lị và nhiệt
độ lị, nó là một khâu qn tính bậc nhất có trễ, trong đó khâu trễ e  s đã được xấp xỉ bằng một
khâu quán tính bậc nhất theo xấp xỉ Taylor.
Ngoài ra, [10] xét đối với trường hợp khâu qn tính bậc nhất có trễ có thời gian trễ ()
là khá nhỏ so với hằng số thời gian (T) của nó, cụ thể là tỷ số T/ thỏa mãn điều kiện T/ > 10
[7]. Sau khi đưa thêm vào khâu chuyển đổi, hàm điều khiển tối ưu cần tìm là p*(t) chính là

cơng suất cung cấp cho lị chứ khơng phải là nhiệt độ lị v(t). Như vậy, dù hàm điều khiển tối
ưu p*(t) có dạng bang-bang (dạng xung vng) tức là có dạng biến thiên nhảy cấp thì hồn tồn
có thể thực hiện được vì qn tính của các phần tử điện là rất nhỏ so với các phần tử nhiệt.
Nội dung luận án [10] đã giải quyết được một số vấn đề chính như sau:
- Xét với cơng nghệ gia nhiệt cho các lị nung phôi cán được cung cấp năng lượng bằng
việc đốt nguyên liệu là dầu nặng FO. Việc điều chỉnh công suất cung cấp cho lò là điều chỉnh
lưu lượng dầu để phối hợp với lượng khơng khí trong q trình đốt.
- Xét với đối tượng có trễ nhỏ, cụ thể là đối tượng có tỷ số T/ thỏa mãn điều kiện
T/  10 [7], khâu trễ e  s được thay thế bằng khâu quán tính bậc nhất theo xấp xỉ Taylor.
- Đã giải được bài toán điều khiển tối ưu cho hệ với tham số phân bố, có trễ khi xét các
hệ số a, , k của lò và vật nung là hằng số.
- Chưa đề cập đến phần phi tuyến, cụ thể là chưa giải bài toán điều khiển tối ưu khi xét
các hệ số a, , k là phi tuyến (thực tế các hệ số này luôn thay đổi theo nhiệt độ của mơi trường
khơng khí trong lị nung, tức là chúng có tính phi tuyến).
1.3. Những vấn đề cần tiếp tục nghiên cứu về điều khiển tối ưu cho hệ với tham số phân
bố, có trễ, phi tuyến và hướng nghiên cứu của luận án
Một số vấn đề tồn tại cần được tiếp tục nghiên cứu hoàn thiện:
Cho đến thời điểm này, tác giả vẫn chưa tìm thấy được nhiều cơng trình nghiên cứu về
điều khiển tối ưu cho hệ với tham số phân bố, có trễ, phi tuyến, đặc biệt hệ này được áp dụng
cho hệ thống truyền nhiệt một phía trong lị điện trở để điều khiển nhiệt độ cho vật nung có
dạng tấm phẳng theo tiêu chuẩn nung chính xác nhất (áp dụng cho một số công nghệ như ủ vật
liệu từ, tôi ram nhiệt luyện các chi tiết cơ khí,…). Ngồi ra, hiện nay cũng chưa có nhiều cơng
trình khoa học ở trong và ngồi nước đưa ra một cách chính xác biểu thức tốn học mơ tả các
hệ số a, , k trong phương trình truyền nhiệt là phi tuyến, các hệ số này chủ yếu được xác định
gần đúng thông qua thực nghiệm.
Hướng nghiên cứu mới của luận án là:
- Thành lập 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


5

- Nghiên cứu đối tượng có trễ khá lớn, cụ thể là đối tượng có tỷ số T/  thỏa mãn điều kiện
6  T/ < 10 [7]. Thay thế khâu trễ e  s trong khâu quán tính bậc nhất, có trễ (lị điện trở)
bằng phép xấp xỉ Pade bậc một.
- Phân tích tính phi tuyến của hệ số truyền tĩnh k của lò theo nhiệt độ bằng cách nhận dạng lị
điện trở thực tế.
- Sau đó tiến hành 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ố. Hệ này được áp dụng cho hệ thống truyền nhiệt một phía trong lị điện
trở để điều khiển nhiệt độ cho vật nung có dạng tấm phẳng (xét với vật dầy) theo tiêu chuẩn
nung chính xác nhất.
- Lập các chương trình tính và xây dựng các thuật tốn, sau đó tiến hành mơ phỏng trên
Matlab để kiểm tra lời giải của bài toán điều khiển tối ưu.
- Tiến hành thí nghiệm trên mơ hình hệ thống thực là lị điện trở và vật nung để kiểm tra các kết
quả mô phỏng.
1.4. Kết luận chương 1
Chương 1 của luận án đã tập trung nghiên cứu 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; đề cập và phân tích các cơng trình, các bài báo của các tác
giả trong và ngoài nước xung quanh vấn đề này. Căn cứ vào việc phân tích, tổng hợp và nghiên
cứu kỹ lưỡng những vấn đề mà các nhà nghiên cứu đã đưa ra, cuối chương tác giả đã đưa ra
những vấn đề cần tiếp tục nghiên cứu về điều khiển tối ưu cho hệ với tham số phân bố, có trễ,
phi tuyến và hướng nghiên cứu mới của luận á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
2.1. Thành lập bài toán điều khiển tối ưu
2.1.1. Mơ hình đối tượng
Q trình đốt nóng một phía cho vật nung có dạng tấm phẳng trong lị điện trở được mơ
tả bằng phương trình vi phân đạo hàm riêng [10], [79].
a

2
�
q ( x, t ) �
q ( x, t )

2
�
t
�
x

(2.1)

trong đó: q(x,t) là phân bố nhiệt độ trong vật nung, phụ thuộc vào tọa độ không gian x với
(0  x  L) và thời gian t với (0  t  tf ), L là bề dầy của vật (m), tf là thời gian nung cho phép
(s), a là hệ số dẫn nhiệt độ (m2/s).
Các điều kiện đầu và điều kiện biên [10,79]:
q ( x, 0)  q0 ( x)  const
(2.2)


�
q ( x, t )
�
x

x 0

   q (0, t )  v (t ) 

(2.3)


6
�
q( x, t )
�
x

0

(2.4)

xL

Với  là hệ số trao đổi nhiệt giữa mơi trường khơng khí trong lị và vật (w/m 2.độ); v(t) là
nhiệt độ của mơi trường khơng khí trong lò (0C); q(0,t) là phân bố nhiệt độ tại bề mặt vật;
q0(x) là hàm phân bố nhiệt độ ban đầu của vật (hằng số, coi như bằng nhiệt độ môi trường);
 là hệ số dẫn nhiệt của vật (W/m.độ).
Nhiệt độ của mơi trường khơng khí trong lị v(t) là đại lượng trung gian được điều khiển
bởi đầu vào là điện áp cung cấp u(t), phân bố nhiệt độ trong vật q(x,t) được điều khiển thông qua
nhiệt độ của môi trường khơng khí trong lị v(t), nhiệt độ v(t) này lại được điều khiển bởi điện áp
u(t). Như vậy, thực chất sự phân bố nhiệt trong vật nung q(x,t) sẽ phụ thuộc vào điện áp cung cấp
u(t).
Quan hệ giữa điện áp cung cấp cho lò u(t) và nhiệt độ lò v(t) thường gặp là một khâu
qn tính bậc nhất, có trễ theo phương trình [6,8,10,79]:
T .v&(t )  v(t )  k .u (t   )
(2.5)
trong đó: T là hằng số thời gian của lò (s),  là thời gian trễ của lò (s), k là hệ số truyền tĩnh của

lò (hằng số), u(t) là điện áp cung cấp cho lò (hàm điều khiển của hệ thống).
Tuy nhiên, trong biểu thức (2.5), lúc này k là hệ số thay đổi phụ thuộc vào nhiệt độ
trong lò, tức là k là một hàm số theo nhiệt độ v, khi đó hệ số truyền tĩnh có thể được biễu qua
phương trình: k  k (v) , do đó k là một hệ số phi tuyến.
Thực tế qua việc nhận dạng lò điện trở thì k thay đổi khá nhiều, ví dụ như trong lò điện
trở với dải nhiệt độ thay đổi từ 0-5000C. ( Việc này sẽ được chứng minh ở phần sau).
Lúc đó quan hệ giữa u(t) và v(t) có thể được biểu diễn theo phương trình:
.v&
(t )  v (t )  k .u (t   )
(2.6)
với k là hệ số truyền tĩnh phi tuyến của lò.
Như vậy, đối tượng điều khiển (lò điện trở - vật nung) được mơ tả bằng phương trình vi
phân đạo hàm riêng bậc hai-dạng parabolic (2.1) với các điều kiện đầu và điều kiện biên (2.2),
(2.3), (2.4) kết hợp với phương trình vi phân thường, có trễ, phi tuyến (2.6). Có thể thấy đây là
một dạng bài tốn điển hình của một hệ với tham số phân bố, có trễ, phi tuyến.
Tuy nhiên, khi hệ số k là phi tuyến thì rất khó tìm được lời giải và khơng áp dụng được
phép biến đổi Laplace. Vì vậy, luận án sẽ thực hiện tuyến tính hóa hệ số k thành N giá trị là:
k1 , k2 , k3 ,..., k N . Giả thiết, k1 , k 2 , k3 ,..., k N là các hằng số.
2.1.2. Phiếm hàm mục tiêu
Bài toán điều khiển tối ưu được đặt ra như sau: tìm một hàm điều khiển u(t), với
0  t  tf sao cho cực tiểu hoá sai lệch nhiệt độ giữa phân bố nhiệt độ mong muốn q*(x) với
nhiệt độ thực của vật tại thời điểm t = tf cho trước q(x,tf ), tức là hàm mục tiêu:
L

�

2

Jc  �
q * ( x )  q ( x, t f ) �

�
�dx  min
0

(2.7)


7

Trong hàm mục tiêu cần cực tiểu thì q*(x) là phân bố nhiệt độ mong muốn (cho trước),
còn q(x,tf) là hàm chưa biết. Rõ ràng hàm q(x,tf) là giá trị của hàm q(x,t) tại thời điểm t=tf, được
hiểu là cuối quá trình gia nhiệt đảm bảo sự đồng đều nhiệt độ nhất trong tồn bộ vật nung. Bài
tốn loại này được gọi là bài tốn nung chính xác nhất. Hàm q(x,t) là nghiệm của phương trình
vi phân (2.1) với các điều kiện đầu và điều kiện biên (2.2), (2.3), (2.4), hàm q(x,t) sau khi tính
được chắc chắn phải phụ thuộc vào điện áp cung cấp cho lò u(t).
2.1.3. Điều kiện ràng buộc
2.1.4. Các bước giải
Quá trình tìm lời giải tối ưu gồm hai bước sau:
- Bước 1: Tìm quan hệ giữa phân bố nhiệt độ trong vật nung q(x,t) và tín hiệu điều khiển
điện áp u(t). Đây chính là việc giải phương trình truyền nhiệt (phương trình đạo hàm riêng
Parabolic) với điều kiện biên loại 3 (cho biết quy luật trao đổi nhiệt giữa bề mặt của vật với
môi trường xung quanh và nhiệt độ của môi trường xung quanh) kết hợp với phương trình vi
phân thường có trễ, phi tuyến (quan hệ giữa u(t) và v(t))
- Bước 2: Phân bố nhiệt độ trong vật nung q(x,t) tìm được sẽ phụ thuộc vào hàm điều
khiển u(t). Thay q(x,t) tìm được ở bước 1 vào phiếm hàm mục tiêu (2.7), sau đó dùng phương
pháp số để thay cho việc cần cực tiểu một phiếm hàm thành việc cực tiểu một hàm nhiều biến
để tìm ra nghiệm tối ưu u*(t).
2.2. Giới thiệu phương pháp xấp xỉ Pade
2.2.1. Đặt vấn đề
2.2.2. Phương pháp xấp xỉ Pade

Xét một đối tượng có trễ dạng e  s được khai triển thành chuỗi lũy thừa:

e

 s


 s 

�
!
�

(2.17)

 0

- Với r = 1, ta có xấp xỉ Pade bậc một (Pade 1):
2  s
e  s �
2  s

(2.21)

- Với r = 2, ta có xấp xỉ Pade bậc hai (Pade 2):
12  6 s   2 s 2
e  s �
12  6 s   2 s 2

(2.22)

- Với r = 3,4,...ta có xấp xỉ Pade với bậc cao hơn.

với r là số bậc cần thay thế;  là thời gian trễ của đối tượng.
2.3. Phương pháp tính gần đúng tích phân xác định
2.4. Nhận dạng mơ hình lị điện trở
2.4.1. Mơ hình lị điện trở
Input

Lị điện trở

Điện áp
Hình 2.5. Mơ hình lị điện trở
2.4.2. Hàm truyền lò điện trở

Output
Nhiệt độ


8

Theo Ziegler – Nichols thì mơ hình lị điện trở có thể được biểu diễn dưới dạng hàm
truyền là một khâu qn tính bậc nhất có trễ như sau [6]:
W( s) 

V (s)
k .e  s

U ( s ) (s  1)

(2.35)

trong đó: T là hằng số thời gian của lò (s);  là thời gian trễ của lò (s); k là hệ số khuếch đại
(hệ số truyền tĩnh) của lò; V(s) là nhiệt độ của lò (0C); U(s) là điện áp đặt vào lị (V).
Nhận dạng mơ hình lị điện trở: Sơ đồ khối hệ thống thí nghiệm để nhận dạng như hình 2.7

Hình 2.7. Sơ đồ hệ thống thu thập dữ liệu.
Đặt vào lò một điện áp dạng bước nhảy: u(t) = 220.1(t) ta thu được ở đầu ra đáp ứng
nhiệt độ trong lị như hình 2.8:

Hình 2.8. Đáp ứng nhiệt độ của lò với u = 220.1(t)
Nhận xét: Từ kết quả hình 2.8, ta thấy: Đặc tính của lị có dạng của một khâu qn tính bậc
nhất có trễ, các hằng số thời gian của lị được xác định như sau: T1200 (s); τ 130 (s)
- Hệ số truyền tĩnh k của lò: k  v / U �500 / 220 �2, 27
f

f

trong đó v là nhiệt độ đặt, U là điện áp cung cấp cho lò
- Khi nhiệt độ lị v(t) bằng vf thì k  const , tuy nhiên v(t) phụ thuộc vào thời gian t, khi v(t) thay
đổi từ nhiệt độ môi trường v0 đến nhiệt độ đặt vf thì k cũng thay đổi phụ thuộc vào v(t), tức là
k �const . Để xác định chính xác hệ số k tại mỗi thời điểm t là rất khó khăn
Mục tiêu của luận án là tìm lời giải cho bài tốn (2.1), (2.2), (2.3), (2.4) và (2.6) khi
xét hệ số truyền tĩnh

k

của lò là phi tuyến và coi các hệ số a,  của vật là hằng số. Để phân

tích sự thay đổi hệ số k theo nhiệt độ v(t), về lý thuyết có thể thực hiện như sau:

9

- Giữ nhiệt độ đặt vf = const, gọi các khoảng thay đổi điện áp là Δu (V), các khoảng thay đổi
của nhiệt độ lò là Δv (0C). Đặt vào lò các điện áp thay đổi theo dạng bậc thang, sau một
khoảng thời gian Δt, điện áp tăng lên một lượng Δu cho tới điện áp u = 220 (V), khi đó hệ số
truyền tĩnh k của lị ứng với mỗi khoảng thời gian Δt có thể được tính:
(i  1, 2, 3,..., N )

ki  vi / ui

(2.36)

Từ (2.36), ta thấy với nhiệt độ đặt vf và điện áp hiệu dụng 220V, ứng với mỗi cặp giá trị
(Δu, Δv) ta sẽ có N giá trị k , do đó khi áp dụng để tìm lời cho bài tốn thì khối lượng tính tốn
sẽ rất lớn. Để đơn giản hóa lời giải cho bài toán, tác giả đã thực hiện như sau:
Không cung cấp trực tiếp điện áp 220V mà đặt vào lò các điện áp thay đổi theo dạng bậc
thang, sau mỗi khoảng thời gian Δt = 4500 (s), nếu điện áp tăng lên một lượng Δu thì nhiệt độ
lò sẽ tăng tương ứng một lượng là Δv cho tới điện áp U = 220 (V) thì nhiệt độ lò sẽ đạt tới
nhiệt độ đặt vf =500 (0C), thời gian thí nghiệm là t = 13500 (s). Thực tế ta sẽ có vơ số cặp giá
trị (Δu, Δv) tương ứng sẽ có nhiều giá trị k , qua nhận dạng lị điện trở thực tế, nhận thấy k có
thể chỉ cần 3, 4 hoặc 5 giá trị đã đảm bảo độ chính xác (nhiệt độ đầu ra sẽ đạt nhiệt độ đặt là
vf 5000C). Do đó trong khoảng nhiệt độ đặt vf, để đơn giản hóa lời giải cho bài toán tối ưu về
sau, luận án chỉ xét với 3 giá trị của hệ số k , đó là k   k1 , k2 , k3  . Ta có đường cong thực
N h ie t d o (o C )

nghiệm như hình 2.9.
400
200

D ie n a p ( V )

0

0

2000

4000

0

2000

4000

6000

8000

10000

12000

6000

8000

10000

12000

200
150
100
50
0

Thoi gian (s)

Hình 2.9. Kết quả thực nghiệm nhận dạng mơ hình lị điện trở
Nhận xét: Từ đường đặc tính thực nghiệm hình 2.9 và cơng thức (2.36), ta thấy hệ số truyền
tĩnh k của lị khơng phải hằng số mà thay đổi phụ thuộc vào nhiệt độ trong lò. Việc chia hệ số
truyền tĩnh k của lò ra làm 3 giá trị là kết quả sau khi nhận dạng lò điện trở thực tế. Để xác định
các hệ số k như hình 2.9, ta có sơ đồ hình 2.10.


10

Hình 2.10. Đáp ứng nhiệt độ của lị để xác định các Δv .
Từ hình 2.9 và hình 2.10, ta xác định được các khoảng thay đổi nhiệt độ Δv và các
khoảng thay đổi điện áp Δu, từ đó tính được các hệ số truyền tĩnh k của lò ứng với mỗi khoảng
thời gian Δt như sau:
Bảng 2.3. Bảng xác định hệ số truyền tĩnh k
STT

Δv (0C)

Δu (V)

k �v u

1

350

185

1,8

2

100

30

3,3

3

25

5

5

Từ hình 2.10 và bảng 2.3, ta thấy hệ số truyền tĩnh k của lò sẽ tăng khi nhiệt độ trong
lò tăng lên theo thời gian. Khi nhiệt độ lò thay đổi từ nhiệt độ môi trường đến khoảng 500 ( 0C)
thì giá trị của k thay đổi khá lớn.

Như vậy, qua việc phân tích ở trên, ta thấy hệ số truyền tĩnh k thay đổi phụ thuộc vào
nhiệt độ lò v(t) mà nhiệt độ lò v(t) lại thay đổi theo thời gian t tức là hệ số truyền tĩnh k cũng
thay đổi theo thời gian t . Trong khoảng nhiệt độ cho trước, chính sự tuyến tính hóa hệ số k ra
làm 3 giá trị là k1 ; k2 ; k3 sẽ làm cho lời giải của bài toán trở nên đơn giản hơn mà độ chính xác
có thể chấp nhận được (kết quả của lời giải sẽ được chứng minh ở nội dung chương 3)..
2.5. Lời giải của bài toán tối ưu
Luận án đề xuất phương án giải bài toán tối ưu cho hệ trên như sau: chia khoảng thời
gian nung vật từ 0÷tf ra làm 3 khoảng thời gian bằng nhau Δt1 = Δt2 = Δt3 và gọi:
+ Δt1 = 0÷t1 ứng với khoảng nhiệt độ lị thay đổi từ v0÷v1 là Δv1, ta có hệ số k1
+ Δt2 = t1÷t2 ứng với khoảng nhiệt độ lị thay đổi từ v1÷v2 là Δv2, ta có hệ số k2
+ Δt3 = t2÷tf ứng với khoảng nhiệt độ lị thay đổi từ v2÷vf là Δv3, ta có hệ số k3
với: t  t / 3; t  2t / 3; tf là thời nung cho phép; v0 là nhiệt độ môi trường, vf là nhiệt độ xác
1
f
2
f
lập (nhiệt độ yêu cầu). Cụ thể xét trong khoảng thời gian Δt1 =0÷t1 như sau:


11

2.5.1. Tìm quan hệ giữa q1(x,t) và tín hiệu điều khiển u1(t)
Để tìm quan hệ giữa q1(x,t) và u1(t), ta dùng phép biến đổi Laplace thuận đối với tham
số thời gian t (khi áp dụng phép biến đổi Laplace với tham số thời gian t thì phương trình vi
phân đạo hàm riêng (2.1) đã được đưa về phương trình vi phân thường đối với biến x), sau đó
dùng phương pháp số để đưa ra lời giải cho quá trình truyền nhiệt.
Để giải phương trình đạo hàm riêng (2.1) với các điều kiện đầu và điều kiện biên (2.2),
(2.3), (2.4), áp dụng phép biến đổi Laplace đối với tham số thời gian t, được phương trình:
2
�

Q1 ( x, s)
a
 sQ1 ( x, s )
�
x2

(2.38)

2
�
Q1 ( x, s ) s
 Q1 ( x, s )
�
x2
a

hoặc

(2.39)

trong đó: Q1 ( x, s )  L q1 ( x, t )
Sau khi biến đổi các điều kiện biên (2.3), (2.4), ta được:
�
Q ( x, s )
 1
   Q1 (0, s )  V1 ( s ) 
�
x
x 0

�
Q1 ( x, s )
�
x

(2.40)

0

(2.41)

xL

Từ phương trình (2.6), dùng phép biến đổi Laplace, trong đó ta coi khâu trễ e -s được thay thế
gần đúng bằng một khâu xấp xỉ Pade bậc một, ta được:
s
1
2
(Ts  1)V1 ( s )  k1.U1 ( s ).e  s �k1.U1 ( s )
(2.42)
s
1
2
Trong đó: V1 (s)  L v1 (t ) ; U1 ( s )  L u1 (t )
Nghiệm tổng quát của (2.39) là:

Q1 ( x, s )  A1 ( s ).e

s
.x

a

 B1 ( s ).e



s
.x
a

(2.43)

Giải phương trình (2.43), cuối cùng ta được :

Q1 ( x, s )   .V1 ( s).

e
s �
. . �
e
a �
�

s
.L
a



s

.L
a

s
.x
a

s
.L
a



e
.e
s
�
.L �
a
e
e
�  . �
�
�
� �
.e

s
.x
a

s
.L
a

e

s
.L
a

�
�
�
�

(2.58)

Từ (2.42), ta có:

s
2
V1 ( s )  k1.U1 ( s )
s�
1 �
 Ts  1 �
�
� 2 �
1

Thay (2.59) vào (2.58), ta được:

(2.59)


12
�
� s �
�
k1. �
1
e
�
�
� 2 �
�
Q1 ( x, s )  U1 ( s )
�
�
� s . L
�
� s �
�
Ts

1
1

e a e

�

�
�
�
� 2 �
�
�
�
�

s
.( L  x )
a

s
.L
a

e

s
.( L  x ) �
a
�

�
�

s
�
� . a


�
�

�
.�
e
�
�

s
.L
a

e

s
.L
a

�
�
�
�
�
�
�
�
�
�

(2.61)

Đặt:

G1 ( x, s ) 

�
� s �
�
k1.�
1
e
�
�
� 2 �
�
�
�
�
�
� s �
�
Ts

1
1

e


�
�
�
�
� 2 �
�
�
�
�

s
.L
a

s
.( L  x )
a

e

s
.( L  x )
a

�
�
�
�

s

�
 . a . �
e
 �
�
�
�

s
.L �
a �

e

s
.L
a

�

e

s
�
.L �
�
a �

(2.62)

�
�
�
�
�

Từ (2.61) và (2.62), suy ra: Q1 ( x, s )  G1 (x, s).U1 ( s)
(2.63)
Như vậy ta đã xây dựng được mối quan hệ giữa tín hiệu điều khiển điện áp U1(s) và
phân bố nhiệt độ Q1(x,s) dưới dạng toán tử.
Từ (2.63), theo định lý về tích chập [8,10,12,79], ta có: q1(x,t) = g1(x,t)* u1(t)
Vậy ta có thể viết:

q1 ( x, t ) 

�

g ( x, ).u (t   ) d
�
1

1

(2.64)

�

q1 ( x, t ) 

hoặc:

�

g ( x, t   ).u ( ) d
�
1

1

(2.65)

�

trong đó:

1

g1 ( x, t )  L

 G1 ( x, s)

(2.66)

Vì vậy, nếu ta biết được hàm g1(x,t) ta sẽ tính được phân bố nhiệt q1(x,t) từ hàm điều
1
khiển u1(t). Ta sẽ đi tìm hàm g1 ( x, t )  L  G1 ( x, s) . Để tìm hàm gốc g1(x,t), ta áp dụng

cơng thức biến đổi ngược [8,10,14,79]. Sau khi biến đổi, cuối cùng ta được:
 Hàm g1(x,t) theo Pade 1:






�k
�
k1.k02 2   k02 .cos � 0 ( L  x) �
�a
� .e  k02t 
g1 ( x, t ) 
� �k0 L �  k0
�
�k L �
2   k02 �
cos � �
sin �0 �
�
� �a� a �a�
�



�k
�
2k1.k12 .cos � 1  L  x  �
2
�a
�


.e  k1 t 
� �k1L �  k1
�
k L�
�
1  Tk12 �
cos � �

sin �1 �
�
�a�
� �a� a
�






13





�
i 2



�
�
2 k1 2   . i 2 cos � i ( L  x) �
�a
�

�

.e

2

 i t

�   L � i .L � L � i .L �
�
 1  T  i 2 2   . i 2 �
sin �

c
os
�
� �
�
 . i a � a � a
� a �
�
�







(2.94)


Để tìm phân bố nhiệt độ q2(x,t) với t nằm trong khoảng Δt2 =t1÷t2 và phân bố nhiệt độ
q3(x,t) với t nằm trong khoảng Δt3 =t2÷tf , ta cũng biến đổi tương tự như trường hợp tìm q1(x,t),
cuối cùng ta cũng được kết quả như sau:
 Hàm g2(x,t) theo Pade 1:





�k
�
k2 .k02 2   k02 .cos � 0 ( L  x ) �
�a
� .e k02t 
g 2 ( x, t ) 
� �k0 L �  k0
�
�k L �
2   k02 �
cos � �
sin �0 �
�

� �a� a �a�
�





�k
�
2k2 .k12 .cos � 1  L  x  �
2
�a
�

.e  k1 t 
� �k1L �  k1
�
k L�
�
1  Tk12 �
cos � �

sin �1 �
�
�a�
� �a� a
�









�
i 2



�
�
2 k2 2   . i 2 cos � i ( L  x) �
�a
�

�

.e

�   L � i .L � L � i .L �
�
 1  T  i 2 2   . i 2 �
sin �
� cos � �
�
 . i a � a � a
� a �
�

�







2

 i t

(2.95)


 Hàm g3(x,t) theo Pade 1:





�k
�
k3.k02 2   k02 .cos � 0 ( L  x) �
�a
� .e  k02t 
g3 ( x, t ) 
� �k0 L �  k0
�
�k L �

2   k02 �
cos � �
sin �0 �
�
�a�
� �a� a
�





�k
�
2k3 .k12 .cos � 1  L  x  �
2
�a
�

.e  k1 t 
� �
�
k L �  k1
k L�
�
1  Tk12 �
cos �1 �

sin �1 �
�

�a�
� �a�  a
�





2 k3 2   . i

�



�
i 2





 1  T i 2

2


�
 cos �
� ( L  x) �
�a

�
i

.e

�
 L
� .L � L � .L �
sin �
  2   .  �
�
� cos �
�
�
 . a
� a �a
� a �
2

i

i

i

� i
�


2

 i t

(2.96)


14

k0  1 / T ; k1 

2 /  ; k1; k2 ; k3 là các hệ số truyền tĩnh của lò ứng với các khoảng thời

gian Δt1; Δt2; Δt3. Trường hợp khâu trễ e s được thay thế bằng phép xấp xỉ Taylor thì phương
trình (2.42) trở thành:
U ( s)
(Ts  1)V1 ( s)  k1.U1 ( s ).e  s �k1 1
(2.42)’
 s 1
Để tìm các hàm gμ(x,t) (   1 �3 ) ứng với 3 trường hợp của hệ số k , ta biến đổi tương tự như
trường hợp Pade 1, cuối cùng ta cũng được kết quả các hàm gμ(x,t) theo khai triển Taylor.
Trong các biểu thức (2.94), (2.95) và (2.96) các i được tính từ cơng thức:
 i  i a / L ; với i là nghiệm của phương trình:  .tg   L /   Bi ; Bi là hệ số BIO của

vật liệu;  là hệ số truyền nhiệt từ khơng gian lị vào vật (w/m2.độ);  là hệ số dẫn nhiệt của vật
cần gia nhiệt (w/m.độ); L là bề dày vật nung (m); a là hệ số dẫn nhiệt độ (m2/s);  là thời gian
trễ của lò (s); T là hằng số thời gian của lò (s)
2.5.2. Tìm lời giải cho hàm phân bố trường nhiệt độ q(x,t)
Tại mỗi thời điểm t (0 ≤ t ≤ tf ), hàm q(x,t) được tính tương ứng với 3 trường hợp:
+ Nếu 0 ≤ t ≤ t1 thì:
t

�

q ( x, t )  g1 ( x, t   ).u ( ) d

(2.116)

0

+ Nếu t1 ≤ t ≤ t2 thì:
t1

t

�

�

q ( x, t )  g1 ( x, t   ).u ( )d  g 2 ( x, t   ).u ( )d
0

(2.117)

t1

+ Nếu t2 ≤ t ≤ tf thì:
t1

t2

0

t1

�

t

�

�

q ( x, t )  g1 ( x, t   ).u ( )d  g 2 ( x, t   ).u ( )d  g 3 ( x, t   ).u ( )d

(2.118)

t2

o Kết luận: Ta đã giải được một hệ thống gồm phương trình vi phân đạo hàm riêng dạng
Parabolic với điều kiện biên loại 3 (quan hệ giữa v(t) và q(x,t)) kết hợp với phương trình vi
phân thường có trễ, phi tuyến (quan hệ giữa u(t) và v(t)). Như vậy, nếu chưa quan tâm tới bài
tốn tối ưu thì ta có thể tính được trường nhiệt độ trong vật nung khi biết điện áp cung cấp cho
lị (bài tốn biết vỏ tìm lõi).
Trường hợp tổng quát: Quan hệ giữa điện áp cung cấp cho lò u(t) và phân bố trường
nhiệt độ trong vật nung q(x,t) được tính theo các cơng thức (2.116), (2.117) và (2.118) tương
ứng với 3 miền thời gian phân chia theo hệ số ki (i=1,2,3).
2.5.3. Lời giải bài toán điều khiển tối ưu
2.5.3.1. Đặt bài tốn
Sau khi tìm được quan hệ giữa q(x,t) và u(t) dưới dạng phương trình tích phân như ở
mục 2.5.2, bài toán được đặt ra: Hãy xác định hàm điều khiển tối ưu u*(t) với (0  t  tf ) sao

cho làm cực tiểu phiếm hàm mục tiêu:
L

�

2

Jc  �
q * ( x )  q ( x, t f ) �
�
�dx
0

(2.119)


15
*

trong đó q (x) là phân bố nhiệt độ cho trước còn q(x,tf ) là phân bố nhiệt độ trong vật nung tại
thời điểm cuối quá trình nung t = tf.
Thay t = tf vào công thức (2.118) sẽ được hàm q(x,tf ):
t1

t2

tf

0

t1

t2

�

�

q( x, t f )  g1 ( x, t f   ).u ( )d  g 2 ( x, t f   ).u ( )d 

g ( x, t
�
3

f

  ).u ( )d

(2.120)

trong đó các hàm g1(x,t), g2(x,t) và g3(x,t) được tính từ các cơng thức (2.94), (2.95) và (2.96),
với tf là thời gian nung cho phép tính bằng giây (s).
Thay (2.120) vào (2.119) sẽ có dạng thức của Jc.
2.5.3.2. Tìm tín hiệu điều khiển tối ưu u*(t) bằng phương pháp số
Để tìm u*(t) ta phải cực tiểu hoá phiếm hàm (2.121):
L

2

�

Jc  �
q * ( x )  q ( x, t f ) �
�
�dx 
0

2

tf
t1
t2
�
�
�
�
�
�dx
�
 q * ( x)  g1 ( x, t f   ).u1 ( )d  g 2 ( x, t f   ).u 2 ( )d  g 3 ( x, t f   ).u3 ( )d �
�
�
�
�
0�
0
t1
t2
�
�

�
L

�

�

�

�

(2.121)
Trước hết ta dùng phương pháp tích phân số [10,11,13,15], áp dụng cơng thức Simson
đối với tích phân vế trái của phiếm hàm (2.121). Khoảng không gian là bề dày tấm từ 0 đến L
ta chia làm n phần bằng nhau (n là một số chẵn).
Lúc đó ta có thể biểu thị hàm mục tiêu Jc như sau:
n

2

J c [u*]  L �i �
q *( xi )  q ( xi , t f ) �
�
�

(2.122)

i 0

trong đó: i là các trọng số gán cho giá trị của hàm dưới dấu tích phân tại điểm xi . Các giá trị xi

và trọng số i là biết trước với mỗi cơng thức tích phân. Nếu dùng cơng thức Simson, các giá
trị của xi và i được xác định như sau [10,79]:

�xi  iL / n
�
 0   n  1 / 3n
�
với i  0,1,..., n và n là một số chẵn
�






4
/
3
n
1
3
n

1
�
�
 2   4   n  2  2 / 3n
�

(2.123)

Do q(xi,tf ) trong (2.122) được xác định theo (2.120) nên để tính J c [u*] ta áp dụng một
lần nữa cơng thức tích phân số Simson và áp dụng tương tự đối với vế phải của (2.121).
Khoảng thời gian từ 0 đến tf được chia ra ba khoảng thời gian bằng nhau là 0÷t1; t1÷t2 và t2÷tf ,
trong đó: - Khoảng thời gian từ 0 đến t1 ta chia ra thành m1 khoảng bằng nhau.
- Khoảng thời gian từ t1 đến t2 ta chia ra thành m2 khoảng bằng nhau.
- Khoảng thời gian từ t2 đến tf ta chia ra thành m3 khoảng bằng nhau.
(với m1 , m2 , m3 cũng là một số chẵn).
Khi đó giá trị của q(xi,tf ) được tính :
m1

m2

j1  0

j2  0

q ( xi , t f )  t1 � j1 g1 ( xi , t f   j1 ).u1 ( j1 )  (t 2  t1 ) � j2 g 2 ( xi , t f   j2 ).u2 ( j2 ) 


16
m3

 (t f  t2 ) � j3 g 3 ( xi , t f   j3 ).u3 ( j3 )

(2.124)

j3  0

trong

đó:

các

giá

trị

 j1  j1t1 / m1
�
�
 0   m1  1/ 3m1
�
�
1  3   m1 1  4 / 3m1
�
�
 2   4   m1  2  2 / 3m1
�

của

 j1 ;  j2 ;  j3

và

 j1 ;  j2 ;  j3 được

 j2  j2 (t2  t1 ) / m2

�
�
 0   m2  1/ 3m2
�
�
1  3   m2 1  4 / 3m2
�
�
 2   4   m2 2  2 / 3m2
�

xác

định

như

 j3  j3 (t f  t2 ) / m3
�
�
 0   m3  1/ 3m3
�
�
1  3   m3 1  4 / 3m3
�
�
 2   4   m3  2  2 / 3m3
�

sau:

(2.125)

với j1  0,1,..., m1 ; j2  0,1,..., m2 ; j3  0,1,..., m3 .Đặt: c1ij1  t1. j1 .g1 ( xi , t f   j1 );

c2ij2  (t2  t1 ). j2 .g 2 ( xi , t f   j2 ); c3ij3  (t f  t2 ). j3 .g3 ( xi , t f   j3 );
u1 ( j1 )  u j1 ; u2 ( j2 )  u j2 ; u3 ( j3 )  u j3 ; u j1  u j2  u j3  u j ; q * ( xi )  qi*

(2.126)

um1  um1 1 là điểm nối điện áp điều khiển tại thời điểm t1 .
um1  m2 1  um1  m2  2 là điểm nối điện áp điều khiển tại thời điểm t2.

thay (2.124); (2.125) và (2.126) vào (2.122), ta được :
2

m3
m2
�* �m1
�
�
J c [u * ]  L �i �
qi  �
c
.
u

c
.
u


c
.
u
�
�
�
1ij1
j1
2ij2
j2
3ij3
j3 �
��
�
i 0
j

0
j

0
j

0
2
3
�
� �1
�

n

(2.127)

Ràng buộc của hàm điều khiển (giới hạn điện áp cung cấp cho lò) được viết là:
U1  uj  U2 (j = 0,1,2,…,m ) với m = m1 + m2 + m3
(2.128)
trong đó: U1 là giới hạn dưới điện áp, U2 là giới hạn trên điện áp
Như vậy, bài toán được đặt ra là hãy tìm cực tiểu của hàm (2.127) với mj+1 biến uj tuân theo
ràng buộc (2.128). Bài toán trở thành bài tốn quy hoạch bậc hai [8,10,79]. Bài tốn này có thể tìm
nghiệm đúng bằng phương pháp số sau một số hữu hạn phép lặp.
2.6. Tính tốn các giới hạn khi giải bài tốn nung chính xác nhất.
2.7. Tính tốn nhiệt độ lò v(t) và sự phân bố nhiệt độ trong vật q(x,t)
2.7.1. Đặt vấn đề
2.7.2. Tính tốn nhiệt độ lị v(t)
Nhiệt độ lị v(t) được tính như sau:
k .(l1  l2  l3 ).u ( j  1)  v( j  1).[T  (l1  l2  l3 )]
T
với j  0,1, 2...m; l1  t1 / m1 ; l2  (t2  t1 ) / m2 ; l3  (t f  t2 ) / m3 ;
v( j ) 

(2.170)

t1  t f / 3; t2  2t f / 3 ; tf là thời gian nung cho phép (s), k  kmax  k3 ;

m1 , m2 , m3 là số khoảng thời gian tương ứng với khoảng thời gian t1 ; t2 ;t3.
T là hằng số thời gian của lò (s).
Như vậy khi đã biết u*(t) ta có thể tính được v(t) từ phương trình (2.170).
2.7.3. Tính tốn phân bố nhiệt độ trong vật nung q(x,t)
Để tính q(x,t) khi biết u*(t) ta cũng dùng phương pháp số [8,10,13,15].

Phân bố nhiệt độ trong toàn bộ vật nung trong khoảng thời gian từ 0  tf được tính như sau:


17
j1

j2

 0

 0

q ( xi , t ) �j1l1 � .g1 ( xi , j1l1    ).u1 (  )  j2l2 � .g 2 ( xi , j2l2    ).u2 (  ) 
j3

 j3l3 � .g3 ( xi , j3l3    ).u3 (  )

(2.83)

 0

với t  0 �t f ; t1  t f / 3 ; t2  2t f / 3 .
2.8. Kết luận chương 2:
Chương 2 là nội dung trọng tâm (đóng góp chính thứ nhất) của luận án, chương này đã
giải quyết được một số vấn đề sau:
- Thành lập 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
- Nhận dạng mơ hình lị điện trở và phân tích tính phi tuyến của hệ số truyền tĩnh k của lò.
- Đã thay thế khâu trễ e s trong khâu qn tính bậc nhất, có trễ (lị điện trở) bằng phép xấp xỉ
Pade 1 khi đối tượng điều khiển có tỷ số T/ thỏa mãn điều kiện 6  T/ < 10 [7].
- Đã đưa ra được lời giải tường minh cho bài tốn tìm trường nhiệt độ của vật nung (dạng tấm

phẳng) khi biết điện áp cung cấp cho lò - xét cả trong trường hợp hệ số truyền tĩnh k của lò là
phi tuyến và thời gian trễ của lị lớn.
- Đã tìm ra lời giải cho bài toán điều khiển tối ưu hệ với tham số phân bố, có trễ, phi tuyến
bằng phương pháp số – Hệ này ứng dụng cho bài tốn nung chính xác nhất trong q trình
gia nhiệt.
- Đã tính tốn các điều kiện giới hạn của bài toán tối ưu cũng như đưa ra các thuật tốn để tính
nhiệt độ lị v(t) và phân bố nhiệt độ trong vật nung q(x,t).
CHƯƠNG 3
CÁC CHƯƠNG TRÌNH TÍNH TỐN VÀ CÁC KẾT QUẢ MƠ PHỎNG
3.1. Đặt vấn đề
3.2. Các chương trình tính
3.2.1 Chương trình tính các giá trị  i
3.2.2 Tính giá trị hàm g(x,t)
3.2.3 Chương trình tính hàm g1(x,t- )
3.2.4. Chương trình giải bài tốn tối ưu
3.3. Các kết quả mơ phỏng chính
3.3.1- Mơ phỏng với mẫu Samot
 = 60 (w/m2 độ);  = 0.955 (w/mđộ); a = 4.86.10-7 (m2/s); L = 0.03 (m); T = 1200 (s);

 = 130 (s); k1  1,8; k2  3,3; k3  5 ; q* = 300 0C; tf = 4200 (s); U1 = 125 V; U2 = 205 V.
Bi   .L/   60.0,03 / 0.955 �1,9 ; 6 �T /  �1200 / 130 �9, 2  10

Đây là vật dầy vì có hệ số Bi >> 0,5. Chọn số lớp không gian là n=4, số khoảng thời gian
m1 =m2=m3=m= 16. Sau khi chạy chương trình ta được kết quả như hình 3.2


18
U*(t)

v(t)

q(x,t)

q(x,tf)=q*
q*=q*=

Hình 3.2: Thực hiện chế độ nung tối ưu với mẫu Samot (sai số e  0)
3.3.2- Mô phỏng với mẫu Diatomite
 = 60 (w/m2 độ);  = 0.2 (w/mđộ); a = 3.6.10-7 (m2/s); L = 0.04 (m); T = 1200 (s);

 = 130 (s); k1  1,8; k2  3,3; k3  5 ; q* = 400 0C; tf = 4500 (s); U1 = 125 V; U2 = 220 V.
Bi   .L/   60.0,04 / 0.2 �12 ; 6 �T /  �1200 / 130 �9, 2  10

Đây cũng là vật dầy vì có Bi >> 0,5. Chọn số lớp không gian là n=10, số khoảng thời gian
m1 =m2=m3=m= 100.
Sau khi chạy chương trình ta được kết quả như hình 3.15
U*(t)

v(t)
q(x,t)

q*

q(x,tf)

Hình 3.15: Thực hiện chế độ nung tối ưu với mẫu Diatomite (sai số e  0.003)
Nhận xét: Trên hình 3.2 và hình 3.15 thì U*(t) là tín hiệu điều khiển tối ưu (điện áp tối ưu);
v(t) là nhiệt độ lò; q(x,t) là trường nhiệt độ trong vật nung (đồ thị biểu diễn nhiệt độ của bề mặt,
nhiệt độ các lớp bên trong và nhiệt độ lớp dưới cùng của vật), e là sai số của hàm mục tiêu Jc.
Tại thời điểm cuối của quá trình nung t = tf , phân bố nhiệt độ tại các lớp trong vật nung q(x,tf )
đều xấp xỉ nhiệt độ đặt q*.

Như vậy lời giải tối ưu đã được kiểm chứng.
So sánh nghiệm tối ưu khi thay thế khâu trễ e  s bằng phép xấp xỉ Taylor và phép
xấp xỉ Pade bậc một (phần mở rộng của luận án).

Kết quả mô phỏng với mẫu Samot được thể hiện ở hình 3.3 và hình 3.4.


19
U*(t)

v(t)
q(x,t)

q*
q(x,tf)

Hình 3.3: Thực hiện chế độ nung tối ưu với mẫu Samot với n=6, m1=m2=m3=40
Sai số e = 3.9104.e-08 (Theo Taylor)

U*(t)

v(t)
q(x,t)

q(x,tf)

q*

Hình 3.4: Thực hiện chế độ nung tối ưu với mẫu Samot với n=6, m1=m2=m3=40
Sai số e =5.8559.e-10

(Theo Pade 1)
Nhận xét: Kết quả mơ phỏng trên hình 3.3 và hình 3.4 ta thấy, khi khâu trễ e  s được
thay thế gần đúng bằng phép xấp xỉ Pade 1 thì sai số của hàm mục tiêu Jc (e =5.8559.e-10) nhỏ
hơn so với trường hợp khi khâu trễ được thay thế gần đúng bằng phép xấp xỉ Taylor với
(e = 3.9104.e-08), nghĩa là bài tốn tối ưu có độ chính xác cao hơn.
3.4. Kết luận chương 3
Như vậy, thông qua các kết quả mơ phỏng có thể khẳng định rằng việc giả thiết tuyến tính
hóa hệ số truyền tĩnh k trong khoảng nhiệt độ cho trước vf ra thành 3 giá trị là k1 ; k2 ; k3 và khi
áp dụng phương pháp số để giải bài toán điều khiển tối ưu như ở chương 2 thì độ chính xác của


20

lời giải có thể chấp nhận được, tức là sai lệch (e) giữa phân bố nhiệt độ tại các lớp trong vật nung
q(x,t) so với phân bố nhiệt độ cho trước q* tại thời điểm cuối quá trình nung t=tf là rất nhỏ, đặc
biệt với đối tượng nhiệt thì việc sai số nhiệt độ như trên hoàn toàn chấp nhận được. Nội dung
trong chương 3 là đóng góp chính thứ hai của luận án.
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
4.1. Giới thiệu mơ hình hệ thống thí nghiệm
Thí nghiệm nhằm mục đích kiểm tra tính đúng đắn của các thuật tốn, chương trình tối ưu
đã được tính tốn bằng mơ phỏng và để kiểm tra độ chính xác mơ hình tốn học của đối tượng
cũng như các thơng số vật lý của lò nhiệt và vật nung.
0220V
Sơ đồ khối hệ thống thí nghiệm như hình 4.2:
Cặp nhiệt đo nhiệt
độ lị

Card NI

USB 6008

uđk 05 V
Mạch pha
xung

Matlab
Simulink
(PC)

Card
NI
USB
6008

BBĐ AC/AC
một pha

Lò

Khuếch
đại

1

010 V

Khuếch
đại

Buồng đốt

2

Vật nung
Cặp nhiệt đo
nhiệt độ vật

Hình 4.2. Sơ đồ khối hệ thống thí nghiệm
4.2. Q trình thí nghiệm thực
Với lị để tiến hành thí nghiệm, qua nhận dạng lại lị điện trở như đã trình bày ở
mục 2.4 cho thấy lị là một khâu qn tính bậc nhất, có trễ, với các hệ số được xác định:
k1 �1.8 , k2 �3.3 , k3 �5 là các hệ số truyền tĩnh của lò ứng với 3 khoảng nhiệt độ v1 ;

v2 ; v3 . T1200 (s) là hằng số thời gian của lò;  130 (s) là thời gian trễ của lò

Trong quá trình thí nghiệm, ta đo nhiệt độ trong khơng gian lò và tại ba điểm của vật: điểm thứ
nhất là lớp bề mặt vật (ký hiệu là q m(t)); điểm thứ hai là lớp ở tâm vật (ký hiệu q t(t)); điểm thứ
ba là lớp dưới cùng của vật (ký hiệu qc(t). Cả bốn tín hiệu đo này được đưa qua bộ khuếch đại
để cho điện áp ra từ 0 đến 10V, sau đó đưa vào máy tính thơng qua card NI USB 6008.
Các tham số
đầu vào
(lị, vật)

Chương trình *
u (t)
tối ưu

Vật nung
Lị điện v(t)

q(x,t)
(đối tượng với tham
trở
số phân bố)

Hình 4.10: Sơ đồ khối mạch thí nghiệm điều khiển tối ưu
Trong sơ đồ khối hình 4.10 thì các tham số của lị và vật đã biết cịn chương trình tối ưu đã
được tính tốn từ lời giải của bài tốn tối ưu, chương trình này tạo ra điện áp tối ưu đặt vào lị.
Đường điện áp tối ưu u*(t) có dạng xung vng, tiếp theo là lị điện trở (đối tượng có trễ, phi
tuyến), cuối cùng là vật nung (đối tượng với tham số phân bố). Từ sơ đồ khối hình 4.10, ta có
sơ đồ mạch thí nghiệm thực như hình 4.11.


21

Hình 4.11: Sơ đồ mạch thí nghiệm thực
4.3. Một số kết quả thí nghiệm
4.3.1. Thí nghiệm với mẫu Samot
4.3.1.1. Thí nghiệm 1. Tiến hành thí nghiệm với một mẫu Samot có các thơng số sau:
L = 0,03 (m); a = 4,86.10-7 (m2/s);  = 0.955 (W/m.độ); q* = 3000C; tf =4200 (s)
U1 = 125 (V); U2 =205 (V). Kết quả thí nghiệm thực như hình 4.12
U*(t)

q*

v(t)
qm(t)
qc(t)

qt(t)

Hình 4.12: Kết quả thí nghiệm với mẫu Samot (q*=3000C)
4.3.1.2. Thí nghiệm 2.
Giữ nguyên các thông số của vật nung, thời gian nung; thay đổi thông số sau:
*
q = 4000C; U1 = 140 (V); U2 =220 (V).
Kết quả thí nghiệm như hình 4.13


22
U*(t)

q*

v(t)
qm(t)
qt(t)

qc(t)

Hình 4.13: Kết quả thí nghiệm với mẫu Samot (q*=4000C)
Trên hình 4.12 và hình 4.13 thì U*(t) là điện áp tối ưu đặt lên thanh dây đốt, q* là nhiệt
độ đặt, v(t) là nhiệt độ trong khơng gian lị; qm(t), qt(t) và qc(t) tương ứng cũng là nhiệt độ lớp
bề mặt, nhiệt độ ở tâm và nhiệt độ lớp dưới cùng của vật nung
Kết quả thí nghiệm cũng cho thấy, cuối quá trình nung (tại tf = 4200s) thì nhiệt độ tại lớp
bề mặt, lớp ở giữa và lớp dưới cùng đều xấp xỉ nhiệt độ đặt là 3000C và 4000C.
4.3.2. Thí nghiệm với mẫu Diatomite
Để kiểm tra tính đa dạng của chương trình tối ưu, ta tiếp tục tiến hành thí nghiệm với
một mẫu Diatomite có độ dầy lớn hơn mẫu Samot. Trong trường hợp này ta cũng thí nghiệm
với nhiệt độ đặt là 3000C và 4000C.

4.3.2.1. Thí nghiệm 1. Thí nghiệm với mẫu Diatomite có các thơng số sau:
L = 0,04 (m); q* = 3000C; tf = 4500 (s); U1 =125 (V); U2 =205 (V).
Kết quả thí nghiệm như hình 4.14
U*(t)

q*

v(t)
qm(t)
qc(t)

qt(t)

Hình 4.14: Kết quả thí nghiệm với mẫu Diatomite (q*=3000C)