第6章叠加法对一般动力荷载的反应-1

高等结构动力学第六章高等结构动力学对一般动力荷载的反应---叠加法高等结构动力学§6.1无阻尼体系的Duhamel积分§6.2无阻尼体系Duhamel积分的数值计算§6.3有阻尼体系的反应第六章对一般动力荷载的反应—叠加法高等结构动力学图6-1（无阻尼）Duhamel积分的推导结构承受图6-1所示任意的一般性荷裁p(t)，求解体系的反应。§6.1无阻尼体系的Duhamel积分§6.1无阻尼体系的Duhamel积分高等结构动力学求解的思路：第六章中所述的短持续时间冲击作用下计算结构反应的近似分析方法，可以用作推导在一般动力荷载下计算反应的基础；t=τ时的荷载强度P（τ），在dτ时间间隔内作用结构上，将产生荷载冲量p(τ)dτ。用方程(6-19)来计算结构对这个冲量的反应。无阻尼体系的Duhamel积分§6.1无阻尼体系的Duhamel积分高等结构动力学对于有限持续时间的冲量来说，这种方法是近似的，但当荷载的持续时间趋于零时，它却是精确的;因此，在微分意义下dτ期间，荷载p(τ)所产生的反应恰为(对于tτ)()()sin()pddvttm（6-1）dv(t)项表示在tτ时程范围内微分冲击p(τ)dτ引起的自由振动反应的贡献，不是时间间隔dt内v的变化。§6.1无阻尼体系的Duhamel积分高等结构动力学图6-1（无阻尼）Duhamel积分的推导整个荷载时程可以视作由一系列连续的短脉冲所组成，每一个脉冲将产生一个如方程（6-1）所示的微分反应。对于这个线性弹性体来说，则总反应可将在荷载时程所产生的全部微分反应相加而获得，亦即对方程（6-1）进行积分：tdtpmtv0)(sin)(1)(（6-2）§6.1无阻尼体系的Duhamel积分高等结构动力学对问题的再认识数学上的积分问题；力学上，对冲量引起的自由振动影响的迭加原理应用。这样，Duhamel积分的应用条件是什么呢？§6.1无阻尼体系的Duhamel积分高等结构动力学方程（6-2）一般称作无阻尼体系的Duhamel积分。它可用来计算任意形式的动力荷载p(t)作用下无阻尼单自由度体系的反应。但在荷载变化很不规则时，计算必须利用数值积分来进行。（6-3）方程（6-2）也可写成如下形式tdthptv0)()()(式中新的符号的定义如下)(sin1)(tmth（6-4）方程（6-3）称作卷积积分。计算任意荷裁作用下结构的反应时，利用这个积分可以获得整个时间域的反应。函数h(t-τ)一般称作单位脉冲反应（这是对无阻尼体系定义的），因为它表示，在t=τ时，在一个单位大小的脉冲作用下，结构的反应。§6.1无阻尼体系的Duhamel积分高等结构动力学不用说，在方程（6-2）中，必须假定在荷载开始作用时（t=0时），结构处于静止状态。对于其它特定的初始条件：v(0)≠0和(0)≠0，这个解还必须加上一个附加自由振动反应。因此，一般情况下（6-5）tdtpmtvtvtv0)(sin)(1cos)0(sin)0()(§6.1无阻尼体系的Duhamel积分高等结构动力学如果所施加的荷载函数是可积的，结构的动力反应可以用前面(6-2)或(6-5)式的积分来计算。在许多实际情况中，因为荷载仅由试验数据给出，此时反应就必须用数值的方法才能十算。注意到在这样的分析中，三角恒等式，和把方程(6-2)写成如下形式（在零初始假定条件下）将是有用的：（6-10）或ttdpmtdpmttv00sin)(1coscos)(1sin)(ttBttAtvcos)(sin)()(其中ttdpmtBdpmtA00sin)(1)(cos)(1)(（6-11）§6.2无阻尼体系Duhamel积分的数值计算§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学因此，要进行Duhamel积分的数值积分，必须进行积分和的数值计算。例如，苜先考察这两个积分中的第一个，其被积函数绘于图7-2中。为了数值计算的方便，计算等时间增量△τ的函数值，函数的相继值用相应的角标来区分。可将这些纵坐标乘以适当的加权系数然后相加而得到近似的积分值。以数学式子表达，即（6-8*）其中，而表示数值求和过程，其具体形式依赖于所利用的近以积分的序号。对于三种基本的近似方法，其求和方法如下进行：AttmdymtA)(1)(1)(0简单求和法(ζ=1):12101)(NAyyyyt（6-9a*）§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学图6-2Duhamel积分数值求和法的公式化梯形法则(ζ=2):（6-9b*）NNAyyyyyt12102222)(§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学simpson法则(ζ=3):（6-9c*）NNAyyyyyt12103424)(对于simpson法则来说，式中必须是偶数。利用式(6-9*)中任何一个式子和式(6-8*)，可得到所考虑的特定时间t的积分的近似值。但是，一般来说，需要求整个反应过程而不只是求某特定时间的位移，换句话说，必须计算一系列相继时刻t1，t2，…的反应，其中两相邻时间之间的间隔为（如用simpson法则寸为2）。为了得到全部的反应过程，把方程（6-9*）的求和表示成增量形式更为方便。/Nt§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学simpson法则(ζ=3):（6-10c*）]cos)()(cos)([)()(22ttpttpttAA33()(2)(2)cos(2)4()cos()cos()AAttpttpttptt简单求和法(ζ=1):)(cos)()()(11ttpttAA（6-10a*）梯形法则(ζ=2):（6-10b*）§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学简单求和0121()N=1,2,3,(6-12a)NNAyyyym梯形法则0121(222)N=1,2,3,(6-12b)2NNNAyyyyymSimpson法则0121(424)N=2,4,6,(6-12c)3NNNAyyyyym利用其中任意一个式子，都可以直接由多需N的任意具体值来获得。但是，通常是要求整个反应的时间历程，为此必须计算直到获得所期望反应时间历程的一系列N的值。为此，应用这些式子的如下递归形式更有效：NANA无阻尼体系§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学11()N=1,2,3,(6-13a)NNNAAym简单求和梯形法则11()N=1,2,3,(6-13b)2NNNNAAyymSimpson法则221(4)N=2,4,6,(6-13c)3NNNNNAAyyym这里00A§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学（6-11*）项可用完全同样的方法来计算，即BtmtB)(1)(其中可以用与方程（6-10*）一样的表达式来计算，但用正弦函数代替余弦函数，将方程（6-8*）和（6-11*）代入方程（6-10）则可导得无阻尼体系最后的反应方程：（6-12*）BAttttmtvcos)(sin)(1)(§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学例题E6-1为了说明Duhamel积分的数值计算，我们来计算一个承受冲击波的水塔的动力反应。结构和冲击波荷载的理想化模型示于图E6-1.图E6-1承受冲击波荷载的水塔§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学对于这个体系，振动频率和周期为：270032.33096.6kgW20.209TNNAyyyyyt12103424)(弧度/秒秒•在数值积分中所用的时间增量为秒，其相应的自由振动角增量为弧度（较长的增量大概也可获得同样满意的结果）。计算的公式为：•其中•在这个无阻尼分析中，应用Simpson法则求和，公式如下：0.005300.0050.153因此。（6-7*）计算过程如下：()vtttdpmtBdpmtA00sin)(1)(cos)(1)(（6-9*）()()sin()cosvtAttBtt（6-6*）§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学初值的选取：由图6-1可知，当时，，当时，根据三角形比例关系，依次计算出相应的，从而得出相应的和的值。02(0)0p（）0.005,0.010,......()19.32,38.64......psintcost0()os05c01pt当时，sint=0,cost=1(3)(4)（）乘数=1（6）（根据公式（6-9*）确定）乘数=4（6）（根据公式（6-9*）确定）()cos1(7)(010pt乘数）0.005()cos19.320.98919.095pt当时，p()=19.32,sint=0.150,cost=0.(3)(4),98（9）50.010()cos38.640.95536.9pt当时，p()=38.64,sint=0.295,cost=0.(3)(4),95（5）()cos4(19.09476.4(7)pt乘数）§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学乘数=1（6）（根据公式（6-9*）确定）(7)()cos1(36.9136.9pt乘数）（三项之和）076.436.9113.(8)3A0113.31133(9).A（两项之和）0()p当时，sint=0(10)乘数=1（11）（根据公式（6-9*）确定）0.005()0.152.9p当时，sint=19.32(10)()sin1((12010)pt乘数）乘数=4（11）（根据公式（6-9*）确定）(12)()sin4(2.9411.6pt乘数）0.010()p当时，sint=38.640.295=11.4(10)乘数=1（11）（根据公式（6-9*）确定）(12)()sin1(11.4111.4pt乘数）§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学011.611.423.0(13)B（三项之和）023.(14)023.0B（两项之和）0,05001sinAt当时（）cos01.000106Bt（）根据公式（6-10）sinco(7)s0010AtBt51.85210Gm()(sin(18)cos)0vGAtBt英尺K=2700千磅/英尺()()0sfkv千磅(19)以时为例同理进行计算：0.020(0.020)77.28p千磅sinsin(0.6)0.564tcoscos(0.6)0.826t§6.2无阻尼体系Duhamel积分的数值计算高等结构动力学()cos()(0.020)cos(0.6)77.280.82663.8ptp乘数=1(0.020)cos(0.6)163.8p()cos()4284.4pt千磅()cos()148.0pt千磅当时0.025当时0.03063.8284.448.0396.2A422.4396.2818.6A()sin()(0.020)sin(0.6)77.280.56443.6ptp乘数=1(0.020)