ansys优化方法简介以及实例.pdf

该【ansys优化方法简介以及实例 】是由【小屁孩】上传分享，文档一共【9】页，该文档可以免费在线阅读，需要了解更多关于【ansys优化方法简介以及实例 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。 : .
ansys优化方法简介以及实例--第1页
拓扑优化理论及在 ANSYS软件中的实现
一．拓扑优化概论：
连续体结构的拓扑优化设计是继结构的尺寸优化设计和形状优化设计之后，在结构优化领域
出现的一种富有挑战性的研究方向，它是一种比尺寸优化和形状优化更高层次的优化方法，
也是结构优化问题中最为复杂的一类问题。 拓扑优化处于结构的概念设计阶段 ，其优化结果
是一切后续设计的基础。因而在初始设计阶段需要确定结构的最佳拓扑形式。 拓扑优化的目
的是寻求结构的刚度在设计空间最佳的分布形式， 或在设计域空间需求结构最佳的传力路线，
以优化结构的某些性能或减轻结构的重量。
目前对于拓扑优化的研究主要集中在以下几个方面：
结构拓扑描述方式和材料插值模型；
拓扑优化中结构拓扑描述方式和材料的插值模型非常重要， 是一切后续拓扑优化工作的基础。
常用的拓扑描述方式和材料插值模型有均匀化方法、 密度法、 变厚度法和拓扑函数描述法等。
拓扑优化求解数值算法，新型优化算法在拓扑优化中的应用；
拓扑优化的数值计算方法主要包括有限元法和无网格法，基于成熟的有限元理论的拓扑优化
格式简单，便于实现，但在优化过程中常因网格的重分和细化导致计算困难，结构中常出现
中间密度材料、棋盘格现象和网格依赖性等问题。无网格法是今年发展的一种新型数值求解
技术，摆脱了有限元繁琐的网格生成过程，从理论上看比有限元法拥有更广阔的应用前景，
但目前尚处于发展和完善中。
拓扑优化的特点是： 设计变量多，计算规模大，目标函数和约束函数一般为设计变量的非线
性、非单调函数。 目前应用于连续体结构拓扑优化计算的优化算法主要包括两类：优化准则
法和序列凸规划法。
去除优化过程中数值计算不稳定的方法，优化结果的提取和重构；
拓扑优化中经常出现的数值计算问题有： 多孔材料、棋盘格现象、网格依赖性和局部极值问
题。优化结果的提取和应用主要考虑的是如何将优化的结果转化为可用的 CAD模型问题，实
现 CAE和 CAD之间的数据共享和交流。
随着拓扑优化理论研究的不断深入，拓扑优化在航空和汽车领域已开始得到初步的应用，主
要是通过拓扑优化获得结构的最初拓扑形式，并在最初拓扑形式的基础上进行相关的后续优
化设计。解决的问题范围包括：线弹性静态结构优化问题、动力优化问题及非线性等复杂情
况下的优化问题。
二．ANSYS中拓扑优化相关理论及应用：
ANSYS中拓扑优化技术采用的是均匀化方法，具有优化准则法和序列凸规划两种优化算法。
1．均匀化方法:
均匀化方法是一种经典的拓扑优化方法，有着严密的数学和力学理论基础。连续体结构拓扑
优化的均匀化方法（ Homogenization Method for Optimization ）是 Bendsoe 等人于 1988
年提出的。其基本思想是：在组成拓扑结构的材料中引入细观结构，以宏观解结构单元模型
对设计区域进行有限元离散划分，用周期性细观结构来描述宏观单元，优化过程中 以细观结
构的几何尺寸作为设计变量 ，把弹性模量、材料密度等参量表示成细观结构几何尺寸变量的
函数。以细观结构的消长实现材料的增减。并产生介于由中间尺寸细观结构组成的复合材料，
从而实现结构拓扑优化模型与尺寸优化模型的统一。它将复杂的拓扑优化问题挂靠在低层次
的尺寸优化变量问题上来求解，但求解过程中均匀化弹性张量计算非常复杂，且微单元的最
佳形状和方向难以确定， 结构响应函数的密度求解复杂，优化变量过多， 计算效率低等缺点，
主要用于拓扑优化理论方面的研究。
2．优化算法：
目前应用于连续体结构拓扑优化计算的优化算法主要包括两类：优化准则法和序列凸规划法。
ansys优化方法简介以及实例--第1页 : .
ansys优化方法简介以及实例--第2页
优化准则法收敛速度快，计算过程不使用导数信息，但其一般适用于 单约束条件下的问题优
化，且不同的优化问题需要推导不同的优化准则。序列凸规划算法包括序列线性规划方法、
序列二次规划算法以及移动近似算法， 其中移动近似算法 （MMA）是目前使用最为广泛的算法
之一，能够广泛应用于多约束情况，其计算过程中要使用前一步或多步的计算信息。
3．应用：
ANSYS中的拓扑优化工具可用来解决以下问题：
（1） 体积约束下的最大刚度设计：以柔顺度为目标函数，体积为约束函数；
（2） 刚度约束下的最小体积优化：以体积为目标函数，刚度为约束函数；
（3） 体积约束下的最大动刚度设计：以 n 阶自振频率为目标函数 ,体积为约束函数；
（4） 以上多种工况的组合优化问题；


