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

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

Damage detection by the topology design formulation using modal parameters

The feasibility of using the topology design method for structural damage identification is investigated for the first time. The finite element model of an undamaged structure and some point‐frequency response functions of a damaged structure are assumed to be available. To carry out the feasibility study, the topology optimization formulation suitable for structural damage detection is newly set up, where both resonances and anti‐resonances are used as the damage indication modal parameters. An idea to progressively reduce the candidate damaged elements is also developed to improve the accuracy and efficiency of the proposed method. Copyright © 2006 John Wiley & Sons, Ltd.     

Structural shape and topology optimization based on level-set modelling and the element-propagating method

This article introduces the element-propagating method to structural shape and topology optimization. Structural optimization based on the conventional level-set method needs to solve several partial differential equations. By the insertion and deletion of basic material elements around the geometric boundary, the element-propagating method can avoid solving the partial differential equations and realize the dynamic updating of the material region. This approach also places no restrictions on the signed distance function and the Courant–Friedrichs–Lewy condition for numerical stability. At the same time, in order to suppress the dependence on the design initialization for the 2D structural optimization problem, the strain energy density is taken as a criterion to generate new holes in the material region. The coupled algorithm of the element-propagating method and the method for generating new holes makes the structural optimization more robust. Numerical examples demonstrate that the proposed approach greatly improves numerical efficiency, compared with the conventional level-set method for structural topology optimization.     

桁架结构智能布局优化设计

结构的布局优化由于涉及尺寸、形状和拓扑三个层次的综合设计而成为优化问题中的难点,结合桁架结构提出了一个基于多个初始基结构的布局优化方法。以智能生成的、型式多样合理的基结构代替传统模型中的单一基结构,然后从不同基结构下的拓扑优化结果中找出最优设计。在克服传统基结构法有可能限制求解空间而丢失最优解这一局限性的同时,将形状和拓扑优化设计有效分离,降低了求解的难度,并且结合拓扑变化法,实现了桁架结构从选型生成、分析计算到优化设计的一体化智能设计过程。算例表明:利用该文提出的方法进行桁架结构的最优布局设计是可靠有效的。     

结构损伤与荷载的快速同步识别方法研究

本文仅以损伤因子为优化变量,提出一种结构损伤和荷载同步识别的方法。首先通过时域荷载识别的方法将未知荷载转化为损伤因子的函数,将近似荷载作用下的结构响应和实测响应的平方距离作为目标函数,从而降低了需要识别未知参数的数目;然后在目标函数的计算过程中,利用虚拟变形法（VDM）可进行结构快速重分析的思想,快速构造给定损伤因子下系统的脉冲响应,避免每步迭代重新集装系统矩阵,并通过荷载形函数方法进一步提高荷载识别的效率;最后利用二次多项式插值近似结构每个时刻的响应方法和推导对应目标函数的梯度表达式来提高优化搜索的速度。本文利用刚架模型进行数值模拟,准确识别了结构中柱子单元刚度损伤、附加质量以及梁上的未知移动荷载,并通过一个悬臂梁试验进一步验证所提出方法的准确性和可行性。     

Topology optimization of energy absorbing structures with maximum damage constraint

A novel density‐based topology optimization framework for plastic energy absorbing structural designs with maximum damage constraint is proposed. This framework enables topologies to absorb large amount of energy via plastic work before failure occurs. To account for the plasticity and damage during the energy absorption, a coupled elastoplastic ductile damage model is incorporated with topology optimization. Appropriate material interpolation schemes are proposed to relax the damage in the low‐density regions while still ensuring the convergence of Newton‐Raphson solution process in the nonlinear finite element analyses. An effective method for obtaining path‐dependent sensitivities of the plastic work and maximum damage via adjoint method is presented, and the sensitivities are verified by the central difference method. The effectiveness of the proposed methodology is demonstrated through a series of numerical examples. Copyright © 2017 John Wiley & Sons, Ltd.     

