一类非耗散显式时间积分方法.docVIP

一类非耗散的显式时间积分方法 　　摘要： 针对大型复杂结构的动力学问题，为了得到高效的计算方法，以泰勒公式为基础，通过调整参数改善算法的稳定性，得到了一类具有非耗散特性且以加速度为基本变量的显式数值积分方法。对提出的方法进行了稳定性和精度分析，并通过多个算例对提出的方法、蛙跳式中心差分法、Newmarkβ法和精确解进行比较。结果表明，无阻尼情况下，所提出的方法的稳定性条件与中心差分法相同；提出的方法的振幅衰减率为0，具有非耗散特性，其周期误差约为隐式Newmarkβ法的一半；所提出的方法给出了蛙跳式中心差分法和无耗散特性的翟方法的统一格式，并可以衍生出更多的显式时间积分方法。 　　关键词： 显式方法； 非耗散； 差分格式； 精度； 稳定性 　　中图分类号： O241.8； O313 文献标志码： A 文章编号：2015 　　DOI：10.16385/ki.issn2015.03.014 　　引 言 　　结构动力分析的数值解法可以分为模态叠加法和直接积分法，其中直接积分法按照求解方式还可以分为隐式法和显式法。隐式法有Newmarkβ法[1]和Wilsonθ法[2]等，通常应用于低频响应占主导地位的长时间结构响应问题或准线性问题。显式法中最具代表性的是中心差分法[3]，通常应用于高频响应占主导地位的非线性问题或瞬态问题。对于大型非线性的动力学计算问题，显式方法的应用更加广泛。随着结构动力响应分析向复杂非线性的方向发展，计算量和计算代价显著提高，计算方法的速度和精度越来越受到重视。 　　文献[4]提出了一种具有三阶精度的显式差分方法，适用于有阻尼系统的动力响应分析，高阶精度需要增加中间变量而使整体计算量增加[5]。隐式算法也可以转换为显式算法，一般以位移为第一求解变量[6]。文献[7]提出了Newmark更新精细积分法，具有良好的稳定性，但要求阻尼矩阵为对角矩阵，所以只能属于半显式的方法。很多新开发的算法一般都是基于泰勒展开或加权残余法[8]，应用泰勒公式展开并截除高阶小量会引入一定的误差，而加权残余法也同样存在一定的误差，所以在考虑计算速度的同时要尽量提高算法精度。 　　现有的研究大多以位移为基本变量，速度和加速度通过对位移的求导获得，而以加速度为基本变量的方法很少受到重视。本文以加速度为基本变量，速度和位移通过对加速度的近似积分获得，在泰勒展开的基础上，通过调整积分参数改善算法的稳定性，提出了一类无耗散特性的显式时间积分方法。文中提出的方法的稳定性条件与中心差分法相同，其周期误差仅为Newmarkβ法的一半。文中提出的显式时间积分方法统一了蛙跳式中心差分法和翟方法（=φ=0.5），并可以衍生出更多的无耗散特性的显式积分格式。 　　1 显式方法 　　显式是指增量步结束时的状态仅依赖于该增量步开始之前的位移、速度和加速度，计算时不需要进行迭代或矩阵求逆。结构的动力学运动方程为M（tn）+C（tn）+Ku（tn）=F（tn） 　　（1）式中 M，C，K分别为振动系统的质量矩阵、阻尼矩阵、刚度矩阵；，，u，F分别为相应时刻的加速度、速度、位移和外力矢量。 　　由稳定性原理可知，λi为矩阵A的第i个特征值，当矩阵A的谱半径ρ（A）=max|λi| 　　如图1所示，在无阻尼情况下，α1时，谱半径存在小于1的区间，积分格式是收敛的，存在算法阻尼，即振幅会随着时间增长而不断衰减；当α=1且β≠γ时，积分格式存在起始点不为0的稳定区间，这意味着步长取某值时算法是稳定的，但进一步减小时间步长时算法可能不稳定；只有在α=1并且β=γ时，通过分析高次方程的解，无论γ取何值，谱半径在Ω∈（0，2）都是恒等于1的，振幅既不衰减也不发散而保持恒定，此时算法具有非耗散特性。图2显示了阻尼对加速度显式法的稳定性影响，随着阻尼比的增大，算法稳定区间逐渐减小。谱半径越小，表示算法收敛越快。 　　3.2 精度分析 　　积分逼近算子矩阵具有一对共轭复根，可以用这对共轭复根表示出振动方程的解λ1，2=a±ib=e（-ξ±i） 　　（16）式中 =arctan（|b/a|），ξ=-ln（a2+b2）/（2），a和b为实数。令=Δt，对于欠阻尼情况，式（13）的数值解可以表示为u（tn）=e-ξtn（c1costn+c2sintn） 　　（17）式中 c1和c2为常数，tn=nΔt，与无阻尼自由振动精确解相比，ξ表示算法阻尼比，表示算法圆频率。 　　振幅衰减率AD反映了幅值误差，周期延长率TD反映了周期误差，见式（18）和（19），因此算法的精确度可以由振幅衰减率和周期延长率来度量。AD=1-e-2πξ 　　（19） 在无阻尼情况下，将加速度显式法、Newmarkβ法和中心差分法