[发明专利]汽轮发电机组碰摩故障的弯扭耦合振动特性分析方法

权利要求书

1.一种汽轮发电机组碰摩故障的弯扭耦合振动特性分析方法，其特征是所述方法包括：

步骤1：设定初始时刻t₀和步进长度Δt；

步骤2：采集时刻t₀机组各节点的位移和速度，计算时刻t₀机组各节点的加速度;

步骤3：令θ＝0；

步骤4：根据时刻t₀+θΔt机组各节点的位移、速度和加速度，计算时刻t₀+θΔt机组各轴段的外力/外力矩；

步骤5：将时刻t₀+θΔt机组各轴段的外力/外力矩分别作为时刻t₀+(θ+1)Δt机组各轴段的临时外力/临时外力矩；

步骤6：根据时刻t₀+(θ+1)Δt机组各轴段的临时外力/临时外力矩，利用基于Riccati变换的轴系弯扭耦合振动增量传递矩阵计算时刻t₀+(θ+1)Δt机组各节点的位移增量；

步骤7：根据时刻t₀+(θ+1)Δt机组各节点的位移增量，利用改进的wilson-θ增量表达式计算时刻t₀+(θ+1)Δt的位移、速度和加速度；

步骤8：根据时刻t₀+(θ+1)Δt的位移、速度和加速度，计算时刻t₀+(θ+1)Δt机组各轴段的外力/外力矩；

步骤9：判断是否满足设定条件，如果满足设定条件，则执行步骤11；否则，执行步骤10；

所述设定条件为：时刻t₀+(θ+1)Δt机组各轴段的外力/外力矩与时刻t₀+(θ+1)Δt机组各轴段的临时外力/临时外力矩的差值绝对值小于设定值；

步骤10：令时刻t₀+(θ+1)Δt机组各轴段的临时外力/临时外力矩等于时刻t₀+(θ+1)Δt机组各轴段的外力/外力矩，返回步骤6；

步骤11：令θ=θ+1，返回步骤5。

2.根据权利要求1所述的方法，其特征是所述计算时刻t₀+θΔt机组各轴段的外力/外力矩具体为，当轴段为碰摩轴段时，计算时刻t₀+θΔt机组该轴段的摩擦力矩；当轴段为轴颈轴段时，计算时刻t₀+θΔt机组该轴段的油膜力；当轴段为动叶栅所在轴段时，计算时刻t₀+θΔt机组动叶栅受到的蒸汽弯曲力矩；当轴段为发电机绕组轴段时；计算时刻t₀+θΔt机组的发电机绕组轴段受到的电磁力矩。