如何利用ANSYS进行拓扑优化(转)

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

2、彳1利用ANSYS进行拓扑优化的过程在ANSYS中，执行优化，通常分为以下 6个步骤：1.1、定义需要求解的结构问题要定义对于结构进行优化分析，定义结构的物理特性必不可少，例如, 结构的杨氏模量、泊松比（其值在 0.10.4之间）、密度等相关的结构特性方面时需Sr彳 的信息，以供结构计算能够正常执行下去。1.2、选择合理的优化单元类型在ANSYS中，不是所有的单元类型都可以执行优化的，必须满足如下的规疋：(1) 2D平面单元：PLANE82单元和PLANE183单元；(2) 3D实体单元：SOLID92单元和SOLID95单元；(3) 壳单元：SHELL93单元。上述单元的特性在帮助文件中有详

3、细的说明，同时对于2D单元，应使用平面应力或者轴对称的单元选项。1.3、指定优化和非优化的区域在ANSYS中规定，单元类型编号为1的单元，才执行优化计算；否则，就 不执行优化计算。例如，对于结构分析中，对于不能去除的部分区域将单元类型 编号设定为2,就可以不执行优化计算，请见下面的代码片段：Et,1,solid92Et,2,solid92Type,1Vsel,s, num,1,2时磊忖呎Vmesh,allType,2Vsel,s ,nu m,3Vmesh,all说明：上述代码片段定义相同的单元类型（solid92），但编号分别为1和2, 并将单元类型编号1利用网格划分分配给了 1#体和2#体，

4、从而对其进行优化 计算；而单元编号为2利用网格划分分配给了 3#体，从而不执行优化计算。1.4、定义载荷步或者需要提取的频率1.4.1线性结构静态分析对于结构优化而言，其总是在特定的载荷（或者载荷步），约束和目标下进 行的，在优化分析的过程中，必须执行线性结构静态分析，才能获得需要的优化 之后的形状。在ANSYS中，可以对单步载荷或者多步载荷执行优化分析，当然, 单步载荷是最简单的了。然而，对于某个特定载荷步，必须使用 LSWRITE载荷 步存储命令将载荷步预存起来，再用 LSSOLVE命令进行求解。先看看下面的代码片段：D,10,all,0,20,1Nsel,s,loc,y,0Sf,Alls

5、elLswrite,1Ddel,Sfdel,Nsel,s,loc,x,0,1D,all,all,0F,212,fx,Lswrite,2Lswrite,3Fi nishTocomp,mcomp,multiple,3Tovar,mcomp,objTovar,volume,c on ,10TodefToloop,20说明：该代码片段首先定义了 3个载荷步，并利用LSWRITE命令将载荷步 预存；之后利用Tocomp命令定义优化任务目标名称 mcomp,并将体积减少10% 作为优化的约束条件，之后用 Todef初始化优化过程，最后利用 Toloop命令执 行优化计算，最大计算次数20次。相关命令：TO

6、COMP、TOVAR、TODEF、TOEXE、TOLOOP 和简要说明。(1) TOCOMP :定义结构优化任务目标。(如何理解COMPLIANCE : Complianee本意是一致性，统一性，在结构优化分析中，特别是对于多个载荷步，需要在多个载荷步之间取得一致性的结果，才能满足结构优化分析的目标)(2) TOVAR :定义优化变量，可以是目标变量，也可以是约束变量等；TODEF :定义优化的初始化条件或者收敛准则;(4)TOEXE :执行单次优化计算；(5)TOLOOP :批量执行多次优化分析计算。说说TOEXE和TOLOOP之间的区别：TOEXE执行单次优化分析计算， 其本身不执行结构分

7、析过程，因此，在利用 TOEXE命令执行优化计算之前，需 要利用SOLVE或者LSSOLVE命令先执行结构静态分析计算；而TOLOOP是- 个执行优化计算的宏命令，其中包含了 SOLVE和LSSOLVE等命令，因此在上 述代码片段中没有出现 SOLVE或者LSSOLVE命令。就使用的便利性而言，利用TOLOOP命令可能更方便，但是利用 TOEXE命令用户可以创建自己的优化 时需Sr彳宏命令，各有所长，主要是看用户如何使用这两个命令了。1.4.2模态优化分析在执行模态优化分析之前，需要使用 TOFREQ和TOVAR定义模态分析任 务名称和优化分析变量，这点和上述的结构优化分析过程类似。值得注意的

8、是， 在利用MXPAND指定模态频率分析计算中，需要将单元计算设定为“YES”这 样才能利用TOEXE命令正常执行优化计算过程，当然在 MXPAND之前还要使 用modopt命令设定模态计算的选项；而利用 TOLOOP命令执行模态优化分析 计算过程中，由于TOLOOP是编制好的宏命令，不需要使用MXPAND等命令， 请见下面的代码片段：/soluti onAn type,modalFi nishTofreq,mfreq,reciprocal,3Tovar,mfreq,objTovar,volume,c on ,50Todef,1.0d-4Toloop,20说明：首先定义模态求解作为分析类型，之

