文档库 最新最全的文档下载
当前位置:文档库 › 拓扑优化

拓扑优化

拓扑优化
拓扑优化

拓扑优化研究方法综述

结构拓扑优化是近20年来从结构优化研究中派生出来的新分支,它在计算

结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化

的工程应用研究还很不成熟,在国外处在发展的初期,尤其在国内尚属于起步

阶段。1904年Michell在桁架理论中首次提出了拓扑优化的概念。自1964年Dorn等人提出基结构法,将数值方法引入拓扑优化领域,拓扑优化研究开始活跃。20世纪80年代初,程耿东和N.Olhoff在弹性板的最优厚度分布研究中首

次将最优拓扑问题转化为尺寸优化问题,他们开创性的工作引起了众多学者的

研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化

设计,开创了连续体结构拓扑优化设计研究的新局面。1993年XieYM和StevenGP提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002年罗鹰等提出三角网格进化法,该方法在优化过程中

实现了退化和进化的统一,提高了优化效率。

通常把结构优化按设计变量的类型划分成三个层次:结构尺寸优化、形状优化和拓扑优化。尺寸优化和形状优化已得到充分的发展,但它们存

在着不能变更结构拓扑的缺陷。在这样的背景下,人们开始研究拓扑优化。拓

扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻

求最优材料的分布问题。寻求一个最佳的拓扑结构形式有两种基本的原理:一

种是退化原理,另一种是进化原理。退化原理的基本思想是在优化前将结构所

有可能杆单元或所有材料都加上,然后构造适当的优化模型,通过一定的优化

方法逐步删减那些不必要的结构元素,直至最终得到一个最优化的拓扑结构形式。进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化,

它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。

退化法即传统的拓扑优化方法,一般通过求目标函数导数的零点或一

系列迭代计算过程求最优的拓扑结构。目前常用于拓扑优化的退化法有基结构

方法、均匀化方法、变密度法、变厚度法等。

进化法是一类全局寻优方法,目前常用于拓扑优化的进化法主要有遗

传算法、模拟退火算法和渐进结构优化法等。

什么是拓扑优化?

拓扑优化是指形状优化,有时也称为外型优化。拓扑优化的目标是寻找承

受单载荷或多载荷的物体的最佳材料分配方案。这种方案在拓扑优化中表现为“最大刚度”设计。

与传统的优化设计不同的是,拓扑优化不需要给出参数和优化变量的定义。目标函数、状态变量和设计变量(参见“优化设计”一章)都是预定义好的。用

户只需要给出结构的参数(材料特性、模型、载荷等)和要省去的材料百分比。

拓扑优化的目标——目标函数——是在满足结构的约束(V)情况下减少

结构的变形能。减小结构的变形能相当于提高结构的刚度。

下面是从振动论坛的回帖,有帮助的:

========================================================= ===============

求助:结构动力学优化设计(拓扑优化)

veasha 发表于: 2009-2-27 11:44 来源: 振动资讯

我的研究课题就是结构动力学优化设计,动力特性和动力响应优化都有涉及。但是组里没有人做过,没有任何相关的背景资料,我想知道在动力学优化设计领域常用的软件、书籍等信息,另外,拓扑优化适合用来做动力学的优化设计吗?ANSYS 适合分析这类问题吗?

希望大家能给予帮助。谢谢~!

最新回复

lq12131010 at 2009-2-27 13:17:59

这些问题我略略研究过。现在也在做一些东西。你是哪个学校啊?

建议没学过优化的,先看看程先生的书。

学过优化,就看bendsoe的书 Olhoff的论文。

软件 nastran hyperwork.。

veasha at 2009-2-27 16:24:41

我是北航的。没有学过优化,程先生是指?因为组里基本都是用ANSYS,缺点就

是优化算法太少了,就两种,而且参数化建模好像比较麻烦。你说的

nastran hyperwork我们组没有人用过。。。我现在的设想是自己写优化算法,但优化过程仍是用ANSYS,这样可行吗?

谢谢!

mjhzhjg at 2009-3-12 13:14:07

拓扑优化应个主要是指拓扑结构发生变化,这方面的软件还是

tosca hyperwork optishape软件要厉害一些。动力学优化设计用ansys就可以了,当然要做二维的拓扑优化,ansys也够了,我前面说的几个软件在三维方面要强一些。

至于书籍,如果你不准备系统学习,只是用必须的话,一般ansys的书都有动力学优化的内容,看看就行了

基于拓扑优化的车身结构研究---经典