基于不完备实测模态数据的结构损伤识别方法研究

在传统的基于模型修正的损伤识别中,由于实测模态信息有限而待识别参数过多,往往导致损伤识别方程出现较大误差,从而限制了该方法在复杂结构中的应用。为解决这一问题,对结构的自由度进行了分解,将损伤结构中模态振型的未测量部分表达为已测量到的模态振型、模态频率以及结构其它参数的函数。将损伤视为结构单元刚度的减小,利用完好结构的计算模态数据以及损伤结构扩充后的实测模态数据,建立了结构的损伤识别方程。运用信赖域优化算法对具有双重约束条件的目标函数进行最小化,识别出了结构各单元的刚度损伤参数。通过两个损伤识别数值仿真算例及实验验证,结果表明,在测点数量有限及测试噪声等不利因素影响下,所提方法只需运用少量的实测模态信息,即可实现结构损伤位置及程度的准确识别,同时算法具有较好的鲁棒性。     

Structural damage detection using time domain responses and an optimization method

An efficient method is proposed to determine the location and severity of structural damage using time domain responses and an optimization method. The time domain responses utilized here are the nodal accelerations measured at the limited points of a structure subjected to an impulse load. The nodal accelerations of the structure are obtained by Newmark time integration method. Firstly, using nodal accelerations extracted for the damaged structure and an analytical model of the structure, an objective function is defined for optimization. Then, the optimization-based damaged detection problem is solved via a differential evolution algorithm for finding the location and severity of damage. In order to assess the accuracy of the proposed method, four numerical examples are considered. Simulation results reveal the efficiency of the method for properly identifying damage with considering measurement noise.     

Simultaneous partial topology and size optimization of a wing structure using ant colony and gradient based methods

This article presents a methodology and process for a combined wing configuration partial topology and structure size optimization. It is aimed at achieving a minimum structural weight by optimizing the structure layout and structural component size simultaneously. This design optimization process contains two types of design variables and hence was divided into two sub-problems. One is structure layout topology to obtain an optimal number and location of spars with discrete integer design variables. Another is component size optimization with continuous design variables in the structure FE model. A multi city-layer ant colony optimization (MCLACO) method is proposed and applied to the topology sub-problem. A gradient based optimization method (GBOM) built in the MSC.NASTRAN SOL-200 module was employed in the component size optimization sub-problem. For each selected layout of the wing structure, a size optimization process is performed to obtain the optimum result and feedback to the layout topology process. The numerical example shows that the proposed MCLACO method and a combination with the GBOM are effective for solving such a wing structure optimization problem. The results also indicate that significant structural weight saving can be achieved.     

Interval analysis of mode shapes to identify damage in beam structures

Various damage detection methods have been proposed by several researchers in the past few decades. Amongst them, the efficiency of mode shapes in detecting damage has been demonstrated by many researchers when further processed. In most cases, the processing involves expansion or reduction of the mode shape data. However, vital information that are damage-prints are often lost during processing of the mode shape data. In addition, most of these processes involve long and complex computation, thus, leading to inaccurate damage identification. In this study, a simple and fast damage identification technique is proposed to identify damage in beam structures. Interval analysis is applied to the mode shapes of a beam structure in the damaged and undamaged states. The interval situations of each of the beam's segment via mode shape are derived to obtain the upper and lower bounds and the derived bounds are compared. To establish a relationship for identify the damaged point, a possibility of damage existence is defined for each segment of the beam structure. The mode shape increment is defined as the increase in the mode shape value. Furthermore, a damage measure index that provide enhance damage information is obtained as the product of the possibility of damage existence and mode shape increment. A numerical model of a simply supported steel beam is applied to demonstrate this method by imposing damage through thickness reduction of elements in segments. In addition, a parametric analysis is carried out to evaluate noise effect by considering varying damage severities and different noise levels. The results showed that this method is simple and provides considerable accurate results.     

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.     

