20XX_20XX版高中物理第1章机械振动章末总结学案沪科版.docx

《20XX_20XX版高中物理第1章机械振动章末总结学案沪科版.docx》由会员分享，可在线阅读，更多相关《20XX_20XX版高中物理第1章机械振动章末总结学案沪科版.docx（53页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、20XX_20XX版高中物理第1章机械振动章末总结学案沪科版2022_2022版中学物理第1章机械振动章末总结学案沪科版 本文关键词：振动，中学物理，学案，机械，_20222022_2022版中学物理第1章机械振动章末总结学案沪科版 本文简介：第1章机械振动章末总结一、简谐运动的图像及作用简谐运动的图像描述了振动质点的位移随时间改变的规律从图像中可以确定位移、速度、加速度、动能和势能等物理量以及它们的改变规律例1(多选)如图1所示，下列说法正确的是()图1A振动图像上的A、B两点振动物体的速度相同B在t0.1s和t0.3s时，2022_2022版中学物理第1章机械振动章末总结学案沪科版 本文内

2、容：第1章机械振动章末总结一、简谐运动的图像及作用简谐运动的图像描述了振动质点的位移随时间改变的规律从图像中可以确定位移、速度、加速度、动能和势能等物理量以及它们的改变规律例1(多选)如图1所示，下列说法正确的是()图1A振动图像上的A、B两点振动物体的速度相同B在t0.1s和t0.3s时，质点的加速度大小相等，方向相反C振动图像上A、B两点的速度大小相等，方向相反D质点在t0.2s和t0.3s时的动能相等答案BC解析A、B两点位移相同，速度大小相等，但方向相反，因此A错，C对在t0.1s和t0.3s时，质点离开平衡位置的位移最大，方向相反，由Fkx，a可知B对t0.2s时，物体通过平衡位置，

3、速度最大，动能最大，而t0.3s时，物体在最大位移处，速度为零，动能最小，故D错二、简谐运动的周期性和对称性1周期性做简谐运动的物体在完成一次全振动后，再次振动时则是重复上一个全振动的形式，所以做简谐运动的物体具有周期性做简谐运动的物体经过同一位置可以对应不同的时刻，物体的位移相同而速度可能等大反向，这样就形成简谐运动的多解问题2对称性(1)速率的对称性：系统在关于平衡位置对称的两位置具有相等的速率(2)加速度和回复力的对称性：系统在关于平衡位置对称的两位置具有等大反向的加速度和回复力(3)时间的对称性：系统通过关于平衡位置对称的两段位移的时间相等振动过程中通过随意两点A、B的时间与逆向通过的

4、时间相等例2物体做简谐运动，通过A点时的速度为v，经过1s后物体第一次以相同速度v通过B点，再经过1s物体紧接着又通过B点，已知物体在2s内所走过的总路程为12cm，则该简谐运动的周期和振幅分别是多大？答案T4s，A6cm或Ts，A2cm解析物体通过A点和B点时的速度大小相等，A、B两点肯定关于平衡位置O对称依题意作出物体的振动路径草图如图甲、乙所示，在图甲中物体从A向右运动到B，即图中从1运动到2，时间为1s，从2运动到3，又经过1s，从1到3共经验了0.5T，即0.5T2s，T4s，2A12cm，A6cm.在图乙中，物体从A先向左运动，当物体第一次以相同的速度通过B点时，即图中从1运动到2

5、时，时间为1s，从2运动到3，又经过1s，同样A、B两点关于O点对称，从图中可以看出从1运动到3共经验了1.5T，即1.5T2s，Ts，1.54A12cm，A2cm.针对训练(多选)一个弹簧振子的振幅是A，若在t时间内物体运动的路程是s，则下列关系中可能正确的是(包括肯定正确的)()At2T，s8ABt，s2ACt，s2ADt，sA答案ABD三、单摆周期公式的应用单摆的周期公式T2是在当单摆的最大偏角不大于5，单摆的振动是简谐运动的条件下才适用的，单摆的周期与振幅无关，与质量也无关，只与摆长和重力加速度有关另外由公式的变形式g还可以测重力加速度例3一个摆长为l1的单摆，在地面上做简谐运动，周期

