连续体结构拓扑优化方法评述_夏天翔

第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.360docs.net/doc/7715176331.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)等。从拓

扑优化方法的基本思路来看,可以将它们分为两类:(1)改变优化对象的材料特性,主要包括均匀化方法和变密度法。均匀化方法将结构变成多孔材料,而变密度法改变了结构的密度。(2)改变优化对象的几何形状,主要包括变厚度法、ICM 法、ESO 法、水平集法。变厚度法改变平面单元的厚度,而ICM 法和ESO 法通过删除或者增加单元。近年来,拓扑优化在飞机设计领域应用日益广泛[10-14]。本文首先介绍了以上几种优化方法的基本模型、改进模型和应用,然后用这些方法对三个典型算例进行了优化和分析比较,并进行了讨论。

1 拓扑优化方法回顾

1.1 均匀化方法

1988年,Bendsoe 等人提出了均匀化方法,虽

然最初此方法在作者看来是一种形状优化方法。其基本思想是在拓扑优化中引入所谓微结构,该方法以空孔尺寸为设计变量,以空孔尺寸的消长实现微结构的增删,从而改变结构拓扑。均匀化方法含有矩形或正方形孔的正方形单元,如图1

所示。

1

o

? ?

图1 四种微结构形式F ig.1 Four kinds of micr ostructure

Bendso e 给出了以结构柔度最小为目标的拓扑优化模型 find a i (i =1,2,,,n) min C

