微分求积时间单元方法.pdfVIP

第九部分 结构动力学 微分求积时间单元方法 1 郭静 ，邢誉峰 （北京航空航天大学固体力学研究所，北京市海淀区学院路37 号，北京100191） 摘要：基于微分求积法则，提出了一种求解动力学常微分方程的高效高精度微分求积时间单元方法 （DQTEM ）。给出了DQTEM 施加初始位移和初始速度的方法，其结果相当于构造了C1 时间单元。与递 推格式的直接积分方法不同，对于考虑的时间域通常只需用一个微分求积时间单元。与RK 法和Newmark 法相比，用少量时间结点的DQTEM 结果就与精确解吻合。稳定性分析表明，DQTEM 通常是条件稳定的。 关键词：微分求积法则；直接积分方法；时间单元；相位误差；数值耗散 中图分类号：O241.81; O321 Differential Quadrature Time Element Method 1 Guo Jing , Xing Yufeng (The Solid Mechanics Research Center, Beihang University, Beijing, 100191, P.R. China) Abstract: Based on differential quadrature rule, an accurate and efficient differential quadrature time element method (DQTEM) is proposed for solving dynamical ordinary differential equation, whose numerical dissipation and phase error is much smaller than the conventional direct integration method. The method of imposing initial displacements and velocities of DQTEM is given; subsequently the C1 time element is constructed. On the contrary to the recursive direct integration method, one differential quadrature time element is usually enough for the whole time domain under consideration. Comparison with the RK method and the Newmark method, DQTEM solutions with an evidently smaller number of sampling points agree extremely well with the exact solutions. Stability analysis indicates that DQTEM is usually conditionally stable. Key words: differential quadrature rule; direct integration method; time element; phase error; numerical dissipation 1 引言 动力学方程描述了系统的运动学特性和动力学特性。最常用的动力学算法是基于有限 差分的直接积分方法[1] 。经典直接积分方法，如RK 方法、中心差分方法、Wilson  法和 Newmark 法等，均存在累积的数值耗散或相位误差，因此其长期跟踪能力不尽人意。虽然 可以通过减小时间步长达到降低耗散减小相位误差以提高计算精度的目的，但这要以提高 计算量