结构动力学使用中心差分法计算单自由度体系动力反应的MATLAB程序_Image_Marked.pdfVIP

中心差分法计算单自由度体系动力反映的报告 前言 基于叠加原理的时域积分法与频域积分法一样，都假设结构在在全部反应过程中都是线性的。而时域 逐步积分法只是假设结构本构关系在一个微小的时间步距内是线性的，相当于分段直线来逼近实际的曲线。 时域逐步积分法是结构动力问题中研究并应用广泛的课题。中心差分法是一种目前发展的一系列结构动力 反应分析的时域逐步积分法的一种，时域逐步积分法还包括分段解析法、平均常加速度法、线性加速度法、 - 和 ‒ 法等。 中心差分法(central difference method)原理 中心差分法的基本思路 将运动方程中的速度向量和加速度向量用位移的某种组合来表示，将微分方程组的求解问题转化为代 数方程组的求解问题，并在时间区间内求得每个微小时间区间的递推公式，进而求得整个时程的反应。 中心差分法是一种显示的积分法，它基于用有限差分代替位移对时间的求导（即速度和加速度）。如 ∆ = ∆ ∆ 果采用等时间步长， （ 为常数），则速度与加速度的中心差分近似为 + + 1 ‒ 1 = 2∆ (1) ‒ 2 + + 1 ‒ 1 = 2 ∆ (2) 用 表示位移，离散时间点的运动为： = () ， = () ， = () ( = 0,1,2…) 体系的运动方程为 () + () + () = () (3) + 1 将速度和加速度的差分近似公式(1)和(2)代入(3) 中得出在 时刻的运动方程，将方程整理得到 由 ‒ 1 和 表示的两步法的运动方程(4) ： 2 (∆2 + 2∆) + 1 = ‒ ( ‒ ∆2) ‒ (∆2 ‒ 2∆) ‒ 1 (4) + 1 这样就