BÀI TẬP LỚN TÍNH TOÁN ĐỘNG CƠ ĐỐT TRONG Đề tài Động cơ 1.5L DOHC VTEC Turbo trên xe Honda CRV 2020

Nội dung thuyết minh gồm: - Nhiệm vụ bài tập lớn được GV hướng dẫn thông qua - Phần số liệu tính toán nhiệt, động học và động lực học có kèm theo các đồ thị - Kết luận của bài tập lớn...

BÀI TẬP LỚN TÍNH TOÁN ĐỘNG CƠ ĐỐT TRONG

Honda CRV 2020

GVHD:

PGS.T

S LÝ VĨNH ĐẠT

TP.HCM Tháng 10 năm 2022

SVTH : Trần Nguyễn Nhật Trường 20145738 Nguyễn Thanh Trường 20145737 Đào Phú Quý 20145371

TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH

KHOA CƠ KHÍ ĐỘNG LỰC

BÀI TẬP LỚN TÍNH TOÁN ĐỘNG CƠ ĐỐT TRONG

Honda CRV 2020

GVHD:

PGS.T

S LÝ VĨNH ĐẠT

SINH VIÊN THỰC HIỆN MSSV ĐÁNH GIÁ HOÀN THÀNH

TRƯỜNG ĐẠI HỌC SPKT TP.HCM CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM

KHOA CƠ KHÍ ĐỘNG LỰC Độc Lập- Tự Do- Hạnh Phúc

BỘ MÔN ĐỘNG CƠ

NHIỆM VỤ BÀI TẬP LỚN MÔN TÍNH TOÁN ĐỘNG CƠ ĐỐT TRONG

1 Số liệu ban đầu

2 Nội dung thuyết minh

2.1 Tính toán nhiệt và xây dựng giản đồ công chỉ thị động cơ;

2.2 Tính toán động lực học cơ cấu piston – trục khuỷu – thanh truyền.

3 Nội dung bản vẽ

3.1 Bản vẽ đồ thị công chỉ thị P-V;

3.2 Bản vẽ đồ thị P – φ, PJ, P1;

3.3 Bản vẽ đồ thị quãng đường Sp, vận tốc Vp, gia tốc Jp của piston.

Ngày giao nhiệm vụ:

Ngày hoàn thành nhiệm vụ:

GIÁO VIÊN HƯỚNG DẪN CHÍNH

(Ký ghi họ tên)

PGS.TS Lý Vĩnh Đạt

LỚP 201453A NHÓM

5

THỜI GIAN THỰC HIỆN

 Tiến độ thực hiện:

1 Tính toán nhiệt: 3 tuần

- Chọn các thông số đầu vào: 1 tuần

- Tính toán nhiệt: 1 tuần

- Vẽ đồ thị công chỉ thị: 1 tuần

2 Tính toán động học, động lực học và vẽ bản vẽ về động học và động lực học: 3

tuần

- Tính toán động học và động lực học: 1 tuần

- Vẽ bản vẽ về động học và động lực học: 2 tuần

3 Viết thuyết minh và chuẩn bị báo cáo: 2 tuần

 Nội dung qui định:

1 Bản vẽ đồ thị công chỉ thị, động lực học,

2 Thuyết minh viết trên khổ giấy A4 (in và file) Nội dung thuyết minh gồm:

- Nhiệm vụ bài tập lớn được GV hướng dẫn thông qua

- Phần số liệu tính toán nhiệt, động học và động lực học có kèm theo các đồ thị

- Kết luận của bài tập lớn

CÁC THÔNG SỐ ĐỘNG CƠ 1.5L DOHC VTEC TURBO TRÊN XE

HONDA CRV 2020

STT Tên thôngsố Ký hiệu Giá trị Đơn vị Ghi chú

11

Góc đóngmuộn

MỤC LỤC

CHƯƠNG 1 TÍNH TOÁN NHIỆT 1

1.1 Trình tự tính toán 1

1.1.1 Các số liệu ban đầu

1.1.2 Các thông số cần chọn

1.2 Tính toán các quá trình công tác 4

1.2.1 Quá trình nạp

1.2.2 Quá trình nén

1.2.3 Quá trình cháy

1.2.4 Quá trình giãn nở

1.2.5 Tính toán các thông số đặc trưng của chu trình

1.2.6 Vẽ đồ thị công chỉ thị

CHƯƠNG 2: TÍNH TOÁN ĐỘNG HỌC VÀ ĐỘNG LỰC HỌC CỦA CƠ CẤU TRỤC KHUỶU – THANH TRUYỀN 19

2.1 Động học của piston 19

