机械振数值分析.ppt

机械振动数值分析;参考书目 有限单元法，王勖成，清华大学出版社 计算动力学，王天舒，清华大学出版社 有限元分析的概念与应用，Robet D.Cook，西安交通大学出版社 ;求解域的离散 —— 将连续体求解域离散为若干个子域，通过子域边界协调条件构造连续体； 利用近似函数求解 —— 每个子域内用近似函数分片表示全域内待定的未知变量； 原问题的等效 —— 利用变分或加权余量法，建立基本变量代数方程组或常微分方程组。;复杂几何构型的适应性； 各种物理问题的可应用性； 建立于严格理论基础上的可靠性； 适合计算机实现的高效性。;单元类型和形式； 有限元法的理论基础和离散格式； 有限元方程的求解方法； 有限元分析程序开发。;标准化 —— 任意复杂问题 → 模块化分解，单元建模 → 有限种类模块化单元 规范化 —— 几何建模 → 力学建模 → 求解 → 后处理分析 通用化 —— 形成标准模块化程序 应用规模化、普及性 —— 求解问题规模庞大，易于为工程技术人员掌握;连续性假设 —— 变形体内部处处连续 均匀性假设 —— 变形体内部物质分配均匀 各向同性假设 —— 物质在各方向上特性相同 线弹性假设 —— 变形与外力作用的关系为线性 小变形假设 —— 变形量远小于物体本身尺寸;亨掺厂护遣桨肥诈脾篓筹敞淀臆补娇椿寄蛹勉瘪店咸浩谬汰饵飞三魏皋挚机械振动数值分析机械振动数值分析;运动微分方程;边界条件;变形体域内任意一点在任意时刻均满足运动微分方程。 变形体边界上任意一点在任意时刻均满足边界条件。;残余力方程;当上式对任意权函数均满足，则称为式（7）为微分方程（1）和边界条件（4）的等效积分形式。;权函数可以选N个函数的线性组合，即;配点法 —— 取Dirac函数为权函数 子域法 —— 权函数在N个子域内取1，在子域外取零，即 ;最小二乘法 —— 调整近似函数中的参数，使余量均方和最小，即 ;迦辽金法 —— 取试探函数为权函数 迦辽金法的特点 ;例1 用各种加权余量法计算图示弹性基础梁的挠度;弹性基础梁的基本微分方程和边界条件为;取试探解为无弹性基础时的精确解;残差方程可写为; 达朗伯—拉格朗日原理; 阶连续性函数; 哈密顿原理;将式（21）代入式（20），得普遍意义下的哈密顿原理;考虑粘滞力后;考虑到应力应变关系，外力虚功可改写为;例2 用哈密顿原理推导受均布动载荷的等截面悬臂梁的振动微分方程;系统总势能可表示为;对系统势能取变分;代入哈密顿方程得;拉格朗日乘子法; 拉格朗日乘子法修正泛函的变分;由此得; 罚函数法修正泛函的变分;罚函数法的特点 附加条件近似满足； 不增加方程阶数； 罚参数相当于刚度系数； 解的精度与罚参数有关。;例3 如图所示弹簧系统，在第一个自由度上施加外力，在第二个自由度上施加强制位移。用拉格朗日乘子法和罚函数法计算第一个自由度的位移;将强制位移边界代入，可解得精确解;可解得; 广义变分原理; 广义变分原理; 广义变分原理;习题2 如图所示弹簧系统，在第一个和第三个自由度上施加外力，在第二个自由度上施加强制位移。用拉格朗日乘子法和罚函数法计算第一个和第三个自由度的位移，并计算第二个自由度上的作用力;几何关系 —— 应力应变关系 —— ; 杆单元的刚度矩阵和质量矩阵;将式（2）代入式（1），得;记;2结点杆单元的刚度矩阵和质量矩阵如式（10）所示;考虑在平面中变形的杆单元，如图所示;则单元刚度矩阵、质量矩阵和载荷向量转换到总体坐标系中为;考虑如图所示杆系结构;总体坐标系下的单元刚度矩阵和单元质量矩阵为; 平面杆单元例子; 有限元离散的一般流程; 形函数的确定; 形函数的确定; 形函数的确定; 形函数的确定; 单元内部积分; 单元内部积分; 单元内部积分; 单元内部积分; 约束条件的处理; 约束条件的处理; 约束条件的处理; 等参变换;等参变换;等参变换;等参变换后的单元刚度矩阵、单元质量矩阵格式;第三章 大型系统特征值问题; 系统运动方程; 特征向量正交性; 特征向量正交性;刚度矩阵为半正定矩阵 质量矩阵为半正定矩阵 ; Rayleigh商; Rayleigh商; 误差估计; 误差估计; 误差估计; 逆迭代法; 逆迭代法; 逆迭代法; 逆迭代法; 逆迭代法; 逆迭代法; 逆迭代法; 带移轴的逆迭代法; Rayleigh商迭代法; 变换法; 变换法; 雅可比变换法; 雅可比变换法; 雅可比变换法; 雅可比变换法; 雅可比变换法; 雅可比变换法; 雅可比变换法; 雅可比变换法; 广义雅可比变换法; 广义雅可比变换法; 广义雅可比变换法; 广义雅可比变换法; 广义雅可比变换法;