6、为T1，已知地球质量为M1，半径为R1.另一摆长为l2的单摆，在质量为M2，半径为R2的星球表面做简谐运动，周期为T2.若T12T2，l14l2，M14M2，则地球半径与星球半径之比R1R2为()A21B23C12D32答案A解析在地球表面单摆的周期T12.在星球表面单摆的周期T22.又因为g，Gg由联立得.1一质点做简谐运动的图像如图2所示，下列说法正确的是()图2A质点振动频率是4HzB在10s内质点经过的路程是20cmC第4s末质点的速度为零D在t1s和t3s两时刻，质点位移大小相等，方向相同答案B解析振动图像表示质点在不同时刻相对平衡位置的位移，由题图可看出，质点运动的周期T4s，其频

7、率f0.25Hz，A错误；10s内质点运动了T，其运动路程为sT4A42cm20cm，B正确；第4s末质点在平衡位置，其速度最大，C错；t1s和t3s两时刻，由题图可看出，位移大小相等，方向相反，D错2(多选)关于简谐运动的周期，下列说法正确的是()A间隔一个周期的两个时刻，物体的振动状况完全相同B间隔半个周期奇数倍的两个时刻，物体的速度和加速度可能同时相同C半个周期内物体动能的改变肯定为零D一个周期内物体势能的改变肯定为零答案ACD解析依据周期的意义知，物体完成一次全振动，全部的物理量都复原到初始状态，所以A、D正确；当间隔半个周期的奇数倍时，全部的矢量都变得大小相等、方向相反，故B选项错误

8、；由于间隔半个周期各矢量大小相等，所以物体的动能必定相等，没有改变，所以C也正确3站在升降机里的人发觉，升降机中摇摆的单摆周期变大，以下说法正确的是()A升降机可能加速上升B升降机肯定加速上升C升降机可能加速下降D升降机肯定加速下降答案C解析由单摆周期公式T2知，周期变大，则等效重力加速度g变小，故升降机的加速度方向向下，故可能加速下降，也可能减速上升，故选项C正确4(多选)一弹簧振子沿x轴振动，平衡位置在坐标原点t0时振子的位移x0.1m；ts时x0.1m；t4s时x0.1m该振子的振幅和周期可能为()A0.1m，sB0.1m,8sC0.2m，sD0.2m,8s答案ACD解析若振幅A0.1m

9、，Ts，则s为半个周期，从0.1m处运动到0.1m处，符合运动实际，4sss为一个周期，正好返回0.1m处，所以A对；若A0.1m，T8s，s只是T的，不行能由负的最大位移处运动到正的最大位移处，所以B错；若A0.2m，Ts，则s，振子可以由0.1m处运动到对称位置，4sssT，振子可以由0.1m处返回0.1m处，所以C对；若A0.2m，T8s，则s2，而sin，即时间内，振子可以从平衡位置运动到0.1m处，再经s又恰好能由0.1m处运动到0.2m处后，再返回0.1m处，所以D对5一个摆长为2m的单摆，在地球上某处振动时，测得完成100次全振动所用的时间为284s.(1)求当地的重力加速度g；

10、(2)把该单摆拿到月球上去，已知月球上的重力加速度是1.60m/s2，则该单摆的振动周期是多少？答案(1)9.78m/s2(2)7.02s解析(1)周期Ts2.84s由周期公式T2得gm/s29.78m/s2.(2)T223.14s7.02s.篇2：桥梁结构振动与限制分析探讨结题报告桥梁结构振动与限制分析探讨结题报告 本文关键词：分析探讨，振动，桥梁，结题，结构桥梁结构振动与限制分析探讨结题报告 本文简介：第五期SRTP结题报告桥梁结构振动与限制分析探讨第五期SRTP结题报告桥梁结构振动与限制分析探讨指导老师：沈火明教授项目组成员：张徐20224968李恒20224977张冰清20225379

