下载文档原格式
如何利用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,solid92Et,2,solid92……Type,1Vsel,s,num,,1,2Vmesh,all……Type,2Vsel,s,num,,3Vmesh,all……说明：上述代码片段定义相同的单元类型（solid92），但编号分别为1和2，并将单元类型编号1利用网格划分分配给了1＃体和2＃体，从而对其进行优化计算；而单元编号为2利用网格划分分配给了3＃体，从而不执行优化计算。
![一种基于MATLAB与ANSYS的联合优化设计方法[发明专利]](https://img.taocdn.com/s1/m/bf4de0e602020740bf1e9b52.png)
专利名称：一种基于MATLAB与ANSYS的联合优化设计方法专利类型：发明专利
发明人：周金宇,邱爽,曹清林,邱睿,刘浏
申请号：CN201710470799.4
申请日：20170620
公开号：CN107391796A
公开日：
20171124
专利内容由知识产权出版社提供
摘要：本发明提供了一种基于MATLAB与ANSYS的联合优化设计方法，该方法基于MATLAB和ANSYS的运行机制，并利用MATLAB和ANSYS各自的优势，把MATLAB作为主程序，在主程序中设置结构外形尺寸、材料参数、迭代条件及单元退化等，通过MATLAB主程序自动调用ANSYS分析的数据结果，实现了优化设计运算迭代的自动化；本发明能够使MATLAB在编程、运算、数据处理和图形处理方面的优势与ANSYS强大的有限元分析方面的权威性有机融合，可以有效的提高运算效率，实现了优化设计运算迭代的自动化，为结构优化设计提供新的设计途径。
申请人：江苏理工学院
地址：213001 江苏省常州市中吴大道1801号
国籍：CN
代理机构：常州佰业腾飞专利代理事务所(普通合伙)
代理人：刘娟娟
更多信息请下载全文后查看。
ANSYS教程二---拓扑优化拓扑优化什么是拓扑优化？拓扑优化是指形状优化，有时也称为外型优化。
拓扑优化的目标是寻找承受单载荷或多载荷的物体的最佳材料分配方案。
这种方案在拓扑优化中表现为“最大刚度”设计。
与传统的优化设计不同的是，拓扑优化不需要给出参数和优化变量的定义。
目标函数、状态变量和设计变量（参见“优化设计”一章）都是预定义好的。
用户只需要给出结构的参数（材料特性、模型、载荷等）和要省去的材料百分比。
拓扑优化的目标――目标函数――是在满足结构的约束（V）情况下减少结构的变形能。
减小结构的变形能相当于提高结构的刚度。
这个技术通过使用设计变量（i）给每个有限元的单元赋予内部伪密度来实现。
这些伪密度用PLNSOL，TOPO命令来绘出。
例如，给定V=60表示在给定载荷并满足最大刚度准则要求的情况下省去60%的材料。
图2-1表示满足约束和载荷要求的拓扑优化结果。
图2-1a表示载荷和边界条件，图2-2b表示以密度云图形式绘制的拓扑结果。
图2-1 体积减少60%的拓扑优化示例如何做拓扑优化拓扑优化包括如下主要步骤：1．定义拓扑优化问题。
2．选择单元类型。
3．指定要优化和不优化的区域。
4．定义和控制载荷工况。
5．定义和控制优化过程。
6．查看结果。
拓扑优化的细节在下面给出。
关于批处理方式和图形菜单方式不同的做法也同样提及。
定义拓扑优化问题定义拓扑优化问题同定义其他线性，弹性结构问题做法一样。
用户需要定义材料特性（杨氏模量和泊松比），选择合适的单元类型生成有限元模型，施加载ANSYS教程二---拓扑优化荷和边界条件做单载荷步或多载荷步分析。
参见“ANSYS Analysis Procedures Guides”第一、二章。
选择单元类型拓扑优化功能可以使用二维平面单元，三维块单元和壳单元。
要使用这个功能，模型中只能有下列单元类型：二维实体单元：SOLID2和__三维实体单元：__和__壳单元：__二维单元用于平面应力问题。
88-51-FB-5B-3A-66现代设计理论和方法大作业----基于MATLAB和ANSYS的悬臂梁拓扑优化指导老师：刘志刚项目组长：薛亚波项目成员：机自66学院：机械工程学院基于MATLAB和ANSYS的悬臂梁拓扑优化一、计划和任务安排表1.1学习时间安排表：1.2任务分配表：二、项目背景介绍及问题描述2.1项目背景及意义：2.1.1工程背景及基本原理：通常把结构优化按设计变量的类型划分成三个层次：结构尺寸优化、形状优化和拓扑优化。
尺寸优化和形状优化已得到充分的发展，但它们存在着不能变更拓扑结构的缺陷。
在这样的背景下，人们开始研究拓扑优化。
拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料的分布问题。
寻求一个最佳的拓扑结构形式有两种基本的原理：一种是退化原理，另一种是进化原理。
退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上，然后构造适当的优化模型，通过一定的优化方法逐步删减那些不必要的结构元素，直至最终得到一个最优化的拓扑结构形式。
进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化，它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。
2.1.2 本文研究意义：目前，结构优化大部分集中在尺寸设计变量(如板厚、杆的剖面积及管梁的直径)。
拓扑结构优化较尺寸优化复杂,但对于有些问题拓扑结构优化比尺寸优化有效,悬臂梁是其中的例子之一。
本文讨论悬臂梁的拓扑优化问题,围绕这一问题,怎样使结构具有最大刚度的设计占有相当重要的地位;怎样优化结构的形状使材料的分布，更加合理从而达到使结构具有最大刚度的目的是本文要研究的问题。
2.2研究现状2.2.1 理论研究现状：结构拓扑优化是近20年来从结构优化研究中派生出来的新分支，它在计算结构力学中已经被认为是最富挑战性的一类研究工作。
目前有关结构拓扑优化的工程应用研究还很不成熟，在国外处在发展的初期，尤其在国内尚属于起步阶段。
ansys workbench拓扑优化实例
拓扑优化（Topology Optimization）是一种结构优化的方法，通过在给定的设计空间内寻找最佳材料分布来实现结构的轻量化和性能优化。
在ANSYS Workbench 中，可以通过拓扑优化模块来进行这种优化分析。
以下是一个简单的ANSYS Workbench 拓扑优化实例的步骤：
1. 创建几何模型：首先，在ANSYS Workbench 中创建一个几何模型，可以是一个零件或者一个装配体。
2. 设定材料属性：为模型中的材料定义材料属性，包括弹性模量、泊松比等。
3. 设定加载和约束条件：定义模型的加载和约束条件，包括受力点、受力大小和方向，以及约束点和约束类型。
4. 添加拓扑优化模块：在Workbench 中选择拓扑优化模块，并将模型导入到拓扑优化模块中。
5. 设定优化目标和约束条件：定义优化的目标，比如最小化结构的重量或者最大化结构的刚度，同时设置一些约束条件，比如最大应力、最大变形等。
6. 设定优化参数：定义拓扑优化的参数，比如拓扑优化的迭代次数、网格分辨率等。
7. 运行优化分析：运行拓扑优化分析，软件会根据设定的优化目标和约束条件，在给定的设计空间内寻找最佳的材料分布。
8. 分析结果：分析优化结果，查看优化后的结构形态和性能指标，根据需要对设计进行进一步的调整和优化。
请注意，以上步骤是一个简化的示例，实际的拓扑优化分析可能涉及到更多的参数设定和分析步骤。
在实际应用中，建议根据具体的工程问题和软件版本进行详细的参数设定和分析。
如何利用ANSYS进行拓扑优化前言就目前而言，利用有限元进行优化主要分成两个阶段：（1）进行拓扑优化，明确零件最佳的外形、刚度、体积，或者合理的固有频率，主要目的是确定优化的方向；（2）进行尺寸优化，主要目的是确定优化后的的零件具体尺寸值，通常是在完成拓扑优化之后，再执行尺寸优化。
在ANSYS中，利用拓扑优化，可以完成以下两个目的：（1）在特定载荷和约束的条件下，确定零件的最佳外形，或者最小的体积（或者质量）；（2）利用拓扑优化，使零件达到需要的固有频率，避免在使用过程中产生共振等不利影响。
本文主要就在ANSYS环境中如何执行拓扑优化进行说明。
1、利用ANSYS进行拓扑优化的过程在ANSYS中，执行优化，通常分为以下6个步骤：定义需要求解的结构问题选择合理的优化单元类型设定优化和非优化的区域定义载荷步或者需要提取的频率对优化过程进行定义和控制计算并查看结果1.1、定义需要求解的结构问题对于结构进行优化分析，定义结构的物理特性必不可少，例如，需要定义结构的杨氏模量、泊松比（其值在0.1～0.4之间）、密度等相关的结构特性方面的信息，以供结构计算能够正常执行下去。
1.2、选择合理的优化单元类型在ANSYS中，不是所有的单元类型都可以执行优化的，必须满足如下的规定：（1）2D平面单元：PLANE82单元和PLANE183单元；（2）3D实体单元：SOLID92单元和SOLID95单元；（3）壳单元：SHELL93单元。
上述单元的特性在帮助文件中有详细的说明，同时对于2D单元，应使用平面应力或者轴对称的单元选项。
1.3、指定优化和非优化的区域在ANSYS中规定，单元类型编号为1的单元，才执行优化计算；否则，就不执行优化计算。
例如，对于结构分析中，对于不能去除的部分区域将单元类型编号设定为≥2，就可以不执行优化计算，请见下面的代码片段：…………Et,1,solid92Et,2,solid92……Type,1Vsel,s,num,,1,2Vmesh,all……Type,2Vsel,s,num,,3Vmesh,all…………说明：上述代码片段定义相同的单元类型（solid92），但编号分别为1和2，并将单元类型编号1利用网格划分分配给了1＃体和2＃体，从而对其进行优化计算；而单元编号为2利用网格划分分配给了3＃体，从而不执行优化计算。
88-51-FB-5B-3A-66现代设计理论和方法大作业----基于MATLAB和ANSYS的悬臂梁拓扑优化指导老师：志刚项目组长：薛亚波项目成员：机自66学院：机械工程学院基于MATLAB和ANSYS的悬臂梁拓扑优化一、计划和任务安排表1.1学习时间安排表：1.2任务分配表：二、项目背景介绍及问题描述2.1项目背景及意义：2.1.1工程背景及基本原理：通常把结构优化按设计变量的类型划分成三个层次：结构尺寸优化、形状优化和拓扑优化。
尺寸优化和形状优化已得到充分的发展，但它们存在着不能变更拓扑结构的缺陷。
在这样的背景下，人们开始研究拓扑优化。
拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域寻求最优材料的分布问题。
寻求一个最佳的拓扑结构形式有两种基本的原理：一种是退化原理，另一种是进化原理。
退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上，然后构造适当的优化模型，通过一定的优化方法逐步删减那些不必要的结构元素，直至最终得到一个最优化的拓扑结构形式。
进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化，它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。
2.1.2 本文研究意义：目前，结构优化大部分集中在尺寸设计变量(如板厚、杆的剖面积及管梁的直径)。
拓扑结构优化较尺寸优化复杂,但对于有些问题拓扑结构优化比尺寸优化有效,悬臂梁是其中的例子之一。
本文讨论悬臂梁的拓扑优化问题,围绕这一问题,怎样使结构具有最大刚度的设计占有相当重要的地位;怎样优化结构的形状使材料的分布，更加合理从而达到使结构具有最大刚度的目的是本文要研究的问题。
2.2研究现状2.2.1 理论研究现状：结构拓扑优化是近20年来从结构优化研究中派生出来的新分支，它在计算结构力学中已经被认为是最富挑战性的一类研究工作。
目前有关结构拓扑优化的工程应用研究还很不成熟，在国外处在发展的初期，尤其在国尚属于起步阶段。
1904年Michell在桁架理论中首次提出了拓扑优化的概念。
自1964年Dorn等人提出基结构法，将数值方法引入拓扑优化领域，拓扑优化研究开始活跃。
20世纪80年代初，程耿东和N.Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题，他们开创性的工作引起了众多学者的研究兴趣。
1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计，开创了连续体结构拓扑优化设计研究的新局面。
1993年XieYM和StevenGP 提出了渐进结构优化法。
1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。
2002年罗鹰等提出三角网格进化法，该方法在优化过程中实现了退化和进化的统一，提高了优化效率。
2.2.2 应用研究现状：在前人提出的重要理论基础上，后人也将其跟其他现代设计的方法相结合，衍生出了其他一些拓扑结构优化方法：如与可靠性相结合的情况下，MAUTE等应用变密度法并结合可靠性分析对一微机电系统进行了基于可靠性的拓扑优化设计，PAPADRAKAKIS等将遗传算法应用于具有可靠性约束的桁架结构拓扑优化设计中，国学者马洪波也对基于遗传算法的结构可靠性优化问题进行了讨论。
华南理工大学机械工程学院欧阳高飞等对基于水平集方法的结构可靠性拓扑优化进行了研究。
2.3研究目标：2.3.1 设计目标：设计一悬臂梁在自由端受到竖直向下集中力P 的作用，要求保持悬臂梁长度尺寸不变。
2.3.2 优化目标：在满足以上条件情况下对原实体悬臂梁做结构拓扑优化设计，使其结构刚度最大、重量最轻。
（优化的结果应该使原设计区域产生孔洞，使结构拓扑发生变化。
）优化目标分析：要使梁满足以上条件情况必须同时兼顾重量和刚度问题。
重量太轻不能满足刚度问题。
刚度太大又必须足够的材料以满足其要求。
将两个因素进行耦合分析，我们可以得出的结论是必须通过结构层次来进行优化。
连续体结构拓扑优化较成熟的方法主要有均匀化方法、变密度方法和渐进结构优化法，通过综合比较我们最终选用变密度法。
2.3.3实现方法：本文在对题目进行分析后，决定从两个独立的方向来分别研究，以此来互相验证，保证结果的可靠性。
实验中要使用ansys 和matlab 两种软件分别独立进行优化分析，然后将两种结果来进行比较。
2.3.4实践目标：通过这次的作业我们期望通过搜索资料，团体讨论，分组作业的方式，以期望达到对拓扑优化的方法和过程有个具体的了解和学习。
另外对分析软件的应用能够达到一个新的高度。
这些不仅能使我们现在的知识体系得到充实和优化，而且也是我们今后人生的财富2.4研究容2.4.1问题描述 如图所示，悬臂梁在自由端受到竖直向下集中力P 的作用，要求在保持悬臂梁长度尺寸不变的情况下对原实体悬臂梁做结构拓扑优化设计，优化目标是使结构刚度最大、重量最轻。
（优化的结果应该使原设计区域产生孔洞，使结构拓扑发生变化。
）原实体悬臂梁横截面为矩形，矩形宽度为0.2m ，悬臂梁的材料为45钢，密度ρ为7.8×10-6kg/mm3，弹性模量E=2×105MPa ，许用切应力[τ]=60MPa 。
竖直向下集中力为N 1045.1P ⨯=。
图悬臂梁受集中力作用（单位：m）2.4.2问题转化本研究中，要求在保持悬臂梁在结构尺寸不变的情况下，使悬臂梁的刚度最大、质量最轻。
一般悬臂梁都为三维立体结构，载荷是在竖直方向上加载的，所以在垂直纸面方向上，并且在悬臂梁的对称平面上加载，根据材料力学知识可知，所加载在和只对加载区周围较小围产生影响，其余部分与加载形式无关。
其受力与如图所示平面的手里相同，故可以将问题抽象为平面问题来处理。
当占有区域S的一个结构上作用有强度为P ( S )的载荷时,使结构具有最大刚度的问题和使结构所受外力功W 具有最小值的问题是等价的。
结构所受的外力功W 与结构变形能C之间的关系为W = 2C，所以我们可以将该问题处理为求在一定的约束条件下，该悬臂梁的的最小变形能最小的问题。
2.4.3数学模型：该问题中，要求同时满足刚度最大，质量最轻，这两个变量若同时改变，则问题复杂度太大，并亲切可能导致问题不可求解。
所以我们采用在确定的质量下，来讨论刚度最大的问题。
由于对特定的材料，其质量和体积有一定的关系，并且我们采用去除法的思想来建立模型的，故我们可以采用给优化后的体积与优化前的体积比赋确定的值，来达到在给定质量条件下满足刚度最大的问题。
其数学模型如下：注：其中C（x）为结构变形能，U为结构变形总位移矩阵，K为结构u为单元位移向量，o k为单元刚度，总刚度矩阵，N为划分单元总数，e（由于划分单元的时候，我们采用等分举行单元，所以个单元的刚度可用V0为未经一个常量来处理）V(x)是拓扑结构优化过程中变化着的体积，过优化前悬臂梁的体积。
F为结构所受的总载荷。
x为悬臂梁的相对密度。
2.4.3模型分析求解：该问题的优化方法有很多种，常用的有如下方法：Optimality Criteria(OC) methods,（优化准则方法）Sequential Linear Programming (SLP) methods（序列线性规划法）Method of Moving Asymptotes (MMA bySvanberg 1987)等为了简化问题的复杂度，此处我们采用standard OC-method.方法来实现。
在处理过程中，关于设计变量相对密度x每一步的更新，我们采用在1995年提出的如下算法来实现：注：其中m是移动限制量，即相对密度变化的最大步长，更新条件中，加上最小边界条件判断，是防止相对密度变为0，使问题出现奇异；加上相对密度最大边界条件判断，是为了防止出现相对密度大于1，是问题是去物理意义。
是数值阻尼系数。
B e可以由以下条件来确定：其中λ为拉格朗日乘数，由双向切片算法来确定。
目标函数对单元相对密度的变化率，即单元敏感度如下确定：为了确保该问题的拓扑结构优化解的存在，在处理技术上，我们得对该问题的求解过程加些限制。
此处我们采用过滤技术来进行处理。
我们采用网格独立性滤波器来对单元敏感度的计算来进行改进，其改进算法如下：H f∧为卷积因子（即权重因子），其由以下表达式来确定：其中采用该算法对单元敏感度的更新进行改进，是为了在更新单元e的敏感度r m in的临域的单元来进行均匀化，防止各乡邻单元的时，将距其距离不超过敏感度变化过大。
注：其中dist(e,f)定义为单元体e的中心到单元体f中心之间的距r m in为过滤大小。
在优化过程中，我们采用改进的敏感度表达式离，（5）来代替表达式（4）来对（3）是进行计算。
三、实验设计及结果分析：3.1技术路线：在本实验中，由于我们提出了两种独立的方法，故我们将小组分为两部分;小组一：该组成员采用Ansys软件来对该悬臂梁模型进行建模分析，并进行拓扑结构优化。
小组二：该组成员采用上面的数学模型，用Matlab原来来实现整个数学模型的算法。
并对处理结果进行动态显示，以便让我们清晰地看到整个拓扑结构优化的过程。
3.2试验结果及分析：3.2.1 Ansys分析的拓扑结构优化结果：3.2.1.2 建立几何模型：3.2.1.2.1 几何参数设置如下：3.2.1.2.1在Ansys/workbench中建立如下几何模型：3.2.1.2网格化：将几何模型直接调入有限元环境。
在workbench里建立的模型可以无缝导入到有限元分析环境里。
采用自动网格划分工具，得到的模型如下：其中，网格单元和节点数设置如下：StatisticsNodes 节点数目3120Elements 单元数561目选择优化类型为“shape optimization”——形状优化在上截面施加15000N的力3.2.1.4求解：StructuralYoung's Modulus 氏模量 2.e+011 PaPoisson's Ratio 0.3Density 密度7850. kg/m³Thermal Expansion 1.2e-005 1/°C Tensile Yield Strength 2.5e+008 PaCompressive Yield Strength 2.5e+008 PaTensile Ultimate Strength 4.6e+008 PaCompressive Ultimate Strength 0. PaThermalThermal Conductivity 60.5 W/m·°C Specific Heat 434. J/kg·°C ElectromagneticsRelative Permeability 10000 Resistivity 1.7e-007 Ohm·mScopeGeometry All BodiesDefinitionTarget Reduction 50. %3.2.1.4.4 其他两种结果：定义质量减少率为80%ScopeGeometry All BodiesDefinitionTarget Reduction 80. %拓扑结果：定义质量减少率为60%ScopeGeometry All BodiesDefinitionTarget Reduction 60. %60%的拓扑结果如下3.2.1.4试验结果分析及结论：从拓扑结构我们可以看到，在不同的材料去除率下，总的趋势都是边界面梁，其中包含着等强度设计的思想，这种拓扑结果较为符合现实情况，在工程上具有很强的实用性。
