多杆复合摆的运动研究

  • 格式:doc
  • 大小:572.50 KB
  • 文档页数:8

具有质量的多杆复合摆运动研究 宫华胜 机械茅以升班 20090873 摘要:本文首先建立了双杆复合摆的拉格朗日运动方程,通过拉普拉斯变换方程求解出其运动的解析方程,运用“科学计算与模拟平台”模拟出了在摆副比较小的条件下的运动,并在此基础上推广出了更多杆的复合摆的微幅运动方程,并对三杆复合摆进行了求解模拟。

关键词:多杆复合摆 欧拉-拉格朗日方程 拉普拉斯变换方程 前言: 在很多多杆复合摆的研究中,大多是求解到摆角的微分方程便截止,用其它软件直接模拟,虽同样可看到良好的模拟效果,但并不能让人了解到其本质,即方程本身的变化规律。本文在前人研究的基础上,对多杆复合摆进行了更深一层的研究,求解出了每个摆随时间变化的解析解,还推广到了两个以上的复合摆的运动,并给出了运动方程及求解过程,最终得到多摆运动的解析方程,以便更好地用“科学计算与模拟平台”进行模拟。

一、双杆复合摆 1、 双杆复合摆运动方程的建立 双杆摆运动示意图如图一,其中每根杆的长度、质量都相等,分别为都为l、m,第一与第二根杆与竖直方向

的夹角分别为记为1和2,重力加速度g。

_______OA杆绕O轴的转动惯量

21

1

3Jml

_______OA杆动能

22211111126k

EJml

_______AB 杆在xOy面做平面运动,其绕质心的转动惯量 图一 双复合摆

22

1

12Jml

在任意位置,_______AB 杆质心坐标

2121sinsin2xll ,2121coscos2yll

动能为

12l

lmg

mgy

x 2222222211()22kC

EmxyJ

222221212221

11[44cos()]824mlml

重力势能为 2121coscos2pEmglmgl 则系统的总动能为 12kkkEEE

222121221

1[43cos()]6ml

2221212

1(43)6ml

由于双摆为在平衡位置的微幅振动,故可以考虑21cos()1。 _______OA杆的质心坐标为

111sin2xl ,111cos2yl

重力势能为 111cos2pEmgl 则系统的重力势能为 12pppEEE

121(3coscos)2mgl

综上可得双摆系统的拉格朗日函数为 kpLEE

222121212

11[43](3coscos)62mlmgl

则可求得 21211(83)6Lml , 212

1

1()(83)6dLmldt

 , 113sin2Lmgl 由拉格朗日方程 11()0dLLdt得 22121

413

322mlmlmgl ----------------(1)

同理,对于2可求得 22122

111

232mlmlmgl ----------------(2)

2、 双杆复合摆运动方程的求解 本文采用拉普拉斯积分变换来求解1()t和2()t的具体解析式,可设1()t和2()t拉氏变换为1[()]()LtXs和2[()]()LtYs,并且 2111

2222

(0)(0)(0)(0)sXssYs



由于本文研究的是多杆复合摆的微振动下的运动,但是为了更好的模拟及展示其运动过程,在此将其运动摆角放大,虽然摆角放大,但是在模拟环境下摆仍然会按照微幅摆动的运动规律运动。所以,规定初始值为

1(0)4 , 1(0)0 ; 2(0)4 , 2(0)0

令gkl,则将初始值和变换公式带入(1)、(2)式中得, 2222

11(83)3453(2)4skXsYssXskYs





求解可得: 34223422

733()474227745()474227sksXssksksksYssksk







()Xs和()Ys分别有四个一级极点且分别相等,为

16737si ,26737si , 36737si ,46737si 则再由拉普拉斯逆变换 1()Re()e (0)knstsskftsFst



其中,Re[()]stsXse是()stXse的留数 ,ks为其一级级数, n为一级级数,k为第k个级数。

从而可以得到 11()[()]tLXs, 12()[()]tLYs

即1212733()42884istssskteksk+22273342884istssskeksk+32273342884istssskeksk +42273342884istssskeksk 21121733cos()22884skstksk+2332

3

733cos()22884skstksk

1222745()42884istssskteksk+22274542884istssskeksk+32274542884istssskeksk

+42274542884istssskeksk 21121745cos()22884skstksk+2332

3

745cos()22884skstksk

综上可得初始值为1(0)4 , 1(0)0 ; 2(0)4 , 2(0)0时的1()t和2()t

关于时间的解析式: 22

31

11322

13

22

31

21322

13

733733()cos()cos()2288422884745745()cos()cos()2288422884sksktststkskksksksktststkskksk











同样,若规定不同的初始值,会得到不同的1()t和2()t。 模拟效果如图: 1(0)4 , 1(0)0 ; 2(0)4 , 2(0)0

图二 双复合摆DTP模拟 二、多杆复合摆运动 1、多杆复合摆方程的建立 多杆系统的总动能: 12kkkknEEEE

2222112111[()]222nkkCCkk

kJmxyJ

 (2)n

总势能: 12ppppnEEEE

11212[()()]nmgyyyyyy

112211cos[cos(cos)]22nkkikimglmgl

 (2)n

同样,可得拉格朗日函数: kpLEE

且由带入拉格朗日方程中,可得n个摆关于与竖直方向的夹角的微分方程组: 0AB

其中,A,B为常矩阵,12,Tn,12,Tn





2、多杆复合摆运动方程的求解 首先规定初始值,的初始值为c,的初始值为d,其中,12,Tn







,