BÀI TẬP VÀ BÀI TẬP LỚN CƠ SỞ ĐIỀU KHIỂN QUÁ TRÌNH

Phần III: Hiệu chỉnh bộ PID Mỗi sinh viên tự làm riêng Cho hệ thống có sơ đồ như sau: + Hàm truyền của bộ PID: 1 i T K Trong đó: KPID là hệ số khuyếch đại của bộ điều khiển Ti là thờ

PHABN 1

BÀI TẬP VÀ BÀI TẬP LỚN

Phần I: Khảo sát các đặc tính của khâu động học cơ bản

(mỗi sinh viên tự là riêng)

1) Khâu quán tính bậc nhất

Cho bảng số liệu thí nghiệm thu được từ

1 )

(











Ts

K s G x K y dt

dy

T

Cho các tham số K=20; T=50 và K=20; T=100

Nhiệm vụ:

- Viết chương trình trên Matlab

- Khảo sát các đặc tính trong miền thời gian: hàm quá độ h(t), hàm trọng lượng g(t)

- Xác định các tham số K và T trên đồ thị

2) Khâu bậc hai

1 2 )

( )

( ) (

2

2

2













Ts s

T

K s

G t

Kx t y dt

dy T

dt

y

d

T





Cho các tham số K=20, T=10,zeta =0, 0.25, 0.5, 0.75, 1

Nhiệm vụ:

- Viết chương trình trên Matlab và Simulink

- Khảo sát các đặc tính trong miền thời gian h(t),w(t)

- Nhận xét ảnh hưởng của độ suy giảm zeta đến đặc tính quá độ của khâu bậc hai

Hướng dẫn dùng Matlab:

Khởi động MATLAB thấy hiện ra dấu nhắc >> ta sẽ đánh các câu lệnh sau:

>>num=[b0 b1 bm-1 bm]; % Định nghĩa tử số nếu hệ số nào không có % gán bằng 0

>>den=[a0 a1 an-1 an]; % Định nghĩa mẫu số nếu hệ số nào không có % gán bằng 0

>>w=tf(num,den) % Định nghĩa hàm truyền đạt w

>>step(w) % Vẽ hàm quá độ h(t)

>>impulse(w) % Vẽ hàm quá độ xung k(t)

>> hold on : giữ hình vẽ hiện tại trong cửa sổ Figure Lệnh này hữu ích khi ta cần vẽ nhiều biểu đồ trong cùng một cửa sổ Figure Sau khi vẽ xong biểu đồ thứ nhất, ta gõ lệnh hold on để giữ lại hình vẽ sau đó vẽ tiếp các biểu đồ khác Các biểu đồ lúc sau sẽ vẽ đè lên biểu đồ thứ nhất trong cùng một cửa sổ Figure này Nếu không muốn giữ hình nữa, ta

gõ lệnh hold off

>>grid on : kẻ lưới trên cửa sổ Figure Nếu không muốn kẻ lưới , ta gõ lệnh grid off

PHABN 2

Phần II: Khảo sát đặc tính của hệ thống

(mỗi sinh viên tự làm riêng)

1) Tìm hàm truyền tương đương của hệ thống

Bằng cách sử dụng các lệnh cơ bản conv, tf, series, parallel, feedback, tìm biểu

thức hàm truyền tương đương G(s) của hệ thống sau:

Hướng dẫn: Bước đầu tiên nhập hàm truyền cho các khối G1, G2, dùng lệnh tf Sau

đó tuỳ theo cấu trúc mắc nối tiếp, song song hay phản hồi mà ta gõ lệnh series ( hoặc dấu

*) , parallel ( hay dấu +), feedback tương ứng để thực hiện việc kết nối các khối với nhau Trong báo cáo trình bày rõ trình tự thực hiện các lệnh này

Ví dụ:

>> G1=tf([1 1],conv([1 3],[1 5])) % Nhập hàm truyền G1

>> G3=tf(1 ,[1 0]) % Nhập hàm truyền G3

