基于自然邻近无网格法的形状优化方法-中国力学学会.docVIP

基于自然邻近无网格法的形状优化方法 王凯，周慎杰，聂志峰，孔胜利 （山东大学南校区机械工程学院，山东济南，250061） 摘要：[2]。 无网格法的理论特点特别适合移动边界及大变形这类问题，学者们纷纷将其应用于形状优化[1~3]。在形状优化中，无网格法的优势不仅在于不需要网格重构，而且可以获得更为准确的灵敏度信息，能够提高形状优化的效率与精度。本文把自然邻近无网格法中控制方程的局部积分“弱”形式[4~6]应用于连续型物质导数法灵敏度分析中，使用直接微分法导出了基于子域积分的连续形式的灵敏度分析公式，采用自然邻近无网格法对其离散求解。使用该灵敏度分析方法，结合约束变尺度序列二次规划法构造了基于自然邻近无网格法的形状优化方法。讨论了数值实施中设计边界的参数化方式、设计速度计算及分析模型更新等问题。算例表明，优化过程不需要网格重构，具有较高的收敛速度，使用较少的设计变量就可以获得良好的优化效果。 2 自然邻近无网格法 自然邻近无网格法通过定义在局部子域上“弱”形式来满足平衡方程和边界条件[]。弹性力学边值问题的积分形式 (1) 式中是权函数是局部子域通常由三部分组成[5]，内部边界Ls、位于整体自由边界上的Lt，以及位于整体本质边界上的Lu。利用高斯散度定理，函数使在Ls上的积分为零： (2) 上式中考虑几何关系与本构关系，将位移作为基本变量，并以non-Sibsonian