滚动信息2

摘要：拓扑优化能够在给定设计空间内实现材料最优分布，但传统方法忽略实际工况中的应力约束与制造约束。本文将有限元分析嵌入拓扑优化迭代流程，建立应力约束下的渐进结构优化模型，提出敏度过滤与边界平滑处理策略，形成一套面向机械结构设计的工程化拓扑优化方法。

关键词：拓扑优化；有限元分析；机械结构

一、引言

机械结构轻量化设计对提升装备性能具有重要意义，拓扑优化作为一种概念设计阶段的结构形态生成技术受到广泛关注。然而传统拓扑优化结果往往存在应力集中、锯齿边界等问题，直接用于工程设计尚有差距。将有限元分析作为力学性能评估手段嵌入优化迭代过程，能够有效解决上述问题。

二、有限元分析与拓扑优化的理论基础

（一）拓扑优化的基本数学模型与求解思路

拓扑优化以连续体结构为对象，在设计域内寻找材料的最优分布形式。变密度法将每个有限单元的相对密度作为设计变量，单元弹性模量与密度之间通过插值函数建立关联。优化目标通常设定为结构柔顺度最小化，即刚度最大化，同时满足体积约束。灵敏度分析计算目标函数对设计变量的导数，为优化迭代提供方向信息。求解过程中采用优化准则法或移动渐近线法更新设计变量，每次更新后单元密度发生变化，结构拓扑逐步演化。

（二）有限元分析在拓扑优化中的核心作用

有限元分析为拓扑优化提供每一步迭代所需的力学响应数据，没有准确的应力、位移和应变信息，优化算法无法判断当前拓扑构型是否满足设计要求。在每次密度更新后，需要重新建立有限元模型，根据当前单元密度重新计算单元刚度矩阵并组装总体刚度矩阵，然后求解平衡方程得到节点位移，再通过本构关系计算单元应力和应变能密度。这些力学量直接用于灵敏度分析，例如柔顺度灵敏度表达式包含单元应变能密度项。在应力约束拓扑优化中，单元应力值决定了该单元是否超过许用应力，从而影响该单元的保留或删除决策。

（三）拓扑优化结果应力集中的产生机理

传统基于柔顺度最小化的拓扑优化只控制整体刚度，不对局部应力水平施加约束，导致优化结果中出现明显的应力集中现象。应力集中通常发生在结构几何突变处，例如尖锐转角、细颈区域或材料突然减少的位置。从力学角度分析，当拓扑优化算法为了降低整体柔顺度而大量移除低应变能密度的单元时，剩余材料需要承受更大的载荷传递路径，局部截面积变小必然导致该处应力升高。敏度过滤技术在消除数值不稳定的同时会模糊材料边界，使优化边界呈现锯齿状，这些微小凹凸进一步加剧了应力集中。优化过程中单元删除是不可逆的，一旦某个单元的密度降至接近零，后续迭代中很难再恢复，这意味着早期决策错误的应力集中点无法被修正。

三、基于有限元分析的拓扑优化改进方法

（一）嵌入应力约束的渐进结构拓扑优化模型构建

在渐进结构优化法的删除准则中引入局部应力约束条件，使材料删除决策基于单元应变能密度与应力许用比值的综合指标。每次有限元分析完成后计算每个单元的冯·米塞斯等效应力，并与平均应力比较，两个指标加权组合成综合灵敏度数，权重系数动态调整。优化初期以应变能密度为主形成主要传力路径，中期增加应力权重消除集中点，后期以应力约束为主导，对高应力区域保留或补充材料。应力严重超限的单元不予删除，必要时进行密度恢复，解决不可逆删除问题。采用移动渐近线算法处理全局收敛控制，确保应力峰值逐步下降并稳定在许用范围内，迭代直至满足体积分数和应力约束双收敛准则。

（二）基于敏度过滤的拓扑边界应力平滑处理策略

将敏度过滤技术扩展应用于边界应力平滑。传统固定半径圆形邻域会模糊材料边界形成锯齿状过渡带。改进策略采用各向异性过滤半径，沿主应力方向加大过滤长度，垂直方向保持较小长度，使边界沿应力流线自然平滑。高应力区域周围的单元获得正向灵敏度增量倾向于保留材料，低应力凹角区域获得负向增量加速删除。每次有限元分析后对节点应力进行局部二次平均，将节点应力值回传替代高斯积分点应力值，消除单元畸变引起的应力伪峰值。对于边界上的孤立材料岛，自动识别并施加惩罚因子，使其在下一迭代中被优先删除。经过多轮处理后边界呈现清晰台阶状轮廓，应力分布趋于均匀。

（三）面向制造约束的拓扑优化结果后处理技术

拓扑优化结果需后处理以满足机械加工工艺要求。以有限元网格为操作对象，首先提取密度场等值线，生成初始几何边界轮廓。采用B样条曲线拟合技术将折线轮廓转化为光滑曲线边界，控制偏差不超过一个单元尺寸。进行特征识别与简化，删除小直径孔洞，加宽或填充狭长通道，将直角改为圆角，圆角半径根据相邻区域应力水平自动计算。对于内部加强筋，分析应力分布判断是否将直线筋改为曲线筋以降低应力集中。几何修复后重新划分网格并进行校核分析，若局部应力超标则回退部分光顺操作保留原始材料。最终输出可用于CAD建模的标准格式文件，包含边界曲线和实体模型。

（四）拓扑优化与有限元分析的迭代耦合实现流程

拓扑优化与有限元分析构成双向数据交换的闭环流程。初始阶段建立三维有限元模型，定义材料属性、边界条件和载荷工况，将设计空间离散化为网格，每个单元对应一个设计变量。主循环中调用有限元求解器进行线性静力分析，计算位移和应力，采用稀疏矩阵存储与预处理共轭梯度法提高效率。提取应变能密度和冯·米塞斯应力传递到优化模块，优化模块自动调整应力约束松弛系数。更新设计变量后重新计算单元弹性模量与总体刚度矩阵，对密度变化小的相邻迭代步采用不完全更新策略。收敛判断同时检查体积分数变化率、目标函数变化率和应力峰值波动幅度，连续多次低于门限时终止迭代，实现力学评估与优化决策的深度耦合。

四、结语

将有限元分析系统性地嵌入拓扑优化流程，建立了应力约束下的结构渐进优化方法，有效抑制了传统拓扑优化中常见的应力集中问题。提出的敏度过滤边界平滑策略和面向制造的后处理技术提升了优化结果的可制造性。

参考文献：

[1]周莹.连续体结构拓扑优化应力约束处理方法研究进展[J].机械工程学报，2020,56(08)：112-123.

[2]龙凯.基于变密度法的结构拓扑优化及其在机械设计中的应用[J].机械设计与研究，2022,38(03)：45-51.

[3]杜义贤.考虑应力约束的渐进结构拓扑优化方法及其算法实现[J].中国机械工程，2024,35(02)：267-275.

杨帆