基于NURBS插值的三维渐进结构优化方法

提出一种基于等几何控制点密度变量的三维双向渐进结构拓扑优化方法。在当前列式下,高阶NURBS基函数被同时用于CAD模型中NURBS实体片的几何场、位移场和温度场以及密度场插值,实现了几何模型、分析模型和优化模型的有效统一,确保了位移场、温度场及密度场的高阶连续性;详细推导了基于等几何控制点密度变量的三维渐进结构法模型及其灵敏度分析列式;最后几个典型的数值算例,包括最小柔顺性、热传导优化问题及三维结构自由振动的基频最大化问题,验证了本文方法的有效性。     

基于无单元Galerkin方法的受迫振动下的连续体结构拓扑优化

基于无单元Galerkin法(EFG)对受迫振动下的连续体结构进行拓扑优化设计。选取节点的相对密度作为设计变量,以动柔度最小化为目标函数,基于带惩罚的各向同性固体材料模型（SIMP）建立了受迫振动下的连续体结构拓扑优化模型,采用伴随法求解得到目标函数的敏度分析公式,利用优化准则法对优化模型进行求解。通过经典的二维连续体结构拓扑优化算例证明了该方法的可行性和有效性。     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

基于无网格数值求解技术的二维连续体结构拓扑优化设计

将无网格径向点插值法(RPIM)引入到连续体结构拓扑优化设计中。在优化问题中,选取节点的相对密度作为设计变量,以结构的柔度最小化作为目标函数,基于带惩罚的各向同性固体材料模型(SIMP)建立了结构拓扑优化的数学模型,推导了目标函数和体积约束的灵敏度,利用优化准则法进行求解。算例表明了应用无网格径向点插值法进行结构拓扑优化设计的可行性和有效性,同时表明选取节点的相对密度作为设计变量可以有效地克服拓扑优化中的棋盘格现象。     

基于扩展多尺度有限元法的含液闭孔材料拓扑优化

基于扩展多尺度有限元方法,提出了含液闭孔结构多尺度拓扑优化方法。该多尺度优化方法旨在研究含液闭孔胞元布局对整体含液闭孔结构力学性能的影响。首先针对含液闭孔结构的整体结构柔顺性问题,采用类似SIMP模型对结构的宏观粗网格等效刚度阵进行插值,建立含液闭孔结构柔顺性的拓扑优化列式;其次,针对含液闭孔材料能够利用胞体内部液体腔体积增量产生变形的特性,提出含液闭孔材料柔性机构的概念,并以结构指定位置方向输出位移为目标,建立液体体积膨胀作用下的含液闭孔柔性机构多尺度拓扑优化数学模型。本文基于自主软件平台SiPESC完成了程序研发,并通过数值算例验证了所提出的拓扑优化方法的有效性。     

基于无网格Galerkin方法的整体式柔性机构的多准则拓扑优化设计

基于一种新型的数值计算技术—无网格伽辽金法,提出了一种整体式柔性机构拓扑优化设计的新方法.利用移动最小二乘形函数和伽辽金的弱变分形式建立弹性问题的控制方程,用Lagrange乘子法增强本质边界条件.在优化问题中,同时综合考虑机构的柔性和结构的刚度要求,用折衷规划法建立了柔性机构拓扑优化的多准则优化模型,这样对于非凸的优化问题也能保证搜索到Pareto解集所有的有效解.基于SIMP密度函数惩罚模型和优化准则法,建立了一种设计变量的显式迭代格式.运用经典算例证明了文中方法的正确性和有效性.     

Sensitivity analysis for fixed-grid shape optimization by using oblique boundary curve approximation

The remesh-free property is the most attractive feature of the various versions of fixed-grid-based shape optimization methods. When the design boundary curves do not pass through the predetermined analysis grids, however, the element stiffness as well as the stress along the curves may be computed inaccurately. Even with the popular area-fraction-based stiffness evaluation approach, the whole optimization process may become quite inefficient in such a case. As an efficient alternative approach, we considered a stiffness matrix evaluation method based on the boundary curve approximation by piecewise oblique curves which can cross several elements. The main contribution of this work is the analytic derivation of the shape sensitivity for the discretized system by the fixed-grid method. Since the force term in the sensitivity equation is associated only with the elements crossed by the design boundary curve, we only need the design velocities of the intersecting points between the curve and the fixed mesh. The present results obtained for two-dimensional elasticity and Poisson’s problems are valid for both the single-scale standard fixed-grid method and the multiscale fictitious domain-based interpolation wavelet-Galerkin method.     

考虑结构自重的基于NURBS插值的3D拓扑描述函数法

