i. Mục đích, đối tƣợng và phạm vi nghiên cứu của đề tài
3.2.1. Giới thiệu phần mềm DX6
Phần mềm DX6 là phần mềm QHTN được sử dụng rộng rãi trong các ngành khoa học kỹ thuật như hoá học, vật liệu, cơ khí v.v. nhằm xây dựng và giải các bài toán thực nghiệm. Phần mềm được được xây dựng dựa trên những lý thuyết cơ bản của xác suất thống kê và lý thuyết QHTN.
Phần mềm là công cụ trực quan giúp người nghiên cứu giảm bớt được số lượng các mẫu thí nghiệm, từ đó giúp giảm thời gian cũng như kinh phí trong quá trình thiết kế và thực hiện thí nghiệm.
DX6 cho phép xác định mô hình toán học trong đó có sự tương tác giữa các biến đến các hàm mục tiêu bằng phương pháp đáp ứng bề mặt trên mô hình không gian ba chiều, đưa ra kết luận mô hình có đạt yêu cầu hay không thông qua tính toán các thông số xác suất thống kê mà các phần mềm QHTN khác như MS.excel hay IRRISTAT không làm được. DX6 cho phép phân tích tối đa 999 hàm mục tiêu cùng với số biến lên đến 12. Đặc biệt phần mềm đưa ra các dự báo xác định được điều kiện tối ưu cho cả bài toán một mục tiêu và đa mục tiêu giúp người thực hiện có thể tìm được đáp án phù hợp với mục tiêu thí nghiệm của mình.
62
3.2.2. Các bƣớc thực hiện cơ bản trên phần mềm DX6
Phần mềm DX6 được thiết kế gần gũi với người sử dụng thông qua thao tác lựa chọn các thanh công cụ được thực hiện bằng chuột. Để bắt đầu chương trình người sử dụng nhấn đúp vào biểu tượng của phần mềm, xuất hiện giao diện chính của phần mềm như trên Hình 3.3.
Nhấn vào File trên thanh công cụ, lựa chọn New Design để bắt đầu chương trình mới hoặc Open Dedesign để mở một chương trình trong thư viện của người sử dụng. Với số lượng yếu tố đầu vào biết trước (number of factors) dóng theo cột tương ứng để chọn phương pháp quy hoạch và từ đó dóng theo hàng ngang để xác định số thử nghiệm cơ bản phải thực hiện (Experiments) như Hình 3.4. Sau khi đã lựa chọn, tiếp theo bấm nút
continue ở góc phải màn hình để tiếp tục.
Các bước được thực hiện tuần tự theo lưu đồ thể hiện ở Hình 3.5, phương pháp sử dụng phần mềm giới thiệu ở phụ lục 4.
Hình 3.3. Màn hình giao diện chính của phần mềm DX6
63
3.2.3. Phân tích kết quả
Sau khi hoàn thành thiết kế thí nghiệm. Để phân tích các kết quả tính toán sử dụng
File, Open Design trên thanh công cụ và chọn
Rsm.dx6 từ thư mục dữ liệu chương trình phần mềm. Bấm vào Open để tải dữ liệu.
Kích chuột vào Block chọn Display Point Type để xem phân tích tại từng điểm thí nghiệm, như Hình 3.6 .
Để phân tích các kết quả tính toán bấm vào mục Analysis và yếu tố đầu ra cần phân tích. Bấm vào ANOVA để kiểm tra sự phù hợp của mô hình đã lựa chọn, như ở Hình 3.7.
Hình 3.7. Kiểm tra sự phù hợp của mô hình theo chuẩn Fisher
Hình 3.6. Lựa chọn hiển thị dạng điểm
Hình 3.5. Các bước thực hiện cơ bản trên phần mềm DX6
Đặt tên và thứ nguyên của các yếu tố đầu vào Khai báo các biến đầu ra
Lựa chọn phương pháp thiết kế thí nghiệm Các điểm thử nghiệm
Lựa chọn các thông số quy hoạch Lựa chọn mô hình toán học
64
Với bài toán thiết lập ở trên từ Hình 3.7 nhận thấy giá trị của hệ số Fisher theo tính toán là Ftn = 25,95 và mô hình toán học như lựa chọn trên Hình 3.7 là phù hợp (Significant).
Dạng phương trình hồi quy phụ thuộc các yếu tố đầu vào được thể hiện trên Hình 3.8.
Để khảo sát giá trị của yếu tố đầu ra phụ thuộc các yếu tố đầu vào bấm vào Model Graphs, kết quả được thể hiện theo dạng đường đồng mức như hoặc dạng không gian như Hình 3.9.
3.2.4. Giải bài toán tối ƣu hoá trên phần mềm DX6
Để giải bài toán tối ưu nhấn chuột vào mục Optimization, kích vào Numerical để lựa chọn mục tiêu tối ưu (giá trị lớn nhất, nhỏ nhất…), khoảng khảo sát của các yếu tố đầu vào. Để biết kết quả tính toán nhấn vào Solutions, kết quả tính toán được thể hiện như trên Hình 3.10.
Hình 3.8. Dạng phương trình hồi quy tìm được
65
Giá trị tối ưu của yếu tố đầu ra được thể hiện là giá trị thực còn giá trị của các yếu tố đầu vào là các giá trị đã được mã hóa, để chuyển từ giá trị mã hóa này sang giá trị thực thì
có thể áp dụng công thức 3.4.
3.3. Thực hiện tối ƣu tham số điều chỉnh động cơ
3.3.1. Ảnh hƣởng của các tham số điều chỉnh tới các tính năng kinh tế kỹ thuật của động cơ sử dụng hệ thống CR
Động cơ diesel nói chung cũng như động cơ diesel sử dụng HTNL CR nói riêng có nhiều tham số điều chỉnh như pf, s, , Tlm, Tbt, Tkn...
3.3.1.1. Ảnh hưởng của nhiệt độ nước làm mát
Hệ thống làm mát có tác dụng tản nhiệt khỏi các chi tiết, giữ cho nhiệt độ của các chi tiết không vượt quá giá trị cho phép và do đó bảo đảm điều kiện làm việc bình thường của động cơ.
Nhiệt độ nước làm mát được xác định thông qua cân bằng nhiệt của động cơ. Việc nâng cao nhiệt độ trong hệ thống làm mát sẽ giúp giảm lượng nhiệt truyền cho nước làm mát do đó góp phần nâng cao hiệu suất có ích của động cơ. Tuy nhiên việc nâng cao nhiệt độ trong hệ thống làm mát chỉ tới một giới hạn nào đó đảm bảo sức bền, độ cứng vững và tuổi thọ các chi tiết và không để xẩy ra kích nổ trong động cơ xăng. Theo nhiều nghiên cứu cho thấy cứ tăng nhiệt độ nước làm mát lên 100C sẽ làm giảm nhiệt lượng truyền cho nước làm mát 4÷5% tuy nhiên khi tăng nhiệt độ nước làm mát cũng sẽ làm tăng nhiệt lượng do khí thải mang đi [3].
Việc tăng nhiệt độ nước làm mát cũng làm tăng nhiệt độ thành xilanh và nhiệt độ trung bình của màng dầu giữa piston và xilanh. Tăng nhiệt độ thành xilanh sẽ làm giảm hệ số nạp, theo [64] nếu tăng nhiệt độ nước làm mát từ 60 lên 900C thì hệ số nạp sẽ giảm 2÷3%. Ngoài ra, khi đó tăng nhiệt độ lớp dầu giữa piston và xilanh tăng sẽ làm giảm công ma sát do giảm độ nhớt của dầu. Cứ tăng nhiệt độ nước làm mát 100C thì công ma sát sẽ giảm 3,5÷10% [64].
66
Do tác động qua lại của những yếu tố kể trên dẫn tới để động cơ làm việc tốt, phát huy được các tính năng kinh tế kỹ thuật thì mỗi động cơ sẽ có vùng nhiệt độ làm việc của nước làm mát tối ưu. Đối với các động cơ làm mát bằng nước nhiệt độ này khoảng 75÷850C [3].
3.3.1.2. Ảnh hưởng của nhiệt độ dầu bôi trơn và nhiệt độ khí nạp
Khi nhiệt độ dầu bôi trơn thấp, độ nhớt cao dẫn đến lực cản ma sát lớn làm tăng công tổn hao cơ giới trong động cơ. Tuy nhiên, nhiệt độ dầu bôi trơn quá cao sẽ làm giảm độ nhớt của dầu dẫn tới có thể phá vỡ màng dầu bôi trơn làm tăng mài mòn các chi tiết.
Nhiệt độ khí nạp cũng là một thông số ảnh hưởng đến các chỉ tiêu kinh tế kỹ thuật của động cơ. Khi nhiệt độ khí nạp tăng sẽ làm cho mật độ không khí bị giảm. Hình 3.11 cho thấy ảnh hưởng của nhiệt độ khí nạp lớn tới suất tiêu thụ nhiên liệu của động cơ…
Hình 3.11. Ảnh hưởng của nhiệt độ khí nạp lớn tới suất tiêu thụ nhiên liệu của động cơ [64]
a. Tại tốc độ 1500v/ph và mô mem khác nhau b. Tại mô mem (50 N.m) và tốc độ khác nhau
a)
67
Hình 3.12. Ảnh hưởng của các tham số điều chỉnh tới các tính năng kinh tế kỹ thuật của động cơ [65]
3.3.1.3. Ảnh hưởng của hệ số dư lượng không khí
Trên các xe đời mới áp dụng công nghệ phun nhiên liệu CR, để đáp ứng được các qui định ngày càng nghiêm ngặt về lượng khí độc hại phát ra, người ta bố trí thêm cảm biến ô xy trên đường xả để xác định nồng độ ô xy trong khí xả, từ đó quyết định đến lượng phun liệu. Như vậy hệ số dư lượng không khí ảnh hưởng đến mức tiêu hao nhiên liệu, sự phát thải của động cơ.
3.3.1.4. Ảnh hưởng của áp suất phun và góc phun sớm
Hình 3.12 cho thấy ảnh hưởng của một số tham số điều chỉnh cơ bản tới tính năng kinh tế kỹ thuật của động cơ trong đó việc tăng áp suất phun có thể cải thiện công suất, giảm tiêu thụ nhiên liệu và phát thải độc hại [65]. Trong khi tối ưu góc phun sớm không chỉ cải thiện các tính năng này mà nó còn có khả năng giúp giảm ồn cho động cơ. Ngoài ra bằng cách thay đổi quy luật cung cấp nhiên liệu như tăng số lần phun trong một chu trình có thể giúp giảm phát thải độc hại và ồn.
Cùng với sự phát triển của công nghệ điện – điện tử và công nghệ chế tạo, áp suất phun cũng như số lần phun trong một chu trình ngày một tăng, áp suất phun cao sẽ làm tăng chiều dài tia phun và tốc độ nhiên liệu giúp cho nhiên liệu phun tơi hơn như thể hiện trên Hình 3.13 [65] và được phân bố khắp không gian buồng cháy đảm bảo nhiên liệu bay hơi nhanh và hòa trộn đồng đều do đó sẽ giúp quá trình cháy được cải thiện. Tuy nhiên áp suất cao cũng làm tăng tốc độ cháy của hỗn hợp, làm tăng rung động và ồn của động cơ. Ngoài ra để tăng áp suất phun thì công dẫn động bơm cao áp cũng tăng nên sẽ làm tăng tổn hao cơ giới của động cơ.
Qua phân tích và khảo sát như trên thì hai tham số điều chỉnh φs,pf là hai tham số cơ bản quan trọng nhất của động cơ bởi nó quyết định tới cả công suất, tiêu hao nhiên liệu, sự phát thải và độ ồn của động cơ. Còn các thông số khác có ảnh hưởng nhưng không nhiều. Vì vậy chọn hai tham số φs,pf được điều chỉnh để khảo sát được xem như là tham số tiêu biểu. Hơn nữa, nếu khảo sát tất các tham số, sẽ phải mất rất nhiều thời gian cũng như số
68
lượng mẫu thí nghiệm tăng lên nhiều lần nằm ngoài khả năng của nghiên cứu sinh. Do đó trong quá trình thí nghiệm các tham số như , Tlm, Tbt, Tkn được giữ ổn định ở chế độ hợp lý đối với từng loại động cơ (cụ thể theo số liệu khuyến cáo của nhà sản xuất).
3.3.2. Quy trình tối ƣu các tham số điều chỉnh
Động cơ sử dụng hệ thống nhiên liệu CR có nhiều tham số điều chỉnh như áp suất phun, góc phun sớm, số lần phun trong một chu trình, hệ số dư lượng không khí λ... Việc điều chỉnh chính xác các tham số này sẽ nâng cao tính năng kinh tế và kỹ thuật của động cơ.
Công việc tối ưu các tham số điều chỉnh được thực hiện trên các băng thử hiện đại ở các chế độ làm việc và các hàm mục tiêu khác nhau thông qua thực hiện nhiều ma trận thử nghiệm, trong đó mỗi ma trận thử nghiệm tương ứng với một chế độ làm việc của động cơ. Việc xác định bộ tham số (φs, pf) tối ưu được tiến hành bằng thực nghiệm. Ứng với mỗi chế độ làm việc, giá trị (φs, pf) tối ưu được lựa chọn theo mục tiêu tính kinh tế hoặc tính hiệu quả cao nhất. Do hai thông số φs và pf ảnh hưởng đồng thời đến đặc tính làm việc của động cơ [65]. Cho nên việc xác định hai thông số này cần phải được tiến hành đồng thời.
Theo đặc tính điều chỉnh, mô men có ích (Me) theo pf và φs luôn có xu hướng đạt cực trị ở một giá trị nhất định [2]. Xét về bài toán QHTN, khi sử dụng QHTG cấp I thì giảm lần thử nghiệm mà vẫn tìm ra hàm hồi quy (quan hệ biến vào và ra). Nếu không đạt được tiêu chuẩn của bài toán QHTN thì phải thực hiện tiếp bài toán QHTG cấp II và như vậy số thử nghiệm sẽ tăng lên. Nhưng nếu so sánh hàm hồi quy thì QHTG cấp II nó thể hiện quan hệ chính xác hơn. Ngoài ra do giới hạn nghiên cứu của luận án, với hai thông số đầu vào nên việc thực hiện số điểm thử nghiệm theo lý thuyết QHTN là không nhiều, cùng với theo nhiều tài liệu của nghiên cứu khác thì có thể bỏ qua bài toán QHTG cấp I. Do đó nghiên cứu này đã lựa chọn phương pháp QHTG cấp II mà bỏ qua tuần tự đi từ phương pháp QHTG cấp I.
69
3.3.3. Tiến hành bài toán quy hoạch trực giao cấp II
Cácbước thực hiện bài toán QHTN được thể hiện theo lưu đồ thuật toán được thể hiện trên hình 3.14, với 2 thông số đầu vào là thời điểm phun và áp suất phun và một thông số đầu ra là mô men.
Hình 3.14. Lưu đồ thuật toán các bước thực hiện bài toán QHTN TG cấp II
Phương trình hồi quy sẽ có dạng như sau:
y = b0 + b1x1 + b2x2 + b12 x1x2 + b11x12 + b22x22 (3.33) Ma trận thử nghiệm được xây dựng theo phương án thực nghiệm bậc 2 của Box – Wilson. Số thử nghiệm được tính theo công thức.
N = 2k + 2k +N0 (3.34)
Trong đó:
Bắt đầu
- Xác định miền biến thiên, tâm qui hoạch - Mã hóa
Chọn dạng phương trình hồi quy
Thực hiện N thí nghiệm
Tính toán các hệ số bj
Giả bài toán tối ưu cực trị
Kết thúc Đúng Loại bỏ hệ số bj Sai KT có nghĩa của các hệ số bj ti >0,05(N0 - 1) KT có nghĩa của PTHQ F > Ftt Đúng Sai
70
- 2k là số thí nghiệm của QHTG cấp I với các Zj = Zjmin hoặc Zj = Zjmax.
- 2k là số thí nghiệm tại các điểm “sao”: xj = +TG hoặc xj = -TG. - N0 là thí nghiệm tại tâm Zj = 0
j
Z
Sau khi tính toán các hệ số bj, sự có nghĩa của các hệ số hồi quy bj được kiểm tra bằng cách áp dụng tiêu chuẩn Student [32], thực hiện các thí nghiệm tại tâm quy hoạch hoặc sử dụng các thí nghiệm song song. Mục đích của bước này là loại bỏ các hệ số bj nhỏ ít ảnh hưởng tới yếu tố đầu ra nhằm đơn giản hóa phương trình hồi quy. Sau khi kiểm tra sự có nghĩa của các hệ số bj, bước tiếp theo tiến hành kiểm tra sự có nghĩa của phương trình hồi quytheo tiêu chuẩn Fisher (F), các bước kiểm tra được trình bày trong [29].
Để giải bài toán tối ưu tìm điểm cực trịcó thể lựa chọn phương pháp giải bài toán tối ưu theo phương pháp leo dốc hoặc luân phiên từng biến giải bài toán phỏng định [29]. Trong đề tài luận án này kết quả giải tối ưu được thực hiện theo phương pháp leo dốc thông qua phần mềm DX6.
3.4. Kết luận chƣơng 3
Nêu lên các phương pháp QHTN trong việc nghiên cứu thực nghiệm. Phương pháp QHTN trực giao trong nghiên cứu thực nghiệm nhiều yếu tố cho độ tin cậy cao khi xây dựng các mô hình toán học thực nghiệm, cho phép phân tích ảnh hưởng đồng thời của nhiều yếu tố đầu vào tới một hay nhiều yếu tố đầu ra của một quá trình công nghệ. Từ đó chọn ra phương pháp QHTN trực giao cấp II xây dựng bộ dữ liệu điều chỉnh của động cơ.
Trong thực tế, các bài toán quy hoạch thực nghiệm thường có kích thước lớn, do đó khó có thể tính toán bằng tay. Vì vậy cần chương trình hoá thuật toán bằng các ngôn ngữ lập trình hoặc bằng các phần mềm QHTN chuyên dụng. Trong chương đã giới thiệu áp dụng phần mềm QHTN DX6 cho quá trình thiết kế thí nghiệm cũng như việc giải bài toán tối ưu tìm bộ dữ liệu chuẩn.
Động cơ diesel sử dụng HTNL CR sử dụng trên phương tiện cơ giới có nhiều tham số điều chỉnh. Cùng với sự phân tích như ở trên, luận án giới hạn hai tham số điều chỉnh φs,pf