连续体结构的变密度拓扑优化方法研究

为了实现使连续体结构的体积约束和柔顺度最小的拓扑优化及解决采用经典变密度法引起的结构优化结果存在如灰度单元、棋盘格等数值不稳定问题，提出了一种新的拓扑优化方法。首先，采用改进的固体各向同性材料惩罚法作为材料插值方案，建立结构拓扑优化模型；其次，通过引入基于高斯权重函数的敏度过滤法和设计新灰度单元抑制算子来解决数值不稳定问题；最后，借助优化准则法求解优化模型。通过算例分析可知：新策略可以改进拓扑优化方法；新的拓扑优化方法具有收敛速度较快、能更好地获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等优势。研究结果为其他连续体结构的拓扑优化研究提供了新思路。     

结构拓扑优化中敏度过滤技术研究

为了抑制连续体结构拓扑优化结果中的棋盘格和灰度单元问题,借鉴粒子群优化算法中粒子状态的更新方法,提出一种改进的敏度更新技术.以结构的柔度最小为优化目标,构建了基于固体各项同性微惩罚结构的结构拓扑优化模型,根据结构的力学响应分析,采用优化准则法进行设计变量更新,进行载荷作用下二维连续体结构的拓扑优化设计,得到了材料在设计域内的最优分布.通过与已有敏度过滤技术的对比分析,验证了文中方法的正确性和有效性．     

自由阻尼层结构阻尼材料配置优化的拓扑敏度法

提出阻尼胞单元和阻尼拓扑敏度等概念，建立了基于阻尼拓扑敏度综合评价的阻尼材料拓扑优化准则，并用于自由阻尼层结构振动控制中阻尼材料的配置优化。建立待控结构阻尼材料布局的拓扑基结构，计算各单元的阻尼拓扑敏度。再建立考虑重量目标及结构频响峰值约束的阻尼材料配置拓扑优化模型。根据所提出的阻尼材料拓扑优化准则，求解上述配置优化问题，确定阈值和各单元拓扑值。并用若干典型结构算例，验证所提出方法的正确性，讨论了阻尼材料布局拓扑基结构的规模与优化效率的关系。     

基于多分辨率-多边形单元建模策略的多材料结构动刚度拓扑优化方法

利用多边形有限单元的高精度求解优势,融合多分辨率拓扑优化方法,实现粗糙位移网格条件下的高分辨率构型设计,由此提出多材料结构动刚度问题的拓扑优化方法。将多边形单元(位移场求解单元)劈分为精细的小单元,构造设计变量与密度变量的重叠网格,形成多分辨率-多边形单元的优化建模策略;以平均动柔度最小化为目标和多材料的体积占比为约束,建立多材料结构的动力学拓扑优化模型,通过HHT-α方法求解结构动响应,采用伴随变量法推导目标函数和约束的灵敏度表达式,利用基于敏度分离技术的ZPR设计变量更新方案构建多区域体积约束问题的优化迭代格式;通过典型数值算例分析优化方法的可行性和动态载荷作用时间对优化结果的影响机制。     

基于复合模态阻尼比的阻尼板渐进法减振优化

研究阻尼板黏弹减振优化方法,为阻尼板设计奠定一定基础.建立了阻尼板运动方程及其求解模型,推导出模态阻尼比敏度算式;构建了以多模态复合模态阻尼比最大、黏弹材料用量及频率变动最小为约束的阻尼板优化模型,采用拓扑渐进法求解优化模型;引入绕单元敏度均化技术解决优化迭代中的棋盘格问题;在引入虚拟单元及逻辑删除技术基础上,编程实现了渐进优化算法;对任意四边形板进行了优化仿真,实现了该板复合阻尼比最大化设计;引入阻尼比体积密度评析拓扑优化效能.研究表明,在阻尼板设计时若实现其多模态阻尼比最大化,则能使板获得良好减振效果.     

考虑泊松效应的材料/结构一体化设计方法

