CN108009322B - 高速机械开关多场耦合仿真计算方法及其计算装置 - Google Patents

高速机械开关多场耦合仿真计算方法及其计算装置 Download PDF

Info

Publication number
CN108009322B
CN108009322B CN201711131430.7A CN201711131430A CN108009322B CN 108009322 B CN108009322 B CN 108009322B CN 201711131430 A CN201711131430 A CN 201711131430A CN 108009322 B CN108009322 B CN 108009322B
Authority
CN
China
Prior art keywords
mechanical switch
speed mechanical
coil
speed
repulsive
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.)
Active
Application number
CN201711131430.7A
Other languages
English (en)
Other versions
CN108009322A (zh
Inventor
李志兵
田宇
田阳
董恩源
徐晓东
朱玉
刘北阳
王浩
颜湘莲
刘焱
曹德新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
Dalian University of Technology
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology, State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI filed Critical Dalian University of Technology
Priority to CN201711131430.7A priority Critical patent/CN108009322B/zh
Publication of CN108009322A publication Critical patent/CN108009322A/zh
Application granted granted Critical
Publication of CN108009322B publication Critical patent/CN108009322B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Abstract

本发明提供了一种高速机械开关多场耦合仿真计算方法及其计算装置,其中,该方法包括如下步骤:对高速机械开关进行电磁场仿真分析,以确定高速机械开关的瞬态线圈电流和斥力金属盘的电磁斥力;根据瞬态线圈电流对高速机械开关进行热电耦合场仿真分析,以获得线圈的温升曲线;根据电磁斥力对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得高速机械开关的位移曲线。本发明中通过对高速机械开关进行电磁场仿真分析、热电耦合场仿真分析和分闸反弹运动场耦合仿真分析,可以全面反映高速机械开关运动、触头碰撞和线圈温升情况,从而可以全面准确反映高速机械开关的机械性能,进而为混合式直流断路器快速机械开关优化设计提供了理论指导。

Description

高速机械开关多场耦合仿真计算方法及其计算装置
技术领域
本发明涉及断路器仿真计算技术领域,具体而言,涉及一种高速机械开关多场耦合仿真计算方法及其计算装置。
背景技术
高速机械开关是混合式直流断路器关键部件,其性能主要由其操动机构决定,目前主要采用电磁斥力机构,普遍存在出力时间短、操动距离小的缺点,因此多被应用于中低压电压等级。对于高电压等级高速机械开关,国内外普遍采用多个中低压等级高速机械开关单元串联同步操作结构形式。因此,高速机械开关的机械特性及稳定性显得尤为重要。
高速机械开关动作过程涉及电磁场、温度场及运动场,而目前对快速开关的研究多基于电磁场单场进行仿真设计,并不能全面准确反映机械性能。
发明内容
鉴于此,本发明提出了一种高速机械开关多场耦合仿真计算方法及其计算装置,旨在解决目前对高速机械开关的研究多基于电磁场单场进行仿真设计,导致的不能全面准确反映机械性能的问题。
一个方面,本发明提出了一种高速机械开关多场耦合仿真计算方法,该方法包括如下步骤:对高速机械开关进行电磁场仿真分析,以确定高速机械开关的瞬态线圈电流和斥力金属盘的电磁斥力;根据瞬态线圈电流对高速机械开关进行热电耦合场仿真分析,以获得线圈的温升曲线;根据电磁斥力对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得高速机械开关的位移曲线。
进一步地,上述高速机械开关多场耦合仿真计算方法中,上述确定电磁斥力的方法为:根据方程组
Figure BDA0001469637890000021
Figure BDA0001469637890000022
Figure BDA0001469637890000023
确定电磁斥力Fc,其中公式(3)的初始条件为v|t=0=v0=0,x|t=0x0=0;上式中,ψq、ψg均为电磁系统全磁链,t为高速机械开关动作时间,Lq为斥力线圈自感,iq为脉冲电流,M为斥力线圈与斥力金属盘的互感,ig为感应电流,Uc为电容电压,Rq为斥力线圈电阻,Lg为斥力金属盘的自感,Rg为斥力金属盘的电阻,Ff为斥力金属盘受到的运动反力,x为斥力金属盘的位移,v1f为高速机械开关运动部件碰撞后初速度,v2f为高速机械开关运动部件碰撞后终速度,v1i为高速机械开关运动部件碰撞前初速度,v2i为高速机械开关运动部件碰撞前终速度,e为碰撞恢复系数。
进一步地,上述高速机械开关多场耦合仿真计算方法中,上述获得线圈的温升曲线的方法为:根据公式
Figure BDA0001469637890000024
其中初始条件为
Figure BDA0001469637890000025
确定线圈随时间的温升;上式中,τ为斥力线圈温升,T为热时间常数,t为高速机械开关动作时间,e为碰撞恢复系数,c为线圈的材质的比热容,m为线圈的质量,KT为线圈的散热系数,A为线圈的散热面积,P为线圈的发热功率;根据线圈随时间的温升确定温升曲线。
进一步地,上述高速机械开关多场耦合仿真计算方法中,上述获得高速机械开关的位移曲线的方法为:将电磁斥力导入预先建立的计算模型中;对高速机械开关的部件连接部分添加相关类型运动副,并设置部件的参数;对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得位移曲线。
进一步地,上述高速机械开关多场耦合仿真计算方法中,电磁斥力通过AKISPL样条差值函数和DM距离函数导入计算模型中。
本发明通过对高速机械开关进行电磁场仿真分析、热电耦合场仿真分析和分闸反弹运动场耦合仿真分析,可以全面反映高速机械开关运动、触头碰撞和线圈温升情况,从而可以全面准确反映高速机械开关的机械性能,以得到对高速机械开关最优化的设计,进而为混合式直流断路器快速机械开关优化设计提供了理论指导。
另一方面,本发明还提出了一种高速机械开关多场耦合仿真计算装置,该装置包括:电磁场分析模块,用于对高速机械开关进行电磁场仿真分析,以确定高速机械开关的瞬态线圈电流和斥力金属盘的电磁斥力;瞬态热分析模块,用于根据瞬态线圈电流对高速机械开关进行热电耦合场仿真分析,以获得线圈的温升曲线;动力学分析模块,用于根据电磁斥力对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得高速机械开关的位移曲线。
进一步地,上述高速机械开关多场耦合仿真计算装置中,电磁场分析模块包括:电磁斥力确定子模块,用于根据方程组
Figure BDA0001469637890000031
Figure BDA0001469637890000032
Figure BDA0001469637890000033
确定电磁斥力Fc,其中公式(3)的初始条件为v|t=0=v0=0,x|t=0x0=0;上式中,ψq、ψg均为电磁系统全磁链,t为高速机械开关动作时间,Lq为斥力线圈自感,iq为脉冲电流,M为斥力线圈与斥力金属盘的互感,ig为感应电流,Uc为电容电压,Rq为斥力线圈电阻,Lg为斥力金属盘的自感,Rg为斥力金属盘的电阻,Ff为斥力金属盘受到的运动反力,x为斥力金属盘的位移,v1f为高速机械开关运动部件碰撞后初速度,v2f为高速机械开关运动部件碰撞后终速度,v1i高速机械开关运动部件碰撞前初速度,v2i为高速机械开关运动部件碰撞前终速度,e为碰撞恢复系数。
进一步地,上述高速机械开关多场耦合仿真计算装置中,瞬态热分析模块包括:线圈温升确定子模块,用于根据公式
Figure BDA0001469637890000041
其中初始条件为
Figure BDA0001469637890000042
确定线圈随时间的温升;上式中,τ为斥力线圈温升,T为热时间常数,t为高速机械开关动作时间,e为碰撞恢复系数,c为线圈的材质的比热容,m为线圈的质量,KT为线圈的散热系数,A为线圈的散热面积,P为线圈的发热功率;温升曲线确定子模块,用于根据线圈随时间的温升确定温升曲线。
进一步地,上述高速机械开关多场耦合仿真计算装置中,运动学分析模块包括:导入子模块,用于将电磁斥力导入预先建立的计算模型中;设置子模块,用于对高速机械开关的部件连接部分添加相关类型运动副,并设置部件的参数;仿真分析子模块,用于对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得位移曲线。
进一步地,上述高速机械开关多场耦合仿真计算装置中,电磁斥力通过AKISPL样条差值函数和DM距离函数导入计算模型中。
本发明通过对高速机械开关进行电磁场仿真分析、热电耦合场仿真分析和分闸反弹运动场耦合仿真分析,可以全面反映高速机械开关运动、触头碰撞和线圈温升情况,从而可以全面准确反映高速机械开关的机械性能,以得到对高速机械开关最优化的设计,进而为混合式直流断路器快速机械开关优化设计提供了理论指导。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1为本发明实施例提供的高速机械开关多场耦合仿真计算方法的流程图;
图2为本发明实施例提供的线圈匝数对高速机械开关出力影响的示意图;
图3为本发明实施例提供的线圈匝数对线圈电流影响的示意图;
图4为本发明实施例提供的线圈匝数对高速机械开关运动速度影响的示意图;
图5为本发明实施例提供的线圈匝数对位移影响的示意图;
图6为本发明实施例提供的斥力盘材质对高速机械开关出力影响的示意图;
图7为本发明实施例提供的斥力盘材质对线圈电流影响的示意图;
图8为本发明实施例提供的斥力盘材质对高速机械开关运动速度影响的示意图;
图9为本发明实施例提供的斥力盘材质对位移影响的示意图;
图10为本发明实施例提供的线圈与斥力盘初始间距对高速机械开关出力影响的示意图;
图11为本发明实施例提供的线圈与斥力盘初始间距对线圈电流影响的示意图;
图12为本发明实施例提供的线圈与斥力盘初始间距对高速机械开关运动速度影响的示意图;
图13为本发明实施例提供的线圈与斥力盘初始间距对位移影响的示意图;
图14为本发明实施例提供的磁感线回路对高速机械开关出力影响的示意图;
图15为本发明实施例提供的磁感线回路对线圈电流影响的示意图;
图16为本发明实施例提供的磁感线回路对高速机械开关运动速度影响的示意图;
图17为本发明实施例提供的磁感线回路对位移影响的示意图;
图18为本发明实施例提供的不同匝数的线圈的温升曲线示意图;
图19为本发明实施例提供的不同匝数的线圈对位移影响的示意图;
图20为本发明实施例提供的高速机械开关等效电路示意图;
图21为本发明实施例提供的获得线圈的升温曲线的方法的流程图;
图22为本发明实施例提供的获得高速机械开关的位移曲线的方法的流程图;
图23为本发明实施例提供的高速机械开关多场耦合仿真计算装置的结构框图;
图24为本发明实施例提供的高速机械开关多场耦合仿真计算装置中,电磁场分析模块的结构框图;
图25为本发明实施例提供的高速机械开关多场耦合仿真计算装置中,瞬态热分析模块的结构框图;
图26为本发明实施例提供的高速机械开关多场耦合仿真计算装置中,运动学分析模块的结构框图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
方法实施例:
参见图1,图1为本实施例提供的高速机械开关多场耦合仿真计算方法的流程图。如图所示,该方法包括如下步骤:
步骤S110,对高速机械开关进行电磁场仿真分析,以确定高速机械开关的瞬态线圈电流和斥力金属盘的电磁斥力。
具体而言,以成熟通用的ANSYS有限元分析软件中低频电磁场分析模块(Maxwell)、瞬态热分析模块(Transient thermal)以及多体动力学分析软件(Adams)为基础,利用C++语言编写高速机械开关多场耦合仿真分析程序,此C++程序包括3个子程序:C++调用Maxwell计算电磁斥力、脉冲电流、速度及位移、线圈损耗程序;关键参数导入Transient thermal温升计算程序;关键参数导Adams碰撞反弹计算程序。在给定高速机械开关的电容电压、位移初值的情况下,首先使用Maxwell电磁场分析软件对高速机械开关机械特性有明显影响的关键因素进行仿真优化设计,以确定高速机械开关的瞬态线圈电流、斥力金属盘的电磁斥力、斥力金属盘的速度、位移和线圈损耗。关键因素包括分闸线圈匝数、斥力盘材质、线圈与斥力盘初始间距、磁感线回路变化等影响因素,即分别取不同线圈匝数、不同斥力盘材质、不同工作间距及不同的空间磁路进行电磁场的仿真分析,研究上述因素对高速机械开关出力、运动速度、线圈电流及位移的影响,仿真结果可以参见图2至图17所示。
步骤S120,根据瞬态线圈电流对高速机械开关进行热电耦合场仿真分析,以获得线圈的温升曲线。
具体而言,将电磁场仿真分析得到的通入线圈的瞬态电流波形导入AnsysWorkbench软件中进行热电耦合场的仿真分析,研究线圈匝数及高速机械开关工作频率对线圈温升的影响,仿真得到高速机械开关完成单次分合操作时,不同匝数线圈的温升曲线,如图18所示。
步骤S130,根据电磁斥力对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得高速机械开关的位移曲线。
具体而言,利用外编程序调用maxwell电磁场仿真计算的不同斥力线圈匝数下电磁斥力结果,开展不同线圈匝数下分闸反弹仿真分析,进而得到高速机械开关的位移曲线,如图19所示。通过对高速机械开关电磁场、温度场和运动场多场耦合仿真分析,获得影响高速机械开关动态特性影响因素及其变化规律,进而可以得出其中主要影响因素为线圈匝数,结合多场耦合分析结果可知,斥力线圈为20匝时,高速机械开关分闸时间、允许温升及分闸反弹等关键技术指标满足技术要求。
可以看出,本实施例中,通过对高速机械开关进行电磁场仿真分析、热电耦合场仿真分析和分闸反弹运动场耦合仿真分析,可以全面反映高速机械开关运动、触头碰撞和线圈温升情况,从而可以全面准确反映高速机械开关机械性能,以得到对高速机械开关最优化的设计,进而为混合式直流断路器快速机械开关优化设计提供了理论指导。
高速机械开关动作过程包含着电、磁、热和机械参量等多种场域变量的变化及相互影响,其中电磁参量耦合是激励源,热及机械参量变化是受控源,因此需从激励源电路维度建立斥力机构等效电路,进而建立全面的多物理场耦合微分方程组。其电磁斥力机构等效电路可以参见图20。
根据方程组
Figure BDA0001469637890000081
Figure BDA0001469637890000082
Figure BDA0001469637890000083
确定电磁斥力Fc,其中公式(3)的初始条件为v|t=0=v0=0,x|t=0x0=0;
上式中,ψq、ψg均为电磁系统全磁链,t为高速机械开关的动作时间,Lq为斥力线圈的自感,iq为脉冲电流,M为斥力线圈与斥力金属盘的互感,ig为感应电流,Uc为电容电压,Rq为线圈的电阻,Lg为斥力金属盘的自感,Rg为斥力金属盘的电阻,Ff为斥力金属盘受到的运动反力,x为斥力金属盘的位移,v1f为高速机械开关运动部件碰撞后初速度,v2f为高速机械开关运动部件碰撞后终速度,v1i为高速机械开关运动部件碰撞前初速度,v2i为高速机械开关运动部件碰撞前终速度,e为碰撞恢复系数。
参见图21,图21为本实施例提供的获得线圈的升温曲线的方法的流程图。如图所示,该方法可以包括:
步骤S2110,根据公式
Figure BDA0001469637890000091
其中初始条件为
Figure BDA0001469637890000092
确定线圈随时间的温升;上式中,τ为斥力线圈的温升,T为热时间常数,t为高速机械开关对的动作时间,e为碰撞恢复系数,c为线圈的材质的比热容,m为线圈的质量,KT为线圈的散热系数,A为线圈的散热面积,P为线圈的发热功率;
具体而言,一般情况下线圈都是采用铜线,所以c可以为铜的比热容。
步骤S2120,根据线圈随时间的温升确定温升曲线。
具体而言,根据线圈的不同时间的温度,可以绘制线圈随时间的温升曲线。
本实施例中,通过公式可以计算得到线圈随时间的温升,进而可以确定温升曲线,通过观察温升曲线,可以更清晰的了解到线圈匝数对温升的影响。
参见图22,图22为本实施例提供的获得高速机械开关的位移曲线的方法的流程图。如图所示,该方法可以包括:
步骤S2210,将电磁斥力导入预先建立的计算模型中。
具体而言,通过AKISPL样条差值函数和DM距离函数,将电磁斥力导入预先建立的Adams计算模型中。
步骤S2220,对高速机械开关的部件连接部分添加相关类型运动副,并设置部件的参数。
具体而言,高速机械开关可以为双稳态弹簧机构,对其部件连接部分添加相关类型运动副,并设置部件的参数,使其能够按照设计要求进行分闸碰撞动作过程仿真分析。
步骤S2230,对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得位移曲线。
本实施例中,通过将电磁斥力导入预先建立的计算模型中,并对高速机械开关进行分闸反弹运动场耦合仿真分析,可以获得位移曲线,通过位移曲线可以清晰地了解到线圈匝数对位移的影响。
综上,本实施例通过对高速机械开关进行电磁场仿真分析、热电耦合场仿真分析和分闸反弹运动场耦合仿真分析,可以全面反映高速机械开关运动、触头碰撞和线圈温升情况,从而可以全面准确反映高速机械开关的机械性能,以得到对高速机械开关最优化的设计,进而为混合式直流断路器快速机械开关优化设计提供了理论指导。
装置实施例:
参见图23,图23为本实施例提供的高速机械开关多场耦合仿真计算装置的结构框图。如图所示,该装置包括:电磁场分析模块100、瞬态热分析模块200和动力学分析模块300。其中,电磁场分析模块用于对高速机械开关进行电磁场仿真分析,以确定高速机械开关的瞬态线圈电流和斥力金属盘的电磁斥力;瞬态热分析模块用于根据瞬态线圈电流对高速机械开关进行热电耦合场仿真分析,以获得线圈的温升曲线;动力学分析模块用于根据电磁斥力对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得高速机械开关的位移曲线。其中,该装置的具体实施过程参见上述方法实施例中的说明即可,本实施例在此不再赘述。
可以看出,本实施例中,通过对高速机械开关进行电磁场仿真分析、热电耦合场仿真分析和分闸反弹运动场耦合仿真分析,可以全面反映高速机械开关运动、触头碰撞和线圈温升情况,从而可以全面准确反映高速机械开关的机械性能,以得到对高速机械开关最优化的设计,进而为混合式直流断路器快速机械开关优化设计提供了理论指导。
参见图24,图24为本实施例提供的高速机械开关多场耦合仿真计算装置中,电磁场分析模块的结构框图。如图所示,电磁场分析模块100包括:电磁斥力确定子模块110,用于根据方程组
Figure BDA0001469637890000111
Figure BDA0001469637890000112
Figure BDA0001469637890000113
确定电磁斥力Fc,其中公式(3)的初始条件为vt=0=v0=0,xt=0x0=0;
上式中,ψq、ψg均为电磁系统全磁链,t为斥力机构动作时间,Lq为斥力线圈自感,iq为脉冲电流,M为斥力线圈与斥力金属盘的互感,ig为感应电流,Uc为电容电压,Rq为线圈电阻,Lg为斥力金属盘的自感,Rg为斥力金属盘的电阻,Ff为斥力金属盘受到的运动反力,x为斥力金属盘的位移,v1f为高速机械开关运动部件碰撞后初速度,v2f为高速机械开关运动部件碰撞后终速度,v1i高速机械开关运动部件碰撞前初速度,v2i为高速机械开关运动部件碰撞前终速度,e为碰撞恢复系数。其中,该装置的具体实施过程参见上述方法实施例中的说明即可,本实施例在此不再赘述。
参见图25,图25为本实施例提供的高速机械开关多场耦合仿真计算装置中,瞬态热分析模块的结构框图。如图所示,瞬态热分析模块200包括:线圈温升确定子模块210和温升曲线确定子模块220。其中,线圈温升确定子模块210用于根据公式
Figure BDA0001469637890000121
其中初始条件为
Figure BDA0001469637890000122
确定线圈随时间的温升;上式中,τ为斥力线圈温升,T为热时间常数,t为斥力机构动作时间,e为碰撞恢复系数,c为线圈的材质的比热容,m为线圈的质量,KT为线圈的散热系数,A为线圈的散热面积,P为线圈的发热功率;温升曲线确定子模块220用于根据线圈随时间的温升确定温升曲线。其中,该装置的具体实施过程参见上述方法实施例中的说明即可,本实施例在此不再赘述。
本实施例中,通过公式可以计算得到线圈随时间的温升,进而可以确定温升曲线,通过观察温升曲线,可以更清晰的了解到线圈匝数对温升的影响。
参见图26,图26为本实施例提供的高速机械开关多场耦合仿真计算装置中,运动学分析模块的结构框图。如图所示,运动学分析模块300包括:导入子模块310、设置子模块320和仿真分析子模块330。其中,导入子模块用于将电磁斥力导入预先建立的计算模型中。设置子模块用于对高速机械开关的部件连接部分添加相关类型运动副,并设置部件的参数。仿真分析子模块用于对高速机械开关进行分闸反弹运动场耦合仿真分析,以获得位移曲线。其中,该装置的具体实施过程参见上述方法实施例中的说明即可,本实施例在此不再赘述。
本实施例中,通过将电磁斥力导入预先建立的计算模型中,并对高速机械开关进行分闸反弹运动场耦合仿真分析,可以获得位移曲线,通过位移曲线可以清晰地了解到线圈匝数对位移的影响。
综上,本实施例通过对高速机械开关进行电磁场仿真分析、热电耦合场仿真分析和分闸反弹运动场耦合仿真分析,可以全面反映高速机械开关运动、触头碰撞和线圈温升情况,从而可以全面准确反映高速机械开关的机械性能,以得到对高速机械开关最优化的设计,进而为混合式直流断路器快速机械开关优化设计提供了理论指导。
需要说明的是,高速机械开关多场耦合仿真计算方法和高速机械开关多场耦合仿真计算装置原理相同,相似之处可相互参照。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (4)

1.一种高速机械开关多场耦合仿真计算方法,其特征在于,包括如下步骤:
对高速机械开关进行电磁场仿真分析,以确定所述高速机械开关的瞬态线圈电流和斥力金属盘的电磁斥力;
根据所述瞬态线圈电流对所述高速机械开关进行热电耦合场仿真分析,以获得线圈的温升曲线;
根据所述电磁斥力对所述高速机械开关进行分闸反弹运动场耦合仿真分析,以获得所述高速机械开关的位移曲线;
上述确定所述电磁斥力的方法为:
根据方程组
Figure FDA0003268623270000011
Figure FDA0003268623270000012
Figure FDA0003268623270000013
确定所述电磁斥力Fc,其中公式(3)的初始条件为vt=0=v0=0,xt=0x0=0;上式中,ψq、ψg均为电磁系统全磁链,t为所述高速机械开关的动作时间,Lq为斥力线圈的自感,iq为脉冲电流,M为所述斥力线圈与所述斥力金属盘的互感,ig为感应电流,Uc为电容电压,Rq为所述斥力线圈的电阻,Lg为所述斥力金属盘的自感,Rg为所述斥力金属盘的电阻,Ff为所述斥力金属盘受到的运动反力,x为所述斥力金属盘的位移,v1f为高速机械开关运动部件碰撞后初速度,v2f为高速机械开关运动部件碰撞后终速度,v1i为高速机械开关运动部件碰撞前初速度,v2i为高速机械开关运动部件碰撞前终速度,e为碰撞恢复系数;
上述获得线圈的温升曲线的方法为:
根据公式
Figure FDA0003268623270000014
其中初始条件为
Figure FDA0003268623270000015
确定所述线圈随时间的温升;上式中,τ为斥力线圈的温升,T为热时间常数,t为所述高速机械开关的动作时间,e为碰撞恢复系数,c为所述线圈的材质的比热容,m为所述线圈的质量,KT为所述线圈的散热系数,A为所述线圈的散热面积,P为所述线圈的发热功率;
根据所述线圈随时间的温升确定所述温升曲线。
2.根据权利要求1所述的高速机械开关多场耦合仿真计算方法,其特征在于,上述获得所述高速机械开关的位移曲线的方法为:
将所述电磁斥力导入预先建立的计算模型中;
对所述高速机械开关的部件连接部分添加相关类型运动副,并设置所述部件的参数;
对所述高速机械开关进行分闸反弹运动场耦合仿真分析,以获得所述位移曲线;
所述电磁斥力通过AKISPL样条差值函数和DM距离函数导入所述计算模型中。
3.一种高速机械开关多场耦合仿真计算装置,利用权利要求1至2中任一项所述的高速机械开关多场耦合仿真计算方法,其特征在于,包括:
电磁场分析模块,用于对高速机械开关进行电磁场仿真分析,以确定所述高速机械开关的瞬态线圈电流和斥力金属盘的电磁斥力;
瞬态热分析模块,用于根据所述瞬态线圈电流对所述高速机械开关进行热电耦合场仿真分析,以获得线圈的温升曲线;
动力学分析模块,用于根据所述电磁斥力对所述高速机械开关进行分闸反弹运动场耦合仿真分析,以获得所述高速机械开关的位移曲线;
所述电磁场分析模块包括:
电磁斥力确定子模块,用于根据方程组
Figure FDA0003268623270000031
Figure FDA0003268623270000032
Figure FDA0003268623270000033
确定所述电磁斥力Fc,其中公式(3)的初始条件为vt=0=v0=0,xt=0x0=0;
上式中,ψq、ψg均为电磁系统全磁链,t为所述高速机械开关动作时间,Lq为斥力线圈的自感,iq为脉冲电流,M为所述斥力线圈与所述斥力金属盘的互感,ig为感应电流,Uc为电容电压,Rq为所述斥力线圈的电阻,Lg为所述斥力金属盘的自感,Rg为所述斥力金属盘的电阻,Ff为所述斥力金属盘受到的运动反力,x为所述斥力金属盘的位移,v1f为高速机械开关运动部件碰撞后初速度,v2f为高速机械开关运动部件碰撞后终速度,v1i为高速机械开关运动部件碰撞前初速度,v2i为高速机械开关运动部件碰撞前终速度,e为碰撞恢复系数;
所述瞬态热分析模块包括:
线圈温升确定子模块,用于根据公式
Figure FDA0003268623270000034
其中初始条件为
Figure FDA0003268623270000035
确定所述线圈随时间的温升;上式中,τ为斥力线圈的温升,T为热时间常数,t为所述高速机械开关的动作时间,e为碰撞恢复系数,c为所述线圈的材质的比热容,m为所述线圈的质量,KT为所述线圈的散热系数,A为所述线圈的散热面积,P为所述线圈的发热功率;
温升曲线确定子模块,用于根据所述线圈随时间的温升确定所述温升曲线。
4.根据权利要求3所述的高速机械开关多场耦合仿真计算装置,其特征在于,所述动力学分析模块包括:
导入子模块,用于将所述电磁斥力导入预先建立的计算模型中;
设置子模块,用于对所述高速机械开关的部件连接部分添加相关类型运动副,并设置所述部件的参数;
仿真分析子模块,用于对所述高速机械开关进行分闸反弹运动场耦合仿真分析,以获得所述位移曲线;
所述电磁斥力通过AKISPL样条差值函数和DM距离函数导入所述计算模型中。
CN201711131430.7A 2017-11-15 2017-11-15 高速机械开关多场耦合仿真计算方法及其计算装置 Active CN108009322B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711131430.7A CN108009322B (zh) 2017-11-15 2017-11-15 高速机械开关多场耦合仿真计算方法及其计算装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711131430.7A CN108009322B (zh) 2017-11-15 2017-11-15 高速机械开关多场耦合仿真计算方法及其计算装置

Publications (2)

Publication Number Publication Date
CN108009322A CN108009322A (zh) 2018-05-08
CN108009322B true CN108009322B (zh) 2021-11-16

Family

ID=62052415

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711131430.7A Active CN108009322B (zh) 2017-11-15 2017-11-15 高速机械开关多场耦合仿真计算方法及其计算装置

Country Status (1)

Country Link
CN (1) CN108009322B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104636554A (zh) * 2015-02-06 2015-05-20 福州大学 智能交直流电磁电器新型设计方法
CN106250626A (zh) * 2016-08-01 2016-12-21 华北电力大学 一种基于l‑i曲线的特高压变压器空载直流偏磁快速计算方法
CN106468889A (zh) * 2016-09-06 2017-03-01 平高集团有限公司 一种快速机械开关半实物仿真平台
CN106484928A (zh) * 2015-08-26 2017-03-08 北京卫星环境工程研究所 基于多软件联合的开关电源电热耦合仿真方法
CN106951664A (zh) * 2017-04-26 2017-07-14 福州大学 一种接触器的自动建模及硬件在环仿真方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITTO20130574A1 (it) * 2013-07-09 2015-01-10 St Microelectronics Srl Procedimento per la progettazione automatica di un circuito elettronico, relativo sistema e prodotto informatico

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104636554A (zh) * 2015-02-06 2015-05-20 福州大学 智能交直流电磁电器新型设计方法
CN106484928A (zh) * 2015-08-26 2017-03-08 北京卫星环境工程研究所 基于多软件联合的开关电源电热耦合仿真方法
CN106250626A (zh) * 2016-08-01 2016-12-21 华北电力大学 一种基于l‑i曲线的特高压变压器空载直流偏磁快速计算方法
CN106468889A (zh) * 2016-09-06 2017-03-01 平高集团有限公司 一种快速机械开关半实物仿真平台
CN106951664A (zh) * 2017-04-26 2017-07-14 福州大学 一种接触器的自动建模及硬件在环仿真方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
高速机械开关单元多场耦合仿真优化设计;田阳 等;《电网技术》;20180430;第42卷(第4期);第1273-1280页 *
高速电磁斥力机构运动特性影响因素研究;卢彦辉 等;《高压电器》;20150816;第51卷(第8期);第0029-0034页 *

Also Published As

Publication number Publication date
CN108009322A (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
CN108416169B (zh) 一种接触器触头系统带载多物理场耦合仿真优化设计系统
Li et al. 3D Numerical simulation method of electromagnetic forming for low conductive metals with a driver
Peng et al. Active damping of ultrafast mechanical switches for hybrid AC and DC circuit breakers
CN101164125A (zh) 电磁促动器上的无传感器位置识别
Fang et al. A novel flux-weakening control strategy for permanent-magnet actuator of vacuum circuit breaker
CN111079299B (zh) 容性负载下电磁继电器闭合弹跳电接触力学特性计算方法
CN108009322B (zh) 高速机械开关多场耦合仿真计算方法及其计算装置
Vilchis-Rodriguez et al. Modelling thomson coils with axis-symmetric problems: practical accuracy considerations
Li et al. An enhanced analytical calculation of the phase inductance of switched reluctance machines
Riba et al. Dynamic model for AC and DC contactors–simulation and experimental validation
Zhu et al. Vibration characteristics of thomson coil actuator based on simulation and experiments
Jiang et al. Optimization design of a permanent magnet actuator for 126-kV vacuum circuit breaker
Lee et al. Reduced modeling of eddy current-driven electromechanical system using conductor segmentation and circuit parameters extracted by FEA
CN104170041B (zh) 电路开关
CN115828812A (zh) 一种计及位移耦合的变压器短路电动力计算方法
Hou et al. Research on fast opening process in electromagnetic repulsion mechanism
CN107229788B (zh) 智能接触器虚拟样机设计平台
Zeng et al. Unsymmetrical bistable multimagnetic circuit permanent magnet actuator for high‐voltage circuit breaker application: analysis, design, and dynamic simulation
Joo et al. Dynamic analysis of linear electromagnetic solenoid for Electric Vehicle Relay
CN111695258B (zh) 一种电磁继电器动态特性仿真动能注入仿真方法
CN110232211B (zh) 一种考虑结构热场作用的接触器弹跳特性计算方法
Craciun et al. Multi-domain simulation and analysis of electromagnetically actuated reclosers
He et al. Multiphysics Modeling of the Dynamic Response of a Vacuum Automatic Circuit Recloser System
Mosallanejad Analyses of the effect of length difference between plunger and winding on tubular linear reluctance motor performance
Zhai et al. Study on optimization of electromagnetic relay's reaction torque characteristics based on adjusted parameters

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