>>G13=G1+G3 % Tính hàm truyền tương đương của G1, G3

% hoặc G13=parallel(G1,G3)

>>G21=feedback(G2,H1) % Tính hàm truyền tương đương của G2, H1

Tiếp tục tính tương tự cho các khối còn lại

2) Tìm hàm truyền tương đương của hệ thống

Cho hệ thống kín có cấu trúc như hình vẽ:

Với K= 8; K=17.564411; K=20

- Viết chương trình xác định hàm truyền đạt của hệ thống khi thay đổi K trong ba

trường hợp đã cho?

- Khảo sát các đặc tính trong miền thời gian h(t),w(t) của hệ với mỗi K tương ứng?

PHABN 3

- Nhận xét so sánh 3 trường hợp trên?

Phần III: Hiệu chỉnh bộ PID

(Mỗi sinh viên tự làm riêng)

Cho hệ thống có sơ đồ như sau:

+ Hàm truyền của bộ PID:

1

i

T K

Trong đó: KPID là hệ số khuyếch đại của bộ điều khiển

Ti là thời gian tích phân

Td là thời gian vi phân

+ Hàm truyền của đối tượng

Các thông số ban đầu của bộ điều khiển PID và đối tượng điều khiển được cho trong bảng sau:

Số cuối cùng

Nhiệm vụ:

PHABN 4

a Viết chương trình mô phỏng cho hệ thống sử dụng Simulink

b Vẽ quá trình quá độ với các thông số ban đầu, tính độ quá điều chỉnh, thời gian quá độ, sai lệch tĩnh

c Hiệu chỉnh bộ PID để nâng cao chất lượng của hệ thống

Phần IV: Ứng dụng SIMULINK để tổng hợp hệ thống điều khiểntự động

(làm theo nhóm – Bài tập lớn)

1 MỤC ĐÍCH :

SIMULINK là một công cụ rất mạnh của Matlab để xây dựng các mô hình một cách trực quan và dễ hiểu Để mô tả hay xây dựng hệ thống ta chỉ cần liên kết các khối có sẵn trong thư viện của SIMULINK lại với nhau Sau đó, tiến hành mô phỏng hệ thống để xem xét ảnh hưởng của bộ điều khiển đến đáp ứng quá độ của hệ thống và đánh giá chất lượng hệ thống

2 CHUẨN BỊ :

Để thực hiện các yêu cầu trong bài thí nghiệm này, sinh viên cần phải chuẩn bị kỹ và hiểu

rõ các khối cơ bản cần thiết trong thư viện của SIMULINK Sau khi khởi động Matlab, ta

gõ lệnh simulink hoặc nhấn vào nút simulink trên thanh công cụ thì cửa sổ SIMULINK hiện ra: 2 thư viện chính áp dụng trong bài thí nghiệm này

2.1 Các khối được sử dụng trong bài thí nghiệm:

a Các khối nguồn – tín hiệu vào (source):

PHABN 5

Khối Step (ở thư viện Simulink \ Sources) có chức năng xuất ra tín hiệu hàm bậc thang Double click vào khối này để cài đặt các thông số:

• Step time : khoảng thời gian đầu ra chuyển sang mức Final value

kể từ lúc bắt đầu mô phỏng Cài đặt giá trị này bằng 0

• Initial value : Giá trị ban đầu Cài đặt bằng 0

• Final value : Giá trị lúc sau Cài đặt theo giá trị ta muốn tác động tới hệ thống Nếu là hàm bậc thang đơn vị thì giá trị này bằng 1

• Sample time : thời gian lấy mẫu Cài đặt bằng 0

Khối Signal Generator (ở thư viện Simulink \ Sources) là bộ phát tín hiệu xuất ra các tín hiệu hìng sin, hìng vuông, hìng răng cưa và ngẫu nhiên (cài đặt các dạng hìng này trong mục Wave form)

b Các khối tải – thiết bị khảo sát ngõ ra (sink):

