弹簧质量阻尼系统模型Word文件下载.docx

弹簧质量阻尼系统模型Word文件下载.docx

《弹簧质量阻尼系统模型Word文件下载.docx》由会员分享，可在线阅读，更多相关《弹簧质量阻尼系统模型Word文件下载.docx（13页珍藏版）》请在冰豆网上搜索。

//////////

b2

ki

|pk2

图1.1机械系统图

1.2要求完成的任务

（1）推导传递函数Y（s）/X（s）,X（s）/P（s）,

（2）给定m=0.2g,b2=0.6N*s/m,k1=8N/m,k2=5N/m，以p为输入

u（t）

（3）用Matlab画出开环系统的波特图和奈奎斯特图，并用奈奎斯特判据分析系统的稳定性。

（4）求出开环系统的截止频率、相角裕度和幅值裕度。

（5）对上述任务写出完整的课程设计说明书，说明书中必须进行原理分

析，写清楚分析计算的过程及其比较分析的结果，并包含Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写。

1.3任务分析

由初始条件和要求完成的主要任务，首先对给出的机械系统进行受力分析，列出相关的微分方程，对微分方程做拉普拉斯变换，将初始条件中给定的数据代入，即可得出Y（s）/X（s），X（s）/P（s）两个传递函数。

由于本系统是一个单位负反馈系统，故求出的传递函数即为开环传函。

后在MATLAB中画出开环波特图

和奈奎斯特图，由波特图分析系统的频率特性，并根据奈奎斯特判据判断闭环系统位于右半平面的极点数，由此可以分析出系统的稳定性。

最后再计算出系统的截止频率、相角裕度和幅值裕度，并进一步分析其稳定性能。

2系统分析及传递函数求解

2.1系统受力分析

单自由度有阻尼振系的力学模型如图2-1所示，包括弹簧、质量及阻尼器

以物体的平衡位置0为原点，建立图示坐标轴X。

则物体运动微分方程为

mx=—ex—kx（2-1）

式中：

—ex为阻尼力，负号表示阻尼力方向与速度方向相反。

2-1

（2-2）

将上式写成标准形式，为

mxexkx二0

令p2=—,2n=—,则上式可简化为

mm

x2nxp2=0（2-3）

这就是有阻尼自由振动微分方程。

它的解可取x=est，其中

s是待定常数。

代入（2-1）式得（s2•2ns•p2）est=0,要使所有时间内上式都

能满足，必须s22ns•p2=0，此即微分方程的特征方程，其解为

S,2=_n±

”n2_p2（2-4）

于是微分方程（2-1）的通解为

st_nt,n2_p2tin2_p2t、

+c2e2=e（c1e+c2）

式中待定常数C1与C2决定与振动的初始条件。

振动系统的性质决定于根式

n2二p2是实数、零、还是虚数。

对应的根si与S2可以是不相等的负实根、相等的负实根或复根。

若si与S2为等根时，此时的阻尼系数值称之为临界阻尼系数，记为Cc,即Cc=2mp。

引进一个无量纲的量,称为相对阻尼系数或阻尼比。

当n>

p或>

1根式.n2-p2是实数，称为过阻尼状态，当nvp或<

1，根式

汕2二p2是虚数，称为弱阻尼状态，当n=p,即=1,称为临界阻尼状态。

现

分别讨论三种状态下的运动特性。

1•过阻尼状态

此时>

1，即n2-p2<

n，（b）式中S1及S2均为

负值，则eS1t及eS2t是两根下降的指数曲线，故（2-2）式

图2-2

所表示的是两条指数曲线之和，仍按指数衰减，不是振动。

图3-2所示为C1>

C2,C1<

0时的情况

2•临界阻尼状态

此时■二1（b）式中S1=S2=—n=—p，特征方程的根是重根，方程（2-1）

的另一解将为te—pt，故微分方程（2-1）的通解为

（2-7）

x=（C1+C2t）e—pt

式中等号右边第一项C1e—pt是一根下降的指数曲线，第二项则可应用麦克劳林级

数展开成以下形式:

c2te"

C2

pt/t

C

1/tp■p2t/2!

p3t2/3!

……pntn/n!

（2-8）

从上式看出，当时间t增长时，第二项QteP也趋近于零。

因此（c）式表示的运动也不是振动，也是一个逐渐回到平衡位置的非周期运动。

3.弱阻尼状态

此时p>

n,或<

1。

利用欧拉公式

可将（2-2）式改写为

22

_nt）（2-10）

ip2』2tC2e2卩"

）乂』"

1込；

p2匚n2tDzSin「p

或

x二Ae"

sin（p2-n2t「）（1-11）

令Pd=jp2-n2，贝U

x二Ae^sinpdt）（2-12）

式中A与「为待定常数，决定于初始条件。

设t=0时，X=X0，x=x0，则可求

得

A「x2（X。

nx0）2,-tg'

^^（2-13）

VPdX。

+nx。

将A与「代入（2-4）式，即可求得系统对初始条件的响应，由式（2-13）可知，系统振动已不再是等幅的简谐振动，而是振幅被限制在曲线_Ae之内随时间不断衰减的衰减振动。

如图3-3所示。

X

Xo

r

I

1

ti

L

J

y

图2-3

这种衰减振动的固有圆频率、固有频率和周期分别为

巳=p2一n2=P,^~2（2-14）