s.t .V [V max

(1)

式中:a i 为微结构中正方形孔的边长(如图1所示);C 为结构柔度;V 为结构实际体积;V max 为给

定的结构最大体积。

Bendsoe 还对比分析了的四种微结构形式,结果发现:增加变量可以提高计算精度,但不能明显减小优化结果的柔度;增加孔的旋转角度变量可以在一定程度上减小中间值的存在。因此,为了保证计算精度往往最多为每个单元设置三个变量。由于Bendsoe 使用了数学规划法求解,变量增多会造成求解更加复杂费时。

运用以上的优化模型,Bendsoe 对受单一工况作用的平板进行了拓扑优化设计。之后,Diaz 等人

[15-16]

在均匀化法中引入权重因子,解决了多工

况下的平面结构拓扑优化和频率优化问题。O-l ho ff [17]

基于CAD,将拓扑、形状、尺寸优化进行了

整合。Sigm und [18]

将均匀化方法运用到复合材料设计中,通过对微结构进行拓扑优化得到了有良好力学特性的复合材料。针对上述情况,一些学者对均匀化方法进行了改进。Suzuki 等人[19]

采用形式四的微结构形式,推导出a 、b 、H 的迭代公式,然后使用优化准则法求解,从而使模型求解更加简单。Bendsoe 等人[1]运用平滑处理消除结果中的孔洞,使优化结果更为光滑。

均匀化方法提出之后,连续体拓扑优化得到了迅速发展。该方法理论严谨、直观,在早期连续体拓扑优化中起到了重要作用。近年来,国内外学者

主要将其用于柔性机构拓扑优化设计[20-26]和复合

材料拓扑优化设计[27-29]

中。

1.2 变厚度法

变厚度法以单元厚度为设计变量,通过删除厚

度处于尺寸下限的单元实现结构拓扑的变更。其优化模型如下: find h i (i =1,2,,,n) min V =E n

i=1

h i s

i

s.t.

R i [

R

h i I (0,h 0)(2)

式中:h i 为第i 个单元的厚度;s i 为单元的面积;R i 为第i 个单元的应力;h 0为单元厚度的上限。

基于变厚度法,Tenek 等人[30]对薄壳结构进行

2航空工程进展 第2卷

了拓扑优化的研究;程耿东等人[31]对平面膜结构进行了拓扑优化的研究;周克民等人[32]对平面体拓扑优化进行了总结,提出用变厚度有限单元求解平面连续体拓扑优化问题;王健等人[33]研究了具有应力和厚度约束的平面结构拓扑优化问题。

变厚度法避免了均匀化方法中构造微结构的麻烦,因此可以较为方便地解决平面拓扑优化问题。但是,由于它把拓扑变量挂靠在低层次的单元厚度上,将连续体拓扑优化问题转化为广义尺寸优化问题,因此无法运用于三维结构中。

1.3变密度法

Mlejnek等人[34]根据均匀化方法提出了变密度法。其基本思想:定义取值范围为[0,1]的相对密度L,将优化目标用相对密度L的显性函数表示,然后运用数学规划法或优化准则法求解。

Mlejnek给出了以结构柔度最小为目标,以质量保留百分比为约束的变密度拓扑优化模型

find L i(i=1,2,,,n)

min m ax

1[j[p

C j

s.t.E n

i=1

L i Q0V0i

m0

[A

0[L i[1(3)

式中:L i为第i个单元的相对密度;p为工况总数; Q0为单元原始密度;V0i为第i个单元的体积;m0为所有单元取原始密度时的结构总质量,可定义为结构原始质量;A定义为质量保留百分比,是优化结果的结构质量与结构原始质量的比值。运用以上模型,M lejnek成功解决了多工况下平面和立体结构的拓扑优化问题。

变密度法主要存在以下几点不足:

(1)由于优化过程中不直接删除材料,所以单元的相对密度在[0,1]上均有分布。对于具有中等大小相对密度(即中间密度)的单元,是否删除就变得难以抉择;

(2)以柔度最小为优化目标,在解决含有强度和刚度约束的优化问题时不够方便;

(3)对于优化结果,需要人为设置阀值,将相对密度小于阀值的单元删除。

对于第一点不足,可通过人为给定相对密度与弹性模量的插值函数关系来减少中间密度去解决。根据插值函数的类型分为SIM P模型[35]和RAM P 模型[36]两种。SIM P模型采用幂函数,而RAM P 模型采用有理函数形式。大量工作证明,通过选取合适的插值函数,不仅可以有效减少中间密度的出现,还可以简化优化模型,减少计算时间。对于第二点不足,Lin[37]等人提出了一种体积约束适应算法(AVC),使体积约束在优化过程中可以根据强度条件等的满足情况自行改变。

变密度法理论严密,Bendsoe等人[38]已经证实了其物理意义的存在,变密度法应用实例较多。早期,Yang等人[39]使用变密度法对发动机部件进行了拓扑优化。随着Nastran、ANSYS等商用CAE 软件开发出基于变密度法的拓扑优化模块,设计人员利用变密度法进行拓扑优化的算例日益增多[40-42],如夏天翔[43]完成了复杂机翼接头的拓扑优化。

1.4独立连续映射方法(ICM法)

1998年,隋允康[44]提出了ICM法。其基本思想:使用在[0,1]上连续变化的拓扑变量t i表征第i 号单元的有与无;建立拓扑变量与约束的函数关系,在每步迭代中求出所有单元的t i,然后依据某种删除准则将t i小于阀值的单元从结构中删除;之后对改变后的结构进行下一步迭代分析直到收敛。以结构重量最小为目标的多工况应力约束下的拓扑优化模型如下[45]:

find t i(i=1,2,,n)

min W=E n

i=1

t A i W0i

s.t.R ij[t B i R0i(j=1,2,,p)

0[t i[1(4)

式中:R ij为j工况下第i个单元的应力;p为工况总数;A、B为常数。

之后,隋允康等人[46-49]提出了含位移和屈曲约束的ICM模型;彭细荣等人[50]提出频率约束下的ICM模型;叶红玲等人[51]基于ICM法对三维连续体结构进行了拓扑优化设计;隋允康等人[52-54]根据独立连续拓扑变量概念,建立了桁架和平面膜结构拓扑优化的有无复合体模型;隋允康等人[55]在Nastran及Patran的二次开发环境下,对位移

3

第1期夏天翔等:连续体结构拓扑优化方法评述

约束、频率约束、位移及频率约束、简谐载荷激励下动位移幅值约束等拓扑优化进行了研究,证明了ICM 法在处理静力问题及动力问题的拓扑优化都是可行的。由于该方法对每个单元单独设置变量,隋允康

[3]

提出运用对偶理论将模型进行转换从而

降低求解规模,提高求解效率。目前已有将ICM 法运用于工程的实例

[56-57]

。

1.5 渐进结构优化方法(ESO 法)

1993年,Xie 等人[58]提出了ESO 法。其主要思想:通过逐步将无效或者低效的材料从模型中删除,使结构最终趋于符合一定工程要求的优化结构。

Xie 等人

[59]

提出了解决以结构重量最小为目

标,多工况应力约束下拓扑优化问题的优化步骤:

(1)在给定的荷载和边界条件下,定义设计区域,在其上划分有限元网格;

(2)在每个工况下,对结构分别进行静力分析;

(3)求出每个工况下单元的Vo n M ises 的应力R ij 和最大的单元应力R max j 。若单元的Von M -i ses 应力在所有工况下均满足

R ij

R max j

[RR n

(5)

则认为该单元处于低应力状态,可从结构中删除。其中,RR n 称为第n 步迭代中的删除率。

以上有限元分析和单元删除重复进行,直到上面的不等式无法满足。然后引进另一个参数ER ,将删除率修改为:

RR n +1=RR n +ER n =0,1,2,,

(6)

(4)重复第(2)~(3)步,直到某个工况下的最大应力达到许用应力。

基于以上删除准则的ESO 法可称为传统ESO 法。Zhou 等人[60]提出传统ESO 法在某些情况下会得到非最优解。为解决这一问题,Xie 等人引入敏感数的概念,敏感数表示删除某个单元对约束或目标函数的影响程度。基于敏感数删除准则的ESO 法成功解决了含应力、刚度、屈曲约束和自

然频率最大的多工况拓扑优化问题

[61-65]

。ESO 法

的改进还有很多,如引入结构效率度量标准(PI)的ESO 法[66]

、逆适应方法[67]

、双向ESO 法(BE -SO)[68-69]、与遗传算法相结合的GESO [70]

、与SIM P 模型相结合的ERPM 法[71]

、基于单元材料属性更改的EPCM [72]、基于应变能的ESO 法[73]、含有周长约束的BESO 法[74]

等。其中,最重要的改进方法是Q uerin 等人提出的BESO 法。BESO 法的一次循环中,结构中不仅有单元被删除,某些已经删除的单元也可以复原。通过对比可以发现,ESO 法与BESO 法的优化结果相似,但后者效率更高。

ESO 法原理简单,算法实现方便,物理概念明确,与现有的有限元分析软件可以很好结合。因此,它的应用实例很多[75-79]

。然而,由于ESO 法是一种演化方法,所以它的迭代次数较多,计算效率较低。此外,虽然T anskanen [80]讨论了ESO 的理论基础,但是它的算法收敛性还没有得到证明。1.6 水平集法

水平集的概念由Osher 等人[81]于1988年提出,但直到2000年才由Sethian 等人[82]运用到结构拓扑优化中。其基本思想:将平面或三维结构表达为一个高一维水平集函数5(x ,t)的零水平集5(x ,t)=0,如图2所示。

图2 水平集函数Fig.2 L ev el set function

上图中,左侧的阴影部分8是待优化的平面结构,它内部含有一个孔洞8。其边界#(x ,t)在优化过程中不断变化。基于水平集的思想,建立右侧所示的水平集函数5(x ,t),它表示了一个比8高一维的外形随时间变化的三维实体。它具有这样的性质:5(x ,t)=0所表示的曲面在任何时刻都等于8。

Sethian 研究了以重量最轻为优化目标,具有柔度约束的平面结构拓扑优化问题。通过改变水平集函数使结构在边界应力大的地方向外扩张,即增加材料;在边界应力小的地方向外收缩,即去除

4航空工程进展 第2卷

材料。

之后,Wang等人[83]和Allaire等人[84]对水平集拓扑优化进行了较为系统的研究和扩展。Wang 证明了水平集拓扑优化的关键是计算一个含有速度场的H amilton-Jacobi方程

55(x,t)

5t+?5(x,t)d x

d t

=0(7)

而Allaire引入形状导数,建立了H amilton-Jaco bi方程同形状导数之间的关系,然后建立了与Wang相同数学模型。

以上的水平集法都归为传统水平集法。它有以下三点不足:

(1)对于二维问题时,结构在优化过程中不能自动产生新孔。对于这一现象,H uang等人[76]给出了解释:空间问题不同于平面问题,由于表面的曲面运动,在优化过程中可以产生新的孔洞;

(2)由于传统水平集法是一种边界演化方法,所以它的收敛速度往往较慢;

(3)每步迭代结束时要重新计算水平集函数5(x,t),需要消耗大量计算时间。

为了解决以上传统水平集法的不足,不少学者对其进行了改进[85-92]。

水平集法作为一种较为年轻的拓扑优化方法受到越来越多的重视。该方法原理简单,过程形象直观,优化结果的边界光滑清晰。由于水平集法诞生时间不长,故将其应用实例较少[93-94]。

2经典算例对比

任何一种拓扑优化方法的提出,都需要用算例来验证。下面给出了三个在文献中频繁出现的经典算例,对比了不同拓扑优化方法的表现。由于各拓扑优化方法的改进方法层出不穷,为方便比较,仅选择最具代表性的文献中的算例。

2.1悬臂梁

由于长宽比的不同会导致不同的优化结果,所以悬臂梁可分为两种:(1)短悬臂梁(L[H);

(2)长悬臂梁(L>H)。下面就这两种情况分别进行讨论。

2.1.1短悬臂梁

短悬臂梁结构如图3所示,其中H U2L,其右端面中点处受到竖直方向的集中力载荷作用。

图3短悬臂梁

Fig.3Sho rt cantilever

不同方法下的短悬臂梁优化后的拓扑结构如图4所示[19,31,34,55,68,83]。由于该问题较为简单,所以各方法的优化结果都是与竖直方向夹45b角且相互垂直的桁架结构。

均匀化方法变厚度法变密度法

ICM法BESO法水平集法

图4短悬臂梁优化结果

Fig.4Optim izat ions of short cantilever

2.1.2长悬臂梁

长悬臂结构如图5所示,其中L U2H,其右端面中点处受到竖直方向的集中力作用。

图5长悬臂梁

F ig.5L o ng cantilever

5

第1期夏天翔等:连续体结构拓扑优化方法评述

不同方法给出的长悬臂梁的拓扑结构如图6所示

[19,32,55,76,89,95]

,优化结果的质量保留百分比列

于如表1

所示。

均匀化方法(50步) 变厚度法(6步

)

变密度法(12步) IC M 法(32步

)

BES O 法(33步)

水平集法(30步)

图6 长悬臂梁优化结果

Fig.6 O pt imizatio ns of long cantilever

由上图可以看出,各优化方法给出的拓扑结构是相同的,但是均匀化方法有较严重的棋盘格式现象;水平集法的优化结果的边界很光滑,这是因为它作为一种边界演化方法,可以保证结构在优化过程中一直具有光滑的边界;变厚度法中,由于矩形单元在迭代中可退化为三角单元,所以结果的边界也较为光滑。各方法的迭代步数如图6所示,从中可以发现各方法的收敛速度已有明显不同。2.2 MBB 梁

MBB 梁最早由德国MBB 公司为解决/空中客车0大型客机机舱地板梁的设计问题而提出,问题描述如图7所示。问题要求梁腹板上预留用于

铺设电缆和管道等设备的孔洞。

图7 M BB 梁

F ig.7 M BB cantilever

不同方法的优化结果[17,31,61,83,96,97]如图8所示。

均匀化方法 变厚度法

变密度法 ICM 法

BE SO 法

水平集法

图8 M BB 梁优化结果

Fig.8 O ptimizations of M BB cantilev er

由上图可以看出,使用变密度法、ICM 法和水平集法可以得到边界清晰、规则的结构;而均匀化方法和ESO 法所得的结果有严重的棋盘格式问题;变厚度法得到的拓扑结构不够规则,而且其右下角与其他方法的结果相比缺少一部分材料。2.3 三个经典算例的对比

不同方法优化结果的质量保留百分比如表1所示。

表1 三个典型算例结果的质量保留百分比T able 1 Per centag es o f remaining mass o f abo ve

thr ee t ypical ex amples

单位:%算例短悬臂梁

长悬臂梁

M BB 梁均匀化方法255542变厚度法243660变密度法604040ICM 法283440ESO 法105025水平集法

30

40

30

3 分析讨论

下面就基本思想、优化模型、优化结果、方法的

不足与改进等方面对上述六种优化方法的总结如表2所示。

6航空工程进展 第2卷

表2优化方法总结

T able2Summar y of abov e optimizatio n met ho ds

方法基本思想

优化模型

设计变量优化目标约束条件

优化结果的拓扑特点

均匀化方法优化过程中以空孔尺寸的消长实现微

结构的增删,从而改变结构拓扑。

微结构空孔

尺寸与取向

柔度最小体积约束含有大量孔洞

变厚度法在迭代收敛后,通过删除厚度处于尺寸

下限的单元实现结构拓扑的变更。

平面单元厚度体积最小应力约束只能是平面结构

变密度法在迭代收敛后,通过删除相对密度低于

某一阀值的单元来改变结构拓扑。

单元相对密度柔度最小体积约束

边界呈锯齿状;会出现棋盘格式等

数值不稳定问题

ICM法每步迭代中删除拓扑变量小于阀值的

单元,直到迭代收敛。

单元拓扑变量重量最轻

应力、位移、屈

曲、频率等

边界呈锯齿状;会出现棋盘格式等

数值不稳定问题

ESO法逐步将低效材料从结构中删除,使其趋

于符合一定工程要求的优化结构。

表征单元有无

的状态变量

多种目标

应力、位移、屈

曲、频率等

边界呈锯齿状;会出现棋盘格式等

数值不稳定问题

水平集法通过改变高一维的水平集函数来改变

结构拓扑,直到符合一定工程要求。

表征单元有无

的状态变量

多种目标

柔度、体积、

位移等

边界光滑;对平面结构进行优化时

难以产生新的孔洞

由上表可以看出:

(1)均匀化方法和变密度法的优化目标为柔度最小,这在通常以重量最轻为目标的结构优化设计中显得不够实用。ESO法和水平集法适用于多种优化目标,如重量最轻、固有频率最大、柔度最小等,因此它们的适用范围更广;

(2)变密度法只能应用于平面结构的拓扑优化中,这极大得限制了它的应用。将均匀化方法运用于三维结构设计时,计算量会大幅增加,因此它一般也只应用于平面结构;

(3)由于ESO法和水平集法是边界演化方法,与其他几种方法相比,计算效率较低;

(4)水平集法优化结果的结构边界相当光滑,其优化结果可以直接使用数控加工方法制造。对其经过进一步改进,有可能实现自动化智能设计。而其他方法的优化结果的结构边界都呈现锯齿状。因此,若要将它们的优化结构制造出来,必须首先进行曲线或曲面拟合。

目前,本文提及的各优化方法经过诸多学者的改进,对简单规则的结构都可以给出很好的优化结构,然而它们在处理实际的工程问题时还存在很多不足。结合飞行器设计的具体要求,提出几点连续体拓扑优化的发展方向:

(1)优化中需进一步解决诸如棋盘格式、网格依赖、局部最优解等问题[98-99];

(2)飞行器结构通常受力形式复杂,有限元模型庞大,这样会导致拓扑优化的设计变量过多,计算时间过长。因此,针对连续体结构拓扑优化的快速分析方法,如无网格法[100-102]、并行计算等需要进一步发展;

(3)拓扑优化需进一步考虑诸如可靠性、疲劳寿命、制造成本、几何形状等约束;

(4)复合材料在航空航天结构上的应用日益增多,然而目前的拓扑优化方法都没有考虑复合材料结构。因此,下一步需要发展能有效处理复合材料结构的拓扑优化方法。

由于连续体结构拓扑优化应用领域广泛,潜力巨大,国内外越来越多的学者投入到该领域的研究中。随着连续体结构拓扑优化的成熟,结构自动化智能设计必将成为现实。

参考文献

[1]Bendsoe M P,Kikuchi N.Gen erating optimal topologies in

stru ctural design u sing a homogenization meth od[J].Comput-er M ethods,1988,71:197-224.

[2]Esch enauer H A,Olhoff N.Topology optimization of contin-

u um structures:A r eview[J].American S ociety of M echan-i cal Engin eers,2001,54(4):331-390.

[3]隋允康,叶红玲,杜家政.结构拓扑优化的发展及其模型

转化为独立层次的迫切性[J].工程力学,2005,22: 107-118.

Sui Yunkang,Ye H ongling,Du Jiazhen g.Developm ent of structural topological optimiz ation an d imm inen cy of its m odel trans formation into indep endent level[J].Engineerin g M e-chanics,2001,54(4):331-390.(in Chin ese)

[4]Rozvan y G I N.Aims,scope,meth od s,his tory and unified

terminology of compu ter-aided topology optimization in s truc-tural m echan ics[J].Structure M ultidisciplinary Optimiza-tion,2001,21:90-108.

7

第1期夏天翔等:连续体结构拓扑优化方法评述

[5]郭中泽,张卫红,陈裕泽.结构拓扑优化设计综述[J].机械

设计,2007,24(8):1-5.

Guo Zh ongze,Zhang W eihong,Chen Yuze.An overview on th e topological optimiz ation design of structu res[J].J ournal of M achine Design,2007,24(8):1-5.(in Ch ines e)

[6]周克民,李俊峰,李霞.结构拓扑优化研究方法综述[J].力

学进展,2005,35(1):69-76.

Zhou Kem in,Li Ju nfeng,Li Xia.A review on topology opt-i miz 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,Ch en Liping,Hu ang Yuying,et al.T opological optimiz ation design for continuu m stru ctures[J].Advances in M ech anics,2004,34(4):463-476.(in C hinese)

[8]王振海,张卫红.柔结构拓扑优化设计发展概况[J].机械设

计,2004,21(3):1-4.

W an g Zhenhai,Zhang W eihong.T opological optim ization de-s ign for flex ible s tructure[J].Journal of M achine Des ign, 2004,21(3):1-4.(in C hinese)

[9]张卫红,王敏.拓扑优化技术在汽车工业的应用[J].昆明理

工大学学报(理工版),2005,30(1):77-81.

Zhang Weihong,Wang M in.Applications of topology optim-i zation in the automobile in dustry[J].J ou rnal of Kunm ing U-niversity of Science and T echnology(Science and T echnology), 2005,30(1):77-81.(in Chinese)

[10]邓扬晨,陈华,马明亮,等.基于拓扑优化技术的飞机普通

框设计方法研究[J].强度与环境,2005,32(2):39-45.

Deng Yangchen,Ch en Hua,M a M ingliang,et al.Studies of aircraft frame design b as ed on topology optimization[J].

S tructure&Environment Engin eering,2005,32(2):39-45.

(in Ch ines e)

[11]W aqas S,Fan Y Q,Lu H,Z Qasim.T opology optimization

for manufacturable configuration design of aircraft tail rib [J].CADDM,2009,19(1):94-101.

[12]邓扬晨,刘晓欧,朱继宏.飞机加强框的一种结构拓扑优化

设计方法[J].飞机设计,2004,4:11-16.

Deng Yang chen,Liu Xiaoou,Zhu Jihong.An approach to s tructural topology optimiz ation about aircraft reinforced frames design[J].Aircraft Design,2004,4:11-16.(in

C hinese)

[13]E schenauer H,Becker W,Schumacher A.M ultidisciplinary

s tructural optim ization in aircraft design[R].BM BF-Bun desvorhaben DYNAFLEX.Germ an:1998.

[14]Krog L.T opology Optimization of Aircraft Wing Box Rib s

[R].AIAA-2004-4481,2004.

[15]Diaz A R,Bend soe M P.Shape optim ization of s tructures for

m ultiple loading conditions u sing a homogeniz ation method [J].S tructural Optimiz ation,1992,4:17-22.

[16]Diaz A R.Solu tion to shape an d topology eigen valu e optim-i

zation problem s us ing a homogen ization m ethod[J].In terna-

tional Journal for Num erical M eth ods in Engineer,1992,35

(7):1487-1502.

[17]Olhoff N.On CAD-integrated structu ral topology and design

optim ization com puter meth ods in applied m echanics and En-gin eering[J].Com puter M ethods in Applied M echanics and Engineering,1991,89(1-3):259-279.

[18]Sigmu nd O.M aterials w ith p rescribed constitutive param e-

ters:an invers e h omogeniz ation pr ob lem[J].International Journ al of S olids and Structu res,1994,31(17):2313-2329.

[19]S uzuki K,Kikuchi N.A homogenization method for s hape

and topology optimization[J].Com puter M ethods in Ap-plied M echan ics an d En gineering,1991,93:291-318. [20]张宪民.柔顺机构拓扑优化设计[J].机械工程学报,2003,

11:47-51.

Zhang Xian min.T opology optim ization of compliant mecha-n isms[J].Ch ines e J ournal of M echanical En gineering, 2003,11:47-51.(in Chinese)

[21]王晓明,刘震宇,郭东明.基于均匀化理论的微小型柔性结

构拓扑优化的敏度分析[J].中国机械工程,1999,10(11): 953-955.

Wang Xiaoming,Liu Zhenyu,Guo Dongmin g.Sensitivity a-n aly sis of topology optimization m icro flexible mechanical structur e by homogeniz ation m ethod[J].China M echanical Engineering,1999,10(11):953-955.(in Chin ese)

[22]罗震,郭文德,蒙永立,等.全柔性微型机构的拓扑优化设

计技术研究[J].航空学报.2005.26(5):617-626.

Luo Zh en,Guo W end e,M eng Yongli,et al.Fully compliant micro-mechanism usin g topological optimiz ation approach [J].Acta Aeronautica et As tron autica Sin ica,2005.26(5): 617-626.(in Chinese)

[23]黄杰,葛文杰,杨方.实现机翼前缘形状连续变化柔性机构

的拓扑优化[J].航空学报,2007,28(4):988-992.

Hu ang Jie,Ge W enjie,Yang Fang.Topology optimiz ation of the compliant mechanis m for shape change of airfoil leading edge[J],Acta Aeronautica et Astronau tica Sinica,2007,28

(4):988-992.(in C hinese)

[24]M aute K,Nikbay G W,Farhat C.Conceptual layout of aero

elastic w ing s tructures by top ology optimization[R].AIAA-2002-1480,2002.

[25]Santer M,Pellegrin o S.Topological optimization of com pl-i

ant adaptive w ing s tructure[J].AIAA Journal,2009,47

(3):523-534.

[26]谢先海,廖道训.基于均匀化方法的柔顺机构的设计[J].中

国机械工程,2003,14(11):1264-1266.

Xie Xianhai,Liao Daoxun.On th e design of com plian t mech-anisms usin g h om ogenization meth od[J].China M echanical Engineerin g,2003,14(11):1264-1266.(in Chin ese) [27]袁振,吴长春.复合材料扭转轴截面微结构拓扑优化设计

[J].力学学报,2003,35(1):39-42.

Yuan Zhen,Wu Changchun.Topology opti mization for micro-structures of cross sec tion of composite torsional shafts[J].Acta

8航空工程进展第2卷

M echanica Sinica,2003,35(1):39-42.(in Chinese)

[28]袁振,吴长春.复合材料周期性线弹性微结构的拓扑优化设

计[J].固体力学学报,2003,24(1):40-45.

Yuan Zhen,W u C han gchun.T op ology optimization for per-i

odic linear elastic micros tru ctures of composite materials[J].

Acta M echan ica S olida Sinica,2003,24(1):40-45.(in Ch-i

nese)

[29]任善,李书.多孔材料微结构拓扑优化设计[J].安阳工学院

学报,2008,2:19-23.

Ren S han,Li S hu.T opology optim ization design for micro-s tructures of cellular m aterials[J].J ou rnal of Anyang In st-i

tute of Technology2008,2:19-23.(in Chin ese)

[30]Tenek L H,Hagiw ara I.Optimal rectan gular plate and s hal-

low shell topologies us ing thick nes s distribu tion or homogen-i

z ation[J].C om puter M ethods in Applied M echanics and E n-gineering,1994,115(1-2):111-124.

[31]程耿东,张东旭.受应力约束的平面弹性体的拓扑优化[J].

大连理工大学学报,1995,35(1):1-9.

C heng Gengdong,Zhang Dongxu.Topological optimization

of plan e elastic continu um with s tress constraints[J].J ou r-nal of Dalian U nivers ity of T echnology,1995,35(1):1-9.

(in Ch ines e)

[32]周克民,胡云昌.利用变厚度单元进行平面连续体的拓扑优

化[J].天津城市建设学院学报,2001,7(1):33-35.

Zhou Kemin,H u Yunchang.T opology optimization of plane con tin uum usin g non-uniform thick nes s finite elemen t[J].

J ou rnal of Tianjin In stitu te of Urban Constru ction,2001,7

(1):33-35.(in Chinese)

[33]王健,程耿东.具有应力和厚度约束的平面弹性体结构拓扑

优化设计[J].机械科学与技术,2002,21(5):741-744.

W an g Jian,Chen Gen gdon g.Topology optimiz ation design of the plain elastic continu um w ith stress and thickness con-straints[J].M echanical S cien ce and T echnology,2002,21

(5):741-744.(in Chinese).

[34]M lejnek H P,R Schirrmacher.An engineer.s approach to opt-i

mal material distribution and shape finding computer method in applied mechanic and engineering[J].Computer M ethods in Ap-plied M echanics and Engi neering,1993,106(1-2):1-26. [35]Rozvan y G I N,Kir sch U,Bendsoe M https://www.360docs.net/doc/7715176331.html,yout optimiza-

tion of s tru ctures[J].Ap plied M ech anics Review s,1995,48

(2):411-419.

[36]S tolpe M,Svanberg K.An alternative in terpolation sch eme

for min imum compliance topology optimization[J].Structu r-al M ultidiscipline Optimization,2001,22(2):l16-124. [37]Lin C Y,S heu F M.Adaptive volum e con straint algorith m

for stress limit-based topology optim ization[J].Computer Aid ed Design,2009,41:685-694.

[38]Bends oe M P,Sigmun d O.Topology optimization-theory,

meth ods and applications[M].New York:Springer Pres s,

2003:5-7.

[39]Yang R J,Ch uan g C H.Optimal topology design u sing lin ear

program ming[J].Computers and Structur es,1995,52(2): 265-275.

[40]李玉霞.刮板输送机轨座的拓扑优化分析[D].石家庄:河

北工业大学机械工程学院,2006.

Li Yuxia.T he top ology optim ization analysis for scraper con-veyor track base[D].Shijiazhu ang:School of M echanical Engineering,H ebei University of T echnology,2006.(in Chinese)

[41]胡方海,王智永.履带架三维结构拓扑优化及ANS YS实现

[J].煤矿机械,2009,30(12):74-77.

Hu Fanghai,W an g Zhiyong.Structural topology optimiza-tion of craw ler frame and ANS YS r ealization[J].Coal M in e M achinery,2009,30(12):74-77.(in Ch ines e)

[42]刘秀敏,何斌,沙巍,等.空间相机主反射镜的拓扑优化设

计[J].中国光学与应用光学,2010,3(3):239-244.

Liu Xiumin,He Bin,S ha Wei,et al.Topological light-w eight design of prim ary mirror in s pace cam era[J].Ch ines e Journal of Optics an d Applied Optics,2010,3(3):239-244.

(in C hinese)

[43]夏天翔,姚卫星.机翼加强肋接头拓扑优化设计.民用飞机

设计与研究.已投稿.

Xia Tianxiang,Yao W eixing.T op ology optimization on wing-join ts of the reinforced rib[J].Civil Aircraft Design And Research.Su bmitted.(in C hinese)

[44]Sui Y K,Yang D Q.A n ew method for structural topological

optim ization based on the concept of in dependent con tin uous variables and s mooth model[J].Acta M ech anica Sinica, 1998,18(2):179-185.

[45]杨德庆,隋允康,刘正兴.多工况应力约束下的连续体结构

拓扑优化映射变换解法[J].上海交通大学学报,2000,34

(8):1061-1065.

Yang Deqin g,Sui Yunkang,Liu Zhengxing.M appin g trans-formation method for topology optimization of con tin uum structur es under m ultiple loadin g cases and stres s cons traint [J].Journal of Shanghai Jiaotong University.2000,34(8): 1061-1065.(in Chinese)

[46]隋允康,张学胜,龙连春,等.位移约束集成化处理的连续

体结构拓扑优化[J].固体力学学报,2006,27(1):102-107.

Sui Yunkang,Zhang Xuesheng,Long Li anchun,et al.Topology optimization of continuum structures w i th the integrated approach of displacement constrai nts[J].Acta Mechanica Solida Sinica, 2006,27(1):102-107.(in Chinese)

[47]隋允康,杨德庆,王备.多工况应力和位移约束下连续体结

构拓扑优化[J].力学学报,2000,32(2):171-179.

Sui Yunkang,Yang Deqing,Wang Bei.T opological optim-i z ation of continuum structu re w ith stress an d displacem ent constraints under multiple loading cases[J].Acta M echanica Sinica,2000,32(2):171-179.(in Chinese)

[48]隋允康,彭细荣,叶红玲.应力约束全局化处理的连续体结

构ICM拓扑优化方法[J].工程力学,2006,23(7):1-7.

Sui Yunkang,Peng Xirong,Ye H onglin g.T opology optim-i

9

第1期夏天翔等:连续体结构拓扑优化方法评述

zation of continuum structu re w ith glob aliz ation of stress con-straints by ICM meth od[J].Engin eering M echanics,2006, 23(7):1-7.(in Chinese)

[49]隋允康,边炳传,叶红玲.连续体结构屈曲约束的ICM方法

拓扑优化[J].计算力学学报,2008,25(3):315-351.

Sui Yunkang,Bian Bingchuan,Ye Hongling.Topological opt-i

mization of the continuum structure s ubjected to the buckling constraints with ICM method[J].Chinese Journal of Computa-tional M echanics,2008,25(3):315-351.(in Chinese)

[50]彭细荣,隋允康.有频率禁区的连续体结构拓扑优化[J].固

体力学学报,2007,28(2):145-150.

Peng Xirong,Su i Yunkang.T opological optimization of the continuum structu re w ith non-frequen cy-b and constraints [J].Acta M echanica Solida Sinica,2007,28(2):145-150.

(in C hinese)

[51]叶红玲,隋允康.基于ICM方法三维连续体结构拓扑优化

[J].固体力学学报,2006,27(4):387-393.

Ye H ongling,Sui Yunk ang,ICM-based topological optimiza-tion of three-dimen sional con tin uum structure[J].Acta M e-chan ica Solida Sinica,2006,27(4):387-393.(in Chinese) [52]隋允康,于新.平面膜结构拓扑优化的有无复合体方法[J].

力学学报,2001,33(3):357-364.

Sui Yunkang,Yu Xin,T he exist-null combination method for the topological optimization of plane membrane structure[J].

Acta M echanica Si nica,2001,33(3):357-364.(in Chinese) [53]于新,隋允康.平面连续体多工况应力约束下结构拓扑优化

的有无复合体方法[J].计算力学学报,2001,18(2): 200-204.

Yu Xin,S ui Yunkang.Ex ist-n ull com bination meth od for the topological optimization of plan e continuum stru cture under m ult-i loading cas e w ith s tr ess constraint[J].C hinese J ournal of Computational M echan ics,2001,18(2):200-204.(in Ch ines e)

[54]隋允康,于新,叶宝瑞.应力和位移约束下桁架拓扑优化的

有无复合体方法[J].固体力学学报,2004,25(3):355-359.

S ui Yunk ang,Yu Xin,Ye Baor ui.Exist-null combination m ethod for the topological optimization of truss w ith stres s and displacement cons traint[J]Acta M echanica Solida Sin-i

ca,2004,25(3):355-359.(in Chin ese)

[55]隋允康,彭细荣.结构拓扑优化ICM方法的改善[J].力学

学报,2005,37(2):190-198.

S ui Yun kang,Peng Xirong.Th e im provement for the ICM m ethod of structural topology optimiz ation[J].Acta M e-chan ica Sin ica2005,37(2):190-198.(in Chin ese)

[56]边炳传,隋允康,李善坡.拓扑优化在变压器吊耳设计中的

应用[J].力学与实践,2008,35(5):75-78.

Bian Bingchuan,Sui Yunkang,Li Shanpo.Application of topolo-gical opti mization to the trunnion design for transformers[J].

M echani cs in Engineering,2008,35(5):75-78.(in Chi nese) [57]李芳,张建润,卢熹.基于ICM法的零件拓扑优化方法应用

研究[J].电子机械工程,2008,24(1):40-42.

Li Fang,Zhang J ian run,Lu Xi.Application s tu dy of ICM-based topology optimiz ation[J].Electro-M echan ical Eng-i n eering,2008,24(1):40-42.(in C hinese)

[58]Xie Y M,Steven G P.A s imple evolutionary procedure for

structural optimiz ation[J].C om puters and Structures, 1993,49:885-896.

[59]Xie Y M,Steven G P.Optimal design of m ultiple load cas e

structur es usin g an evolutionary procedure[J].En gineering Computations,1994,11:295-302.

[60]Zh ou M,Roz vany G I N.On the validity of ESO type m eth-

ods in topology optimization[J].Structu ral M ultidisciplin e Optim ization,2001,21:80-83.

[61]Chu D N,Xie Y M,A Hear,et al.Evolutionary structure opt-i

mization for problems w ith s tiffness constrains[J].Fini te Ele-ments i n Analysis and Design,1996,21:239-251.

[62]Xie Y M,Steven G P.E volution ary structur e optimiz ation

for dynamic problem s[J].Computers and S tructures.1996, 58(6):1067-1073.

[63]M anickarajah D,Xie Y M,Steven G P.An evolu tionary

meth od for optimization of plate b uckling resistance[J].F-i

nite Elements in Analysis and Design,1998,29:205-230.

[64]M anickarajah D,Xie Y M,Steven G P.Optimization of col-

umns an d frames again st buckling[J].Computers and S truc-tu res,2000,75:45-54.

[65]Rong J H,Xie Y M,Yang X Y.An improved meth od for

evolutionary s tructural optimization against buckling[J].

Computer s and Structur es,2001,79:253-263.

[66]Lian g Q Q,Xie Y M,Steven G P.Optimal selection of to-

pologies for the minimu m-w eight des ign of continuum s truc-tu res w ith stres s constraints[J].Proceedin gs of the Institu-tion of M echanical Engin eer s,Part C:Journ al of M echanical Engineering Science,1999,213(8):755-762.

[67]Reynolds D,M cconn achie J,Bettess P,et al.Reverse adap-

tivity a n ew evolutionary tool for structu ral optimization[J].

International J ou rnal for Nu merical M ethods in En gineering.

1999,45:529-552.

[68]Querin O M,Young https://www.360docs.net/doc/7715176331.html,putational efficiency and valida-

tion of b-i directional evolutionary structur e optimization[J].

Computer M ethod in Applied M echanics Engineering,2000, 189:559-573.

[69]Yang X Y,Xie Y M,S teven G P,et al.T opology optimiza-

tion for frequ ency usin g an evolutionary method[J].Journal of Structu ral Engin eering,1999,12:1432-1438.

[70]Liu X,Yia W J,Lib Q S,et al.Genetic evolutionary s truc-

tu ral optimization[J].Journal of Constructional Steel Re-search,2008,64:305-311.

[71]Zhu J H,Zhang W H,Qiu K P.B-i directional evolu tionary

top ology optimization u sing element replaceable method[J].

Computational M echanics,2006,40(1):97-109.

[72]郭中泽,陈裕泽,张卫红,等.基于单元材料属性更改的结

构渐进拓扑优化方法[J].机械科学与技术,2006,25(8):

10航空工程进展第2卷

928-931.

Guo Zh on gze,Chen Yuz e,Zhang W eih ong,et al.An evolu-tionary s tructural optimization(ES O)meth od b as ed on ele-m ent's properties changing method[J].M ech anical Science An d T echnology,2006,25(8):928-931.(in C hinese) [73]张爱林,杨海军,姚力.基于应变能的预应力平面实体结构

拓扑优化设计[J].工程力学,2008,25(10):110-116.

Zhang Ailin,Yang Haijun,Yaoli.Topology optimiz ation de-sign of pr estressed plane en tity structures b as ed on strain en-ergy[J].E ngineering M echan ics,2008,25(10):110-116.

(in Ch ines e)

[74]Yang X Y,Xie Y M,Liu J S,et al.Peri m eter control in th e

bidirectional evolutionary optimization meth od[J].Stru ctural and M ultidisciplinary Optimization,2003,24(6):430-440.

[75]Cervera E,Trevelyan J.Evolutionary structural optimization

based on boundary representation of NURBS.Part II:3D algo-

ri thms[J].Computers and Structures,2005,83:1917-1929.

[76]H uang X,Xie Y M.B-i directional evolu tion ary topology op-

timiz ation of continuum structures w ith on e or mu ltiple mate-rials[J].Computational M echanics,2009,43:393-401. [77]王伟.机翼结构拓扑优化的一种新渐进结构优化方法[J].

飞机设计,2007,27(4):19-20.

W an g W ei.A n ew evolu tionary s tructural optimization meth-od and its application in the w ing structure topology optimiza-tion[J].Aircraft Design,2007,27(4):19-20.(in Chin ese) [78]王伟,杨伟,赵美英.大展弦比飞翼结构拓扑、形状与尺寸综

合优化设计[J].机械强度,2008,30(4):596-600.

W ang W ei,Yang Wei,Zhao M eiying.Integrate topology/

shape/size optim ization into h igh aspect-ratio w ing des ign [J].Journal of M echan ical Stren gth,2008,30(4):596-600.(in Chinese)

[79]倪晓宇,易红,汤文成,等.机床床身结构的有限元分析与

优化[J].设计与研究,2005,2:47-50.

Ni Xiaoyu,Yi H ong,Tang W en cheng,et al.Finite element analysis and optim ization of machine tool b ody[J].Des ign and Research,2005,2:47-50.(in C hinese)

[80]T an skanen P.T he evolutionary structu ral optimization

method:theoretical aspects[J].Compu ter M eth od in Ap-

plied M ech anics Engineering,2002,191:5485-5498. [81]Osher S,Seth ian J A.Fronts propagating w ith curvatu re de-

pendent s peed:algorithms based on ham ilton-jacobi for mula-tions[J].J ou rnal of Computational Physics,1988,79:

12-49.

[82]S ethian J A,Wiegman n A.Structural b oundary design via

level s et and immersed interface m ethods[J].Journal of

C om putational Ph ysics,2000,163:489-528.

[83]W an g M Y,Wang X M,Guo D M.A level set method for

s tructural topology optimiz ation[J].Computer M ethods in

A pplied M echanics and Engineering.2003,192:227-246.

[84]Allaire G,Jouve F,T oader A M.Structu ral optimization

u sing sensitivity analysis and a leve-l set method[J].Journal of Com putational Physics,2004,194:363-393.

[85]陈淑玉.连续体拓扑优化设计[D].济南:山东大学机械工

程学院,2006.

Chen S huyu.Topology optimiz ation of continuum structu re

[D].Jinan:Sch ool of M echan ical Engineering,S han Dong

University,2006.

[86]荣见华,梁清泉,杨端生.基于结构拓扑随机变异的水平集

优化方法[J].力学学报,2007,39(6):796-803.

Rong J ian hua,Liang Qingquan,Yang duansh eng.A level set method for structural topology optimization based on to-pology random mutations[J].Chinese Journal of T heoretical and Applied M echan ics,2007,39(6):796-803.(in Ch ines e) [87]罗俊召,陈立平,罗震.基于半隐式格式的水平集法连续体

结构形状和拓扑优化方法[J].固体力学学报,2008,29(2): 175-180.

Luo Junzhao,Chen Liping,Luo Zhen.A level set method for to-pology opti mization usi ng AOS scheme[J].Chi nese Journal of Solid M echanics,2008,29(2):175-180.(in Chinese)

[88]M ei Y L,W an g X M.A level s et meth od for structu ral to-

pology optimization and its applications[J].Advances in En-gin eering Softw are,2004,35:415-441.

[89]王海滨.Level set算法及其应用[D].西安:西北工业大学

理学院,2006.

Wang H aibin.Level s et algorithm an d its application[D].

Xi.an:School of S cien ce,North w estern Polytechnical U n-i

versity,2006.(in Chin ese)

[90]荣见华.一种改进的结构拓扑优化水平集方法[J].力学学

报,2007,39(2):253-260.

Rong Jianhua.An improved level set m ethod for structural top ology optimization[J].Chinese Journal of T heoretical and Applied M echan ics,2007,39(2):253-260.(in Ch ines e) [91]W ang S Y,Lim K M,Wang M Y.An extended level set

m ethod for sh ape and topology optimiz ation[J].Journ al of Compu tational Phys ics,2007,221:395-421.

[92]欧阳高飞,张宪民.结构拓扑优化中水平集方法的快速初始

化研究[J].应用基础与工程科学学报,2008,16(1): 136-143.

Ouyang Gaofei,Zhan g Xianmin.Research on fast initializa-tion for the level s et method in stru ctural topology optimiza-tion[J].Journal of Basic Science and Engineering,2008,16

(1):136-143.(in Chin ese)

[93]胥志刚,林忠钦,来新民,等.面向车身结构轻量化设计的

水平集拓扑优化[J].上海交通大学学报,2007,41(9): 1393-1401.

Xu Zhigang,Lin Zhongqin,Lai Xinmin g,et al.A level set top ology optimization m ethod for auto-b ody structu re light-w eight design[J].Journ al of Shanghai Jiaotong University, 2007,41(9):1393-1401.(in Chinese)(下转第55页)

11

第1期夏天翔等:连续体结构拓扑优化方法评述

10-14.

Zh u Aipin g,M a Ling.Area control am munition[J].Winged M issiles J ou rnal,2009(8):10-14.(in Chinese)

[14]樊燕君.基于DSP的巡飞弹弹载控制系统设汁研究[D].南

京理工大学机械工程学院.2007.

Fan Yanju n.Control system of loitering mu nition s depen ds on DS P[D].School of M ech anical Engineering,Nanjing U-niversity of Science and T echn ology,2007.(in Chinese) [15]吕强,叶正寅,李栋.充气结构机翼的设计和试验研究[J].

飞行力学,2007,25(4):77-80.

L v Qiang,Ye Zhengyin,Li Dong.Design and capability a-nalysis of an aircraft w ith inflatab le w ing[J].Flight Dynam-ics,2007,25(4):77-80.(in C hines e)

[16]宋淼,苏德利,薛宏伟,等.悬浮弹空中姿态稳定技术探讨

[J].兵工自动化,2009,28(11):33-35.

Song M iao,Su Deli,Xue Hongw ei,et al.Study of air atti tude stability technique of suspended ammuni tion[J].Ordnance Indus-try Automation,2009,28(11):33-35.(in Chinese)

[17]程养民.微重力火箭设计中减旋方案等问题[J].固体火箭,

1999,22(1):1-5.

C heng Yan gmin.On the problems of gyration reduction etc

in microgravity rocket design[J].J ou rnal of S olid Rock et T echnology,1999,22(1):1-5.(in Chinese)[18]乐挺,王立新,艾俊强.Z型翼变体飞机的纵向多体动力学

特性[J].航空学报,2010,31(4):679-686.

Le Tin g,Wang Li Xin,Ai Junqiang.Longitudin al multibody dynamic characteristics of Z-w ing morphing aircraft[J].Acta Aeronautica E t Astronautica S inica,2010,31(4):679-686.

(in C hines e)

[19]郭美芳,范宁军.多模式战斗部与起爆技术研究[J].探测与

控制学报,2005,27(1):31-34.

Guo M eifang,Fan Ningjun.The study on a mutibody w arhead and i nitiation technology[J].Journal of Detecti on&Control, 2005,27(1):31-34.(in Chinese)

[20]胡照.炮弹新型战斗部技术发展概况[J].科学技术,2010,

5:18-19.

Hu zhao.T he development of new w arh ead in s hell[J].

Science&T echnology,2010,5:18-19.(in Chinese)

作者简介:

杨延年(1985-),男,硕士研究生。主要研究方向:飞行器设计。

张元明(1964-),女,研究员。主要研究方向:无人机设计和飞行试验研究。

(编辑:赵毓梅)

(上接第11页)

[94]Yamada T,Izui K,Ni shiw aki S,et al.A topology optimization

method based on the level set method incorporating a fictitious in-terface energy[J].Computer Methods in Applied M echanics and Engineering.doi:10.1016/j.cma.2010.05.013.

[95]黄杰.基于变密度法的连续体拓扑优化技术及其应用研究

[D].西安:西北工业大学机电学院,2007.

H uang Jie.Research of theory and applic ation about topology op-

timization of continuum structure based on alterable density method[D].Xi.an:School of M echanical Engineering,North-w estern Polytechnical University,2007.(i n Chinese)

[96]杨志勇.基于变密度法的连续体结构拓扑优化研究[D].西

安:西安电子科技大学机械电子工程学院,2009.

Yang Zhiyong.Res earch on topology optim ization of con tin u-um stru cture based on alterable density m ethod[D].Xi.an: S chool of Electron ic Engineering,Xidian un iversity,2009.

(in Ch ines e)

[97]隋允康,叶红玲,彭细荣,等.连续体结构拓扑优化应力约束

凝聚化的ICM方法[J].力学学报,2007,39(4):554-563.

S ui Yunkang,Ye H on gling,Pen g Xir on g,et al.T he ICM m ethod for continuum structu ral top ology optimization w ith con dens ation of stress cons traints[J].Chin es e Journ al of T heoretical and Applied M echan ics,2007,39(4):554-563.

(in Chin es e)

[98]Dias A R,S igmund O.Checkerboard patterns in layout opt-i

m ization[J].Structu ral Optimization,1995,10(1):40-45.[99]Sigmun d O.Numerical ins tabilities in topology optimiz ation:

a sur vey on procedures dealing w ith checkerboard,m es h de-

pendencies,local minim a[J].Structural Optimization,1998, 16(1):68-75.

[100]Cho S.Kw ak J.T opology design optimization of geom etr-i cally non-lin ear stru ctures usin g mesh free meth od[J].

Com puter M ethods in Applied M echanics and En gineering,

2006,195:5909-5925.

[101]Zh eng J,Long S Y,Li G Y.Th e topology optimiz ation de-sign for continuum structu res based on the elem ent free

Galerkin meth od[J].Engineering Analysis w ith Boundary

Elements,2010,34:666-672.

[102]宋康祖,陆明万,张雄.固体力学中的无网格方法[J].力学进展,2000,30(1):55-65.

Song Kangzu,Lu M ingw an,Zhang Xiong.M eshles s m eth-

od for s olid m ech anics[J].Advances in M echan ics,2000,

30(1):55-65.(in Chin ese)

作者简介:

夏天翔(1987-),男,博士研究生。主要研究方向:结构优化设计、结构疲劳理论。

姚卫星(1957-),男,博导,教授。主要研究方向:结构疲劳理论、结构可靠性、结构优化设计。

(编辑:赵毓梅)

55

第1期杨延年等:一种新型无人飞行器的发展现状与关键技术