11、桥梁结构振动与限制分析探讨一、课题的探讨意义及探讨方法1.1课题的探讨意义桥梁结构的振动是引起桥梁损坏(破坏)桥梁结构振动与限制分析探讨结题报告 本文内容：第五期SRTP结题报告桥梁结构振动与限制分析探讨第五期SRTP结题报告桥梁结构振动与限制分析探讨指导老师：沈火明教授项目组成员：张徐20224968李恒20224977张冰清20225379桥梁结构振动与限制分析探讨一、课题的探讨意义及探讨方法1.1课题的探讨意义桥梁结构的振动是引起桥梁损坏(破坏)的一个重要因素，引起桥梁振动的因素主要有：地震引起的振动、荷载引起的振动及车-桥耦合作用引起的振动。传统的结构强度设计方法通过增加结构物自身抗力

12、来抵挡地震作用，即由结构本身储存和消耗地震能量。但由于人类测震技术的不成熟，尚不能精确估计振动的强度和特性。因此，可能会出现结构不满意平安性的要求而产生平安事故。近年来发展起来的结构限制技术是建筑结构抗震领域内的一个新的探讨热点,它是通过采纳结构振动限制的理论与方法变更结构系统的动力学性能或阻尼耗散性能来增加和改善结构的抗震实力，是一种主动主动的对策。因此，近年来桥梁结构的振动限制倍受学术界、工程界的广泛关注,并获得了长足的进步。结构振动与限制的探讨与应用有着广泛的前景，它的探讨和发展将给结构工程抗震设计带来一张革命，其巨大的经济效益和社会效益已得到证明。1.2本文的探讨思路和方法本文以竖向弯

13、曲振动时桥梁跨中挠度的振动幅度为限制目标，通过分析安装TMD前后桥梁跨中挠度的振动幅度改变量来探讨TMD对桥梁振动的限制效果，并探究TMD的参数优化。对于TMD限制下的车桥耦合系统的振动，已被采纳的数值探讨方法有两类：一类在建立系统耦合方程组的基础上，借助编程语言或数值计算软件（如MATLAB和VB等），利用数值积分方法编程求解耦合方程组；另一类，借助仿真分析软件（如有限元软件和Simulink等）实现对系统的仿真分析。本文将实行第一类方法，编程求解方程组。针对简支梁桥在列车匀速通过时的竖向弯曲振动，本文先建立车桥耦合振动理论模型，利用数值计算方法结合MATLAB软件，编程求解车桥时变系统振动

14、微分方程组，获得列车过桥时桥梁竖向振动位移响应；再建立车桥-TMD耦合振动理论模型，求解获得单个以及多个TMD限制下的桥梁竖向位移响应，分析TMD的限制效果，并探讨TMD参数优化对限制效果的影响。1、理论模型及求解方法2.1车辆-简支梁桥竖向振动模型2.1.1模型简化模型中，车体、转向架、车轮均被认为是刚体，相互间通过弹簧阻尼系统连接；不考虑车轮与钢轨表面粗糙，认为钢轨固结于桥梁上作为桥梁的一部分，不考虑钢轨及轨下结构局部变形造成的影响，于是轮轴的竖向位移等于轮轨接触点（即车轮与桥梁接触点）的竖向位移。由于仅考虑系统的竖向振动，本文中车辆简化为二系弹簧悬挂系统，仅考虑车体沉浮、点头，前后转向架

15、构架沉浮运动，每节车辆四个自由度。桥梁采纳简支欧拉梁模型。系统简化模型及坐标系建立如图1所示。图1系统简化计算模型图中各物理量意义如下：v-行车速度；yc-车体中心竖向沉浮位移；i-第i节车体点头位移；Jc-车体点头惯量；mc-车体质量；cs2-二系悬挂阻尼；ks2-二系悬挂刚度；mt-转向架构架与轮对质量之和；yt-构架中心竖向沉浮位移；cs1-一系悬挂阻尼；ks1-一系悬挂刚度；a-同一节车转向架中心距离；EI-桥梁抗弯刚度；d-前后两节车相邻轮对间距；L-桥梁全长；w-桥梁挠度，以水平位置为坐标起点；mk-TMD质量；kk-TMD刚度；Ck-TMD阻尼；yz-TMD竖向位移；li-第i个

