基于虚拟材料的T-B梁拓扑优化ESO方法

自从Zhou和Rozvany于2001年以T-B梁(Tie-Beam)问题质疑渐进结构优化方法ESO(Evolutionary Structural Optimization)的算法收敛性以来,结构最优化领域的研究者提出了多种算法以求解决这一问题。T-B梁问题仍是当前结构拓扑优化领域的研究热点,因为现有的各种T-B梁解法或是存有不足或是无法获得满意解答。本文在传统ESO算法中,将虚拟材料引入待删除单元,用以检测结构传力路径是否受到破坏,进而确定单元删除的合理性。算例表明,本文算法可有效防止ESO方法求解T-B梁问题时的失效并获得最优解,且本文算法只需在ESO迭代中附加一次检测,不改变ESO方法的迭代进程和寻优能力。     

薄元分解与Laplacian光顺相结合的四面体有限元网格优化方法

提出一种有效的三维实体四面体有限元网格质量优化方法以满足有限元分析对网格质量的要求。对薄元分解法进行改进,改进的薄元分解法更全面地考虑了各种劣质单元类型,能够对三维实体网格剖分中产生的各种类别的孤立劣质单元进行有效的分解;将改进的薄元分解法与Laplacian光顺优化方法相结合以解决某些网格剖分算法如推进波前法和Delaunay三角化方法产生的非孤立劣质单元问题。经过实例检验,本文提出的四面体单元网格优化算法健壮有效、易于实现,能够显著提高最差单元的质量。     

基于OPS和拓扑交换的离心泵网格质量提高算法

针对OPS算法本身固有的无法提高"全边界单元"质量的缺陷,提出了一种结合OPS和拓扑优化的网格质量提高算法。该方法通过判别"全边界单元"中边(或面)受到相邻单元共享的形式,采用不同的拓扑交换策略,且在拓扑交换过程中仅考虑网格的拓扑相关性,并不考虑网格的质量,保证将"全边界单元"转化为OPS算法提高其质量的其他类型边界单元。算例分析表明,本文算法相比Freitag的组合算法在消除极端二面角方面更具优势,相比OPS算法对提高离心泵边界网格及整体网格质量也具有较好的优势。     

一种变位移约束限的结构拓扑优化方法

针对仅有位移约束和重量最小的结构拓扑优化问题, 基于ICM(独立、连续、映射)方法思路,提出了一种变位移约束限的结构拓扑优化方法. 在每一轮子循环迭代求解开始时,为了控制拓扑设计变量的变化量, 形成和引进了新的位移约束限. 另外,建立了单元删除阈值和几轮迭代循环的单元删除策略. 为了确保优化迭代中结构非奇异和方法具有增添单元的功能,在结构孔洞和边界周围引入了一层人工材料单元,并建立了一套有效结构信息到结构最大设计域信息的映射转换方法. 结合对偶求解方法,形成了一种新的连续体结构的拓扑优化方法. 给出的算例表明该方法没有目标函数的振荡现象,且验证了该方法的正确性和有效性. 关键词:拓扑优化; 位移约束; 连续体结构; ICM方法; 渐进结构优化;      

一种基于单元几何变形的网格修匀方法

基于单元几何变形操作提出一种高效的非结构网格质量修匀方法。其基本过程是先对每个单元独立地进行拉伸-收缩操作以优化单元的形状,然后在整个网格中通过对各单元的节点位置进行加权平均来获得改善后的网格。为进一步提高修匀方法对网格质量的优化效果,并使得该方法具备一定的网格调整能力,结合动网格技术提出了对单元进行大范围和较大幅度移动的策略;在修匀过程中还通过适当算法调整单元形心位置和单元尺寸,进一步增强了修匀方法对网格局部进行疏密调节的能力。本文方法可适用于平面和三维非结构网格的质量改善及网格调整。若干算例表明了方法的有效性。     

基于SIMP插值模型的渐进结构优化方法

在传统渐进结构优化算法(ESO)及带惩罚的变密度法(SIMP)的基础上,本文提出将二者相结合的基于SIMP插值模型的渐进结构优化算法.该方法通过缩小传统ESO算法中的进化步长,从而缩小了由于进化步长过大而导致的敏度评估误差,使得ESO算法在合理性及通用性上获得了较大改善.数值算例表明,该方法在保持了常规ESO方法的优点的同时,拥有更高的稳定性和可靠性.     

用反向渐进结构优化方法研究洞室支护优化

