一种大型汽轮发电机负载励磁电流的自动迭代计算方法
【技术领域】
[0001] 本发明涉及一种大型汽轮发电机励磁电流的数值迭代计算方法，可用于大型汽轮 发电机励磁电流的计算研究和发电机产品开发研究中，属于发电机电磁设计技术领域。
【背景技术】
[0002] 励磁电流是发电机的主要运行数据之一，其设计值的准确程度不但影响着励磁系 统的选型与性能参数，更重要的是影响着运行过程中转子绕组的实际温升情况，因此励磁 电流的准确计算成为发电机产品设计过程中的重要环节。
[0003] 在计算额定励磁电流的过程中，非线性所带来的影响是主要的误差来源，尤其是 对于大容量发电机。美国前西屋公司计算机分析方法和试验数据表明，额定电压下总的视 在空载饱和磁势MSNL可能高达磁路磁铁部分饱和磁势的160~180%。我国曾长期使用的 手算方法（按DZ28-63)，由于没有计及电枢反应增加的转子电流所产生的转子漏磁及其对 转子磁路饱和的影响，在大型汽轮发电机中实际的I fN值要比计算值大。
[0004] 西屋公司70年代初编写的一种空载饱和计算程序是将定、转子和气隙分成许多 小块形成平面磁网络，通过迭代取得磁通势平衡，计算出各小块区域的磁场数据和总的励 磁安匝，这实际上就是数值计算方法。西屋公司70年代开始尝试利用数值计算方法来确定 励磁电流，EXCITE程序的开发正是出于此目的，它采用矢量位网络法，计算得出电机工况所 需要的气隙磁通，从而确定励磁电流。多年来此程序也不断得到了改进和完善，如从有限差 分模式到有限元模式的修正等，但受到早期计算机条件的限制和一些算法上的局限性，此 程序在初期存在一定缺陷，计算精度不高，且计算上费时较多。针对四极发电机的励磁电流 计算，取而代之的是FE-EXCITE程序，该程序并非EXCITE程序的升级，而是基于有限元计算 方法开发的另一种程序，但该程序同样存在一定的不足，如发电机定子槽数多于66槽时该 程序就无法计算，因此在目前世界上最大单机容量（1800MW)的Fin5发电机励磁电流计算 过程中，也未能得到应用。Fin5发电机励磁电流是采用德国西门子公司开发的另一种有限 元程序来计算的，但是该程序存在的缺点是不能考虑发电机转子单边偏槽结构的励磁电流 计算。
[0005] 我国目前已规划并正在开展世界最大容量两极1200MW、四极核电2000MW水氢冷 发电机的研发工作，在大容量发电机以及百万千瓦级四极发电机自主化设计过程中，必须 研究新的适用程度更高、能够保证精度的励磁电流计算方法，以满足各种发电机结构（如 偏槽、浅槽、各种转子槽型等）、非线性材料和不同工况下的励磁电流的计算要求。
【发明内容】
[0006] 本发明的目的是提供一种能够适应各种槽型结构，可计算任意工况下的励磁电流 的方法。
[0007] 为了达到上述目的，本发明的技术方案是提供了一种大型汽轮发电机负载励磁电 流的自动迭代计算方法，其特征在于，包括以下步骤：
[0008] 步骤1、输入定子模型参数、转子模型参数及电气参数；
[0009] 步骤2、分别设定：定子电流与q轴的夹角λ的初值λ 〇及负载励磁电流If的初 值 IfO ;
[0010] 步骤3、求取气隙中心的磁密数据，包括以下步骤：
[0011] 步骤3. 1、依据步骤1输入的定子模型参数、转子模型参数及电气参数进行定子建 模及转子建模，建模时转子q轴与X轴重合；
[0012] 步骤3. 2、选用2维电磁单元对建立的定子模型及转子模型进行网格划分，将发电 机本体截面结构模型离散为数值模型，并在气隙中心设立路径线；
[0013] 步骤3. 3、定子线圈载荷加载：
[0014] 将0时刻电流Ia= I "os (-π +λ 〇)作为载荷加载到A相定子线圈上，式中，1"为 运行定子电流幅值，加载在B相定子线圈和C相定子线圈上的电流在加载在A相定子线圈 的电流的基础上相位分别增加120°和240° ;
[0015] 步骤3. 4、转子线圈载荷加载：
[0016] 根据转子不同槽内匝数Ν，将IfOXN加载到对应的转子槽线圈内；
[0017] 步骤3. 5、在定子铁心外圆处设定平行边界条件；
[0018] 步骤3. 6、设定求解方式为2D静态场，求解得到气隙中心的磁密数据；
[0019] 步骤4、依据步骤3得到的气隙中心的磁密数据，计算得到气隙中心基波磁密的幅 值Bni和相位Φ "，气隙磁密基波空间分布为B = B^Xco^P-ilO,式中，β为以q轴为参考 轴的空间电角度变量；
[0020] 步骤5、根据步骤3. 1的建模原则计算出定子A相线圈轴线与X轴的夹角为Φ ZA， A相气隙感应电压为
，式中，kw为定子绕组系数，w为定子绕组每相串 连匝数，Lrff为定子铁心有效长度，τ为气隙极距，f为频率，根据发电机电势相量关系可得 Εδ与1」司的夹角θ = π/2-λ〇-(φ ZA-i])m)，再由公式：
[0023] 计算出端电压U。和功率因数角私，式中，\为定子漏抗；
[0024] 步骤6、设定电压误差限值ε。和相位误差限值％，并制定循环策略如下：
[0025] 1)当U。〉UΝ+ ε u，妗> 、式中，Un为额定端电压，終^为额定功率因数角，此时 希望1]。和矜都减小，需要调整If变小；
[0026] 2)当Uc< U Ν_ ε u，％ >羚+% >此时希望I变大而钱减小，需要调整λ变大；
[0027] 3)当U。〉U Ν+ ε u，灼< 此时希望U。变小而_增大，需要调整λ变小；
[0028] 4)当Uc< U Ν_ ε u，涔 <終^ 一^，此时希望Uc和铃均增大，需要调整If变大。
[0029] 步骤7、将计算的端电压U。和功率因数角魏与额定端电压Un和额定功率因数角 进行比较，判断是否达到预设的误差限值ε。和力〃只有同时满足预设限值才达到迭代 要求，进入步骤8,否则根据步骤6的循环策略返回到步骤2对初值λ O和/或初值IfO进 行修改，再次计算；
[0030] 步骤8、得到最终的端电压U。和功率因数角鈐《
[0031] 优选地，在所述步骤3. 1中，建模时保证转子磁极中心d轴与定子A相线圈轴线重 合，转子q轴与X轴重合。
[0032] 优选地，在所述步骤2中，将夹角λ的取值范围设定为[As。，λΛ]，则λ〇 = Us(]+A J/2 ;将所述负载励磁电流If的取值范围设定为[Ifs，IfJ，则IfO= (Ifs+IfJ/2。
[0033] 本发明的优点是计算方法适用程度高，可以满足各种发电机结构（如偏槽、浅槽、 各种转子槽型等）、非线性材料和不同工况下的励磁电流的计算要求，帮助设计人员在方案 设计阶段准确获取励磁参数，从而指导发电机产品的优化设计；另外，本发明计算过程按功 能进行模块化设计，并通过循环策略的设置实现自动迭代寻解功能，具有功能模块清晰、计 算程序易于重复调用、易于阅读和升级、计算速度快、求解精度高等优点。
【附图说明】
[0034] 图1为本发明的算法流程图。
【具体实施方式】
[0035] 为使本发明更明显易懂，兹以优选实施例，并配合附图作详细说明如下。
[0036] 结合图1，本发明提供了一种大型汽轮发电机负载励磁电流的自动迭代计算方法， 其步骤为：
[0037] 步骤1、模型参数及相关电气参数的输入。
[0038] 按照发电机定、转子本体截面结构及额定工况，输入模型参数及电气参数，作为有 限元建模及迭代计算的输入参数，包括定子模型参数、转子模型参数、额定电压、额定功率、 功率因数、频率等参数，同时还包括工作目录设置等，该一系列参数可通过编制MATlAB界 面程序进行规范化输入，由输入数据模块程序实现。
[0039] 表1为某电厂大型汽轮发电机实际运行参数，在本发明中，编制了界面化的发电 机定、转子结构参数及电气参数输入模块，所有建模数据及电气参数的输入由该模块实现。
[0041] 表1某电厂实测运行参数
[0042] 步骤2、文件预处理。
[0043] 输入、输出数据采用文件形式进行存储和读写，便于不同模块间的数据共享和 传输，也便于计算结果的标准化显示。根据界面化输入中设置的工作目录路径建立相应 的文件夹D :\Excitation_File_Ansys，并将上一步输入的定子模型参数、转子模型参数 以及电气参数以文件的形式分别保存于所建目录下，数据文件分别为Stat〇rPara.tXt、 RotorPara. txt 和 AnPara. txt。该项功能由 Matlab 子程序 PretreatDoc ()完成。
[0044] 步骤3、负载励磁电流计算模块
[0045] 计算发电机在负载状态下的励磁电流大小，采用ANSYS进行气隙磁密数值计算， 由MATLAB进行循环策略的控制，判断是否满足计算精度要求，按循环策略改变迭代参数， 驱动ANSYS进行自动迭代计算，直到满足精度要求为止，具体由RatedPFLoad ()模块程序实 现。具体步骤为：
[0046] 步骤3. 1、定子电流与q轴的夹角λ的初值λ 〇,夹角λ是考虑转子偏槽影响后 的综合夹角，首先设定夹角λ的范围为[λ s。，λ J，该范围能够保证λ的收敛性，取初始 值λ0 = (λ^+λ J/2。在本实施例子中，夹角λ的范围取为[0°，90° ]，则初值λ〇 = 45。。
[0047] 步骤3. 2、设定负载励磁电流If的初值IfO,后续迭代采用二分法，为了迭代过程 的实现，也首先设