基于拓扑优化的车身结构研究 瞿元王洪斌张林波吴沈荣 奇瑞汽车股份有限公司,安徽芜湖,241009 摘要:随着CAE技术的发展,虚拟仿真技术在汽车开发中的作用也愈来愈显著。而前期工程阶段,如何布置出合理的车身骨架架构,一直是个相对空白的地带,也是整车正向开发过程中绕不过的坎。尽管研发工程师根据经验,参照现有车型的结构特点,也能进行车身骨架架构的设定,但总是缺乏有效手段直观地反映不同车型结构布置的特点。本文用拓扑优化的方法,从结构基本特征的角度来审视这一问题,并运用该方法对某SUV车身结构进行研究,获得一些直观性的结论。 关键词:车身,前期工程,拓扑优化 1引言 随着对整车研发过程认识的加深,以及对正向开发过程的探索,在车型开发前期,对车身结构做出更合理的规划显得愈来愈重要。常规的研发思路之一是通过参考已有车型的结构,经过适当的修改,形成新的结构,并用于新车型中。但是对于原始车型的设计思路、结构布置的原因等缺乏系统的理解,或者理解不深,往往在更改过程中产生新的问题。为了部分解决上述问题,本文从结构拓扑优化的角度,对某SUV 车型车身结构的总体布置进行初步探讨,以期加深对结构布置的理解。 2研究方法概述 合理化的车身结构,是满足整车基本性能的重要保障。为了能够实现结构的最优布置,文献[1]使用了拓扑优化工具来布置车身结构。其基本思路是从造型以及车内空间布置出发,建立车身空间的基础网格模型,然后根据一定的工况要求,对基础网格进行拓扑分析,并根据拓扑结果建立梁、板壳模型,并进行多项性能的优化,从而实现车身结构的正向开发。本文借助于该思想,建立研究对象的结构空间包络,并对该包络进行拓扑分析,然后将仿真结果与原始结构进行比较,寻找车身结构中的关键点,推测初始结构可能的布置思想,从而加深对该研究思路的理解。其基本过程如下图所示:

工程结构拓扑优化的理论研究及应用_满宏亮.

提要 本文首先介绍了国内外拓扑优化技术的研究发展现状,讨论了拓扑优化的 原理、方法以及各种拓扑优化算法。其次,着重研究了SIMP 材料插值方法,建立了基于SIMP 理论的连续体结构拓扑优化模型,选取准则优化法对其密 度迭代格式进行了推导;并且利用MATLAB软件编程实现,有效地进行了平面结构的分析和拓扑优化设计。然后,分析了拓扑优化中的数值计算不稳定性现象,研究了能够有效消除拓扑优化中的数值计算不稳定性现象的各种解决方法,并对其进行了比较。最后,利用连续体结构拓扑优化求解理论和算法,使用结构有限元分析软件Hyperworks 对具体工程结构部件进行了拓扑优化设计研究,成功地应用到了实际工程问题中,算例结果表明了该优化方法的有效性和正确性。 关键词:有限元拓扑优化材料插值模型数值计算不稳定性优化 求解算法 Key words: FEA Topology optimization Material Interpolation Model Numerical Calculation Instabilities Optimization Solution Algorithm-i- 目录 第一章绪论 (1) 1.1 前言 (1) 1.2 国内外拓扑优化研究概况 (3) 1.3 本文研究内容及意义 (9) 第二章现代结构拓扑优化理论 (11) 2.1 拓扑的概念 (11) 2.1.1 拓扑学的由来 (11) 2.1.2 拓扑学及拓扑性质 (13) 2.2 结构拓扑优化原理和方法 (16) 2.2.1 拓扑优化的基本原理 (17) 2.2.2 结构拓扑优化设计方法 (17) 2.2.3 拓扑优化设计方法比较 (21) 2.3 拓扑优化设计的优化算法概述 (22) 2.3.1 优化算法分类 (22) 2.3.2 拓扑优化常用算法 (24) 第三章连续体结构拓扑优化的模型建立与求解算法 (27) 3.1 连续体结构拓扑优化设计的模型描述 (29) 3.2 数学模型的有限元离散 (34) 3.2.1 单元应变和应力.........................................34吉林大学硕士研究生学位论文-ii- 3.2.2 单元平衡方程 (35) 3.2.3 连续体结构拓扑优化的数学模型的有限元离散形式 (38) 3.3 基于SIMP 理论的优化准则法 (39) 第四章结构拓扑优化程序实现 (45) 4.1 基于SIMP 理论的优化准则法迭代分析流程 (45) 4.2 优化过程的MA TLAB 编程实现 (47) 4.3 计算实例 (48)

北航拓扑优化程序学习报告

拓扑优化的 99行程序学习报告4月19日2011 《结构优化设计》课程学习报告 任课教师:李书

一、前言: 在最近的结构优化设计课程上学习了O.Sigmund的《A 99 line topology optimization code written in Matlab》一文,对拓扑优化的理论原理与实际的计算机程序实现都有了一定的理解,文章主要是通过拓扑优化的原理来实现对简单结构的静力学问题的优化求解,而编写的代码仅有99行,包括36行的主程序,12行的OC优化准则代码,16行的网格过滤代码和35行的有限元分析代码。 自1988 年丹麦学者Bendsoe与美国学者Kikuchi提出基于均匀化方法的结构拓扑优化设计基本理论以来,均匀化方法应用到具有周期性结构的材料分析中,近几年该方法已经成为分析夹杂、纤维增强复合材料、混凝土材料等效模量,以及材料的细观结构拓扑优化常用的手段之一。其基本思想是在组成拓扑结构的材料中引入微结构,优化过程中以微结构的几何尺寸作为设计变量,以微结构的消长实现其增删,并产生介于由中间尺寸微结构组成的复合材料,从而实现了结构拓扑优化模型与尺寸优化模型的统一。文章就是通过均匀化的基础,结合拓扑结构优化的工程实际,以计算机模拟的方法将拓扑优化的一般过程呈现出来,有助于初涉拓扑优化的读者对拓扑优化有个基础的认识。 二、拓扑优化问题描述 为了简化问题的描述,文中假设设计域是简单的矩形形式,且在进行有限元离散的时候采用正方形单元对其进行离散。这样不仅便于进行单元离散和单元编号,也利于对结构进行几何外形的描述。 一般说来,基于指数逼近法的拓扑优化最小化的问题可作如下描述: 文中采用的对结构材料属性的描述是所谓的“指数逼近法”或者称为SIMP 逼近法,即(Solid Isotropic Material with Penalization带惩罚因子的各项同性材料模型法),该方法是拓扑优化中常用的变密度材料插值模型中最具代表性的一种。

