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

第28卷第3期 V01．28No．3 砒 E Jun．2015 2015年6月 O振眦m动甜 V工h 程∞ 学嚼 腿 巷一g 一类非耗散的显式时间积分方法 杨 超，肖守讷，阳光武，朱 涛，杨 冰 (西南交通大学牵引动力国家重点实验室，四川成都610031) 摘要：针对大型复杂结构的动力学问题，为了得到高效的计算方法，以泰勒公式为基础，通过调整参数改善算法的 稳定性，得到了一类具有非耗散特性且以加速度为基本变量的显式数值积分方法。对提出的方法进行了稳定性和 精度分析，并通过多个算例对提出的方法、蛙跳式中心差分法、Newmark一口法和精确解进行比较。结果表明，无阻 尼情况下，所提出的方法的稳定性条件与中心差分法相同；提出的方法的振幅衰减率为0，具有非耗散特性，其周期 误差约为隐式Newmark-fl法的一半；所提出的方法给出了蛙跳式中心差分法和无耗散特性的翟方法的统一格式， 并可以衍生出更多的显式时间积分方法。 关键词：显式方法；非耗散；差分格式；精度；稳定性 中图分类号：0241．8；0313文献标志码：A 文章编号：1004—4523(2015)03—0441—08 DOI：10．16385／j．cnki．issn．1004—4523．2015．03．014 现有的研究大多以位移为基本变量，速度和加 引 言 速度通过对位移的求导获得，而以加速度为基本变 量的方法很少受到重视。本文以加速度为基本变 结构动力分析的数值解法可以分为模态叠加法 量，速度和位移通过对加速度的近似积分获得，在泰 和直接积分法，其中直接积分法按照求解方式还可 勒展开的基础上，通过调整积分参数改善算法的稳 定性，提出了一类无耗散特性的显式时间积分方法。 以分为隐式法和显式法。隐式法有Newmark-fl 文中提出的方法的稳定性条件与中心差分法相同， 法[1]和Wilson一0法[21等，通常应用于低频响应占主 导地位的长时间结构响应问题或准线性问题。显式 其周期误差仅为Newmark-fl法的一半。文中提出 法中最具代表性的是中心差分法[3]，通常应用于高 的显式时间积分方法统一了蛙跳式中心差分法和翟 频响应占主导地位的非线性问题或瞬态问题。对于 方法(乒一9—0．5)，并可以衍生出更多的无耗散特 大型非线性的动力学计算问题，显式方法的应用更 性的显式积分格式。 加广泛。随着结构动力响应分析向复杂非线性的方 向发展，计算量和计算代价显著提高，计算方法的速 1 显式方法 度和精度越来越受到重视。 文献I-4]提出了一种具有三阶精度的显式差分 显式是指增量步结束时的状态仅依赖于该增量 方法，适用于有阻尼系统的动力响应分析，高阶精度 步开始之前的位移、速度和加速度，计算时不需要进 需要增加中间变量而使整体计算量增加口]。隐式算 行迭代或矩阵求逆。结构的动力学运动方程为 法也可以转换为显式算法，一般以位移为第一求解 磁(≠。)+a；：(f。)+勋(￡。)一F(t。)(1) 变量[6]。文献[7]提出了Newmark更新精细积分式中M，c，K分别为振动系统的质量矩阵、阻尼矩 法，具有良好的稳定性，但要求阻尼矩阵为对角矩 阵、刚度矩阵；五，五，H，F分别为相应时刻的加速度、 阵，所以只能属于半显式的方法。很多新开发的算 速度、位移和外力矢量。 法一般都是基于泰