基于虚载荷变量的无网格法形状优化的研究

将选择施加在"虚结构"控制点上的虚载荷作为形状优化的设计变量,并将它与无网格Galerkin法相结合来开展结构形状优化研究,采用罚函数法来施加边界条件,通过直接微分法建立了结构形状优化的离散型灵敏度分析算法,利用无网格法研究了节点坐标关于设计变量导数的计算。所提出的算法简单明了,它不仅解决了网格的畸变问题,而且简化了优化模型和迭代流程,并可使结构的受力特性得到进一步的改善。最后用2个工程实例验证了所建立的算法,并得到了形状优化结果。     

基于水平集描述的结构拓扑与构件布局联合优化新方法

该文提出了一种结构拓扑与内嵌构件布局联合优化的新颖方法。这种方法突出的特点是利用水平集函数隐式地描述不规则的构件形状,因此可以非常方便地处理构件之间的互不覆盖约束条件。数值算例表明:较之文献中已有的方法,该文算法能够以更小的计算量有效地实现结构拓扑与内嵌构件布局的联合优化。     

Isogeometric topology optimization for continuum structures using density distribution function

This paper will propose a more effective and efficient topology optimization method based on isogeometric analysis, termed as isogeometric topology optimization (ITO), for continuum structures using an enhanced density distribution function (DDF). The construction of the DDF involves two steps. (1)  Smoothness: the Shepard function is firstly utilized to improve the overall smoothness of nodal densities. Each nodal density is assigned to a control point of the geometry. (2) Continuity: the high-order NURBS basis functions are linearly combined with the smoothed nodal densities to construct the DDF for the design domain. The nonnegativity, partition of unity, and restricted bounds [0, 1] of both the Shepard function and NURBS basis functions can guarantee the physical meaning of material densities in the design. A topology optimization formulation to minimize the structural mean compliance is developed based on the DDF and isogeometric analysis to solve structural responses. An integration of the geometry parameterization and numerical analysis can offer the unique benefits for the optimization. Several 2D and 3D numerical examples are performed to demonstrate the effectiveness and efficiency of the proposed ITO method, and the optimized 3D designs are prototyped using the Selective Laser Sintering technique.     

基于人工鱼群算法的结构模型修正与损伤检测

提出结构模型修正结构损伤检测的人工鱼群算法。将结构模型修正与结构损伤检测结构动力学逆问题转化为约束优化数学问题,并尝试用人工鱼群算法求解。介绍人工鱼群算法基本原理,定义关键参数并描述觅食、聚群、追尾及随机等行为;据模型修正原理利用结构损伤前后模态特性数据定义优化问题目标函数;通过两层刚架不同损伤工况数值仿真、三层框架试验数据验证方法的有效性。结果表明,基于人工鱼群算法的结构模型修正与损伤检测方法能有效修正结构有限元模型,在不同噪声水平及各种结构损伤工况下不仅能准确定位结构损伤且能精确识别损伤程度。     

A level set method for shape and topology optimization of large‐displacement compliant mechanisms

A parameterization level set method is presented for structural shape and topology optimization of compliant mechanisms involving large displacements. A level set model is established mathematically as the Hamilton–Jacobi equation to capture the motion of the free boundary of a continuum structure. The structural design boundary is thus described implicitly as the zero level set of a level set scalar function of higher dimension. The radial basis function with compact support is then applied to interpolate the level set function, leading to a relaxation and separation of the temporal and spatial discretizations related to the original partial differential equation. In doing so, the more difficult shape and topology optimization problem is now fully parameterized into a relatively easier size optimization of generalized expansion coefficients. As a result, the optimization is changed into a numerical process of implementing a series of motions of the implicit level set function via an existing efficient convex programming method. With the concept of the shape derivative, the geometrical non‐linearity is included in the rigorous design sensitivity analysis to appropriately capture the large displacements of compliant mechanisms. Several numerical benchmark examples illustrate the effectiveness of the present level set method, in particular, its capability of generating new holes inside the material domain. The proposed method not only retains the favorable features of the implicit free boundary representation but also overcomes several unfavorable numerical considerations relevant to the explicit scheme, the reinitialization procedure, and the velocity extension algorithm in the conventional level set method. Copyright © 2008 John Wiley & Sons, Ltd.     

