拓扑优化算法及其实现111_第1页
拓扑优化算法及其实现111_第2页
拓扑优化算法及其实现111_第3页
拓扑优化算法及其实现111_第4页
拓扑优化算法及其实现111_第5页
已阅读5页,还剩30页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、拓扑优化简介拓扑优化设计流程算例拓扑优化:拓扑优化是结构优化的一种。结构优化可分为尺寸优化、形状优化、形貌优化和拓扑优化。尺寸优化以参数为优化对象,比如板厚、梁的截面宽、长和厚等;形状优化以结构件外形或者孔洞形状为优化对象,比如凸台过渡倒角的形状等;形貌优化是在已有薄板上寻找新的凸台分布,提高局部刚度;拓扑优化以材料分布为优化对象,通过拓扑优化,可以在均匀分布材料的设计空间中找到最佳的分布方案。拓扑优化相对于尺寸优化和形状优化,具有更多的设计自由度,能够获得更大的设计空间,是结构优化最具发展前景的一个方面。 优化设计过程:将区域离散成足够多的子区域,对这些子区域进行结构分析,再按某种优化策略和

2、准则从这些子区域中删除某些单元,用保留下来的单元描述结构的最优拓扑。变密度法Level Set法 (水平集法)ICM(独立映射法)ESO(进化法) 拓扑优化建模方法SIMP( Solid Isotropic Microstructures with Penalization )(固体各向同性惩罚函数法)RAMP( Rational Approximation of Material Properties )(材料属性的理性近似模型)OC法(优化准则法)MMA法(移动渐进线法)SLP(序列线性规划法)SQP(序列二次规划法)优化求解方法拓扑优化简介算例基于99行拓扑优化程序代码top(nelx,

3、nely,volfrac,penal,rmin)312SIMP法+OC法4节点矩形单元TeeeKB DB tdA e e1234目的:求整体位移矩阵B、D、t分别代表什么?eKK单元刚度矩阵整体刚度矩阵KUFUF和 为整体位移与载荷 e ee划分网格数(nelx,nely)整体节点编排:纵向1nely+1nely+22(nely+1)12875643局部整体 e e(1)(2)(3)(4)KUF(有限元基本方程)U 各节点位移矩阵建立优化模型(min&max)目标函数约束函数设计变量10minmin(). .( )01nTepTeoeeCU Fu k uKUFst VV( )()epx

4、设计变量e 本质上,结构拓扑优化是个( 0,10,1)整数规划问题,属于组合优化的范畴。 2 2n n次计算有限元分析才能求得全局最优解,是个指数时间算法,非多项式时间算法,随着单元数量n n的增加,计算量会激增,也即是困扰组合优 化 领 域 的 N PN P 难 题 。10minmin(). .( )01nTepTeoeeCU Fu k uKUFst VV( )()epx设计变量e密度变量的引入: 在工程中,材料的刚度线性依赖材料的密度,即刚度大的材料,密度也大。比如,钢的密度比铝的密度大,因此钢的刚度比铝的刚度大。按照这个朴素的逻辑,就可用单元的密度来代替材料的有无,如下式 peeKKeK

5、 为惩罚的单元刚度矩阵;eK 为真实的单元刚度矩阵;0,1为每个单元的密度,是 到 之间的连续变量;01()p为总大于 的罚因子,一般取 ,大的罚因子可使密度快速的趋于 或者 。1Penalization301(0,1)整数规划问题0,1区间内的单元密度的连续变量优化问题OC法优化求解优化求解变量约束优化的变量约束优化的KTKT条件:条件:定义在闭区间上的一元函数的优化(仅有变量上下限约束)对应于目标函数的拉格朗日函数为:对应于目标函数的拉格朗日函数为:u优化设计准则优化设计准则优化设计准则优化设计准则 上式即为设计变量的迭代准则。由该式可以看出当柔度取得极值时,在整个设计区域内单元的应变能密

6、度是恒定的常数 。由此建立更新设计变量迭代格式:minmin1minmax(,)max(,)max(,)min(1,)min(1,)min(1,)nnneeennnnneeeennneeemifBmBifmBmmifBm101(),pTneeeneepu k uBBvminmin1minmax(,)max(,)max(,)min(1,)min(1,)min(1,)nnneeennnnneeeennneeemifBmBifmBmmifBm1()nTepTeeeeCU Fu k uImagesc(-x)xnew = max(0.001,max(x-move,min(1.,min(x+move,x.

7、*sqrt(-dc./lmid)优化结果:各单元密度组成的矩阵X目的:消除棋盘格效应及网格依赖性目的:消除棋盘格效应及网格依赖性方法:方法:1、高阶单元法(计算量大) 2、周长约束法(周长约束的上限值需要依靠经验来确定,因为局部尺寸和周长边界间没有直接的关系。如果周长约束边界定得太紧,则可能导致没有计算结果,如果定的太松又达不到预期的效果。因此约束边界很难确定,这种情况在三维问题下特别明显) 3、局部梯度约束方法(局部斜率约束属于局部约束,可防止局部细条的形成,从而降低结构拓扑的几何复杂性,但优化结果难以满足全局最佳,并且,这种方法在优化问题中引入了2N(二维)或3N(三维)个额外约束,使计算

8、效率大大降低) 4、网格过滤法(网格过滤方法只需定义一个局部长度尺寸,相对、网格过滤法(网格过滤方法只需定义一个局部长度尺寸,相对较为容易,在约束尺度下的结构变量都被过滤掉。网格过滤方法的优点是较为容易,在约束尺度下的结构变量都被过滤掉。网格过滤方法的优点是不需要在优化问题中加入额外约束,且容易实施。缺点是过滤方法为一种不需要在优化问题中加入额外约束,且容易实施。缺点是过滤方法为一种基于启发式求解规则的方法)基于启发式求解规则的方法)网格过滤法原理:数字图像处理中的降噪技术网格过滤法原理:数字图像处理中的降噪技术 一张图片被离散成有限个像素点,一个像素点代表一个灰度值(256个灰度水平)。假设

9、灰度值是连续变化的。 图像处理的一个常见问题是噪音来自电气传感器噪声,噪声传输错误等。主要表现为离散孤立的像素变化,这些地方与其周围通常出现显著的不同。 解决方法:线性和非线性解决方法:线性和非线性 Pratt,(1991) 证明了非线性技术优于线性技术。但是非线性技术不能应用于拓扑优化,因为非线性方法常会造成拓扑优化问题的不光滑。 线性技术又分为两种:基于傅里叶变换对的方法(频率)和基于卷积的方法。但是傅里叶变换法只能用于常规的矩形网格。 目前 SIMP 方法常用的过滤技术是Sigmund提出的过滤方案,用于修改目标函数的敏度信息:拓扑优化简介OC法拓扑优化设计流程左边界各节点受横向约束右下角节点受纵向约束F(2,1) = -1;fixeddofs = union(1:2:2*(nely+1),2*(nelx+1)*(nely+1);top(60,20,0.5,3,3)6020在Matlab中

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论