内容简介
《结构序列近似建模优化方法与应用》面向大型运载火箭关键承载舱段轻量化设计需求，聚焦大型复杂结构高效优化设计理论方法，系统研究结构序列近似建模优化方法及其在大型运载火箭加筋圆柱壳舱段设计中的应用。在理论研究方面，提出面向高维非规则设计域的混合整数填充采样算法，实现了设计样本点对非规则设计空间的均匀覆盖；发展了面向设计域全局近似的增广径向基函数近似建模方法、序列近似建模方法和多保真度近似建模方法，提升了大型复杂结构非线性力学问题分析效率。在应用研究方面，重点针对大型运载火箭两类典型加筋圆柱壳舱段，研究考虑整体稳定性和桁条局部稳定性的大直径加筋圆柱壳舱段的优化设计方法，并协同考虑结构承载能力和载荷扩散性能的集中力扩散舱段的优化设计方法，获得满足工程应用需求的结构轻量化设计方案。
目录
目录
“空天前沿技术丛书”序一
“空天前沿技术丛书”序二
前言
第1章 绪论 1
1.1 结构优化概述 1
1.1.1 结构优化设计 1
1.1.2 结构优化算法 4
1.2 结构序列近似优化方法研究进展 6
1.2.1 初始采样方法 6
1.2.2 近似建模方法 9
1.2.3 序列加点采样方法 14
1.3 运载火箭典型薄壁加筋舱段优化设计研究进展 16
1.3.1 薄壁均匀加筋圆柱壳优化设计 16
1.3.2 集中力扩散舱段优化设计 18
第2章 高维非规则设计域的混合整数序列填充采样算法 21
2.1 立方体域的混合整数序列填充采样算法 21
2.1.1 考虑混合整数的拉丁超立方均布性准则 21
2.1.2 基于递归演化和概率跃迁的混合整数序列填充采样算法 24
2.1.3 基于递归演化和概率跃迁的低相关性混合整数序列填充采样算法 30
2.2 非规则域的混合整数序列填充采样算法 35
2.2.1 基于样本点局部密度函数改造的非规则域序列填充采样算法 36
2.2.2 基于逐点采样-排列优化的非规则域混合整数序列填充采样算法 46
2.3 不同采样方法在典型加筋圆柱壳优化设计中的应用 51
2.3.1 典型加筋圆柱壳描述 51
2.3.2 面向加筋圆柱壳轻量化的非规则域混合整数序列填充采样 52
2.3.3 不同初始样本点的加筋圆柱壳轻量化对比分析 54
2.4 本章小结 59
第3章 大型复杂结构的增广径向基函数近似建模方法 61
3.1 增广径向基函数近似模型 61
3.2 基于矩估计的ARBF近似建模方法 64
3.2.1 矩估计方法 64
3.2.2 算例测试 69
3.3 基于分片K折交叉验证的ARBF近似建模方法 75
3.3.1 分片K折交叉验证方法 75
3.3.2 算例测试 83
3.4 面向全局近似的ARBF序列近似建模方法 89
3.4.1 辅助模型与*终模型差异性的偏差-方差分解 90
3.4.2 基于非精确Voronoi图解的ARBF序列近似建模算法 92
3.4.3 算例测试 99
3.5 基于ARBF的多保真度近似建模方法 105
3.5.1 多保真度近似模型 105
3.5.2 基于分片K折交叉验证的ARBF多保真度近似建模方法 107
3.5.3 算例验证 109
3.6 本章小结 113
第4章 基于传统优化算法的加筋圆柱壳舱段优化设计 115
4.1 基于近似建模的网格加筋圆柱壳优化设计 115
4.1.1 网格加筋圆柱壳有限元建模 115
4.1.2 网格加筋圆柱壳稳定性分析 116
4.1.3 基于近似模型的网格加筋圆柱壳优化设计 118
4.2 基于模拟退火算法的集中力扩散舱段优化设计 121
4.2.1 基于多层级加筋的集中力扩散舱段多区域联合设计 121
4.2.2 基于静力分析和工程算法的集中力扩散舱段优化模型 129
4.2.3 集中力扩散舱段多区域联合优化设计 134
4.3 基于协同模拟退火算法的集中力扩散舱段优化设计 139
4.3.1 集中力扩散舱段结构分解策略 140
4.3.2 协同模拟退火算法 142
4.3.3 集中力扩散舱段多子域协同优化设计 148
4.4 本章小结 159
第5章 考虑整体稳定性的大直径加筋圆柱壳优化 161
5.1 大直径薄壁加筋圆柱壳稳定性分析及优化问题建模 161
5.1.1 加筋圆柱壳有限元建模 161
5.1.2 加筋圆柱壳稳定性分析 163
5.1.3 大直径加筋圆柱壳轻量化设计问题建模 164
5.2 基于探索/开发竞争并行采样的ARBF序列近似优化算法 165
5.2.1 序列并行采样策略 165
5.2.2 探索/开发竞争采样机制 172
5.2.3 算法整体框架 173
5.3 考虑整体承载性能的重型运载火箭加筋圆柱壳轻量化设计 174
5.3.1 优化问题建模 174
5.3.2 加筋圆柱壳轻量化设计 175
5.4 本章小结 184
第6章 考虑桁条局部稳定性的大直径加筋圆柱壳优化 185
6.1 桁条稳定性对大直径加筋圆柱壳承载性能影响分析 185
6.1.1 大直径加筋圆柱壳失稳机理 185
6.1.2 桁条失稳模式对大直径加筋圆柱壳承载性能影响分析 186
6.2 基于近似模型的大直径加筋圆柱壳参数灵敏度分析 189
6.2.1 基于ARBF近似模型正交分解的全局灵敏度分析方法 189
6.2.2 数值算例分析 197
6.2.3 大直径加筋圆柱壳参数灵敏度分析 198
6.3 基于搜索空间重构的ARBF序列近似优化算法 203
6.3.1 初始采样空间描述及优化问题重构 204
6.3.2 基于改进双精英种群进化的约束空间开发采样策略 205
6.3.3 基于MIVDO算法的约束空间探索采样策略 208
6.3.4 优化算法整体框架 209
6.4 考虑桁条局部稳定性约束的重型运载火箭加筋圆柱壳轻量化设计 210
6.5 本章小结 218
第7章 基于多保真度近似模型的集中力扩散舱段优化 220
7.1 集中力扩散舱段多保真度优化问题描述 220
试读
第1章 绪论
1.1 结构优化概述
如何获取满足使役功能、低成本和轻量化需求的*优结构设计方案，一直是工程师和学者关注的重点。随着有限元方法、数学规划理论和计算机技术的发展，结构优化技术逐渐成为计算力学领域的重要分支，并在航空航天、车辆、船舶等工程领域得到广泛的关注和应用。下面从结构优化设计和结构优化算法两方面对结构优化技术进行概述。
1.1.1 结构优化设计
集计算力学、数学规划、计算机科学与技术，以及其他工程学科于一体的结构优化设计是现代结构设计领域的重要研究方向，可以为获取创新结构*优化设计提供强有力的先进工具。从设计变量的不同类型来看，结构优化设计通常分为三个层次(图1.1)，即拓扑优化、形状优化、尺寸优化。这与工程设计中概念设计、基本设计和详细设计恰好对应。
图1.1 结构优化中的三个设计层次
1.拓扑优化设计
拓扑优化通过优化设计域内的材料分布，使结构在指定载荷及位移边界条件作用下，满足多种设计需求并达到*优性能。拓扑优化是结构优化领域层级*高的设计手段，优化效果*明显，优化结果更新颖，但是同时也更具挑战性，是结构优化领域公认的极具挑战性的课题。
拓扑优化算法主要分为变密度拓扑优化及其衍生方法、边界演化拓扑优化算法两大类。变密度拓扑优化算法采用有限元网格单元作为拓扑优化的基本设计变量，将连续体结构的拓扑优化问题等价转换为结构有限元模型中离散单元的保留和删除问题，进而通过均匀化思想将离散优化问题转化为便于数值求解的连续优化问题。其衍生方法还包括固体各向同性材料/微结构插值惩罚函数方法(solid isotropic macrostructural/material with penalization method， SIMP)[1， 2]、渐进结构优化(evolutionary structural optimization， ESO)方法[3]、*立-连续-映射(independent continuous and mapping， ICM)方法[4-7]等。变密度拓扑优化算法理论完备、便于数值实现，虽然其设计结果具有边界不光滑、应力计算不准确、难加工制造等不足，但是可通过细分单元、使用新的过滤策略、方法演化或后处理手段进行缓解，因此在工程设计应用中具有一定的优势。特别地，Xie等[3， 8-12]提出的渐进结构优化算法具有算法理念简单、收敛速度快和边界清晰等优点，并成功集成于商用软件Ameba[13]，已在建筑工程、机械工程、航空航天等领域的创新设计方面展现出强大的实际工程解决能力[14，15]。
边界演化拓扑优化算法通过控制实边界的演化过程来实现拓扑优化。其优势在于能够得到光滑和清晰的优化结构边界，并且几何信息可以通过函数表示，易于提取*终构型的几何信息，典型方法包括水平集法、特征驱动法、相场法、拓扑导数法和移动可变性组件/孔洞法等。水平集法[16， 17]通过引入水平集函数的等值面隐式描述结构边界的移动、合并，以及空洞的生成，并通过函数演化更新结构的边界形状，但是该方法的隐式拓扑优化框架会制约几何描述能力，难以实现新型结构对特征尺寸、结构连通性、支撑结构布置等方面的精细控制要求[18]。为此，周莹[19]和Zhang等[20， 21]将拓扑优化算法与CAD(computer aided design，计算机辅助设计)特征设计方法相结合，提出特征驱动的结构拓扑优化算法，在实现*优化创新结构力学性能的同时保留工程结构的特征属性。此外，Zhang等[22-25]和Liu等[26，27]原创性地提出基于可动变形组件(moving morphable components，MMC)法的显式拓扑优化新框架，为破解隐式拓扑优化算法面临的结构特征精细控制难题提供了全新思路，可以充分发挥显式与隐式拓扑描述的优势，实现对结构局部特征几何进行显式、灵活控制。
目前，拓扑优化算法在解决线性问题中已表现出巨大的应用潜力，成功解决了诸如热传导、热力耦合、增材制造、振动抑制、线性屈*、稳健性和可靠性设计等问题。然而，工程实际中存在大量强非线性因素，如几何非线性、材料非线性、接触非线性，现有的拓扑优化算法在优化效率、收敛性和准确度方面都还面临不小的挑战[28-30]。
2.形状优化设计
形状优化通过优化结构内部或边界几何轮廓，使整体结构达到性能*优、缓解边界应力集中或质量*轻。通常情况下，在形状优化过程中，结构的拓扑构型不变，即没有新的孔洞生成。形状优化将结构的几何边界作为关键控制变量，因此主要难点在于如何建立结构几何边界的描述方法，并在优化过程中对其进行调整和形状灵敏度分析[31]。
*早的形状优化策略是将有限元网格的边界节点坐标作为设计变量，实现结构边界几何轮廓的调控。然而，该策略存在两个明显不足：一方面，过多的设计变量使优化效率低下；另一方面，边界节点坐标相互*立变化极易导致有限元网格出现过度畸变和锯齿形状边界。为此，多项式和样条函数相继被用于参数化描述结构边界形状，即采用直线、圆弧、样条*线、二次参数*线、二次*面、柱面来描述结构边界，将结构边界形状优
