1.一种考虑多重事件发生的电磁暂态仿真方法，其特征在于，包括如下步骤：步骤1：全局采用隐式梯形积分法计算至t+Δt时刻，搜索是否有事件发生；其中，所述事件为电力电子元件动作；步骤2：若检测到事件发生，通过线性插值算法确定事件发生时间t a ，并修改相关状态变量值，重新形成系统节点导纳矩阵；若检测不到事件发生，则跳转至步骤6；步骤3：在t a 时刻，根据修改后的相关状态变量值和重新形成的节点导纳矩阵，采用半步长后退欧拉法和换路定理，进行初始化过程；步骤4：在t a +Δt/2时刻，搜索是否有其他事件发生，若检测到其他事件发生，返回步骤2，否则执行步骤5；步骤5：在t a 时刻，采用带阻尼的隐式梯形积分法和半步长后退欧拉法进行重新同步化过程；步骤6：计算完成后，继续采用隐式梯形积分法进行下一步长计算；所述带阻尼的隐式梯形积分法为：其中，a为阻尼系数，t为时间，Δt为仿真步长，f(t)为非状态变量，y(t)为状态变量；当a＝0时，积分方法变为纯隐式梯形积分方法；当a＝1时，积分方法变为后退欧拉法；当0<a<1时，积分方法介于梯形积分法和后退欧拉法之间；所述采用带阻尼的隐式梯形积分法和半步长后退欧拉法进行重新同步化过程，具体包括：当0≤t+Δt-t a ≤Δt/2时，仿真步长变为Δt′，根据以下公式计算仿真步长Δt′以及阻尼系数a；Δt′＝(t+Δt-t a )+Δt/2；其中，t a 为事件发生的时刻；根据所述阻尼系数a，按照所述仿真步长Δt′进行仿真，然后进行一次半步长后退欧拉法，完成同步过程；若当Δt/2≤t+Δt-t a ≤Δt，仿真步长变为Δt′，根据以下公式计算仿真步长Δt′以及阻尼系数a；Δt′＝(t+Δt-t a )其中，t a 为事件发生的时刻；根据所述阻尼系数a，按照所述仿真步长Δt′进行仿真，然后进行两次半步长后退欧拉法，完成同步过程。
2.一种考虑多重事件发生的电磁暂态仿真装置，其特征在于，包括搜索模块，用于全局采用隐式梯形积分法计算至t+Δt时刻，搜索是否有事件发生；其中，所述事件为电力电子元件动作；线性插值模块，用于若检测到事件发生，通过线性插值算法确定事件发生时间t a ，并修改相关状态变量值，重新形成系统节点导纳矩阵；若检测不到事件发生，则跳转至执行下一步长计算模块；初始化模块，用于在t a 时刻，根据修改后的相关状态变量值和重新形成的节点导纳矩阵，采用半步长后退欧拉法和换路定理，进行初始化过程；所述搜索模块，还用于在t a +Δt/2时刻，搜索是否有其他事件发生，若检测到其他事件发生，返回执行所述线性插值模块，否则执行重新同步化过程模块；重新同步化过程模块，用于在t a 时刻，采用带阻尼的隐式梯形积分法和半步长后退欧拉法进行重新同步化过程；下一步长计算模块，用于计算完成后，继续采用隐式梯形积分法进行下一步长计算；所述带阻尼的隐式梯形积分法为：其中，a为阻尼系数，t为时间，Δt为仿真步长，f(t)为非状态变量，y(t)为状态变量；其中，当a＝0时，积分方法变为纯隐式梯形积分方法；当a＝1时，积分方法变为后退欧拉法；当0<a<1时，积分方法介于梯形积分法和后退欧拉法之间；所述采用带阻尼的隐式梯形积分法和半步长后退欧拉法进行重新同步化过程，具体包括：当0≤t+Δt-t a ≤Δt/2时，仿真步长变为Δt′，根据以下公式计算仿真步长Δt′以及阻尼系数a；Δt′＝(t+Δt-t a )+Δt/2；其中，t a 为事件发生的时刻；根据所述阻尼系数a，按照所述仿真步长Δt′进行仿真，然后进行一次半步长后退欧拉法，完成同步过程；当Δt/2≤t+Δt-t a ≤Δt时，仿真步长变为Δt′，根据以下公式计算仿真步长Δt′以及阻尼系数a；Δt′＝(t+Δt-t a )其中，t a 为事件发生的时刻，根据所述阻尼系数a，按照所述仿真步长Δt′进行仿真，然后进行两次半步长后退欧拉法，完成同步过程。
3.一种设备，其特征在于，包括处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1所述的考虑多重事件发生的电磁暂态仿真方法。
4.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如权利要求1所述的考虑多重事件发生的电磁暂态仿真方法。