16、轮对与第一个轮对之间的距离；2.1.2车辆系统振动微分方程第i节车振动方程车体沉浮运动：(2-1)车体点头运动：（2-2）构架沉浮运动：(2-3)(2-4)对车厢整体可得（2-5）（2-6）其中p2i-1，p2i分别为第i节车厢前后两个车轮与桥梁之间的作用力。2.1.3简支梁桥振动方程本为采纳简支欧拉梁模型，不考虑桥梁阻尼时，振动方程为：（2-7)其中，A-桥梁单位长度的质量；F(x,t)-t时刻x处作用在梁上的外力，包括桥梁自重和轮轨相互作用力，即（2-8）其中，N表示车辆节数；表示dirac函数；将式（2-8）代入式（2-7）得，（2-9）为求解方程（2-9），利用分别变量法设（2-10）

17、其中为简支梁的振型函数，n为模态截断数，Tj(t)形态振幅函数。将式（2-10）代入式（2-9），各项自0到L积分，利用振型函数的正交性与dirac函数的性质，并令,得j=1,2,.,n(2-11)将式（2-5）与式（2-6）代入式(2-11),消去变量Pi，得j=1,2,.,n(2-12)为简化表达，令，整理得，j=1,2,.,n(2-13)这样，式（2-1）式（2-4）与式(2-13)一起组成车桥时变系统耦合振动微分方程组，方程组共有（4N+n)个方程，以yci,yt2i-1,yt2i,i,Tj(t)共（4N+n）个未知量为求解变量，结合初值条件，可利用动力学连续数值积分方法联合求解。2.

18、2车桥TMD耦合模型（以跨中悬挂单个TMD为例）设TMD质量为mk，弹簧刚度为kk，阻尼为ck，悬挂位置为跨中，以静平衡位置为坐标起点，振动位移为yz，则TMD的运动方程为（2-14）简支梁跨中悬挂TMD时，车辆系统振动方程不受影响，简支梁所受外力应考虑梁与TMD之间的相互作用，即通过TMD弹簧与阻尼器传递的力，于是，式（2-8）变为（2-15）其中，为dirac函数。此时，式（2-13）变为j=1,2,.,n（2-16）这样，式（2-1）式（2-4）与式(2-14)、式（2-16）一起组成车桥TMD时变系统耦合振动微分方程组，方程组共有（4m+n+1）方程，以yci,yt2i-1,yt2i,

19、i,Tj(t),yz共（4m+n+1）个未知量为求解变量，结合初值条件，可利用动力学连续数值积分方法联合求解。2.3求解步骤与方法（以单个TMD限制系统为例）车桥TMD系统耦合方程组写成矩阵形式为：（2-17）其中，x-以yci,yt2i-1,yt2i,i,Tj(t),yz构成的未知向量，即分块表示为：；M,C,K-总体等效质量,阻尼,刚度矩阵;f-等效载荷向量；2.3.1系数矩阵a)质量矩阵分块表示其中,；,I表示单位矩阵；mz=mk;,;其中，限于篇幅，hij(t)简写为hij.b)阻尼矩阵分块表示为,;，c)刚度矩阵分块表示为,;，d)载荷向量分块表示为;,.2.3.2求解方法本文采纳N

20、ewmark-法，利用MATLAB软件编程求解矩阵方程组（2-17），获得系统的位移响应。三、实际算例依据（二）中建立的计算模型与求解方法，本文给出以下实际算例。实例探讨一辆10节编组的列车通过一座简支梁桥时，引起的桥梁振动，并利用TMD限制桥梁振动。列车采纳德国ICE动车和拖车，前后2节动车中间8节拖车编组，桥梁为一全长32m简支梁桥。详细车辆与桥梁参数见附录。3.1桥梁静挠度依据材料力学中简支梁在均布压力下跨中挠度的计算公式，本文算例中的桥梁在自重下跨中挠度为：3.2实施限制前，列车过桥时引起的跨中振动响应如图2，给出了列车以100km/h车速匀速通过时，简支梁跨中挠度的时程曲线。t=0s