在许多如大坝、桥梁等大型土木工程结构中,结构的自重是初始设计阶段必须考虑的重要载荷之一,因此研究自重载荷作用下的结构拓扑优化设计问题具有十分重要的意义.针对考虑自重载荷作用的拓扑优化问题所面临的主要困难,总结了现有处理考虑自重载荷的拓扑优化问题的三类主要方法;提出一种基于非均匀有理B样条(non-uniform rational B-splines,NURBS)基函数插值的拓扑描述函数方法,基于此方法研究了考虑设计依赖自重载荷作用的2D/3D结构优化设计问题.在列式下,高阶NURBS基函数被同时用于三维NURBS实体片中的几何场、位移场及设计变量场插值,实现了几何模型、分析模型和优化模型的有效统一,确保了位移场及设计变量场的高阶连续性;详细推导了基于NURBS基函数插值的考虑自重载荷作用的三维结构拓扑优化模型及其灵敏度列式,并采用移动渐进线方法(method of moving asymptotes,MMA)进行了优化求解;多个算例验证了方法的有效性和稳定性,结果表明,优化迭代过程稳健,收敛快,能够有效地克服自重载荷作用下连续体结构拓扑优化中经常遇到的低密度区域材料的寄生效应及目标函数的非单调性等问题.     

Evolutionary topology optimization for acoustic-structure interaction problems using a mixed u/p formulation

In this work, we present a topology optimization method for acoustic-structure interaction problems, which combines bi-directional evolutionary structural optimization (BESO) with a mixed displacement-pressure (u/p) formulation as an effective and straightforward design method for a multi-physics system involving acoustic-structure interactions. Due to the binary characteristics of the BESO and the multi-physics modeling approach of the mixed formulation, the proposed optimization procedure could benefit from high computational efficiency and high-quality design in acoustic-structure interaction problems. Several topology optimization problems for vibro-acoustic systems are carried out, in order to demonstrate the effectiveness of the presented method.     

连续体结构的拓扑优化设计

对基于有限元数值求解技术的连续体结构的拓扑优化设计技术进行了综述. 利用密度-刚度插值格式和优化准则方法, 以结构的柔度最小化作为优化的目标函数, 论述并建立线弹性结构的静力学拓扑优化设计的数学模型和设计变量显示的迭代格式; 基于数学规划方法中的一种凸规划方法----移动渐近线方法和密度方法, 以结构的频率最大化作为优化的目标函数, 论述并建立了特征值问题拓扑优化设计的数学模型和设计变量隐式的更新方法. 对多目标拓扑优化问题、柔性机构的拓扑优化问题以及多物理场拓扑优化设计问题进行了讨论. 对优化结构中出现的棋盘格式和网格依赖性等数值计算问题进行剖析和讨论, 介绍和分析了目前解决数值计算问题常见的方法, 在此基础上对边界扩散现象进行了讨论. 给出了连续体结构拓扑优化设计的程序流程, 并用Matlab程序实现了算法, 通过几个典型的算例证明所综述方法的有效性.      

积分背景网格对自适应EFG法拓扑优化的影响研究

采用应力能量范数作为误差指标,探讨了EFG法中积分背景网格对计算精度的影响,得到了合理划分背景网格的建议;建立了以节点密度为设计变量、以最小化柔度为优化日标的拓扑优化模型。采用以节点密度值为加点判据的自适应规则加点方案,开展了连续体结构的拓扑优化研究,该加点方案能有效地减少设计变量的个数,探讨了背景网格对拓扑优化结果的影响。算例结果表明,采用合适的背景网格不仅能进一步减少设计变量的个数,而且能够改善拓扑优化结果的光滑性,使计算效率和精度得到提高。     

Optimum shape and topology design using the boundary element method

A procedure is developed for simultaneous shape and topology design optimization of linear elastic two-dimensional continuum structures. An intuitive approach is presented to treat such topological problems whereby material is eliminated from within the structure (by introducing holes at regions of low stress) through a sequence of shape optimization processes. A mathematical programming technique coupled with the boundary element (BE) method of response and sensitivity analyses enables the optimal positioning of these holes plus optimization of the overall structural shape. The analytical derivative BE formulation is explained together with the use of appropriate design velocity fields, and example problems are solved to demonstrate the optimization procedure.     

传热结构的多目标拓扑优化设计研究

深入分析了传热结构多目标拓扑优化设计中的几个关键问题。提出了基于结构柔度最小化和结构散热弱度最小化的多目标拓扑优化设计方法,建立了传热结构的多目标拓扑优化设计模型,推导了传热结构多目标拓扑优化中用于迭代分析求解的优化准则算法和敏度分析方程。通过数值计算验证了理论和算法的有效性。     

