《2022年2022年连续体结构拓扑优化方法评述 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年连续体结构拓扑优化方法评述 .pdf（12页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第 2 卷 第 1 期2011 年 2 月航空工程进展ADV AN CES IN AERON AU TICA L SCIENCE A ND EN GIN EERIN GVol 2N o 1Feb2011收稿日期 : 2010 12 01;修回日期 : 2011 01 20基金项目 : 教育部长江学者创新团队项目( Irt0906)通信作者 : 姚卫星 , wxyao nuaa. edu. cn文章编号 :1674 8190( 2011) 01001 12连续体结构拓扑优化方法评述夏天翔 , 姚卫星(南京航空 航天大学 飞行器先进设计技术国防重点学科实验室, 南京210016)摘要: 连续体结构
2、拓扑优化在优化中能产生新的构型,对实现自动化 智能结构设 计具有重要 意义。目前 , 连续体结构拓扑优化方法主要有: 均匀化方法、变厚度法、变密度法、渐进结构优化方法、 水平集法、独立连 续映射方法。本文首先系统回顾了以上方法的发展历程, 介绍了它们的研究现状。其次, 通过对比以 上拓扑优 化方法对若干典型算例的优化结果, 表明以上方法都有较好的减重效果。最后 , 对以 上方法 进行了 总结 , 列出了 它们的优缺点和发展方向。关键词 : 拓扑优化 ; 均匀化方法 ; 变厚 度法 ;变密度法 ; 渐进结构优化方法; 水平集法 ; 独立连续映射方法中图分类号 : V 211. 7文献标识码 :
3、AA Survey of Topology Optimizationof Continuum StructureXia Tianx iang, Yao Weix ing( Key Laboratory of Fundamental Science for National Defenseadvanced Design Technology of FlightV ehicle, Nanjing University of Aeronautics and Astronautics, N anjing 210016, China)Abstract: A s the topology optim iz
4、ation of continuum structure can producenew configurations during the optimization, it is significant for automatic structure design. At present, the most commonly used topology optimization methods of continuum structur eare: the homogenization method, variablethickness method, variable density met
5、ho d, evolutio nary str uctur al optimizatio n method, level set method,independent continuous mappingmethod. Firstly,the development processes of above metho ds are systematically reviewed, their current r esearch is briefly intro ducedin this paper. Then, these methods are comparedand discussedthr
6、ough a number oftypical examples. The typical examples show that all of above methods have greatabilities to reduceweight. Finally , the advantages, disadvantages and development directio ns of above methods are discussed.Key words: topology optimization;homogenization metho d; variable thickness me
7、thod; variable density method;evolutionar y structure optimizationmethod; level set method; independent continuo us mapping method0引言按照设计变量的不同, 结构优化可分为以下三个层次 : 尺寸优化、 形状优化和拓扑优化。结构拓扑优化能在给定的外载荷和边界条件下, 通过改变结构拓扑使结构在满足约束的前提下性能达到最优。与尺寸优化、 形状优化相比, 结构拓扑优化的经济效果更为明显, 在优化中能产生新的构型, 是结构实现自动化智能设计所必不可少的。按照优化对象的性质,
8、拓扑优化可分为离散体拓扑优化和连续体拓扑优化两种。连续体拓扑优化与离散体拓扑优化相比 , 在 应用范围更广的同时, 模型描述困难 , 设计变量多 , 计算量大。在过去很长一段时间里, 连续体拓扑优化发展得十分缓慢, 直到 1988 年 Bendsoe等人 1提出均匀化方法之后 , 它才得到了迅速发展。目前, 国内外学者对结构拓扑优化问题已经进行了大量研究 2 9。目前最常用的连续体拓扑优化方法有均匀化方法、变厚度法、 变密度法、渐进结构优化方法(ESO) 、 水平集法( Level set) 、 独立连续映射方法(ICM) 等。从拓名师资料总结 - - -精品资料欢迎下载 - - - - -
9、- - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 12 页 - - - - - - - - - 扑优化方法的基本思路来看 , 可以 将它们分为两类:( 1) 改变优化对象的材料特性, 主要包括均匀化方法和变密度法。均匀化方法将结构变成多孔材料, 而变密度法改变了结构的密度。( 2) 改变优化对象的几何形状, 主要包括变厚度法、ICM 法、ESO 法、水平集法。变厚度法改变平面单元的厚度,而 ICM 法和 ESO 法通过删除或者增加单元。近年来 , 拓扑优化在飞机设计领域应用日益广泛 10 14。本文首先介绍了以上几种优化方法的基本模型
10、、改进模型和应用, 然后用这些方法对三个典型算例进行了优化和分析比较, 并进行了讨论。1拓扑优化方法回顾1. 1均匀化方法1988 年 , Bendsoe 等人提出了均匀化方法, 虽然最初此方法在作者看来是一种形状优化方法。其基本思想是在拓扑优化中引入所谓微结构, 该方法以空孔尺寸为设计变量, 以空孔尺寸的消长实现微结构的增删 , 从而改变结构拓扑。均匀化方法含有矩形或正方形孔的正方形单元, 如图 1 所示。!?图 1四种微结构形式Fig. 1Four kinds of micr ostructureBendsoe 给出了以结构柔度最小为目标的拓扑优化模型find ai( i= 1, 2, #
11、 , n)minCs. t. V ? Vmax( 1)式中 : ai为微结构 中正方形孔的边长 ( 如图 1 所示) ; C 为结构柔度 ; V 为结构实际体积; Vmax为给定的结构最大体积。Bendsoe 还对比分析了的四种微结构形式, 结果发现 : 增加变量可以提高计算精度,但不能明显减小优化结果的柔度;增加孔的旋转角度变量可以在一定程度上减小中间值的存在。因此,为了保证计算精度往往最多为每个单元设置三个变量。由于 Bendsoe 使用了数学规划法求解, 变量增多会造成求解更加复杂费时。运用以上的优化模型, Bendsoe 对受单一工况作用的平板进行了拓扑优化设计。之后, Diaz 等人
12、 15 16在均匀化法中引入权重因子, 解决了多 工况下的平面结构拓扑优化和频率优化问题。 Olhoff 17基于 CAD, 将拓扑、形状、尺寸优 化进行了整合。 Sigmund 18将均匀化方法运用到复合材料设计中 , 通过对微结构进行拓扑优化得到了有良好力学特性的复合材料。针对上述情况,一些学者对均匀化方法进行了改进。Suzuki 等人 19采用形式四的微结构形式 , 推导出 a、 b、 的迭代公式 , 然后使用优化准则法求解, 从而使模型求解更加简单。Bendsoe 等人 1运用平滑处理消除结果中的孔洞,使优化结果更为光滑。均匀化方法提出之后, 连续体拓扑优化得到了迅速发展。该方法理论严
13、谨、直观, 在早期连续体拓扑优化中起到了重要作用。近年来,国内外学者主要将其用于柔性机构拓扑优化设计 20 26和复合材料拓扑优化设计 27 29中。1. 2变厚度法变厚度法以单元厚度为设计变量,通过删除厚度处于尺寸下限的单元实现结构拓扑的变更。其优化模型如下 :find hi( i =1, 2, # , n)min V =%ni=1hisis.t.i?hi& (0, h0)( 2)式中 : hi为第 i 个单元的厚度 ; si为单元的面积 ;i为第 i 个单元的应力 ; h0为单元厚度的上限。基于变厚度法 , Tenek 等人 30对薄壳结构进行2航空工程进展第 2 卷名师资料总结 - -
14、-精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 12 页 - - - - - - - - - 了拓扑优化的研究 ; 程耿东等人31对平面膜结构进行了拓扑优化的研究;周克民等人 32对平面体拓扑优化进行了总结 , 提出用变厚度有限单元求解平面连续体拓扑优化问题;王健等人33研究了具有应力和厚度约束的平面结构拓扑优化问题。变厚度法避免了均匀化方法中构造微结构的麻烦, 因此可以较为方便地解决平面拓扑优化问题。但是 , 由于它把拓扑变量挂靠在低层次的单元厚度上 , 将连续体拓扑优化问题转化为广义尺寸优
15、化问题 , 因此无法运用于三维结构中。1. 3变密度法Mlejnek 等人 34根据均匀化方法提出了变密度法。其基本思想: 定义取值范围为 0, 1 的相对密度, 将优化目标用相对密度的显性函数表示,然后运用数学规划法或优化准则法求解。Mlejnek 给出了以结构柔度最小为目标, 以质量保留百分比为约束的变密度拓扑优化模型findi( i = 1, 2, # , n)min max1 ? j ? pCjs. t.%ni=1i!0V0im0? ?0 ?i? 1( 3)式中:i为第 i 个单元的相对密度; p 为工况总数 ;!0为单元原始密度; V0i为第 i 个单元的体积; m0为所有单元取原始
16、密度时的结构总质量, 可定义为结构原始质量; ? 定义为质量保留百分比, 是优化结果的结构质量与结构原始质量的比值。运用以上模型 , M lejnek 成功解决了多工况下平面和立体结构的拓扑优化问题。变密度法主要存在以下几点不足:( 1) 由于优化过程中不直接删除材料, 所以单元的相对密度在 0, 1 上均有分布。对于具有中等大小相对密度 ( 即中间密度 )的单元 , 是否删除就变得难以抉择 ;( 2) 以柔度最小为优化目标, 在解决含有强度和刚度约束的优化问题时不够方便;( 3) 对于优化结果, 需要人为设置阀值, 将相对密度小于阀值的单元删除。对于第一点不足,可通过人为给定相对密度与弹性模
17、量的插值函数关系来减少中间密度去解决。根据插值函数的类型分为SIM P 模型 35和 RAM P模型 36两种。 SIM P 模 型采用幂函数, 而 RAM P模型采用有理函数形式。大量工作证明,通过选取合适的插值函数 , 不仅可以有效减少中间密度的出现, 还可以简化优化模型, 减少计算时间。对于第二点不足 , Lin 37等人提出了一种体积约束适应算法( AVC) , 使体积约束在优化过程中可以根据强度条件等的满足情况自行改变。变密度法理论严密, Bendsoe 等人 38已经证实了其物理意义的存在,变密度法应用实例较多。早期, Yang 等人 39使用变密度法对发动机部件进行了拓扑优化。随
18、着Nastran 、 ANSYS等商用CAE软件开发出基于变密度法的拓扑优化模块, 设计人员 利 用 变 密 度法 进 行 拓 扑 优 化 的 算 例 日 益 增多 40 42,如夏天翔 43完成了复杂机翼接头的拓扑优化。1. 4独立连续映射方法(ICM 法)1998 年, 隋允康 44提出了 ICM 法。其基本思想: 使用在 0, 1 上连续变化的拓扑变量ti表征第 i号单元的有与无; 建立拓扑变量与约束的函数关系, 在每步迭代中求出所有单元的ti,然后依据某种删除准则将ti小于阀值的单元从结构中删除; 之后对改变后的结构进行下一步迭代分析直到收敛。以结构重量最小为目标的多工况应力约束下的拓
19、扑优化模型如下 45:find ti( i =1, 2, # n)min W =%ni=1t?iW0is. t.ij? t#i0i( j =1, 2, # p)0 ? ti? 1( 4)式中 :ij为 j 工况下第i 个单元的应力; p 为工况总数 ; ?、 #为常数。之后 , 隋允 康等人 46 49提出 了含位移和屈曲约束的 ICM 模型 ; 彭细荣等人 50提出频率约束下的 ICM 模型 ;叶红玲等人 51基于 ICM 法对三维连续体结构进行了拓扑优化设计 ; 隋允康等人 52 54根据独立连续拓扑变量概念, 建立了桁架和平面膜结构拓扑优化的有无复合体模型; 隋允康等人 55在 Nast
20、ran 及 Patran 的二次开发环境下, 对 位移3第 1 期夏天翔等 : 连续体结构拓扑优化方法评述名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 12 页 - - - - - - - - - 约束、 频率约束、位移及频率约束、 简谐载荷激励下动位移幅值约束等拓扑优化进行了研究, 证明了ICM 法在处理静力问题及动力问题的拓扑优化都是可行的。由于该方法对每个单元单独设置变量,隋允康 3提出运用对偶理论将模型进行转换从而降低求解规模, 提高求解效率。目前已有将ICM法
21、运用于工程的实例 56 57。1. 5渐进结构优化方法(ESO 法)1993 年 , Xie 等人 58提出了 ESO 法。其主要思想: 通过逐步将无效或者低效的材料从模型中删除, 使 结 构最 终 趋于 符 合一 定 工程 要求 的 优化结构。Xie 等人 59提出了解决以结构重量最小为目标,多工况应力约束下拓扑优化问题的优化步骤:( 1) 在给定的荷载和边界条件下, 定义设计区域,在其上划分有限元网格;( 2)在 每个 工 况 下, 对 结 构 分 别 进 行静 力分析 ;( 3) 求出每个工况下单元的Vo n M ises 的应力ij和最大的单元应力max j。若单元的Von M ise
22、s 应力在所有工况下均满足ijmaxj? RRn( 5)则认为该单元处于低应力状态, 可从结构中删除。其中 , RRn称为第 n 步迭代中的删除率。以上有限元分析和单元删除重复进行, 直到上面的不等式无法满足。然后引进另一个参数ER ,将删除率修改为 :RRn+ 1=RRn+ ER n =0, 1, 2, #( 6)( 4) 重复第 ( 2) ( 3)步, 直到某个工况下的最大应力达到许用应力。基于以 上 删 除 准 则 的 ESO 法 可 称 为 传 统ESO 法。 Zhou 等人 60提出传统ESO 法在某些情况下会得到非最优解。为解决这一问题, Xie 等人引入敏感数的概念, 敏感数表示
23、删除某个单元对约束或目标函数的影响程度。基于敏感数删除准则的 ESO 法成功解决了含应力、 刚度、 屈曲约束和自然频率最大的多工况拓扑优化问题 61 65。ESO 法的改进还有很多 , 如引入结构效率度量标准(PI) 的ESO 法 66、 逆 适 应 方 法 67、 双 向 ESO 法 ( BESO) 68 69、与 遗 传 算 法 相 结 合 的 GESO 70、与SIM P 模型相结合的ERPM 法 71、 基于单元材料属性更改的EPCM 72、 基于应变能的ESO 法 73、含有周长约束的BESO 法 74等。其中 , 最重要的改进方法是Q uerin 等人提出的BESO 法。 BESO
24、法的一次循环中 , 结构中不仅有单元被删除, 某些已经删除的单元也可以复原。通过对比可以发现,ESO 法与 BESO 法的优化结果相似, 但 后者效率更高。ESO 法原理简单 , 算法实现方便 , 物理概念明确, 与现有的有限元分析软件可以很好结合。因此,它的应用实例很多 75 79。然而 , 由于 ESO 法是一种演化方法 , 所以它的迭代次数较多, 计算效率较低。此外 , 虽然 T anskanen 80讨论了ESO 的理论基础 , 但是它的算法收敛性还没有得到证明。1. 6水平集法水平集的概念由 Osher 等人 81于 1988 年提出, 但直到 2000 年才由 Sethian 等人
25、 82运用到结构拓扑优化中。其基本思想: 将平面或三维结构表达为一个高一维水平集函数?(x, t) 的零水平集?( x, t) = 0, 如图 2 所示。图 2水平集函数Fig. 2Level set function上图中 , 左侧的阴影部分%是待优 化的平面结构 , 它内部含有一个孔洞%。其边界&( x, t) 在优化过程中不断变化。基于水平集的思想, 建立右侧所示的水平集函数?(x, t) , 它表示了一个比%高一维的外形随时间变化的三维实体。它具有这样的性质 : ?( x,t) =0 所表示的曲面在任何时刻都等于 %。Sethian 研究了以重量最轻为优化目标, 具有柔度约束的平面结构
26、拓扑优化问题。通过改变水平集函数使结构在边界应力大的地方向外扩张, 即增加材料 ; 在边界应力小的地方向外收缩, 即去除4航空工程进展第 2 卷名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 12 页 - - - - - - - - - 材料。之后 , Wang 等人 83和 Allaire等人 84对水平集拓扑优化进行了较为系统的研究和扩展。Wang证明了水平集拓扑优化的关键是计算一个含有速度场的 H amilton Jacobi 方程?( x, t)t+ ? ?(x,
27、 t)dxdt= 0( 7)而 Allaire引 入形状 导数 , 建立 了 H amiltonJacobi 方程同形状导数之间的关系,然后建立了与Wang 相同数学模型。以上的水平集法都归为传统水平集法。它有以下三点不足 :( 1) 对于二维问题时 , 结构在优化过程中不能自动产生新孔。对于这一现象, H uang 等人 76给出了解释 : 空间问题不同于平面问题, 由于表面的曲面运动 , 在优化过程中可以产生新的孔洞;( 2) 由于传统水平集法是一种边界演化方法,所以它的收敛速度往往较慢;( 3) 每步迭代结束时要重新计算水平集函数?( x,t) , 需要消耗大量计算时间。为了解决以上传统
28、水平集法的不足, 不少学者对其进行了改进 85 92。水平集法作为一种较为年轻的拓扑优化方法受到越来越多的重视。该方法原理简单, 过程形象直观, 优化结果的边界光滑清晰。由于水平集法诞生时间不长 , 故将其应用实例较少 93 94。2经典算例对比任何一种拓扑优化方法的提出, 都需要用算例来验证。下面给出了三个在文献中频繁出现的经典算例 , 对比了不同拓扑优化方法的表现。由于各拓扑优化方法的改进方法层出不穷, 为方便比较,仅选择最具代表性的文献中的算例。2. 1悬臂梁由于长宽比的不同会导致不同的优化结果, 所以悬臂 梁 可分 为 两种 : ( 1)短 悬 臂 梁 ( L ? H ) ;(2) 长
29、悬臂梁 ( L H ) 。下面就这两种情况分别进行讨论。2. 1. 1短悬臂梁短悬臂梁结构如图3 所示 , 其中 H ? 2L , 其右端面中点处受到竖直方向的集中力载荷作用。图 3短悬臂梁Fig. 3Short cantilever不同方法下的短悬臂梁优化后的拓扑结构如图 4 所示 19, 31, 34, 55, 68, 83。由于该问题较为简单, 所以各方法的优化结果都是与竖直方向夹45( 角且相互垂直的桁架结构。均匀化方法变厚度法变密度法ICM 法BESO 法水平集法图 4短悬臂梁优化结果Fig. 4Optim izations of short cantilever2. 1. 2长悬臂
30、梁长悬臂结构如图5 所示 , 其中 L ? 2H ,其右端面中点处受到竖直方向的集中力作用。图 5长悬臂梁Fig. 5L ong cantilever5第 1 期夏天翔等 : 连续体结构拓扑优化方法评述名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 12 页 - - - - - - - - - 不同方法给出的长悬臂梁的拓扑结构如图6所示 19, 32 , 55, 76, 89, 95, 优化结果的质量保留百分比列于如表 1 所示。均匀化方法( 50 步)变厚度法 ( 6
31、步)变密度法 ( 12 步)ICM 法( 32步)BES O法( 33 步)水平集法 ( 30 步)图 6长 悬臂梁优化结果Fig. 6Optimizations of long cantilever由上图可以看出,各优化方法给出的拓扑结构是相同的 , 但是均匀化方法有较严重的棋盘格式现象;水平集法的优化结果的边界很光滑, 这是因为它作为一种边界演化方法, 可以保证结构在优化过程中一直具有光滑的边界; 变厚度法中, 由于矩形单元在迭代中可退化为三角单元, 所以结果的边界也较为光滑。各方法的迭代步数如图6 所示 , 从中可以发现各方法的收敛速度已有明显不同。2. 2MBB 梁MBB 梁最早由德国
32、 MBB 公司为 解决 ) 空中客车?大型客机机舱地板梁的设计问题而提出, 问题描述如图7 所示。问题要求梁腹板上预留用于铺设电缆和管道等设备的孔洞。图 7M BB 梁Fig. 7M BB cantilever不同 方 法 的 优 化 结 果 17, 31, 61, 83, 96 , 97如 图 8所示。均匀化方法变厚度法变密度法ICM 法BESO法水平集法图 8MBB 梁优化结果Fig. 8Optimizations of MBB cantilev er由上图可以看出, 使用变密度法、 ICM 法和水平集法可以得到边界清晰、规则的结构; 而均匀化方法和 ESO 法所得的 结果有严重的棋盘格式
33、 问题; 变厚度法得到的拓扑结构不够规则, 而且其右下角与其他方法的结果相比缺少一部分材料。2. 3三个经典算例的对比不同方法优化结果的质量保留百分比如表1所示。表 1三个典型算例结果的质量保留百分比Table 1Percentagesof remaining mass of abovethr eetypical examples单位 : %算例短悬臂梁长悬臂梁MBB 梁均匀化方法255542变厚度法243660变密度法604040ICM 法283440ESO 法105025水平集法3040303分析讨论下面就基本思想、 优化模型、优化结果、方法的不足与改进等方面对上述六种优化方法的总结如表
34、2 所示。6航空工程进展第 2 卷名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 12 页 - - - - - - - - - 表 2优化方法总结Table 2Summary of above optimizatio n methods方法基本思想优化模型设计变量优化目标约束条件优化结果的拓扑特点均匀化方法优化过程中以空 孔尺 寸的 消长 实 现微结构的增删, 从而改变结构拓扑。微结构空孔尺寸与取向柔度最小体积约束含有大量孔洞变厚度法在迭代收敛后, 通过删 除厚度处于尺寸
35、下限的单元实现结构拓扑的变更。平面单元厚度体积最小应力约束只能是平面结构变密度法在迭代收敛后, 通过删 除相对密度低于某一阀值的单元来改变结构拓扑。单元相对密度柔度最小体积约束边界呈锯齿状; 会出现棋盘格式等数值不稳定问题ICM 法每步迭代中删除 拓扑 变量 小于 阀 值的单元 , 直到迭代收敛。单元拓扑变量重量最轻应力、位移、 屈曲、 频率等边界呈锯齿状; 会出现棋盘格式等数值不稳定问题ESO 法逐步将低效材料从结构中删除 , 使其趋于符合一定工程要求的优化结构。表征单元有无的状态变量多种目标应力、位移、 屈曲、 频率等边界呈锯齿状; 会出现棋盘格式等数值不稳定问题水平集法通过改变高一维 的
36、水 平集 函数 来 改变结构拓扑 , 直到符合一定工程要求。表征单元有无的状态变量多种目标柔度、 体积、位移等边界光滑 ; 对平面结构进行优化时难以产生新的孔洞由上表可以看出:( 1) 均匀化方法和变密度法的优化目标为柔度最小 , 这在通常以重量最轻为目标的结构优化设计中显得不够实用。ESO 法和水平集法适用于多种优化目标 , 如重量最轻、固有频率最大、 柔度最小等,因此它们的适用范围更广;( 2) 变密度法只能应用于平面结构的拓扑优化中, 这极大得限制了它的应用。将均匀化方法运用于三维结构设计时, 计算量会大幅增加, 因此它一般也只应用于平面结构;( 3) 由于 ESO 法和水 平集法 是边
37、界 演化方法,与其他几种方法相比,计算效率较低;( 4) 水平集法优化结果的结构边界相当光滑,其优化结果可以直接使用数控加工方法制造。对其经过进一步改进, 有可能实现自动化智能设计。而其他方法的优化结果的结构边界都呈现锯齿状。因此, 若要将它们的优化结构制造出来, 必须首先进行曲线或曲面拟合。目前 , 本文提及的各优化方法经过诸多学者的改进, 对简单规则的结构都可以给出很好的优化结构,然而它们在处理实际的工程问题时还存在很多不足。结合飞行器设计的具体要求, 提出几点连续体拓扑优化的发展方向:( 1) 优化中需进一步解决诸如棋盘格式、网格依赖、 局部最优解等问题 98 99;( 2) 飞行器结构
38、通常受力形式复杂, 有限元模型庞大 , 这样会导致拓扑优化的设计变量过多, 计算时间过长。因此, 针对连续体结构拓扑优化的快速分析方法, 如无网格法 100 102、 并行计算等需要进一步发展 ;( 3) 拓扑优化需进一步考虑诸如可靠性、疲劳寿命、 制造成本、几何形状等约束;( 4) 复合材料在航空航天结构上的应用日益增多 , 然而目前的拓扑优化方法都没有考虑复合材料结构。因此 , 下一步需要发展能有效处理复合材料结构的拓扑优化方法。由于连续体结构拓扑优化应用领域广泛, 潜力巨大 , 国内外越来越多的学者投入到该领域的研究中。随着连续体结构拓扑优化的成熟,结构自动化智能设计必将成为现实。参考文
39、献 1Bendsoe MP, KikuchiN.Generating optimaltopologiesinstru ctural design using a homogenizationmeth od J . Computer M ethods,1988, 71:197 224. 2Eschenauer H A, OlhoffN.Topologyoptimizationof continu um structures:A r eview J .American Society of Mechanical Engin eers, 2001, 54( 4) :331 390. 3隋允康 ,叶红
40、玲 ,杜 家政 . 结 构拓 扑优 化 的发 展 及其 模 型转 化 为 独 立 层 次 的 迫 切 性 J .工 程 力 学,2005 ,22:107 118.Sui Yunkang, Ye H ongling ,Du Jiazheng.Developm ent ofstructuraltopologicaloptimiz ation and imm inen cy of its m odeltrans formationintoindependent level J .Engineerin g M echanics, 2001,54( 4) : 331 390. ( in Chin ese
41、) 4Rozvan y G I N.Aims,scope, meth ods, historyand unifiedterminologyof computer aided topologyoptimizationin structural m echanics J .StructureM ultidisciplinaryOptimization,2001, 21: 90 108.7第 1 期夏天翔等 : 连续体结构拓扑优化方法评述名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7
42、页，共 12 页 - - - - - - - - - 5郭中泽 ,张卫红 , 陈裕泽 . 结构 拓扑优化设计综述 J . 机械设计 ,2007, 24( 8) : 1 5.Guo Zhongze, Zhang W eihong,Chen Yuze.An overview onth e topologicaloptimiz ation design of structu res J . Journalof M achine Design , 2007, 24( 8) : 1 5. ( in Chinese) 6周克民 ,李俊峰 , 李霞 . 结构拓 扑优化研究方法综述 J . 力学进展 ,200
43、5, 35( 1) : 69 76.Zhou Kem in, Li Junfeng,Li Xia.A review on topologyoptimiz ation of structu res J . Advances in M echanics, 2005, 35( 1) : 69 76. ( in Chin ese) 7罗震 ,陈立平 , 黄玉盈 , 等. 连续体结构的拓扑优化设计 J .力学进展 , 2004, 34( 4) : 463 476.Lu o Zhen,Chen Liping,Hu ang Yuying,et al.T opologicaloptimiz ation des
44、ign for continuu m stru ctures J . Advances inM echanics, 2004, 34( 4) :463 476. ( in Chinese) 8王振海 ,张卫红 . 柔结构拓扑优化设计发展概况 J . 机 械设计,2004, 21( 3) : 1 4.W ang Zhenhai,Zhang W eihong.T opological optim ization design for flex ible structure J .Journalof M achine Design,2004,21( 3) : 1 4.( in Chinese) 9张
45、卫红 ,王敏 . 拓扑优化技术在汽车工业的应用 J . 昆 明理工大学学报( 理工版 ) ,2005, 30( 1) : 77 81.Zhang Weihong, Wang Min.Applicationsof topologyoptim ization in the automobilein dustry J . Journal of Kunm ing Universityof Scienceand T echnology( Science and T echnology) ,2005,30( 1) : 77 81. ( in Chinese) 10邓扬晨 ,陈华 , 马明亮 , 等. 基于
46、拓 扑优化技术的飞机 普通框设计方法研究 J . 强度与环境, 2005,32( 2) : 39 45.Deng Yangchen,Chen Hua, M aM ingliang,et al. Studies ofaircraftframedesign b ased on topologyoptimization J .Structure & EnvironmentEngin eering, 2005, 32( 2) : 39 45.( in Chinese) 11W aqas S, Fan Y Q, Lu H,Z Qasim.T opology optimizationformanufact
47、urableconfigurationdesign of aircrafttailrib J . CADDM , 2009, 19( 1) :94 101. 12邓扬晨 ,刘晓欧 , 朱继宏 . 飞机加强框的一种结构拓扑优化设计方法 J . 飞机设计 , 2004,4: 11 16.Deng Yang chen, LiuXiaoou, ZhuJihong.Anapproach tostructuraltopologyoptimiz ationaboutaircraftreinforcedframes design J .AircraftDesign,2004,4:11 16.( inChine
48、se) 13Eschenauer H , Becker W,Schumacher A.M ultidisciplinarystructuraloptim izationinaircraftdesign R .BM BFBun desvorhaben DYNAFLEX.Germ an: 1998. 14 Krog L.T opology Optimizationof AircraftWingBox Rib s R . AIAA2004 4481, 2004. 15 Diaz A R,BendsoeM P. Shape optim ization of structures formultiple
49、loadingconditionsu sing a homogeniz ation method J . StructuralOptimiz ation,1992, 4:17 22. 16Diaz A R.Solution to shape and topologyeigenvalu e optim ization problem s using a homogenization m ethod J.InternationalJournal for Num erical M ethods in Engineer,1992, 35( 7) : 1487 1502. 17Olhoff N.On C
50、AD integratedstructu ral topology and designoptim ization com puter meth ods in applied mechanics and Engin eering J . Com puter Methods in AppliedMechanics andEngineering,1991, 89(1 3) :259 279. 18Sigmund O.Materialswith prescribed constitutiveparameters:an invers e homogeniz ation pr oblem J . Int