Sec15_复特征值分析.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * S15-* NAS122, Section 15, August 2010 Copyright? 2010 MSC.Software Corporation S15-* NAS122, Section 15, August 2010 Copyright? 2010 MSC.Software Corporation 第15章 复特征值分析 复特征值分析 复特征值分析用于 颤振分析 计算结构阻尼模态 评估系统的稳定性 转子系统(如转子临界转速) 刹车稳定性(如刹车啸叫分析) …… 质量和刚度矩阵可以是非对称的，也可能包含复数 复特征值分析(续) 具有粘滞阻尼系数的自由运动方程如下 假设位移向量{u(t)}和模态向量{y}之间有如下关系 将方程(2)代入(1)得到 (3) (1) (2) 复特征值分析(续) 方程(3) 描述了复特征值问题 相应的特征方程的根 li 可能是实数，可能是虚数，或者是复数 相应的特征向量可能是复数形式的 复特征值分析(续) 粘滞阻尼矩阵[B]和在频率响应分析中使用的相同 直接法(SOL 107)用d集(即d-set，物理自由度和附加自由度)的质量阵、阻尼阵和刚度阵来求解特征方程 模态法(SOL 110)用h集(即h-set，模态自由度和附加自由度)的质量阵、阻尼阵和刚度阵来求解特征方程 复特征值分析(续) MD Nastran有3种复特征值求解方法： INV – Inverse power逆幂法 HESS – Upper Hessenberg上海森伯格法 CLAN – Complex Lanczos复兰索斯法 复特征值分析(续) INV逆幂法 与实特征值分析中的实INV方法相似 [M]可以是奇异阵 用户必须在复平面中指定搜索区域 对于相当规模的模型，计算花销比较大，不容易得到结果 没有HESS法可靠 不推荐使用该方法 复特征值分析(续) HESS上海森伯格法 与实特征值分析中的GIV方法类似 要求[M]非奇异 对于大自由度问题计算时间耗费很大。因此，除了小自由度问题，通常推荐使用模态法。 求解两类标准方程形式： For [B] = [0], For [B] = [0], 复特征值分析(续) HESS上海森伯格法(续) 用标准方程形式 For [B] = [0] (续) 复特征值分析(续) CLAN复兰索斯法 与实特征值分析中的实兰索斯法类似 是变换法和跟踪法的混合方法 没有类似实特征值计算中的斯图姆序列检查(Sturm Sequence Check) 移频点可能需要指定 用实特征值分析作为参考 复特征值分析(续) 在Patran设置复特征值分析 复特征值分析(续) 执行控制段 直接复特征值分析为SOL 107 模态复特征值分析为SOL 110 复特征值分析(续) 工况控制段 复特征值分析(续) 模型数据段 复特征值分析(续) 模型数据段(续) 复特征值分析(续) 模型数据段(续) 复特征值分析(续) 模型数据段(续) 复特征值分析(续) 模型数据段(续) 复特征值分析(续) 模型数据段(续) 案例分析：复特征值分析 计算振动打桩机复频率和振型 M2 1.5 lbf*sec2/inch M3 3.0 lbf*sec2/inch K1 12500.0 lbf/inch K2 50000.0 lbf/inch C1 30.0 lbf*sec/inch Exciter Pile Ground K1 K2 C1 M2 M3 3 2 1 1.0 1.0 案例分析：复特征值分析(续) 振动打桩机模型 模型文件pile_driver.bdf 案例分析：复特征值分析(续) 振动打桩机(续) 案例分析：复特征值分析(续) 复特征值 案例分析：复特征值分析(续) 复特征向量 案例分析：复特征值分析(续) 复特征向量(续) 案例分析：复特征值分析(续) 注意粘滞阻尼系统的根是复共轭对( a + b ) 稳定系统要求a 0 通用建模建议 复特征值计算时间可能会大大多于实特征值计算 不同于实特征值计算，复特征值计算没有类似斯图姆序列检查的机制 移频之后可能会丢根 通用建模建议续) 求解方法选择 小问