结构力学15动力学(1)ppt课件

1、二、单自在度体系自在振动微分方程的解二、单自在度体系自在振动微分方程的解)(mk=w02yy=+w& &)(0akyym=+LL& &) sin()(aw+=tatysincos)(00www+=tvtyty)0(020=yCyycossin)(21ww+=tCtCty)0(010w=vCvy &y(t)ty0y0y(t )tv0/v0/TtaaT/其中是沿质点振动方向的构造柔度系数，它表示在质点上沿振动方向加单位荷载使质点沿振动方向所产生的位移。 k使质点沿振动方向发生单位位移时，须在质点上沿振动方向施加的力。 st=W在质点上沿振动方向施加数值为W的

2、荷载时质点沿振动方向所产生的位移。计算时可根据体系的详细情况，视、 k、 st 三参数中哪一个最便于计算来选用。自振周期计算公式：自振周期计算公式：圆频率计算公式：圆频率计算公式：一些重要性质：一些重要性质：1 1自振周期与且只与构造的质量和构造的刚度有关，与外界自振周期与且只与构造的质量和构造的刚度有关，与外界的干扰要素无关。干扰力只影响振幅的干扰要素无关。干扰力只影响振幅a a。2 2自振周期与质量的平方根成正比，质量越大，周期越大自振周期与质量的平方根成正比，质量越大，周期越大频率越小；自振周期与刚度的平方根成反比，刚度越大，频率越小；自振周期与刚度的平方根成反比，刚度越大，周期越小频率

3、越大；要改动构造的自振周期，只需从改动周期越小频率越大；要改动构造的自振周期，只需从改动构造的质量或刚度着手。构造的质量或刚度着手。3 3两个外形类似的构造，假设周期相差悬殊，那么动力性能两个外形类似的构造，假设周期相差悬殊，那么动力性能相差很大。反之，两个外形看来并不一样的构造，假设其自振相差很大。反之，两个外形看来并不一样的构造，假设其自振周期相近，那么在动荷载作用下的动力性能根本一致。周期相近，那么在动荷载作用下的动力性能根本一致。stgWgmmk=w1gkmTst=22例例4 4、图示三根单跨梁，、图示三根单跨梁，EIEI为常数，在梁中点有集中质量为常数，在梁中点有集中质量m m，不思

4、索梁的质量，试比较三者的自振频率。不思索梁的质量，试比较三者的自振频率。l/2l/2l/2l/2l/2l/2mmm解：解：1 1求求EIl4831=P=13l/165l/32P=1l/2EIlllllEIl7687)325216322(61321=EIl768732=EIl19233=311481mlEIm=w32277681mlEIm=w3331921mlEIm=w据此可得：据此可得：1? 2 ? 3= 1 ? 1.512 ? 2 构造约束越强构造约束越强, ,其刚度越大其刚度越大, ,刚度越大刚度越大, ,其自振动频率也越大。其自振动频率也越大。1例例5 5、求图示构造的自振圆频率。、求图

5、示构造的自振圆频率。解法解法1 1：求：求 k k=1/hMBA=kh = MBCklhmIEIBAClhEIlEI33=lmhEImk23=w23lhEIk =1h解法解法2 2：求：求 EIlhhlhEI332212=21131mlhEIm=w例例6 6、求图示构造的自振频率。、求图示构造的自振频率。lEImk1k11k11k33lEI解：求解：求 k3113lEIkk+=mkmklEI+=3311w对于静定构造普通计算柔度系数方便。对于静定构造普通计算柔度系数方便。假设让振动体系沿振动方向发生单位位移时，一切刚节点假设让振动体系沿振动方向发生单位位移时，一切刚节点都不能发生转动如横梁刚度

6、为都不能发生转动如横梁刚度为刚架刚架) )计算刚度系数方便。计算刚度系数方便。312lEI一端铰结的杆的侧移刚度为：一端铰结的杆的侧移刚度为：33lEI两端刚结的杆的侧移刚度为：两端刚结的杆的侧移刚度为：五、阻尼对自在振动的影响五、阻尼对自在振动的影响忽略阻尼影响时所得结果忽略阻尼影响时所得结果 能不能能不能 反映实践构造的振动规律。反映实践构造的振动规律。大体上大体上忽略阻尼的振动规律思索阻尼的振动规律构造的自振频率是构造的固有特性，与外因无关。构造的自振频率是构造的固有特性，与外因无关。简谐荷载作用下有能够出现共振。简谐荷载作用下有能够出现共振。自在振动的振幅永不衰减。自在振动的振幅永不衰

