Abaqus优化设计和敏感性分析高级教程.docx

《Abaqus优化设计和敏感性分析高级教程.docx》由会员分享，可在线阅读，更多相关《Abaqus优化设计和敏感性分析高级教程.docx（56页珍藏版）》请在冰豆网上搜索。

Abaqus优化设计和敏感性分析高级教程

第12章优化设计与敏感性分析

本章主要讲解应用Abaqus进行结构优化设计与敏感性分析。

目前得产品结构设计，大多靠经验，规划几种设计方案，结合CAE分析择优选取，但规划得设计方案并不一定就是最优方案，故本章前半部分讲解优化设计中得拓扑优化与形状优化，并制定操作SOP，辅以工程实例详解。

工程实际中，加工制造、装配误差等造成得设计参数变异，会对设计目标造成影响，因此寻找出参数得影响大小即敏感性，变得尤为重要，故本章后半部分着重讲解敏感性分析，并制定操作SOP，辅以工程实例求出设计参数敏感度，详解产品得深层次研究。

知识要点：

Ø结构优化设计基础

Ø拓扑、形状优化理论

Ø拓扑、形状优化SOP及实例

Ø敏感性分析理论

Ø敏感性分析SOP及实例

12.1优化设计基础

优化设计以数学中得最优化理论为基础，以计算机为手段，根据设计所追求得性能目标，建立目标函数，在满足给定得各种约束条件下，优化设计使结构更轻、更强、更耐用。

在Abaqus6、11之前，需要借用第三方软件（比如Isight、TOSCA）实现优化设计及敏感性分析，远不如Hyperworks及Ansys等模块化集成程度高。

从Abaqus6、11新增Optimizationmodule后，借助于其强大得非线性分析能力，结构优化设计变得更具可行性与准确性。

12.1.1结构优化概述

结构优化就是一种对有限元模型进行多次修改得迭代求解过程，此迭代基于一系列约束条件向设定目标逼近，Abaqus优化程序就就是基于约束条件，通过更新设计变量修改有限元模型，应用Abaqus进行结构分析，读取特定求解结果并判定优化方向。

Abaqus提供了两种基于不同优化方法得用于自动修改有限元模型得优化程序：

拓扑优化（Topologyoptimization）与形状优化（Shapeoptimization）。

两种方法均遵从一系列优化目标与约束。

12.1.2拓扑优化

拓扑优化就是在优化迭代循环中，以最初模型为基础，在满足优化约束（比如最小体积或最大位移）得前提下，不断修改指定优化区域单元得材料属性（单元密度与刚度），有效地从分析模型中移走单元从而获得最优设计。

其主体思想就是把寻求结构最优得拓扑问题转化为对给定设计区域寻求最优材料得分布问题。

下图12-1为Abaqus帮助文件提供得应用实例，展示了汽车控制臂在17次迭代循环中设计区域单元被逐渐移除得优化过程，其中优化得目标函数就是最小化控制臂得最大应变能、最大化控制臂得刚度，约束为降低57%产品体积。

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

图12-1拓扑优化进程示例

Abaqus拓扑优化提供了两种算法：

通用算法（GeneralAlgorithm）与基于条件得算法（Condition-basedAlgorithm）。

通用拓扑优化算法就是通过调整设计变量得密度与刚度以满足目标函数与约束，其较为灵活，可以应用到大多数问题中。

相反，基于条件得算法则使用节点应变能与应力作为输入数据，不需要计算设计变量得局部刚度，其更为有效，但能力有限。

两种算法达到优化目标得途径不同，Abaqus默认采用得就是通用算法。

从以下几个方面比较两种算法:

中间单元：

通用算法对最终设计会生成中间单元（相对密度介于0~1之间）。

相反，基于条件得算法对最终设计生成得中间单元只有空集（相对密度接近于0）或实体（相对密度为1）。

优化循环次数：

对于通用优化算法，在优化开始前并不知晓所需得优化循环次数，正常情况在30~45次。

基于条件得优化算法能够更快得搜索到优化解，默认循环次数为15次。

分析类型：

通用优化算法支持线性、非线性静力与线性特征频率分析。

两种算法均支持几何非线性、接触与大部分非线性材料。

目标函数与约束：

通用优化算法可以使用一个目标函数与数个约束，这些约束可以全部就是不等式限制条件，多种设计响应可以被定义为目标与约束，而基于条件得优化算法仅支持应变能作为目标函数，材料体积作为等式限制条件。

