基于子结构的宏微结构协同优化方法

第３３卷第２３期
中国机械工程
V o l ．３３㊀N o ．２３２０２２年１２月
C H I N A M E C HA N I C A LE N G I N E E R I N G
p p
．２８５１Ｇ２８５８基于子结构的宏微结构协同优化方法
吴紫俊１㊀肖人彬２
１．武汉纺织大学数字化纺织装备湖北省重点实验室,武汉,４３００２０
２．华中科技大学人工智能与自动化学院,武汉,４３００７４
摘要:利用子结构自由度的缩减实现结构宏微观的一体化设计,需要预先建立不同子结构构型的样本.宏观结构优化时只能在预定的微观构型中插值,限制了结构宏微观的优化空间.结合子结构自由度缩减与反求理论,针对周期性结构的设计提出了结构的宏微观协同设计方法.利用固体各向同性材料惩罚(S I M P )方法和子结构方法,构建了基于子结构的周期性结构宏微观协同设计框架;分析了宏观结构与微观周期性单胞的材料体积变化关联关系,并给出了两个尺度上的惩罚因子㊁变量过滤半径等优化参数的匹配.在此基础上,结合传统的结构柔度最小优化模型,计算了设计变量在宏观结构㊁微观周期性结构的灵敏度值,实现了宏微结构的协同设计.通过悬臂梁结构的设计,验证了所提方法的有效性
.
关键词:协同设计;周期性结构;宏微观优化设计;子结构中图分类号:T P ３９１
D O I :１０．３９６９/j ．i s s n ．１００４ １３２X．２０２２．２３．００９开放科学(资源服务)标识码(O S I D )
:AS u b s t r u c t u r e Ｇb a s e dC o Ｇo p
t i m i z a t i o n M e t h o d f o rM a c r o Ｇm i c r o S t r u c t u r e s WUZ i j
u n １㊀X I A O R e n b i n ２
１．H u b e iK e y L a b o r a t o r y o fD i g i t a lT e x t i l eE q u i p m e n t ,W u h a nT e x t i l eU n i v e r s i t y ,W u h a n ,４３００２０２．S c h o o l o fA r t i f i c i a l I n t e l l i g e n c e a n dA u t o m a t i o n ,H u a z h o n g U n i v e r s i t y o
f S c i e n c e a n d T e c h n o l o
g y
,W u h a n ,４３００７４A b s t r a c t :T h e i n t e g r a t e d s t r u c t u r a lm a c r o Ｇm i c r o s c o p i c d e s i g nw a s a c h i e v e d b y u s i n g
t h e r e d u c t i o n o f s u b s t r u c t u r a l d e g r e e s o f f r e e d o m ,w h i c h r e q u i r e d t h e p r e Ｇe s t a b l i s h m e n t o f s a m p l e s o f d i f f e r e n t s u b Ｇs t r u c t u r a l c o n f i g u r a t i o n s ．T h em a c r o s c o p i c s t r u c t u r e o p t i m i z a t i o nm i g h t o n l y b e i n t e r p
o l a t e d i n p r e d e Ｇf i n e dm i c r o s c o p i c c o n f i g u r a t i o n s ,w h i c h r e s t r i c t e d t h e o p t i m i z a t i o ns p
a c eo f s t r u c t u r a lm a c r oa n dm i Ｇc r o s c o p i c ．A m a c r o Ｇm i c r o c o l l a
b o r a t i v em e t h o dw a s p r o p o s e d f o r t h ed e s i g no f p e r i o d i
c s t r u c t u r e sb y c o m b i n i n g s u b s t r u c t u r e
d
e g r e e o
f f r e e d o mr e d u c t i o n a n d i n v e r s e t h e o r y ．U s i n
g t
h e s o l
i d i s o t r o p
i cm a Ｇt e r i a lw i t h p e n a l i z a t i o n (S I M P )m e t h o d a n d t h e s u b s t r u c t u r em e t h o d ,a s u b s t r u c t u r e Ｇb a s e dm a c r o Ｇm i Ｇ
c r o c o l l a b o r a t i v e f r a m e w o r k f o r p e r i o
d i c s t r u c t u r
e sw a s c o n s t r u c t e d ．T h e c o r r e l a t i o nb e t w e e n t h em a Ｇ
t e r i a l v o l u m e v a r i a t i o no fm a c r o s t r u c t u r e s a n d m i c r o s c o p i c p e r i o d i c c e l lw a s a n a l y
z e d ,a n d t h em a t c Ｇh i n g o fo p t i m i z a t i o n p a r a m e t e r s w a s g i v e ns u c ha s p e n a l t y f
a c t o ra n dv a r i a
b l ef i l t e rr a d i u sa tt w o s
c a l e s ．T h e n ,t h e s e n s i t i v i t y v a l u e so f t h e
d
e s i g nv a r i a b l e s i nt h em a c r o s c o p
i cs t r u c t u r e sa n d m i c r o Ｇs c o p i c p e r i o d i c s t r u c t u r e sw e r e c a l c u l a t e d i n c o m b i n a t i o nw i t h t h e t r a d i t i o n a l s t r u c t u r a l f l e x i b i l i t y m
i n Ｇi m u mo p t i m i z a t i o nm o d e l ,a n dt h ec o Ｇd e s i g
no fm a c r oa n d m i c r os t r u c t u r e sw a s r e a l i z e d ．T h ee f f e c Ｇt i v e n e s s o f t h e p r o p o s e dm e t h o dw a s v e r i f i e db y t h e d e s i g
no f t h e c a n t i l e v e r b e a ms t r u c t u r e s ．K e y w
o r d s :c o l l a b o r a t i v e d e s i g n ;p e r i o d i c s t r u c t u r e ;m a c r o Ｇm i c r o s c o p i c o p t i m i z a t i o nd e s i g n ;s u b Ｇs t r u c t u r e
收稿日期:２０２１１２０６
基金项目:国家自然科学基金(５１８７５２２０);湖北省教育厅科学技术研究项目(D ２０２２１７０１
)０㊀引言
基于子结构方法的宏微结构一体化优化方法中,利用子结构自由度缩减实现了结构宏微观的一体化设计,但需要预先建立不同子结构构型的
样本[１]
,宏观结构优化时只能在预定的微观构型中进行插值,即限制了结构宏微观的优化空间[
２]
.宏微结构的协同优化设计是从宏微多个尺度寻找材料的最佳分布位置,组合成不同构型的微观构型和宏观结构,使得材料潜力充分发挥[
３
].由于结构设计引入了微观尺度,因而扩大了结构
设计空间[４
],赋予了结构更多的功能,如最小频率响应的结构[５Ｇ６
]㊁具有良好可制造性的宏观结构等[
７
].目前,周期性子结构的优化方法主要有两
类[８Ｇ９]:基于均匀化方法的周期性结构设计[１０Ｇ１１
]与
基于宏观结构周期性约束条件的设计方法.在基
１５８２ Copyright ©博看网. All Rights Reserved.
于均匀化方法[１２Ｇ１３]的周期性结构设计方法中,周期性位移边界条件施加在单胞上[１４],并假设结构的应力和应变在宏观结构上是周期性变化的[１５Ｇ１６].而基于周期性约束的方法中,假设宏观结构中任意两个周期性结构中的单胞具有相同的单元密度和相等的灵敏度值.这两类方法相比,尽管都能获得相近的优化构型,但基于周期性约束条件的设计方法的计算效率要低于均匀化方法.
作为结构拓扑优化的热点分支之一,周期性结构设计得到了许多关注.L I等[１７Ｇ１８]提出基于维度缩减的代理模型来设计周期性结构,降低了结构频率带隙计算的计算资源消耗.Z H E N G 等[１９]基于不确定性负载的高斯分布提出了周期结构鲁棒拓扑优化框架,研究了周期结构设计空间小的结构对不确定性载荷的敏感度.H E 等[２０]针对不规则有限元网格提出了一种周期性结构优化方法,用非结构设计点来清晰地描述周期结构的拓扑.R I V A等[２１]提出了一种基于波在细胞界面反射和传输的纯解析方法,可判断所提出的解决方案的最优性.C H E N等[２２]提出了一种用于结构和蜂窝材料拓扑并行设计的移动等值面阈值(m o v i n g i s oＧs u r f a c e t h r e s h o l d,M I S T)方法,利用均匀化理论推导出了宏观结构和微观周期性单元两个维度的物理响应函数.C H E N G 等[２３]研究了周期性均质多孔材料板和均匀加筋实心板的最大面外屈曲载荷设计,并基于渐进均匀化方法和C A I方法建立了结构宏观和微观设计变量与灵敏度分析的关系.D A H L B E R G 等[２４]引入周期性结构单元以获得具有更高对称性的结构,提出利用高对称性设计导频结构,大大减小了低频电波传输时的漏波现象.
赵清海等[２５]基于周期性约束提出了多材料结构稳态热传导拓扑优化设计方法,并结合O r d e r e dＧR AM P材料插值方法获得了周期性热传导拓扑构型.焦洪宇等[２６]利用固体各向同性材料惩罚(s o l i d i s o t r o p i cm a t e r i a lw i t h p e n a l i z aＧt i o n,S I M P)方法,通过增加周期性约束条件,提出了各向同性微结构材料的周期性拓扑优化方法,并在其随后的研究中[２７],对循环对称结构㊁梯形结构等进行了周期性布局优化研究.杜义贤等[２８]利用能量均匀化方法建立了基于宏观力学性能的微观点阵结构的优化模型,以周期性单胞为研究对象,将材料用量和力学方程等作为约束条件,获得了边界清晰的周期性点阵结构;并在其后的研究中[１０],结合负泊松比表征吸能特性,提
出了兼具吸能和承载特性的周期性宏细观结构多尺度优化方法.
在周期性结构的多尺度设计中,为在设计域中设计微结构构型的同时兼顾微结构在宏观结构中的最优分布[９,２９],本文提出了基于子结构法的周期性结构的宏微观设计方法,并利用子结构法中的自由度凝聚来实现周期性单胞的矩阵规模缩减,提高其计算效率的同时,也可利用子结构的反求来完成周期性单胞结构与宏观构型的协同设计.
与此同时,文献[１]和文献[２]中提出的基于自定义微结构构型的宏微结构一体化设计方法A R S P(a p p r o x i m a t i o n o fr e d u c e d s u b s t r u c t u r e w i t h p e n a l i z a t i o n),由于限定了微结构构型的设计范围,其宏观结构与自定义微结构构型之间存在兼容匹配问题.本文为突破预定义微结构构型的设计限制,利用子结构的自由度缩减与反求来实现两个尺度上结构的协同设计,有助于实现大规模个性化设计[３０].
１㊀子结构与优化插值模型
１．１㊀子结构凝聚与反求
在优化子结构中,单个子结构的本构方程可表示为[３１]
K b b K b i
K i b K i i
é
ë
êêù
û
úú
a b
a i
é
ë
êêù
û
úú＝
P b
P i
é
ë
êêù
û
úú(１)其中,下标b㊁i分别表示子结构边界节点和内部节点;a b和a i为对应节点的位移向量;子矩阵K b b为边界节点的刚度矩阵;K i i为子结构内部节点的刚度矩阵;K T b i＝K i b为子结构中边界节点和内部节点的耦合子刚度矩阵;P b㊁P i分别为子结构边界节点和内部节点的载荷列阵.
由式(１)的第二式可得子结构内部节点位移向量:
a i＝K－１i i(P i－K i
b a b)(２)㊀㊀将式(２)代入式(１)的第一式,可得方程:
(K b b－K b i K－１i i K i b)a b＝P b－K b i K－１i i P i(３)㊀㊀式(３)可简单表示为
K∗b b a b＝P∗b(４)
K∗b b＝K b b－K b i K－１i i K i b(５)此时K∗b b即为子结构凝聚后的刚度矩阵.
本文中,将利用式(５)进行子结构凝聚组建宏观本构方程的过程称为子结构凝聚,通过凝聚后的本构方程进行结构设计称为宏观结构设计;利用式(２)计算子结构内部自由度的位移的过程称为子结构反求,通过子结构反求获得对应的周期性单胞的过程称为微观结构设计[３２].
２５８２
中国机械工程第３３卷第２３期２０２２年１２月上半月Copyright©博看网. All Rights Reserved.
１．２㊀优化插值模型的构建
在周期性宏观结构设计中,假设宏观结构由若干个子结构组成,如图１所示
.
图１㊀周期性子结构与宏观结构
F i g
．１㊀T h e p e r i o d i c s u b s t r u c t u r e a n dm a c r o s t r u c t u r e 根据子结构可知,宏观结构设计域Ω的计算可直接利用凝聚后的刚度矩阵表示:
K s U s ＝F s
(６
)K s ＝ʏ
Ω
K ∗
b b d
Ω其中,K s 为子结构凝聚后所形成的整体刚度矩
阵;U s ㊁F s 为凝聚后子结构交界面上节点的位移向量及载荷列阵.
针对宏观结构设计,设计变量设置为子结构所含的材料体积,由于宏观结构由周期性结构组成,故只需要判断子结构材料的含量.此时宏观子结构材料含量设计变量的取值为(０,１],对应的刚度矩阵可表示为
K ＝ρp
K ∗b b
(７
)其中,ρ为设计变量;
p 为ρ的惩罚因子.子结构密度决定了周期性子结构构型的材料含量.微观子结构构型设计中,设计密度为子结构中各个单元的材料含量,利用式(２)可获得子结构的内部节点值,通过优化变量更新即可获得微观结构的构型.假设子结构内部单元的刚度矩阵为K ０,它与子结构优化密度变量之间的关系可表示为
K b ＝(ρ
b )q
K ０(８
)其中,ρ
b 为该子结构内部单元的材料相对含量;q 为ρb 的惩罚因子.
２㊀基于周期性微观单胞的宏微观结构优化
２．１㊀宏微观优化模型定义
以结构柔度最小为优化目标,实现宏微观结构的协同优化:
f i n d X (ρ)
m i n ㊀c (ρ)＝U T
K U
s ．t ．K U ＝F ㊀㊀㊀V ɤυ㊀㊀㊀０＜ρm i n ɤρi ɤ１
㊀㊀i ＝１,２, ,N üþýï
ïï
ï
ï
ï(９)式中,X 为宏观结构中所有子结构的密度集合;c 为宏观结构柔度设计目标函数;V 为宏观结构的材料含量总量;
K ㊁U ㊁F 分别为整体刚度矩阵㊁位移场和边界约束;υ为体积约束;ρi 为各子结构设计变量;ρm i n 是为避免奇异问题设定的最小密度,ρm
i n ＝０．００１.基于有限元子结构方法,优化问题(式(９
))在宏观层面的优化可表示为
f i n d ㊀ρ
m i n ㊀c (ρ)＝U T K U ＝ðN
i ＝１
U T b b ,i K ∗
b b ,i U b b ,i s ．t ．K U ＝F ㊀㊀㊀V ɤυ
㊀㊀㊀０＜ρm i n ɤρi ɤ１
㊀㊀i ＝１,２, ,N üþýï
ïïïï
ïïï(
１０)式中,ρ为每个子结构的材料相对含量;
U b b ,i ㊁K ∗
b b ,i 分别为第i 个子结构凝聚后的边界自由度对应的位移场和刚度矩阵;N 为宏观设计域子结构总数.
微观结构是构型相同的周期性结构,优化模型可表示为
f i n d ㊀ρ
b m i n ㊀
c b (ρb )＝ðm
j ＝１
U T
j ,i K ０U j ,i s ．t ．K ∗
b b ,i U b b ,i ＝F b b ,i ㊀㊀㊀V b ɤυN /ðN
i ＝１
s g n (ρ
i ＞ρm i n )㊀㊀㊀ρj ɪ[ρm
i n ,１]㊀㊀j ＝１,２, ,m üþý
ï
ïïïïïïïï(１１
)s g
n (x )＝１㊀ρi ＞ρm i n ０㊀ρ
i ɤρm i n {
式中,c b 为子结构柔度优化目标;K ０㊁U j ,i 分别为第i 个子结构中单元的单元刚度矩阵和位移场;F b b ,i 为当前子结构凝聚后的边界自由度对应的位移场;m 为当前子结构中单元总数;ρb 为微结构构型的设计变量;V b 为当前子结构的体积约束,此时该约束是变动的,与子结构在宏观结构的分布数量相关.
通过上述宏微观的优化模型可知,结合S I M P 方法计算,
宏观结构的拓扑构型由宏观子结构的材料含量确定.而在微观结构设计中,微观周期性结构设计的材料约束不是固定值,而是根据宏观结构的子结构的材料密度分布进行更新.当宏观结构中周期性子结构分布在整个设计域时,微观结构的材料含量约束与宏观材料约束相同;当周期性结构分布在宏观设计域的子域时,微观结构的材料含量约束将增大.２．２㊀宏微观两尺度的灵敏度分析宏观结构由构型相同的周期性结构组成,假
设对式(１０)的目标函数计算密度导数:∂c (ρ)∂ρi ＝－∂(U T ρp
K U )∂ρ
＝－p ρp －
１i U T b b ,i K ∗b b ,i U b b ,i (１２
)㊀㊀微观结构的灵敏度可表示为
３５８２ 基于子结构的宏微结构协同优化方法
吴紫俊㊀肖人彬Copyright ©博看网. All Rights Reserved.
∂c b (ρ)
∂ρ
j ＝－∂ðm
j ＝１
U T
j ,i ρq b K ０U j ,i ∂ρ
j ＝
－q (ρ
j )q
－１
ðm
j ＝１
U
T
j ,
i K ０U j ,i
(１３
)㊀㊀微观上灵敏度与宏观子结构相关,
结构优化过程中,在微观结构设计变量增加宏观的惩罚因子项可加快其收敛速度.
在基于子结构的优化方法中,两个维度上的
变量更新策略均采用传统的优化准则法(o p
t i Ｇm a l i t y c
r i t e r i am e t h o d ,O C )[３３]
:ρ
n e w
i ＝m
a x (ρm i n ,ρi －ε)㊀㊀㊀㊀㊀㊀㊀㊀㊀ρi D ηi ɤm a x (ρm i n ,ρ
i －ε)ρi D ηi ㊀㊀m a x (ρm i n ,ρi －ε)ɤρi D ηi ɤm i n (１,ρi ＋ε)m i n (１,ρi ＋ε)㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀m i n (１,ρ
i ＋ε)ɤρi D ηi ìîíïïïïïïïï(１４
)㊀㊀在上述更新策略中,ε是变量迭代步长.η＝０．５,而优化条件D i 可由子结构的目标函数和体积导数定义:
D i ＝－∂c ∂ρi (κ∂V ∂ρi )－１
(１５
)式中,κ为拉格朗日因子;∂V /∂ρi ＝v i 为每个子结构的体积约束.
２．３㊀周期性结构的宏微观拓扑优化算法
在该宏微观结构协同优化中,通过子结构自由度凝聚和反求计算,可以获得设计域上所有子结构边界节点和内部节点的位移场.通过式(１２
)和式(１３)计算对应的宏观结构和微观结构的灵敏度.由于微观周期性结构的构型会影响宏观子结构凝聚的刚度矩阵,因此在进行变量更新时,应先更新周期性结构的设计变量,然后更新宏观结构设计变量.对应的优化设计流程如图２所示,关键算法如下所示.
D e g i n M a c r o S t r u c t u r e ()ʊ宏观结构设计㊀w h i l e (判断是否收敛)㊀{
㊀㊀s u b s t r c u t u r e M a c r o (网格划分参数);ʊ根据优化参数获得凝聚的子结构有限元模型㊀㊀c a l c M a c r o M o d e l (边界条件);ʊ计算宏观设计域的有限元模型
㊀㊀f o r １:N
㊀㊀㊀i S e n s i t i v i t y M a c r o ();ʊ计算宏观结构每个子结构的设计变量灵敏度
㊀㊀u p d a t i n g M a c r o ();ʊ根据灵敏度值更新设计变量㊀㊀D e g i n M i c r o S t r u c t u r e ();ʊ微观结构设计㊀}
D e g i n M i c r o S t r u c t u r e ()ʊ微观结构设计㊀w h i l e
(判断是否收敛)㊀{
㊀㊀f o r １:N㊀{
㊀㊀㊀c a l c s u b s t r c u t u r e M i c r o (子结构节点编号);ʊ反求子结构内部节点位移场㊀㊀㊀i S e n s i t i v i t y
M i c r o ();ʊ计算宏观结构每个子结构的设计变量灵敏度㊀㊀㊀}㊀㊀u p d a t i n g
M i c r o ();ʊ根据灵敏度值更新微结构设计变量㊀
}
图２㊀结构设计流程
F i g ．２㊀T h e p r o c e s s o f s t r u c t u r a l d e s i g
n ㊀㊀该协同优化计算框架中,
每一步迭代均会判断宏微观结构是否收敛,而微观结构优化过程中,其材料约束受宏观子结构分布限制,因此,微观结构优化中需要从宏观结构中更新微观材料约束.
３㊀周期性结构优化实例
３．１㊀悬臂梁结构的设计(
实例１)以悬臂梁的计算为例,如图３所示,设计域大
小为L １ˑL ２＝２ˑ１,其左侧固定,右上角承受
F ＝１N 的垂直向下压力.所选用材料的弹性模量为E ０＝１,泊松比υ０＝０．３,体积约束为０．４.N x 与N y 分别是水平与垂直方向上的子结构数目
.
图３㊀优化模型(实例１
)F i g ．３㊀O p
t i m i z a t i o nm o d e l (e x a m p l e １)
４５８２ 中国机械工程第３３卷第２３期２０２２年１２月上半月
Copyright ©博看网. All Rights Reserved.
设计域划分为３２０ˑ１６０个四边形平面单元,
并选择不同大小的子结构进行宏观结构凝聚.宏观结构设计中,惩罚因子p 设置为２,过滤半径设置为与单元边长相等;微观结构设计中,惩罚因子
q 设置为３,
过滤半径设置为１．１.具体参数设置如表２所示.基于S I M P 方法,不考虑宏微观协同设计时,该模型的优化构型和柔度值c 如图４
所示.
表２㊀宏微观优化设计参数
T a b ．２㊀M a c r o Ｇm i c r o s c o p i c o p t i m i z a t i o nd e s i g
n p a r a m e t e r s
N x ˑN y
宏观
微观
惩罚因子p 过滤半径惩罚因子q 过滤半径子结构大小(m ˑm )２ˑ１２１３１．１１６０ˑ１６０８ˑ４
２１３１．１４０ˑ４０１６ˑ８２１３１．１２０ˑ２０３２ˑ１６
２
１
３
１．１
１０ˑ１０
图４㊀基于传统S I M P 方法的优化设计结构
F i g ．４㊀T h e o p t i m i z a t i o nd e s i g
n s t r u c t u r e u n d e r t r a d i t i o n a l S I M Pm e t h o d
㊀㊀在不同子结构大小下的优化结果如图５所
示.从优化结果可看出,设计域划分的子结构越多则结构杆系特征越明显,随着子结构数量的增多,其宏观结构的柔度值先增大后减小.
㊀㊀(a )N x ˑN y ＝２ˑ１,㊀㊀㊀(b )N x ˑN y ＝８ˑ４,㊀㊀㊀㊀c ＝７９．７２㊀㊀㊀㊀㊀㊀㊀㊀c ＝９０．０３
㊀(c )N x ˑN y ＝１６ˑ８,㊀㊀㊀(d )N x ˑN y ＝３２ˑ１６,㊀㊀㊀c ＝８６．４６㊀㊀㊀㊀㊀㊀㊀㊀㊀c ＝６７．４１
图５㊀不同子结构大小的优化结果(实例１
)F i g ．５㊀T h e o p
t i m i z a t i o n r e s u l t sw i t hd i f f e r e n t s i z e s o f s u b s t r u c t u r e (e x a m p
l e １)微观结构优化过程中,随着宏观结构的杆系
特征越来越明显,材料逐渐向微结构中集中,使得微观结构周期性子结构中的材料增多,微结构的杆系特征越来越不明显.
随着设计域中的子结构数量增多,通过式
(１１
)可得出宏观结构的材料会逐渐集中到各个子结构中,使得周期子结构单胞中的材料逐步增多.对比图５c 与图５d ,宏观结构的杆系特征越明显,材料在子结构单胞中的集中程度越大.由此可以推导出,在利用所提出的周期性子结构设计方法进行结构设计时,需要平衡宏观结构划分和周期子结构单胞大小.
３．２㊀固端梁的结构优化设计(
实例２)如图６所示的固端梁结构优化问题中,设计域大小为L １ˑL ２＝２ˑ１,其左右两侧固定,设计域中心受F ＝１N 的载荷.宏观结构划分为N x ˑ
N y ＝８ˑ４个子结构,子结构大小m ˑm 分别设置为１０ˑ１０㊁２０ˑ２０㊁３０ˑ３０㊁４０ˑ４０,体积约束为０．３.宏微观优化参数同上一实例
.
图６㊀优化模型(实例２
)F i g ．６㊀T h e d e s i g
nm o d e l (e x a m p l e ２)在不同子结构大小下,该优化模型的设计结果如图７所示.从优化构型来看,宏观结构基本相似,其结构柔度值随着子结构的增大先增大后减小.从微观构型来看,无论子结构尺寸的大小,均可得到相似的微观构型,随着子结构尺寸的增大,构型越来越清晰,其周期性子结构如图８所示.
㊀㊀㊀(a )m ˑm ＝１０ˑ１０,㊀㊀㊀㊀(b )m ˑm ＝２０ˑ２０,㊀㊀㊀㊀㊀㊀c ＝１２．９７㊀㊀㊀㊀㊀㊀㊀㊀㊀c ＝１３．５５
㊀㊀㊀(c )m ˑm ＝３０ˑ３０,㊀㊀㊀㊀(d )m ˑm ＝４０ˑ４０,㊀㊀㊀㊀㊀㊀c ＝１３．５９㊀㊀㊀㊀㊀㊀㊀㊀㊀c ＝１０．７４
图７㊀不同子结构大小的优化结果(实例２
)F i g ．７㊀T h e o p
t i m i z a t i o n r e s u l t sw i t hd i f f e r e n t s i z e s o f s u b s t r u c t u r e (e x a m p
l e ２)从优化目标函数分析,本文所提的周期性结
构设计方法所得结构柔度值均小于A R S P 方
法[１
]所得结构柔度值.相比于A R S P 方法,本文方法由于扩展了子结构内部的设计空间,可以得
５５８２ 基于子结构的宏微结构协同优化方法
吴紫俊㊀肖人彬Copyright ©博看网. All Rights Reserved.
(a )m ˑm ＝１０ˑ１０㊀㊀㊀(b )m ˑm ＝２０ˑ２０
(c )m ˑm ＝３０ˑ３０㊀㊀㊀(d )m ˑm ＝４０ˑ４０
图８㊀周期性微观结构构型F i g
．８㊀P e r i o d i cm i c r o s t r u c t u r e 到柔度更小的优化构型.
在宏微观协同优化过程中,由于微观结构优化嵌入到了宏观结构优化中,为了得到清晰的微观构型,微观优化变量迭代步长应小于宏观优化步长,避免迭代步长太大使得微观构型变化剧烈,从而引起宏观优化时的收敛速度缓慢的问题.图７d 的优化收敛过程如图９所示.在迭代计算中,
最后一步的柔度值有较大程度的减小,这是由于在该迭代步的前一步已经得到了最优构型,即该迭代步不属于优化过程的迭代步,而是为了获得优化构型的实际柔度值,通过有限元方法进行结构的直接计算,即惩罚因子设置为１时获得的结构柔度值.通过S I M P 方法,在整个优化过程中宏观材料含量不变,其材料在微观周期性单胞构型和宏观周期性结构分布之间进行匹配,以同时获得宏微观构型
.
图９㊀图７d 的收敛过程F i g ．９㊀E v o l u t i o nh i s t o r y o f f i g
．７d 对比上述两个实例,为了避免由于宏观子结构数增多而导致的周期性子结构大小减小,进而导致材料在周期性子结构中过度集中,难以获得具有明显几何特征的周期性子结构,在利用所提方法进行周期性宏微观结构设计时,尽量采用固定周期性子结构单胞大小的形式进行优化计算,既可获得清晰的宏观构型,又可避免由于材料集中引起的周期子结构单胞结构特征不清晰的问题.
４㊀结论
本文提出了基于子结构法的周期性微观构型的宏微观协同设计方法,根据子结构凝聚和反求构建对应的宏微观优化模型.基于S I M P 方法的优化框架,利用子结构的自由度凝聚实现结构的宏观构型设计,利用子结构内部自由度反求来实现周期性微观构型设计,同时给出了该优化方法的灵敏度计算和变量更新迭代方法.通过两个实
例验证了所提方法的合理性.
相比于传统基于均匀化的宏微观协同优化方法,基于子结构的宏微观协同优化方法可利用子结构及其对应超单元实现宏观结构的快速迭代和微观结构的再设计.由于子结构是宏观设计域的一部分,网格㊁节点等是一一对应的,故解决了均匀化方法的尺度分离问题.同时通过子结构自由度凝聚提高了宏观结构的设计效率.
尽管本文基于子结构实现了周期性结构的宏微观协同设计,但仍然存在以下几个问题,需要进一步研究.一是优化效率需要进一步研究.通过构建子优化问题实现微观周期性结构的设计,子结构划分的网格数越多,则微结构构型越精细,但也制约着整个优化问题的计算效率.二是宏微观设计参数的匹配仍需要进一步研究.较大的迭代步长会加快结构设计的收敛速度,本文中宏观结构与微结构通过两个关联的优化问题分别设计,较小的子结构设计变量迭代步长将影响宏观结构的收敛,较大的宏观设计变量迭代步长可能会让子结构得不到合适的微观构型.
参考文献:
[１]㊀WUZ i j u n ,X I AL i a n g ,WA N GS h u t i n g
,e t a l ．T o Ｇp o l o g y O p t i m i z a t i o no f H i e r a r c h i c a lL a t t i c eS t r u c Ｇt u r e sw i t hS u b s t r u c t u r i n g [J ]．C o m p u t e rM e t h o d s i n A p p l i e d M e c h a n i c sa n d E n g i n e e r i n g
,２０１９,３４５:６０２Ｇ６１７．
[２]㊀WU Z i j
u n ,F A N F e i ,X I A O R e n b i n ,e ta l ．T h e S u b s t r u c t u r i n g Ｇb a s e d T o p o l o g y O p t i m i z a t i o n f o r M a x i m i z i n g t h e F i r s t E i g
e n v a l u e o
f H i e r a r c h i c a l L a t t i c eS t r u c t u r e [J ]．I n t e r n a t i o n a l J o u r n a l f o rN u Ｇ
m e r i c a lM e t h o d si n E n g i n e e r i n g
,２０２０,１２１(１３):２９６４Ｇ２９７８．
[３]㊀S I V A P U R A M R ,D U N N I N G P D ,K I M H A．
S i m u l t a n e o u s M a t e r i a la n dS t r u c t u r a lO p
t i m i z a t i o n
６５８２ 中国机械工程第３３卷第２３期２０２２年１２月上半月
Copyright ©博看网. All Rights Reserved.
b y M u l t i s
c a l eT o p o l o g y O p t i m i z a t i o n[J]．S t r u c t u r a l
a n d M u l t i d i s c i p l i n a r y O p t i m i z a t i o n,２０１６,５４:１２６７Ｇ
１２８１．
[４]㊀董小虎,丁晓红．基于自适应成长法的周期性加筋结构拓扑优化设计方法[J]．中国机械工程,２０１８,２９(１７):２０４５Ｇ２０５１．
D O N G X i a o h u,D I N G X i a o h o n g．T o p o l o g y O p t iＧ
m u m D e s i g n M e t h o df o rP e r i o d i cS t i f f e n e rS t r u cＧ
t u r e sB a s e do nA d a p t i v eG r o w t h M e t h o d[J]．C h i n a
M e c h a n i c a l E n g i n e e r i n g,２０１８,２９(１７):２０４５Ｇ２０５１．[５]㊀V I C E N T E W M,Z U O Z H,P A V A N E L L O R,e t
a l．C o n c u r r e n tT o p o l o g y O p t i m i z a t i o nf o r M i n i m iＧ
z i n g F r e q u e n c y R e s p o n s e s o fT w oＧl e v e lH i e r a r c h i c a l
S t r u c t u r e s[J]．C o m p u t e r M e t h o d s i n A p p l i e d M eＧ
c h a n i c s a n dE n g i n e e r i n g,２０１６,３０１:１１６Ｇ１３６．[６]㊀张横,丁晓红,沈磊,等．考虑连接性的三明治阻尼复合结构拓扑优化设计[J]．中国机械工程,２０２１,３２(２０):２４０３Ｇ２４１０．
Z HA N G H e n g,D I N GX i a o h o n g,S H E NL e i,e t a l．T o p o l o g y O p t i m i z a t i o no f S a n d w i c hD a m p i n g C o mＧ
p o s i t e S t r u c t u r ew i t hC o n n e c t i v e S t i f f n e s sP h a s e[J]．
C h i n aM e c h a n i c a l E n g i n e e r i n g,２０２１,３２(２０):２４０３Ｇ
２４１０．
[７]㊀WUJ,A A G E N,W E S T E R MA N N R,e t a l．I n f i l l O p t i m i z a t i o nf o r A d d i t i v e M a n u f a c t u r i n gＧa p p r o a cＧ
h i n g B o n eＧl i k eP o r o u sS t r u c t u r e s[J]．I E E E T r a n sＧ
a c t i o n so n V i s u a l i z a t i o n a n d C o m p u t e r G r a p h i c s,
２０１８,２４(２):１１２７Ｇ１１４０．
[８]㊀X I A Q i,S H IT i e l i n．AC a s c a d i cM u l t i l e v e lO p t i m iＧz a t i o nA l g o r i t h mf o r t h eD e s i g n o f C o m p o s i t e S t r u cＧ
t u r e sw i t hC u r v i l i n e a r F i b e r B a s e d o nS h e p a r d I n t e rＧ
p o l a t i o n[J]．C o m p o s i t eS t r u c t u r e s,２０１８,１８８:２０９Ｇ
２１９．
[９]㊀F UJ u n j i a n,X I AL i a n g,G A OL i a n g,e t a l．T o p o l oＧ
g y O p t i m i z a t i o n o fP e r i o d i c S t r u c t u r e s w i t h S u bＧ
s t r u c t u r i n g[J]．J o u r n a l o f M e c h a n i c a l D e s i g n,２０１９,１４１(７):０７１４０３．
[１０]㊀杜义贤,尹鹏,李荣,等．兼具吸能和承载特性的梯度结构宏细观跨尺度拓扑优化设计[J]．机械工
程学报,２０２０,５６(７):１８５Ｇ１９４．
D U Y i x i a n,Y I NP e n g,L IR o n g,e t a l．M a c r o a n d
M i c r oT r a n sＧs c a l eT o p o l o g i c a l O p t i m i z a t i o nD e s i g n
o fG r a d i e n tS t r u c t u r e w i t h B o t h E n e r g y A b s o r pＧ
t i o na n dL o a dＧb e a r i n g C h a r a c t e r i s t i c s[J]．J o u r n a l
o fM e c h a n i c a l E n g i n e e r i n g,２０２０,５６(７):１８５Ｇ１９４．[１１]㊀廖中源,王英俊,王书亭．基于拓扑优化的变密度点阵结构体优化设计方法[J]．机械工程学报,
２０１９,５５(８):６５Ｇ７２．
L I A OZ h o n g y u a n,WA N G Y i n g j u n,WA N GS h uＧ
t i n g．G r a d e dＧd e n s i t y L a t t i c e S t r u c t u r e O p t i m i z aＧ
t i o nD e s i g n B a s e do n T o p o l o g y O p t i m i z a t i o n[J]．
J o u r n a l o fM e c h a n i c a lE n g i n e e r i n g,２０１９,５５(８):
６５Ｇ７２．
[１２]㊀B E N D SØE M P,K I K U C H IN．G e n e r a t i n g O p t iＧm a lT o p o l o g i e s i nS t r u c t u r a l D e s i g nU s i n g aH o mＧ
o g e n i z a t i o n M e t h o d[J]．C o m p u t e rM e t h o d s i nA pＧ
p l i e d M e c h a n i c sa n d E n g i n e e r i n g,１９８８,７１(２):
１９７Ｇ２２４．
[１３]㊀B E N D SØE M P,S I GMU N D O．M a t e r i a l I n t e r p oＧl a t i o nS c h e m e s i nT o p o l o g y O p t i m i z a t i o n[J]．A rＧ
c h i v eo f A p p l i e
d M
e c h a n i c s,１９９９,６９(９):６３５Ｇ
６５４．
[１４]㊀S I GMU N D O．M a t e r i a l sw i t hP r e s c r i b e dC o n s t i t uＧt i v eP a r a m e t e r s:a nI n v e r s eH o m o g e n i z a t i o nP r o bＧ
l e m[J]．I n t e r n a t i o n a l J o u r n a l o fS o l i d sa n dS t r u cＧ
t u r e s,１９９４,３１(１７):２３１３Ｇ２３２９．
[１５]㊀Z HUJ i h o n g,Z HA N G W e i h o n g,X I A L i a n g．T oＧp o l o g y O p t i m i z a t i o n i n A i r c r a f t a n d A e r o s p a c e
S t r u c t u r e sD e s i g n[J]．A r c h i v e so fC o m p u t a t i o n a l
M e t h o d s i nE n g i n e e r i n g,２０１６,２３(４):５９５Ｇ６２２．[１６]㊀B I E G L E R L T,L A N G Y i d o n g,L I N W e i j i e．M u l t iＧs c a l eO p t i m i z a t i o nf o rP r o c e s sS y s t e m sE nＧ
g i n e e r i n g[J]．C o m p u t e r s&C h e m i c a l E n g i n e e r i n g,
２０１４,６０:１７Ｇ３０．
[１７]㊀L IM i n,C H E N GZ h i b a o,J I A G a o f e n g,e t a l．D iＧm e n s i o nR e d u c t i o na n dS u r r o g a t eB a s e dT o p o l o g y
O p t i m i z a t i o no fP e r i o d i cS t r u c t u r e s[J]．C o m p o s i t e
S t r u c t u r e s,２０１９,２２９:１１１３８５．
[１８]㊀L I M i n,J I A G a o f e n g,C H E N G Z h i b a o,e ta l．
G e n e r a t i v eA d v e r s a r i a lN e t w o r k G u i d e d T o p o l o g y
O p t i m i z a t i o n o f P e r i o d i c S t r u c t u r e s v i a S u b s e t
S i m u l a t i o n[J]．C o m p o s i t eS t r u c t u r e s,２０２１,２６０:
１１３２５４．
[１９]㊀Z H E N G Y o n g f e n g,X I A O M i,G A OL i a n g,e t a l．R o b u s tT o p o l o g y O p t i m i z a t i o nf o rP e r i o d i cS t r u cＧ
t u r e sb y C o m b i n i n g S e n s i t i v i t y A v e r a g i n g w i t ha
S e m i a n a l y t i c a l M e t h o d[J]．I n t e r n a t i o n a lJ o u r n a l
f o rN u m e r i c a lM e t h o d s i nE n
g i n e e r i n g,２０１９,１１７
(５):４７５Ｇ４９７．
[２０]㊀H E G u a n g q i a n g,HU A N G X i a o d o n g,WA N G
H u,e t a l．T o p o l o g y O p t i m i z a t i o n o f P e r i o d i c
S t r u c t u r e s U s i n g B E S O B a s e d o n U n s t r u c t u r e d
D e s i g nP o i n t s[J]．S t r u c t u r a l a n d M u l t i d i s c i p l i n a r y
O p t i m i z a t i o n,２０１６,５３:２７１Ｇ２７５．
[２１]㊀R I V A E,C A Z Z U L A N IG,B E L L O N IE,e ta l．
A nO p t i m a lM e t h o d f o rP e r i o d i c S t r u c t u r e sD e s i g n
[C]ʊA S M E２０１７C o n f e r e n c e o nS m a r tM a t e r i a l s,
A d a p t i v e S t r u c t u r e s a n d I n t e l l i g e n t S y s t e m s．
S n o w b i r d,U t a h,２０１７:S MA S I S２０１７Ｇ３８３７．[２２]㊀C H E N W e n j i o n g,T O N G L i y o n g,L I U S h u t i a n．
C o n c u r r e n tT o p o l o g y
D e s i g no fS t r u c t u r e a n d M aＧ
t e r i a l U s i n g a T w oＧs c a l e T o p o l o g y O p t i m i z a t i o n
７５８２
基于子结构的宏微结构协同优化方法 吴紫俊㊀肖人彬
Copyright©博看网. All Rights Reserved.
[J]．C o m p u t e r s&S t r u c t u r e s,２０１７,１７８:１１９Ｇ
１２８．
[２３]㊀C H E N GG e n g d o n g,X UL i a n g．T w oＧs c a l eT o p o l oＧ
g y D e s i g n O p t i m i z a t i o n o f S t i f f e n e d o r P o r o u s
P l a t eS u b j e c t t o O u tＧo fＧp l a n eB u c k l i n g C o n s t r a i n t
[J]．S t r u c t u r a l a n dM u l t i d i s c i p l i n a r y O p t i m i z a t i o n,
２０１６,５４(５):１２８３Ｇ１２９６．
[２４]㊀D A H L B E R G O,M I T C H E L LＧT HOMA S R C, Q U E V E D OＧT E R U E L O．R e d u c i n g t h eD i s p e r s i o n
o fP e r i o d i cS t r u c t u r e sw i t hT w i s t a n dP o l a rG l i d e
S y mm e t r i e s[J]．S c i e n t i f i c R e p o r t s,２０１７,７:
１０１３６．
[２５]㊀赵清海,张洪信,华青松,等．周期性多材料结构稳态热传导拓扑优化设计[J]．工程力学,２０１９,
３６(３):２４７Ｇ２５６．
Z HA O Q i n g h a i,Z HA N G H o n g x i n,HU A Q i n gＧ
s o n g,e t a l．M u l t iＧm a t e r i a lT o p o l o g y O p t i m i z a t i o n
o fS t e a d yＧs t a t e H e a tC o n d u c t i o n S t r u c t u r eu n d e r
P e r i o d i c C o n s t r a i n t[J]．E n g i n e e r i n g M e c h a n i c s,
２０１９,３６(３):２４７Ｇ２５６．
[２６]㊀焦洪宇,李英,胡顺安,等．基于导重法的结构类周期性布局优化方法研究[J]．机械工程学报,
２０２０,５６(１３):２１８Ｇ２３０．
J I A O H o n g y u,L IY i n g,HUS h u n a n,e t a l．S t u d y
o f S t r u c t u r a l P e r i o d i cＧl i k e L a y o u t O p t i m i z a t i o n
B a s e do nG u i d eＧw e i g h tM e t h o d[J]．J o u r n a l o fM eＧ
c h a n i c a l E n g i n e e r i n g,２０２０,５６(１３):２１８Ｇ２３０．[２７]㊀焦洪宇,周奇才,李文军,等．基于变密度法的周期性拓扑优化[J]．机械工程学报,２０１３,４９(１３):
１３２Ｇ１３８．
J I A O H o n g y u,Z HO U Q i c a i,L I W e n j u n,e ta l．
P e r i o d i c T o p o l o g y O p t i m i z a t i o n U s i n g V a r i a b l e
D e n s i t y M e t h o d[J]．J o u r n a lo f M e c h a n i c a l
E n g iＧ
n e e r i n g,２０１３,４９(１３):１３２Ｇ１３８．
[２８]㊀杜义贤,李涵钊,田启华,等．基于能量均匀化的高剪切强度周期性点阵结构拓扑优化[J]．机械工
程学报,２０１７,５３(１８):１５２Ｇ１６０．
D U Y i x i a n,L IH a n z h a o,T I A N Q i h u a,e t a l．T oＧ
p o l o g y O p t i m i z a t i o no fP e r i o d i cL a t t i c eS t r u c t u r e
w i t h H i g h S h e a r S t r e n g t h U s i n g E n e r g yＧb a s e d
H o m o g e n i z a t i o n[J]．J o u r n a lo f M e c h a n i c a lE n g iＧ
n e e r i n g,２０１７,５３(１８):１５２Ｇ１６０．
[２９]㊀熊敏,丁晓红,季懿栋,等．树状分支传热结构层次生长优化设计技术[J]．中国机械工程,２０１９,
３０(２２):２６６８Ｇ２６７４．
X I O N G M i n,D I N G X i a o h o n g,J IY i d o n g,e ta l．
H i e r a r c h y G r o w t h M e t h o df o rO p t i m a lD e s i g no f
B r a n c h i n g H e a tT r a n s f e r S t r u c t u r e s[J]．
C h i n aM eＧ
c h a n i c a l E n g i n e e r i n g,２０１９,３０(２２):２６６８Ｇ２６７４．[３０]㊀肖人彬,赖荣燊,李仁旺．从大规模定制化设计到大规模个性化设计[J]．南昌工程学院学报,２０２１,
４０(１):１Ｇ１２．
X I A O R e n b i n,L A IR o n g s h e n,L i R e n w a n g．F r o m
D e s i g n f o rM a s sC u s t o m i z a t i o n t oD e s i g n f o rM a s s
P e r s o n a l i z a t i o n[J]．J o u r n a lo fN a n c h a n g I n s t i t u t e
o fT e c h n o l o g y,２０２１,４０(１):１Ｇ１２．
[３１]㊀李昂,刘初升．基于超单元技术的大型复杂结构的拓扑优化设计[J]．中国机械工程,２０１７,２８(２０):
２４６７Ｇ２４７４．
L I A n g,L I U C h u s h e n g．T o p o l o g y O p t i m i z a t i o n
D e s i g n s o f L a r g e a n dC o m p l e xS t r u c t u r e sB a s e d o n
S u p e rE l e m e n t T e c h n i q u e[J]．C h i n a M e c h a n i c a l
E n g i n e e r i n g,２０１７,２８(２０):２４６７Ｇ２４７４．
[３２]㊀付君健,孙鹏飞,杜义贤,等．基于子结构法的多层级结构拓扑优化[J]．中国机械工程,２０２１,３２
(１６):１９３７Ｇ１９４４．
F UJ u n j i a n,S U N P e n g f e i,D U Y i x i a n,e t a l．H iＧ
e r a r c h i c a lS t r u c t u r eT o p o l o g y O p t i m i z a t i o nB a s e d
o nS u b s t r u c t u r eM e t h o d[J]．C h i n aM e c h a n i c a l E nＧ
g i n e e r i n g,２０２１,３２(１６):１９３７Ｇ１９４４．
[３３]㊀B E N D SØE M P,S I GMU N DO．T o p o l o g y O p t i m iＧz a t i o n:T h e o r y,M e t h o d sa n d A p p l i c a t i o n s[M]．
B e r l i n:S p r i n g e r,２００３．
(编辑㊀王旻玥)
作者简介:吴紫俊,男,１９８５年生,讲师㊁博士.研究方向为结构优化与分析,C A D/C AM技术.EＧm a i l:z j w u＠w t u．e d u．c n.肖人彬(通信作者),男,１９６５年生,教授㊁博士研究生导师.研究方向为智能设计㊁复杂产品创新设计.EＧm a i l:r b x i a o＠h u s t．e d u．c n.
８５８２
中国机械工程第３３卷第２３期２０２２年１２月上半月Copyright©博看网. All Rights Reserved.。