【摘 要】
:
在诸多高新技术领域对于超轻质、多功能等先进结构需求的强劲推动下,结构拓扑优化近年来发展迅猛。经典的拓扑优化方法以有限元模型的密度(像素)为设计变量(密度取0或1分别表示孔洞或材料),因此拓扑优化的数学本质是大规模非线性0-1离散整数规划。为了回避组合复杂性、维度爆炸等困难,传统的方法是将0或1的离散变量松弛为0到1的连续变量,允许中间密度(灰度)单元。中间密度(灰度)单元尽管可以使用材料插值格式进
论文部分内容阅读
在诸多高新技术领域对于超轻质、多功能等先进结构需求的强劲推动下,结构拓扑优化近年来发展迅猛。经典的拓扑优化方法以有限元模型的密度(像素)为设计变量(密度取0或1分别表示孔洞或材料),因此拓扑优化的数学本质是大规模非线性0-1离散整数规划。为了回避组合复杂性、维度爆炸等困难,传统的方法是将0或1的离散变量松弛为0到1的连续变量,允许中间密度(灰度)单元。中间密度(灰度)单元尽管可以使用材料插值格式进行惩罚但是依然存在,因此需要复杂的后处理获得可制造的设计。灰度单元也不利于结构几何/拓扑信息的提取与控制,也会破坏物理模型的真实性从而带来分析困难。事实上,如何消除非0-1中间密度破解灰度难题是拓扑优化领域长期关注的重点挑战性课题。直接求解大规模非线性整数规划可从根本上破解这一难题,然而该方案长久以来被认为是不可能的。本博士论文旨在打破这一瓶颈,从大规模离散变量拓扑优化的数学本源出发,通过研究序列近似整数规划算法的理论、方法和程序实现,构造出通用高性能的离散变量拓扑优化新方法。本论文的主要研究内容为下列四个方面。(一)提出了基于离散变量灵敏度和序列近似整数规划(SAIP)的离散变量拓扑优化新方法。回答了序列近似整数规划中一系列理论问题,成功构造出求解大规模可分离整数规划子问题的正则松弛算法,其计算效率与连续变量的梯度类算法相当,并且该方法具有全局收敛性和渐近强对偶性质。提出了两种适合于离散变量的运动极限策略。该方法成功求解了四种以柔顺性为目标的不同约束不同物理场的拓扑优化问题,优化结果表明该方法统一于数学规划框架下,可高效求解具有多重非线性约束及大量局部约束的离散变量拓扑优化问题。这些理论推导和算法应用充分说明了序列近似整数规划可以有效破解灰度难题,突破已有离散变量方法局限性。(二)如何处理非单调(灵敏度符号不定)的目标或约束函数是离散变量方法求解复杂非柔顺性拓扑优化问题的阻碍。为此,本文构造了基于信赖域的序列近似整数规划方法,即SAIP-TR方法。信赖域约束可直接控制离散设计变量的改变量从而直接控制近似模型的精度,将信赖域约束引入序列近似整数规划可以建立离散变量方法的理性运动极限策略。本文指出了离散变量下信赖域约束是线性约束,构造了基于merit函数的动态调整信赖域半径方案。数值算例表明SAIP-TR方法可以求解以柔顺机构为典型代表的非单调目标函数拓扑优化问题。然后,在黑白拓扑描述下构造了表示结构材料相和空白相的最小尺寸几何约束,并使用SAIP-TR方法成功求解几何约束以消除柔顺机构的点铰和应力集中。最后,SAIP-TR方法成功求解多工况结构响应约束下的无材料用量约束问题。综上所述,SAIP-TR方法极大地扩展了离散变量拓扑优化方法的适用范围并且具有优良的程序通用性。(三)对于离散变量拓扑优化给出的清晰黑白拓扑结构,给出了仅利用节点密度和节点邻域信息的可编程的Euler-Poincaré定理,成功计算了拓扑优化结果的拓扑不变量,即欧拉数和贝蒂数。进一步说明拓扑优化结果的贝蒂数与结构复杂性(孔洞个数)直接相关。结合可编程的Euler-Poincaré定理和三种策略(线性松弛结构复杂性控制约束,骨架约束和控制分比下降)给出了可自由控制结构复杂性的可求解列式。最后在序列近似整数规划框架下使用正则松弛算法成功求解了不同初始孔洞个数,不同目标孔洞个数的二维三维结构的复杂性控制问题。这些研究结果解决了复杂结构拓扑不变量的可计算性问题,构建了计算拓扑学与结构拓扑优化之间的桥梁。(四)为了应用本文的离散变量方法求解大规模拓扑优化问题,提出了将物理多尺度方法和数值多尺度结合的多分辨率离散变量拓扑优化方法。该方法首先利用扩展多尺度有限元计算出粗网格内的数值形函数,然后利用数值形函数构造出表示粗网格与细网格位移转换关系的投影矩阵,最后结合投影矩阵与数值代数中的多重网格方法得到了基于扩展多尺度有限元的多重网格共轭梯度迭代求解器。该求解器避免了传统多分辨率方法的孤岛结构(QR模式),使得本文的离散变量方法可以求解近千万量级自由度的大规模三维拓扑优化问题,获得了具有丰富细节的多分辨率设计。
其他文献
海洋天然气水合物储量巨大并广泛分布于深海沉积层中,具有极高的开发价值。而水合物的开采可能会导致储层失稳,进而引发一系列严重地质灾害,因此,安全保障将是水合物商业化开采的首要前提。储层的各项环境因素是复杂多变的,其中地层温压环境控制着水合物稳定带,并导致了不同区域、不同埋深处的水合物沉积物力学特性的差异性;在不同区域水合物的开采所导致的沉积物强度弱化以及储层变形也将是各异的,给开采工艺方法的设计带来
残余位移是结构震后可修复性评价的重要指标,合理分析并准确估计结构的震后残余位移对于地震损伤控制及性能评估具有重要意义。SDOF体系的残余位移分析是实现结构震后可修复性能控制的基础,现阶段与残余位移相关的研究主要是针对SDOF体系开展的,这些研究在时程分析时通常采用了常阻尼系数假定,事实上,这一阻尼假定并不适用于结构固有阻尼耗能机制的模拟,基于该假定可能会导致不切实际的分析结果。切线刚度正比阻尼已经
在现代造船模式下,数字化技术应用成为船舶产业的研究热点,为船舶设计制造及其自动化带来了新的发展模式和技术支持。数字化造船的实施,以降低生产成本、提高造船效率和质量为目标,成为我国船舶企业在国际船舶市场中提高核心竞争力的重要手段。船舶分段建造精度管理与控制作为数字化造船的关键技术之一,是关于建造精度问题的研究,贯穿于船舶生产全生命周期。全面加强建造全过程精度管理与控制,是保证船体建造质量和快速合拢的
天然气水合物是一种清洁高效的新型能源,其因储量巨大、分布广泛等特点,被认为是21世纪最有开采潜力的接替性能源,是目前世界各国力争的未来能源战略制高点。然而,天然气水合物开采过程中引起的水合物分解,会导致储层胶结结构弱化或消失,进而引起储层强度、刚度等发生巨大变化,可能诱发地层沉降、变形以及井筒失稳等工程地质灾害。阐明水合物沉积物变形微观响应机理,是实现天然气水合物资源安全开采的重要前提。本文以此为
皮肤黑色素瘤是一种由表皮基底层黑素细胞高度病变所导致的恶性肿瘤,其恶性程度高、转移早,已成为皮肤癌中最致命的一种疾病。藤黄酸(Gambogic acid,GA)是从藤黄树胶中提取的一种具有抗肿瘤活性的化合物,其对恶性黑色素瘤有很强的抑制作用。GA的作用靶点多,能通过不同信号途径抑制黑色素瘤活性,且能与多种化疗药物联用达到协同抗肿瘤的效果,其在癌症的预防和治疗中有很好的应用前景。但由于GA的肝首过效
目的:青少年抑郁障碍呈逐年增加的趋势,为家庭和社会带来了巨大的经济压力和负担,是目前重点关注的问题。女性青少年的抑郁障碍发生率更高,且抑郁障碍常常伴随着违纪、攻击、社交退缩、厌学、欺凌等一系列行为问题。而目前对于女性青少年的内化问题研究较多,但是对于女性青少年抑郁障碍患者伴发的行为问题,尤其是外化行为问题的研究较少。另外,青少年的既往经历、父母的教养方式、学校因素与青少年抑郁障碍密切相关,但对于女
近几十年来,高超声速飞行器已经成为各个经济政治大国的军事研究热点。在高超声速飞行器的各种技术当中,热防护技术有着极为重要的地位。一方面,由于飞行器的速度达到了5马赫数以上,飞行器外表面会受到极为严重的气动加热,需要更加耐热且坚实可靠的外层热防护系统;另一方面,飞行器要达到如此高的速度,需要有足够大推力的超燃冲压发动机为其提供动力,这类发动机内部燃料燃烧时也会产生极高的温度,同样需要一套热防护系统和
参数化在计算几何、计算机辅助几何设计(Computer Aided Geometric Design,简称CAGD)、计算机图形学(Computer Graphics,简称CG)、有限元分析(Finite Element Analysis,简称FEA)和等几何分析(Isogeometric Analysis,简称IGA)等领域中应用广泛.实际上,曲线/曲面/体的参数化就是将参数域映射为曲线/曲面/
实际工程中存在着各种难以避免的不确定性,对不确定性问题量化分析的实际需求,使得处理不确定性问题的新方法变得越来越重要,也使“不确定性分析的方法在数值建模技术领域中获得契机”。本文以粘弹性、瞬态非线性热传导和拉压不同模量问题为研究对象,开展模糊不确定性正/反问题的数值建模与求解方法研究。这三个问题的共同特点是确定性问题的求解过程复杂、计算量大、影响计算精度和计算效率的因素多。由于不确定分析过程中一般
近些年来,包括智能手机在内的许多成像设备已经普及人们的工作和生活,并沿伸出丰富的实际应用。受图像捕获过程中的复杂成像环境、图像存储与显示条件等因素的影响,获得的图像信息会退化。例如,有雾的图像,有雨的图像,低分辨率的图像或者破损的图像等等。图像复原的主要目标通过合适的算法,从退化图像复原相应的高质量图像。图像复原涵盖计算机视觉和图像处理中许多经典问题,这些问题在过去的几十年中引起了很多研究与关注。