CN108875255B - 基于电动汽车实际行驶工况的永磁驱动电机温升分析方法 - Google Patents
基于电动汽车实际行驶工况的永磁驱动电机温升分析方法 Download PDFInfo
- Publication number
- CN108875255B CN108875255B CN201810721926.8A CN201810721926A CN108875255B CN 108875255 B CN108875255 B CN 108875255B CN 201810721926 A CN201810721926 A CN 201810721926A CN 108875255 B CN108875255 B CN 108875255B
- Authority
- CN
- China
- Prior art keywords
- permanent magnet
- driving motor
- actual
- electric automobile
- motor
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 44
- 238000004088 simulation Methods 0.000 claims abstract description 30
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 32
- 238000004804 winding Methods 0.000 claims description 30
- 238000000034 method Methods 0.000 claims description 22
- 230000005291 magnetic effect Effects 0.000 claims description 16
- 229910052742 iron Inorganic materials 0.000 claims description 15
- 238000001816 cooling Methods 0.000 claims description 13
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 12
- 239000010949 copper Substances 0.000 claims description 12
- 229910052802 copper Inorganic materials 0.000 claims description 12
- 230000033228 biological regulation Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 230000017525 heat dissipation Effects 0.000 claims description 5
- 239000003302 ferromagnetic material Substances 0.000 claims description 4
- 230000002093 peripheral effect Effects 0.000 claims description 4
- 238000005096 rolling process Methods 0.000 claims description 4
- 239000000110 cooling liquid Substances 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 239000011162 core material Substances 0.000 claims 2
- 239000011810 insulating material Substances 0.000 abstract description 4
- 230000008859 change Effects 0.000 description 6
- 238000011217 control strategy Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000000295 fuel oil Substances 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 238000003475 lamination Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005347 demagnetization Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/72—Electric energy management in electromobility
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
本发明公开了一种基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,包括以下步骤:S1、将实际电动汽车的电机模型嵌入到Advisor仿真软件中,并在给定该电动汽车行驶工况的条件下,记录电机的实际转矩和转速;S2、根据电机的实际转矩和转速,计算电动汽车行驶时永磁驱动电机的实际工作电流;S3、将永磁驱动电机的实际工作电流作为Ansys仿真模型的输入,分析永磁驱动电机的温升。本发明提供的基于电动汽车实际行驶工况的永磁驱动电机温升计算方法,能够更准确地计算电动汽车实际行驶时驱动电机的温升分布,可以更好地发挥绝缘材料性能,降低系统成本,提高电动汽车驱动电机运行的稳定性和可靠性。
Description
技术领域
本发明属于电机温升计算技术领域,具体涉及一种基于电动汽车实际行驶工况的永磁驱动电机温升分析方法。
背景技术
20世纪70年代以来,全球能源危机与环境问题日益严重使得电动汽车行业快速发展。目前国内车用燃油过度消耗,石油对外依存度不断提高,环境问题日益加剧,而然后汽车尾气排放是主要污染源之一。在减少环境污染、环境能源危机的前提下,电动汽车因其污染小、节约能源的优点,处于重要的战略发展机遇期。
电动汽车是指以车载电源为动力,用电机驱动车轮行驶,负荷道路交通、安全法规各项要求的车辆。相较于传统燃油汽车,纯电动汽车以电池代替燃油,以驱动电机代替内燃机,从而将电能转化为机械能。目前,纯电动汽车用驱动电机主要有直流电机、永磁同步电机、交流异步电机和开关磁阻电机。而电动汽车驱动电机的温升计算时电机设计阶段最重要的校验手段之一,对电机运行是的等效电阻,电感和效率等参数有显著影响,一直是电力行业广受关注的热点问题。
电动汽车作为一种交通运输工具,要求其驱动电机具有比工业牵引电机更高的功率密度和效率,准确计算驱动电机温升对电动汽车而言尤为重要,过高的温升计算结果会使电机设计尺寸和重量过大,不仅会侵占电动汽车有限的空间,也会缩短电动汽车的续驶里程,而过低的温升计算结果可能导致电机实际温升超过绝缘材料和永磁铁的最高允许温度,致使绕组绝缘的过早失效和永磁体退磁,严重时甚至会威胁到人员的生命安全。
当前对电动汽车驱动电机温升问题的研究,多采用工业牵引电机温升的分析方法,基于额定电流和冷却调节进行计算。而电动汽车驱动电流是按行驶需求不断变化的,采用该类方法计算电机温升必然会产生较大的误差。大量的统计数据表明,行驶工况能够较好地描述电动汽车行驶规律,因此采用行驶工况对电动汽车驱动电机温升进行分析具有良好的现实意义。在现有技术中,一种方法为利用整车工况计算驱动电机的转矩和转速,然后利用有限元仿真软件计算损耗,最后基于热网络分析电动机温升,但该方法需要实验测试来确定电机内损耗分布,很难在设计阶段对电机性能进行校核;另外一种方法中,利用有限元计算了UDDS工况下永磁驱动电机的损耗,进而给出了电机温升结果,但其无法与目前电动汽车广泛采用的矢量控制等先进控制策略直接接口,在某种程度上阻碍了该方法在工程实践中的实际应用。
发明内容
针对现有技术中的上述不足,本发明提供的基于电动汽车实际行驶工况的永磁驱动电机温升计算方法解决了现有技术中无法在实际行驶工况下对永磁驱动电机进行温升分析的问题。
为了达到上述发明目的,本发明采用的技术方案为:基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,包括以下步骤:
S1、将实际电动汽车的电机模型嵌入到Advisor仿真软件中,并在给定该电动汽车行驶工况的条件下,记录电机的实际转矩和转速;
S2、根据电机的实际转矩和转速,计算电动汽车行驶时永磁驱动电机的实际工作电流;
S3、将永磁驱动电机的实际工作电流作为Ansys仿真模型的输入,得到电动汽车实际行驶时永磁驱动电机的温升随其实际工作电流变化的规律。
进一步地,所述步骤S1具体为:
S11、根据永磁驱动电机所装配的实际电动汽车的整车参数,确定Advisor整车仿真模型中的整车参数;
S12、采用理论计算方法计算实际永磁驱动电机输出特性曲线,得到不同转速对应的电机转矩;
S13、根据转速和转矩数据,查找电机手册数据,给定不同转速下电机损耗,并调整Advisor整车仿真模型中的所对应的损耗数据;
S14、在Advisor软件中,给定电动汽车行驶工况,运行该行驶工况若干周期,记录电动汽车永磁电动机实际转矩和转速。
进一步地,所述步骤S12中,额定转速下采用恒转矩调速,输出额定转矩,额定转速以上采用横定功率调速,输出恒功率,得到永磁驱动电机输出特性曲线。
进一步地,其特征在于,所述步骤S2具体为:
S21、计算电机电磁转矩标幺值,进而计算永磁驱动电机d轴和q轴电流的标幺值;
S22、根据永磁驱动电机d轴和q轴电流的标幺值,分别计算d轴和q轴旋转电流的瞬时值,并将其变换为三相电流瞬时值;
S23、利用数理统计方法对所得永磁驱动电机三相电流瞬时值进行统计,得到永磁驱动电机实际工作电流的变化曲线。
进一步地,其特征在于,所述步骤S21中:
所述电机电磁转矩标幺值Tb的计算公式为:
Tb=TE/TB
式中,TB为电磁转矩基值,且TB=nPψFiB;
TE为实际电磁转矩,且TE=nP(ψFiQ+(LD-LQ)iDiQ);
nP为极对数;
ΨF为永磁铁体空间磁链;
LD、LQ分别为永磁驱动电机的直、交轴等效电感;
所述永磁驱动电机d轴电流的标幺值id的计算公式为:
id=iD/iB
所述永磁驱动电机q轴电流的标幺值iq的计算公式为:
iq=iQ/iB
式中,iB为电流基值,且iB=ψF/LQ-LD;
iD、iQ分别为永磁驱动电机的直、交轴电流。
进一步地,其特征在于,所述步骤S22具体为:
所述d轴旋转电流的瞬时值iQd的计算公式为:
iQd=iqiB
所述q轴旋转电流的瞬时值iQq的计算公式为:
iQq=idiB
将所述d轴和q轴旋转电流的瞬时值变换为三相电流瞬时值的方法具体为:
通过park逆变换将旋转的d轴和q轴的选择电流瞬时值转变为静止的α、β轴电流瞬时值,然后通过clark逆变换将静止的α、β轴电流瞬时值变换为三相电流瞬时值;
所述三相电流瞬时值为电动汽车永磁驱动电机三相输入电流的瞬时值。
进一步地,其特征在于,所述步骤S3具体为:
确定永磁驱动电机的结构参数和冷却方式,将永磁驱动电机的损耗作为热源,将永磁驱动电机的实际工作电流作为Ansys仿真模型的输入,在给定的若干运行工况下运行,在电机温升达到稳定时,得到电动汽车实际行驶时永磁驱动电机的温升随其实际工作电流变化的规律。
进一步地,其特征在于,所述永磁驱动电机的结构参数包括绕组匝数尺寸、永磁铁的形状位置、永磁铁的槽数、永磁铁的槽型尺寸、永磁铁的铁心叠片材料尺寸;
所述冷却方式为包括各部分材料的散热系数、冷却液流量、冷却通道的内径、外径和圆心角。
进一步地,其特征在于,所述永磁驱动电机的损耗包括铜损、铁耗和轴承摩擦损耗;
所述铜损为各相绕组上产生的损耗,所述铜损PCu计算公式为:
式中,I为绕组相电流;
R为绕组电阻;
m为相数;
所述铁耗为交变磁场穿过定子铁芯时,产生的涡流损耗和磁滞损耗,所述铁耗Pfe的计算公式为:
式中,Pfe为单位质量的铁耗;
ω为角频率;
Bp磁密峰值;
γ1、γ2、ε均是铁磁材料有关的常数;
所述轴承摩擦损耗Pf为电动汽车内置式永磁驱动电机滚动轴承的摩擦损耗,所述轴承摩擦损耗的计算公式为:
式中,Q为轴承载荷;
d为滚珠中心处的直径;
v为滚珠中心的圆周速度。
本发明的有益效果为:本发明提供的基于电动汽车实际行驶工况的永磁驱动电机温升计算方法,能够更准确地计算电动汽车实际行驶时驱动电机的温升分布,可以更好地发挥绝缘材料性能,降低系统成本,提高电动汽车驱动电机运行的稳定性和可靠性。
附图说明
图1为本发明提供的实施例中基于电动汽车实际行驶工况的永磁驱动电机温升分析方法实现流程图。
图2为本发明提供的实施例中计算电动汽车行驶时永磁驱动电机的实际工作电流方法流程图。
图3为本发明提供的实施例中分析永磁驱动电机的温升方法流程图。
图4为本发明提供的实施例中Advisor中电动汽车驱动电机的电流逆解耦模型结构图。
图5为本发明提供的实施例中在市区工况UDC中电动汽车电机的转矩和瞬时工作电流对应图。
图6为本发明提供的实施例中市郊工况EUDC中电动汽车电机转矩和瞬时工作电流对应图。
图7为本发明提供的实施例中UDC工况下电机定子槽温度分布图。
图8为本发明提供的实施例中EUDC工况下电机定子槽温度分布图。
图9为本发明提供的实施例中UDC工况下定子绕组仿真和实测结果对比图。
图10为本发明提供的实施例中EUDC工况下定子绕组温升仿真和实测结果对比图。
图11为本发明提供的实施例中UDC工况下永磁铁温升仿真和实测结果对比图。
图12为本发明提供的实施例中EUDC工况下永磁铁温升仿真和实测结果对比图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,包括以下步骤:
S1、将实际电动汽车的电机模型嵌入到Advisor仿真软件中,并在给定该电动汽车行驶工况的条件下,记录电机的实际转矩和转速。
如图2所示,上述步骤S1具体为:
S11、根据永磁驱动电机所装配的实际电动汽车的整车参数,确定Advisor整车仿真模型中的整车参数。
分析永磁驱动电机所装配实际电动汽车的整车参数,在Advisor整车仿真模型中选择和该实际车型相近的参数,并修改Advisor整车仿真模型中的整车参数(如整车重量、迎风面积等)使其与实际车型的参数一致。
S12、采用理论计算方法计算实际永磁驱动电机输出特性曲线,得到不同转速对应的电机转矩。
上述步骤中,在额定转速下以采用恒转速的调速方式,在额定转速以上时采用恒功率的调速方式,计算实际永磁驱动电机输出特性曲线,即得到关于不同转速对应的电机转矩。
S13、根据转速和转矩数据,查找电机手册数据,给定不同转速下电机损耗,并调整Advisor整车仿真模型中的所对应的损耗数据。
同时将所得转矩与转速数据带入到电动汽车Advisor整车仿真模型的m文件中,并按照实际电机的参数修改m文件中的峰值转矩,最高转速等数据。
S14、在Advisor软件中,给定电动汽车行驶工况,运行该行驶工况若干周期,记录电动汽车永磁电动机实际转矩和转速。
S2、根据电机的实际转矩和转速,计算电动汽车行驶时永磁驱动电机的实际工作电流;
如图3所示,上述步骤S2具体为:
S21、计算电机电磁转矩标幺值,进而计算永磁驱动电机d轴和q轴电流的标幺值;
上述步骤S21中:
电机电磁转矩标幺值Tb的计算公式为:
Tb=TE/TB (1)
式中,TB为转矩基值,且TB=nPψFiB;
TE为实际电磁转矩,且TE=nP(ψFiQ+(LD-LQ)iDiQ);
nP为极对数;
ΨF为永磁铁体空间磁链;
LD、LQ分别为永磁驱动电机的直、交轴等效电感;
永磁驱动电机d轴电流的标幺值id的计算公式为:
id=iD/iB (2)
所述永磁驱动电机q轴电流的标幺值iq的计算公式为:
iq=iQ/iB (3)
式中,iB为电流基值,且iB=ψF/LQ-LD;
iD、iQ分别为永磁驱动电机的直、交轴电流。
根据上述公式中的标幺值定义,得到永磁驱动电机电磁转矩的标幺值计算式为:
Tb=iq(1-id) (4)
为提高电动汽车驱动电机的输出性能,多数电动汽车采用最大转矩电流比控制策略对永磁驱动电机进行驱动,最大转矩电流比控制策略下永磁电机的实际驱动电流的d、q轴分量满足:
式中,is=i/iB,i定子电流的瞬时值;
基于拉格朗日极值定理,得到方程:
对上述方程(6)进行求解,并消去λ,得到:
求解方程(7)中中iq的反函数,由于公式(7)中第二行方程中id的反函数不存在),因此需要借助其他方程,将iq反函数求解结果带入到方程Tb-(iq(1-id))中,进而给出id对于Tb和iq反函数,并搭建相应的Simulink电流逆解耦模型,即可计算Advisor仿真模型中电动汽车永磁驱动电机的d轴和q轴电流的标幺值。
S22、根据永磁驱动电机d轴和q轴电流的标幺值,分别计算d轴和q轴旋转电流的瞬时值,并将其变换为三相电流瞬时值;
上述步骤S22具体为:
d轴旋转电流的瞬时值iQd的计算公式为:
iQd=iqiB (8)
q轴旋转电流的瞬时值iQq的计算公式为:
iQq=idiB (9)
将d轴和q轴旋转电流的瞬时值变换为三相电流瞬时值的方法具体为:
通过park逆变换将旋转的d轴和q轴的选择电流瞬时值转变为静止的α、β轴电流瞬时值,然后通过clark逆变换将静止的α、β轴电流瞬时值变换为三相电流瞬时值;
三相电流瞬时值为电动汽车永磁驱动电机三相输入电流的瞬时值。
S23、利用数理统计方法对所得永磁驱动电机三相电流瞬时值进行统计,得到永磁驱动电机实际工作电流的变化曲线。
S3、将永磁驱动电机的实际工作电流作为Ansys仿真模型的输入,得到电动汽车实际行驶时用词驱动电机的温升随实际工作电流变化的规律。
上述步骤S3具体为:
确定永磁驱动电机的结构参数和冷却方式,将永磁驱动电机的损耗作为热源,将永磁驱动电机的实际工作电流作为Ansys仿真模型的输入,在给定的若干运行工况下运行,在电机温升达到稳定时,得到电动汽车实际行驶时永磁驱动电机的温升随其实际工作电流变化的规律。
其中,永磁驱动电机的结构参数包括绕组匝数尺寸、永磁铁的形状位置、永磁铁的槽数、永磁铁的槽型尺寸、永磁铁的铁心叠片材料尺寸等等;
冷却方式为包括各部分材料的散热系数、冷却液流量、冷却通道的内径、外径和圆心角等等。
其中,永磁驱动电机的损耗包括铜损、铁耗和轴承摩擦损耗;
铜损为各相绕组上产生的损耗,铜损PCu计算公式为:
式中,I为绕组相电流;
R为绕组电阻;
m为相数;
铁耗为交变磁场穿过定子铁芯时,产生的涡流损耗和磁滞损耗,铁耗Pfe的计算公式为:
式中,Pfe为单位质量的铁耗;
ω为角频率;
Bp磁密峰值;
γ1、γ2、ε均是铁磁材料有关的常数;
轴承摩擦损耗Pf为电动汽车内置式永磁驱动电机滚动轴承的摩擦损耗,轴承摩擦损耗的计算公式为:
式中,Q为轴承载荷;
d为滚珠中心处的直径;
v为滚珠中心的圆周速度。
需要说明的是,在电动汽车处于加速或减速工况下,为满足实际行驶需求,永磁驱动电机的输入电流在该加速的时段内,电流呈现上升或下降的直线,在加速度不变的情况下,电流瞬时值曲线的斜率不变,上述方法中取该时段内的电流平均值作为Ansys中永磁驱动电机模型的输入。
在本发明的一个实施例中,提供了将本发明提供的方法的应用实例:
本实例以内置式永磁同步电机作为温升的研究对象,采用纯电动仿真模型,整车及内置式永磁同步电机相关参数如表1所示:
表1电机及整车参数表
在输出转矩小于额定转矩时,永磁驱动电机采用恒转矩的调速方式,在输出转矩大于或等于额定转矩时,永磁驱动电机采用恒功率的调速方式,其转速—转矩参数如表2所示:
表2电机转矩和转速
将表1和表2的参数嵌入到Advisor软件中,对实际工况进行仿真分析,本实例中的仿真工况为新欧洲循环工况中的市区典型工况UDC和市郊典型工况EUDC,将表2所示的内置式永磁驱动电机的输出特性嵌入到表1的电动汽车模型中,对两种工况循环运行若干周期后的仿真结果表明,22kW永磁电机能够较好地实现电动汽车的动力性能,证明了所选电机功率和电动机车型的匹配是合理的。
基于拉格朗日定理的电流逆解模型求解永磁驱动电机实际工作电流的过程为:
永磁电机的电磁转矩T满足:
T=p(ψFiQ+LDiDiQ-LQiDiQ) (13)
式中,p为极对数;
ΨF为永磁铁空间磁链;
LD、LQ分别为直、交轴等效电感;
iD、iQ分别为直、交轴电流。
为简化计算过程,采用标幺值进行计算,选取转矩和电流的基值IB、TB分别为
IB=ψF/LQ-LD (14)
TB=pψFiB (15)
则式(13)的标幺化形式为:
Tnor=iq(1-id) (16)
由于最大转矩电流比控制策略在永磁电机驱动中具有较好的优越性;本实例在此策略下建立内置式永磁电机的工况电流进行逆解耦模型,其d、q轴电流分量满足
式中,inor=i/IB,i是定子电流的瞬时值。
最大转矩电流比控制的基本思想是利用最小的电流下产生最大的转矩。为求解电流最小值,综合考虑式(16)、(17)并构造拉格朗日极值函数
式中,λ为拉格朗日算子,分别对式(6)计算id,iq和δ的偏微分,可得
Tnor=(id(id-1)3)0.5 (20)
利用式(20)结合式(16),即可给出Advisor中电动汽车驱动电机的电流逆解耦模型,如图4所示。
图4中,电流基准值iB可通过电机给定数据得到,基于图4的电流逆解耦模型,以新欧洲驱动循环中的市区工况UDC和市郊工况EUDC为例,分析其电动汽车电机的转矩和瞬时工作电流分别如图5和图6所示。
图5和图6中上半部分为两种工况下驱动电机的实际输出转矩Te,下半部分为两种工况下驱动电机的瞬时工作电流,为正弦电流波形,为行文简洁,仅给出了一个工况内的波形。
结合图4可以看出,电机的工况电流随着转矩和转速的变化而变化,并且,在汽车加速或减速时,所需的驱动转矩和电流也较大。行驶速度不变时,转矩和电流也保持不变,表明了所给电流逆解耦模型的有效性。
在确定驱动电机温升随工作电流变化时,首先进行永磁电机热源的计算:
(1)铜损:绕组铜损耗包括各相绕组上产生的损耗:
式中,I为绕组相电流,R为绕组电阻,m相数。
(2)铁耗:交变磁场穿过定子铁心时,会产生涡流损耗和磁滞损耗。涡流损耗和磁滞损耗统称为铁芯损耗,一般将其统一计算为:
式中,Pfe为单位质量的铁耗,w为角频率,Bp磁密峰值,γ1、γ2、ε均是铁磁材料有关的常数
(3)轴承摩擦损耗
电动汽车内置式永磁驱动电机滚动轴承的摩擦损耗可通过下式计算:
式中,Q为轴承载荷(N),d为滚珠(或滚柱)中心处的直径(m),v为滚珠中心的圆周速度(m/s)。
为简化利用瞬态工作电流分析驱动电机温升问题,根据电机的实际结构和传热学理论做出以下假设:
(1)电机绕组,定转子内的损耗都是均匀分布的,并且电机的外部冷却条件相同;
(2)忽略温差造成的定子绕组电阻率的变化。
根据上述假设,以及相应的导热微分方程和边界条件,可得此各向异性介质模型的二维瞬态温度场的边值问题为:
式中:λ为导热系数;ρ为材料密度;cp为定压比热容;qV为已知热流密度或热流密度函数;Tf为冷却介质温度;α为边界上的散热系数;S为求解区域的边界,且方向为逆时针方向;T为环境温度;t为时间;x、y、z为坐标系坐标。
上式的等价变分问题为:
基于式(26)的导热微分方程和内置式永磁电机损耗计算式,考虑到永磁电机的机壳外表面周围介质的温度以及该面与周围介质之间的热交换系数均为已知,满足第三类边界条件,将Advisor模型中抽取的电机实际电流作为Ansys模型的输入,可以给出得电动汽车内置式永磁驱动电机内温度场分布的定解。
确定永磁驱动电机的参数中:电机采用往复水冷方式,水流量为1L/min,扇形内径为0.215m,外径为0.225m,对应的圆心角为7.8度。
分别运行10个UDC和EUDC循环工况后,驱动电机二维温度场分布分别如图7、8所示。
图7是在UDC驱动工况下的驱动电机温度场分布,图8是在驱动EUDC工况下电动汽车驱动内置式永磁电机温度场的分布。驱动电机的温度变化情况仿真时间分别为2000s和4500s.
对比图7~8可以看出,在UDC和EUDC工况下内置式永磁电机的槽部温度不同,这是由于不同行驶工况下电动汽车对驱动电机的需求不同,导致不同的工况下,槽部温度存在着较大差异。
定子绕组电流随时间变化较大,定子铜耗以及槽部的发热密度变化显著,槽部温度分布从绕组中心向四周呈辐射状降低。而定子铁心铁耗和发热密度变化相对较小,定子热阻发热密度较高,产生的热量通过绝缘传递到定子铁心,再经由外部水冷系统带走。为此,电机最高温度出现在定子绕组上,定子铁心的温度低于绕组内部温度。
从图7和图8可以进一步看出,在永磁体的不同部位其温度分布不同,这是因为转子侧损耗主要为集中在转子表面附件的铁心附件损耗,相对定子侧损耗大小和发热密度均大大减小,热量分别向气隙和转子铁心传递,为此从转子铁心外表面向内部的温度逐渐降低。
另外,在同一工况下,不同时刻内置式永磁驱动电机的槽部温度不同,这是由于在不同时刻的电动汽车的行驶状态不同,导致驱动电机的定子电流不同,因此,不同时刻驱动电机槽部的温度是变化的。
为对所提基于工况的电动汽车驱动电机温升分析方法进行验证,搭建自行设计的具有表1和表2所给定参数的内置式永磁电机实验台架,将图5和图6所示的电流作为变频器输出的电机驱动电流,利用电动汽车测试系统,通过埋在定子绕组和磁钢附近的Pt100传感器检测温度。
分别连续运行20个UDC和EUDC工况,得到热平衡时内置式永磁电机的动态温升,不同工况下的定子绕组和永磁铁温升的仿真和实测结果分别如图9-12所示。
为对比清晰,图9-图12仅给出了达到热平衡时一个周期内的暂态温升。可以看出,内置式永磁电机定子绕组和永磁铁仿真和实验结果的一致性较好。实测的温升结果略低于仿真计算结果,原因可能是仿真温升分析的散热系数和热阻等参数计算过于理想,同时PT100传感器埋设位置无法和仿真计算位置具有一定的偏差。
在同一行驶工况下,由于实际电机驱动电流的周期性变化,电动汽车驱动电机温升也呈现有规律的周期变化。本实例的结果表明,UDC和EUDC工况下,定子绕组温升的变化范围均约为30℃,变化幅度较大。永磁体温升的变化范围均约为20℃,变化幅度略小于定子绕组。
在两种行驶工况下电动汽车驱动电机温升不同。其中,EUDC工况的温升略高于UDC,结合图5和图6可知,EUDC工况下电动汽车驱动电机的工作电流较大,导致在该工况运行下电机绕组和永磁铁温升高于UDC工况。
对图9-12结果进一步计算可知,在UDC工况下,所给仿真分析方法和实测温升的最大绝对误差为4.8℃,EUDC工况下的最大绝对误差为3.2℃。两种工况下仿真和实测温升的相对误差均小于2%,证明了基于工况电动汽车温升分析方法具有较好的精度。
本发明的有益效果为:本发明提供的基于电动汽车实际行驶工况的永磁驱动电机温升计算方法,能够更准确地计算电动汽车实际行驶时驱动电机的温升分布,可以更好地发挥绝缘材料性能,降低系统成本,提高电动汽车驱动电机运行的稳定性和可靠性。
Claims (6)
1.基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,其特征在于,包括以下步骤:
S1、将实际电动汽车的电机模型嵌入到Advisor仿真软件中,并在给定该电动汽车行驶工况的条件下,记录电机的实际转矩和转速;
S2、根据电机的实际转矩和转速,计算电动汽车行驶时永磁驱动电机的实际工作电流;
S3、将永磁驱动电机的实际工作电流作为Ansys仿真模型的输入,得到电动汽车实际行驶时永磁驱动电机的温升随其实际工作电流变化的规律;
所述步骤S1具体为:
S11、根据永磁驱动电机所装配的实际电动汽车的整车参数,确定Advisor整车仿真模型中的整车参数;
S12、采用理论计算方法计算实际永磁驱动电机输出特性曲线,得到不同转速对应的电机转矩;
S13、根据转速和转矩数据,查找电机手册数据,给定不同转速下电机损耗,并调整Advisor整车仿真模型中的所对应的损耗数据;
S14、在Advisor软件中,给定电动汽车行驶工况,运行该行驶工况若干周期,记录电动汽车永磁电动机实际转矩和转速;
所述步骤S12中,额定转速下采用恒转矩调速,输出额定转矩,额定转速以上采用恒定功率调速,输出恒功率,得到永磁驱动电机输出特性曲线;
所述步骤S2具体为:
S21、计算电机电磁转矩标幺值,进而计算永磁驱动电机d轴和q轴电流的标幺值;
S22、根据永磁驱动电机d轴和q轴电流的标幺值,分别计算d轴和q轴旋转电流的瞬时值,并将其变换为三相电流瞬时值;
S23、利用数理统计方法对所得永磁驱动电机三相电流瞬时值进行统计,得到永磁驱动电机实际工作电流的变化曲线。
2.根据权利要求1所述的基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,其特征在于,所述步骤S21中:
所述电机电磁转矩标幺值Tb的计算公式为:
Tb=TE/TB
式中,TB为电磁转矩基值,且TB=nPψFiB;
TE为实际电磁转矩,且TE=nP(ψFiQ+(LD-LQ)iDiQ);
nP为极对数;
ΨF为永磁铁体空间磁链;
LD、LQ分别为永磁驱动电机的直、交轴等效电感;
所述永磁驱动电机d轴电流的标幺值id的计算公式为:
id=iD/iB
所述永磁驱动电机q轴电流的标幺值iq的计算公式为:
iq=iQ/iB
式中,iB为电流基值,且iB=ψF/LQ-LD;
iD、iQ分别为永磁驱动电机的直、交轴电流。
3.根据权利要求1所述的基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,其特征在于,所述步骤S22具体为:
所述d轴旋转电流的瞬时值iQd的计算公式为:
iQd=iqiB
所述q轴旋转电流的瞬时值iQq的计算公式为:
iQq=idiB
将所述d轴和q轴旋转电流的瞬时值变换为三相电流瞬时值的方法具体为:
通过park逆变换将旋转的d轴和q轴的选择电流瞬时值转变为静止的α、β轴电流瞬时值,然后通过clark逆变换将静止的α、β轴电流瞬时值变换为三相电流瞬时值;
所述三相电流瞬时值为电动汽车永磁驱动电机三相输入电流的瞬时值。
4.根据权利要求1所述的基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,其特征在于,所述步骤S3具体为:
确定永磁驱动电机的结构参数和冷却方式,将永磁驱动电机的损耗作为热源,将永磁驱动电机的实际工作电流作为Ansys仿真模型的输入,在给定的若干运行工况下运行,在电机温升达到稳定时,得到电动汽车实际行驶时永磁驱动电机的温升随其实际工作电流变化的规律。
5.根据权利要求4所述的基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,其特征在于,所述永磁驱动电机的结构参数包括绕组匝数尺寸、永磁铁的形状位置、永磁铁的槽数、永磁铁的槽型尺寸、永磁铁的铁心叠片材料尺寸;
所述冷却方式为包括各部分材料的散热系数、冷却液流量、冷却通道的内径、外径和圆心角。
6.根据权利要求4所述的基于电动汽车实际行驶工况的永磁驱动电机温升分析方法,其特征在于,所述永磁驱动电机的损耗包括铜损、铁耗和轴承摩擦损耗;
所述铜损为各相绕组上产生的损耗,所述铜损PCu计算公式为:
式中,I为绕组相电流;
R为绕组电阻;
m为相数;
所述铁耗为交变磁场穿过定子铁芯时,产生的涡流损耗和磁滞损耗,所述铁耗Pfe的计算公式为:
式中,Pfe为单位质量的铁耗;
ω为角频率;
Bp磁密峰值;
γ1、γ2、ε均是铁磁材料有关的常数;
所述轴承摩擦损耗Pf为电动汽车内置式永磁驱动电机滚动轴承的摩擦损耗,所述轴承摩擦损耗的计算公式为:
式中,Q为轴承载荷;
d为滚珠中心处的直径;
v为滚珠中心的圆周速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810721926.8A CN108875255B (zh) | 2018-07-04 | 2018-07-04 | 基于电动汽车实际行驶工况的永磁驱动电机温升分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810721926.8A CN108875255B (zh) | 2018-07-04 | 2018-07-04 | 基于电动汽车实际行驶工况的永磁驱动电机温升分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108875255A CN108875255A (zh) | 2018-11-23 |
CN108875255B true CN108875255B (zh) | 2022-08-16 |
Family
ID=64298945
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810721926.8A Active CN108875255B (zh) | 2018-07-04 | 2018-07-04 | 基于电动汽车实际行驶工况的永磁驱动电机温升分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108875255B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110442944B (zh) * | 2019-07-29 | 2023-03-31 | 江苏大学 | 一种基于多工况运行的变漏磁永磁同步电机的建模方法 |
CN113394922B (zh) * | 2020-03-13 | 2022-06-14 | 中车永济电机有限公司 | 永磁电机调整方法、装置、设备以及存储介质 |
CN115422780A (zh) * | 2022-10-12 | 2022-12-02 | 中车长春轨道客车股份有限公司 | 牵引系统的温度仿真方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103308325A (zh) * | 2013-06-26 | 2013-09-18 | 东莞中山大学研究院 | 电动汽车驱动系统半实物仿真平台 |
CN104007358A (zh) * | 2014-06-11 | 2014-08-27 | 安徽江淮汽车股份有限公司 | 电动汽车永磁同步电机驱动系统短路故障诊断方法及系统 |
CN104777825A (zh) * | 2014-12-29 | 2015-07-15 | 吉林大学 | 轮毂电机系统转矩控制方法及转矩控制方法的检验台架 |
CN105844026A (zh) * | 2016-03-25 | 2016-08-10 | 东南大学 | 基于田口方法确定永磁直线电机工作电流与时间的方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002006076A1 (de) * | 2000-07-17 | 2002-01-24 | Abm Greiffenberger Antriebstechnik Gmbh | Verfahren zur sensorlosen antriebsregelung eines elektrofahrzeugs sowie danach arbeitende antriebsregelung |
US20040007997A1 (en) * | 2002-07-11 | 2004-01-15 | Visteon Global Technologies, Inc. | Vector control system for permanent magnet sychronous machines using parameter scheduling table |
CN101769797A (zh) * | 2009-01-06 | 2010-07-07 | 李虎 | 一种用于预测永磁同步电动机中永磁体温度的温升分析方法 |
CN102331719B (zh) * | 2011-07-11 | 2013-04-17 | 联合汽车电子有限公司 | 基于热模型的混合动力汽车动态降额控制方法 |
CN106240341B (zh) * | 2016-08-05 | 2018-10-09 | 武汉理工大学 | 一种电动汽车用永磁同步电机冷却系统及其控制方法 |
CN106599439A (zh) * | 2016-12-08 | 2017-04-26 | 重庆大学 | 面向能耗的纯电动汽车双电机动力系统参数优化匹配方法 |
-
2018
- 2018-07-04 CN CN201810721926.8A patent/CN108875255B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103308325A (zh) * | 2013-06-26 | 2013-09-18 | 东莞中山大学研究院 | 电动汽车驱动系统半实物仿真平台 |
CN104007358A (zh) * | 2014-06-11 | 2014-08-27 | 安徽江淮汽车股份有限公司 | 电动汽车永磁同步电机驱动系统短路故障诊断方法及系统 |
CN104777825A (zh) * | 2014-12-29 | 2015-07-15 | 吉林大学 | 轮毂电机系统转矩控制方法及转矩控制方法的检验台架 |
CN105844026A (zh) * | 2016-03-25 | 2016-08-10 | 东南大学 | 基于田口方法确定永磁直线电机工作电流与时间的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108875255A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Performance analysis of an axial magnetic-field-modulated brushless double-rotor machine for hybrid electric vehicles | |
CN108875255B (zh) | 基于电动汽车实际行驶工况的永磁驱动电机温升分析方法 | |
CN104331573A (zh) | 无刷复合结构电机系统的优化设计方法 | |
Si et al. | Analysis of temperature field for a surface-mounted and interior permanent magnet synchronous motor adopting magnetic-thermal coupling method | |
Chen et al. | Calculation analysis of thermal loss and temperature field of in-wheel motor in micro-electric vehicle | |
Villan et al. | Experimental comparison between induction and synchronous reluctance motor-drives | |
Zhang et al. | A cast copper rotor induction motor for small commercial EV traction: Electromagnetic design, analysis, and experimental tests | |
Cheng et al. | Research on induction motor for mini electric vehicles | |
Qinglong et al. | Indirect field oriented control technology for asynchronous motor of electric vehicle | |
An et al. | Calculation model of armature reaction magnetic field of interior permanent magnet synchronous motor with segmented skewed poles | |
Li et al. | Optimization of external rotor surface permanent magnet machines based on efficiency map over a target driving cycle | |
Sun et al. | Analysis of torque ripple and fault-tolerant capability for a 16/10 segmented switched reluctance motor in HEV applications | |
Liu et al. | Kriging assisted on-line torque calculation for brushless DC motors used in electric vehicles | |
JP2013119295A (ja) | ハイブリッド自動車用モータの性能解析方法 | |
Ba et al. | A generalized per-phase equivalent circuit model of the PMSM with predictable core loss | |
CN106934174B (zh) | 一种电力推进船舶实验平台的建模及稳定性分析方法 | |
CN106844923A (zh) | 一种多相永磁电机的参数化设计方法 | |
CN106646230A (zh) | 一种高速电机测试加载系统及采用该系统实现的负载功率调节方法 | |
Song et al. | Efficiency map calculation for surface-mounted permanent-magnet in-wheel motor based on design parameters and control strategy | |
Kitzberger et al. | Influence of hysteresis and eddy current losses on electric drive energy balance in driving cycle operation | |
Zhu | The key technologies for powertrain system of intelligent vehicles based on switched reluctance motors | |
Melo et al. | Evaluation of an energy loss-minimization algorithm for EVs based on induction motor | |
Rossi et al. | Design and thermal assessment of a high performance electric motor for racing applications | |
Rezaeealam | Permanent Magnet Tubular Generator with Quasi-Halbach Array for Free-Piston Generator System | |
Finken et al. | Design rules for energy efficient IPM motors in HEV applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240101 Address after: Room 10, 07, 3rd Floor, Hainan Ruize Office Building, Yazhou Bay Science and Technology City, Yazhou District, Sanya City, Hainan Province, 572000 Patentee after: Hainan Longsheng Technology Trading Co.,Ltd. Address before: 150022 No. 2468 Puyuan Road, Songbei District, Harbin City, Heilongjiang Province Patentee before: HEILONGJIANG University OF SCIENCE AND TECHNOLOGY |