21、时，列车起先上桥，约9s时，列车刚好完全离开桥梁。图2100km/h车速下跨中振动响应从图中可以看出，当列车以100km/h的速度通过时，跨中最大挠度发生在列车刚上桥后，为60.4mm，相对于跨中静挠度27.9mm超出一倍以上，振幅为32.5mm，这对于桥梁结构本身和列车的平安都是有害的也是危急的，因此须要对其进行限制。图中，t=9s后列车离开桥梁，跨中位移并没有衰减，而是接着作以静挠度为均值，以列车刚离开桥瞬间振幅为幅值的简谐振动。这是因为算例中的简支梁没有考虑桥梁本身的阻尼，当列车离开桥梁后，桥梁自由振动，此后跨中响应仅确定于其初值条件即车离开桥时跨中的振动状况，由于不存在阻尼这种振动将持

22、续不会衰减，但这仅是算例中假设前提下的状况。为探讨桥梁跨中位置在列车通过时相对于静平衡位置的振动，将图3.1中的时程曲线沿纵坐标向下平移27.9mm，得到图3如下：图3跨中振幅时程曲线3.3实施限制前，跨中振幅随车速的改变状况本为计算了列车以60km/h-200km/h不同车速通过桥梁时，跨中振幅的改变状况，如图4所示。图4桥梁的振幅随车速的改变曲线从图中可以看出，随着列车过桥车速的增加，简支梁跨中振幅基本呈增加趋势，个别车速下振幅出现局部峰值，车速为185km/h时振幅甚至接近42mm。3.4利用单个TMD限制桥梁在列车通过时的振动本文以车速为100km/h为例，探讨单个TMD对桥梁振动的限

23、制作用。计算中的TMD质量、阻尼和刚度参数由DenHartog参数调整公式（式3-1）给出。，（3-1）其中，为TMD质量比，cz和cc分别为TMD的阻尼系数和临界阻尼系数。如图5所示，给出了车速100km/h，跨中悬挂质量比为0.08%的单个TMD时，桥梁跨中振幅的时程曲线。图5单个TMD限制下的跨中振动曲线从图中可以看出，跨中悬挂质量比为0.08%的单个TMD时，跨中振幅为30.6mm。相对于限制前的32.5mm振幅，限制效果为5.85%。若变更TMD质量比，将获得不同的限制效果。于是，以限制效果最大为限制目标，可获得最佳TMD质量比参数。同时，图中看出当列车离开桥梁后，跨中位移呈现向静平

24、衡位置衰减趋势，这是由于模型中考虑了TMD阻尼的作用。3.5单个TMD限制的最佳质量比如图6所示，本文给出了车速为100km/h时，单个TMD的限制效果随TMD质量比改变曲线。图6单个TMD的限制效果随TMD质量比改变曲线从上图可以看出质量比=0.05%时，获得最佳限制效果为24.3%，此时跨中振幅为24.6mm，最大挠度52.5mm，该TMD减振效果显著。3.6实现MTMD对桥梁振动的限制MTMD形式多样，涉及的参数也多如：TMD的悬挂位置，各TMD的质量，阻尼，刚度，TMD相互间的频率间隔等。本文算例仅探讨频率呈线性分布，等间距悬挂的5个TMD的限制作用。图7给出了等质量比，频率间线性分布

25、的5个TMD等间距悬挂时，跨中振动的时程曲线，其中，MTMD参数为质量比=0.001%，频率间隔df=3Hz。图7MTMD限制下的跨中振幅曲线从上图可以看出，5个TMD限制时，跨中振幅为24.2mm，限制效果为25.54%，相对于最佳质量比的单个TMD限制效果虽然只增加了1.24%，但总质量却只有单个TMD的1/10，这特别有利于降低悬挂TMD对桥梁静挠度的影响。同时，假如调整MTMD的参数，可以获得不同的限制效果。于是，以限制效果最大为目标，以质量比和频率间隔为优化变量，依据优化理论，利用二维优化搜寻方法可以获得MTMD的最佳参数。由于涉及的程序较大，本文受条件限制未能利用该方法探讨MTMD

26、的最佳参数，只是给出了保持总质量比的前提下限制效果随频率间隔的改变状况（如图8）。图8MTMD限制效果随频率间隔的改变曲线从上图可以看出，限制总质量为0.05%，当频率间隔较低时，随着df增加，限制效果呈增加趋势但改变不明显；当df=11.5Hz时，限制效果取最大值；当频率间隔高于11.5HZ时，限制效果明显快速衰减。所以取频率间隔为11.5Hz，可以获得最佳限制效果为26.1%。3.7算例结论3.13.6的算例中，本文运用（二）中建立的模型以及求解方法，实现了TMD和MTMD对简支梁桥在列车匀速通过时的振动响应的限制。四、总结本文探讨了列车过桥时引起的车桥耦合振动，建立了车辆简支梁桥与车辆简

