Tổng quan nghiên cứu

Trong các ngành công nghiệp luyện kim, chế tạo máy và gia công nhiệt luyện, quá trình gia nhiệt kim loại chiếm khoảng 15% đến 25% tổng chi phí năng lượng của toàn bộ dây chuyền sản xuất. Chất lượng cơ tính của sản phẩm như độ bền kéo, giới hạn chảy, độ dẻo dai và độ cứng phụ thuộc trực tiếp vào sự phân bố trường nhiệt độ bên trong lòng phôi thép. Tuy nhiên, việc đặt cảm biến nhiệt độ vào sâu bên trong lõi phôi trong điều kiện môi trường lò nung khắc nghiệt từ 650°C đến 1200°C là bất khả thi về mặt kỹ thuật và phá hủy cấu trúc vật liệu. Bài toán "Biết vỏ tìm lõi" – tức là chỉ dựa vào nhiệt độ bề mặt ngoài và nhiệt độ không gian lò để xác định và điều chỉnh chính xác nhiệt độ tại các lớp sâu bên trong phôi – là một thách thức công nghệ then chốt cần được giải quyết.

Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình toán học giải tích dưới dạng hàm truyền đạt để quan sát trường nhiệt độ phân bố theo tiết diện phôi thép tấm đa lớp, từ đó thiết kế bộ điều khiển nhiệt độ kinh điển đạt chất lượng tối ưu mà không cần cấy cảm biến đo trực tiếp trong phôi. Phạm vi nghiên cứu tập trung vào đối tượng phôi thép tấm phẳng có chiều dài 0,5 mét, chiều rộng 0,2 mét, chiều dày 0,05 mét (50 mm), khối lượng riêng 7800 kg/m³ và hệ số dẫn nhiệt 55 W/m.K trong lò điện trở gia nhiệt một phía. Nghiên cứu mang ý nghĩa thực tiễn lớn khi giúp hệ thống đạt độ quá điều chỉnh 0%, triệt tiêu hoàn toàn dao động và rút ngắn thời gian quá độ từ 1175 giây xuống còn 950 giây (cải thiện 19,1%), nâng cao hiệu suất truyền nhiệt và tiết kiệm đáng kể điện năng tiêu thụ trong công nghiệp gia nhiệt.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng định luật truyền nhiệt Fourier và phương trình vi phân đạo hàm riêng mô tả quá trình dẫn nhiệt đẳng hướng trong vật rắn ba chiều. Bằng việc áp dụng định lý Ostrogradsky để chuyển đổi tích phân mặt kín sang tích phân thể tích, phương trình cân bằng nhiệt vi phân thuần nhất được thiết lập cho vật liệu đồng chất với hệ số khuếch tán nhiệt xác định. Đối với vật thể mỏng có tiêu chuẩn Biot nhỏ hơn 0,25, phôi kim loại được xem như một khối nhiệt dung đồng nhất một lớp. Đối với các phôi tấm dày, lý thuyết chia phôi thành $n$ lớp đẳng nhiệt liên tiếp được áp dụng, kết hợp tương đương nhiệt học - điện học trong đó điện trở nhiệt $R$ và nhiệt dung $C$ của từng lớp được xác định từ hệ số dẫn nhiệt 55 W/m.K, hệ số truyền nhiệt đối lưu 335 W/m²K và nhiệt dung riêng 460 J/kg.K.

Để thiết kế hệ thống điều khiển tự động, hai lý thuyết tổng hợp điều khiển kinh điển được áp dụng chuyên sâu:

  1. Tiêu chuẩn tối ưu độ lớn (Tiêu chuẩn phẳng - Modulus Optimum): Triệt tiêu các hằng số thời gian lớn trong mạch hở ($T_n = 556$ giây, $T_u = 500$ giây) để đưa hàm truyền hệ kín về dạng bậc hai chuẩn với hệ số suy giảm tối ưu 0,707.
  2. Phương pháp đa thức đặc trưng có hệ số suy giảm thay đổi được: Sử dụng đa thức chuẩn bậc cao với hệ số suy giảm $\alpha = 2,8$ (vượt ngưỡng giới hạn ổn định $\alpha > 1,6$) kết hợp bộ biến đổi công suất Tiristor có hàm truyền quán tính bậc nhất với hằng số thời gian 0,57 giây, đảm bảo khử triệt để sai lệch tĩnh và triệt tiêu dao động quá độ.

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

Nghiên cứu áp dụng phương pháp giải tích kết hợp mô hình hóa toán học và thực nghiệm mô phỏng trên phần mềm MATLAB/Simulink. Thay vì sử dụng phương pháp sai phân lưới phức tạp đòi hỏi tài nguyên tính toán lớn trong thời gian thực, nghiên cứu lựa chọn phương pháp hàm truyền liên tục nhờ khả năng tích hợp trực tiếp vào các cấu trúc điều khiển tự động công nghiệp.

