47应用拓扑优化技术降低柴油机的噪声---方华

应用拓扑优化技术降低柴油机的噪声

方华1孔令彬1 宫传刚2

（1-吉林大学汽车学院，2-道依茨一汽（大连）柴油机有限公司。）

应用拓扑优化技术降低柴油机的噪声Reducing Diesel Engine Noise with Topology

Optimization

方华1孔令彬1 宫传刚2

(1-吉林大学汽车学院，2-道依茨一汽（大连）柴油机有限公司。）

摘要：本文应用Altair HyperWorks 的OptiStruct优化模块，采用拓扑优化技术，对某柴油机正时齿罩进行了优化，得到了最佳的质量分布，据此进行了结构改进，使该正时齿罩的模态振型避开了正时轮系啮合激励频段，进而降低辐射噪声。同时利用SYSNOISE软件进行了辐射声功率的计算，得到了理论上的验证。

关键词：拓扑优化，啮合频率，模态，噪声

Abstract: In this paper, the optimization of timing gear cover was done with Altair HyperWorks, the optimal mass distribution was gained, so the structural modification could be carried out. The modes of timing gear cover have avoided the frequency range of gear train excitation, and consequently the radiated noise would be reduced. At the same time, the radiated sound power was calculated by SYSNOISE software, the theoretic verification was done.

Keywords: T opology Optimization, engaging frequency, modal analysis noise

1 概述

某柴油机前端声强（1/3倍频程）测试表明，2000 Hz中心频率的辐射噪声最大，位置如图1（a）中的圆圈区域。而此时正时轮系的啮合频率为1813 ±53Hz，恰好与2000Hz 的辐射噪声相对应。为此应使罩盖的模态固有频率避开这些频段，尽可能减少透射噪声以及共振的可能性。

2 有限元建模及模态分析

正时齿轮罩的三维实体模型及有限元模型如图1所示，有限元网格采用六面体网格划分，单元类型为solid45，共有节点74278个,单元56145个。材料参数取自材料测试结果，弹性模量61GPa，泊松比0.3，密度2.68Kg/cm3。

基于此模型进行了有限元自由模态分析，并与实验模态的结果进行了对比，各阶固有振型误差小于6％（图表略），验证了有限元模型的正确性。

约束正时齿轮罩与机体的结合面及螺栓孔，进行工作状态下的模态分析，计算模型及结果如图2所示。由图可知，原型在1800~2000Hz左右的模态有三阶，且位置与声强测试的结果相吻合。为此要进行结构优化，减少或消除该频率范围内的模态。

(

(a) (b)

图1 正时齿轮罩盖的实体(a)与有限元模型(b)

Figure 1 Entity of timing gear cover (a) and its FE model (b)

第一阶振型（1271Hz）第二阶振型（1848Hz）

第三阶振型（2041Hz）第四阶振型（2165Hz）

图2 原型的工作状态模态分析结果

Figure 2 Results of working modal analysis on the original

3 拓扑优化分析及优化结构的确定

基于以上分析，同时考虑实际的安装等因素，建立了罩盖的优化分析模型，如图3(a)，

其中优化区域如图3(b)所示。

(a) (b)

图3 优化分析模型(a)与优化区域(b)

Figure 3 Optimization analysis model (a) and optimization area

设置设计变量为优化区域；响应为前两阶固有频率及体积比；目标函数为体积比最小；

约束条件为第一，二阶模态固有频率。利用OptiStruct优化模块，采用拓扑优化技术进行了

分析，得到如图4所示的优化结果。

图

4 优化区域拓扑图

Figure 4 Topology map of optimization area

根据优化区域拓扑图，结合零件的加工及装配约束，对原始模型进行结构修改，得到正时齿罩的最终优化模型，如图5所示，优化模型共有节点75764个,单元57670个。并进行了工作模态分析，结果表明，优化模型位于1800~2000Hz左右的固有模态已消除，2000Hz 以上的模态频率值提高了300Hz左右，基本避开了中心频率为2000Hz的共振频段。

图5 优化后模型的工作状态模态分析结果

Figure 5 Results of working modal analysis on the optimized timing gear cover

4 辐射声功率的计算

为了考察改进前后的罩盖的辐射噪声情况，首先进行有限元的响应分析，得到了表面的响应位移，为辐射声场的分析提供了输入边界。又建立了图6所示的边界元模型，在正时齿轮罩盖与机体连接的内表面添加一个对称平面来形成闭合的空间，该平面为刚性平面，其法向速度为零，将正时齿轮罩盖的辐射噪声完全反射。图中的平面网格为域点网格，距离正时齿轮罩盖模型的外表面为1m，用于考察1m噪声情况。选取1000~3000Hz频域进行边界元分析。

图7为边界元分析得到的辐射声功率的频谱，红线为原型，蓝线为优化后。结果表明，优化模型与原始模型相比，在1800-2200Hz的范围内，输出声功率明显降低，并且其他峰值也比原型的最高峰值低3dB以上，总输出声功率降低了2.4dB。

（a）(b)

图6 原型边界元模型（a ）和优化结构边界元模型（a ）

Figure 6 Original BEM model (a) and optimization BEM model (b)

图7 输出声功率频谱

Figure 7 Output sound power spectrum

5 结论

利用Hyperworks 软件的建模及优化功能模块，对正时齿轮罩盖进行优化，可以达到很好的效果。优化后的正时齿轮罩盖模态分析表明，结构的固有频率已经远离了齿频激振频域，且辐射声功率降低了2.4dB ，达到了预期的效果。

本文的研究过程，代表了低噪声结构零部件的开发流程，即包括试验分析研究，模拟仿真分析，结构优化设计等。对零部件的声学优化具有一定的指导意义。

6 参考文献

1. Laurent Moulin Simulation Study of the Internal Mechanical Excitations to Reduce the Noise Emitted by a Heavy Duty Engine SAE 2003-01-1727

2． 于开平 周传月 谭惠丰等编著 HyperMesh 从入门到精通 科学出版社 2005.5 3． 周传月 腾万秀 张俊堂 主编 工程有限元与优化分析应用实例教程 科学出版社 2005.6

工程结构拓扑优化的理论研究及应用_满宏亮.

提要 本文首先介绍了国内外拓扑优化技术的研究发展现状，讨论了拓扑优化的 原理、方法以及各种拓扑优化算法。其次，着重研究了SIMP 材料插值方法，建立了基于SIMP 理论的连续体结构拓扑优化模型，选取准则优化法对其密 度迭代格式进行了推导；并且利用MATLAB软件编程实现，有效地进行了平面结构的分析和拓扑优化设计。然后，分析了拓扑优化中的数值计算不稳定性现象，研究了能够有效消除拓扑优化中的数值计算不稳定性现象的各种解决方法，并对其进行了比较。最后，利用连续体结构拓扑优化求解理论和算法，使用结构有限元分析软件Hyperworks 对具体工程结构部件进行了拓扑优化设计研究，成功地应用到了实际工程问题中，算例结果表明了该优化方法的有效性和正确性。 关键词：有限元拓扑优化材料插值模型数值计算不稳定性优化 求解算法 Key words: FEA Topology optimization Material Interpolation Model Numerical Calculation Instabilities Optimization Solution Algorithm-i- 目录 第一章绪论 (1) 1.1 前言 (1) 1.2 国内外拓扑优化研究概况 (3) 1.3 本文研究内容及意义 (9) 第二章现代结构拓扑优化理论 (11) 2.1 拓扑的概念 (11) 2.1.1 拓扑学的由来 (11) 2.1.2 拓扑学及拓扑性质 (13) 2.2 结构拓扑优化原理和方法 (16) 2.2.1 拓扑优化的基本原理 (17) 2.2.2 结构拓扑优化设计方法 (17) 2.2.3 拓扑优化设计方法比较 (21) 2.3 拓扑优化设计的优化算法概述 (22) 2.3.1 优化算法分类 (22) 2.3.2 拓扑优化常用算法 (24) 第三章连续体结构拓扑优化的模型建立与求解算法 (27) 3.1 连续体结构拓扑优化设计的模型描述 (29) 3.2 数学模型的有限元离散 (34) 3.2.1 单元应变和应力.........................................34吉林大学硕士研究生学位论文-ii- 3.2.2 单元平衡方程 (35) 3.2.3 连续体结构拓扑优化的数学模型的有限元离散形式 (38) 3.3 基于SIMP 理论的优化准则法 (39) 第四章结构拓扑优化程序实现 (45) 4.1 基于SIMP 理论的优化准则法迭代分析流程 (45) 4.2 优化过程的MA TLAB 编程实现 (47) 4.3 计算实例 (48)

Abaqus中Topology和Shape优化指南

Abaqus中Topology和Shape优化指南

目录 1. 优化模块界面......................................................................................................- 1 - 2. 专业术语..............................................................................................................- 1 - 3.定义拓扑优化Task(general optimization和condition-based optimization).......- 2 - 3.1 General Optimization 参数设置.................................................................- 3 - 3.1.1 Basic选项参数..................................................................................- 3 - 3.1.2 Density选项参数..............................................................................- 4 - 3.1.3 Perturbation选项参数.......................................................................- 5 - 3.1.4 Advanced选项参数...........................................................................- 5 - 3.2 Condition-based topology Optimization 参数设置....................................- 6 - 3.2.1 Basic选项参数..................................................................................- 7 - 3.2.2 Advanced选项参数...........................................................................- 7 - 4 定义Shape Optimization Task方法....................................................................- 8 - 4.1 Basic选项参数............................................................................................- 8 - 4.2 Mesh Smoothing Quality选项参数............................................................- 9 - 4.3 Mesh Smoothing Quality选项参数..........................................................- 11 - 5 定义design response变量方法.........................................................................- 13 - 5.1 单个design response定义方法...............................................................- 14 - 5.2 combined design response定义方法........................................................- 15 - 5.3 design response使用注意事项.................................................................- 17 - 5.3.1 定义design response的操作.........................................................- 17 - 5.3.2 condition-based topology optimization的design response............- 18 - 5.3.3 general topology optimization的design response..........................- 18 - 5.3.4 design response for shape optimization...........................................- 21 - 6 定义objective function方法..............................................................................- 22 - 6.1 目标函数定义...........................................................................................- 23 - 6.2 目标函数的运算.......................................................................................- 23 - 6.2.1 min运算..........................................................................................- 23 - 6.2.2 max运算..........................................................................................- 24 -

北航拓扑优化程序学习报告

拓扑优化的 99行程序学习报告4月19日2011 《结构优化设计》课程学习报告 任课教师：李书

一、前言： 在最近的结构优化设计课程上学习了O.Sigmund的《A 99 line topology optimization code written in Matlab》一文，对拓扑优化的理论原理与实际的计算机程序实现都有了一定的理解，文章主要是通过拓扑优化的原理来实现对简单结构的静力学问题的优化求解，而编写的代码仅有99行，包括36行的主程序，12行的OC优化准则代码，16行的网格过滤代码和35行的有限元分析代码。 自1988 年丹麦学者Bendsoe与美国学者Kikuchi提出基于均匀化方法的结构拓扑优化设计基本理论以来,均匀化方法应用到具有周期性结构的材料分析中,近几年该方法已经成为分析夹杂、纤维增强复合材料、混凝土材料等效模量,以及材料的细观结构拓扑优化常用的手段之一。其基本思想是在组成拓扑结构的材料中引入微结构，优化过程中以微结构的几何尺寸作为设计变量，以微结构的消长实现其增删，并产生介于由中间尺寸微结构组成的复合材料，从而实现了结构拓扑优化模型与尺寸优化模型的统一。文章就是通过均匀化的基础，结合拓扑结构优化的工程实际，以计算机模拟的方法将拓扑优化的一般过程呈现出来，有助于初涉拓扑优化的读者对拓扑优化有个基础的认识。 二、拓扑优化问题描述 为了简化问题的描述，文中假设设计域是简单的矩形形式，且在进行有限元离散的时候采用正方形单元对其进行离散。这样不仅便于进行单元离散和单元编号，也利于对结构进行几何外形的描述。 一般说来，基于指数逼近法的拓扑优化最小化的问题可作如下描述： 文中采用的对结构材料属性的描述是所谓的“指数逼近法”或者称为SIMP 逼近法，即（Solid Isotropic Material with Penalization带惩罚因子的各项同性材料模型法），该方法是拓扑优化中常用的变密度材料插值模型中最具代表性的一种。

ABAQUS拓扑优化手册

ABAQUS拓扑优化分析手册/用户手册 分析手册： 13. Optimization Techniques优化技术 13.1 结构优化：概述 13.1.1 概述 ABAQUS结构优化是一个帮助用户精细化设计的迭代模块。结构优化设计能够使得结构组件轻量化，并满足刚度和耐久性要求。ABAQUS提供了两种优化方法——拓扑优化和形状优化。拓扑优化（Topology optimization）通过分析过程中不断修改最初模型中指定优化区域的单元材料性质，有效地从分析的模型中移走/增加单元而获得最优的设计目标。形状优化（Shape optimization）则是在分析中对指定的优化区域不断移动表面节点从而达到减小局部应力集中的优化目标。拓扑优化和形状优化均遵从一系列优化目标和约束。 最优化方法（Optimization）是一个通过自动化程序增加设计者在经验和直觉从而缩短研发过程的工具。想要优化模型，必须知道如何去优化，仅仅说要减小应力或者增大特征值是不够，做优化必须有更专门的描述。比方说，想要降低在两种不同载荷工况下的最大节点力，类似的还有，想要最大化前五阶特征值之和。这种最优化的目标称之为目标函数（Object Function）。另外，在优化过程中可以同时强制限定某些状态参量。例如，可以指定某节点的位移不超过一定的数值。这些强制性的指定措施叫做约束（Constraint）。 ABAQUS/CAE可以创建模型然后定义、配置和执行结构优化。更多信息请参考用户手册第十八章。 13.1.2 术语（Terminology） 设计区域（Design area）: 设计区域即模型需要优化的区域。这个区域可以是整个模型，也可以是模型的一部分或者数部分。一定的边界条件、载荷及人为约束下，拓扑优化通过增加/删除区域中单元的材料达到最优化设计，而形状优化通过移动区域内节点来达到优化的目的。 设计变量（Design variables）：设计变量即优化设计中需要改变的参数。拓扑优化中，设计区域中单元密度是设计变量，ABAQUS/CAE优化分析模块在其优化迭代过程中改变单元密度并将其耦合到刚度矩阵之中。实际上，拓扑优化将模型中单元移除的方法是将单元的质量和刚度充分变小从而使其不再参与整体结构响应。对于形状优化而言，设计变量是指设计区域内表面节点位移。优化时，ABAQUS或者将节点位置向外移动或者向内移动，抑或不移动。在此过程中，约束会影响表面节点移动的多少及其方向。优化仅仅直接修改边缘处的节点，而边缘内侧的节点位移通过边缘处节点插值得到。 设计循环（Design cycle）: 优化分析是一种不断更新设计变量的迭代过程，执行ABAQUS进行模型修改、查看结果以及确定是否达到优化目的。其中每次迭代叫做一个设计循环。 优化任务（Optimization task）：一次优化任务包含优化的定义，比如设计响应、目标、限制条件和几何约束。 设计响应（Design responses）: 优化分析的输入量称之为设计响应。设计响应可以直接从ABAQUS的结果输出文件.odb中读取，比如刚度、应力、特征频率及位移等。或者ABAQUS 从结果文件中计算得到模型的设计响应，例如质心、重量、相对位移等。一个设计响应与模型紧密相关，然而，设计响应存在一定的范围，例如区域内的最大应力或者模型体积。另外，设计响应也与特点的分析步和载荷状况有关。 目标函数（Objective functions）: 目标函数决定了优化的目标。一个目标函数是从设计响应中萃取的一定范围内的值，如最大位移和最大应力。一个目标函数可以用多个设计响应

ANSYS拓扑优化原理讲解以及实例操作

拓扑优化是指形状优化，有时也称为外型优化。 拓扑优化的目标是寻找承受单载荷或多载荷的物体的最佳材料分配方案。这种方案在拓扑优化中表现为“最大刚度”设计。与传统的优化设计不同的是，拓扑优化不需要给出参数和优化变量的定义。目标函数、状态变量和设计变量（参见“优化设计”一章）都是预定义好的。用户只需要给出结构的参数（材料特性、模型、载荷等）和要省去的材料百分比。给每个有限元的单元赋予内部伪密度来实现。这些伪密度用PLNSOL ，TOPO 命令来绘出。拓扑优化的目标——目标函数——是在满足结构的约束（V ）情况下减少结构的变形能。减小结构的变形能相当于提高结构的刚度。这个技术通过使用设计变量。 结构拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料分布的问题。通过拓扑优化分析，设计人员可以全面了解产品的结构和功能特征，可以有针对性地对总体结构和具体结构进行设计。特别在产品设计初期，仅凭经验和想象进行零部件的设计是不够的。只有在适当的约束条件下，充分利用拓扑优化技术进行分析，并结合丰富的设计经验，才能设计出满足最佳技术条件和工艺条件的产品。连续体结构拓扑优化的最大优点是能在不知道结构拓扑形状的前提下，根据已知边界条件和载荷条件确定出较合理的结构形式，它不涉及具体结构尺寸设计，但可以提出最佳设计方案。拓扑优化技术可以为设计人员提供全新的设计和最优的材料分布方案。拓扑优化基于概念设计的思想，作为结果的设计空间需要被反馈给设计人员并做出适当的修改。最优的设计往往比概念设计的方案结构更轻，而性能更佳。经过设计人员修改过的设计方案可以再经过形状和尺寸优化得到更好的方案。 5.1.2优化拓扑的数学模型 优化拓扑的数学解释可以转换为寻求最优解的过程，对于他的描述是：给定系统描述和目标函数，选取一组设计变量及其范围，求设计变量的值，使得目标函数最小（或者最大）。一种典型的数学表达式为： ()()()12,,0,,0min ,g x x v g x x v f x v ?=??≤???? 式中，x -系统的状态变量；12g g 、-一等式和不等式的结束方程；(),f x v -目标函数；v -设计变量。 注：在上述方程中，x 作为系统的状态变量，并不是独立的变量，它是由设计变量得出的，并且与设计变量相关。 优化拓扑所要进行的数学运算目标就是，求取合适的设计变量v ，并使得目标函数值最小。 5.2基于ANSYS 的优化拓扑的一般过程 （进行内容排版修改） 在ANSYS 中，进行优化拓扑，一般分为6个步骤。具体流程见图5-1：

如何利用ANSYS进行拓扑优化

如何利用ANSYS进行拓扑优化 前言 就目前而言，利用有限元进行优化主要分成两个阶段： （1）进行拓扑优化，明确零件最佳的外形、刚度、体积，或者合理的固有频率，主要目的是确定优化的方向； （2）进行尺寸优化，主要目的是确定优化后的的零件具体尺寸值，通常是在完成拓扑优化之后，再执行尺寸优化。 在ANSYS中，利用拓扑优化，可以完成以下两个目的： （1）在特定载荷和约束的条件下，确定零件的最佳外形，或者最小的体积（或者质量）； （2）利用拓扑优化，使零件达到需要的固有频率，避免在使用过程中产生共振等不利影响。 本文主要就在ANSYS环境中如何执行拓扑优化进行说明。

1、利用ANSYS进行拓扑优化的过程 在ANSYS中，执行优化，通常分为以下6个步骤： 、定义需要求解的结构问题 对于结构进行优化分析，定义结构的物理特性必不可少，例如，需要定义结构的杨氏模量、泊松比（其值在～之间）、密度等相关的结构特性方面的信息，以供结构计算能够正常执行下去。

、选择合理的优化单元类型 在ANSYS中，不是所有的单元类型都可以执行优化的，必须满足如下的规定： （1）2D平面单元：PLANE82单元和PLANE183单元； （2）3D实体单元：SOLID92单元和SOLID95单元； （3）壳单元：SHELL93单元。 上述单元的特性在帮助文件中有详细的说明，同时对于2D单元，应使用平面应力或者轴对称的单元选项。 、指定优化和非优化的区域 在ANSYS中规定，单元类型编号为1的单元，才执行优化计算；否则，就不执行优化计算。例如，对于结构分析中，对于不能去除的部分区域将单元类型编号设定为≥2，就可以不执行优化计算，请见下面的代码片段：…… …… Et,1,solid92 Et,2,solid92 …… Type,1 Vsel,s,num,,1,2 Vmesh,all …… Type,2 Vsel,s,num,,3 Vmesh,all ……

连续体结构拓扑优化方法及存在问题分析

编号：SY-AQ-00556 ( 安全管理) 单位：_____________________ 审批：_____________________ 日期：_____________________ WORD文档/ A4打印/ 可编辑 连续体结构拓扑优化方法及存 在问题分析 Topology optimization method of continuum structure and analysis of existing problems

连续体结构拓扑优化方法及存在问 题分析 导语：进行安全管理的目的是预防、消灭事故，防止或消除事故伤害，保护劳动者的安全与健康。在安全管理的四项主要内容中，虽然都是为了达到安全管理的目的，但是对生产因素状态的控制，与安全管理目的关系更直接，显得更为突出。 文章深入分析国内外连续体结构拓扑优化的研究现状，介绍了拓扑优化方法的发展及实现过程中存在的问题。对比分析了均匀化方法，渐进结构优化法，变密度法的优缺点。研究了连续体结构拓扑优化过程中产生数值不稳定现象的原因，重点讨论了灰度单元，棋盘格式，网格依赖性的数值不稳定现象，并针对每一种数值不稳定现象提出了相应的解决办法。 结构拓扑优化设计的主要对象是连续体结构，1981年程耿东和Olhof在研究中指出：为了得到实心弹性薄板材料分布的全局最优解，必须扩大设计空间，得到由无限细肋增强的板设计。此研究被认为是近现代连续体结构拓扑优化的先驱。 目前，国内外学者对结构拓扑优化问题进行了大量研究，这些

研究大多数建立在有限元法结构分析的基础上，但由于有限元法中单元网格的存在，结构拓扑优化过程中常常出现如灰度单元，网格依赖性和棋盘格等数值不稳定的现象。本文介绍了几种连续体结构拓扑优化方法及每种方法存在的问题，并提出了相应的解决办法。 1.拓扑优化方法 连续体结构拓扑优化开始于1988年Bendoe和Kikuchi提出的均匀化方法，此后许多学者相继提出了渐进结构优化方法、变密度法等拓扑优化数学建模方法。 1.1.均匀化方法 均匀化方法即在设计区域内构造周期性分布的微结构，这些微结构是由同一种各向同性材料实体和孔洞复合而成。采用有限元方法进行分析，在每个单元内构造不同尺寸的微结构，微结构的尺寸和方向为拓扑优化设计变量。1988年Bendsoe研究发现，通过在结构中引入具有空洞微结构的材料模型，将困难的拓扑设计问题转换为相对简单的材料微结构尺寸优化问题。 很多学者发展了均匀化方法，Suzhk进行了基于均匀化方法结

拓扑优化经典99行程序解读

3188-1-1.html Sigmund教授所编写的top优化经典99行程序，可以说是我们拓扑优化研究的基础； 每一个新手入门都会要读懂这个程序，才能去扩展，去创新； 99行程序也有好多个版本，用于求解各种问题，如刚度设计、柔顺机构、热耦合问题，但基本思路大同小异； 本文拟对其中的一个版本进行解读，愿能对新手有点小小的帮助。 不详之处，还请论坛内高手多指点 读懂了该程序，只能说是略懂拓扑优化理论了， 我手里就有一些水平集源程序是成千上万行，虽然在99行的基础上成熟了很多，但依然还有很多的发展空间。 源程序如下： %%%% A 99 LINE TOPOLOGY OPTIMIZATION CODE BY OLE SIGMUND, JANUARY 2000 %%% %%%% CODE MODIFIED FOR INCREASED SPEED, September 2002, BY OLE SIGMUND %%% function top(nelx,nely,volfrac,penal,rmin); nelx=80; nely=20; volfrac=0.4; penal=3; rmin=2; % INITIALIZE x(1:nely,1:nelx) = volfrac; loop = 0; change = 1.; % START ITERATION while change > 0.01 loop = loop + 1; xold = x; % FE-ANAL YSIS [U]=FE(nelx,nely,x,penal); % OBJECTIVE FUNCTION AND SENSITIVITY ANAL YSIS [KE] = lk; c = 0.; for ely = 1:nely for elx = 1:nelx n1 = (nely+1)*(elx-1)+ely; n2 = (nely+1)* elx +ely; Ue = U([2*n1-1;2*n1; 2*n2-1;2*n2; 2*n2+1;2*n2+2; 2*n1+1;2*n1+2],1); c = c + x(ely,elx)^penal*Ue'*KE*Ue; dc(ely,elx) = -penal*x(ely,elx)^(penal-1)*Ue'*KE*Ue; end end

连续体结构拓扑优化方法及存在问题分析(最新版)

( 安全管理 ) 单位：_________________________ 姓名：_________________________ 日期：_________________________ 精品文档 / Word文档 / 文字可改 连续体结构拓扑优化方法及存在问题分析(最新版) Safety management is an important part of production management. Safety and production are in the implementation process

连续体结构拓扑优化方法及存在问题分析 (最新版) 文章深入分析国内外连续体结构拓扑优化的研究现状，介绍了拓扑优化方法的发展及实现过程中存在的问题。对比分析了均匀化方法，渐进结构优化法，变密度法的优缺点。研究了连续体结构拓扑优化过程中产生数值不稳定现象的原因，重点讨论了灰度单元，棋盘格式，网格依赖性的数值不稳定现象，并针对每一种数值不稳定现象提出了相应的解决办法。 结构拓扑优化设计的主要对象是连续体结构，1981年程耿东和Olhof在研究中指出：为了得到实心弹性薄板材料分布的全局最优解，必须扩大设计空间，得到由无限细肋增强的板设计。此研究被认为是近现代连续体结构拓扑优化的先驱。 目前，国内外学者对结构拓扑优化问题进行了大量研究，这些

研究大多数建立在有限元法结构分析的基础上，但由于有限元法中单元网格的存在，结构拓扑优化过程中常常出现如灰度单元，网格依赖性和棋盘格等数值不稳定的现象。本文介绍了几种连续体结构拓扑优化方法及每种方法存在的问题，并提出了相应的解决办法。 1.拓扑优化方法 连续体结构拓扑优化开始于1988年Bendoe和Kikuchi提出的均匀化方法，此后许多学者相继提出了渐进结构优化方法、变密度法等拓扑优化数学建模方法。 1.1.均匀化方法 均匀化方法即在设计区域内构造周期性分布的微结构，这些微结构是由同一种各向同性材料实体和孔洞复合而成。采用有限元方法进行分析，在每个单元内构造不同尺寸的微结构，微结构的尺寸和方向为拓扑优化设计变量。1988年Bendsoe研究发现，通过在结构中引入具有空洞微结构的材料模型，将困难的拓扑设计问题转换为相对简单的材料微结构尺寸优化问题。 很多学者发展了均匀化方法，Suzhk进行了基于均匀化方法结构

abaqus常用技巧总结

a b a q u s常用技巧总结-标准化文件发布号：（9556-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

Abaqus常用技巧总结 本手册是由simwe等论坛上精华帖以及本人下载的其他资料整理，由于很多资料搜集已经很久，而且时间有限，都没有注明原作者，也没有最资料进行分类整理，见谅。如需要，请PM给我。 Shelly31 Python.tzy@https://www.wendangku.net/doc/307972536.html, 2007.8.1 建议阅读方式：

目录 ABAQUS常用技巧总结 (2) 目录 (3) 1.对TIME INCREMENT的根本理解 (5) 2.ABAQUS 请问 MOMENT的加载 (5) 3.ABAQUS计算时C盘的临时文件太大了，怎么改目录？ (6) 4.CAE中如何加预应力 (6) 5.HYPERMESH里面看到ABAQUS分析的结果 (6) 6.X-Y PLOTS (6) 7.把上一次的分析结果作为下一次分析的初始条件该怎么做 (7) 8.材料方向与增量步 (8) 9.多个INP文件如何实现批处理 (9) 10.关于ABAQUS的任务管理 (10) 11.关于数据的输入输出 (12) 12.后处理积分 (12) 13.接触分析激活杀死 (13) 14.利用QUEUE的功能由本地机器向远程UNIX机器提交ABAQUS作业的方法[精华] (14) 15.利用命令进行计算时如何设置调用内存量 (17) 16.清华大学BBS的ABAQUS精华 (17) 17.请问怎么实现双曲线 (55)

18.取消坐标系等的显示 (56) 19.如何在计算中修改材料特性 (57) 20.输出计算过程中的总质量和总刚度矩阵 (60) 21.先张预应力： (61) 22.用户子程序的使用 (61) 23.怎样设定用双CPU机器进行ABAQUS计算 (61) 24.中途停止正在运算的JOB (62) 25.自适应网格技术 (62) 26.ABAQUS计算与内存 (63) 27.质量缩放 (64) 28.ABAQUS多处理器进行并行计算的效果研究 (79) 29.YAHOO讨论组摘录--CONTACT+OVERCLOSURE (81) 30.原创：无限元建立方法，希望得到加分 (95) 31.[分享]ABAQUS 使用问答 (102) 32.[转帖]ABAQUS6.4导入外来模型的几点小经验！ (122) 33.ABAQUS的多图层绘图 (125) 34.子结构 (125) 35.如何在不同的分析步改变材料的参数 (126) 36.模型的重启动分析－RESTART (127) 37. ABAQUS的单位心得 (128)

拓扑优化技术

拓扑优化技术 第1节基本知识 一、拓扑优化的概念 拓扑优化是指形状优化，有时也称为外型优化。拓扑优化的目标是寻找承受单载荷或多载荷的物体的最佳材料分配方案。这种方案在拓扑优化中表现为“最大刚度”设计。 与传统的优化设计不同的是，拓扑优化不需要给出参数和优化变量的定义。目标函数、状态变量和设计变量都是预定义好的。用户只需要给出结构的参数（材料特性、模型、载荷等）和要省去的材料百分比。 拓扑优化的目标—目标函数—是在满足结构的约束（V）情况下减少结构的变形能。减小结构的变形能相当于提高结构的刚度。这个技术通过使用设计变量（ i）给每个有限元的单元赋予内部伪密度来实现。这些伪密度用PLNSOL，TOPO命令来绘出。 ANSYS提供的拓扑优化技术主要用于确定系统的最佳几何形状，其原理是系统材料发挥最大利用率，同时确保系统的整体刚度（静力分析）、自振频率（模态分析）在满足工程要求的条件下获得极大或极小值。 拓扑优化应用场合：线性静力分析和模态分析。 拓扑优化原理：满足结构体积缩减量的条件下使目标函数结构柔量能量（the enery of structure compliance—SCOMP）的极小化。结构柔量能量极小化就是要求结构刚度的最大化。 例如，给定V=60表示在给定载荷并满足最大刚度准则要求的情况下省去60%的材料。图19-1表示满足约束和载荷要求的拓扑优化结果。图19-1a表示载荷和边界条件，图19-b 表示以密度云图形式绘制的拓扑结果。 图19-1 体积减少60%的拓扑优化示例 二、拓扑优化的基本过程 拓扑优化的基本步骤如下：

1．定义结构问题定义材料弹性模量、泊松系数、材料密度。 2．选择单元类型拓扑优化功能中的模型只能采用下列单元类型： ● 二维实体单元：Plane2和Plane82，用于平面应力问题和轴对称问题。 ● 三维实体单元：Solid92、Solid95。 ● 壳单元：SHELL93。 3．指定优化和不优化区域ANSYS只对单元类型编号为1的单元网格部分进行拓扑优 化，而对单元类型编号大于1的单元网格部分不进行拓扑优化，因此，拓扑优化时要确保进行拓扑优化区域单元类型编号为1，而不进行拓扑优化区域单元类型编号大于1即可。 4．定义并控制载荷工况或频率提取可以在单个载荷工况和多个载荷工况下做拓扑优化，单载荷工况是最简便的。 要在几个独立的载荷工况中得到优化结果时，必须用到写载荷工况和求解功能。在定义完每个载荷工况后，要用LSWRITE命令将数据写入文件，然后用LSSOLVE命令求解载荷工况的集合。 5．定义和控制优化过程拓扑优化过程包括定义优化参数和进行拓扑优化两个部分。用户可以用两种方式运行拓扑优化：控制并执行每一次迭代或自动进行多次迭代。 ANSYS有三个命令定义和执行拓扑优化：TOPDEF，TOPEXE和TOPITER。TOPDEF 命令定义要省去材料的量，要处理载荷工况的数目，收敛的公差；TOPEXE命令执行一次优化迭代；TOPITER命令执行多次优化迭代。 （1）定义优化参数首先要定义优化参数。用户要定义要省去材料的百分比，要处理载荷工况的数目，收敛的公差。 命令：TOPDEF GUI：Main Menu>Solution>Solve>Topological opt 注:本步所定义的内容并不存入ANSYS数据库中，因此在下一个拓扑优化中要重新使用TOPDEF命令。 （2）执行单次迭代定义好优化参数以后，可以执行一次迭代。迭代后用户可以查看收敛情况并绘出或列出当前的拓扑优化结果。可以继续做迭代直到满足要求为止。如果是在GUI方式下执行，在Topological Optimization 对话框（ITER域）中选择一次迭代。 命令：TOPEXE GUI：Main Menu>Solution>Solve>Topological opt TOPEXE的主要优点是用户可以设计自己的迭代宏进行自动优化循环和绘图。在下一节，可以看到TOPITER命令是一个ANSYS的宏，用来执行多次优化迭代。 （3）自动执行多次迭代 在定义好优化参数以后，用户可以自动执行多次迭代。在迭代完成以后，可以查看收敛情况并绘出或列出当前拓扑形状。如果需要的话，可以继续执行求解和迭代。TOPITER 命令实际是一个ANSYS的宏，可以拷贝和定制。

拓扑优化

一种新的优化方法——拓扑优化。是一种以多种使用条件为目标优化参数的优化方式，可以提高零件的真正使用效益，更加准确的反映了设计的优化过程。 优化设计可以在很大程度上改善和提高铸造件、锻造件和冲压件的性能，并减轻产品重量。然而，优化设计特别是拓扑优化很少应用在实际工程中。一方面是因为工程问题的复杂性和高度非线性，拓扑优化技术目前还无法实现这些系统优化问题，但更重要的是一门新的技术和方法很难取代人们已经习惯多年的思维模式和工作方式。 工程设计人员需要有更系统、更科学的设计思想和方法，以达到提高产品开发效率、节约原材料、降低成本及提高产品质量的目的，结构优化设计则是实现这些目的较佳手段[1]。由于设计变量类型的不同，结构优化设计可以分为由易到难的四个不同层次：尺寸优化、形状优化、形貌优化和拓扑布局优化。由于拓扑优化设计的难度较大，被公认为是当前结构优化领域内最具有挑战性的课题之一。但是在工程应用中，拓扑优化可以提供概念性设计方案，取得的经济效益比尺寸优化、形状优化更大，因此，拓扑优化技术对工程设计人员更具吸引力，已经成为当今结构优化设计研究的一个热点。 发动机运转期间，主轴承座承受多种载荷，这些载荷包括：螺栓预紧载荷、轴瓦过盈载荷及曲轴动载荷等。目前，主轴承座的主要评价指标是结构的强度、刚度是否满足设计需求。在明确主轴承座承载情况和设计要求的前提下，作者对某大马力发动机原有主轴承座进行了最大爆发压力工况下的有限元分析。分析模型及主轴承座轴瓦径向变形量见图1（a）、图1 （b）和图1（c）。通过主轴承座的强度分析和动态疲劳安全系数分析可以得知：主轴承座的动态疲劳安全系数为1.843，远远大于安全系数阀值1，所以主轴承座的强度足以满足设计需求。而从图1（b）可以得知轴瓦在变形后水平方向径向减小0.0739mm ，已经接近曲轴、轴瓦径向间隙最小值0.079mm，这容易导致曲轴与轴瓦间缺少油膜润滑，形成干摩擦，最终导致曲轴磨损加剧，发动机动载荷增加，甚至机毁人亡的悲剧；另外从图1（c）可以得知轴瓦在变形后上下方向径向增加0.0971mm ，小于轴瓦径向变形许可值0.147mm 。所以，根据有限元分析结果可以判断：主轴承座在水平方向的刚度不足够，应该改进现有结构，提高其刚度性能。

连续体结构拓扑优化方法评述_夏天翔

第2卷第1期2011年2月航空工程进展 A DV A N CES IN A ERON A U T ICA L SCIEN CE A N D EN GIN EERIN G Vo l 12N o 11Feb 1 2011 收稿日期:2010-12-01; 修回日期:2011-01-20基金项目:教育部长江学者创新团队项目(Irt0906)通信作者:姚卫星,w xyao@https://www.wendangku.net/doc/307972536.html, 文章编号:1674-8190(2011)01-001-12 连续体结构拓扑优化方法评述 夏天翔,姚卫星 (南京航空航天大学飞行器先进设计技术国防重点学科实验室,南京 210016) 摘 要:连续体结构拓扑优化在优化中能产生新的构型,对实现自动化智能结构设计具有重要意义。目前,连续体结构拓扑优化方法主要有:均匀化方法、变厚度法、变密度法、渐进结构优化方法、水平集法、独立连续映射方法。本文首先系统回顾了以上方法的发展历程,介绍了它们的研究现状。其次,通过对比以上拓扑优化方法对若干典型算例的优化结果,表明以上方法都有较好的减重效果。最后,对以上方法进行了总结,列出了它们的优缺点和发展方向。 关键词:拓扑优化;均匀化方法;变厚度法;变密度法;渐进结构优化方法;水平集法;独立连续映射方法中图分类号:V 211.7 文献标识码:A A Survey of Topology Optimization of Continuum Stru cture Xia Tianx iang ,Yao Weix ing (K ey L abor ator y of F undamental Science fo r N atio nal Defense -adv anced Design T echno lo gy of F lig ht V ehicle,Nanjing U niver sity o f A eronautics and A st ronautics,N anjing 210016,China) Abstract:A s the to po log y optim izat ion o f continuum structure can pr oduce new config ur atio ns during the optim-i zatio n,it is significant for automatic str ucture design.A t present,the most commo nly used t opolo gy o ptimiza -t ion methods of continuum st ructur e ar e:the ho mog enization method,var iable t hickness method,v ariable dens-i t y metho d,evo lutio nar y str uctur al o pt imizatio n met ho d,lev el set metho d,independent co ntinuous mapping method.Firstly,the develo pment pro cesses of above metho ds ar e sy stematically review ed,their cur rent r e -sear ch is br iefly intro duced in this paper.T hen,these methods ar e com par ed and discussed t hr ough a number of typical ex amples.T he typical ex amples show that all of above methods have gr eat abilities to r educe w eig ht.F-i nally ,the adv ant ag es,disadv ant ag es and dev elo pment directio ns of abov e metho ds ar e discussed. Key words:to po lo gy o ptimization;homog enizat ion metho d;va riable thickness method;var iable density method;evolutionar y structure optimization metho d;lev el set method;independent continuo us mapping method 0 引言 按照设计变量的不同,结构优化可分为以下三个层次:尺寸优化、形状优化和拓扑优化。结构拓 扑优化能在给定的外载荷和边界条件下,通过改变结构拓扑使结构在满足约束的前提下性能达到最优。与尺寸优化、形状优化相比,结构拓扑优化的经济效果更为明显,在优化中能产生新的构型,是 结构实现自动化智能设计所必不可少的。 按照优化对象的性质,拓扑优化可分为离散体拓扑优化和连续体拓扑优化两种。连续体拓扑优化与离散体拓扑优化相比,在应用范围更广的同 时,模型描述困难,设计变量多,计算量大。在过去很长一段时间里,连续体拓扑优化发展得十分缓慢,直到1988年Bendso e 等人[1] 提出均匀化方法之后,它才得到了迅速发展。目前,国内外学者对结构拓扑优化问题已经进行了大量研究[2-9]。目前最常用的连续体拓扑优化方法有均匀化方法、变厚 度法、变密度法、渐进结构优化方法(ESO)、水平集法(Level set)、独立连续映射方法(ICM)等。从拓

结构拓扑优化的发展现状及未来

结构拓扑优化的发展现状及未来 王超 中国北方车辆研究所一、历史及发展概况 结构拓扑优化是近20年来从结构优化研究中派生出来的新分支，它在计算结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化的工程应用研究还很不成熟，在国外处在发展的初期,尤其在国内尚属于起步阶段。1904 年Michell在桁架理论中首次提出了拓扑优化的概念。自1964 年Dorn等人提出基结构法,将数值方法引入拓扑优化领域，拓扑优化研究开始活跃。20 世纪80 年代初，程耿东和N. Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题，他们开创性的工作引起了众多学者的研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计，开创了连续体结构拓扑优化设计研究的新局面。1993年和提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002 年罗鹰等提出三角网格进化法，该方法在优化过程中实现了退化和进化的统一，提高了优化效率。 二、拓扑优化的工程背景及基本原理 通常把结构优化按设计变量的类型划分成三个层次：结构尺寸优化、形状优化和拓扑优化。尺寸优化和形状优化已得到充分的发展，但它们存在着不能变更结构拓扑的缺陷。在这样的背景下，人们开始研究拓扑优化。拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料的分布问题。寻求一个最佳的拓扑结构形式有两种基本的原理：一种是退化原理，另一种是进化原理。退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上，然后构造适当的优化模型，通过一定的优化方法逐步删减那些不必要的结构元素,直至最终得到一个最优化的拓扑结构形式。进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化，它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。 三、结构拓扑优化设计方法 目前常使用的拓扑优化设计方法可以分为两大类：退化法和进化法。 退化法即传统的拓扑优化方法，一般通过求目标函数导数的零点或一系列迭代计算过程求最优的拓扑结构。目前常用于拓扑优化的退化法有基结构方法、均匀化方法、变密度法、变厚度法等。 基结构方法（GSA）的思路是假定对于给定的桁架节点，在每两个节点之间用杆件连结起来得到的结构称为基结构。按照某种规则或约束,将一些不必要的杆件从基本结构中删除，认为最终剩下的构件决定了结构的最佳拓扑。基结构方法更适合于桁架和框架结构的拓扑优化。基结构法是在有限的子空间内寻优，容易丢失最优解，另外还存在组合爆炸、解的奇异性等问题。 均匀化方法（HA）引入微结构的单胞，通过优化计算确定其材料密度分布，并由此得出最优的拓扑结构。均匀化方法主要应用于连续体的拓扑优化设计，它不仅能用于应力约束和位移约束，也能用于频率约束。目前用均匀化方法来进行拓扑优化设计的有一般弹性问题、热传导问题、周期渐进可展曲面问题、非线性热弹性问题、振动问题和骨改造问题等。 变密度法是一种比较流行的力学建模方式，与采用尺寸变量相比，它更能反映拓