1、电气传动控制系统课程设计一、引言MATLAB 作为一个强大的数学及仿真软件，在科研与工程中被广泛使用。对于我们自动化系的学生而言，不论是专业发展、学术科研还是今后参加工作，认真学习 MATLAB 都是有很大必要的。利用 MATLAB/Simulink 验证“直流电动机转速/电流双闭环 PID 控制方案”可以熟悉 MATLAB 以及 Simulink 的使用方法，并掌握利用 MATLAB 分析控制系统性能的技巧。二、实验原理与建模1.系统建模(1) 额定励磁下的直流电动机的动态数学模型图 1 给出了额定励磁下他励直流电机的等效电路，其中电枢回路电阻 R 和电感 L 包含整流装置内阻和平波电抗器电
2、阻与电感在内，规定的正方向如图所示。图 1 直流电动机等效电路由图 1 可列出微分方程如下：（主电路，假定电流连续）0ddIURLEt（额定励磁下的感应电动势）eCn（牛顿动力学定律，忽略粘性摩擦）2375eLGDTdt（额定励磁下的电磁转矩）emI定义下列时间常数：电枢回路电磁时间常数，单位为 s；lLTR电力拖动系统机电时间常数，单位为 s；2375meGDC代入微分方程，并整理后得： 0()dddlIUERITtmdLTEIRt式中， 负载电流。/dLmITC在零初始条件下，取等式两侧得拉氏变换，得电压与电流间的传递函数1（1）0()/dlIsRUETs电流与电动势间的传递函数为（2）(
3、)ddLmsIIs0()dUs+_ ()E1/lRTs()dI_+s()dLmRTs()Ea) b)mRTs1/lRTs1eC+_+0()dUs()dIL()Ensc)图 2 额定励磁下直流电动机的动态结构图a) 式（1）的结构图 b)式（2）的结构图c)整个直流电动机的动态结构图(2) 晶闸管触发和整流装置的动态数学模型在分析系统时我们往往把它们当作一个环节来看待。这一环节的输入量是触发电路的控制电压Uct，输出量是理想空载整流电压 Ud0。把它们之间的放大系数 Ks 看成常数，晶闸管触发与整流装置可以看成是一个具有纯滞后的放大环节，其滞后作用是由晶闸管装置的失控时间引起的。下面列出不同整流
4、电路的平均失控时间：表 1 各种整流电路的平均失控时间（f=50Hz）整流电路形式 平均失控时间 Ts/ms单相半波 10单相桥式（全波） 5三相全波 3.33三相桥式，六相半波 1.67用单位阶跃函数来表示滞后，则晶闸管触发和整流装置的输入输出关系为 01()dsctsUKT按拉氏变换的位移定理，则传递函数为（3）0()sTdctUsKe由于式（3）中含有指数函数 ，它使系统成为非最小相位系统，分析和设计都比较麻烦。为sTe2了简化，先将 按台劳级数展开，则式（3）变成sTe0 23() 1!ssTd sct s sUsKKeTT考虑到 Ts 很小，忽略其高次项，则晶闸管触发和整流装置的传递
5、函数可近似成一阶惯性环节（4）0()1dsctUKsT其结构图如图 3 所示。()ctUs 0()dUssTKe ()cts 0()dUs1sKTa) b)图 3 晶闸管触发和整流装置的动态结构图a) 准确的结构图 b)近似的结构图(3)比例放大器、测速发电机和电流互感器的动态数学模型比例放大器、测速发电机和电流互感器的响应都可以认为是瞬时的，因此它们的放大系数也就是它们的传递函数，即（5）()ctpnUsK（6）()s（7）()idUI(4)双闭环控制直流电动机调速系统的动态数学模型根据以上分析，可得双闭环控制系统的动态结构图如下31/lRTsms1eCsK()ACRW()ASRsnUi c
6、tU 0dU dLI dI nnU iU_ _ _图 4 双闭环控制系统的动态结构图2. 实验系统参数系统中采用三相桥式晶闸管整流装置，基本参数如下：直流电动机：220V，13.6A，1480r/min， =0.131V/（r/min ） ，允许过载倍数 =1.5。eC晶闸管装置： 。76sK电枢回路总电阻：R=6.58。时间常数： =0.018s， =0.25s。lTm反馈系数：=0.00337V/（r/min） ，=0.4V/A。反馈滤波时间常数： =0.005s， =0.005s。oionT3.PID 调节器参数设计设计多闭环控制系统的一般原则是：从内环开始，一环一环地逐步向外扩展。在这
7、里是：先从电流环入手，首先设计好电流调节器，然后把整个电流环看作是转速调节系统中的一个环节，再设计转速调节器。双闭环控制系统的动态结构图绘于图 5，它增加了滤波环节，包括电流滤波、转速滤波和两个给定滤波环节。其中 Toi 为电流反馈滤波时间常数，T on 为转速反馈滤波时间常数41/lRTsms1eCsK()ACRW()ASRsnUiU ctU 0dU dLI dI n_ _1oT1onTsoi1oTs电流环图 5 双闭环控制系统的动态结构图(1) 电流调节器的设计对于电力拖动控制系统，电流环通常按典型型系统来设计。要把内环校正成典型型系统，显然应该采用 PI 调节器，其传递函数可以写成（8）
8、1()iACRsWK式中 Ki电流调节器的比例系数；电流调节器的超前时间常数。i为了让调节器零点对消掉控制对象的大时间常数（极点） ，选择（9）ilT一般情况下，希望超调量 %5%时，取阻尼比 =0.707， ，得：0.5IiKT， （ ） （10）12Iiisoi又因为 （11）isIKR得到 （12）0.52il liIssisiRTTRKK(2) 转速调节器的设计对于电力拖动控制系统，转速环通常希望具有良好的抗扰性能，因此我们要把转速环校正成典型型系统。要把转速环校正成典型型系统，ASR 也应该采用 PI 调节器，其传递函数为5（13）1()nASRsWsK式中 Kn电流调节器的比例系数
9、；电流调节器的超前时间常数。n转速开环增益 （14）nNemKRCT按照典型型系统的参数选择方法， （ ） （15）nnhT2nionT（16）21NnhK考虑到式（14）和（15） ，得到 ASR 的比例系数（17）(1)2emnnhCTR一般以选择 h=5 为好所以：， （18）5nnT2650NnKT经过如上设计,得到的双闭环控制系统从理论上讲有如下动态性能：电动机起动过程中电流的超调量为 4.3%，转速的超调量为 8.3%.(3) ACR 和 ASR 的理论设计及结果电流环的设计电流环的设计具体设计步骤如下：a，确定时间常数整流装置滞后时间常数 Ts按表 1，三相桥式电路的平均失控时间
10、 Ts=0.00167s。电流滤波时间常数 Toi=0.005s。电流环小时间常数 取 。i0.167.05.67sisoib，选择电流调节器结构电流调节器选择 PI 型，其传递函数为（19）1()iACRsWK6c，选择电流调节器参数ACR 超前时间常数： 。0.18silTACR 的比例系数为（20）0.186.574.90.297iiIsRKd，校验近似条件由电流环截止频率，晶闸管装置传递函数近似条件，忽略反电势对电流环影响的条件，小时间常数近似处理条件等考虑得电流调节器传递函数为（21）0.180.18().2962ACRssWs转速环的设计具体设计步骤如下：a，确定时间常数按小时间常
11、数近似处理，取 。20.134.05.1834snionTb，选择转速调节器结构由于设计要求无静差，转速调节器必须含有积分环节；又根据动态要求，应按典型型系统设计转速环。故 ASR 选用 PI 调节器，其传递函数为（22）1()nASRsWsKc，选择转速调节器参数按典型型系统最佳参数的原则，取 h=5，则 ASR 的超前时间常数为50.1834s0.917snnT转速开环增益 2-2216/s6.7.4NnhK于是，ASR 的比例系数为 (1)0.13.2519.32576804emnnhCTRd，校验近似条件从转速环截止频率，电流环传递函数简化条件，小时间常数近似处理条件等考虑得：转速调节
12、器传递函数为7（23）0.9170.917()35ASRssWsASR 输出限幅值的确定当 ASR 输出达到限幅值 U*im，转速外环呈开环状态，转速的变化对系统不再产生影响。双闭环系统变成一个电流无静差的单闭环系统。稳态时（24）*imddUII式中，最大电流 Idm 是由设计者选定的，取决于电机的过载能力和拖动系统允许的最大加速度。在这里，我们选取 Idm=20A，那么 ASR 输出限幅值为（25）*0.428VimdmUI三、实验设计与分析1.系统总体建立根据以上参数，采用 Simulink 建模如下：部分模块数据说明：饱和非线性模块：上限为 8，下限为-8；模拟结果如下：8从图 13
13、中可以清楚地看出，输出转速有很大的超调，最大可达 83.3%，这是我们所不能接受的。实践表明：应用这些工程设计方法来设计电流调节器参数，其实际电流特性与预期的比较接近。但是，由于这两种设计方法从理论上来讲都只适用于零初始条件下对线性控制系统的设计，因此，对于含有非线性环节的可控硅调速系统来说，理论和实际的矛盾比较突出。在电机起动过程的大部分时间内，转速器处于饱和限幅状态，转速环相当于开环，系统表现为恒值电流调节的单环系统。因而转速的动态响应一定有超调，只是在转速超调后，转速调节器退出饱和，才真正发挥线性调节的作用。从另一个角度上看，在转速调节器起着饱和的非线性控制作用，只有这样，才能保证内环的
14、恒值调节。所以可以看出，上述的很大的转速超调是因为我们用了零初始条件下线性控制系统的工程设计方法设计了具有非线性环节的速度环参数的结果。因此，速度调节器的设计参数与实际调试结果相差比较大，使系统对负载扰动引起的动态速降（升）缺乏有效的抑制能力，存在起动和制动过程中超调量大，突加（减）负载时，动态速降（升）大等缺点。所以，我们对 ACR 和 ASR 的参数进行整定，特别是速度控制器的参数。我们就对其作出了适当的调整，将速度控制器的传递函数改成 ，将电流调节器的传递函数改为 。0.813s0.1867s修改之后系统如图：9模拟结果如下：2.电流环跟随性能仿真实验如上文所述：电流环的作用就是保持电枢
15、电流在动态过程中不超过允许值，在突加控制作用时不希望有超调，或者超调量越小越好。这就需要我们对电流环的跟随性能加以分析。将电流环从系统中分离出来（将电枢电压对电流环影响看成是扰动） ，电流环的模型如图 17 所示。10通过如下命令可以得到电流环的 bode 图和 nyquist 图以及电流环的单位阶跃响应。 num,den=linmod(current_loop); sys=tf(num,den); margin(sys); mag,phase,w=bode(sys); gm,pm,wcg,wcp=margin(mag,phase,w); Nyquist(sys); axis equal; S
16、tep(sys)得到如下图像：11我们还可以得到以下的数据:gm =4.3078pm =48.4499wcg =12345.6682wcp =163.7923剪切频率 c=163.7923rad/s；相角相对裕度 =48.4499；-穿越频率 g=345.6682rad/s 幅值相对裕度 Lh=20lg（4.3078）=12.6851dB从图 18 与 19 种可以看出我们设计的电流环控制器是正确的，电流环是稳定的，根据剪切频率就可以看出电流的响应很快，即跟随性很好。从图 20 中可以更直接的看到这一点。在图 20 中还可以看出电流环的超调量很小（3.6%）与过渡过程时间很短（0.07s） 。
17、3.转速环抗扰性能仿真将仿真时间设置为 1.0s 到 1.5s，再对系统模型作如下修改：则可以得到该时间内 ASR、ACR 以及 I 的图线：13四、计算机仿真与结果处理由图可见，系统地工作过程可概括为如下几点：（1）ASR 从起动到稳速运行的过程中经历了两个状态，即饱和限幅输出与线性调节状态；（2）ACR 从起动到稳速运行的过程中制工作在一种状态，即线性调节状态；（3）该系统对于起动特性来说，已达到预期目的；14（4）对于系统性能指标来说，起动过程中电流的超调量为 5.3%，转速的超调量为 21.3%。这与理论最佳设计有一定差距，尤其是转速超调量略高一些。抗扰性能分析实验中我们选取 Start time=0.0，Stop time=5.0，仿真时间从 0s 到 5.0s。扰动加入的时间均为3.5s。取 Idl=12A，得到输出结果如图：由图像可见，系统的抗干扰能力较强，只是恢复时间稍长一点。