Bộ dữ liệu nghiên cứu bao gồm 12 kịch bản mô phỏng toàn diện được phân loại theo 4 cấu trúc hình học của phôi tấm (1 lớp, 2 lớp, 3 lớp và 4 lớp). Cỡ mẫu thử nghiệm bao quát trạng thái danh định và các trạng thái biến thiên thông số nhiệt vật lý trong biên độ $\pm 15%$ nhằm kiểm tra tính bền vững của hệ thống điều khiển. Phương pháp chọn mẫu thử nghiệm được thiết kế theo nguyên tắc quét tham số có trọng số, tập trung vào sự thay đổi của hệ số dẫn nhiệt, nhiệt dung và kích thước phôi. Thời gian quan sát cho mỗi mẫu mô phỏng kéo dài liên tục 4500 giây với chu kỳ lấy mẫu 0,1 giây. Lý do lựa chọn phương pháp mô phỏng số giải tích là nhằm đánh giá chính xác động học phân bố nhiệt độ tại từng lớp phôi mà các phương pháp đo lường vật lý truyền thống không thể thực hiện được do giới hạn phá hủy của que đo nhiệt điện.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và phân tích động học hệ thống điều khiển trường nhiệt độ trong phôi tấm đã ghi nhận các kết quả định lượng nổi bật sau:

Thứ nhất, mô hình quan sát nhiệt độ dạng hàm truyền phân lớp phản ánh chính xác động học truyền nhiệt từ bề mặt vào sâu trong lõi. Với mô hình 1 lớp, hàm truyền quán tính có hằng số thời gian $T = 536,41$ giây; khi chia thành 3 lớp, hằng số thời gian của lớp trong cùng đạt 17,05 giây. Đáp ứng nhiệt độ tại các lớp diễn ra tuần tự, lớp ngoài tiếp nhận nhiệt lượng tức thời trong khi lớp thứ ba có độ trễ pha tự nhiên khoảng 120 giây trước khi tăng tốc đồng đều.

Thứ hai, bộ điều khiển PID thiết kế theo tiêu chuẩn phẳng cho hệ thống gia nhiệt phôi 3 lớp đạt thời gian quá độ danh định là 1175 giây, sai số xác lập bằng 0 và độ quá điều chỉnh đạt 0%. Khi các thông số phôi thay đổi $\pm 15%$, thời gian quá độ dao động ổn định trong khoảng từ 1060 giây (giảm 9,8%) đến 1300 giây (tăng 10,6%) mà không hề xuất hiện hiện tượng dao động hay mất ổn định.

Thứ ba, bộ điều khiển PID thiết kế theo phương pháp đa thức đặc trưng với hệ số suy giảm $\alpha = 2,8$ mang lại chất lượng vượt trội. Thời gian quá độ danh định được rút ngắn xuống còn 950 giây, nhanh hơn 225 giây (tương đương cải thiện tốc độ đáp ứng 19,15% so với tiêu chuẩn phẳng). Trong các điều kiện thông số phôi biến động mạnh, thời gian quá độ đạt giá trị cực tiểu là 760 giây và cực đại là 1075 giây, độ quá điều chỉnh luôn duy trì tuyệt đối ở mức 0%.

Thứ tư, đối với các mô hình 1 lớp, 2 lớp và 4 lớp, phương pháp đa thức đặc trưng đều duy trì thời gian xác lập dưới 970 giây (cụ thể: 1 lớp đạt 970 giây, 2 lớp đạt 960 giây, 4 lớp đạt 940 giây), chứng minh tính thích ứng cao với độ dày chia nhỏ của vật liệu.

Thảo luận kết quả

Các phát hiện nghiên cứu chứng minh rằng việc lựa chọn hệ số suy giảm $\alpha = 2,8$ trong phương pháp đa thức đặc trưng đã giải quyết triệt để sự mâu thuẫn giữa tốc độ tác động và độ quá điều chỉnh. Về mặt lý thuyết, khi $\alpha < 1,5$, hệ thống sẽ xuất hiện dao động mạnh và độ quá điều chỉnh có thể vượt quá 40%; tuy nhiên với $\alpha = 2,8$, nghiệm của phương trình đặc trưng được phân bố sâu vào nửa trái mặt phẳng phức, tạo ra đáp ứng quá độ đơn điệu không vọt lố.

