ANSYS NCORPORATED ANSYS结构分析指南(下) 结构动力学 v0205
ANSYS动力学分析指南 谱分析指南 第一章模态分析 1.1模态分析的定义及其应用 ….-1 1.2模态分析中用到的命令 .1-1 1.3模态提取方法 ..1-1 1.4矩阵缩减技术和主自由度选择准则 .1-5 1.5模态分析过程 ..1-7 1.6建 .1-7 1.7加载及求解 .1-7 1.8扩展模 ..-12 1.9观察结……… .1-14 1.10有预应力模态分析 .1-16 1.11大变形预应力模态分析 .1-17 1.12循环对称结构的模态分析 .1-17 1.13模态分析实例 1-22 第二章谐响应分析 2.1谐响应分析的定义与应用 .2-1 2.2谐响应分析中用到的命令 2-1 2.3三种求解方法 ..2-1 析¨ 2.4完全法谐响应分析 ..2-3 2.5缩减法谐响应分析 2-11 2.6模态叠加法谐响应分析 ….-6 2.7有预应力的完全法谐响应分………….-19 2.8谐响应分析实例 ….-20 第三章瞬态动力学分析 3.1瞬态动力学分析的定义 .3-1 3.2学习瞬态动力学的预备工作 ..3-1 3.3三种求解方法… 3-2 3.4完全法瞬态动力学分析 ..3-3 3.5模态叠加法瞬态动力分析 .3-15 3.6缩减法瞬态动力学分析过程 ..3-21 3.7有预应力瞬态动力学分析 .3-28 3.8瞬态分析的关键技术细节 .3-29 3.9瞬态动力学分析实例 .3-33 第四章谱分析 4.1谱分析的定义 .- 4.2什么是谱 4-1 4.3谱分析使用的命令 4.4单点响应谱(SPRS)分析步骤 .4-3 4.5随机振动(PSD)分析步骤 …….- 第1页
NSYS动力学分析指南 谱分析指南 4.6随机振动分析结果应用 4-17 4.7DDAM(动力设计分析方法)谱分析. 4-20 4.8多点响应谱(MPRS)分析 4-20 4.9谱分析的实例(GUI命令流和批处理) 4-21
ANSYS动力学分析指南 模态分析 第一章模态分析 §1.1模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振 型,它们是承受动态载荷结构设计中的重要参数。同时,也可以作为其它动力学分析问题的起点 例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分 析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转 的涡轮叶片等的模态分析,后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态 分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性,如塑性和接触(间隙)单元, 即使定义了也将被忽略。 ANSYS提供了七种模态提取方法,它们分别是子空间法、分块 Lanczos法 PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存 在阻尼。后面将详细介绍模态提取方法。 §1.2模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样,无论进行何种类型的分 析,均可从用户图形界面(GUI)上选择等效于命令的菜单选项来建模和求解问题 后面的“模态分析实例(命令流或批处理方式)”将给出进行该实例模态分析时要输入的命令(手 工或以批处理方式运行 ANSYS时)。而“模态分析实例(GUI方式)”则给出了以从 ANSYS GU 中选择菜单选项方式进行同一实例分析的步骤。(要想了解如何使用命令和GUI选项建模,请参阅 )。>中有更详细的按字母顺序列出的 ANSYS 命令说明。 §1.3模态提取方法 典型的无阻尼模态分析求解的基本方程是经典的特征值问题 K{o,}=aM{} 其中: [K]=刚度矩阵, }=第i阶模态的振型向量(特征向量) O2=第i阶模态的固有频率(ao2是特征值) 质量矩阵 有许多数值方法可用于求解上面的方程。 ANSYS提供了7种方法模态提取方法,下面分别进行 讨论。 第1页
ANSYS动力学分析指南 模态分析 1.分块 Lanczos法 2.子空间( Subspace)法 3. Power Dynamics法 4.缩减( Reduced/ Householder)法 5.非对称( Unsymmetric)法 6.阻尼(Damp)法(阻尼法求解的是另一个方程,参见>中关于此法的 详细信息) 7.QR阻尼法(QR阻尼法求解的是另一个方程,参见>中关于此法的详 细信息) 注意一阻尼法和对称法在 ANSYS/Professional中不可用 前四种方法(分块 Lanczos法、子空间法、 PowerDynamics法和缩减法)是最常用的模态提取 方法。下表比较了这四种模态提取方法,并分别对每一种方法进行了简要描述 对称系统特征值求解法表 模态提取法 适用范围 内存要求存贮要求 缺省提取方法 分块 用于提取大模型的多阶模态(40阶以上) Lanczos法|建议在模型中包含形状较差的实体和壳单元时采用此法 最适合于由壳或壳与实体组成的模型 速度快,但要求比子空间法内存多50% 用于提取大模型的少数阶模态(40阶以下 子空间法适合于较好的实体及壳单元组成的模型 可用内存有限时该法运行良好 用于提取大模型的少数阶模态(20阶以下 Power适合于100K以上自由度模型的特征值快速求解 Dynamics对于网格较粗的模型只能得到频率近似值 复频情况时可能遗漏模态 用于提取小到中等模型(小于10K自由度)的所有模态 缩减法 选取合适主自由度时可获取大模型的少数阶(40阶以下)模态, 此时频率计算的精度取决于主自由度的选取。 §1.3.1分块 Lanczos法 分块 Lanczos法特征值求解器是却省求解器,它采用 Lanczos算法,是用一组向量来实现 Lanczos 递归计算。这种方法和子空间法一样精确,但速度更快。无论 EQSLV命令指定过何种求解器进行求 解,分块 Lanczos法都将自动采用稀疏矩阵方程求解器。 计算某系统特征值谱所包含一定范围的固有频率时,采用分块 Lanczos法方法提取模态特别有 效。计算时,求解从频率谱中间位置到髙频端范围内的固有频率时的求解收敛速度和求解低阶频率 时基本上一样快。因此,当采用频移频率( FREQB)来提取从 FREQB(起始频率)的n阶模态时 该法提取大于 FREQB的n阶模态和提取n阶低频模态的速度基本相同。 第2页
ANSYS动力学分析指南 模态分析 §1.3.2子空间法 子空间法使用子空间迭代技术,它内部使用广义 Jacobi迭代算法。由于该方法采用完整的] 和[M]矩阵,因此精度很高,但是计算速度比缩减法慢。这种方法经常用于对计算精度要求高,但 无法选择主自由度(DOF)的情形 做模态分析时如果模型包含大量的约束方程,使用子空间法提取模态应当采用波前( front)求 解器,不要采用JCG求解器;或者是使用分块 Lanczos法提取模态。当你的分析中存在大量的约束 方程时,如果采用JCG求解器组集内部单元刚度,致使计算要求有很大的内存才能进行下去 §1.3.3 PowerDynamics法 Power Dynamics法内部采用子空间迭代计算,但采用PCG迭代求解器。这种方法明显地比子空 间法和分块 Lanczos法快。但是,如果模型中包含形状较差的单元或病态矩阵时可能出现不收敛问 题。该法特别适用于求解超大模型(大于10000°个自由度)的起始少数阶模态。谱分析不要使用 该方法提取模态。 PowerDynamics法不进行Sum序列检查(即不检查模态遗漏问题),这可能影响有多个重复频 率问题的解。此法总是采用集中质量近似算法,即自动釆用集中质量矩阵( LUMPM,ON。 注意如果用 Powerdynamics法求解含刚体运动的模型的模态,则一定要用RGD命令或选择等效 的GUI途径 yEE-Main Menu> Solution Analysis Options E Main Menu> Preprocessor>-Loads-> Analysis Options )o §1.3.4缩减法 缩减法采用HBI算法( Householder-二分逆迭代)来计算特征值和特征向量。由于该方法采用 个较小的自由度子集即主自由度(DOF)来计算,因此计算速度更快。主自由度(DOF)导致计 算过程中会形成精确的[K]矩阵和近似的M]矩阵(通常会有一些质量损失)。因此,计算结果的精 度将取决于质量阵[M]的近似程度,近似程度又取决于主自由度的数目和位置 §1.3.5非对称法 非对称法也采用完整的[k]和[M]矩阵,适用于刚度和质量矩阵为非对称的问题(例如声学中 流体结构耦合问题)。此法采用 Lanczos算法,如果系统是非保守的(例如轴安装在轴承上),这种 算法将解得复数特征值和特征向量。特征值的实部表示固有频率,虚部是系统稳定性的量度一负值 表示系统是稳定的,而正值表示系统是不稳定的。该方法不进行Sum序列检查,因此有可能遗漏 些高频端模态。 第3页
ANSYS动力学分析指南 模态分析 §1.3.6阻尼法 阻尼法用于阻尼不能被忽略的问题,如转子动力学研究。该法使用完整矩阵([K]、[M]及阻 阵[C])。阻尼法采用 lanczos算法并计算得到复数特征值和特征向量(如下所述)。此法不能用 Sum序列检査。因此,有可能遗漏所提取频率的一些高频端模态。 §1.3.5.1阻尼法特征值的实部和虚部 特征值的虚部代表系统的稳态角频率。特征值的实部σ代表系统的稳定性。如果σ小于零, 系统的位移幅度将按EXP(σ)指数规律递减。如果σ大于零,位移幅度将按指数规律递増。(或者换 句话说,负的σ表示按指数规律递减的稳定响应;正的σ则表示按指数规律递增的不稳定响应。) 如果不存在阻尼,特征值的实部将为零 ANSYS报告的特征值结果实际上是被2丌除过的。这样给出的频率是以H(周/秒)为单位的。 报告的特征值虚部=0)(27) 报告的特征值实部=/(2) §1.3.5.2阻尼法特征向量的实部和虚部 在有阻尼系统中,不同节点上的响应可能存在相位差。对任何节点,幅值应是特征向量实部和 虚部分量的矢量和。 §1.3.7QR阻尼法 QR阻尼法同时具有分块 Lanczos法与复 Hessenberg法的优点,最关键的思想是,以线性合并无 阻尼系统少量数目的特征向量近似表示前几阶复阻尼特征值。采用实特征值求解(分块 Lanczos法) 无阻尼振型之后,运动方程将转化到模态坐标系。然后,采用QR阻尼法,一个相对较小的特征值 问题就可以在特征子空间中求解出来了 该方法能够很好地求解大阻尼系统模态解,阻尼可以是任意阻尼类型,即无论是比例阻尼或非 比例阻尼。由于该方法的计算精度取决于提取的模态数目,所以建议提取足够多的基频模态,特别 是阻尼较大的系统更应当如此,这样才能保证得到好的计算结果。该方法不建议用于提取临界阻尼 或过阻尼系统的模态。该方法输出实部和虚部特征值(频率),但仅仅输岀实特征向量(模态振型)。 参见CE方法的详细内容,掌握使用QR阻尼法 MODOPT命令)处理约束方程(E)的技术 约束方程(CE)方法 约束方程 Cakey 处理方法 应用范围 直接消去法模型中只有少量约束方程时使用。例如,在一个1000由度问 题中,只有大约1,000个约束方程。一旦约束方程太多,该方法需 要的内存极高。此时,建议使用拉格朗日乘子法( Cakey=1或2)。 拉格朗日乘模型中存在大量约束方程时使用。例如,在一个1000度问 子法 题中,具有1,000以上的约束方程。特别注意,当使用 CEINTE CERIG或 CYCSOI命令创建约束方程时,一条命令就可以生成多 个约束方程。此时,建议使用拉格朗日乘子法。 Cakey=1:" Quick Solution"是一个快速处理方法,占用CPU时间 接近于直接消去法。但是,提取较高阶频率值一般是实际值的1 2%。当高阶频率比低阶频率高出二次或更高次的数量级时,就会 第4页
ANSY动力学分析指南 模态分析 出现这种误差 Cakey=0:" Accurate Solution"是一个严密精确的方法。但是,占 用CPU的时间大致是" Quick Solution"的两倍 §1.4矩阵缩减技术和主自由度选择准则 下面介绍如何矩阵缩减技术以及选择主自由度(DOF)的基本准则 §1.4.1矩阵缩减技术 矩阵缩减是通过缩减模型矩阵的大小以实现快速、简便的分析过程的方法。它主要用于动力学 分析,如模态分析、谐响应分析和瞬态动力学分析。矩阵缩减也用于子结构分析中以生成超单元。 矩阵缩减允许按照静力学分析那样建立一个详细的模型,而仅将“有动力学特征”部分用于动 力学分析。可以通过辨识定义为主自由度的关键自由度来选择模型的“有动力学特征”部分,但必 须注意,主自由度应足以描述系统的动力学行为。 ANSYS程序根据主自由度(DOF)来计算缩减矩 阵和缩减自由度(DOF)解,然后通过执行扩展处理将解扩展到完整的自由度(DOF)集上。矩阵 缩减的主要优点是,计算缩减解可以大大节省CPU时间,大问题的动力学分析时更是如此。 ANSYS程序釆用的矩阵缩减基础理论是 Guyan缩减法计算缩减矩阵。此法的一个关键假设是 对于较低的频率,从自由度(被缩减掉的自由度(DOF))上的惯性力和从主自由度传递过来的弹性 力相比是可以忽略的。因此,结构的总质量只分配到主自由度(DOF)上。最终结果是缩减的刚度 矩阵是精确的,而缩减的质量和阻尼矩阵是近似的。关于如何计算缩减矩阵的详细内容参见 > §1.4.2人工选择主自由度的准则 选择主自由度是缩减法分析中很重要的一步。缩减质量矩阵的精度(求解精确)将取决于主自 由度的位置和数目。对于给定的问题,可以选择多种不同的主自由度集,在所多种情形下都可以得 到能够接受的结果 用命令M和MGEN来选择主自由度,也可用 TOTAL命令让程序在求解过程中选择主自由度。 建议两种方式兼用:自己选择少量主自由度,同时让 ANSYS程序选择一些自由度。这样,程序将 弥补那些可能被遗漏的模态。 下面是选择主自由度的基本准则 主自由度的总数至少应是感兴趣的模态数的两倍。 把预计结构或部件要振动的方向选为主自由度。 例如对于平板问题,应至少在法向上选择几个主自由度(见图la)。如果在一个方向上的运动 会引起另一个方向上的大运动时,应在两个方向上都选择主自由度(见图b)。 11
ANSYS动力学分析指南 模态分析 图1(a)平板可能有的法向主自由度 (b)Ⅹ方向运动引起Y方向运动 3.在相对较大的质量或较大转动惯量但相对较低刚度的位置选择主自由度(见图2)。凸肩或“松 散”连接的结构是这种位置的实例。相反地,不要选择质量相对较小或有较高刚度(如靠近约 束处的自由度(DOF))的位置作为主自由度。 质量单元 ROTX 图2应选择主自由度的位置:(a)大转动惯量(b)大质量 4.如果最关注的是弯曲模态,则可以忽略转动和“拉伸”自由度 5.如果要选的自由度属于一个耦合约束集,则只须选中耦合集中第一个(首要的)自由度 6.在施加力或非零位移的位置选择主自由度 7.对于轴对称壳模型( SHELL51或 SHELL61),选择模型中的平行于或接近平行于中心线部分的 所有节点的全局U自由度为主自由度,这样就可以避免主自由度间的振荡运动(见图3)。如 果运动基本上是平行于中心线,这条建议可以放宽。对于MODE≥2的轴对称周期单元,应将 其UX、UZ自由度都选择为主自由度 图3在轴对称壳模型中选择主自由度 检查主自由度集的有效性的最好方法是用两倍(或一半)数目的主自由度再次进行分析然后 较结果。另一种方法是观察在模态分析解中输出的缩减质量分布。缩减质量最起码在运动的主要方 向上的分量应该占结构整个质量的10%~15%。 §14.3程序选择主自由度的要点 如果让 ANSYS程序选择主自由度(命令[ TOTAL]),选出的主自由度的分布将取决于求解时单 元被处理的顺序。例如,程序将按单元是从左到右还是从右到左被处理的而选择出不同的主自由度 集。然而,这种差异通常在结果中只会产生无关紧要的差别 对于有统一的大小和特征的网格(如平板),主自由度通常不会是统一的。在这种情况下,应当 用命令M和MGEN人为地指定一些主自由度。在质量分布不规则的结构中也应做同样的处理,因 第6页
ANSYS动力学分析指南 模态分析 为程序选出的主自由度可能集中在高质量区。 §1.5模态分析过程 模态分析过程由四个主要步骤组成 建模; 2.加载及求解 3.扩展模态 4.观察结果。 下面分别展开进行详细讨论: §1.6建模 主要完成下列工作:首先指定工作名和分析标题,然后在前处理器(PREP7)中定义单元类型 单元实常数、材料性质以及几何模型。 ANSYS的《建模和网格指南》中对这些工作有更详细的说明。 注意以下两点: 在模态分析中只有线性行为是有效的。如果指定了非线性单元,它们将被当作是线性的 例如,如果分析中包含了接触单元,则系统取其初始状态的刚度值并且不再改变此刚度值 材料性质可以是线性的,各向同性的或正交各向异性的,恒定的或和温度相关的。在模态 分析中必须指定杨氏模量EX(或某种形式的刚度)和密度DENS(或某种形式的质量) 而非线性特性将被忽略。 §1.7加载及求解 主要完成下列工作:首先定义分析类型、指定分析设置、定义载荷和边界条件和指定加载过程 设置,然后进行固有频率的有限元求解。在得到初始解后,再对模态进行扩展,以供査看。扩展模 态将在下一节“扩展模态”中进行详细说明。 §1.7.1进入 ANSYS求解器 命令:/SOLU GUI: Main Menu>Solution §1.72指定分析类型和分析选项 ANSYS提供的用于模态分析的选项如下表所示,表中的每一个选项都将在随后详细解释 分析类型和分析选项 命令 GUI选择途径 New Analys ANTYPE Main Menu>Solutions-Analysis Type-New Analysis Analysis Type: Modal ANTYPE Main Menu> Solution>-Analysis Type- New Analysis >Modal MODOPT Main Menu>Solution> Analysis Options Number of Modes to Extract MODOPT ain Menu> Solution>Analysis Options