基于节点密度的自适应EFG法拓扑优化方法研究

针对无网格法拓扑优化中设计变量庞大、计算效率低的问题,提出了基于节点密度的自适应EFG拓扑优化方法。该方法以节点相对密度为设计变量,建立了以结构柔度最小为优化目标的自适应EFG拓扑优化模型,给出了拓扑优化中自适应节点加密准则、节点加密方法和新增节点设计变量更新规则。算例结果表明,采用本文提出的自适应EFG拓扑优化方法进行结构拓扑优化时,可以减少结构分析和优化的设计变量数目,提高优化效率。     

Finite Element Error Analysis Approach for Three-Dimensional Incompressible Viscous Fluid Flow Analysis

In our previous research, the modified Galerkin method was proposed as one of the most efficient methods for the analyses of convection-diffusion problems and two-dimensional viscous fluid flow problems. In this modified Galerkin method, the inertia term is considered explicitly, so only the symmetrical matrixes appear. Then an artificial viscosity is introduced through an error analysis approach to improve its accuracy and stability. In this paper, we proposed a new finite element formulation for three-dimensional incompressible viscous fluid flow analysis. This formulation (‘MS’ algorithm and ‘MSR’ algorithm) is based on the modified Galerkin method coupled with the Semi-Implicit Method for Pressure-Linked Equations. The cubic cavity flow problems were investigated for the Reynolds number of 400, 1,000, 2,000 and 3,200 using non-uniform meshes. Finally, we confirmed the effectiveness of our proposed method through the comparison with other research works.     

弹性力学问题的局部Petrov—Galerkin方法

提出了弹性力学平面问题的局部Petrov-Galerkin方法,这是一种真正的无网格方法。这种方法采和移动最小二乘近似函数作为试函数,并且采用移动最小二乘近似函数的权函数作为加权残值法加权函数;同时这种方法只包含中心在所考虑点处的规则局部区域上以及局部边界上的积分,所得系统矩阵是一个带状稀疏矩阵,该方法可以容易推广到求解非线性问题以及非均匀介质的力学问题。还计算了两个弹性力学平面问题的例子,给出了位移和能量的索波列夫模及其相对误差。所得计算结果证明：该方法是一种具有收敛快、精度高、简便有效的通用方法;在工程中具有广阔的应用前景。     

考虑微结构连接性的双尺度结构自然频率拓扑优化

多孔材料因具有轻量化、高孔隙率和减振/散热等优良多物理特性,在航空航天等领域具有广阔应用前景。采用拓扑优化方法对含多种多孔材料的结构进行结构与材料微结构构型一体化设计,有助于获得具有优良力学性能的结构设计。然而,传统逆均匀化微结构设计方法无法确保不同多孔材料微结构之间的连接性,设计结果不具备可制造性。本文面向含多种多孔材料的双尺度结构基频最大化设计问题,考虑不同微结构之间的连接性,协同设计多孔材料的微结构构型及其在宏观尺度下的布局。采用均匀化方法计算多孔材料的宏观等效力学性能,通过对不同多孔材料微结构单胞的边界区域采用相同的拓扑描述确保双尺度优化过程中任意空间排布下不同微结构的连接性,并通过优化算法确定微结构间的连接形式及微结构拓扑。在宏观尺度,提出结合离散材料插值模型和RAMP插值模型RAMP (Rational Approximation of Material Properties)的多孔材料各向异性宏观等效刚度及质量插值模型,获得清晰的多孔材料宏观尺度布局并减轻优化过程中伪振动模态的影响。建立以双尺度结构基频最大化为目标,以材料用量为约束的优化列式,推导灵敏度表达式,并基于梯度优化算法求解双尺度结构拓扑优化问题。数值算例表明,采用本文优化方法能够有效确保基频最大化双尺度结构设计中不同多孔材料微结构之间的连接性,增强优化设计结果的可制造性。     

Modified Galerkin Method for Unsteady Two-Dimensional Viscous Fluid Flow Analysis

The Modified Galerkin Method (MGM) has been proposed as one of the most efficient methods for two-dimensional convection-diffusion equations. In the MGM, the non-symmetric matrices, which are derived from the convection term in the Galerkin formulation, are not used, and an artificial diffusion is introduced through an error analysis approach to improve its discretization accuracy in both time and space directions. In this study, the MGM is applied for two-dimensional viscous fluid flow analysis, and the driven cavity flow problems are solved up to Reynolds number of 10,000 using the vorticity-stream function formulation and non-uniform meshes. The results show the effectiveness of MGM.