式中P、f、T是无阻尼自由振动的固有圆频率、固有频率和周期。

由上可见，阻尼对自由振动的影响有两个方面：

一方面是阻尼使自由振动的周期增大、频率减小，但在一般工程问题中n都比P小得多，属于小阻尼的情况。

例=n/p=0.05时,fd=0.9990f,Td=1.00125T;

而在=0.20时,fd=0.98f,Td=1.02T,所以在阻尼比较小时，阻尼对系统的固有频率和周期的影响可以略去不计，即可

以近似地认为有阻尼自由振动的频率和周期与无阻尼自由振动的频率和周期相等。

另一方面，阻尼对于系统振动振幅的影响非常显著，阻尼使振幅随着时间不

断衰减，其顺次各个振幅是：

t=t1时，A1=Ae-nt1;

t=t1+Td时，A2=Ae"

（t1Td）;

t=t1+2Td时，A3=Ae'

（t12Td），…..。

而相邻两振幅之比是个常数。

即

（2-16）

=Aj/Aj1=enTd

式中n称为减幅系数或振幅衰减率，n称为衰减系数，n越大表示阻尼越大，振幅衰减也越快。

当=0.05时，n=1.37,A2=Ai/1.37=0.73Ai,每一个周期内振幅减少27%,振幅按几何级数衰减，经过10次振动后，振幅将减小到初值的4.3%可见，衰减是非常显著的。

在工程上，通常取（2-6）式的自然对数以避免取指数的不便，即

二Ln（Aj/Aj.J二nTd（2-17）

式中S称为对数减幅或对数衰减率。

将Td=2二/.p2-n2代入，得

:

=2n/、p2-n2=2二/.J-2（2-18）

当'

<

1时，

2n（2-19）

因为任意两个相邻的振幅之比是一个常数enTd，即

Ai/A2=A?

/A3^Aj/Ai-=Aj/Ajr=e^d=e

故有

A/Aji=（A/A2）（A2/Aj）……（Aj/AjJ=ej「

因此对数减幅S也可表达为

1A

Ln—（2-20）

jA（j1）

此外，根据（3-6）式，可以用实测法来求得系统的阻尼系数。

因为

IA1Ajc1A

LnnTd=nLnLn—

AjiTdAji2mTdAj1

（2-21）

-呱丄

TdAj卡

所以只要实测得出衰减振动的周期Td及相邻两次振幅Aj和Aj+i，即可计算出系统的阻尼系数C。

根据弹簧和阻尼器的特性可得以下关系式：

Fki（t）=kix（t），Fk2（t）=k2[x（t）—y（t）],Fb2（t）=b2dy（t）/dt

设不加p（t）时，质量块处于平衡状态，此时x=0,y=0,即x（0）=0,y（0）=0,

根据受力平衡方程，在不计重力时，可得出以下方程：

（2-22）

k2[x（t）-y（t）]=b2dy（t）/dt

又根据牛顿第二定律，有方程:

md2x（t）/dt2=p（t）—Fki（t）—Fk2（t）—Fb2（t）（2-23）

2.2传递函数求解

（1）求Y（s）/X（s）:

对式（2-1）进行拉普拉斯变换，得：

k2X（s）—k2Y（s）=b2*sY（s），化简得传递函

数：

（2）求X（s）/P（s）:

一2

对式（2-2）进行拉普拉斯变换，得：

msX（s）=P（s）—k1X（s）—2k2[X（s）—Y（s）],

并将式（2-3）代入可解得传递函数：

32

X（s）/P（s）=（b2s+k2）/[mb2s+mk2s+b2（k1+2k2）s+k1k2]（2-25）

已知条件为：

给定m=0.2g,b2=0.6N.s/m,/=8N/m,k2=5N/m,设p（t）

是输入u（t）的阶跃力。

（2-26）

将所给参数代入传递函数式（2-3）和式（2-4）中，可求得具体的传递函数如下：

Y（s）/X（s）=5/（0.6s+5）

X（s）/P（s）=（0.6s+5）/（1.2*10A-4s3+10A-3s2+10.8s+40）（2-27）

2.3系统开环传递函数的求解

（1）对于丫（s）/X（s）:

由微分方程丫（s）/X（s）=5/（0.6s+5）可画出单位负反馈系统方框结构图如下:

故开环传递函数为：

G（S）=5/（0.6s+5）

（2）对于X（s）/P（s）：

由微分方程ms2X（s）=P（s）—k1X（s）—2k2[X（s）—Y（s）]及丫（s）/X（s）=k2/（b2s+k2）可

画出系统方框结构图如下：

P（s）

X（s）

故开环传递G（s）=

3.用MATLAB对系统作开环频域分析

3.1开环系统波特图

（

画波特图时采用的MATLAB语句如下：

>

num=[5];

den=（[0.6,5]）;

margin（num,den）%画系统的开环对数幅频、相频特性运

行结果如图3-1

aweJiagrem

Gm=kif,Pm«

.130deg（at0rad/sl

£

3寻呈UECiz

Q

10

1D10'

Frequency（rad/s）

CTaJs&

SEUa.

图3-1Y（s）/X（s）的开环波特图

（2）对于X（s）/P（s）:

G（s）=

num=[0.6,5];

den=（[]）;

行结果如图3-2所示：

3.2开环系统奈奎斯特图及稳定性判断

（1）对于丫（s）