连续系统的振动之集中质量法、假设模态法、模态综合法和有限元法.pptVIP

质量矩阵 M 也可直接利用单元质量矩阵组集而成 方法：将单元质量矩阵 me1、me2 和 me3 的各个元素统一按 qi (i=1,2,3) 的下标重新编号，放入 M 中与编号相对应的行和列中 连续系统的振动 / 有限元法 单元质量矩阵： 和广义坐标 相对应的质量矩阵： 连续系统的振动 / 有限元法 全系统的势能： 也可组集得到： 连续系统的振动 / 有限元法 当杆上有常值轴向力作用时，三根杆的广义外力阵为： 系统的广义力阵： 作用力的总虚功： 与广义坐标 q 对应的广义力阵 也可将Fe1、Fe2 和 Fe3 的各个元素统一按 qi (i=1,2,3) 的下标重新编号，放入 Q 中与编号相对应的行和列中 连续系统的振动 / 有限元法 用广义坐标阵 q 表示的广义质量阵、广义刚度阵和广义外力阵： 用广义坐标阵 q 表示的全系统的动力学方程： 连续系统的振动 / 有限元法 梁的弯曲振动 单元质量矩阵和刚度矩阵的求解 将梁划分为多个单元 取出其中一个单元进行分析 ：单元形函数 单元长 l，节点横向位移 u1(t)、u3(t)，节点截面转角u2(t)、u4(t) 单元 x 位置截面的横向位移： 单元节点坐标： 选为均质梁在端点常值位移作用下的静挠度曲线： 连续系统的振动 / 有限元法 解得： 单元形函数应满足的边界条件： 1 连续系统的振动 / 有限元法 代入： 形函数列阵： 单元的动能： ：材料密度 ：截面积 单元质量矩阵 为常数时 连续系统的振动 / 有限元法 单元的势能： ：抗弯刚度 单元刚度矩阵 为常数时 假设梁上有分布外力，其虚功： ：与节点坐标ue 对应的单元广义力列阵 对于均布载荷， f 为常数，有： 连续系统的振动 / 有限元法 全系统的动力学方程 以上对单元所作的分析必须进行综合，以扩展到总体结构 以一个例子进行说明 梁划分为两个单元 单元质量矩阵： 单元刚度矩阵： 单元坐标阵 连续系统的振动 / 有限元法 全部节点坐标列阵： 节点坐标约束条件： 只有四个独立 定义独立的广义坐标： 广义坐标列阵： 节点坐标与广义坐标之间的关系： 全系统动能： 与q对应的质量阵 连续系统的振动 / 有限元法 因此，有： 单元质量矩阵： 与杆相同，M 也可直接将 me1、me2 的各个元素统一按 qi (i=1,2,3,4)的下标重新编号，放入对应的行列组集而成 同理可得与广义坐标 q 所对应的刚度矩阵 连续系统的振动 / 有限元法 系统的广义力阵： 作用力的总虚功： 与广义坐标 q 对应的广义力阵 也可将Fe1、Fe2 的各个元素统一按 qi (i=1,2,3) 的下标重新编号，放入 Q 中与编号相对应的行和列中 假设梁上有简谐变化的均布载荷 与节点坐标ue 对应的单元广义力列阵： 节点坐标与广义坐标之间的关系： 与节点坐标U对应的广义力阵 连续系统的振动 / 有限元法 用广义坐标阵 q 表示的广义质量阵、广义刚度阵和广义外力阵： 用广义坐标阵 q 表示的全系统的动力学方程： 连续系统的振动 / 有限元法 梁的稳态响应： 外力写成分布力形式： Ma 强迫振动方程： 广义力： 广义力列阵： 连续系统的振动 / 假设模态法 离散化强迫振动方程： 令： 坐标变换： 梁的稳态响应： 求得 得 代入梁的稳态响应方程中得解 连续系统的振动 / 假设模态法 假设模态法 动力学方程 瑞利法 里兹法 连续系统的瑞利法是基于能量法的假设模态法，是多自由度系统的瑞利法的推广 以梁的弯曲振动为例 假设梁以某阶模态函数作频率为 的自由振动： 设系统为保守系统，机械能守恒 即 引入系统的参考动能： 连续系统的振动 / 假设模态法 定义瑞利商： 参考动能： 当 为准确的第 i 阶模态函数时，瑞利商即为相应的特征值，即第 i 阶固有频率 若 是试函数，它满足梁的几何边界条件，但不满足动力学方程，则瑞利商是一个依赖于 的标量 试函数 越接近某阶真实模态，瑞利商越接近该阶固有频率 与多自由度系统相同，瑞利商大于基频 实际计算时可选择梁的静变形函数，或选择条件相近的梁的精确解作为试函数 连续系统的振动 / 假设模态法 若梁上存在集中质量和弹性支撑 则最大势能和参考动能相应改为： 连续系统的振动 / 假设模态法 例：等截面悬臂梁 端部有一集中质量 用瑞利法估计基频 解： 选择等截面悬臂梁在均布载荷下的静挠度曲线作为试函数： 选择端部集中质量作用下的静挠度曲线作为试函数： 因集中质量大于梁的分布质量，选用后一种试函数好 连续系统的振动 / 假设模态法 假设模态法 动力学方程 瑞利法 里兹法 里兹法是瑞利法的改进 瑞利法使用单个试函数，而里