Khối Mux (ở thư viện Simulink \ Signals Routing) là bộ ghép kênh nhiều ngõ vào 1 ngõ ra, từ ngõ ra này ta đưa vào Scope để xem nhiều tín hiệu trên cùng một cửa sổ Double click vào khối này để thay đổi số kênh đầu vào (trong mục Number of inputs)

Khối Scope (ở thư viện Simulink \ Sinks) là cửa sổ xem các tín hiệu theo thời gian, tỉ lệ xích của các trục được điều chỉnh tự động để quan sát tín hiệu một cách đầy đủ

Khối XY Graph dùng để xem tương quan 2 tín hiệu trong hệ thống (quan sát mặt phẳng pha)

c.Các khối xử lý – khối động học :

Khối Sum (ở thư viện Simulink \ Math Operations) là bộ tổng (cộng hay trừ) các tín hiệu, thường dùng để lấy hiệu số của tín hiệu đặt với tín hiệu phản hồi Double click để thay đổi dấu của bộ tổng

Khối Gain (ở thư viện Simulink \ Math Operations) là bộ tỉ lệ Tín hiệu sau khi qua khối này sẽ được nhân với giá trị Gain Double click để thay đổi giá trị độ lợi Gain

Khối Transfer Fcn (ở thư viện Simulink \ Continuous) là hàm truyền của hệ tuyến tính Double click để thay đổi bậc và các hệ số của hàm truyền Cài đặt các thông số:

- Numerator : các hệ số của đa thức tử số

- Denominator : các hệ số của đa thức mẫu số

PHABN 6

Khối Relay (ở thư viện Simulink \ Discontinuities) là bộ điều khiển rơle

2 vị trí có trễ (còn gọi là bộ điều khiển ON-OFF) Các thông số : _ Switch on point : nếu tín hiệu đầu vào lớn hơn giá trị này thì ngõ ra của khối Relay lên mức ‘on’

- Switch off point : nếu tín hiệu đầu vào nhỏ hơn giá trị này thì ngõ

ra của khối Relay xuống mức ‘off’

- Output when on : giá trị của ngõ ra khi ở mức ‘on’

- Output when off : giá trị của ngõ ra khi ở mức ‘off’

Nếu tín hiệu đầu vào nằm trong khoảng (Switch on point, Switch off point) thì giá trị ngõ ra giữ nguyên không đổi

Khối PID controller (ở thư viện Simulink Extras \ Additional

Linear) là bộ điều khiển PID với hàm truyền ( )

KP : hệ số tỉ lệ (proportional term) KI: hệ số tích phân (integral term) KD: hệ số vi phân (derivative term)

Khối Saturation (ở thư viện Simulink \ Discontinuities) là một khâu bão hòa Các thông số cài đặt:

- Upper limit : giới hạn trên Nếu giá trị đầu vào lớn hơn Upper limit thì ngõ ra luôn bằng giá trị Upper limit

- Lower limit : giới hạn dưới Nếu giá trị đầu vào nhỏ hơn Lower limit thì ngõ ra luôn bằng giá trị Lower limit

Khâu bão hoà dùng để thể hiện giới hạn biên độ của các tín hiệu trong thực tế như : áp ra cực đại của bộ điều khiển đặt vào đối tượng, áp nguồn…

2.2 Các bước tiến hành để xây dựng một ứng dụng mới trong SIMULINK:

- Sau khi khởi động Matlab, gõ lệnh simulink hoặc nhấn vào nút simulink trên thanh công cụ thì cửa sổ SIMULINK hiện ra (cần chở 1 chút)

- Trong cửa sổ SIMULINK, vào menu File / New để mở cửa sổ cho một ứng dụng mới Kích chuột vào các thư viện đã giới thiệu ở mục II.1 để chọn khối cần tìm Kích chuột trái vào khối này, sau đó kéo và thả vào cửa sổ ứng dụng vừa mới tạo ra Double click vào khối này để cài đặt và thay đổi các thông số

