结构动力学大作业.docx

《结构动力学大作业.docx》由会员分享，可在线阅读，更多相关《结构动力学大作业.docx（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、结构动力学大作业目录1.力插值法(1)1.1分段常数插值法(1)1.2分段线性插值法(4)2.加速度插值法(7)2.1常加速度法(7)2.2线加速度法(9)附录(12)分段常数插值法源程序(12)分段线性插值法源程序(12)常加速度法源程序(13)线加速度法源程序(13)1.力插值法力插值法对构造的外荷载进行插值，分为分段常数插值法和分段线性插值法，这两种方法均适用于线性构造的动力反响计算。1.1分段常数插值法图1-1为一个单自由度无阻尼系统，构造的刚度为k，质量为m，位移为y(t)，施加的外力为P(t)。图1-2为矩形脉冲荷载的示意图，图中td表示作用的时间，P0表示脉冲荷载的大小。图1-1

2、单自由度无阻尼系统示意图图1-2矩形脉冲荷载示意图对于一个知足静止初始条件的无阻尼单自由度体系来讲，当施加一个td时间的矩形脉冲荷载，此时构造在td时间内的位移反响能够用杜哈梅积分得到：0()sin()2(1cos)(1cos)(0)tststdPyttdmtytyttT=-=-=-?(1-1)假如构造本身有初始的位移和速度，那么叠加上构造自由振动的部分，构造的位移反响为：02()cossin(1cos)(0)stdytytyttyttT=+-(1-2)图1-3分段常数插值法微段示意图对于施加于构造任意大小的力，将其划分为t的微段，每一段的荷载都为一个常数(每段相当于一个矩形的脉冲荷载)，如图

3、1-3所示，则将每一段的位移和速度写成增量的形式为：1costsint(1cost)iiiiyPyyk+=?+?+-?(1-3)i+1/sintcostsintiiiyPyyk=-?+?+?(1-4)程序流程图如下i+1costsint(1cost)iiiyPyyk=?+?+-?i+1/sintcostsintiiiyPyyk=-?+?+?图1-4分段常数插值法流程图根据流程图能够编写相应的算法，利用MATLAB进行编程，程序源代码见附录。为了验证程序的正确性，本文选取的下面的例题进行验证。对于一个单自由度的无阻尼构造，当其遭到一个周期荷载时，其构造响应分为稳态解和瞬态解，由于没有阻尼的影响，

4、其瞬态解并不会衰减，其理论表达式为：021()()(sinsin)1pxtttk=-(1-5)式中，()xt为位移响应，0p为鼓励，k为刚度，为荷载频率与固有振动频率之比，为荷载频率，为构造固有频率。现令0p为1，k为1，则为1，取为2/3。程序求得的解与解析解比照如图1-5所示由于理论解与程序基本重合，所以将理论解乘以-1，方便比拟：位移y时间ta)位移速度v时间tb)速度图1-5分段常数插值法结果验证由图1-5可知理论解与程序算得的解基本重合，能够验证程序的准确性。1.2分段线性插值法与分段常数插值法不同，分段线性插值法将每一微段的力当成一个线性的直线，对于每一个微段，可看成一个矩形和一个

5、三角形脉冲的叠加。图1-6为分段线性插值微段示意图。图1-6分段线性插值法微段示意图对于无阻尼的体系，后一个时间步的位移和速度可由前一个时间步相应的值求得：11cossin(1cos)(1sin)iiiiiyPPyyttttkkt+?=?+?+-?+-?(1-6)11/sincossin(1cos)iiiiiyPPyyttttkkt+?=-?+?+?+-?(1-7)分段线性插值法的流程图如图1-7所示，与分段常数插值法仅仅是迭代的方式有所不一样。11cossin(1cos)(1sin)iiiiiyPPyyttttkkt+?=?+?+-?+-?11/sincossin(1cos)iiiiiyPP

6、yyttttkkt+?=-?+?+?+-?图1-7分段线性插值法流程图程序源代码见附录，同样利用1.1节的算例进行验证，所得的结果如图1-8所示。2.加速度插值法加速度插值法也叫逐步积分法，其对加速度进行插值，可分为常加速度法和线加速度法。2.1常加速度法图2-1常加速度法微段示意图对于一个单自由度构造，其运动方程为：my()y()()()tctkytpt+=(2-1)将式1-1转变为增量方程：mycykyp?+?+?=?(2-2)在通过逐步积分，将时间转化为一系列微小的时间段t?，如图2-1所示，现令111()()2iiiiytyyttt+=+4)将式(2-3)、(2-4)代入式(2-1)，

7、即：1111myiiiicykyp+=(2-5)此时，式(2-5)中只要1iy+为未知变量，可直接求出1iy+，之后再利用式(2-3)、(2-4)，可求出ti+1时刻的速度与位移。算法的流程图如下所示：将荷载作用的时间划分为0243yt=?tt?*4(2)2iiiimPPcymyt?=?+?22iiiyyyt?=?-?24()2iiiiyyytyt?=?-?-?图2-2常加速度法流程图算例验证的结果如下列图所示，讲明了该程序的正确性。由于需要对加速度进行插值，此处增加了加速度验证。位移y时间ta)位移速度v时间tb)速度加速度a时间tc)加速度图2-3常加速法结果验证2.2线加速度法线加速度法与常加速度法原理类似，其速度与位移的增量方程与常加速度法对应的方程略有不同，图2-4为线加速度法微段示意图。a)位移b)速度图2-4线加速度法微段示意图关于详细原理不在赘述，下列图为线加速度法的流程图。