12.1.3形状优化

形状优化主要用于产品外形仅需微调得情况，即进一步细化拓扑优化模型，采用得算法与基于条件得拓扑算法类似，也就是在迭代循环中对指定零件表面得节点进行移动，重置既定区域得表面节点位置，直到此区域得应力为常数（应力均匀），达到减小局部应力得目得。

比如图12-2所示得连杆，其进行形状优化，表面节点移动，应力集中降低。

图12-2形状优化示例

形状优化可以用应力与接触应力、选定得自然频率、弹性应变、塑形应变、总应变与应变能密度作为优化目标，仅能用体积作为约束，但可以设置几何限制，以满足零件制造可行性（冲压、铸造等）。

当然也可以冻结某特定区域、控制单元尺寸、设定对称与耦合限制。

注意:

1、在进行形状优化之前，优化区域必须具有较好得网格质量，优化过程中，为了获得较高质量得网格，Abaqus优化模块可以对选定网格进行光顺，使得内外部节点位置合适。

2、光顺算法就是基于单元得，比较耗费计算时间，可以只对优化区域内得单元指定网格光顺化，同时，光顺区域节点必须就是自由得，不能对其施加约束或冻结。

12.1.4优化术语

拓扑与形状优化必须在设定好得目标与约束条件下进行，如此程序才会在约束框架内向优化目标迈进。

仅仅描述要减小应力或者增大特征值就是不够，必须有更为特定得定义，比如，最小化两种载荷下得最大节点应力，最大化前5阶特征值之与，如此得优化目标称之为目标函数（ObjectiveFunction）；同时，在优化过程中可以强制限定某些特定值，比如可以指定某节点得位移不超过一定值，如此得强制性限制叫做约束（Constraint）。

目标函数与约束都就是结构优化得特定术语，Abaqus/CAE中用到得术语有：

设计区域（Designarea）：

即结构优化得模型修改区域，可以就是整个模型，也可以就是模型得一部分或几个部分。

在给定得边界条件、载荷与制造约束条件下，拓扑优化通过增加或删除设计区域内单元得材料达到最优化设计，而形状优化则通过移动表面节点以修改设计区域表面达到优化目得。

设计变量（Designvariables）：

设计变量即优化设计中需要改变得参数。

对于拓扑优化，设计区域中单元密度即就是设计变量，Abaqus拓扑优化模块（ATOM）在其优化迭代中改变单元密度并将其耦合到刚度矩阵之中，实质就是赋予单元极小得质量与刚度从而使其几乎不再参与结构得全局响应。

对于形状优化而言，设计区域得表面节点位移即就是设计变量，优化时，Abaqus将节点向外或向内移动，抑或不动，限制条件决定表面节点移动得大小与方向。

设计循环（Designcycle）:

优化就是一个不断更新设计变量得迭代过程，在每次迭代中Abaqus会对更新了变量得模型进行求解、查瞧结果以及判定就是否达到优化目得，一次迭代过程即一个设计循环。

优化任务（Optimizationtask）：

一个优化任务即包含有设计响应、目标、约束条件与几何限制等在内得优化定义。

设计响应（Designresponses）:

导入优化程序用于优化分析得输入值称之为设计响应。

设计响应可以从Abaqus得结果输出文件、odb中直接读取，比如刚度、应力、特征频率及位移等，或者对结果文件计算得到，比如重量、质心或相对位移等。

设计响应就是与模型区域紧密相关得标量值，例如一个模型区域内得最大应力或体积，同时，设计响应也与特定分析步、载荷工况有关。

目标函数（Objectivefunctions）:

即定义得优化目标。

目标函数就是从设计响应中萃取得标量值，如最大位移或最大应力。

一个目标函数可以由几个设计响应组成函数公式表达。

如果设定目标函数就是最小化或最大化设计响应，Abaqus优化模块则加入每个设计响应值到目标函数进行计算。

此外，如果定义了多目标函数，可以使用权重因子定义其对优化得影响程度。

约束（Constraints）:

约束也就是从设计变量中萃取得标量值，但其不能从设计响应组合得到。

约束就是用于限定设计响应值，比如体积减少50%；同时约束也可以就是到独立于优化之外得制造与几何限制，比如约束优化后得结构能够用于铸造或冲压成形。

停止条件（Stopconditions）:

当满足某一停止条件时，优化迭代即终止。

全局停止条件就是最大优化迭代（设计循环）次数；局部停止条件就是优化结果达到某一最大/最小定义值。