2.2 Động lực học của cơ cấu trục khuỷu –thanh truyền 19

2.2.1 Lực khí thể

2.2.2 Lực quán tính của các chi tiết chuyển động

BẢNG SỐ LIỆU 23

CODE MATLAP 27

TÀI LIỆU THAM KHẢO 31

DANH MỤC HÌNH ẢNH

Hình 1.1: Đồ thị công chỉ thị P-V 19

Hình 2.1: Đồ thị P−φ ,PJ , P 1 22

DANH MỤC BẢNG BIỂU

Bảng 1: Số liệu tính toán nhiệt động cơ 15Bảng 2: Trị số áp suất của MCCT của quá trình nén và giãn nở tính toán 23Bảng 3: Kết quả tính toán động lực học cơ cấu khuỷu trục - thanh truyền 24

CHƯƠNG 1 TÍNH TOÁN NHIỆT 1.1 Trình tự tính toán

1.1.1 Các số liệu ban đầu

Loại động cơ : Xăng, tăng áp

Nhiệt độ không khí nạp mới (To)

Miền Nam nước ta thuộc khi vực nhiệt đới, nhiệt độ trung bình trongngày có thể chọn là tkk = 29𝑜C cho khu vực miền Nam, do đó:

Áp suất khí nạp trước xupap nạp (P k)

Với động cơ tăng áp, ta có thể chọn Áp suất khí nạp trước xupap nạp P k

là:

P k >P0=0,15 MPa

Nhiệt độ khí nạp trước xuppap nạp (Tk)

Đối với động cơ bốn kỳ có két nước làm mát trung gian:

Áp suất cuối quá trình nạp (P a)

Đối với động cơ tăng áp, ta có thể chọn áp suất cuối quá trình nạp P a

Độ tăng nhiệt độ khi nạp mới (ΔT)

Trị số độ tăng nhiệt độ khi nạp mới ΔT căn cứ vào số liệu thực nghiệm.Đối với động cơ xăng, ta chọn trong phạm vi:

ΔT =0÷20℃=20℃

Hệ số nạp thêm (λ1 )

Hệ số nạp thêm λ1 biểu thị sự tương quan lượng tăng tương đối của hỗnhợp khí công tác sau khi nạp thêm so với lượng khí công tác chiếm chỗ ở thểtích Va Thông thường ta có thể chọn trong phạm vi sau:

Chọn hệ số quét buồng cháy (λ¿¿ 2) ¿

Với động cơ tăng áp (có quét buồng cháy), ta chọn hệ số quét buồng cháy

λ2:

Chọn hệ số hiệu đính tỷ nhiệt (λ t)

Hệ số hiệu đính tỷ nhiệt λ t được chọn theo hệ số dư lượng không khí α

đối với động cơ Xăng có α=0,85−0,95 Ta chọn

α=0,9; λ t=1,15

Hệ số lợi dụng nhiệt tại điểm Z (ξ Z )

Hệ số lợi dụng nhiệt tại điểm Z (ξZ) phụ thuộc vào chu trình công tác củađộng cơ Với loại động cơ Xăng ta thường chọn:

ξ z =0,75 ÷ 0,92=0,8

Hệ số lợi dụng nhiệt tại điểm b (ξ b)

Hệ số lợi dụng nhiệt tại điểm b ξ b tùy thuộc vào loại động cơ Xăng hayDiesel Với các loại động cơ Xăng ta chọn:

Hệ số lượng dư không khí (α)

Hệ số lượng dư không khí α tùy thuộc vào loại động cơ Đối với loạiđộng cơ Xăng ta chọn:

α=0,85÷ 0,95=0,9

Hệ số điền đầy đồ thị công (φ d)

Hệ số điền đầy đồ thị công φ d phụ thuộc vào loại động cơ Đối với động

cơ Xăng thường chọn:

Tỷ số tăng áp

Đối với động cơ xăng trị số λ ta chọn

γ r=0,0273

Nhiệt độ cuối quá trình nạp (T a)

Hệ số nạp η v được xác định theo công thức

Tỷ số nén đa biến trung bình (n1 )

Tỷ số nén đa biến trung bình n l được xác định bằng cách giải phươngtrình:

a '

v+b

' v

2 .T a (ε n1 −1 +1)Thay số vào công thức tính n1 ta được:

Nhiệt độ cuối quá trình nén (T¿¿c)¿

Nhiệt độ cuối quá trình nén (T¿¿c)¿ được xác định theo công thức:

T c =T a .ε n1 −1

Thay số vào công thức tínhT c ta được:

T c=361,1.10,6 1,37−1=864,95 K

1.2.3 Quá trình cháy

