本发明涉及继电器产品设计技术领域，特别涉及一种基于变化贡献率的电磁继电器分步容差优化方法。
背景技术：
电磁继电器产品设计包括总体设计、参数设计和容差设计。总体设计完成系统功能实现，参数设计通过改变中心值提高产品抗干扰性，容差设计通过控制设计公差提高产品抗干扰性。中心值的改变不仅会影响抗干扰性而且会改变产品性能，分析复杂；容差设计不会改变中心值，应用范围更广泛。因此，采用容差设计方法提高电磁继电器批次产品的质量一致性意义重大。
容差优化是基于正交试验方法确定各因素对输出结果的贡献率，然后依据贡献率对容差进行反分配得到满足要求的容差方案。传统的容差优化方法依据各因素一、二次贡献率及输出特性改变量确定各因素容差的变化范围，这种直接确定新的容差方案有三个缺点：分配方案并不唯一，需要手动调节各因素容差范围以满足输出特性要求；分配方式只根据原有公差的贡献率，没有考虑缩短容差后原先不显著的因素可能变得显著的情况，分配出的方案并非最优；分配后的方案并没有验证，不能保证新的容差方案一定满足要求。
技术实现要素：
本发明目的是为了解决目前传统的容差优化方法需要手动调节各因素容差范围等问题，提供了一种基于变化贡献率的电磁继电器分步容差优化方法。
为达上述目的，本发明提出一种基于变化贡献率的电磁继电器分步容差优化方法，包括以下步骤：
s1：根据电磁继电器产品质量一致性要求确定容差优化的标准；根据关键因素与输出特性之间的关系以及实际加工能力，确定所述关键因素的公差范围；
s2：根据步骤s1中确定的所述公差范围确定所述关键因素的水平数；根据所 述关键因素的数量和所述水平数设计正交试验，利用快速计算模型计算每组所述正交试验的输出结果；根据所述输出结果求取各所述关键因素的贡献率；
s3：将所述贡献率按照一定容差步长减小容差，获得所述关键因素的新公差范围；按照所述新公差范围进行正交试验，再次计算各所述关键因素的贡献率；
s4：基于蒙特卡罗方法，在步骤s3确定的所述新公差范围内产生虚拟样本，应用快速计算模型计算得到所述虚拟样本的输出特性分布，求出对应的新公差范围的输出特性波动范围；
s5：检验步骤s4得到的输出特性波动范围是否满足容差优化标准，如果不满足，则重复步骤s3-s4的工作，直到所述输出特性波动范围满足容差优化标准，或者各所述关键因素的新公差范围已达到实际生产过程中设备的加工精度下限，此时各关键因素的所述新公差范围即为分步容差优化的最终结果。
根据本发明提出的优化方法，其中步骤s3中，在容差优化初始阶段设置的所述容差步长比在容差优化末尾阶段设置的容差步长大。
根据本发明提出的优化方法，其中，所述步骤s3还包括，在容差步长计算中引入符合所述容差步长的设置的加速系数，然后结合贡献率得到所述容差步长的计算公式：
ii+1＝αρimi
式中ii+1为第i+1步的步长；ρi为第i步计算的贡献率；mi为第i步的容差极限值；α为加速系数，按照下式计算
式中ηi表示第i步计算的输出特性波动范围与容差优化标准的比值；
第i+1步的容差极限值mi+1为：
式中ε为实际加工能力下限值；[-mi+1,mi+1]即为新的容差范围。
与现有技术相比，本发明具有以下有益效果：
本发明方法基于变化贡献率的电磁继电器分步容差优化方法，使得传统容差分配方法中对于缩短容差后原先不显著的因素可能变得显著的情况能够得以 考虑，保证了容差分配方案最优；另外，本发明容差优化速度高，能够自动调节各因素容差范围以满足输出特性质量一致性要求。
本发明用于在电磁继电器产品的设计阶段，根据产品的质量一致性要求以及实际加工能力，利用正交试验设计计算容差分配后的变化了的贡献率，引入容差步长，通过多次迭代，得到能够满足产品质量一致性要求与加工能力要求的最优容差分配方案。
附图说明
图1为本发明基于变化贡献率的电磁继电器分步容差优化方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
具体实施方式一：
下面结合图1说明本实施方式，本实施方式所述基于变化贡献率的电磁继电器分步容差优化方法，它包括以下步骤：
s1：根据电磁继电器产品质量一致性要求确定容差优化的标准；根据关键因素与输出特性之间的关系以及实际加工能力，确定关键因素的公差范围；
s2：根据新确定的公差范围确定关键因素的水平数；根据关键因素数量和水平数设计正交试验，利用快速计算模型计算每组试验方案的输出结果；根据正交试验结果求取各因素的贡献率；
s3：根据各因素贡献率按照一定容差步长减小容差，获得关键因素的新公差范围；按照新公差范围进行正交试验，计算各因素的贡献率；
s4：基于蒙特卡罗方法，在步骤s3确定的关键因素新公差范围内产生虚拟样本，应用快速计算模型计算得到虚拟样本的输出特性分布，求出对应关键因素新公差范围的输出特性的波动范围；
s5：检验步骤s4得到的输出特性波动范围是否满足容差优化标准，如果不满足，则重复步骤s3-步骤s5的工作，直到输出特性波动范围满足容差优化标准或者各关键因素的公差范围已达到实际生产过程中设备的加工精度下限，此 时各关键因素的公差范围即分步容差优化的最终结果。
具体实施方式二：
下面结合图1说明本实施方式，本实施方式对实施方式一作进一步说明，步骤s1中根据电磁继电器产品具体输出特性参数(如吸合电压、释放电压等)的质量一致性要求确定输出特性参数容差优化的标准，即输出特性参数上下波动的范围；然后根据关键因素(即设计参数)与输出特性之间的关系，确定关键因素的公差范围，即如果关键因素变化对输出特性影响小，则可分配较大的公差范围，如果关键因素变化对输出特性影响较大，则可分配较小的公差范围，分配的公差都应在加工能力许可范围之内。
在步骤s2中，确定关键因素的水平数为3水平，中心值、中心值+公差上限、中心值-公差下限各占1个水平；根据关键因素数量和水平数选取正交表，安排试验方案，通过电磁继电器动态特性快速计算模型对每组试验方案进行计算，求得输出特性值yi，然后代入下式计算贡献率：
式中n为试验方案总数；ρil、ρiq分别为关键因素xi对应输出特性y的一次项贡献率和二次项贡献率；ρe为误差项贡献率；sil、siq分别为关键因素xi对应输出特性y的一次项波动平方和和二次项波动平方和；se为误差项平方和；ve为误差方差；sil、siq、se、ve的计算公式如下：
式中t1i、t2i、t3i为试验结果中，分别对应关键因素xi三个水平的部分和；dfe为误差项自由度。如果sil和siq中有比ve小的项存在，需要将其归并到se中，并且每归并一项，dfe需要加1；st为输出特性的总偏差平方和，按照下式计算：
式中
步骤s3所述的容差步长在容差优化初始阶段可以设置较大以加快优化速度的目的，在容差优化的末尾阶段容差步长应该设置较小以达到结果最优的目的。因此，在容差步长计算中引入能够符合上述要求的加速系数，然后结合贡献率 得到容差步长的计算公式：
ii+1＝αρimi
式中ii+1为第i+1步的步长；ρi为第i步计算的贡献率；mi为第i步的容差极限值；α为加速系数，按照下式计算：
式中ηi表示第i步计算的输出特性波动范围与容差优化标准的比值。
第i+1步的容差(或公差)极限值为：
式中mi+1为第i+1步的容差极限值；ε为实际加工能力下限值。[-mi+1,mi+1]即为新的容差(或公差)范围。
步骤s4中，在各关键因素公差范围内基于蒙特卡罗思想产生符合正态分布的各关键因素虚拟样本100个，然后进行随机组合得到100个电磁继电器虚拟样本，通过继电器动态特性快速计算模型计算得到虚拟样本的输出特性值yk，代入标准偏差计算公式即可求出能够反映输出特性波动范围的标准偏差。
本领域普通技术人员可以理解：附图只是一个实施例的示意图，附图中的模块或流程并不一定是实施本发明所必须的。
本领域普通技术人员可以理解：实施例中的装置中的模块可以按照实施例描述分布于实施例的装置中，也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。