7、减。自在振动的振幅逐渐衰减。自在振动的振幅逐渐衰减。共振时的振幅趋于无穷大。共振时的振幅趋于无穷大。共振时的振幅较大但为有限值。共振时的振幅较大但为有限值。f产生阻尼的缘由：构造与支承之间的外摩擦；资料之间的内摩产生阻尼的缘由：构造与支承之间的外摩擦；资料之间的内摩擦；周围介质的阻力。擦；周围介质的阻力。f阻尼力确实定：总与质点速度反向阻尼力确实定：总与质点速度反向; ;大小与质点速度有如下关系：大小与质点速度有如下关系：f 与质点速度成反比比较常用，称为粘滞阻尼。与质点速度成反比比较常用，称为粘滞阻尼。f 与质点速度平方成反比如质点在流体中运动遭到的阻与质点速度平方成反比如质点在流体中运动遭

8、到的阻力。力。f 与质点速度无关如摩擦力。与质点速度无关如摩擦力。f粘滞阻尼力的分析比较简单，粘滞阻尼力的分析比较简单，( (由于由于R(t)=R(t)=Cy ).Cy ).其他阻尼力也可化为等效粘滞阻尼力来分析。其他阻尼力也可化为等效粘滞阻尼力来分析。f振动模型振动模型ykyym& &kmyc&y有阻尼的自在振动，动平衡方程：有阻尼的自在振动，动平衡方程：)14.15(0LLL& &=+kyycymwwmcmk2,=( 阻尼比)1(2=wl0222=+wwll)( =ltCety设解为：设解为：特征方程为：特征方程为：1)1低阻尼情况21wwwwl=r

9、ri其中tCtCeyrrtwwwsincos21+=+=tyvtyeyrrrtwwwwwsincos000c)16.15(022LLL& &=+yyyww令令低阻尼体系的自振圆频率低阻尼体系的自振圆频率000220020)()sin(yvytgyvyataeyrrrtwwawwaww+=+=+=ae-ttyty低阻尼低阻尼y- t曲线曲线无阻尼无阻尼y- t曲线曲线阻尼对自振频率的影响阻尼对自振频率的影响. .=而随www,12r 当0.2，那么存在0.96r/1。在工程构造问题中，假设0.011 强阻尼：不出现振动，实践问题不常见。2)=1临界阻尼情况)1(2=wl=w lt

10、etCCyw+=)(21tetvtyyww+=)1 (00)16.15(022LL& &=+yyywwtyy0000vtg=这条曲线仍具有衰减性，这条曲线仍具有衰减性，但不具有动摇性。但不具有动摇性。m15-3 单自在度体系的受迫振动 受迫振动强迫振动：构造在动力荷载作用下的振动。受迫振动强迫振动：构造在动力荷载作用下的振动。ky(t)ymkyym& &P(t )mP(t )P(t )弹性力弹性力ky、惯性力、惯性力ym& &和荷载和荷载P(t)之间的平衡方程为之间的平衡方程为:)()(atPkyymLL& &=+)24.15()

