2025-11-19T11:34:14.421060

Isogeometric Topology Optimization Based on Topological Derivatives

Teixeira, Krenn, Gangl et al.
Topology optimization is a valuable tool in engineering, facilitating the design of optimized structures. However, topological changes often require a remeshing step, which can become challenging. In this work, we propose an isogeometric approach to topology optimization driven by topological derivatives. The combination of a level-set method together with an immersed isogeometric framework allows seamless geometry updates without the necessity of remeshing. At the same time, topological derivatives provide topological modifications without the need to define initial holes [7]. We investigate the influence of higher-degree basis functions in both the level-set representation and the approximation of the solution. Two numerical examples demonstrate the proposed approach, showing that employing higher-degree basis functions for approximating the solution improves accuracy, while linear basis functions remain sufficient for the level-set function representation.
academic

Isogeometric Topology Optimization Based on Topological Derivatives

基本信息

  • 论文ID: 2509.09236
  • 标题: Isogeometric Topology Optimization Based on Topological Derivatives
  • 作者: Guilherme Henrique TEIXEIRA (Graz University of Technology), Nepomuk KRENN (Johann Radon Institute), Peter GANGL (Johann Radon Institute), Benjamin MARUSSIG (Graz University of Technology)
  • 分类: math.NA cs.CE cs.NA math.OC
  • 发表时间: 2025年10月15日 (arXiv版本v3)
  • 论文链接: https://arxiv.org/abs/2509.09236

摘要

拓扑优化是工程中的重要工具，用于设计优化结构。然而，拓扑变化通常需要重新网格划分，这可能变得具有挑战性。本文提出了一种基于拓扑导数驱动的等几何拓扑优化方法。水平集方法与浸入式等几何框架的结合允许无缝几何更新而无需重新网格划分。同时，拓扑导数提供拓扑修改而无需定义初始孔洞。研究了高阶基函数在水平集表示和解近似中的影响。两个数值算例验证了所提方法，表明采用高阶基函数近似解可提高精度，而线性基函数对水平集函数表示仍然足够。

研究背景与动机

问题定义

本研究要解决的核心问题是传统拓扑优化中的网格重构难题。在结构优化过程中，设计域的拓扑变化（如孔洞的产生、合并或消失）通常需要重新生成计算网格，这不仅计算成本高昂，还可能导致数值不稳定。

重要性

拓扑优化作为结构设计的重要工具，广泛应用于航空航天、汽车制造、建筑工程等领域。解决网格重构问题对于提高优化效率、确保数值稳定性具有重要意义。

现有方法局限性

  1. 传统水平集方法：虽然可以避免显式网格重构，但难以自动生成新的孔洞，依赖于初始几何的孔洞配置
  2. 密度方法：虽然避免了网格问题，但可能产生灰色单元和棋盘格模式
  3. 形状导数方法：只能处理边界变化，无法处理拓扑变化

研究动机

作者提出将等几何分析(IGA)、水平集方法和拓扑导数相结合，旨在：

  • 消除网格重构需求
  • 自动生成新孔洞而无需预定义
  • 利用B样条的高阶光滑性提高精度
  • 简化几何更新过程

核心贡献

  1. 创新方法框架：首次将Amstutz和Andrä的纯拓扑导数驱动方法与等几何分析相结合
  2. 无网格重构优化：通过浸入式等几何框架实现拓扑变化过程中的无缝几何更新
  3. 高阶基函数研究：系统研究了不同多项式阶数对水平集表示和解近似的影响
  4. 数值验证：通过经典算例验证方法有效性，证明线性基函数足以表示水平集而高阶基函数改善解精度

方法详解

任务定义

输入：设计域D、边界条件、载荷、材料属性 输出：最优材料分布Ω ⊂ D 目标：最小化柔度同时满足体积约束 约束：线性弹性方程、边界条件

模型架构

1. 线性弹性问题formulation

在固定域D上求解：

-∇ · (αΩσ(u)) = 0 in D
u = 0 on ΓD
σ(u) · n = τ on ΓN
σ(u) · n = 0 on Γ0

其中材料参数αΩ定义为：

αΩ = {
  αin   in Ω (材料区域)
  αout  in D\Ω (空洞区域，αout ≪ 1)
}

2. 水平集表示

域Ω通过水平集函数ϕ隐式表示：

ϕ(x) < 0 ⟺ x ∈ Ω
ϕ(x) = 0 ⟺ x ∈ ∂Ω  
ϕ(x) > 0 ⟺ x ∈ D\Ω

水平集函数采用B样条离散化：

ϕ(ξ,η) = Σᵢⱼ Bᵢ,d(ξ)Bⱼ,d(η)cᵢⱼ

3. 拓扑导数计算

对于线性弹性问题，拓扑导数的解析表达式为：

dℐ(Ω)(x₀) = {
  -3αin((αout-αin)/(2αout+αin))σ(u):ε(u) - l  if x₀ ∈ Ω
  -3αout((αin-αout)/(2αin+αout))σ(u):ε(u) + l  if x₀ ∈ D\Ω
}

4. 水平集更新算法

采用球面线性插值更新水平集：

ϕⁱ⁺¹ = (1/sin θᵢ)[sin((1-κᵢ)θᵢ)ϕⁱ + sin(κᵢθᵢ)gΩᵢ]

其中θᵢ是当前水平集与广义拓扑导数gΩᵢ之间的L²角度。

技术创新点

1. 浸入式等几何方法

  • 使用固定的B样条背景网格，避免网格重构
  • 通过材料参数αΩ区分材料和空洞区域
  • 保持几何表示的高阶光滑性

