tính toán nhiệt ứng dụng matlab

clc clear all % LAP TRINH TINH TOAN NHIET DONG CO DIESEL disp('LAP TRINH TINH TOAN NHIET DONG CO DIESEL') % programmer: aaaaaaaaaaaaaaaaaa disp('programmer: aaaaaaaaaaa') %======================================================= % SO LIEU CHO TRUOC disp('SO LIEU CHO TRUOC') %======================================================= % Cong suat dinh muc Ne disp('1 Cong suat dinh muc: Ne [kW]') Ne = 86 % [kW] % Ty so nen e disp('2 Ty so nen: e') e = 19.8 % So vong quay dinh muc n disp('3 So vong quay dinh muc: n [v/p]') n = 3500 % [v/p] % Duong kinh xy lanh: D disp('4 Duong kinh xy lanh: D [mm]') D = 93 % [mm] % Hanh trinh pistong: S disp('5 Hanh trinh pistong: S [mm]') S = 92 % So xy lanh: i disp('6 So xy lanh: i') i=4 % So ky: To % [mm] disp('7 So ky: To') To = % Goc mo som xupap nap: phi1 disp('8 Goc mo som xupap nap: phi1 = 10 do') % Goc dong muon xupap nap: phi2 disp('9 Goc dong muon xupap nap: phi2 = 24 do') % 10 Goc mo som xupap thai: phi3 disp('10 Goc mo som xupap thai: phi3 = 61 do') % 11 Goc dong muon xupap thai: phi4 disp('11 Goc dong muon xupap thai: phi4 = do') % 12 Buong chay ngan cach disp('12 Buong chay ngan cach') % 13 So xupap nap: disp('13 So xupap nap: 2') % 14 So xupap thai: disp('14 So xupap thai: 1') % 15 Dong co tang ap disp('15 Dong co tang ap') % 16 Ap suat nhien lieu bat dau phun: 114-121 bar disp('16 Ap suat nhien lieu bat dau phun: 114-121 bar') % 17 Toc cam chung: 700 +- 20 (vong/phut) disp('17 Toc cam chung: 700 +- 20 (vong/phut)') % 18 Ap suat dau boi tron 4.02 - 4.8 bar disp('18 Ap suat dau boi tron 4.02 - 4.8 bar') %======================================================= % CAC THONG SO CHON disp('CAC THONG SO CHON') %======================================================= % Tinh toc trung binh cua Pistong: Vtb [NLDCDT T99] disp('1 Tinh toc trung binh cua Pistong: Vtb') Vtb = (S*10^-3*n) / 30 % [m/s] % Loai dong co: if Vtb < disp(' Loai dong co: Dong co thap toc ') elseif and(6 = thi: % mcv2 = ((19.876+(1.634/anpha))+((1/2)*(427.86 + (187.36/anpha))*10^-5)*T)/(1+gamar) % Dat E = E = 19.876+(1.634/anpha) % Dat F = F = (1/2)*(427.86 + (187.36/anpha))*(10^-5) % Neu 0.7 x1 Tz = x2 else Tz = x1 end % 5.Ap suat tai diem z: Pz [NLDCDT T174] disp('5.Ap suat tai diem z: Pz') Pz = lamda*Pc %[MPa] % 6.Ty so gian no som: ro disp('6.Ty so gian no som: ro') ro = betaz*(Tz/(lamda*Tc)) % 7.Ty so gian no sau: sma [NLDCDT T182](5-93) % Trong qua trinh gian no, ta co: disp('7.Ty so gian no sau: sma') sma = e/ro %=========================================================== % IV QUA TRINH GIAN NO disp('IV QUA TRINH GIAN NO') % 1.Chi so gian no da bien trung binh: n2 [NLDCDT T184](5-101) disp('Chi so gian no da bien trung binh: n2') %n2 = input('Chon gia tri n2 trong(1.15-1.25),chon n2 = ') n2 = 1.285 % Tb la nhiet tri tai diem b va duoc xac dinh theo cong thuc: Tb = Tz/(20^(n2-1)); % Thay n2 vao hai ve cua phuong trinh tren va so sanh neu sai so giua ve % < 0.2% la dat yeu cau: % Ve trai: O = n2 - 1; % Ve phai: U = 8.314/((((entab-entaz)*Qh)/(M1*(1+gamar)*beta*(Tz-Tb)))+avz2+(bvz2*(Tz+Tb))); y = abs(O-U)*100 disp('Xem cai y o tren no nho hon 0,2% suy gia tri n2 chon la dung') % 2.Ap suat cuoi qua trinh gian no: Pb [NLDCDT T182](5-94) disp('2.Ap suat cuoi qua trinh gian no:: Pb') % Ap suat cuoi qua trinh gian no duoc xac dinh theo cong thuc: Pb = Pz/(sma^n2) % 3.Nhiet cuoi qua trinh gian no: Tb %[MPa] [NLDCDT T182](5-95) disp('3.Nhiet cuoi qua trinh gian no: Tb') Tb = Tz / (e^(n2-1)) %[K] % 4.Kiem tra nhiet sot: Trt disp('4.Kiem tra nhiet sot: Trt') Trt = Tb*((Pr/Pb)^((m-1)/m)) %[K] % Kiem tra: tentaTr = ((abs(Trt-Tr))/Trt)*100 if tentaTr = 4, % duong knh xylanh 90 den 120 mm: % Pm = 0.09 + 0.0120* Cm [MPa] % Doi voi dong co Diezel bon ky 12 xylanh D >= 150 mm, va dong co ky 12 % xylanh cong suat lon: % Pm = 0.03 + 0.0120* Cm [MPa] % Doi voi dong co Diezel o to may keo co buong chay du bi: % Pm = 0.015 + 0.0156* Cm % Doi voi dong co Diezel bon ky ,i >= xy lanh buong chay ngan cach kieu % xoay loc, D = 90 den 120 mm disp('Buong chay ngan cach loai xoay loc') Cm = (92/1000)*3500/30; Pm = 0.09 + (0.0138*Cm) % 6.Ap suat co ich trung binh: pe [NLDCDT T202](6-33) disp('6.Ap suat co ich trung binh: pe') % Do la cong co ich cua 1m^3 the tich cong tac xylanh chu trinh: pe = Pi - Pm %[MPa] % 7.Hieu suat co gioi: nm [NLDCDT T202](6-36) disp('7.Hieu suat co gioi: nm') nm = pe / Pi % 8.Suat tieu hao nhien lieu co ich: ge disp('8.Suat tieu hao nhien lieu co ich: ge') ge = gi / nm %[g/ml.h] % 9.Hieu suat co ich: ne disp('9.Hieu suat co ich: ne') ne = ni * nm % 10.Kiem nghiem duong kinh xy lanh disp('10.Kiem nghiem duong kinh xy lanh') % Tinh the tich Vh : Vhe = (30*To*Ne)/(pe*i*n) % Duong kinh xylanh tinh toan: Dtt = sqrt((4*Vhe)/(pi*S)) % So sanh: TentaD = Dtt*100 - D*100 disp('Xem cai TentaD o tren no nho hon 0,1% la duoc') if TentaD < (1/10) disp(' Duong kinh xylanh tinh toan duoc la dung') else disp(' Duong kinh xylanh tinh toan duoc la sai') end % VI.VE VA HIEU DINH DO THI CONG disp('VI.VE VA HIEU DINH DO THI CONG') disp('Tinh the tich Vc [cm^3]') Vc = (Vh / (e - 1))*1000 %[cm^3] i = [ro 10 11 12 13 14 15 16 17 18 19 19.8]; % Dat Vx = i*Vc , i = ro : e Vx = i.*Vc ; % Qua trinh nen Px = Pc*(1./(i.^n1)); % Qua trinh gian no px = Pz*(ro.^n2)*(1./(i.^n2)); % Ve thi cong P-V figure PV = plot(Vx,Px,Vx,px*0.85); hold on title (' Do Thi Cong P-V') xlabel (' V (cm^3)') ylabel (' P (MPa)') legend('Duong nen','Duong chay gian no') % - END