ANSYS拓扑优化原理讲解以及实例操作

拓扑优化是指形状优化,有时也称为外型优化。 拓扑优化的目标是寻找承受单载荷或多载荷的物体的最佳材料分配方案。这种方案在拓扑优化中表现为“最大刚度”设计。与传统的优化设计不同的是,拓扑优化不需要给出参数和优化变量的定义。目标函数、状态变量和设计变量(参见“优化设计”一章)都是预定义好的。用户只需要给出结构的参数(材料特性、模型、载荷等)和要省去的材料百分比。给每个有限元的单元赋予内部伪密度来实现。这些伪密度用PLNSOL ,TOPO 命令来绘出。拓扑优化的目标——目标函数——是在满足结构的约束(V )情况下减少结构的变形能。减小结构的变形能相当于提高结构的刚度。这个技术通过使用设计变量。 结构拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料分布的问题。通过拓扑优化分析,设计人员可以全面了解产品的结构和功能特征,可以有针对性地对总体结构和具体结构进行设计。特别在产品设计初期,仅凭经验和想象进行零部件的设计是不够的。只有在适当的约束条件下,充分利用拓扑优化技术进行分析,并结合丰富的设计经验,才能设计出满足最佳技术条件和工艺条件的产品。连续体结构拓扑优化的最大优点是能在不知道结构拓扑形状的前提下,根据已知边界条件和载荷条件确定出较合理的结构形式,它不涉及具体结构尺寸设计,但可以提出最佳设计方案。拓扑优化技术可以为设计人员提供全新的设计和最优的材料分布方案。拓扑优化基于概念设计的思想,作为结果的设计空间需要被反馈给设计人员并做出适当的修改。最优的设计往往比概念设计的方案结构更轻,而性能更佳。经过设计人员修改过的设计方案可以再经过形状和尺寸优化得到更好的方案。 优化拓扑的数学模型 优化拓扑的数学解释可以转换为寻求最优解的过程,对于他的描述是:给定系统描述和目标函数,选取一组设计变量及其范围,求设计变量的值,使得目标函数最小(或者最大)。一种典型的数学表达式为: ()()()12,,0 ,,0min ,g x x v g x x v f x v ?=??≤?? ?? && 式中,x -系统的状态变量;12g g 、-一等式和不等式的结束方程;(),f x v -目标函数; v -设计变量。 注:在上述方程中,x 作为系统的状态变量,并不是独立的变量,它是由设计变量得 出的,并且与设计变量相关。 优化拓扑所要进行的数学运算目标就是,求取合适的设计变量v ,并使得目标函数值最小。 基于ANSYS 的优化拓扑的一般过程 (进行内容排版修改) 在ANSYS 中,进行优化拓扑,一般分为6个步骤。具体流程见图5-1:

结构拓扑优化的发展现状及未来

结构拓扑优化的发展现状及未来 王超 中国北方车辆研究所一、历史及发展概况 结构拓扑优化是近20年来从结构优化研究中派生出来的新分支,它在计算结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化的工程应用研究还很不成熟,在国外处在发展的初期,尤其在国内尚属于起步阶段。1904 年Michell在桁架理论中首次提出了拓扑优化的概念。自1964 年Dorn等人提出基结构法,将数值方法引入拓扑优化领域,拓扑优化研究开始活跃。20 世纪80 年代初,程耿东和N. Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题,他们开创性的工作引起了众多学者的研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计,开创了连续体结构拓扑优化设计研究的新局面。1993年和提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002 年罗鹰等提出三角网格进化法,该方法在优化过程中实现了退化和进化的统一,提高了优化效率。 二、拓扑优化的工程背景及基本原理 通常把结构优化按设计变量的类型划分成三个层次:结构尺寸优化、形状优化和拓扑优化。尺寸优化和形状优化已得到充分的发展,但它们存在着不能变更结构拓扑的缺陷。在这样的背景下,人们开始研究拓扑优化。拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料的分布问题。寻求一个最佳的拓扑结构形式有两种基本的原理:一种是退化原理,另一种是进化原理。退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上,然后构造适当的优化模型,通过一定的优化方法逐步删减那些不必要的结构元素,直至最终得到一个最优化的拓扑结构形式。进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化,它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。 三、结构拓扑优化设计方法 目前常使用的拓扑优化设计方法可以分为两大类:退化法和进化法。 退化法即传统的拓扑优化方法,一般通过求目标函数导数的零点或一系列迭代计算过程求最优的拓扑结构。目前常用于拓扑优化的退化法有基结构方法、均匀化方法、变密度法、变厚度法等。 基结构方法(GSA)的思路是假定对于给定的桁架节点,在每两个节点之间用杆件连结起来得到的结构称为基结构。按照某种规则或约束,将一些不必要的杆件从基本结构中删除,认为最终剩下的构件决定了结构的最佳拓扑。基结构方法更适合于桁架和框架结构的拓扑优化。基结构法是在有限的子空间内寻优,容易丢失最优解,另外还存在组合爆炸、解的奇异性等问题。 均匀化方法(HA)引入微结构的单胞,通过优化计算确定其材料密度分布,并由此得出最优的拓扑结构。均匀化方法主要应用于连续体的拓扑优化设计,它不仅能用于应力约束和位移约束,也能用于频率约束。目前用均匀化方法来进行拓扑优化设计的有一般弹性问题、热传导问题、周期渐进可展曲面问题、非线性热弹性问题、振动问题和骨改造问题等。 变密度法是一种比较流行的力学建模方式,与采用尺寸变量相比,它更能反映拓