27、支梁桥TMD耦合振动微分方程，推导了方程的求解过程。探讨了TMD对桥梁振动的限制作用，利用数值计算探讨了TMD的参数优化对限制效果的影响。在建立理论模型的基础上，运用该模型计算了实际车桥耦合振动响应，并实现了TMD对桥梁的限制作用。但本文中的模型还非常粗糙，与实际状况差距较大，若要指导实际工程应用尚需更加符合实际的细化模型。同时，本文没有验证数值计算结果的牢靠性。参考文献1肖盐平.桥梁振动限制的初步探讨.西南交通高校探讨生学位论文，2022，12.2王均刚,马汝建.TMD振动限制结构的发展及应用.济南高校学报(自然科学版),2022，04.3杨宜谦,张煅,周宏业,孙宁.用调频质量阻尼器抑制铁路

28、桥梁竖向共振的探讨.中国铁道科学,1998年3月.4郭文华,路萍.TMD对高速列车通过简支箱梁桥时的振动限制探讨.振动与冲击,2022,12.5许磊平,刘伟庆等.TMD抑制桥梁振动仿真分析.防灾减灾工程学报.2022,6.6黄维平,强士中.大跨度悬索桥的双向TMD振动限制.桥梁建设,2000年第4期7张红艳,白长青,许庆余.多自由度困难结构的TMD调谐减震限制探讨.应用力学学报,2022,12.8张晶,王志强.调谐质量阻尼器(TMD)在大跨斜拉桥减震限制中的应用.工业建筑,20229肖艳平,沈火明,叶献辉.两种桥梁振动限制方法的对比分析.噪声与振动限制,2022年8月10李小珍,强士中.列车-

29、桥梁耦合振动探讨的现状与发展趋势.铁道学报，2002年10月。11郝超，强士中.移动荷载作用下桥梁的振动限制.国外桥梁，1999.12左一舟.关于列车-轨道(桥梁)时变系统空间振动方程的建立及其求解.铁道科学与工程学报,2022年2月13晋智斌,强士中,李小珍.高速列车-桥梁竖向随机振动的时域分析方法.地震工程与工程振动，2022,06.附录一：实际算例中车桥参数表桥梁参数全长L=32m弯曲刚度EI=5.18e10单位长度质量m=1.08e4车辆参数表动车参数拖车参数车体质量/kg5.88e44.55e4构架质量/kg6.45e33.87e3车体的点头刚度/kg*m23.089e62.391e

30、6二系阻尼系数/kg/s9e40.292e5一系阻尼系数/kg/s3e40.219e5二系弹簧刚度/N/m1.52e60.324e6一系弹簧刚度/N/m3.418e62.82e7轮对间隔/m11.4617附录二：数值计算程序代码%用纽马克法计算结构的动力学运动方程，返回值为简支梁跨中挠度%单个TMD限制functionw_mid=bz2t1(vv,t,mu)globalLvax1ax2dax1dax2;rn=5;n=rn;%所考虑桥梁振型的前n阶nn=n+41;g=9.81;%重力加速度-m/s2EI=5.18e10;m=10.8e3;L=32;v=vv/3.6;%列车移动荷载的速度-km/h

31、(-m/s)Ms1=5.88e4;%车体质量-kgMp1=6.45e3;%构架质量（含轮对质量）-kgIs1=3.089e6;%车体的点头刚度-kg*m2Cs1=0.9e5;%二系阻尼系数-kg/sCp1=0.3e5;%一系阻尼系数-kg/sKs1=1.52e6;%二系弹簧刚度-N/mKp1=2.418e6;%一系弹簧刚度-N/max1=11.46;%轮对间隔-mMs2=4.55e4;%车体质量-kgMp2=3.87e3;%构架质量（含轮对质量）-kgIs2=2.391e6;%车体的点头刚度-kg*m2Cs2=0.292e5;%二系阻尼系数-kg/sCp2=0.219e5;%一系阻尼系数-kg

