中心差分法的基本理论与程序设计.pdfVIP

中心差分法的基本理论与程序设计 1 程序设计的目的与意义 该程序通过用 C 语言（部分 C++语言）编写了有限元中用于求解动力学问 题的中心差分法，巩固和掌握了中心差分法的基本概念，提高了实际动手能力， 并通过实际编程实现了中心差分法在求解某些动力学问题中的运用， 加深了对该 方法的理解和掌握。 2 程序功能及特点 该程序采用 C 语言（部分 C++语言）实现了用于求解动力学问题的中心差 分法，可以求解得到运动方程的解答，包括位移，速度和加速度。计算简便且在 算法稳定的条件下，精度较高。 3 中心差分法的基本理论 在动力学问题中，系统的有限元求解方程（运动方程）如下所示： Ma t Ca t Ka t Q t 式中， a t 是系统结点位移向量， a t 和 a t 分别是系统的结点加速度向 量和结点速度向量， M ,C , K 和 Q t 分别是系统的质量矩阵、阻尼矩阵、刚度矩 阵和结点载荷向量，并分别由各自的单元矩阵和向量集成。 与静力学分析相比，在动力分析中， 由于惯性力和阻尼力出现在平衡方程中， 因此引入了质量矩阵和阻尼矩阵， 最后得到的求解方程不是代数方程组， 而是常 微分方程组。常微分方程的求解方法可以分为两类， 即直接积分法和振型叠加法。 中心差分法属于直接积分法， 其对运动方程不进行方程形式的变换而直接进 行逐步数值积分。通常的直接积分是基于两个概念，一是将在求解域 0 t T 内 的任何时刻 t 都应满足运动方程的要求，代之仅在一定条件下近似地满足运动方 程，例如可以仅在相隔 t 的离散的时间点满足运动方程； 二是在一定数目的 t 区 域内，假设位移 a 、速度 a 、加速度 a 的函数形式。 中心差分法的基本思路是用有限差分代替位移对时间的求导， 将运动方程中 的速度和加速度用位移的某种组合表示， 然后将常微分方程组的求解问题转换为 代数方程组的求解问题， 并假设在每个小的时间间隔内满足运动方程， 则可以求 得每个时间间隔的递推公式，进而求得整个时程的反应。 在中心差分法中，加速度和速度可以用位移表示，即： 1 a a 2 a a 2 t t t t t t 1 a at t at t 2 t 时间 t t 的位移解答 at t ，可由时间 t 的运动方程应得到满足，即由下式： Ma Ca Ka Q t t t t 而得到。为此将加速度和速度的表达式代入上式中， 即可得到中心差分法的递推 公式： 1 1