基于Matlab的单闭环直流调速系统课程设计

基于Matlab的单闭环直流调速系统课程设计基于Matlab的单闭环直流调速系统2/37目录一、摘要……………………………………………………………………2二、总体方案设计………………………………………………………31、控制原理2、控制结构图三、参数计算………………………………………………………………51、静态参数设计计算2、动态参数设计计算四、稳定性分析……………………………………………………………81、基于经典自控理论得分析2、利用MATLAB辅助分析A、利用根轨迹分析B、在频域内分析奈氏曲线：bode图利用单输入单输出仿真工具箱分析用Simulink仿真五、系统校正………………………………………………………………141、系统校正的工具2、调节器的选择3、校正环节的设计4、限流装置的选择基于Matlab的单闭环直流调速系统3/37六、系统验证………………………………………………………………151、分析系统的各项指标2、单位阶跃响应3、Simulink仿真系统验证系统运行情况七、心得体会………………………………………………………………20八、参考文献………………………………………………………………20一、摘要运动控制课是后续于自动控制原理课的课程，是更加接近本专业实现应用的一门课程。直流电动机具有良好的起、制动性能，宜于在大范围内平滑调速，在许多需要调速和快速正反向的电力拖动领域中得到了广泛的应用。由于直流拖动控制系统在理论上和实践上都比较成熟，而且从控制的角度来看，它又是交流拖动控制系统的基础。所以加深直流电机控制原理理解有很重要的意义。本设计首先进行总体系统设计，然后确定各个参数，当明确了系统传函之后，再进行稳定性分析，在稳定的基础上，进行整定以达到设计要求。另外，设计过程中还要以matlab为工具，以求简明直观而方便快捷的设计过程。基于Matlab的单闭环直流调速系统4/37二、总体方案设计1、控制原理根据设计要求，所设计的系统应为单闭环直流调速系统，选定转速为反馈量，采用变电压调节方式，实现对直流电机的无极平滑调速。所以，设计如下的原理图：基于Matlab的单闭环直流调速系统5/37图1、单闭环直流调速系统原理图转速用与电动机同轴相连的测速电机产生的正比于转速的电压信号反馈到输入端，基于Matlab的单闭环直流调速系统6/37再与给定值比较，经放大环节产生控制电压，再通过电力电子变换器来调节电机回路电流，达到控制电机转速的目的。这里，电压放大环节采用集成电路运算放大器实现，主电路用晶闸管可控整流器调节对电机的电源供给。所以，更具体的原理图如下：图2、单闭环直流调速系统具体原理图2、控制结构图有了原理图之后，把各环节的静态参数用自控原理中的结构图表示，就得到了系统的稳态结构框图。基于Matlab的单闭环直流调速系统7/37图3、单闭环直流调速系统稳态结构框图同理，用各环节的输入输出特性，即各环节的传递函数，表示成结构图形式，就得到了系统的动态结构框图。由所学的相关课程知：放大环节可以看成纯比例环节，电力电子变换环节是一个时间常数很小的滞后环节，这里把它看作一阶惯性环节，而额定励磁下的直流电动机是一个二阶线性环节。所以，可以得到如下的框图：基于Matlab的单闭环直流调速系统8/37图4、单闭环直流调速系统动态结构框图三、参数计算设计完系统框图，就可以用已知的传递函数结合设计要求中给定的参数进行对系统静态和动态两套参数的计算。以便于后续步骤利用经典控制论对系统的分析。为了方便以下的计算，每个参数都采用统一的符号，这里先列出设计要求中给出的参数及大小：电动机：PN=10kwUN=220vIdN=55AnN=1000rpmRa=0.5Ω晶闸管整流装置：二次线电压E2l=230vKs=44基于Matlab的单闭环直流调速系统9/37NadNNnCeRIU主回路总电阻：R=1Ω测速发电机：PNc=23.1kwUNc=110vIdN=0.21AnNc=1900rpm系统运动部分：飞轮矩GD2=10Nm2电枢回路总电感量：要求在主回路电流为额定值10﹪时，电流仍连续生产机械：D=10s≤5﹪1、静态参数设计计算A、空载到额定负载的速降ΔnN由公式：（其中D，s已知）得：ΔnN≤5.26rpmB、系统开环放大倍数K由公式：（由公式可算出Ce=0.1925Vmin/r）得：K=53.3C、测速反馈环节放大系数a)1(NNsnsnD)1(edKCRIn基于Matlab的单闭环直流调速系统10/37设：测速发电机电动势系数=UNc/nNc=0.0579Vmin/r测速发电机输出电位器RP2分压系数a2根据经验，人为选定a2=0.2则a=Ceca2=0.01158注：1、a2正确性的验证：反馈通道的输出用于和给定比较，参照图3的标注，Un略小于Un※即可，当a2=0.2时，Un=11.58v满足要求（图1中，3为-，2为+，7要求+，也可验证）2、RP2的选择主要从功耗方面考虑，以不致过热为原则。D、运算放大器的放大系数Kp由公式（其中α即a）Kp≥20.14取Kp=21（若向小方向取，可能影响快速性，由于后加限幅电路，略大无妨）espCKKK基于Matlab的单闭环直流调速系统11/37此处的近似，使k由53.3变为55.582、动态参数设计计算在经典控制论中，动态分析基于确定系统的传函，所以要求出传函并根据已知求的传函中的未知参数，再用劳斯判据得出系统稳定性，在稳定的基础上再加校正以优化系统，使稳、准、快指标平衡在要求范围内的值上。由图4，得系统开环传函其中,Ts晶闸管装置滞后时间常数Tm机电时间常数Tl：电磁时间常数○1主电路电感值L根据要求在主回路电流为额定值10﹪时，电流仍连续。)1)(1()(m2mssTsTTsTKsWl基于Matlab的单闭环直流调速系统12/37结合抑制电流脉动的措施中关于L的讨论，得：公式：其中，整流变压器副边额定相电压（二次相电压）得：L=0.017H○2其他未知参数计算电磁时间常数机电时间常数对于三相桥式整流电路，晶闸管装置的滞后mind2693.0IULV8.1323230322lUUs017.00.1017.0RLTls075.01925.0301925.03750.110375me2mCCRGDT基于Matlab的单闭环直流调速系统13/37时间常数为Ts=0.00167s系统传函为：四、稳定性分析稳定是系统首要的条件，一切的分析只有建立在稳定的基础上才有意义。1、基于经典自控理论得分析根据系统闭环传函111)(1)1()1)(1(/1)1)(1(/)(sm2sm3smespm2msespm2msespcsKTTsKTTTsKTTTKCKKsTsTTsTCKKsTsTTsTCKKsWlllll））（（）（1s075.0s001275.01s00167.058.55sW2基于Matlab的单闭环直流调速系统14/37特征方程应用三阶系统的劳斯-赫尔维茨判据，系统稳定的充分必要条件是代入整理得：或所以：把所得参数代入就是说，k小于49.4系统才稳定。0111)(1sm2sm3smsKTTsKTTTsKTTTll0000030213210aaaaaaaa，，，，ssms)1())((TTKTTTTll0111)(smsmsmKTTTKTTKTTTlls2ssm)(TTTTTTKll4.49)(s2ssmTTTTTTKll基于Matlab的单闭环直流调速系统15/37但是，按稳态调速性能指标要求计算出的要K≥53.3它们是矛盾的。所以，当前的系统是不满足要求的。2、利用MATLAB辅助分析设计使用MATLAB版本：7.2.0.232（R2006a）【SimulinkV7.2ControlSystemToolboxV7.0】A、利用根轨迹分析系统传函可化为w=26123577/((s+20.4)*(s+38.5)*(s+600)){此处对原传函稍作近似，把原传函s2系数由0.001275近似为0.001275，即可得到本传函}○1化成LTI标准型传函程序：k=26123577;d=conv(conv([120.4],[138.5]),[1600]);------------------卷积函数基于Matlab的单闭环直流调速系统16/37s1=tf(k,d)---------------------------------------------------------化成LIT标准型函数运行结果：○2求零极点：程序：num=[2.612e007];den=[1658.93.613e0042.659e007];[z,p,k]=tf2zp(num,den)--------------------------------------基于Matlab的单闭环直流调速系统17/37----解出零极点运行结果:显然，在虚轴右半面有极点系统不稳定○3绘制根轨迹:程序：k=26123577;d=conv(conv([120.4],[138.5]),[1600]);基于Matlab的单闭环直流调速系统18/37s1=tf(k,d);rlocus(s1)-----------------------------------------------------------------绘制根轨迹函数运行结果:上图右下方极点放大图从根轨迹图上很直观地看到位于在虚轴右半基于Matlab的单闭环直流调速系统19/37面有极点，系统不稳定。B、在频域内分析○1奈氏曲线：程序：k=26123577;d=conv(conv([120.4],[138.5]),[1600]);s1=tf(k,d);figure(1);nyquist(s1)基于Matlab的单闭环直流调速系统20/37运行结果：○2bode图程序：k=26123577;d=conv(conv([120.4],[138.5]),[1600]);基于Matlab的单闭环直流调速系统21/37s1=tf(k,d);figure(1);bode(s1);margin(s1);grid运行结果：基于Matlab的单闭环直流调速系统22/37所得参数：○3利用单输入单输出仿真工具箱分析rltool和sisotool命令都能调用单输入单输出仿真程序，只是rltool仅通过根轨迹来分析，sisotool还包括bode图。菜单DesignsEditcompensator可修改增益，修改的值为当前增益的倍数（图5）。图5基于Matlab的单闭环直流调速系统23/37图6K=21时的阶跃响应：图6菜单tools下可以直接调出输出波形，可见系统振荡（图7）。图7基于Matlab的单闭环直流调速系统24/37○4用Simulink仿真图8Simulink仿真连接图输出示波器波形：很明显系统振荡（由于没加限幅，电流早已过大，电机已毁，实际中是不存在的）基于Matlab的单闭环直流调速系统25/37图9、当前系统输出曲线五、系统校正为了满足要求，还保证系统的稳定性，一般采用加调节器校正的方法来整定系统。1、系统校正的工具在设计校正装置时，主要的研究工具是Bode图，即开环对数频率特性的渐近线。它的绘制方法简便，可以确切地提供稳定性和基于Matlab的单闭环直流调速系统26/37稳定裕度的信息，而且还能大致衡量闭环系统稳态和动态的性能。因此，Bode图是自动控制系统设计和应用中普遍使用的方法。伯德图与系统性能的关系中频段以-20dB/dec的斜率穿越0dB，而且这一斜率覆盖足够的频带宽度，则系统的稳定性好；截止频率（或称剪切频率）越高，则系统的快速性越好；低频段的斜率陡、增益高，说明系统的稳态精度高；高频段衰减越快，即高频特性负分贝值越低，说明系统抗高频噪声干扰的能力越强。基于Matlab的单闭环直流调速系统2