- Có thể nhân số lượng các khối bằng cách dùng chức năng Copy và Paste Kích chuột trái nối các ngõ vào / ra của các khối để hình thành sơ đồ hệ thống

- Có thể dời một hoặc nhiều khối từ vị trí này đến vị trí khác bằng cách nhấp chuột để chọn các khối đó và kéo đến vị trí mới Dùng phím Delete để xóa các phần không cần thiết hay bị sai khi chọn

- Có thể viết chú thích trong cửa sổ ứng dụng bằng cách double click vào một vị trí trống

và gõ câu chú thích vào Vào menu Format / Font để thay đổi kiểu chữ

PHABN 7

- Như vậy, mô hình hệ thống đã xây dựng xong Bây giờ tiến hành mô phỏng hệ thống bằng cách vào menu Simulation / Simulation Parameters để cài đặt các thông số mô phỏng

Cửa sổ Simulation Parameters hiện ra như sau:

- Start time : thời điểm bắt đầu mô phỏng Mặc định chọn bằng 0

- Stop time : thời điểm kết thúc mô phỏng Giá trị này chọn theo đặc tính của hệ thống Nếu hệ thống có thời hằng lớn thì giá trị Stop time cũng phải lớn để quan sát hết thời gian quá độ của hệ thống

- Các thông số còn lại chọn mặc định như ở hình kế bên

- Chạy mô phỏng bằng cách vào menu Simulation / Start Khi thời gian mô phỏng bằng giá trị Stop time thì quá trình mô phỏng dừng lại Trong quá trình mô phỏng, nếu ta muốn dừng nửa chừng thì vào menu Simulation / Stop

3 Thiết kế hệ thống điều khiển

Một caloriphe được sử dụng để đốt nóng không khí trước khi đưa vào buồng sấy bằng

cách sử dụng chất tải nhiệt là hơi nước từ nồi hơi công nghiệp có lưu lượng là Gs, áp suất

không khí ngoài trời có nhiệt độ t kk1 , độ ẩm tương đối φ 1 , lưu lượng G kk1 Không khí sau

khi ra khỏi caloriphe t kk2 cần đảm bảo nhiệt độ được ổn định ở giá trị chủ đạo t sp

Yêu cầu:

1) Vẽ sơ đồ quy trình công nghệ? xác định mục đích điều khiển từ đó xác định các biến quá trình?

2) Các sách lược điều khiển có thể áp dụng? Hãy lựa chọn một sách lược hợp lý nhất nhằm đạt được yêu câu đặt ra? Biểu diễn sơ đồ khối thể hiện nguyên lý điều khiển và lưu

đồ điều khiển P&ID tương ứng theo sách lược đã chọn?

3) Thiết kế hệ thống điều khiển phản hồi với bộ điều khiển ON//OFF, P, PI và PID để điều khiển quá trình này So sánh và đánh giá chất lượng điều khiển ứng với mỗi bộ điều khiển trên?

Bảng số liệu thí nghiệm dưới đây được ghi nhận khi theo dõi diễn biến của nhiệt độ theo thời gian tính từ lúc khởi động hệ thống (mở nhanh van cấp dòng hơi vào thiết bị từ trạng thái đóng hoàn toàn sang trạng thái mở hoàn toàn) đến khi hệ đạt ổn định

PHABN 8

NHÓM: N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 N11 N12 N13 N14 N15 N16 N17 N18 N19 N20 N21 N22 N23 N24 N25

Thời gian,

phút

0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0

0.5 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0

1 30.0 30.0 30.2 30.2 30.2 30.0 30.0 30.0 30.0 30.0 30.5 30.5 30.5 30.5 30.5 30.5 30.5 30.5 30.5 30.2 30.2 30.2 30.2 30.2 30.2

1.5 30.5 30.5 30.5 30.5 30.5 30.1 30.1 30.1 30.1 30.1 30.5 30.5 30.5 30.5 30.5 30.5 30.5 30.5 30.5 30.2 30.2 30.2 30.2 30.2 30.2