洞室开挖支护的优化是在铁路隧道、矿山开挖、水电站建设中经常面对的课题,特别是研究洞室开挖支护拓扑优化有着重要的意义。本文利用近年来发展的渐进优化方法(ESO)的概念,提出了反向ESO方法的迭代格式,建立了洞室开挖支护优化的准则及其敏感度,研究了均匀地基不同地应力条件使得开挖卸荷产生的总应变能最小的最优支护拓扑。结果证明了反向ESO方法的适用性,最优支护拓扑大致成椭圆或者圆形,其主轴方向与最大主应力主轴方向相同。     

基于位移和应力灵敏度的结构拓扑优化设计

针对仅有应力约束以及含位移和应力约束的重量最小结构拓扑优化问题,基于ICM(独立、连续、映射)方法和渐进结构优化方法的思路,提出了一种结构拓扑优化方法. 在优化迭代循环的每一轮子循环迭代求解开始时,通过形成和引进新的位移和应力约束限, 自动构建设计变量移动限.将结构应力约束归并为几个最可能的有效应力约束,大大减少了应力灵敏度的分析量. 另外,建立了单元删除阈值和几轮迭代循环的单元删除策略. 为了确保优化迭代中结构非奇异和方法具有增添单元的功能,在结构孔洞和边界周围引入了一层人工材料单元,构建了其等效的优化模型. 结合位移和应力灵敏度分析,形成了一种新的连续体结构的拓扑优化方法. 给出的算例验证了该方法的正确性和有效性.      

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

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

基于一种动态删除率的ESO方法

渐进结构优化方法 ESO(Evolutionary Structural Optimization)的基本思想是基于单元灵敏度,通过把无效或低效的单元逐步从结构中删除,从而得到优化的拓扑结构。经过20年的发展,其算法结构、理论研究以及实际工程应用领域已经取得了大量成果。在原始ESO算法中删除单元的数目是由固定删除率RR的取值决定的,设计者无法预期每一个迭代步删除单元的数量。对于一个初始满设计区域,合理的删除策略为随着迭代的进行,随着结构应力分布越均布,应该逐渐减少单元删除数量。基于此,本文构造了一种动态删除率,使得随着迭代的进行单元的删除数量逐渐减少,相对于前人构造的动态删除率更为简单明了,人为控制参数更少,并且通过算例证明该方法相比原始删除策略具有更好的优化效果。     

Collocated discrete least‐squares (CDLS) meshless method: Error estimate and adaptive refinement

Meshless methods are new approaches for solving partial differential equations. The main characteristic of all these methods is that they do not require the traditional mesh to construct a numerical formulation. They require node generation instead of mesh generation. In other words, there is no pre‐specified connectivity or relationships among the nodes. This characteristic make these methods powerful. For example, an adaptive process which requires high computational effort in mesh‐dependent methods can be very economically solved with meshless methods. In this paper, a posteriori error estimate and adaptive refinement strategy is developed in conjunction with the collocated discrete least‐squares (CDLS) meshless method. For this, an error estimate is first developed for a CDLS meshless method. The proposed error estimator is shown to be naturally related to the least‐squares functional, providing a suitable posterior measure of the error in the solution. A mesh moving strategy is then used to displace the nodal points such that the errors are evenly distributed in the solution domain. Efficiency and effectiveness of the proposed error estimator and adaptive refinement process are tested against two hyperbolic benchmark problems, one with shocked and the other with low gradient smooth solutions. These experiments show that the proposed adaptive process is capable of producing stable and accurate results for the difficult problems considered. Copyright © 2007 John Wiley & Sons, Ltd.     

惯性稳定平台扩张状态观测器/PD复合控制

为提高惯性稳定平台控制系统的稳定精度,在常规PID控制的基础上提出了一种扩张状态观测器与PID相结合的复合控制算法。利用扩张状态观测器将惯性稳定平台的各种内部扰动和外部扰动都视为总和扰动并观测出来,然后通过PID控制器进行误差反馈控制,从而提高控制系统的扰动抑制能力与稳定精度。以Lu Gre摩擦模型加入控制模型进行仿真分析,并通过北航自研的惯性稳定平台进行实验验证。结果表明:扩张状态观测器/PD复合控制方法具有高的扰动抑制能力,可显著提高稳定平台稳定精度。相比常规PID方法,扩张状态观测器/PD复合控制使横滚框和俯仰框的稳定精度分别提高了33.23%和55.01%。     

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

Versatile anisotropic mesh adaptation methodology applied to pure quantity of interest error estimator. Steady,laminar incompressible flow

We introduce a new flexible mesh adaptation approach to efficiently compute a quantity of interest by the finite element method. Efficiently, we mean that the method provides an evaluation of that quantity up to a predetermined accuracy at a lower computational cost than other classical methods. The central pillar of the method is our scalar error estimator based on sensitivities of the quantity of interest to the residuals. These sensitivities result from the computation of a continuous adjoint problem. The mesh adaptation strategy can drive anisotropic mesh adaptation from a general scalar error contribution of each element. The full potential of our error estimator is then reached. The proposed method is validated by evaluating the lift, the drag, and the hydraulic losses on a 2D benchmark case: the flow around a cylinder at a Reynolds number of 20.     

