基于变密度法的平模台振加强筋布局优化

针对混凝土预制构件水平成型振动台(简称平模台振)振动效果改进的问题,提出了通过优化振动台面加强筋布局来提高其振动效果的方法.基于变密度拓扑优化法,建立以加强筋单元相对密度为设计变量,振动台面柔度最小为目标的优化模型,并进行灵敏度分析,采用最优准则法求解优化模型;根据拓扑优化结果,对振动台面加强筋进行布局优化,获得了合理...     

基于拓扑与形状优化的小型化金属天线设计方法

天线小型化设计需要基于先进的设计方法,基于拓扑优化的设计往往存在灰度单元,因此设计结果无法直接应用,需要进一步规整设计。而对于电磁金属结构,粗糙的规整方法会引起结构性能的很大变化以致偏离最优结果。提出一种拓扑优化和形状优化相结合的方法,用于金属天线结构的小型化设计。该方法通过拓扑优化获得金属天线结构的概念构型,进而利用形状优化对概念构型进行边界规整和精细化设计。形状优化方法采用多控制点贝塞尔曲线描述拓扑概念构型,通过贝塞尔曲线控制点的移动实现天线构型的调控。给出了贝塞尔曲线控制点的设置原则,基于拓扑优化得到场量分布结果,利用较少的贝塞尔曲线控制点实现天线拓扑构型结构特征的有效调控。该方法可以获得无灰度单元残留的拓扑结果,同时可有效避免密度阈值规整方法中天线性能改变的问题,并且获得的拓扑构型边界光滑。数值算例表明拓扑优化和形状优化相结合方法的有效性。此外,该方法可拓展到其他类型电磁器件的优化设计中。     

基于应力及其灵敏度的结构拓扑渐进优化方法

在导出应力灵敏度的基础上,建立了一种改进的基于应力及其灵敏度的结构拓扑双方向渐进优化算法．该方法是对传统ESO和BESO方法的改进和完善．算例表明该方法能较大程度减少解的振荡状态数,获得更佳的拓扑结构．具有概念上的简洁性和应用上的有效性。     

热结构瞬态响应的耦合灵敏度分析方法与优化设计

研究结构瞬态热变形和热应力的灵敏度分析方法及其优化设计,灵敏度计算给出了直接法和伴随法两种算法.考虑了温度场的耦合作用,在直接法中需要计算温度场对设计变量的导数,在伴随法中需要计算热载荷对温度场的导数.数值算例验证了该方法的精度.伴随法在应用程序中的实现,为大型结构优化提供了高效率的灵敏度计算方法.     

一种新的形状灵敏度分析方法——具有两类变量的伴随方程法

本文将所考虑的问题视为具有两类独立变量的力学系统,通过建立具有两类变量的伴随系统方程,得到了定义在变化边界上的目标或约束泛函的敏度分析公式,由此建立了完全边界型的形状优化方法。     

基于应变能与频率灵敏度的静动力双目标ESO

当前在运用渐进结构优化(ESO)时,大多仅设定了单一的静力或动力目标,难以满足工程结构设计的需求。为此,将单目标优化常用的应变能灵敏度和频率灵敏度进行无量纲处理,再与多目标优化理论结合,开发出静动力双目标ESO。通过多个不同边界条件的深受弯构件数值算例,证实了新方法的运行稳定性和普遍适用性,同时还得到了静力优化与动力优化间的权重系数比取值建议。有限元对比分析结果表明,该新方法相较于传统的单目标优化,能够兼顾结构的静动力性能,使结构耗材减少但静力刚度基本维持,同时材料利用率和一阶固有频率还能不断提升。     

基于目标函数的渐进结构频率优化二次灵敏度及其应用

渐进结构频率优化中的二次灵敏度反映了结构变化后单元灵敏度的变化情况,利用二次灵敏度的二次删除法能够在一定程度上解决标准渐进结构频率优化方法中效率和精度的矛盾,是一种可行的折衷算法.本文将二次灵敏度的表达式进行了进一步推广,给出了基于目标函数的二次灵敏度的表达式,然后结合算例分析了利用基于目标函数二次灵敏度的二次删除法的性能,可以看出,推广的二次灵敏度以及利用该二次灵敏度的二次删除法能够应用于复杂条件下的结构频率优化问题.     

形状设计灵敏度分析的改进的再生核质点法

