Nghiên cứu đánh giá ảnh hưởng của quá trình tăng tốc tới khả năng hoạt động của xe điện

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - TRỊNH QUANG HUY NGHIÊN CỨU ĐÁNH GIÁ ẢNH HƯỞNG CỦA QUÁ TRÌNH TĂNG TỐC TỚI KHẢ NĂNG HOẠT ĐỘNG CỦA XE ĐIỆN Chuyên ngành: Kỹ thuật khí động lực LUẬN VĂN THẠC SĨ KỸ THUẬT CHUYÊN NGÀNH: KỸ THUẬT CƠ KHÍ ĐỘNG LỰC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS ĐÀM HOÀNG PHÚC Hà Nội - Năm 2017 MỤC LỤC LỜI CAM ĐOAN .1 MỞ ĐẦU .2 DANH MỤC CÁC KÝ TỰ VÀ KÝ HIỆU VIẾT TẮT TRONG LUẬN VĂN C I T ệ 11 L ệ 12 T ụ ệ ế Hệ ệ ệ 19 131N ă ợ 132Đ ô ô ệ 19 ệ 31 14 N C ởVệ N ế ệ : 33 II: Xâ Mô ô ỏ ô ỏ ệ .35 ọ o ọ ệ .37 ỏ ệ 38 ứ 2.2.2 Các côn ô ỏ ắ ệ 39 : 39 234 T ệ ô ỏ ỏ ắ ệ 2.4.1 Dò ă 242 N ợ ắ Mô .35 ệ 35 ọ ỏ ỏ ệ : .38 2 Đặ Mô ọ ọ 212 P Mô ô ọ 2.1 P â 2 Mô ứ : 45 ô ô ệ ô ỏ ỏ .45 .46 o 2.5 T ậ o ệ : 45 P ầ ệ 46 ô ỏ : .48 T ậ o 48 252P ầ ô C ỏ 31P III: Mô ỏ 50 k ảo .51 k ảo sát .51 I ợng khảo sát: 51 3.2 L a chọ 3.3 Khảo sát khả ă o c khác 54 ò ện chu trình mẫu v i gia t c khác 55 3.3.2 Khảo sát khả ă o ng c a xe v i k = 58 3.3.3 Khảo sát khả ă o ng c a xe v i k=0,8 .58 3.3.4 Khảo sát khả ă o ng c a xe v i k= 0,6 59 3.3.5 Khảo sát khả ă o ng c a xe v i k=0,4 .59 3.3.6 Khảo sát khả ă o ng c a xe v i k=0,2 .60 3.3.1 Khảo ứ ă ng c a xe v 3.3.7 Phân tích kết khảo 3.4 Khảo sát khả ă o 3.4.1 Mơ hình quy luậ ã ợc v i hệ s k khác nhau: 60 ng c a xe ch y t khác 61 i vận t c: 61 3.4.2 Các chu trình mẫu: 61 3.4.3 Khảo ã ợc v i vận t c 20km/h .63 3.4.4 Khảo ã ợc v i vận t c 30km/h .64 3.4.5.Khảo ã ợc v i vận t c 40km/h 64 3.4.6 Khảo ã ợc v i vận t c 50km/h .65 3.4.7 Khảo ã ợc v i vận t c 60km/h .65 3.4.8 Phân tích kết khảo ã ợc v i vận t c khác nhau: 66 3.5 Phân tích kết luận kết khảo sát .66 Kế ậ .68 II DANH MỤC HÌNH ẢNH 40 ă Hình 1.2: L H 3: X ô ô Hình 1.4:Mứ ứ ô ô ện i-M EV ện c a Mitsubishi Motors ợ ờng phát thải khí CO2 a lo t ện ă Hình 1.5:Hiệu su ợ o a loa ện .10 Hình 1.6:Mứ H phát thả 7: P â b k oả o a lo ầ o ứ ô ô ệ ện 10 Ho Kỳ ă 2009 .11 Hình 1.8: Tesla – Kế ả ợ o ỗ ệ Mỹ kỳ Obama .11 Hinh 9: X H ệ ỡ 10: X b ỏ b ế N U .12 ệ ụ ụ Hình 1.11: Angkor EV 2013 – Xe ă Hình 1.12: Chỉ s s dụ H 13: S H 14: M H 15: N ệ ợ ợng c ệ T ọ T ợ Hả 14 ả Campuchia .14 ô ô ầ ợ ă ậ k ẩ ệ ện .15 oVệ N .17 ệ .19 ý ó ọ ếb o ệ F C ( :Wk ) .21 H 16: M H 17: C H 18: S ọ ệ o ếb o ụ ụ ệ ệ F C ô ô ệ 21 22 ệ o ụ a Maxwell Technology 23 Hình 1.19 : M H 20: P â H 21: Q Hình 1.22: Biế ệ ả o ắ ệ o ô ic ô ỏ ắ ậ 24 ă ă ợ 26 ợ ợng phóng ắc quy v inhiệ III o ắ 26 .28 ện 31 Hình 1.23: Phân lo H 24: Đ H 1: S H 22S H 3: Đặ H 4: M H 5: M H 6:M ắ H 7: Dị Hình 2.8: S H 1: X o ệ 31 ụng lên xe 35 ệ o ế b 36 ệ .38 ả .40 ợ ầ ả 42 ô ô ậ o 40 ệ ó ă ợ 45 49 ện GM EV1 General Motor sản xu t 52 Hình 3.2: Khả ă H ặ ă c c a GM EV1 .53 3: Đồ th vận t c c a GM EV1 ch y theo quy luật mẫu v i gia t c khác 54 H 4: Đồ th vận t c c a GM EV1 ch y trong1 chu trình mẫu v i gia t c khác 55 H 5:C ò ện phóng m t chu trình mẫu v i gia t c khác 56 Hình 3.6: S i c a I max theo gia t c m t chu kỳ 57 Hình 3.7: Kết khảo sát v i a =2,39 m/s2 58 Hình 3.8: Kết khảo sát v i a =1,91 m/s2 58 Hình 3.9: Kết khảo sát v i a =1,43 m/s2 59 Hình 3.10: Kết khảo sát v i a =0,96 m/s2 59 Hình 3.11: Kết khảo sát v i a =0,49m/s2 60 Hình 3.12: Kết khảo sát v i gia t c khác 60 ã ờng 5km v i vận t c khác 62 ò ện chu trình mẫu 63 Hình 3.13: Chu trình mẫu v Hình 3.14: Khảo Hình 3.15: Kết khảo ã ợc v i V=20 km/h 63 Hình 3.16: Kết khảo ã ợc v i V=30 km/h 64 IV Hình 3.17: Kết khảo ã ợc v i V=40 km/h 64 Hình 3.18: Kết khảo ã ợc v i V=50 km/h 65 Hình 3.19: Kết khảo ã ợc v i V= 60 km/h 65 Hình 3.20: T ng hợp kết khảo v i vận t c khác 66 V LỜI CAM ĐOAN Tô b o o ậ ă ú ỡ ã N C ợ b bè, ă ký ệ , kế ô ứ o ầ Hồng Phúc, â ệ o, o ệ ẫ k o Đ N ậ ọ B ậ ệ ả o ô N Hệ ă ọ TS Đàm K o H N o o T Đ ù ù ợ ọ B K o H ă T ả Trịnh Quang Huy MỞ ĐẦU C ú ả kệ ặ ệ ặ ó bằ o ĩ o T o ô b ô ậ k , b , ợ , â ễ , o ế ợ k oả 8% Đứ ), k ă V ệ (E ă ô ậ ô ặ bệ T ệ ( ũ ă ệ ă k ả ă ụ óý o ĩ â o ứ … Đó ế o ụ ă ợ o, ó k ả ă ệ giao ả ợ ý o ô o ă ế k ả ă ũ ó ă k ả ă ó ụ ọ ợ o ế k ó o , ó o ă , ỏ, k o ó ế ệ k ả ă ù ợ ó ệ ả 5685 k ệ ợ b ứ ệ ( ầ ợ ậ ệ ệ 1980-2011 Cùng ệ ứ ệ ệ ó ứ o ) o ) K ả ă ọ , â ế ệ V …X ó o o ệ ệ ệ ,X ả ả k , o ù 2565 ợ , ặ ả, ế o b ò :K ô ă o o ợ ứ ệ , ọ ắ ệ ó ă ó (IEA 2012) ế ợ ă ứ kệ ô o o ẽ giả ả ă o ệ CO2 o 2,6% ô ệ ả, ặ bệ o ă ứ o ò o ă ợ , ậ N C 1: T C 2: Xâ C 3: Mơ ă b o ó3 : ệ ô ỏ ứ ô ỏ k ảo Kế L ậ V ỗ k ỏ ế ô N â ô t ầ bả P ú â , ậ b ứ ụ N â ú ọ ,Vệ ỡ ô o ậ bế â ả Vệ C k Đ ã ế ò ỏ ị Tơ ã o ó , mong ã ệ ă ậ ợ â ắ ế ế ả ầ o o Đ ọ ậ ó ý ầ ẫ ô ã ặ , ô o b ọ T ă k ô â o TS Đ o Ho ệ o ậ ô ô ô Đ ọ B k o H ứ H N , N … … ă ệ Trịnh Quang Huy ă 2017 DANH MỤC CÁC KÝ TỰ VÀ KÝ HIỆU VIẾT TẮT TRONG LUẬN VĂN TT Ký tự Giải thích ký tự Đơn vị SPbatt Mật độ công suất ắc quy W/kg SEbatt Mật độ lượng ắc quy Wh/ kg Edis Năng lượng xả Joule Echg Năng lượng nạp Joule EV Xe ô tô điện HEV Xe tơ hybrid điện ICV Ơ tơ trang bị động đốt ESD Năng lượng tổn thất tự phóng ắc quy AC Điện xoay chiều 10 DC Điện chiều 11 BLDC 12 m Khối lượng kg 13 a Gia tốc m/s2 14 Voc Điện áp hở mạch V 15 A Diện tích cản gió m2 16 r Bán kính bánh xe hiệu dụng m 17 Frr Lực cản lăn N Động điện chiều không chổi than Hệ số cản lăn 18 19 Wh Fad Lực cản khơng khí N 20 Mật độ khơng khí Ω 21 Hệ số cản khí động học Ah 22 Vận tốc dài xe m/s 23 Vận tốc gió m/s 24 Lực cản dốc N 25 Fla Lực cản quán tính N 26 Fte Lực kéo động truyền đến bánh xe N 27 Fωa Lực cản quán tính khối lượng chuyển động quay N Tài liêụ tham khảo [1] “Electric Vehicle Technology Explained”, James Larminie Oxford Brookes University, Oxford, Uk; John Lowry - Acenti Designs Ltd, UK; 2003 [2] “Kết cấu ô tô ” Nguyễn Khắc Trai, Nguyễn Trọng Hoan, Hồ Hữu Hải, Phạm Huy Hường, Nguyễn Văn Chưởng, Trịnh Minh Hoàng, Nhà xuất Bách khoa Hà Nội,2009 Phụ lục 1.Công thức mô ác quy chì –axit hở mạch: open_circuit_voltage_LA.m function E_oc=open_circuit_voltage_LA(DoD,N) % find the open circuit voltage of lead acid if DoD < error('Depth of discharge 1') end % see aqu 2.10 in text E_oc= (2.15-((2.15-2.00)*DoD))*N; Khảo sát khả tăng performmance_simulation_GM.m tốc xe GM % GMEV1 ***************** % Simulates the WOT test of the GM EV1 electric car t=linspace(0,15,151); % to 15 seconds, in 0.1 sec steps v=zeros(1,151); % 151 readings of velocity dT=0.1; % 0.1 second time step % In this case there are three phases to acceleration, as % explained in the text for n=1:150 if v(n)30.5 % Controller stop any more speed increase v(n+1)=v(n); else % Equ 7.22 v(n+1)=v(n)+dT*((62.1/v(n))-0.046-(0.000137*(v(n)^2))); end; end; v=v*3.6; % Multiply all v values by 3.6 to convert to kph plot(t,v); xlabel('Time/seconds'); ylabel('velocity/kph'); title('Full power (WOT) acceleration of GM EV1 electric car'); Mô chu kỳ xe: % % % % % % % % one_cylce_oto.m *********************** ONE CYCLE This script file performs one cycle, of any drive cycle of N points with any vehicle and for lead acid or NiCad batteries All the appropriate variables must be set by the calling program ********************** % V_ECE47; V_GM; %V_array; N=length(V) V=V./3.6; % % % first we set up the vehicle data mass=1540; % vehicle mass+ two 70 kg passengers area=1.8; % frontal area in square metres Cd= 0.19; % Drag ceofficient Gratio = 37; % Gearing ratio, = G/r G_eff= 0.95; % Transmission effi Regen_ratio = 0.5; % this set the proportion of the braking that is done regeneratively using the motor bat_type='LA'; % lead acid battery Nocells=156; % 26 of cell (12 volt) batteries Capacity=60; % 60 Ah batteries This is assumed to be the 10 hour rate capacity k=1.12; Pac=250; % Peukert % average power of accessories % % % These are the constant for the motor efficiency % % equation, 7.23 kc= 0.3; ki=0.01; kw=0.000005; ConL=600; % % some constants which are calculated Frr= 0.0048*mass*9.8; % 7.1 Rin= (0.022/Capacity)*Nocells; % 2.2 Rin=Rin+0.05; PeuCap= ((Capacity/10)^k)*10; % 2.12 % Set up arrays for storing data for battery DoD_end= zeros(1,100); CR_end=zeros(1,100); D_end=zeros(1,100); % We now need similar arrays for use within each cycle DoD=zeros(1,N); % Depth of discharge, as in chap CR=zeros(1,N); % charge removes from battery, Peukert D=zeros(1,N); % record of distance traveled in km % Pte(1)=0; % Pmot_in(1)=0; % no power into motor for C = 2:N accel=V(C)-V(C-1); Fad=0.5*1.25*area*Cd*V(C)^2; % Equ 7.2 Fhc=0; % Equ 7.3 Fla=1.05*mass*accel; % The mass is increased modestly to compensate for % the fact that we have excluded the moment of inertia Pte(C)= (Frr+Fad+Fhc+Fla)*V(C); % Equ 7.9 & 7.23 omega = Gratio*V(C-1); if omega ==0 % stationary Pte(C)=0; Pmot_in(C)=0; % no power into motor Torque=0; eff_mot=0.5; % Dummy value, to make sure not zero elseif omega>0 % moving if Pte(C)=0 Pmot_out(C) = Pte(C)/G_eff; % Motor power > shaft power elseif Pte(C)0 % Now use equ 7.23 eff_mot=(Torque*omega)/((Torque*omega)+((Torque^2)*kc)+(omega*ki)+(omega^ 3)*kw+600); elseif Torque= Pmot_in(C) = Pmot_out(C)/eff_mot; % Equ 7.23 elseif Pmot_out(C) < Pmot_in(C) = Pmot_out(C)*eff_mot; end; end; Pbat(C) = Pmot_in(C-1) + Pac; % equ 7.26 if bat_type == 'NC' E(C)=open_circuit_voltage_NC(DoD(C-1),Nocells); elseif bat_type == 'LA' E(C)=open_circuit_voltage_LA(DoD(C-1),Nocells); else error('Invalid battery type'); end; if Pbat(C) > % use equ 2.20 I(C) =(E(C)-((E(C)*E(C))-(4*Rin*Pbat(C)))^0.5)/(2*Rin); CR(C)= CR(C-1)+ ((I(C)^k)/3600); % equ 2.18 elseif Pbat(C)==0 I(C)=0; elseif Pbat(C)1 DoD(C)=1; end; % Since we are taking one second time intervals, the distance traveled in % metres is same as velocity Divide by 1000 for km D(C)=D(C-1)+(V(C)/1000); XDATA(C)=C; % See section 7.4.4 YDATA(C)=eff_mot; % of these two arrays end; figure(2) plot(I); % figure(3) % plot(t,E); % figure(4) % plot(t,Pbat) % figure(5) % plot(t,Pmot_in) % Now return to calling program Mô khả hoạt động xe : range_simulation_oto.m % Siluation of the GM EV1 running the SFUDS % driving cycle this simulation is for range % measurement The run continues until the battery depth of discharge >90% % V_ECE47; % Get the velocity values ECE47, they are in %V_GM; % an arrays V %thu_vantoc2, %a5, V_array; N=length(V); % find out how many readings % Divide all velocities by 3.6, to conver to m/sec V=V./3.6; % first we set up the vehicle data mass=1540; % vehicle mass+ two 70 kg passengers area=1.8; % frontal area in square metres Cd= 0.19; % Drag ceofficient Gratio = 37; % Gearing ratio, = G/r G_eff= 0.95; % Transmission effi Regen_ratio = 0.5; % this set the proportion of the braking that is done regeneratively using the motor bat_type='LA'; % lead acid battery Nocells=156; % 26 of cell (12 volt) batteries Capacity=60; % 60 Ah batteries This is assumed to be the 10 hour rate capacity k=1.12; % Peukert Pac=250; % average power of accessories % These are the constant for the motor efficiency % equation, 7.23 kc= 0.3; ki=0.01; kw=0.000005; ConL=600; % some constants which are calculated Frr= 0.0048*mass*9.8; % 7.1 Rin= (0.022/Capacity)*Nocells; % 2.2 Rin=Rin+0.05; PeuCap= ((Capacity/10)^k)*10; % 2.12 % Set up arrays for storing data for battery DoD_end= zeros(1,100); CR_end=zeros(1,100); D_end=zeros(1,100); % We now need similar DoD=zeros(1,N); % CR=zeros(1,N); % D=zeros(1,N); % arrays for use within each cycle Depth of discharge, as in chap charge removes from battery, Peukert record of distance traveled in km CY=1; % CY controls the outer loop, and counts the number % of cycle completed We want to keep cycling till the battery flat This % we define as being more than 90% discharged That is, DoD_end >0.9 % We also use the variable XX to monitor the dischager, and stop the loop % going too far DD=0; while DD=19.8 vel(n+1)=vel(n)+a*dT*(62.1/vel(n)-0.046-(0.000137*(vel(n)^2))); end; d(n+1)=d(n) + 0.1*vel(n); % Compute distance traveled end; vel=vel*3.6; % Multiply by 3.6 to convert m/sec to kph % figure; plot(t,vel); % % % % %axis([0 150 70]); xlabel('Time/seconds'); ylabel('Velocity/kph'); title(' acceleration of electric GM EV1'); %figure; %plot(t,d); axis([0 30 140]); Khảo sát khả hoạt động xe với gia tốc khác nhau: two_cyles.m; speedup.m 6.1 two_cyles.m n=zeros(1,311); vel=zeros(1,311); dT=1; t=linspace(1,310,311); a=2.39; for n=1:15 if vel(n) 16.7 vel(n+1)=vel(n); end; end for n=16:136 vel(n+1)=vel(n); end decel=(vel(137)-0)/15; for n=137:151 vel(n) = vel(n-1) - decel; end; for n= 151:156 vel(n)=0 end for n=156:171 if vel(n) 16.7 vel(n+1)=vel(n); end; end for n=171:291 vel(n+1)=vel(n); end for n=291:306 vel(n) = vel(n-1) - decel; end; for n= 306:311 vel(n)=0; end vel=vel*3.6; % Multiply by 3.6 to convert m/sec to kp plot(t,vel) 6.2 speedup.m %Creating driving patern for analysing current in accelerating phase %acceleration must be lower than 2,39m/s2 n=zeros(1,151); % step = second, cycles last in 135 seconds V=zeros(1,151); dT=1; t=linspace(1,150,151); a=2.39; % Maximum acceleration of GMEV1 k=1; % investigation constants k= full power % k= 0.8 80% full power for n=1:136 % speed up the car to 60 km/h then remains this speed in 135 s if V(n)= 16.7 V(n+1)=V(n); end; end; decel=(V(136)-0)/15; for n=137:151 % slow down the car to km/h in 15 seconds V(n) = V(n-1) - decel; end; V=V*3.6; % Multiply by 3.6 to convert m/sec to km/h plot(V) xlabel('Time/seconds'); ylabel('Velocity/kph'); title('one cycle'); Khảo sát cường độ dòng điện với gia tốc khác nhau: % % % % % % % % Tính I 1cycle v?i V=60 km/h ONE CYCLE This script file performs one cycle, of any drive cycle of N points with any vehicle and for lead acid or NiCad batteries All the appropriate variables must be set by the calling program ********************** % V_ECE47; %V_GM; %speedup02; %v60_test; a_test %V_array; N=length(V); V=V./3.6; % % % first we set up the vehicle data mass=1540; % vehicle mass+ two 70 kg passengers area=1.8; % frontal area in square metres Cd= 0.19; % Drag ceofficient I_a_test.m Gratio = 37; % Gearing ratio, = G/r G_eff= 0.95; % Transmission effi Regen_ratio = 0.5; % this set the proportion of the braking that is done regeneratively using the motor bat_type='LA'; % lead acid battery Nocells=156; % 26 of cell (12 volt) batteries Capacity=60; % 60 Ah batteries This is assumed to be the 10 hour rate capacity k=1.12; % Peukert Pac=250; % average power of accessories % % % These are the constant for the motor efficiency % % equation, 7.23 kc= 0.3; ki=0.01; kw=0.000005; ConL=600; % % some constants which are calculated Frr= 0.0048*mass*9.8; % 7.1 Rin= (0.022/Capacity)*Nocells; % 2.2 Rin=Rin+0.05; PeuCap= ((Capacity/10)^k)*10; % 2.12 % Set up arrays for storing data for battery DoD_end= zeros(1,100); CR_end=zeros(1,100); D_end=zeros(1,100); % We now need similar arrays for use within each cycle DoD=zeros(1,N); % Depth of discharge, as in chap CR=zeros(1,N); % charge removes from battery, Peukert D=zeros(1,N); % record of distance traveled in km % Pte(1)=0; % Pmot_in(1)=0; % no power into motor for C = 2:N accel=V(C)-V(C-1); Fad=0.5*1.25*area*Cd*V(C)^2; % Equ 7.2 Fhc=0; % Equ 7.3 Fla=1.05*mass*accel; % The mass is increased modestly to compensate for % the fact that we have excluded the moment of inertia Pte(C)= (Frr+Fad+Fhc+Fla)*V(C); % Equ 7.9 & 7.23 omega = Gratio*V(C-1); if omega ==0 % stationary Pte(C)=0; Pmot_in(C)=0; % no power into motor Torque=0; eff_mot=0.5; % Dummy value, to make sure not zero elseif omega>0 % moving if Pte(C)=0 Pmot_out(C) = Pte(C)/G_eff; % Motor power > shaft power elseif Pte(C)0 % Now use equ 7.23 eff_mot=(Torque*omega)/((Torque*omega)+((Torque^2)*kc)+(omega*ki)+(omega^ 3)*kw+600); elseif Torque= Pmot_in(C) = Pmot_out(C)/eff_mot; % Equ 7.23 elseif Pmot_out(C) < Pmot_in(C) = Pmot_out(C)*eff_mot; end; end; Pbat(C) = Pmot_in(C-1) + Pac; % equ 7.26 if bat_type == 'NC' E(C)=open_circuit_voltage_NC(DoD(C-1),Nocells); elseif bat_type == 'LA' E(C)=open_circuit_voltage_LA(DoD(C-1),Nocells); else error('Invalid battery type'); end; if Pbat(C) > % use equ 2.20 I(C) =(E(C)-((E(C)*E(C))-(4*Rin*Pbat(C)))^0.5)/(2*Rin); CR(C)= CR(C-1)+ ((I(C)^k)/3600); % equ 2.18 elseif Pbat(C)==0 I(C)=0; elseif Pbat(C)1 DoD(C)=1; end; % Since we are taking one second time intervals, the distance traveled in % metres is same as velocity Divide by 1000 for km D(C)=D(C-1)+(V(C)/1000); XDATA(C)=C; % See section 7.4.4 YDATA(C)=eff_mot; % of these two arrays end; %figure(2) 11 plot(I);%axis([0 950 -100 400]); xlabel('Time/seconds'); ylabel('I/ampere'); title('discharge current'); % figure(3) % plot(t,E); % figure(4) % plot(t,Pbat) % figure(5) % plot(t,Pmot_in) % Now return to calling program Mô hoạt động xe với vận tốc khác nhau: v-test.m %GM EV1 accelaeration %This file is for constructing the accelaeration profile for GM EV1 t=linspace(0,380,381); % to 50 s, in 0.1 s steps vel=zeros(1,381); % 501 readings of velocity acc=zeros(1,381); d=zeros(1,381);% Array for storing distance traveled dT=0.1; % 0.1 second time step a= 1; for n= 1:380 % Now follow equations 7.17 & 7.18 if vel(n) 35.8 % Controller stop any more speed increase vel(n+1)=vel(n); else vel(n)>=19.8; vel(n+1)=vel(n)+a*dT*(62.1/vel(n)-0.046-(0.000137*(vel(n)^2))); % vel(n+1) = vel(n) + dT*a; end; d(n+1)=d(n) + 0.1*vel(n); % Compute distance traveled end; vel=vel*3.6; % Multiply by 3.6 to convert m/sec to kph figure; plot(t,vel); axis([0 400 200]); xlabel('Time/seconds'); ylabel('Velocity/kph'); title(' acceleration of electric GM EV1'); %figure; %plot(t,d); axis([0 30 140]); Khảo sat cường độ dòng điện với vận tôc khác nhau: % % % % % % Tính I 1cycle v?i V=60 km/h ONE CYCLE This script file performs one cycle, of any drive cycle of N points with any vehicle and for lead acid or NiCad batteries All the appropriate variables must be set 12 Iv60_test.m % by the calling program % ********************** % V_ECE47; %V_GM; %speedup02; v60_test; %a_test %V_array; N=length(V); V=V./3.6; % % % first we set up the vehicle data mass=1540; % vehicle mass+ two 70 kg passengers area=1.8; % frontal area in square metres Cd= 0.19; % Drag ceofficient Gratio = 37; % Gearing ratio, = G/r G_eff= 0.95; % Transmission effi Regen_ratio = 0.5; % this set the proportion of the braking that is done regeneratively using the motor bat_type='LA'; % lead acid battery Nocells=156; % 26 of cell (12 volt) batteries Capacity=60; % 60 Ah batteries This is assumed to be the 10 hour rate capacity k=1.12; % Peukert Pac=250; % average power of accessories % % % These are the constant for the motor efficiency % % equation, 7.23 kc= 0.3; ki=0.01; kw=0.000005; ConL=600; % % some constants which are calculated Frr= 0.0048*mass*9.8; % 7.1 Rin= (0.022/Capacity)*Nocells; % 2.2 Rin=Rin+0.05; PeuCap= ((Capacity/10)^k)*10; % 2.12 % Set up arrays for storing data for battery DoD_end= zeros(1,100); CR_end=zeros(1,100); D_end=zeros(1,100); % We now need similar arrays for use within each cycle DoD=zeros(1,N); % Depth of discharge, as in chap CR=zeros(1,N); % charge removes from battery, Peukert D=zeros(1,N); % record of distance traveled in km % Pte(1)=0; % Pmot_in(1)=0; % no power into motor for C = 2:N accel=V(C)-V(C-1); Fad=0.5*1.25*area*Cd*V(C)^2; % Equ 7.2 Fhc=0; % Equ 7.3 Fla=1.05*mass*accel; % The mass is increased modestly to compensate for % the fact that we have excluded the moment of inertia Pte(C)= (Frr+Fad+Fhc+Fla)*V(C); % Equ 7.9 & 7.23 13 omega = Gratio*V(C-1); if omega ==0 % stationary Pte(C)=0; Pmot_in(C)=0; % no power into motor Torque=0; eff_mot=0.5; % Dummy value, to make sure not zero elseif omega>0 % moving if Pte(C)=0 Pmot_out(C) = Pte(C)/G_eff; % Motor power > shaft power elseif Pte(C)0 % Now use equ 7.23 eff_mot=(Torque*omega)/((Torque*omega)+((Torque^2)*kc)+(omega*ki)+(omega^ 3)*kw+600); elseif Torque= Pmot_in(C) = Pmot_out(C)/eff_mot; % Equ 7.23 elseif Pmot_out(C) < Pmot_in(C) = Pmot_out(C)*eff_mot; end; end; Pbat(C) = Pmot_in(C-1) + Pac; % equ 7.26 if bat_type == 'NC' E(C)=open_circuit_voltage_NC(DoD(C-1),Nocells); elseif bat_type == 'LA' E(C)=open_circuit_voltage_LA(DoD(C-1),Nocells); else error('Invalid battery type'); end; if Pbat(C) > % use equ 2.20 I(C) =(E(C)-((E(C)*E(C))-(4*Rin*Pbat(C)))^0.5)/(2*Rin); CR(C)= CR(C-1)+ ((I(C)^k)/3600); % equ 2.18 elseif Pbat(C)==0 I(C)=0; elseif Pbat(C)1 DoD(C)=1; end; % Since we are taking one second time intervals, the distance traveled in % metres is same as velocity Divide by 1000 for km D(C)=D(C-1)+(V(C)/1000); XDATA(C)=C; % See section 7.4.4 YDATA(C)=eff_mot; % of these two arrays end; %figure(2) plot(I,'r');axis([0 950 -100 400]); xlabel('Time/seconds'); ylabel('I/ampere'); title('discharge current'); % figure(3) % plot(t,E); % figure(4) % plot(t,Pbat) % figure(5) % plot(t,Pmot_in) % Now return to calling program 15 ... .58 3.3.4 Khảo sát khả ă o ng c a xe v i k= 0,6 59 3.3.5 Khảo sát khả ă o ng c a xe v i k=0,4 .59 3.3.6 Khảo sát khả ă o ng c a xe v i k=0,2 .60 3.3.1 Khảo ứ ă ng c a xe v 3.3.7... ợng khảo sát: 51 3.2 L a chọ 3.3 Khảo sát khả ă o c khác 54 ị ện chu trình mẫu v i gia t c khác 55 3.3.2 Khảo sát khả ă o ng c a xe v i k = 58 3.3.3 Khảo sát khả ă o ng c a xe v... kg Edis Năng lượng xả Joule Echg Năng lượng nạp Joule EV Xe ô tô điện HEV Xe ô tô hybrid điện ICV Ơ tơ trang bị động đốt ESD Năng lượng tổn thất tự phóng ắc quy AC Điện xoay chiều 10 DC Điện chiều