9、后利用tofreq设定频率优化作为优化任务，名称为mfreq (其中reciprocal表示多阶模态频率，本例中为前3阶 模态)，并利用命令tovar设定了体积减少50%作为优化的约束条件，之后设定 优化的收敛准则，最后利用toloop命令执行最多20次优化迭代计算。在本例中， 因为toloop是一个宏命令，所以并没有显式的使用 modopt和mxpand命令，如 果使用toexe命令，则用户必须显式的使用 modopt和mxpand命令。1.5、对优化过程进行定义和控制在ANSYS中执行优化过程有2种方式：(1) 采用solve和toexe命令相结合的机制，一步一步的执行优化计算，直到满足用

10、户的优化目标和约束条件为止。(2) 采用toloop宏命令执行优化计算，可执行多次优化分析计算。说明：采用第二种方式操作更简单，推荐采用。总体上来说，执行优化计算过程大致上分为4个步骤：(1) 定义优化函数(或者方式和任务)；(2) 定义优化目标或者约束条件；(3) 初始化优化过程；(4) 执行优化计算(可以是单次计算，也可以自动批量计算)。以下分别进行描述：1.5.1、定义优化函数在ANSYS中，主要采用以下两种方式定义优化函数:时磊忖呎（1）利用TOCOMP定义结构优化；（2）利用TOFREQ定义频率优化。说明：可以利用TOLIST命令例举出所有定义了的优化函数。1.5.2、定义优化目标或

11、者约束条件定义优化目标或者约束条件主要是利用TOVAR命令设定，看看下面的代码片段：Tocomp，comp，single 1Tovar， comp，objTovar，Volume，con,，25Todef， 1.0d-4Toloop， 10， 1说明：首先利用tocomp命令定义了一个单步（single）的结果优化分析任务， 并将任务名称（comp）作为优化任务，并设定体积减少25%作为优化约束条件， 之后初始化优化过程，并利用toloop命令批量执行优化计算，当然，可以使用tostat 命令查询tovar、todef和totype命令的设定状态（totype命令设定优化类型，具 体说明可以参