ABAQUS拓扑优化手册

ABAQUS拓扑优化分析手册/用户手册 分析手册: 13. Optimization Techniques优化技术 13.1 结构优化:概述 13.1.1 概述 ABAQUS结构优化是一个帮助用户精细化设计的迭代模块。结构优化设计能够使得结构组件轻量化,并满足刚度和耐久性要求。ABAQUS提供了两种优化方法——拓扑优化和形状优化。拓扑优化(Topology optimization)通过分析过程中不断修改最初模型中指定优化区域的单元材料性质,有效地从分析的模型中移走/增加单元而获得最优的设计目标。形状优化(Shape optimization)则是在分析中对指定的优化区域不断移动表面节点从而达到减小局部应力集中的优化目标。拓扑优化和形状优化均遵从一系列优化目标和约束。 最优化方法(Optimization)是一个通过自动化程序增加设计者在经验和直觉从而缩短研发过程的工具。想要优化模型,必须知道如何去优化,仅仅说要减小应力或者增大特征值是不够,做优化必须有更专门的描述。比方说,想要降低在两种不同载荷工况下的最大节点力,类似的还有,想要最大化前五阶特征值之和。这种最优化的目标称之为目标函数(Object Function)。另外,在优化过程中可以同时强制限定某些状态参量。例如,可以指定某节点的位移不超过一定的数值。这些强制性的指定措施叫做约束(Constraint)。 ABAQUS/CAE可以创建模型然后定义、配置和执行结构优化。更多信息请参考用户手册第十八章。 13.1.2 术语(Terminology) 设计区域(Design area): 设计区域即模型需要优化的区域。这个区域可以是整个模型,也可以是模型的一部分或者数部分。一定的边界条件、载荷及人为约束下,拓扑优化通过增加/删除区域中单元的材料达到最优化设计,而形状优化通过移动区域内节点来达到优化的目的。 设计变量(Design variables):设计变量即优化设计中需要改变的参数。拓扑优化中,设计区域中单元密度是设计变量,ABAQUS/CAE优化分析模块在其优化迭代过程中改变单元密度并将其耦合到刚度矩阵之中。实际上,拓扑优化将模型中单元移除的方法是将单元的质量和刚度充分变小从而使其不再参与整体结构响应。对于形状优化而言,设计变量是指设计区域内表面节点位移。优化时,ABAQUS或者将节点位置向外移动或者向内移动,抑或不移动。在此过程中,约束会影响表面节点移动的多少及其方向。优化仅仅直接修改边缘处的节点,而边缘内侧的节点位移通过边缘处节点插值得到。 设计循环(Design cycle): 优化分析是一种不断更新设计变量的迭代过程,执行ABAQUS进行模型修改、查看结果以及确定是否达到优化目的。其中每次迭代叫做一个设计循环。 优化任务(Optimization task):一次优化任务包含优化的定义,比如设计响应、目标、限制条件和几何约束。 设计响应(Design responses): 优化分析的输入量称之为设计响应。设计响应可以直接从ABAQUS的结果输出文件.odb中读取,比如刚度、应力、特征频率及位移等。或者ABAQUS 从结果文件中计算得到模型的设计响应,例如质心、重量、相对位移等。一个设计响应与模型紧密相关,然而,设计响应存在一定的范围,例如区域内的最大应力或者模型体积。另外,设计响应也与特点的分析步和载荷状况有关。 目标函数(Objective functions): 目标函数决定了优化的目标。一个目标函数是从设计响应中萃取的一定范围内的值,如最大位移和最大应力。一个目标函数可以用多个设计响应

如何利用ANSYS进行拓扑优化

如何利用ANSYS进行拓扑优化 前言 就目前而言,利用有限元进行优化主要分成两个阶段: (1)进行拓扑优化,明确零件最佳的外形、刚度、体积,或者合理的固有频率,主要目的是确定优化的方向; (2)进行尺寸优化,主要目的是确定优化后的的零件具体尺寸值,通常是在完成拓扑优化之后,再执行尺寸优化。 在ANSYS中,利用拓扑优化,可以完成以下两个目的: (1)在特定载荷和约束的条件下,确定零件的最佳外形,或者最小的体积(或者质量); (2)利用拓扑优化,使零件达到需要的固有频率,避免在使用过程中产生共振等不利影响。 本文主要就在ANSYS环境中如何执行拓扑优化进行说明。

1、利用ANSYS进行拓扑优化的过程 在ANSYS中,执行优化,通常分为以下6个步骤: 、定义需要求解的结构问题 对于结构进行优化分析,定义结构的物理特性必不可少,例如,需要定义结构的杨氏模量、泊松比(其值在~之间)、密度等相关的结构特性方面的信息,以供结构计算能够正常执行下去。