Dữ liệu mô phỏng có thể được biểu diễn trực quan qua biểu đồ đường đáp ứng thời gian (từ 0 đến 4500 giây) so sánh giữa nhiệt độ đặt và nhiệt độ thực tế của từng lớp, kết hợp bảng thống kê so sánh các chỉ tiêu chất lượng (thời gian quá độ, độ quá điều chỉnh, sai số tĩnh). So với các nghiên cứu gia nhiệt truyền thống chỉ khống chế nhiệt độ không gian lò dẫn đến hiện tượng chênh lệch nhiệt độ bề mặt và lõi phôi lên tới 80°C - 100°C gây nứt nẻ ứng suất nhiệt, phương pháp điều khiển theo mô hình hàm truyền phân lớp giúp kiểm soát độ dốc nhiệt giữa các lớp luôn nằm trong giới hạn an toàn dưới 15°C trong suốt quá trình nâng nhiệt.

Đồng thời, sự tích hợp khâu lọc quán tính đầu vào với hằng số thời gian $T_{loc} = 317,6$ giây đã triệt tiêu hoàn toàn tác động của các điểm không trong hàm truyền kín, ngăn ngừa triệt để các xung nhiệt đột ngột có thể làm suy giảm tuổi thọ của dây đốt lò điện trở.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu lý thuyết và kiểm chứng mô phỏng, bốn giải pháp công nghệ trọng tâm được đề xuất nhằm tối ưu hóa quy trình nhiệt luyện công nghiệp:

Thứ nhất, chuyển đổi cấu hình điều khiển lò nung từ chế độ điều khiển nhiệt độ không gian lò sang thuật toán PID đa thức đặc trưng ($\alpha = 2,8$) tích hợp bộ quan sát hàm truyền. Kỹ sư tự động hóa tại nhà máy cần hoàn thành việc lập trình khối thuật toán trên PLC hoặc hệ thống điều khiển phân tán trong vòng 3 tháng đầu, hướng tới mục tiêu rút ngắn 15% đến 20% chu kỳ gia nhiệt mỗi mẻ phôi.

Thứ hai, chuẩn hóa mô hình phân lớp phôi tấm với số lớp $n \ge 3$ cho các chi tiết thép có chiều dày trên 0,05 mét. Phòng kỹ thuật và công nghệ vật liệu chịu trách nhiệm xây dựng thư viện thông số nhiệt vật lý (hệ số dẫn nhiệt, nhiệt dung riêng) tương ứng với từng mác thép cụ thể trong thời gian 6 tháng, nhằm loại bỏ hoàn toàn nhu cầu sử dụng que đo nhiệt phá hủy bên trong phôi, giảm 100% chi phí cảm biến tiêu hao.

Thứ ba, tích hợp khâu tiền lọc tín hiệu đặt với hằng số thời gian lọc thích nghi bằng 4 lần hằng số tích phân để khử triệt để tác động quá điều chỉnh đối với các lò nung công suất lớn trên 100 kW. Đội ngũ bảo trì và vận hành hệ thống thực hiện cài đặt và hiệu chuẩn trong thời hạn 2 tuần trước khi đưa lò vào vận hành thương mại.

Thứ tư, thiết lập quy trình tự động bù suy giảm công suất do già hóa dây đốt điện trở định kỳ 6 tháng một lần. Phòng cơ điện nhà máy cần đo đạc sự biến thiên điện trở thuần của thanh gia nhiệt để cập nhật kịp thời hệ số khuếch đại của bộ biến đổi Tiristor, giữ sai số kiểm soát nhiệt độ toàn hệ thống luôn dưới mức 1,5%.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho bốn nhóm đối tượng chuyên môn:

  1. Kỹ sư điều khiển và tự động hóa: Nắm vững phương pháp thiết kế bộ quan sát trạng thái dạng hàm truyền liên tục và kỹ thuật tổng hợp bộ điều khiển PID theo tiêu chuẩn phẳng và đa thức đặc trưng để ứng dụng trực tiếp cho các hệ thống lò điện trở công nghiệp từ 650°C đến 1200°C.
  2. Kỹ sư công nghệ luyện kim và nhiệt luyện: Hiểu rõ cơ chế truyền nhiệt và quy luật phân bố nhiệt độ theo tiết diện lá thép, từ đó thiết lập giản đồ nhiệt nung tối ưu, đảm bảo độ đồng đều tổ chức tế vi pha thép và hạn chế tối đa ứng suất dư gây nứt phôi tấm dày 50 mm.
  3. Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật điều khiển: Sử dụng luận văn như tài liệu tham khảo chuẩn mực về mô hình hóa các đối tượng có thông số phân bố thành các khâu thông số tập trung tương đương và giải pháp giải quyết bài toán "Biết vỏ tìm lõi".
  4. Giám đốc kỹ thuật và nhà quản lý sản xuất tại các nhà máy cơ khí, cán thép: Áp dụng giải pháp điều khiển tiên tiến để tối ưu hóa hiệu suất năng lượng, cắt giảm từ 10% đến 15% lượng điện tiêu thụ trong các công đoạn ủ, tôi, ram và nâng cao tuổi thọ thiết bị lò nhiệt.

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

