[发明专利]基于受扰轨迹动态特性的数值积分提前终止方法

权利要求书

1.基于受扰轨迹动态特性的数值积分提前终止方法，其特征在于，所述方法包括：根据受扰轨迹信息，结合包括：失稳模式易变度、受扰轨迹非同调度、最大转子角间隙、同向相邻摆次转子角间隙信息相似度及振荡衰减度在内的动态特性，判断是否提前终止数值积分；

所述方法具体包括如下步骤：

对测试算例全集中的每个算例分别执行下述步骤：

步骤A：求取算例的失稳模式易变度K；

步骤B：采用小步长逐步积分方法向前积分一步，直至到达主导映象首摆最远点或动态鞍点；

步骤C：若采用小步长逐步积分已达主导映象首摆最远点，则根据下述判定条件判定是否提前终止积分：

1)失稳模式易变度K值小于第一阈值ε₁；

2)主导映象机械功率P_m大于电磁功率分量P_c；

若同时满足上述两个判定条件则提前终止积分；否则，执行下一步；

步骤D：若采用小步长逐步积分已达主导映象首摆动态鞍点，则根据下述判定条件判定是否提前终止积分：

1)失稳模式易变度K值小于第一阈值ε₁；

若满足上述判定条件则提前终止积分；否则，执行下一步；

步骤E：基于首摆暂态过程受扰轨迹，求取首摆非同调度μ₁；

步骤F：若采用小步长逐步积分已达主导映象首摆最远点，则根据下述判定条件判定是否提前终止积分：

1)失稳模式易变度K值小于第二阈值ε₂；

2)主导映象机械功率P_m大于电磁功率分量P_c；

3)首摆非同调度μ₁值小于第三阈值ε₃；

若同时满足上述三个判定条件则提前终止积分；否则，执行下一步；

步骤G：若采用小步长逐步积分已达主导映象首摆动态鞍点，则根据下述判定条件判定是否提前终止积分：

1)失稳模式易变度K值小于第二阈值ε₂；

2)首摆非同调度μ₁值小于第三阈值ε₃；

若同时满足上述两个判定条件则提前终止积分；否则，执行下一步；

步骤H：采用小步长逐步积分方法向前积分一步，直至最大转子角间隙大于预设的第四阈值ε₄，或到达主导映象第三摆最远点或动态鞍点；

步骤I：若最大转子角间隙已大于第四阈值ε₄，则提前终止积分；否则，执行下一步；其中：δ_ij表示机组i转子角同机组j转子角的差值，即：δ_ij＝δ_i-δ_j；i表示领前群中第i个机组；j表示余下群中第j个机组；S表示领前群；A表示余下群；

步骤J：若采用小步长逐步积分已达主导映象第三摆最远点或动态鞍点，求取前三摆非同调度μ_1-3；若前三摆非同调度μ_1-3值小于第五阈值ε₅，则提前终止积分；否则，执行下一步；

步骤K：采用小步长逐步积分方法向前积分一步，直至最大转子角间隙大于预设的第四阈值ε₄，或到达主导映象下一摆最远点；

步骤L：若最大转子角间隙已大于第四阈值ε₄，则提前终止积分；否则，执行下一步；

步骤M：若采用小步长逐步积分已达主导映象下一摆最远点，从该摆次始，向后溯洄四个摆次，若已小于第四摆，执行步骤K，否则执行下一步；

步骤N：根据步骤M中向后溯洄四个摆次的暂态过程，求取正向及反向相邻摆次转子角间隙信息相似度，若相似度小于第六阈值ε₆，则提前终止积分；否则，执行下一步；

步骤O：若从该摆次始，向后溯洄六个摆次，若已小于第四摆，则执行步骤K，否则执行下一步；

步骤P：根据步骤O中向后溯洄的六个摆次的暂态过程，求取连续两个正向及反向相邻摆次振荡衰减度，若振荡衰减度小于第七阈值ε₇，则提前终止积分；否则，执行步骤K。

2.根据权利要求1所述的基于受扰轨迹动态特性的数值积分提前终止方法，其特征在于，所述失稳模式易变度K采用大步长泰勒级数展开快速求取。