基于物质导数概念和直接微分法,将再生核质点法应用于形状设计灵敏度分析（DSA）中。导出了基于无网格近似的灵敏度方程,特别强调了在考虑形状函数关于设计变量的物质导数时无网格方法与有限元法的不同。通过对RKPM形状函数及其物质导数进行矩式显式表述,提高了无网格方法的计算效率。对两个二维线弹性问题进行了位移灵敏度和应力灵敏度分析,计算结果与解析解吻合的很好;同时通过对通常的RKPM和改进的RKPM计算耗时的比较,显示了该方法不仅有效,而且可以显著地提高计算效率。     

基于虚荷载变量的形状优化和灵敏度分析

基于选择施加在结构“控制点”上的虚荷载作为优化设计变量,针对一种新的承受约束的形状优化数值方法进行了研究。借助于节点位移与虚荷载之间的线性关系,提出了一种新的计算灵敏度系数的解析方法。利用节点移动速度域概念构造了优化新形状产生的计算公式,以结构中节点的最大应力最小化作为优化目标,通过控制网格结点的最大位移量,较好地解决了单元网格在形状优化中的扭曲问题。对三个不同的实例成功地完成了形状优化。     

Robust design in aerodynamics using third‐order sensitivity analysis based on discrete adjoint. Application to quasi‐1D flows

In this paper, the second‐order second moment approach, coupled with an adjoint‐based steepest descent algorithm, for the solution of the so‐called robust design problem in aerodynamics is proposed. Because the objective function for the robust design problem comprises first‐order and second‐order sensitivity derivatives with respect to the environmental parameters, the application of a gradient‐based method , which requires the sensitivities of this function with respect to the design variables, calls for the computation of third‐order mixed derivatives. To compute these derivatives with the minimum CPU cost, a combination of the direct differentiation and the discrete adjoint variable method is proposed. This is presented for the first time in the relevant literature and is the most efficient among other possible schemes on condition that the design variables are much more than the environmental ones; this is definitely true in most engineering design problems. The proposed approach was used for the robust design of a duct, assuming a quasi‐1D flow model; the coordinates of the Bézier control points parameterizing the duct shape are used as design variables, whereas the outlet Mach number and the Darcy–Weisbach friction coefficient are used as environmental ones. The extension to 2D and 3D flow problems, after developing the corresponding direct differentiation and adjoint variable methods and software, is straightforward. Copyright © 2011 John Wiley & Sons, Ltd.     

Sensitivity analysis for 3-D elastoplastic structural reliability

By taking the elastoplastic effect of structural material into account and based on 3-D elastoplastic stochastic finite element method, methods for sensitivity analysis with respect to both the distribution parameters of random variables and parameters in the limit state function are suggested. In the incremental iterative calculation, the sub-increment changingK, method and the corresponding formulas for accelerating convergence are used. The sensitivity of 3-D structural system reliability with respect to random variables is also studied.     

THE SATELLITE STRUCTURE TOPOLOGY OPTIMIZATION BASED ON HOMOGENIZATION METHOD AND ITS SIZE SENSITIVITY ANALYSIS

With the development of satellite structure technology, more and more design parameters will affect its structural performance. It is desirable to obtain an optimal structure design with a minimum weight, including optimal configuration and sizes. The present paper aims to describe an optimization analysis for a satellite structure, including topology optimization and size optimization. Based on the homogenization method, the topology optimization is carried out for the main supporting frame of service module under given constraints and load conditions, and then the sensitivity analysis is made of 15 structural size parameters of the whole satellite and the optimal sizes are obtained. The numerical result shows that the present optimization design method is very effective.     

Second-order sensitivity of eigenpairs in multiple parameter structures

This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenvectors are transformed into multiple parameter forms,and the second-order perturbation sensitivity matrices of eigenvalues and eigenvectors are developed.With these formulations,the efficient methods based on the second-order Taylor expansion and second-order perturbation are obtained to estimate changes of eigenvalues and eigenvectors when the design parameters are changed. The presented method avoids direct differential operation,and thus reduces difficulty for computing the second-order sensitivity matrices of eigenpairs.A numerical example is given to demonstrate application and accuracy of the proposed method.     

强迫谐振动下连续体结构拓扑优化

应用结构拓扑优化ICM(独立连续映射)方法,对强迫谐振动下结构拓扑优化问题建立了以重量极小为目标,位移幅值为约束的优化模型.位移幅值采用一阶泰勒展式近似,由于拓扑优化中设计变量数目通常很多,对强迫谐振动位移幅值的敏度分析推导了伴随法公式,使得一次敏度分析可以计算出对所有设计变量的偏导数,克服了采用直接法敏度分析中一次只能计算出对一个设计变量的偏导数的不足.算例表明用伴随法分析敏度在结构拓扑优化中可以大幅提高计算效率,ICM方法采用独立于截面及形状参数的拓扑优化设计变量更清晰地反映了拓扑优化的本质.     

