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）: 目标函数决定了优化的目标。

一个目标函数是从设计响应中萃取的一定范围内的值，如最大位移和最大应力。

一个目标函数可以用多个设计响应来公式表示。

如果设定目标函数最小化或者最大化设计响应，ABAQUS拓扑优化模块则通过增加每个设计响应值代入目标函数进行计算。

另外，如果有多个目标函数，可以试用权重因子定义每个目标函数的影响程度。

约束（Constraints）: 约束亦是从设计变量中萃取的一定范围的数值。

然而，一个约束不能由设计响应集合而来。

约束限定了设计响应，比如可以指定体积必须降低45%或者某个区域的位移不能超过1mm。

约束也可以指定制造跟优化无关的制造或者几何约束，比如轴承面的直径不能改变。

停止条件（Stop conditions）: 全局停止条件决定了优化的最大迭代次数。

局部停止条件在局部最大/最小达成之后指定优化应该停止。

13.1.3 ABAQUS/CAE结构优化步骤下面的步骤需要合并到ABAQUS/CAE模型结构优化设计中：1）创建需要优化的ABAQUS模型。

2）创建一个优化任务。

3）创建设计响应。

4）利用设计响应创建目标函数和约束。

5）创建优化进程，提交分析。

基于优化任务的定义及优化程序，ABAQUS/CAE拓扑优化模块进行迭代运算：1）准备设计变量（单元密度或者表面节点位置）。

2）更新ABAQUS有限元模型。

3）执行ABAQUS/Standard分析。

这些迭代或者设计循环不会停止，除非：1）最大迭代数达到2）指定的停止条件达到。

下图描述了ABAQUS优化分析的过程：13.1.4 拓扑优化拓扑优化开始于包含指定条件（例如边界条件和载荷）的初始设计开始。

优化分析过程在符合优化约束（比如最小体积或者最大位移）的前提下改变初始设计区域的单元密度和刚度从而确定结构新的材料分布方式。

下图展示了汽车控制臂在17次设计循环中拓扑优化的过程，其中优化的目标函数是试图最小化控制臂的最大应变能，最大化控制臂的刚度。

约束为降低57%的产品体积。

优化过程中，控制臂中部的部分单元不断被移除。

ABAQUS可以应用如下目标到拓扑优化过程中：1）应变能（结构刚度的度量值）2）特征频率3）内力和支反力4）重量和体积5）重心6）惯性矩。

可以应用其他相同约束变量到拓扑优化分析中。

另外，拓扑优化同样可以考虑标准产品制造过程。

例如铸造和冲压。

可以冻结指定区域、应用数量尺寸、对称性及耦合约束。

拓扑优化的例子在ABAQUS Example Problems Manual的Section11.1.1中。

13.1.5 基于密度（一般）VS刚度的拓扑优化拓扑优化支持两种算法——一般算法比较灵活，可以应用到大多数问题中；基于刚度的算法，更为有效，当时应用能力有限。

ABAQUS默认采用一般算法，但是当创建优化任务的时候可以进行优化算法的选择。

每种优化算法达到优化目标的过程是不同的。

算法一般拓扑优化算法在满足目标函数和约束前提下使得密度及刚度较好地匹配设计变量，可参考Bendsøe and Sigmund (2003)，文献中有一般算法的部分描述。

相反的，刚度拓扑优化作为一种更为有效的算法，使用了应变能和节点力作为输入量而且并不需要计算设计变量的局部刚度。

基于刚度的优化算法是德国卡尔斯鲁厄大学（U.Karlsruhe）的Bakhtiary (1996)提出的。

中间密度单元一般算法导致了最终设计中间单元的产生（相对密度处于0~1之间的单元）。

相反地，基于刚度的优化算法会产生不是空材料即是实体（0或者1）。

优化设计循环步优化设计需要的循环步在优化计算开始的时候是并不知晓的，但是一般来说，这个步数处于30~45之间。

基于刚度的优化算法能够更快地达到优化算法的解（默认15步）。

分析类型一般优化算法支持线性静力、非线性静力及线性模态分析。

