第6章：动力有限元法.ppt

STDU DYNAMICS OF STRUCTURES 结构动力学 Dynamics of Structures 武兰河 QQ:281017372 TelChapter 6 Finite Element Method §6-1 有限元法的基本概念 §6-3 梁的弯曲振动 §6-2 杆的纵向振动 § 6.1 有限元法的基本概念 有限元法(The finite element method)是20世纪五六十年代发展起来的一种数值近似方法，最初用于求解飞机的应力问题 特点:兼有集中质量法和假设模态法的优点 方法:将复杂结构分割为有限个单元(element)，单元的端点称为节点(node)，将节点的位移作为广义坐标，并将单元的质量和刚度集中到节点上 每个单元作为弹性体，单元内各点的位移用节点位移的插值函数来表示，这种插值函数实际就是单元的假设模态 二. 有限元法的主要特点 1.以简单逼近复杂 2.采用矩阵表达，便于编制计算机程序 3.比较容易处理边界条件，可解决复杂域问题 4.适应性强，应用范围广泛 一. 有限元法的基本步骤 1.求解区域离散化 2.选择插值函数 3.分析单元的力学特性（建立单元的动力学方程） 4.集合各单元的平衡方程，建立结构的平衡方程 5.求解结构的平衡方程 § 6.2 杆的纵向振动 一. 结构的离散化 将杆件划分为多个单元 取出其中一个长度为l的单元 设单元两端的节点位移为 和 用一个向量表示 相应的两端杆端力向量为 与多自由度系统类似，该单元的运动方程可以写成 二. 位移插值 该单元内x点处的位移表示为 称为形函数(插值函数） 通常取一个节点坐标有单位位移而其余节点位移为零时单元的静变形函数为形函数 单元的假设模态 显然有 代入式 构成单元内部连续的位移场 写成矩阵形式 其中 形函数矩阵 单元位移列矩阵 三. 单元刚度矩阵和质量矩阵 因为轴力单元只有轴向变形，其微段的势能为 积分可以得到该单元的势能 将位移插值代入此式，整理后得到 其中 单元刚度矩阵 单元刚度矩阵的各元素为 若单元的材料常数和横截面积为常数,则有 类似地，单元微段的动能为 积分可以得到该单元的动能 单元质量矩阵 将位移插值代入此式，整理后得到 其中 单元刚度矩阵的各元素为 若单元的材料常数和横截面积为常数,则有 显然，单元刚度矩阵和质量矩阵均是对称的 设单元上作用有分布力 计算其对虚位移 的虚功 与节点位移相对应的单元广义力向量 若单元上作用有分布力 为常数 则等效结点力为 如此可将单元上的力等效转化到结点上 四. 单元动力平衡方程的建立 由Lagrange方程，可以建立单元的动力平衡方程 五. 系统的动力平衡方程 以图所示的变截面杆的纵向振动为例来说明综合的一般规律 将杆划分为三个单元 各单元的质量矩阵和刚度矩阵分别为 各单元的节点位移为 组合为全部节点的位移列阵 各节点位移之间有以下的约束关系 因此6个位移分量只有3个是独立的 定义以下的独立广义结点位移 并令 结构的节点位移列阵 则单元的节点位移与结构节点位移有如下关系 其中 计算整个系统的动能 其中 整个系统的质量矩阵 系统的质量矩阵也可以直接由单元质量矩阵组集而成 做法: 将各单元质量矩阵 的各元素按照 的下标编号,放入与编号相对应的行和列 然后将各个位置的元素相叠加即得到 矩阵 对号入座 得到 类似地,整个系统的势能 其中 整个系统的刚度矩阵 系统的刚度矩阵也可以直接由单元刚度矩阵组集而成 当杆上有常值力f作用时 各单元的等效结点力为 将3个杆单元的广义力组合为系统的广义力F 作用力的总虚功为 与广义位移q对应的广义力列阵 实际计算时,也可以由对号入座的方法,将各单元的结点力列阵元素按总体编码组集而得 将整个系统的动能、势能和虚功代入拉格郎日方程，得系统的动力学方程 以上对杆的纵向振动的计算方法也完全适用于轴的扭转振动等同一类型的振动。对于梁的弯曲振动，步骤完全相同，只是插值函数不同，未节点位移数目不同而已，有兴趣的同学请参照课本自学，由于学时关系，这里就不讲了 得到动力平衡方程之后，其他求解过程不须详述 § 6.3 梁的弯曲振动 一. 结构的离散化 将梁划分为多个单元 取出其中一个长度为l的单元 设单元两端的节点位移为 用一个向量表示 相应的杆端力向量为 其中 节点的横向位移 节点的截面转角 二. 位移插值 该单元内x点处的位移表示为 称为单元的形函数(插值函数） 单元的假设模态 通常取一个节点坐标有单位位移而其余节点位移为零时单元的静变形函数为形函数 单元形函数应满足的边界条件： 满足此边界条件的形函数为： 各