ansys的几种优化方法_命令流

  • 格式:txt
  • 大小:11.18 KB
  • 文档页数:5
/PREP7
ET,1,PLANE42
MP,EX,1,6.89E10
MP,NUXY,1,0.3
K,1
K,5,254/H
KFILL! 在第1至第5个关键点之间生成2,3,4关键点
K,6,,TK16
K,7,63.5/H,TK27
K,8,127/H,TK38
K,9,190.5/H,TK49
opanl,SCRATCH !指定分析文件
opvar,TVOL,Obj,,,1/H !定义优化目标函数,收敛误差为1
opvar,STRS,sv,,206E6 ! 定义状态变量即优化过程的约束条件
opvar,DEFL,sv,,12.5/H
*GET,DEFL,SORT,,MAX !取出最大的位移值赋给defl
*STATUS,PARM !显示当前参数变量的状态值
DEFL=ABS(DEFL)
DIF1=TK16-TK27 !设置参数值,以保证曲线的光滑性
DIF2=TK27-TK38
DIF3=TK38-TK49
FINISH
/OPT
2.优化算法:
目前应用于连续体结构拓扑优化计算的优化算法主要包括两类:优化准则法和序列凸规划法。
优化准则法收敛速度快,计算过程不使用导数信息,但其一般适用于单约束条件下的问题优化,且不同的优化问题需要推导不同的优化准则。序列凸规划算法包括序列线性规划方法、序列二次规划算法以及移动近似算法,其中移动近似算法(MMA)是目前使用最为广泛的算法之一,能够广泛应用于多约束情况,其计算过程中要使用前一步或多步的计算信息。
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 !取出结构总体积
ANSYS中拓扑优化技术采用的是均匀化方法,具有优化准则法和序列凸规划两种优化算法。
1.均匀化方法:
均匀化方法是一种经典的拓扑优化方法,有着严密的数学和力学理论基础。连续体结构拓扑优化的均匀化方法(Homogenization Method for Optimization)是Bendsoe等人于1988年提出的。其基本思想是:在组成拓扑结构的材料中引入细观结构,以宏观解结构单元模型对设计区域进行有限元离散划分,用周期性细观结构来描述宏观单元,优化过程中以细观结构的几何尺寸作为设计变量,把弹性模量、材料密度等参量表示成细观结构几何尺寸变量的函数。以细观结构的消长实现材料的增减。并产生介于由中间尺寸细观结构组成的复合材料,从而实现结构拓扑优化模型与尺寸优化模型的统一。它将复杂的拓扑优化问题挂靠在低层次的尺寸优化变量问题上来求解,但求解过程中均匀化弹性张量计算非常复杂,且微单元的最佳形状和方向难以确定,结构响应函数的密度求解复杂,优化变量过多,计算效率低等缺点,主要用于拓扑优化理论方面的研究。
形状优化实例
FINISH
/CLEAR,START
/TRIAD,OFF !关闭整体直角坐标系的三角符号
H=1000 !设置比例尺,采用iso
TK16=6.35/H !设置参数变量并附初值
TK27=6,35/H
TK38=6.35/H
TK49=6.35/H
TK50=4/H
拓扑优化的数值计算方法主要包括有限元法和无网格法,基于成熟的有限元理论的拓扑优化格式简单,便于实现,但在优化过程中常因网格的重分和细化导致计算困难,结构中常出现中间密度材料、棋盘格现象和网格依赖性等问题。无网格法是今年发展的一种新型数值求解技术,摆脱了有限元繁琐的网格生成过程,从理论上看比有限元法拥有更广阔的应用前景,但目前尚处于发展和完善中。
CSYS,0 !设置坐标系为整体坐标系
KSEL,,LOC,Y,0 !选择y=0的关键点
LSTR,KPNEXT(0),KPNEXT(KPNEXT(0)) !对选择的关键点进行连线,即生成底线
KSEL,ALL
CSYS,1
KSEL,,LOC,Y,60 !选择theta=60°的关键点
LSTR,KPNEXT(0).KPNEXT(KPNEXT(0)) !由选择的关键点生成线,即左侧边线
拓扑优化理论及在ANSYS软件中的实现
拓扑优化理论及在ANSYS软件中的实现
拓扑优化理论及在ANSYS软件中的实现
一.拓扑优化概论:
连续体结构的拓扑优化设计是继结构的尺寸优化设计和形状优化设计之后,在结构优化领域出现的一种富有挑战性的研究方向,它是一种比尺寸优化和形状优化更高层次的优化方法,也是结构优化问题中最为复杂的一类问题。拓扑优化处于结构的概念设计阶段,其优化结果是一切后续设计的基础。因而在初始设计阶段需要确定结构的最佳拓扑形式。拓扑优化的目的是寻求结构的刚度在设计空间最佳的分布形式,或在设计域空间需求结构最佳的传力路线,以优化结构的某些性能或减轻结构的重量。
*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
LSEL,ALL
A,1,2,7,6
*REPEAT,4,1,1,1,1 !重复上述命令,共生成4个面
ESIZE,,4
AMESH,ALL
FINISH
/SOLU
NSEL,S,LOC,Y
DSYM,SYMM,X !对选择的节点施加x方向的对称约束
NSEL,S,LOC,X
DSYM,ASYM,Y !对选择的节点施加x方向的反对称约束
拓扑优化的特点是:设计变量多,计算规模大,目标函数和约束函数一般为设计变量的非线性、非单调函数。目前应用于连续体结构拓扑优化计算的优化算法主要包括两类:优化准则法和序列凸规划法。
去除优化过程中数值计算不稳定的方法,优化结果的提取和重构;
拓扑优化中经常出现的数值计算问题有:多孔材料、棋盘格现象、网格依赖性和局部极值问题。优化结果的提取和应用主要考虑的是如何将优化的结果转化为可用的CAD模型问题,实现CAE和CAD之间的数据共享和交流。
L3=LSNEXT(L2)
LFILLET,L1,L2,FIL !对内控的线进行倒角
LFILLET,L2,L3,FIL
LFILLET,L3,L1,FIL
LSEL,AL !选择所有的线
WPROT,,90
LSBW,ALL !所有线由工作平面分割
WPROT,,,60
LSBW,ALL
CSYS,1 !设置坐标系为柱坐标系
NSEL,ALL
FK,10,FX,66725*4 !在10号关键点施加集中载荷,实现弯矩
DK,1,ALL,0 !在1号关键点施加全约束
SOLVE
FINISH
/POST
SET,LAST
ETABLE,EVOL,VOLU !建立单元表,并取出每个单元的体积
PRNSOL,S,PRIN !列出节点的主应力
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
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) !得到下一个编号最小的线的编号
opvar,TK4,9dv,4/H,7/H
opsave,INITIAL,opt !保存所有的优化数据到文件optp里
optype,SUBP !设置优化方式,子问题逼近算法
OPSUBP,30 !指定迭代次数
opexe ! 执行优化
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
TVOL=TVOL*2 !由于分析时只计算了结构的一半,总体积要乘2

下载文档原格式

  / 5