2 30.7 30.7 30.7 30.7 30.7 30.1 30.1 30.1 30.1 30.1 31.0 31.0 31.0 31.0 31.0 31.0 31.0 31.0 31.0 30.5 30.5 30.5 30.5 30.5 30.5

2.5 31.0 31.0 31.0 31.0 31.0 30.2 30.2 31.0 31.0 31.0 33.0 33.0 34.0 33.0 33.0 32.0 32.0 31.0 31.0 30.5 30.5 30.5 30.5 30.5 30.5

3 31.5 31.5 31.0 31.0 31.5 30.3 30.5 32.0 32.0 33.0 36.0 36.0 38.0 35.0 35.0 33.0 33.0 32.0 32.0 31.0 31.0 31.0 31.0 31.0 31.0

3.5 32.0 32.0 32.0 32.0 32.0 30.6 31.0 34.6 36.1 37.7 39.2 40.7 42.3 38.9 52.7 36.7 34.0 34.0 34.0 31.0 31.0 31.0 32.5 32.5 32.5

4 32.5 33.0 32.0 33.0 33.0 31.0 31.5 38.5 41.3 44.2 47.0 49.8 52.7 46.8 54.8 42.6 40.5 38.4 36.3 31.5 31.5 31.5 33.0 33.0 34.0

4.5 34.0 35.0 33.5 35.0 35.5 31.5 32.5 41.8 45.7 49.7 53.6 57.5 61.5 53.8 56.1 47.8 44.9 41.9 38.9 32.9 33.8 34.8 35.7 36.7 37.6

5 36.7 38.0 35.3 39.3 40.6 32.7 34.0 44.6 49.5 54.3 59.2 64.1 68.9 60.0 57.0 52.5 48.8 45.0 41.3 35.4 37.3 39.1 40.9 42.7 44.5

6 42.4 44.9 39.9 47.4 49.9 35.0 37.5 47.0 52.6 58.3 63.9 69.6 75.2 65.6 57.6 56.7 52.2 47.8 43.3 37.8 40.4 43.0 45.6 48.1 50.7

7 47.4 50.9 43.9 54.4 57.9 37.0 40.5 49.0 55.3 61.6 67.9 74.2 80.6 70.5 58.0 60.4 55.3 50.3 45.2 39.9 43.2 46.5 49.8 53.1 56.4

8 51.8 56.1 47.4 60.5 64.8 38.7 43.1 50.7 57.5 64.4 71.3 78.2 85.1 74.9 58.2 63.7 58.1 52.4 46.8 41.8 45.7 49.7 53.6 57.5 61.5

9 55.5 60.6 50.4 65.7 70.8 40.2 45.3 52.1 59.5 66.8 74.2 81.5 88.9 78.8 58.5 66.6 60.5 54.4 48.3 43.5 48.0 52.6 57.1 61.6 66.1

10 58.8 64.5 53.0 70.3 76.1 41.5 47.3 53.3 61.1 68.8 76.6 84.4 92.1 82.3 58.6 69.2 62.7 56.1 49.6 45.1 50.1 55.2 60.2 65.2 70.3

11 61.6 67.9 55.3 74.2 80.6 42.6 49.0 54.3 62.4 70.6 78.7 86.8 94.9 85.3 58.7 71.5 64.6 57.7 50.7 46.5 52.0 57.5 63.0 68.5 74.1

12 64.1 70.9 57.2 77.7 84.5 43.6 50.4 55.2 63.6 72.0 80.4 88.8 97.2 88.1 58.8 73.6 66.3 59.0 51.8 47.8 53.7 59.7 65.6 71.5 77.5

13 66.2 73.4 58.9 80.6 87.9 44.5 51.7 55.9 64.6 73.2 81.9 90.5 99.2 90.5 58.9 75.4 67.8 60.3 52.7 49.0 55.3 61.6 67.9 74.2 80.6

