一种自锁定的复合材料预变形舱门结构优化设计方法

一种自锁定的复合材料预变形舱门结构优化设计方法
【专利摘要】本发明涉及一种自锁定的复合材料预变形舱门结构优化设计方法，在飞行器预变形舱门设计中充分考虑预变形曲线形状和铺层形式对门结构的影响，将预变形设计引入到传统的飞行器门结构设计中，发挥门结构预加载后预紧力的锁紧作用，在减轻结构重量的同时保证了可靠性与使用性，同时，采用复合材料设计思想，结合预变形设计与复合材料铺层设计，利用全局优化算法得到全局最优的复合材料预变形门结构设计方案。
【专利说明】
-种自锁定的复合材料预变形舱口结构优化设计方法
技术领域
[0001] 本发明设及复合材料结构优化设计方法领域，特别设及一种飞行器预变形设备舱 口结构优化设计与复合材料设计方法；
【背景技术】
[0002] 在航空航天设计中，保证设计要求的前提下使结构重量最轻一直是设计人员不懈 的追求。随着复合材料的兴起，越来越多的结构采用复合材料取代W往传统的金属材料。与 传统金属材料的均质性不同，复合材料的铺层设计、各向异性等独特的设计方法对设计人 员提出了新的挑战。
[0003] 飞行器设备舱口是航空结构中典型的结构之一，也是保证飞行器使用安全和机内 工作人员安全的重要部件。在使用过程中，飞行器设备舱口主要受到内外压差及气动载荷 的作用，需要具由足够的强度、刚度与可靠性W保持结构的使用性能与使用安全。如果设备 舱不能保证强度、刚度与可靠性，将对飞行器的气动特性与气密性造成严重影响，甚至直接 危及机内人员安全。
[0004] 一方面，是航空航天领域迫切的减重需求不断敦促飞行器结构设计抛弃不必要的 重量，另一方面，是为避免复杂的锁定机构带来的可靠性下降问题。故摸弃锁机构，采用预 变形口取代，舱口关闭后，预变形结构受挤压产生预应力能够有效地减小舱口在受到外部 载荷与惯性载荷作用下的弹性变形，起到锁紧舱口的作用。

【发明内容】

