机电系统动力学第七章第1234节.ppt

* 振动分析的方法很多，数值仿真方法是进行振动分析的最直接的一类方法，它们可以应用于包括非线性振动在内的各种振动问题，这类方法是用以研究动态响应有效手段之一。 第七章 振动的仿真 从数学的观点来看，数值仿真方法是解微分方程边值问题和初值问题的逐步方法。在结构动力学响应计算方面采用实用有效的数值仿真方法，可以对系统在任意激励下的动态响应进行分析。 数值仿真方法的特点： ?在时间域内对响应的时间历程进行离散，把运动微分方程分为各离散时刻的方程； ?将某时刻的速度和加速度用相邻时刻的各位移的线性组合表示，于是系统的运动微分方程就化为一个由位移组成的某离散时刻的代数方程组； ?对耦合的系统运动微分方程进行逐步数值积分，从而求出在一系列离散时刻上的响应值。 ? 这种数值仿真方法称为逐步积分法(或直接积分法)。 目前用于求解多自由度线性振动系统的常用方法有： ?中心差分法； ?侯博特(Houbolt)法； ?威尔逊(Wilson-?)法； ?纽马克(Newmark-?)法等数值方法。 ? 对于高频分量和低频分量混合的问题，若采用无条件稳定的解法，可以提高计算效率。 ◆中心差分法是直接积分法的一种。 ◆它是将系统的运动微分方程在时间域内离散，化成对时间的差分格式，然后根据初始条件，利用逐步积分求出在一系列离散时刻上的响应值。 7.1 中心差分法 离散系统的运动微分方程为 (7.1-1) 式中M，C，K分别为系统的质量矩阵，阻尼矩阵和刚度矩阵； , ，x分别表示系统的加速度，速度和位移；R(t)是外力向量。 以确定在时刻 的近似解。 假定在t=0时，位移、速度和加速度分别为已知的x0, ， 。为了求解系统的运动微分方程在时间区间[0，T]的解，把时间全程T划分为n等份，即： 在中心差分法中，是按中心差分将速度和加速度向量离散化为 (7.1-2) 由以上两式可见，t时刻的速度和加速度是以相邻时刻的位移表示的。考虑在t时刻的动力方程为 (7.1-5) 将式(7.1-3)和式(7.1-4)代入式(7.1-5)，有 (7.1-6) 式中 (7.1-7) (7.1-8) ◆求解方程式(7.1-6)，可得xt+?t。 ◆但由式(7.1-8)可以看出，为求xt+?t必须使用xt和xt-?t的值。 ◆可见开始计算时，即t=0时，要计算x?t的值，就需要未知的x-?t值。 ◆因此应该有一个起始技术，因而这种算法不是自起步的。 ◆由于x0， , 是已知的，由t=0时的式(7.1-3)和式(7.1-4)可得 (7.1-9) 表7.1-1 中心差分法的计算机实施格式 A. 初始计算 1. 形成质量矩阵M，阻尼矩阵C和刚度矩阵K。 2. 给出初始值 x0， ， 。 3. 选择时间步长△t，△t?△tcr，并计算积分常数： ， ， ， 。 4. 计算 。 5. 形成有效刚度矩阵 ： 6. 对 作三角分解： B. 关于每一时间增量计算 1. 计算t时刻的有效载荷 2. 计算t＋△t时刻的位移 3. 如果需要，计算t时刻的加速度和速度 ◆中心差分法是一种显式积分方法。 ◆需要特别指出的是，使用中心差分法必须考虑积分的时间步长△t不能大于临界值△tcr，即 式中Tn为离散系统的最小周期。 (7.1-10) ◆如果不满足式(7.1-10)，数值解将出现发散现象，表明这种算法不是无条件稳定的。 ◆侯博特(Houbolt)法是Houbolt为研究飞机振动所提出的方法。 ◆该方法是以三级位移插值为基础的，通过四点的位移建立三次式，用两个向后差分公式表示在时刻t+△t的速度和加速度，即 7.2 侯博特法 (7.2-1) (7.2-2) 于是，在t＋△t时刻的动力方程为 将加速度和速度的向量形式(7.2-1)式和式(7.2-2)代入式(7.2-3)，即得关于xt+?t的求解方程 式中 (7.2-3) (7.2-4) (7.2-5) (7.2-6) ◆由式(7.2-4)-式(7.2-6)可以看出，要计算xt+?t时刻的解，必须使用前三步的位移xt, xt-?t和xt-2?t。 ◆由于该方法不是自起步的，要用其它方法由x0， ， 起步，例如可用中心差分法求出xt和x2?t后，才能