第4章 动力分析降阶.doc

第四章 动力分析降阶 3.1 概述 1 定义 动力分析降阶：指将给定的动力学数学模型用具有较自由度动力学模型替代 2 动力学降阶的原因 数学模型太大，不降阶无法求解 数学模型比实际要求的太详细，完全求解代价太大 动力学降阶比建立一个单独的、较小动力学模型分析更精确，并且更加经济 3.2 MSC/NASTRAN中使用的降阶方法 Guyan降阶法（静态缩聚） 广义动力降阶法(GDR) 模态降阶法 分量模态综合法（(superelement) 3.2.1 静态缩聚 （1）理论 若为结构的未约束（自由）坐标集合，分为 其中， 将刚度矩阵分为A-Set 和O-Set, 得到 如果P0 = 0, 求解u0得到 其中， 因此得到 动力学方程为 矩阵Maa, Baa, Kaa和Pa的维数比原方程相应维数小 （2）Nastran卡片 或者 或者 OMIT, OMIT1 注：（1） 指定A集合用ASET卡片, 指定O集合用OMIT卡片； （2）同一节点的同一自由度不能同时指定为A集合与O集合 （3）静态缩聚求解控制 （4）特点 A集合选择与经验有关 分析结果精度与A集合选择有关，因此与用户经验有关 局部模态可能丢失 高阶模态出现较大误差 A集合选择有时困难 3.2.2 模态降阶 Nastran求解线性动力问题的两种方法 直接法：根据A-SET中的坐标直接求解 模态法：在模态坐标下进行求解（H-SET）（A集合的运动方程以模态坐标表示，成为 H集合） (1)模态降阶理论 将质量矩阵和刚度矩阵分为A-SET和O-SET, 并进行模态分析得到模态矩阵 并进行模态变换 如果为质量标准化模态矩阵，则得到 （2）模态降阶求解控制 3.2.3 例子 问题：使用静态缩聚方法对平板进行模态分析 （1）有限元网格 （2）载荷边界条件 注：O为要保留的结点 (3) 输入文件 ID SEMINAR, PROB2 SOL 103 TIME 10 CEND TITLE = REDUCTION PROCEDURES, NORMAL MODES EXAMPLE SUBTITLE = USING STATIC REDUCTION ECHO = UNSORTED SUBCASE 1 SUBTITLE=USING GIVENS METHOD = 1 SPC = 1 VECTOR=ALL BEGIN BULK EIGR, 1, AGIV, , , , 5 PARAM, COUPMASS, 1 PARAM, WTMASS, 0.00259 INCLUDE ’plate.bdf’ $ $ SELECT A-SET, STATIC REDUCTION IS DONE AUTOMATICALLY $ ASET1, 345, 3, 5, 7, 9, 11 ASET1, 345, 25, 27, 29, 31, 33 ASET1, 345, 47, 49, 51, 53, 55 ENDDATA （4）部分结果文件 CQUAD4 12 1 13 14 25 24 CQUAD4 13 1 14 15 26 25 CQUAD4 14 1 15 16 27 26 CQUAD4 15 1 16 17 28 27 CQUAD4 16 1 17 18 29 28 CQUAD4 17 1 18 19 30 29 CQUAD4 18 1 19 20 31 30 CQUAD4 19 1 20 21 32 31 CQUAD4 20 1 21 22 33 32 CQUAD4 21 1 23 24 35 34 CQUAD4 22 1 24 25 36 35 CQUAD4 23 1 25 26 37 36 CQUAD4 24 1 26 27 38 37 CQUAD4 25 1 27 28 39 38 CQUAD4 26 1 28 29 40 39 CQUAD4 27 1 29 30 41 40 CQUAD4 28 1 30 31 42 41 CQUAD4 29 1 31 32 43 42 $ $ plate.bdf $ $ geometric input file for plate model $ PSHELL 1 1 .1 1 1 CQUAD4 1 1 1 2 1312 CQUAD4 2 1 2 3 1413 CQUAD4 3 1 3 4 1514 CQUAD4 4 1 4 5 1615 CQUAD4 5 1 5 6 1716 CQUAD4 6 1 6 7 1817 CQUAD4 7 1 7 8 1918 CQUAD4 8 1 8 9 2019 CQUAD4 9 1 9 102120 CQUAD4 10 1