32、/sKs2=0.324e6;%二系弹簧刚度-N/mKp2=28.2e6;%一系弹簧刚度-N/max2=17;%轮对间隔-mdax1=3.67+3.5;dax2=3.67*2;c1=EI/m/L4*pi4;%桥梁自由振动频率-1/sc2=pi/L;%c2*v车辆行进谐振频率-1/sc3=2*Mp1/m/L;c4=Ms1/m/L;c5=2*Is1/m/L/ax1;c6=(2*Mp1+Ms1)*g/m/L;c7=2*Mp2/m/L;c8=Ms2/m/L;c9=2*Is2/m/L/ax2;c10=(2*Mp2+Ms2)*g/m/L;l=zeros(20,1);%第i个轮对与以一个轮对之间的距离li-m

33、l(2)=ax1;l(3)=ax1+dax1;forj=4:18ifrem(j,2)=0%推断奇偶l(j)=l(j-1)+ax2;elsel(j)=l(j-1)+dax2;endendl(19)=l(18)+dax1;l(20)=l(19)+ax1;k0=zeros(1,n);nt=5;f0=sqrt(c1);%基频Mk=mu*m*L;Kk=m*L*mu*c1/(1+mu)2;Ck=2*Mk*sqrt(3*mu/8/(1+mu)3)*f0;m10=ones(1,n);m11=diag(m10);m21=zeros(41,n);m20=Mp1,Mp1,Ms1,Is1,.Mp2,Mp2,Ms2,I

34、s2,.Mp2,Mp2,Ms2,Is2,.Mp2,Mp2,Ms2,Is2,.Mp2,Mp2,Ms2,Is2,.Mp2,Mp2,Ms2,Is2,.Mp2,Mp2,Ms2,Is2,.Mp2,Mp2,Ms2,Is2,.Mp2,Mp2,Ms2,Is2,.Mp1,Mp1,Ms1,Is1,Mk;m22=diag(m20);c11=zeros(n,n);c12=zeros(n,41);c221=Cp1+Cs1,0,-Cs1,0.5*ax1*Cs1;0,Cp1+Cs1,-Cs1,-0.5*ax1*Cs1;-Cs1,-Cs1,2*Cs1,0;0.5*ax1*Cs1,-0.5*ax1*Cs1,0,0.5*ax1*

35、ax1*Cs1;c222=Cp2+Cs2,0,-Cs2,0.5*ax2*Cs2;0,Cp2+Cs2,-Cs2,-0.5*ax2*Cs2;-Cs2,-Cs2,2*Cs2,0;0.5*ax2*Cs2,-0.5*ax2*Cs2,0,0.5*ax2*ax2*Cs2;c220=zeros(4,4);c40=zeros(4,1);c42=zeros(1,40);c22=c221,c220,c220,c220,c220,c220,c220,c220,c220,c220,c40;c220,c222,c220,c220,c220,c220,c220,c220,c220,c220,c40;c220,c220,c2

36、22,c220,c220,c220,c220,c220,c220,c220,c40;c220,c220,c220,c222,c220,c220,c220,c220,c220,c220,c40;c220,c220,c220,c220,c222,c220,c220,c220,c220,c220,c40;c220,c220,c220,c220,c220,c222,c220,c220,c220,c220,c40;c220,c220,c220,c220,c220,c220,c222,c220,c220,c220,c40;c220,c220,c220,c220,c220,c220,c220,c222,c2

37、20,c220,c40;c220,c220,c220,c220,c220,c220,c220,c220,c222,c220,c40;c220,c220,c220,c220,c220,c220,c220,c220,c220,c221,c40;c42,Ck;k11=zeros(n,n);ii=1:n;Pkk=sin(0.5*ii*pi);k11=(ii.4*c1);k11=diag(k11);k12=zeros(n,41);k221=Kp1+Ks1,0,-Ks1,0.5*ax1*Ks1;0,Kp1+Ks1,-Ks1,-0.5*ax1*Ks1;-Ks1,-Ks1,2*Ks1,0;0.5*ax1*Ks

