东华大学：《理论力学》课程教学资源（PPT课件讲稿）第06章 多自由度体系的微振动

§63多自由度保守体系的谐振子系统 考虑理想稳定约束的俐力体系 、拉格朗日函数 在平衡位置q1=0(j=1,2,…s)附近势能的泰勒展开 aU 02U U=U.+ ∑ q;+ ∑ qqk+高阶项 了(aq 2(aq, aqk aU 平衡位置: 02U =0,记 =常数 0q) ag, dqk) →U=∑qq12,T=∑t1q/2

§6.3 多自由度保守体系的谐振子系统 U u q q /2 , T t q q / 2 q q U 0 u q U : q q q q U 2 1 q q U U U q 0 (j 1,2, s ) j,k jk j k j,k jk j k o j k 2 jk o j j,k j k o j k 2 j j o j o j      = = =            = =           +            +           = + = =    平衡位置 ， 记 常 数 高阶项 在平衡位置 附近势能的泰勒展开： 一、拉格朗日函数 考虑理想稳定约束的保守力体系

L=T-U=2(kqk-uk*)/2 二、特征频率和振动程 d OLOL ∑(t1q+u1q)=0(j=,2…s) dt(aq;aq 设解的形式为:q1=A;cos(ot+φ)(j=1,2…s) 代入运动方程得线性差程组: ∑(-0t+u1)A=0(j=1,2… 非零解的充要条件是Ak系数行列式为零,即 02t1+u1=0(j=1,2…s)∈本征方程 本征值:o,O2,g…,(正的实数;

: ( ) t u 0 j 1 2 s A ( t u )A 0 j 1 2 s q A cos( t ) j 1 2 s (t q u q ) 0 j 1 2 s q L q L dt d L T U (t q q u q q )/ 2 1 2 3 s jk jk 2 k k jk jk k 2 j j k jk k jk k j j j,k jk j k jk j k 本征值 ， ， ， 正的实数； （ ， ） 本征方程 非零解的充要条件是 系数行列式为零，即 （ ， ） 代入运动方程得线性方程组： 设解的形式为： （ ， ） （ ， ） 二、特征频率和振动方程     −  + = =  − + = = =  +  = = + = =   −           = − = −            

线性方程组:∑(02t1k+u1)Ak=0(j=1,2…s) 本征方程:0t+u=0(j=12…s) 特解:q= A! cos(ot+q) 将o1代入线性方程组,把A看作已知,可得 AW,A A,…A 2A0(1=1,2…s) 这些都是常数,共有s(s-1)个 通解:q1=∑Acos(aot+q)(j=1,2…s A共有s2个,但有s(s-1)个比值,因此A独立的只有 s2-s(s-1)=s个,再加上s个相角q,q2,…φ,共有2s 个待定常数,它们可由始条件决定

个待定常数，它们可由初始条件决定。 个，再加上 个相角 ， ， ，共有 共 有 个，但有 个比值，因 此 独立的只有 通 解 （ ， ） 这 些 都是常数，共有 个 。 ， ， （ ） 将 代入线性方程组，把 看作已知，可得 特解： 本征方程 （ ， ） 线性方程组： （ ， ） s s(s-1) s s 2s A s s(s-1) A : q A cos( t ) j 1 2 s s(s-1) A A A A A A l 1,2 s ω A q A cos( t ) : t u 0 j 1 2 s ( t u )A 0 j 1 2 s 1 2 s 2 (l) j (l) 2 j s l 1 l l (l) j j (l) j (l) 1 (l) s (l) s (l) 1 (l) 3 (l) 3 (l) 1 (l) 2 (l) 2 (l) l 1 l l (l) j (l) j jk jk 2 k jk jk k 2 − =    =  +  =  =  =  =  = =  +  −  + = = − + = =   =       

例:求振动方程见图) 解:选取广义坐标θ1,θ2 T=ml(61+02)/2 0 重=mgl(1-cos61) +mgl(1-cos 02) 2 ≈mgl(01+02)/2 mg mg U弹=k(x2-x1)2/2≈k12(2-01)2/2 k(61-26162+62)/2 U=(mgl+kl)1/2-kl6162/2 k021/2+(mgl+kl2)02/2

mg mg l l y1 y2 k θ1 θ2 kl / 2 (mgl kl ) / 2 U (mgl kl ) / 2 kl / 2 kl ( 2 )/ 2 U k(x x ) / 2 kl ( ) / 2 mgl( )/ 2 mgl(1 cos ) U mgl(1 cos ) T ml ( )/ 2 : : ( ) 2 2 2 2 1 2 1 2 2 2 1 2 2 1 2 2 2 1 2 2 2 1 2 2 2 1 2 2 2 1 2 1 2 2 2 1 2 1 2 −   + +  = +  −   =  −   +  = −   −    +  + −  = −  =  +    弹 重 解 选取广义坐标 ， 。 例 求振动方程见 图  

T=m12(02+02)/2 U=(mgl+k2)1/2-k612/2 kI2O1/2+(mgl+k2)02/2 mgl +kl-o'ml kIt 0 - KIZ mgl+kI-@ ml →(mgl+k2)-02ml2=±k2 01=g/I g ①2=g/l+2k/m Q2=g/1+2K/m

     = +  =      = +  =   + −  =  = − + −  + −  − −   + +  = +  −   =  +  g / l 2k / m g / l g / l 2k / m g / l (mgl kl ) ml kl 0 kl mgl kl ml mgl kl ml kl kl / 2 (mgl kl ) / 2 U (mgl kl ) / 2 kl / 2 T ml ( )/ 2 2 1 2 2 2 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 2 1 2 2 2 1 2 2 2 2 1 2  