Damage identification in framed structures using natural frequencies

A novel procedure for damage identification of framed structures is proposed, where both the location and the extent of structural damage in framed structures can be correctly determined using only a limited number of measured natural frequencies. No knowledge of the modal shapes of the damaged structure is required. On the basis of the characteristic equations for the original and the damaged structure, a set of equations is generated. Two computational techniques, the direct iteration (DI) technique and the Gauss–Newton least-squares (GNLS) technique, are utilized to determine structural damage from the derived equations. Finally, different numerical examples are used to demonstrate the effectiveness of the proposed method. © 1997 John Wiley & Sons, Ltd.     

Robust concurrent topology optimization of multiscale structure under single or multiple uncertain load cases

Concurrent topology optimization of macrostructure and material microstructure has attracted significant interest in recent years. However, most of the existing works assumed deterministic load conditions, thus the obtained design might have poor performance in practice when uncertainties exist. Therefore, it is necessary to take uncertainty into account in structural design. This article proposes an efficient method for robust concurrent topology optimization of multiscale structure under single or multiple load cases. The weighted sum of the mean and standard deviation of the structural compliance is minimized and constraints are imposed to both the volume fractions of macrostructure and microstructure. The effective properties of the microstructure are calculated via the homogenization method. An efficient sensitivity analysis method is proposed based on the superposition principle and orthogonal similarity transformation of real symmetric matrices. To further reduce the computational cost, an efficient decoupled sensitivity analysis method for microscale design variables is proposed. The bidirectional evolutionary structural optimization method is employed to obtain black and white designs for both macrostructure and microstructure. Several two-dimensional and three-dimensional numerical examples are presented to demonstrate the effectiveness of the proposed approach and the effects of load uncertainty on the optimal design of both macrostructure and microstructure.     

Topology optimization of particle‐matrix composites for optimal fracture resistance taking into account interfacial damage

This paper presents a topology optimization framework for optimizing the fracture resistance of two‐phase composites considering interfacial damage interacting with crack propagation through a redistribution of the inclusions phase. A phase field method for fracture capable of describing interactions between bulk brittle fracture and interfacial damage is adopted within a diffuse approximation of discontinuities. This formulation avoids the burden of remeshing problem during crack propagation and is well adapted to topology optimization purpose. Efficient design sensitivity analysis is performed by using the adjoint method, and the optimization problem is solved by an extended bidirectional evolutionary structural optimization method. The sensitivity formulation accounts for the whole fracturing process involving crack nucleation, propagation, and interaction, either from the interfaces and then through the solid phases, or the opposite. The spatial distribution of material phases are optimally designed using the extended bidirectional evolutionary structural optimization method to improve the fractural resistance. We demonstrate through several examples that the fracture resistance of the composite can be significantly increased at constant volume fraction of inclusions by the topology optimization process.     

基于性能的钢框架结构失效模式识别及优化

该文定义了损伤和滞回耗能两个性能指标,以增量动力分析(IDA)方法为基础,提出基于性能的钢框架结构失效模式识别方法,并以性能指标为目标函数,以构件截面尺寸作为变量,建立钢框架结构失效模式多目标优化方法。在多条地震波作用下,对一个20层benchmark钢框架结构进行了失效模式识别与优化分析,结果表明,以损伤和滞回耗能作为评价指标的基于性能的钢框架结构失效模式识别方法,能有效识别最不利地震作用下的结构失效模式,基于性能的失效模式多目标优化方法能够显著提高结构整体的抗震性能。