用拓扑优化设计多体分比蜂窝状结构（英文）

提出用拓扑优化及增材制造法设计同一材料制作的蜂窝状结构,并把设计可归结为在不同的设计域内使用不同的体分比的拓扑优化问题。首先推导了多体分比约束条件下的蜂窝状结构拓扑优化问题;然后提出了基于多阈值的等值面法及与商用有限元分析软件对接的算法。给出了有三个体积约束的最小柔顺性问题的数值算例。计算结果表明,对于具有三个体积约束的拓扑优化设计问题,存在一个最佳的体积约束组合;使用等量的材料,与两个体积约束的最佳拓扑设计相比,用最佳体积约束组合所得到的优化结构的柔顺性还要小10%。     

考虑回转性能的三维结构拓扑优化设计

以考虑运动学性能的回转运动体材料最优分布问题为背景,提出并求解了转动惯量约束下的三维结构拓扑优化问题.常规的结构拓扑优化通常构造为给定材料体积约束下的结构最小柔顺性问题.本文指出,引入结构转动惯量约束对于回转结构的拓扑优化具有重要意义.本文将重量约束下的优化准则法计算格式拓广到单个转动惯量约束下的最小柔顺性问题,而对于包含转动惯量约束的多约束问题采用数学规划法求解.本文给出的数值算例验证了所提出算法的正确性,并且表明引入转动惯量约束后可得到与常规拓扑优化问题不同的结构布局形式.     

多工况下结构鲁棒性拓扑优化设计

针对工程中存在多个随机不确定工况载荷作用的情况, 将鲁棒性设计思想引入到连续体结构拓扑优化设计, 发展和完善不确定性优化理论和计算方法. 基于概率模型和SIMP方法,提出以结构柔顺度标准差最小化为目标、具有体积约束的连续体鲁棒性拓扑优化数学模型.通过对目标函数及其灵敏度计算公式的推导, 采用数学规划法实现优化问题的求解. 数值算例验证了所提优化模型的正确性及算法的有效性, 并通过与确定性优化结果的比较,证明鲁棒性拓扑优化能够给出结构柔顺性变异更小的材料分布.      

具有体积约束的连续体结构拓扑优化仿生方法

基于Wolff法则的连续体拓扑优化仿生方法是将待优化的结构被看作是一块遵从Wolff法则生长的骨骼,把寻找结构最优拓扑的过程比拟为骨骼的生长过程.为了解决具有非应变约束的连续体拓扑优化问题,引入浮动参考应变区间.该文以具有体积约束的最小柔顺性连续体为例,分析了浮动参考应变区间法的可行性,即在生长过程中,设计域内各点处的材料更新由当前参考应变区间确定,而参考应变区间的更新由当前结构体积率与约束体积率确定.数值结果表明该方法可行.     

基于多相材料的稳态热传导结构轻量化设计

在多相材料的结构拓扑优化问题中,通常给定各相材料体积约束或材料总重量约束作为材料的控制用量.在结构轻量化设计的实际工程背景下,以结构总重量最小化为目标的优化模型具有明确的工程意义.针对含多相材料的稳态传热结构拓扑优化问题,提出了以结构总重量最小化为目标和给定热柔顺度为约束的多工况连续体结构拓扑优化建模方法.遵循独立连续映射建模方式,采用两类独立拓扑变量分别表征单元热传导矩阵和单元重量状态.推导了热柔顺度和总重量对设计变量的敏度,基于一阶和二阶泰勒展开得到各自的近似表达式.通过求解偏微分方程,实现了约束函数一次项过滤,消除了棋盘格现象和网格依赖性问题,并保证了约束方程在过滤后严格成立.建立的近似优化模型具有二次函数形式的目标函数和一次函数形式的约束函数.基于对偶序列二次规划方法对优化模型进行求解直至收敛.通过四个三维结构数值算例分析对比了热柔顺度约束限值、不同材料混合及多工况、多约束条件对优化结果的影响.数值算例结果表明,本文提出的优化方法在基于多相材料的多工况稳态热传导结构轻量化设计中具有可行性和有效性.     

连接组合结构协同动力学拓扑优化设计

