Đo thay đổi đầu vào để tính trước phần nhiệt cần thêm
Một dòng sản phẩm đi qua thiết bị gia nhiệt có thể tăng lưu lượng trước khi cảm biến nhiệt độ ra thấy lạnh đi. Nếu đo được lưu lượng đầu vào, có thể tính phần nhiệt bổ sung ngay khi nhiễu xuất hiện. Đó là ứng dụng feedforward cho nhiễu đo được; bộ điều khiển phản hồi nhiệt độ vẫn xử lý phần sai lệch còn lại.
Cooper và Houtz mô tả chuỗi đo nhiễu, dự đoán ảnh hưởng rồi thực hiện tác động bù, đồng thời nhấn mạnh phép đo phải tạo đủ thời gian cho tác động đến biến cần giữ. [1] Việc lắp thêm một đồng hồ lưu lượng chỉ mới cung cấp thông tin. Cần quan hệ năng lượng, động học và giới hạn chấp hành để biến thông tin đó thành lệnh đúng.
Cân bằng năng lượng cho một bồn gia nhiệt giả định
Xét bồn khuấy trộn hoàn toàn, khối lượng giữ trong bồn M không đổi, dòng vào bằng dòng ra m. Chất lỏng có nhiệt dung c không đổi; nhiệt độ trong bồn và ra đều bằng T. Dòng vào có T_i, môi trường T_a, hệ số thất thoát nhiệt gộp H đơn vị W/K. Công suất điện P đi vào bộ gia nhiệt với hiệu suất chuyển thành nhiệt cho bồn η.
Mc dT/dt = mc(T_i−T)+ηP−H(T−T_a).
M dùng kg, m dùng kg/s, c dùng J/(kg·K), P dùng W, t dùng s. H là độ dẫn nhiệt tổn thất, không phải diện tích hay một hằng số thời gian. Phương trình giả định không sôi, không phản ứng và không có biến thiên tồn lượng gây thêm tích trữ.
Muốn giữ T=T_sp ở trạng thái ổn định, lệnh bù từ mô hình là P_ff=[mc(T_sp−T_i)+H(T_sp−T_a)]/η. Đây là lệnh công suất điện; không tự đồng nhất với phần trăm mở van hơi hay phần trăm dòng điều khiển.
Bước thay đổi lưu lượng và nhiệt thiếu nếu không bù
Cho M=200 kg, c=4.000 J/(kg·K), m_0=2 kg/s, T_i=T_a=25°C, T_sp=65°C, H=100 W/K và η=0,95. Nhiệt mang theo dòng là 320 kW, thất thoát 4 kW, nên P_0=341,052632 kW.
Khi lưu lượng tăng lên 3 kg/s và các đại lượng còn lại giữ nguyên, nhiệt cần vào dòng thành 480 kW. Lệnh mới là P_1=509,473684 kW, tăng 168,421053 kW điện, tương ứng 160 kW nhiệt hữu ích. Không nhân toàn bộ P_0 với 1,5 vì phần thất thoát 4 kW không tăng theo lưu lượng trong mô hình này.
Nếu công suất vẫn giữ P_0, nhiệt độ ổn định mới là T_ss=[mcT_i+ηP_0+HT_a]/(mc+H)=51,776860°C. Khi bù đúng và nhiệt cấp thay đổi tức thời, vế phải tại T=65°C bằng không nên nhiệt độ có thể giữ nguyên trong mô hình. Điều kiện tức thời sẽ được bỏ ở phần tiếp theo.
Sai số phép đo truyền trực tiếp vào phần bù
Tại m=3 kg/s, nếu cảm biến nhiệt đầu vào đọc cao hơn thực tế 1 K, thuật toán cho rằng cần ít nhiệt hơn. Nhiệt cấp bị thiếu mc×1=12 kW; sai lệch nhiệt độ ổn định là −12.000/(12.000+100)=−0,991736 K khi chưa có feedback trim.
Nếu đồng hồ lưu lượng đọc thấp 5%, tức 2,85 thay cho 3 kg/s, nhiệt cấp bị thiếu 0,15×4.000×40=24 kW. Sai lệch ổn định là −1,983471 K. Dấu của sai lệch là một kiểm tra quan trọng: đọc lưu lượng thấp trong nhiệm vụ gia nhiệt phải dẫn đến bù thiếu, không phải bù thừa.
Phép tính này cũng cho thấy vì sao không nên cộng cùng một lưu lượng đầy đủ hai lần. Nếu bộ phản hồi đang có bias cung cấp P_0, feedforward có thể dùng phần thay đổi ΔP_ff=P_ff−P_0. Nếu feedforward cung cấp toàn bộ P_ff, phần phản hồi phải được định nghĩa là hiệu chỉnh quanh mức đó. Hai cấu trúc đều có thể dùng, nhưng cần nhất quán khi chuyển chế độ.
Bù đúng trạng thái cuối vẫn có thể tạo một hõm nhiệt
Giữ bước lưu lượng từ 2 lên 3 kg/s, nhưng giả định nhiệt thực cấp vào bồn có trễ bậc nhất τ_a=15 s sau lệnh. Nhiệt hữu ích khi đó là q(t)=484.000−160.000exp(−t/15) W. Lệnh điện đã tăng đúng mức, nhưng nhiệt thực chưa tăng ngay.
Đặt θ=T−65 và τ_p=Mc/(mc+H)=66,115702 s. Với θ(0)=0, phương trình sau bước là Mc dθ/dt+(mc+H)θ=−160.000exp(−t/15). Nghiệm:
θ(t)=−[160.000/(Mc)] × [exp(−t/τ_p)−exp(−t/τ_a)]/(1/τ_a−1/τ_p).
Nhiệt độ thấp nhất ở t=ln(τ_p/τ_a)/(1/τ_a−1/τ_p)=28,779746 s, với θ=−1,941224 K. Tại 60 s, sai lệch vẫn −1,494776 K trước khi trở về không. Đây là đáp ứng của bồn và chấp hành đã nêu, chưa có vòng phản hồi nhiệt độ. Giá trị cuối đúng không chứng minh chất lượng quá độ tốt.
Từ bù tĩnh sang bù động theo hai đường truyền
Trong mô hình tuyến tính quanh một điểm làm việc, gọi y là sai lệch nhiệt độ, u là sai lệch lệnh thao tác và d là nhiễu đo được. Viết y=G_p(s)u+G_d(s)d. Nếu u=F(s)d, triệt tiêu lý tưởng cần F=−G_d/G_p. MathWorks trình bày cấu trúc này cùng việc ghép với phản hồi. [2]
Xét một bộ mô hình độc lập để minh họa động học, không phải mô hình bồn vừa tính: G_p=0,12exp(−8s)/(40s+1), đơn vị K/kW; G_d=−5exp(−20s)/(25s+1), đơn vị K/(kg/s). Khi đó F=41,666667(40s+1)/(25s+1)exp(−12s), đơn vị kW/(kg/s).
Với bước nhiễu d=0,2 kg/s, lệnh bù bắt đầu sau 12 s ở 13,333333 kW rồi giảm theo u(t)=8,333333[1+0,6exp{−(t−12)/25}] kW. Lệnh cuối là 8,333333 kW. Phần nhô đầu là tác dụng của lead–lag, không phải tích phân sai số nhiệt độ.
Nếu chỉ dùng gain tĩnh 41,666667 và vẫn đặt đúng trễ 12 s, sau thời điểm t=20 s nhiệt độ sai lệch là y=exp[−(t−20)/25]−exp[−(t−20)/40] K. Nó đạt −0,171329 K tại t=51,333575 s rồi về không. Mô hình bù động lý tưởng triệt tiêu toàn bộ đáp ứng này, nhưng chỉ trong điều kiện mô hình đúng và chấp hành thực hiện được lệnh.
Điều kiện nhân quả không thể bỏ bằng phép chia hàm truyền
Với hai đường bậc nhất có trễ θ_p và θ_d, F chứa exp[−(θ_d−θ_p)s]. Nếu θ_d nhỏ hơn θ_p, hệ số trở thành một phép đi trước thời gian, đòi hỏi giá trị nhiễu tương lai. NPTEL nêu rõ trường hợp đó không thể thực hiện trực tiếp bằng bộ điều khiển chỉ biết tín hiệu hiện tại. [3]
Trong mô hình minh họa, nếu trễ đường nhiễu giảm từ 20 xuống 5 s trong khi đường thao tác vẫn 8 s, phép bù lý tưởng cần biết trước nhiễu 3 s. Có thể xem xét đo xa hơn phía thượng lưu, giảm trễ chấp hành hoặc chấp nhận bù không hoàn hảo. Không thay exp(+3s) bằng một khối trễ âm rồi coi bài toán triển khai đã được giải.
Giới hạn công suất và vai trò phần phản hồi
Giả sử riêng bộ gia nhiệt bồn bị giới hạn ở 450 kW điện. Sau bước lưu lượng, dù thuật toán yêu cầu 509,473684 kW, mức nhiệt độ ổn định cao nhất trong mô hình là 25+0,95×450.000/12.100=60,330579°C. Không bộ bù nào tạo được công suất chưa có; phần feedback cần thấy giới hạn để tránh tích lũy lệnh vô ích.
Với feedback u=C(r−y)+Fd, đường nhiễu thành y/d=(G_d+G_pF)/(1+G_pC) khi r=0. Feedforward giảm phần nhiễu đã mô tả, feedback xử lý sai mô hình và nhiễu chưa đo trong khả năng vòng. Khi thử nghiệm, cần ghi đồng thời nhiễu thô, nhiễu đã lọc, phần bù, phần phản hồi, tổng lệnh, công suất thực và nhiệt độ. Nhờ đó mới phân biệt được bù sai gain, lệch thời gian, nhiễu cảm biến hay giới hạn thiết bị.
Nguồn tham khảo
[1] The Feed Forward Controller — Doug Cooper and Allen Houtz, Control Guru / Control Station Inc., 9 April 2015. Truy cập ngày 09/11/2026.
[2] Temperature Control in Heat Exchanger — MathWorks, Control System Toolbox Documentation, R2026b. Truy cập ngày 09/11/2026.
[3] Module 5: Other control strategies — Lecture 1: Feedforward control, page 6 — NPTEL, course 103103037. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một dòng sản phẩm đi qua thiết bị gia nhiệt có thể tăng lưu lượng trước khi cảm biến nhiệt độ ra thấy lạnh đi. Nếu đo được lưu lượng đầu vào, có thể tính phần nhiệt bổ sung ngay khi nhiễu xuất hiện. Đó là ứng dụng feedforward cho nhiễu đo được; bộ điều khiển phản hồi nhiệt độ vẫn xử lý phần sai lệch còn lại.
Cooper và Houtz mô tả chuỗi đo nhiễu, dự đoán ảnh hưởng rồi thực hiện tác động bù, đồng thời nhấn mạnh phép đo phải tạo đủ thời gian cho tác động đến biến cần giữ. [1] Việc lắp thêm một đồng hồ lưu lượng chỉ mới cung cấp thông tin. Cần quan hệ năng lượng, động học và giới hạn chấp hành để biến thông tin đó thành lệnh đúng.
Cân bằng năng lượng cho một bồn gia nhiệt giả định
Xét bồn khuấy trộn hoàn toàn, khối lượng giữ trong bồn M không đổi, dòng vào bằng dòng ra m. Chất lỏng có nhiệt dung c không đổi; nhiệt độ trong bồn và ra đều bằng T. Dòng vào có T_i, môi trường T_a, hệ số thất thoát nhiệt gộp H đơn vị W/K. Công suất điện P đi vào bộ gia nhiệt với hiệu suất chuyển thành nhiệt cho bồn η.
Mc dT/dt = mc(T_i−T)+ηP−H(T−T_a).
M dùng kg, m dùng kg/s, c dùng J/(kg·K), P dùng W, t dùng s. H là độ dẫn nhiệt tổn thất, không phải diện tích hay một hằng số thời gian. Phương trình giả định không sôi, không phản ứng và không có biến thiên tồn lượng gây thêm tích trữ.
Muốn giữ T=T_sp ở trạng thái ổn định, lệnh bù từ mô hình là P_ff=[mc(T_sp−T_i)+H(T_sp−T_a)]/η. Đây là lệnh công suất điện; không tự đồng nhất với phần trăm mở van hơi hay phần trăm dòng điều khiển.
Bước thay đổi lưu lượng và nhiệt thiếu nếu không bù
Cho M=200 kg, c=4.000 J/(kg·K), m_0=2 kg/s, T_i=T_a=25°C, T_sp=65°C, H=100 W/K và η=0,95. Nhiệt mang theo dòng là 320 kW, thất thoát 4 kW, nên P_0=341,052632 kW.
Khi lưu lượng tăng lên 3 kg/s và các đại lượng còn lại giữ nguyên, nhiệt cần vào dòng thành 480 kW. Lệnh mới là P_1=509,473684 kW, tăng 168,421053 kW điện, tương ứng 160 kW nhiệt hữu ích. Không nhân toàn bộ P_0 với 1,5 vì phần thất thoát 4 kW không tăng theo lưu lượng trong mô hình này.
Nếu công suất vẫn giữ P_0, nhiệt độ ổn định mới là T_ss=[mcT_i+ηP_0+HT_a]/(mc+H)=51,776860°C. Khi bù đúng và nhiệt cấp thay đổi tức thời, vế phải tại T=65°C bằng không nên nhiệt độ có thể giữ nguyên trong mô hình. Điều kiện tức thời sẽ được bỏ ở phần tiếp theo.
Sai số phép đo truyền trực tiếp vào phần bù
Tại m=3 kg/s, nếu cảm biến nhiệt đầu vào đọc cao hơn thực tế 1 K, thuật toán cho rằng cần ít nhiệt hơn. Nhiệt cấp bị thiếu mc×1=12 kW; sai lệch nhiệt độ ổn định là −12.000/(12.000+100)=−0,991736 K khi chưa có feedback trim.
Nếu đồng hồ lưu lượng đọc thấp 5%, tức 2,85 thay cho 3 kg/s, nhiệt cấp bị thiếu 0,15×4.000×40=24 kW. Sai lệch ổn định là −1,983471 K. Dấu của sai lệch là một kiểm tra quan trọng: đọc lưu lượng thấp trong nhiệm vụ gia nhiệt phải dẫn đến bù thiếu, không phải bù thừa.
Phép tính này cũng cho thấy vì sao không nên cộng cùng một lưu lượng đầy đủ hai lần. Nếu bộ phản hồi đang có bias cung cấp P_0, feedforward có thể dùng phần thay đổi ΔP_ff=P_ff−P_0. Nếu feedforward cung cấp toàn bộ P_ff, phần phản hồi phải được định nghĩa là hiệu chỉnh quanh mức đó. Hai cấu trúc đều có thể dùng, nhưng cần nhất quán khi chuyển chế độ.
Bù đúng trạng thái cuối vẫn có thể tạo một hõm nhiệt
Giữ bước lưu lượng từ 2 lên 3 kg/s, nhưng giả định nhiệt thực cấp vào bồn có trễ bậc nhất τ_a=15 s sau lệnh. Nhiệt hữu ích khi đó là q(t)=484.000−160.000exp(−t/15) W. Lệnh điện đã tăng đúng mức, nhưng nhiệt thực chưa tăng ngay.
Đặt θ=T−65 và τ_p=Mc/(mc+H)=66,115702 s. Với θ(0)=0, phương trình sau bước là Mc dθ/dt+(mc+H)θ=−160.000exp(−t/15). Nghiệm:
θ(t)=−[160.000/(Mc)] × [exp(−t/τ_p)−exp(−t/τ_a)]/(1/τ_a−1/τ_p).
Nhiệt độ thấp nhất ở t=ln(τ_p/τ_a)/(1/τ_a−1/τ_p)=28,779746 s, với θ=−1,941224 K. Tại 60 s, sai lệch vẫn −1,494776 K trước khi trở về không. Đây là đáp ứng của bồn và chấp hành đã nêu, chưa có vòng phản hồi nhiệt độ. Giá trị cuối đúng không chứng minh chất lượng quá độ tốt.
Từ bù tĩnh sang bù động theo hai đường truyền
Trong mô hình tuyến tính quanh một điểm làm việc, gọi y là sai lệch nhiệt độ, u là sai lệch lệnh thao tác và d là nhiễu đo được. Viết y=G_p(s)u+G_d(s)d. Nếu u=F(s)d, triệt tiêu lý tưởng cần F=−G_d/G_p. MathWorks trình bày cấu trúc này cùng việc ghép với phản hồi. [2]
Xét một bộ mô hình độc lập để minh họa động học, không phải mô hình bồn vừa tính: G_p=0,12exp(−8s)/(40s+1), đơn vị K/kW; G_d=−5exp(−20s)/(25s+1), đơn vị K/(kg/s). Khi đó F=41,666667(40s+1)/(25s+1)exp(−12s), đơn vị kW/(kg/s).
Với bước nhiễu d=0,2 kg/s, lệnh bù bắt đầu sau 12 s ở 13,333333 kW rồi giảm theo u(t)=8,333333[1+0,6exp{−(t−12)/25}] kW. Lệnh cuối là 8,333333 kW. Phần nhô đầu là tác dụng của lead–lag, không phải tích phân sai số nhiệt độ.
Nếu chỉ dùng gain tĩnh 41,666667 và vẫn đặt đúng trễ 12 s, sau thời điểm t=20 s nhiệt độ sai lệch là y=exp[−(t−20)/25]−exp[−(t−20)/40] K. Nó đạt −0,171329 K tại t=51,333575 s rồi về không. Mô hình bù động lý tưởng triệt tiêu toàn bộ đáp ứng này, nhưng chỉ trong điều kiện mô hình đúng và chấp hành thực hiện được lệnh.
Điều kiện nhân quả không thể bỏ bằng phép chia hàm truyền
Với hai đường bậc nhất có trễ θ_p và θ_d, F chứa exp[−(θ_d−θ_p)s]. Nếu θ_d nhỏ hơn θ_p, hệ số trở thành một phép đi trước thời gian, đòi hỏi giá trị nhiễu tương lai. NPTEL nêu rõ trường hợp đó không thể thực hiện trực tiếp bằng bộ điều khiển chỉ biết tín hiệu hiện tại. [3]
Trong mô hình minh họa, nếu trễ đường nhiễu giảm từ 20 xuống 5 s trong khi đường thao tác vẫn 8 s, phép bù lý tưởng cần biết trước nhiễu 3 s. Có thể xem xét đo xa hơn phía thượng lưu, giảm trễ chấp hành hoặc chấp nhận bù không hoàn hảo. Không thay exp(+3s) bằng một khối trễ âm rồi coi bài toán triển khai đã được giải.
Giới hạn công suất và vai trò phần phản hồi
Giả sử riêng bộ gia nhiệt bồn bị giới hạn ở 450 kW điện. Sau bước lưu lượng, dù thuật toán yêu cầu 509,473684 kW, mức nhiệt độ ổn định cao nhất trong mô hình là 25+0,95×450.000/12.100=60,330579°C. Không bộ bù nào tạo được công suất chưa có; phần feedback cần thấy giới hạn để tránh tích lũy lệnh vô ích.
Với feedback u=C(r−y)+Fd, đường nhiễu thành y/d=(G_d+G_pF)/(1+G_pC) khi r=0. Feedforward giảm phần nhiễu đã mô tả, feedback xử lý sai mô hình và nhiễu chưa đo trong khả năng vòng. Khi thử nghiệm, cần ghi đồng thời nhiễu thô, nhiễu đã lọc, phần bù, phần phản hồi, tổng lệnh, công suất thực và nhiệt độ. Nhờ đó mới phân biệt được bù sai gain, lệch thời gian, nhiễu cảm biến hay giới hạn thiết bị.
Nguồn tham khảo
[1] The Feed Forward Controller — Doug Cooper and Allen Houtz, Control Guru / Control Station Inc., 9 April 2015. Truy cập ngày 09/11/2026.
[2] Temperature Control in Heat Exchanger — MathWorks, Control System Toolbox Documentation, R2026b. Truy cập ngày 09/11/2026.
[3] Module 5: Other control strategies — Lecture 1: Feedforward control, page 6 — NPTEL, course 103103037. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
