结构动力学大作业2.docx

结构动力学大作业班级： 学号：姓名： 目 录1.Wilson-θ法原理简介22.Wilson-θ程序验算32.1△t的影响42.2θ的影响53.非线性问题求解54.附录8Wilson-θ法源程序8Wilson-θ法原理简介图1-1 Wilson-θ法示意图Wilson-θ法是基于对加速度a的插值近似得到的，图1-1为Wilson-θ法的原理示意图。推导由t时刻的状态求t+△t时刻的状态的递推公式： (1-1)对τ积分可得速度与位移的表达式如下： (1-2) (1-3)其中τ=θt，由式(1-2)、(1-3)可以解出： (1-4) (1-5)将式(1-4)、(1-5)带入运动方程： (1-6) (1-7)注意到此时的式子为{}和上一个时刻、、以及t+θ△t时刻的荷载相关，可以运用迭代的思想来求解，下图给出线弹性条件下Wilson-θ法的流程图：图1-2 Wilson-θ法流程图Wilson-θ程序验算对线弹性条件下的Wilson-θ法进行MATLAB编程，源代码见附录。选取如下算例进行验证。对于一个单自由度的无阻尼结构，当其受到一个周期荷载时，其结构响应分为稳态解和瞬态解，由于没有阻尼的影响，其瞬态解并不会衰减，其理论表达式为： (2-1)式中，为位移响应，为激励，为刚度，为荷载频率与固有振动频率之比，为荷载频率，为结构固有频率。现令为1，为1，则为1，取为2/3。程序求得的解与解析