12、见帮助文件）。同样，也可以利用如下的命令格式删除已经定义好 的优化任务名称或者约束条件:Tovar, Reference_Name DEL可以定义的优化目标和约束条件的组合方式见下表:序号优化类型优化目标约束条件1结构优化(TOCOMP )单步优化(sin gle)体积(VOLUME )2多步优化(Multiple)体积(VOLUME )3频率优化(TOFREQ )单阶频率优化(Sin gle)体积(VOLUME )4加权平均频率优化(Weighted Mea n)体积(VOLUME )5一致平均频率优化(Reciprocal Mean)体积(VOLUME )6欧几里德范数频率优化(Eucli

13、dea n norm)体积(VOLUME )7体积(VOLUME )单步结构载荷8体积(VOLUME )多部结构载荷1.5.3、初始化优化过程初始化优化过程主要是利用 TOTYPE和TODEF命令分别设定优化计算执 行的方式和终止计算(或者计算收敛的精度)，具体说明如下：(1) TOTYPE :设定优化计算执行方式，主要有优化准则方式(OptimalityCriteria)和顺序凸状程序方式(Sequential Convex Programming OC 方式 可以定义仅将体积作为约束条件，而 SCP方式必须定义优化计算的目标。(2) TODEF :设定优化计算的终止方式或者收敛精度。说明：

14、优化计算的设定参数并不保存到数据库中，如果用户利用RESUME恢复数据库，贝U需要重新禾I用 TOCOMP、TOFREQ、TOVAR、TOTYPE、TODEF等命令进行重新设定。时磊忖呎1.5.4、执行优化计算执行优化计算分为单次优化计算和批量多次自动优化计算两种方式，以下分别说明：（1）利用SOLVE TOEXE方式，看看下面的代码片段：TOCOMP，COMP，Single 1TOVAR，COMP，OBJTOVAR，VOLUME，CON，25TOTYPE， OCTODEF，1.0d-4/SOLUTIONSOLVETOEXEFINISH/POST1PLNSOL， TOPO*GET，TOPSTA

15、T，TOPO,，CONV*STAT，TOPSTAT/SOLUTIONSOLVETOEXEFINISH/P0ST1说明：首先利用TOCOMP、TOVAR等命令设定好优化定义，之后进入求 解模块，利用SOLVE TOEXE命令组合求解单次优化计算，之后进入后处理模 块，显示拓扑优化的结果；再进入求解模块，再一次执行优化计算，直到满足用 户的计算目标为止。再看看一个利用单步优化计算进行频率优化的代码片段：TOFREQ，FREQ1, SINGLE，1TOVAR，FREQ1，OBJTOVAR，VOLUME，CON，25TOTYPE， SCPTODEF，1.0d-4/SOLUTIONANTYPE， MOD

16、ALMODOPT，LANB，1MXPAND，1,，YESSOLVETOEXEFINISHTOPLOT, 0*GET , TOPSTAT, TOPO, CONV*STAT, TOPSTAT/SOLUTIONSOLVETOEXEFINISHTOPLOT, 0说明：采用单步优化频率计算过程类似于单步结构优化计算过程，不再赘述。（2）自动批量多次优化计算该方式是ANSYS中推荐采用的优化计算方式，主要命令是 TOLOOP，其 本身是一个宏命令，看看下面的代码片段，并要注意和SOLVE TOEXE方式的 区别：LSWRITELSWRITELSWRITETOCOMP , MCOMP, MULTIPLE ,

17、 3TOVAR , MCOMP, OBJTOVAR , VOLUME , CON , 30TODEF , 0.001TOLOOP , 20, 1说明：本段程序代码首先定义 3个载荷步，之后定义了优化分析的设定等，最后利用TOLOOP命令执行最多20次的优化计算迭代。1.6、查看优化计算的结果当优化计算完毕，优化的结果存放在结果文件中（*.RST），用户可以使用如下的命令格式查看优化计算的结果：（1）对于节点计算的结果，用PLNSOL，TOPO或者PRNSOL，TOPO显示虚密度（虚密度是指去除和保留部分的指标，是一个 01之间的值，虚密 度为0,表示可以去除；虚密度为 1,表示应该保留），当然

18、也可以使用 TOPLOT，0命令显示节点计算结果。（2） 对于单元计算的结果，用 PLESOL，TOPO或者PRESOL，TOPO显示虚密度，当然也可以使用TOPLOT，1命令显示单元计算结果。（3） 另外，也可以使用TOGRAPH和TOPRINT显示优化计算过程，看看下面的代码片段：/P0ST1TOPLOT, 1PLNSOL, TOPOTOGRAPH, OBJTOGRAPH , CON , VOLUMETOPRINT , OBJTOPRINT , CON(4) 也可以利用ANSYS中的单元表格功能进行结果显示，看看下面的代码片段：ETABLE，EDENS，TOPOPLETAB， EDENSP

19、RETAB， EDENSESEL，S, ETAB，EDENS，0.9, 1.0EPLOT(5) 如何查看优化计算过程是否收敛和迭代的次数？看看下面的代码片段:*GET，TOPCV，TOPO , CONV*GET，TITER，TOPO，ITER*GET，TOBJ，TOPO，ITER 1, TOHO*GET，TCON，TOPO，ITER 1，TOHC，1*STAT，说明：变量TOPCV = 1表示优化计算收敛，TITER变量表示优化计算的次 数，TOBJ变量表示优化计算最终的目标值，TCON变量表示约束条件值。(6) 对于频率优化计算，特别是对于加权平均、一致性平均和欧几里德范数，可以通过如下的代

20、码片段命令获得优化计算后的频率：*GET，FREQ1, MODE，1，FREQ*GET，FREQ2，MODE，2,FREQGET , FREQ3, MODE , 3, FREQ时需Sr彳说明：可以进入后处理模块（/POST1）和利用SET和PLDISP命令获得对 应的振形。时磊5说-2、拓扑优化有关的命令简要说明总体说明：拓扑优化的命令均以 TO开头，共计12个，现简单说明如下：(1) TOCOMP, REFName, Type, NUMLC , LCARR命令功能：定义单步或者多部一致性结果拓扑优化任务。参数解释：REFName：优化任务名称，最多8个字符；Type:优化载荷步类型，可以是

21、Single (单步)或者Multiply (多步);NUMLC :对于单步载荷，指定对哪个载荷步进行优化计算；对于多步载荷，指 定需要考虑的载荷步总数；LCARR：对于多步载荷，利用该参数可以定义多个载荷步之间权重，需要将各个载荷步的权重先定义在一个数组中，再用权重数组的方式加以引用。(2) TOFREQ， REFName，Type, NFreq， Frqarr，TargVal命令功能：定义单个或者多个频率拓扑优化任务。参数解释：REFName：优化任务名称，最多8个字符；Type:可以是 Single Weighted Reciproca、Euclidean四种优化形式；NFreq：对于T

22、ype= Single,该参数表示频率的标志值； 对于 Weighted ReciprocalEuclidean这三种方式，该参数表示频率的总体阶数；Frqarr：对于Type= Weighted Reciprocal该参数表示加权数组，用权重数组 的方式加以引用；对于Type= Euclidean,该参数表示目标频率；TargVal:对于Type= Reciprocal该参数用于定义目标频率。（3） TOVAR , REFName, Type, Lower, Upper, Boundtype 命令功能：定义优化变量和类型。参数解释：REFName：优化任务名称，最多 8个字符，它是利用 TOCOMP或者TOFREQ命令定义的；Type:优化变量类型，可以是 OBJ、CON和DEL ;Lower:当Type= CON时，优化变量的下限；Upper：当Type= CON时，优化变量的上限；Boundtype：当Type= CON时，表示优化变量的类型，可以是 PERCENT （百分 数）或者ACTUAL （实际值）。命令缺省：TOVAR，VOLUME，OBJ（4）TOTYPE，Type命令功能：定义拓扑优化求解类型。参数解释：Type:可以是OC （优化准则方式，默认值），或者SCP （顺序凸度程序方式）。 在