工程实际结构通常是由多个部件组合而成, 且各部件通过连接构件传递彼此间的载荷和振动能量. 连接构件的布局设计与约束状况对整个结构的拓扑构型、动态性能以及承载能力等均有较大的影响. 本文研究连接组合结构构型与部件间连接构件布局的协同动力学优化问题, 使整体结构在简谐激励作用下动柔顺度达到最小. 以弹簧连接单元模拟连接构件的约束及承载状况, 将承力构件材料的相对密度与弹簧连接单元的相对刚度同时作为设计变量. 在材料体积约束以及连接构件数量约束的条件下, 采用基于梯度的优化算法开展组合结构的拓扑构型与连接构件的布局协同优化设计. 通过与无连接约束构件的单体式结构拓扑优化结果进行对比, 展示了组合结构拓扑构型的变化, 以及连接约束的布局设计对整体材料分布和结构动力性能的影响. 数值结果表明, 虽然组合结构协同优化设计的动柔顺度总是大于单体式结构的结果, 但结构固有频率的变化却具有一定的偶然性, 即可提供更加优越的结构构型与连接布局设计.      

狭长结构拓扑优化

通常的拓扑优化是在给定区域内，通过设计材料分布实现结构拓扑形式优化。对于设计区域的长和宽相近的平面问题，现行的方法可得到清晰的拓扑。但是，狭长结构的设计域具有大的长宽比。为了保证基结构包含足够多的拓扑形式，宽度方向要求有一定量的有限元分割，从而导致整体网格数和设计变量多、问题求解困难。本文提出了通过基本结构拼装的狭长结构拓扑优化方法，建立了以最小平均柔顺性密度为目标、同时设计材料分布和设计域几何尺度的基本结构的拓扑优化问题的数学提法和求解方法。利用所提出的问题提法和求解方法，设计了狭长悬臂梁的拓扑形式，讨论了危险截面的弯矩与剪力的相对值以及材料体积约束对拓扑形式的影响。数值结果表明，不同的弯矩与剪力的相对数值对应不同的拓扑形式，随着相对数值的增加，梁的拓扑形式由类桁架结构逐渐变成竖直立板加强的框架式结构。     

基于可行域调整策略的多约束材料优化设计

研究基于可行域调整策略的多约束结构材料优化问题。建立了以结构柔顺度为目标函数,考虑体积、位移、结构自振频率为约束条件的结构材料优化模型和目标函数及约束函数近似式,并采用可行域调整策略和对偶算法求解。典型结构材料优化设计结果表明:可行域调整策略使优化过程平稳,获得的序列细观结构拓扑构型清晰;相比于单一体积约束结构拓扑优化求解,考虑多约束结构拓扑优化能获得更好的黑白占优解。     

面向应力约束的独立连续映射方法

大多数已有的拓扑优化研究为系统刚度最大化设计,尤其以体积比约束下的静态柔顺度最小化问题为典型.从工程角度出发,结构强度设计至关重要.以往的应力研究表明,应力约束拓扑优化存在着奇异性、约束数目庞大、高度非线性特性等诸多数值困难.为了实现应力约束下的拓扑优化设计,采用归一化p范数应力指标以减少单元应力约束数目.遵循独立连续映射建模方式,引入密度变量的倒变量函数作为设计变量.推导了应力约束函数和体积目标函数对设计变量的敏度,并基于一阶和二阶泰勒近似得到各自的显式表达式.通过构造的系列二次规划子问题,原拓扑优化问题采用序列二次规划算法高效求解.二维数值算例考察了结构刚度和强度设计结果的异同,以及不同应力约束上限值对应力约束拓扑优化结果的影响.通过提出方法与传统变密度法结果的比较,说明提出的独立连续映射方法在应力约束下具有可行性和有效性.优化结果也表明了考虑应力约束的连续体拓扑优化具有必要性.      

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

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

基于改进的双向渐进结构优化法的应力约束拓扑优化

工程结构设计时经常需要限制最大名义应力,以避免发生断裂或疲劳破坏,一个有效的策略是采用拓扑优化方法. 常规的双向渐进结构优化法(bi-evolutionary structural optimization, BESO)不能有效求解应力约束拓扑优化问题,为此本文提出一种改进的双向渐进结构优化方法,处理体积和应力约束下的最小柔顺性问题. 引入基于K-S函数的全局应力度量,以减小大量局部应力约束引起的计算代价. 采用拉格朗日乘子法将应力约束函数引入到目标函数,然后由二分法确定合适的拉格朗日乘子的值使得应力约束得到满足. 而且,详细推导了基于BESO方法的应力约束拓扑优化模型及其灵敏度列式,最后通过三个典型拓扑优化算例验证改进方法的有效性. 为展示考虑应力约束的优点,将应力约束设计与传统的基于刚度的设计进行了比较. 结果表明, 改进的BESO方法优化迭代过程稳健,获得了边界灰度单元很少的清晰的拓扑构型,并实现了有效降低应力集中效应的设计.      