Lượng không khí lý thuyết cần để đốt cháy 1kg nhiên liệu (M o¿

Lượng không khí lý thuyết cần để đốt cháy 1kg nhiên liệu Mo được tínhtheo công thức :

Lượng khí nạp mới thực tế nạp vào xylanh (M1 )

Đối với động cơ xăng lượng khí nạp mới được tính theo công thức:

M1=α M0+ 1μ

n.l (kmolkk /kg.nl) Chọn μ n l =112 kg/kmol

Hệ số biến đổi phân tử khí lý thuyết (β¿¿ 0) ¿

Ta có β0 được xác định bởi công thức:

β0=M2

Hệ số biến đổi phân tử khí thực tế (β)

Ta có hệ số biến đổi phân từ khí thực tế β xác định theo công thức:

β=1+ β 1+γ0−1

r

Thay số ta được

β=1+ 1,078−11+0,0273=1,076

Hệ số biến đổi phân tử khí tại điểm (Z)

Hệ số biến đổi phân tử khí tại điểm Z được xác định theo công thức:

β z=1+β 1+γ0−1

r ξ ξ z b

Thay số vào ta được

β z=1+ 1,078−11+0,0273 0,80,9=1,0675

Tổn thất nhiệt do cháy không hoàn toàn (∆Q H)

Vì α <1 nên ∆ Q H được tính theo công thức:

∆ Q H=120 10 3.(1−α) M0 (KJ/kg.nl)

Thay số vào ta được:

∆ Q H=120 10 3.( 1−0,9 ) 0,516=6192 (KJ/kg.nl)

Tỷ nhiệt mol đẳng tính trung bình của môi chất tại điểm (Z)

Tỷ nhiệt mol đẳng tích trung bình của môi chất tại điểm Z được tính theo

công thức:

mc}} = {{M} rsub {2} left ({X} rsub {z} + {{γ} rsub {r}} over {{β} rsub {0}} right ) overline {{mc} rsub {v} rsup {vz +M1(1−X Z).mc v ¿

M2

(

β0

)

=a} rsub {vz} + {{1} over {2} b vz .T z

mc}} = {{β} rsub {0} left ({X} rsub {z} + {{γ} rsub {r}} over {{β} rsub {0}} right ) overline {{mc} rsub {v} rsup {vz +(1−X Z).mc v ¿

Nhiệt độ cuối quá trình cháy (T¿¿z)¿

Nhiệt độ cuối quá trình cháy đối với động cơ xăng được tính bằng cáchgiải phương trình sau:

Đối với động cơ xăng áp suất cuối quá trình cháy được tính bằng công

Đối với động cơ xăng: δ= ε ρ=10,6

Xác định chỉ số dãn nở đa biến trung bình (n¿¿ 2) ¿

Ở nhiệt độ từ 1200÷2600K, sai khác của tỷ nhiệt không lớn lắm, do đó ta

Nhiệt độ cuối quá trình giãn nở (T¿¿b)¿

Đối với động cơ xăng:

T b= T z

ε(n2 −1 ) =2662,3166

10,6 1,264−1 =1427,5093(K )

Áp suất cuối quá trình giãn nở (P b)

Đối với động cơ xăng

P b=P z

ε n2 = 11,2701 10,6 1,2640=0,57 (MPa)

Kiểm nghiệm nhiệt độ khí sót (T r)

1.2.5 Tính toán các thông số đặc trưng của chu trình

Áp suất chỉ thị trung bình tính toán (P i '¿

Với động cơ Xăng, áp suất chỉ thị trung bình p i ' được xác định theo côngthức:

λ P=P z

P c= 11,2701

Thay vào công thức ta được:

Là tỷ số giữa phần nhiệt lượng chuyển thành công mà ta thu được vànhiệt lượng mà nhiên liệu tỏa ra khi đốt cháy 1kg nhiên liệu dạng lỏng hay 1 m3

Tính suất tiêu hao nhiên liệu chỉ thị (g i)

Suất tiêu hao nhiên liệu chỉ thị g i được tính theo công thức:

3 )

Thể tích toàn bộ:

Bảng 1: Số liệu tính toán nhiệt động cơ

- Điểm a: điểm cuối quá trình nạp, áp suất P a, thể tích V a

V c= V h ε−1= 0,37510,6−1=0,039(dm3)

𝑃𝑥𝑛, 𝑉𝑥𝑛 là áp suất và thể tích tại một điểm bất kì trên đường congnén

P xn =P a .

(

V xn

)

Bằng cách cho các giá trị 𝑉𝑥𝑛 đi từ 𝑉𝑎 đến 𝑉𝑐 ta lần lượt xác định được cácgiá trị 𝑃𝑥𝑛, kết quả tính toán được ghi trong bảng

Dựng đường cong giãn nở

Trong quá trình giãn nở, khí cháy được giãn nở theo chỉ số đa biến n2 Tương tự ta có:

Điểm r’: x r ' =R

(

)

(

)

V r '=x r ' ⋅ π ⋅B2

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

The tich V (lit)

0 1 2 3 4 5 6 7 8 9

CHƯƠNG 2: TÍNH TOÁN ĐỘNG HỌC VÀ ĐỘNG LỰC HỌC CỦA CƠ

CẤU TRỤC KHUỶU – THANH TRUYỀN 2.1 Động học của piston

2.2 Động lực học của cơ cấu trục khuỷu –thanh truyền

2.2.1 Lực khí thể

Từ đồ thị công chỉ thị, tiến hành vẽ đồ thị lực khí thể Pkt dựa trên nhữnggiá trị áp suất đã có ở đồ thị công chỉ thị P – V và vẽ lại theo góc quay trụckhuỷu α

Quá trình nạp: α = [4°, 180°]

p kt = p a

Quá trình nén: α = [180°, 360° - θ ] = [180°, 330°]

2.2.2 Lực quán tính của các chi tiết chuyển động

 Khối lượng cơ cấu khuỷu trục – thanh truyền

Khối lượng của nhóm piston (khối lượng của các chi tiết chuyển độngtịnh tiến): Dựa vào bảng 2 – 1 giáo trình Động cơ đốt trong 2 trang 24 và đườngkính xilanh ta chọn được:

mnp = 123,9(kg/m2) (piston hợp kim nhôm)

Khối lượng của khuỷu trục (các chi tiết chuyển động quay): Dựa vàobảng 2-1 giáo trình Động cơ đốt trong 2 trang 24 và đường kính xilanh ta chọnđược:

mk = 165,2(kg/m2) (trục khuỷu gang đúc)

Khối lượng nhóm thanh truyền

Dựa vào bảng 2 – 1 giáo trình Động cơ đốt trong 2 trang 24 và đườngkính xilanh ta chọn được:

mtt = 165,2 (kg/m2)

Để đơn giản cho việc tính toán và sai số của nó cũng không đáng kế nên

ta chọn phương pháp dùng khối lượng thay thế Khối lượng thay thế được tínhtheo công thức:

m A =m tt ⋅

(

k

)

L

Do khó khăn trong việc tìm tài liệu về cách bố trí khối lượng thay thế nên

ta chọn theo công thức kinh nghiệm:

với α đi theo từng quá trình tương tự như đối với lực khí thể p kt

 Lực quán tính (lực ly tâm) của khối lượng chuyển động quay

p k =−m r .R ω2

(

m2

)

 Lực tổng cộng p 1 : là lực tổng hợp của lực khí thể và lực quán tính được tính theo công thức:

p1= p kt + p j

Hình 2.1: Đồ thị P−φ ,P J ,P1

BẢNG SỐ LIỆU

Bảng 2: Trị số áp suất của MCCT của quá trình nén và giãn nở tính toán

Quá trình nén (V a −V c) Quá trình giãn nở ¿¿)

Bảng 3: Kết quả tính toán động lực học cơ cấu khuỷu trục - thanh truyền

x2 = R.*((1-cosd(a2))+(lambda/4).*(1- cosd(2.*a2)));

V2 = x2*Fp + Vc;

P2 = Pa.*(Va./V2).^n1;

% qua trinh chay - gian no

% hieu chinh doan c'-c"

a2hc = linspace (340,360,100);

x2hc = R.*((1-cosd(a2hc))+(lambda/4).*(1- cosd(2.*a2hc))); V2hc = x2hc*Fp + Vc;

P1b = Pz.*(Vz./V1b).^n2;

% xac dinh diem b'

%xac dinh diem b"

100 %% do thi P-phi Pj P1

%% xuat bang excel

%%% lap bang ap suat MCCT

%%% lap bang gia tri ket qua tinh toan dong luc hoc

TÀI LIỆU THAM KHẢO

1 Giáo trình Động cơ Đốt trong 1, Nguyễn Văn Trạng, Đại học Sư phạm Kỹthuật Tp HCM, 2005

2 Giáo trình Động cơ Đốt trong 2, Nguyễn Văn Trạng, Đại học Sư phạm Kỹthuật Tp HCM, 2005

3 Lý Vĩnh Đạt, Tính toán nhiệt và động lực học động cơ đốt trong, Đại học sưphạm kỹ thuật thành phố Hồ Chí Minh, 2016