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

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

三杆类桁架材料模型多工况最小柔度优化

研究了多工况结构柔度最小化方法。提出了3杆类桁架连续体材料模型。推导了该材料的刚度矩阵及其导数。通过优化杆件分布场得到优化的类桁架连续体。克服了目前普遍采用单元的"有"和"无"表示结构拓扑的轮廓粗糙、锯齿状边界问题。结点位置的杆件密度和方向作为设计变量,杆件在单元内的密度和方向通过结点位置的数值插值得到,并且在单元内连续变化。由于没有抑制中间密度,完全不存在数值不稳定问题。类桁架连续体由于与杆系结构有明确的对应关系,可以合理地转化杆系结构。选择杆件分布场中的部分杆件可以形成杆系结构。如果再进一步作尺寸和形状优化就可以得到最终的拓扑优化结构。     

基于密度-敏度层次更新策略的三维连续体结构拓扑优化

针对三维连续体结构静力学拓扑优化中灰度单元过多的问题,以结构的柔度最小为目标函数,运用基于SIMP插值方法建立的静力学拓扑优化数学模型,提出一种密度-敏度层次更新策略.以数控插齿机床身为例,运用该更新策略,采用MATLAB编程进行三维拓扑优化.优化结果中灰度单元明显减少,说明了密度-敏度层次更新策略的合理性.     

基于等效静态载荷法的多相材料结构动态拓扑优化设计

为了满足多相材料结构动态性能要求,提出一种基于等效静态载荷法的多相材料结构动态拓扑优化设计方法。采用序列固体各向同性料插值模型惩罚刚度矩阵和质量矩阵,以多个动载荷作用的多相材料结构总动态柔顺度最小化为优化目标,以结构质量和成本为约束条件,构建多相材料结构动态拓扑优化模型,利用等效静态载荷法将多相材料结构动态拓扑优化问题转化为多工况下的多相材料结构静态拓扑优化问题,以降低灵敏度分析的复杂程度;采用移动渐近线算法求解多相材料结构动态拓扑优化问题。数值算例结果表明所提方法是有效性的,与传统的拓扑优化方法相比,基于等效静态载荷法的多相材料结构动态拓扑优化求解时间节省了75%,大大提高了计算效率;与单相材料拓扑优化结果相比,多相材料结构动态拓扑优化设计获得的结构具有更好的动态性能。     

空间独立插值下的节点ICM拓扑优化方法

为了克服连续体结构拓扑优化中的数值不稳定性现象,在独立连续映射法中,采用节点拓扑变量描述单元的有无。单元内任意一点的弹性模量和单元体积建模基于空间独立的拓扑变量场,基于低阶单元形函数和改进过滤方式下的插值函数,得到了各类不同空间插值下的ICM法。通过二维拓扑优化数值算例比较了各类方法下的优化结果。结果表明,混合插值模式下的ICM法能得到较为理想的拓扑优化结果。     

基于物质点描述的双向渐进式拓扑优化方法

:为了克服连续体结构拓扑优化中的数值不稳定现象,定义了表征物质点及其领域有无的物质点拓扑变量,提出基于物质点描述的双向渐进式拓扑优化方法.基于过滤法构造拓扑变量场的插值函数,从而在拓扑优化模型中自然消除了棋盘格现象.为适用于不同单元类型和网格离散形式等,重新定义了灵敏度密度.通过二维数值算例对理论方法进行验证.结果表明:方法在连续体结构拓扑优化设计中具有可行性和有效性.     

考虑混合约束的柔顺机构拓扑优化设计

为了满足制造工艺和静强度要求,提出一种综合考虑最小尺寸控制和应力约束的柔顺机构混合约束拓扑优化设计方法。采用改进的固体各向同性材料插值模型描述材料分布,利用多相映射方法同时控制实相和空相材料结构的最小尺寸,采用最大近似函数P范数求解机构的最大应力,以机构的输出位移最大化作为目标函数,综合考虑最小特征尺寸控制和应力约束建立柔顺机构混合约束拓扑优化数学模型,利用移动渐近算法求解柔顺机构混合约束拓扑优化问题。数值算例结果表明,混合约束拓扑优化获得的柔顺机构能够同时满足最小尺寸制造约束和静强度要求,机构的von Mises等效应力分布更加均匀。     

考虑混合约束的柔顺机构拓扑优化设计

为了满足制造工艺和静强度要求,提出一种综合考虑最小尺寸控制和应力约束的柔顺机构混合约束拓扑优化设计方法。采用改进的固体各向同性材料插值模型描述材料分布,利用多相映射方法同时控制实相和空相材料结构的最小尺寸,采用最大近似函数P范数求解机构的最大应力,以机构的输出位移最大化作为目标函数,综合考虑最小特征尺寸控制和应力约束建立柔顺机构混合约束拓扑优化数学模型,利用移动渐近算法求解柔顺机构混合约束拓扑优化问题。数值算例结果表明,混合约束拓扑优化获得的柔顺机构能够同时满足最小尺寸制造约束和静强度要求,机构的von Mises等效应力分布更加均匀。     

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

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