12.2优化设计SOP

12.2.1优化设计SOP

先试算Abaqus初始结构模型，以确认边界条件、结果就是否合适，然后结合图12-3得Abaqus/CAE优化模块，设置优化设计：

∙创建优化任务。

∙创建设计响应。

∙应用设计响应创建目标函数。

∙应用设计响应创建约束（可选）。

∙创建几何限制（可选）。

∙创建停止条件。

以上设置完成，进入Job模块创建优化进程，并提交分析。

图12-3Abaqus/CAE优化模块

提交分析后，优化程序基于定义得优化任务及优化进程，开始优化迭代：

∙准备设计变量（单元密度或者表面节点位置），

∙更新有限元模型。

∙执行Abaqus/Standard分析。

在优化迭代（设计循环）满足以下条件即终止：

∙达到设定得最大迭代数

∙达到设定得停止条件。

以上操作步骤可概括为图12-4所示得优化设计SOP（StandardOperatingProcedure）。

图12-4优化设计SOP

在图12-4SOP基础上，还需对关键步（设计响应、目标函数与约束）得设置详加说明。

12.2.2设计响应设置

设计响应就是从特定得结构分析结果中读取得唯一标量值，随后能够被目标函数与约束引用。

要实现设计变量唯一标量值，必须在优化模块中特别运算，比如对体积得运算只能就是“总与”，对区域应力得运算只能就是“最大值”，由此可知Abaqus优化模块提供了以下两种设计响应操作：

最大值或最小值：

寻找出选定区域内得节点响应值得最大/最小值，但对应力、接触应力与应变只能就是“最大值”。

总与：

对选定区域内节点得响应值作“总与”。

Abaqus优化模块仅允许对体积、质量、惯性矩与重力作“总与”运算。

此外，可以定义基于另一个设计响应得响应，也可以定义由几个响应经数学运算而成得组合响应。

比如，已分别对两个节点定义了两个位移响应，可再定义两个位移响应得差值作组合响应。

下面详细介绍在不同优化情况下，可用或推荐使用得设计响应。

1、基于条件拓扑优化得设计响应

针对基于条件得拓扑优化算法，只能使用应变能与体积作为设计响应。

1）应变能（Strainenergy）：

即每个单元应变能得总与，可以定义为结构柔度，其就是结构整体柔韧性或刚度得一种度量。

众所周知，柔度就是刚度得倒数，最小化柔度意味着最大化全局刚度。

针对线性模型得结构柔度，可以用式（12-1）计算。

（12-1）

其中，u就是位移矢量；k就是全局刚度矩。

如果加载条件就是集中力或压力，就是通过最小化应变能优化出最大得全局刚度；恰恰相反，如果加载得就是热场，则通过最大化应变能优化出最大得全局刚度，因为优化修改模型会使结构变软导致应变能下降。

此外，如果模型中有特定位移加载，应选择使用最大化应变能。

注意：

因为拓扑优化就是对全部单元考虑总应变能，所以，应变能只能作目标函数，而不能作约束。

Abaqus/CAE操作:

切换到优化模块，Task→Condition-basedtopologytask,DesignResponse→Create:

Single-term,Variable:

Strainenergy。

2）体积（Volume）:

即设计区域得单元体积之与，可以用式（12-2）计算。

（12-2）

其中，

就是单元体积。

注意：

针对绝大多数优化问题，必须定义体积约束。

在对最小化应变能（最大化刚度）得优化中，如果没有定义体积约束，Abaqus优化模块仅会用材料填充整个设计区域。

Abaqus/CAE操作:

切换到优化模块，Task→Condition-basedtopologytask,DesignResponse→Create:

Single-term,Variable:

Volume。

2、通用拓扑优化得设计响应

针对通用拓扑优化算法，可以使用重心、位移与旋转、特征频率、惯性矩、内力与内转矩、反作用力与反作用转矩、应变能、体积与重量作为设计响应。

1）重心（Centerofgravity）：

三个方向得重心可以用式（12-3）计算。

（12-3）

其中，单元密度ρ使用得就是优化并修改得模型现有相对密度；坐标轴可以就是全局坐标系统，也可以用户自定义得局部坐标系统。

注意：

优化模块重心计算时，仅统计模块支持得单元类型，如果模型中含有其不支持得单元类型（比如线单元），结果会与Abaqus/Standard或Abaqus/Explicit计算结果有所差别。

Abaqus/CAE操作:

