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)上。最终结果是缩减的刚度 矩阵是精确的,而缩减的质量和阻尼矩阵是近似的。关于如何计算缩减矩阵的详细内容参见 << ANSYS理论参考手册>> §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