动态系统建模仿真 实验报告.docx

动态系统建模仿真实验报告实验二，实验四姓 名 学 号 实验二 直流电动机-负载建模及仿真实验1实验内容在运动控制系统中电机带动负载转动，电机-负载成为系统的被控对象。本实验项目要求根据电机工作原理及动力学方程，建立模型并仿真。2实验目的掌握直流电动机-负载的模型的建立方法；3实验器材（1）硬件：PC机。（2）工具软件：操作系统：Windows系列；软件工具：MATLAB及simulink。4实验原理在很多应用场合中，直流电动机的输出轴直接与负载轴相连，转动部件固定在负载轴上，即为常见的电机直接驱动负载形式。如果不考虑传动轴在转动过程中的弹性形变，即把传动轴的刚度看作无穷大，就可以在系统设计过程中，将执行电机和负载视为一个整体对象，这样被控对象的模型就可以用如图2.1所示的框图来表示。其中表示控制电压；，分别表示电机的电枢电压，电枢电感和电枢电阻；为电机的转动惯量，为负载的转动惯量，包括由电机驱动的转动体、轴承内圈、转动轴、轴套、速度测量元件、角度测量元件以及被测试件折合到电机轴上的转动惯量等；、分别表示电机和负载的粘性阻尼系数；为电机的电磁力矩系数；为电机的反电势系数；为电机-负载的转角，为电机-负载的角速度。在这一实验中，认为电机与负载的转角是相同的，并考虑了电机及负载转动中产生的粘滞阻尼力矩，所以其电压方程、力矩方程变为如下形式 (2.1)由方程组（2.1）可以得到相应的结构框图如图1所示。图1直流电动机-负载数学模型结构框图5实验要求：（1）建立从到的传递函数模型，求其频率特性，并与项目1中的电机频率特性进行对比。（2）分别取（Dm+DL）1=0.1（Dm+DL）和（Dm+DL）2=0.01（Dm+DL），编制MATLAB或simulink程序，比较阻尼系数不同时电机-负载模型的频率特性。（3）分别取JL1=0.1JL和JL=10JL，编制MATLAB或simulink程序，比较电机-负载模型的频率特性。实验所需具体参数如下表。参数备注电枢电阻Ra4.80电枢电感La21mH力矩系数km46.32N.m/A反电势系数ke55.3V/(rad/s)电机转动惯量Jm0.5Kg.m2电枢部分电机阻尼系数Dm40Nm/(rad/s)负载转动惯量JL25 Kg.m2折合到转动轴上Dm+DL270Nm/(rad/s)6实验步骤与结果（1） 建立从到的传递函数模型，由电机-负载系统框图，可以得到以下关系式：代入数据有。首先用simulink建立电机-负载模型，如图2 图2 电机负载仿真框图实验一中电机A的系统的传递函数为：根据这两个传递函数将其闭环频率特性曲线绘制出来，如下图：图3 电机A与电机负载频率特性曲线由波特图中可以看到，此系统的频率特性由其主导极点-35.3433所决定，近似于以-35.3433为极点的一阶惯性环节的频率特性，系统的带宽也大致在35rad/s左右。而对于电机A，无负载的情况下频率特性曲线为标准的二阶振荡环节的频率特性。由此可见，加入负载后的频率特性加入负载后闭环幅频曲线的谐振峰值减小，相应的时域里面超调量减小；带宽变小，调整时间变长，快速性降低。取（Dm+DL）1=0.1（Dm+DL），设为sys2，传递函数为：取（Dm+DL）2=0.01（Dm+DL），设为sys3，传递函数为：分别画出波特图，为方便比较，将sys2与sys3放在一起，我们可以看到，随着阻尼系数的减小，系统的带宽有所降低，使系统的快速性降低。但是sys2和sys3的曲线几乎重合，说明了只改变电机和负载的粘性阻尼系数并不能对系统有较大的影响。由于Dm+DL=270Nm/(rad/s)，乘以很小的系数以后，此项在参数中的比例很小，甚至可以忽略其影响。图4 阻尼系数不同时电机-负载模型的频率特性图取JL1=0.1JL，设为sys4，传递函数为：取JL1=0.1JL，设为sys5，传递函数为：画出sys4和sys5的频率特性曲线图5 负载转动惯量不同时电机-负载模型的频率特性图由此可以看出：当负载减小到一定程度时，系统会变为二阶振荡环节，从绘制的频率特性曲线上可以看出系统的带宽大大提高，系统快速性相应提高。随着负载转动惯量的增大，主导极点将快速右移，从绘制的频率特性曲线上可以看出系统的带宽降低，系统调整时间变长，时间常数增大，系统快速性相应降低，动态性能变差。实验四 位置伺服系统建模及仿真实验1实验内容本项目在建立电机-负载这一被控对象模型后，加入系统的控制器，组成一个闭环控制系统，通过仿真检验系统性能。2实验目的了解位置伺服系统的组成及工作原理；了解不同控制策略对系统性能的影响。3实验器材（1）硬件：PC机。（2）工具软件：操作系统：Windows系列；软件工具：MATLAB及simulink。4实验原理图1是一个以直流电机为驱动元件的位置伺服系统的方块图，Gc（s）为控制器，uf为与作用于转动轴上的摩擦力矩相对应的电压值。对于位置伺服控制系统，控制器的输出并不是直接驱动电机，而是经过D/A转换及功率放大后驱动电机带动负载运动。控制的目标，是使由位置传感器及测量装置给出位置反馈信号跟踪指令信号。实际的控制对象中包含D/A、功率放大器、电机、负载、位置传感器及测量装置等环节，在本实验项目中，将各环节的模型适当简化，得到广义被控对象为如下形式： (4.10)其中J为等效转动惯量，B为等效阻尼系数。图1位置伺服系统方块图5实验要求：（1）采用PID控制器对系统进行仿真，求出负载转角的响应曲线。要求考虑摩擦力矩、控制器输出饱和等非线性因素的影响。（2）采用PID+干扰观测器（DOB）控制器对系统进行仿真，求出负载转角的响应曲线，并与PID控制的响应曲线进行比较。仿真时要求考虑摩擦力矩、控制器输出饱和等非线性因素的影响。6实验步骤与结果（1）系统仿真框图如图3.6所示，图中的theta为负载转角，u为控制器输出，Kp、Kv分别为比例和微分系数。该仿真软件包含四个程序，具体见下表。名称功能备注Pd_with_friction.mdl主运行程序在simulink环境下运行Fmdl1.m仿真摩擦力矩模型用matlab编写Para_pdwithfriction.m加载仿真所需参数用matlab编写，双击图3.6中的按钮“double click here for load parameters”则可运行该程序。Postpdwithfriction.m绘制重要参数的曲线用matlab编写，双击图3.6中的按钮“double click here for plot”则可运行该程序。仿真实验的操作顺序为：第一步，设置输入指令。系统的位置指令为正弦信号，在仿真之前应设置其幅值和频率。应当注意，matlab中正弦信号的幅值是以rad为单位的，而频率是以rad/s为单位的。第二步，双击图2中的按钮“double click here for load parameters”加载对象模型及控制器参数。第三步，利用simulink的simulationparameter菜单设置仿真参数，令仿真步长为变步长，选则ode45算法。第四步，点击simulink的按钮进行仿真。第五步，双击图2中的按钮“double click here for plot”绘制曲线。也可以将变量输出至工作空间（matlab workspace）自编程序画图。在simulink中搭建的以PD控制器控制电机-负载模型的系统框图如下：图2 PD控制的simulink仿真框图其中ts=0.001s，Kp=500，Kv=10，输入电压u幅值限制：-10V+10V，死区电压：-0.5V+0.5V。输入为幅值为5弧度，频率为1rad/s的正弦信号时，电机-负载模型的输出转角跟踪正弦参考输入的效果如下图所示图3 PD控制转角的跟踪效果由上图可知，PD控制器对此正弦输入信号的跟踪效果一般，无论是幅值还是相位都存在一定的偏差。 （2）干扰观测器的仿真框图如图4所示： 图4 干扰器simulink仿真框图其主要作用是抑制摩擦力矩的影响，消除低速时的死区。取干扰器的离散域传递函数为：图5 PD+干扰器控制的simulink仿真结构图系统仿真框图如上图所示，图中的theta为负载转角，u为控制器输出，Kp、Kv分别为比例和微分系数，Q（z）和G2（z）组成干扰观测器。该仿真软件包含四个程序，具体见下表。名称功能备注Pddob.mdl主运行程序在simulink环境下运行Fmdl1.m仿真摩擦力矩模型用matlab编写Paradob.m加载仿真所需参数用matlab编写，双击图3.6中的按钮“double click here for load parameters”则可运行该程序。Postdob.m绘制重要参数的曲线用matlab编写，双击图3.6中的按钮“double click here for plot”则可运行该程序。仿真实验的操作顺序为：第一步，给定的仿真框图中没有加入控制器输出饱和环节，要求学生参照PID控制的仿真框图加入饱和环节，范围为 -10V +10V。第二步，设置输入指令。系统的位置指令为正弦信号，在仿真之前应设置其幅值和频率。第三步，双击图5中的按钮“double click here for load parameters”加载对象模型及控制器参数。第四步，利用simulink的simulationparameter菜单设置仿真参数，令仿真步长为变步长，选则ode45算法。第五步，点击simulink的按钮进行仿真。第六步，双击图5中的按钮“double click here for plot”绘制曲线。也可以将变量输出至工作空间（matlab workspace）自编程序画图。当系统的参考