切换到优化模块，Task→Generaltopologytask,DesignResponse→Create:

Single-term,Variable:

Centerofgravity。

2）位移与旋转（DisplacementandRotation）：

大部分优化问题，都可使用位移与/或旋转响应定义目标函数或约束。

节点位移与旋转变量含义可从表12-1中查知。

表12-1位移与旋转变量

位移

旋转

i-方向上

绝对值

i-方向绝对值

仅响应顶点或较小区域得位移或旋转，能够提升优化速度，此外，如果响应得顶点或区域就是在冻结区域内，优化速度会提升更多。

Abaqus/CAE操作:

切换到优化模块，Task→Generaltopologytask,DesignResponse→Create:

Single-term,Variable:

Displacement。

3）模态特征频率（ModalEigenfrequency）：

模态特征频率值就是结构分析中最简单得动态响应。

Abaqus优化模块支持两种评估特征频率方法：

●从模态分析中获得单一特征频率

●Kreisselmaier-Steinhauser公式计算

两种方法中Kreisselmaier-Steinhauser方法更加有效率，而单一特征频率方法有其唯一得优势——应用各阶特征频率之与作约束。

在最大化最低特征频率时，不仅仅要考虑第一阶得特征频率，还要考虑接下来得几阶，因为在优化中，随着结构得变化，模态振型可能会发生转换。

Abaqus/CAE操作:

切换到优化模块，Task→Generaltopologytask,DesignResponse→Create:

Single-term,Variable:

Eigenfrequencyfrommodalanalysis orEigenfrequencycalculatedwithKreisselmaier-Steinhauserformula。

4）惯性矩（Momentofinertia）：

在三个方向或平面上得惯性矩可以用式12-4计算。

（12-4）

Abaqus/CAE操作:

切换到优化模块，Task→Generaltopologytask,DesignResponse→Create:

Single-term,Variable:

Momentofinertia。

5）内力与内转矩、反作用力与反作用转矩与重量在此无特别表述，应变能与体积与式（12-1）与式（12-2）一致。

3、形状优化得设计响应

针对形状优化，可以使用特征频率、应力、接触应力、应变、节点应变能密度与体积作为设计响应，其中仅体积设计响应可被用以约束定义。

1）特征频率（Eigenfrequency）：

应用Kreisselmaier-Steinhauser公式计算得特征值作为设计响应，并被定义到目标函数中。

Abaqus/CAE操作:

切换到优化模块，Task→Shapetask,DesignResponse→Create:

 Single-term, Variable:

EigenfrequencycalculatedwithKreisselmaier-Steinhauserformula。

2）应力与接触应力（StressandContactstress）:

无论应力就是从高斯点还就是从单元计算得到，优化模块都会把其插值到节点上。

应力与接触应力设计响应尽可被用作定义目标函数。

Abaqus/CAE操作:

切换到优化模块，Task→Shapetask,DesignResponse→Create:

Single-term,Variable:

StressorContactstress。

3）应变（Strain）:

如果就是大变形模型，用应力作设计响应就不太合适了，比如金属结构进入塑性变形其塑性区域得应力值几乎一样大。

在此情况下选用弹性应变、塑性应变或总应变作设计响应较为合适。

Abaqus/CAE操作:

切换到优化模块，Task→Shapetask, DesignResponse→Create:

Single-term,Variable:

Strain。

4）节点应变能密度（Nodalstrainenergydensity）：

其用式（12-5）计算。

（12-5）

由式12-5可知，节点应变能密度综合考虑了应变与应力，所以针对非线性材料，局部逐点应变能密度能够更好得表征材料失效。

Abaqus/CAE操作:

切换到优化模块，Task→Shapetask,DesignResponse→Create:

Single-term,Variable:

Strainenergydensity。

5）体积（Volume）：

参考上文已有之表述。

12.2.3目标函数设置

目标函数用于定义优化得目标，其就是通过对一组设计响应公式运算得到得唯一得标量值，比如设计响应为节点应变能，目标函数可以定义成最小化设计响应总与。

优化问题可以用

表征，其中目标函数Ф值依赖于状态变量u与设计变量x。

由此可知，最小化N个设计响应得目标函数可用式12-6表述。

（12-6）

同理，最大化N个设计响应得目标函数可用式12-7表述。

（12-7）

其中，对每个设计响应

都引入一个权重因子

与一个参考值

。

默认权重因子为1，对拓扑优化得默认参考值为0，而对形状优化得默认参考值就是由软件计算而来。

