论文部分内容阅读
Level Set方法是在研究曲线以曲率相关的速度演化时提出来的,是描述曲线(或曲面)演化过程的一种新方法。该方法可用于计算和分析依赖于时间、位置、界面的几何性质以及外部物理特性的速度场作用下的界面的并发运动问题(如界面的融合、交叉、断裂等),在图形图像处理、晶体生长、火焰燃烧、计算流体的数值模拟以及计算机视觉和半导体加工等运动界面追踪领域得到了广泛的应用,已经成为处理运动界面追踪问题的主流方法。本文在求解Hamilton-Jacobi方程的粘性解基础上,结合Godunov单调差分格式、本质无震荡插值(ENO)格式、加权本质无震荡插值(WENO)格式和TVD Runge-Kutta积分,给出了Level Set方程的高分辨率的数值求解算法,并成功处理了多维标量双曲守恒律方程及方程组的激波追踪问题以及非反应激波和多介质流中的爆轰间断的追踪问题。在此基础上构造了窄带Level Set算法,并通过修正局部区域边界附近的Level Set函数值,有效地消除了局部区域边界条件的伪振荡问题,但对运动界面的演化趋势不会造成任何影响。在局部求解Level Set函数时,引入了正则化函数c(φ),避免了计算时考虑局部区域边界的问题。论文将这种局部Level Set方法应用于多物质材料界面变形运动上,并与全局Level Set算法进行了计算时间的比较,充分显示了局部Level Set方法的快速性,并揭示应用Level Set方法进行界面追踪时每一步策略的物理和几何本质。另一方面,论文给出了一种基于Level Set方法的结构拓扑优化设计方法,将结构拓扑构形集中在拓扑边界的处理上。通过引入Level Set函数,将结构边界作为零等值线嵌入到高一维的Level Set函数中。这种方法处理拓扑变形非常灵活,并且非常易于描述结构边界的形状,自然地达到结构拓扑优化设计域中孔洞自然融合的目的,是一种渐近型的结构优化方法。在优化迭代过程中,节点上的Level Set函数值被看作是设计变量,并且为了便于微分计算能够正常进行,对Heaviside函数及其导函数做了光滑化处理,选取了一种正则化的Heaviside函数,避免了均匀化方法中出现的振荡现象(棋盘格式),可以平稳地收敛到最终结果。在此算法基础上,论文又提出了一种基于Level Set算法的局部结构拓扑优化设计方法,在灵敏度分析过程中引入窄带思想,以达到减小计算量,进而加快计算速度的目的。通过对数值算例的分析对比,指出了这种局部优化方法在降低计算量方面的局限性,并提出了一种新的修正的结构拓扑优化方法,以达到大规模减少计算量的目的。本文研究工作主要集中在应用偏微分方程数值求解上,发展了一种局部的Level Set方法,并将这种算法用于追踪界面的演化变形过程,构建了包含应力、应变等信息的灵敏度函数,从而实现了结构的拓扑优化设计。数值算例说明了本文方法确实是快速稳健的,具有重要的理论意义与工程应用价值。