积分格式对数值模拟波动特性的影分析.docVIP

积分格式对数值模拟波动特性的影响分析 唐 晖1 李小军1 亓兴军2 1 中国地震局工程力学研究所，哈尔滨，150080 2 山东建筑大学土木工程学院，济南，250101 提要 本文通过分析时空离散波动数值模拟中振动传递函数的特性，探讨了显式积分格式对波动传播规律的影响。结果表明，所研究的显式积分格式应用于波动的有限元模拟时，波动的截止频率只受空间离散网格尺寸控制，而与时间离散步距无关。证明了在波动的有限元模拟中，所研究的显式积分格式的引入不影响离散网格中计算波动的截止频率。 关键词：地震波动，有限元模拟，显式积分格式，截止频率，寄生振荡 THE INFLUENCE OF AN EXPLICIT INTEGRATION FORMULA ON PROPERTY OF NUMERICAL SIMULATION OF WAVE PROPAGATION Tang Hui1 Li Xiaojun1 Qi Xingjun2 1 Institute of Engineering Mechanics, China Earthquake Administration, Harbin, 150080 2 School of Civil Engineering, Shandong Jianzhu University, Jinan,250101 Abstract The property of vibration transfer function is analyzed in this paper when wave is simulated numerically with time and space discretization, the effect of the integration formula on wave propagation is discussed. As the result shows, when the explicit integration formula analyzed in the paper is applied to finite element simulation of wave motion, the cut-off frequency of wave motion is only controlled by the size of space discretization and is independent of size of time discretization step. It is approved that the application of the explicit integration formula in finite element simulation of wave motion has no effect on the cut-off frequency of wave propagation. Key word: wave propagation, finite element simulation, explicit integration formula, cut-off frequency, parasitic oscillation 1 引言 波动问题的解法大致可以分为解析解法和数值解法，后者相对于前者在求解较为复杂的波动问题时更为有效，因而具有广泛的应用前景，其中以有限元法和有限差分法为主。文献[1]给出了一种将有限元法与有限差分法相结合的类似于差分法的有限元解法，这种方法不仅保持了有限差分的计算简便和便于编程等的优点，还利用有限元的优点消除了有限差分方法在处理非均匀介质和非规则边界时所遇到的困难。在应用该有限元法对波动问题进行求解的过程中，用有限元离散模型代替原连续体以及对时间域的离散，并结合中心差分格式，使得弹性波的传播规律发生变化[2,3]。然是，中心差分法在任意阻尼的情况下蜕变为隐式积分格式，就失去了显式积分格式计算量方面的优势。文献[4]给出了一种任意阻尼的情况均为显式格式的逐步积分方法。对于这一显式积分格式的计算精度、稳定性以及能耗特性，在一些文章中作了较为详细的分析[5，6,7]，而该显式积分格式应用于波动的有限元模拟时，对波动传播规律的影响尚未有具体的研究成果。为此，本文将对该问题进行具体研究。 2 显式积分方法 对于动力方程： (1) 本文要研究积分方法求解的积分格式[4]为 (2)