An efficient method for computing local quantities of interest in elasticity based on finite element error estimation

We present an efficient finite element method for computing the engineering quantities of interest that are linear functionals of displacement in elasticity based on a posteriori error estimate. The accuracy of quantities is greatly improved by adding the approximate cross inner product of errors in the primal and dual problems, which is calculated with an inexpensive gradient recovery type error estimate, to the quantities obtained from the finite element solution. With less CPU time, the accuracy of the improved quantities obtained with the proposed method on the coarse finite element mesh is similar to that of the quantities obtained from the finite element solutions on the finer mesh. Three quantities related to the local displacement, local stress and stress intensity factor are computed with the proposed method to verify its efficiency.     

Goal‐oriented space–time adaptivity in the finite element Galerkin method for the computation of nonstationary incompressible flow

This paper presents a general strategy for designing adaptive space–time finite element discretizations of the nonstationary Navier–Stokes equations. The underlying framework is that of the dual weighted residual method for goal‐oriented a posteriori error estimation and automatic mesh adaptation. In this approach, the error in the approximation of certain quantities of physical interest, such as the drag coefficient, is estimated in terms of local residuals of the computed solution multiplied by sensitivity factors, which are obtained by numerically solving an associated dual problem. In the resulting local error indicators, the effects of spatial and temporal discretization are separated, which allows for the simultaneous adjustment of time step and spatial mesh size. The efficiency of the proposed method for the construction of economical meshes and the quantitative assessment of the error is illustrated by several test examples. Copyright © 2012 John Wiley & Sons, Ltd.     

Three-dimensional adaptive finite element computations and applications to non-Newtonian fluids

An adaptive strategy for the finite element solution of three-dimensional viscous flow problems is defined and implemented. The solution strategy is based on an advancing front mesh generator making use of binary data structures for fast geometrical data handling. The error is estimated a posteriori with a residual-type bound. The error estimate is shown to exhibit proper convergence for tetrahedral elements. Its combination with the mesh generator and an interpolation scheme for unstructured meshes is shown to generate adaptive meshes and to reduce the solution cost for a given error level, as illustrated by the isothermal flow of a shear-thinning fluid.     

Accurate and robust adaptive mesh refinement for aerodynamic simulation with multi‐block structured curvilinear mesh

A multi‐block curvilinear mesh‐based adaptive mesh refinement (AMR) method is developed to satisfy the competing objectives of improving accuracy and reducing cost. Body‐fitted curvilinear mesh‐based AMR is used to capture flow details of various length scales. A series of efforts are made to guarantee the accuracy and robustness of the AMR system. A physics‐based refinement function is proposed, which is proved to be able to detect both shock wave and vortical flow. The curvilinear mesh is refined with cubic interpolation, which guarantees the aspect ratio and smoothness. Furthermore, to enable its application in complex configurations, a sub‐block‐based refinement strategy is developed to avoid generating invalid mesh, which is the consequence of non‐smooth mesh lines or singular geometry features. A newfound problem of smaller wall distance, which negatively affects the stability and is never reported in the literature, is also discussed in detail, and an improved strategy is proposed. Together with the high‐accuracy numerical scheme, a multi‐block curvilinear mesh‐based AMR system is developed. With a series of test cases, the current method is verified to be accurate and robust and be able to automatically capture the flow details at great cost saving compared with the global refinement. Copyright © 2015 John Wiley & Sons, Ltd.     

黏弹性人工边界等效荷载计算的改进方法

黏弹性人工边界在场地地震反应和结构-地基动力相互作用等问题的计算中已得到了广泛的应用.地震波在黏弹性人工边界中的输入是通过将地震波转化为作用于人工边界处的等效载荷来实现的.计算等效节点载荷的常规方法默认边界节点对应区域内的应力为均布力,但实际上该节点对应区域内的应力分布通常是不均匀的.本文在有限元方法结合黏弹性局部人工边界的显式时域波动方法的基础上,建立了无限域散射问题地震波等效载荷计算的一种改进方法.该方法采用细化网格与应力积分相结合的方法计算人工边界等效节点力,有效地降低了人工边界上等效节点力的计算误差.以不同角度入射地震波的二维算例为例,算例给出的波场位移云图和节点位移时程曲线验证了本文方法的有效性,其计算精度与网格尺寸和地震波入射角度密切相关,且网格越小、入射角度越小,计算精度越高.对于相同的网格尺寸,本文采用方法的计算精度明显高于常规方法,尤其是对于斜入射问题优势更为明显.