iv Các nội dung chính trong đề tài
4.5. Tiến hành thử nghiệm và kết quả
4.5.1. Xây dựng đường đặc tính ngồi
Đặc tính ngồi của động cơ được lựa chọn khảo sát là đặc tính ngồi sử dụng, được
79
định nghĩa là đường đặc tính tiến sát tới đường đặc tính khĩi đen và được xác định thơng qua việc đo độ khĩi bằng thiết bị Smokemeter [67].
Trong luận án để xác định điểm thử nghiệm nằm trên đường đặc tính ngồi thì tại điểm tâm quy hoạch ( 0 0
,
s pf
) tiến hành điều chỉnh lượng nhiên liệu phun Gnl sao cho độ khĩi đo được đạt tới giá trị đặt trước FSN = 9,5 [67]. Tại mỗi ma trận thử nghiệm giá trị Gnl được giữ cố định.
Trên cơ sở các thơng số kỹ thuật của động cơ, chọn dải tốc độ khảo sát từ 1000 đến 3000 v/ph, chia dải tốc độ này thành 10 phần đều nhau. Quá trình xây dựng đặc tính được xuất phát từ điểm 1 ứng với tốc độ là 3000 v/ph hướng đến điểm số 2 ứng với tốc độ là 1000 v/ph như thể hiện trên Hình 4.10 và
bộ giá trị (φs, pf) từ các cơng
trình nghiên cứu đã thực hiện trên băng thử. Giá trị φs và pf tối ưu tại các tốc độ khác nhau
được xác định theo lưu đồ thuật tốn thể hiện trên Hình 4.11.
Hình 4.11. Thuật tốn xác định bộ số liệu φs và pf tại đặc tính ngồi
Hình 4.10. Giới hạn vùng làm việc khi khảo sát đường đặc tính ngồi
Bắt đầu
k:=0
Chọn (0s, p0f) tại điểm nmax
Chọn dải (s, pf) khảo sát tại điểm nmax - k.n Thực hiện bài tốn quy hoạch thực
nghiệm
Xác định được (s, pf) tối ưu tại điểm nmax- k.n với hàm mục tiêu là Mmax
k = 10 Kết thúc
Đúng k:= k+1
(s, pf) tối ưu tại điểm nmax- (k-1).n được chọn là (0s, p0f) tại điểm nmax- k.n
80
4.5.1.1. Xây dựng thử nghiệm tại tốc độ 3000 v/ph
Quy trình thử nghiệm tại tốc độ 3000 v/ph được thiết kế theo phương pháp QHTG cấp
II. Giá trị (φs, pf) tối ưu được xác định bằng cách tiến hành tuần tự theo các bước như mục
4.4.1. Thơng qua các thơng số kỹ thuật của động cơ, băng thử và các kết quả nghiên cứu tại Phịng thí nghiệm Động cơ đốt trong. Miền khảo sát áp suất phun và gĩc phun sớm được lựa chọn dựa vào đường đặc tính khĩi đen của động cơ cũng như qua thực hiện khảo nghiệm tại phịng thí nghiệm để tìm miền làm việc chủ yếu và ổn định của động cơ. Đồng thời kết hợp các loại tài liệu và đề tài đã nghiên cứu trước đĩ để xác định miền khảo sát như sau: 0 0 s 400 bar p 800 bar 8 32 f
Như vậy tâm quy hoạch ( 0 0 ,
s pf
) = (200; 600 bar). Để thống nhất ký hiệu và thuận tiện trong quá trình thiết lập thử nghiệm, tiến hành mã hĩa với x1 là φs, x2 là pf.
Với số lượng yếu tố đầu vào là 2, áp dụng cơng thức (3.30) số thử nghiệm cần thực hiện N = 10 thể hiện như bảng 4.3 với 4 thử nghiệm ở mức trên và mức dưới, 4 thử nghiệm ở mức “sao” và 2 thử nghiệm ở mức trung tâm [18, 29, 31]. Trong bảng 4.3 cũng như các bảng của các điểm thử nghiệm tại các tốc độ ở các vị trí xét cĩ số liệu của hai biến mã hĩa X1 và X2 là ± 1,41, theo nguyên tắc QHTN ngồi việc xét ở vùng ± 1 cịn phải xét điểm lân cận để tăng thêm độ chính xác của bộ số liệu, điểm lân cận chính bằng (±) căn bậc hai của số biến đầu vào.
Bảng 4.3. Các điểm thử nghiệm tại tốc độ 3000 v/ph
TT x1 x2 X1 X2 Me (N.m) 1 8 400 -1 -1 25,5 2 32 400 1 -1 26,6 3 8 800 -1 1 27,5 4 32 800 1 1 26,4 5 3 600 -1,41 0 23,6 6 37 600 1,41 0 23,2 7 20 320 0 -1,41 28 8 20 880 0 1,41 27,6 9 20 600 0 0 28,5 10 20 600 0 0 28,5
81
Việc tính tốn tiếp theo cĩ thể tính tốn bằng tay như đã được trình bày trong mục 3.1.3.2. Tuy nhiên để rút ngắn thời gian tính tốn trong luận án tác giả đã sử dụng phần mềm quy hoạch thực nghiệm DX6. Kết quả tính tốn các giá trị bj bằng phần mềm được thể hiện như trên bảng 4.4.
Bảng 4.4. Các giá trị bj tại tốc độ 3000 v/ph
b0 b1 b2 b11 b22 b12
28,95 -0,071 0,154 -2,55 -0,35 -0,55
Theo yêu cầu của bài tốn QHTN bước tiếp theo cần kiểm tra sự cĩ nghĩa của các hệ số bj theo tiêu chuẩn Student nhằm đơn giản hĩa phương trình hồi quy bằng cách loại bỏ các hệ số bj cĩ mức độ ý nghĩa nhỏ. Tuy nhiên việc này cĩ thể dẫn tới bỏ đi các hệ số cĩ ảnh hưởng lớn tới việc giải bài tốn tối ưu (các hệ số b11, b22), do đĩ theo hàm mục tiêu tại
đặc tính ngồi là xác định giá trị của (φs, pf) mà tại đĩ mơ men đạt giá trị lớn nhất nên bước
kiểm tra sự cĩ nghĩa của các hệ số bj sẽ được bỏ qua.
Như vậy hàm số mơ men phụ thuộc áp suất phun và gĩc phun sớm tại tốc độ n = 3000 v/ph được thể hiện như cơng thức sau:
Me = 28,95 – 0,071*X1 + 0,154*X2 – 2,55*X12 -0,35*X22 – 0,55*X1*X2 (4.1) Kết quả tính tốn hệ số Fisher trên phần mềm cho thấy Ftn=25,95 trong khi đĩ tra bảng tiêu chuẩn Fisher (phụ lục 3) với mức ý nghĩa p= 0,05, bậc tự do f1 = N-L = 4, f2 = n0-1 = 1 (N là số điểm thực nghiệm, L là số các hệ số bj cĩ nghĩa, n0 là số thử nghiệm tại tâm quy hoạch) F = 224,6. So sánh nhận thấy Ftn < F do đĩ mơ hình tốn đã chọn phù hợp với thực nghiệm.
Từ phương trình 4.1 giá trị Me phụ thuộc vào pf và φs được khảo sát trên phần mềm DX6, kết quả được thể hiện như trên Hình 4.12 và Hình 4.13.
Kết quả khảo sát cho thấy giá trị mơmen lớn nhất là Memax = 28,97 Nm tại điểm (X1, X2) = (-0,04; 0,25) quy đổi từ giá trị mã hĩa sang giá trị thực (x1, x2) = (19,520; 650 bar).
82
4.5.1.2. Xây dựng thử nghiệm tại điểm nmax – Δn
Các bước xác định bộ thơng số tối ưu (φs, pf) tại điểm 2800 v/ph được tiến hành tương
tự như điểm 3000 v/ph. Tuy nhiên cĩ một số điểm cần lưu ý như sau:
- Giá trị Gnl được xác định để đạt được độ khĩi theo yêu cầu, giá trị (φs, pf) tối ưu xác định ở điểm số 3000 v/ph được chọn là tâm quy hoạch ( 0 0
,
s pf
) tại điểm 2800 v/ph.
- Do cùng nằm trên đặc tính ngồi và bước tốc độ (Δn) khơng quá lớn nên dải áp suất phun và gĩc phun sớm tại điểm 2800 v/ph được lựa chọn hẹp hơn, chọn dải áp suất và gĩc phun sớm khảo sát như sau:
0 0 s 500 bar p 800 bar 12 27, 04 f
Bước tiếp theo, tương tự như tại tốc độ 3000 v/ph, tiến hành xây dựng ma trận thí nghiệm với 10 thử nghiệm như bảng 4.5. Các giá trị bj tại điểm 2800 v/ph sau khi tính tốn được thể hiện trong bảng 4.6.
Như vậy giá trị của mơ men phụ thuộc vào áp suất và gĩc phun sớm được thể hiện thơng qua cơng thức sau:
Me = 29,46 – 0,072*X1 – 0,071*X2 – 0,55*X12 -0,22*X22 – 0,05*X1*X2. (4.2) Theo kết quả tính tốn từ phần mềm Ftn = 2,69 <F = 224,6. Như vậy mơ hình tốn đã chọn phù hợp với thực nghiệm.
83
Bảng 4.5. Các điểm thử nghiệm tại tốc độ 2800 v/ph
TT x1 x2 X1 X2 Me (N.m) 1 12 500 -1 -1 28,9 2 27,04 500 1 -1 28,5 3 12 800 -1 1 29 4 27,04 800 1 1 28,4 5 8,92 650 -1,41 0 28,2 6 30,12 650 1,41 0 28,5 7 19,52 438,5 0 -1,41 29,2 8 19,52 861,5 0 1,41 28,8 9 19,52 650 0 0 29,45 10 19,52 650 0 0 29,45 Bảng 4.6. Các giá trị bj tại tốc độ 2800 v/ph b0 b1 b2 b11 b22 b12 29,46 -0,072 -0,071 -0,55 -0,22 -0,05
Kết quả khảo sát được thể hiện trên Hình 4.14 và Hình 4.15. Giá trị mơ men lớn nhất là 29.47 tại điểm (X1, X2) = (-0,06; - 0,15) quy đổi từ giá trị mã hĩa sang giá trị thực (x1, x2) = (19,070; 627,5 bar).
84
4.5.1.3. Xây dựng thử nghiệm tại điểm nmax – k.Δn
Với biến k chạy từ 2 đến 10, việc thiết kế thí nghiệm được tiến hành tương tự như đối
với điểm 3000 v/ph và 2800 v/ph, với các giá trị (φs, pf) tối ưu tìm được tại điểm nmax – (k-
1).Δn sẽ được chọn làm tâm quy hoạch cho điểm nmax – k.Δn.
Sau khi tính tốn trên phần mềm DX6, kết quả thể hiện mối quan hệ giữa mơ men với gĩc phun sớm, áp suất phun cũng như bộ thơng số gĩc phun sớm và áp suất phun tối ưu theo tốc độ động cơ ở đường đặc tính ngồi được tổng hợp theo bảng 4.7, bảng 4.8.
Bảng 4.7. Giá trị bj ở các tốc độ khác nhau trên đường đặc tính ngồi
n (v/ph) b0 b1 b2 b11 b22 b12 1000 27,32 -0,93 -0,25 -1 -0,3 0,038 1200 28,36 -0,89 -0,067 -0,97 -0,47 0,28 1400 29,34 -0,47 -0,37 -0,62 -0,9 -0,05 1600 30,38 -0,084 -0,066 -0,41 -0,14 -0,05 1800 31,19 -0,049 -0,06 -0,47 -0,3 -0,2 2000 31,61 -0,035 -0,003 -0,5 -0,33 -0,2 2200 31,33 -0,043 -0,016 -0,48 -0,088 0 2400 30,29 -0,051 -0,0017 -0,45 -0,062 -0,012 2600 30,18 -0,07 -0,024 -0,42 -0,16 -0,075 2800 29,46 -0,072 -0,071 -0,55 -0,22 -0,05 3000 28,95 -0,071 0,15 -2,55 -0,35 -0,55
85
Bảng 4.8. Giá trị mơ men lớn nhất ở đường đặc tính ngồi
n (v/ph) φs (độ) pf (bar) Me (N.m) 1000 10,44 542,29 27,59 1200 12,61 575,2 28,45 1400 15,14 588,2 29,46 1600 16,98 598,5 30,38 1800 17,47 608,1 31,19 2000 17,64 615,7 31,61 2200 17,87 614,8 31,33 2400 18,11 621,8 30,59 2600 18,5 622,6 30,18 2800 19,07 627,5 29,58 3000 19,52 650 28,97
Từ Hình 4.16 cho thấy giá trị mơ men lớn nhất trên đường đặc tính ngồi đạt được tại
tốc độ 2000 v/ph. Giá trị (φs, pf) tối ưu cĩ xu hướng tăng khi tốc độ động cơ tăng. Trong đĩ giá trị của φs thay đổi trong khoảng từ 100 đến 200, giá trị của pf tương ứng thay đổi trong
khoảng từ 540 bar đến 650 bar.
86