为实现含有不同泊松比组分复合材料的优化设计,并考虑宏观结构及复杂的边界条件,提出了考虑泊松效应的材料/结构一体化设计方法,其显著特征在于不同组分材料中引入了泊松比插值,假设宏观结构由周期性排列的复合材料组成,复合材料含两种各向同性且泊松比不同的组分材料,以静态问题中柔顺度最小化或动态问题中特征值最大化为目标以及宏微观体积比为约束建立了拓扑优化模型。采用均匀化理论预测了复合材料等效性能,推导了目标函数对宏微观密度变量的敏度表达式。分别采用密度过滤和敏度过滤来消除宏微观拓扑优化中的不稳定性现象。采用优化准则法分别更新宏观、微观密度变量,考察了微观体积比和组分材料泊松比参数对优化结果的影响。三维数值算例结果表明所提出的一体化方法具有可行性和优越性。     

基于对数型Heaviside近似函数作为过滤函数的动力响应结构拓扑优化ICM方法

应用ICM(Independent Continuous and Mapping)方法, 建立了以重量最小为目标函数, 以连续频率带或离散点频率的简谐激励下的响应振幅为约束的拓扑优化模型. 引入了对数型Heaviside近似函数作为过滤函数, 并做了敏度分析, 利用对偶二次规划进行优化模型的求解, 并运用敏度过滤的方法处理动力响应数值不稳定的问题. 数值算例比较了利用对数型函数和幂函数作为过滤函数时对拓扑结构的影响, 结果显示利用对数型函数较幂函数结构优化迭代次数更少, 收敛更快.     

基于优化准则法的自由阻尼材料拓扑优化

基于各向正交惩罚材料密度法法,建立了以自由阻尼结构模态阻尼比最大化为目标,以阻尼材料体积为约束条件,以阻尼材料单元相对密度为设计变量的拓扑优化数学模型。推导了模态阻尼比对阻尼单元相对密度的灵敏度和设计变量的更新准则,基于优化准则算法用MSC.Nastran的Direct Matrix Abstraction Program语言编制了通用的阻尼材料拓扑优化程序。以一铝板-声腔耦合系统为对象,以结构一阶弯曲模态阻尼比最大为目标,利用该拓扑优化程序和有限元分析方法对铝板上的阻尼材料进行优化布置,并用试验进行了验证。把铝板全阻尼处理和拓扑优化后部分阻尼处理的降噪效果进行对比,在仿真环境中,参考点处声压分别下降了110.6Pa和107.7Pa,在实验条件下,参考点处声压分别下降了22.47Pa和20.91Pa,从而验证了优化方法的有效性。     

考虑重叠过滤及稳定性约束的桁架拓扑优化方法

从工程实际的角度来说,一般不允许结构形式中包含重叠单元。目前桁架拓扑优化的基结构法在选定基节点的情况下一般不建重叠单元,这导致可行域缩小,使优化不能找到更优解,人为增删杆件缺乏科学依据。针对该问题,该文对重叠给出准确的数学描述,建立包含重叠杆的基结构,利用Heaviside函数将拓扑变量连续化处理,使之在优化过程中可以获取目标函数、约束函数的敏度信息,同时考虑在拓扑优化中加入基频约束以避免出现机构,并加入稳定性约束防止出现压杆失稳,通过优化模型实现重叠过滤。最后通过两个案例计算证明可以找到更优解,验证了该方法的有效性。     

考虑设计相关载荷的稳态热传导拓扑优化设计

针对稳态热传导问题,以结构散热弱度最小为目标,建立了连续体传热结构的拓扑优化模型和方法,给出了相应的算例。优化方法中分别建立了设计相关载荷和非相关载荷的灵敏度列式,采用Rational Approximation of Material Properties (RAMP)方法对材料密度进行惩罚,利用优化准则法控制设计目标与材料分布,以敏度过滤技术抑制棋盘格效应。算例的结果直观显示了设计相关载荷和非设计相关载荷以及复合载荷对结构拓扑构型的影响规律,表明了该文考虑设计相关载荷的稳态热传导结构拓扑优化方法的合理性。     

Explicit level‐set‐based topology optimization using an exact Heaviside function and consistent sensitivity analysis

