磁悬浮小球仿真报告

  • 格式: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