同步电机三相短路电流和转矩计算
同步发电机发生三相短路电流分析1
发电机简介
发电机通常由定子、转子、端盖及轴承等部件
组成。由轴承及端盖将发电机的定子，转子连接 组装起来，使转子能在定子中旋转，做切割磁力 线的运动，从而产生感应电势，通过接线端子引 出，接在回路中，便产生电流。
汽轮发电机是与汽轮机配套的发电机 。其转 速通常为3000转/分（频率为50赫）或3600转/分 （频率为60赫）。高速汽轮发电机为了减少因离 心力而产生的机械应力以及降低风磨耗，转子直 径一般较小，长度较大（即细长转子）。这种细 长转子使大型高速汽轮发电机的转子尺寸受到限 制。20世纪70年代以后，汽轮发电机的最大容量 达130～150万千瓦。
主磁通交链三相磁链的表达 ψb |0 |=ψ0 cos(θ0-120°)
式为：Βιβλιοθήκη ψa0=ψ0cos(θ0+ω0t)
ψc |0|=ψ0cos(θ0+120°)
ψb0=ψ0 cos(θ0+ω0t-120°)
ψc0=ψ0cos(θ0+ω0t+120°)
三相短路电流的磁链： ψai=ψa|0|-ψa0= ψ0cosθ0 -ψ0cos(ω0t)
事故
题 课题简介
同步发电机是电力系统中最重要和最复杂的元件，
由多个具有 电磁耦合关系的绕组构成。同步发电机三相
突然短路时，定子绕组中会产生很大的冲击电流，其峰 值可达额定电流的 10倍以上，从而将在电机内部产生很 大的电磁力和电磁转矩，如果设计和制造时未加充分考 虑，会使定子绕组端部受到损伤，或使转轴发生有害的 变形，还可以破坏电网的稳定和正常运行。因此，虽然 突然短路的瞬态过程时间很短，却受到设计和运行人员 的密切关注。了解短路后电流变化情况至关重要。
三相同步电机电磁计算公式
150 31.91495584
空载励磁电流Ifo= （130）励磁绕组线规a*b= 励磁绕组导线截面积qf= （132）第n层线圈平均匝长度lfn= Qm
Wm
rm 第n层线圈n= lcf=
(134)Rf(75。)=
0.6479688 0.004138
3.19748665 0.36780445 0.04597556 0.40488548
(172)直流分量时间常数Ta= 控制励磁持续短路电流倍数fko= 额定励磁持续短路电流倍数fkN= 冲击短路电流倍数fy= （176）整步功率Pr=
9有效材料
定子绕组铜重Gcu1= 励磁绕组铜重Gcu2= 定子硅钢片Gfe=
6.482295371 5.325
9.11855881 91.5988 9.16 28 0.5 2 71.98
1.016634338 3.5
43.44247714 6.946167169
5 1061.609472 42.46437888 0.736361089 0.809997198 3.039066562 0.155765707
磁极铁芯净长度lfem=
75 400 1500 50 0.8
3 135.3204388
2 36.92
26 20.4204 20.4204
24 0 0 0.96 23.04 24 22.8
铁芯计算长度li
24.2
最小气隙δ=
0.1
最大气隙δm=
0.15
定子绕组
（20）每极每相槽数q=
4
定子槽数Z1=
0.03280029 0.33303111 1.63964654 7.47053142 147.946378
三相交流系统短路电流计算
三相交流系统短路电流计算武汉华能阳光电气三相交流系统短路电流计算武汉华能阳光电气有限公司三相交流系统短路电流计算第一部分：短路电流计算1. 2. 3.图15 一台发电机馈电短路示意图4.短路电流峰值ip武汉华能阳光电气三相交流系统短路电流计算武汉华能阳光电气有限公司短路电流峰值ip的计算见第9.1.1.2条，发电机的电阻和电抗应用校正后的值，即KGRG和KGI〃d。
5. 对称开断电流Ib对称开断电流Ib是电流I〃k在经一短暂时间衰减后达到的值，用系数μ表示衰减常数，即Ib=μI〃k(46)μ与tmin和I〃k/IrG比值有关，可根据I〃k/IrG比值和选择的tmin /IrG对tmin 0.02s, =0.84+0.26e-0.26IkG/对tmin 0.52s, =0.71+0.51e-0.30IkGI对tmin 010.s, =0.62+0.72e-0.32对tmin 0.025s, =0.56+0.94e(47)上式中，μ步调相机(tmin1.6倍额定负载下的励磁电压)。
式中的I〃kG(I〃kG)和IkG应归算到同一电压下的值。
计算电动机的μ值时，式中I〃kG/I〃rM当I〃kG/IkG≤2，式(47)中μ值取μ=1。
μ值也可按图16曲线间对应的μ值，可用线性插值求取。
图16tmin＞0.1s的低压发电武汉华能阳光电气三相交流系统短路电流计算武汉华能阳光电气有限公司图16 计算开断电流Ib用系数6. 稳态短路电流Ik稳态短路电流Ik与铁芯饱和度和电网中开关状态有关，因此对Ik的计算精度比对称短路电流初始值I〃k计算精度要低。
这里给出的方法是足以估算由一台发电机或同步电机分别供电的短路电流的上、下限值。
Imax=λ式中，IrG为发电机额定电流，系数λ倒数。
a.稳态短路电流的最小值Ikminb.稳态短路电流的最大值(上限)ImaxIkmin=λλmin根据图minI(49)max根据图max IrG(48)17或18求得，xdsat()17或图18励磁。
三相短路实用计算
0.93
E8
E6 // E7
E6 X13 X 12
E7 X12 X13
1.04 1.9 0.93 0.68 0.68 1.9
1.01
13
a 1.9
11 0.03
E7 0.93
k 14 1.93
E8 1.01
3.
起始次暂态电流计算 由变压器T3方面提供的电流为
变压器T1： X6=0.105×100/31.5＝0.33
变压器T2： X7=0.105×100/20＝0.53
变压器T3： X8=0.105×100/7.5＝1.4
线路L1： X9=0.4×60×100/1152＝0.18
线路L2： X10=0.4×20×100/1152＝0.06
线路L3： X11=0.4×10×100/1152＝0.03
解：1. 取 SB 100MVA,UB Uav ，各元件电抗的标幺值计算如下：
发电机： X1=0.12×100/60＝0.2 调相机： X2=0.2×100/5＝4 负荷LD1： X3=0.35×100/30＝1.17 负荷LD2 ： X4=0.35×100/18＝1.95 负荷LD3 ： X5=0.35×100/6＝5.83
Uc
Ua
I
L2
X7
X10
0.75
0.096 0.53
C 工频周期分量，其幅值将从起始次暂态电流逐渐衰 减至稳态值
C 非周期分量和倍频周期分量，它们将逐渐衰减至零
• 短路电流计算一般指起始次暂态电流或稳态短 路电流计算；而其它任意时刻短路电流工频周
期分量有效值计算工程上采用运算曲线方法。
6.3 同步发电机突然三相短路的物理过程及短路电流分析
6.3 同步发电机突然三相短路的物理过程及短路电流分析6.3.1 同步发电机在空载情况下突然三相短路的物理过程上一节讨论了无限大电源供电电路发生三相对称短路的情况。
实际上电力系统发生短路故障时，大多数情况下作为电源的同步发电机不能看成无限大容量，其内部也存在暂态过程，因而不能保持其端电压和频率不变。
所以一般在分析和计算电力系统短路时，必须计及同步发电机的暂态过程。
由于发电机转子的惯量较大，在分析短路电流时可以近似地认为发电机转子保持同步转速，只考虑发电机的电磁暂态过程。
同步发电机稳态对称运行时，电枢磁势的大小不随时间而变化，在空间以同步速度旋转，由于它与转子没有相对运动，因而不会在转子绕组中感应出电流。
但是在发电机端突然三相短路时，定子电流在数值上将急剧变化。
由于电感回路的电流不能突变，定子绕组中必然有其它自由电流分量产生，从而引起电枢反应磁通变化。
这个变化又影响到转子，在转子绕组中感生出电流，而这个电流又进一步影响定子电流的变化。
定子和转子绕组电流的互相影响是同步电机突然短路暂态过程区别于稳态短路的显著特点，同时这种定、转子间的互相影响也使暂态过程变得相当复杂。
图6-6 凸极式同步发电机示意图图6-6为凸极同步发电机的示意图。
定子三相绕组分别用绕组，，表示，绕组的中心轴，，轴线彼此相差120o。
转子极中心线用轴表示，称为纵轴或直轴；极间轴线用轴表示，称为横轴或交轴。
转子逆时针旋转为正方向，轴超前轴90o。
励磁绕组的轴线与轴重合。
阻尼绕组用两个互相正交的短接绕组等效，轴线与轴重合的称为阻尼绕组，轴线与轴重合的称为阻尼绕组。
定子各相绕组轴线的正方向作为各绕组磁链的正方向，各相绕组中正方向电流产生的磁链的方向与绕组轴线的正方向相反，即定子绕组中正电流产生负磁通。
励磁绕组及轴阻尼绕组磁链的正方向与轴正方向一致，轴阻尼绕组磁链的正方向与轴正方向一致，转子绕组中正向电流产生的磁链与轴线的正方向相同，即在转子方面，正电流产生正磁通。
三相同步电机电磁计算公式（精）
三相同步电机电磁计算公式（精）给定区额定功率PN=75额定电压UN=400额定转速nN=1500额定频率f=50额定功率因数cos θ=0.8额定相数m=3额定电流IN=135.3204388定⼦计算区极对数P=2通风道数nK=0通风道宽度bK=0定⼦叠压系数Kfet=0.96定⼦铁芯净长度Lfet=23.04磁极铁芯总长度lm=24磁极铁芯净长度lfem=22.8线负荷A=437.3640556发热参数Aj=3679.473134(39定⼦齿距ts=1.7017 ts1=1.734425 ts2=1.930775 ts1/3=1.778058333 (40定⼦齿宽度bt1=0.834425 bt2=0.910775定⼦齿计算宽度bts=0.859875定⼦槽深hs=2.26定⼦齿计算⾼度hts‘=1.82磁路计算定⼦轭⾼度hjs=3.2定⼦轭计算⾼度hjs’=3.37定⼦轭磁路长度ljs=13.175085极弧系数αp’=0.7 (47极靴宽度bp=13.42824128磁极偏⼼距H=0.33121825极靴圆弧半径Rp=12.56878175极靴边缘⾼度hp'=0.25 (51极靴中⼼⾼度hp=2.19358252初取漏磁系数ζ‘=1.048970637磁极宽度bm=7.668255488转⼦轭内径Dir=9转⼦轭外径Djr=14磁极中⼼⾼度hm=3.70641748磁极侧⾼度hm‘=3.768404852转⼦轭⾼度hjr=2.5 (59转⼦轭计算⾼度hjr‘=4转⼦轭磁路长度ljr=3.927转⼦轭轴向长度lr=24.3磁极与轭间的残隙δ2=0.0088实际极弧系数=αp=0.693730948⽓隙⽐δm/δ=1.5最⼩⽓隙⽐极距δ/τ=0.004897064 (66基波磁场幅度系数α1=1.1151三次谐波磁场幅度系数α3=0.00646758磁场分部系数fd=0.711265508磁场波形系数fb=1.108747511直轴电枢反应磁场幅度系数Ad1=0.8532交轴电枢反应磁场幅度系数Aq1=0.33884电枢磁动势直轴折算系数Kad=0.765133172电枢磁动势交轴折算系数Kaq=0.303865124 (74定⼦卡⽒系数K δ1=1.113452078阻尼笼卡⽒系数K δ2=1.030852295卡⽒系数Kδ=1.147804629 (77空载每极总磁通θ=0.025001113斜槽系数Ksk=0.997146644⽓隙磁密最⼤值B δ=0.711292106定⼦视在磁密Bts‘=1.478524606定⼦轭磁密Bjs=1.609967011 (82⽓隙磁压降Fδ=653.1394979定⼦齿磁压降Fts=10.738定⼦轭磁压降Fjs=82.2125304⽓隙,定⼦齿,轭磁压降之和F δtj=746.0900283 (86计算漏磁⼏何尺⼨Υ1=0.409973094Υ2=0.554905146 am=3.100659528 ap=3.206420626hpm=1.54572168 (87磁极压板厚d‘=0.6磁极压板宽b’=8.2磁极计算长度lm‘=25.2极靴漏磁导Λp=6.89576E-07极⾝漏磁导Λm=8.91309E-07磁极漏磁导Λ=1.58088E-06 (93每极漏磁通θζ=0.001179482漏磁系数ζ=1.047177195磁极磁通θm=0.026180596磁极极⾝截⾯积Sm=184.6762251 (97极⾝磁密Bm=1.417648411转⼦轭磁密Bjr=1.346738462残隙处磁密B ζ2=1.422563654极⾝磁压降Fm=58.56139619转⼦轭磁压降Fjr=67.5444残隙磁压降F ζ2=100.1484812空载每极磁压降Ffo=972.3443057(104定⼦线圈尺⼨αc=0.685397076ηy=20.5257745 lF=13.25667083 lE=8.391213665 lB 30 (105线圈半匝平均长度lca=56.51334166定⼦绕组相电阻(75。
第三章 电力系统三项短路电流的使用计算
近似计算2：
假设条件：
所有发电机的电势为1，相角为 0，即 E 10 不计电阻、电纳、变压器非标准变比。 不计负荷（空载状态）或负荷用等值电抗表示。 短路电路连接到内阻抗为零的恒定电势源上
起始次暂态电流和冲击电流的 实用计算
没有给出系统信息
X S*
IB IS
有阻尼绕组 jxd
jxd 无阻尼绕组
E
E
三、起始次暂态电流和冲击电流的实用计算 1. 起始次暂态电流的计算
•起始次暂态电流：短路电流周期分量（基频分量） 的初值。
•静止元件的次暂态参数与稳态参数相同。
•发电机：用次暂态电势 E 和次暂态电抗 X d
表示。
E G 0 U G 0 jX dIG 0
三、起始次暂态电流和冲击电流的实用计算 1. 起始次暂态电流的计算
（3）短路电流使用计算步骤
较精确计算步骤
绘制电力系统等值电路图 进行潮流计算 计算发电机电势 给定短路点，对短路点进行网络简化 计算短路点电流 由短路点电流推算非短路点电流、电压。
例题
三、起始次暂态电流和冲击电流的实用计算 1. 起始次暂态电流的计算
电力系统三相短路的实用计算
三、起始次暂态电流和冲击电流的实用计算 1. 起始次暂态电流的计算
（1）同步发电机的模型
ia
Eq xd
cos(t
0 )
Ed xq
sin(t
0 )
I cos(t 0－）
ia
Eq|0| xd
当cos(xtd
0
)xq(时Exqd|0|
Exqd|0I| )cos(x1td0E)qe|0|Ttd E(qE|0x|qd|0| ExE|dx0q|d|0|
浅谈三相交流系统中短路电流计算
浅谈三相交流系统中短路电流计算摘要：在三相交流系统中，对电气元件选型，对电气元件整定计算，单相分断保护校验等均会涉及到安装点的短路电流计算，短路电流计算的正确与否，关系到所选设备的安全使用，以及该设备所在线路的甚至是系统的安全运行。
关键词：电流；电路；计算1计算短路电流的目的1.1目的1.1.1电气主接线比选。
1.1.2选择导体和电器。
1.1.3确定中性点接地方式。
1.1.4计算软导线的短路摇摆。
1.1.5确定分裂导线间隔棒的间距。
1.1.6验算接地装置的接触电压和跨步电应。
1.1.7选择继电保护装置和进行整定计算。
1.1.8为采取限制短路电流的措施提供依据（实用法计算时）。
1.1.9通常情况下，应计算两种不同幅值的短路电流:1)最大短路电流，用于选择电气设备的容量或额定值以校验电气设备的动稳定、热稳定及分断能力，整定继电保护装置;2)最小短路电流，用于选择熔断器、设定保护定值或作为校验继电保护装置灵敏度和校验感应电动机启动的依据。
2短路电流计算2.1短路电流计算的条件2.1.1短路前三相系统是正常运行情况下的接线方式，不考虑仅在切换过程中短时出现的接线方式。
.2.1.2设定短路回路各元件的磁路系统为不饱和状态，即认为各元件的感抗为一常数。
若电网电压在6kV以上时，除电缆线路应考虑电阻外，网络阻抗一般可视为纯电抗(略去电阻);若短路电路中总电阻RΣ大于总电抗XΣ的1/3，则应计人其有效电阻。
2.1.3电路电容和变压器的励磁电流略去不计。
2.1.4在短路持续时间内，短路相数不变，如三相短路保持三相短路，单相接地短路保持单相接地短路。
2.1.5电力系统中所有发电机电势相角都认为相同(大多数情况下相角很接近)。
2.1.6对于同类型的发电机，当它们对短路点的电气距离比较接近时，则假定它们的超瞬态电势的大小和变化规律相同。
因此，可以用超瞬态网络(发电机来代表)进行网络化简，并将这些发电机合并成一台等值发电机。
同步电机三相短路电流和转矩计算
同步电机三相短路电流与电磁转矩计算编写 佘名寰本文就是按照陈珩教授所著得‘同步电机运行基本理论与计算机算法’一书介绍得算法与例题计算同步电机得三相短路电流。
计算程序用MATLAB 语言编写,计算结果与书中结果基本一致。
本文可供电力系统电气技术人员与大专院校电力专业学生参考。
1． 计算方法1、1初始数据计算由短路前得机端电压u [0], 定子绕组电流i [o], 与功率因数角φ[0] 求得短路前得功率角δ0=tan −1u [0]sinφ[0]+x q i [0]u [0]cosφ[0]+ri [0]−φ[0]从而得u [0], i [0] 得正、交轴分量u d[0]=u [0]sin δ0 u q[0]=u [0]cos δ0 i d[0]=i [0]sin(δ0+φ[0]) i q[0]=i [0]cos(δ0+φ[0])短路前得空载电势就是E q[0]=u q[0]+ri q[0]+x d i d[0]励磁电流为i f[0]= E q[0]/x af式中 x d 为同步电机正轴同步电抗x q 同步电机交轴同步电抗 x af 定子绕组与劢磁绕组间得互感电抗r 定子绕组电阻1、2电流变化量得状态空间方程式同步电机突然短路时各绕组电流得变化量 ∆i d ∆i q ∆i f ∆i D ∆i Q 得计算可运用以派克分量表示得状态空间方程式[ ∆u d ∆u q ∆u f 00]=[−x d x af x aD −x qx aQ −x af x f x fD −x aDx fDx D−x aQx Q ] [ ∆i d ∆i q ∆i f ∆i D ∆i Q ] +[−r x q−x aQ −x d−rx af x aD r fr Dr Q ] [ ∆i d∆i q ∆i f ∆i D ∆i Q ]方程中各下标变量得含义为d---纵轴,q---横轴,f----励磁绕组,D---纵轴阻尼绕组,Q---横轴阻尼绕组,a---定子绕组 上式可简化为∆u dq0=X dq0(3)∆I dq0+Z dq0(3)∆I dq0化作电流变化量得常系数一阶微分方程组形式∆I dq0=−X dq0（3）−1Z dq0(3)∆I dq0+X dq0（3）−1∆u dq0在三相短路时若励磁电压不可调,则∆u dq0=[−u d [0] −u q [0] 0 0 0 ]t由于电流不能突变,t=0瞬间电流变化量得初值∆i dq0 0=[ 0 0 0 0 0 ]t将电压变化量与电流变化量得初值代入微分方程,用数值计算得龙格---库塔法即可求出 t=0+Δh 时刻得各电流变化量,反复计算则可求得各个时刻得∆i dq0 ,叠加短路前绕组电流i dq0 [0]=[ i d [0] i q [0] i f [0] 00]t可得短路时电流全量i dq0=[ i d i q i f i D i Q ]t用派克逆变换可得定子三相电流,以a 相为例i a =i d cos (t +θ0)−i q sin (t +θ0)θ0 短路t=0时转子位置角2.、同步电机三相短路电流计算例题与程序电机参数r=0、005, r f =0、000656,r D =0、00151, r Q =0、00159 x d =1,0, x q =0、60, x f =1、03, x D =0、95, x Q =0、70 x af =0、85, x aD =0、85, x fD =0、85, x aQ =0、45原始运行条件为额定负载U [0]=1, i [0]=1, φ[0]=0、5548 (单位为弧度,相对于cos φ=0、8) 短路时得转子位置角 θ0=3、1416 三相短路计算程序:CMSHORT3、M % part 1ra=0、005;rf=、000656;rzd=、00151;rzq=、00159;xd=1、0;xq=、60;xf=1、03;xzd=、95;xzq=、70;xaf=、85;xazd=、85;xfzd=、85;xazq=、450;u0=1、0;i0=1、0;phas=、5548;cita0=3、1416;p=31、4160;h=、5236; x1=[-1、0,0、0,0、85,0、85,0、0;0、0,-0、60,0、0,0、0,0、45;-0、85,0、0,1、03,0、85,0、0;-0、85,0、0,0、85,0、95,0、0;0、0,-0、45,0、0,0、0,0、70];z1=[-0、005,0、6,0、0,0、0,-0、45;-1、0,-0、005,0、85,0、85,0、0;0、0,0、0,0、000656,0、0,0、0;0、0,0、0,0、0,0、00151,0、0;0、0,0、0,0、0,0、0,0、00159];g0=(u0*sin(phas)+xq*i0)/(u0*cos(phas)+ra*i0);g0=atan(g0)-phas;ud0=u0*sin(g0);uq0=u0*cos(g0);di0=i0*sin(g0+phas);qi0=i0*cos(g0+phas);eq0=uq0+ra*qi0+xd*di0;fi0=eq0/xaf;du=[-ud0,-uq0,0、0,0、0,0、0];x2=inv(x1);z2=-x2*z1;i1=x2*du';y=[0、0,0、0,0、0,0、0,0、0];% part 2t=0、0 ;for i=1:5b(i)=y(i);enddy=z2*(y)、'+i1;tt(1)=t;di(1)=y(1)+di0;qi(1)=y(2)+qi0;fi(1)=y(3)+fi0;zdi(1)=y(4);zqi(1)=y(5);ai(1)=di(1)*cos(t+cita0)-qi(1)*sin(t+cita0);tm(1)=qi(1)*(-xd*di(1)+xaf*fi(1)+xazd*zdi(1))-di(1)*(-xq*qi(1)+ xazq*zqi(1));% while (tt(n)<p)% k=0;for n=2:10a(1)=h/2;a(2)=a(1);a(3)=h;a(4)=h;for k=1:3 for i=1:5c(i)=b(i)+a(k)*dy(i); y(i)=y(i)+a(k+1)*dy(i)/3、0; enddy=z2*(c)、'+i1; end for i=1:5y(i)=y(i)+h*dy(i)/6、0; end t=t+h; for i=1:5 b(i)=y(i); enddy=z2*(y)、'+i1; % k=k+1; % part 3 tt(n)=t; di(n)=y(1)+di0; qi(n)=y(2)+qi0; fi(n)=y(3)+fi0; zdi(n)=y(4); zqi(n)=y(5);ai(n)=di(n)*cos(t+cita0)-qi(n)*sin(t+cita0);tm(n)=qi(n)*(-xd*di(n)+xaf*fi(n)+xazd*zdi(n))-di(n)*(-xq*qi(n)+ xazq*zqi(n)); end plot(tt,ai);xlabel('t(rad)');ylabel('ia');title('3-phase fault of synchronous machine'); grid3.程序说明:X1为状态空间方程得 X dq0 矩阵 Z1 为状态空间方程得 Z dq0 矩阵 X2=X dq0（3）−1I1=X dq0（3）−1∆u dq0程序 % part 1 第一部分输入初始数据,建立状态空间方程;程序 % part 2 第二部分用定步长四阶龙格---库塔法解一阶微分方程组,计算 t n+1=t n +h 时刻各绕组电流得变化量。
三相短路电流计算
短路电流次暂态值 I ′′ = Id = 9.16 = 2.92kA X*KΣ 3.14
短路电流冲击值 ish = 2K sh I ′′ = 2.55 × 2.92 = 7.45kA
次暂态短路功率 S ′′ = 3U av I ′′ = 3 × 6.3 × 2.92 = 31.86MVA
4.4.2 有限容量电源系统的三相短路电流计算
1．有效值的计算---运算曲线法 有限容量电源系统发生三相短路后，其母线电压不再保持恒定，短路电流周期分量也随 之发生变化。如果已知短路后某一时刻发电机的电势，则短路电流周期分量相应时刻的有效 值可按下式求取
I pt =
Et 3X KΣ
（4.4.9）
式中： Et ——短路后 t 时刻的发电机的电势；
（4.4.4）
式中： I d ——计算点所在电压级的基准电流， I d =
Sd 3U d
X *KΣ ——短路回路总电抗的标幺值， X *KΣ = X KΣ
3I d Ud
于是可得：
2.短路电流冲击值
I ′′ =I ∞=
Ip
=
Id X *KΣ
（4.4.5）
由式（4.2.9）和（4.2.10）可求得短路电流冲击值和短路冲击电流有效值
X KΣ ——短路回路总电抗。
但是同步发电机突然短路时，电势随时间变化的规律是很复杂的，用上式计算比较困难。 电力部门根据国产同步发电机参数和容量配置情况，用概率统计的方法分别制定了汽轮发电 机和水轮发电机的短路电流运算曲线。利用运算曲线可以方便地查出三相短路电流周期分量
的有效值 I*pt ，因此在实际工程计算中，通常采用“运算曲线”来求解三相短路电流周期分
Id1 =
Sd = 3U av
暂态分析-三相短路电流计算
§3-4 三相短路电流周期分量起始值的 计算机算法
一.等值网络 二.用节点阻抗矩阵的计算方法 三.用节点导纳矩阵的计算方法
§3-4 三相短路电流周期分量起始值的计算机算法 在实际电力系统中进行短路电流周期分量起始值 的计算时，由于系统结构复杂，一般可采用计算 机计算。 在上机计算前需要完成两部分工作，一是根据计 算原理确定计算所用数学模型和计算方法；二是 根据所选定的数学模型和计算方法编制计算程 序。 这里主要介绍基本的数学模型和计算方法。
异步电动机的次暂态电势 E 0可由正常运行方 • 式计算而得，设正常时电动机端电压为 U 0 ，吸收 • 的电流为 I 0 ， 则：
• "
E 0 = U 0 − jx I
"
• "
•
•
0
（3－5）
由于异步电动机电阻较大，因而非周期电流分量 衰减较快。考虑到此因素，在计算短路冲击电流时 " 虽然仍应用公式 iM = K M I m ，但一般将冲击系数 KM取得较小，如容量为1000kw以上的异步电动机取 KM=1.7～1.8。 在实用计算中只对于短路点附近，显著供给短 " " 路电流的大容量电动机，才按上述方法以 E 0 、 x 作 " 为电动机的等值参数计算 I 。
实际上，x” 和异步电动机启动时的电抗相等。 启动瞬间，转子尚未转动，定子绕组和短接的鼠 笼绕组相应于一个副边短接的双绕组变压器，等 值电路与图3－1完全相同，故x” 即启动电抗，可 直接由异步电动机启动电流求得。即：
1 x = x st = I st
"
（3－4）
式中xst为电动机启动电抗标么值，Ist为启动电 流标么值，其值一般为4～7，故x” 可近似取 0.2。
三相短路分析及短路电流计算
X f ?
14
b、短路电流衰减的时间常数？
• 在超暂态过程中，只有D绕组电流存在衰减， 衰减时间常数为？ T=L/R, L=? R=? 定子、励磁绕组均短路时D绕组的时间常数 ：
TD X D Td rD
15
超暂态过程结束时刻d绕组电流值Id
• iD＝0（阻尼绕组可忽略）
d＝0，f绕组磁链不变
第二章 同步发电机的数学模型 及机端三相短路分析
（回顾）
第十六讲 三相短路分析及短路电流计算
1
问题
1、如何将短路电流计算结果与派克方程结合来分 析短路过程？
2、什么是发电机的超暂态过程、暂态过程？ 3、超暂态电抗、暂态电抗、同步电抗?大小关系？ 4、哪些绕组短路瞬间磁链不突变？ 5、短路电流计算时如何等值？ 6、为什么要计算0时刻短路电流？短路容量？
时间常数 Td0>>TD
结论：三相短路后励磁绕组中电流衰减比阻
尼绕组中电流衰减慢得多！
8
五、短路电流变化过程的假设
转子绕组直流电流（d绕组短路电流直流量） 的衰减分两个阶段：
1、超暂态过程
励磁绕组中直流电流不衰减，而D绕组中直流电流 衰减引起d绕组直流分量衰减；
2、暂态过程
D绕组中电流已衰减为零，即忽略阻尼绕组，励 磁绕组中直流电流衰减，引起d绕组直流分量衰 减到稳态。
Xad Xdl
20
b、短路电流衰减的时间常数？
无阻尼绕组，定子绕组短路时励磁绕组的时间常数为：
Xf Td
rf
X f Xad 2
Xd
fr
Xf rf
Xd Xad 2
Xf
Xd
Td0
Xd
Xd
21
第三章：电力系统三相短路实用计算
E ＿
''
＋
1
x '' d1
xL1
E ＿
''
＋
2
xd'' 2
+ xL2
U f |0|
x '' d1 xL1
xd'' 2
xL2 U f |0|
正常分量
故障分量
采用
E'' |0|
1
和忽略负荷的近似后
I
'' f
1
x '' d1
xL1
1
x'' d2
x '' L2
或者应用叠加原理，直接由故障分量求的
G
G
S LD1
L1 L2
S LD 2
f (3)
K
S LD 3
SLD1 SLD 2 SLD 3 为负荷
短路发生在 K 点
发生三相短路后的等效电路图
_
+ E1''
x '' d1
_
+ E2''
xd'' 2
xL1
零点电势等效为
xL2
U f |0|
U f |0|
上图可以等效 故障后网络=正常分量＋故障分量
SB
30 103
1650A
3U B 3 10.5
k (3) 115kV
50km
xd
xd
U S
2 N
N
U
2 B
xd 0.2
三相同步电机的转速公式
三相同步电机的转速公式
三相同步电机的转速公式是：
N = 120 * f / P
其中，N表示电机的转速（单位为转每分钟），f表示电源频率（单位为赫兹），P表示电动机的极数。
拓展：
除了转速公式，三相同步电机还有以下几个重要的公式：
1.电枢电流公式：
I_s = V_s / X_s
其中，I_s表示电枢电流（单位为安培），V_s表示电机的输入电压（单位为伏特），X_s表示电机的同步电抗（单位为欧姆）。
2.转矩公式：
T = k * I_s * I_a * sin(δ)
其中，T表示电机输出的转矩（单位为牛顿·米），k表示电机的系数，I_a表示电机的方向电流（单位为安培），δ表示电机的电击转角（单位为弧度）。
3.转差角公式：
δ = arccos((E_s - E_a) / (E_s))
其中，E_s表示电机的同步电势（单位为伏特），E_a表示电机的电枢电势（单位为伏特）。
以上这些公式可以用来计算、理解三相同步电机的性能特征和工作原理。
同步电机突然三相短路的物理分析
(
T'd
) cos(ωt + a0 ) (
Ta
E'q0 1 1 -t + + e 2 x'd xq
(
Ta
)
E'q0 1 1 -t cosa0 + e 2 x'd xq
) cos(2ωt + a0 )
Ta
i f = i f[0] +
(xd - x'd )E'q0 xad x'd
具有何种变化规律(自由电流、强制电流) ? (3) 计及电阻，确定各电流分量的衰减规律。
假设：同步电机是理想同步机——转速不变，参数用标幺值表示。
2
三、无阻尼绕组的G的f(3) 物理分析
分析方法与思路：先不计绕组电阻，磁链守恒——各绕组电流分量及初值；后计 电阻，确定各电流分量的衰减！
= ψ a ψ 0 cos(α 0 + ωt )
E′ ψ dω ψ0 q0 i dω = =cosωt = cosωt x′ x′ x′ d d d xd - x'd E'q0 xad Δi fω = i dω = - cos ω t xf x x' ad d ψ qω ψ 0 E'q0 i qω = = sinωt = sinωt xq xq xq
对于超导体闭合回路，任何扰动后都将具有“维持所环磁链 永久不变”的特性； 对于实际有电阻的线圈回路，在任何扰动瞬间都将维持其所 环磁链不突变——楞次定则！
分析方法：
(1) 电机存在多个互感耦合的绕组→电阻相对较小→首先作超导回路对待； (2) 基于磁链守恒原则，确定在突然短路暂态过程中将有哪些电流分量?
电力系统三相短路的分析计算
d q
X did X qiq
X adi f
f X adid X f i f
X X
d f
X ad X ad
Xa X f
∴
d q
( Xa X ad )id X qiq
X adi f
X aid
X ad (i f
id )
电力系统中某一处发生短路和 断相故障的情况
两个以上简单故障的组合
电力系统 短路故障
电力系统 断相故障
1.三相对称短路 2.单相接地短路 3.两相短路 4.两相接地短路
1.断一相故障 2.断两相故障
属
于
在各种短路故障中，单相
不
接地占大多数（65%），三
对
相短路的机会最少（5%）.
称
但三相短路的短路电流最
当电源距短路点的电气距离较远时，内阻抗相对于外阻抗要小得多，且 由短路而引起的电源送出功率的变化远小于电源的容量，这时认为电源的 电压幅值和频率都不发生变化。这样，可以将该电源视为恒定电势源或无 限大容量电源。
二、恒定电势源的三相短路
2.1、三相短路的暂态过程
短路前，系统中的A相电压和电流分别为
e Em sin(t ) i Im sin(t )
因此，在有限容量系统突然发生三相短路时，短路电流的 初值将大大超过稳态短路电流。
实际电机的绕组中都存在电阻，励磁绕组中的直流分量将衰减至零。 与该分量对应的定子电流中的自由分量也将逐步衰减，定子电流最终为 稳态短路电流。
三、同步发电机三相短路的暂态过程
同步发电机的暂态电势和暂态电抗
为了便于描述同步电机突然短路的暂态过程,需要确定一个短路瞬间不突 变的电势—交轴暂态电势 (通常以暂态后电势 代替)。
同步电机突然三相短路的物理分析
理论基础：回路磁链守恒原则。
对于超导体闭合回路，任何扰动后都将具有“维持所环磁链 永久不变”的特性；
对于实际有电阻的线圈回路，在任何扰动瞬间都将维持其所 环磁链不突变——楞次定则！
分析方法：
(1) 电机存在多个互感耦合的绕组→电阻相对较小→首先作超导回路对待； (2) 基于磁链守恒原则，确定在突然短路暂态过程中将有哪些电流分量?
E'q0 x'd
-
Eq[0] xd
注意： ψq0=0 → 定子基频电流q轴分量=0
9
5-4 二、不计衰减时的短路电流
2、定子iap、i2ω、转子Δifω
Δifω
令：ψf=0、 ψa= ψa0(b、c类似)
xσf
t=0 磁平衡等值电路
idω xσa
xad ψdω ↑
iqω
xσa
xaq
ψ
↑
qω
a0 0 cos0
二、不计衰减时的短路电流
1、定子基频交流和 f-f 直流电流
设：空载短路——iω、if[0]、Δifa共同作用，
保持各相定子绕组磁链为0、f-f 初值ψf0
if[0] + Δifa
xσf ↑Ψf0
id xσa
xad
id
i fa
xad x f
xd
xad xf
id
f 0 Eq0
xd xd xd
id
↑ψd
xσf xad x'd
E'q x'd
= =
xad xf
xσa
ψf
= σf
xad xσf
+ xσf xad xσf + xad
电力系统三相短路实用算法
3 电力系统三相短路的实用计算①起始次暂态电流I"(短路电流基频交流分量的初始值)、冲击电流（短路电流最大瞬时值）、短路电流最大有效值、短路容量；（用于效验断路器开断电流、继电保护整定、电气设备动稳定效验）；②采用运算曲线法近似计算电网三相短路暂态过程中，任一时刻短路电流（交流分量的有效值）3.1交流电流初始值的计算一、计算近似假设（各个元件次暂态参数的获取）1）发电机①电抗：用x d"；②电动势：用E"(近似认为短路前后瞬间保持不变)相量表示：E0"=U0+jI0x d"标量表示：E0"≈U0+jI0x d"sinφ|0|其中：I|0|=P|0|−jQ|0|U0③近似计算中可取E"=1.05~1.08④不计负荷影响时（短路前空载），E"=1，且同相位。
⑤当电源远离短路点，可将发电机看作恒定电压源，取其额定电压U N。
2）线路、变压器① 并联支路：忽略线路对地电容、变压器励磁回路； ② 高压输电线路：仅考虑线路电抗，忽略电阻； ③低压输电线路或电缆：近似用阻抗模值z = 2+x 2 ④变压器变比：不考虑实际变比，用平均电压比。
3) 一般负荷①不考虑负荷（即短路前空载）：基于负荷电流远小于短路电流。
②考虑负荷：恒定阻抗负荷：z i =U i|0|2P i|0|−jQ i|0|综合负荷：E "=0.8，x "=0.35远离短路点的负荷：略去不计或x "=0.354) 短路点附近的大型异步（同步）电动机负荷：①正常运行时，异步电动机的转差率很小（2%~5%），可作同步机看待。
则根据短路瞬间磁链守恒原理，可用与转子绕组总磁链成正比的E "、x "(为启动电抗)表示。
如短路瞬间的机端电压小于E "，则考虑到送短路电流，当作发电机看待。
E "、x "的确定：x "=1I st =14～7=0.14～0.25，近似x "≅0.2E 0 "≈U 0 −jI 0 x "sin φ|0|，近似E 0 "≅0.9（I "≅0.45）②如短路瞬间的机端电压大于E "，当作综合负荷看待。
三相交流系统短路电流计算(GBT 15544—1995 )
中华人民共和国国家标准三相交流系统短路电流计算GB/T15544—1995Short-circuit current calculation in three-phase a.c.systems 国家技术监督局1995-04-06批准1996-01-01实施本标准等效采用IEC909(1988)《三相交流系统短路电流计算》(以下简称《909标准》)。
第一篇概述1主题内容与适用范围1.1主题内容本标准规定了用等效电压源法计算三相交流系统短路电流，并提出了计算中采用的校正系数的求取方法及推荐值。
1.2适用范围本标准适用于标称电压380V～220kV，频率50Hz的三相交流系统的短路电流计算。
本标准不适用于受控条件(短路试验站)下人为短路和飞机、船舶用电气设备的短路计算。
本标准主要作为进出口设备及对外工程投标使用，在国内工程计算中逐步推广采用。
2引用标准GB156—93额定电压GB2900.1—92电工术语基本术语GB2900.25—94电工术语旋转电机3术语3.1短路short-circuit通过一个比较低的电阻或阻抗，偶然地或有意地对正常电路中不同电压下的两个或几个点之间的连接。
3.2短路电流short-circuit current在电路中，由于故障或不正确连接造成短路而产生的过电流。
注：需区别流过短路点和电网支路中的短路电流。
3.3预期(可达到的)短路电流prospective(available)short-circuit current电源不变，将短路点用阻抗可忽略的理想连接代替时，流过短路点的电流。
注：假设三相短路电流是由于三相同时短路而产生的。
由于三相不在同一瞬间短路，在短路电流中可能出现较大的非周期分量的研究不属于本标准范围。
3.4对称短路电流symmetrical short-circuit current不计非周期分量时的预期(可达到的)短路电流对称交流分量的有效值。
3.5对称短路电流初始值initial symmetrical short-circuit current系统非故障元件的阻抗保持为短路前瞬间值时的预期(可达到的)短路电流的对称交流分量有效值(见图1和图12)。
同步发电机突然三相短路的物理过程及短路电流分析
磁链匝链定子三相绕组的磁链随着"λ'的变化而变化，因此6。
3同步发电机突然三相短路的物理过程及短路电流分析6.3.1同步发电机在空载情况下突然三相短路的物理过程上一节讨论了无限大电源供电电路发生三相对称短路的情况。
实际上电力系统发生短路故障时，大多数 情况下作为电源的同步发电机不能看成无限大容量,其内部也存在暂态过程，因而不能保持其端电压和频率 不变.所以一般在分析和计算电力系统短路时，必须计及同步发电机的暂态过程。
由于发电机转子的惯量较 大，在分析短路电流时可以近似地认为发电机转子保持同步转速，只考虑发电机的电磁暂态过程。
同步发电机稳态对称运行时，电枢磁势的大小不随时间而变化，在空间以同步速度旋转,由于它与转子 没有相对运动，因而不会在转子绕组中感应岀电流。
但是在发电机端突然三相短路时，定子电流在数值上将 急剧变化。
由于电感回路的电流不能突变,定子绕组中必然有其它自由电流分量产生，从而引起电枢反应磁 通变化。
这个变化又影响到转子，在转子绕组中感生岀电流，而这个电流又进一步影响定子电流的变化。
定 子和转子绕组电流的互相影响是同步电机突然短路暂态过程区别于稳态短路的显著特点，同时这种定、转子 间的互相影响也使暂态过程变得相当复杂。
图6—6凸极式同步发电机示意图图6-6为凸极同步发电机的示意图。
定子三相绕组分别用绕组二 上,一 一，L —表示,绕组的中心轴'，-，-轴线彼此相差120°。
转子极中心线用「轴表示，称为纵轴或直轴;极间轴线用 「轴表示， 称为横轴或交轴。
转子逆时针旋转为正方向,「轴超前「轴90°。
励磁绕组㈡的轴线与•轴重合。
阻尼绕组 用两个互相正交的短接绕组等效,轴线与 「轴重合的称为 Wr 阻尼绕组，轴线与£轴重合的称为少；阻尼绕组。
定子各相绕组轴线的正方向作为各绕组磁链的正方向,各相绕组中正方向电流产生的磁链的方向与绕组 轴线的正方向相反，即定子绕组中正电流产生负磁通。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
同步电机三相短路电流和电磁转矩计算编写佘名寰本文是按照珩教授所著的‘同步电机运行基本理论与计算机算法’一书介绍的算法和例题计算同步电机的三相短路电流。
计算程序用MATLAB语言编写，计算结果与书中结果基本一致。
本文可供电力系统电气技术人员和大专院校电力专业学生参考。
1．计算方法1.1初始数据计算由短路前的机端电压u[0], 定子绕组电流i[o], 和功率因数角φ[0] 求得短路前的功率角−t[t]δ0=ttt−t t[t]tttt[t]+t t t[t]t[t]tttt[t]+tt[t]从而得u[0], i[0]的正、交轴分量u d[0]=u[0]sinδ0u q[0]=u[0]cosδ0i d[0]=i[0]sin(δ0+φ[0])i q[0]=i[0]cos(δ0+φ[0])短路前的空载电势是E q[0]=u q[0]+ri q[0]+x d i d[0]励磁电流为i f[0]= E q[0]/x af式中x d为同步电机正轴同步电抗x q同步电机交轴同步电抗x af定子绕组与劢磁绕组间的互感电抗r 定子绕组电阻1.2电流变化量的状态空间方程式同步电机突然短路时各绕组电流的变化量 ∆i d ∆i q ∆i f ∆i D ∆i Q的计算可运用以派克分量表示的状态空间方程式[∆t t∆t t ∆t t t t ]=[−t t t ttt tt −t tt tt −t tt t t t tt −t tt t ttt t−t ttt t ] [∆t t ∆t t ∆t t ∆t t ∆t t]+[−tt t −t tt −t t −tt tt t tt t tt tt t ] [∆t t∆t t ∆t t ∆t t ∆t t ]方程中各下标变量的含义为d---纵轴，q---横轴，f----励磁绕组，D---纵轴阻尼绕组，Q---横轴阻尼绕组，a---定子绕组上式可简化为∆u dq0=X dq0(3)∆I dq0+Z dq0(3)∆I dq0化作电流变化量的常系数一阶微分方程组形式∆I dq0=−X dq0（3）−1Z dq0(3)∆I dq0+X dq0（3）−1∆u dq0在三相短路时若励磁电压不可调，则∆u dq0=[−u d [0] −u q [0] 0 0 0 ]t由于电流不能突变，t=0瞬间电流变化量的初值∆i dq0 0=[ 0 0 0 0 0 ]t将电压变化量和电流变化量的初值代入微分方程，用数值计算的龙格---库塔法即可求出 t=0+Δh 时刻的各电流变化量，反复计算则可求得各个时刻的∆i dq0 ，叠加短路前绕组电流i dq0 [0]=[ i d [0] i q [0] i f [0] 00]t可得短路时电流全量i dq0=[ i d i q i f i D i Q ]t用派克逆变换可得定子三相电流，以a 相为例i a =i d cos (t +θ0)−i q sin (t +θ0)θ0 短路t=0时转子位置角2．.同步电机三相短路电流计算例题与程序 电机参数r=0.005, r f =0.000656，r D =0.00151, r Q =0.00159x d =1,0, x q =0.60, x f =1.03, x D =0.95, x Q =0.70 x af =0.85, x aD =0.85, x fD =0.85, x aQ =0.45原始运行条件为额定负载U [0]=1, i [0]=1, φ[0]=0.5548 (单位为弧度，相对于cos φ=0.8) 短路时的转子位置角 θ0=3.1416三相短路计算程序：CMSHORT3.M % part 1ra=0.005;rf=.000656;rzd=.00151;rzq=.00159;xd=1.0;xq=.60;xf=1.03;xzd=.95;xzq=.70;xaf=.85;xazd=.85;xfzd=.85;xazq=.450; u0=1.0;i0=1.0;phas=.5548;cita0=3.1416;p=31.4160;h=.5236; x1=[-1.0,0.0,0.85,0.85,0.0;0.0,-0.60,0.0,0.0,0.45; -0.85,0.0,1.03,0.85,0.0; -0.85,0.0,0.85,0.95,0.0;0.0,-0.45,0.0,0.0,0.70]; z1=[-0.005,0.6,0.0,0.0,-0.45;-1.0,-0.005,0.85,0.85,0.0; 0.0,0.0,0.000656,0.0,0.0; 0.0,0.0,0.0,0.00151,0.0;0.0,0.0,0.0,0.0,0.00159];g0=(u0*sin(phas)+xq*i0)/(u0*cos(phas)+ra*i0); g0=atan(g0)-phas; ud0=u0*sin(g0); uq0=u0*cos(g0); di0=i0*sin(g0+phas); qi0=i0*cos(g0+phas); eq0=uq0+ra*qi0+xd*di0; fi0=eq0/xaf;du=[-ud0,-uq0,0.0,0.0,0.0]; x2=inv(x1); z2=-x2*z1; i1=x2*du';y=[0.0,0.0,0.0,0.0,0.0]; % part 2 t=0.0 ; for i=1:5 b(i)=y(i); enddy=z2*(y).'+i1;tt(1)=t;di(1)=y(1)+di0;qi(1)=y(2)+qi0;fi(1)=y(3)+fi0;zdi(1)=y(4);zqi(1)=y(5);ai(1)=di(1)*cos(t+cita0)-qi(1)*sin(t+cita0);tm(1)=qi(1)*(-xd*di(1)+xaf*fi(1)+xazd*zdi(1))-di(1)*(-xq*qi(1)+xazq*zqi(1));% while (tt(n)<p)% k=0;for n=2:10a(1)=h/2;a(2)=a(1);a(3)=h;a(4)=h;for k=1:3for i=1:5c(i)=b(i)+a(k)*dy(i);y(i)=y(i)+a(k+1)*dy(i)/3.0;enddy=z2*(c).'+i1;endfor i=1:5y(i)=y(i)+h*dy(i)/6.0;endt=t+h;for i=1:5b(i)=y(i);enddy=z2*(y).'+i1;% k=k+1;% part 3tt(n)=t;di(n)=y(1)+di0;qi(n)=y(2)+qi0;fi(n)=y(3)+fi0;zdi(n)=y(4);zqi(n)=y(5);ai(n)=di(n)*cos(t+cita0)-qi(n)*sin(t+cita0);tm(n)=qi(n)*(-xd*di(n)+xaf*fi(n)+xazd*zdi(n))-di(n)*(-xq*qi(n)+ xazq*zqi(n)); endplot(tt,ai);xlabel('t(rad)');ylabel('ia');title('3-phase fault of synchronous machine'); grid3．程序说明：X1为状态空间方程的 X dq0 矩阵 Z1 为状态空间方程的 Z dq0 矩阵 X2=X dq0（3）−1I1=X dq0（3）−1∆u dq0程序 % part 1 第一部分输入初始数据，建立状态空间方程； 程序 % part 2 第二部分用定步长四阶龙格---库塔法解一阶微分方程组，计算 t n+1=t n +h 时刻各绕组电流的变化量。
时间步长 h 用标么值，单位为弧度,h 取1/6 π 定步长四阶龙格---库塔法解一阶微分方程组的基本公式：y 1’=f 1(t,y 1,y 2,…,y m ), y 1(t 0)=y 10 y 2’=f 2(t,y 1,y 2,…,y m ), y 2(t 0)=y 20 …….Y m ’=f m (t,y 1,y 2,…,y m ), y m (t 0)=y m0微分方程组由 t j 积分一步到 t j+1=t j +hY i,j+1=y ij +(k 1i +2k 2i +2k 3i +k 4i )h/6, i=1,2,…,m k 1i =f i (t j ,y 1j ,y 2j ,…,y mj ), i=1,2,…,mk 2i =f i (t j +h/2,y 1j +k 11h/2,…,y mj +k 1m h/2), i=1,2,…,m k 3i =f i (t j +h/2,y 1j +k 21h/2,…,y mj +k 2m h/2), i=1,2,…,m k 4i =f i (t j +h,y 1j +k 31h,…,y mj +k 3m h), i=1,2,…,m程序 % part 3 第三部分输出计算结果，tm 为转矩Tm =i q ψd +i d ψq ψd =-x d i d +x af i f +x aD i D ψq =-x q i q +x aQ i Q本程序计算时间周期为t=0—10X1/6π,绘画取20π4.程序计算结果：T=0--10X1/6π,取10个点纵轴短路电流didi =Columns 1 through 70.6714 2.3222 4.2319 6.0811 8.8684 9.9698 9.7469Columns 8 through 108.0545 5.8005 3.3833A相短路电流aiai =Columns 1 through 7-0.4470 -0.0735 0.9780 2.3802 5.4181 8.2065 9.0178Columns 8 through 108.4739 5.0433 3.4802横轴短路电流qiqi =Columns 1 through 70.7368 1.0078 2.6660 2.4421 1.2145 0.3252 -1.0709Columns 8 through 10-2.6902 -3.5734 -3.5201励磁绕组电流fifi =Columns 1 through 72.2884 2.69513.09844.0924 4.35965.3790 4.7447Columns 8 through 104.1245 3.3933 3.3557纵轴阻尼绕组电流zdizdi =Columns 1 through 70 0.1670 2.7617 3.8854 4.8734 5.9673 5.8786 Columns 8 through 104.0548 3.8983 1.6384横轴阻力尼绕组电流zqizqi =Columns 1 through 70 0.8486 1.6528 1.3091 0.7257 -0.7101 -1.0671 Columns 8 through 10-2.2395 -2.6311 -2.6190电磁转矩tmtm =Columns 1 through 70.1161 2.5244 4.4641 5.5277 4.3285 1.9010 -1.7246Columns 8 through 10-3.4925 -4.8964 -4.4002三相短路电流和转矩波形图A相短路电流图（故障前额定负载，θ0=3.1416）转矩010203040506070t(rad)ia3-phase fault of synchronous machineA 相短路电流（短路前空载，转子角 θ0=π）转矩010203040506070-6-4-2246t(rad)t m10203040506070t(rad)i a3-phase fault of synchronous machineA 相短路电流（短路前空载，转子角 θ0=π/2）转矩10203040506070-5-4-3-2-1012345t(rad)t m10203040506070-5-4-3-2-1012345t(rad)i a3-phase fault of synchronous machine参考文献：1．同步电机运行基本理论与计算机算法 箸者 珩 水利电力 2．FORTRAN 常用算法程序集 徐士良 编 清华大学3．MATLAB 电力系统设计与分析 吴天明等编箸 国防工业10203040506070-5-4-3-2-1012345t(rad)t m。