11、(2LL& &mtPyy=+w一、简谐荷载：一、简谐荷载：tmFtAwsinsin)(22=+tmFtAtAwsinsinsin22=+tAysin=mtFyywsin2=+& &)(22w+=mFAtytmFystwwwsin)1 (1sin)1 (22222=wFmFyst=2单自在度体系强迫单自在度体系强迫振动的微分方程振动的微分方程特解：特解：最大静位移最大静位移yst是把荷载幅值当作静荷载作用时构造所产生是把荷载幅值当作静荷载作用时构造所产生的位移。的位移。tyystwsin1122=特解可写为：特解可写为：通解可写为：通解可写为：tytCtCystw

12、wwsin11cossin2221+=设设t=0时的初始位移和初始速度均为零，那么：时的初始位移和初始速度均为零，那么：0,12221=CyCstww)sin(sin1122ttyystwww=过渡阶段：振动开场两种振动同时存在的阶段；过渡阶段：振动开场两种振动同时存在的阶段；平稳阶段：后来只按荷载频率振动的阶段。由于阻尼的存在平稳阶段：后来只按荷载频率振动的阶段。由于阻尼的存在按自振频率振动按荷载频率振动平稳阶段：平稳阶段：tyystwsin1122=最大动位移振幅为：最大动位移振幅为：22max11w=styy22max11w=styy动力系数动力系数为：为：1023123w重要的特性：重

13、要的特性：当当/0时时,1，荷载变化得，荷载变化得很慢，可当作静荷载处置。很慢，可当作静荷载处置。当当0 / 1，并且随，并且随/的增大而增大。的增大而增大。当当/ 1时时,。即当荷载频。即当荷载频率接近于自振频率时，振幅会无率接近于自振频率时，振幅会无限增大。称为限增大。称为“共振。通常把共振。通常把0.75 / 1时时,的绝对值随的绝对值随/f的增大而减小。当的增大而减小。当很大时，荷载变化很快，构造来不及反响。很大时，荷载变化很快，构造来不及反响。 当动荷载作用在单自在度体系的质点上时，由于体系上各当动荷载作用在单自在度体系的质点上时，由于体系上各截面的内力、位移都与质点处的位移成正比，

14、故各截面的最大截面的内力、位移都与质点处的位移成正比，故各截面的最大动内力和最大动位移可采用一致的动力系数，只需将干扰力幅动内力和最大动位移可采用一致的动力系数，只需将干扰力幅值乘以动力系数按静力方法来计算即可。值乘以动力系数按静力方法来计算即可。 例：知例：知m=300kg，EI=90105N.m2 ,k=48EI/l3 ,P=20kN,=80s-1 求梁中点的位移幅值及最大动力弯矩。求梁中点的位移幅值及最大动力弯矩。2mEImkPsint2m解：解：1)1)求求kEIl212148321+=+=EIlEIlEIl192519248333=+=1316.13451921=smlEImw2)2

