《ansys优化方法简介以与实例》由会员分享，可在线阅读，更多相关《ansys优化方法简介以与实例（9页珍藏版）》请在金锄头文库上搜索。
1、拓扑优化理论及在ANSYS软件中的实现一拓扑优化概论：连续体结构的拓扑优化设计是继结构的尺寸优化设计和形状优化设计之后，在结构优化领域出现的一种富有挑战性的研究方向，它是一种比尺寸优化和形状优化更高层次的优化方法，也是结构优化问题中最为复杂的一类问题。拓扑优化处于结构的概念设计阶段，其优化结果是一切后续设计的基础。因而在初始设计阶段需要确定结构的最佳拓扑形式。拓扑优化的目的是寻求结构的刚度在设计空间最佳的分布形式，或在设计域空间需求结构最佳的传力路线，以优化结构的某些性能或减轻结构的重量。目前对于拓扑优化的研究主要集中在以下几个方面：结构拓扑描述方式和材料插值模型；拓扑优化中结构拓扑描述方式和
2、材料的插值模型非常重要，是一切后续拓扑优化工作的基础。常用的拓扑描述方式和材料插值模型有均匀化方法、密度法、变厚度法和拓扑函数描述法等。拓扑优化求解数值算法，新型优化算法在拓扑优化中的应用；拓扑优化的数值计算方法主要包括有限元法和无网格法，基于成熟的有限元理论的拓扑优化格式简单，便于实现，但在优化过程中常因网格的重分和细化导致计算困难，结构中常出现中间密度材料、棋盘格现象和网格依赖性等问题。无网格法是今年发展的一种新型数值求解技术，摆脱了有限元繁琐的网格生成过程，从理论上看比有限元法拥有更广阔的应用前景，但目前尚处于发展和完善中。拓扑优化的特点是：设计变量多，计算规模大，目标函数和约束函数一般
3、为设计变量的非线性、非单调函数。目前应用于连续体结构拓扑优化计算的优化算法主要包括两类：优化准则法和序列凸规划法。去除优化过程中数值计算不稳定的方法，优化结果的提取和重构； 拓扑优化中经常出现的数值计算问题有：多孔材料、棋盘格现象、网格依赖性和局部极值问题。优化结果的提取和应用主要考虑的是如何将优化的结果转化为可用的CAD模型问题，实现CAE和CAD之间的数据共享和交流。随着拓扑优化理论研究的不断深入，拓扑优化在航空和汽车领域已开始得到初步的应用，主要是通过拓扑优化获得结构的最初拓扑形式，并在最初拓扑形式的基础上进行相关的后续优化设计。解决的问题围包括：线弹性静态结构优化问题、动力优化问题及非
4、线性等复杂情况下的优化问题。二ANSYS中拓扑优化相关理论及应用：ANSYS中拓扑优化技术采用的是均匀化方法，具有优化准则法和序列凸规划两种优化算法。1均匀化方法:均匀化方法是一种经典的拓扑优化方法，有着严密的数学和力学理论基础。连续体结构拓扑优化的均匀化方法（Homogenization Method for Optimization）是Bendsoe等人于1988年提出的。其基本思想是：在组成拓扑结构的材料中引入细观结构，以宏观解结构单元模型对设计区域进行有限元离散划分，用周期性细观结构来描述宏观单元，优化过程中以细观结构的几何尺寸作为设计变量，把弹性模量、材料密度等参量表示成细观结构几何
5、尺寸变量的函数。以细观结构的消长实现材料的增减。并产生介于由中间尺寸细观结构组成的复合材料，从而实现结构拓扑优化模型与尺寸优化模型的统一。它将复杂的拓扑优化问题挂靠在低层次的尺寸优化变量问题上来求解，但求解过程中均匀化弹性量计算非常复杂，且微单元的最佳形状和方向难以确定，结构响应函数的密度求解复杂，优化变量过多，计算效率低等缺点，主要用于拓扑优化理论方面的研究。2优化算法：目前应用于连续体结构拓扑优化计算的优化算法主要包括两类：优化准则法和序列凸规划法。优化准则法收敛速度快，计算过程不使用导数信息，但其一般适用于单约束条件下的问题优化，且不同的优化问题需要推导不同的优化准则。序列凸规划算法包括
6、序列线性规划方法、序列二次规划算法以及移动近似算法，其中移动近似算法（MMA）是目前使用最为广泛的算法之一，能够广泛应用于多约束情况，其计算过程中要使用前一步或多步的计算信息。3应用：ANSYS中的拓扑优化工具可用来解决以下问题：（1） 体积约束下的最大刚度设计：以柔顺度为目标函数，体积为约束函数；（2） 刚度约束下的最小体积优化：以体积为目标函数，刚度为约束函数；（3） 体积约束下的最大动刚度设计：以n阶自振频率为目标函数,体积为约束函数；（4） 以上多种工况的组合优化问题；形状优化实例FINISH/CLEAR,START/TRIAD,OFF !关闭整体直角坐标系的三角符号H=1000 !设
7、置比例尺，采用isoTK16=6.35/H !设置参数变量并附初值TK27=6,35/HTK38=6.35/HTK49=6.35/HTK50=4/H/PREP7ET,1,PLANE42MP,EX,1,6.89E10MP,NUXY,1,0.3K,1K,5,254/HKFILL!在第1至第5个关键点之间生成2,3,4关键点K,6,TK16K,7,63.5/H,TK27K,8,127/H,TK38K,9,190.5/H,TK49K,10,254/H,TK50SPLINE,6,7,8,9,10L,1,6*REPEAT,5,1,1 !重复l命令，关键点编号自动加1，分别在2,7、3,8、4,9、5,10
8、之间共生成4条线其中包含了、命令已生成的线，共5条LSEL,S,LINE,5,9 !选择上述生成的5条线LESIZE,ALL,1 !指定线在划分网格前的等分数为1LSEL,ALLA,1,2,7,6*REPEAT,4,1,1,1,1 !重复上述命令，共生成4个面ESIZE,4AMESH,ALLFINISH/SOLUNSEL,S,LOC,YDSYM,SYMM,X !对选择的节点施加x方向的对称约束NSEL,S,LOC,XDSYM,ASYM,Y !对选择的节点施加x方向的反对称约束NSEL,ALLFK,10,FX,66725*4 !在10号关键点施加集中载荷，实现弯矩DK,1,ALL,0 !在1号关
9、键点施加全约束SOLVEFINISH/POSTSET,LASTETABLE,EVOL,VOLU !建立单元表，并取出每个单元的体积PRNSOL,S,PRIN !列出节点的主应力NSORT,S,1NSEL,U,LOC,X,0,230/H !选择介于0到230/H的节点*GET,STRS,SORT,MAX !取出最大的应力值并赋给strsNSEL,ALLSSUM !体积相加*GET,VTOT,SSUM,ITEM,VOLU !取出结构总体积TVOL=TVOL*2 !由于分析时只计算了结构的一半，总体积要乘2NSEL,U,LOC,X,250/H,265/H !选择介于250/H到265/H的节点PRN
10、SOL,U,Y !列表显示出所选节点在y向的位移值NSORT,U,Y,1 !位移值升序排列PRNSOL,U,Y !列出排序后的结果*GET,DEFL,SORT,MAX !取出最大的位移值赋给defl*STATUS,PARM !显示当前参数变量的状态值DEFL=ABS(DEFL)DIF1=TK16-TK27 !设置参数值，以保证曲线的光滑性DIF2=TK27-TK38DIF3=TK38-TK49FINISH/OPTopanl,SCRATCH !指定分析文件opvar,TVOL,Obj,1/H !定义优化目标函数,收敛误差为1opvar,STRS,sv,206E6 ! 定义状态变量即优化过程的约束
11、条件opvar,DEFL,sv,12.5/H opvar,DIF1,sv,1/H opvar,DIF2,sv,1/H opvar,DIF3,sv,1/H opvar,TK16,dv,4/H,7/H !定义第一个设计变量及其变化围opvar,TK27,dv,4/H,7/H opvar,TK38,dv,4/H,7/H opvar,TK4,9dv,4/H,7/H opsave,INITIAL,opt !保存所有的优化数据到文件optp里optype,SUBP !设置优化方式,子问题逼近算法OPSUBP,30 !指定迭代次数opexe ! 执行优化PARSAV,RSET1 !将参数的值输出到文件中且文
12、件名为rset1oplist,all,1 !列表显示所有序列/AXLAB,Y,TVOL !设置曲线输出时y轴的说明plvaropt,TVOL !绘图显示目标函数随优化次数的变化规律FINISH/POST1PLNSOL,U,SUM,0,1 !显示结构在优化后的总位移分布PLNSOL,S,EQV,0,1 !显示结构在优化后的应力分布FINISH尺寸优化实例FINISH/CLEAR,START/TRIAD,OFF !关闭整体直角坐标系的三角符号 !CREATE,HEXPLATE生成一个宏文件，名为hexplate*AFUN,DEG !设置三角函数的角度单位为INRAD=200*COS(30)-20
13、!指定半径TI=30 !壁厚FIL=10 !倒角半径/PREP7CYL4,-200,INRAD,-30,INRAD+T1,30 !生成3个部分圆环面CYL4,200*COS(60),200*SIN(60),INRAD,-90,INRAD+T1,150CYL4,200*COS(60),200*SIN(-60),INRAD,90,INRAD+T1,150AADD,ALL !面相加ADELE,ALL !仅删除面，线保留LSEL,RADIUS,INRAD+T1 !选择零件控的线L1=LSNEXT(0) !得到编号最小的线的编号L2=LSNEXT(L1) !得到下一个编号最小的线的编号L3=LSNEXT
14、(L2) LFILLET,L1,L2,FIL !对控的线进行倒角LFILLET,L2,L3,FILLFILLET,L3,L1,FILLSEL,AL !选择所有的线WPROT,90 LSBW,ALL !所有线由工作平面分割WPROT,60LSBW,ALLCSYS,1 !设置坐标系为柱坐标系LSEL,U,LOC,Y,0,60 !选择介于0到60围之外的所有线LDELE,ALL,1 !删除所选择的线LSEL,ALL !选择所有的线KSLL !选择依附线上的关键点KSEL,INVE !选择不依附于线的关键点KDELE,ALL !删除不依附于线的关键点KSEL,ALL !选择所有的关键点CSYS,0 !设置坐标系为整体坐标系KSEL,LOC,Y,0 !选择y=0的关键点LSTR,KPNEXT(0),KPNEXT(KPNEXT(0) !对选择的关键点进行连线，即生成底线KSEL,ALLCSYS,1KSEL,LOC,Y,60 !选择theta=60的关键点LSTR,KPNEXT(0).KPNEXT(KPNEXT(0) !由选择的关键点生成线，即左侧边线KSEL,ALLLSEL,LOC,X,100 !选择x=100
