✅作者简介热爱科研的Matlab仿真开发者擅长毕业设计辅导、数学建模、数据处理、建模仿真、程序设计、完整代码获取、论文复现及科研仿真。 关注我领取海量matlab电子书和数学建模资料个人信条格物致知,完整Matlab代码获取及仿真咨询内容私信。 内容介绍一、PHEV 的发展背景与优势能源与环境驱动随着全球能源需求的增长和环境问题的日益突出传统燃油汽车面临着严峻的挑战。石油资源的有限性以及燃烧化石燃料所带来的环境污染促使汽车行业寻求更可持续的动力解决方案。插电式混合动力汽车PHEV作为一种过渡性技术结合了传统燃油发动机和电动驱动系统的优势既能利用电力实现零排放行驶又能在电量不足时依靠燃油发动机继续行驶有效缓解了消费者对续航里程的焦虑。政策与市场推动各国政府纷纷出台相关政策鼓励发展新能源汽车PHEV 因其较低的排放和较好的燃油经济性受到了政策的支持。同时消费者对环保和节能汽车的需求不断增加使得 PHEV 在市场上具有较大的发展潜力。二、PHEV 动力系统组件及原理发动机工作原理PHEV 中的发动机与传统燃油汽车发动机类似通常为内燃机常见的有汽油发动机或柴油发动机。以汽油发动机为例它通过将汽油与空气混合后在气缸内点火燃烧产生高温高压气体推动活塞运动进而带动曲轴旋转将化学能转化为机械能。作用发动机在 PHEV 中主要用于在电池电量不足时提供动力确保车辆能够持续行驶。此外在一些工况下发动机还可以与电机协同工作共同为车辆提供动力以满足更高的功率需求例如在加速或爬坡时。电机工作原理电机是 PHEV 实现电动驱动的关键组件常见的有永磁同步电机或异步电机。以永磁同步电机为例当定子绕组通入三相交流电时会产生一个旋转磁场该磁场与转子中的永磁体相互作用使转子跟随旋转磁场转动从而将电能转化为机械能。作用电机可以在纯电动模式下单独驱动车辆实现零排放行驶。它具有响应速度快、扭矩大的特点能够为车辆提供良好的加速性能。此外在车辆制动过程中电机还可以作为发电机使用将车辆的动能转化为电能回收到电池组中实现能量的回收利用提高车辆的能源效率。电池组工作原理电池组是 PHEV 储存电能的部件目前常用的是锂离子电池。锂离子电池通过锂离子在正负极之间的嵌入和脱嵌来实现充放电过程。在充电过程中锂离子从正极脱出经过电解液嵌入负极放电时锂离子则从负极脱出回到正极。作用电池组为电机提供电能支持车辆的纯电动行驶。其容量和性能直接影响车辆的纯电续航里程和动力性能。此外电池组还需要具备良好的充放电效率、安全性和耐久性以满足车辆的日常使用需求。三、PHEV 动力系统建模组件建模为了准确分析和设计 PHEV 动力系统需要对各个组件进行建模。发动机建模通常基于其热力学和动力学原理考虑进气、燃烧、排气等过程以及发动机的扭矩、功率输出特性与转速、负荷之间的关系。电机建模则基于电磁学原理考虑电机的电磁转矩、转速、电流、电压等参数之间的关系以及电机的效率特性。电池组建模需要考虑电池的电化学特性如电池的充放电容量、内阻、开路电压与 SOCState of Charge荷电状态的关系等。系统集成建模在对各个组件建模的基础上还需要将它们集成起来建立 PHEV 动力系统的整体模型。这个模型需要考虑各个组件之间的相互作用和能量流。例如发动机与电机之间的功率分配电池组与电机之间的电能传递以及车辆行驶过程中的能量消耗和回收等。通过系统集成建模可以对 PHEV 在不同工况下的性能进行模拟和分析为动力系统的设计和优化提供依据。四、PHEV 动力系统设计目标设定PHEV 动力系统设计的目标通常包括提高车辆的燃油经济性、降低排放、增加纯电续航里程以及保证车辆的动力性能。这些目标之间往往相互制约需要在设计过程中进行权衡。例如增加电池组容量可以提高纯电续航里程但会增加车辆的重量和成本可能对燃油经济性和动力性能产生一定影响。参数优化根据设定的目标需要对动力系统的各个参数进行优化。这包括发动机的排量、压缩比、燃烧策略电机的功率、扭矩特性、效率以及电池组的容量、电压等级等。通过优化这些参数可以使动力系统在不同工况下都能达到较好的性能表现。例如合理选择发动机和电机的功率匹配可以使车辆在城市拥堵工况下更多地使用纯电模式而在高速行驶时发动机和电机协同工作提高能源利用效率。控制策略设计PHEV 动力系统的控制策略对于实现其性能目标至关重要。控制策略需要根据车辆的行驶工况如车速、加速度、负载等、电池 SOC 以及驾驶员的需求实时调整发动机、电机和电池组之间的工作状态实现动力的合理分配和能量的高效管理。常见的控制策略有基于规则的控制策略和优化算法驱动的控制策略。基于规则的控制策略根据预先设定的规则来决定发动机和电机的工作模式例如当电池 SOC 高于某一阈值且车速较低时采用纯电模式当电池 SOC 低于某一阈值或需要较大功率时发动机启动。优化算法驱动的控制策略则通过数学优化算法如动态规划、模型预测控制等在满足车辆行驶需求的前提下最小化燃油消耗或排放。⛳️ 运行结果 部分代码clear allclose allrun(Drive_cycles);%storing various energy consumption results for UDDS 505[UDDS_505_Net_tractive_energy,UDDS_505_Fuel_energy_WhpKm,...UDDS_505_Total_energy_consumption_WhpKm,...UDDS_505_Total_energy_consumption_mpgge,...UDDS_505_WTW_PEU,UDDS_505_WTW_GHG,...UDDS_505_Range_km,Total_mass]energy_cunsumption_results(UDDS_505);%storing various energy consumption results for HWFET[HWFET_Net_tractive_energy,HWFET_Fuel_energy_WhpKm,...HWFET_Total_energy_consumption_WhpKm,...HWFET_Total_energy_consumption_mpgge,...HWFET_WTW_PEU,HWFET_WTW_GHG,...HWFET_Range_km,Total_mass]energy_cunsumption_results(HWFET);%storing various energy consumption results for US06 city[US06_city_Net_tractive_energy,US06_city_Fuel_energy_WhpKm,...US06_city_Total_energy_consumption_WhpKm,...US06_city_Total_energy_consumption_mpgge,...US06_city_WTW_PEU,US06_city_WTW_GHG,...US06_city_Range_km,Total_mass]energy_cunsumption_results(US06_city);%storing various energy consumption results for US06 highway[US06_highway_Net_tractive_energy,US06_highway_Fuel_energy_WhpKm,...US06_highway_Total_energy_consumption_WhpKm,...US06_highway_Total_energy_consumption_mpgge,...US06_highway_WTW_PEU,US06_highway_WTW_GHG,...US06_highway_Range_km,Total_mass]energy_cunsumption_results(US06_highway);%creating tableenergy_cnsmp{Net tractive energy Wh/km;Fuel energy_Wh/km;...Total energy consumption Wh/km;...Total energy consumption mpgge;...WTW PEU Wh PE/km;WTW GHG emission g/km;Range km;Mass};udds505[UDDS_505_Net_tractive_energy;UDDS_505_Fuel_energy_WhpKm;...UDDS_505_Total_energy_consumption_WhpKm;...UDDS_505_Total_energy_consumption_mpgge;...UDDS_505_WTW_PEU;UDDS_505_WTW_GHG;...UDDS_505_Range_km;Total_mass];hwfet[HWFET_Net_tractive_energy;HWFET_Fuel_energy_WhpKm;...HWFET_Total_energy_consumption_WhpKm;...HWFET_Total_energy_consumption_mpgge;...HWFET_WTW_PEU;HWFET_WTW_GHG;...HWFET_Range_km;Total_mass];us06city[US06_city_Net_tractive_energy;US06_city_Fuel_energy_WhpKm;...US06_city_Total_energy_consumption_WhpKm;...US06_city_Total_energy_consumption_mpgge;...US06_city_WTW_PEU;US06_city_WTW_GHG;...US06_city_Range_km;Total_mass];us06highway[US06_highway_Net_tractive_energy;US06_highway_Fuel_energy_WhpKm;...US06_highway_Total_energy_consumption_WhpKm;...US06_highway_Total_energy_consumption_mpgge;...US06_highway_WTW_PEU;US06_highway_WTW_GHG;...US06_highway_Range_km;Total_mass];%weighted resultsweighted[(0.29*UDDS_505_Net_tractive_energy0.12*HWFET_Net_tractive_energy...0.14*US06_city_Net_tractive_energy0.45*US06_highway_Net_tractive_energy);...(0.29*UDDS_505_Fuel_energy_WhpKm0.12*HWFET_Fuel_energy_WhpKm...0.14*US06_city_Fuel_energy_WhpKm0.45*US06_highway_Fuel_energy_WhpKm);...(0.29*UDDS_505_Total_energy_consumption_WhpKm0.12*HWFET_Total_energy_consumption_WhpKm...0.14*US06_city_Total_energy_consumption_WhpKm0.45*US06_highway_Total_energy_consumption_WhpKm);...(0.29*UDDS_505_Total_energy_consumption_mpgge0.12*HWFET_Total_energy_consumption_mpgge...0.14*US06_city_Total_energy_consumption_mpgge0.45*US06_highway_Total_energy_consumption_mpgge);...(0.29*UDDS_505_WTW_PEU0.12*HWFET_WTW_PEU...0.14*US06_city_WTW_PEU0.45*US06_highway_WTW_PEU);...(0.29*UDDS_505_WTW_GHG0.12*HWFET_WTW_GHG...0.14*US06_city_WTW_GHG0.45*US06_highway_WTW_GHG);...(0.29*UDDS_505_Range_km0.12*HWFET_Range_km...0.14*US06_city_Range_km0.45*US06_highway_Range_km);Total_mass];table(Drive_Cycle,Max_Velocity_kph,Avg_Velocity_kph,Cycle_Time_s,Cycle_Distance_km)%printing table in command windowtable(energy_cnsmp,udds505,hwfet,us06city,us06highway,weighted)City_mpg0.33* US06_city_Total_energy_consumption_mpgge0.67*UDDS_505_Total_energy_consumption_mpggeHighway_mpg0.78* US06_highway_Total_energy_consumption_mpgge0.22* HWFET_Total_energy_consumption_mpggeCombined0.43*City_mpg0.57*Highway_mpg%local function for repeated calculationsfunction [Net_tractive_energy,Fuel_energy_WhpKm,...Total_energy_consumption_WhpKm,Total_energy_consumption_mpgge,...WTW_PEU,WTW_GHG,Range_km,Total_mass]energy_cunsumption_results(v)run(Glider_specs.m);run(Transmission_Honda.m);run(Engine_R16A.m);run(Willans_engine_model_Honda);run(Drive_cycles);run(Fuel_properties);run(EC_factors);%conventional vehiclemmm_enginem_fuelm_transmission;%added mass% pm_acc(acc_load/omega_max)*(4*pi)*(1/Vd)*(1/10^5);%accessory loadCY1;for c1:length(Gear_ratio)V_max(c)omega_max*r/Gear_ratio(c);endDzeros(1,length(v));fuel_tank37.85;%10 gallons of E10while fuel_tank0for n2:length(v)if v(n)0if v(n)v(n-1)Fw(n)cr*m*g0.5*rho*CdAf*v(n)^2m*(v(n)-v(n-1));Tw(n)Fw(n)*r;Pw(n)Fw(n)*v(n);omega_w(n)v(n)/r;if ((v(n)0) || v(n)V_max(1))G(n)Gear_ratio(1);elseif ((v(n)V_max(1)) v(n)V_max(2))G(n)Gear_ratio(2);elseif ((v(n)V_max(2)) v(n)V_max(3))G(n)Gear_ratio(3);elseif ((v(n)V_max(3)) v(n)V_max(4))G(n)Gear_ratio(4);elseif ((v(n)V_max(4)) v(n)V_max(5))G(n)Gear_ratio(5);elseif ((v(n)V_max(5)) v(n)V_max(6))G(n)Gear_ratio(6);endPe(n)Pw(n)/eff_dr;%with power loss consideration, can account for other power cunsumption sinks like auxiliariesPe(n)Pe(n)acc_load;omega_e(n)G(n)*omega_w(n);Te(n)Pe(n)/omega_e(n);%Te(n)Tw(n)/G(n);pme(n)(4*pi*Te(n))/(Vd);%pascalpme(n)pme(n);cm(n)S*omega_e(n)/pi;e(n)e00e01*cm(n)e02*cm(n)^2;ploss(n)ploss0ploss1*cm(n)ploss2*cm(n)^2;%barploss(n)ploss(n)*10^5;%pascaleff(n)e(n)/(1(ploss(n)/pme(n)));pma(n)pme(n)/eff(n);%or %pma(n)(pme(n)ploss(n))/e(n);%pascalmf(pma(n)*omega_e(n)*Vd)/(4*pi*Hlv);fuel_tankfuel_tank-mf;elseif v(n)v(n-1)Ft_b(n)abs(cr*m*g0.5*rho*CdAf*v(n)^2m*(v(n)-v(n-1)));Pt_b(n)Ft_b(n)*v(n);clutch(n)1;endelseT_idle(n)1;endD(n)D(n-1)(v(n)/1000);if fuel_tank0break;endendD_EOC(CY)D(n);Pw_EOC(CY)trapz(Pw)/3600;CYCY1;endrange_milessum(D_EOC)*0.621371;%milesmiles_per_gallonrange_miles/10;%mpgmpggefuel_tank_liter10*3.78541;%fuel amount in litersfuel_tank_energy9500*fuel_tank_liter; %WhFuel_energy_WhpMilefuel_tank_energy/range_miles;%Wh/mileFuel_energy_kWhpMileFuel_energy_WhpMile/1000;%kWh/mile%tableRange_kmsum(D_EOC);%kmNet_tractive_energytrapz(Pw_EOC)/Range_km;Fuel_energy_WhpKmfuel_tank_energy/Range_km;%Wh/kmTotal_energy_consumption_WhpKmFuel_energy_WhpKm;%Wh/kmTotal_energy_consumption_mpggemiles_per_gallon;%mpggeWTW_PEUFuel_energy_kWhpMile*PEU_WTW_factor_E10;%Wh PE/mileWTW_GHGFuel_energy_kWhpMile*GHG_WTW_factor_E10;%g GHG/mileTotal_massm;end 参考文献 往期回顾可以关注主页点击搜索