机械振动数值分析.ppt

机械振动数值分析 主讲：何 欢 博士 工作单位：振动工程研究所 电话 办公室：18楼602 Email: hehuan@nuaa.edu.cn 求解域的离散 —— 将连续体求解域离散为若干个子域，通过子域边界协调条件构造连续体； 利用近似函数求解 —— 每个子域内用近似函数分片表示全域内待定的未知变量； 原问题的等效 —— 利用变分或加权余量法，建立基本变量代数方程组或常微分方程组。 复杂几何构型的适应性； 各种物理问题的可应用性； 建立于严格理论基础上的可靠性； 适合计算机实现的高效性。 单元类型和形式； 有限元法的理论基础和离散格式； 有限元方程的求解方法； 有限元分析程序开发。 标准化 —— 任意复杂问题 → 模块化分解，单元建模 → 有限种类模块化单元 规范化 —— 几何建模 → 力学建模 → 求解 → 后处理分析 通用化 —— 形成标准模块化程序 应用规模化、普及性 —— 求解问题规模庞大，易于为工程技术人员掌握 连续性假设 —— 变形体内部处处连续 均匀性假设 —— 变形体内部物质分配均匀 各向同性假设 —— 物质在各方向上特性相同 线弹性假设 —— 变形与外力作用的关系为线性 小变形假设 —— 变形量远小于物体本身尺寸 变形体域内任意一点在任意时刻均满足运动微分方程。 变形体边界上任意一点在任意时刻均满足边界条件。 配点法 —— 取Dirac函数为权函数 子域法 —— 权函数在N个子域内取1，在子域外取零，即 最小二乘法 —— 调整近似函数中的参数，使余量均方和最小，即 迦辽金法 —— 取试探函数为权函数 迦辽金法的特点 拉格朗日乘子法 几何关系 —— 应力应变关系 —— 几何体的离散逼近； 自由度选择与形函数的确定； 单元应变场的表示； 单元应力场的表示； 哈密顿原理或达朗伯-拉格朗日原理； 系统泛函变分； 单元内部积分； 单元坐标系到总体坐标系的转换； 总体刚度矩阵和质量矩阵的组装； 约束条件的处理或引入。 形函数的阶数与单元形式和自由度数有关 结点上的插值函数满足 单元中任意点处形函数和应有 选取多项式时，常数项和一次项要完备 为表达统一起见，可采用参数变换 单点积分 两点积分 直接代入法 特点 方程阶数降低 对角元素置1法 特点 适用于零位移情况 对角元素置大数法 特点 相当于罚函数法 局部坐标系向物理坐标系的映射 类似于构造形函数，可构造物理坐标和局部坐标之 间的关系。仍然以三角形平面单元为例 需要注意的地方 对面积分来说 刚度矩阵为半正定矩阵 质量矩阵为半正定矩阵 Rayleigh商的上下界 物理意义 —— 假设模态相当于增加了约束，使得系统刚度提高 基本假定 刚度矩阵非奇异 迭代方程 考虑到任意向量都可由特征向量表示为 综合式（26）和式（27），得 从式（30）中容易看出，当 考虑Rayleigh商 结合式（27）、（28）和（29）可将式（33）化为 逆迭代法的一般步骤 给定初始迭代向量； 若已计算出前i-1阶模态，根据已得到的特征向量进行正交化处理； 逆迭代法的一般步骤 解方程； 计算特征值 几个要点 正交化处理； 归一化处理； 当刚度矩阵奇异时的处理； 系统包含近频时的处理 考虑特征值问题 Rayleigh商迭代法 —— 以Rayleigh商作为移轴量 两类应用广泛的变换法 雅可比变换法 Householder-QR法 中心差分法 将位移函数按Taylor级数展开，得 由式（13）可得 中心差分法 将式（14）代入任意时刻的运动方程 得 其中 中心差分法 起始步的计算 定义积分常数 中心差分法 计算步骤 1、初始计算 1） 形成刚度矩阵 和质量矩阵 和阻尼矩阵 ； 2） 给定初始条件 和 ，并计算 ； 3） 给定时间步长，计算积分常数； 4） 计算起步位移： ； 5） 计算有效质量矩阵： ； 6） 对有效质量矩阵进行三角分解： ； 中心差分法 2、对每一个时间步 1） 计算时刻 的有效载荷 ； 2） 计算 时刻的位移： ； 3） 如有需要，计算时刻