航空结构大规模并行分析与优化应用

开篇：润墨网以专业的文秘视角，为您筛选了一篇航空结构大规模并行分析与优化应用范文，如需获取更多写作素材，在线客服老师一对一协助。欢迎您的阅读与分享！

摘要：

针对飞机设计精细化数值分析模型自由度已经达到亿级，对高性能计算的要求也越来越高的问题，围绕大规模并行计算环境下结构分析和优化的若干关键问题，研究满足高性能计算体系特点的区域分解并行算法、超大规模结构变量敏度高效求解和结构非线性振动特性求解等关键技术.对国产CAE软件HAJIF进行并行化改造，初步实现基于最大航程的气动结构综合优化设计和基于精细化模型的复合材料机翼综合优化设计.HAJIF的计算效率和精度得到明显提高.

关键词：

区域分解算法； 超大规模结构变量； 解析敏度； HAJIF

中图分类号： V221

文献标志码： B

0 引 言

进入21世纪以来，以高性能计算机为基础的计算科学得到长足的发展，在航空航天、地震预测、核爆模拟等许多工业领域得到应用，高性能计算水平的高低已经成为衡量企业、研究机构科技水平的重要标志，更是检验一个国家高科技能力的试金石.高性能计算是未来高精度数值仿真和数值试验的基本要求，其基础是并行计算机硬件与其上运行的高性能计算软件，两者缺一不可.随着天河系列、曙光系列等高性能计算机的发展，我国在高性能计算的硬件研究和制造方面已经达到国际先进水平，但是在软件方面，仍缺乏自主研发的高性能计算系统.从飞行器结构强度分析领域来看，目前国内可解决飞行器综合优化设计问题的系统性软件平台还处于空白.另一方面，由于CAE软件市场长期由国外垄断，授权费用昂贵，战略性、高水准的工业工程数值模拟软件发展受到限制.

1 高性能并行结构分析与优化的若干关键技术

现代工业设计的产品复杂性越来越高，对高性能计算的要求也越来越高，例如我国自主研发的某型飞机，其精细化模型自由度已经达到数亿量级.这对高性能计算能力，包括并行计算、超大规模高效数值算法等提出更高的要求，涉及到以下几类关键技术.

1.1 满足高性能计算体系特点的区域分解并行算法

并行计算方法的发展一般遵循2个途径：重构和分而治之.按照算法设计特点的不同，现有的结构并行分析技术也归为2类：一类称为算法级的并行有限分析技术，根据有限元分析过程中不同环节的计算特点，通过重构计算序列来实施并行求解；另一类称为子结构级的并行有限元分析技术，是指将一个大结构分解成多个子结构，然后根据分而治之的策略进行并行求解.区域分解法是后者的典型代表，其算法是高度并行的，即计算的主要步骤在各子域内独立进行，同时允许在不同子域中选用不同的数学模型，以便整体模型更适合于工程物理实际，具有很好的工程适用性.[13]

采用区域分解算法进行结构计算时，依照大规模并行计算特点，先按全机大部件（机翼、机身、中央翼、垂尾、平尾等）分别建立结构的大部件分析模型.每个大部件分别进行自动区域划分，按单元的不同计算量分别考虑不同的权数，形成静态负载平衡的若干子结构（与高性能计算机的计算节点数匹配）.在每个子结构内部将子结构全部自由度缩聚到子结构的边界点上，形成以子结构边界点为未知量的子结构刚度矩阵和载荷向量矩阵，按自由度形成全局界面方程.用分布式预处理共轭梯度（Distributed Preconditioned Conjugate Gradient， DPCG）法求解界面方程，为提高并行效率，按结构的受载工况数量将界面方程划分为若干个单工况的界面方程进行并行求解.最后，将边界点的位移回代到每个子结构计算子结构内部自由度.[46]求解流程见图1.

1.2 超大规模结构变量敏度高效求解技术

经过几十年的发展，参数优化依然是航空结构优化设计的最重要手段.随着飞机综合性能和设计生产自动化水平越来越高，超大规模变量优化成为飞机设计的关键技术.[3]采用自然网格进行全机有限元模拟时，模型节点数一般为数千至上万，而采用精细网格进行全机有限元模拟时，节点数可达数十万至上千万，精细化模型自然带来结构设计变量的大幅增长，当前国际先进优化设计水平的设计变量为10万个.[78]

敏度求解的常用方法有3种：有限差分法、半解析法和解析法.按照约束性质的不同，超大规模（10万以上）结构变量的敏度求解方法分为2类：对数量较少的全局约束（位移、模态）采用直接解析法进行求解，对数量众多的局部约束（应力、应变）采用伴随矩阵法等半解析法进行求解.敏度求解过程中需要对刚度矩阵求逆，流程见图2.

但是，实际工程结构刚度矩阵往往规模太大或奇异性较强，不能直接求逆，虚载荷法可有效解决该问题.在虚载荷计算中提前建立约束及实位移、虚位移关系表，将实、虚位移挑选过程从敏度计算循环过程中分离，按边界条件计算虚载荷，可有效提高敏度计算效率.同时，在结构微分刚度矩阵组装中建立设计变量与关联的属性卡及单元编号关系表，避免微分刚度计算时反复查找单元和属性.

1.3 结构非线性振动特性求解技术

在变形较大的情况下，大展弦比机翼位移与应变关系的线性假设不再适用，此时小应变假设不能正确反映其物理本质，必须采用大应变理论.在求解结构发生大变形时的切线刚度矩阵时，基于更新的拉格朗日列式的非线性方程的求解方法NewtonRaphson法，由Gaussia消去法求解增量位移的线性方程组，最K误差向量由Jacobian法形成结构的切线刚度矩阵（几何大变形的贡献）.在得到结构大变形下的切线刚度矩阵后，根据求解结构的不同特点，分别采用逆幂法、子空间迭代和FEAST法进行非线性振动特性的计算，求解流程见图3.分别计算载荷包线内所有工况下的变形及其非线性振动特性，形成非线性振动特性库，供耦合分析调用.[910]

2 高性能并行结构分析与优化应用

2.1 结构分析优化平台简介

航空结构强度分析与优化系统HAJIF是中航工业强度所研制推出的目前国内航空领域功能最为全面的大型CAE软件系统.该系统以强度试验数据库为支撑，提供飞行器结构基础分析、优化设计、气动弹性分析、热分析等功能[11]，见图4.系统提供图形前后置功能，集成自动快速建模、切面刚度计算、细节强度分析等航空特色模块.系统采用先进的开放式、可扩充的软件架构，可方便地为用户提供个性化定制开发和服务.