、选择合理的优化单元类型 在ANSYS中,不是所有的单元类型都可以执行优化的,必须满足如下的规定: (1)2D平面单元:PLANE82单元和PLANE183单元; (2)3D实体单元:SOLID92单元和SOLID95单元; (3)壳单元:SHELL93单元。 上述单元的特性在帮助文件中有详细的说明,同时对于2D单元,应使用平面应力或者轴对称的单元选项。 、指定优化和非优化的区域 在ANSYS中规定,单元类型编号为1的单元,才执行优化计算;否则,就不执行优化计算。例如,对于结构分析中,对于不能去除的部分区域将单元类型编号设定为≥2,就可以不执行优化计算,请见下面的代码片段:…… …… Et,1,solid92 Et,2,solid92 …… Type,1 Vsel,s,num,,1,2 Vmesh,all …… Type,2 Vsel,s,num,,3 Vmesh,all ……

机械结构拓扑优化设计研究现状及其发展趋势

机械结构拓扑优化设计研究现状及其发展趋势 发表时间:2018-12-27T16:17:28.400Z 来源:《河南电力》2018年13期作者:谢进芳 [导读] 机械产品应用范围相对较广,为确保机械产品在我国日常生活及企业从生产中得到有效应用,实施优化设计十分必要。 (广东科立工业技术股份有限公司广东省佛山市 528000) 摘要:随着现代科学技术的发展,市场产品竞争也越来越激烈,产品品种的换代速度加快,产品的复杂性在不断增加。所以产品生产正在以小批量、多品种的生产方式取代过去的单一品种大批量生产方式。而这种生产方式,肯定会缩短产品的生产周期,产品的成本也会降低,产品提高市场的占有率和竞争力也会提高。所以在机械结构设计中采用优化设计是满足市场竞争的需要。 关键词:机械结构拓扑;现状;发展趋势 引言 机械产品应用范围相对较广,为确保机械产品在我国日常生活及企业从生产中得到有效应用,实施优化设计十分必要。目前我国已经针对机械结构优化设计进行了研究,并取得一定成果,主要表现在船舶行业、焊工航天以及汽车行业等。机械结构的优化设计可有效提高其产品性能并增加其自身市场竞争力,对其市场发展起重要作用。 1.机械结构优化设计 随着科学技术的发展,机械产品更新换代的速度越来越快。过去,机械产品主要是大批量生产,产品相对单一。目前采用的是小批量加工方式,以保证产品的多样性。为了保证生产企业的利润,必须在保证质量的前提下,缩短生产周期,降低生产成本。优化设计能够达到上述目标,在一定程度上缩短了生产时间,降低了成本,有效地抢占了市场。机械结构优化设计已广泛应用于造船、运输、航空航天、冶金、纺织、建筑等领域。 机械结构优化设计流程主要包括:(1)针对所优化机械产品尽心目标函数优化设计,可确保机械产品相关技术指标符合优化要求。(2)设计机械产品优化函数变量,变量设计包括机械产品长度、厚度以及弧度等相关结构参数。(3)对机械产品优化设计约束条件进行设定,对计算过程中各项变量浮动范围进行限定。(4)通过以上步骤得出多种优化设计方案,分别对不同方案进行评价,根据机械结构优化设计需求选择最佳方案实施。 2.机械结构拓扑优化设计常用方法 (1)均匀化方法 常用的连续结构拓扑优化设计方法主要有均匀化方法、变密度方法、水平集方法以及进化结构优化方法等。 均匀化方法属于材料描述方式,基本思想是将微结构模型引入结构拓扑优化设计领域,以微结构的单胞尺寸参数为设计变量,根据单胞尺寸的变化实现微结构的增删,优化实体与孔的分布形成带孔洞的板,达到结构拓扑优化的目的。优化过程:①设计区域的划分;②确定设计变量;③进行拓扑优化设计;④以不同的微结构形式的分布显示连续结构的形状和拓扑状态。 图1 微结构单胞示意图 微结构的划分形式通常有空孔、实体和开孔 3种,空孔是指没有材料的微结构,其孔的尺寸为 1;实体是指具有各向同性材料的微结构,其孔的尺寸为 0;开孔是指具有正交各向异性材料的微结构,其孔的尺寸介于 0~1 且可变化。设计区域划分为空孔、实体和开孔的微结构形式。简单的二维微结构单胞示意图如图 1 所示。微结构上孔的尺寸和方位角是设计变量,其中孔的尺寸是微结构材料主方向,它可以由坐标转换矩阵体现在材料的有效弹性模量上,通过微结构的密度与有效弹性模量之间的关系曲线,把设计变量与结构各处的形态联结起来。在结构拓扑优化设计过程中,微结构中孔的尺寸和在 0~1 的变化区域就可使各微结构在空孔与实体之间变化,这样就可用连续变量对结构优化设计问题进行描述。 均匀化结构拓扑优化方法涉及的设计变量非常多,用的较多的优化算法是准则优化算法。 (2)变密度方法 变密度方法式是引入一种假想的密度在 0~1可变的材料,采用材料的密度作为优化设计变量,实现结构的拓扑变化;材料弹性模量等物理参数与材料密度间的关系也是人为假定的;这样不但将结构的拓扑优化问题转换为材料的最优分布问题,还可使优化结果尽可能具有非 0 即 1 的密度分布。变密度结构拓扑优化方法与采用尺寸变量相比,它更能反映拓扑优化的本质特征。因此,在实际工程的结构优化设计中大多采用变密度方法来解决结构优化问题。变密度结构拓扑优化方法常用的插值模型是固体各向同性惩罚微结构模型(SIMP)。由于变密度结构拓扑优化方法更能反映拓扑优化的本质特征,且概念简单、设计变量数目少,简化了计算求解过程,因此,变密度结构拓扑优化方法成为目前最常用的、也是用的最多的结构优化设计方法。 3.机械结构优化的应用趋势 随着优化方法的不断发展和完善,结构优化设计也逐渐发展起来。近年来,在结构优化算法方面,由于结构优化设计中变量较多,结构优化设计往往采用接近实际情况的复杂结构模型来模拟一些大型结构系统。因此,新的准则优化方法备受关注,但如何为一些特殊结构