This paper presents a level‐set‐based topology optimization method based on numerically consistent sensitivity analysis. The proposed method uses a direct steepest‐descent update of the design variables in a level‐set method; the level‐set nodal values. An exact Heaviside formulation is used to relate the level‐set function to element densities. The level‐set function is not required to be a signed‐distance function, and reinitialization is not necessary. Using this approach, level‐set‐based topology optimization problems can be solved consistently and multiple constraints treated simultaneously. The proposed method leads to more insight in the nature of level‐set‐based topology optimization problems. The level‐set‐based design parametrization can describe gray areas and numerical hinges. Consistency causes results to contain these numerical artifacts. We demonstrate that alternative parameterizations, level‐set‐based or density‐based regularization can be used to avoid artifacts in the final results. The effectiveness of the proposed method is demonstrated using several benchmark problems. The capability to treat multiple constraints shows the potential of the method. Furthermore, due to the consistency, the optimizer can run into local minima; a fundamental difficulty of level‐set‐based topology optimization. More advanced optimization strategies and more efficient optimizers may increase the performance in the future. Copyright © 2012 John Wiley & Sons, Ltd.     

Reliability‐based topology optimization of structures under stress constraints

In this paper, we propose an approach for reliability‐based design optimization where a structure of minimum weight subject to reliability constraints on the effective stresses is sought. The reliability‐based topology optimization problem is formulated by using the performance measure approach, and the sequential optimization and reliability assessment method is employed. This strategy allows for decoupling the reliability‐based topology optimization problem into 2 steps, namely, deterministic topology optimization and reliability analysis. In particular, the deterministic structural optimization problem subject to stress constraints is addressed with an efficient methodology based on the topological derivative concept together with a level‐set domain representation method. The resulting algorithm is applied to some benchmark problems, showing the effectiveness of the proposed approach.     

A new hole insertion method for level set based structural topology optimization

Structural shape and topology optimization using level set functions is becoming increasingly popular. However, traditional methods do not naturally allow for new hole creation and solutions can be dependent on the initial design. Various methods have been proposed that enable new hole insertion; however, the link between hole insertion and boundary optimization can be unclear. The new method presented in this paper utilizes a secondary level set function that represents a pseudo third dimension in two‐dimensional problems to facilitate new hole insertion. The update of the secondary function is connected to the primary level set function forming a meaningful link between boundary optimization and hole creation. The performance of the method is investigated to identify suitable parameters that produce good solutions for a range of problems. Copyright © 2012 John Wiley & Sons, Ltd.     

多根非线性方程组求解的探路者灰狼算法

针对传统算法在多根非线性方程组求解时依赖初始值的选定,求解个数不完全,求解精度不高的问题,提出了一种结合探路者算法的灰狼优化算法 (PGWO)。由于灰狼优化算法存在后期收敛速度慢等问题,结合了探路者算法,根据探路者中跟随者的更新机制对灰狼个体的位置进行改变,进而平衡算法的全局搜索和局部搜索能力。通过 9 组多根非线性方程组的仿真实验结果和其他群智能算法进行比较,实验结果表明 PGWO 算法提高了多根非线性方程组求解的精度,在求解个数上得到明显提升,进而说明了算法的有效性。     

基于等几何裁剪分析的拓扑与形状集成优化

提出了将设计和分析、拓扑与形状优化集成的思想,探索了基于等几何裁剪分析的拓扑与形状集成优化设计算法,该方法统一了结构优化的计算机辅助设计、计算机辅助工程分析和优化设计的模型,基于B样条的等几何裁剪分析既能准确表达几何形状,又可以用裁剪面分析方便处理任意复杂拓扑优化问题,由裁剪选择标准确定合理的拓扑结构变动方向,结构变动时无需重新划分网格,设计结果突破初始设计空间的限制,还可方便优化形状。建立了等几何裁剪灵敏度分析的计算方法,给出了等几何裁剪分析拓扑与形状集成优化算法,通过典型实例表明所用方法的正确性和有效性。     

Topology optimization with incompressible materials under small and finite deformations using mixed u/p elements

