1、M1-1,目录,第一章: 动力学绪论 M1-4第一节 动力学分析概述 M1-6 第二节 动力学分析类型 M1-9第三节 基本概念和术语 M1-14第四节 动力学分析实例 M1-33第二章: 模态分析 M2-1第一节 模态分析概述 M2-3第二节 模态分析术语和概念 M2-5第三节 模态分析步骤 M2-14第四节 模态分析实例 M2-33第五节 有预应力的模态分析 M2-34第六节 模态的循环对称性问题 M2-42第三章: 谐分析 M3-1第一节 谐分析概述 M3-3第二节术语和概念 M3-6第三节谐分析步骤 M3-12第四节谐分析实例 M3-33,M1-2,目录 （接上页）,第四章: 瞬态动力
2、学分析 M4-1第一节 瞬态动力学分析概述 M4-3第二节 瞬态动力学分析术语和概念 M4-5第三节 瞬态动力学分析步骤 M4-15第四节 瞬态动力学分析实例 M4-49第五章: 谱分析 M5-1第一节 谱分析概述 M5-3第二节 响应谱分析 M5-5第三节 随即震动分析 M5-28第六章: 模态叠加 M6-1第一节 什么是模态叠加? M6-3第二节 模态叠加步骤 M6-5第三节 模态叠加实例 M6-23,M1-3,第一章,动力学绪论,M1-4,第一章: 动力学绪论,第一节: 动力学分析的定义和目的 第二节: 动力学分析的不同类型 第三节: 基本概念和术语 第四节: 动力学分析的一个实例,M1
3、-5,动力学 第一节: 定义和目的,什么是动力学分析? 动力学分析是用来确定惯性（质量效应）和阻尼起着重要作用时结构或构件动力学特性的技术。 “动力学特性” 可能指的是下面的一种或几种类型: 振动特性 - （结构振动方式和振动频率） 随时间变化载荷的效应（例如：对结构位移和应力的效应） 周期（振动）或随机载荷的效应,M1-6,动力学 定义和目的（接上页）,静力分析也许能确保一个结构可以承受稳定载荷的条件，但这些还远远不够，尤其在载荷随时间变化时更是如此。 著名的美国塔科马海峡吊桥（Galloping Gertie） 在 1940年11月7日，也就是在它刚建成4个月后，受到风速为42英里/小时的
4、平稳载荷时发生了倒塌。,M1-7,动力学 定义和目的（接上页）,动力学分析通常分析下列物理现象: 振动 - 如由于旋转机械引起的振动 冲击 - 如汽车碰撞，锤击 交变作用力 - 如各种曲轴以及其它回转机械等 地震载荷 - 如地震，冲击波等 随机振动 - 如火箭发射，道路运输等上述每一种情况都由一个特定的动力学分析类型来处理,M1-8,动力学 第二节: 动力学分析类型,请看下面的一些例子: 在工作中，汽车尾气排气管装配体的固有频率与发动机的固有频率相同时，就可能会被震散。那么，怎样才能避免这种结果呢? 受应力（或离心力）作用的涡轮叶片会表现出不同的动力学特性，如何解释这种现象呢?答案：进行 模态
5、分析 来确定结构的振动特性,M1-9,动力学 动力学分析类型（接上页）,汽车防撞挡板应能承受得住低速冲击 一个网球排框架应该设计得能承受网球的冲击，但会稍稍发生弯曲解决办法 ：进行 瞬态动力学分析 来计算结构对随时间变化载荷的响应,M1-10,动力学 动力学分析类型（接上页）,回转机器对轴承和支撑结构施加稳态的、交变的作用力，这些作用力随着旋转速度的不同会引起不同的偏转和应力解决办法 ： 进行谐分析来确定结构对稳态简谐载荷的响应,M1-11,位于地震多发区的房屋框架和桥梁应该设计应当能够承受地震载荷要求.解决办法：进行谱分析来确定结构对地震载荷 的影响,Courtesy: US Geologi
6、cal Survey,动力学 动力学分析类型（接上页）,M1-12,太空船和飞机的部件必须能够承受持续一段时间的变频率随机载荷。解决办法 ：进行随机振动分析来确定结构对随机震动的影响,Courtesy: NASA,动力学 动力学分析类型（接上页）,M1-13,动力学 第三节: 基本概念和术语,讨论的问题: 通用运动方程 求解方法 建模要考虑的因素 质量矩阵 阻尼,M1-14,动力学 - 基本概念和术语 运动方程,通用运动方程如下：,不同分析类型是对这个方程的不同形式进行求解 模态分析：设定F（t）为零 ，而矩阵 C 通常被忽略； 谐响应分析：假设F（t） 和 u（t） 都为谐函数，例如 Xsi
7、n（wt），其中，X 是振幅， w 是单位为弧度/秒的频率； 瞬间动态分析：方程保持上述的形式。,M1-15,动力学 - 基本概念和术语 运动方程(接上页),其中: M = 结构质量矩阵 C = 结构阻尼矩阵 K = 结构刚度矩阵 F = 随时间变化的载荷函数 u = 节点位移矢量 = 节点速度矢量 = 节点加速度矢量,M1-16,动力学 -基本概念和术语 求解方法,如何求解通用运动方程 ? 两种主要方法： 模态叠加法 直接积分法模态叠加法 按自然频率和模态将完全耦合的通用运动方程转化为一组独立的非耦合方程 可以用来处理瞬态动力学分析和谐响应分析 详见第六章,M1-17,动力学 - 基本概念和
8、术语 求解方法 （接上页）,直接积分法 直接求解运动方程 在谐响应分析中，因为载荷和响应都假定为谐函数，所以运动方程是以干扰力频率的函数而不是时间的函数的形式写出并求解的 对于瞬态动力学，运动方程保持为时间的函数，并且可以通过显式或隐式的方法求解,M1-18,动力学 - 基本概念和术语 求解方法 （接上页）,显式求解方法 也称为闭式求解法或预测求解法 积分时间步 Dt 必须很小，但求解速度很快（没有收敛问题） 可用于波的传播，冲击载荷和高度非线性问题 ANSYS-LS/DYNA 就是使用这种方法，此处不作介绍,显式求解法 也可成为开式求解法或修正求解法 积分时间步 Dt 可以较大，但方程求解时
9、间较长（因为有收敛问题） 除了 Dt 必须很小的问题以外，对大多数问题都是有效的 ANSYS 使用 Newmark 时间积分方法,M1-19,动力学 - 基本概念和术语 求解方法 （接上页）,显式方法 当前时间点的位移 ut 由包含时间点t-1 的方程推导出来 有条件稳定: 如果Dt 超过结构最小周期的确定百分数,计算位移和速度将无限增加隐式方法 当前时间点的位移 ut 由包含时间点 t 的方程推导出来 无条件稳定: Dt的大小仅仅受精度条件控制, 无稳定性。,M1-20,动力学 - 基本概念和术语 建模要考虑的问题,几何形状和网格划分 材料性质 各种非线性几何形状和网格划分: 一般同于静态分
10、析要考虑的问题 要包括能充分描绘模型几何形状所必须的详细资料 在关心应力结果的区域应进行详细的网格划分，在仅关心位移结果的时候，粗糙的网格划分可能就足够了,M1-21,动力学 - 基本概念和术语 建模要考虑的问题（接上页）,材料性质: 需要定义杨氏模量和密度 请记住要使用一致的单位 当使用英制单位时，对于密度，要定义质量密度而不是重力密度： 质量密度=重力密度（lb/in3） / g （in/sec2） 钢的密度 = 0283/386 = 73 x 10-4 lb-sec2/in4,M1-22,动力学 - 基本概念和术语 建模要考虑的问题（接上页）,非线性 （大变形，接触，塑性等等）： 仅在完
11、全瞬态动力学分析中允许使用。 在所有其它动力学类型中（如模态分析、谐波分析、谱分析以及简化的模态叠加瞬态分析等） ，非线性问题均被忽略，也就是说最初的非线性状态将在整个非线性求解过程中一直保持不变。,M1-23,1,2,BEAM3,动力学 - 基本概念和术语 质量矩阵,对于动力学分析需要质量矩阵 M，并且这个质量矩阵是按每个单元的密度以单元计算出来的。有两种类型的质量矩阵 M: 一致质量矩阵 和集中质量矩阵， 对于2-D 梁单元BEAM3，其质量分布矩阵和集中质量矩阵如下所示：,M1-24,动力学 - 基本概念和术语 质量矩阵（接上页）,一致质量矩阵 通过单元形函数计算出来； 是大多数单元的缺
12、省选项； 某些单元有一种称为简化质量矩阵 的特殊形式的质量矩阵，其中对应于转动自由度的各元素均被置零。集中质量矩阵 质量被单元各节点所平分，非对角线元素均为零； 通过分析选项来激活。,M1-25,动力学 - 基本概念和术语 质量矩阵（接上页）,应当采用哪种质量矩阵? 对大多数分析来说，一致质量矩阵为缺省设定； 若结构在一个方向的尺寸与另两个方向相比很小时，可采用简化质量矩阵（如果可能得到的话）或集中质量矩阵例如细长的梁或很薄的壳； 集中质量矩阵可用于波的传播问题。,M1-26,动力学 - 基本概念和术语 阻尼,什么是阻尼? 阻尼是一种能量耗散机制，它使振动随时间减弱并最终停止 阻尼的数值主要取
13、决于材料、运动速度和振动频率 阻尼可分类如下： 粘性阻尼 滞后或固体阻尼 库仑或干摩擦阻尼,M1-27,动力学 - 基本概念和术语 阻尼（接上页）,粘性阻尼 粘性阻尼一般物体在液体中运动时发生 由于阻尼力与速度成正比，因此在动力学分析中要考虑粘性阻尼 比例常数 c 称作阻尼常数 通常用 阻尼比 x （阻尼常数 c 对临界阻尼常数 cc*的比值）来量化表示 临界阻尼定义为出现振荡和非振荡行为之间的阻尼的极值, 此处阻尼比 = 10*对一个质量为 m ，频率为 w的单自由度弹簧质量系统, cc = 2mw,注意: 阻尼比 x = 对于螺栓或铆钉链接结构为2%到15%,M1-28,动力学 - 基本概
14、念和术语 阻尼（接上页）,滞后和固体阻尼 是材料的固有特性 在动力学分析中应该考虑 认识还不是很透彻，因此很难定量的确定库仑或干摩擦阻尼 物体在干表面上滑动时产生的阻尼 阻尼力与垂直于表面的力成正比 比例常数 m 就是摩擦系数 动力学分析中一般不予考虑,M1-29,动力学 - 基本概念和术语 阻尼（接上页）,ANSYS 允许上述所有三种形式的阻尼 通过规定阻尼比x, Rayleigh阻尼常数 a （后面将进行讨论），或定义带有阻尼矩阵的单元，可将粘性阻尼纳入考虑 通过规定另一种Rayleigh 阻尼常数 b （后面将进行讨论）可将滞后或固体阻尼纳入考虑通过规定带有摩擦性能的接触表面单元和间隙单
15、元，可将库仑阻尼纳入考虑，（此处不进行讨论，可参见ANSYS 结构分析指南）,M1-30,动力学 - 基本概念和术语 阻尼（接上页）,Rayleigh 阻尼常数a 和 b 用作矩阵 M 和 K 的乘子来计算 C:C = aM + bKa/2w + bw/2 = x此处 w 是频率, x 是阻尼比 在不能定义阻尼比 x时，需使用这两个阻尼常数 a 是粘度阻尼分量， b 是滞后或固体或刚度阻尼分量,M1-31,动力学 - 基本概念和术语 阻尼（接上页）,a 阻尼 亦可称作质量阻尼 只有当粘度阻尼是主要因素时才规定此值，如在进行各种水下物体、减震器或承受风阻力物体的分析时 如果忽略b 阻尼，a 可通
16、过已知值x（阻尼比） 和已知频率w来计算： a = 2xw 因为只允许有一个a值,所以要选用最主要的响应频率来计算 a,a=3,1,2,05,M1-32,动力学 - 基本概念和术语 阻尼（接上页）,b阻尼 亦可称作结构或刚度阻尼 是大多数材料的固有特性 b阻尼对每一个材料进行规定（作为材料性质DAMP），或作为一个单一的总值 如果忽略a 阻尼， b可以通过已知的x（阻尼比）和已知频率w来计算： b = 2x/w选用最主要的响应频率来计算b,b=0004,0003,0001,0002,M1-33,动力学 - 基本概念和术语 阻尼（接上页）,定义a 和 b 阻尼： 使用方程a/2w + bw/2 = x 因为有两个未知数，所以近似的假设alpha 和beta 阻尼的总和在频率范围w1 至w2 之间是一个长阻尼比x 这将给出两个联立方程，从而可以计算出a 和 bx = a/2w1 + bw1/2x = a/2w2 + bw2/2,a+b,b,a,w1,w2,M1-34,动力学 第四节：实例介绍,在实例中，你可运行 “Galloping Gertie” （塔可马吊桥）的动力学分析实例只须遵循动力学实例附刊中的说明主要目的是向初学者介绍典型动力学分析的步骤，每一步具体含义参见本指南的后面的介绍资料。,