14 68.0 75.6 60.4 83.2 90.8 45.2 52.8 56.6 65.4 74.3 83.1 92.0 100.8 92.7 59.0 77.0 69.2 61.3 53.5 50.0 56.7 63.4 70.0 76.7 83.4

15 69.6 77.5 61.7 85.5 93.4 45.8 53.8 57.1 66.1 75.2 84.2 93.2 102.2 94.6 59.0 78.4 70.4 62.3 54.2 51.0 58.0 64.9 71.9 78.9 85.9

16 71.0 79.2 62.8 87.4 95.6 46.4 54.6 57.5 66.7 75.9 85.1 94.3 103.4 96.3 59.0 79.7 71.4 63.2 54.9 51.8 59.1 66.4 73.6 80.9 88.2

17 72.2 80.6 63.8 89.1 97.5 46.9 55.3 57.9 67.2 76.5 85.8 95.1 104.4 97.8 59.1 80.9 72.4 63.9 55.4 52.6 60.1 67.7 75.2 82.7 90.3

18 73.2 81.9 64.6 90.5 99.2 47.3 55.9 58.2 67.6 77.1 86.5 95.9 105.3 99.2 59.1 81.9 73.2 64.6 55.9 53.3 61.1 68.8 76.6 84.4 92.1

19 74.1 83.0 65.3 91.8 100.6 47.7 56.5 58.5 68.0 77.5 87.0 96.5 106.0 100.4 59.1 82.8 74.0 65.2 56.4 53.9 61.9 69.9 77.9 85.9 93.8

20 74.9 83.9 65.9 92.9 101.9 48.0 56.9 58.7 68.3 77.9 87.5 97.0 106.6 101.4 59.1 83.6 74.7 65.7 56.8 54.5 62.7 70.9 79.0 87.2 95.4

21 75.6 84.7 66.5 93.8 102.9 48.2 57.4 58.9 68.6 78.2 87.9 97.5 107.1 102.4 59.2 84.3 75.2 66.2 57.1 55.0 63.4 71.7 80.1 88.4 96.8

22 76.2 85.4 66.9 94.7 103.9 48.5 57.7 59.1 68.8 78.5 88.2 97.9 107.6 103.2 59.2 84.9 75.8 66.6 57.5 55.5 64.0 72.5 81.0 89.5 98.0

23 76.7 86.0 67.3 95.4 104.7 48.7 58.0 59.2 69.0 78.7 88.5 98.2 108.0 104.0 59.2 85.5 76.2 67.0 57.7 55.9 64.6 73.2 81.9 90.5 99.2

24 77.1 86.6 67.7 96.0 105.4 48.9 58.3 59.4 69.1 78.9 88.7 98.5 108.3 104.7 59.2 86.0 76.7 67.3 58.0 56.3 65.1 73.9 82.7 91.4 100.2

25 77.5 87.0 68.0 96.5 106.0 49.0 58.5 59.5 69.3 79.1 88.9 98.7 108.5 105.2 59.2 86.4 77.0 67.6 58.2 56.7 65.6 74.5 83.4 92.2 101.1

26 77.8 87.4 68.3 97.0 106.5 49.1 58.7 59.5 69.4 79.2 89.1 98.9 108.8 105.8 59.2 86.8 77.4 67.9 58.4 57.0 66.0 75.0 84.0 93.0 102.0

27 78.1 87.8 68.5 97.4 107.0 49.3 58.9 59.6 69.5 79.3 89.2 99.1 109.0 106.2 59.2 87.2 77.7 68.1 58.6 57.3 66.4 75.5 84.6 93.6 102.7

28 78.4 88.1 68.7 97.7 107.4 49.4 59.0 59.7 69.6 79.4 89.3 99.2 109.1 106.7 59.2 87.5 77.9 68.3 58.7 57.5 66.7 75.9 85.1 94.3 103.4

29 78.6 88.3 68.9 98.0 107.8 49.4 59.2 59.7 69.6 79.5 89.4 99.3 109.2 107.0 59.2 87.8 78.1 68.5 58.9 57.8 67.0 76.3 85.5 94.8 104.1

