FEM结构动力学问题有限元重点分析.ppt

瞬态分析- 术语和概念求解方法 系统的动力响应（振型叠加法） 系统的动力响应（振型叠加法） 系统的动力响应（振型叠加法） 系统的动力响应（振型叠加法） 系统的动力响应（振型叠加法） 在振型叠加法中，将系统位移转换到以固有振型为基向量的空间，这对系统的性质并无影响，而是以求解广义特征值为代价，得到n个单自由度的运动方程。 振型叠加法中对于n个自由度系统运动方程的积分，比对联立方程组的直接积分节省计算时间。需要注意的是(1) 有限元高频成分对于系统的实际影响较小，所以通常求前p个特征解的响应；(2) 有限元高阶特征解与实际情况相差较大，因此，用有限元求高阶特征解的意义不大，此时，需要采用其它方法。 对于非线性系统通常采用直接积分法，此时，无法采用振型叠加法。 对于一般的n个自由度系统，如果只积分p个单自由度系统的运动方程，即使积分是精确的，最后得到的系统响应也会由于忽略高阶振型而引入误差。 系统的动力响应 系统的动力响应(cont.) 中心差分法是显式算法； 中心差分法是条件稳定算法； 中心差分法比较适合冲击、爆炸等载荷引起的波传播问题； 不适合结构动力学问题。结构的动力响应通常低频成分是主要的，从计算精度考虑，允许采用较大的时间步长。此时需要采用无条件稳定的隐式算法，如Newmark算法。 系统的动力响应(cont.) Newmark计算结果 中心差分法 Newmark方法 Newmark方法求解步骤 Newmark方法是隐式算法。 当d≥0.5,a≥0.25(0.54+d)2时，算法是无条件稳定的。 例题： (1) 中心差分法 解： 中心差分计算结果 (2) Newmark方法 * 有限元方法与应用结构动力学问题有限元 结构动力学：研究结构在动力荷载作用下的动力反应。 地震荷载 风荷载：Tacoma大桥风毁录像 动力荷载：荷载的大小、方向、作用位置 随时间而变化。 输入 input 输出 Output 结构体系 静力响应 静荷载 位移 内力 应力 刚度、约束杆件尺寸 截面特性 大小 方向 作用点 结构体系 动力响应 输入 input 输出 Output 动荷载 动位移 加速度 速度 动应力 动力系数 随时间变化 质量、刚度阻尼、约束频率、振型 大小 方向 作用点 时间变化 数值 时间函数 结构动力体系 动载荷（又称动力分析） 固有特性分析 响应分析 固 有 频 率 振 型 位 移 响 应 速 度 响 应 加 速 度 响 应 动 应 变 动 应 力 固有特性：是一组模态参数构成，它由结构本身（质量与刚度分布）决定，而与外部载荷无关，但决定了结构对动载荷的响应； 响应分析：是计算结构对给定动载荷的各种响应特性。 模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。 固有特性分析就是对模态参数进行计算，其目的一是避免结构出现共振和有害的振型，二是为响应分析提供必要依据。 结构动力学问题的有限元法的实质就是将一个弹性连续体的振动问题，离散为一个以有限个节点位移为广义坐标的多自由度系统的振动问题。其基本原理和分析方法类同静力学的有限元法，按杆梁、薄板等不同结构进行分析。不同的是，应用振动理论建立动力学方程时，在单元分析中除需形成刚度矩阵外，还需形成质量矩阵，阻尼矩阵；在整体分析中，不仅求动力响应，还有求解特征值问题（结构振动的固有频率及相应的振动型（或模态）） 从以上步骤可以看出，和静力分析相比，在动力分析中，由于惯性力和阻尼力出现在平衡方程中，因此引入了质量矩阵和阻尼矩 阵，最后得到求解方程不是代数方程组，而是常微分方程组。其 它的计算步骤和静力分析是完全相同的。 关于二阶常微分方程组的解法有两类：直接积分法和振型叠加法。 直接积分法是直接对运动方程积分。而振型叠加法是首先求解一无阻尼的自由振动方程，然后用解得的特征向量，即固有振型对 运动方程式进行变换。 动力学有限元法的特点 一、载荷特点 结构所受的载荷是随时间变化的动载荷。 这是与静力分析的一个根本区别。 二、位移特点 1、节点位移{q}不仅是坐标的函数，而且也是时间的函数。仍以节点位移{q}作为基本未知量。 2、节点具有速度 加速度。 3、利用节点位移插值表示单元内任一点的位移,一般仍采用与静力分析相同的形函数 [N]。当单元数量较多时，上述插值可以得到较好的插值精度。 4、在线弹性条件下，单元内的应变和应力与节点