拓扑优化学习报告_北理工_王路.pdf

《拓扑优化学习报告_北理工_王路.pdf》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、北京理工大学车辆工程王路基于99行程序的拓扑优化学习报告（一）背景和前言随着汽车工业的飞速发展以及日益突出的能源问题，汽车工业面临的挑战以及竞争环境也越来越激烈，对汽车产品提出了降低其制造成本及燃油经济性的新要求。在提高汽车安全性、减少汽车排放和解决能源消耗的背景下，提出了汽车轻量化技术。实现汽车轻量化的途径包括三个方面：结构优化技术、新型材料和先进性制造工艺。其中，我们所讨论的是结构优化技术，其中结构优化设计分为三个层次：尺寸优化（SizeOptimization）、形状优化（ShapeOptimization）和拓扑优化（TopologyOptimiz

2、ation）。本文我们基于99行matlab程序初步学习拓扑优化技术中的理论和优化方法。拓扑优化技术指的是在给定的设计空间内寻求最佳的材料布局，同时在满足平衡方程、物理关系、几何关系和边界约束条件下使得结构达到某种性能最优的应用技术。拓扑优化的理论研究最早可以追溯到Michel提出的桁架理论，连续体结构的拓扑优化由于描述和数值计算得困难，发展一直相对缓慢，直到Bendsoe和Kikuchi在1988年提出的均匀化方法之后才得到迅速的发展，其基本思想是在组成拓扑结构中引入微结构，通过微结构的几何参数作为设计变量，通过微结构的增加和删减实现结构的拓扑形状的改变

3、，实现拓扑优化和尺寸优化的统一。在微结构的基础上，我们介绍变密度法的应用，变密度法是在均匀化方法的基础上产生的，把材料引入微结构代之以密度在0~1之间变化的假想材料，把密度作为设计变量，从而实现材料的删减，因其模型简单、计算变量相对较少成为目前广泛采用的方法。根据不同的插值模式，变密度法又有不同的插值模型：SIMP法（SolidIsotropicMaterialwithPenalization）、Hashin-Shtrikman法，以及RAMP法（RationalApproximationofMaterialProperties）。（二）拓扑优化问题的描述

4、（1）材料插值模型变密度理论的材料插值模型通过引入中间密度单元，将离散型问题转化成连续型优化问题，而实际上，中间密度单元是无法存在和制造，因此要尽量避免中间密度单元的产生，减少中间密度单元的数目，此时就需要对设计变量中出现的中间密度值进行惩罚。北京理工大学车辆工程王路目前材料插值模型主要分为两类，SIMP法和RAMP法，基于SIMP格式的材料插值模型为：pE(x)E(x)(EE),x0,1（1）imini0mini其中：E(x)——插值以后的弹性模量iE——实体部分材料的弹性模量0E——孔洞部分材料的弹性模量minx——单元相对密度，取值为1时

5、表示有材料，为0时表示无材料即孔洞ip——惩罚因子在本文中，在本篇文章中，讨论的SIMP的表达为：pE(x)(x)E（2）i,ji,j0其中：xi,j为第i个子域内第j个单元的相对密度。（2）数学模型根据在体积或质量约束下求最小柔度，即最大刚度来建立优化模型，基于变密度理论的SIMP法的周期性拓扑优化问题的数学模型可以表达为：TfindX(x1,1,x1,2,x1,3,,xi,j)Ri1,2,,mj1,2,,nmnmnTTpTminC(X)FUUKUui,jki,jui,j(xi,j)ui,jk0ui,ji1j

6、1i1j1s.tKUFmnVfV0xi,ji,ji1j10xminxi,jxmax1其中：X——单元相对密度矢量C——结构的柔度F——载荷矢量U——位移矢量北京理工大学车辆工程王路k——结构刚度矩阵ui,j——单元位移矢量ki,j——单元刚度矩阵k——初始段元刚度矩阵0——单元体积i,jV——优化后的体积f——保留的体积分数V——初始体积0x——设计变量的取值下限minx——设计变量的取值上限maxn——子域内单元的个数（三）程序解读99行程序主要包括主函数、OC优化准则、网格过滤、有限元分析、灵敏度分析5个部分，其流程图

7、为：北京理工大学车辆工程王路设计域的离散化和程序参数的初始化主程序调用有限元分析求解子程序位移U灵敏度分析子程序dc否网格过滤子程序dc(new)OC优化准则优化设计变量x(new)是否满足精度要求是结束循环图1.99行程序的流程图具体程序分析：总体程序为：%%%%A99LINETOPOLOGYOPTIMIZATIONCODEBYOLESIGMUND,JANUARY2000%%%%%%%CODEMODIFIEDFORINCREASEDSPEED,September2002,BYOLESIGMUND%%%functiontop(nelx,nely,volfr

8、ac,penal,rmin);nelx=80;nely=20;vo