CN104112036A - 混联式混合动力电动汽车的仿真方法 - Google Patents

混联式混合动力电动汽车的仿真方法 Download PDF

Info

Publication number
CN104112036A
CN104112036A CN201410259855.6A CN201410259855A CN104112036A CN 104112036 A CN104112036 A CN 104112036A CN 201410259855 A CN201410259855 A CN 201410259855A CN 104112036 A CN104112036 A CN 104112036A
Authority
CN
China
Prior art keywords
motor
energy
model
power
parallel hybrid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410259855.6A
Other languages
English (en)
Other versions
CN104112036B (zh
Inventor
杜云峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan University of Arts and Science
Original Assignee
Hunan University of Arts and Science
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hunan University of Arts and Science filed Critical Hunan University of Arts and Science
Priority to CN201410259855.6A priority Critical patent/CN104112036B/zh
Publication of CN104112036A publication Critical patent/CN104112036A/zh
Application granted granted Critical
Publication of CN104112036B publication Critical patent/CN104112036B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

一种混联式混合动力电动汽车系统的仿真方法,包括首先建立由发动机、发电机、电动机、储能装置(蓄电池)、动力分配装置、制动能量回收系统、能量管理系统组成的数学仿真模型,再对该系统进行系统性能分析和优化设计,最后得出仿真系统的各项性能指标;其仿真方法包括如下步骤:1.确定系统组成,2.建立系统模型,3.进行系统仿真;该方法可以实现能量利用率的最优管理,使有限的燃油发挥最大的功效;使混联式混合动力电动汽车可以在再生制动状态下将制动能量回馈给电池,实现电动汽车能量的回收,从而大量减少混联式混合动力电动汽车的油耗和排放,达到节能和环保的目的,还可对能量进行最优利用,迅速、平稳地进行加、减速,具有良好的动力性能。

Description

混联式混合动力电动汽车的仿真方法
技术领域
本发明涉及一种混联式混合动力电动汽车系统的仿真方法,属于能源技术领域。 
背景技术
现代水力发电、核能发电、风力发电和太阳能的利用,为人们提供了巨大的能源。各种高能蓄电池和高效率电机不断地出现,使人们把目光转向了电动汽车,电动车辆可以摆脱对石油燃料的依赖和实现“零污染”或“超低污染”的排放,为电动汽车的“东山再起”创造了有利条件。因此,各国和各大汽车公司都重视了电动汽车的研究、开发和试制; 
混联式混合动力电动汽车的主要特征是同时具有串联式和并联式两种驱动方式。混联式混合动力电动汽车是在综合了串联式混合动力电动汽车和并联式混合动力电动汽车的结构特点而组成的,既可以在串联混合模式下工作,也可以在并联混合模式下工作。它与串联式相比,主要增加了机械动力传递路线,而与并联式相比它增加了电力驱动传输路线。即同时兼顾了串联式和并联式的优点,结构上保证了在各种复杂工况下可灵活采取串联方式或并联方式,以实现热效率最高、污染量排放最低的目标。 
发明内容
本发明的目的是提供一种混联式混合动力电动汽车系统的仿真方法,用以实现混联式混合动力电动汽车的运行,使电动车辆可以摆脱对石油燃料的依赖,实现“零污染”或“超低污染”的排放,为电动汽车的“东山再起”创造有利条件;使混联式混合动力电动汽车可以实现能量的最优利用,从而实现低油耗、低排放和节能的目标。 
为实现发明目的,本发明采用的技术方案是:根据混联式混合动力电动汽车系统的工作原理,建立系统的仿真模型,运用混联式驱动系统和能量管理系统的控制策略,对系统进行仿真,得到电动汽车行驶速度图、发动机输出功率图、发电机转速图、电池荷电状态图和电动机转速图。为了进行系统性能分析 和优化设计,建立由发动机、发电机、电动机、储能装置(蓄电池)、动力分配装置、制动能量回收系统、能量管理系统等组成的仿真模型,并在MATLAB/SIMULINK平台上对此模型进行仿真,其仿真方法包括如下步骤:1.确定系统组成:2.建立系统模型,3.进行系统仿真。 
其中1:所述确定系统组成,是确定本系统由发动机、发电机、电动机、储能装置即蓄电池、动力分配装置、制动能量回收系统、能量管理系统组成。 
其中2:所述建立系统模型是根据系统工作原理,建立混联式混合动力电动汽车系统各个组成部分的模型。以下主要介绍发动机模型、发电机模型、电动机模型、储能装置即蓄电池模型、动力分配装置模型、制动能量回收系统模型、能量管理系统模型: 
其中:A.发动机模型: 
考虑发动机单独驱动的情况,汽车行驶时的功率平衡方程式为 
P e = 1 η T ( Gfu a 3600 + Giu a 3600 + C D Au a 3 76140 + δ mu a 3600 du dt ) - - - ( 1 )
式中,Pe为发动机输出功率;ηT为传动系效率;G为汽车重量;f为滚动阻力系数;i为汽车行驶路面的坡度;ua为汽车行驶速度;CD为空气阻力系数;A为迎风面积;δ为旋转质量换算系数;m为汽车质量; 
B.发电机和电动机模型: 
发电机和电动机均为永磁同步电机,两个模块内部模型相同,只是在能量管理系统的控制下其运行方式不同。永磁同步电机转子是永久磁铁,磁通有正弦波和梯形波两种模式,本文采用正弦波控制永磁同步电机。正弦波控制永磁同步电机模型假定定子磁通是正弦分布的,因此产生的感应电动势也是正弦的,因为一般永磁同步电机气隙比较大,因此模型没有考虑定子磁路的饱和和铁损。永磁同步电机模型是建立在二相旋转坐标系(dq坐标系)上的,所有参数也折算到旋转坐标系上。电机方程如下: 
d dt i d = 1 L d u d - R L d i d + L q L d p ω r i q - - - ( 2 )
d dt i q = 1 L q u q - R L q i q - L d L q p ω r i d - λp ω r L q - - - ( 3 )
Te=1.5p[λiq+(Ld-Lq)idiq]     (4) 
式中,Lq、Ld为q、d轴电感;R为定子绕组电阻;iq、id为q、d轴电流分量;uq、ud为q、d轴电压分量;ωr为转子角速度;p为极对数;Te为电磁转矩;λ为永磁转子对定子的磁通系数。 
永磁同步电机的机械方程: 
d dt ω r = 1 J ( T e - F ω r - T m ) - - - ( 5 )
dθ dt = ω r - - - ( 6 )
式中,J为转子和负载的转动惯量;F为转子和负载的摩擦系数;θ为转子位置角;Tm为轴上机械转矩; 
C.储能装置即蓄电池模型: 
本系统的蓄电池采用KiBaM模型,该模型可以实时反映蓄电池容量与充、放电电流的关系。设在仿真步长Δt内,系统额定电压为U,负载需求功率为Pd,系统剩余功率为Pe,则有: 
充电电流:Ic=Pe/U 
最大充电电流Icmax: 
Icmax=[-kcqmax+kq10e-kΔt+q0kc(1-e-kΔt)]/[1-e-kΔt+c(kΔt-1+e-kΔt)]     (7) 
放电电流:Id=Pd/U 
最大放电电流Idmax: 
Idmax=[kq10e-kΔt+q0kc(1-e-kΔt)]/[1-e-kΔt+c(kΔt-1+e-kΔt)]     (8) 
式中,c表示可用电荷容量与总容量的比值;q10表示在时间步长Δt起始时刻可用电荷容量;q0表示在时间步长Δt起始时刻总电荷容量;k表示比率常数;qmax表示最大容量; 
D.动力分配装置模型: 
在混联式混合动力传动结构中,用于实现能量分流和综合的动力分配装置是一个行星齿轮机构,其中行星架与发动机的输出轴相连,齿圈与电动机的转轴相连,同时也与输出齿轮相连,而太阳齿轮轴发出的动力驱动发电机发电,中间与一离合器相连,必要时锁死太阳轮,使行星齿轮机构以一定的传动比工作。发电机的转子刚性联接在发动机的输出轴上; 
通过对行星机构的变速比和受力分析,可以得到如下方程组: 
ω 1 + k p ω 2 - ( 1 + k p ) ω 3 = 0 T 3 = ( 1 + k p ) η S T 1 = ( 1 + k p ) η R k p T 2 - - - ( 9 )
式中,kp为齿数比,kp=z2/z1;z1为太阳轮齿数;z2为齿圈齿数;ηS为由太阳轮到行星架的效率;ηR为由齿圈到行星架的效率;T1、T2、T3分别为太阳轮、齿圈和行星架的转矩;ω1、ω2、ω3分别为太阳轮、齿圈和行星架的角速度; 
作用在驱动轮上的转矩Tt是由发动机产生的转矩经过动力分配装置后传至车轮上的,因此驱动力为: 
F t = T tq i 0 η T r - - - ( 10 )
式中,Ttq为作用在齿圈上的总转矩;i0为主减速器速比;ηT为传动系效率;r为车轮半径; 
E.制动能量回收系统模型: 
制动能量回收系统是混合动力电动汽车的重要组成部分之一,又称再生制动系统,是指汽车在制动或下坡时将储存于车身上的势能和动能,通过电动机转化为电能,并储存于储能装置中的过程; 
再生制动系统的结构与原理如图1所示,由驱动轮、主减速器、变速器、电动机、AC/DC转换器、DC/DC转换器、能量储存系统以及控制器组成; 
汽车在制动或滑行过程中,根据驾驶员的制动意图,由制动控制器计算得到汽车需要的总制动力。再根据一定的制动力分配控制策略得到电动机应该提供的电动机再生制动力。电动机控制器计算需要的电动机电枢中的制动电流,通过一定的控制方法使电动机跟踪需要的制动电流,从而较准确地提供再生制动力矩,在电动机的电枢中产生的电流经AC/DC转换器整流再经DC/DC转换器反充到储能装置中保存起来; 
F.能量管理系统模型: 
能量管理系统随时向驾驶员提供有用的信息,如充电时,显示充电状态;运行时,显示电池的电压、电流、温度、电解液的密度、剩余电量以及车速等,使驾驶员心中有数,操作正确。上述信息也随时提供给计算机,用来分析诊断故障、保持工作正常或保持最优运行状态,如通过电动机控制器避免电池过度放电,通过平衡单元保持各电池的温度均衡和电量均衡。 
其中3:所述进行系统仿真是在MATLAB/SIMULINK中,对混联式混合动力电动汽车系统进行建模仿真。电动汽车的整车主要参数为:汽车质量约1325kg,满载1600kg;旋转质量换算系数为1.2;迎风面积为2.57m2;空气阻力系数为0.3;滚动阻力系数为0.018;车轮半径为0.3m; 
根据系统工作原理,建立混联式混合动力电动汽车系统的各个子模块。以下主要介绍发动机、发电机、电动机、蓄电池、动力分配装置模块的主要参数: 
(1)发动机模块: 
发动机模块的主要参数为:最大功率为57kW,最大转速为6000rpm,排量为1.5L; 
(2)发电机模块: 
发电机模块的主要参数为:最大功率为30kW,最大转速为13000rpm,定子绕组电阻为0.05Ω,q、d轴电感均为0.000635H,极对数为1; 
(3)电动机模块: 
电动机模块的主要参数为:最大功率为50kW,最大转速为6000rpm,定子绕组电阻为0.0065Ω,q、d轴电感分别为0.002057H和0.001598H,极对数为4; 
(4)蓄电池模块: 
蓄电池模块的主要参数为:最大功率为21kW,额定电压为200V,额定容量为6.5Ah,初始荷电状态为41.53%,满充电电压为235.6V,额定放电电流为1.3A,内阻为0.30769Ω; 
(5)动力分配装置模块: 
动力分配模块的主要参数为:太阳轮齿数为29,齿圈齿数为75,主减速器速比为4.133。 
将系统的各个子模块组合在一起,建立混联式混合动力电动汽车的整体模型,对系统进行仿真,仿真时间为8s。图2、3、4、5、6分别为电动汽车行驶速度图、发动机输出功率图、发电机转速图、电池荷电状态图和电动机转速图。 
在0-1.1s期间,油门踏板状态为80%,对电动汽车进行加速。此时混联式混合动力电动汽车所需功率较小,混联式混合动力电动汽车只由电动机驱动,所需能量由电池提供,混联式混合动力电动汽车以纯电动工作模式运行。在图2中,混联式混合动力电动汽车的速度从0km/h变为12km/h。在图3中,发动机输出功率为0kW。发电机转速如图4所示。在图5中,电池荷电状态从41.53%变为41.36%,电池处于放电状态。在图6中,电动机转速从0rpm变为400rpm。 
在1.1-2s期间,油门踏板状态保持80%不变,对电动汽车继续进行加速。此时混联式混合动力电动汽车所需功率较大,能量管理系统启动混合运行模式, 所需能量由发动机和电池共同提供。在图2中,混联式混合动力电动汽车的速度从12km/h变为22km/h。在图3中,发动机输出功率首先从0kW变为16kW,然后再从16kW变为27kW。发电机转速如图4所示。在图5中,电池荷电状态从41.36%变为41.34%,电池处于放电状态。在图6中,电动机转速从400rpm变为800rpm。 
在2-4s期间,油门踏板状态为20%,对电动汽车继续进行加速。在图2中,混联式混合动力电动汽车的速度从22km/h变为27.5km/h。当t=2s时,发动机不能立刻减少输出功率,在能量管理系统的控制下,电池停止向电动机供电,而是开始由发电机对其进行充电。当t=2.2s时,发电机停止工作,所需能量只由电池提供。发动机输出功率如图3所示。发电机转速如图4所示。电池荷电状态如图5所示,在2-2.2s期间,电池处于充电状态,在2.2-4s期间,电池处于放电状态。在图6中,电动机转速从800rpm变为1000rpm。 
在4-6s期间,油门踏板状态为75%,对电动汽车继续进行加速。在图2中,混联式混合动力电动汽车的速度从27.5km/h变为47.5km/h。当t=4.2s时,发动机重新启动,向混联式混合动力电动汽车提供能量。发动机输出功率如图3所示。发电机转速如图4所示。在图5中,电池荷电状态从41.09%变为40.78%,电池处于放电状态。在图6中,电动机转速从1000rpm变为1720rpm。 
在6-8s期间,油门踏板状态为-60%(即刹车踏板状态为60%),对电动汽车进行减速(即进行刹车)。在图2中,混联式混合动力电动汽车的速度首先从47.5km/h变为48km/h,然后再从48km/h变为40km/h。混联式混合动力电动汽车进入再生制动状态,发电机停止工作,而电动机开始作发电机运行,将制动能量回馈给电池,电池处于充电状态,将能量储存在电池中,从而实现电动汽车能量的回收。发动机输出功率如图3所示。发电机转速如图4所示。在图5中,电池荷电状态从40.78%变为41.45%。在图6中,电动机转速首先从1720rpm变为1740rpm,然后再从1740rpm变为1450rpm。 
与现有技术相比,本发明有如下特点和进步:由于本混联式混合动力电动汽车系统的仿真方法,建立了由发动机、发电机、电动机、储能装置(蓄电池)、动力分配装置、制动能量回收系统、能量管理系统等组成的混联式混合动力电动汽车系统的数学仿真模型,然后在MATLAB软件中应用SIMULINK对系统进行了建模仿真。因此该方法可以实现能量利用率的最优管理,使有限的燃油发挥最大的功效。混联式混合动力电动汽车可以在再生制动状态下将制动能量回馈给电池,实现电动汽车能量的回收,从而大量减少混联式混合动力电动汽车的油耗和排放,达到节能和环保的目的。另外,在能量管理系统的控制下,电动机辅助发动机工作,混联式混合动力电动汽车可以迅速平稳地响应加、减速的 要求。能量管理系统根据电池的荷电状态,对电池进行充、放电等管理,保持电池的电量平衡。混联式混合动力电动汽车能够对能量进行最优利用,可以迅速、平稳地进行加、减速,具有良好的动力性能。混联式混合动力电动汽车在再生制动状态下,能够回收电动汽车的制动能量,可以获得较高的燃油经济性,使发动机的排污大为降低。 
附图说明
下面结合附图和实施例对本本发明进一步说明: 
图1是本发明的再生制动系统结构与原理图; 
图2是本发明的电动汽车行驶速度图; 
图3是本发明系统运行时的发动机输出功率图; 
图4是本发明系统运行时的发电机转速图; 
图5是本发明系统运行时的电池荷电状态图; 
图6是本发明系统运行时的电动机转速图。 
具体实施方式
参照附图:设计出的混联式混合动力电动汽车系统的仿真方法包括如下步骤:1.确定系统组成:所述混联式混合动力电动汽车系统,由发动机、发电机、电动机、储能装置(蓄电池)、动力分配装置、制动能量回收系统、能量管理系统等组成;发动机输出的一部分功率通过机械传动装置直接输送到驱动桥,另一部分则驱动发电机发电。发电机发出的电能输送给蓄电池或电动机,电动机产生的驱动力矩通过动力复合装置传送给驱动桥。混联式驱动系统的控制策略是:当汽车低速行驶时,驱动系统主要以串联方式工作;当汽车高速稳定行驶时,驱动系统则以并联工作方式为主。混联式混合动力结构采用行星齿轮机构作为动力分配装置。将发动机、发电机和电动机通过一个行星齿轮装置连接起来,动力从发动机输出到与其相连的行星架,行星架将一部分转矩传送到发电机,另一部分传送到传动轴,同时发电机也可以驱动电动机来驱动传动轴。这种机构有两个自由度,可以自由地控制两个不同的速度。此时车辆并不是串联式或并联式,而是两种驱动形式同时存在。这种形式充分利用了两种驱动形式 的优点。能量管理系统根据汽车运行工况,控制各部分协调工作。能量管理系统在满足汽车动力性能的前提下,最优分配电动机、发动机、电池等动力部件的功率输出,从而实现能量的最优分配,减少汽车的油耗和排放。另外,能量管理系统还要考虑电池的荷电状态平衡,对电池进行充、放电等管理,以延长电池的寿命,降低汽车的维护成本; 
2.建立系统模型:根据系统工作原理,建立混联式混合动力电动汽车系统各个组成部分的模型。以下主要介绍发动机模型、发电机模型、电动机模型、储能装置(蓄电池)模型、动力分配装置模型、制动能量回收系统模型、能量管理系统模型: 
其中:A.发动机: 
考虑发动机单独驱动的情况,汽车行驶时的功率平衡方程式为 
P e = 1 η T ( Gfu a 3600 + Giu a 3600 + C D Au a 3 76140 + δ mu a 3600 du dt ) - - - ( 1 )
式中,Pe为发动机输出功率;ηT为传动系效率;G为汽车重量;f为滚动阻力系数;i为汽车行驶路面的坡度;ua为汽车行驶速度;CD为空气阻力系数;A为迎风面积;δ为旋转质量换算系数;m为汽车质量。 
B.发电机和电动机: 
发电机和电动机均为永磁同步电机,两个模块内部模型相同,只是在能量管理系统的控制下其运行方式不同。永磁同步电机转子是永久磁铁,磁通有正弦波和梯形波两种模式,本文采用正弦波控制永磁同步电机。正弦波控制永磁同步电机模型假定定子磁通是正弦分布的,因此产生的感应电动势也是正弦的,因为一般永磁同步电机气隙比较大,因此模型没有考虑定子磁路的饱和和铁损。永磁同步电机模型是建立在二相旋转坐标系(dq坐标系)上的,所有参数也折算到旋转坐标系上。电机方程如下: 
d dt i d = 1 L d u d - R L d i d + L q L d p ω r i q - - - ( 2 )
d dt i q = 1 L q u q - R L q i q - L d L q p ω r i d - λp ω r L q - - - ( 3 )
Te=1.5p[λiq+(Ld-Lq)idiq]     (4) 
式中,Lq、Ld为q、d轴电感;R为定子绕组电阻;iq、id为q、d轴电流分量;uq、ud为q、d轴电压分量;ωr为转子角速度;p为极对数;Te为电磁转矩;λ为永磁转子对定子的磁通系数。 
永磁同步电机的机械方程: 
d dt ω r = 1 J ( T e - F ω r - T m ) - - - ( 5 )
dθ dt = ω r - - - ( 6 )
式中,J为转子和负载的转动惯量;F为转子和负载的摩擦系数;θ为转子位置角;Tm为轴上机械转矩。 
C.储能装置(蓄电池): 
本文的蓄电池采用KiBaM模型,该模型可以实时反映蓄电池容量与充、放电电流的关系。设在仿真步长Δt内,系统额定电压为U,负载需求功率为Pd,系统剩余功率为Pe,则有: 
充电电流:Ic=Pe/U 
最大充电电流Icmax: 
Icmax=[-kcqmax+kq10e-kΔt+q0kc(1-e-kΔt)]/[1-e-kΔt+c(kΔt-1+e-kΔt)]     (7) 
放电电流:id=Pd/U 
最大放电电流Idmax: 
Idmax=[kq10e-kΔt+q0kc(1-e-kΔt)]/[1-e-kΔt+c(kΔt-1+e-kΔt)]     (8) 
式中,c表示可用电荷容量与总容量的比值;q10表示在时间步长Δt起始时刻可用电荷容量;q0表示在时间步长Δt起始时刻总电荷容量;k表示比率常数;qmax表示最大容量。 
D.动力分配装置: 
在混联式混合动力传动结构中,用于实现能量分流和综合的动力分配装置是一个行星齿轮机构,其中行星架与发动机的输出轴相连,齿圈与电动机的转轴相连,同时也与输出齿轮相连,而太阳齿轮轴发出的动力驱动发电机发电,中间与一离合器相连,必要时锁死太阳轮,使行星齿轮机构以一定的传动比工作。发电机的转子刚性联接在发动机的输出轴上; 
通过对行星机构的变速比和受力分析,可以得到如下方程组: 
ω 1 + k p ω 2 - ( 1 + k p ) ω 3 = 0 T 3 = ( 1 + k p ) η S T 1 = ( 1 + k p ) η R k p T 2 - - - ( 9 )
式中,kp为齿数比,kp=z2/z1;z1为太阳轮齿数;z2为齿圈齿数;ηS为由太阳轮到行星架的效率;ηR为由齿圈到行星架的效率;T1、T2、T3分别为太阳轮、齿圈和行星架的转矩;ω1、ω2、ω3分别为太阳轮、齿圈和行星架的角速度; 
作用在驱动轮上的转矩Tt是由发动机产生的转矩经过动力分配装置后传至车轮上的,因此驱动力为: 
F t = T tq i 0 η T r - - - ( 10 )
式中,Ttq为作用在齿圈上的总转矩;i0为主减速器速比;ηT为传动系效率;r为车轮半径。 
E.制动能量回收系统: 
制动能量回收系统是混合动力电动汽车的重要组成部分之一,又称再生制动系统,是指汽车在制动或下坡时将储存于车身上的势能和动能,通过电动机转化为电能,并储存于储能装置中的过程; 
再生制动系统的结构与原理如图1所示,由驱动轮、主减速器、变速器、电动机、AC/DC转换器、DC/DC转换器、能量储存系统以及控制器组成; 
汽车在制动或滑行过程中,根据驾驶员的制动意图,由制动控制器计算得到汽车需要的总制动力。再根据一定的制动力分配控制策略得到电动机应该提供的电动机再生制动力。电动机控制器计算需要的电动机电枢中的制动电流,通过一定的控制方法使电动机跟踪需要的制动电流,从而较准确地提供再生制动力矩,在电动机的电枢中产生的电流经AC/DC转换器整流再经DC/DC转换器反充到储能装置中保存起来。 
F.能量管理系统: 
能量管理系统随时向驾驶员提供有用的信息,如充电时,显示充电状态;运行时,显示电池的电压、电流、温度、电解液的密度、剩余电量以及车速等,使驾驶员心中有数,操作正确。上述信息也随时提供给计算机,用来分析诊断故障、保持工作正常或保持最优运行状态,如通过电动机控制器避免电池过度放电,通过平衡单元保持各电池的温度均衡和电量均衡等。 
3.进行系统仿真: 
在MATLAB/SIMULINK中,对混联式混合动力电动汽车系统进行建模仿真。电动汽车的整车主要参数为:汽车质量约1325kg,满载1600kg;旋转质量换算系数为1.2;迎风面积为2.57m2;空气阻力系数为0.3;滚动阻力系数为0.018;车轮半径为0.3m。 
根据系统工作原理,建立混联式混合动力电动汽车系统的各个子模块。以下主要介绍发动机、发电机、电动机、蓄电池、动力分配装置模块的主要参数: 
(1)发动机模块: 
发动机模块的主要参数为:最大功率为57kW,最大转速为6000rpm,排量为1.5L; 
(2)发电机模块: 
发电机模块的主要参数为:最大功率为30kW,最大转速为13000rpm,定子绕组电阻为0.05Ω,q、d轴电感均为0.000635H,极对数为1; 
(3)电动机模块: 
电动机模块的主要参数为:最大功率为50kW,最大转速为6000rpm,定子绕组电阻为0.0065Ω,q、d轴电感分别为0.002057H和0.001598H,极对数为4; 
(4)蓄电池模块: 
蓄电池模块的主要参数为:最大功率为21kW,额定电压为200V,额定容量为6.5Ah,初始荷电状态为41.53%,满充电电压为235.6V,额定放电电流为1.3A,内阻为0.30769Ω; 
(5)动力分配装置模块: 
动力分配模块的主要参数为:太阳轮齿数为29,齿圈齿数为75,主减速器速比为4.133。 
将系统的各个子模块组合在一起,建立混联式混合动力电动汽车的整体模型,对系统进行仿真,仿真时间为8s。图2、3、4、5、6分别为电动汽车行驶速度图、发动机输出功率图、发电机转速图、电池荷电状态图和电动机转速图。 
在0-1.1s期间,油门踏板状态为80%,对电动汽车进行加速。此时混联式混合动力电动汽车所需功率较小,混联式混合动力电动汽车只由电动机驱动,所需能量由电池提供,混联式混合动力电动汽车以纯电动工作模式运行。在图2中,混联式混合动力电动汽车的速度从0km/h变为12km/h。在图3中,发动机输出功率为0kW。发电机转速如图4所示。在图5中,电池荷电状态从41.53%变为41.36%,电池处于放电状态。在图6中,电动机转速从0rpm变为400rpm。 
在1.1-2s期间,油门踏板状态保持80%不变,对电动汽车继续进行加速。此时混联式混合动力电动汽车所需功率较大,能量管理系统启动混合运行模式,所需能量由发动机和电池共同提供。在图2中,混联式混合动力电动汽车的速度从12km/h变为22km/h。在图3中,发动机输出功率首先从0kW变为16kW,然后再从16kW变为27kW。发电机转速如图4所示。在图5中,电池荷电状态从41.36%变为41.34%,电池处于放电状态。在图6中,电动机转速从400rpm变为800rpm。 
在2-4s期间,油门踏板状态为20%,对电动汽车继续进行加速。在图2中,混联式混合动力电动汽车的速度从22km/h变为27.5km/h。当t=2s时,发动机不能立刻减少输出功率,在能量管理系统的控制下,电池停止向电动机供电,而是开始由发电机对其进行充电。当t=2.2s时,发电机停止工作,所需能量只由电池提供。发动机输出功率如图3所示。发电机转速如图4所示。电池荷电状态如图5所示,在2-2.2s期间,电池处于充电状态,在2.2-4s期间,电池处于放电状态。在图6中,电动机转速从800rpm变为1000rpm。 
在4-6s期间,油门踏板状态为75%,对电动汽车继续进行加速。在图2中,混联式混合动力电动汽车的速度从27.5km/h变为47.5km/h。当t=4.2s时,发动机重新启动,向混联式混合动力电动汽车提供能量。发动机输出功率如图3所示。发电机转速如图4所示。在图5中,电池荷电状态从41.09%变为40.78%,电池处于放电状态。在图6中,电动机转速从1000rpm变为1720rpm。 
在6-8s期间,油门踏板状态为-60%(即刹车踏板状态为60%),对电动汽车进行减速(即进行刹车)。在图2中,混联式混合动力电动汽车的速度首先从47.5km/h变为48km/h,然后再从48km/h变为40km/h。混联式混合动力电动汽车进入再生制动状态,发电机停止工作,而电动机开始作发电机运行,将制动能量回馈给电池,电池处于充电状态,将能量储存在电池中,从而实现电动汽车能量的回收。发动机输出功率如图3所示。发电机转速如图4所示。在图5中,电池荷电状态从40.78%变为41.45%。在图6中,电动机转速首先从1720rpm变为1740rpm,然后再从1740rpm变为1450rpm。 

Claims (3)

1.一种混联式混合动力电动汽车系统的仿真方法,其特征是:根据混联式混合动力电动汽车系统工作原理,先确定系统组成,再载建立系统的仿真模型,然后运用混联式驱动系统和能量管理系统的控制策略,对系统进行仿真,得到电动汽车行驶速度图、发动机输出功率图、发电机转速图、电池荷电状态图和电动机转速图;为了进行系统性能分析和优化设计,建立由发动机、发电机、电动机、储能装置(蓄电池)、动力分配装置、制动能量回收系统、能量管理系统组成的仿真模型,并在MATLAB/SIMULINK平台上对此模型进行仿真,其仿真方法包括如下步骤:1.确定系统组成:2.建立系统模型,3.进行系统仿真。 
2.根据权利要求1所述混联式混合动力电动汽车系统的仿真方法,其特征是:所述建立系统模型是根据系统工作原理,建立混联式混合动力电动汽车系统各个组成部分的模型,主要是建立发动机模型、发电机模型、电动机模型、储能装置即蓄电池模型、动力分配装置模型、制动能量回收系统模型和能量管理系统模型: 
其中:A.发动机模型: 
考虑发动机单独驱动的情况,汽车行驶时的功率平衡方程式为: 
式中,Pe为发动机输出功率;ηT为传动系效率;G为汽车重量;f为滚动阻力系数;i为汽车行驶路面的坡度;ua为汽车行驶速度;CD为空气阻力系数;A为迎风面积;δ为旋转质量换算系数;m为汽车质量; 
B.发电机和电动机模型: 
发电机和电动机均为永磁同步电机,两个模块内部模型相同,只是在能量管理系统的控制下其运行方式不同;永磁同步电机转子是永久磁铁,磁通有正弦波和梯形波两种模式,本文采用正弦波控制永磁同步电机;正弦波控制永磁同步电机模型假定定子磁通是正弦分布的,因此产生的感应电动势也是正弦的,因为一般永磁同步电机气隙比较大,因此模型没有考虑定子磁路的饱和和铁损;永磁同步电机模型是建立在二相旋转坐标系(dq坐标系)上的,所有参数也折算到旋转坐标系上;电机方程如下: 
Te=1.5p[λiq+(Ld-Lq)idiq]     (4) 
式中,Lq、Ld为q、d轴电感;R为定子绕组电阻;iq、id为q、d轴电流分量;uq、ud为q、d轴电压分量;ωr为转子角速度;p为极对数;Te为电磁转矩;λ为永磁转子对定子的磁通系数; 
永磁同步电机的机械方程: 
式中,J为转子和负载的转动惯量;F为转子和负载的摩擦系数;θ为转子位置角;Tm为轴上机械转矩; 
C.储能装置即蓄电池模型: 
本系统的蓄电池采用KiBaM模型,该模型可以实时反映蓄电池容量与充、放电电流的关系;设在仿真步长Δt内,系统额定电压为U,负载需求功率为Pd,系统剩余功率为Pe,则有: 
充电电流:Ic=Pe/U, 
最大充电电流Icmax: 
Icmax=[-kcqmax+Kq10e-kΔt+q0kc(1-e-kΔt)]/[1-e-kΔt+c(kΔt-1+e-kΔt)]     (7) 
放电电流:Id=Pd/U, 
最大放电电流Idmax: 
Idmax=[kq10e-kΔt+q0kc(1-e-kΔt)]/[1-e-kΔt+c(kΔt-1+e-kΔt)]     (8) 
式中,c表示可用电荷容量与总容量的比值;q10表示在时间步长Δt起始时刻可用电荷容量;q0表示在时间步长Δt起始时刻总电荷容量;k表示比率常数;qmax表示最大容量; 
D.动力分配装置模型: 
在混联式混合动力传动结构中,用于实现能量分流和综合的动力分配装置是一个行星齿轮机构,其中行星架与发动机的输出轴相连,齿圈与电动机的转轴相连,同时也与输出齿轮相连,而太阳齿轮轴发出的动力驱动发电机发电,中间与一离合器相连,必要时锁死太阳轮,使行星齿轮机构以一定的传动比工作;发电机的转子刚性联接在发动机的输出轴上; 
通过对行星机构的变速比和受力分析,可以得到如下方程组: 
式中,kp为齿数比,kp=z2/z1;z1为太阳轮齿数;z2为齿圈齿数;ηS为由太阳轮到行星架的效率;ηR为由齿圈到行星架的效率;T1、T2、T3分别为太阳轮、齿圈和行星架的转矩;ω1、ω2、ω3分别为太阳轮、齿圈和行星架的角速度; 
作用在驱动轮上的转矩Tt是由发动机产生的转矩经过动力分配装置后传至车轮上的,因此驱动力为: 
式中,Ttq为作用在齿圈上的总转矩;i0为主减速器速比;ηT为传动系效率;r为车轮半径; 
E.制动能量回收系统: 
制动能量回收系统是混合动力电动汽车的重要组成部分之一,又称再生制动系统,是指汽车在制动或下坡时将储存于车身上的势能和动能,通过电动机转化为电能,并储存于储能装置中的过程; 
再生制动系统的结构与原理如图1所示,由驱动轮、主减速器、变速器、电动机、AC/DC转换器、DC/DC转换器、能量储存系统以及控制器组成; 
汽车在制动或滑行过程中,根据驾驶员的制动意图,由制动控制器计算得到汽车需要的总制动力。再根据一定的制动力分配控制策略得到电动机应该提供的电动机再生制动力。电动机控制器计算需要的电动机电枢中的制动电流,通过一定的控制方法使电动机跟踪需要的制动电流,从而较准确地提供再生制动力矩,在电动机的电枢中产生的电流经AC/DC转换器整流再经DC/DC转换器反充到储能装置中保存起来; 
F.能量管理系统模块: 
能量管理系统随时向驾驶员提供有用的信息,如充电时,显示充电状态;运行时,显示电池的电压、电流、温度、电解液的密度、剩余电量以及车速等,使驾驶员心中有数,操作正确。上述信息也随时提供给计算机,用来分析诊断故障、保持工作正常或保持最优运行状 态,如通过电动机控制器避免电池过度放电,通过平衡单元保持各电池的温度均衡和电量均衡。 
3.根据权利要求1所述混联式混合动力电动汽车系统的仿真方法,其特征是:所述进行系统仿真是在MATLAB/SIMULINK中,对混联式混合动力电动汽车系统进行建模仿真;电动汽车的整车主要参数为:汽车质量约1325kg,满载1600kg;旋转质量换算系数为1.2;迎风面积为2.57m2;空气阻力系数为0.3;滚动阻力系数为0.018;车轮半径为0.3m; 
根据系统工作原理,建立混联式混合动力电动汽车系统的各个子模块;以下分别介绍发动机、发电机、电动机、蓄电池、动力分配装置模块的主要参数: 
(1)发动机模块: 
发动机模块的主要参数为:最大功率为57kW,最大转速为6000rpm,排量为1.5L; 
(2)发电机模块: 
发电机模块的主要参数为:最大功率为30kW,最大转速为13000rpm,定子绕组电阻为0.05Ω,q、d轴电感均为0.000635H,极对数为1; 
(3)电动机模块: 
电动机模块的主要参数为:最大功率为50kW,最大转速为6000rpm,定子绕组电阻为0.0065Ω,q、d轴电感分别为0.002057H和0.001598H,极对数为4; 
(4)蓄电池模块: 
蓄电池模块的主要参数为:最大功率为21kW,额定电压为200V,额定容量为6.5Ah,初始荷电状态为41.53%,满充电电压为235.6V,额定放电电流为1.3A,内阻为0.30769Ω; 
(5)动力分配装置模块: 
动力分配模块的主要参数为:太阳轮齿数为29,齿圈齿数为75,主减速器速比为4.133。 
CN201410259855.6A 2014-06-12 2014-06-12 混联式混合动力电动汽车的仿真方法 Expired - Fee Related CN104112036B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410259855.6A CN104112036B (zh) 2014-06-12 2014-06-12 混联式混合动力电动汽车的仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410259855.6A CN104112036B (zh) 2014-06-12 2014-06-12 混联式混合动力电动汽车的仿真方法

Publications (2)

Publication Number Publication Date
CN104112036A true CN104112036A (zh) 2014-10-22
CN104112036B CN104112036B (zh) 2017-03-01

Family

ID=51708824

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410259855.6A Expired - Fee Related CN104112036B (zh) 2014-06-12 2014-06-12 混联式混合动力电动汽车的仿真方法

Country Status (1)

Country Link
CN (1) CN104112036B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104614986A (zh) * 2014-12-08 2015-05-13 天津大学 基于工况自学习的混合动力总成优化标定实验系统及方法
CN104698864A (zh) * 2015-02-26 2015-06-10 安徽江淮汽车股份有限公司 湿式双离合器变速器的离合器模型控制方法及系统
CN105292112A (zh) * 2015-10-29 2016-02-03 清华大学 一种混合动力车辆的结构优化设计方法
CN105446316A (zh) * 2015-11-12 2016-03-30 东方电气(乐山)新能源设备有限公司 一种纯电动汽车整车控制器模型在环测试装置及其测试方法
CN105608273A (zh) * 2015-12-23 2016-05-25 华晨汽车集团控股有限公司 一种基于cfd软件优化动力电池包温度分布的方法
CN105868487A (zh) * 2016-04-11 2016-08-17 中国建筑科学研究院 一种基于能耗设备模型的节能专家系统
CN105893648A (zh) * 2015-01-26 2016-08-24 德尔福派克电气系统有限公司 混合动力和电动汽车高压电气导线线径选择方法
CN106054669A (zh) * 2016-06-22 2016-10-26 北京理工大学 一种单轴并联混合动力汽车仿真平台及仿真方法
CN106407587A (zh) * 2016-09-29 2017-02-15 肇庆学院 混合动力货车防怠速系统的建立方法
CN107144793A (zh) * 2017-05-03 2017-09-08 山东大学 一种动力电池的分数阶KiBaM模型参数辨识方法及系统
CN107491571A (zh) * 2016-08-11 2017-12-19 宝沃汽车(中国)有限公司 一种整车性能仿真的方法及系统
CN107656207A (zh) * 2017-08-17 2018-02-02 同济大学 一种车用燃料电池发动机在动力总成回路的动态测试方法
CN107848423A (zh) * 2015-07-29 2018-03-27 日产自动车株式会社 电动车辆的控制装置以及电动车辆的控制方法
WO2018153116A1 (zh) * 2017-02-21 2018-08-30 山东大学 计及非线性容量特性的分数阶KiBaM电池模型及参数辨识方法
CN108957359A (zh) * 2018-08-15 2018-12-07 合肥工业大学 一种用于整车环境下复合电源系统瞬时效率的测试方法
CN110194179A (zh) * 2019-06-26 2019-09-03 重庆大学 一种串联式混合动力电动汽车动力模式的确定系统
CN112906296A (zh) * 2021-02-02 2021-06-04 武汉理工大学 混合动力汽车全服役期能量优化方法、系统和存储介质
CN116666787A (zh) * 2023-07-28 2023-08-29 中国标准化研究院 一种锂电池智能管理系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101882174A (zh) * 2010-05-31 2010-11-10 北汽福田汽车股份有限公司 用于优化混合动力汽车的动力系统参数的方法
CN102467598A (zh) * 2010-11-10 2012-05-23 台达电子工业股份有限公司 车辆性能的仿真系统
CN102692326B (zh) * 2012-06-06 2015-04-29 奇瑞汽车股份有限公司 一种模拟混合动力汽车性能仿真试验装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101882174A (zh) * 2010-05-31 2010-11-10 北汽福田汽车股份有限公司 用于优化混合动力汽车的动力系统参数的方法
CN102467598A (zh) * 2010-11-10 2012-05-23 台达电子工业股份有限公司 车辆性能的仿真系统
CN102692326B (zh) * 2012-06-06 2015-04-29 奇瑞汽车股份有限公司 一种模拟混合动力汽车性能仿真试验装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
白羽鹤 等: "基于Matlab/Simulink的混联式混合动力电动汽车建模与仿真", 《计算机与现代化》 *
龚贤武 等: "混联型混合动力汽车建模及控制策略研究", 《郑州大学学报(工学版)》 *

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104614986B (zh) * 2014-12-08 2017-04-05 天津大学 基于工况自学习的混合动力总成优化标定实验系统及方法
CN104614986A (zh) * 2014-12-08 2015-05-13 天津大学 基于工况自学习的混合动力总成优化标定实验系统及方法
CN105893648A (zh) * 2015-01-26 2016-08-24 德尔福派克电气系统有限公司 混合动力和电动汽车高压电气导线线径选择方法
CN105893648B (zh) * 2015-01-26 2020-10-16 安波福电气系统有限公司 混合动力和电动汽车高压电气导线线径选择方法
CN104698864B (zh) * 2015-02-26 2017-03-01 安徽江淮汽车集团股份有限公司 湿式双离合器变速器的离合器模型控制方法及系统
CN104698864A (zh) * 2015-02-26 2015-06-10 安徽江淮汽车股份有限公司 湿式双离合器变速器的离合器模型控制方法及系统
CN107848423A (zh) * 2015-07-29 2018-03-27 日产自动车株式会社 电动车辆的控制装置以及电动车辆的控制方法
CN107848423B (zh) * 2015-07-29 2020-08-18 日产自动车株式会社 电动车辆的控制装置以及电动车辆的控制方法
CN105292112A (zh) * 2015-10-29 2016-02-03 清华大学 一种混合动力车辆的结构优化设计方法
CN105446316A (zh) * 2015-11-12 2016-03-30 东方电气(乐山)新能源设备有限公司 一种纯电动汽车整车控制器模型在环测试装置及其测试方法
CN105608273A (zh) * 2015-12-23 2016-05-25 华晨汽车集团控股有限公司 一种基于cfd软件优化动力电池包温度分布的方法
CN105608273B (zh) * 2015-12-23 2018-08-07 华晨汽车集团控股有限公司 一种基于cfd软件优化动力电池包温度分布的系统
CN105868487A (zh) * 2016-04-11 2016-08-17 中国建筑科学研究院 一种基于能耗设备模型的节能专家系统
CN105868487B (zh) * 2016-04-11 2019-04-16 中国建筑科学研究院 一种基于能耗设备模型的节能专家系统
CN106054669A (zh) * 2016-06-22 2016-10-26 北京理工大学 一种单轴并联混合动力汽车仿真平台及仿真方法
CN106054669B (zh) * 2016-06-22 2019-05-24 北京理工大学 一种单轴并联混合动力汽车仿真平台及仿真方法
CN107491571A (zh) * 2016-08-11 2017-12-19 宝沃汽车(中国)有限公司 一种整车性能仿真的方法及系统
CN106407587B (zh) * 2016-09-29 2019-05-07 肇庆学院 混合动力货车防怠速系统的建立方法
CN106407587A (zh) * 2016-09-29 2017-02-15 肇庆学院 混合动力货车防怠速系统的建立方法
WO2018153116A1 (zh) * 2017-02-21 2018-08-30 山东大学 计及非线性容量特性的分数阶KiBaM电池模型及参数辨识方法
US11526639B2 (en) 2017-02-21 2022-12-13 Shandong University Fractional-order KiBaM battery model considering nonlinear capacity characteristics and parameter identification method
CN107144793A (zh) * 2017-05-03 2017-09-08 山东大学 一种动力电池的分数阶KiBaM模型参数辨识方法及系统
CN107144793B (zh) * 2017-05-03 2019-05-21 山东大学 一种动力电池的分数阶KiBaM模型参数辨识方法及系统
CN107656207A (zh) * 2017-08-17 2018-02-02 同济大学 一种车用燃料电池发动机在动力总成回路的动态测试方法
CN107656207B (zh) * 2017-08-17 2021-02-02 同济大学 一种车用燃料电池发动机在动力总成回路的动态测试方法
CN108957359A (zh) * 2018-08-15 2018-12-07 合肥工业大学 一种用于整车环境下复合电源系统瞬时效率的测试方法
CN110194179A (zh) * 2019-06-26 2019-09-03 重庆大学 一种串联式混合动力电动汽车动力模式的确定系统
CN112906296A (zh) * 2021-02-02 2021-06-04 武汉理工大学 混合动力汽车全服役期能量优化方法、系统和存储介质
CN112906296B (zh) * 2021-02-02 2022-05-10 武汉理工大学 混合动力汽车全服役期能量优化方法、系统和存储介质
CN116666787A (zh) * 2023-07-28 2023-08-29 中国标准化研究院 一种锂电池智能管理系统
CN116666787B (zh) * 2023-07-28 2023-10-17 中国标准化研究院 一种锂电池智能管理系统

Also Published As

Publication number Publication date
CN104112036B (zh) 2017-03-01

Similar Documents

Publication Publication Date Title
CN104112036A (zh) 混联式混合动力电动汽车的仿真方法
CN106427599B (zh) 电动车辆的控制装置
Mansour et al. Dynamic modeling of the electro-mechanical configuration of the Toyota Hybrid System series/parallel power train
US10906385B2 (en) Hybrid vehicle powertrains with flywheel energy storage systems
EP2125414B1 (en) Hybrid vehicle auxiliary equipment energy management
CN103909922A (zh) 串联式混合动力汽车的整车控制策略
CN104309490A (zh) 电动汽车制动能量回收装置及方法
CN103631149A (zh) 一种增程式电动汽车续驶里程的仿真系统及其仿真方法
CN103338999A (zh) 用于混合动力车辆的驱动控制装置和驱动控制方法及混合动力车辆
CN104290591A (zh) 一种混联式混合动力拖拉机动力系统及其控制方法
CN201506357U (zh) 混合动力输出功率平衡装置
CN102390248A (zh) 单离合器混合动力系统
CN115879275A (zh) 混合动力电动汽车性能仿真系统、方法、设备及存储介质
CN206202005U (zh) 一种用于纯电动车的动力系统
Sarathkumar et al. Modelling and simulation of electric vehicle drive through SAEJ227 & EUDC cycles
CN106526351B (zh) 混合动力专用车辆再生辅助动力试验台及其测试方法
CN201002506Y (zh) 混合动力汽车用的动力耦合装置
Gokce et al. Modeling and Simulation of a Serial–Parallel Hybrid Electrical Vehicle
CN106042893B (zh) 对蓄电装置的充放电进行控制的混合动力车辆
CN202200824U (zh) 单离合器混合动力系统
Zhou et al. Research on the energy management of composite energy storage system in electric vehicles
CN205239137U (zh) 一种内燃机改为油电混合动力装置
Krishna et al. Plug in Hybrid Vehicle-A Simulation Study
Dhand Design of electric vehicle propulsion system incorporating flywheel energy storage
Jung et al. Voltage control of electric vehicle permanent magnet synchronous motor for maximizing energy recovery rate

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170301

Termination date: 20180612