形状优化实例
FINISH
/CLEAR,START
/TRIAD,OFF ! 关闭整体直角坐标系的三角符号
H=1000 ! 设置比例尺，采用 iso
TK16=H !设置参数变量并附初值
TK27=6,35/H
TK38=H
TK49=H
TK50=4/H
/PREP7
ET,1,PLANE42
MP,EX,1,
MP,NUXY,1,
K,1
K,5,254/H
KFILL! 在第 1 至第 5 个关键点之间生成 2,3,4 关键点
K,6,,TK16
K,7,H,TK27
K,8,127/H,TK38
K,9,H,TK49
K,10,254/H,TK50
SPLINE,6,7,8,9,10
L,1,6
*REPEAT,5,1,1 ! 重复 l 命令，关键点编号自动加 1，分别在 2,7 、3,8 、4,9 、5,10 之间共
生成 4 条线其中包含了、命令已生成的线，共 5 条
LSEL,S,LINE,,5,9 ! 选择上述生成的 5 条线
LESIZE,ALL,,,1 ! 指定线在划分网格前的等分数为 1
LSEL,ALL
A,1,2,7,6
*REPEAT,4,1,1,1,1 ! 重复上述命令，共生成 4 个面
ESIZE,,4
ansys优化方法简介以及实例--第2页 : .
ansys优化方法简介以及实例--第3页
AMESH,ALL
FINISH
/SOLU
NSEL,S,LOC,Y
DSYM,SYMM,X !对选择的节点施加 x 方向的对称约束
NSEL,S,LOC,X
DSYM,ASYM,Y !对选择的节点施加 x 方向的反对称约束
NSEL,ALL
FK,10,FX,66725*4 ! 在 10 号关键点施加集中载荷，实现弯矩
DK,1,ALL,0 ! 在 1 号关键点施加全约束
SOLVE
FINISH
/POST
SET,LAST
ETABLE,EVOL,VOLU ! 建立单元表，并取出每个单元的体积
PRNSOL,S,PRIN ! 列出节点的主应力
NSORT,S,1
NSEL,U,LOC,X,0,230/H ! 选择介于 0 到 230/H 的节点
*GET,STRS,SORT,,MAX ! 取出最大的应力值并赋给 strs
NSEL,ALL
SSUM ! 体积相加
*GET,VTOT,SSUM,,ITEM,VOLU !取出结构总体积
TVOL=TVOL*2 !由于分析时只计算了结构的一半，总体积要乘 2
NSEL,U,LOC,X,250/H,265/H ! 选择介于 250/H 到 265/H 的节点
PRNSOL,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-TK38
DIF3=TK38-TK49
FINISH
/OPT
opanl,SCRATCH ! 指定分析文件
opvar,TVOL,Obj,,,1/H ! 定义优化目标函数 ,收敛误差为 1
opvar,STRS,sv,,206E6 ! 定义状态变量即优化过程的约束条件
opvar,DEFL,sv,,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
ansys优化方法简介以及实例--第3页 : .
ansys优化方法简介以及实例--第4页
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 ! 将参数的值输出到文件中且文件名为 rset1
oplist,all,,1 ! 列表显示所有序列
/AXLAB,Y,TVOL !设置曲线输出时 y 轴的说明
plvaropt,TVOL ! 绘图显示目标函数随优化次数的变化规律
FINISH

/POST1
PLNSOL,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 ! 指定内半径
TI=30 ! 壁厚
FIL=10 ! 倒角半径
/PREP7
CYL4,-200,,INRAD,-30,INRAD+T1,30 ! 生成 3 个部分圆环面
CYL4,200*COS(60),200*SIN(60),INRAD,-90,INRAD+T1,150
CYL4,200*COS(60),200*SIN(-60),INRAD,90,INRAD+T1,150
AADD,ALL !面相加
ADELE,ALL !仅删除面，线保留
LSEL,,RADIUS,,INRAD+T1 ! 选择零件内控的线
L1=LSNEXT(0) !得到编号最小的线的编号
L2=LSNEXT(L1) !得到下一个编号最小的线的编号
L3=LSNEXT(L2)
LFILLET,L1,L2,FIL ! 对内控的线进行倒角
LFILLET,L2,L3,FIL