拓扑优化

结构拓扑优化设计现状及前景 目前, 最优化设计理论和方法在机械结构设计中得到了深入的研究和广泛的应用。所谓优化设计就是根据具体的实际问题建立其优化设计的数学模型, 并采用一定的最优化方法寻找既满足约束条件又使目标函数最优的设计方案。根据优化问题的初始设计条件, 目前结构优化技术有四大领域: 1) 尺寸优化; 2) 形状优化; 3) 拓扑与布局优化; 4) 结构类型优化。结构尺寸优化是在结构的拓扑确定的前提下, 首先用少量尺寸对结构的某些变动进行表达, 如桁架各单元的横截面尺寸、某些节点位置的变动等, 然后在此基础上建立基于这些尺寸参数的数学模型并采用优化方法对该模型进行求解得到最优的尺寸参数。在尺寸优化设计中, 不改变结构的拓扑形态和边界形状, 只是对特定的尺寸进行调整, 相当于在设计初始条件中就增加了拓扑形态的约束。而结构最初始的拓扑形态和边界形状必须由设计者根据经验或实验确定, 而不能保证这些最初的设计是最优的, 所以最后得到的并不是全局最优的结果。结构形状优化是指在给定的结构拓扑前提下, 通过调整结构内外边界形状来改善结构的性能。以轴对称零件的圆角过渡形状设计的例子。形状设计对边界形状的改变没有约束,和尺寸优化相比其初始的条件得到了一定的放宽,应用的范围也得到了进一步的扩展。拓扑优化设计是在给定材料品质和设计域内,通过优化设计方法可得到满足约束条件又使目标函数最优的结构布局形式及构件尺寸。拓扑设计的初始约束条件更少, 设计者只需要提出设计域而不需要知道具体的结构拓扑形态。拓扑设计方法是一种创新性

的设计方法, 能为我们提供一些新颖的结构拓扑。目前, 拓扑设计理论在柔性受力结构、MEMS 器件及其它柔性微操作机构的设计中得到了广泛的研究。 结构拓扑优化的发展概况 结构拓扑优化包括离散结构的拓扑优化和连续变量结构的拓扑优化。近10 年来, 结构拓扑优化设计虽然取得了一些进展, 但大部分是针对连续变量的, 关于离散变量的研究为数甚少。由于离散变量优化的目标函数和约束函数是不连续、不可微的, 可行域退化为不连通的可行集, 所以难度远大于连续变量优化问题。在离散结构中, 桁架在工程中的应用较为广泛, 由于其重要性, 也由于其分析比较简单, 桁架结构的拓扑优化在文献中研究得最多. 结构拓扑优化的历史可以追溯到1904 年Michell提出的桁架理论, 但这一理论只能用于单工况并依赖于选择适当的应变场, 不能应用于工程实际。1964 年Dorn、Gomory、Greenberg 等人提出基结构法( ground structure approach) , 将数值方法引入该领域, 此后拓扑优化的研究重新活跃起来, 陆续有一些解析和数值方面的理论被 提出来。所谓基结构就是一个由结构节点、荷载作用点和支承点组成的节点集合, 集合中所有节点之间用杆件相连的结构。该方法的基本思路是: 从基结构的模型出发, 应用优化算法( 数学规划法或准则法) , 按照某种规划或约束, 将一些不必要的杆件从基结构中删除, 例如截面积达到零或下限的杆件将被删掉, 并认为最终剩下的杆件 决定了结构的最优拓扑。因此应用基结构, 可以将桁架拓扑优化当作

拓扑优化技术

