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

如何利用ANSYS进行拓扑优化前言就目前而言，利用有限元进行优化主要分成两个阶段：（1） 进行拓扑优化，明确零件最佳的外形、刚度、体积，或者合理的固有频率，主要目的是确定优化的方向；（2） 进行尺寸优化，主要目的是确定优化后的的零件具体尺寸值，通常是在完成拓扑优化之后，再执行尺寸优化。在ANSYS中，利用拓扑优化，可以完成以下两个目的：（1） 在特定载荷和约束的条件下，确定零件的最佳外形，或者最小的体积（或者质量）；（2） 利用拓扑优化，使零件达到需要的固有频率，避免在使用过程中产生共振等不利影响。本文主要就在ANSYS环境中如何执行拓扑优化进行说明。欢迎下载1、利用ANSYS进行拓扑优化的过程在ANSYS中，执行优化，通常分为以下6个步骤：定义需要求解的结构问题选择合理的优化单元类型设定优化和非优化的区域定义载荷步或者需要提取的频率对优化过程进行定义和控制计算并查看结果1.1、定义需要求解的结构问题对于结构进行优化分析，定义结构的物理特性必不可少，例如，需要定义结构的杨氏模量、泊松比（其值在0.10.4之间）、密度等相关的结构特性方面的信息，以供结构计算能够正常执行下去。1.2、选择合理的优化单元类型在ANSYS中，不是所有的单元类型都可以执行优化的，必须满足如下的规定：（1）2D平面单元：PLANE82单元和PLANE183单元；（2）3D实体单元：SOLID92单元和SOLID95单元；（3）壳单元：SHELL93单元。上述单元的特性在帮助文件中有详细的说明，同时对于2D单元，应使用平面应力或者轴对称的单元选项。1.3、指定优化和非优化的区域在ANSYS中规定，单元类型编号为1的单元，才执行优化计算；否则，就不执行优化计算。例如，对于结构分析中，对于不能去除的部分区域将单元类型编号设定为2，就可以不执行优化计算，请见下面的代码片段：Et,1,solid92Et,2,solid92Type,1Vsel,s,num,1,2Vmesh,allType,2Vsel,s,num,3Vmesh,all说明：上述代码片段定义相同的单元类型（solid92），但编号分别为1和2，并将单元类型编号1利用网格划分分配给了1体和2体，从而对其进行优化计算；而单元编号为2利用网格划分分配给了3体，从而不执行优化计算。1.4、定义载荷步或者需要提取的频率1.4.1 线性结构静态分析对于结构优化而言，其总是在特定的载荷（或者载荷步），约束和目标下进行的，在优化分析的过程中，必须执行线性结构静态分析，才能获得需要的优化之后的形状。在ANSYS中，可以对单步载荷或者多步载荷执行优化分析，当然，单步载荷是最简单的了。然而，对于某个特定载荷步，必须使用LSWRITE载荷步存储命令将载荷步预存起来，再用LSSOLVE命令进行求解。先看看下面的代码片段：D,10,all,0,20,1Nsel,s,loc,y,0Sf,AllselLswrite,1Ddel,Sfdel,Nsel,s,loc,x,0,1D,all,all,0F,212,fx,Lswrite,2Lswrite,3FinishTocomp,mcomp,multiple,3Tovar,mcomp,objTovar,volume,con,10TodefToloop,20说明：该代码片段首先定义了3个载荷步，并利用LSWRITE命令将载荷步预存；之后利用Tocomp命令定义优化任务目标名称mcomp，并将体积减少10作为优化的约束条件，之后用Todef初始化优化过程，最后利用Toloop命令执行优化计算，最大计算次数20次。相关命令：TOCOMP、TOVAR、TODEF、TOEXE、TOLOOP和简要说明。（1） TOCOMP：定义结构优化任务目标。（如何理解COMPLIANCE：Compliance本意是一致性，统一性，在结构优化分析中，特别是对于多个载荷步，需要在多个载荷步之间取得一致性的结果，才能满足结构优化分析的目标）（2） TOVAR：定义优化变量，可以是目标变量，也可以是约束变量等；（3） TODEF：定义优化的初始化条件或者收敛准则；（4） TOEXE：执行单次优化计算；（5） TOLOOP：批量执行多次优化分析计算。说说TOEXE和TOLOOP之间的区别：TOEXE执行单次优化分析计算，其本身不执行结构分析过程，因此，在利用TOEXE命令执行优化计算之前，需要利用SOLVE或者LSSOLVE命令先执行结构静态分析计算；而TOLOOP是一个执行优化计算的宏命令，其中包含了SOLVE和LSSOLVE等命令，因此在上述代码片段中没有出现SOLVE或者LSSOLVE命令。就使用的便利性而言，利用TOLOOP命令可能更方便，但是利用TOEXE命令用户可以创建自己的优化宏命令，各有所长，主要是看用户如何使用这两个命令了。1.4.2 模态优化分析在执行模态优化分析之前，需要使用TOFREQ和TOVAR定义模态分析任务名称和优化分析变量，这点和上述的结构优化分析过程类似。值得注意的是，在利用MXPAND指定模态频率分析计算中，需要将单元计算设定为“YES”，这样才能利用TOEXE命令正常执行优化计算过程，当然在MXPAND之前还要使用modopt命令设定模态计算的选项；而利用TOLOOP命令执行模态优化分析计算过程中，由于TOLOOP是编制好的宏命令，不需要使用MXPAND等命令，请见下面的代码片段：/solutionAntype,modalFinishTofreq,mfreq,reciprocal,3Tovar,mfreq,objTovar,volume,con,50Todef,1.0d-4Toloop,20说明：首先定义模态求解作为分析类型，之后利用tofreq设定频率优化作为优化任务，名称为mfreq（其中reciprocal表示多阶模态频率，本例中为前3阶模态），并利用命令tovar设定了体积减少50作为优化的约束条件，之后设定优化的收敛准则，最后利用toloop命令执行最多20次优化迭代计算。在本例中，因为toloop是一个宏命令，所以并没有显式的使用modopt和mxpand命令，如果使用toexe命令，则用户必须显式的使用modopt和mxpand命令。1.5、对优化过程进行定义和控制在ANSYS中执行优化过程有2种方式：（1） 采用solve和toexe 命令相结合的机制，一步一步的执行优化计算，直到满足用户的优化目标和约束条件为止。（2） 采用toloop宏命令执行优化计算，可执行多次优化分析计算。说明：采用第二种方式操作更简单，推荐采用。总体上来说，执行优化计算过程大致上分为4个步骤：（1） 定义优化函数（或者方式和任务）；（2） 定义优化目标或者约束条件；（3） 初始化优化过程；（4） 执行优化计算（可以是单次计算，也可以自动批量计算）。以下分别进行描述：1.5.1、定义优化函数在ANSYS中，主要采用以下两种方式定义优化函数：（1） 利用TOCOMP定义结构优化；（2） 利用TOFREQ定义频率优化。说明：可以利用TOLIST命令例举出所有定义了的优化函数。1.5.2、定义优化目标或者约束条件定义优化目标或者约束条件主要是利用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命令设定优化类型，具体说明可以参见帮助文件）。同样，也可以利用如下的命令格式删除已经定义好的优化任务名称或者约束条件：Tovar，Reference_Name，DEL可以定义的优化目标和约束条件的组合方式见下表：序号优化类型优化目标约束条件1结构优化（TOCOMP）单步优化（single）体积（VOLUME）2多步优化（Multiple）体积（VOLUME）3频率优化（TOFREQ）单阶频率优化（Single）体积（VOLUME）4加权平均频率优化（Weighted Mean）体积（VOLUME）5一致平均频率优化（Reciprocal Mean）体积（VOLUME）6欧几里德范数频率优化（Euclidean norm）体积（VOLUME）7体积（VOLUME）单步结构载荷8体积（VOLUME）多部结构载荷1.5.3、初始化优化过程初始化优化过程主要是利用TOTYPE和TODEF命令分别设定优化计算执行的方式和终止计算（或者计算收敛的精度），具体说明如下：（1） TOTYPE：设定优化计算执行方式，主要有优化准则方式（Optimality Criteria）和顺序凸状程序方式（Sequential Convex Programming），OC方式可以定义仅将体积作为约束条件，而SCP方式必须定义优化计算的目标。（2） TODEF：设定优化计算的终止方式或者收敛精度。说明：优化计算的设定参数并不保存到数据库中，如果用户利用RESUME恢复数据库，则需要重新利用TOCOMP、TOFREQ、TOVAR、TOTYPE、TODEF等命令进行重新设定。1.5.4、执行优化计算执行优化计算分为单次优化计算和批量多次自动优化计算两种方式，以下分别说明：（1） 利用SOLVETOEXE方式，看看下面的代码片段：TOCOMP，COMP，Single，1TOVAR，COMP，OBJTOVAR，VOLUME，CON，25TOTYPE，OCTODEF，1.0d-4/SOLUTIONSOLVETOEXEFINISH/POST1PLNSOL，TOPO*GET，TOPSTAT，TOPO，CONV*STAT，TOPSTAT/SOLUTIONSOLVETOEXEFINISH/POST1说明：首先利用TOCOMP、TOVAR等命令设定好优化定义，之后进入求解模块，利用SOLVETOEXE命令组合求解单次优化计算，之后进入后处理模块，显示拓扑优化的结果；再进入求解模块，再一次执行优化计算，直到满足用户的计算目标为止。再看看一个利用单步优化计算进行频率优化的代码片段：TOFREQ，FREQ1，SINGLE，1TOVAR，FREQ1，OBJTOVAR，VOLUME，CON，25TOTYPE，SCPTODEF，1.0d-4/SOLUTIONANTYPE，MODALMODOPT，LANB，1MXPAND，1，YESSOLVETOEXEFINISHTOPLOT，0*GET，TOPSTAT，TOPO，CONV*STAT，TOPSTAT/SOLUTIONSOLVETOEXEFINISHTOPLOT，0说明：采用单步优化频率计算过程类似于单步结构优化计算过程，不再赘述。（2） 自动批量多次优化计算该方式是ANSYS中推荐采用的优化计算方式，主要命令是TOLOOP，其本身是一个宏命令，看看下面的代码片段，并要注意和SOLVETOEXE方式的区别：LSWRITELSWRITELSWRITETOCOMP，MCOMP，MULTIPLE，3TOVAR，MCOMP，OBJTOVAR，VOLUME，CON，30TODEF，0.001TOLOOP，20，1说明：本段程序代码首先定义3个载荷步，之后定义了优化分析的设定等，最后利用TOLOOP命令执行最多20次的优化计算迭代。1.6、查看优化计算的结果当优化计算完毕，优化的结果存放在结果文件中（*.RST），用户可以使用如下的命令格式查看优化计算的结果：（1） 对于节点计算的结果，用PLNSOL，TOPO或者PRNSOL，TOPO显示虚密度（虚密度是指去除和保留部分的指标，是一个01之间的值，虚密度为0，表示可以去除；虚密度为1，表示应该保留），当然也可以使用TOPLOT，0命令显示节点计算结果。（2） 对于单元计算的结果，用PLESOL，TOPO或者PRESOL，TOPO显示虚密度，当然也可以使用TOPLOT，1命令显示单元计算结果。（3） 另外，也可以使用TOGRAPH和TOPRINT显示优化计算过程，看看下面的代码片段：/POST1TOPLOT，1PLNSOL，TOPOTOGRAPH，OBJTOGRAPH，CON，VOLUMETOPRINT，OBJTOPRINT，CON（4） 也可以利用ANSYS中的单元表格功能进行结果显示，看看下面的代码片段：ETABLE，EDENS，TOPOPLETAB，EDENSPRETAB，EDENSESEL，S，ETAB，EDENS，0.9，1.0EPLOT（5） 如何查看优化计算过程是否收敛和迭代的次数？看看下面的代码片段：*GET，TOPCV，TOPO，CONV*GET，TITER，TOPO，ITER*GET，TOBJ，TOPO，ITER1，TOHO*GET，TCON，TOPO，ITER1，TOHC，1*STAT，说明：变量TOPCV1表示优化计算收敛，TITER变量表示优化计算的次数，TOBJ变量表示优化计算最终的目标值，TCON变量表示约束条件值。（6） 对于频率优化计算，特别是对于加权平均、一致性平均和欧几里德范数，可以通过如下的代码片段命令获得优化计算后的频率：*GET，FREQ1，MODE，1，FREQ*GET，FREQ2，MODE，2，FREQ*GET，FREQ3，MODE，3，FREQ说明：可以进入后处理模块（/POST1）和利用SET和PLDISP命令获得对应的振形。2、拓扑优化有关的命令简要说明总体说明：拓扑优化的命令均以TO开头，共计12个，现简单说明如下：（1）TOCOMP，REFName，Type，NUMLC，LCARR命令功能：定义单步或者多部一致性结果拓扑优化任务。参数解释：REFName：优化任务名称，最多8个字符；Type：优化载荷步类型，可以是Single（单步）或者Multiply（多步）；NUMLC：对于单步载荷，指定对哪个载荷步进行优化计算；对于多步载荷，指定需要考虑的载荷步总数；LCARR：对于多步载荷，利用该参数可以定义多个载荷步之间权重，需要将各个载荷步的权重先定义在一个数组中，再用权重数组的方式加以引用。（2）TOFREQ，REFName，Type，NFreq，Frqarr，TargVal命令功能：定义单个或者多个频率拓扑优化任务。参数解释：REFName：优化任务名称，最多8个字符；Type：可以是Single、Weighted、Reciprocal、Euclidean四种优化形式；NFreq：对于TypeSingle，该参数表示频率的标志值；对于Weighted、Reciprocal、Euclidean这三种方式，该参数表示频率的总体阶数；Frqarr：对于TypeWeighted、Reciprocal，该参数表示加权数组，用权重数组的方式加以引用；对于TypeEuclidean，该参数表示目标频率；TargVal：对于TypeReciprocal，该参数用于定义目标频率。（3） TOVAR，REFName，Type，Lower，Upper，Boundtype命令功能：定义优化变量和类型。参数解释：REFName：优化任务名称，最多8个字符，它是利用TOCOMP或者TOFREQ命令定义的；Type：优化变量类型，可以是OBJ、CON和DEL；Lower：当TypeCON时，优化变量的下限；Upper：当TypeCON时，优化变量的上限；Boundtype：当TypeCON时，表示优化变量的类型，可以是PERCENT（百分数）或者ACTUAL（实际值）。命令缺省：TOVAR，VOLUME，OBJ（4） TOTYPE，Type命令功能：定义拓扑优化求解类型。参数解释：Type：可以是OC（优化准则方式，默认值），