磁悬浮小球仿真报告
- 格式:doc
- 大小:594.00 KB
- 文档页数:19
0 / 19下载文档可编辑 磁悬浮小球控制仿真报告 一.仿真要求 采用根轨迹和频域法仿真磁悬浮小球系统 二.系统建模 磁悬浮系统方程可以由下面的方程描述: 22dx(t)
mF(i,x)mg dt动力学方程
2iF(i,x)K()x 电学力学关联方程
(,)Fixmg0 边界方程 ()()1
diUtRitL
dt
电学方程
对2xiKxiF)(),(泰勒展开: )x-)(xx,(iF)i-)(ix,(iF)x,F(ix)F(i,000x000i00 )x-(xK)i-(iK)x,F(ix)F(i,0x0i00 平衡点小球电磁力和重力平衡,有 (,)Fixmg0
|,00i00iixxF(i,x)F(i,x)i;|,00x00iixx
F(i,x)F(i,x)
x
对2iF(i,x)K()x求偏导数得:
20xx00302KiKF(i,x)x 0ii00202KiKF(i,x)x 1 / 19下载文档可编辑
此系统的方程式如下: xx2Kiix2Ki)x-(xK)i-(iKdtxdm30202000x0i22
拉普拉斯变换后得: )()()(sxmx2Kisimx2Kissx30202002 由边界方程 )2020xiK(mg 代入得系统的开环传递函数: 200x(s)-1=i(s)as-b 定义系统对象的输入量为控制电压inU ,系统对象输出量为x所反映出来的输出电压为outU,则该系统控制对象的模型可写为: outssa2ina00
U(s)Kx(s)-(K/K)G(s)===
U(s)Ki(s)as-b
0000
0
iia=, b=
2gx
特征方程为:2
00as-b=0
解得系统的开环极点为:0
00
b2g
s==ax
取系统状态变量分别为1out2out x=u,x=u&
系统的状态空间表示法如下: 2 / 19下载文档可编辑
•11
ins
•
2200a
0 1 0xx
=+u2g2g?K
0-x
xxi?K
12
1xxx0 1y
代入实际参数,可以得到
in2
1
21
U124990xx0098010xx
••
..
系统的状态方程可以写为
•CXYBUAXXin
故YUin间的传递函数为
BA)(sICsUsYsG1Tin0)()()(
将以上参数值代入有
5250300.0311s77.8421sG20.)(
三.根轨迹法仿真 根据系统模型,采用根轨迹法设计一个控制器 对于传递函数5250300.0311s77.8421sG20.)(的系统,设计控制器,使得校正后系统达到以下指标:调整时间2s(%2)0ts.;最大超调量%10Mp;稳态误差=2%; 3 / 19下载文档可编辑
步骤如下: 1) 确定闭环期望极点ds的位置,由最大超调量
21pMe10(/)%
可以得到: =0.591;近似取0.6;
由cos();可以得到:=0.938(弧度)。
性能指标与根轨迹关系图 又有:2s04tns.
可以得到:3383nω=.
,于是可以得到期望的闭环极点为:
3383.(cos()jsin()) 2) 未校正系统的根轨迹不通过闭环期望极点,因此需要对
系统进行超前校正,设控制器为: cccccTs+1ks-zG(s)=k=αTs+1αs-p )(1
3) 计算超前校正装置应提供的相角,已知期望的闭环主导极点和系统原来的极点的相角和为:
-1-1d
33.8321Sin(θ)33.8321Sin(θ)G(s)=-tan-tan
33.8321Cos(θ)-31.329133.8321Cos(θ)+31.3291
=-3.803
4 / 19下载文档可编辑
因此校正装置提供的相角为:31438030661.(.). 4) 设计超前校正装置,已知:=0.938
对于最大的值的角度可由下式计算得到:)(2
1
磁悬浮根轨迹计算图 所以有:10.7712()
按最佳确定法作图规则,在上图中画出相应的直线,求出超前校正装置的零点和极点,分别为: cz2376.,cp4805.
校正后系统的开环传递函数为:
c
co
k(s+23.76)2502.96
G(s)G(s)=αs+48.05s+31.33s-31.33
5) 由幅值条件cdodG(s)G(s)=1,得0495.;0308ck. 6) 系统的校正后开环传递函数
co
(s+23.76)2502.96G(s)G(s)=0.622
s+48.05s+31.33s-31.33 5 / 19下载文档可编辑
校正后系统的根轨迹如下图所示: 校正后的根轨迹图 从图中可以看出,系统的三条根轨迹都有位于左半平面的部分,选取适当的K就可以稳定系统。系统的阶跃响应如下所示: 6 / 19下载文档可编辑
根轨迹校正后的阶跃响应 可以看出,系统有较好的稳定性,但系统存一定的稳态误差,并且误差过大,为使系统瞬态响应满足要求,可以采用直接对系统增加零点和极点的方法:
系统阶跃响应如下: 7 / 19下载文档可编辑
超前滞后控制器下的系统阶跃响应 针对磁悬浮系统,利用simulink搭建仿真结构框图如下:
simulink仿真框图 8 / 19下载文档可编辑
分别连接开关,利用示波器可以看到两种控制器下的效果图: 超前滞后控制器下的系统阶跃响应 根轨迹校正后的阶跃响应 附仿真程序如下: clear; >> t=0:0.01:1; >> s=tf('s'); G1=0.622*(s+23.76)/(s+48.05)*2502.96/((s+31.33)*(s-31.33)); 9 / 19下载文档可编辑
G2=0.991*(s+23.424)/(s+48.8648)*2502.96/((s+31.33)*(s-31.33)) ; G3=1*(s+23.424)*(s+1)/((s+48.8648)*(s+0.3))*2502.96/((s+31.33)*(s-31.33)); figure(1); rlocus(G1); G22=G2/(1+G2); G33=G3/(1+G3); figure(2); step(G22,t) grid on; figure(3); step(G33,t) grid on; 四、频率法仿真 依系统模型,采用频率法设计一个超前校正控制器, 单位负反馈系统,其开环传递函数为:
5250300.0311s77.8421sG20.)( 10 / 19下载文档可编辑
设计控制器Gc(s),使得系统的静态位置误差常数为2%,相位裕量为50°,增益裕量等于或大于10分贝。 控制器设计如下: 1) 选择控制器,由系统的Bode图
系统oG的Bode图 可以看出,给系统增加一个超前校正就可以满足设计要求,设超前校正装置为: c
cc
Ts+1ks+1/TGS=k=αTs+1αs+1/Tα
已校正系统具有开环传递函数
coc2
Τs+177.8421G(S)G(S)=kαΤs+10.0311s-30.5250
2) 根据稳态误差要求计算增益公式
pcoC2
s?0s?0
Ts+1/77.8421K=limG(s)G(s)==limkαΤs+10.0311s-30.5250