拓扑优化技术 第1节基本知识 一、拓扑优化的概念 拓扑优化是指形状优化,有时也称为外型优化。拓扑优化的目标是寻找承受单载荷或多载荷的物体的最佳材料分配方案。这种方案在拓扑优化中表现为“最大刚度”设计。 与传统的优化设计不同的是,拓扑优化不需要给出参数和优化变量的定义。目标函数、状态变量和设计变量都是预定义好的。用户只需要给出结构的参数(材料特性、模型、载荷等)和要省去的材料百分比。 拓扑优化的目标—目标函数—是在满足结构的约束(V)情况下减少结构的变形能。减小结构的变形能相当于提高结构的刚度。这个技术通过使用设计变量( i)给每个有限元的单元赋予内部伪密度来实现。这些伪密度用PLNSOL,TOPO命令来绘出。 ANSYS提供的拓扑优化技术主要用于确定系统的最佳几何形状,其原理是系统材料发挥最大利用率,同时确保系统的整体刚度(静力分析)、自振频率(模态分析)在满足工程要求的条件下获得极大或极小值。 拓扑优化应用场合:线性静力分析和模态分析。 拓扑优化原理:满足结构体积缩减量的条件下使目标函数结构柔量能量(the enery of structure compliance—SCOMP)的极小化。结构柔量能量极小化就是要求结构刚度的最大化。 例如,给定V=60表示在给定载荷并满足最大刚度准则要求的情况下省去60%的材料。图19-1表示满足约束和载荷要求的拓扑优化结果。图19-1a表示载荷和边界条件,图19-b 表示以密度云图形式绘制的拓扑结果。 图19-1 体积减少60%的拓扑优化示例 二、拓扑优化的基本过程 拓扑优化的基本步骤如下:

1.定义结构问题定义材料弹性模量、泊松系数、材料密度。 2.选择单元类型拓扑优化功能中的模型只能采用下列单元类型: ● 二维实体单元:Plane2和Plane82,用于平面应力问题和轴对称问题。 ● 三维实体单元:Solid92、Solid95。 ● 壳单元:SHELL93。 3.指定优化和不优化区域ANSYS只对单元类型编号为1的单元网格部分进行拓扑优 化,而对单元类型编号大于1的单元网格部分不进行拓扑优化,因此,拓扑优化时要确保进行拓扑优化区域单元类型编号为1,而不进行拓扑优化区域单元类型编号大于1即可。 4.定义并控制载荷工况或频率提取可以在单个载荷工况和多个载荷工况下做拓扑优化,单载荷工况是最简便的。 要在几个独立的载荷工况中得到优化结果时,必须用到写载荷工况和求解功能。在定义完每个载荷工况后,要用LSWRITE命令将数据写入文件,然后用LSSOLVE命令求解载荷工况的集合。 5.定义和控制优化过程拓扑优化过程包括定义优化参数和进行拓扑优化两个部分。用户可以用两种方式运行拓扑优化:控制并执行每一次迭代或自动进行多次迭代。 ANSYS有三个命令定义和执行拓扑优化:TOPDEF,TOPEXE和TOPITER。TOPDEF 命令定义要省去材料的量,要处理载荷工况的数目,收敛的公差;TOPEXE命令执行一次优化迭代;TOPITER命令执行多次优化迭代。 (1)定义优化参数首先要定义优化参数。用户要定义要省去材料的百分比,要处理载荷工况的数目,收敛的公差。 命令:TOPDEF GUI:Main Menu>Solution>Solve>Topological opt 注:本步所定义的内容并不存入ANSYS数据库中,因此在下一个拓扑优化中要重新使用TOPDEF命令。 (2)执行单次迭代定义好优化参数以后,可以执行一次迭代。迭代后用户可以查看收敛情况并绘出或列出当前的拓扑优化结果。可以继续做迭代直到满足要求为止。如果是在GUI方式下执行,在Topological Optimization 对话框(ITER域)中选择一次迭代。 命令:TOPEXE GUI:Main Menu>Solution>Solve>Topological opt TOPEXE的主要优点是用户可以设计自己的迭代宏进行自动优化循环和绘图。在下一节,可以看到TOPITER命令是一个ANSYS的宏,用来执行多次优化迭代。 (3)自动执行多次迭代 在定义好优化参数以后,用户可以自动执行多次迭代。在迭代完成以后,可以查看收敛情况并绘出或列出当前拓扑形状。如果需要的话,可以继续执行求解和迭代。TOPITER 命令实际是一个ANSYS的宏,可以拷贝和定制。

拓扑优化经典99行程序解读

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

连续体结构拓扑优化方法及存在问题分析(最新版)

( 安全管理 ) 单位:_________________________ 姓名:_________________________ 日期:_________________________ 精品文档 / 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进行了基于均匀化方法结构

结构拓扑优化的发展现状及未来说课讲解

结构拓扑优化的发展现状及未来

结构拓扑优化的发展现状及未来 王超 中国北方车辆研究所一、历史及发展概况 结构拓扑优化是近20年来从结构优化研究中派生出来的新分支,它在计算结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化的工程应用研究还很不成熟,在国外处在发展的初期,尤其在国内尚属于起步阶段。1904 年Michell在桁架理论中首次提出了拓扑优化的概念。自1964 年Dorn等人提出基结构法,将数值方法引入拓扑优化领域,拓扑优化研究开始活跃。20 世纪80 年代初,程耿东和N. Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题,他们开创性的工作引起了众多学者的研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计,开创了连续体结构拓扑优化设计研究的新局面。1993年Xie.Y.M和Steven.G.P 提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002 年罗鹰等提出三角网格进化法,该方法在优化过程中实现了退化和进化的统一,提高了优化效率。 二、拓扑优化的工程背景及基本原理 通常把结构优化按设计变量的类型划分成三个层次:结构尺寸优化、形状优化和拓扑优化。尺寸优化和形状优化已得到充分的发展,但它们存在着不能变更结构拓扑的缺陷。在这样的背景下,人们开始研究拓扑优化。拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料的分布问题。寻求一个最佳的拓扑结构形式有两种基本的原理:一种是退化