Design of buckling constrained multiphase material structures using continuum topology optimization

This study contributes to a possibility of evaluating composite structures configuration such as steel and concrete using buckling and volume constraints based on multi-material topology optimization. A Jacobi active-phase algorithm is used to generate multiphase topology optimization. It provides a rational solution appropriated to the topology optimizer, Method of Moving Asymptotes due to the conflict in updating the design variables. A modified material interpolation scheme solving spurious buckling modes problem which occurs in the multi-material topology optimization process is given and discussed. An investigation of buckling constraint parameter is described. It allows a single-objective minimum compliance topology optimization to obtain two objectives of maximizing both structure stiffness and first buckling load factor. The optimal changing topologies of single material structure and multi-material structure corresponding to different buckling constraints are presented. Numerical examples of compression-only structures and compression-tension structures considering structural instability are performed using both single material and multiple materials to verify the efficiency and superiority of the present method.

     

离散变量结构拓扑优化设计研究

研究了离散变量结构拓扑优化设计的若干问题，讨论了离散型优化模型的合理性，提出截面设计变量的离散程度和全局约束影响最优拓扑，是优化中不可忽视的因素，文中还提出了一种解离散变量桁架，刚架结构拓扑优化的启发式算法。     

Material microstructure optimization for linear elastodynamic energy wave management

We describe a systematic approach to design material microstructures to achieve desired energy propagation in a two-phase composite plate. To generate a well-posed topology optimization problem we use the relaxation approach which requires homogenization theory to relate the macroscopic material properties to the microstructure, here a sequentially ranked laminate. We introduce an algorithm whereby the laminate layer volume fractions and orientations are optimized at each material point. To resolve numerical instabilities associated with the dynamic simulation and constrained optimization problem, we filter the laminate parameters. This also has the effect of generating smoothly varying microstructures which are easier to manufacture. To demonstrate our algorithm we design microstructure layouts for tailored energy propagation, i.e. energy focus, energy redirection, energy dispersion and energy spread.     

一种基于应力的双方向结构拓扑优化算法

基于传统的渐进结构优化方法，提出了一种基于应力的双方向渐进拓扑优化算法。该方法是对传统方法和目前的双方向法的改进和完善。算例表明该方法能避免目前的双方向法具有的解的振荡问题，具有更高的可靠性，能获得更佳的拓扑结构。     

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

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

Studies on optimal topology design of structures with discrete variables

Some problems in the optimal topology design of structures with discrete variables are studied in this paper. The problem of a model of discrete optimization is discussed and a neglected fact that discrete optimum design may be controlled by the discreteness of sizing variables and global constraints is pointed out. A heuristic algorithm for solving discrete topology optimization problems of trusses and frames is proposed.     

考虑最小尺寸精确控制的SIMP和MMC混合拓扑优化方法

拓扑优化作为一种先进设计方法, 已被成功用于多个学科领域优化问题求解, 但从拓扑优化结果到其工程应用之间仍存在诸多阻碍, 如在结构设计中存在难以制造的小孔或边界裂缝和单铰链连接等. 在拓扑优化设计阶段考虑结构最小尺寸控制是解决上述问题的一种有效手段. 在最小尺寸控制的结构拓扑优化方法中, 通用性较强的固体各向同性材料惩罚法SIMP优化结果边界模糊不光滑, 包含精确几何信息的移动变形组件法MMC对初始布局具有较强依赖性. 本文提出一种考虑最小尺寸精确控制的SIMP和MMC混合拓扑优化方法. 所提方法继承了二者优势, 避免了各自缺点. 在该方法中, 首先采用活跃轮廓算法ACWE获取SIMP输出的拓扑结构边界轮廓数据, 提出了SIMP优化结果到MMC组件初始布局的映射方法. 其次, 通过引入组件的3个长度变量, 建立了半圆形末端的多变形组件拓扑描述函数模型. 最后, 以组件厚度变量为约束, 构建了考虑结构最小尺寸控制的拓扑优化模型. 采用最小柔度问题和柔性机构问题验证了所提方法的有效性. 数值结果表明, 所提方法在无需额外约束的条件下, 仅通过组件厚度变量下限设置, 可实现整体结构的最小尺寸精确控制, 并获得了具有全局光滑的拓扑结构边界.