CN112982544A - 一种基于kpca的地下铲运机故障诊断方法 - Google Patents
一种基于kpca的地下铲运机故障诊断方法 Download PDFInfo
- Publication number
- CN112982544A CN112982544A CN202110215155.7A CN202110215155A CN112982544A CN 112982544 A CN112982544 A CN 112982544A CN 202110215155 A CN202110215155 A CN 202110215155A CN 112982544 A CN112982544 A CN 112982544A
- Authority
- CN
- China
- Prior art keywords
- data
- matrix
- statistics
- spe
- kpca
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000003745 diagnosis Methods 0.000 title claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 20
- 239000013598 vector Substances 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 44
- 239000010705 motor oil Substances 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 239000010720 hydraulic oil Substances 0.000 claims description 5
- 239000003921 oil Substances 0.000 claims description 5
- 230000000717 retained effect Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 239000002826 coolant Substances 0.000 claims description 2
- 239000000446 fuel Substances 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 abstract description 2
- 229910052802 copper Inorganic materials 0.000 abstract description 2
- 239000010949 copper Substances 0.000 abstract description 2
- 238000012795 verification Methods 0.000 abstract description 2
- 230000002159 abnormal effect Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000110 cooling liquid Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02F—DREDGING; SOIL-SHIFTING
- E02F9/00—Component parts of dredgers or soil-shifting machines, not restricted to one of the kinds covered by groups E02F3/00 - E02F7/00
- E02F9/26—Indicating devices
- E02F9/267—Diagnosing or detecting failure of vehicles
- E02F9/268—Diagnosing or detecting failure of vehicles with failure correction follow-up actions
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02F—DREDGING; SOIL-SHIFTING
- E02F9/00—Component parts of dredgers or soil-shifting machines, not restricted to one of the kinds covered by groups E02F3/00 - E02F7/00
- E02F9/20—Drives; Control devices
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Civil Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明提供一种基于KPCA的地下铲运机故障诊断方法,属于地下铲运机故障诊断技术领域。该方法利用地下铲运机实际生产过程中产生的历史数据搭建KPCA模型,通过对特征值分解,进一步求得特征向量来建立主元子空间与残差子空间的T2与SPE统计量,并利用建模数据与两个统计量的分布得到各自的控制限,并与当前检测数据的统计量比较,判断系统是否出现异常,完成故障检测,计算各参量对统计量的贡献值,得到贡献率图,定位故障变量。本发明采用谦比希铜矿的地下铲运机实际生产数据进行案例验证,试验结果表明本发明所提出的方法是十分有效的。
Description
技术领域
本发明涉及地下铲运机故障诊断技术领域,特别是指一种基于KPCA的地下铲运机故障诊断方法。
背景技术
矿山装备地下铲运机运用于矿井内出矿作业,工作环境恶劣,工况复杂,受颠簸,粉尘,潮湿等影响容易发生故障,致使装备长时间停机维修甚至堵塞巷道,常造成大量的经济损失。因此,地下铲运机的故障诊断十分重要,对地下铲运机故障的及时检测和准确定位能大幅度减少维修时间,提高装备工作效率。
对于故障诊断问题,传统方法根据原理不同,分为基于解析模型的故障诊断方法、基于知识的故障诊断方法,这两种方法偏重于局部遥信量信息,难以满足故障诊断算法对全局性与快速性的要求,并缺乏对铲运机变量间耦合性强且具有强非线性特点的考虑。同时随着系统结构日趋复杂,导致此类算法所需的完备知识库和样本更加难以准确获取,故障诊断精度难以保证。地下铲运机故障特征关联复杂,变量间具有很强的非线性,且噪声等干扰情况实时存在,采用上述算法可能会引起故障检测发生错误。
发明内容
本发明要解决的技术问题是提供一种基于KPCA的地下铲运机故障诊断方法。
该方法首先调用地下铲运机历史正常运行状态数据对模型进行训练,确定T2统计量和SPE统计量的置信限,然后再标准化处理采集到的实时数据,计算特征空间中的T2统计量和残差空间中的SPE统计量是否超过建模过程中计算得到的置信限,若上述统计量超过置信限,则进一步计算测试数据对统计量的贡献,判断故障发生位置,完成故障定位。
具体包括步骤如下:
(1)调用地下铲运机历史正常运行状态数据对KPCA铲运机故障检测模型进行训练,确定T2统计量和SPE统计量的置信限;其中,KPCA铲运机故障检测模型通过历史正常运行数据及KPCA算法构建得到;
(2)标准化处理采集到的铲运机实际运行时的数据,计算该数据通过分解得到的主元矩阵空间,即特征空间中的T2统计量和该数据通过分解得到的残差矩阵空间,即残差空间中的SPE统计量是否超过建模过程中计算得到的置信限;
(3)若步骤(2)所得统计量超过置信限,则进一步计算测试数据对统计量的贡献,判断故障发生位置,完成故障定位。
其中,步骤(1)中运行状态数据包括发动机油温、发动机油压、发动机进气歧管压力、发动机进气歧管温度、发动机冷却液温度、发动机转速、发动机扭矩、发动机燃油率、发动机负载、大臂与铲斗压力、转向泵压力、变速箱油温、变速箱油压、液压油温度、前轴制动压力、后轴制动压力、制动液压油温度、制动回路气压、右前轮胎压力、油门踏板位置、车速共21维参数。
步骤(1)中确定T2统计量和SPE统计量的置信限的具体步骤如下:
S1:选取地下铲运机正常状况下的数据,通过计算各变量的均值和标准差对数据进行标准化处理,
特征空间下的数据协方差矩阵为
协方差矩阵的特征方程可表示为
Cv=λv
其中,协方差矩阵C的特征值λ>0,特征向量v(v∈Rn)为单位向量;
由上式所获得的最大特征值λ所对应的特征向量v成为特征空间F中第一个主成分(PC),特征向量vk可表示为:
将特征向量vk代入内积运算式并化简得
nλα=Kα
上述,nλ为核矩阵K的特征值,α为特征值nλ对应的特征向量,σ为选定的参数;
S4:通过nλα=Kα计算特征值nλ,规范化特征向量αk使得<αk,αk>=1/λk;
其中,nλ为核矩阵K的特征值,α为特征值nλ对应的特征向量;
其中,ti为v在特征空间上的投影,αi,j为K的第i个特征值对αi的第j个系数,x为原始空间的输入向量;
S6:由上述过程,建立KPCA模型;
S7:计算正常运行数据的T2和SPE统计量并确定置信限;
T2统计量:T2=[t1,t2,…,tk]Λ-1[t1,t2,…,tk]T,其中k表示保留的主元个数,Λ-1表示与保留的主元相对应的特征值对角矩阵的逆矩阵;
步骤(2)中当T2计量满足T2>Tα 2或SPE统计量满足SPE>SPEα至少有一个成立时,表明系统可能发生故障,反之则为正常运行的状态。其中,T为实时数据T统计量,Tα为T统计量的置信限,SPE为实时数据SPE统计量,SPEα为SPE统计量置信限。
步骤(3)中具体步骤包括:
S1:计算T2统计量的贡献
用CT2表示T2对vi的偏导数,进一步刻画各个变量对T2的贡献程度:
其中第二项和第四项对ytest不影响,作为常数项处理,变换成如下公式,为第i个变量对T2统计量的贡献,当某一个变量对应的CT2为最大值时,表明相应变量为故障变量:
本发明的上述技术方案的有益效果如下:
上述方案中,采用谦比希铜矿的地下铲运机实际生产数据进行案例验证,试验结果表明本发明所提出的方法是十分有效的。
附图说明
图1为本发明方法的流程图;
图2为本发明实施例提供的KPCA主元贡献率为85%时的诊断结果图;
图3为本发明实施例提供的KPCA主元贡献率为95%时的诊断结果图;
图4为本发明实施例提供的用于故障变量定位的贡献图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明提供一种基于KPCA的地下铲运机故障诊断方法。
如图1所示,本方法首先调用地下铲运机历史正常运行状态数据对模型进行训练,确定T2统计量和SPE统计量的置信限,然后再标准化处理采集到的实时数据,计算特征空间中的T2统计量和残差空间中的SPE统计量是否超过建模过程中计算得到的置信限,若上述统计量超过置信限,则进一步计算测试数据对统计量的贡献,判断故障发生位置,完成故障定位。
下面结合具体实施例予以说明。
实施例1
基于KPCA算法的地下铲运机故障模型在建模过程中,需要采集地下铲运机在稳定运行条件下的多变量测量数据,如发动机油温、发动机油压、发动机冷却液温度、变速箱油温、变速箱油压、液压油温等由相应传感器获取的信息。这些采集的测量数据构成数据矩阵X(x1,x2,…,xn)∈Rn×m,n表示样本数量,m表示特征数量。
建立基于KPCA算法的地下铲运机故障诊断模型,需要经过以下三个阶段:
步骤1:模型训练。调用地下铲运机历史正常运行状态数据对模型进行训练,确定T2统计量和SPE统计量的置信限。本实施例选取地下铲运机正常运行时的4000个数据点作为训练数据,如表1所示:
表1:训练数据
步骤1.1:选取地下铲运机正常状况下的数据,通过计算各变量的均值和标准差对数据进行标准化处理。
步骤1.4:通过nλα=Kα计算特征值nλ,规范化特征向量αk使得<αk,αk>=1/λk;
步骤1.6:建立KPCA模型;
步骤1.7:计算正常运行数据的T2和SPE统计量并确定置信限。
T2统计量:T2=[t1,t2,…,tk]Λ-1[t1,t2,…,tk]T,其中k表示保留的主元个数,Λ-1表示与保留的主元相对应的特征值对角矩阵的逆矩阵。
步骤2:在线监测,标准化处理采集到的实时数据,计算特征空间中的T2统计量和残差空间中的SPE统计量是否超过建模过程中计算得到的置信限。
本实施例将地下铲运机运行过程中采集的出现柴油污染导致牵引力低故障的实际数据作为测试数据,如表2所示,该测试数据包括100组故障数据,共500组数据,其中故障出现在第200-300组数据。
表2:测试数据
本实施例将主元贡献率选择为85%、90%分别进行如上实验,所得结果如图2及图3所示。由该结果分析可知此方法可准确识别出地下铲运机故障节点,即此方法对地下铲运机故障诊断具有很好的适用性,但准确的进行故障定位还需要进行进一步分析。
步骤3:故障变量定位。如上述步骤2所得统计量在200-300组数据点超过置信限,则进一步计算测试数据对统计量的贡献,判断故障发生位置,完成故障定位。
步骤3.1:T2统计量的贡献
用CT2表示T2对vi的偏导数,进一步刻画各个变量对T2的贡献程度。
其中第二项和第四项对ytest不影响,作为常数项处理,变换成如下公式,为第i个变量对T2统计量的贡献,当某一个变量对应的CT2为最大值时,表明相应变量为故障变量。
通过以上分析,得到该故障的故障变量如图4所示,由该贡献图可知该故障的故障变量为变量3,其中变量3对应为发动机油压,符合实际机理分析结果,由此该故障诊断方法对地下铲运机具有较好的适用性。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种基于KPCA的地下铲运机故障诊断方法,其特征在于:包括步骤如下:
(1)调用地下铲运机历史正常运行状态数据对KPCA铲运机故障检测模型进行训练,确定T2统计量和SPE统计量的置信限;其中,KPCA铲运机故障检测模型通过历史正常运行数据及KPCA算法构建得到;
(2)标准化处理采集到的铲运机实际运行时的数据,计算该数据通过分解得到的主元矩阵空间,即特征空间中的T2统计量和该数据通过分解得到的残差矩阵空间,即残差空间中的SPE统计量是否超过建模过程中计算得到的置信限;
(3)若步骤(2)所得统计量超过置信限,则进一步计算测试数据对统计量的贡献,判断故障发生位置,完成故障定位。
2.根据权利要求1所述的基于KPCA的地下铲运机故障诊断方法,其特征在于:所述步骤(1)中运行状态数据包括发动机油温、发动机油压、发动机进气歧管压力、发动机进气歧管温度、发动机冷却液温度、发动机转速、发动机扭矩、发动机燃油率、发动机负载、大臂与铲斗压力、转向泵压力、变速箱油温、变速箱油压、液压油温度、前轴制动压力、后轴制动压力、制动液压油温度、制动回路气压、右前轮胎压力、油门踏板位置、车速共21维参数。
3.根据权利要求1所述的基于KPCA的地下铲运机故障诊断方法,其特征在于:所述步骤(1)中确定T2统计量和SPE统计量的置信限的具体步骤如下:
S1:选取地下铲运机正常状况下的数据,通过计算各变量的均值和标准差对数据进行标准化处理,
特征空间下的数据协方差矩阵为
协方差矩阵的特征方程可表示为
Cv=λv
其中,协方差矩阵C的特征值λ>0,特征向量v,(v∈Rn为单位向量;
由上式所获得的最大特征值λ所对应的特征向量v成为特征空间F中第一个主成分PC,特征向量vk表示为:
将特征向量vk代入内积运算式并化简得
nλα=Kα
上述,nλ为核矩阵K的特征值,α为特征值nλ对应的特征向量,σ为选定的参数;
S4:通过nλα=Kα计算特征值nλ,规范化特征向量αk使得<αk,αk>=1/λk;
其中,nλ为核矩阵K的特征值,α为特征值nλ对应的特征向量;
其中,ti为v在特征空间上的投影,αi,j为K的第i个特征值对αi的第j个系数,x为原始空间的输入向量;
S6:由上述过程,建立KPCA模型;
S7:计算正常运行数据的T2和SPE统计量并确定置信限;
T2统计量:T2=[t1,t2,…,tk]Λ-1[t1,t2,…,tk]T,其中k表示保留的主元个数,Λ-1表示与保留的主元相对应的特征值对角矩阵的逆矩阵;
4.根据权利要求1所述的基于KPCA的地下铲运机故障诊断方法,其特征在于:所述步骤(2)中当T2计量满足T2>Tα 2或SPE统计量满足SPE>SPEα至少有一个成立时,表明系统可能发生故障,反之则为正常运行的状态,其中,T为实时数据T统计量,Tα为T统计量的置信限,SPE为实时数据SPE统计量,SPEα为SPE统计量置信限。
5.根据权利要求1所述的基于KPCA的地下铲运机故障诊断方法,其特征在于:所述步骤(3)中具体步骤包括:
S1:计算T2统计量的贡献
用CT2表示T2对vi的偏导数,进一步刻画各个变量对T2的贡献程度:
其中第二项和第四项对ytest不影响,作为常数项处理,变换成如下公式,为第i个变量对T2统计量的贡献,当某一个变量对应的CT2为最大值时,表明相应变量为故障变量:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110215155.7A CN112982544B (zh) | 2021-02-24 | 2021-02-24 | 一种基于kpca的地下铲运机故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110215155.7A CN112982544B (zh) | 2021-02-24 | 2021-02-24 | 一种基于kpca的地下铲运机故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112982544A true CN112982544A (zh) | 2021-06-18 |
CN112982544B CN112982544B (zh) | 2022-03-29 |
Family
ID=76350959
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110215155.7A Active CN112982544B (zh) | 2021-02-24 | 2021-02-24 | 一种基于kpca的地下铲运机故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112982544B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120232756A1 (en) * | 2011-03-03 | 2012-09-13 | Eaton Corporation | Fault detection, isolation and reconfiguration systems and methods for controlling electrohydraulic systems used in construction equipment |
US20130035910A1 (en) * | 2010-09-29 | 2013-02-07 | Yingwei Zhang | Continuous annealing process fault detection method based on recursive kernel principal component analysis |
CN103422525A (zh) * | 2013-08-24 | 2013-12-04 | 烟台兴业机械股份有限公司 | 地下铲运机的故障智能诊断监视系统 |
-
2021
- 2021-02-24 CN CN202110215155.7A patent/CN112982544B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130035910A1 (en) * | 2010-09-29 | 2013-02-07 | Yingwei Zhang | Continuous annealing process fault detection method based on recursive kernel principal component analysis |
US20120232756A1 (en) * | 2011-03-03 | 2012-09-13 | Eaton Corporation | Fault detection, isolation and reconfiguration systems and methods for controlling electrohydraulic systems used in construction equipment |
CN103422525A (zh) * | 2013-08-24 | 2013-12-04 | 烟台兴业机械股份有限公司 | 地下铲运机的故障智能诊断监视系统 |
Non-Patent Citations (1)
Title |
---|
董佳: "基于主元分析的多故障状态监测与诊断方法研究", 《信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112982544B (zh) | 2022-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110298455B (zh) | 一种基于多变量估计预测的机械设备故障智能预警方法 | |
CN109710983B (zh) | 一种基于关键性能指标的柴油机气缸分层故障诊断方法 | |
CN104714537B (zh) | 一种基于联合相对变化分析和自回归模型的故障预测方法 | |
CN104699077B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 | |
CN104777830B (zh) | 一种基于kpca混合模型的多工况过程监控方法 | |
CN112685910B (zh) | 一种基于混合预测模型的复杂装备电源组故障预测方法 | |
CN102736546A (zh) | 一种流程工业复杂机电系统的状态监测装置及方法 | |
CN104793606B (zh) | 基于改进的kpca和隐马尔科夫模型的工业故障诊断方法 | |
CN108549908B (zh) | 基于多采样概率核主成分模型的化工过程故障检测方法 | |
CN107451392A (zh) | 一种含有多个相关退化过程的剩余寿命预测方法 | |
CN109187060B (zh) | 列车速度传感器信号异常检测及轴抱死故障诊断方法 | |
CN106907927A (zh) | 一种核灵活流形嵌入电熔镁炉故障监测方法 | |
CN108830006B (zh) | 基于线性评价因子的线性-非线性工业过程故障检测方法 | |
CN113569990A (zh) | 面向强噪声干扰环境的演艺装备故障诊断模型构建方法 | |
CN112982544B (zh) | 一种基于kpca的地下铲运机故障诊断方法 | |
CN112149054B (zh) | 基于时序扩展的正交邻域保持嵌入模型的构建与应用 | |
CN109001492A (zh) | 一种列车轴速度传感器故障诊断方法 | |
CN116894475A (zh) | 一种基于神经网络的燃气轮机状态预警方法 | |
CN111832442A (zh) | 从海量桥梁动态应变数据中自动化分离温度应变成分的方法 | |
CN116933003A (zh) | 基于DaNet的无人机发动机的剩余使用寿命预测方法 | |
CN114995338A (zh) | 一种基于规范变量分析与js散度融合的工业过程微小故障检测方法 | |
CN114371677A (zh) | 基于谱半径-区间主成分分析的工业过程状态监测方法 | |
CN110288724B (zh) | 一种基于小波函数主元分析的批次过程监测方法 | |
CN114200914A (zh) | 一种基于mw-occa的质量相关早期故障检测方法 | |
CN114266401B (zh) | 一种底盘发动机故障预测方法及系统 |
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 |