1.5. .Lựa chọn mơ hình tính toán
2.4. Mơ hình nghiên cứu ảnh hưởng của phụ tải nhiệt đến sự tương tác
Từ mơ hình tương tác và mơ hình tính tốn trường nhiệt độ, biến dạng nhiệt, ta rút ra được các hệ phương trình sử dụng để nghiên cứu ảnh hưởng của phụ tải nhiệt như sau:
ÕT = «v 2T = 1 Õ2T + -- ^T +Õ2T Õ 2T A Õ 2 Õt c p ^ Õx2 Õy2 ự h } = AT [ [ 4 0 0 0 ] T A = D - d A - D - d' = [ d + D . A D - T0 ) ] - [ [ + [.Ap [T - T0)] (2.61)
Hệ phương trình (2.61) biểu diễn sự cân bằng nhiệt trong một vật thể, sử dụng để tính tốn trường nhiệt độ của pít tơng và xi lanh (phương trình thứ nhất); biểu diễn biến dạng nhiệt của vật thể, sử dụng để tính tốn biến dạng nhiệt sau khi xác định được trường nhiệt độ của cặp pít tơng - xi lanh (phương trình thứ hai); xác định khe hở nhiệt giữa pít tơng và xi lanh ở trạng thái nguội (phương trình thứ ba); xác định khe hở nhiệt giữa pít tơng và xi lanh tại mỗi chế độ phụ tải nhiệt khác nhau (phương trình thứ tư). Sau khi tính được trường nhiệt độ và trường biến dạng nhiệt, từ phương trình này sẽ xác định được khe hở nhiệt phụ thuộc vào phụ tải và hành trình của pít tơng. Giá trị này được sử dụng làm giá trị đầu vào để tính ảnh hưởng của phụ tải nhiệt đến chuyển động phụ của pít tơng.
t* = 3 r = N ( I OI + ỉcmiH ) A.l .m, I + lcm1H Ú N x = t2P xX m x 6m, 2m, - y t ỉ pỗ xO mx 4 L + y Ụ ị N c t 2M mc 5 m c _|_ t ỉ l Pờ X m x r 6 I 2I,O, 2 I r Pk.- A.t (2.62) >
Hệ phương trình (2.62) sử dụng để tính chuyển động phụ của pít tơng: phương trình thứ nhất dùng để xác định thời gian chuyển động ngang của pít tơng trong khe hở phụ thuộc vào thông số khe hở nhiệt A; phương trình thứ hai để tính tốn góc xoay của pít tơng trong khe hở phụ thuộc vào thông số khe hở nhiệt A; phương trình thứ ba là phương trình chuyển vị ngang của pít tơng, phương trình này cũng phụ thuộc vào thơng số khe hở nhiệt A, từ phương trình này sẽ tính được vận tốc chuyển động ngang của pít tơng trong khe hở phụ thuộc vào khe hở nhiệt A tại các góc quay trục khuỷu tương ứng. Giá trị vận tốc ngang này là đầu vào để tính tốn lực va đập của pít tơng lên xi lanh.
V 1 N £ = kN'3 = —— v0 s i n o pt H—-—r (1 - c o s o pt ) - o p mj co2 1 ị ---- --- J N vd ( t ) sin0p (t - t1 )dt1 - (2.63) m1o p 0 - ẳ — T X ‘- Ị N rj( t ) e ~- ( '- s i n o ’j t ~ t1)d t1 i- PLFL< Ị X ( x ) d x 0 0
Phương trình (2.63) sử dụng để tính lực va đập của pít tơng lên xi lanh. Ảnh hưởng của phụ tải nhiệt lên lực va đập được thể hiện qua vận tốc ban đầu vo (vận tốc chuyển động ngang của pít tơng tương ứng tại thời điểm xảy ra va đập).
Mơ hình tổng thể trình bày ở hệ phương trình (2.61), (2.62) và (2.63) có tính tổng qt hơn các mơ hình đã xét như là các mơ hình thành phần, vì đây là mơ hình có khe hở và có tương tác, ngồi ra cịn xét thêm ảnh hưởng của phụ tải nhiệt thông qua trị số khe hở nhiệt tại các chế độ công tác của động cơ.
2.5. Lựa chọn phương pháp tính và phần mềm tính tốn
2.5.1. Các phương pháp tính tốn trường nhiệt độ và lựa chọn phương pháp tính pháp tính
2.5.1.1. Phân loại bài tốn truyền nhiệt
Việc phân loại các phương pháp tính tốn trường nhiệt độ có liên hệ mật thiết với dạng biểu diễn toán học của bài tốn truyền nhiệt. Ngồi ra theo dấu hiệu chung có thể phân chia chúng thành các phương pháp phân tích chính xác và gần đúng, phương pháp số.
Căn cứ mục tiêu cuối cùng, bài toán truyền nhiệt được chia theo các dạng: - Bài toán trực tiếp: Cho phương trình vi phân cơ bản, hoặc cho hệ phương trình vi phân với điều kiện biên xác định trường nhiệt độ.
- Bài toán nghịch: Cho biết trường nhiệt độ và phương trình vi phân, xác định điều kiện hoặc phân bố nhiệt ban đầu.
- Bài toán thụ động: Cho trước trường nhiệt độ và điều kiện biên cần xác định các hệ số của phương trình vi phân.
- Bài toán cảm ứng: Cho biết sự phân bố nhiệt độ, điều kiện biên và lí tính của vật thể, phải chuẩn hóa hoặc thiết lập phương trình vi phân của quá trình.
Nếu phân loại bài toán truyền nhiệt từ mối quan hệ giữa nhiệt độ và thời gian thì chúng ta chia làm hai loại:
- Bài toán tĩnh (ổn định).
- Bài tốn động (khơng ổn định).
Khi mà hệ số của phương trình vi phân và của điều kiện biên khơng phụ thuộc vào nhiệt độ thì bài tốn gọi là tuyến tính, trường hợp ngược lại là phi tuyến.
2.5.1.2. Các phương pháp tính tốn trường nhiệt độ và lựa chọn phương pháp tính
Tính trường nhiệt độ của các kết cấu có thể dùng các phương pháp sau: các phương pháp phân tích chính xác và phân tích gần đúng, các phương pháp số, các phương pháp đồ thị, các phương pháp tương tự. Các phương pháp trên có thể được áp dụng độc lập hoặc hợp thành các tổ hợp khác nhau.
Việc sử dụng các phương pháp phân tích giới hạn bởi dạng miền và điều kiện biên, bởi thế phương pháp phân chia biến chỉ phù hợp cho vật thể có hình dáng đơn giản. Phương pháp thế nhiệt chỉ có thể sử dụng trong trường hợp điều kiện đồng nhất.
Phương pháp giải bài toán phi tuyến được nghiên cứu ít hơn vì nó rất phức tạp, phổ biến nhất vẫn là bước chuyển từ bài toán phi tuyến sang bài tốn tuyến tính. Khó khăn nhất vẫn là việc giải bài tốn phi tuyến, khơng ổn định với sự bố trí 2, 3 số đo cho vật thể nhiều thành phần. Hiện nay, phương pháp duy nhất có thể giải bài toán trên là phương pháp số, đương nhiên phương pháp số có thể giải được các bài tốn tuyến tính.
Trong các phương pháp số để giải bài tốn biên thì phổ biến nhất trong kỹ thuật là phương pháp sai phân hữu hạn. Những năm gần đây thì phương pháp phần tử hữu hạn (PTHH) rất hay được sử dụng trong các ngành kỹ thuật vì có nhiều ưu điểm. Phương pháp PTHH kế thừa tư tưởng của phương pháp
sai phân hữu hạn và phương pháp xấp xỉ hàm, với sự hỗ trợ của phương pháp này có thể giải được các bài toán nhiệt, kết cấu, điện từ, mô phỏng, tối ưu hóa... Phương pháp PTHH là cơ sở của lĩnh vực mơ phỏng hóa trong các bài tốn thiết kế, kiểm nghiệm. Thơng qua sự phát triển của khoa học kỹ thuật đồ họa trên máy tính người ta có thể mơ phỏng các hoạt động của kết cấu. Điều này cho phép giảm chi phí và thời gian thực hiện các thí nghiệm theo phương pháp truyền thống.
Phương pháp mô phỏng toán học (phương pháp tương tự) có cơ sở là mơ phỏng tương tự các q trình vật lý khác nhau theo bản chất: lan tỏa nhiệt, sự phân bố thế trong mạch điện, dòng chất lỏng, sự khuếch tán và các hiện tượng khác bằng các phương trình vi phân.
Các phương pháp đồ thị sử dụng để giải các bài toán tuyến tính với cấu trúc đơn giản.
Trường hợp bài tốn động (không ổn định) với chi tiết của thiết bị năng lượng thì phần lớn các trường hợp đó là phi tuyến. Điều đó có liên quan đến sự thay đổi tính chất nhiệt - vật lí của vật liệu, cấu trúc vật liệu vào nhiệt độ. Càng đặc biệt quan trọng với các chi tiết như xi lanh, pít tơng, nắp máy của động cơ đốt trong, rơ to tuabin khí, nơi mà nhiệt độ tăng quá 500oC. Trong những trường hợp mà điều kiện biên phụ thuộc vào nhiệt độ bề mặt thì được xếp vào bài tốn phi tuyến loại II.
Việc giải bài toán bù truyền nhiệt ở dạng trực tiếp hoặc dạng ẩn có thể sử dụng phương pháp đồng dạng. Lí thuyết đồng dạng cho phép biến đổi các đại lượng vật lí thơng thường thành biến tổ hợp, qua đó làm giảm số lượng biến. Ngoài ra lí thuyết đồng dạng cịn cho phép nhận các phương trình đồng dạng chứa sự phụ thuộc suy rộng, trên cơ sở nghiên cứu thực nghiệm các hiện tượng vật lí và dùng phương pháp phân tích hoặc dùng phương pháp số để giải các mơ hình tốn học của chúng cũng như xác định giới hạn sử dụng các mối liên hệ đó.
Trong bài tốn truyền nhiệt, lí thuyết đồng dạng được sử dụng để nhận các lời giải ở dạng đại lượng không thứ nguyên. Sử dụng lí thuyết đồng dạng trong phương pháp số giải các bài tốn biên khi khơng thể tìm thấy sự phụ thuộc giữa các biến ở dạng trực tiếp sẽ tiết kiệm được rất nhiều thời gian. Trong trường hợp này bằng việc thiết lập mối liên hệ giữa các biến có thể làm giảm các phương án tính tốn.
Qua các phương pháp tính tốn trường nhiệt độ của các chi tiết động cơ đốt trong, ta thấy phương pháp PTHH là phương pháp có thể giải được bài tốn phi tuyến khơng ổn định, cho vật thể nhiều thành phần, có kết cấu phức tạp [8], [17], [19]. "
Trường nhiệt độ của chi tiết xác định bằng lí thuyết PTHH được thực hiện theo các nội dung sau:
Định luật nhiệt động I [17]:
r dT , „ , , , ế „ ế ,
Ệẳ (2.64)
Pc
õt■ + {v}T {L}T^ + {L}T {q} = q
Định luật Fourier vê thông lượng nhiệt [17]:
ịq}= -[D]ị L}T trong đó: [ D] = Kxx 0 0 0 Kyy 0 0 0 Kzz - Ma trận các hệ số dẫn nhiệt. Trường nhiệt độ được xác định bởi [17]:
T = {n}t {Te}
trong đó:
T = T(x, y, z, t) - Trường nhiệt độ cần xác định. {N} = {N(x, y, z)} - Ma trận hàm dạng.
{Te} = {Te(x, y, z, t)} - Véc tơ nhiệt độ tại các nút của phần tử. Hệ phương trình xác định trường nhiệt độ [17]:
[ c i l í í M t ^ M C l + M H r . H a í + í s í + í Q ? } Gradient nhiệt được xác định như sau [17]:
ÕT ÕT ÕT õ x Õy Õx
Véc tơ mật độ dịng nhiệt được tính theo gradient nhiệt [17]:
ị a } = { L } T = ịq] = - [ D]ịa} = - [ D][ B]ịT } (2.65) (2.66) (2.67) (2.68) (2.69) Mật độ dòng nhiệt qua bê mặt đối lưu [17]:
ịq Ỵ = hf (T - Tb ) (2.70) Việc dùng phương pháp PTHH để giải bài tốn tính tốn trường nhiệt độ và biến dạng nhiệt của cặp xi lanh và pít tơng là phù hợp vì xi lanh và pít
tơng phải chịu phụ tải nhiệt lớn trong quá trình làm việc. Bên cạnh đó, với sự phát triển của máy tính và các phần mêm chuyên dụng, việc sử dụng phương pháp PTHH để giải các bài tốn truyền nhiệt có nhiêu thuận lợi. Các phần mêm chuyên dụng có thể sử dụng để giải bài toán truyền nhiệt của động cơ đốt trong bằng phương pháp PTHH như phần mêm ANSYS, FEMAB, SAP v.v... Thơng qua kết quả tính được, người ta có thể mơ phỏng các hoạt động của kết cấu. Điêu này cho phép giảm chi phí và thời gian thực hiện các thí nghiệm theo phương pháp truyên thống. Phần mêm ANSYS do hãng ANSYS của Mỹ phát triển. Đây là một bộ phần mêm hoàn chỉnh dựa trên phân tích PTHH để xác định các ứng xử của hệ vật lí khi chịu tác động của các loại tải trọng khác nhau. Phần mêm ANSYS được sử dụng rộng rãi trên thế giới để giải quyết các bài tốn thiết kế, mơ phỏng tối ưu kết cấu và các quá trình truyên nhiệt, dòng chảy, tĩnh điện, điện từ v . v . cũng như tương tác giữa các mơi trường và hệ vật lí khác nhau. Trong luận án này, tác giả sử dụng phương pháp PTHH cùng phần mêm ANSYS Workbench để tính tốn trường nhiệt độ và biến dạng nhiệt của cặp pít tơng, xi lanh.
2.5.2. Các phương pháp tính tốn tương tác của cặp p ít tơng - xi lanh
Để tính tốn tương tác của cặp pít tơng - xi lanh, trên cơ sở lí thuyết đã xây dựng, có thể dùng phương pháp giải tích trun thống hoặc phương pháp số. Việc sử dụng giải tích để tính tốn sẽ cho cái nhìn tổng quát, rõ ràng, tuy nhiên khi đó việc tính tốn sẽ gặp nhiêu khó khăn, phức tạp nhất là khi khảo sát nhiêu trường hợp, nhiêu yếu tố ảnh hưởng. Sử dụng phương pháp số với sự hỗ trợ của máy tính và các phần mêm chuyên dụng là hướng đi phù hợp hiện nay, trong đó có bài tốn tương tác giữa các vật thể.
Trong bài toán tương tác của cặp pít tơng - xi lanh cần giải quyết hai vấn đê sau: Tìm vận tốc chuyển động phụ của pít tơng khi chuyển động ngang trong khe hở giữa pít tơng và xi lanh và góc quay của pít tơng trong khe hở đó; tìm lực va đập của pít tơng với xi lanh.
Để giải quyết vấn đê thứ nhất, do vận tốc chuyển động phụ thuộc vào lực ngang N nên có thể sử dụng phương pháp ma trận và phần mêm Matlab để giải. Trong luận án này sử dụng phần mêm Matlab để xác định chuyển động phụ.
Vấn đê thứ hai, do sử dụng mơ hình va đập để xác định lực tương tác nên có thể sử dụng các phần mêm tính tốn va đập để giải. Có thể sử dụng phần mêm như LS-DYNA, AUTODYN v . v .
Trong luận án này, sau khi nghiên cứu phương pháp giải bài toán va chạm bằng phần tử hữu hạn [19], [46] đã quyết định lựa chọn phần mêm ANSYS Workbench LS-Dyna để giải quyết bài tốn va chạm giữa pít tơng và xi lanh. Kết luận chương 2
Trên cơ sở mơ hình tính tốn sự tương tác của cặp pít tơng - xi lanh của của Nikishin V.N., lí thuyết biến dạng đàn hồi và lí thuyết truyên nhiệt và biến dạng nhiệt của vật thể đã xây dựng được mơ hình lí thuyết để xác định sự ảnh hưởng của phụ tải nhiệt đến sự tương tác của cặp pít tơng - xi lanh động cơ diesel lai máy phát tàu thủy, cụ thể như sau:
- Khe hở nhiệt là đại lượng phản ánh phụ tải nhiệt của cặp pít tơng - xi lanh động cơ đốt trong.
- Sự tương tác của cặp pít tơng - xi lanh được thể hiện qua hai yếu tố: Thứ nhất là chuyển động phụ của pít tơng trong khe hở giữa pít tơng và xi lanh; Thứ hai là lực va đập của pít tơng lên xi lanh trong khe hở nêu trên.
- Qua kết quả tính tốn trường nhiệt độ và biến dạng nhiệt của pít tơng và xi lanh tại các chế độ công tác của động cơ sẽ xác định được khe hở nhiệt giữa chúng. Đây là cơ sở để tính tốn xác định ảnh hưởng của phụ tải nhiệt đến sự tương tác của cặp pít tơng - xi lanh theo mơ hình đã chọn.
Việc xây dựng mơ hình tương tác của cặp chi tiết pít tơng - xi lanh có xét đến ảnh hưởng của phụ tải nhiệt (thông qua khe hở nhiệt) là đóng góp mới của luận án.
Từ mơ hình lí thuyết đã xây dựng lựa chọn các phần mêm phù hợp để tính tốn là kết hợp phần mêm Matlab để tính tốn chuyển động phụ của pít tơng, phần mêm ANSYS Workbench LS-Dyna để tính tốn va chạm của pít tơng với thành xi lanh; phần mêm ANSYS Workbench Steady - State Thermal và ANSYS Workbench Static Structural để tính tốn nhiệt, biến dạng nhiệt và khe hở nhiệt của cặp pít tơng - xi lanh.
Chương 3. KHẢO SÁT ẢNH HƯỞNG CỦA PHỤ TẢI NHIỆT TỚI SỰ TƯƠNG TÁC CỦA CẶP PÍT TƠNG - XI LANH
ĐỘNG CƠ 6H 12/14
3.1. Đặc điểm kết cấu của cặp pít tơng - xi lanh động cơ 6H 12/14
3.1.1. Giới thiệu sơ bộ về động cơ 6 H 12/14
Động cơ 6 ^ 12/14 là động cơ diesel 4 kì, 6 xi lanh do Liên Xô (cũ) chế tạo để lai máy phát điện và được sử dụng trên các tàu quân sự kiểu 205, 206 của Hải quân nhân dân Việt Nam.
Động cơ có loại có tăng áp, có loại khơng có tăng áp. Cơng suất động cơ từ (59 - 66) kW cho loại không tăng áp, từ (84,5 - 132) kW cho loại có tăng áp.
Thơng số chính của động cơ có tăng áp [5]
- Số xi lanh : 6 - Đường kính xi lanh : 120 mm - Hành trình pít tơng: 140 mm - Thứ tự nổ của xi lanh : 1 - 5 - 3 - 6 - 2 - 4 - Tỉ số nén: 14 - Tỉ số tăng áp: 1,4 - Góc mở sớm xupap nạp (độ GQTK): 45±8 - Góc đóng muộn xupap nạp (độ GQTK): 45±8 - Góc mở sớm xupap thải (độ GQTK): 45±8 - Góc đóng muộn xupap thải (độ GQTK): 45±8