利用Matlab仿真平台设计单闭环直流调速系统

1 大连工业大学课程设计（论文）任务书 专 业 自动化 班 级 11 级 2 班 题 目利用 Matlab 仿真平台设计单闭环直流调速系统 子 题 指导教师唐树森答疑教师唐树森 设计时间20142014 年 11 月 24 日至 20142014 年 12 月 1 日共 1 1 周 设计要求 设计（论文）的任务和基本要求，包括设计任务、查阅文献、方案设计、说明书（计算、 图纸、撰写内容及规范等） 、工作量等内容。 1、 技术要求： 电动机额定值为：，电枢 1 min/1000,55,220,10rAVKW5 . 0 a R 晶闸管装置，三相桥式可控整流电路，整流变压器 Y/Y 连接，二次线电 2 压。44,230 2 sl KVE 、主回路总电阻 3 MV 0 . 1R 测速发电机：永磁式，型，额定值为 4 110/231ZYS min/1900,21. 0 ,110, 1 . 23rAVW 生产机械要求调速范围 5 %5,10sD 系统运动部分的飞轮力矩，按时计算电 4 22 10NmGD dnomd II%10 min 枢回路总电感量：，VVEE l 8 . 1323/2303/ 22 2、设计要求： 阐述单闭环直流调速系统的基本构成和工作原理 1 分析所设计系统的静态性能指标和动态性能指标 2 根据动态性能指标设计校正装置 3 设计出系统的 Simulink 仿真模型，验证所设计系统的性能。 4 给出所设计系统的性能指标：上升时间、超调量、调节时间 5 r t% p 、最大启动电流、稳态误差。 s t maxd I ss e 指导教师签字： 系（教研室）主任签字： 2014 年 12 月 1 日 2 目目 录录 一、摘要 3 二、总体方案设计 4 1、控制原理 2、控制结构图 三、参数计算6 1、静态参数设计计算 2、动态参数设计计算 四、稳定性分析9 1、基于经典自控理论得分析 2、利用 MATLAB 辅助分析 A、利用根轨迹分析 B、在频域内分析 奈氏曲线： bode 图 利用单输入单输出仿真工具箱分析 用 Simulink 仿真 五、系统校正15 1、系统校正的工具 2、调节器的选择 3、校正环节的设计 4、限流装置的选择 六、系统验证19 1、分析系统的各项指标 2、单位阶跃响应 3、Simulink 仿真系统验证系统运行情况 七、心得体会21 八、参考文献21 3 一、摘一、摘 要要 运动控制课是后续于自动控制原理课的课程，是更加接近本专业实现应用的 一门课程。 直流电动机具有良好的起、制动性能，宜于在大范围内平滑调速，在许多需 要调速和快速正反向的电力拖动领域中得到了广泛的应用。由于直流拖动控制系 统在理论上和实践上都比较成熟，而且从控制的角度来看，它又是交流拖动控制 系统的基础。所以加深直流电机控制原理理解有很重要的意义。 本设计首先进行总体系统设计，然后确定各个参数，当明确了系统传函之后， 再进行稳定性分析，在稳定的基础上，进行整定以达到设计要求。 另外，设计过程中还要以 matlab 为工具，以求简明直观而方便快捷的设计 过程。 二、总体方案设计二、总体方案设计 1、控制原理、控制原理 4 根据设计要求，所设计的系统应为单闭环直流调速系统，选定转速为反馈量，采用 变电压调节方式，实现对直流电机的无极平滑调速。 所以，设计如下的原理图： 图 1、单闭环直流调速系统原理图 转速用与电动机同轴相连的测速电机产生的正比于转速的电压信号反馈到输入端， 再与给定值比较，经放大环节产生控制电压，再通过电力电子变换器来调节电机回路电 流，达到控制电机转速的目的。 这里，电压放大环节采用集成电路运算放大器实现，主电路用晶闸管可控整流器调 节对电机的电源供给。 所以，更具体的原理图如下： 5 图 2、单闭环直流调速系统具体原理图 2、控制结构图、控制结构图 有了原理图之后，把各环节的静态参数用自控原理中的结构图表示，就得到了系统 的稳态结构框图。 图 3、单闭环直流调速系统稳态结构框图 同理，用各环节的输入输出特性，即各环节的传递函数，表示成结构图形式，就得 到了系统的动态结构框图。 由所学的相关课程知：放大环节可以看成纯比例环节，电力电子变换环节是一个时间常 数很小的滞后环节，这里把它看作一阶惯性环节，而额定励磁下的直流电动机是一个二 6 阶线性环节。 所以，可以得到如下的框图： 图 4、单闭环直流调速系统动态结构框图 三、参数计算三、参数计算 设计完系统框图，就可以用已知的传递函数结合设计要求中给定的参数进行对系统 静态和动态两套参数的计算。以便于后续步骤利用经典控制论对系统的分析。 为了方便以下的计算，每个参数都采用统一的符号，这里先列出设计要求中给出的 参数及大小： 电动机：PN=10kw UN=220v IdN=55A nN=1000rpm Ra=0.5 晶闸管整流装置：二次线电压 E2l=230v Ks=44 主回路总电阻：R=1 测速发电机：PNc=23.1kw UNc=110v IdN=0.21A nNc=1900rpm 系统运动部分：飞轮矩 GD2=10Nm2 电枢回路总电感量：要求在主回路电流为额定值 10时，电流仍连续 生产机械：D=10 s5 1、静态参数设计计算、静态参数设计计算 7 N a dN N n Ce RIU A、空载到额定负载的速降 nN 由公式： （其中 D，s 已知） 得：nN5.26rpm B、系统开环放大倍数 K 由公式： （由公式 可算出 Ce =0.1925Vmin/r） 得：K=53.3 C、测速反馈环节放大系数 a 设：测速发电机电动势系数= UNc/ nNc=0.0579 Vmin/r 当 a2=0.2 时，U 测速发电机输出电位器 RP2分压系数 a2 根据经验，人为选定 a2=0.2 则 a=Cec a2=0.01158 注：1、a2正确性的验证： 反馈通道的输出用于和给定比较，参照图 3 的标注，Un略小于 Un即可， n=11.58v 满足要求（图 1 中，3 为-， 2 为+ ，7 要求+，也可验证） 2、RP2的选择主要从功耗方面考虑，以不致过热为原则。 D、运算放大器的放大系数 Kp 由公式 （其中 即 a） Kp 20.14 取 Kp=21 （若向小方向取，可能影响快速性，由于后加限幅电路，略大无妨） 此处的近似，使 k 由 53.3 变为 55.58 2、动态参数设计计算、动态参数设计计算 )1 ( N N sn sn D e sp C KK K )1 ( e d KC RI n 8 在经典控制论中，动态分析基于确定系统的传函，所以要求出传函并根据已知求的 传函中的未知参数，再用劳斯判据得出系统稳定性，在稳定的基础上再加校正以优化系 统，使稳、准、快指标平衡在要求范围内的值上。 由图 4，得系统开环传函 其中,Ts晶闸管装置滞后时间常数 Tm机电时间常数 Tl：电磁时间常数 主电路电感值 L 1 根据要求在主回路电流为额定值 10时，电流仍连续。 结合抑制电流脉动的措施中关于 L 的讨论，得： 公式： 其中，整流变压器副边额定相电压（二次相电压） 得：L=0.017H 其他未知参数计算 2 电磁时间常数 机电时间常数 对于三相桥式整流电路，晶闸管装置的滞后时间常数为 Ts = 0.00167 s 系统传函为：系统传函为： ) 1)(1( )( m 2 ms sTsTTsT K sW l mind 2 693 . 0 I U L V 8 .132 3 230 3 2 2 l U U s 017. 0 0 . 1 017 . 0 R L Tl s 075 . 0 1925 . 0 30 1925 . 0 375 0 . 110 375 me 2 m CC RGD T ）（ ）（ 1s075 . 0 s001275 . 0 1s00167 . 0 58.55 sW 2 9 四、稳定性分析四、稳定性分析 稳定是系统首要的条件，一切的分析只有建立在稳定的基础上才有意义。 1、基于经典自控理论得分析、基于经典自控理论得分析 根据系统闭环传函 特征方程 应用三阶系统的劳斯-赫尔维茨判据，系统稳定的充分必要条件是 代入整理得： 或 所以： 把所得参数代入 就是说，k 小于 49.4 系统才稳定。 但是，按稳态调速性能指标要求计算出的要 K 53.3 它们是矛盾的。 所以，当前的系统是不满足要求的。 01 11 )( 1 sm 2 sm 3 sm s K TT s K TTT s K TTT ll 1 11 )( 1 )1 ( ) 1)(1( / 1 ) 1)(1( / )( sm 2 sm 3 sm e sp m 2 ms esp m 2 ms esp c s K TT s K TTT s K TTT KC KK sTsTTsT CKK sTsTTsT CKK sW l l l l l 00000 30213210 aaaaaaaa， ssms )1 ()(TTKTTTT ll 0 111 )( smsmsm K TTT K TT K TTT ll s 2 ssm )( TT TTTT K l l 4 . 49 )( s 2 ssm TT TTTT K l l 10 2、利用、利用 MATLAB 辅助分析辅助分析 设计使用 MATLAB 版本：7.2.0.232 （R2006a） 【Simulink V7.2 Control System Toolbox V7.0】 A、利用根轨迹分析、利用根轨迹分析 系统传函可化为 w=26123577/(s+20.4)*(s+38.5)*(s+600) 此处对原传函稍作近似，把原传函 s2系数由 0.001275 近似为 0.001275，即可得到本传函 化成 LTI 标准型传函 1 程序： k=26123577; d=conv(conv(1 20.4,1 38.5),1 600);-卷积函数 s1=tf(k,d) -化成 LIT 标准型函数 运行结果： 求零极点： 2 程序：num=2.612e007; den=1 658.9 3.613e004 2.659e007; z,p,k=tf2zp(num,den) -解出零极点 运行结果: 显然，在虚轴右半面有极点 系统不稳定 11 绘制根轨迹: 3 程序： k=26123577; d=conv(conv(1 20.4,1 38.5),1 600); s1=tf(k,d); rlocus(s1) -