Shape optimization approaches to free‐surface problems

Three different reformulations of a free‐surface problem as shape optimization problems are considered. These give rise to three different cost functionals that apparently have not been exploited in literature. The shape derivatives of the cost functionals are explicitly determined. The gradient information is combined with the boundary variation method in a preconditioned steepest descent algorithm to solve the shape optimization problems. Numerical results that compare the performance of the proposed cost functionals are presented. Copyright © 2014 John Wiley & Sons, Ltd.     

Application of dual kriging to structural shape optimization based on the boundary contour method

Summary The paper presents an approach in which the coupling of dual kriging and the boundary contour method (BCM) is applied to structural shape optimization problems in mechanical engineering design. The problems consist of optimizing the shape of an elastic body, which requires minimizing an objective function subject to some given constraints, such as those of displacement, stress or manufacturing. The originality of the present work is involved with the use of two novel methods that are combined here to solve structural shape optimization problems. The first one, called dual kriging, is a general, versatile interpolation and geometric modeling tool. The second one is a new variant of the boundary element method (BEM), called the BCM, which achieves a further reduction in dimensionality of analysis problems. Based on the advantages of these two methods, the coupling approach presented here is expected to offer an effective as well as a straightforward manner for solving shape optimal design problems. Received 18 December 1997; accepted for publication 21 April 1998     

基于显式几何更新算法的地下管道自动形状优化

基于计算力学中的结构优化思想,应用一种新型的显式几何更新算法,自行编制C++程序,实现地下管道形状设计的自动优化。管道内的流体假设为牛顿不可压缩流,并考虑惯性项。优化区域主要为管道竖直方向和水平方向的过渡段。形状优化的设计变量是几何边界的有限元节点坐标,优化目标是实现流体黏性能耗散的最小化。优化过程基于形状梯度,即通过形状敏感度分析来求解目标函数相对于设计变量的偏导数。所使用的显式几何更新算法既可以通过网格清晰描述形状,也可以大范围地自动更新网格。详细介绍了地下管道自动形状优化过程的关键步骤。通过数值算例探讨了不同注入速度、密度和黏度对其最优形状的影响。     

悬式绝缘子的形状优化

本文将带自适应运动极限的序列线性规划法应用于悬式绝缘子瓷头的形状优化问题中。建立了优化的数学模型;给出了相应的算法;分析了优化计算的结果;为悬式绝缘子的形状优化设计提供了可行的有效手段。     

双向固定网格渐进结构优化方法

近年来发展的渐进结构优化方法是一种有前途的结构拓扑和形状优化方法.本文在渐进结构优化方法的框架内建立了统一敏感度的概念,并基于固定网格有限元技术,发展了一种新的增加材料技术,提出了双向固定网格渐进结构优化方法.将该方法应用于复合材料壳结构开孔形状优化,以孔周等Tsai-Hill强度值作为优化目标,可以得到合理的最优解,证明了双向固定网格渐进结构优化方法的适用性.不同的初始点能得到几乎相同的最优解,展示了本文方法良好的全局最优性.     

考虑动荷载响应的结构形状优化设计研究

以产生相同位移场为基础,将动荷载转化为一系列的等效静荷载,然后将这些等效静荷载作为多个载 荷工况进行有限元优化设计分析. 通过几个算例的验证,表明这种基于动荷载等效转换的算 法能够反映动荷载对结构的动力响应,对具有大规模自由度的结构形状优化问题是有效的.     

Optimal Shape of Column with Own Weight: Bi and Single Modal Optimization

The optimality conditions, via Pontryagin’s maximum principle, in the case of bimodal optimization of columns are derived. When these conditions are applied to the stability of a compressed column with own weight, the problem of determining the optimal cross–sectional area function is reduced to the solution of a nonlinear boundary value problem. Two specific problems are analyzed in detail. In Problem 1, that is new, the shape of a heavy compressed column with clamped ends stable against buckling and having minimal volume is determined. In Problem 2, formulated by Keller and Niordson, optimal shape of a vertical column with one end clamped the other end free is determined.     

两类变量综合处理的结构形状优化设计方法

