多杆复合摆的运动研究
- 格式: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
动能为
12l
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
,