第五章有限元动力学基本原理说课.ppt

四、机械结构动力响应的计算 1.振型叠加法 于是，振动方程解耦为 依次用求解常微分方程的方法可以得到解： 再代回 可以得到问题的解。 四、机械结构动力响应的计算 1.振型叠加法 通常情况下，高阶振型对动力响应影响较小，因此，只取最低的3～5阶（≤10）阶振型就可以得到满意精度的动力响应。 2.直接积分法 对于动力学方程 假设已知其初始条件 如果将求解时间T等分成nT个时间区间 ，并能通过前若干个时刻的解来确定下一时刻的解，就可获得问题的解，直接积分法可以解决这一问题。 四、机械结构动力响应的计算 2.直接积分法 直接积分中对动力学方程是逐步地进行数值积分的，“直接”的意思是指：进行数值积分前没有进行把方程变为另一种形式的变换。 直接积分有两个假设，一个是动力学方程的解只在相隔 的一些离散时间区间上满足方程，而不要求在任意时刻都满足方程；另一个是假设位移速度和加速度在每一个时间区间 内按一定的规律变化。 直接积分法有中心差分法、Houbolt法、Wilson- 法、Newmark和龙格-库塔法。 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 对于加速度 上式的误差为 高阶小量 对于速度 方程在t时刻为 将速度和加速度的差分格式代入方程，可以得到 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 显然，求解 ，需要2个初始条件 和 我们有初始条件 ，根据 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 可以解得 注意，有 于是，可以求得方程的解。归纳起来后，中心差分法的计算步骤为2大步，9小步。 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 第一步：初始计算 1).形成 2).计算初始值 3).选取时间步长 ,并计算积分常数 4).计算 5).形成有效质量阵: 6).做三角分解: 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 第二步：对每一时间步长 7).计算在时刻t的有效载荷 8).求解在时刻 的位移 9).如果需要，计算时刻 的速度和加速度 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 说明： 1). 的取值 是结构最高频率对应的周期，是最小周期。 2).求解方程是显式差分，对于对角线质量阵和无阻尼振动求解很方便，可以直接求解，不需进行三角分解。 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 例题：一个二自由度振动系统，其动力方程为： 已知该系统的自由振动周期为 ，试用中心差分法求解步长为 和 时，方程的解。 解：取12个步长的系统响应，假设 计算 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 考虑 ，有 因而有 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 对每个时刻步长求解方程 2.60 2.54 2.78 3.37 4.17 4.90 5.26 5.02 4.14 2.83 1.45 0.39 δ2 1.02 2.04 2.77 3.07 2.91 2.40 1.70 1.02 0.49 0.17 0.03 0 δ1 12△t 11△t 10△t 9△t 8△t 7△t 6△t 5△t 4△t 3△t 2△t △t 时间 四、机械结构动力响应的计算 2.直接积分法 ●Houbolt法 使用新的差分格式 误差同样为 的高阶小量 方程在t+ 时刻为 四、机械结构动力响应的计算 2.直接积分法 ●Houbolt法 将速度和加速度的差分格式代入方程，可以得到 显然，求解 ，需要3个初始条件 、 和 一般有初始条件 ， 可用特殊的方法获得。 四、机械结构动力响应的计算 2.直接积分法 ●Houbolt法 归纳后，Houbolt法的计算步骤为： 第一步：初始计算 1).形成 2).计算初始值 3).选取时间步长 ,并计算积分常数 四、机械结构动力响应的计算 2.直接积分法 ●Houbolt法 4).使用特殊方法计算： 5).形成有效刚度阵: 6).做三角分解: 第二步：对每一时间步长 7).计算在时刻t+ 的有效载荷 8).求解在时刻 的位移 四、机械结构动力响应的计算 2.直接积分法 9).如果需要，计算时刻 的速度和加速度