本发明是一种采用变阶变步长3S‑DIRK算法的电磁暂态仿真方法，其特点是，具有4个分算法，在正常仿真计算时采用算法A(步骤0)，当检测到开关动作需要插值时，采用算法D进入(步骤1)，并将算法切换到算法B(步骤2)，等到数值振荡消除后采用算法C(步骤3)同步后切换到算法A(步骤0)，若在步骤2阶段发生插值时，则进入步骤1插值，4个分算法具有不同的优点，可以用在电磁暂态的不同计算情形，本发明在切换算法时可以保证元件的等值导纳不变，整个仿真过程中计算精度不低于2阶，且故障期间具有L稳定可以消除数值振荡，并且可以变步长计算，能够应用于大规模电网的电磁暂态仿真计算。