《数控工作台直线运动单元控制系统的建模与仿真分析》由会员分享，可在线阅读，更多相关《数控工作台直线运动单元控制系统的建模与仿真分析（24页珍藏版）》请在金锄头文库上搜索。
1、精选优质文档-倾情为你奉上 机电控制工程 数控工作台直线运动单元控制系统建模与仿真分析学 号 姓 名： 班 级： 指导老师： 日期: 一、题目介绍1实践题目 数控工作台单自由度直线运动单元速度开闭环控制系统建模与仿真分析2 实践目的1）、结合自动控制原理，掌握机电控制系统建模、仿真分析方法和技能； 2)、学习使用MATLAB软件Simulink工具箱构建控制系统的数学模型，绘制时域、频域曲线；3 实践任务1）建立如图（1）所示的数控工作台的直线运动单元速度控制系统数学模型，以给定电压为输入、以实际丝杠转速为输出，求出系统开环传递函数；参考给定的相关数据表1，确定关键参数，进行相应简化处理后进行
2、MATLAB/Simulink仿真分析，分析结构参数对系统性能的影响，并判断稳定性；比较matlab仿真分析结果与直线运动单元的实际运行结果，进行模型验证。2）建立如图（2）所示的数控工作台直线运动单元的速度闭环的数学模型，以给定电机转速为输入、以实际电机轴转速为输出，求出系统闭环传递函数；参考给定的相关数据表1，确定关键参数，进行相应简化处理后进行MATLAB仿真分析，分析结构参数对系统性能的影响，并判断稳定性；比较matlab仿真分析结果与直线运动单元的实际运行结果，进行模型验证。图 （1）速度开环系统图 （2）速度闭环系统表1工作台及电机参数额定电压24V反电动势常数0.0802v/rp
3、m齿轮减速比1电压放大Apm2.4电机电阻1.18欧转矩常数0.08048Nm/A电机电感1.37mH电机转子转动惯量丝杠导程5mm等效阻尼系数（参考）0.0015丝杠直径14mm速度放大增益Kw 暂取20丝杠长度360mm丝杠密度工作台质量块质量15kg4实验步骤（1）分别就图（1）与图（2）两个系统按建模步骤写出建模过程； （2）画出动态结构图； （3）图（1）以给定电压为输入、以实际丝杠转速为输出，求出系统开环传递函数； （4）图（2）以给定电机转速为输入、以实际电机轴转速为输出， 求出系统闭环传递函数； （5）采用MATLAB 对速度控制系统进行仿真分析，包括时域和频域分析， 分析结构
4、参数对系统性能的影响，并判断稳定性； （6）比较matlab 仿真与XY 工作台的实际运行效果，验证模型。二、直线运动单元的开环系统模型及仿真 1、速度开环系统建模 （1） 根据克希霍夫定律，电枢回路电压平衡方程为：（2）一般电磁转矩与电枢电流成正比，即： 其中为转矩常数 。（3）电动机轴上的转距平衡方程：电磁转矩用以驱动负载并克服摩擦力矩，其中：为电动机轴上的总转动惯量, 包括转子转动惯量及丝杠的转动惯量和工作台折算到丝杠的转动惯量。 （4）电磁感应反电动势为：上面四个方程中各个参数的含义及取值见下表2： 表2 参数含义及取值符号含义计算过程及结果说明输入电压放大倍数2.4见表1输入电压由外
5、界驱动电压决定输入量电机电枢电感1.37mH见表1电枢电流中间变量，无需求出精确值中间变量电机电阻1.18欧见表1电枢感应反电动势中间变量，无需求出精确值中间变量电机转矩中间变量，无需求出精确值中间变量转矩常数0.08048Nm/A见表1电机转子、丝杠和工作太折合到电机转轴上的转动惯量由表1参数推导得出电机转轴角速度输出量电动机和负载折合到电动机轴上的粘性摩擦系数0.0015见表1负载转矩在允许范围内由外界负载决定输入量电枢反电动势常数0.0802v/rpm见表1电枢转动惯量见表1:丝杠转动惯量工作台折合到转轴上的转动惯量L丝杠导程5mm见表1D丝杠直径14mm见表1h丝杠长度360mm见表1
6、工作台质量块质量15kg见表1丝杠密度见表1速度放大增益暂取20见表1（5）将上述方程进行拉氏变换后得到：当不计负载时，该开环系统的职能框图为： 开环系统的传递函数为： 2、使用Simulink进行仿真 当不计外加负载时，在Simulink中构建模型如下： 取输入电压为幅值24V，频率1Hz的正弦波，用Matlab对输出转速的仿真效果如下图所示：3. 在matlab程序输入窗口建立系统模型num = 0.1932; %传递函数分子den = 1.51e-7,1.33e-4,0.0082; %传递函数分母sys = tf(num,den); %建立传递函数figure(1); % 第1张图 st
7、ep(sys,g); % 阶跃响应，时域图，图线显示为绿色grid on; %网格线figure(2); % 第2张图 bode(sys,g); % 频域分析，伯德图，绿线表示grid on; %网格线figure(3); %第3张图 nyquist(sys,m); %奈氏图，红线表示grid on; %网格线figure(4); %第4张图 nichols(sys); %尼柯尔斯grid on; %网格线figure(5); %第5张图margin(sys); %求系统的稳定裕度grid on; %网格线仿真得到如下图线; 1）、阶跃响应2）、伯德图3）、奈奎斯特图4）、尼柯尔斯图5）、计算
8、稳定裕度分析：从matlab仿真得到的图线可知，该系统稳定，无超调，响应迅速，调节时间为0.05s左右，幅值和相位稳定裕度都非常大，符合系统设计的要求。4、改变系统结构参数对系统性能的影响（1）改变输入电压放大倍数 对系统性能的影响 1）、 变大：例如增大为原来的2倍 ，即取4.8，开环系统的传递函数为： M文件原始代码如下：num = 0.3863; %传递函数分子den = 1.51e-7,1.33e-4,0.0082; %传递函数分母sys = tf(num,den); %建立传递函数figure(1); % 第1张图 step(sys,g); % 阶跃响应，时域图，图线显示为绿色gri
9、d on; %网格线figure(2); % 第2张图 margin(sys); %求系统的稳定裕度grid on; %网格线此时的阶跃响应如下：稳定裕度计算图： 由上两图分析可知，当输入电压放大倍数增大时，系统的相应速度比之前快了，但是稳定裕度减小了，也就是在增加快速性的同时牺牲了稳定性。 2）、 减小：例如减小为1 开环系统的传递函数为： M文件原始代码如下：num = 0.0966; %传递函数分子den = 1.51e-7,1.33e-4,0.0082; %传递函数分母sys = tf(num,den); %建立传递函数figure(1); % 第1张图 step(sys,g); %
10、阶跃响应，时域图，图线显示为绿色grid on; %网格线figure(2); % 第2张图 margin(sys); %求系统的稳定裕度grid on; %网格线此时的阶跃响应图如下：稳定裕度计算图： 由上两图分析可知，当输入电压放大倍数减小时，系统的相应速度比之前慢了，但是稳定裕度增大了，也就是当减小电压放大倍数时，系统的快速性变差了，当变得更加稳定了。（2）阻尼系数 的影响 当取0.0005时 开环系统的传递函数为： 当取0.0025时 开环系统的传递函数为： M文件代码：num1 = 0.1932; %传递函数分子den1= 1.51e-7,1.31e-4,0.0082; %传递函数分母sys1 = tf(num1,den1); %建立传递函数num2 = 0.1932; %传递函数分子den2= 1.51e-7,1.32e-4,0.0082; %传递函数分母sys2= tf(num2,den2); %建立传递函数num3= 0.1932; %传递函数分子den3 = 1.51e-7,1.33e-4,0.0082; %传递函数分母sys3 = tf(num3,den3);