原理,另一种是进化原理。退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上,然后构造适当的优化模型,通过一定的优化方法逐步删减那些不必要的结构元素,直至最终得到一个最优化的拓扑结构形式。进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化,它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。 三、结构拓扑优化设计方法 目前常使用的拓扑优化设计方法可以分为两大类:退化法和进化法。 退化法即传统的拓扑优化方法,一般通过求目标函数导数的零点或一系列迭代计算过程求最优的拓扑结构。目前常用于拓扑优化的退化法有基结构方法、均匀化方法、变密度法、变厚度法等。 基结构方法(GSA)的思路是假定对于给定的桁架节点,在每两个节点之间用杆件连结起来得到的结构称为基结构。按照某种规则或约束,将一些不必要的杆件从基本结构中删除,认为最终剩下的构件决定了结构的最佳拓扑。基结构方法更适合于桁架和框架结构的拓扑优化。基结构法是在有限的子空间内寻优,容易丢失最优解,另外还存在组合爆炸、解的奇异性等问题。 均匀化方法(HA)引入微结构的单胞,通过优化计算确定其材料密度分布,并由此得出最优的拓扑结构。均匀化方法主要应用于连续体的拓扑优化设计,它不仅能用于应力约束和位移约束,也能用于频率约束。目前用均匀化方法来进行拓扑优化设计的有一般弹性问题、热传导问题、周期渐进可展曲面问题、非线性热弹性问题、振动问题和骨改造问题等。 变密度法是一种比较流行的力学建模方式,与采用尺寸变量相比,它更能反映拓扑优化的本质特征。通常,单元密度与弹性模量之间的关系采用人为给出的

连续体结构拓扑优化方法及存在问题分析

编号: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进行了基于均匀化方法结

拓扑优化简介

拓扑优化 什么是拓扑优化? 拓扑优化是指形状优化,有时也称为外型优化。拓扑优化的目标是寻找承受单载荷或多载荷的物体的最佳材料分配方案。这种方案在拓扑优化中表现为“最大刚度”设计。 与传统的优化设计不同的是,拓扑优化不需要给出参数和优化变量的定义。目标函数、状态变量和设计变量(参见“优化设计”一章)都是预定义好的。用户只需要给出结构的参数(材料特性、模型、载荷等)和要省去的材料百分比。 拓扑优化的目标——目标函数——是在满足结构的约束(V)情况下减少结构的变形能。减小结构的变形能相当于提高结构的刚度。这个技术通过使用设计变量( i)给每个有限元的单元赋予内部伪密度来实现。这些伪密度用PLNSOL,TOPO命令来绘出。 例如,给定V=60表示在给定载荷并满足最大刚度准则要求的情况下省去60%的材料。图2-1表示满足约束和载荷要求的拓扑优化结果。图2-1a表示载荷和边界条件,图2-2b表示以密度云图形式绘制的拓扑结果。 图2-1 体积减少60%的拓扑优化示例 如何做拓扑优化 拓扑优化包括如下主要步骤: 1.定义拓扑优化问题。 2.选择单元类型。 3.指定要优化和不优化的区域。 4.定义和控制载荷工况。 5.定义和控制优化过程。 6.查看结果。 拓扑优化的细节在下面给出。关于批处理方式和图形菜单方式 不同的做法也同样提及。 定义拓扑优化问题 定义拓扑优化问题同定义其他线性,弹性结构问题做法一样。用户需要定义材料特性(杨氏模量和泊松比),选择合适的单元类型生成有限元模型,施加载

荷和边界条件做单载荷步或多载荷步分析。参见“ANSYS Analysis Procedures Guides”第一、二章。 选择单元类型 拓扑优化功能可以使用二维平面单元,三维块单元和壳单元。要使用这个功能,模型中只能有下列单元类型: 二维实体单元:SOLID2和SOLID82 三维实体单元:SOLID92和SOLID95 壳单元:SHELL93 二维单元用于平面应力问题。 指定要优化和不优化的区域 只有单元类型号为1的单元才能做拓扑优化。可以使用这种限制控制模型优化和不优化的部分。例如,如果要保留接近圆孔部分或支架部分的材料,将这部分单元类型号指定为2或更大即可: … ET,1,SOLID92 ET,2,SOLID92 … TYPE,1 VSEL,S,NUM,,1,,2 !用这些单元划分的实体将被优化 VMESH,ALL TYPE,2 VSEL,S,NUM,,3 !用这些单元划分的实体将保持原状 VMESH,ALL … 用户可以使用ANSYS的选择和修改命令控制单元划分和类型号定义。 定义和控制载荷工况 可以在单个载荷工况和多个载荷工况下做拓扑优化。单载荷工况是最简便的。 要在几个独立的载荷工况中得到优化结果时,必须用到写载荷工况和求解功能。在定义完每个载荷工况后,要用LSWRITE命令将数据写入文件,然后用LSSOLVE命令求解载荷工况的集合。 例如,下面的输入演示如何将三个载荷工况联合做一个拓扑优化分析。 … D,10,ALL,0,,20,1 !定义第一个载荷工况的约束和载荷 NSEL,S,LOC,Y,0 SF, ALLSEL LSWRITE,1 !写第一个载荷工况 DDEL, SFDEL, NSEL,S,LOC,X,0,1 D,ALL,ALL,0

相关文档
相关文档 最新文档