[发明专利]电磁继电器静、动态特性快速计算方法

权利要求书

1.一种电磁继电器静、动态特性快速计算方法，其特征在于：它包括修正系数的定义及求取、静态吸力特性的快速求解和动态特性的快速求解，其特征在于：

所述的修正系数的定义是指同一电磁系统模型的有限元法计算结果与磁路法计算结果的比值；

所述的修正系数的求取步骤如下：

步骤一：对尺寸参数的任意取值组合(x′₁，x′₂，…，x′_n)，根据公式：

η(x1,x2,...,xn)=fFEM(x1,x2,...,xn)fMEC(x1,x2,...,xn)]]>

式中：

η(x₁，x₂，…，x_n)——尺寸参数取值为x₁，x₂，…，x_n时的修正系数；

f_FEM(x₁，x₂，…，x_n)——尺寸参数取值为x₁，x₂，…，x_n时电磁系统的有限元法计算结果；

f_MEC(x₁，x₂，…，x_n)——尺寸参数取值为x₁，x₂，…，x_n时电磁系统的磁路法计算结果；

分别采用磁路法和有限元法求电磁系统在尺寸参数取值点及其邻域内点处的修正系数其中或或或

xn*≤xn′<xnr;]]>

步骤二：由公式：∂η∂x1≈η(x1r,x2*,...xn*)-η(x1*,x2*,...xn*)x1r-x1*∂η∂x2≈η(x1*,x2r,...xn*)-η(x1*,x2*,...xn*)x2r-x2*...∂η∂xn≈η(x1*,x2*,...xnr)-η(x1*,x2*,...xn*)xnr-xn*]]>

式中：

——分别为修正系数关于尺寸参数x₁，x₂和x_n的偏导数；

——分别为尺寸参数取值组合为时的修正系数；

求得各尺寸参数在点处的偏导数下角标i＝1，2，…，n；

步骤三：由公式：dη′|(x1*,x2*,...,xn*)=∂η∂x1(x1′-x1*)+∂η∂x2(x2′-x2*)+...+∂η∂xn(xn′-xn*)]]>

式中：

——尺寸参数取值点(x′₁，x′₂…，x′_n)关于取值点的修正系数增量；

——分别为修正系数关于尺寸参数x₁，x₂和x_n的偏导数；

求得修正系数η(x′₁，x′₂，…，x′_n)相对于修正系数的增量

步骤四：由公式：η(x1′,x2′,...,xn′)=η(x1*,x2*,...,xn*)+dη′|(x1*,x2*,...,,xn*)]]>

式中：

η(x′₁，x′₂…，x′_n)——尺寸参数取值为x′₁，x′₂…，x′_n时的修正系数；

——尺寸参数取值为时的修正系数；

——尺寸参数取值点(x′₁，x′₂，…，x′_n)关于取值点的修正系数增量；

求得尺寸参数的任意取值组合(x′₁，x′₂，…，x′_n)对应的磁路模型结果和有限元模型计算结果对比得出的修正系数；

所述的静态吸力特性的快速求解步骤如下：

步骤一：列出电磁系统尺寸参数x₁，x₂，...，x_n及其在各自取值范围[a_i，b_i]内的m_r+2个取值点x_i0，x_i1，…，其中该步骤中下角标i＝1，2，...，n，下角标r＝1，2，...，n；

步骤二：对尺寸参数的取值点x_ij及，其中下角标i＝1，2，...，n，下角标i＝0，1，...，m_i，分别采用磁路和有限元法对电磁系统模型进行静态特性求解，以求得其对应的修正系数η_Fijk及其中，下角标F表示是吸力的修正系数；下角标i＝1，2，...，n；下角标j＝0，1，...，m_i；下角标k＝1，2，...，K；K为求解静态特性时衔铁位移计算点的个数；

步骤三：对尺寸参数的任意取值组合(x′₁，x′₂，…，x′_n)，a_i≤x′_i≤b_i，先将各取值点x′_i，与其对应的参数x_i的取值点进行比较，找到距离x′_i最近的除之外的x_i的取值点x_it  ，根据公式求得对应的偏导数近似值然后用式和式即可求得取值组合(x′₁，x′₂，…，x′_n)下的η′_Fk(x′₁，x′₂，…，x′_n)，再根据磁路模型求得该组合下的静态特性各计算点的吸力并利用公式即可完成静态特性的快速准确求解，该步骤中，下角标i＝1，2，...，n，下角标t取0，1，...，j_i-1，j_i，...，m_i其中一个值，下角标k＝1，2，...，K；

所述的动态特性的快速求解步骤如下：

步骤一：列出电磁系统尺寸参数x₁，x₂，...，x_n及其在各自取值范围[a_i，b_i]内的m_r+2个取值点x_i0，x_i1，...，其中该步骤中，下角标i＝1，2，...，n，下角标r＝1，2，...，n；

步骤二：对尺寸参数的取值点x_ij及其中i＝1，2，...，n，下角标j＝0，1，2，...，m_i，分别采用磁路和有限元法对不同线圈电流不同衔铁位移下的电磁系统模型进行衔铁吸力和线圈磁链的求解，其中，电流值均匀分布在0A到线圈电流稳态值之间，衔铁位移计算点均匀分布在吸合和释放位置之间，以求得其对应的吸力和磁链修正系数η_Fijpq、η_ψijpq和其中，下角标F表示是吸力的修正系数；下角标ψ表示是磁链的修正系数；下角标i＝1，2，...，n；下角标j＝0，1，2，...，m_i；下角标p＝1，2，...，P；下角标q＝1，2，...，Q；P为电流在0A到线圈电流稳态值之间的取值个数，Q为衔铁位移在吸合和释放位置之间的取值个数；

步骤三：对尺寸参数的任意取值组合(x′₁，x′₂，…，x′_n)，a_i≤x′_i≤b_i，求其对应电磁系统模型在不同线圈电流不同衔铁位移下的衔铁吸力和线圈磁链的修正系数η_Fpq(x′₁，x′₂…，x′_n)和η_ψpq(x′₁，x′₂，…，x′_n)该步骤中，下角标i＝1，2，...，n，下角标p＝1，2，...，P；下角标q＝1，2，...，Q；

步骤四：根据步骤三求得的尺寸参数任意取值组合(x′₁，x′₂，…，x′_n)对应的电磁系统在不同线圈电流不同衔铁位移下的衔铁吸力和线圈磁链的修正系数，建立该电磁系统磁路模型求解不同线圈电流不同衔铁位移下的吸力和磁链，并利用公式：求得该取值组合下电磁系统在不同线圈电流不同衔铁位移下修正后的吸力和磁链；

步骤五：根据已求得的尺寸参数任意取值组合(x′₁，x′₂，…，x′_n)对应的电磁系统在不同线圈电流不同衔铁位移下的修正后的衔铁吸力和线圈磁链，采用4阶龙格库塔方法求解电磁继电器吸合过程状态

方程组：dψdt=U-Ri(ψ,s)dvdt=Fs(ψ,s)-Ff(s)mdsdt=vψ|t=0=0,v|t=0=0,s|t=0=s0]]>

式中：

ψ——线圈磁链；t——时间；U——线圈电压；R——线圈电阻；

i——线圈电流；s——衔铁位移；v——衔铁速度；F_x——衔铁

受到的电磁吸力；F_f——衔铁受到的机械反力；m——衔铁的质

量；s_o——衔铁初始位移；

完成电磁继电器动态特性的快速求解。