本文针对截面变量为离散变量为连续变量的结构优化问题提出了一种优化设计的方法,首先将单元内力作一阶近似,利用凝聚函数多约束问题转化了单约束问题。在解过程中,把定义在连续区间上的形状变量看成是在一些离散以值的离工用变量,然后将两类变量统一考虑并利用相对差商法求解。将该算法应用于几个经典的结构优化算例,运算结果显示了该方法是可行的,优化结果也比较满意。     

Numerical solution of steady free‐surface flows by the adjoint optimal shape design method

Numerical solution of flows that are partially bounded by a freely moving boundary is of great importance in practical applications such as ship hydrodynamics. Free‐boundary problems can be reformulated into optimal shape design problems, which can in principle be solved efficiently by the adjoint method. In this work we investigate the suitability of the adjoint shape optimization method for solving steady free‐surface flows. The asymptotic convergence behaviour of the method is determined for free‐surface flows in 2D and 3D. It is shown that the convergence behaviour depends sensitively on the occurrence of critical modes. The convergence behaviour is moreover shown to be mesh‐width independent, provided that proper preconditioning is applied. Numerical results are presented for 2D flow over an obstacle in a channel. The observed convergence behaviour is indeed mesh‐width independent and conform the derived asymptotic estimates. Copyright © 2003 John Wiley & Sons, Ltd.     

Optimization of the shapes of obstacles in jet-separation flow

The model of an ideal incompressible fluid is used to study the solvability of optimal control problems for the shape of a nozzle which discharges free-boundary fluid flow with and without accounting for gravity (internal aerodynamics) and shape optimization problems for an obstacle with jet separation (external aerodynamics). The qualitative properties of such flows are studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 30–39, May–June, 2007.     

虚拟层合单元在平面梁和拱形状优化中的应用

本文在虚拟层合单元的基础上,为适应形状优化时边界变化的特点,构造了曲线分层的虚拟层合梁单元,这样,求解优化问题时仅把特定的层厚作为设计变量,无须改变原始数据结构,有限元网格只在优化初始生成一次即可,具有单元数少、设计变量少的优点,算例表明此法计算简便,结果可靠,具有良好的工程应用前景。     

响应面方法在二维连续体形状优化中的应用

简单介绍了响应面方法的实质,并将其应用于二维连续体结构形状优化中,结合有限元方法,近似拟合出了状态约束(应力和位移)对设计变量的显函数,代替了复杂的解析灵敏度分析,以MSC/Patran为平台采用PCL语言开发了二维连续体形状优化的程序．文中算例通过与差分法的比较表明了这一方法的有效性和可靠性．     

Sensitivity of a general class of shape functionals to topological changes

The topological derivative represents the first term of the asymptotic expansion of a given shape functional with respect to the small parameter which measures the size of singular domain perturbations. The topological derivative has been successfully applied in the treatment of problems such as topology optimization, inverse analysis and image processing. In this paper, the calculation of the topological derivative for a general class of shape functionals is presented. In particular, we evaluate the topological derivative of a modified energy shape functional associated to the steady-state heat conduction problem, considering the nucleation of a small circular inclusion as the topological perturbation. Several methods were proposed to calculate the topological derivative. In this paper, the so-called topological-shape sensitivity method is extended to deal with a modified adjoint method, leading to an alternative approach to calculate the topological derivative based on shape sensitivity analysis together with a modified Lagrangian method. Since we are dealing with a general class of shape functionals, which are not necessarily associated to the energy, we will show that this new approach simplifies the most delicate step of the topological derivative calculation, namely, the asymptotic analysis of the adjoint state.     

Backward uncertainty propagation in shape optimization

