[物理]机械振动数值分析.ppt

中心差分法 将式（14）代入任意时刻的运动方程 得 其中 中心差分法 起始步的计算 定义积分常数 中心差分法 计算步骤 1、初始计算 1） 形成刚度矩阵 和质量矩阵 和阻尼矩阵 ； 2） 给定初始条件 和 ，并计算 ； 3） 给定时间步长，计算积分常数； 4） 计算起步位移： ； 5） 计算有效质量矩阵： ； 6） 对有效质量矩阵进行三角分解： ； 中心差分法 2、对每一个时间步 1） 计算时刻 的有效载荷 ； 2） 计算 时刻的位移： ； 3） 如有需要，计算时刻 的加速度和速度 中心差分法 中心差分法的数值稳定性 对任意时刻的运动方程 考虑其中心差分格式，有 将式（21）写为递推矩阵格式有 中心差分法 其中 对无阻尼系统，由式（23）计算出递推矩阵特征值为 中心差分法 根据收敛条件，可知 从而得 对有阻尼系统，中心差分法的临界时间步长为 中心差分法 由于对任意一个解耦后的运动方程来说均应满足式（27），因此，满足收敛条件的时间步长由系统中的最高阶模态对应的固有频率决定。容易看出，当满足式（25）时，谱半径为 1 ，此时中心差分法没有数值阻尼。 Newmark法 在位移和速度的Taylor级数中只保留一阶导数项，则可得 若在式（28）中将导数的平均值代替 时刻的导数值，则可将式（28）改写为 Newmark法 将式（29）中的第二式代入第一式得 从式（29）中可以看出，加速度近似取平均加速度，因此该方法又称为平均加速度法。若在积分过程中对加速度取不同近似格式，可得到各种不同的积分格式。 Newmark引入如下速度、位移关系 未知量 Newmark法 由式（31）中的第二式可解出 再将式（32）代入式（31）中的第一式可解出 考虑 时刻的系统运动方程 Newmark法 将式（32）和式（33）代入式（34），得 简记为 Newmark法 令 并记 得 计算步骤 1、初始计算 1） 形成刚度矩阵 和质量矩阵 和阻尼矩阵 ； 2） 给定初始条件 和 ，并计算 ； 3） 给定时间步长，计算积分常数； 4） 形成有效刚度矩阵 ； 5） 对有效刚度矩阵进行三角分解 ； Newmark法 2、对每一个时间步 1） 计算 时刻的有效载荷 ； 2） 计算 时刻的位移： ； 3） 计算 时刻的加速度和速度 Newmark法 Newmark法的数值稳定性 仍然讨论解耦后的无阻尼模态运动方程 将Newmark法格式的加速度近似表达式代入，得 利用式（31），并考虑到如下关系 Newmark法 得 从式（42）中可联立解得 由于时间的任意性，可将上式改写为 Newmark法 将上式代入式（40），得 写成递推矩阵格式，有 Newmark法 其中 将上式代入式（40），得 Newmark法 递推矩阵的特征值为 要求Newmark法稳定，即要求 Newmark法 由式（49）的第一式和第二式可分别得到 从式（51）可解出，Newmark法无条件稳定的条件为 Newmark法 若不满足式（52）所示条件，则算法是条件稳定的，要得到稳定的解，要求计算时间步长小于临界时间步长，即 雅可比变换法 经第一次循环，对非对角元（1，2）消零，得 由式（52）可计算得 雅可比变换法 类似的，分别对非对角元（1，3）和（2，3）消零，得 雅可比变换法 雅可比变换法的一般步骤 给定当前循环的阀值，一般取第m次迭代的阀值为10-2m； 对上三角块的非对角元按式（48）计算，若计算结果大于给定的阀值，则按式（47）进行旋转变换； 检查式（49）式是否满足，若满足则由迭代最终步骤计算出的矩阵给出特征值，计算特征向量，否则转向第1步，开始新一轮循环。 广义雅可比变换法 对雅可比变换矩阵进行修改，得转换矩