动力学自然单元法的谐波激励下的连续体结构拓扑优化

自然单元法是一种基于Voronoi图构造形函数的无网格方法,根据自然单元法的优点,提出了动力学自然单元法频率激励载荷下连续体的结构拓扑优化计算。采用各向同性固体微结构惩罚(SIMP)模型,将节点相对密度作为设计变量,建立以动柔度最小为目标函数,频率激励载荷作用下的拓扑优化模型。采用伴随分析法进行灵敏度分析并利用优化准则法对优化模型进行求解。通过数值算例计算,不仅得到了无棋盘格现象的优化结果,而且相比其它无网格方法提高了计算效率,说明该方法具有可行性和优越性。     

任意负泊松比超材料结构设计的功能基元拓扑优化法

基于拓扑优化方法,提出了设计具有任意负泊松比超材料及结构的功能基元拓扑优化法。针对功能基元的不同初始拓扑基结构,包括矩形和三角形初始拓扑基结构,以指定的负泊松比值作为约束条件,以功能基元柔顺度最大化为目标函数,建立了任意负泊松比超材料拓扑优化模型并求解。提取拓扑优化得到的功能基元最优构型,经周期性分布从而形成负泊松比结构。建立优化得到的超材料结构有限元模型,验算了功能基元的泊松比,计算分析了该超材料试件的静、动力学特性。结果表明,该负泊松比效应超材料试件具有较好的承载能力,且在中低频段有较好的减振效果。     

A novel algorithm using an orthotropic material model for topology optimization

This article presents a novel algorithm for topology optimization using an orthotropic material model. Based on the virtual work principle, mathematical formulations for effective orthotropic material properties of an element containing two materials are derived. An algorithm is developed for structural topology optimization using four orthotropic material properties, instead of one density or area ratio, in each element as design variables. As an illustrative example, minimum compliance problems for linear and nonlinear structures are solved using the present algorithm in conjunction with the moving iso-surface threshold method. The present numerical results reveal that: (1) chequerboards and single-node connections are not present even without filtering; (2) final topologies do not contain large grey areas even using a unity penalty factor; and (3) the well-known numerical issues caused by low-density material when considering geometric nonlinearity are resolved by eliminating low-density elements in finite element analyses.     

Topological evolutionary computing in the optimal design of 2D and 3D structures†

An application of evolutionary algorithms and the finite-element method to the topology optimization of 2D structures (plane stress, bending plates, and shells) and 3D structures is described. The basis of the topological evolutionary optimization is the direct control of the density material distribution (or thickness for 2D structures) by the evolutionary algorithm. The structures are optimized for stress, mass, and compliance criteria. The numerical examples demonstrate that this method is an effective technique for solving problems in computer-aided optimal design.     

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.     

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

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

Bidirectional evolutionary structural optimization for nonlinear structures under dynamic loads

This study aims to develop efficient numerical optimization methods for finding the optimal topology of nonlinear structures under dynamic loads. The numerical models are developed using the bidirectional evolutionary structural optimization method for stiffness maximization problems with mass constraints. The mathematical formulation of topology optimization approach is developed based on the element virtual strain energy as the design variable and minimization of compliance as the objective function. The suitability of the proposed method for topology optimization of nonlinear structures is demonstrated through a series of two- and three-dimensional benchmark designs. Several issues relating to the nonlinear structures subjected to dynamic loads such as material, geometric, and contact nonlinearities are addressed in the examples. It is shown that the proposed approach generates more reliable designs for nonlinear structures.     

考虑结构稳定性的变密度拓扑优化方法

将稳定性问题引入传统变密度法中,可实现包含稳定性约束的平面模型结构拓扑优化。以单元相对密度为设计变量,结构柔度最小为目标函数,结构体积和失稳载荷因子为约束条件建立优化问题数学模型,提出了一种考虑结构稳定性的变密度拓扑优化方法。通过分析结构柔度、体积、失稳载荷因子对设计变量的灵敏度,并基于拉格朗日乘子法和Kuhn-Tucker条件,推导了优化问题的迭代准则。同时,利用基于约束条件的泰勒展开式求解优化准则中的拉格朗日乘子。通过推导平面四节点四边形单元几何刚度矩阵的显式表达式,得到了优化准则中的几何应变能。最后,通过算例对提出的方法进行了验证,并与不考虑稳定性的传统变密度拓扑优化方法进行对比,结果表明该方法能显著提高拓扑优化结果的稳定性。研究结果对细长受压结构的优化设计有重要指导意义,对结构的稳定性设计有一定参考价值。