We aim at quantifying the impact of state uncertainties in shape optimization. This provides confidence bounds for the optimal solution. The approach is presented for inverse designs where the target is assumed uncertain. No sampling of a large dimensional space is necessary, and the approach uses what is already available in a deterministic gradient‐based inversion algorithm. Our proposal is based on the introduction of directional quantile‐based extreme scenarios knowing the probability density function of the target data. We use these scenarios to define a matrix having the structure of the covariance matrix of the optimization parameters. We compare this construction to another one using the gradient of the functional by an adjoint method. The paper goes beyond inverse design and shows how to apply the method to general optimization problems. The ingredients of the paper are illustrated on a model problem with the Burgers equation and on the optimization of the shape of an aircraft. Overall, the computational complexity is comparable with the deterministic case. Copyright © 2015 John Wiley & Sons, Ltd.     

3-D Aerodynamic Shape Design Using Hessians Based on Incomplete Sensitivities

We investigate the practicability of an optimization algorithm based on the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method for 3D shape design problems, using approximate sensitivities of objective functions, where the contribution of the partial derivatives of the flow state with respect to the control variables is neglected. Therefore, it is worthwhile to investigate how an optimization method based on the Hessian behaves in this context. Indeed, the Hessian should be far from its real value if the gradient approximation is wrong. The optimization methodology is characterized by an unstructured CAD-free framework for shape and mesh deformations, an automatic differentiation of programs for the computation of the gradient of the cost function, and an unstructured flow solver. The redesign of transonic and supersonic wings has been considered and the performance of the BFGS method has been analyzed in comparison with a steepest descent method. Taking into account that a line search is too expensive to be carried out in such problems, a step size proportional to the gradient modulus has been employed for updating the control variables. Numerical results show that the BFGS method does not suffer from the approximation used in the evaluation of sensitivities, and leads to an effective improvement of the efficiency of the optimization methodology. These results can be then considered an a posteriori justification for incomplete sensitivities.     

结构形状优化与智能模型化

阻碍形状优化软件广泛应用的障碍之一是依据自然设计度量描述和建立设计模型、分析模型、优化模型以及实现三个模型之间的转换.本文称这一困难为结构形状优优设计软件的适用性.本文提出了一种方法用来动态地确定平面连续体结构形状优化过程中的边界,应用基于人工智能方法的启发式规则与技术,自动生成由设计单元法表示的几何模型,也就是将一个结构自动剖分成若干个大的四边形映射单元.这些大单元对于进一步的网络生成是必要的,同时也是向全自动计算机辅助形状优化系统前进的重要一步.     

基于解析解和新准则的深埋隧道断面形状优化

针对深埋隧道断面形状优化问题,考虑复杂受力下洞周可能存在拉应力情况,提出基于摩尔-库伦屈服准则的新的优化准则,基于解析解用模拟退火法获得洞周最小破坏范围的最优断面形状．优化模型中,以断面映射函数系数为优化的设计变量,利用混合罚函数使优化孔形满足一定的约束条件,设立孔边最大有效应力（与屈服准则相关）最小的优化准则．针对存在内压的深埋隧道断面形状优化问题,对不同内压、不同侧压系数、不同宽高比限制等情况进行优化计算,并对比两种优化准则在该情况下优化结果的异同．结果表明,对压应力主导优化的情况,新准则与原准则没有差别;对可能出现拉应力的复杂应力情况,新准则有更好的适应性,可以获得更优断面形状．     

Development and assessment of an intrusive polynomial chaos expansion-based continuous adjoint method for shape optimization under uncertainties

This article contributes to the development of methods for shape optimization under uncertainties, associated with the flow conditions, based on intrusive Polynomial Chaos Expansion (iPCE) and continuous adjoint. The iPCE to the Navier–Stokes equations for laminar flows of incompressible fluids is developed to compute statistical moments of the Quantity of Interest which are, then, compared with those obtained through the Monte Carlo method. The optimization is carried out using a continuous adjoint-enabled, gradient-based loop. Two different formulations for the continuous adjoint to the iPCE PDEs are derived, programmed, and verified. Intrusive PCE methods for the computation of the statistical moments require mathematical development, derivation of a new system of governing equations and their numerical solution. The development is presented for a chaos order of two and two uncertain variables and can be used as a guide to those willing to extend this development to a different set of uncertain variables or chaos order. The developed method and software, programmed in OpenFOAM, is applied to two optimization problems pertaining to the flow around isolated airfoils with uncertain farfield conditions.