(结构动力学)对简谐和周期荷载的反应

1、结构动力学（2010）结构动力学 第四章 单自由度体系对 简谐和周期荷载的反应单自由度体系对简谐荷载作用下的反应是结构动力学中的一个经典内容。不仅工程中实际存在这种形式的荷载，而且简谐荷载作用下单自由度体系的解提供了了解结构动力特性和用于分析更复杂荷载作用反应的手段和方法。4.1 无阻尼体系的简谐振动运动方程： 其中：p0 简谐荷载的幅值； 简谐荷载的圆频率。 初始条件 ： tpkuumsin0 )0(,)0(00uuuutt4.1 无阻尼体系的简谐振动运动方程是带有初值条件的二阶常微分方程，全解=齐次方程的通解+特解通解对应的方程是一个自由振动方程，其解uc为无阻尼自由振动： c - com

2、plementary tpkuumsin0 mktBtAtunnnc/sincos)(4.1 无阻尼体系的简谐振动特解满足运动方程的解，记为up(t) ，是由动荷载p0sint直接引起的振动解。设特解为： 其中，/n频率比，外荷载的激振频率与结构自振频率之比 ； pparticular tpkuumsin0 tDtCtupcossin)(0,)/(1120DkpCntpkuumsin0 4.1 无阻尼体系的简谐振动全解=通解+特解 待定系数A、B由初值条件确定 tkptBtAtututunnnpcsin)/(11sincos)()()(2020)/(1/)0()0(nnnkpuBuA)0()0

