基于拓扑优化的车身结构研究
瞿元王洪斌张林波吴沈荣
奇瑞汽车股份有限公司，安徽芜湖，241009
摘要：随着CAE技术的发展，虚拟仿真技术在汽车开发中的作用也愈来愈显著。而前期工程阶段，如何布置出合理的车身骨架架构，一直是个相对空白的地带，也是整车正向开发过程中绕不过的坎。尽管研发工程师根据经验，参照现有车型的结构特点，也能进行车身骨架架构的设定，但总是缺乏有效手段直观地反映不同车型结构布置的特点。本文用拓扑优化的方法，从结构基本特征的角度来审视这一问题，并运用该方法对某SUV车身结构进行研究，获得一些直观性的结论。
关键词：车身,前期工程,拓扑优化
1引言
随着对整车研发过程认识的加深，以及对正向开发过程的探索，在车型开发前期，对车身结构做出更合理的规划显得愈来愈重要。常规的研发思路之一是通过参考已有车型的结构，经过适当的修改，形成新的结构，并用于新车型中。但是对于原始车型的设计思路、结构布置的原因等缺乏系统的理解，或者理解不深，往往在更改过程中产生新的问题。为了部分解决上述问题，本文从结构拓扑优化的角度，对某SUV 车型车身结构的总体布置进行初步探讨，以期加深对结构布置的理解。
2研究方法概述
合理化的车身结构，是满足整车基本性能的重要保障。为了能够实现结构的最优布置，文献[1]使用了拓扑优化工具来布置车身结构。其基本思路是从造型以及车内空间布置出发，建立车身空间的基础网格模型，然后根据一定的工况要求，对基础网格进行拓扑分析，并根据拓扑结果建立梁、板壳模型，并进行多项性能的优化，从而实现车身结构的正向开发。本文借助于该思想，建立研究对象的结构空间包络，并对该包络进行拓扑分析，然后将仿真结果与原始结构进行比较，寻找车身结构中的关键点，推测初始结构可能的布置思想，从而加深对该研究思路的理解。其基本过程如下图所示：
图1 研究思路
3某SUV车身结构研究
本文选取的研究对象是某比较受欢迎的城市SUV，通过对其结构的研究，有利于了解其总体结构的布置原理。本文采用上文提到的方法，对该车型结构进行剖析。
3.1 研究对象及结构包络
图2显示为该SUV车型的结构模型，经过适当处理，产生图3的结构包络模型，可以用该结构包络进行后面的结构拓扑分析。该结构包络模型完全包含了原结构的布置空间，材料为普通钢材。考虑到整车在前后碰撞过程中需要纵向刚度比较好，将前后纵梁以及后轮罩直接设定为非设计空间，其余网格均作为设计空间。
图2 某SUV结构有限元模型
图3 某SUV结构包络模型
3.2 工况
车身在实际使用过程中承受非常复杂的载荷，这些载荷对车身的影响各不相同，有的影响局部，有的影响整个车身。在实际研发过程中，不可能对所有可能的工况进行考察，而且，不同的设计阶段，考察的指标也不相同。在概念设计阶段，更重要的是保证车身的总体结构刚度，避免后期产生较大变更，导致项目延期或者增加较多的开发成本。本文主要考察某SUV车型结构布置特点，因此，主要考虑NVH以及碰撞两个方面的工况。其具体考察工况如下表1所示，4个NVH工况，主要考察整体刚度以及前后端的弯曲性能；4个碰撞方面的工况，主要考察车身承受不同方向的撞击。
表1 主要考察工况[1]
对于NVH类的工况，其导致的车身变形都是非常小的，所以，在拓扑优化过程中，不需要考虑结构失效问题。而Safety工况对车身的影响都是大变形、非线性的，还有接触力存在。考虑到前后纵梁作为非设计空间，可以看做是刚度很硬的结构，这样对车身其他部位而言，Safety工况的影响将限定在线性范围内。因此，在进行车身拓扑优化的时候，上述工况均作为小变形来处理。
为了更好地反映不同工况对车身结构布置的影响，对以上8种工况组合成4种研究方案，分别为（1）Safety工况单独考察；（2）NVH工况单独考察；（3）NVH和Safety同时考察；（4）先Safety工况然后NVH工况。如表1所示，8种工况中，每个工况在不同考察方案中的权重比率一致，比如对第一种方案，表1中5-8号工况的权重比例为：15:10:10:5，具体权重根据上述比例分配来设定，其余类同。
3.3 优化目标与约束
对于一个优化问题而言，设计变量、约束与优化目标是其主要元素。对于本文的优化问题而言，其设计变量为单元密度，约束分别是体积分数和对称约束，而优化的目标为考察工况下车身结构的柔度最小化，其中，体积分数是指剔除初始非设计体积的当前总体积与初始总体积的比值，而对称约束是指相对于车身XZ平面，车身结构左右对称。以上优化问题可简化描述成如下形式[2]：
目标：最小化
，i表示设计变量数，j表示考察工况数，
为第j个工况下的柔度，为第j个工况的权重；
约束：体积分数，车身对称约束；
设计变量：，其中表示第i个单元的单元密度。
3.4优化结果与结构分析
为了能够获取上述不同方案下的结构拓扑布置，采用OPTISTRUCT软件对上述结构进行优化计算[3]，并将结果分别与该SUV的结构布置进行比较。对于车身结构而言，一般可分成前舱、顶盖、地板、侧围、后围等结构部分。考虑到部分结构模块的复杂性以及碰撞分析的非线性，所以在下面对比过程中，主要考察对象为顶盖、下部车身、侧围、防火墙这四个部分。
3.4.1 顶盖
对于车顶盖结构而言，四种方案中，方案2-4涉及NVH工况，且扭转工况所占比重较大，为了获得较高的抗扭转能力，主要结构成斜网状分布，结果如图4所示。方案1为安全工况，载荷以纵、横分布为主，所以顶盖材料主要呈横向（Y向）分布，顶部横梁分布于B柱后侧③以及后背门上部②。方案2为完全NVH工况，除了前顶横梁①以及后部横梁外②，中间梁均为网状，而前后顶横梁①②对于车身抗扭能力具有重要作用。方案3中同时考虑了安全性能与NVH要求，前后横梁依然是重要结构，同时在B柱上部后侧同方案1一样，存在一根横向布置梁③。方案4中，将方案1（安全工况）中的主要横梁设定为
非设计空间，然后进行NVH工况优化，从结果看，顶盖材料分布与方案3相近。从4个方案的结果可以看出，对于车身顶盖而言，前后横梁①②以及B柱上部横梁③对车身主要结构的贡献相对比较大，在结构布置的时候，需要重点关注。其结果同原结构比较来看，前后横梁相互对应，其他中间横梁对整车刚度的作用比较小，其主要作用表现为对顶盖外钣金的支撑以及抗雪压等。
①前顶横梁；②后顶横梁；③B柱上横梁；
图 4 顶盖结构拓扑图
3.4.2 下部车身
下部车身共考察9个部位梁的分布，如下图5所示。方案1主要体现了①②④⑧⑨五个横梁，横梁⑤⑥较原始结构有偏移，拆分成了3个横梁。从SUV的地板纵梁布置看，其结构形式明显不同于一般贯通式布置，其纵梁②与整车X向呈一定角度，且中止于前地板中横梁⑤，该纵梁的布置在方案1中得到较显著提示，其拓扑材料的分布与SUV的该纵梁分布相同，同样止于前地板中横梁⑤附近。后地板前横梁⑦对碰撞安全的影响不明显，拓扑结果中在该位置没有材料布置，而横梁⑧与⑨的作用显著。
方案2中主要考察NVH性能，从结果看，前地板材料以中通道为中心连接梁，呈网状分布，说明该部位的Y向横梁的作用相对弱化，只有如结果所示的布置才能最大化扭转刚度；在后地板部位，横梁⑧⑨位置清晰，其中横梁⑧前部材料呈X状分布，前点部位对应原型结构的横梁⑦，后点部位连接横梁⑧，说明这两个横梁对刚度的贡献比较大，尽管横梁⑦与所给的拓扑材料分布有差别。该布置方式也符合一般的认知，从车身实际情况看，该部位也是设计重点关注部位之一，车身整体扭转刚度的提升、路试开裂问题的解决，都需要在这里开展更多的工作。
①防火墙；②前地板侧纵梁；③前地板；
④前地板前横梁；⑤前地板中横梁；⑥前地板后横梁；
⑦后地板前横梁；⑧后地板中横梁；⑨后地板后横梁；
图 5 下部车身结构拓扑图
方案3综合考虑的Safety与NVH工况，从结果可以看出，梁①-⑤、⑧⑨均可以较清晰地体现，梁⑥在方案中有所体现，不过位置较原型结构而言，偏靠前一些。而梁⑦与纵梁搭接处形成了X型交叉结构，与方案2相同，说明，该处位置主要体现了刚度方面的要求。对于前排驾乘舱地板部位而言，交叉型的梁布置结构有利于提升整体的扭转刚度，对比到原型结构可以看出，原型结构中前地板纵梁布置不同于一般的垂直布置，呈现一定角度的交错分布结构，同时兼顾了方案1中的前地板纵梁的分布，从拓扑结果来看，该SUV的前地板梁的布置考虑比较综合，兼顾多个性能。
方案4将方案1中部分梁结构设定为非设计空间，然后按照NVH工况进行优化分析。从结果看，①和⑨是必要的整体结构，②在该方案中也有一定体现，但是梁的长度相对较短，截止于梁④。方案中梁③较原结构有变化，斜向连接于①④之间。后地板部位梁的布置同方案2和3，均会形成X结构。
从上述方案来看，SUV的地板结构布置考虑性能比较综合，对于该尺寸的车身结构布置而言，前后纵梁的不贯通设计也符合拓扑结构的布置。梁④和⑦、⑧是地板横梁中非常重要的结构，梁②-③、⑤-⑥的具体布置，需要根据实际需要，进行适当的调整，保证车身Safety和NVH的性能综合平衡。
3.4.3侧围
对于侧围结构而言，从下图6比较看，4种组合方案对侧围结构的影响主要体现在车身后部三角窗处，方案1（安全工况）的布置与其他三种方案略有差别，材料分布上只有三角窗前、后梁①②以及与轮罩的连接件④。其他方案中，除了原型中的梁①-④的布置比较清晰，梁③下部的材料也有所布置，从而可以提供比较好的NVH性能，这在方案1中没有得到反映，说明这些材料分布对NVH性能而言，具有更多正向作用。
①三角窗前梁；②三角窗后梁；③三角窗下梁；
④三角窗与轮罩连接件；
图 6 侧围结构拓扑图
3.4.4 防火墙
对于防火墙结构而言，主要承担碰撞以及扭转载荷。从图7所示拓扑结果来看，防火墙的材料呈网状分布，其中方案1能比较清晰地看到梁①-③的分布，其他方案中，只有梁②③比较清晰可见。
①防火墙横梁；②纵梁；③流水槽横梁；
图7防火墙结构拓扑图
4 总结
通过材料布置的拓扑分析，可以看出，该SUV车身原结构在布置上比较符合较佳的材料分布。初步获得下列直观性结论：
（1）对车顶盖结构而言，前后顶部横梁对车身具有比较重要的影响，尤其是后顶部横梁；B柱上部顶横梁对顶部Safety具有重要影响；
（2）对于下部车身而言，除防火墙部位以及后围横梁外，驾乘舱下部的座椅横梁以及后轮罩部位的横梁，无论是Safety工况还是NVH工况，都比较重要；
（3）地板纵梁的布置可以是非贯通式的，其走势上也并非一定要与整车X向一致，可以与门槛梁的走向成一定角度；地板梁一定程度的网状布置，对提升抗扭性能有正向作用；
（4）后侧围三角窗部位结构对安全和NVH都比较重要；
（5）对于防火墙布置而言，网状的结构对NVH性能有正向作用。
本文所述方法在有效性方面还需要更多的研究例证。而且，对于Safety工况进行了线性化处理，所得结果是否符合碰撞安全的设计理念还需要进一步研究。
5 参考文献
[1] Carl Reed, Application of Optistruct optimization to body in white design, Jaguar Cars L imited, 2002
[2] 周传月等主编，工程有限元与优化分析应用实例教程，科学出版社，2005.6
[3] Altair, HyperMesh User’s Guide, 2009
Body Structure Investigation with Topology O
ptimization
Qu Yuan Wang Hongbing Zhang Linbo Wu Shenrong
Abstract: With development of CAE technology, virtual simulation becomes more and more import ant. During the early stage of vehicle design, there is lack of information for how to layout an id eal body framework. To accomplish such a task is also critical for vehicles program development. Although, by experiences, engineers may be able to design body framework based on structure features of existing vehicles, there is lack of effective and intuitive methods to describe the struct ural characteristics for various vehicles. In this paper, the topology optimization is performed to in vestigate the body structure of some SUV. The numerical results are compared to the vehicle str ucture. Certain similar structural features are found.
Key words: Body Early Engineering Topology
基于拓扑优化的车身结构研究 瞿元王洪斌张林波吴沈荣 奇瑞汽车股份有限公司，安徽芜湖，241009 摘要：随着CAE技术的发展，虚拟仿真技术在汽车开发中的作用也愈来愈显著。而前期工程阶段，如何布置出合理的车身骨架架构，一直是个相对空白的地带，也是整车正向开发过程中绕不过的坎。尽管研发工程师根据经验，参照现有车型的结构特点，也能进行车身骨架架构的设定，但总是缺乏有效手段直观地反映不同车型结构布置的特点。本文用拓扑优化的方法，从结构基本特征的角度来审视这一问题，并运用该方法对某SUV车身结构进行研究，获得一些直观性的结论。 关键词：车身,前期工程,拓扑优化 1引言 随着对整车研发过程认识的加深，以及对正向开发过程的探索，在车型开发前期，对车身结构做出更合理的规划显得愈来愈重要。常规的研发思路之一是通过参考已有车型的结构，经过适当的修改，形成新的结构，并用于新车型中。但是对于原始车型的设计思路、结构布置的原因等缺乏系统的理解，或者理解不深，往往在更改过程中产生新的问题。为了部分解决上述问题，本文从结构拓扑优化的角度，对某SUV 车型车身结构的总体布置进行初步探讨，以期加深对结构布置的理解。 2研究方法概述 合理化的车身结构，是满足整车基本性能的重要保障。为了能够实现结构的最优布置，文献[1]使用了拓扑优化工具来布置车身结构。其基本思路是从造型以及车内空间布置出发，建立车身空间的基础网格模型，然后根据一定的工况要求，对基础网格进行拓扑分析，并根据拓扑结果建立梁、板壳模型，并进行多项性能的优化，从而实现车身结构的正向开发。本文借助于该思想，建立研究对象的结构空间包络，并对该包络进行拓扑分析，然后将仿真结果与原始结构进行比较，寻找车身结构中的关键点，推测初始结构可能的布置思想，从而加深对该研究思路的理解。其基本过程如下图所示：
多工况应力约束下连续体结构拓扑优化设计ΞTOPOLOG Y OPTIMIZATION DESIGN OF THE CONTINUUM STRUCTURE FOR MU L TIPL E LOADING CON DITIONS WITH STRESS CONSTRAINTS 王 健ΞΞ (山东理工大学交通与车辆工程学院,淄博255012) 程耿东 (大连理工大学工程力学研究所,大连116024) WAN G Jian (Traffic and Vehicle Engineering School,Shandong Univer sity o f Technology,Zibo255012,China) CHEN G Gengdong (Research Institute o f Engineering Mechanics,Dalian Univer sity o f Technology,Dalian116024,China) 摘要 建立多工况应力约束条件下连续体结构拓扑优化的数学模型,给出求解方法。采用包络法处理大量的应力约束,用改进的满应力法进行求解,方法简单、实用。提出的分层优化技术能使最优结构更为清晰。分层优化方法的基本思想是按载荷大小分为几个层次,后面层次的拓扑优化以前面层次得到的最优拓扑为基础,通过逐层优化,最终得到最优结构。分层优化时主要考虑属于本层载荷的影响,避免大小载荷混在一起,最优拓扑模糊不清的问题。为解决各层优化单元厚度相差太大,易造成结构刚度矩阵奇异的问题,提出对相应参数的调整方法。算例表明该方法是有效的。 关键词 结构拓扑优化 应力约束 连续体结构 满应力法 分层优化技术 中图分类号 T B114.3 T B115 Abstract The mathematical m odel of topology optimization design of the continuum structure for multiple loading conditions with stress constraints are presented in the paper,and the s olving method is als o given.The problem is s olved by m odified fully stress method combined with a bundle method to deal with plentiful stress constraints,both the method are sim ple and practical.The multilevel opti2 mization technique is proposed in this paper to make clearer optimal topology of structures.The main idea of the multilevel optimization method is to partition the load cases into several levels according to their magnitude.In every level,we mainly consider the in fluence of the loads belonged to this level.In this way,we av oid the blending of various loads and the dim topological structure.T o s olve the prob2 lem that the single structure stiffness matrix caused by the too big dispersion of element thickness between different levels,it proposed the adjustive method of relevant parameters.Numeral com putations show that the method is effective and efficient. K ey w ords Structure topology optimization;Stress constraints;Continuum structure;Fully stress method;Multilevel optimization technique Correspondent:WANG Jian,E2mail:wangjian0721@https://www.doczj.com/doc/343308439.html,,Fax:+86253322313164 The project supported by the Natural Science F oundation of Shandong Province,China(N o.Y96F03085). Manuscript received20010920,in revised form20011225. 1 引言 在多工况、多约束情况下,结构的最优拓扑往往是超静定的,必须考虑变形协调条件,其数学模型是一个非线性规划问题。文献[1～4]是离散结构拓扑优化方面成功采用非线性规划方法求解的范例。连续体结构拓扑优化方面也有考虑多工况情况的文章发表[5,6],但这方面的工作不多,且没有研究应力约束问题。实际工程结构多半在多种工况下工作,应力约束是最基本的约束条件,所以研究多工况应力约束下连续体结构的拓扑优化问题是非常必要的。 多工况下受到应力约束的结构拓扑优化问题的数学模型可以描述为式(1),用数学规划法求解时自然将其作为一个多约束问题来处理;连续体结构拓扑优化的设计变量很多,采用文献[1～4]中的数学规划方法求解意味着将有浩大的计算工作量,因此一般采用准则法———满应力法解决。用满应力法求解多工况问题时往往使用包络法处理大量的应力约束[7]。包络法的基本思想是把每一个应力约束先单独地考虑,求出在这个应力约束下改进后的新设计变量,然后对每一个设计变量,在所有的值中挑出最大的作为新的设计。这种方法可以保证应力约束条件满足,并且也易于将 机械强度 Journal of Mechanical Strength2003,25(1):055～057 Ξ ΞΞ王 健,男,1962年7月生,山东省济南市长清县人,汉族。山东理工大学交通与车辆工程学院院长,教授,博士,长期从事结构优化研究,发表相关论文20余篇。 20010920收到初稿,20011225收到修改稿。山东省自然科学基金资助项目(Y96F03085)。
结构拓扑优化的发展现状及未来 王超 中国北方车辆研究所一、历史及发展概况 结构拓扑优化是近20年来从结构优化研究中派生出来的新分支，它在计算结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化的工程应用研究还很不成熟，在国外处在发展的初期,尤其在国内尚属于起步阶段。1904 年Michell在桁架理论中首次提出了拓扑优化的概念。自1964 年Dorn等人提出基结构法,将数值方法引入拓扑优化领域，拓扑优化研究开始活跃。20 世纪80 年代初，程耿东和N. Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题，他们开创性的工作引起了众多学者的研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计，开创了连续体结构拓扑优化设计研究的新局面。1993年和提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002 年罗鹰等提出三角网格进化法，该方法在优化过程中实现了退化和进化的统一，提高了优化效率。 二、拓扑优化的工程背景及基本原理 通常把结构优化按设计变量的类型划分成三个层次：结构尺寸优化、形状优化和拓扑优化。尺寸优化和形状优化已得到充分的发展，但它们存在着不能变更结构拓扑的缺陷。在这样的背景下，人们开始研究拓扑优化。拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料的分布问题。寻求一个最佳的拓扑结构形式有两种基本的原理：一种是退化原理，另一种是进化原理。退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上，然后构造适当的优化模型，通过一定的优化方法逐步删减那些不必要的结构元素,直至最终得到一个最优化的拓扑结构形式。进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化，它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。 三、结构拓扑优化设计方法 目前常使用的拓扑优化设计方法可以分为两大类：退化法和进化法。 退化法即传统的拓扑优化方法，一般通过求目标函数导数的零点或一系列迭代计算过程求最优的拓扑结构。目前常用于拓扑优化的退化法有基结构方法、均匀化方法、变密度法、变厚度法等。 基结构方法（GSA）的思路是假定对于给定的桁架节点，在每两个节点之间用杆件连结起来得到的结构称为基结构。按照某种规则或约束,将一些不必要的杆件从基本结构中删除，认为最终剩下的构件决定了结构的最佳拓扑。基结构方法更适合于桁架和框架结构的拓扑优化。基结构法是在有限的子空间内寻优，容易丢失最优解，另外还存在组合爆炸、解的奇异性等问题。 均匀化方法（HA）引入微结构的单胞，通过优化计算确定其材料密度分布，并由此得出最优的拓扑结构。均匀化方法主要应用于连续体的拓扑优化设计，它不仅能用于应力约束和位移约束，也能用于频率约束。目前用均匀化方法来进行拓扑优化设计的有一般弹性问题、热传导问题、周期渐进可展曲面问题、非线性热弹性问题、振动问题和骨改造问题等。 变密度法是一种比较流行的力学建模方式，与采用尺寸变量相比，它更能反映拓
结构拓扑优化的组合准则及应用 丁繁繁* 郭兴文 (河海大学工程力学系,江苏,南京,210098) 摘要：本文研究了拓扑相关荷载作用下连续体结构拓扑优化设计问题，探讨了ESO 方法中单独应用最大拉应变准则或主应力准则来删除单元的问题,提出了基于主压应力删除准则与最大拉应变删除准则的组合优化删除准则,给出了组合准则的迭代步骤.依据所提准则与迭代步骤, 应用Ansys 分析软件对一受拓扑相关径向均布荷载作用的连续体进行了拓扑优化设计,获得了相应的最优拓扑结构,算例表明,本文提出的组合优化法可以消除单一应力删除准则在优化过程中出现的迭代波动问题,能加快拓扑优化的收敛速度. 关键词:拓扑优化, 拓扑相关荷载, 主应力准则, 最大拉应变准则，组合准则 1.前言 结构拓扑优化设计是目前结构优化设计领域最赋有挑战性的研究课题,近十几年来，随着科学技术的进步, 结构拓扑优化设计得到了迅速的发展. 有关结构拓扑优化设计的最新发展,文献以综述的形式作了详细的叙述.连续体结构拓扑优化方法主要有均匀化法、两相法、内力法、变厚度法、变密度法、人工材料、渐进结构优化法及线性规划法等。其中渐进结构优化法(简称ESO)是通过一定的删除准则，将无效或低效的材料逐步去掉,结构将逐渐趋于优化。该方法可采用已有的有限元分析软件,通过迭代过程在计算机上实现,该法的通用性很好。 ESO 法最早是由澳大利亚华裔学者谢忆民于1993年提出来的。随后得到了荣见华等人的发展，成功应用于包含应力、位移（刚度）、临界应力和动力学约束的众多结构拓扑优化领域。基于主应力的ESO 法考虑了实际材料在拉、压应力方面的特性差异，特别适用于一些拉压性质明显的建筑类型，例如桥梁工程，从而改进了ESO 法的工程适用性。 ]4~1[]5[目前,连续体结构拓扑优化研究主要集中在荷载作用位置及作用方向不变情况下的结构拓扑优化问题，而对于荷载作用位置变动情况下的连续体结构拓扑优化研究刚刚起步. ]6[本文研究了荷载位置随拓扑变化而变化作用下的连续体结构拓扑优化问题，该连续体结构是一混凝土受压结构。优化过程中在进行尝试使用不同删除准则的基础上,提出了基于主压应力删除准则与最大拉应变删除准则的组合优化删除准则.依据提出的组合优化删除准则, 应用Ansys 分析软件对一受径向均布荷载作用简支的矩形初始构型进行了拓扑优化设计, 获得了相应的最优拓扑结构,算例表明,本文提出的组合优化法可以消除单一应力删除准则https://www.doczj.com/doc/343308439.html,
3188-1-1.html Sigmund教授所编写的top优化经典99行程序，可以说是我们拓扑优化研究的基础； 每一个新手入门都会要读懂这个程序，才能去扩展，去创新； 99行程序也有好多个版本，用于求解各种问题，如刚度设计、柔顺机构、热耦合问题，但基本思路大同小异； 本文拟对其中的一个版本进行解读，愿能对新手有点小小的帮助。 不详之处，还请论坛内高手多指点 读懂了该程序，只能说是略懂拓扑优化理论了， 我手里就有一些水平集源程序是成千上万行，虽然在99行的基础上成熟了很多，但依然还有很多的发展空间。 源程序如下： %%%% A 99 LINE TOPOLOGY OPTIMIZATION CODE BY OLE SIGMUND, JANUARY 2000 %%% %%%% CODE MODIFIED FOR INCREASED SPEED, September 2002, BY OLE SIGMUND %%% function top(nelx,nely,volfrac,penal,rmin); nelx=80; nely=20; volfrac=0.4; penal=3; rmin=2; % INITIALIZE x(1:nely,1:nelx) = volfrac; loop = 0; change = 1.; % START ITERATION while change > 0.01 loop = loop + 1; xold = x; % FE-ANAL YSIS [U]=FE(nelx,nely,x,penal); % OBJECTIVE FUNCTION AND SENSITIVITY ANAL YSIS [KE] = lk; c = 0.; for ely = 1:nely for elx = 1:nelx n1 = (nely+1)*(elx-1)+ely; n2 = (nely+1)* elx +ely; Ue = U([2*n1-1;2*n1; 2*n2-1;2*n2; 2*n2+1;2*n2+2; 2*n1+1;2*n1+2],1); c = c + x(ely,elx)^penal*Ue'*KE*Ue; dc(ely,elx) = -penal*x(ely,elx)^(penal-1)*Ue'*KE*Ue; end end
升降辊床连杆摇臂结构拓扑优化设计 升降辊床作为一种新型输送设备，具有高速、稳定、易于维护等优点，在各汽车焊装车间得到了广泛应用。文章对辊床连杆摇臂结构进行动力学分析，在此基础上针对摇臂结构进行结构拓扑优化，改善机构应力应变并提升疲劳寿命。 标签：升降辊床;摇臂结构;有限元;拓扑优化;疲劳寿命 引言 “冲压、焊装、油漆和总装”被称为当代汽车制造的四大工艺[1]，在上汽大众仪征工厂焊装车间，焊接工艺种类多达8至10种，用来转运车身的工艺生产线多达12条，拥有德国KUKA自动化机械臂800多台，工艺过程极其复杂，工位数量繁多。基于曲柄连杆摇臂结构的Siemens高速输送升降辊床的大量应用，极大地提高了生产节拍，使生产线实现了柔性生产，产能得到大幅度提高[2]。 1 辊床结构及动力学分析 本文以西门子公司11-0908-1200系列升降辊床为研究对象，主要参数如表1所示。升降辊床主要由底座、升降机構、水平输送辊床和控制系统四大部分组成，实现其升降功能的是一个典型的多连杆机构，并可拆分为两个四杆机构，即前半部分为曲柄连杆摇臂机构[3-4]，后半部分为平行四杆机构，因此，在运动学分析计算中可以忽略后半部分的平行四杆机构，仅分析前半部分的曲柄连杆摇臂机构[5]（图1）。 为了解曲柄连杆摇臂机构在其运动周期内各构件的受力情况，在Adams软件中创建升降辊床曲柄连杆摇臂动力学仿真模型，施加辊床框架及雪橇、车身的重力负载为13000N，直接作用在前后摇臂上，受力方向始终竖直向下，经求解，后摇臂受到来自连杆的峰值拉力为17588N，在升降辊床从低位向高位运行过程中，摇臂克服负载力并将其向上举升，拉力从峰值开始逐渐降低为0N。 2 辊床有限元仿真分析 对辊床连杆结构进行有限元分析。摇臂的制造原材料为Q235B，建立摇臂模型并导入到ANSYS软件中，网格划分后共得到47478个节点、19295个单元。连杆与后摇臂相连的铰接转动副-单孔摇臂关节轴承处，其转动副处最大受力为17588N，选取此瞬态时刻，对后摇臂进行静力学分析，施加负载、约束后进行计算，得到其应力、应变分析结果情况如图2所示。 通过分析发现，在主轴中部轴颈与曲柄连接处是应力集中最严重的部位，从有限元分析结果可以看出，最大应力为91.36MPa，虽然小于摇臂材料的屈服强度235MPa，但这些应力集中部位极易出现疲劳裂纹，直至机械失效损坏，该分析结果与摇臂在实际生产作业中发生的断裂故障一致。
结构拓扑优化设计现状及前景 目前, 最优化设计理论和方法在机械结构设计中得到了深入的研究和广泛的应用。所谓优化设计就是根据具体的实际问题建立其优化设计的数学模型, 并采用一定的最优化方法寻找既满足约束条件又使目标函数最优的设计方案。根据优化问题的初始设计条件, 目前结构优化技术有四大领域: 1) 尺寸优化; 2) 形状优化; 3) 拓扑与布局优化; 4) 结构类型优化。结构尺寸优化是在结构的拓扑确定的前提下, 首先用少量尺寸对结构的某些变动进行表达, 如桁架各单元的横截面尺寸、某些节点位置的变动等, 然后在此基础上建立基于这些尺寸参数的数学模型并采用优化方法对该模型进行求解得到最优的尺寸参数。在尺寸优化设计中, 不改变结构的拓扑形态和边界形状, 只是对特定的尺寸进行调整, 相当于在设计初始条件中就增加了拓扑形态的约束。而结构最初始的拓扑形态和边界形状必须由设计者根据经验或实验确定, 而不能保证这些最初的设计是最优的, 所以最后得到的并不是全局最优的结果。结构形状优化是指在给定的结构拓扑前提下, 通过调整结构内外边界形状来改善结构的性能。以轴对称零件的圆角过渡形状设计的例子。形状设计对边界形状的改变没有约束,和尺寸优化相比其初始的条件得到了一定的放宽,应用的范围也得到了进一步的扩展。拓扑优化设计是在给定材料品质和设计域内,通过优化设计方法可得到满足约束条件又使目标函数最优的结构布局形式及构件尺寸。拓扑设计的初始约束条件更少, 设计者只需要提出设计域而不需要知道具体的结构拓扑形态。拓扑设计方法是一种创新性
的设计方法, 能为我们提供一些新颖的结构拓扑。目前, 拓扑设计理论在柔性受力结构、MEMS 器件及其它柔性微操作机构的设计中得到了广泛的研究。 结构拓扑优化的发展概况 结构拓扑优化包括离散结构的拓扑优化和连续变量结构的拓扑优化。近10 年来, 结构拓扑优化设计虽然取得了一些进展, 但大部分是针对连续变量的, 关于离散变量的研究为数甚少。由于离散变量优化的目标函数和约束函数是不连续、不可微的, 可行域退化为不连通的可行集, 所以难度远大于连续变量优化问题。在离散结构中, 桁架在工程中的应用较为广泛, 由于其重要性, 也由于其分析比较简单, 桁架结构的拓扑优化在文献中研究得最多. 结构拓扑优化的历史可以追溯到1904 年Michell提出的桁架理论, 但这一理论只能用于单工况并依赖于选择适当的应变场, 不能应用于工程实际。1964 年Dorn、Gomory、Greenberg 等人提出基结构法( ground structure approach) , 将数值方法引入该领域, 此后拓扑优化的研究重新活跃起来, 陆续有一些解析和数值方面的理论被 提出来。所谓基结构就是一个由结构节点、荷载作用点和支承点组成的节点集合, 集合中所有节点之间用杆件相连的结构。该方法的基本思路是: 从基结构的模型出发, 应用优化算法( 数学规划法或准则法) , 按照某种规划或约束, 将一些不必要的杆件从基结构中删除, 例如截面积达到零或下限的杆件将被删掉, 并认为最终剩下的杆件 决定了结构的最优拓扑。因此应用基结构, 可以将桁架拓扑优化当作
文章编号:1004Ο8820(2003)02Ο0138206 具有多种约束的连续体结构拓扑优化 江允正,王子辉,初明进 (烟台大学土木工程系,山东烟台264005) 摘要:对于具有多种约束条件的连续体结构的拓扑优化设计,本文提出一种通用优化方 法:首先用优化方法确定微孔或称为基点的位置,然后再扩大微孔并确定其边界.文中对 于具有应力和位移约束的几个平面问题进行拓扑优化,计算结果十分令人满意. 关键词:结构拓扑优化;结构优化;连续体; 中图分类号:TP391.72 文献标识码:A 近年来,Bendsoe 和K ikuchi [1]等广泛采用连续体拓扑优化的均匀方法.首先从连续介质中人为地引进某一形式的微结构,例如周期性分布的微孔洞;然后用以数学中扰动理论为基础的均匀化方法这一数学工具建立材料的宏观弹性性质和微结构尺寸的关系,连续介质的拓扑优化就转化为决定微结构尺寸最优分布的尺寸优化问题,可以采用成熟的尺寸优化算法.迄今为止的均匀化方法还不能给出带有微观结构的材料的宏观许用应力和微结构尺寸的关系,因此到目前为止均匀优化方法可以求解的拓扑优化问题还很有限.均匀化方法的另一缺点是求得的最终设计可能具有很不清晰的拓扑,即结构中有的区域是相对密度介于0和1之间的多孔介质;文献[2]提出修改的满应力法来求解受应力约束的平面弹性体的拓扑优化问题,也仅能考虑应力约束问题;文献[3]提出统一骨架与连续体的结构拓扑优化的ICM 理论与方法.这些方法,基本上都采用有限元法进行结构分析,为了使边界光滑,不得不划分很细的单元,对于一般平面问题,单元数目都在数千个之上,计算效率低.总之,拓扑优化是最具挑战性而又困难的问题,优化方法仍然处在发展初期.这一领域迫切需要取得进展,开发通用的算法仍是挑战. 如上所述,采用均匀方法时,首先从连续介质中人为地引进某一形式的微结构,例如周期性分布的微孔洞.我们认为微孔洞的数量和位置应该用优化方法确定.并称这种微孔的中心叫做删除区的基点.然后扩大微孔,用优化方法确定孔的边界.于是,连续体结构的拓扑优化,可以归结为确定删除区的基点位置及其边界的问题. 1 方 法 对于一个二维连续体,当给定外载和支承位置时,满足应力、位移等各种约束条件下的结构最优拓扑问题,都可以按如下步骤来求解: 收稿日期:2002-12-17 作者简介:江允正(1942-),男,湖南衡阳人,教授,主要从事结构优化方向教学与研究工作. 第16卷第2期 烟台大学学报(自然科学与工程版)Vol.16No.22003年4月Journal of Y antai University (Natural Science and Engineering Edition ) Apr.2003
( 安全管理 ) 单位：_________________________ 姓名：_________________________ 日期：_________________________ 精品文档 / Word文档 / 文字可改 连续体结构拓扑优化方法及存在问题分析(最新版) Safety management is an important part of production management. Safety and production are in the implementation process
连续体结构拓扑优化方法及存在问题分析 (最新版) 文章深入分析国内外连续体结构拓扑优化的研究现状，介绍了拓扑优化方法的发展及实现过程中存在的问题。对比分析了均匀化方法，渐进结构优化法，变密度法的优缺点。研究了连续体结构拓扑优化过程中产生数值不稳定现象的原因，重点讨论了灰度单元，棋盘格式，网格依赖性的数值不稳定现象，并针对每一种数值不稳定现象提出了相应的解决办法。 结构拓扑优化设计的主要对象是连续体结构，1981年程耿东和Olhof在研究中指出：为了得到实心弹性薄板材料分布的全局最优解，必须扩大设计空间，得到由无限细肋增强的板设计。此研究被认为是近现代连续体结构拓扑优化的先驱。 目前，国内外学者对结构拓扑优化问题进行了大量研究，这些
研究大多数建立在有限元法结构分析的基础上，但由于有限元法中单元网格的存在，结构拓扑优化过程中常常出现如灰度单元，网格依赖性和棋盘格等数值不稳定的现象。本文介绍了几种连续体结构拓扑优化方法及每种方法存在的问题，并提出了相应的解决办法。 1.拓扑优化方法 连续体结构拓扑优化开始于1988年Bendoe和Kikuchi提出的均匀化方法，此后许多学者相继提出了渐进结构优化方法、变密度法等拓扑优化数学建模方法。 1.1.均匀化方法 均匀化方法即在设计区域内构造周期性分布的微结构，这些微结构是由同一种各向同性材料实体和孔洞复合而成。采用有限元方法进行分析，在每个单元内构造不同尺寸的微结构，微结构的尺寸和方向为拓扑优化设计变量。1988年Bendsoe研究发现，通过在结构中引入具有空洞微结构的材料模型，将困难的拓扑设计问题转换为相对简单的材料微结构尺寸优化问题。 很多学者发展了均匀化方法，Suzhk进行了基于均匀化方法结构
第19卷 第3期应用力学学报Vol.19 No.3 2002年9月CHINESE JOURNAL OF APPL IE D MECHANICS Sep.2002 文章编号:100024939(2002)0320050204 结构拓扑优化设计的三角网格进化法Ξ 罗 鹰 段宝岩 (西安电子科技大学 西安 710072) 摘要:针对进化式拓扑优化方法的不足,提出了一种基于遗传算法的新型进化式拓扑优化方法—三角网格进化法,该方法不仅能够同时进行拓扑、形状与截面变量优化设计,而且在优化过程中实现了退化和进化的统一,提高了优化效率。另外本文还首次对结构类型变量进行了优化计算,取得了有益的结果。最后几个数值算例证明了本方法的可行性和有效性。 关键词:拓扑优化;进化法;类型优化;遗传算法 中图分类号:039TB121 文献标识码: A 1 引 言 工程结构拓扑优化方法可分为两类:退化法和进化法。退化法又可进一步分为基结构方法(ground structural approach)[1]和均匀化方法(ho2 mogenization method)[2],退化法的基本思想是在优化前将结构所有可能杆单元(对基结构方法而言)或所有材料(对均匀化方法而言)都加上,然后构造适当的优化模型,通过一定的优化方法逐步删减那些不必要的结构元素(杆单元及节点)或材料,直至最终得到一个最优化的拓扑结构形式。当然,在删减的同时也可能伴随着少量结构元素的再加入。进化法[3～6]正好与退化法相反,它是从另一个途径考虑问题。根据给定的固定节点与载荷,首先给出简单拓扑结构形式,然后通过一定的优化策略不断增加结构元素,直到获得最优的拓扑结构。K irsch[5,6]曾对此类方法进行过分析与展望,并且由William在1995年提出了自然生长方法[3],Mc Keown在1998年又提出了节点增加方法[4]。它们的不足之处在于,优化过程中,只有结构元素(包括杆单元和节点)的增加而不能够删减。另外,根据目前所掌握的文献看,结构类型变量优化还未被问津。本文利用遗传算法(G A)将结构类型也作为一类设计变量,对它进行了数学优化计算的尝试。 2 优化模型 本文讨论的是结构的整体优化问题,设计目标是使结构整体重量最轻(或体积最小),而约束条件包括应力约束以及各节点坐标位移约束。设计变量包括结构类型、拓扑、可动节点坐标以及单元截面积四种参数。由于遗传算法(G A)[5,7,8,9]不能直接处理结构优化中各设计变量,而必须将它们转换成遗传空间中由基因个体排列组成的染色体或个体。为此,引入以下几组参数: 211 结构类型参数αi 杆系结构的类型不仅有桁架、刚架(梁)结构,还有杆、梁组合结构(即结构中既有杆单元又有梁单元)。为此引入参数αi(i=1,2,…,N)分别代表结构中各单元的类型。其中,N表示结构单元数。其数学表达式为: α i = 0 单元i为杆单元 1 单元i为梁单元 (i=1,2,…,N) (1)结构的总刚度方程为: Ξ基金项目:国家自然科学基金项目(95635150) 来稿日期:2001202220 修回日期:2002202227第一作者简介:罗鹰,男,1970年生,西安电子科技大学机电工程学院博士生;研究方向:面向工程的广义优化1
编号：SY-AQ-00556 ( 安全管理) 单位：_____________________ 审批：_____________________ 日期：_____________________ WORD文档/ A4打印/ 可编辑 连续体结构拓扑优化方法及存 在问题分析 Topology optimization method of continuum structure and analysis of existing problems
连续体结构拓扑优化方法及存在问 题分析 导语：进行安全管理的目的是预防、消灭事故，防止或消除事故伤害，保护劳动者的安全与健康。在安全管理的四项主要内容中，虽然都是为了达到安全管理的目的，但是对生产因素状态的控制，与安全管理目的关系更直接，显得更为突出。 文章深入分析国内外连续体结构拓扑优化的研究现状，介绍了拓扑优化方法的发展及实现过程中存在的问题。对比分析了均匀化方法，渐进结构优化法，变密度法的优缺点。研究了连续体结构拓扑优化过程中产生数值不稳定现象的原因，重点讨论了灰度单元，棋盘格式，网格依赖性的数值不稳定现象，并针对每一种数值不稳定现象提出了相应的解决办法。 结构拓扑优化设计的主要对象是连续体结构，1981年程耿东和Olhof在研究中指出：为了得到实心弹性薄板材料分布的全局最优解，必须扩大设计空间，得到由无限细肋增强的板设计。此研究被认为是近现代连续体结构拓扑优化的先驱。 目前，国内外学者对结构拓扑优化问题进行了大量研究，这些
研究大多数建立在有限元法结构分析的基础上，但由于有限元法中单元网格的存在，结构拓扑优化过程中常常出现如灰度单元，网格依赖性和棋盘格等数值不稳定的现象。本文介绍了几种连续体结构拓扑优化方法及每种方法存在的问题，并提出了相应的解决办法。 1.拓扑优化方法 连续体结构拓扑优化开始于1988年Bendoe和Kikuchi提出的均匀化方法，此后许多学者相继提出了渐进结构优化方法、变密度法等拓扑优化数学建模方法。 1.1.均匀化方法 均匀化方法即在设计区域内构造周期性分布的微结构，这些微结构是由同一种各向同性材料实体和孔洞复合而成。采用有限元方法进行分析，在每个单元内构造不同尺寸的微结构，微结构的尺寸和方向为拓扑优化设计变量。1988年Bendsoe研究发现，通过在结构中引入具有空洞微结构的材料模型，将困难的拓扑设计问题转换为相对简单的材料微结构尺寸优化问题。 很多学者发展了均匀化方法，Suzhk进行了基于均匀化方法结
某商用车白车身结构疲劳寿命分析与优化设计 李明1李源2陈斌3 （1湖南工业大学机械工程学院，湖南株洲，412008；2国防科学技术大学指挥军官基础教育学院，湖南长沙，410072；3 湖南大学汽车车身先进设计制造国家重点实验室，湖南长沙，410082） 摘要：本文基于应力分析结果，采用有效的疲劳寿命预估方法，利用专业耐久性疲劳寿命分析系统MSC．Fatigue对该型商用车白车身进行S-N全寿命分析，得其疲劳寿命分布与危险点的寿命值。采用 结构优化、合理选材等方法，提高白车身结构的疲劳寿命。 关键词：白车身；有限元；静态分析；疲劳寿命分析；优化 Body-in-white Fatigue Analysis and Optimization Design of the Commercial Vehicle LI Ming1, LI Yuan2, CHEN Bin3 （1 School of Mechanical Engineering , Hunan University of Technology, Zhuzhou, Hunan 412008, China; 2 College of Basic Education for Officers, National University of Defense Technology, Changsha, Hunan 410072, China；3 State Key Laboratory of Advanced Design and Manufacture for Vehicle Body, Hunan University, Changsha, Hunan 410082,China） Abstract:Based on the results of stress analysis, this paper took the effective way of the fatigue life estimating, used the professional durability fatigue life analysis system MSC. Fatigue, and the S-N life-cycle analysis of the certain type of commercial vehicle body-in-white finite element model, got the distribution of fatigue life and the fatigue life value of the danger points. Finally, by the structural optimization and material selection, writer improved the fatigue life of white body structure. Keywords: B ody-in-white structure, FEM, Static analysis, Fatigue lifetime analysis; Optimization 0 前言 在车身结构疲劳领域的国内研究中，1994年，江苏理工大学陈龙在建立了车辆驾驶室疲劳强度计算的力学和数学模型基础上，提出了车辆驾驶室疲劳强度研究方法[1]。2001 年，清华大学孙凌玉[2]等首次计算机模拟了汽车随机振动过程。2002年，上海汇众汽车制造有限公司王成龙[3]等应用FATIGUE 软件的分析，结合疲劳台架试验，探讨了疲劳强度理论在汽车产品零部件疲劳寿命计算中的应用，提出了提高零部件疲劳强度的方法。2004年，同济大学汽车学院靳晓雄[4]等人提到进行零部件疲劳寿命预估，
结构优化设计的综述与发展 摘要：结构优化设计，就是在计算机技术等高科技手段的支持下，为了提升机械产品的性能、工作效率，延长机械产品的工作寿命，对机械产品的尺寸、形状、拓扑结构和动态性能进行优化的过程。这是机械行业发展的必然要求，也是信息时代的必然要求。结构优化设计，必须在保证机械产品满足工作需要的前提下，通过科学的计算来实行。文章将简单对结构优化设计的发展状况进行介绍，列举几种优化设计方法，以及讨论未来优化的发展情况。 关键词：结构优化设计发展优化设计方法 1 结构优化设计 结构优化简单来说就是在满足一定的约束条件下，通过改变结构的设计参数，以达到节约原材料或提高结构性能的目的。结构优化设计通常是指在给定结构外形，给定结构各元件的材料和相关载荷及整个结构的强度、刚度、工艺等要求的条件下，对结构进行整体和元件优化设计。结构优化设计一般由设计变量、约束条件和目标函数三要素组成。评价设计优、劣的标准，在优化设计中称为目标函数；结构设计中以变量形式参与的称为设计变量；设计时应遵守的几何、刚度、强度、稳定性等条件称为约束条件，而设计变量、约束函数与目标函数一起构成了优化设计的数学模型。结构优化的目的是让设计的结构利用材料更经济、受力分布更合理。 结构优化设计根据设计变量选取的不同可以分为截面（尺寸）优化、形状优化、拓扑优化三个层次。尺寸优化是选取结构元件的几何尺寸作为设计变量，例如，杆元截面积、板元的厚度等等[1]。而形状优化是选取结构的内部形状或者是节点位置作为设计变量。拓扑优化就是选取结构元件的有无作为设计变量，为0-1型逻辑型设计变量。 2 结构优化设计研究概况与现状 结构优化设计最早可以追溯到17世纪，伽利略和伯努利对弯曲梁的研究从而引发了变截面粱形状优化的问题。后来Maxwell和Michell提出了单载荷仅有应力约束条件下最小重量桁架结构布局的基本理论，为系统地分析结构优化理论作出了重大的贡献。然而长期以来，由于缺乏高速可靠的计算手段和理论，结构优化设计一直无法获取较大发展。 到上世纪六十年代，有限元技术借助于计算机技术，得到了极大的发展。1960年Schmit在求解多种载荷情况下弹性结构的最小重量问题时，首次在结构优化中引入入数学规划理论，并与有限元方法结合应用，形成了全新的结构优化思想，标志着现代结构优化技术的开始[2]。 1973年Zienkiewicz和Campbell[3]在解决水坝的形状优化问题时，首次以节点坐标作为设计变量，在结构分析方面使用了等参元，在优化方法上使用了序列线性规划的方法。其后，众多的学者在此基础上，逐渐发展形成了使用边界形状参数化方法描述连续
第2卷第1期2011年2月航空工程进展 A DV A N CES IN A ERON A U T ICA L SCIEN CE A N D EN GIN EERIN G Vo l 12N o 11Feb 1 2011 收稿日期:2010-12-01; 修回日期:2011-01-20基金项目:教育部长江学者创新团队项目(Irt0906)通信作者:姚卫星,w xyao@https://www.doczj.com/doc/343308439.html, 文章编号:1674-8190(2011)01-001-12 连续体结构拓扑优化方法评述 夏天翔,姚卫星 (南京航空航天大学飞行器先进设计技术国防重点学科实验室,南京 210016) 摘 要:连续体结构拓扑优化在优化中能产生新的构型,对实现自动化智能结构设计具有重要意义。目前,连续体结构拓扑优化方法主要有:均匀化方法、变厚度法、变密度法、渐进结构优化方法、水平集法、独立连续映射方法。本文首先系统回顾了以上方法的发展历程,介绍了它们的研究现状。其次,通过对比以上拓扑优化方法对若干典型算例的优化结果,表明以上方法都有较好的减重效果。最后,对以上方法进行了总结,列出了它们的优缺点和发展方向。 关键词:拓扑优化;均匀化方法;变厚度法;变密度法;渐进结构优化方法;水平集法;独立连续映射方法中图分类号:V 211.7 文献标识码:A A Survey of Topology Optimization of Continuum Stru cture Xia Tianx iang ,Yao Weix ing (K ey L abor ator y of F undamental Science fo r N atio nal Defense -adv anced Design T echno lo gy of F lig ht V ehicle,Nanjing U niver sity o f A eronautics and A st ronautics,N anjing 210016,China) Abstract:A s the to po log y optim izat ion o f continuum structure can pr oduce new config ur atio ns during the optim-i zatio n,it is significant for automatic str ucture design.A t present,the most commo nly used t opolo gy o ptimiza -t ion methods of continuum st ructur e ar e:the ho mog enization method,var iable t hickness method,v ariable dens-i t y metho d,evo lutio nar y str uctur al o pt imizatio n met ho d,lev el set metho d,independent co ntinuous mapping method.Firstly,the develo pment pro cesses of above metho ds ar e sy stematically review ed,their cur rent r e -sear ch is br iefly intro duced in this paper.T hen,these methods ar e com par ed and discussed t hr ough a number of typical ex amples.T he typical ex amples show that all of above methods have gr eat abilities to r educe w eig ht.F-i nally ,the adv ant ag es,disadv ant ag es and dev elo pment directio ns of abov e metho ds ar e discussed. Key words:to po lo gy o ptimization;homog enizat ion metho d;va riable thickness method;var iable density method;evolutionar y structure optimization metho d;lev el set method;independent continuo us mapping method 0 引言 按照设计变量的不同,结构优化可分为以下三个层次:尺寸优化、形状优化和拓扑优化。结构拓 扑优化能在给定的外载荷和边界条件下,通过改变结构拓扑使结构在满足约束的前提下性能达到最优。与尺寸优化、形状优化相比,结构拓扑优化的经济效果更为明显,在优化中能产生新的构型,是 结构实现自动化智能设计所必不可少的。 按照优化对象的性质,拓扑优化可分为离散体拓扑优化和连续体拓扑优化两种。连续体拓扑优化与离散体拓扑优化相比,在应用范围更广的同 时,模型描述困难,设计变量多,计算量大。在过去很长一段时间里,连续体拓扑优化发展得十分缓慢,直到1988年Bendso e 等人[1] 提出均匀化方法之后,它才得到了迅速发展。目前,国内外学者对结构拓扑优化问题已经进行了大量研究[2-9]。目前最常用的连续体拓扑优化方法有均匀化方法、变厚 度法、变密度法、渐进结构优化方法(ESO)、水平集法(Level set)、独立连续映射方法(ICM)等。从拓