结构动力学的有限单元法.pptVIP

结构动力学的有限单元法 结构动力学的有限单元法 一、有限元法简介 二、采用有限元法分析动力学问题 三、结构动力学方程 四、方程的特征值及振型的正交性 五、结构动力学响应的求解方法 六、有限元法求解动力学问题实例 一、有限元法简介 有限元法的基本思想： 采用单元将实际上连续的弹性物体离散化，并且各单元彼此之间以节点联结，单元上选取简单的函数组合作为位移模型，利用弹性力学的变分原理原理来获得单元的运动方程组。然后，按照一定的规则把所有单元的运动方程组集合起来，经适当的边界条件处理，便得到整个物体的总体运动方程组。最后，选择适当的方法来积分总体运动方程组，问题的解答将在物体的各离散点上给出。 有限元法分析的基本过程： 连续体的离散化 选取单元的位移模型 建立单元的运动方程 经整体集合建立总体运动方程 求解运动方程 计算结构的应力应变或其它动力学特性 （1）连续体的离散化 用一些假想的面或线将物体所占据的空间区域分割成一系列的子区域，每个子区域就叫做单元或元素。各单元彼此之间仅在有限的指定点(称为节点)处相互连接。 三个方面： a)连续体本身的离散化； b)作用于连续体上的力系离散化 c)边界条件的离散化 连续体的离散化需要注意的几点： 1）单元的密度 （单元越密、节点越多，精度越高，计算量越大。机械结构有限元动力学分析，计算固有频率和振型，可以考虑将网格划分得粗一些。） 2）采用疏密不等的网格划分方法。 3）在采用三角形单元时，要尽量使每个三角形单元边长不要相差过大，不要出现过尖、过钝的内角，避免计算结果出现大的误差。 4) 要尽量将单元的节点和单元边界设置在几何形状、材料特性和载荷发生突变处。 （2）选取单元的位移模型 为了确定单元内任意点M在坐标轴方向上的位移，必须假设一个位移函数，以使点M的位移由单元节点位移通过位移插值函数来获得。 多项式便于微分和积分，增加多项式的阶数可以改善结果的精度，所以位移插值函数大多选用完全的多项式形式： 例：三节点三角形平面单元 （每个节点两个位移分量，三个节点共六个自由度，所以多项式包含六个待定系数，每个位移分量各取三项） 例：平面梁单元 平面梁单元形状函数矩阵N 为了使有限元法的解在单元尺寸逐步趋小时能够收敛于正确解，单元位移函数需要满足的条件： A、位移模式中必须包括反应刚体位移的常数项。 B、位移模式中必须包括反应常应变的线性位移项。 C、位移模式中必须能保证单元之间的连续性。 二、采用有限元法分析动力学问题 动力学问题的有限元法也同结构静力学问题一样，要把物体离散为有限个数的单元体。在考虑单元特性时，物体所受到的载荷还要考虑单元的惯性力和阻尼力等因素。 一般地，单元的刚度矩阵，质量矩阵和阻尼矩阵可以分别表示为： 三、结构的动力学方程 单元的刚度矩阵、质量矩阵和阻尼矩阵，要用来形成整体的刚度矩阵、质量矩阵和阻尼矩阵。在不考虑体积力的条件下，整个结构的动力学方程为： 当f = 0, c = 0时得到自由振动时的无阻尼动力方程 常系数线性齐次常微分方程组，令其解为： 四、方程的特征值及振型的正交性 代入方程得到： 得到齐次的线性代数方程组 N阶自由度系统的自由振动方程应有n个固有频率i=1,2,3,…n 全部频率从小到大按顺序排列起来得到频率向量 四、方程的特征值及振型的正交性 四、方程的特征值及振型的正交性 第n个规格化振型形式 按照 的大小顺序排列起来就得到振型矩阵。 四、方程的特征值及振型的正交性 关于振型的正交性 四、方程的特征值及振型的正交性 关于振型的正交性 五、结构的动力学响应求解方法 振型叠加法 将n阶自由度系统的动力方程，经过振型模态矩阵变换，转化为互相不耦合的n个单自由度问题，进行逐个求解后，然后再叠加得到动力响应的结果。 n个自由度的结构，在激励力p(t)的作用下的动力响应可以表示为各阶主振型的线性叠加。 五、结构的动力学方程求解方法 逐步积分法（直接积分法） 对于有复杂激振力或非比例阻尼情况下，将时间离散化，逐步求出每个时间间隔 上的状态向量（位移、速度和加速度），最后获得的状态向量为结构系统的动力响应解。 六、有限元法求解动力学问题实例 梁直径d=10mm,长度2000mm，材料为钢，E=2.1X105MPa，密度7.8X103kg/m3。采用有限元法计算梁的固有频率。 将梁分成4个单元，单元彼此之间以节点相连 建立单元的刚度矩阵和质量矩阵 根据相邻单元在连接节点在具有相同