C3同步发电机的基本方程电力系统湖南大学

C3同步发电机的基本方程湖南大学·电力系统分析课程Contents本章内容框架同步发电机的数学模型与分析方法，从物理建模到坐标变换再到等效电路分析。01同步发电机的物理模型与理想化条件02同步发电机的原始方程03派克变换与dq0坐标系04dq0坐标系下的基本方程05等效电路与相量图分析06转子运动方程与模型阶数选择CHAPTER01同步发电机的物理模型与理想化条件从实际电机结构出发，建立可用于数学建模的理想化模型ElectromagneticModel同步发电机的绕组结构同步发电机的电磁模型由定子三相绕组（a、b、c）和转子绕组（励磁绕组f、d轴阻尼绕组D、q轴阻尼绕组Q）共同构成，各绕组的空间位置关系决定了互感系数的周期性时变特征。同步发电机定子绕组与转子绕组空间结构示意01定子三相绕组a、b、c在空间上互差120°电角度对称分布，各相绕组的轴线方向即为该相磁动势的作用方向02转子励磁绕组f沿d轴方向布置，产生主磁场的直流励磁磁动势，是同步发电机建立气隙磁场的核心来源03d轴阻尼绕组D和q轴阻尼绕组Q分别等效转子在直轴和交轴方向的涡流阻尼效应，反映转子铁芯和槽楔中的感应电流04定转子绕组之间存在随转子位置周期性变化的互感关系，这是同步发电机方程呈现时变系数的根本原因ASSUMPTIONS建立数学模型的理想化假设同步发电机基本方程的建立依赖于磁路线性、气隙均匀和正弦感应电势三大核心假设，这些假设将复杂的电磁场问题简化为可用常系数或周期性时变系数微分方程描述的集总参数电路问题。磁路线性假设忽略铁芯磁饱和、磁滞和涡流效应，认为磁导率为常数，各绕组电感参数与电流大小无关叠加原理成立，各绕组磁链可线性叠加，大幅简化多绕组耦合电磁关系线性叠加气隙与结构假设定转子槽口和通风沟不影响电感参数，气隙磁导仅与空间角度有关凸极效应通过直轴与交轴磁阻差异体现，隐极机可简化为两轴参数相等均匀气隙电势波形假设空载且转子恒速旋转时，定子绕组感应的空载电势为时间的正弦函数气隙磁密沿圆周按正弦分布，互感系数可用余弦函数精确描述正弦分布同步发电机·基本方程参考正方向的选取规则正方向约定是同步发电机方程推导的基础性前提，直接决定电压方程和磁链方程中各项的正负号，采用电动机惯例时需确保电流、电压、磁链三者的正方向内在自洽。dq旋转坐标系d轴取转子励磁绕组中心轴方向，q轴沿转子旋转方向超前d轴90°电角度，构成dq旋转坐标系。该坐标系随转子同步旋转，是分析凸极效应和电磁转矩的关键数学工具。90°电角度磁链正方向各绕组电流正方向使其产生的磁链与绕组轴线正方向一致，遵循右手螺旋定则。正向电流产生正向磁链，确保电感矩阵中互感系数为正值，简化方程表达。RIGHT-HANDRULE定子电流与电压定子电流以中性点流向机端为正，定子电压以电流流出端为正，即机端电压等于感应电势减去阻抗压降。此约定与电网分析中的惯例保持一致，便于互联系统计算。U=E−IZ转子励磁与阻尼转子励磁电压以提供正向励磁电流的方向为正，阻尼绕组端电压恒为零，视为短路绕组。阻尼绕组在暂态过程中感应电流，抑制转子振荡，其短路特性是稳定分析的重要边界条件。Ud=0CHAPTER02同步发电机的原始方程在abc静止坐标系下建立电压方程与磁链方程，揭示时变系数的本质MotorModel·Voltage定子电压方程定子三相电压方程由法拉第电磁感应定律和欧姆定律导出，表现为磁链时间导数与电阻压降的叠加关系，其核心困难在于磁链本身包含随转子位置周期性变化的互感系数。01由电磁感应定律e=−dψ/dt，定子各相端电压等于感应电势减去电阻压降：ua=−dψa/dt−ia·RaInductionLaw02三相电压方程写成矩阵形式，Rs为定子三相电阻对角矩阵：uabc=−dψabc/dt−Rs·iabcMatrixForm03磁链ψ包含自感磁链和互感磁链，而定转子互感随转子角度θ周期性变化，导致方程系数为时变函数。这种时变特性使得直接求解三相电压方程变得极为复杂，需要引入坐标变换来简化计算。Time-Varying04时变系数使得方程组求解极为困难，这是引入坐标变换（Park变换）的根本动因。通过将三相静止坐标系转换为两相旋转坐标系，可将时变系数转换为常数系数，从而大大简化电机模型的分析和控制设计。ParkTransformCHAPTER03·VOLTAGEEQUATIONS转子电压方程转子三个绕组的电压方程中，励磁绕组具有外加直流励磁电压，而d轴和q轴阻尼绕组均为短路绕组、端电压恒为零，转子方程与定子方程联立构成六绕组完整电压方程组。01FIELD励磁绕组电压方程uf=dψf/dt+if·Rf，其中uf为外加直流励磁电压，Rf为励磁绕组电阻uf直流励磁02D-AXISd轴阻尼绕组（短路）端电压恒为零：0=dψD/dt+iD·RD，RD为d轴阻尼绕组等效电阻0V短路绕组03Q-AXISq轴阻尼绕组（短路）同样为短路绕组：0=dψQ/dt+iQ·RQ，RQ为q轴阻尼绕组等效电阻0V短路绕组04SYSTEM定转子联立完整方程组转子方程与定子方程联立，构成6×6阶完整电压方程组，以分块矩阵形式清晰表达定转子耦合关系6×6分块矩阵MAGNETICFLUXEQUATION磁链方程的矩阵形式六个绕组的磁链方程构成6×6电感矩阵，对角线元素为各绕组自感，非对角线元素为绕组间互感；其中定转子间互感是转子位置角θ的周期函数，导致整个电感矩阵为时变矩阵。OVERVIEW矩阵形式ψ=L·iψ和i均为6维列向量，L为6×6电感矩阵，包含36个电感参数6×6电感矩阵SELFIND.定子自感Laa、Lbb、Lcc含恒定分量和二次谐波分量，是θ的函数Laa=Ls0+Ls2·cos2θMUTUALIND.定子相间互感Mab、Mbc、Mca含θ函数项，反映三相空间耦合关系Mab=−Ms0+Ls2·cos(2θ+60°)S-RCOUPLING定转子互感Maf、MaD、MaQ随转子旋转周期性变化，cosθ与sinθ交替Maf=Maf·cosθInductanceAnalysis电感参数的时变特征分析同步发电机电感参数的时变性源于转子旋转导致的定转子相对位置周期性变化，凸极效应使定子自感含2θ谐波、定转子互感含θ基波分量，这是方程求解困难的根本原因。01定子自感含恒定分量和二次谐波分量：Laa=Ls0+Ls2·cos2θ，其中Ls0=(Ld+Lq)/2，Ls2=(Ld−Lq)/2cos2θ02定子相间互感同样含θ函数项：Mab=−Ms0+Ls2·cos(2θ+π/3)，反映120°空间分布带来的耦合调制120°03定转子互感随θ按正弦规律变化：Maf=Maf_max·cosθ，MaQ=−MaQ_max·sinθ，当轴线重合时耦合最强sin/cos04转子绕组互感与自感与θ无关：如MfD、Lff、LDD，因转子绕组随转子一起旋转，相对位置固定ConstantEQUATIONANALYSIS原始方程的局限性与坐标变换的必要性abc坐标系下的同步发电机方程组因电感矩阵含时变系数而成为变系数微分方程组，既无法获得解析解也难以直接数值求解，必须通过坐标变换将时变系统转化为常系数系统。0112阶变系数微分方程组原始方程组包含6个电压方程和6个磁链方程，联立后形成12阶变系数微分方程组，求解难度极大12阶02电感矩阵的θ依赖电感矩阵中约一半元素是转子位置角θ的函数（含cosθ、cos2θ项），而θ(t)由转子运动方程决定cosθ03线性工具失效变系数微分方程组无法直接应用拉普拉斯变换等线性系统分析工具，暂态响应的解析求解几乎不可能Laplace04Park变换：从时变到常系数Park变换将定子三相量投影到随转子同步旋转的dq0坐标系，消去电感矩阵中的θ依赖项，使方程变为常系数dq0CHAPTER03派克变换与dq0坐标系将定子三相量从静止坐标系映射到旋转坐标系，消除时变系数MotorTheory·CoordinateTransformationPark变换的物理意义Park变换的物理本质是观察者视角的切换——从静止的定子坐标系切换到随转子同步旋转的坐标系，使定子三相合成磁动势在dq轴上的投影变为直流量。旋转磁动势定子三相电流产生的合成磁动势在空间中旋转，当转速等于同步速时，在转子上观察该磁动势是静止的ω=ωs直轴分量id反映定子磁动势沿转子励磁方向的投影，决定电枢反应的去磁或助磁效应去磁/助磁交轴分量iq反映定子磁动势沿交轴方向的投影，与转子磁场正交，主要影响电磁转矩输出电磁转矩零序分量i0三相对称运行时为零，仅在不对称故障或谐波分析中需要考虑i₀=0PARKTRANSFORMPark变换矩阵的数学表达Park变换矩阵P通过三角函数将定子三相量投影到旋转的dq0轴上，变换矩阵的元素是转子位置角θ的三角函数，系数2/3保证等幅值特性，变换前后功率保持不变。d轴变换行fd=(2/3)[fa·cosθ+fb·cos(θ−2π/3)+fc·cos(θ+2π/3)]将三相量投影到d轴方向，cosθ表示d轴与a轴夹角的余弦投影系数，2/3系数确保变换前后幅值相等。cosθ·2/3q轴变换行fq=−(2/3)[fa·sinθ+fb·sin(θ−2π/3)+fc·sin(θ+2π/3)]将三相量投影到q轴方向，sinθ表示q轴与a轴夹角的正弦投影系数，负号保证q轴超前d轴90°的正向旋转。sinθ·2/3零序变换行f0=(1/3)(fa+fb+fc)反映三相不对称程度，三相对称时f₀恒为零。1/3系数保证零序分量与原始三相量具有相同的物理量纲，便于故障分析。f₀=0（对称）逆变换矩阵P⁻¹可将dq0量还原为abc三相量，变换与逆变换均为可逆线性变换，保证信息无损。逆变换矩阵由正变换矩阵转置并调整系数得到，满足P·P⁻¹=I单位矩阵条件。[fafbfc]ᵀ=P⁻¹·[fdfqf0]ᵀP⁻¹=(3/2)·PᵀMATHEMATICALPROPERTIESPark变换的关键数学性质Park变换是可逆线性变换、功率守恒变换，其核心价值在于将定转子间随θ时变的互感全部转化为与θ无关的常系数，使电感矩阵从时变矩阵变为常数矩阵。可逆性P·P⁻¹=I，变换前后信息完全等价，可通过逆变换精确还原abc三相量，保证分析的完备性。逆变换矩阵与正变换互为转置，确保数学上的严格对称。P·P⁻¹=I常系数化变换后定子等效d、q绕组与转子绕组间的互感不再含θ，电感矩阵由时变退化为常数矩阵。这一特性是同步电机建模的关键突破，使微分方程求解成为可能。θ→CONST功率守恒变换前后系统功率平衡不变，三相功率之和等于dq0分量功率之和（含系数修正）。瞬时功率与平均功率均保持守恒，满足能量守恒定律。ΣP=CONST等幅值特性系数2/3保证三相对称稳态下dq分量幅值等于相量幅值，便于工程计算和物理量理解。该特性使标幺值系统构建更加直观，简化电机参数整定。2/3COEFFInductanceMatrixPark变换后的电感矩阵Park变换将6×6时变电感矩阵转化为常系数矩阵，定子侧表现为直轴同步电感Ld、交轴同步电感Lq和零序电感L0三个独立参数，定转子互感全部变为与θ无关的常数。01直轴同步电感Ld定子d轴自感变为Ld=Ls0+Ms0+3Ls2/2，综合了漏感、互感恒定分量和凸极效应分量。Ld=Ls0+Ms0+3Ls2/202交轴同步电感Lq定子q轴自感变为Lq=Ls0+Ms0−3Ls2/2，凸极机中Ld>Lq，隐极机中Ld≈Lq。Lq=Ls0+Ms0−3Ls2/203定转子互感常数Mad定子与转子d轴绕组互感变为常数Mad=√(3/2)·Maf_max，与θ完全无关，反映d轴方向的恒定耦合。Mad=√(3/2)·Maf04矩阵分块解耦变换后电感矩阵呈分块结构：d轴系统(d,f,D)与q轴系统(q,Q)解耦，零序系统独立，大幅简化分析。d·f·D⊥q·Q⊥0CHAPTER04dq0坐标系下的基本方程建立常系数微分方程组，形成完整的同步发电机数学模型MotorControl·VoltageEquationdq0坐标系下的定子电压方程Park变换后的定子电压方程在磁链变化率项之外引入了旋转电势项（-ωψq和+ωψd），这是坐标变换的数学结果，物理上对应转子旋转导致d轴和q轴磁链交叉感应的电动势。01d轴电压方程ud=dψd/dt+Rs·id−ω·ψq−ω·ψq为旋转电势，反映q轴磁链随转子旋转在d轴上的感应电动势。该项源于坐标系旋转，是dq变换引入的耦合项，体现了旋转磁场对静止绕组的相对运动效应。−ωψq02q轴电压方程uq=dψq/dt+Rs·iq+ω·ψd+ω·ψd为旋转电势，反映d轴磁链在q轴上的交叉感应。与d轴方程对称，符号为正源于右手定则和dq坐标系的正方向定义，共同构成旋转电势的完整表达。+ωψd03零序电压方程u0=dψ0/dt+Rs·i0零序分量与dq轴完全解耦，方程形式最简单，不含旋转电势项。在三相平衡对称工况下恒为零，仅在绕组不对称、电源不平衡或存在中性点连接时才需考虑。Decoupled04符号约定与惯例采用电动机惯例时，旋转电势项的符号取决于正方向约定。发电机惯例下各项符号会相应改变，电流方向定义反转。分析时必须首先明确所采用的惯例体系，避免符号混淆导致计算错误。Motor/GenFLUXEQUATIONSdq0坐标系下的磁链方程变换后的磁链方程呈清晰的轴解耦结构：d轴磁链仅与d轴电流、励磁电流和d轴阻尼电流有关，q轴磁链仅与q轴电流和q轴阻尼电流有关，零序磁链独立，三组方程可分别求解。d轴磁链Ld·Madψd=Ld·id+Mad·if+Mad·iDLd为直轴同步电感，Mad为定子与d轴转子绕组的恒定互感q轴磁链Lq·Maqψq=Lq·iq+Maq·iQLq为交轴同步电感，Maq为定子与q轴阻尼绕组的恒定互感零序磁链L₀ψ0=L0·i0L0为定子零序电感，仅取决于定子漏感分量，与转子参数完全无关轴解耦特性3+2+1d轴与q轴系统之间无互感耦合，6绕组问题分解为独立的d轴子系统（3绕组）、q轴子系统（2绕组）和零序分量，三组方程可分别独立求解EquationsSummarydq0坐标系下完整方程组汇总dq0坐标系下的同步发电机模型由定子电压方程（含旋转电势项）、磁链方程（常系数轴解耦矩阵）、转子电压方程和转子运动方程共同构成，是一组可求解的常系数非线性微分方程组。定子电压方程ud=dψd/dt+Rs·id−ω·ψqd-axis·Cross-couplinguq=dψq/dt+Rs·iq+ω·ψdq-axis·Cross-couplingu0=dψ0/dt+Rs·i0Zero-seq·Decoupled磁链方程与转子方程ψd=Ld·id+Mad·if+Mad·iDψq=Lq·iq+Maq·iQψ0=L0·i0FluxLinkage·DecoupledAxesuf=dψf/dt+Rf·if0=dψD/dt+RD·iD0=dψQ/dt+RQ·iQRotorVoltage·DamperWindingsTj·dω/dt=Tm−Te其中Te为电磁转矩SwingEquation·TmvsTeSteady-StateAnalysis同步稳态下的方程简化同步稳态运行时dq坐标系中各量均为直流量，磁链变化率项dψ/dt为零，微分方程退化为代数方程组，结合忽略定子电阻的工程近似后可直接用于潮流计算和稳态分析。01稳态条件下dψd/dt=dψq/dt=0，电压方程退化为：ud=Rs·id−ωψq，uq=Rs·iq+ωψddψ/dt=002同步转速下ω=1(pu)，进一步忽略Rs（Rs远小于同步电抗Xd），简化为：ud≈−Xq·iq，uq≈Eq−Xd·idω=1pu03简化后的代数方程可直接写成复数形式：U=Eq−jXd·Id（隐极机），便于与电力系统网络方程联立求解复数形式04稳态简化模型是电力系统潮流计算的基础，发电机在潮流中以PV节点或平衡节点形式参与计算PV节点CHAPTER05等效电路与相量图分析将数学方程转化为电路模型和几何图形，建立工程直觉EQUIVALENTCIRCUITd轴等效电路d轴等效电路将定子d轴绕组、励磁绕组和d轴阻尼绕组的电磁耦合用三个互感耦合的支路表示，公共互感Mad反映三个绕组共享的d轴主磁通路径，是分析直轴暂态和次暂态过程的基础。定子d轴支路电阻Rs串联漏感Lσ再串联公共互感Mad，端电压为ud，电流为idRs→Lσ→Mad励磁绕组支路电阻Rf串联漏感Lσf再串联公共互感Mad，外加直流电压uf，电流为ifRf→Lσf→Madd轴阻尼绕组支路电阻RD串联漏感LσD再串联公共互感Mad，端电压为零（短路），电流为iDRD→LσD→Mad稳态等效励磁效应等效为内电势Eq，d轴电路简化为Eq串联Xd（直轴同步电抗），便于与网络方程联立Eq+XdEQUIVALENTCIRCUITq轴等效电路q轴等效电路由定子q轴绕组和q轴阻尼绕组两支路通过互感Maq耦合构成，比d轴少一个励磁支路；稳态下简化为交轴同步电抗Xq，凸极机中Xq<Xd反映交轴气隙磁阻更大。定子q轴支路电阻Rs串联漏感Lσ再串联公共互感Maq，端电压为uq，电流为iq，构成q轴等效电路的定子侧主支路uq·iqq轴阻尼绕组支路电阻RQ串联漏感LσQ再串联公共互感Maq，端电压为零（短路绕组），电流为iQ，用于阻尼动态过程中的振荡RQ+LσQ无励磁绕组简化q轴无励磁绕组，结构比d轴简单，稳态下等效为uq=Xq·id+Rs·iq，忽略Rs后简化为uq≈Xq·iduq≈Xq·id凸极机参数特征凸极机Xd与Xq的差异程度（Xd/Xq比值）是区分水轮机和汽轮发电机的重要参数特征，水轮机该比值通常更大Xd/XqPHASORDIAGRAM同步发电机相量图相量图将dq坐标系下的电压方程转化为直观的几何关系，通过内电势Eq、端电压U、电流I及其dq分量之间的矢量运算，清晰展示功率角δ、功率因数和电压调节的工程含义。01q轴为参考方向，端电压U分解为Ud=U·sinδ和Uq=U·cosδ，其中δ为功率角（Eq超前U的角度）02电流I分解为Id=I·sin(δ+φ)和Iq=I·cos(δ+φ)，φ为功率因数角，反映有功和无功的分配比例03凸极机电压方程：U=Eq−jXd·Id−jXq·Iq，在相量图上表现为Eq到U的矢量减法路径04功角特性：P=(Eq·U/Xd)·sinδ+U²·(Xd−Xq)/(2XdXq)·sin2δ，第一项为励磁功率，第二项为磁阻功率同步发电机相量图·dq坐标系下电压与电流的矢量关系ELECTROMAGNETICTORQUE电磁转矩表达式及其物理组成同步发电机的电磁转矩由励磁转矩和磁阻转矩两部分组成，前者由励磁磁场与q轴电流的相互作用产生，后者由凸极效应导致的直轴与交轴磁阻差异产生，隐极机中磁阻转矩为零。01电磁转矩基本表达式：Te=ψd·iq−ψq·id，由d轴和q轴的磁链与交叉轴电流的乘积之差构成Te=ψd·iq−ψq·id02展开包含励磁分量Mad·if·iq、磁阻分量(Ld−Lq)·id·iq及阻尼绕组分量Mad·iD·iq−Maq·iQ·id励磁+磁阻+阻尼03励磁转矩Mad·if·iq是输出有功的主力，与励磁电流和q轴电流成正比，可通过调节励磁控制有功主力04磁阻转矩(Ld−Lq)·id·iq源于凸极效应，使凸极同步电机在失磁后仍能维持一定的异步运行能力凸极效应Chapter06转子运动方程与模型阶数选择建立机电耦合的转子动力学方程，讨论不同应用场景下的模型简化策略DYNAMICMODEL·COREEQUATION转子运动方程（摇摆方程）转子运动方程描述了原动机机械转矩与电磁转矩之间的动态平衡关系，通过惯性时间常数Tj将机械系统的转子动力学与电气系统的电磁过程耦合，是暂态稳定分析的核心方程。01基本方程：J·dωm/dt=Tm-Te-D·Δωm，J为转动惯量，D为综合阻尼系数，Δωm为转速偏差J·dω/dt02标幺值形式：Tj·dΔω/dt=Pm-Pe-D·Δω，Tj=2H为惯性时间常数，典型值火电4-8s、水电2-4s4–8s03功角方程：dδ/dt=ω0·Δω（标幺值下），将转速偏差与功角变化联系起来，构成二阶微分方程组dδ/dt04转子运动方程与电气方程联立形成机电暂态模型，是电力系统暂态稳定计算（如等面积定则）的理论基础