1. Trang chủ
  2. » Luận Văn - Báo Cáo

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

89 19 0

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

THÔNG TIN TÀI LIỆU

Thông tin cơ bản

Định dạng
Số trang 89
Dung lượng 3,15 MB

Nội dung

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

Ngày đăng: 13/12/2021, 23:33

HÌNH ẢNH LIÊN QUAN

Hình 1.2: Lộ trình hơn 40 năm nghiên cứu ô tô điện của Mitsubishi Motors - 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
Hình 1.2 Lộ trình hơn 40 năm nghiên cứu ô tô điện của Mitsubishi Motors (Trang 14)
Hình 1.3: Xe ô tô điện i-MiEV được đưa ra thị trường - 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
Hình 1.3 Xe ô tô điện i-MiEV được đưa ra thị trường (Trang 14)
Hình 1.6:Mức độ phát thải tương đương của các loại động cơ đốt trong và điện - 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
Hình 1.6 Mức độ phát thải tương đương của các loại động cơ đốt trong và điện (Trang 16)
Hình 1.7: Phân bổ khoản đầu tư cho nghiên cứu ô tô điện tại Hoa Kỳ từ năm 2009  - 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
Hình 1.7 Phân bổ khoản đầu tư cho nghiên cứu ô tô điện tại Hoa Kỳ từ năm 2009 (Trang 17)
Hình 1.10: Xe bus điện sử dụng siêu tụ tại Thượng Hải - 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
Hình 1.10 Xe bus điện sử dụng siêu tụ tại Thượng Hải (Trang 20)
Hình 1.12: Chỉ số sử dụng năng lượng của xe ô tô động cơ xăng và điện - 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
Hình 1.12 Chỉ số sử dụng năng lượng của xe ô tô động cơ xăng và điện (Trang 21)
Hình 1.13: Siêu xe điện Tesla đầu tiên được nhập khẩu vào Việt Nam - 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
Hình 1.13 Siêu xe điện Tesla đầu tiên được nhập khẩu vào Việt Nam (Trang 23)
Hình 1.14: Minh họa hệ thống động lực trên xe điện - 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
Hình 1.14 Minh họa hệ thống động lực trên xe điện (Trang 25)
Hình 1.15: Nguyên lý hóa học của tế bào nhiên liệu Fuel Cell (nguồn: Wikipedia).  - 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
Hình 1.15 Nguyên lý hóa học của tế bào nhiên liệu Fuel Cell (nguồn: Wikipedia). (Trang 27)
Hình 1.17: Cấu tạo siêu tụ điện. - 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
Hình 1.17 Cấu tạo siêu tụ điện (Trang 28)
Hình 1.18: Siêu tụ điện và module tụ lớn nhất trên thị trường của Maxwell Technology.  - 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
Hình 1.18 Siêu tụ điện và module tụ lớn nhất trên thị trường của Maxwell Technology. (Trang 29)
Hình 1.19 :Mạch điện đơn giản mô phỏng ắc quy. - 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
Hình 1.19 Mạch điện đơn giản mô phỏng ắc quy (Trang 30)
Hình 1.21: Quan hệ giữa công suất và năng lượng của các loại ắc quy - 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
Hình 1.21 Quan hệ giữa công suất và năng lượng của các loại ắc quy (Trang 32)
Hình 1.23: Phân loại động cơ điện - 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
Hình 1.23 Phân loại động cơ điện (Trang 37)
Hình 2.1: Sơ đồ các lực tác dụng lên xe - 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
Hình 2.1 Sơ đồ các lực tác dụng lên xe (Trang 41)
b Lực cả nh ng hí: là lc do không khí tác dụng lên mặc ca xe khi xe ch y. L c cả  k ô   k    ó          ợc v i chi u chuy     ng - 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 Lực cả nh ng hí: là lc do không khí tác dụng lên mặc ca xe khi xe ch y. L c cả k ô k ó ợc v i chi u chuy ng (Trang 42)
Hình 2.6:Mạch tương đương ắc quy nối với tải - 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
Hình 2.6 Mạch tương đương ắc quy nối với tải (Trang 48)
Hình 2.7: Dòng lượng trên xe ô tô điện có sự thu hồi năng lượng phanh - 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
Hình 2.7 Dòng lượng trên xe ô tô điện có sự thu hồi năng lượng phanh (Trang 51)
Hình 3.2: Khả năng tăng tốc của GMEV1 - 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
Hình 3.2 Khả năng tăng tốc của GMEV1 (Trang 59)
Hình 3.3: Đồ thị vận tốc của GMEV1 chạy theo quy luật mẫu với các gia tốc khác nhau  - 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
Hình 3.3 Đồ thị vận tốc của GMEV1 chạy theo quy luật mẫu với các gia tốc khác nhau (Trang 60)
Hình 3.7: Kết quả khảo sát vớ ia =2,39 m/s2 - 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
Hình 3.7 Kết quả khảo sát vớ ia =2,39 m/s2 (Trang 64)
Hình 3.8: Kết quả khảo sát vớ ia =1,91 m/s2 - 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
Hình 3.8 Kết quả khảo sát vớ ia =1,91 m/s2 (Trang 64)
Hình 3.9: Kết quả khảo sát vớ ia =1,43 m/s2 - 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
Hình 3.9 Kết quả khảo sát vớ ia =1,43 m/s2 (Trang 65)
.Hình 3.10: Kết quả khảo sát vớ ia =0,96 m/s2 - 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
Hình 3.10 Kết quả khảo sát vớ ia =0,96 m/s2 (Trang 65)
Hình 3.12: Kết quả khảo sát với các gia tốc khác nhau - 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
Hình 3.12 Kết quả khảo sát với các gia tốc khác nhau (Trang 66)
Hình 3.11: Kết quả khảo sát vớ ia =0,49m/s2 - 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
Hình 3.11 Kết quả khảo sát vớ ia =0,49m/s2 (Trang 66)
Hình 3.13: Chu trình mẫu với quãng đường 5km với các vận tốc khác nhau - 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
Hình 3.13 Chu trình mẫu với quãng đường 5km với các vận tốc khác nhau (Trang 68)
Hình 3.15: Kết quả khảo sát quãng đường đi được với V=20 km/h - 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
Hình 3.15 Kết quả khảo sát quãng đường đi được với V=20 km/h (Trang 69)
Hình 3.14: Khảo sát cường độ dòng điện trong các chu trình mẫu - 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
Hình 3.14 Khảo sát cường độ dòng điện trong các chu trình mẫu (Trang 69)
Hình 3.19: Kết quả khảo sát quãng đường đi được với V=60 km/h - 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
Hình 3.19 Kết quả khảo sát quãng đường đi được với V=60 km/h (Trang 71)

TÀI LIỆU CÙNG NGƯỜI DÙNG

TÀI LIỆU LIÊN QUAN

w