两种算法均支持几何非线性和接触，以及很多非线性材料。

此外，静力拓扑优化中可以指定位移。

然而，模态分析不可指定位移。

拓扑优化支持符合材料分析，个别复合材料的碾压分析是不支持的，如不能指定复合材料纤维走向等。

目标函数和约束一般优化算法可以使用一个目标函数和数个约束，这些约束可以全部是不等式。

多种设计响应可以被定义为目标和约束，例如应变能、位移和旋转、支反力及内力、特征频率和材料的体积及重量。

基于刚度的拓扑优化算法更为有效，然而适用性比较差，仅支持应变能（一种刚度的两度）作为目标函数，材料体积作为约束方程。

13.1.6 形状优化形状优化采用了跟基于刚度的拓扑优化算法类似的算法。

形状优化一般是对表面节点进行较小的调整以减小局部应力集中。

形状优化用于产品外形需要微调的情况。

形状优化试图重置既定区域的表面节点位置直到此区域的应力成为常数（应力均匀）。

下图是连杆形状优化以减小局部应力集中的例子：形状优化支持一下目标：1）应力和接触应力2）自然频率3）弹性、塑性、全应变和应变能密度形状优化只能应用体积约束，另外，可以使用一定数量的制造几何限制条件使提出的设计能够继续铸造或者冲压过程。

也可以冻结某特定区域、应用数量尺寸、对称性及耦合限制等。

形状优化的网格光顺形状优化过程中，ABAQUS/CAE拓扑优化模块修改模型表面。

如果拓扑优化模块，只对表面节点进行位置调整而不对表层内节点进行调整，单元将会发生扭曲现象。

这样，ABAQUS 的优化结果将变得不可信。

为了获得相当质量的表面单元，ABAQUS/CAE拓扑优化模块可以对既定区域进行网格光顺，从而使得内外节点位置关系更加和谐。

需要注意的是，在进行优化之前，指定区域必须具有较好的网格质量。

ABAQUS/CAE拓扑优化模块可以光顺标准连续单元，比如三角形单元、四边形单元和四面体单元。

其他单元类型将不会被光顺化。

可以指定光顺单元的倾斜角度或者单元质量。

网格光顺化是一个比较耗费计算量的过程。

光顺算法是基于单元的，计算时间也会随着区域内单元密度及自由度情况响应增加。

因此，可以只对优化区域内的单元指定网格光顺化。

网格光顺化的区域节点必须是自由的，不能应用于有约束的节点或者区域。

网格光顺化的设定具体参见ABAQUS/CAE用户手册。

网格光顺化可以应用到优化区域也可以应用到非优化区域，然而优化区域必须包含于网格光顺化指定区域。

ABAQUS/CAE在网格光顺化中会自动将所有表面单元进行约束，从而使优化结果得到保存。

然而可以选择设计区域附近的表层节点进行光顺化移动，从而使得优化区域和非优化设计区域能够更加连续。

默认的网格光顺化采用强制拉普拉斯网格光顺化算法。

但是，如果模型比较小，比如小于1000个节点，可以采用局部梯度化网格光顺算法。

在每次运算中，局部梯度化网格光顺算法进行质量不好的单元的辨别，然后通过移动节点对其进行光顺。

局部梯度化光顺通常产生具有最佳外形的单元。

对于较大模型，局部梯度化网格光顺算法一般会在最佳网格质量形成之前停止，这样只有最坏的单元质量得到了光顺。

13.2 优化模型13.2.1设计响应13.2.1.1概述一个设计响应：1)是一个单范围值，例如结构的体积2)通过ABAQUS/CAE读取odb文件的结果和模型文件来计算3)与目标函数和约束相关联（比如，可以创建一个目标函数和约束去最小化节点位移或者强制降低结构重量50%以上）4)必须在一定的分析过程基础上（比如，在对既定区域选择最大化最小特征频率则需要进行提取特征频率的分析）尽管已有一些限制应用到了响应上，还必须在ABAQUS/CAE拓扑优化模块中通过操作来或者设计响应的范围值。