38、1,-0.5*ax1*Ks1,0,0.5*ax1*ax1*Ks1;k222=Kp2+Ks2,0,-Ks2,0.5*ax2*Ks2;0,Kp2+Ks2,-Ks2,-0.5*ax2*Ks2;-Ks2,-Ks2,2*Ks2,0;0.5*ax2*Ks2,-0.5*ax2*Ks2,0,0.5*ax2*ax2*Ks2;k220=zeros(4,4);k40=zeros(4,1);k42=zeros(1,40);k22=k221,k220,k220,k220,k220,k220,k220,k220,k220,k220,k40;k220,k222,k220,k220,k220,k220,k220,k220,k

39、220,k220,k40;k220,k220,k222,k220,k220,k220,k220,k220,k220,k220,k40;k220,k220,k220,k222,k220,k220,k220,k220,k220,k220,k40;k220,k220,k220,k220,k222,k220,k220,k220,k220,k220,k40;k220,k220,k220,k220,k220,k222,k220,k220,k220,k220,k40;k220,k220,k220,k220,k220,k220,k222,k220,k220,k220,k40;k220,k220,k220,k2

40、20,k220,k220,k220,k222,k220,k220,k40;k220,k220,k220,k220,k220,k220,k220,k220,k222,k220,k40;k220,k220,k220,k220,k220,k220,k220,k220,k220,k221,k40;k42,Kk;Y0=zeros(nn,1);%位移初始条件Y10=zeros(nn,1);%速度初始条件Y20=zeros(nn,1);%加速度初始条件X_mid=sin(pi/2),sin(pi),sin(3*pi/2),sin(2*pi),sin(5*pi/2);%前5阶振型在跨中取值ift=0retur

41、nelsedt=0.005;%积分步长nm=t/dt;Gama=0.50;Beta=0.25;a0=1/Beta/dt2;a1=Gama/Beta/dt;a2=1/Beta/dt;a3=0.5/Beta-1;a4=Gama/Beta-1;a5=dt/2*(Gama/Beta-2);a6=(1-Gama)*dt;a7=Gama*dt;w_mid=zeros(nm+1,1);%跨中挠度diracl=zeros(1,20);fork=0:nmt=k*dt;form=1:20%diraclidiracl(m)=diraci(l(m)/v,(L+l(m)/v,t);endii=1:n;mt=2*Mk/m

42、/L*sin(ii*pi/2);ctt=-Ck*sin(ii*pi/2);ktt=-Kk*sin(ii*pi/2);P011=sin(ii*v*t*c2);P022=sin(ii*(v*t-ax1)*c2);P033=sin(ii*(v*t-ax1-dax1)*c2);P044=sin(ii*(v*t-ax1-dax1-ax2)*c2);P055=sin(ii*(v*t-ax1-dax1-ax2-dax2)*c2);P066=sin(ii*(v*t-ax1-dax1-ax2-dax2-ax2)*c2);P077=sin(ii*(v*t-ax1-dax1-ax2-dax2-ax2-dax2)*c

43、2);P088=sin(ii*(v*t-ax1-dax1-ax2-dax2-ax2-dax2-ax2)*c2);P099=sin(ii*(v*t-ax1-dax1-ax2-dax2-ax2-dax2-ax2-dax2)*c2);P01010=sin(ii*(v*t-ax1-dax1-ax2-dax2-ax2-dax2-ax2-dax2-ax2)*c2);P01111=sin(ii*(v*t-ax1-dax1-4*ax2-4*dax2)*c2);P01212=sin(ii*(v*t-ax1-dax1-5*ax2-4*dax2)*c2);P01313=sin(ii*(v*t-ax1-dax1-5*ax2-5*dax2)*c2);P01414=sin(ii*(v*t-ax1-dax1-6*ax2-5*ax2)*c2);P01515=sin(ii*(v*t-ax1-dax1-6*ax2-6*dax2)*c2);P01616=sin(ii*(v*t-ax1-dax1-7*ax2-6*dax2)*c2);P01717=sin(ii*(v*t-ax1-dax1-7*ax2-7*dax2)*c2);P01818=sin(ii*(v*t-ax1-dax1-8*ax2-7*dax2)*c2);P01