一类非耗散的显式数值积分方法-振动工程.DOCVIP

一类非耗散的显式积分方法β法和精确解进行比较。结果表明，无阻尼情况下，文中提出的方法的稳定性条件与中心差分法相同；提出的方法的振幅衰减率为0，具有非耗散特性，其周期误差约为隐式Newmark-β法的一半；文中提出的方法给出了蛙跳式中心差分法和无耗散特性的翟方法的统一格式，并可以衍生出更多的显式时间积分方法。( 关键词：非耗散，显式方法，差分格式，精度，稳定性 中图分类号：O241.8；O313 文献标识码：A 文章编号： 引 言 结构动力分析的数值解法可以分为模态叠加法和直接积分法，其中直接积分法按照求解方式还可以分为隐式法和显式法。隐式法有Newmark-β法[1]和Wilson-θ法[2]等，通常应用于低频响应占主导地位的长时间结构响应问题或准线性问题。显式法中最具代表性的是中心差分法[3]，通常应用于高频响应占主导地位的非线性问题或瞬态问题。对于大型非线性的动力学计算问题，显式方法的应用更加广泛。随着结构动力响应分析向复杂非线性的方向发展，计算量和计算代价显著提高，计算方法的速度和精度越来越受到重视。 文献[4]提出了一种具有三阶精度的显式差分方法，适用于有阻尼系统的动力响应分析，高阶精度需要增加中间变量而使整体计算量增加[5]。隐式算法也可以转换为显式算法，一般以位移为第一求解变量[6]。文献[7]提出了Newmark更新精细积分法，具有良好的稳定性，但要求阻尼矩阵为对角矩阵，所以只能属于半显式的方法。很多新开发的算法一般都是基于泰勒展开或加权残余法[8]，应用泰勒公式展开并截除高阶小量会引入一定的误差，而加权残余法也同样存在一定的误差，所以在考虑计算速度的同时要尽量提高算法精度。 现有的研究大多以位移为基本变量，速度和加速度通过对位移的求导获得，而以加速度为基本变量的方法很少受到重视。本文以加速度为基本变量，速度和位移通过对加速度的近似积分获得，在泰勒展开的基础上，通过调整积分参数改善算法的稳定性，提出了一类无耗散特性的显式时间积分方法。文中提出的方法的稳定性条件与中心差分法相同，其周期误差仅为Newmark-β法的一半。文中提出的显式时间积分方法统一了蛙跳式中心差分法和翟方法()，并可以衍生出更多的无耗散特性的显式积分格式。 1.显式方法 显式是指增量步结束时的状态仅依赖于该增量步开始之前的位移、速度和加速度，计算时不需要进行迭代或矩阵求逆。结构的动力学运动方程为： (1) 式中，M、C、K分别为振动系统的质量矩阵、阻尼矩阵、刚度矩阵；、、、分别为相应时刻的加速度、速度、位移和外力矢量。 1.1中心差分法 文献[9,10]中通常所描述的中心差分法(CDM)一般都是半显式格式的，在阻尼矩阵为非对角的情况下会退化为隐式格式。中心差分法对运动方程直接进行积分，其主要公式如下： (2) 将式(2)带入式(1)，得到以位移为第一求解变量的运动方程 (3) 中心差分法先利用式(3)求得tn+1时刻的位移，然后根据式(2)进行两次求导得到tn时刻的速度和加速度。这种半显式的中心差分法是以位移为基本变量，不存在算法阻尼，具有二阶精度。 1.2蛙跳式中心差分法 (4) (5) (6) 在获得tn时刻的加速度以后，由式(5)得到tn+1/2时刻的速度，最后根据式(6)得到tn+1时刻的位移。位移、速度和加速度是在不同时刻交错求解的。需要指出的是，蛙跳格式以加速度为基本变量，没有直接给出整数时刻的速度。式(4)中的速度项是用来代替实际的速度，因为在每个循环中，同一时刻的加速度要比速度优先获得。蛙跳式中心差分法是显式有限元软件LS-DYNA的最常用的积分方法[11]。 1.3 翟方法 (7) 式中，为积分参数，一般取为0.5；在得到tn+1时刻的位移和速度后，就可以由式(8)得到该时刻的加速度。 (8) 翟方法具有二阶精度，在积分参数都为0.5时不存在算法阻尼，与Newmark-β法具有类似的算法阻尼特性。该显式方法也是以加速度为基本变量，只要其质量阵是对角阵或经对角化而不论阻尼阵如何复杂，积分过程中都无须联立求解耦合方程组，避免了组装等效刚度矩阵和矩阵求逆，这就大大节省了计算内存和时间，因而用于分析工程中大型非线性结构动力问题时，可以显著地提高计算速度和效率。 2.基于加速度的显式积分方法 (9) 假设加速度的导数满足向后差分公式，则有 (10) 把式(10)带入式(9)，并略去高阶项，得到 (11) 将式(11)中值不为1的系数以变量代替 (12) 式(12)中，推荐取α=1, β=γ, γ∈[0,1]，γ的取值尽量接近0.5，此时该算法仅有唯一的可变参量γ，具有非耗散特性，即不存在算法阻尼，这个结果是通过下文的稳定性分析得