2. 纯拓扑导数驱动

  • 不求解Hamilton-Jacobi方程，仅依赖拓扑导数
  • 自动生成新孔洞，无需预定义初始孔洞
  • 简化优化过程，提高计算效率

3. 切割单元处理

  • 使用隐式几何积分库精确计算切割比例
  • 采用标准高斯积分避免小区域积分不稳定
  • 实施平均滤波平滑材料属性跳跃

实验设置

数值算例

选择两个经典拓扑优化基准问题：

  1. 悬臂梁问题：2×1矩形域，左端固定，右端点载荷
  2. 四分之一环问题：内径1外径2的四分之一环，底部固定，顶部点载荷

参数设置

  • 网格：128×128背景网格
  • 材料参数：αin = 1, αout = 10⁻⁴
  • 体积控制：l = 5
  • 滤波参数：γ = 10⁻⁴
  • 材料属性：E = 1, ν = 1/3
  • 收敛条件：εθ = 1, 最大迭代200次

对比设置

研究两种配置：

  1. 相同阶数：水平集和解采用相同多项式阶数(d = p)
  2. 混合阶数：线性水平集(d = 1)与高阶解(p > 1)

评价指标

  • 目标函数值J（柔度）
  • 收敛角度θ
  • 相对面积Ai/A₀
  • 收敛迭代次数

实验结果

悬臂梁问题主要结果

收敛性能

  • 高阶优势：p ≥ 2比p = 1收敛更快，中间阶段目标函数下降更迅速
  • 一致收敛：所有配置收敛到相似最优解(J ≈ 6.97-7.00)
  • 迭代效率：p = 2和p = 4需要相同迭代次数，p = 3在高阶水平集时需要29步vs线性水平集24步

水平集阶数影响

  • 线性水平集(d = 1)产生与高阶水平集(d = p)相似的最终形状
  • 奇偶阶数差异：奇数阶在Greville点共享导致平均化效应，产生轻微差异

四分之一环问题主要结果

局部最优现象

  • 收敛差异：p = 1,2的低阶配置在较高角度停止(θ = 4.72, 9.33)
  • 最优解：p = 3(d = 1)达到最低目标函数值4.12，p = 3(d = 3)为4.13
  • 效率平衡：p = 2(d = 2)在129次迭代达到4.15，效率较高

形状多样性

不同配置产生不同拓扑结构，体现了拓扑优化中局部最优的固有挑战。

关键发现

  1. 高阶解近似有益：提高收敛速度和解精度
  2. 线性水平集充分：d = 1足以准确表示几何边界
  3. 网格无关性：在128×128精细网格下，高阶优势有限
  4. 滤波重要性：平均滤波和PDE滤波对结果质量关键

相关工作

拓扑优化方法分类

  1. 密度方法：SIMP、RAMP等，可能产生灰色单元
  2. 水平集方法：避免灰色单元但需要初始孔洞
  3. 拓扑导数方法：可自动生成孔洞，本文采用此类

等几何拓扑优化

  • 早期工作：Dedè等(2012)相场模型，Shojaee等(2012)径向基函数
  • B样条方法：Jahangiry和Tavakkoli(2017)，Khatibinia等(2020)
  • 本文创新：首次结合纯拓扑导数与等几何方法

浸入式方法

广泛应用于流体力学、固体力学等领域，本文将其引入拓扑优化实现无网格重构。

结论与讨论

主要结论

  1. 成功开发了基于拓扑导数的浸入式等几何拓扑优化框架
  2. 高阶基函数改善解近似但线性基函数足以表示水平集
  3. 方法避免了网格重构并可自动生成孔洞
  4. 在经典基准问题上验证了方法有效性

局限性

  1. 局部最优问题：不同参数配置可能收敛到不同局部最优解
  2. 材料跳跃处理：高阶基函数对材料-空洞界面跳跃的表示优势不明显
  3. 参数敏感性：滤波参数γ等对结果有重要影响
  4. 二维限制：当前仅在2D问题中验证

未来方向

  1. 深入研究材料界面跳跃的高阶表示
  2. 扩展到三维问题
  3. 多材料拓扑优化
  4. 应力约束等更复杂约束条件

深度评价

优点

  1. 方法创新性强：将三种先进技术有机结合，解决实际工程问题
  2. 理论基础扎实：拓扑导数推导严谨，数学公式完整
  3. 实验设计合理：系统研究不同阶数配置，对比全面
  4. 工程实用性高：避免网格重构，适合实际应用

不足

  1. 三维扩展缺失：仅在2D验证，工程应用多为3D
  2. 参数指导不足：对滤波参数等关键参数选择缺乏深入分析
  3. 计算效率分析缺失：未与传统方法进行计算成本对比
  4. 复杂几何验证不足：算例相对简单，复杂工程结构验证缺乏

影响力

  1. 学术贡献：为拓扑优化领域提供新的技术路线
  2. 工程价值：有望应用于航空航天、汽车等领域的结构设计
  3. 可复现性：基于开源GeoPDEs代码，便于复现和扩展

适用场景

  1. 连续体结构优化：特别适合需要频繁拓扑变化的设计问题
  2. 概念设计阶段：快速探索不同拓扑配置
  3. 高精度要求应用：利用B样条高阶光滑性的场景
  4. 自动化设计：无需人工定义初始孔洞的自动优化

参考文献

论文引用了34篇重要文献，涵盖拓扑优化、等几何分析、水平集方法等核心领域的经典和前沿工作，为研究提供了坚实的理论基础。关键参考文献包括Amstutz和Andrä(2006)的拓扑导数方法、Hughes等(2005)的等几何分析基础理论等。