1. 引言
随着当前汽车行业的迅猛发展与汽车保有量的显著提升[1],其对能源资源的消耗及环境质量的负面影响日益加剧。鉴于此,汽车行业正逐步转型,致力于新能源汽车领域的深入探索。国内外汽车制造商纷纷投身于纯电动汽车、混合动力汽车以及氢燃料电池汽车等多种类型新能源汽车的研发之中。近年来,电动汽车总体市场占有率目前较高,赵礼辉等[2]提出了车辆使用年限主要集中在6~12年之间,车辆年行驶里程与期望使用年限间存在较好的幂函数关系,期望使用年限随年行驶里程的增加而减小,95%用户全寿命周期行驶里程小于70 × 104 km。
新能源汽车与传统燃油车相比,驱动电机系统代替了燃油机,使得新能源汽车具有清洁污染、效率高、噪音小等优点。此外,新能源汽车的动力由电驱动系统提供,电驱动系统中电机和减速器的集成化使得其体积更小、效率更高、运行平稳。
电驱动系统主要是由驱动电机、电机控制器与减速器这三大子系统组成的。对于电驱动系统驱动方式主要包括分布式驱动和集中式驱动。针对电驱动系统各部件的定义可以依据GB/T 18488.1-2015标准的介绍[3]。电驱动系统实际上是复杂的机电一体化系统,包含不同的零部件,不同的零部件具有不同的失效模式,关联多种物理载荷特征。对于电驱动系统可靠性的研究,首先需要明确各部件的失效机理,这一标准[4]旨在通过全面的可靠性试验方法,为电动汽车制造商提供系统的测试依据,确保电动汽车驱动电机系统的质量和可靠性,从而提高消费者对电动汽车的信任度。
结合电驱动系统零部件失效机理构造表征用户风险程度的特征参数,采用多维特征融合方法进行参数降维,建立合理的目标函数,然后采用非监督学习算法实现典型用户聚类,分析不同类别的用户与不同部件的关联,建立失效风险综合评价指标,开展不同用户等级同一部件损伤风险评估,开展同一用户等级不同部件损伤风险评估,定义用户损伤等级。
基于用户原始载荷时域数据,通过计算电驱动系统不同零部件的伪损伤,进而计算出不同零部件的损伤强度,结合聚类后用户风险等级与部件损伤关联结果,设计以损伤强度信号为输入,用户风险等级为输出的机器学习模型,实现云端新增用户的失效风险等级快速辨识。
2. 电驱动系统核心部件失效机理
电驱动系统中不同零部件的失效模式和失效机理也不相同,针对不同的零部件,有不同的受迷宫模型,针对不同的寿命损伤模型,构建不同的伪损伤计算方法,各核心部件其受到的载荷一方面承受汽车在长期服役过程中路面传递的激励载荷[5],另一方面承受电驱动系统在运行过程中各部件受到多物理场(力、热、电、磁)载荷的影响[6]。对于用户来说,不同用户的驾驶习惯也是完全不同的,即使地域环境与道路工况类型是相同的情况下,其载荷特征与损伤效果也会具有很大的差异性。因此这也会对电驱动系统可靠性提出了更高的要求,基于实际用户载荷分析电驱动系统各部件失效载荷与失效机理可为电驱动系统零部件设计提供参考依据,并指导产品可靠性试验验证。
电驱动系统中包含多种部件,各部件关联载荷也是不同的[7],并且不同部件的失效模式和失效机理也不相同。对于电驱动系统失效部件来说,主要包括机械结构部件(轴、齿轮、轴承、磁轭、壳体等);电子电气部件(控制单元、功率器件、电容、逆变器、绕组等)。对于轴类部件来说,失效模式通常为机械疲劳失效,其失效主导载荷来源于扭矩的变幅加载历程[8]。对于齿轮部件,其失效主要为接触疲劳损伤与弯曲疲劳损伤,在用户运行工况下,不同扭矩之间的交变也会对齿轮造成损伤,而转速的不同,也会影响着齿轮转动的圈数。对于轴承类部件来说,失效模式主要为接触疲劳,当轴承要在高速运行工况下进行运行的时候,会引起温度升高,进而加速轴承失效。
目前许多学者经过大量试验与对不同部件的损伤研究,针对不同的部件建立不同的寿命预测模型,最为典型的是Palmgren-Miner准则[9]。假设材料发生疲劳失效时[10],所吸收的能量为W,在不同载荷等级
作用到失效时的循环数为
,载荷
作用
次时材料所吸收的能量
,那么对于不同的载荷等级i,有:
(1)
当吸收的能量值达到极限时,材料发生失效,则有:
(2)
式(2)中,
个载荷等级下总的疲劳损伤为
,某一载荷下造成的疲劳损伤
与载荷循环比成线性关系。
在工程应用上,对于不同载荷谱作用于同一零部件时,采用伪损伤的方法可以来反映不同载荷谱对零件的损伤效果。伪损伤计算是以应力与寿命的曲线为基础,利用miner损伤准则,计算在不同载荷等级下所对应的损伤,进行累加损伤。S-N曲线的表达式为:
(3)
式(3)中,
为应力幅值,
为试件在幅值
作用下的疲劳寿命,
为常数、与材料有关,
为疲劳强度指数。
3. 风险用户等级定义
针对电驱动系统中各零部件的复杂失效机理与多样化的失效模式,本研究构建了一套详尽的特征参数体系,旨在精确量化并表征用户所面临的风险级别。为了显著提升数据处理效率并优化模型预测精度,我们引入了多维特征融合技术与主成分分析(PCA)方法[11],对原始高维数据进行有效降维处理。在此基础上,我们设计了合理的目标函数,并采纳了非监督学习领域的经典算法——Kmeans聚类算法,对用户群体进行了深入的聚类分析,以期实现典型用户特征的精准划分。
基于聚类分析的结果,我们进一步对不同等级的用户群体进行了细致的剖析。在同一等级的用户类群内部,我们评估了不同零部件的失效风险等级,以揭示其内部差异与共性特征。同时,我们也对比分析了不同等级用户在同一部件上的失效风险等级,以探究风险等级随用户群体变化的规律。
3.1. 特征参数构造
在探究风险用户等级划分的深层次过程中,特征参数的遴选与用户驾驶行为的关联性显得尤为重要。驾驶行为的多样性,涵盖操作习惯、行驶路径的选择偏好以及车辆维护保养的重视程度等多个维度,均会对机械部件的损耗程度产生直接或间接的影响。本研究致力于构建一个既科学又全面的特征参数体系框架,旨在实现对用户风险等级更为精准的量化评估。具体而言,我们选取了零部件轴的磨损状况、齿轮的损伤程度以及轴承的损坏情况作为核心特征参数,这些参数深刻而全面地映射出电驱动系统的整体损伤面貌,如表1所示,展示了部分用户的关于各部件的损伤。
Table 1. Feature parameter construction
表1. 特征参数构造
轴 |
齿轮 |
轴承 |
5.4E−07 |
8.8E−10 |
2.3E−08 |
1.2E−06 |
3.6E−09 |
3.9E−08 |
4.8E−07 |
6.4E−10 |
2.6E−08 |
4.7E−07 |
9.6E−10 |
2.4E−08 |
6.9E−07 |
1.6E−09 |
3.2E−08 |
5.3E−07 |
1.0E−09 |
2.5E−08 |
6.1E−07 |
9.8E−10 |
2.6E−08 |
3.2. 主成分分析
成分分析(PCA)是一种广泛应用的降维方法,它通过线性变换将原始的高维数据映射到一个低维空间中,同时尽可能保留数据中的主要变化模式。每个主成分都是原始变量的线性组合,它们按照方差大小进行排序,每一个主成分都是由原始变量的线性组合得到的结果,当我们通过选择前几个主成分的时候,我们可以在降低数据维度的同时也能保留大部分的有效信息。
(一) 主成分基本原理
主成分基本原理的基于数据的方差最大化思想。方差代表了数据的离散程度,方差越大,说明数据在该维度上的变化越丰富,所包含的信息也就越多[12]。因此,PCA通过找到一个正交变换矩阵,将原始数据投影到一个新的坐标系中,使得投影后的数据在每个维度上的方差尽可能大。这些新的维度就是主成分,它们之间互不相关,且能够尽可能多地保留原始数据的信息。
具体来说,PCA的实现过程包括以下几个步骤:
数据标准化:这一步是将原始数据进行标准化的处理,当标准化以后,原始的数据都会在0~1之间。这将会消除量纲和数值大小对分析结果的影响。
计算协方差矩阵:计算标准化后数据的协方差矩阵,反映各变量之间的相关性。
特征值分解:将计算得到的协方差矩阵进行特征值分解,分解后便会得到特征值和特征向量。
选择主成分:根据特征值的大小选择前k个主成分。
转换数据:将原始数据转换到由前k个主成分构成的新坐标系中,得到降维后的数据。
(二) 主成分分析结果
将本文构建所有的用户数据,依据轴、齿轮、轴承这三个核心部件的失效机理,与寿命模型,分别计算轴、齿轮、轴承的损伤,这三维特征参数进行主成分分析,依据累积信息贡献率大于90%提取出前两个主成分,如表2所示,该表展示了部分用户的主成分分析结果,在该表中,每一列代表一维主成分,每一行则代表一个用户。
Table 2. Principal component analysis results
表2. 主成分分析结果
主成分一 |
主成分二 |
−1.17 |
0.36 |
0.66 |
−1.05 |
0.81 |
0.84 |
−0.37 |
−0.17 |
−0.45 |
−0.20 |
−0.16 |
1.26 |
−0.37 |
−0.11 |
−0.60 |
0.66 |
0.28 |
0.18 |
3.3. 聚类分析
本文采用非监督学习领域的K-means聚类算法,旨在不明确预设用户分类标准的情况下,对用户群体进行典型用户划分。具体而言,我们首先设定聚类数目k值,即预设将用户划分为k类。在算法初始阶段,随机选定k个初始类簇中心。随后,依据点与点之间的欧氏距离(或曼哈顿距离、马氏距离),将每个用户(即数据点)分配到距离最近的类簇中心所代表的类别中。这一过程不断迭代进行,直至类簇中心点的位置变化趋于稳定,或达到预设的迭代次数上限。整个K-means聚类算法的具体实施流程如图1所示,通过此流程,我们能够有效地将用户群体划分为k个具有显著特征的典型用户类别。
Figure 1. K-means algorithm flow
图1. K-means算法流程
聚类中输入的样本数据为降维后的主成分得分矩阵,假设样本数据为矩阵
,其中
为样本个数,
为对应的特征维度。初始化
个聚类中心
,计算每个样本到每个聚类中心的欧式距离,其公式如下:
(4)
式(4)中,
表示第
个样本,
;
表示第
个聚类中心,
;
代表第
个样本的第
个特征,
;
代表第
个聚类中心的第
个属性。
将所有的距离通过算法得出后,通过比较每个样本到每个聚类中心之间的距离,进而将样本分配到距离其最近的聚类中心的簇中,便得到
个类簇。
3.4. 聚类数目确定
对于聚类数目
的选取,肘部法则(Elbow method)是随着
值增大时,计算损失函数下降的拐点,其损伤函数为簇内误差平方和(SSE),簇内的样本点到聚类中心的距离反映了簇内的凝聚度,则对于簇内误差平方和计算公式如下:
(5)
式(5)中,
是为第
类簇,
是
中的样本点,
是
的聚类中心点。
当聚类数目为2到7时,分别计算误差平方和SSE结果如图2所示。当聚类数目为5时,近似位于曲线拐点,原因是当聚类数目小于5时,损失函数SSE下降速度较快,各类别之间容易区分;当聚类数目大于5时,损伤函数SSE下降速度较慢,即使增加一个类别,SSE值下降很小,各类簇间不易区分。因此取聚类数目为5。
Figure 2. Error squared and SSE
图2. 误差平方和SSE
3.5. 聚类结果分析
通过K-means聚类分析以后,3703名用户被分为五类,将每类的用户数目进行统计,等级一有220名用户,等级二有410名用户,等级三有780名用户,等级四有532名用户,等级五有1761名用户,结果如表3所示。
Table 3. Number of users result
表3. 用户数目结果
等级一 |
等级二 |
等级三 |
等级四 |
等级五 |
220 |
410 |
780 |
532 |
1761 |
3.6. 同一等级类内不同部件失效风险等级评估
通过开展对同一用户等级类内不同部件失效风险等级评估,可以得到每类用户与各部件失效风险相关性,如图3所示,等级一的用户主要为齿轮的失效风险等级较高,等级二的用户主要为轴的失效风险等级较高,等级三的用户主要为轴承、齿轮的失效风险等级较高,等级四的用户对各个部件的损伤差异都比较小,等级五的用户对轴与轴承的损伤都比较大。
Figure 3. Assessment of the failure risk level of different components within the same class
图3. 同一等级类内不同部件失效风险等级评估
3.7. 不同等级同一部件失效风险等级评估
通过开展对不同等级用户同一部件失效风险等级评估,可以揭示各部件与各等级用户之间的失效风险相关性,如图4所示,等级一的用户主要为轴承失效风险等级最高,等级三的用户主要为轴、齿轮、绕组、IGBT的失效风险等级最高。
Figure 4. Assessment of the failure risk level of the same component at different levels
图4. 不同等级同一部件失效风险等级评估
4. 云端用户失效风险等级辨识
4.1. 随机森林算法模型
在机器学习中,随机森林作为一种监督式学习算法,广泛应用于分类和回归问题。它以决策树为基础进行构建,并通过集成多个决策树来形成一个强大的模型。随机森林并不依赖于整个训练数据集来构建单一的、大型的决策树。相反,它采用多个小的决策树,这些决策树是基于训练数据的不同子集和随机选择的特征属性来构建的。随后,随机森林将这些决策树的结果进行组合,以形成一个更为准确和稳健的模型。随机森林中的每个子集都是通过随机抽样的样本和随机选择的特征来创建的。这种随机性降低了决策树对训练数据的过度依赖,有助于减少过拟合的风险。通过这种方式,随机森林能够提升模型的泛化能力,使其在处理新数据时表现更为出色。
在决策树中,每个节点表示一个特征的测试,每个分支表示一个测试结果。决策树的目标是通过一系列的特征测试,将数据分成尽可能纯的子集。这个过程是通过递归地选择最优特征进行分割来实现的,最优特征的选择通常基于某种度量标准,如信息增益,信息增益越大,数据集的纯度提升的越多,不确定性减少的越多,信息熵用于度量一个随机变量不确定性的大小,是度量信息含量的一个指标。计算公式如下
(6)
式中,x是数据集,n是数据集中的类别数,pi是第i个类别在数据集中出现的概率。
4.2. 用户特征参数输入
风险用户等级定义中有五类用户,将五类用户输入随机森林算法模型进行训练,如表4所示,等级一共220名用户,等级二共410名用户,等级三共780名用户,等级四共532名用户,等级五共1761名用户,一共3703名用户。
Table 4. Number of users result
表4. 用户数目结果
等级一 |
等级二 |
等级三 |
等级四 |
等级五 |
220 |
410 |
780 |
532 |
1761 |
4.3. 用户风险等级辨识结果
3703名用户核心部件损伤作为输入,5个用户类别作为输出,其中3000名用户作为训练集,703名用户作为测试集,如图5所示,预测集的准确率达到100%,如图6所示,测试集的准确率达到96.9697%。
Figure 5. Prediction set accuracy
图5. 预测集准确率
Figure 6. Test set accuracy
图6. 测试集准确率
5. 结论
本研究通过对电动汽车云端用户数据的分析,深入探讨了电驱动系统中轴、齿轮、轴承等关键部件的失效机理与寿命模型,结合用户驾驶行为,评估了不同用户群体中各个部件的损伤情况与失效风险。研究的主要结论包括:
1) 失效机理与寿命模型的构建:基于对电驱动系统中轴、齿轮、轴承部件失效机理的分析,本研究建立了相应的寿命模型。通过这些模型,我们能够准确预测各部件的损伤程度,并为后续的维护与风险管理提供依据。
2) 风险等级评估与分类:本研究提出了基于不同等级用户的失效风险评估模型。通过对同一等级用户内不同部件的失效风险进行评估,发现即使在同一用户类别内,不同部件的风险等级也存在差异。此外,针对不同等级用户的同一部件失效风险等级评估,有助于精确识别高风险用户,进一步提升风险管理的精度。
3) 风险用户等级的定义与应用:研究定义了“风险用户”这一概念,并基于风险用户等级训练了随机森林网络模型。该模型能够有效识别潜在的高风险用户,并为电动汽车的风险管理提供量化的支持。通过这一方法,企业能够更有针对性地进行预防性维护与维修,减少意外故障的发生。
通过以上结论可以看出,本文的研究为电动汽车行业提供了新的思路和方法,尤其是在提高电驱动系统可靠性、优化维护策略、降低运营成本等方面具有重要意义。未来的研究可以进一步深入探讨不同驾驶环境、驾驶模式以及电动汽车使用场景对电驱动系统可靠性的影响,推动智能化和个性化维护服务的实现。