基于Kriging模型的电磁继电器快速计算方法技术

本发明专利技术目的在于提供一种基于近似模型的电磁继电器快速计算方法，属于继电器产品性能分析技术领域。首先，基于插值思想构建反映电磁继电器输出特性与过程变量的自定义函数作为Kriging模型的基函数；其次，通过拉丁超立方抽样方法建立Kriging模型的误差函数，从而构建基于Kriging方法的电磁系统近似模型；之后，基于变形能法建立反力特性计算近似模型；最后，基于数值方法求解继电器动态特性方程组，完成电磁继电器动态特性的快速计算。本发明专利技术提出的电磁继电器快速计算方法兼顾了计算精度及计算速度，可用于电磁继电器结构优化及稳健性设计等分析领域中。

Fast calculation method of electromagnetic relay based on Kriging model

The invention aims at providing a fast calculation method of electromagnetic relay based on an approximate model, belonging to the technical field of relay product performance analysis. First of all, the idea of constructing a custom function interpolation to reflect the output characteristics of electromagnetic relay and process variables as the basis function based on Kriging model; secondly, the error function by Latin hypercube sampling method to establish the Kriging model, so as to construct the electromagnetic system model based on Kriging method; after the establishment of anti deformation method and approximate calculation model based on force characteristics; finally, numerical methods for solving equations based on the dynamic characteristics of the relay, the fast calculation of dynamic characteristics of the electromagnetic relay. The electromagnetic relay rapid calculation method provided by the invention takes into account the calculation accuracy and the calculation speed, and can be used in the analysis field of electromagnetic relay structure optimization and robustness design, and the like.