1. Bài toán "Biết vỏ tìm lõi" trong luận văn giải quyết thách thức kỹ thuật gì? Bài toán giải quyết việc xác định chính xác nhiệt độ bên trong lõi phôi thép dày 0,05 mét mà không cần khoan đặt cảm biến nhiệt phá hủy cấu trúc. Bằng cách thiết lập hàm truyền đạt từ nhiệt độ mặt ngoài với hệ số dẫn nhiệt 55 W/m.K, hệ thống tự động suy biến trường nhiệt độ bên trong với độ chính xác cao.

2. Tại sao phương pháp đa thức đặc trưng lại ưu việt hơn tiêu chuẩn phẳng trong điều khiển nhiệt độ phôi? Phương pháp đa thức đặc trưng với hệ số $\alpha = 2,8$ rút ngắn thời gian quá độ từ 1175 giây xuống 950 giây (nhanh hơn 19,15%) so với tiêu chuẩn phẳng. Phương pháp này phân bố vị trí các điểm cực của hệ kín tối ưu hơn, giúp triệt tiêu hoàn toàn hiện tượng vọt lố và giữ độ quá điều chỉnh ở mức 0%.

3. Phôi thép tấm nên được chia thành bao nhiêu lớp khi xây dựng bộ quan sát? Đối với phôi thép tấm dày 50 mm, việc chia thành 3 lớp hoặc 4 lớp là tối ưu nhất. Mô hình 3 lớp cho phép kiểm soát chi tiết nhiệt độ từng phần dày 0,0167 mét, vừa đảm bảo độ chính xác vi phân nhiệt học vừa giữ bậc hàm truyền ở mức bậc 3, thuận lợi cho việc xử lý thời gian thực trên vi điều khiển.

4. Khi thông số nhiệt vật lý của phôi thay đổi $\pm 15%$, hệ thống có bị mất ổn định không? Hệ thống hoàn toàn không bị mất ổn định. Kết quả mô phỏng cho thấy khi thông số phôi biến thiên $\pm 15%$, thời gian quá độ chỉ dao động trong khoảng từ 760 giây đến 1075 giây, không hề xuất hiện dao động cộng hưởng hay sai lệch tĩnh, chứng minh tính kháng nhiễu và độ bền vững rất cao.

5. Mô hình hàm truyền này có thể áp dụng cho các dải nhiệt độ lò trên 1200°C không? Hoàn toàn có thể áp dụng. Khi vận hành ở dải nhiệt độ cao trên 1200°C, kỹ sư chỉ cần tái hiệu chuẩn hệ số truyền nhiệt bức xạ - đối lưu $\alpha = 335$ W/m²K và nhiệt dung riêng $c = 460$ J/kg.K tương ứng với điểm làm việc nhiệt độ mới của lò điện trở.

Kết luận

Nghiên cứu đã hoàn thành toàn diện các mục tiêu đề ra với những đóng góp khoa học và ứng dụng thực tiễn nổi bật:

  • Xây dựng thành công mô hình toán học giải tích dạng hàm truyền đạt tổng quát mô tả trường nhiệt độ phôi tấm $n$ lớp từ phương trình vi phân truyền nhiệt Fourier.
  • Giải quyết triệt để bài toán "Biết vỏ tìm lõi", cho phép quan sát liên tục nhiệt độ bên trong lòng phôi thép dày 0,05 mét trong suốt chu kỳ gia nhiệt 4500 giây.
  • Tổng hợp thành công bộ điều khiển PID theo tiêu chuẩn phẳng bù hai hằng số thời gian trội 556 giây và 500 giây, loại bỏ sai lệch tĩnh xác lập.
  • Ứng dụng xuất sắc phương pháp đa thức đặc trưng ($\alpha = 2,8$), rút ngắn thời gian quá độ xuống 950 giây (cải thiện 19,1% tốc độ đáp ứng) và duy trì độ vọt lố 0%.
  • Kiểm chứng độ bền vững của thuật toán qua 12 kịch bản mô phỏng MATLAB/Simulink dưới các điều kiện biến đổi tham số phôi $\pm 15%$.

Kế hoạch phát triển tiếp theo là triển khai thử nghiệm thuật toán trên hệ thống nhúng điều khiển lò nhiệt công nghiệp thực tế trong vòng 6 đến 12 tháng tới. Các doanh nghiệp sản xuất cơ khí và nhiệt luyện quan tâm có thể tiếp nhận chuyển giao công nghệ để tối ưu hóa quy trình nung và nâng cao chất lượng sản phẩm phôi thép.