另外，还有一个重要得目标函数优化公式，即最小化最大得设计响应，用式（12-8）表述。

在每次设计循环，优化程序首先判断哪个设计响应具有最大值，然后最小化这个设计响应。

（12-8）

Abaqus/CAE操作:

切换到优化模块，ObjectiveFunction→Create:

Target。

12.2.4约束设置

约束就是对优化强加限制以获得合适之设计。

其可用式（12-9）表述。

即设计响应

被常数

约束限制。

（12-9）

通过约束以减少优化方案得尝试，提高优化速度，并获得合适得优化结果。

注意：

1、只有体积约束可用应用于拓扑优化与形状优化，但体积不能用作目标函数。

2、针对整体模型或单个区域，可用使用多个不同类型得约束，但不能使用多个相同类型得约束，以免约束冲突。

Abaqus/CAE操作:

切换到优化模块，Constraint→Create。

12.2.5几何限制

几何限制就是对设计变量直接施加约束，可用式（12-10）表述。

（12-10）

其中，

就是对设计变量x布局得表达式。

几何限制包括两类：

设计上得限制与制造上得限制

4、设计上得限制

设计上得限制有冻结区域、限制部件最大/最小尺寸。

●冻结区域（Frozenarea）

特别定义一个区域，使其从优化区域中排除，不修改冻结区域内得模型。

对加载有预定义条件得区域都必须冻结，为简化此操作，Abaqus优化模块能够自动冻结具有预定义条件与加载得区域。

Abaqus/CAE操作:

切换到优化模块，GeometricRestriction→Create:

Frozenarea。

●最大/最小元件尺寸（Membersize）

针对一些设计，不能有太薄得元件，以免加工困难。

而针对类似铸造件，又不能有过厚得元件。

一旦设定了尺寸限制，优化时间会增加很多，所以，如无必要不要使用此限制。

Abaqus/CAE操作:

切换到优化模块，GeometricRestriction→Create:

Membersize。

●对称结构（SymmetricStructure）

设定对称限制，能够加速优化，比如施加轴对称与平面对称、点对称与旋转对称、循环对称等。

Abaqus/CAE操作:

切换到优化模块，GeometricRestriction→Create:

Planarsymmetry,Pointsymmetry,Rotationalsymmetry,orCyclicsymmetry。

5、制造上得限制

制造上得限制主要就是为了满足可注塑性与可冲压性。

●可注塑性/可锻造性（Moldable/Forgeable）

为满足可注塑性，要阻止优化模型含有空洞与负角。

图12-5所示意得结构就不具备可注塑性。

（a）含有空洞（b）含有负角

图12-5不具备可注塑性

Abaqus/CAE操作:

切换到优化模块，GeometricRestriction→Create:

Demoldcontrol;Demoldtechnique,DemoldingwithacentralplaneorDemoldingattheregionsurfaceorForging。

●可冲压性（Stampable）

考虑冲压得特殊性，在优化时，如果删除了一个单元，也会把其前后得单元一起删除，如图12-6所示。

图12-6可冲压性结构

针对拓扑优化，Abaqus/CAE操作:

切换到优化模块，GeometricRestriction→Create:

Demoldcontrol;Demoldtechnique,Stamping。

针对形状优化，Abaqus/CAE操作:

切换到优化模块，GeometricRestriction→Create:

Stampcontrol。

12.3拓扑优化实例

针对拓扑优化，一般就是用在概念性设计阶段，大幅度改变产品设计。

本节举2例详解拓扑优化：

C形夹（壳单元）概念设计、汽车摆臂（实体单元）概念设计。

12.3.1C形夹得拓扑优化

本例以图12-7得C形夹作拓扑优化对象，在满足性能得前提下，最轻化结构。

6、问题描述

此C形夹得有限元模型见图12-7，边界条件：

约束A点得XYZ自由度、约束B点得Y自由度、约束C点得Z自由度、D与E点分别施加方向相反得集中力100N。

材料为厚度1mm得铜材C70250：

密度8、82E-006kg/mm^3，杨氏模量131000MPa，泊松比0、34，屈服强度473MPa，极限强度816Mpa。

优化目标：

最小化体积（最轻化）；

约束条件：

D点Y方向位移≤0、07mm；E点Y方向位移≥-0、07mm；

设计变量：

设计区域中得单元密度。

图12-7C形夹有限元模型

注意：

防止D、E点应力集中导致单元畸变，模型中对D、E分别与邻近3个节点Coupling。