参考结构动力学.pptx

结构动力学

（2003秋）

;结构动力学

第五章

单自由度体系对任意荷载的反应;5.3数值计算方法—时域逐步积分法

一般时域逐步积分法的构造

从前面介绍的分段解析法和中心差分法给出的逐步积分公式可以发现，所谓的时域逐步积分方法就是构造出根据某一时刻及其以前时刻的运动，推算下一时刻运动的递推计算公式。

具体情况可表述为，设体系在ti及ti以前时刻的运动已知，求ti＋1时刻的运动（ti＝iΔt）。

体系在ti＋1时刻的运动包括：位移、速度和加速度，需要有三个方程（条件）求这三个量。因此，除体系的运动方程外，还需补充两个方程（条件）。;5.3数值计算方法—时域逐步积分法

一般时域逐步积分法的构造

两个补充方程可以通过对运动状态的假设得到。

例如可以假设在ti和ti＋1时刻，即Δt时间段内，体系的加速度为常数a，则积分(不定积分)得到体系的速度和位移为：

其中，τ为由ti时刻起算的局部时间坐标，c1和c2为积分常数，可由τ＝0时的初值条件

确定，最后得：;5.3数值计算方法—时域逐步积分法

一般时域逐步积分法的构造

当τ＝Δt，即t＝ti＋1时刻，体系得运动状态为：

而ti和ti＋1时间段内的常加速度a=(üi+1+üi)/2，则得到：

再加上ti＋1时刻的运动方程：

可以求得ti＋1时刻的位移、速度和加速度。;5.3数值计算方法—时域逐步积分法

一般时域逐步积分法的构造

以上方法也称为平均加速度法，即假设加速度为ti和ti＋1时间段内的平均值：

也可以假设加速度a为其它形式的变化规律，例如为线性变化：

则采用同样的分析步骤可以得到线性加速度法的时域逐步积分公式。;5.3数值计算方法—时域逐步积分法

一般时域逐步积分法的构造

平均加速度法和线性加速度法的基本假设和补充公式;5.3数值计算方法—时域逐步积分法

3、Newmark—β法

Newmark—β同样将时间离散化，运动方程仅要求在离散的时间点上满足。假设在ti时刻的运动均已求得,然后计算ti+1时刻的运动。

与中心差分法不同的是，它不是用差分对ti时刻的运动方程展开，得到外推计算ui+1的公式，而是通过对加速度的假设，以ti时刻的运动量为初始值，通过积分得到计算ti+1时刻的运动公式。

与平均加速度法和线性加速度法不同的是，它用不同的加速度假设条件给出速度和位移的计算公式。;5.3数值计算方法—时域逐步积分法

3、Newmark—β法

Newmark—β法假设在时

间段[ti，ti+1]内，加速度为

一常量，记为a，经过简单

积分计算可以得到速度、

位移与a之间的关系式：;5.3数值计算方法—时域逐步积分法

3、Newmark—β法

;5.3数值计算方法—时域逐步积分法

3、Newmark—β法

;5.3数值计算方法—时域逐步积分法

Newmark—β法的时域逐步积分公式

;5.3数值计算方法—时域逐步积分法

3、Newmark—β法

在Newmark—β法中，控制参数β和γ的取值影响着算法的精度和稳定性，可以证明，只有当γ取1/2时，这个方法才具有二阶精度，因此一般均取：

γ=1/2,0≤β≤1/4

Newmark—β法的稳定性条件：

当γ=1/2,β=1/4时，Δt≤∞，即成为无条件稳定的。;3、Newmark—β法

通过对Newmark—β法中控制参数β取不同的值也可以得到其它时域逐步积分方法。下表给出了β取不同值时Newmark—β法所对应的逐步积分法。;3、Newmark—β法

在动力问题研究中，常采用γ=1/2，β=1/4的所谓无条件稳定的Newmark—β法，实际上就是平均（常）加速度方法。

Newmark—β法为单步法不需要格外处理计算的“起步”问题，属于自起步方法。;Wilson—θ法

在时域逐步积分法发展的早期，Wilson—θ法曾得到广泛应用。Wilson—θ法是基于线性加速度法基础之上发展的。当参数θ1.37时，方法是无条件稳定的。初略分析，采用了线性加速度假设比平均常加速度法更精确，而且是无条件稳定的，应是一种优秀的逐步积分法。但随着对数值算法特性研究的深入，发现Wilson—θ法存在一系列弊病。而且对于一些强冲击问题，Wilson—θ法可能无法完成计算。;5.3数值计算方法—时域逐步积分法

不同方法的计算精度(单自由度体系无阻尼自由振动,Δt/Tn=0.1);5.3数值计算方法—时域逐步积分法

不同方法的振幅衰减AD(Amplitueddecay)和周期延长PE(Periodelongation);5.3