自动多重子结构法AMLS振型向量快速计算方案.pdfVIP

北京力学会第18届学术年会论文集：参评青年优秀论文 自动多重子结构法AMLS的振型向量快速计算方案事 Heinrich 尹家聪1’2 Voss2陈璞1 (1北京大学工学院力学与空天技术系，100871； 2Instituteof of 1071) Numerical Simulation，HamburgUniversityTechnology，Germany，2 摘要：针对超大规模有限元快速模态分析的自动多重子结构法(AMLS)，提出并实现了一种 基于AMLS剖分树先根遍历和块高斯消去逆变换的快速特征向量计算方案。与已知的计算 方案相比，此方案无需显式组装所谓的变换映射矩阵，可将振型向量及其他动力响应的计 算效率提高至少一到两个数量级． 关键词：自动多重子结构法，模态分析，特征值，振型向量 一、引言 Multi．1evel 自动多重子结构法(Automated 种新型模态分析算法【l，2]。AMLS通过图剖分将结构的刚度、质量矩阵按其稀疏特性， 代数地划分为多个层次和成千上万个子结构；并由子结构的模态综合出整体结构的模 态，显著地缩减了计算量并保持了相当的精度。AMLS具备一次求解几千阶模态的能力： 有报告称AMLS可在工作站上求得290万阶问题的四万多子结构模态，最终综合出824 个整体结构的低阶模态；与在超级计算机上运行的MSC．Nastran相比，AMLS的特征值 误差在3％以下，求解时间仅为后者的一半[2】。目前，AMLS已被集成于MSC．Nastran、 Abaqus等多个商用有限元软件中，并广泛应用于汽车工业的频响和振动．噪声耦合分析。 然而，作为一种模态分析算法，AMLS的振型向量计算尚未完善。为此，本文首先 完成了文献f11中建议的一种基于显式组装变换映射矩阵的振型向量计算方案：然后，提 出并实现了一种基于块高斯消去逆变换的振型向量快速计算方案，并对比了计算效率。 二、自动多重子结构法 设初始结构共被划分为Ⅳ个子结构，并按AMLS剖分树的“后根遍历”顺序排列 为{l，2…．，册。那么在数学意义上，AMLS的计算流程为： 步骤一：利用图剖分工具将待求解的振动方程 Kx=2Mx (1) 的刚度矩阵K、质量矩阵M进行自由度重排，以减少块高斯消去过程中产生的填充元。 步骤二：对重排后的刚度、质量矩阵按剖分树的“后根遍历”顺序执行块高斯消去 被化为块对角阵；而质量矩阵丽=U1MU的分块特性不变，但大量子块变稠密。 步骤三：忽略块高斯消去后质量矩阵丽中的非对角耦合项，求解子结构模态： K，，ZI=MffzfAf，(f-1，．．．，Ⅳ) (2) 其中A。=diag[2},彳熙…，彳：哪)】和z。=[z％z(2’，．．．，z：州’】分别为特征值和振型向量矩阵。 ’国家自然科学基金项目；留学基金委公派联合培养博士生资助项目(2010601199) I．53 北京力学会第18届学术年会论文集：参评青年优秀论文 步骤四：利用N个子结构组成的振型矩阵Z=[z，，Z：一．，ZⅣ】，将(文，fit)降阶为两个 步骤五：求解降阶特征值问题K。Y=gM。Y。解出的q个∥即为方程(1)式的“AMLS 近似特征值”；而V=uZY为“AMLS近似振型向量”，Y=【J，I，Y2，．．．，％】。 三、两种振型向量计算方案 实际计算程序中，步骤一至步骤四均是按照AMLS剖分树后根遍历的顺序“交错执 行”的。记之，=diag[1l．．．，Is-I,Z。，I州，．．．IⅣ】，那么V=UZY的实际数学表达式应为： V；U1幺。U：之：…UN