15、)求求552. 11122=wmEIlPPy35333max1075. 51090192451020552. 11925=3)3)求求ymaxymax， Mmax MmaxmkNlPM.04.31420552. 141)(41max=例、一简支梁例、一简支梁I28b，惯性矩，惯性矩I=7480cm4，截面系数，截面系数W=534cm3，E=2.1104kN/cm2。在跨度中点有电动机分量。在跨度中点有电动机分量Q=35kN，转速，转速n=500r/min。由于具有偏心，转动时产生离心力。由于具有偏心，转动时产生离心力P=10kN，P的竖向分量为的竖向分量为Psint。忽略梁的质量，试求强迫振动

16、。忽略梁的质量，试求强迫振动的动力系数和最大挠度和最大正应力。梁长的动力系数和最大挠度和最大正应力。梁长l=4m解：解：1求自振频率和荷载频率求自振频率和荷载频率 SQlEIg13434 .57400359807480101 . 24848=Sn13 .526050014. 32602=2 2求动力系数求动力系数88. 54 .573 .5211112222=wEIPlEIQlystst484833max+=+=+=+=WlPQWPlWQl4)(44maxstg=w175.6MPa 必需特别留意，这种处置方法只适用于单自在度体系在质点上受干扰力作用的情况。对于干扰力不作用于质点的单自在度体系，

17、以及多自在度体系，均不能采用这一方法。I22b3570cm4357039.739.71.35对于本例，采用较小的截面的梁既可防止共振，又能获得较好的经济效益。325149.2设体系在设体系在t=0时静止，时静止，然后有瞬时冲量然后有瞬时冲量S作用。作用。二、普通荷载二、普通荷载普通荷载作用下的动力反响可利用瞬时冲量的普通荷载作用下的动力反响可利用瞬时冲量的动力反响来推导动力反响来推导1、瞬时冲量的动力反响、瞬时冲量的动力反响P(t)tP瞬时冲量瞬时冲量S S引起的振动可视为引起的振动可视为由初始条件引起的自在振动。由初始条件引起的自在振动。由动量定理：由动量定理：mtPmSv=000=yt c

18、ossin)(00www+=tvtytytmStywwsin)( =t ttttmSty=wwsin)()(sinww=tmScossin)(00www+=tvtytycossin)(00www+=tvtytycossin)(00www+=tvtytycossin)(00www+=tvtytytPSmv=002、恣意荷载、恣意荷载P(t)的动力反响的动力反响P(t)tdPdS)(=时辰的微分冲量对时辰的微分冲量对t瞬时瞬时(t )引起的动力反响引起的动力反响:)(sin)(ww=tmdPdy初始静止形状的单自在初始静止形状的单自在度体系在恣意荷载作用度体系在恣意荷载作用下的位移公式下的位移公式

19、: :wwdtPmtyt)(sin)(1)(0=(Duhamel 积分积分)15.29初始位移初始位移y0y0和初始速度和初始速度v0v0不为零在恣意荷载作用下的位移公式不为零在恣意荷载作用下的位移公式: :wwwwdtPmtvtytyt)(sin)(1sincos)(000+=t3、几种典型荷载的动力反响、几种典型荷载的动力反响1 1突加荷载突加荷载 =0,0, 0)(0tPttP当当P(t)tPwwdtPmtyt)(sin)(1)(0=wwdtPmtyt)(sin1)(00=)cos1 ()cos1 (20tytmPstwww=yst=P0=P0 /m2ysty(t)t023质点围绕静力平

20、衡质点围绕静力平衡位置作简谐振动位置作简谐振动2)(max=styty2 2短时荷载短时荷载 =ututPttP, 00,0, 0)(0P(t)tPu阶段阶段(0tu)：无荷载，体系以：无荷载，体系以t=u时辰的位移时辰的位移 和速度和速度为初始条件作自在振动。为初始条件作自在振动。)cos1 ()(uyuystw=uyuvstwwsin)(=sincos )(00www+=tvtyty)(sinsin)(cos)cos1 ()(utuyutuytystst+=wwww)cos)(costutystww=或者直接由或者直接由Duhamel积分作积分作wwdtPmtyt)(sin)(1)(0=w

21、wdtPmtyu)(sin1)(00=)cos)(cos20tutmPwww=)2(sin2sin2utuyst=ww另解：短时荷载可以为由两个突加荷载叠加而成。另解：短时荷载可以为由两个突加荷载叠加而成。P(t)tPP(t)tPuP(t)tPu)cos1 ()(tytystw=)(cos1 ()(utytyst=w当当0t u)cos1 ()(tytystw=)(cos1 (utystw)cos)(costutystww=)2(sin2sin2utuyst=wwysty(t)t023T最大动反响最大动反响1当 u T/2 最大动位移发生在阶段)cos1 ()(tytystw=styy2max

22、=2当u T/2 最大动位移发生在阶段 =2)2(sin2sin2)(utuytyst=ww2sin2maxuyystw=2sin2uw=21, 221,sin2TuTuTu当当Tu1/611/22动力系数反响谱动力系数反响谱(与与T和和u之间的关系曲线之间的关系曲线)3 3线性渐增荷载线性渐增荷载 =rrrttPttttPtP当当,0,)(00P(t)tP0tr这种荷载引起的动力反响同样可由这种荷载引起的动力反响同样可由DuhamelDuhamel积分来求解积分来求解: :=rrrstrrstttttttytttttyty当当,)(sinsin11,sin)(wwwww 对于这种线性渐增荷载

23、对于这种线性渐增荷载, ,其动力反响与升载时间的长短有其动力反响与升载时间的长短有很大关系。其动力系数的反响谱如下：很大关系。其动力系数的反响谱如下：01.02.03.04.0TtrtrP0动力系数反响谱动力系数反响谱动力系数动力系数介于介于1 1与与2 2之间。之间。假设升载很短，假设升载很短，trT/4,tr4T,tr4T,那么那么接近于接近于1,1,即相当于静荷载情况。即相当于静荷载情况。常取外包虚线作为设计的根据。常取外包虚线作为设计的根据。三、有阻尼的强迫振动三、有阻尼的强迫振动单独由单独由v0v0引起的自引起的自在振动：在振动：瞬时冲量瞬时冲量d

24、s=Pdt=mv0ds=Pdt=mv0所引起的振动，可视为所引起的振动，可视为以以v0=Pdt/mv0=Pdt/m，y0=0y0=0为初为初始条件的自在振动：始条件的自在振动：tveyrrtwwwsin0=tmPdteyrrtwwwsin=将荷载将荷载P(t)P(t)的加载过程的加载过程 看作一系列瞬时冲量：看作一系列瞬时冲量：)(sin)()(www=temdPdyrtr总反响总反响wwwdtemPtyrttr)(sin)()()(0=+tyvtyerrrtwwwwwsincos000P(t)tddPdS)(=tEI=m例、图示一单层建筑物的计算简图。屋盖系统和柱子的质量均集例、图示一单层建

25、筑物的计算简图。屋盖系统和柱子的质量均集中在横梁处合计为中在横梁处合计为m9.8kN ，加一程度力，加一程度力P=9.8kNP=9.8kN，测得侧移，测得侧移A0=0.5cmA0=0.5cm，然后忽然卸载使构造发生程度自在振动。在测得周期然后忽然卸载使构造发生程度自在振动。在测得周期T=1.5s T=1.5s 及一及一个周期后的侧移个周期后的侧移A1=0.4cmA1=0.4cm。求构造的阻尼比。求构造的阻尼比和阻尼系数和阻尼系数c c。解：解：0335. 04 . 05 . 0ln21ln211=+kkyymNAPk/10196005. 0108 . 9430=1189. 45 . 122=s

26、Tw=wk2=wmc 2=wwm22cmsNmsN/2 .332/33220189. 4101960355. 024=1突加荷载P0)sin(cos1 )(20ttemPtyrrrtwwwwww=低阻尼低阻尼y- t曲线曲线无阻尼无阻尼y- t曲线曲线ysty(t)t02345y(t)t02345静力平衡位置具有阻尼的体系在具有阻尼的体系在突加荷载作用下，突加荷载作用下，最初所引起的最大最初所引起的最大位移接近于静位移位移接近于静位移yst=P0/m2的两倍，的两倍，然后逐渐衰减，最然后逐渐衰减，最后停留在静力平衡后停留在静力平衡位置。位置。2简谐荷载P(t)=Fsint)34.15(sin2

27、2LL& &tmFyyyww=+设特解为：设特解为：y=Asin t +Bcos t代入代入15-34得：得：222222222222224)(2,4)(wwwwww+=+=mFBmFAsincos21tCtCeyrrtwww+=+Asin t +Bcos t 齐次解加特解得到通解：齐次解加特解得到通解：自在振动，因阻尼作用，逐渐衰减、消逝。纯强迫振动，平稳振动，振幅和周期不随时间而变化。结论：在简谐荷载作用下，无论能否计入阻尼的作用，纯结论：在简谐荷载作用下，无论能否计入阻尼的作用，纯 强迫振动部分总是稳定的周期运动，称为平稳振动。强迫振动部分总是稳定的周期运动，称为平稳振动。y=Asin t +Bcos t =yPsin(t ) (15-35a)2122222222)(1)(2,4121wwaww=+=+=tgyBAystP振幅：振幅：yp，最大静力位移：最大静力位移：yst=F/k=F/m2=stPyy2122222241+ww=stPyy2122222241+ww动力系数动力系数与频率比与频率比/和阻尼比和阻尼比有关有关4.03.02.01.001.02.03.0/=0=0.1=0.2=0.3=0.5=1.0几点留意：几点留意：随随增大增大曲线渐趋平缓，曲线渐趋平缓， 特别是在特别是在/=1附近附近的的