[0095] 下面结合附图以及具体实施方式进一步说明本发明。
[0096] 如图1所示的系统构架图可知,本发明一种基于不确定优化的飞行器结构功能一体化设计平台系统,由气动外形设计优化模块、热防护系统(前舱)设计优化模块、热防护系统(机翼)设计优化模块、复合材料冷结构设计优化模块、机翼防热承载一体化结构优化模块、机翼气动/结构/热耦合优化分析模块、基础数据库模块及用于集成上述模块的系统框架程序组成。整体实现过程如下:
[0097] (一)气动外形设计优化模块:用于实现气动外形确定性分析及优化、不确定性分析及优化的功能;对于确定性分析及优化部分,根据给定的外形参数进行外形参数化建模,生成参数化模型;将几何文件导入有限元分析软件中进行网格划分,生成网格文件;对网格文件进行气动力/热计算,得到计算结果文件;将结果文件导入结果可视化软件进行数据后处理;最终以全弹道总加热量最小为优化目标,以升阻比为约束条件开展气动外形的确定性优化;对于不确定性分析及优化部分,根据给定的气动外形参数以及不确定参数对不确定参数进行定量化表征与参数化建模,生成参数化模型;将几何文件导入有限元分析软件中进行表面网格划分,生成网格文件;对网格文件进行气动力/热计算,得到全弹道总加热量;在此基础上,利用不确定性传播分析方法得到总加热量的鲁棒性指标,以满足鲁棒性约束条件下弹道总热量最小为目标开展气动外形的不确定性优化;
[0098] (二)热防护系统(前舱)设计优化模块:用于实现热防护系统(前舱)确定性分析及优化、不确定性分析及优化的功能。对于确定性分析及优化部分,根据给定的外形参数对飞行器前舱进行外形参数化建模,生成参数化模型;将几何文件导入有限元分析软件中进行表面网格划分,生成网格文件;利用已建立的气动热数据库插值得到表面热流密度,实现前舱热分析;根据输出的表面温度计算前舱温度场、表面温度以及红外辐射强度;最终以前舱结构重量和红外辐射特性为优化目标,以结构温度为约束条件开展气动外形的确定性优化;对于不确定性分析及优化部分,根据给定的前舱结构设计参数以及不确定参数,在确定性热分析和红外辐射特性分析的基础上,利用不确定性传播分析方法得到温度场分散性,并计算得到其可靠度;在此基础上,以温度可靠度为约束条件,以前舱结构质量和表面红外辐射强度最小为目标开展前舱结构的不确定性优化;
[0099] (三)热防护系统(机翼)设计优化模块:用于实现热防护系统(机翼)确定性分析及优化、不确定性分析及优化的功能。对于确定性分析及优化部分,根据给定的外形参数对飞行器机翼进行外形参数化建模,生成参数化模型;将几何文件导入有限元分析软件中进行表面网格划分,生成网格文件;利用已建立的气动热数据库插值得到表面热流密度,实现前舱热分析;根据输出的表面温度计算机翼温度场、表面温度以及红外辐射强度;最终以机翼结构重量和红外辐射特性为优化目标,以结构温度为约束条件开展气动外形的确定性优化;对于不确定性分析及优化部分,根据给定的机翼结构设计参数以及不确定参数,在确定性热分析和红外辐射特性分析的基础上,利用不确定性传播分析方法得到温度场分散性,并计算得到其可靠度;在此基础上,以温度可靠度为约束条件,以机翼结构质量和表面红外辐射强度最小为目标开展机翼结构的不确定性优化;
[0100] (四)复合材料冷结构设计优化模块:用于实现复合材料冷结构确定性分析及优化、不确定性分析及优化的功能。对于确定性分析及优化部分,利用已有的网格模型,对复合材料结构中具有相同铺层角度的层进行合并,形成超级层模型;利用有限元分析软件进行有限元模型的参数化建模;对冷结构进行强度和刚度计算分析;最终以轻量化为目标,各铺层角度为优化变量,强度和刚度为约束进行优化设计;对于不确定性分析及优化部分,根据已有的参数化模型以及不确定参数,利用有限元程序实现前舱超级层模型的可靠性和鲁棒性分析;在此基础上以轻量化为目标,各区域铺层层数为变量进行优化设计;
[0101] (五)机翼气动/结构/热耦合优化分析模块:用于实现机翼的气动、结构以及热耦合确定性分析及优化、不确定性分析及优化的功能;对于确定性分析及优化部分,根据给定的复合材料机翼结构参数进行参数化建模,生成参数化模型;将几何文件导入有限元分析软件中进行表面网格划分,生成网格文件;对网格文件进行气动力计算,得到计算结果文件;将结果文件导入有限元分析软件进行结构、气动力以及热耦合分析,得到分析结果;最终以机翼结构重量为优化目标,以机翼强度和刚度为约束条件开展气动外形的确定性优化;对于不确定性分析及优化部分,根据给定的机翼结构参数以及不确定参数对不确定参数进行不确定性量化;在确定性热耦合与力耦合分析的基础上,利用不确定性传播分析方法得到分散性结果,获得强度及刚度可靠性指标;在此基础上,以机翼结构重量作为优化目标,以强度和刚度可靠性指标作为约束条件,开展机翼结构的不确定性优化;
[0102] (六)机翼防热承载一体化结构优化模块:用于实现机翼防热承载一体化结构的确定性分析及优化、不确定性分析及优化的功能;对于确定性分析及优化部分,根据给定的机翼外形参数进行外形参数化建模,生成整体机翼参数化模型以及单包结构参数化模型;将整体机翼以及单包结构几何文件分别导入有限元分析软件中,进行热耦合分析;最终以机翼热量及承载为优化目标,以要求的工况条件为约束条件开展机翼防热承载一体化结构的确定性优化;对于不确定性分析及优化部分,根据给定的机翼外形参数以及不确定参数对不确定参数进行不确定性量化;在确定性热耦合分析的基础上,利用不确定性传播分析方法得到分散性结果,获得机翼热量、承载以及不确定性可靠性指标;在此基础上,以工况条件作为约束条件,开展机翼防热承载一体化结构的不确定性优化;
[0103] (七)基础数据库模块:用于处理气动力/热数据库、一体化模型库、算法库等数据资源;包括气动力气动热数据库的维护与应用模块以及设计模型数据库的维护与应用模块;
[0104] (八)系统框架程序:用于实现各子模块的集成和数据的传递与管理,主要功能包括对各模块确定性分析优化设计和不确定性分析优化设计工程的建立、打开操作及对各子模块的调用,以及对基础数据库的调用;
[0105] 上述各模块的实现过程如下:
[0106] 1.气动外形设计优化模块实现方法
[0107] 该模块的实现过程如图2和图3所示:
[0108] (1)从输入参数文件中解析各参数值,其中输出参数文件,最初在单独运行此模块时由用户对输入参数的选择和设置生成,当此模块被不确定性分析模块或确定性优化模块调用时,不确定性分析模块或确定性优化模块将按需要对文件中的数据进行修改;
[0109] (2)打开CATIA外形参数化建模模块,将参数值写入到快速参数化建模窗口的相应位置,生成参数化模型文件;
[0110] (3)打开网格划分模块,调用有限元分析软件,对参数化模型文件进行网格划分,得到网格文件;
[0111] (4)打开弹道气动力/热计算模块,读取弹道数据文件,同时将各个结构的网格文件进行合并,依次计算结构外形的边界层外缘参数、热流参数;
[0112] (5)打开结果后处理模块,将结果文件导入结果可视化软件;
[0113] (6)对于确定性优化,打开确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,以全弹道总加热量最小为优化目标,以升阻比为约束条件开展气动外形的确定性优化循环;
[0114] (7)对于不确定性优化,打开不确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,利用不确定性传播分析方法得到总加热量的鲁棒性指标,以满足鲁棒性约束条件下弹道总热量最小为目标开展气动外形的不确定性优化。
[0115] 2.热防护系统设计优化模块实现方法
[0116] 该模块的实现过程如图4和图5所示:
[0117] (1)从输入参数文件中解析各参数值,其中输出参数文件,最初在单独运行此模块时由用户对输入参数的选择和设置生成,当此模块被不确定性分析模块或确定性优化模块调用时,不确定性分析模块或确定性优化模块将按需要对文件中的数据进行修改;
[0118] (2)打开CATIA外形参数化建模模块,将参数值写入到快速参数化建模窗口的相应位置,生成参数化模型文件;
[0119] (3)打开气动热数据库模块,浏览已建立的气动热数据库,利用库中对应的图表,插值得到表面热流密度;
[0120] (4)打开网格划分模块,调用有限元分析软件,对参数化模型文件进行网格划分,得到网格文件;
[0121] (5)打开温度场计算机后处理模块,导入之前步骤得到的表面热流密度以及表面网格文件,计算得到对应结构的温度场、表面温度以及红外辐射强度;
[0122] (6)对于确定性优化,打开确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,以结构重量和红外辐射特性为优化目标,以结构温度为约束条件开展气动外形的确定性优化;
[0123] (7)对于不确定性优化,打开不确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,以温度可靠度为约束条件,以结构质量和表面红外辐射强度最小为目标开展机翼结构的不确定性优化。
[0124] 3.复合材料冷结构设计优化模块实现方法
[0125] 该模块的实现过程如图6和图7所示:
[0126] (1)从输入参数文件中解析各参数值,其中输出参数文件,最初在单独运行此模块时由用户对输入参数的选择和设置生成,当此模块被不确定性分析模块或确定性优化模块调用时,不确定性分析模块或确定性优化模块将按需要对文件中的数据进行修改;
[0127] (2)打开有限元参数化建模模块,将参数值写入到快速参数化建模窗口的相应位置,生成参数化模型文件;
[0128] (3)打开刚度/强度计算模块,导入参数化模型文件,对冷结构进行强度和刚度计算分析;
[0129] (4)打开刚度/强度分析及后处理模块,调用有限元分析软件,对强度和刚度计算结果进行分析;
[0130] (5)对于确定性优化,调用确定性优化模块,打开集成计算软件,通过预先设置好的计算与处理顺序,以轻量化为目标,各铺层角度为优化变量,强度和刚度为约束先进行层数优化循环,再进行顺序优化;
[0131] (6)对于不确定性优化,调用不确定性优化模块,打开集成计算软件,通过预先设置好的计算与处理顺序,以轻量化为目标,各区域铺层层数为变量先进行层数优化循环,再进行顺序优化。
[0132] 4.机翼气动/结构/热耦合优化分析模块实现方法
[0133] 该模块的实现过程如图8和图9所示:
[0134] (1)从输入参数文件中解析各参数值,其中输出参数文件,最初在单独运行此模块时由用户对输入参数的选择和设置生成,当此模块被不确定性分析模块或确定性优化模块调用时,不确定性分析模块或确定性优化模块将按需要对文件中的数据进行修改;
[0135] (2)打开CATIA外形参数化建模模块,将参数值写入到快速参数化建模窗口的相应位置,生成参数化模型文件;
[0136] (3)打开有限元自动建模模块,调用有限元分析软件,对参数化模型文件进行网格划分,得到表面网格文件;
[0137] (4)打开网格转换模块,调用有限元分析软件,对得到的表面网格文件进行文件转换,得到网格文件;
[0138] (5)打开气动力/热计算模块,根据得到的网格文件进行气动力和热计算,并调用可视化软件浏览计算结果;
[0139] (6)打开气动/结构/热耦合分析模块,调用有限元分析软件,将前面步骤中的计算结果进行耦合分析;
[0140] (7)对于确定性优化,打开确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,以机翼结构重量为优化目标,以机翼强度和刚度为约束条件开展气动外形的确定性优化循环
[0141] (8)对于不确定性优化,打开不确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,以机翼结构重量作为优化目标,以强度和刚度可靠性指标作为约束条件,开展机翼结构的不确定性优化循环。
[0142] 5.机翼防热承载一体化结构优化模块实现方法
[0143] 该模块的实现过程如图10和图11所示:
[0144] (1)从输入参数文件中解析各参数值,其中输出参数文件,最初在单独运行此模块时由用户对输入参数的选择和设置生成,当此模块被不确定性分析模块或确定性优化模块调用时,不确定性分析模块或确定性优化模块将按需要对文件中的数据进行修改。
[0145] (2)打开CATIA外形参数化建模模块,将参数值写入到快速参数化建模窗口的相应位置,生成参数化模型文件。
[0146] (3)打开热耦合分析模块,调用有限元分析软件,对参数化模型文件进行热耦合分析。
[0147] (4)对于确定性优化,打开确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,以机翼热量及承载为优化目标,以要求的工况条件为约束条件开展机翼防热承载一体化结构的确定性优化循环。
[0148] (5)对于不确定性优化,打开不确定性优化模块,调用集成计算软件,通过预先设置好的计算与处理顺序,以工况条件作为约束条件,开展机翼防热承载一体化结构的不确定性优化循环。
[0149] 6.各个模块确定性优化模块与不确定性优化模块的实现
[0150] 各个模块下确定性优化模块与不确定性优化模块执行优化的过程如图12所示:
[0151] (1)判断是否存在配置文件。如果存在,读取优化配置文件内容在界面上显示;如果不存在,读取输入参数文件信息,并在界面显示。
[0152] (2)用户在界面上设置优化模型,包括选择优化变量,设置目标函数和约束条件,选择优化算法并对参数进行设置等。
[0153] (3)将优化模型信息保存到优化配置文件中。
[0154] (4)执行优化。调用优化软件集成程序,将所需的信息包括优化文件夹路径,传递给该程序。
[0155] (5)通过API将从优化配置文件中读取的优化模型信息传入优化软件模板文件中。
[0156] (6)调用优化设计软件,对修改后的模板文件进行优化计算。优化计算运行结束后,将优化结果信息导出保存到工程路径下,并向优化模块程序发送消息。
[0157] (7)程序读取优化结果信息、显示在界面上,并将结果保存到优化配置文件中。
[0158] 确定性优化模块与不确定性优化模块的区别在于执行优化时,确定性优化调用确定性分析模块的执行程序,不确定性优化则调用不确定性分析模块的执行程序。
[0159] 确定性优化模块的执行优化步骤如下:
[0160] (1)将用户设置的输入参数初始值数据写入确定性分析的输入文件。执行确定性分析,得到计算结果即输出文件,由输出文件中解析出相应数据。
[0161] (2)由优化算法计算出下一个的迭代点。
[0162] (3)将下一个迭代点的数据写入确定性分析的输入文件,执行确定性分析,得到输出文件,由输出文件中解析出相应数据。判断数值是否满足约束,目标是否收敛,若满足约束、目标收敛则输出优化结果,否则计算出下一个迭代点,执行此步至优化结束。
[0163] 不确定性优化模块的执行优化步骤如下:
[0164] (1)将用户设置的输入参数初始值数据写入不确定性分析的输入文件。依据在不确定性分析模块中进行的设置,执行不确定性传播流程,得到各响应的不确定性描述参数,进行指标计算,得到指标数值,生成包括响应值和指标数值的输出文件。由输出文件中解析出相应数据。
[0165] (2)由优化算法计算出下一个的迭代点
[0166] (3)将下一个迭代点的数据写入不确定性分析的输入文件,执行不确定性分析,得到输出文件,由输出文件中解析出相应数据。判断数值是否满足约束,目标是否收敛,若满足约束、目标收敛则输出优化结果,否则计算出下一个迭代点,执行此步至优化结束。
[0167] 7.系统框架程序的实现
[0168] 系统框架程序与各模块之间的数据传递均采用读写db文件的方式,实现方法如下:
[0169] 建立基于不确定性优化的结构功能一体化设计工程:
[0170] (1)根据用户输入的路径和工程名称,在路径下建立以工程名称命名的文件夹。
[0171] (2)在文件夹中建立与各个模块对应的文件夹;
[0172] (3)将各模块运行所需文件的模板拷贝到对应文件夹下,并对文件中的参数进行初始化;
[0173] 各子模块的调用:
[0174] (1)获取当前工程所在文件夹路径和优化系统所在路径;
[0175] (2)当前工程所在文件夹路径,加上当前调用模块对应文件夹的名称,就可获得模块调用所需的文件夹路径,如所需的是某个文件的路径,加上该文件名即可;
[0176] (3)优化系统所在路径加上其中调用模块的位置,获得该模块执行程序路径;
[0177] (4)由调用模块执行程序路径、模块调用所需的文件或文件夹路径及命令参数,生成命令行;
[0178] (5)用命令行方式调用模块。
[0179] 本发明未详细阐述部分属于本领域公知技术。
[0180] 以上所述,仅为本发明部分具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域的人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。