基于MATLAB和ANSYS的悬臂梁拓扑优化

现代设计理论和方法大作业----基于MATLAB和ANSYS的悬臂梁拓扑优化指导老师：刘志刚项目组长：薛亚波项目成员：机自66学院：机械工程学院基于MATLAB和ANSYS的悬臂梁拓扑优化一、计划和任务安排表1.1学习时间安排表：任务时间查阅整理资料5月03日——5月08日学习相关知识和论文5月09日——5月15日设计实验方案5月16日——5月17日学习优化理论5月18日——5月20日撰写论文5月21日——5月22日集体交流学习心得5月23日——5月24日1.2任务分配表：学号姓名任务06011136齐寰宇查阅资料、答辩PPT制作06011137王晓庆查阅资料、答辩PPT制作06011139韩冰查阅资料、研究问题分析06011141康辰龙查阅资料、ansys拓扑模块研究06011143梁伟云查阅资料、项目背景研究06011145卢岗查阅资料ansys拓扑模块研究06011153谢清伟查阅资料ansys拓扑模块研究06011154薛亚波查阅资料、后期工作总结06011155严宏伟查阅资料ansys拓扑模块研究06011156杨志亮查阅资料、项目背景研究06011157张龙查阅资料、答辩PPT制作06012020乔百杰查阅资料、Matlab结构优化06073012邓磊查阅资料、项目背景研究06182047魏立峰查阅资料、研究问题分析06055164杜博查阅资料、研究问题分析06011146牛玉辉查阅资料、研究问题分析06011147朴群星查阅资料、项目背景研究二、项目背景介绍及问题描述2.1项目背景及意义：2.1.1工程背景及基本原理：通常把结构优化按设计变量的类型划分成三个层次：结构尺寸优化、形状优化和拓扑优化。尺寸优化和形状优化已得到充分的发展，但它们存在着不能变更拓扑结构的缺陷。在这样的背景下，人们开始研究拓扑优化。拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料的分布问题。寻求一个最佳的拓扑结构形式有两种基本的原理：一种是退化原理，另一种是进化原理。退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上，然后构造适当的优化模型，通过一定的优化方法逐步删减那些不必要的结构元素，直至最终得到一个最优化的拓扑结构形式。进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化，它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。2.1.2本文研究意义：目前，结构优化大部分集中在尺寸设计变量(如板厚、杆的剖面积及管梁的直径)。拓扑结构优化较尺寸优化复杂,但对于有些问题拓扑结构优化比尺寸优化有效,悬臂梁是其中的例子之一。本文讨论悬臂梁的拓扑优化问题,围绕这一问题,怎样使结构具有最大刚度的设计占有相当重要的地位;怎样优化结构的形状使材料的分布，更加合理从而达到使结构具有最大刚度的目的是本文要研究的问题。2.2研究现状2.2.1理论研究现状：结构拓扑优化是近20年来从结构优化研究中派生出来的新分支，它在计算结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化的工程应用研究还很不成熟，在国外处在发展的初期，尤其在国内尚属于起步阶段。1904年Michell在桁架理论中首次提出了拓扑优化的概念。自1964年Dorn等人提出基结构法，将数值方法引入拓扑优化领域，拓扑优化研究开始活跃。20世纪80年代初，程耿东和N.Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题，他们开创性的工作引起了众多学者的研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计，开创了连续体结构拓扑优化设计研究的新局面。1993年XieYM和StevenGP提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002年罗鹰等提出三角网格进化法，该方法在优化过程中实现了退化和进化的统一，提高了优化效率。2.2.2应用研究现状：在前人提出的重要理论基础上，后人也将其跟其他现代设计的方法相结合，衍生出了其他一些拓扑结构优化方法：如与可靠性相结合的情况下，MAUTE等应用变密度法并结合可靠性分析对一微机电系统进行了基于可靠性的拓扑优化设计，PAPADRAKAKIS等将遗传算法应用于具有可靠性约束的桁架结构拓扑优化设计中，国内学者马洪波也对基于遗传算法的结构可靠性优化问题进行了讨论。华南理工大学机械工程学院欧阳高飞等对基于水平集方法的结构可靠性拓扑优化进行了研究。2.3研究目标：2.3.1设计目标：设计一悬臂梁在自由端受到竖直向下集中力P的作用，要求保持悬臂梁长度尺寸不变。2.3.2优化目标：在满足以上条件情况下对原实体悬臂梁做结构拓扑优化设计，使其结构刚度最大、重量最轻。（优化的结果应该使原设计区域产生孔洞，使结构拓扑发生变化。）优化目标分析：要使梁满足以上条件情况必须同时兼顾重量和刚度问题。重量太轻不能满足刚度问题。刚度太大又必须足够的材料以满足其要求。将两个因素进行耦合分析，我们可以得出的结论是必须通过结构层次来进行优化。连续体结构拓扑优化较成熟的方法主要有均匀化方法、变密度方法和渐进结构优化法，通过综合比较我们最终选用变密度法。2.3.3实现方法：本文在对题目进行分析后，决定从两个独立的方向来分别研究，以此来互相验证，保证结果的可靠性。实验中要使用ansys和matlab两种软件分别独立进行优化分析，然后将两种结果来进行比较。2.3.4实践目标：通过这次的作业我们期望通过搜索资料，团体讨论，分组作业的方式，以期望达到对拓扑优化的方法和过程有个具体的了解和学习。另外对分析软件的应用能够达到一个新的高度。这些不仅能使我们现在的知识体系得到充实和优化，而且也是我们今后人生的财富2.4研究内容2.4.1问题描述如图所示，悬臂梁在自由端受到竖直向下集中力P的作用，要求在保持悬臂梁长度尺寸不变的情况下对原实体悬臂梁做结构拓扑优化设计，优化目标是使结构刚度最大、重量最轻。（优化的结果应该使原设计区域产生孔洞，使结构拓扑发生变化。）原实体悬臂梁横截面为矩形，矩形宽度为0.2m，悬臂梁的材料为45钢，密度ρ为7.8×10-6kg/mm3，弹性模量E=2×105MPa，许用切应力[τ]=60MPa。竖直向下集中力为N1045.1P。图悬臂梁受集中力作用（单位：m）2.4.2问题转化本研究中，要求在保持悬臂梁在结构尺寸不变的情况下，使悬臂梁的刚度最大、质量最轻。一般悬臂梁都为三维立体结构，载荷是在竖直方向上加载的，所以在垂直纸面方向上，并且在悬臂梁的对称平面上加载，根据材料力学知识可知，所加载在和只对加载区周围较小范围内产生影响，其余部分与加载形式无关。其受力与如图所示平面的手里相同，故可以将问题抽象为平面问题来处理。当占有区域S的一个结构上作用有强度为P(S)的载荷时,使结构具有最大刚度的问题和使结构所受外力功W具有最小值的问题是等价的。结构所受的外力功W与结构变形能C之间的关系为W=2C，所以我们可以将该问题处理为求在一定的约束条件下，该悬臂梁的的最小变形能最小的问题。2.4.3数学模型：该问题中，要求同时满足刚度最大，质量最轻，这两个变量若同时改变，则问题复杂度太大，并亲切可能导致问题不可求解。所以我们采用在确定的质量下，来讨论刚度最大的问题。由于对特定的材料，其质量和体积有一定的关系，并且我们采用去除法的思想来建立模型的，故我们可以采用给优化后的体积与优化前的体积比赋确定的值，来达到在给定质量条件下满足刚度最大的问题。其数学模型如下：注：其中C（x）为结构变形能，U为结构变形总位移矩阵，K为结构总刚度矩阵，N为划分单元总数，eu为单元位移向量，ok为单元刚度，（由于划分单元的时候，我们采用等分举行单元，所以个单元的刚度可用一个常量来处理）V(x)是拓扑结构优化过程中变化着的体积，V0为未经过优化前悬臂梁的体积。F为结构所受的总载荷。x为悬臂梁的相对密度。2.4.3模型分析求解：该问题的优化方法有很多种，常用的有如下方法：OptimalityCriteria(OC)methods,（优化准则方法）SequentialLinearProgramming(SLP)methods（序列线性规划法）MethodofMovingAsymptotes(MMAbySvanberg1987)等为了简化问题的复杂度，此处我们采用standardOC-method.方法来实现。在处理过程中，关于设计变量相对密度x每一步的更新，我们采用在1995年提出的如下算法来实现：注：其中m是移动限制量，即相对密度变化的最大步长，更新条件中，加上最小边界条件判断，是防止相对密度变为0，使问题出现奇异；加上相对密度最大边界条件判断，是为了防止出现相对密度大于1，是问题是去物理意义。是数值阻尼系数。Be可以由以下条件来确定：其中为拉格朗日乘数，由双向切片算法来确定。目标函数对单元相对密度的变化率，即单元敏感度如下确定：为了确保该问题的拓扑结构优化解的存在，在处理技术上，我们得对该问题的求解过程加些限制。此处我们采用过滤技术来进行处理。我们采用网格独立性滤波器来对单元敏感度的计算来进行改进，其改进算法如下：其中Hf为卷积因子（即权重因子），其由以下表达式来确定：采用该算法对单元敏感度的更新进行改进，是为了在更新单元e的敏感度时，将距其距离不超过rmin的临域内的单元来进行均匀化，防止各乡邻单元的敏感度变化过大。注：其中dist(e,f)定义为单元体e的中心到单元体f中心之间的距离，rmin为过滤大小。在优化过程中，我们采用改进的敏感度表达式（5）来代替表达式（4）来对（3）是进行计算。三、实验设计及结果分析：3.1技术路线：在本实验中，由于我们提出了两种独立的方法，故我们将小组分为两部分;小组一：该组成员采用Ansys软件来对该悬臂梁模型进行建模分析，并进行拓扑结构优化。小组二：该组成员采用上面的数学模型，用Matlab原来来实现整个数学模型的算法。并对处理结果进行动态显示，以便让我们清晰地看到整个拓扑结构优化的过程。3.2试验结果及分析：3.2.1Ansys分析的拓扑结构优化结果：3.2.1.2建立几何模型：3.2.1.2.1几何参数设置如下：3.2.1.2.1在Ansys/workbench中建立如下几何模型：3.2.1.2网格化：将几何模型直接调入有限元环境。在workbench里建立的模型可以无缝导入到有限元分析环境里。采用自动网格划分工具，得到的模型如下：其中，网格单元和节点数设置如下：StatisticsNodes节点数目3120Elements单元数目5613.2.1.3施加载荷和设置边界条件选择优化类型为“shapeoptimization”——形状优化在上截面施加15000N的力3.2.1.4求解：3.2.1.4.1设置材料属性：StructuralYoung'sModulus杨氏模量2.e+011PaPoisson'sRatio0.3Density密度7850.kg/m³ThermalExpansion1.2e-0051/°CTensileYieldStrength2.5e+008PaCompressiveYieldStrength2.5e+008PaTensileUltimateStrength4.6e+008PaCompressiveUltimateStrength0.PaThermalThermalConductivity60.5W/m·°CSpecificHeat434.J/kg·°CElectromagneticsRelativePermeability10000Resistivity1.7e-007Ohm·m3.2.1.4.2定义期望目标值的减少量，即期望减少的质量或体积分数：ScopeGeometryAllBodiesDefinitionTargetReduction50.%3.2.1.4.3拓扑结果如下：3.2.1.4.4其他两种结果：定义质量减少率为80%ScopeGeometryA