This paper focuses on topology optimization utilizing incompressible materials under both small‐ and finite‐deformation kinematics. To avoid the volumetric locking that accompanies incompressibility, linear and nonlinear mixed displacement/pressure (u/p) elements are utilized. A number of material interpolation schemes are compared, and a new scheme interpolating both Young's modulus and Poisson's ratio (E ‐ν interpolation) is proposed. The efficacy of this proposed scheme is demonstrated on a number of examples under both small‐ and finite‐deformation kinematics. Excessive mesh distortions that may occur under finite deformations are dealt with by extending a linear energy interpolation approach to the nonlinear u/p formulation and utilizing an adaptive update strategy. The proposed optimization framework is demonstrated to be effective through a number of representative examples.     

A moving superimposed finite element method for structural topology optimization

Level set methods are becoming an attractive design tool in shape and topology optimization for obtaining efficient and lighter structures. In this paper, a dynamic implicit boundary‐based moving superimposed finite element method (s‐version FEM or S‐FEM) is developed for structural topology optimization using the level set methods, in which the variational interior and exterior boundaries are represented by the zero level set. Both a global mesh and an overlaying local mesh are integrated into the moving S‐FEM analysis model. A relatively coarse fixed Eulerian mesh consisting of bilinear rectangular elements is used as a global mesh. The local mesh consisting of flexible linear triangular elements is constructed to match the dynamic implicit boundary captured from nodal values of the implicit level set function. In numerical integration using the Gauss quadrature rule, the practical difficulty due to the discontinuities is overcome by the coincidence of the global and local meshes. A double mapping technique is developed to perform the numerical integration for the global and coupling matrices of the overlapped elements with two different co‐ordinate systems. An element killing strategy is presented to reduce the total number of degrees of freedom to improve the computational efficiency. A simple constraint handling approach is proposed to perform minimum compliance design with a volume constraint. A physically meaningful and numerically efficient velocity extension method is developed to avoid the complicated PDE solving procedure. The proposed moving S‐FEM is applied to structural topology optimization using the level set methods as an effective tool for the numerical analysis of the linear elasticity topology optimization problems. For the classical elasticity problems in the literature, the present S‐FEM can achieve numerical results in good agreement with those from the theoretical solutions and/or numerical results from the standard FEM. For the minimum compliance topology optimization problems in structural optimization, the present approach significantly outperforms the well‐recognized ‘ersatz material’ approach as expected in the accuracy of the strain field, numerical stability, and representation fidelity at the expense of increased computational time. It is also shown that the present approach is able to produce structures near the theoretical optimum. It is suggested that the present S‐FEM can be a promising tool for shape and topology optimization using the level set methods. Copyright © 2005 John Wiley & Sons, Ltd.     

水下耐压结构拓扑优化设计方法探究

探究了拓扑优化设计方法在水下耐压结构设计中的应用。与固定载荷作用下结构的优化设计相比,此类问题需要正确地确定压力作用面。在拓扑优化过程中,利用变密度法得到的中间结构拓扑实际上可以看成是灰度图。基于此,提出了基于图像分割技术的压力加载面搜索方法,并利用距离正规化水平集方法(DRLSE)检测图像边界。利用数值算例验证了方法的有效性,并研究了静水压力作用下结构的拓扑优化设计问题。在给定材料约束的前提下,研究了不同边界条件下耐压壳体的最小柔顺度及最优结构拓扑形式。优化结果说明了该方法在多球交接耐压壳结构形式优化设计及复杂边界条件下耐压结构新形式探索中的工程应用价值。     

A bidirectional evolutionary structural optimization algorithm for mass minimization with multiple structural constraints

A bidirectional evolutionary structural optimization algorithm is presented, which employs integer linear programming to compute optimal solutions to topology optimization problems with the objective of mass minimization. The objective and constraint functions are linearized using Taylor's first-order approximation, thereby allowing the method to handle all types of constraints without using Lagrange multipliers or sensitivity thresholds. A relaxation of the constraint targets is performed such that only small changes in topology are allowed during a single update, thus ensuring the existence of feasible solutions. A variety of problems are solved, demonstrating the ability of the method to easily handle a number of structural constraints, including compliance, stress, buckling, frequency, and displacement. This is followed by an example with multiple structural constraints and, finally, the method is demonstrated on a wing-box, showing that topology optimization for mass minimization of real-world structures can be considered using the proposed methodology.