




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、同步电机三相短路电流和电磁转矩计算 编写 佘名寰本文是按照陈珩教授所著的同步电机运行基本理论与计算机算法一书介绍的算法和例题计算同步电机的三相短路电流。计算程序用MATLAB语言编写,计算结果与书中结果基本一致。本文可供电力系统电气技术人员和大专院校电力专业学生参考。1 计算方法1.1初始数据计算由短路前的机端电压u0, 定子绕组电流io, 和功率因数角0 求得短路前的功率角0=tan-1u0sin0+xqi0u0cos0+ri0-0从而得u0, i0 的正、交轴分量ud0=u0sin0uq0=u0cos0id0=i0sin(0+0)iq0=i0cos(0+0)短路前的空载电势是Eq0=uq0
2、+riq0+xdid0励磁电流为if0= Eq0/xaf式中 xd 为同步电机正轴同步电抗xq 同步电机交轴同步电抗 xaf 定子绕组与劢磁绕组间的互感电抗 r 定子绕组电阻1.2电流变化量的状态空间方程式同步电机突然短路时各绕组电流的变化量 id iq if iD iQ 的计算可运用以派克分量表示的状态空间方程式uduquf00=-xdxafxaD-xqxaQ-xafxfxfD-xaDxfDxD-xaQxQidiqifiDiQ+-rxq-xaQ-xd-rxafxaDrfrDrQidiqifiDiQ方程中各下标变量的含义为d-纵轴,q-横轴,f-励磁绕组,D-纵轴阻尼绕组,Q-横轴阻尼绕组,a
3、-定子绕组上式可简化为udq0=Xdq0(3)Idq0+Zdq0(3)Idq0化作电流变化量的常系数一阶微分方程组形式Idq0=-Xdq0(3)-1Zdq0(3)Idq0+Xdq0(3)-1udq0在三相短路时若励磁电压不可调,则udq0=-ud0 -uq0 0 0 0 t由于电流不能突变,t=0瞬间电流变化量的初值idq0 0= 0 0 0 0 0 t将电压变化量和电流变化量的初值代入微分方程,用数值计算的龙格-库塔法即可求出t=0+h 时刻的各电流变化量,反复计算则可求得各个时刻的idq0 ,叠加短路前绕组电流idq0 0= id0 iq0 if0 00t可得短路时电流全量idq0= id
4、 iq if iD iQ t用派克逆变换可得定子三相电流,以a相为例 ia=idcost+0-iqsint+0 0 短路t=0时转子位置角2.同步电机三相短路电流计算例题与程序电机参数r=0.005, rf=0.000656,rD=0.00151, rQ=0.00159xd=1,0, xq=0.60, xf=1.03, xD=0.95, xQ=0.70xaf=0.85, xaD=0.85, xfD=0.85, xaQ=0.45原始运行条件为额定负载 U0=1, i0=1, 0=0.5548 (单位为弧度,相对于cos=0.8) 短路时的转子位置角 0=3.1416 三相短路计算程序:CMSHO
5、RT3.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,
6、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;
7、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); end dy=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)*
8、(-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:10 a(1)=h/2; a(2)=a(1); a(3)=h; a(4)=h; for k=1:3 for i=1:5 c(i)=b(i)+a(k)*dy(i); y(i)=y(i)+a(k+1)*dy(i)/3.0; end dy=z2*(c).'+i1; end for i=1:5 y(i)=y(i)+h*dy(i)/6.0; end t=t+h; for i=1:5 b(i)=y(i)
9、; end dy=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); endplot(tt,ai); xlabel('t(rad)');ylabel('
10、ia'); title('3-phase fault of synchronous machine'); grid3程序说明:X1为状态空间方程的 Xdq0 矩阵Z1 为状态空间方程的 Zdq0 矩阵X2=Xdq0(3)-1I1=Xdq0(3)-1udq0程序 % part 1 第一部分输入初始数据,建立状态空间方程;程序 % part 2 第二部分用定步长四阶龙格-库塔法解一阶微分方程组,计算 tn+1=tn+h 时刻各绕组电流的变化量。时间步长 h用标么值,单位为弧度,h 取1/6 定步长四阶龙格-库塔法解一阶微分方程组的基本公式:y1=f1(t,y1,y2,ym)
11、, y1(t0)=y10y2=f2(t,y1,y2,ym), y2(t0)=y20.Ym=fm(t,y1,y2,ym), ym(t0)=ym0微分方程组由 tj 积分一步到 tj+1=tj+h Yi,j+1=yij+(k1i+2k2i+2k3i+k4i)h/6, i=1,2,mk1i=fi(tj,y1j,y2j,ymj), i=1,2,mk2i=fi(tj+h/2,y1j+k11h/2,ymj+k1mh/2), i=1,2,mk3i=fi(tj+h/2,y1j+k21h/2,ymj+k2mh/2), i=1,2,mk4i=fi(tj+h,y1j+k31h,ymj+k3mh), i=1,2,m程
12、序 % part 3 第三部分输出计算结果,tm 为转矩Tm =iqd+idqd=-xdid+xafif+xaDiDq=-xqiq+xaQiQ本程序计算时间周期为t=010X1/6,绘画取204.程序计算结果:T=0-10X1/6,取10个点纵轴短路电流didi = Columns 1 through 7 Columns 8 through 10 8.066892631590545 5.980344929918005 3.734356462313833A相短路电流aiai = Columns 1 through 7 Columns 8 through 10横轴短路电流qiqi = Columns 1 through 7 Columns 8 through 10励磁绕组电流fifi = Columns 1 through 7 Columns 8 through 10纵轴阻尼绕组电流zdizdi = Columns 1 through 7 Columns 8 through 10横轴阻力尼绕组电流zqizqi = Columns 1 through 7 Columns 8 through 10电磁转矩tmtm = Columns 1 through 7 C
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年床头末端装置合作协议书
- 四年级信息技术下册 装饰家园1教学实录 龙教版
- memset使用的注意事项
- 第10课《三峡》教学设计 2024-2025学年统编版语文八年级上册
- 水质提升与水生态修复计划
- 品牌社区建设的实践计划
- 幼儿园日常管理的实施细节计划
- 学校开放日活动安排计划
- 幼儿园学期教学目标探索计划
- 八年级物理上册 第二单元 第3节《我们怎样区分声音(续)》教学实录 (新版)粤教沪版
- 拉萨市2025届高三第一次联考(一模)英语试卷(含答案解析)
- 2024年中考物理压轴题专项训练:电磁继电器核心综合练(原卷版)
- 少先队活动课《民族团结一家亲-同心共筑中国梦》课件
- 2024年河北省公务员考试《行测》真题及答案解析
- 《曲唑酮治疗失眠及其相关抑郁、焦虑的专家共识(2019)》解读
- 第四单元 讲英雄故事继先辈责任(教学设计)-2022-2023学年五年级语文下册大单元教学
- 小儿常见出疹性疾病皮疹图谱和治疗课件
- 《实践论》(原文)毛泽东
- 智能制造系统解决方案“揭榜挂帅”项目申报书
- GB/T 39560.2-2024电子电气产品中某些物质的测定第2部分:拆解、拆分和机械制样
- 语文-山东省济南市2024-2025学年高中三年级九月摸底考试(暨开学考)试题和答案
评论
0/150
提交评论