[0005] 本发明的技术解决问题:克服现有技术的不足，提供一种自锁定的复合材料预变 形舱口结构优化设计方法，该方法在飞行器预变形舱口设计中充分考虑预变形曲线形状和 铺层形式对口结构的影响，将预变形设计引入到传统的飞行器口结构设计中，发挥口结构 预加载后预紧力的锁紧作用，在减轻结构重量的同时保证了可靠性与使用性，同时，采用复 合材料设计思想，结合预变形设计与复合材料铺层设计，利用全局优化算法得到全局最优 的复合材料预变形口结构设计方案。
[0006] 本发明技术解决方法:一种自锁定的复合材料预变形舱口结构优化设计方法，包 括W下步骤：
[0007] (1) W 口结构边框顶点位置为定位点，在口的预变形一侧设置特征点，用于拟合出 符合基本设计需求的拟合曲线，选取特征点坐标为优化设计变量，记为X=(X1，X2,…，Xk); 由于工程需要，作为设计变量的特征点位置有相对应的变化范围，即XG [x^xU]，可W依据 经验或要求给出初始设计变量及各个设计变量的约束范围；
[000引（2)基于上述特征点，带入所设计拟合曲线的方程组中，解出若干曲线上点的坐 标，利用样条曲线连接各定位点建立近似拟合曲线，与其他预变形口形状的控制点一并构 建起口结构几何模型，即读入包含特征点坐标的输入文件到软件中并生成几何模型的过 程，上述过程即为读入过程；
[0009] (3)为将口结构几何模型进一步转化为可用于计算的有限元模型，基于CAE工具软 件平台，对所述口结构几何模型进行材料设置、网格划分、建立接触、复合材料铺层的有限 元建模过程，最终建立的舱口结构模型包括壳单元构成的口结构、梁单元构成的筋结构W 及实体单元构成的纵向止挡=部分，上述过程即为建模过程。将由特征点坐标到建立有限 元模型的过程进行参数化，参数化中所设及的参数包括特征点坐标与复合材料铺层厚度；
[0010] (4)得到有限元模型后，分两步进行仿真计算W模拟结构实际使用中的变形，第一 步为预加载过程，模拟预变形口合上时产生预应力的过程，第二步为实际工况加载过程，根 据读入的载荷数据，模拟预变形口使用中受外载荷变形的过程，两步计算后可得到结果数 据，对结果数据进行下一步的后处理，上述过程即为计算过程；
[0011] (5)利用CAE工具软件的后处理功能，提取上述加载计算的结果数据，对结果数据 进行处理，得到与设计要求直接相关的指标数据，W结果文件的形式输出，上述过程即为输 出过程；
[0012] (6)将整个读入、建模、加载、计算、输出过程编写并整合为AP化自动处理命令流文 件，当各设计变量在给定约束范围内取值时，自动根据特征点位置和复合材料铺层厚度生 成参数化的几何模型与有限元模型，并根据读入的载荷数据进行计算，输出相应的结果文 件，建立AP化自动处理命令流文件是使读入、建模、加载、计算、输出过程可W被CAO工具软 件调用，实现自动运行与处理，为下一步进行优化做准备；
[0013] (7)利用CAO工具软件调用CAE工具软件运行步骤(6)得到的AP化自动处理命令流 文件，并控制（2)中读入的输入文件、（5)输出的结果文件的修改与读入，实现集成优化，建 立相应的设计变量与约束变量，构建CAO工具软件中的优化模型；
[0014] (8)在优化模型中，W特征点坐标、复合材料铺层厚度为优化变量，W结构的强度 即应变、刚度即工作载荷下变形后最大间隙作为约束，W结构重量最轻为优化目标，利用全 局寻优算法即化inter算法进行减重优化，可W实现在保证结构满足设计需要的前提下，对 自锁定的复合材料预变形舱口结构进行减重。
[0015] 本发明与现有技术相比的优点在于:本发明W复合材料预变形口结构形状和铺层 厚度由优化设计变量，W结构质量最小为设计目标，W结构强度、刚度为约束条件，在保证 结构使用要求的基础上实现了结构质量最小化。本发明中的读入、建模、计算、输出过程 AP化命令流都是通过自编语句实现。相对于常规的飞行器舱口结构，本方法所采用的预变 形口结构既能减轻结构重量、省去繁琐而降低可靠性的锁定机构，还可W满足飞行器结构 实际的使用要求。考虑到工程实际中针对不同精度、不同尺寸的口结构对预变形的拟合精 度与拟合要求有所差异，本发明采用的建立预变形口方法可W使用任意高阶的幕函数，从 而拟合任意精度与任意形状的预变形曲线，拓展了本方法的适用范围。在设计中，针对复合 材料运一新兴材料进行设计，使得比起通常的航空侣合金结构，本发明所设及结构具有更 大的轻量化优势。
【附图说明】
[0016] 图1为本发明中复合材料预变形舱口结构参数化有限元模型示意图；
[0017] 图2为本发明中CAO工具软件所建立优化模型示意图；
[0018] 图3为本发明实施例优化过程中质量下降过程示意图；
[0019] 图4为本发明的方法实现流程图。
【具体实施方式】
[0020] 如图所示，本发明自锁定的复合材料预变形舱口结构优化设计方法，包括W下步 骤：
[0021 ] (1) W 口结构边框顶点位置为定位点，在口的预变形一侧设置特征点，用于拟合出 符合基本设计需求的拟合曲线。选取特征点坐标为优化设计变量，记为X=(X1，X2,…，Xk)。 由于工程需要，一般作为设计变量的特征点位置有相对应的变化范围，将设计变量的变化 范围作为约束，即XG [/，/]，通常可W依据经验或要求等条件给出初始设计变量及其约束 范围。
[0022] (2)基于上述特征点，带入所设计拟合曲线方程组中，解出若干曲线上点的坐标， WS次曲线为例，经过定位点（0，0)与定位点（1120，0)，特征点（Xi，y 1)，带入方程y = Cix+ C2X2+C3X3得到如下方程组：
[0023] (I)
[0024]
[0025] (2)
[0026] 解出；
[0027] 饼
[0028] 从而得到所设计拟合曲线的代数表达式，按一定间隔取其上若干点求出坐标作为 定位点，W样条曲线依次连接各个定位点的方式建立近似拟合曲线，用预变形侧的边缘与 其他=个直线边缘，得到口结构几何模型。
[0029] (3)基于上述口结构几何模型，在ANSYS软件中继续进行有限元模型建模，壁板的 材料选用各向异性的T700复合材料，加强筋的材料选用各向同性的7050侣合金。壁板采用 壳单元甜化L181，各加强筋采用梁单元肥AM188,梁单元在两个方向上分别采用工字梁和口 字梁，纵向止挡采用实体单元S0LID185,进而建立有限元模型。
[0030] 在壁板上进行复合材料铺层，铺层采用（45°/-45° /0° /90°)S的对称层合板铺层， 厚度按总厚度均匀分配。对特征点坐标和复合材料铺层厚度分别设置参数，建立参数化有 限元模型，最终建立的模型如图1所示，包括复合材料壁板、=个纵向工字加强筋和五个横 向工资加强筋。
[0031] (4)进行加载计算，共分为两步。第一步利用纵向止挡将预变形口压平，模拟预变 形口合上时产生预应力的过程，第二步读入气动载荷，模拟预变形口实际使用中受到外载 荷变形的过程。
[0032] (5)使用CAE工软件对计算结果进行后处理，利用ANSYS的AP化命令流，提取加载计 算过程的结果数据并进行处理，得到设计指定的约束参数与设计优化的目标参数，即预变 形边缘最大间隙、壁板上复合材料最大第一主应变与最小第=主应变W及预变形口结构总 质量，W结果文件的形式输出；
[0033] (6)将上述的整个读入、建模、加载、计算、输出过程编写并整合为AP化自动处理命 令流文件，通过读入一组符合给定约束范围的设计方案的设计参数，可W自动根据特征点 位置参数和复合材料铺层厚度完成参数化有限元建模过程，并根据读入的载荷文件进行两 步计算，输出相应的结果文件；
[0034] (7)利用Isight软件，调用ANSre程序运行上一步得到的AP化整合命令流文件，并 控制读入文件、输出文件的修改与读入，实现Isight对整个建模计算输出过程的集成，在 Isight中建立相应的优化变量、约束变量与目标函数，构建优化模型，优化模型如图所示， 图中自上而下、自左而右依次为优化控制算法(Optimizationl)、输入文件的修改化dit_ Inputfi 1 e)、ANSYS读入输入文件并运行AP化整合命令流(Go ANSYS)、读入输出文件(ReacL Outputf i 1 e)的优化流程。
[0035] (8)在优化模型中，W特征点坐标参数、复合材料铺层厚度参数为优化变量，W结 果文件中的结构强度(应变）、刚度(工作载荷下口结构变形后最大间隙）为约束条件，W预 变形口结构重量最轻为优化目标，采用IS i ght中的全局最优算法(Po inter算法)进行自锁 定预变形复合材料口结构的减重优化，可W实现在保证结构满足设计需要的前提下，对自 锁定的复合材料预变形舱口结构进行减重，优化过程中质量下降如图所示，横坐标为优化 过程中迭代计算次数、纵坐标为每次迭代计算得到的质量，可W看到随着迭代次数增多、结 构质量不断减轻直至收敛。
[0036] W上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制；其可扩展应 用于飞行器舱口结构优化设计领域，凡采用等同变换或者等效替换而形成的技术方案，均 落在本发明权利保护范围之内。
【主权项】
1. 一种自锁定的复合材料预变形舱门结构优化设计方法，其特征在于实现步骤如下： (1) 以门结构边框顶点位置为定位点，在门的预变形一侧设置特征点，用于拟合出符合 基本设计需求的拟合曲线，选取特征点坐标为优化设计变量，记为X =( Xl，X2，…，Xk);根据 工程需要，作为设计变量的特征点位置有相对应的变化范围，即xe [/，/]，依据经验或要 求给出初始设计变量及各个设计变量的约束范围； (2) 基于上述特征点，带入所设计拟合曲线的方程组中，解出若干曲线上点的坐标，利 用样条曲线连接各定位点建立近似拟合曲线，与其他预变形门形状的控制点一并构建门结 构几何模型，即读入包含特征点坐标的输入文件并生成几何模型的过程，上述操作即为读 入过程； (3) 为将门结构几何模型进一步转化为可用于计算的有限元模型，基于CAE工具软件平 台，对所述门结构几何模型进行材料设置、网格划分、建立接触、复合材料铺层的有限元建 模过程，最终建立的舱门结构模型包括壳单元构成的门结构、梁单元构成的筋结构以及实 体单元构成的纵向止挡三部分，上述操作即为建模过程;将由特征点坐标到建立有限元模 型的过程进行参数化，参数化中所涉及的参数包括特征点坐标与复合材料铺层厚度； (4) 得到有限元模型后，分两步进行仿真计算以模拟结构实际使用中的变形，第一步为 预加载过程，模拟预变形门合上时产生预应力的过程，第二步为实际工况加载过程，根据读 入的载荷数据，模拟预变形门使用中受外载荷变形的过程，两步计算后可得到加载计算的 结果数据，对结果数据进行下一步的后处理，上述操作即为计算过程； (5) 利用CAE工具软件的后处理功能，提取上述加载计算的结果数据，对结果数据进行 处理，得到与设计要求直接相关的指标数据，以结果文件的形式输出，上述操作即为输出过 程； (6) 将整个读入、建模、加载、计算、输出过程编写并整合为APDL自动处理命令流文件， 当各设计变量在给定约束范围内取值时，自动根据特征点位置和复合材料铺层厚度生成参 数化的几何模型与有限元模型，并根据读入的载荷数据进行计算，输出相应的结果文件，建 立APDL自动处理命令流文件是使读入、建模、加载、计算、输出过程通过CAO工具软件调用 CAE软件实现自动运行与处理，为下一步进行优化做准备； (7) 利用CAO工具软件调用CAE工具软件运行步骤(6)得到的APDL自动处理命令流文件， 并控制(2)中读入的输入文件、（5)输出的结果文件的修改与读入，实现集成优化，建立相应 的设计变量与约束变量，构建CAO工具软件中的优化模型； (8) 在优化模型中，以特征点坐标、复合材料铺层厚度为优化变量，以结构的强度即应 变、刚度即工作载荷下变形后最大间隙作为约束，以结构重量最轻为优化目标，利用全局寻 优算法即Pointer算法进行减重优化，实现在保证结构满足设计需要的前提下，对自锁定的 复合材料预变形舱门结构进行减重。2. 根据权利要求1所述的一种自锁定的复合材料预变形舱门结构优化设计方法，其特 征在于:所述步骤(2)中，使用幂函数作为拟合曲线，根据特征点的个数确定幂函数的阶次， 将特征点坐标带入幂函数方程中构成方程组，反解出幂函数各阶系数的唯一解，利用已解 出的幂函数上以一定间隔取点得到定位点，利用样条曲线连接各定位点建立近似拟合曲 线，来近似替代拟合曲线建立模型。3. 根据权利要求1所述的一种自锁定的复合材料预变形舱门结构优化设计方法，其特 征在于:所述步骤(3)中，将门所在面以壳单元划分网格并设定复合材料铺层以建立壳单元 构成的门结构，将面上的线以梁单元划分网格并设置加强筋的几何参数以建立梁单元构成 的筋结构，实体单元构成的纵向止挡部分依据预变形门边缘最高点位置与与门结构自动保 持相对距离关系。4.根据权利要求1所述的一种自锁定的复合材料预变形舱门结构优化设计方法，其特 征在于:所述步骤(4)中，利用纵向止挡压平预变形门模拟预变形门合上时产生预应力的力 学过程，其中门上远离预变形一侧边缘模拟可转动不可移动的轴位置，门上预报形一侧边 缘角点施加单向位移约束，与下压的纵向止挡共同模拟门关闭时合上的相对变形情况。
【文档编号】G06F17/50GK105956342SQ201610464009
【公开日】2016年9月21日
【申请日】2016年6月23日
【发明人】王晓军, 蔡逸如, 王磊, 耿新宇, 张泽晟, 樊维超, 刘鑫
【申请人】北京航空航天大学