【技术实现步骤摘要】
基于Kriging模型的电磁继电器快速计算方法
本专利技术涉及继电器产品性能分析
，特别涉及一种基于Kriging模型的电磁继电器快速计算方法。
技术介绍
电磁继电器作为一种能够对电路实现切换、控制、保护、检测和调节等功能的基本元器件，广泛应用于工业、航空航天、武器装备等领域。应用快速计算方法对其动态特性进行求解，对继电器的设计验证、性能评估以及产品优化具有重要意义。电磁继电器快速计算方法的主要难点在于电磁系统的快速计算。电磁系统的计算涉及电、磁、力多场耦合分析，传统分析方法主要有有限元方法和磁路法。有限元方法计算精度高但时效性差，适用于计算量小的场合；磁路法求解速度快，但结果精度低，适合在初步设计中应用。在结构优化及稳健性设计等分析领域中，往往需要求解速度快并且时效性高的方法，传统方法难以满足要求。因而，快速准确的求解电磁系统是对电磁继电器进行快速计算的前提。随着智能算法的发展，有学者提出采用数学手段建立机电元件计算的近似模型，从而实现电磁特性的快速计算。在常用的近似模型构建方法中，响应面法和移动最小二乘法对于非线性程度较大的函数难以准确求解；基于神经网络的径向基函数法模型参数求解困难，且需要重复计算大量数据，时效性较差；Kriging方法是基于误差的无偏估计确定模型参数，基函数可以更改，使用较灵活，而且求解速度快，适用于计算量大的场合。由于继电器电磁系统设计参数与输出特性的非线性程度较大，比较适合应用Kriging方法求解。
技术实现思路
本专利技术的目的在于提供一种能够提高灵活性及求解速度的基于Kriging模型的电磁继电器快速计算方法。为达上述目的，本专利技术公开一种基于Kriging模型的电磁继电器快速计算方法，其特征在于，包括以下步骤：S1：确定影响电磁系统暂态过程的关键节点；基于插值思想构建所述关键节点的输出特性与其它节点的输出特性之间的函数关系；建立反映电磁继电器输出特性与过程变量关系的自定义插值函数；S2：根据吸力曲线拐点位置选取插值节点(Ui,αj)；在电磁系统关键设计参数公差范围内均匀选取若干参数节点Δxk；应用有限元法计算各关键设计参数在(Ui,αj,Δxk)条件下电磁系统的输出电磁力矩；通过三次样条插值方法得到插值节点处各关键设计参数变化量与电磁力矩变化量之间的关系；应用拉丁超立方抽样在插值节点构成区域内选取采样点，基于量子粒子群算法寻优确定自定义插值函数中的影响系数；S3：以自定义插值函数作为Kriging模型的基函数；应用拉丁超立方抽样在参数公差范围内生成若干随机样本，通过有限元方法及所述基函数分别计算各随机样本的输出特性；将有限元方法及所述基函数这两种方式的计算结果的差作为插值节点建立Kriging模型的误差函数；综合所述基函数及误差函数建立电磁系统的近似模型；S4：根据接触系统反力部件图纸尺寸及尺寸链，确定影响触簧系统反力特性的因素，应用变形能法建立触簧系统反力计算模型；S5：通过电磁系统近似模型和触簧系统反力计算模型分别得到输入参数与电磁特性、机械特性的函数关系；基于数值方法求解电磁继电器动态特性数学模型以完成基于Kriging模型的电磁继电器快速计算，具体计算式如下：式中：U为额定线圈电压；t为持续时间；i为线圈电流；R为线圈电阻；为磁链；J为衔铁转动惯量；ω为衔铁转动角速度；α为衔铁转动角度；T为电磁吸力矩；Tf为反力矩；为线圈磁链初始值；i0为电流初始值；α0为衔铁角位移初始值。2、根据权利要求1所述的电磁继电器快速计算方法，其特征在于，步骤S1中，选取电磁继电器吸力曲线各拐点位置作为影响电磁系统暂态过程的关键节点；基于插值思想构建关键节点输出特性与其它点输出特性的函数关系；根据电磁继电器过程变量与输出特性之间的函数关系，确定包含影响系数的自定义插值函数。3、根据权利要求2所述的电磁继电器快速计算方法，其特征在于，在步骤S2中，将变量空间以关键节点(Ui,αj)划分为若干区域，应用拉丁超立方抽样在区域内选取N个采样点Xi＝[Ui,αi](i＝1,2,…,N)及对应目标函数值Gi，基于量子粒子群算法以误差函数达到最小值为优化目标，确定自定义插值函数中的影响系数。4、根据权利要求3所述的电磁继电器快速计算方法，其特征在于，在步骤S3中，所述随机样本的输出特性包括磁链和电磁吸力矩；将两种方式的计算结果的差作为插值节点得到误差函数Z，关联函数类型为高斯函数；综合所述基函数K及误差函数Z建立电磁系统的近似模型：y＝K+Z。5、根据权利要求4所述的电磁继电器快速计算方法，其特征在于，在步骤S4中，应用变形能法根据反力部件的受力情况计算其等效柔度及变形量，进而计算接触系统提供的反力矩，建立触簧系统反力计算模型。6、根据权利要求5所述的电磁继电器快速计算方法，其特征在于，步骤S5包括：S51：将时间离散化为若干个Δt，通过电磁系统近似模型和触簧系统反力计算模型分别得到输入参数与电磁特性、机械特性的函数关系：上式中为磁链，T为电磁吸力矩，Tf为机械反力矩，x为设计参数，i为线圈电流，α为衔铁转角；S52：根据初始状态α0＝0°、i0＝0A计算出通过如下迭代公式获得下一个Δt时刻的α1、i1、S53：通过反复迭代得到最终结果，完成基于Kriging模型的电磁继电器快速计算。与现有技术相比，本专利技术具有以下优点：本专利技术提出一种基于近似模型的电磁继电器快速计算方法，通过建立基于Kriging方法的电磁系统近似模型，并构建基于变形能法的反力特性计算近似模型，实现了电磁继电器尺寸参数改变的动态特性的高精度快速计算。本专利技术通过对电磁继电器尺寸参数改变的动态特性的高精度快速计算，得到尺寸参数与输出特性的对应关系，为电磁继电器产品的设计验证、性能评估及产品优化奠定了基础。附图说明图1是本专利技术所述基于Kriging模型的电磁继电器快速计算方法的流程图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。具体实施方式一下面结合图1说明本实施方式，本实施方式所述基于Kriging模型的电磁继电器快速计算方法，它包括以下步骤：S1：确定影响电磁系统暂态过程的关键节点；基于插值思想构建关键节点输出特性与其它点输出特性的函数关系；建立反映电磁继电器输出特性与过程变量关系的自定义插值函数；S2：根据吸力曲线拐点位置选取插值节点(Ui,αj)；在电磁系统关键设计参数公差范围内均匀选取若干参数节点Δxk；应用有限元法计算各关键设计参数在(Ui,αj,Δxk)条件下电磁系统的输出电磁力矩；通过三次样条插值方法得到插值节点处各关键设计参数变化量与电磁力矩变化量的关系；应用拉丁超立方抽样在插值节点构成区域内选取采样点，基于量子粒子群算法寻优确定所述自定义插值函数中的影响系数；S3：以所述自定义插值函数作为Kriging模型的基函数；应用拉丁超立方抽样在参数公差范围内生成若干随机样本，通过有限元方法及所述基函数分别计算各随机样本输出特性；将有限元方法及所述基函数这两种方式的计算结果的差本文档来自技高网...

