ABAQUS拓扑优化手册.docx
《ABAQUS拓扑优化手册.docx》由会员分享,可在线阅读,更多相关《ABAQUS拓扑优化手册.docx(25页珍藏版)》请在冰豆网上搜索。
ABAQUS拓扑优化手册
ABAQUS拓扑优化分析手册/用户手册
分析手册:
13.OptimizationTechnique优化技术
13.1结构优化:
概述
13.1.1概述
ABAQUS结构优化是一个帮助用户精细化设计的迭代模块。
结构优化设计能够使得结构
组件轻量化,并满足刚度和耐久性要求。
ABAQUS提供了两种优化方法一一拓扑优化和形状
优化。
拓扑优化(Topologyoptimization)通过分析过程中不断修改最初模型中指定优化区域的单元材料性质,有效地从分析的模型中移走/增加单元而获得最优的设计目标。
形状优化
(Shapeoptimization)则是在分析中对指定的优化区域不断移动表面节点从而达到减小局部应力集中的优化目标。
拓扑优化和形状优化均遵从一系列优化目标和约束。
最优化方法(Optimization)是一个通过自动化程序增加设计者在经验和直觉从而缩短研发过程的工具。
想要优化模型,必须知道如何去优化,仅仅说要减小应力或者增大特征值
是不够,做优化必须有更专门的描述。
比方说,想要降低在两种不同载荷工况下的最大节点
力,类似的还有,想要最大化前五阶特征值之和。
这种最优化的目标称之为目标函数(Object
Function)。
另外,在优化过程中可以同时强制限定某些状态参量。
例如,可以指定某节点的位移不超过一定的数值。
这些强制性的指定措施叫做约束(Constraint)。
ABAQUS/CAE可以创建模型然后定义、配置和执行结构优化。
更多信息请参考用户手册
第十八章。
13.1.2术语(Terminology)
设计区域(Designarea):
设计区域即模型需要优化的区域。
这个区域可以是整个模型,
也可以是模型的一部分或者数部分。
一定的边界条件、载荷及人为约束下,拓扑优化通过增
加/删除区域中单元的材料达到最优化设计,而形状优化通过移动区域内节点来达到优化的目的。
设计变量(Designvariables):
设计变量即优化设计中需要改变的参数。
拓扑优化中,设计区域中单元密度是设计变量,ABAQUS/CAE优化分析模块在其优化迭代过程中改变单
元密度并将其耦合到刚度矩阵之中。
实际上,拓扑优化将模型中单元移除的方法是将单元的
质量和刚度充分变小从而使其不再参与整体结构响应。
对于形状优化而言,设计变量是指设
计区域内表面节点位移。
优化时,ABAQUS或者将节点位置向外移动或者向内移动,抑或
不移动。
在此过程中,约束会影响表面节点移动的多少及其方向。
优化仅仅直接修改边缘处
的节点,而边缘内侧的节点位移通过边缘处节点插值得到。
设计循环(Designcycle):
优化分析是一种不断更新设计变量的迭代过程,执行ABAQUS进行模型修改、查看结果以及确定是否达到优化目的。
其中每次迭代叫做一个设
计循环。
优化任务(Optimizationtask):
—次优化任务包含优化的定义,比如设计响应、目标、
限制条件和几何约束。
设计响应(Designresponses:
优化分析的输入量称之为设计响应。
设计响应可以直接从ABAQUS的结果输出文件.odb中读取,比如刚度、应力、特征频率及位移等。
或者ABAQUS从结果文件中计算得到模型的设计响应,例如质心、重量、相对位移等。
一个设计响应与模
型紧密相关,然而,设计响应存在一定的范围,例如区域内的最大应力或者模型体积。
另外,
设计响应也与特点的分析步和载荷状况有关。
目标函数(Objectivefunctions):
目标函数决定了优化的目标。
一个目标函数是从设计响应中萃取的一定范围内的值,如最大位移和最大应力。
一个目标函数可以用多个设计响应
来公式表示。
如果设定目标函数最小化或者最大化设计响应,ABAQUS拓扑优化模块则通
过增加每个设计响应值代入目标函数进行计算。
另外,如果有多个目标函数,可以试用权重
因子定义每个目标函数的影响程度。
约束(Constraints:
约束亦是从设计变量中萃取的一定范围的数值。
然而,一个约束不能由设计响应集合而来。
约束限定了设计响应,比如可以指定体积必须降低45%或者某
个区域的位移不能超过1mm。
约束也可以指定制造跟优化无关的制造或者几何约束,比如
轴承面的直径不能改变。
停止条件(Stopconditions):
全局停止条件决定了优化的最大迭代次数。
局部停止条
件在局部最大/最小达成之后指定优化应该停止。
13.1.3ABAQUS/CAE结构优化步骤
下面的步骤需要合并到ABAQUS/CAE模型结构优化设计中:
1)创建需要优化的ABAQUS模型。
2)创建一个优化任务。
3)创建设计响应。
4)禾U用设计响应创建目标函数和约束。
5)创建优化进程,提交分析。
基于优化任务的定义及优化程序,ABAQUS/CAE拓扑优化模块进行迭代运算:
1)准备设计变量(单元密度或者表面节点位置)。
2)更新ABAQUS有限元模型。
3)执行ABAQUS/Standard分析。
这些迭代或者设计循环不会停止,除非:
1)最大迭代数达到
2)指定的停止条件达到。
下图描述了ABAQUS优化分析的过程:
13.1.4拓扑优化
拓扑优化开始于包含指定条件(例如边界条件和载荷)的初始设计开始。
优化分析过程在符合优化约束(比如最小体积或者最大位移)的前提下改变初始设计区域的单元密度和刚
度从而确定结构新的材料分布方式。
下图展示了汽车控制臂在17次设计循环中拓扑优化的过程,其中优化的目标函数是试
图最小化控制臂的最大应变能,最大化控制臂的刚度。
约束为降低57%的产品体积。
优化过
程中,控制臂中部的部分单兀不断被移除。
StartAfter5cyclesAfter10cycles
100%volumeB5%volume77%volume
ABAQUS可以应用如下目标到拓扑优化过程中:
1)应变能(结构刚度的度量值)
2)特征频率
3)内力和支反力
4)重量和体积
5)重心
6)惯性矩。
可以应用其他相同约束变量到拓扑优化分析中。
另外,拓扑优化同样可以考虑标准产品制造过程。
例如铸造和冲压。
可以冻结指定区域、应用数量尺寸、对称性及耦合约束。
拓扑优化的例子在ABAQUSExampleProblemsManual的Sectionll.1.1中。
13.1.5基于密度(一般)VS刚度的拓扑优化
拓扑优化支持两种算法一一一般算法比较灵活,可以应用到大多数问题中;基于刚度的
算法,更为有效,当时应用能力有限。
ABAQUS默认采用一般算法,但是当创建优化任务的
时候可以进行优化算法的选择。
每种优化算法达到优化目标的过程是不同的。
算法
一般拓扑优化算法在满足目标函数和约束前提下使得密度及刚度较好地匹配设计变量,可参考Bends?
eandSigmund(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拓扑优化模块中通过操作来
或者设计响应的范围值。
例如,一个体积设计响应可以只用于设计区域的体积之和,计算vonMises应力的设计响应必须是模型中一定区域的应力最大值。
如下设计响应操作为
ABAQUS/CAE拓扑优化模块(以后简称ATOM)所提供:
最小化或者最大化:
既定区域的最小值或者最大值。
ATOM允许的最大值操作为:
应力,接
触应力和应变设计响应。
求和:
既定区域所有值求和值。
ATOM允许体积、重量、惯性矩和重力设计响应的和值优化。
13.2.1.2基于刚度的拓扑优化设计响应
ATOm提供基于刚度的拓扑优化方法以应变能和体积设计响应。
应变能
结构的顺从性是结构挠性及刚性的量度,被定义为所有单元的应变能之和。
顺从性与刚度成反比,因此最小化顺从性计算最大化全局刚度。
如果载荷是集中力或者面力,则必须选择最小化应变能而不是最大化整体刚度;而如果载荷是热应力场,优化就会在使得结构变软的同时,应变能不断减小,因为降低应变能会导致硬化,所以应该一直选择最大化应变能。
另外,如果模型指定了位移,则须最大化应变能。
拓扑优化考虑所有单元的总应变能,因此,如果选择应变能作为目标函数,则须将目标
应用于整个模型。
优化过程中不能使用应变能作为约束。
Abaqus/CAEUsage:
Optimizationmodule:
Task―stiffnesstopologytask,
DesignResponse—Create:
Single-term,VariableStrainenergy
体积
体积是设计区域内所有单元体积之和。
对于大多数优化问题,需要给定一个ABAQUS模
型。
但是如果最小化应变能(最大化刚度)。
如果优化需要最小化应变能,则不可用体积作
为约束。
Abaqus/CAEUsage:
Optimizationmodule:
Task—^stiffnesstopologytask,
DesignResponseCreate:
Single-term,Variable:
Volume
13.2.1.3基于密度的拓扑优化(一般拓扑优化)设计响应
ATOM一般拓扑优化算法支持重心、平动、旋转、特征频率、惯性矩、内力、反力、力矩、应变能、体积和重量设计响应。
重心
在优化分析中将指定区域的重心作为一个设计变量,并且可以选择三个主分量:
r/w(i\
『厂片L
I,
当ATOM计算重心时,单元会以当前相对密度为计算依据。
举例说明,当需要约束Y方向的重心在某个范围内的时候,这即为一个重心优化问题。
设计响应会考虑整个模型或者指定区域的重心位置。
如果使用了局部坐标系,ATOM会根据坐标轴和原点重新计算重心。
如果不选择使用局
部坐标系,ATOM将使用整体坐标系。
壳体及膜区域的重心在ATOM中考虑其厚度后作为三维区域处理的。
ATOM只能对拓扑优化支持的单元类型计算重心。
因此,ATOM计算得到的重心跟ABAQUS/Stardard和
ABAQUS/Explicit不同,例如,ATOM不会考虑线区域。
Abaqus/CAEUsage:
Optimizationmodule:
Task-^generaltopologytask,
DesignResponseCreate:
Single-term,Variable:
Centerofgravity
平动和转动
在大多数优化问题中需要使用位移或者转动去定义目标函数。
如果仅仅应用位移或者转
动到某些制高点或者小区域,则优化性能会得到改善。
另外,如果指定某些区域为一定位移或者冻结某些区域,同样会加快优化的进程。
下表提供了可用的位移和转动变量:
Displacement
Rotation
i-direction
Absolute
Q处+处+M
/於+哼+卑
Absoluteini-direction
Usage:
Optimizationmodule:
Taskgeneraltopologytask,DesignResponseCreate:
Single-term,Variable:
Displacement
模态特征频率分析
模态特征值是结构分析中最简单的动力响应。
特征频率数据在拓扑优化过程中通常的用
途包括如下几种:
最大化最低阶特征频率,最大化选定的特征频率,限定某特征频率高于或
低于一给定值,最大化或最小化某模态下的特征频率,采用带隙(bandgap)优化迫使模态
远离某频率。
Abaqus/CAE拓扑优化模块支持两种方式来评估特征频率:
模态分析中得到的单一特征频率,以及采用Kreisselmaier-Steinhauser方程(?
不懂)。
Kreisselmaier-Steinhauser方程是两种方法中最有效的,并且可以随时使用。
评估单一特征频率的唯一好处就是可以将特征频率的总和作为约束条件,Kreisselmaier-Steinhauser公式
则不能使用。
当尝试最大化最低阶固有频率时,在此推荐除了考虑第一阶特征频率,还至少要考虑之
后两个高点的固有频率。
在优化过程中,很多固有频率根据它们到最低阶频率的距离进行加权一一越接近第一阶固有频率,权重越大。
如果是尝试最大化最小特征频率或者特别的,最
大化不止一个最小特征频率时,应该使用Kreisselmaier-Steinhauser特征值方程。
因为如果
使用了Kreisselmaier-Steinhauser方程来最大化最小特征频率,则不需追踪最小模态,但是应当追踪高阶的模态,因为模态可能已经发生转换。
举例来说,当模型正在被优化时,第一
阶频率被最大化,可能第二个特征模态变成了最小特征频率的模态。
Abaqus/CAEUsage:
Optimizationmodule:
Taskgeneraltopologytask,DesignResponseCreate:
Single-term,Variable:
EigenfrequencyfrommodalanalysisorEigenfrequencycalculatedwithKreisselmaier-Steinhauserformula惯性矩
在优化过程中可使用惯性矩作为设计响应,以最小化关于一选定轴的转动惯量。
在常规
优化中,可使用整个模型或者某一区域的惯性矩之和作为目标函数或者约束条件。
关于惯性矩可以选择在三个坐标轴上或者三个坐标平面内的惯性矩:
4=I卩+z2)(IV.
打=fp(x2+x2)dl
hy=I帖=-]㈣汕、
F-Tj■—■fxi~~,—'[fhF二(t\4G=!
戸=_J㈣氛八1
如果选择了局部坐标系,Abaqus/CAE拓扑优化模块将采用局部坐标系坐标轴方向重新计算重心。
没有指定局部坐标系时,将缺省采用全局坐标系。
当Abaqus/CAE拓扑优化模块计算惯性矩时,将把板桥和膜区域的厚度赋上并做为三维区域考虑。
并且Abaqus/CAE拓扑优化模块只计算那些支持拓扑优化的单元的惯性矩。
所以,
Abaqus/CAE拓扑优化模块计算出的惯性矩可能与Standard和Explicit的值不同。
例如,模型中包括线域,如果对于两个正交的轴,选择其中任何一个作为对称轴的话,对该轴的惯性矩
为0.
Abaqus/CAEUsage:
Optimizationmodule:
Taskgeneraltopologytask,DesignResponseCreate:
Single-term.Variable:
Momentofinertia
内力和弯矩
在常规拓扑优化中,可使用整个模型或某一区域的节点内力和内弯矩,作为优化目标和
约束条件。
Table13.2.1-2列出了可用的节点内力和内弯矩变量
Table13.2.1—Nodalinternalforceandmomentvariablesforageneraltopologyoptimization.
Nodalinternalforce
Nodalinternalmoment
i-direction
EA;Ui
ESi
Absolute
£Ui\
X乩“I
Absoluteini-direction
y«FiI
节点内力节点内弯矩
i-direction
绝对值
绝对值ini-direction
在此不能使用只有节点力和弯矩绝对值的参考坐标系。
在优化中,结构必须在使用的力
的方向上有刚度,否则,在此方向的内力将变为0.
Abaqus/CAEUsage:
Optimizationmodule:
Taskgeneraltopologytask,Design
ResponseCreate:
Single-term,Variable:
InternalforceorInternalmoment
接触反力和力矩
节点反力和力矩只能用于一般算法的拓扑优化设计。
下表列出了可供优化的设计响应。
Nodalreactionforce
Nodalreactionmoment
i-direction
E凤喇
L瓦e冶
Absolute
Z人凹1
工兀丐1
Absoluteini-direction
应
节点反力和力矩不能使用局部坐标系。
结构在力的方向上必须有刚度,否则反力置为零。
Abaqus/CAEUsage:
Optimizationmodule:
Task―^generaltopologytask,
DesignResponse—Create:
Single-term,Variable:
ReactionforceorReactionmoment
应变能
结构的顺从性是结构整体刚度的度量,在线性模型中被定义为所有单元的应变能。
顺从性与
刚度成反比,最小化顺从性等效于最大化整体刚度。
如果载荷为温度载荷,则结构越软,应
变能越低,最小化应变能则变得非常费劲。
另外,当指定模型位移时,需要一直选择最大化
应变能。
应变能只能用于整个模型。
Abaqus/CAEUsage:
Optimizationmodule:
Task—^generaltopologytask,
DesignResponseCreate:
Single-term,Variable:
Strainenergy
体积
体积即为设计区域内所有单元体积之和。
对于大多数问题,需要采用体积约束。
例如,最小化应变能(最大化刚度)而不采用体积约束,则ATOM会简单的
对设计区域进行材料填充。
Abaqus/CAEUsage:
Optimizationmodule:
Task~Jgeneraltopologytask,
DesignResponseCreate:
Single-term,Variable:
Volume
重量
重量是设计区域内所有单元重量之和。
ATOM十算中会不断采用当前密度。
对于多数优
化问题,如果不采用体积约束则须采用重量约束。
ATOM仅对支持的单元类型进行重量优化。
Abaqus/CAEUsage:
Optimizationmodule:
Task-^generaltopologytask,
DesignResponseCreat