3、(00uuuutt4.1 无阻尼体系的简谐振动满足初始条件的解 ： 瞬态反应和稳态反应 tkptkpututunnnnnnsin)/(11sin)/(1/)0(cos)0()(2020瞬态反应 稳态反应 4.1 无阻尼体系的简谐振动稳态反应 ： u0稳态反应的振幅： ust等效静位移，或静位移： Rd动力放大系数： tkptunsin)/(11)(20kpust0200)/(11nkpu20)/(11nstduuR4.1 无阻尼体系的简谐振动无阻尼体系动力放大系数 =0 ，Rd =1 =n，Rd 发生共振 /n2， Rd1 20)/(11nstduuR4.1 无阻尼体系的简谐振动无阻尼体系共振

4、时动力反应时程 共振时(=n)： ttutunnstpcos2)(4.2 有阻尼体系的简谐振动运动方程： 初始条件：利用c=2mn，并将运动方程两边同除m，得到如下形式的运动方程： tpkuucumsin0 )0(,)0(00uuuutttmpuuunnsin202 4.2 有阻尼体系的简谐振动通解uc对应于有阻尼自由振动反应： 特解up可以设为如下形式 : tmpuuunnsin202 )sincos()(tBtAetuDDtcntDtCtupcossin)(0cos)(2sin2)(22022tDCtmpDCnnnn4.2 有阻尼体系的简谐振动 运动方程的全解：u(t)=uc+up ： 2

5、222222)/(2)/(1 /2)/(2)/(1 )/(1nnnstnnnstuDuC0)(1)2()2()(122DCuDCnnstnntDtCtBtAetuDDtncossin)sincos()(4.2 有阻尼体系的简谐振动 图4.3 有初始条件影响的动力反应时程 4.2 有阻尼体系的简谐振动（1）共振反应（=n） 满足零初始条件 运动解：当=0时 ： 与无阻尼时的结果完全相同 2222222)/(2)/(1 /2)/(2)/(1 )/(1nnnstnnnstuDuC2,0stuDCststuBuA2121,21ttteutunDDtstncos)sin1(cos2)(2ttutunns

6、tcos2)(tutBtAetustDDtncos2)sincos()(（1）有阻尼体系的共振反应（=n） 图4.4 有阻尼体系共振反应时程 4.2 有阻尼体系的简谐振动（2）动力放大系数Rd（dynamic magnification factor） 振动的稳态解： u0 稳态振动的振幅 相角，反映体系振动位移与简谐荷载的相位关系 )sin(cossin)(0tutDtCtu)(tan,1220CDDCu212220)/(1)/(2tan)/(2)/(1 1nnnnstuu2222222)/(2 )/(1 /2)/(2 )/(1 )/(1nnnstnnnstuDuC静位移kpust0（2）动

7、力放大系数Rd（dynamic magnification factor） 动力放大系数定义为 ： 2220)/(2)/(1 1nnstuustduuR0222)/(2)/(1 1nndR（2）动力放大系数Rd222)/(2)/(1 1nndR（2）动力放大系数Rd(1) 当21时，1dR，即体系不发生放大反应。 (2) 当21时，22max21)(,121)(峰值ndR。 (3) 当1/n（共振时） ，21dR。 (4) 当2/n时，1dR，对任意 均成立。 4.2 有阻尼体系的简谐振动（3）阻尼体系动力反应与荷载的相位关系 在动力荷载的作用下，有阻尼体系的动力反应（位移）一定要滞后动力荷载

8、一段时间，即存在反应滞后现象。这个滞后的时间即由相角反映，如果滞后时间为t0，则= t0 （t0=/）。由计算的公式可知，滞后的相角与频率比/n和阻尼大小均有关系。 21)/(1)/(2tannn（3）阻尼体系动力反应与荷载的相位关系 右图给出阻尼比=0.2时，相应于不同频率比/n时的外力和位移曲线及滞后相角。相角实际是反映结构体系位移（反应）相应于动力荷载的反应滞后时间，从图中可以发现，频率比越大，即外荷载作用得越快，动力反应的滞后时间越长。 （3）阻尼体系动力反应与荷载的相位关系 （3）阻尼体系动力反应与荷载的相位关系 表 4.1 三种特殊情况时体系振动位移与简谐荷载的相位关系 由n/图判

9、断 物理解释 （根据关系：uufufufICS2, ） 0/n时 0 0则u 和0u ，即cf和0If 则)(tpfS即)(tpku ，u与)(tp相位相同 1/n时 90 90)()(90(),()(0,2相差，则相差与）相同，而与即则则tptuuutputpuctpffffkuumumfcSIsnI n/时 180 180180)()(,，则相差位移反相，所以位移与，而惯性力与则和，则tptpffffIcSI 4.3 振动测量仪器（拾振仪） 测量振动量仪器主要有三种： 加速度计：测量加速时程（强震仪） 位移计：测量位移时程（地震仪） 速度计：测量速度（目前应用逐步多起来） 4.3 振动测量

10、仪器（拾振仪） （1）加速度计（强震仪） 加速度计测量的是加速度在基底加速度作用下仪器质点的运动方程为：仪器基底加速度时程：仪器质点所记录的相对位移u(t)为： )(tumkuucumg tutuggsin)(0 )sin()sin()/(2)/(1 1)(02220tuRkmtkumtugdnng （1）加速度计（强震仪） 为简单起见，仅讨论u(t)的振幅u0 ：通常采用提高加速度计中弹簧刚度的方法来实现提高n的目的。因此，加速度计或强震仪中弹簧刚度比较大，是比较刚性的。 00)(gduRkmu 5 . 0/0nmkn/n5 . 07 . 04.3 振动测量仪器（拾振仪） （2）位移计（地震

11、仪）位移计是用来测量仪器基底的位移量 仪器基底位移时程：在基底位移作用下仪器质点的运动方程为：仪器质点的所记录的相对位移u(t)为： )sin(02tumkuucumg tuuggsin0)sin()()sin()(0202tuRtRukmtugdndg（1）位移计（地震仪） 通常采用降低位移计中弹簧刚度的方法来实现降低n的目的。因此，位移计中弹簧刚度比较小，是比较柔性的。 1/nn5 . 0020)(gdnuRudngRuu200)(位移比4.4 隔振（震） 隔振（震）分两种情况：1）阻止振动的输出。例如，大型机器动力机器振动向地基中的传播；地铁车辆振动传播。力的传递和隔震 2）阻止振动的输

12、入。例如，结构抗震问题中的隔震设计，在振动的结构或地基上安装的精密仪器设备的隔震问题。基底振动的隔离 4.4 隔振（震） 力的传递和隔震 基底振动的隔离 4.4 隔振（震） 1、力的传递和隔震p0 sint机器的不平衡力 机器的转速(角速度)m 机器质量(设为刚性质量块) k、c隔振元件的总刚度和阻尼。fT 从隔振元件传到地基上的力 单质点体系简谐振动问题的解为: uckufffDsT)sin()(tRutudst1、力的传递和隔震传到地基上的力为：作用力fT的最大值为：将ust=p0 /k、c=2mn代入上式得：TR传递率(transmissibility)，是反映隔振效果的量 uckuff

13、fDsT)sin()(tRutudst)cos()sin()(tctkRutfdstT222maxckRufdstT22220max)/(2)/(1 )/(21nnnTpfTR1、力的传递和隔震力的传递率TR 当频率比：传递率：为达到隔振的目的，可采用降低n的办法，减小隔振元件刚度， 或增加仪器质量的方法，提高隔振效果。实际的减震设计方案应在尽量小的刚度和可接受的静位移之间优化选取。阻尼对隔振的影响？ 22220max)/(2)/(1 )/(21nnnTpfTR2n1TR4.4 隔振（震） 2、基底振动的隔离ug(t)基底（地面）的振动位移时程；ut(t)=u(t)+ug(t)质点的绝对运动时

14、程；u(t)相对位移。输入基底运动的位移时程为：单质点体系简谐振动问题的解为: 质点的总位移ut(t)为： tutuggsin)(0)sin()()(02tuRtugdn)sin()/(21)()()(120tRutututundggt2、基底振动的隔离位移的传递率TR为 ：位移的传递率与力的传递率完全相同，说明两种隔振问题是相通和相同的，其隔振设计方法也基本相同 )sin()/(21)()()(120tRutututundggt200)/(21ndgtRuuTR222200)/(2)/(1 )/(21nnngtuuTR222)/(2)/(1 1nndR2、基底振动的隔离对建筑结构的隔震问题与

15、以上讨论的单质点体系隔振问题有类似的地方。例如都是试图通过降低体系自振频率的方法来提高隔震(振)效率。也有不同的地方，建筑结构体系是多自由度体系，其隔震(振)效率的研究更复杂，而且地震动是宽频带的过程，总有与结构自振频率相同的频率成份存在，无法通过避开地震动频率的方法来实现隔震(振)的目地。隔震问题研究已成为一门专门的工程抗震研究领域。加速度传递率为 ：加速度的传递率与位移的传递率相同 00020200gtgtgtuuuuuuTR 2、基底振动的隔离算例1，工程场地竖向加速度为g=0.1g，振动频率为f=10Hz，安放一个重m=50kg的敏感仪器，仪器固定在刚度k=14kN/m，阻尼比=10%

16、的橡胶隔振垫上，问：传递到仪器上的加速度是多少？如果仪器只能承受0.005g的加速度，给出解决方案。2、基底振动的隔离解： 传递到仪器上的加速度是多少？求TR体系自振频率： 频率比：加速度传递率：传递到仪器上的加速度：g=0.1g，f=10Hz，m=50kg，k=14kN/m，=10%sec/73.1610005014radmkn75. 373.16102ngguTRugt009. 01 . 0091. 000 091. 0)(2)(1 )(21222200nnngtuuTR 2、基底振动的隔离解：如果仪器只能承受0.005g的 加速度，给出解决方案。降低体系的自振频率n，即增大/n可以提高隔

17、振效率，由于隔振垫参数不易改变，采用增加附加质量办法降低n，先假设附加质量mb=60kg，体系总质量m=m+mb=110kg。 例题1体系新的n和/n可为：体系新的阻尼比： 57. 528.11102sec/28.11100011014nnrad067. 01 . 028.11)5060(73.165022nnnnmmmmc例题1体系新的传递率：传递到仪器上的加速度：因为方案成功。 如果要求附加质量后，0t=0.005g，则mb应是多少？57. 528.11102n067. 004. 0)57. 5067. 02()57. 51 ()57. 5067. 02(12222RTgguRTugt00

18、4. 01 . 004. 000 ggut005. 0004. 00 算例2 汽车在多跨连续梁上行驶，桥梁跨度均为L=30m，桥面由于长时徐(蠕)变效应而产生15cm的挠度(桥面的中点)。桥面可以用振幅为7.5cm的正弦曲线来近似，汽车可以用一个单质点SDOF体系模拟，如果车重m=1.8t，等效弹簧刚度K=140kN/m，等效阻尼比=40%， 求： 车以80km/h行驶时，汽车的竖向运动ut(t)的振幅u0t 发生共振时汽车的行驶速度(使振幅最大时的速度) 算例2 解： 车以80km/h行驶时，汽车的竖向运动ut(t)的振幅u0t汽车相当于受振幅为ug0=0.075m， 波长为L=30m的简谐

19、运动ug的干扰简谐运动的周期:车辆的固有周期: ssmmvLT35. 1)3600/80000(30/skmTnn71. 01408 . 1222算例2 频率比:振动传递率:汽车竖向运动的振幅： sT35. 1sTn71. 053. 035. 171. 0TTnn3 . 1)53. 04 . 02()53. 01 ()53. 04 . 02(1)2()1 ()2(122222222TRmuuTRuggt0975. 0075. 03 . 13 . 1000算例2 解： 发生共振时汽车的行驶速度(使振幅最大时的速度)如果体系的阻尼比很小，当=n时ut最大，而本问题阻尼比0.4很大，因此使u0t取最

20、大的不一定等于n，此时要采用取极值条件求使u0t最大，即使TR取最大。算例2 使TR取最大值的频率，也使TR2取最大值。当汽车的行驶速度为135km/h时，车辆的振幅达到最大值 vLT/22222)2()1 ()2(1(TR02TR89. 0TTnn798. 089. 071. 0nTTsTn71. 0hkmsmTLv/135/6 .37798. 0304.5 用简谐振动(强迫振动)试验 确定体系的粘性阻尼比 可以用自由振动方法求阻尼比的原因是由于自振衰减的快慢由控制，或说衰减规律可以明显反应出阻尼比的影响。而动力放大系数同样受控制，Rd曲线形状可以反映出的影响，其影响主要有两点： (1)峰值

21、大小， (2)曲线的胖瘦。利用体系对简谐荷载反应的结果也可以得到体系的阻尼比，有两种主要方法：共振放大法和半功率(带宽)法，其原理均是基于对动力放大系数Rd的分析。4.5 用强迫振动试验确定体系的阻尼比 1、共振放大法根据动力放大系数Rd : 当发生共振(/n1)时: 222)/(2)/(1 1nndR21)()(00stnstnduuuuRn)(2)(210nstnduuR4.5 用强迫振动试验确定体系的阻尼比 1、共振放大法由于从动力放大曲线定u0(n)不容易，一般用u0m代替，u0m=max(u0)，则: 用共振放大法确定体系的阻尼比，方法简单。但实际工程中测得的动力放大系数曲线一般以u

22、0图给出，用以上公式计算阻尼比时，还需得到零频时的静位移值ust，实际测量静载位移无论从加载设备和记录(拾振)设备都有一定的困难，即实现动力加荷和测量动力信号的设备不能在零频率时工作。因此工程中往往采用半功率(带宽)法从动力试验中得到阻尼比。 )(2)(210nstnduuRmstduuR0max2)(214.5 用强迫振动试验确定体系的阻尼比 2、半功率带宽法 (半功率点法) 半功率点：动力放大系数Rd上振幅值等于1/2倍最大振幅的点所对应的两个频率点。 记：a和b分别等于半功 率点对应的两个频率。 则阻尼比 可由如下公式计算： nab2ababababffffnabfff2半功率带宽法 (

23、半功率点法) 证明： (4.34) (4.35)nab2 由 Rd可知，Rd的最大值。2max121)(dR。而振幅等于21倍(Rd)max对应的频率满足以下方程： 222212121)/(2)/(1 1nn (a) 对式(a)两边同时取倒数、并开平方，整理后得： 0)1 (81)(21 (2)(22224nn (b) 式(b)是关于(/n)2一元二次方程，可得两个根为： 22212)21 ()(n (c) 式(c)取正号时对应数值较大的根b，负号对应较小的根a。一般的工程结构，阻尼比较小，式(c)中 的平方项可忽略，因此 121n 则对应于半功率点的两个根为： 1,1nanb (d) 由式(

24、d)得到半功率点频率b和a与阻尼比 的关系， 2nab (e) 由此得到式（4.34） 。若再用式(d)得关系2nab，代入式(e)，又得到式（4.35） 。 abab三种阻尼比的测量方法 共介绍了三种测量结构阻尼的方法：对数衰减率法、共振放大法和半功率带宽法，虽然是针对单自由度体系推导的，但这些方法对多自由度体系同样适用。下面对这三种方法给一简要的总结。(1)对数衰减率法 采用自由振动试验，测一阶振型的阻尼比较容易。高阶振型的阻尼比的关键是能激发出按相应振型进行的自由振动。(2)共振放大法 采用强迫振动试验，由于静（零频）荷载下的位移较难确定，应用上存在一定的技术困难，但通过一定数学上的处理

25、还是可用的，例如，利用接近零频的非零频位移通过插值外推得到零频时的位移值。(3)半功率带宽法 采用强迫振动试验，不但能用于单自由度也可用于多自由度体系，对多自由度体系要求共振频率稀疏，即多个自振频率应相隔较远，保证在确定相应于某一自振频率的半功率点时不受相邻频率的影响。4.6 粘性阻尼的能量耗散和等效粘性阻尼 1、粘性阻尼体系的能量耗散SDOF体系在简谐力p(t)=p0sint作用下， 在一个振动循环内的能量耗散记为：ED 阻尼引起的能量耗散，即阻尼力做的功；EI 外力做的功；ES 弹性力做的功；EK 惯性力做的功。在简谐荷载p(t)作用下， SDOF的位移为： )sin()(0tutu1、粘

26、性阻尼体系的能量耗散（1）阻尼引起的能量耗散ED粘性阻尼引起的耗散与振幅u0的平方成正比， 与阻尼比和外荷载的频率成正比。 )sin()(0tutu2020/2020/202/20)(2)cos()(kuucdttucdtucdtuucdufEnDDnnnnkkkmc/2)/(2221、粘性阻尼体系的能量耗散（2）外力做的功EI （IInput) )sin()(0tutu2020)(2kuucEnD2000/2000/20)(2sin)cos()sin()()(kuupdttutpdtutpdutpEnIkpuRndn/)(2)2(sin001、粘性阻尼体系的能量耗散（3）弹性力的功 ES （

27、4）惯性力的功 EK (Kinetic) 可见在简谐振动中的一个循环内，弹性力和惯性力做功均等于零，而由阻尼耗散的能量等于外力做的功。 )sin()(0tutu0)cos()sin()(0/200/20dttutukdtukudufEsS/20002/200)cos()sin()(dttutumdtuumdufEIK 4.6 粘性阻尼的能量耗散和等效粘性阻尼 2、等效粘性阻尼(1) 粘性阻尼是一种理想化的阻尼，具有简单和便于分析计算的优点。(2) 工程中结构的阻尼源于多方面，其特点和数学描述更为复杂，这时可以将复杂的阻尼在一定的意义上等效成粘性阻尼。(3) 一般采用基于能量等效的原则。(4)

28、阻尼耗散能量的大小可以用阻尼力的滞回曲线反映。2、等效粘性阻尼（1）阻尼力的滞回曲线阻尼力的滞回曲线：阻尼力与位移之间的关系曲线，即fDu曲线。 )sin()(0tutu)()sin()cos()(22020200tuuctuuctuctucfD粘性阻尼力滞回曲线2、等效粘性阻尼（1）阻尼力的滞回曲线对粘性阻尼力的滞回曲线整理可以得到：研究滞回曲线的意义：力在一个循环内所做的功等于证明： 滞回曲线所包围的面积。 )(220tuucfD粘性阻尼力的滞回曲线是一椭圆 1)()(2020ucfuuDDnDEkuucucuabS2020002)(椭圆面积：2、等效粘性阻尼（1）阻尼力的滞回曲线抗力曲线

29、：fD+ fsu曲线。fD+ fs有时称为抗力。抗力滞回曲线包围的面积等于阻尼力做的功。 在实际测量时，量测到的量是抗力。2、等效粘性阻尼（2）等效粘性阻尼比确定等效粘性阻尼比的原则：基于能量耗散相等的原理。具体实现方法：在一个振动循环内让等效粘性阻尼做的功等于实际阻尼所做的功。（2）等效粘性阻尼比在一个循环内实际阻尼力作的功： 在一个循环内等效阻尼力作的功： 粘DEDE（2）等效粘性阻尼比202kuEneqD粘DDEE粘20)/(2kuEnDeq202 kuEDeq04SDeqEE20021kuES1/n4.7 滞变阻尼（复阻尼）理论 粘滞阻尼由于其在建立运功方程和求解时的方便性，而在工程中

30、得到广泛应用。但它也存在一个严重的缺陷，即，粘滞阻尼力和能量耗散与激振频率有关。例如在每一振动循环中耗散的能量为：对一结构体系，阻尼比为常数，固定振幅u0，则在每一振动循环中耗散的能量与激振频率成正比，这与结构试验结果不符，试验结果表明，阻尼力或其耗能与频率基本是无关的。为此，人们发展了滞变阻尼理论（hysteretic）。 202kuEnD4.7 滞变阻尼（复阻尼）理论 滞变阻尼（hysteretic）：阻尼力大小与位移幅值成正比 而与速度同相。三种型式的滞变阻尼定义： 克拉夫1981 Clough(克拉夫) 1993 Chopra 1995 其中为滞变阻尼参数。 第一种型式是直接套用滞变阻

31、尼的定义； 第二种是滞变阻尼的复数形式； 第三种是从构造频率无关阻尼的构思出发。 )()()(tututukfD)(tkuifD)(tukfD4.7 滞变阻尼（复阻尼）理论 三种型式定义的滞变阻尼在复数域是完全等价， 例如，假设 u(t)=u0eit，则 fD 均可写成 iKu(t) ， 但在实数域则不尽相同。 共同点是耗能与频率无关，但具体耗能值不同。 第一种： 第三种： 从滞回曲线形状分析，第一种形式滞变阻尼与实际相差太大，不可接受。 第二和第三种形式滞变阻尼的耗能相同。2012 kuED2032kuEEDD滞变阻尼与粘滞阻尼的关系将滞变阻尼滞回耗能关系代入计算等效粘性阻尼比公式 得到：

32、或:当共振时，滞变阻尼参数与阻尼比关系： )/(21n20)/(2kuEnDeq203kuEEDD)(2n24.7 滞变阻尼（复阻尼）理论 滞变阻尼参数与阻尼比关系式：2是在n时取得的，对n时并不成立。有些教科书中没有明确指出这点，有时导致模糊的概念。 图4.15 粘性阻尼与滞变阻尼耗能ED与激振频率的关系4.7 滞变阻尼（复阻尼）理论 用滞变阻尼第二种表达式 fDiKu(t) 时也称为复阻尼，在复阻尼理论中，将阻尼力和弹性恢复力合在一起构成复刚度。由 fDiku(t) 和 fs=ku(t)，可定义复刚度为：在复数形式的简谐荷载作用下质点的运动方程可写成：kik 1tiepukum0 4.7

33、滞变阻尼（复阻尼）理论 复数形式的稳态解可设成： tiUetu)(tiepukum0 tiUeu2 titiepUeKm022220)(1 )(1 nniKPUkik 14.7 滞变阻尼（复阻尼）理论 复数形式的稳态反应为： u(t)是一个复函数，可写成它的模与单位复数积的形式其中，tinneikptu2220)(1 )(1 )()(0)(tieutu2122200)(1tan)(1 1nnkpu4.7 滞变阻尼（复阻尼）理论 当取 时，复阻尼理论的解与粘滞阻尼理论的解完全相同。 2122200)(1tan)(1 1nnkpu)(2n2220)/(2)/(11nnstuu21)/(1)/(2t

34、annn4.7 滞变阻尼（复阻尼）理论 滞变阻尼和粘性阻尼的耗能滞变阻尼的耗能接近实际，而粘性阻尼当外力频率较低时，低估了体系的耗能能力；外力频率较高时，又会过高估计耗能能力。因此，希望通过阻尼比的选取使粘性阻尼理论能正确反映所有频率时体系的耗能是不可能的，一个较为稳妥的方法是使阻尼比的选取能较为正确地反映感兴趣频段内的耗能能力。这可通过设外荷载频率等于感兴趣频率的方法实现。实际的做法是取外荷载频率等于结构自振频率，此时结构的反应最大，是阻尼影响最大的点。由于结构往复试验时，在不同频率下得到的滞回环面积基本相等，因此可以用共振时的公式来定阻尼比，而不考虑实际加荷频率，这样得到的阻尼比对反映共振时的耗能能力相对准确。4.8 单自由度体系对周期荷载的反应 依靠的基础： 依靠已得到的单自由度体系对简谐荷载反应分析结果。在获得简谐荷载作用的结果后，就可以方便地分析单自由度体系对任意周期性荷载的反应，简谐荷载是一种最简单、最具代表性的周期荷载，而任意周期性荷载均可以分解成简谐荷载的代数和。具体实施方法： 利用Fourier级数展开法。将任意的周期荷载p(t)展开成Fourier级数，把任意周期性荷载表示成一系列简谐荷载的叠加，对每一简谐荷载作用下结构的反应可以容易得到其稳态解，再求和，得到结构在任意周期性荷载作用下的反应。限制