【技术保护点】
一种基于Kriging模型的电磁继电器快速计算方法，其特征在于，包括以下步骤：S1：确定影响电磁系统暂态过程的关键节点；基于插值思想构建所述关键节点的输出特性与其它节点的输出特性之间的函数关系；建立反映电磁继电器输出特性与过程变量关系的自定义插值函数；S2：根据吸力曲线拐点位置选取插值节点(Ui,αj)；在电磁系统关键设计参数公差范围内均匀选取若干参数节点Δxk；应用有限元法计算各关键设计参数在(Ui,αj,Δxk)条件下电磁系统的输出电磁力矩；通过三次样条插值方法得到插值节点处各关键设计参数变化量与电磁力矩变化量之间的关系；应用拉丁超立方抽样在插值节点构成区域内选取采样点，基于量子粒子群算法寻优确定自定义插值函数中的影响系数；S3：以自定义插值函数作为Kriging模型的基函数；应用拉丁超立方抽样在参数公差范围内生成若干随机样本，通过有限元方法及所述基函数分别计算各随机样本的输出特性；将有限元方法及所述基函数这两种方式的计算结果的差作为插值节点建立Kriging模型的误差函数；综合所述基函数及误差函数建立电磁系统的近似模型；S4：根据接触系统反力部件图纸尺寸及尺寸链，确定影响触簧系统反力特性的因素，应用变形能法建立触簧系统反力计算模型；S5：通过电磁系统近似模型和触簧系统反力计算模型分别得到输入参数与电磁特性、机械特性的函数关系；基于数值方法求解电磁继电器动态特性数学模型以完成基于Kriging模型的电磁继电器快速计算，具体计算式如下：...

【技术特征摘要】
1.一种基于Kriging模型的电磁继电器快速计算方法，其特征在于，包括以下步骤：S1：确定影响电磁系统暂态过程的关键节点；基于插值思想构建所述关键节点的输出特性与其它节点的输出特性之间的函数关系；建立反映电磁继电器输出特性与过程变量关系的自定义插值函数；S2：根据吸力曲线拐点位置选取插值节点(Ui,αj)；在电磁系统关键设计参数公差范围内均匀选取若干参数节点Δxk；应用有限元法计算各关键设计参数在(Ui,αj,Δxk)条件下电磁系统的输出电磁力矩；通过三次样条插值方法得到插值节点处各关键设计参数变化量与电磁力矩变化量之间的关系；应用拉丁超立方抽样在插值节点构成区域内选取采样点，基于量子粒子群算法寻优确定自定义插值函数中的影响系数；S3：以自定义插值函数作为Kriging模型的基函数；应用拉丁超立方抽样在参数公差范围内生成若干随机样本，通过有限元方法及所述基函数分别计算各随机样本的输出特性；将有限元方法及所述基函数这两种方式的计算结果的差作为插值节点建立Kriging模型的误差函数；综合所述基函数及误差函数建立电磁系统的近似模型；S4：根据接触系统反力部件图纸尺寸及尺寸链，确定影响触簧系统反力特性的因素，应用变形能法建立触簧系统反力计算模型；S5：通过电磁系统近似模型和触簧系统反力计算模型分别得到输入参数与电磁特性、机械特性的函数关系；基于数值方法求解电磁继电器动态特性数学模型以完成基于Kriging模型的电磁继电器快速计算，具体计算式如下：式中：U为额定线圈电压；t为持续时间；i为线圈电流；R为线圈电阻；为磁链；J为衔铁转动惯量；ω为衔铁转动角速度；α为衔铁转动角度；T为电磁吸力矩；Tf为反力矩；为线圈磁链初始值；i0为电流初始值；α0为衔铁角位移初始值。2.根据权...

【专利技术属性】
技术研发人员：叶雪荣，林义刚，董宝旭，王瑛琪，付饶，邓杰，翟国富，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江,23