基于IGA-SIMP法的连续体结构应力约束拓扑优化

建立了一种IGA-SIMP框架下的连续体结构应力约束拓扑优化方法。基于常用的SIMP模型,将非均匀有理B样条（NURBS）函数用于几何建模、结构分析和设计参数化,实现了结构分析和优化设计的集成统一。利用高阶连续的NURBS基函数,等几何分析（IGA）提高了结构应力及其灵敏度的计算精度,增加了拓扑优化结果的可信性。为处理大量局部应力约束,提出了基于稳定转换法修正的P-norm应力约束策略,以克服拓扑优化中的迭代振荡和收敛困难。通过几个典型平面应力问题的拓扑优化算例表明了本文方法的有效性和精确性。应力约束下的体积最小化设计以及体积和应力约束下的柔顺度最小化设计的算例表明,基于稳定转换法修正的约束策略可以抑制应力约束体积最小化设计中的迭代振荡现象,获得稳定收敛的优化解;比较而言,体积和应力约束下的柔顺度最小化设计的迭代过程更加稳健,适合采用精确修正的应力约束策略。     

THE LAYOUT OPTIMIZATION OF STIFFENERS FOR PLATE-SHELL STRUCTURES

The plate-shell structures with stiffeners are widely used in a broad range of engineering structures. This study presents the layout optimization of stiffeners. The minimum weight of stiffeners is taken as the objective function with the global stiffness constraint. In the layout optimization, the stiffeners should be placed at the locations with high strain energy/or stress. Conversely, elements of stiffeners with a small strain energy/or stress are considered to be used inefficiently and can be removed. Thus, to identify the element efficiency so that most inefficiently used elements of stiffeners can be removed, the element sensitivity of the strain energy of stiffeners is introduced, and a search criterion for locations of stiffeners is presented. The layout optimization approach is given for determining which elements of the stiffeners need to be kept or removed. In each iterative design, a high efficiency reanalysis approach is used to reduce the computational effort. The present approach is implemented for the layout optimization of stiffeners for a bunker loaded by the hydrostatic pressure. The numerical results show that the present approach is effective for dealing with layout optimization of stiffeners for plate-shell structures.     

A structural shape optimization system using the two-dimensional boundary contour method

Summary The research recently conducted has demonstrated that the Boundary Contour Method (BCM) is very competitive with the Boundary Element Method (BEM) in linear elasticity Design Sensitivity Analysis (DSA). Design Sensitivity Coefficients (DSCs), required by numerical optimization methods, can be efficiently and accurately obtained by two different approaches using the two-dimensional (2-D) BCM as presented in Refs. [1] and [2]. These approaches originate from the Boundary Integral Equation (BIE). As discussed in [2], the DSCs given by both BIE-based DSA approaches are identical, and thus the users can choose either of them in their applications. In order to show the advantages of this class of DSA in structural shape optimization, an efficient system is developed in which the BCM as well as a BIE-based DSA approach are coupled with a mathematical programming algorithm to solve optimal shape design problems. Numerical examples are presented. Received 20 July 1998; accepted for publication 7 December 1998     

Two forms of continuum shape sensitivity method for fluid–structure interaction problems

Continuum Sensitivity Equation (CSE) methods for deriving and computing derivatives with respect to shape design variables are developed in two forms and compared in their application to fluid–structure interaction (FSI) problems. The local derivative form poses the CSEs in terms of the partial derivatives of the state variables with respect to shape parameters, while the CSEs in total derivative form are posed in terms of the total derivative, also known as the material or substantial derivative. In the literature CSEs are often posed in local form for fluids and total form for solids. The two forms are compared here for the purpose of applying a single form to both fluid and structure domains. The local form, also known as the boundary velocity method, requires design velocity only at the boundaries and interfaces of the domains to pose the CSEs. In contrast, the total form, also known as the domain velocity method, requires the design velocity in the whole domain. The local form requires higher-order spatial derivatives of the analysis solution than the total form, which affects the accuracy of its results. Higher order p-elements are shown to be a remedy to the inaccuracy of local form CSE seen in the literature for finite element solutions. The practicality, accuracy, and efficiency of these two CSE forms are compared based on the implementation and computed derivatives for three examples: a linear Timoshenko beam subject to a tip force, fluid flow around an airfoil, and an airfoil attached to a nonlinear joined beam subject to a gust load.