30 78.8 88.5 69.0 98.3 108.1 49.5 59.3 59.8 69.7 79.6 89.5 99.4 109.4 107.4 59.2 88.0 78.4 68.7 59.0 58.0 67.3 76.6 86.0 95.3 104.6

31 78.9 88.7 69.2 98.5 108.3 49.6 59.4 59.8 69.7 79.7 89.6 99.5 109.5 107.7 59.2 88.2 78.5 68.8 59.1 58.2 67.6 77.0 86.4 95.7 105.1

32 79.1 88.9 69.3 98.7 108.5 49.6 59.5 59.8 69.8 79.7 89.7 99.6 109.5 107.9 59.2 88.4 78.7 69.0 59.2 58.3 67.8 77.2 86.7 96.1 105.6

33 79.2 89.0 69.4 98.9 108.7 49.7 59.5 59.9 69.8 79.8 89.7 99.7 109.6 108.1 59.2 88.6 78.8 69.1 59.3 58.5 68.0 77.5 87.0 96.5 106.0

34 79.3 89.2 69.4 99.0 108.9 49.7 59.6 59.9 69.8 79.8 89.8 99.7 109.7 108.4 59.2 88.8 79.0 69.2 59.4 58.6 68.2 77.7 87.3 96.8 106.4

35 79.4 89.3 69.5 99.2 109.0 49.8 59.6 59.9 69.9 79.8 89.8 99.8 109.7 108.5 59.2 88.9 79.1 69.3 59.5 58.8 68.4 78.0 87.6 97.1 106.7

36 79.5 89.4 69.6 99.3 109.2 49.8 59.7 59.9 69.9 79.9 89.8 99.8 109.8 108.7 59.2 89.0 79.2 69.3 59.5 58.9 68.5 78.2 87.8 97.4 107.0

37 79.6 89.5 69.6 99.4 109.3 49.8 59.7 59.9 69.9 79.9 89.9 99.8 109.8 108.8 59.2 89.1 79.3 69.4 59.6 59.0 68.7 78.3 88.0 97.7 107.3

38 79.6 89.5 69.7 99.5 109.4 49.8 59.8 59.9 69.9 79.9 89.9 99.9 109.8 109.0 59.2 89.2 79.4 69.5 59.6 59.1 68.8 78.5 88.2 97.9 107.6

39 79.7 89.6 69.7 99.5 109.5 49.9 59.8 59.9 69.9 79.9 89.9 99.9 109.9 109.1 59.2 89.3 79.4 69.5 59.7 59.2 68.9 78.6 88.4 98.1 107.8

40 79.7 89.6 69.8 99.6 109.5 49.9 59.8 60.0 69.9 79.9 89.9 99.9 109.9 109.2 59.2 89.4 79.5 69.6 59.7 59.3 69.0 78.8 88.5 98.3 108.0

41 79.7 89.7 69.8 99.6 109.6 49.9 59.8 60.0 69.9 79.9 89.9 99.9 109.9 109.3 59.2 89.5 79.5 69.6 59.7 59.3 69.1 78.9 88.7 98.4 108.2

42 79.8 89.7 69.8 99.7 109.6 49.9 59.9 60.0 70.0 79.9 89.9 99.9 109.9 109.4 59.2 89.5 79.6 69.7 59.8 59.4 69.2 79.0 88.8 98.6 108.4

43 79.8 89.8 69.8 99.7 109.7 49.9 59.9 60.0 70.0 80.0 89.9 99.9 109.9 109.4 59.2 89.6 79.6 69.7 59.8 59.5 69.3 79.1 88.9 98.7 108.5

44 79.8 89.8 69.9 99.8 109.7 49.9 59.9 60.0 70.0 80.0 90.0 99.9 109.9 109.5 59.2 89.6 79.7 69.7 59.8 59.5 69.3 79.2 89.0 98.8 108.7

Nhiệt độ, oC