CN104914846B - 基于自适应滑模观测器的电连接器间歇性失效检测方法 - Google Patents
基于自适应滑模观测器的电连接器间歇性失效检测方法 Download PDFInfo
- Publication number
- CN104914846B CN104914846B CN201510153226.XA CN201510153226A CN104914846B CN 104914846 B CN104914846 B CN 104914846B CN 201510153226 A CN201510153226 A CN 201510153226A CN 104914846 B CN104914846 B CN 104914846B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msup
- mtr
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0224—Process history based detection method, e.g. whereby history implies the availability of large amounts of data
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- External Artificial Organs (AREA)
Abstract
本发明公开了一种基于自适应滑模观测器的电连接器间歇性失效检测方法,包括建立数学模型;基于数学模型设计滑模观测器,在滑模项中对滑模观测器增益变量项设计自适应律以避免故障上界未知对状态估计的影响;对滑模观测器增益常数项求取值范围以使系统能在有限时间内达到并维持滑模运动;利用等效误差注入技术进行故障重构。该故障检测方法只需故障本身有界,对于故障的一阶导数没有要求,对于缓变、突变、间歇故障均可有效检测。当电连接器接触电阻发间歇性突变时可以被有效检测出,通过设置阈值可以判断出电连接器是否失效。本发明能够有效地根据已被故障影响的过程变量的测量值估计出故障值,反映出故障的本来面貌,从而达到间歇性失效检测的目的。
Description
技术领域
本发明公开了一种基于自适应滑模观测器的电连接器间歇性失效检测方法,涉及到乘性间歇性故障的检测技术,属于故障检测技术领域。
背景技术
连接器一般用于实现电信号的传输和控制以及电子与电气设备之间的电连接,连接器一旦失效,不仅会造成其提供的传输功能故障,可能还会对其他与其相连的部件造成影响。连接器的常见失效模式有短路、接触不良、瞬断、绝缘不良、断路、误配线、固定不良、密封不良等。接触不良、瞬断失效可以统称为间歇性失效(Intermittent Failure,IF),是失效模式及影响分析、可靠性以及故障检测和诊断领域研究的难点,如何检测失效、确定失效位置以及对系统的影响仍是一个开放性的问题。
具体到连接器的间歇性失效,主要原因为使用过程中由于环境因素的影响造成的接触对间歇地接触失效。电连接器在贮存和使用中都要受到各种环境应力的影响,当电连接器处于振动、冲击或温度变化等环境时,接触对之间将有微动产生。微动发生时,接触表面微动磨损和微动腐蚀大多同时存在。接触对的间歇性失效主要由电连接器的微动磨损引起,电连接器的微动磨损会导致接触电阻逐渐增加,导致发热增加,进而促使金属接触面氧化加速,随后将伴随接触电阻从几毫欧间歇性地增大减小至几欧姆或几十欧姆,最终引发电气系统回路中发生间歇性的不可复现性的失效。
对于电连接器的间歇性失效可以转换为间歇性故障的检测问题。已有的间歇性故障检测方法大都是传统故障检测方法的改进。故障检测主要分为定性和定量两大方法,其中定量方法主要包括基于模型的方法和基于数据的方法两大类。基于观测器的故障检测方法则是基于解析模型的故障检测方法中使用最广泛的。传统的基于观测器的故障检测方法的基本思想是利用系统的定量模型和测量信号重建某一可测变量,将估计值与测量值之差作为残差,以检测和分离系统故障。检测中人为设定阈值的大小会对被检测系统的灵敏度和鲁棒性产生很大的影响,对于慢变或初始故障的检测能力很低。残差信号本身并不能直接反映故障的真实状态,因此基于残差的故障检测与分离技术是一种间接的故障诊断技术。与之相对应,通过故障重构来进行检测不是通过残差信号得到失效信息,而是利用观测器直接重构失效信号,避免了产生和评价残差信号的复杂性,且通过重构不仅能够实现故障检测和分离,还可以估计其严重程度、辨别其类型。
Wesley G等基于特征分析的方法对控制器与永磁交流电机之间的连接器间歇性失效提出了一种检测算法,通过测量交流电机磁场定向电流,然后对其进行快速傅里叶变换或非抽取离散小波变换,便可以得到相应的变换系数,最后根据系统正常运行下的变换系数设置能量阈值即可实现间歇故障的检测。但该方法具有以下不足:首先,上述故障检测方法需要分别获得系统正常情况和发生故障时的数据。在实际电路系统中,由于技术条件限制有些量难以测量。其次,本算法主要是离线情况进行设计的,不能在线实时检测和诊断间歇故障。再者,电流只是间接地体现了电路中电阻的变化,不能直观地反应故障的类型。Bennett S M等研究了存在间歇性传感器故障时牵引电机输出转矩的估计问题,提出了一种奉献观测器的方法来检测间歇故障,但是,此方法没有考虑间歇故障的持续时间和间隔时间对可检测性的影响,当故障间隔时间较短时不再适用。赵瑾等就乘性故障的检测问题以滑模观测器为基础,利用状态和输出变换方法与奇异值分解方法相结合对系统进行降阶,提出一种鲁棒故障重构观测器;同时给出优化滑模策略,并进行了严格论证,保证对系统不确定性具有鲁棒性以及跟踪系统状态的收敛性;应用等价输出控制概念以及设计的鲁棒故障重构观测器.获取故障信息,实现执行器故障的检测与重构。但文中假设故障上界是已知的,并未考虑乘性故障因与状态及系统输入有关上界未知这一情况。其次,该方法限制条件太多,是就方法寻找适用的对象。再者,文中实际仿真对象并未体现出乘性故障的特点。
基于模型的故障检测方法能够在线实时地进行电连接器失效的检测,对于整个系统安全可靠地运行很有意义。能够及时对失效部件进行更换避免事故的发生。通过电连接器间歇性失效的检测可以对电连接器的使用寿命做出评估,为预防事故的发生提供了条件。通过故障重构可以辨别其类型,从而为系统采取相应的措施有效地消除故障对系统的影响提供更充实的依据。
发明内容
本发明所要解决的技术问题是:针对现有技术的缺陷,提供基于自适应滑模观测器的电连接器间歇性失效检测方法,通过设计滑模观测器进行故障重构可以有效地判断出电连接器是否发生间歇性失效。对滑模控制器增益设计自适应律可以放松一般故障重构中所需的故障有界且上界已知这一条件。本方法即使在存在外界扰动的条件下也可以有效地检测电连接器的间歇性失效,且具有良好的移植性,可用于相似系统的故障检测。
本发明为解决上述技术问题采用以下技术方案:
一种基于自适应滑模观测器的电连接器间歇性失效检测方法,实现步骤包括:
步骤一、对电连接器所在系统建立数学模型;
步骤二、根据步骤一中的数学模型,将电连接器的间歇性失效检测等效于数学模型中的乘性故障检测;
步骤三、针对数学模型中乘性故障上界故障未知的特点,对滑模控制器增益设计自适应律;
步骤四、使用滑模中的等效误差注入技术进行误差重构,结合线性变换Lyapunov函数求解滑模控制器增益中常数项的取值;
步骤五、利用故障重构思想重构故障,根据重构出的接触电阻的曲线设定阈值进行失效判断,实现间歇性故障检测。
作为本发明的进一步优选方案,步骤一中,所述数学模型的具体表达如下:
式中,t是系统运行时间;x(t)∈Rn是系统状态向量;是系统状态向量的一阶导数;u(t)∈Rm是控制输入向量;y(t)∈Rp是测量输出向量;ξ(t)为扰动;M∈Rn×q;A∈Rn ×n;B∈Rn×m;C∈Rp×n;Δ(t,x)∈Rq;Rm、Rn、Rp、Rq分别表示m维、n维、p维、q维实向量空间;Rn×n、Rn×m、Rp×n、Rn×q分别表示n×n维、n×m维、p×n维、n×q维实矩阵;
A、B、C、E的表达式分别为:
C=[0 1];
式中,R为电阻值;ΔA代表系统参数变化的乘性故障;ΔR为电连接器的接触电阻的变化;J为电动机轴上的等效转动惯量;L为电枢回路电感;b为轴上的等效阻尼系数;km为电动机的力矩系数;ke是电动机的电势系数;
fi,i=1,2,…,q是待重构的乘性故障,为已知的与矩阵A和状态x有关的函数。
作为本发明的进一步优选方案,所述步骤三的具体过程为:
设计自适应滑模观测器:
其中,分别为x,y的估计向量,H∈Rn×p为待定的增益矩阵,υ(t)是确保在滑动平面存在独立于ξ(t)的稳定滑模运动的优化滑模策略,其表达式如下所示:
针对故障上界未知的情况,对参数ρ1(t)设计如下所示的自适应律:
其中,ρ0,Γ为大于0的常数。
作为本发明的进一步优选方案,所述步骤四的具体过程为:
设定系统估计误差方程如下所示:
令ε=β-ρ1,构造Lyapunov函数:
V(t)=eTPe+Γ-1ε2;
求取滑模控制器的增益:
取使得滑模运动达到理想的滑动模态的ρ0值,其中,是CT的正交补;
T-1=[ΣT(ΣPTΣT)-1 P-1CT(CP-1CT)-1];
其中,
对状态误差方程做线性变换,即则得到新坐标下的状态误差方程为:
即:
作为本发明的进一步优选方案,所述步骤五的具体过程为:
采用等价输出误差介入原理,当滑模运动到达时,有则系统误差方程变为:
式中,A11=ΣP(A-HC)ΣT(ΣPΣT)-1,取Q0=ΣPΣT;
当计算式:成立时,则A11稳定,得到重构信号为:
作为本发明的进一步优选方案,所述步骤五还包括:采用光滑函数进行信号调整,将不连续的控制输入υ1用连续的方式来代替:
重构故障设定定义重构的乘性故障:
其中,δ'为正常数参数,当时,将精确逼近f。
作为本发明的进一步优选方案,所述参数δ'的取值为,δ'=10-5。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
⑴基于自适应滑模观测器的电连接器间歇性失效检测方法。可以利用观测器直接重构失效信号,避免了产生和评价残差信号的复杂性,降低了故障误报率和漏报率,且可以根据重构信号设置阈值直观地判断出故障类型。
⑵通过滑模的等效误差注入技术进行故障重构不需要故障一阶导数有界,只需故障有界即可,所以对于突变信号也可以有效检测。
⑶通过对滑模观测器增益项设计自适应律可以消除故障上界未知的影响,进一步求取滑模观测器增益常数项可以使滑模运动快速到达并维持在滑模面。
附图说明
图1为本发明实施流程;
图2为含有电连接器的系统实施例电路图;
图3为电连接器的接触电阻变化曲线;
图4为本发明间歇性失效检测方法效果仿真,其中:
图4(a)为状态跟踪曲线,图4(b)重构接触电阻示意图,图4(c)为重构误差示意图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本发明涉及的故障检测技术首先获取系统的数学模型;其次基于解析模型设计滑模观测器,在滑模项中对滑模观测器增益设计自适应律以避免故障上界未知对状态估计的影响。在稳定观测器设计的基础上,对滑模观测器增益常数项规定取值范围,以保证观测器能在有效时间达到滑动模态且维持在滑动模态。从而可以利用等效误差注入技术进行故障重构。该故障检测方法只需故障本身有界,对于故障的一阶导数没有要求,对于缓变、突变、间歇故障均可有效检测。因此当电连接器接触电阻发间歇性突变时可以被有效检测出,从而判断出电连接器是否失效。
下面结合附图对本发明的技术方案做进一步的详细说明:
⑴系统数学模型的建立
图2中,r1为电连接器前端导线电阻,r2为电连接器后端导线电阻,Ra、L分别为电枢回路的电阻和电感,R电为电连接器的接触电阻,e是电动机的反电势,J为电动机轴上的等效转动惯量,b为轴上的等效阻尼系数;根据基尔霍夫定律和牛顿第二定律对图2所示的系统列基本方程组:
式中,ke是电动机的电势系数,km为电动机的力矩系数;Td为扰动力矩;为电机轴的转角θ的一阶导数,即电机轴的角速度;i为电路中的电流;为了将该系统模型转换成便于检测研究的数学模型表达式,选取则根据系统模型(1)建立系统的状态空间表达式为
式中,t是系统运行时间,在仿真中t取为10s;x(t)∈Rn是系统状态向量;是系统状态向量的一阶导数;u(t)∈Rm是控制输入向量;y(t)∈Rp是测量输出向量;ξ(t)为扰动;A∈Rn×n;B∈Rn×m;C∈Rp×n;ΔA∈Rn×n;Rm、Rp分别表示m维和n维实向量空间;Rn×n、Rn×m、Rp×n分别表示n×n维、n×m维、p×n维实矩阵;A、B、C、E、ΔA的表达式如下:
C=[0 1];
式中,R=r1+R电+r2+Ra;ΔA代表系统参数变化的乘性故障;ΔR为电连接器的接触电阻的变化。大部分基于模型的故障检测均是针对加性故障进行的,对系统(2)调整可以重新写成下列形式:
式中M∈Rn×q,Δ(t,x)∈Rq;表达式如下:
其中,fi,i=1,2,…,q是待重构的乘性故障,为已知的与矩阵A和状态x有关的函数,Rq表示q维实向量空间。
对于系统(3),系统可观,因此可以进行观测器的设计。由于rank(CE)≠rank(E),将E拆为E=E1+M;
存在适当维数的正定矩阵P和Q以及矩阵H满足式:
(A-HC)TP+P(A-HC)=-Q (4)
存在非零矩阵G、F满足:PM=CTGT;PE1=CTFT;
扰动ξ(t)范数有界,即||ξ(t)||≤α,且α已知;故障Δ(t)范数有界,即||Δ(t,x)||≤β,且β未知。
(2)自适应滑模观测器的设计
设计如下所示的自适应滑模观测器:
其中,分别为x,y的估计向量,H∈Rn×p为待定的增益矩阵,υ(t)是确保在滑动平面存在独立于ξ(t)的稳定滑模运动的优化滑模策略,其表达式如下所示:
针对故障上界未知的情况,对参数ρ1(t)设计如下所示的自适应律:
其中,ρ0,Γ为大于0的常数,本发明中取Γ为0.01。
由系统(3)和(5)得到系统估计误差方程如下所示:
令ε=β-ρ1,构造如下的Lyapunov函数:
V(t)=eTPe+Γ-1ε2 (9)
则有:
根据式(4)(6)(7)可得:
因为对任意正定矩阵P∈Rn×n和e∈Rn,有λm(P)||e||2≤eTPe≤λM(P)||e||2,其中λm和λM分别表示矩阵P的最小和最大特征值,则:
当,时,所以设计的自适应滑模观测器对系统扰动具有鲁棒性,同时具有一致最终有界的构造误差。
(3)滑模控制器增益的求取
为了求得能使滑模运动能在有限时间内达到理想的滑动模态的ρ0值,取其中是CT的正交补,则T-1=[ΣT(ΣPTΣT)-1 P-1CT(CP-1CT)-1],其中则对状态误差方程做线性变换,即则由式(8)得:
令:
则得到新坐标下的状态误差方程为:
即:
取Lyapunov函数:
V2(t)=ey TWey+Γ-1ε2 (14)
式中,W=(CP-1CT)-1,则:
取则
由于
则
从e1(t)进入Φ的时刻到滑模运动发生的时刻进行积分,有
式中,V2(tΦ)为在t=tΦ时刻V2的初始条件,ts是滑模运动开始的时刻。可以证明达到滑模运动所需的时间ts满足:
因此式(13)所示误差系统将在有限时间内到达滑模面S={e:Ce=0},并保持滑模运动。
(4)故障重构
采用等价输出误差介入原理,当滑模运动到达时,有系统误差方程(13)变为:
式中,A11=ΣP(A-HC)ΣT(ΣPΣT)-1,取Q0=ΣPΣT,结合式(4),下式成立
则A11稳定。
重构信号为:
由式(17)可以看出重构故障信号的精确度与干扰有关,但通常干扰要远远小于故障,所以可以采用式(18)近似故障重构。滑模变结构控制会出现抖振现象,为了削弱抖振,可以采用光滑函数进行信号调整,将不连续的控制输入υ1用一种连续的方式来代替:
其中,δ是一个很小的正值,可以看出只要正确的选择δ值,等效输出注入就可以用式(19)替代,本发明中取δ=0.05。根据式(18)可以重构故障因为设计的自适应滑模观测器可以渐近估计系统状态x,所以定义重构的乘性故障:
其中,δ'为很小的正常数,若则能够较精确的逼近f,本发明中取δ'=10-5。
规定接触电阻超过10Ω时发生间歇性失效;图3中有两段持续不同时间和不同幅值的间歇性失效区。如能精确重构,便能判断出电连接器是否发生间歇性失效。图1描述了本发明的实施流程,图4则是采用本发明的自适应滑模观测器的航空电连接器间歇性失效检测的检测效果。图4(a)为所涉及的自适应滑模观测器在干扰存在条件下的状态跟踪曲线,由图可以看出自适应滑模观测器估计的状态可以很好的跟踪原系统的状态,这说明了本文设计的观测器可以克服故障和干扰的影响。为了达到失效检测的目的,在状态估计收敛的情况下对电连接器的接触电阻进行重构,图4(b)和(c)分别为重构电阻曲线和重构误差曲线,从图中可以看出经过有限时间后,重构信号能够跟踪原信号。在滑模运动达到时其跟踪误差趋向于零,但有一些误差值较大的突变点,将其与原接触电阻变化曲线相比,可以发现这些点正是接触电阻发生突变的点,因为在突变点处状态跟踪误差较大,导致重构信号误差也较大,但很快状态就收敛了。总的来说,利用本文的自适应滑模观测器可以很好地进行失效检测,在有干扰存在的条件下,对突变和缓变信号均能快速地跟踪,从而便于故障类型的判断,进一步可以知道系统部件是否失效。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案内容,依据本发明的技术实质,在本发明的精神和原则之内,对以上实施例所作的任何简单的修改、等同替换与改进等,均仍属于本发明技术方案的保护范围之内。
Claims (7)
1.基于自适应滑模观测器的电连接器间歇性失效检测方法,其特征在于,实现步骤包括:
步骤一、对电连接器所在系统建立数学模型;
步骤二、根据步骤一中的数学模型,将电连接器的间歇性失效检测等效于数学模型中的乘性故障检测;
步骤三、针对数学模型中乘性故障上界未知的特点,对滑模观测器增益设计自适应律;
步骤四、结合线性变换Lyapunov函数求解滑模观测器增益中常数项的取值,以使系统在有限时间内达到并维持滑模运动;
步骤五、利用故障重构思想重构故障,根据重构出的接触电阻的曲线设定阈值进行失效判断,实现间歇性失效检测。
2.如权利要求1所述的基于自适应滑模观测器的电连接器间歇性失效检测方法,
其特征在于,步骤一中,所述数学模型的具体表达如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mi>E</mi>
<mi>&xi;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mi>M</mi>
<mi>&Delta;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mi>C</mi>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中,t是系统运行时间;x(t)∈Rn是系统状态向量;是系统状态向量的一阶导数;u(t)∈Rm是控制输入向量;y(t)∈Rp是测量输出向量;ξ(t)为扰动;M∈Rn×q;A∈Rn×n;B∈Rn×m;C∈Rp×n;Δ(t,x)∈Rq;Rm、Rn、Rp、Rq分别表示m维、n维、p维、q维实向量空间;Rn×n、Rn×m、Rp ×n、Rn×q分别表示n×n维、n×m维、p×n维、n×q维实矩阵;
A、B、C、E的表达式分别为:
C=[0 1];
<mrow>
<mi>M</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<mi>B</mi>
<mo>;</mo>
</mrow>
式中,R为电阻值;ΔA代表系统参数变化的乘性故障;ΔR为电连接器的接触电阻的变化;J为电动机轴上的等效转动惯量;L为电枢回路电感;b为轴上的等效阻尼系数;km为电动机的力矩系数;ke是电动机的电势系数;
fi,i=1,2,…,q是待重构的乘性故障,为已知的与矩阵A和状态x有关的函数。
3.如权利要求2所述的基于自适应滑模观测器的电连接器间歇性失效检测方法,其特征在于,所述步骤三的具体过程为:
设计自适应滑模观测器:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mo>-</mo>
<mi>H</mi>
<mo>(</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<mi>y</mi>
<mo>)</mo>
<mo>+</mo>
<mi>&upsi;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mi>C</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中,分别为x,y的估计向量;H∈Rn×p为待定的增益矩阵;υ(t)是确保在滑动平面存在独立于ξ(t)的稳定滑模运动的优化滑模策略,其表达式如下所示:
<mrow>
<mi>&upsi;</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>P</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>|</mo>
<mo>|</mo>
<mi>G</mi>
<mo>|</mo>
<mo>|</mo>
<mfrac>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>&NotEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
针对故障上界未知的情况,对参数ρ1(t)设计如下所示的自适应律:
<mrow>
<msub>
<mover>
<mi>&rho;</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>&Gamma;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>G</mi>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>;</mo>
</mrow>
式中,ρ0,Γ为大于0的常数;
存在适当维数的正定矩阵P和Q以及矩阵H满足式:
(A-HC)TP+P(A-HC)=-Q (4)
存在非零矩阵G、F满足:PM=CTGT;PE1=C TFT;
扰动ξ(t)范数有界,即||ξ(t)||≤α,且α已知;故障Δ(t)范数有界,即||Δ(t,x)||≤β,且β未知。
4.如权利要求3所述的基于自适应滑模观测器的电连接器间歇性失效检测方法,其特征在于,所述步骤四的具体过程为:
设定系统估计误差方程如下所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>(</mo>
<mi>A</mi>
<mo>-</mo>
<mi>H</mi>
<mi>C</mi>
<mo>)</mo>
<mi>e</mi>
<mo>+</mo>
<mi>&upsi;</mi>
<mo>-</mo>
<mi>E</mi>
<mi>&xi;</mi>
<mo>-</mo>
<mi>M</mi>
<mi>&Delta;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mi>C</mi>
<mi>e</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
令ε=β-ρ1,构造Lyapunov函数:
V(t)=eTPe+Γ-1ε2;
求取滑模控制器的增益:
取使得滑模运动达到理想的滑动模态的ρ0值,其中,是CT的正交补;
T-1=[ΣT(ΣPTΣT)-1P-1CT(CP-1CT)-1];
其中,
对状态误差方程做线性变换,即可得:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>T</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>-</mo>
<mi>H</mi>
<mi>C</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>T</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>&upsi;</mi>
<mo>-</mo>
<mi>E</mi>
<mi>&xi;</mi>
<mo>-</mo>
<mi>M</mi>
<mi>&Delta;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>-</mo>
<mi>H</mi>
<mi>C</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>T</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>T</mi>
<mi>&upsi;</mi>
<mo>-</mo>
<mi>T</mi>
<mi>E</mi>
<mi>&xi;</mi>
<mo>-</mo>
<mi>T</mi>
<mi>M</mi>
<mi>&Delta;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
令:
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>T</mi>
<mi>E</mi>
<mn>1</mn>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>C</mi>
<mo>&perp;</mo>
<mi>T</mi>
</msubsup>
<msup>
<mi>PP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msup>
<mi>F</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msup>
<mi>F</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msup>
<mi>F</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mi>T</mi>
<mi>&upsi;</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>|</mo>
<mo>|</mo>
<mi>G</mi>
<mo>|</mo>
<mo>|</mo>
<mfrac>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&upsi;</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
则得到新坐标下的状态误差方程为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>T</mi>
<mi>&upsi;</mi>
<mo>-</mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mi>&xi;</mi>
<mo>-</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>x</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
即:
5.如权利要求4所述的基于自适应滑模观测器的电连接器间歇性失效检测方法,其特征在于,所述步骤五的具体过程为:
采用等价输出误差介入原理,当滑模运动到达时,有则系统误差方程变为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>11</mn>
</msub>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>13</mn>
</msub>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&upsi;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msup>
<mi>G</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>x</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中,A11=ΣP(A-HC)ΣT(ΣPΣT)-1,取Q0=ΣPΣT;
当计算式:成立时,A11稳定,重构信号为:
<mrow>
<mover>
<mi>&Delta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>C</mi>
<mi>M</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msub>
<mi>&upsi;</mi>
<mn>1</mn>
</msub>
<mo>.</mo>
</mrow>
6.如权利要求5所述的基于自适应滑模观测器的电连接器间歇性失效检测方法,其特征在于,所述步骤五还包括:采用光滑函数进行信号调整,将不连续的控制输入υ1用连续的方式来代替:
<mrow>
<msub>
<mi>&upsi;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>|</mo>
<mo>|</mo>
<mi>G</mi>
<mo>|</mo>
<mo>|</mo>
<mfrac>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>+</mo>
<mi>&delta;</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
重构故障设定定义重构的乘性故障:
其中,δ'为正常数参数,当时,将精确逼近f。
7.如权利要求6所述的基于自适应滑模观测器的电连接器间歇性失效检测方法,其特征在于:所述参数δ'的取值为,δ'=10-5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510153226.XA CN104914846B (zh) | 2015-04-01 | 2015-04-01 | 基于自适应滑模观测器的电连接器间歇性失效检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510153226.XA CN104914846B (zh) | 2015-04-01 | 2015-04-01 | 基于自适应滑模观测器的电连接器间歇性失效检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104914846A CN104914846A (zh) | 2015-09-16 |
CN104914846B true CN104914846B (zh) | 2017-12-22 |
Family
ID=54084002
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510153226.XA Expired - Fee Related CN104914846B (zh) | 2015-04-01 | 2015-04-01 | 基于自适应滑模观测器的电连接器间歇性失效检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104914846B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110058124A (zh) * | 2019-04-25 | 2019-07-26 | 中国石油大学(华东) | 线性离散时滞系统的间歇故障检测方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106843254B (zh) * | 2017-03-08 | 2019-08-09 | 北京航天自动控制研究所 | 一种实时主动重构容错控制方法 |
CN109142961A (zh) * | 2018-08-27 | 2019-01-04 | 中国人民解放军96901部队22分队 | 一种电连接器间歇失效机理研究方法 |
CN110764027B (zh) * | 2019-10-30 | 2021-11-02 | 江西理工大学 | 基于频谱特征变化的电连接器发生间歇故障诊断方法 |
CN111308390A (zh) * | 2020-02-21 | 2020-06-19 | 武汉科技大学 | 一种基于滑模观测器的模块化多电平换流器鲁棒故障检测方法 |
CN115113104A (zh) * | 2022-06-13 | 2022-09-27 | 云南电网有限责任公司电力科学研究院 | 电连接结构故障检测方法、存储介质及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102778889A (zh) * | 2012-07-12 | 2012-11-14 | 南京航空航天大学 | 航天器姿态控制系统的间歇性故障容错分析方法 |
CN103941725A (zh) * | 2014-04-24 | 2014-07-23 | 淮海工学院 | 一种非线性网络控制系统的故障诊断方法 |
JP2014153090A (ja) * | 2013-02-05 | 2014-08-25 | Yazaki Corp | 絶縁状態検出装置 |
CN104458298A (zh) * | 2014-12-09 | 2015-03-25 | 南京航空航天大学 | 基于多模型的高速列车悬架系统多执行器故障检测与隔离方法 |
-
2015
- 2015-04-01 CN CN201510153226.XA patent/CN104914846B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102778889A (zh) * | 2012-07-12 | 2012-11-14 | 南京航空航天大学 | 航天器姿态控制系统的间歇性故障容错分析方法 |
JP2014153090A (ja) * | 2013-02-05 | 2014-08-25 | Yazaki Corp | 絶縁状態検出装置 |
CN103941725A (zh) * | 2014-04-24 | 2014-07-23 | 淮海工学院 | 一种非线性网络控制系统的故障诊断方法 |
CN104458298A (zh) * | 2014-12-09 | 2015-03-25 | 南京航空航天大学 | 基于多模型的高速列车悬架系统多执行器故障检测与隔离方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110058124A (zh) * | 2019-04-25 | 2019-07-26 | 中国石油大学(华东) | 线性离散时滞系统的间歇故障检测方法 |
CN110058124B (zh) * | 2019-04-25 | 2021-07-13 | 中国石油大学(华东) | 线性离散时滞系统的间歇故障检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104914846A (zh) | 2015-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104914846B (zh) | 基于自适应滑模观测器的电连接器间歇性失效检测方法 | |
CN104579816B (zh) | 诊断车辆网络的方法和系统 | |
EP3379270A1 (en) | Method and system for non-destructive lightning strike pre-testing of carbon fiber composite components | |
CN106338385B (zh) | 一种基于奇异谱分解的旋转机械故障诊断方法 | |
US9563722B2 (en) | Sigma algebraic approximants as a diagnostic tool in power networks | |
CN112101554B (zh) | 异常检测方法及装置、设备、计算机可读存储介质 | |
CN107505538B (zh) | 基于线模电压行波幅值衰减特性的半波长输电线路非同步故障定位方法 | |
CN105206888A (zh) | 一种锂离子电池内部温度监测方法 | |
CN104850889A (zh) | 一种飞机旋转作动器驱动装置自适应故障检测、隔离与健康评估的方法 | |
CN107357275A (zh) | 非高斯工业过程故障检测方法及系统 | |
CN103018632B (zh) | 基于费歇信息的小电流接地系统单相接地故障选线方法 | |
CN108614539A (zh) | 机载设备故障诊断与预测模型验证方法 | |
CN103487759A (zh) | 一种混合动力汽车电池soc预测方法 | |
CN104833935A (zh) | 一种针对有功功率方向法定位谐波源局限性的检验方法 | |
CN110543717B (zh) | 一种gis设备温升特性仿真模型的构建方法及装置 | |
CN103827653A (zh) | 检测异常的发生的方法、设备及计算机程序 | |
CN105737881A (zh) | 一种电子设备加速可靠性增长试验方法 | |
CN109284933A (zh) | 一种基于数理统计的电子式互感器状态评估系统及方法 | |
CN111381130A (zh) | 一种计及行波波速的t接线路故障定位方法及系统 | |
CN103487696A (zh) | 一种母排温升一致性评价方法及装置 | |
CN108377209A (zh) | 基于scada的设备故障检测系统和检测方法 | |
CN117669388B (zh) | 故障样本生成方法、装置、计算机介质 | |
Guo et al. | Parameter estimation of a complex load model using phasor measurements | |
Na et al. | Detecting instant of multiple faults on the transmission line and its types using time–frequency analysis | |
CN105373112A (zh) | 一种基于多模型参数估计的舵机故障检测与诊断方法 |
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 |
Granted publication date: 20171222 |
|
CF01 | Termination of patent right due to non-payment of annual fee |