拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识。通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] 。于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] 。如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架。
尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] 。为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] 。Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题。张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法。高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法。廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度。但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法。
为了解决体积约束和柔顺度最小的结构拓扑优化问题及采用变密度法引起的数值不稳定问题,以兼顾改进拓扑优化方法为前提,引入改进的固体各向同性材料惩罚(solid isotropic material with penalization, SIMP)法和基于高斯权重函数的敏度过滤法,并设计新的灰度单元抑制算子,对由变密度法和优化准则法形成的经典拓扑优化方法进行改进,提出了一种新的拓扑优化方法。
1 连续体结构拓扑优化
采用新的拓扑优方法研究连续体结构的优化设计问题,需要确定材料插值方案、建立结构拓扑优化模型、设计灵敏度、确定敏度过滤法、确定更新准则及给出评价指标等几个步骤。
1.1 基于改进SIMP 法的材料插值方案
采用能在0与1之间连续取值的基于SIMP法的材料插值方案可以有效降低结构优化设计的难度。参考Andreassen等[18 ] 的研究成果,目前应用较为普遍的材料插值方案为:
E i = E m i n + x i p E 0 - E m i n , x i ∈ 0,1 (1)
式中:E 0 为实体材料的弹性模量;E min 为无效材料的弹性模量;p 为无效材料的惩罚系数;xi 为拓扑设计变量,i =1,2,…,n;Ei 为当量弹性模量。
通过进一步研究发现,式(1)所示的插值方案对改进拓扑优化方法的效果有限。为了进一步改进拓扑优化方法,采用Zhu等[19 ] 提出的改进SIMP法构建材料插值方案,具体如下:
E i = E m i n + α - 1 α x i p + 1 α x i E 0 , x i ∈ 0,1 (2)
须说明的是,文献[19 ]提出的材料插值方案中并没有E min ,此处加入E min 的目的是规避奇异刚度矩阵的产生和确保xi 的下限值可取0。
1.2 结构拓扑优化模型
建立优化模型是开展结构拓扑优化设计的前提。在式(2)的基础上,参考文献[9 ,18 ],则使结构柔顺度最小的拓扑优化模型为:
f i n d : x = x 1 , x 2 , ⋯ , x n T m i n : c x = F T U = U T K U s . t : K U = F V x V 0 = V f x i ∈ 0,1 (3)
式中: x 为结构拓扑设计变量的矢量;c 为结构柔顺度; K 为结构的刚度矩阵; U 为载荷矢量 F 作用在结构上而产生的位移矢量;V 0 和V 分别为结构设计域的材料体积和结构优化设计后的材料体积;V f 为优化结构的许用材料体积分数。
1.3 灵敏度
采用经典优化准则法求解式(3)所示的优化模型时,须设计目标函数和约束关于拓扑设计变量xi 的灵敏度。根据文献[18 ],基于式(2)和式(3),求解c 和V 关于xi 的灵敏度:
∂ c ∂ x i = - α - 1 α p x i p - 1 + 1 α E 0 u i T k 0 u i ∂ V ∂ x i = 1 (4)
式中: u i 为单元的位移矩阵; k 0 为具有单位弹性模量的刚度矩阵。
1.4 基于高斯权重函数的敏度过滤法
为了有效解决采用变密度法引起的数值不稳定问题及改进拓扑优化方法,选用合适的敏度过滤法极为重要。Cheng等[20 ] 发现,采用基于高斯权重函数的敏度过滤法可以得到轻质且拓扑构型好的优化结构。本文采用基于高斯权重函数的敏度过滤法对c 关于xi 的灵敏度进行过滤处理,具体为:
∂ c ̂ ∂ x i = 1 m a x γ , x i ∑ j ∈ N i H i , j ∑ j ∈ N i H i , j x j ∂ c ∂ x j (5)
式中:γ 为为了规避当xi 取为0时式中分母为0而引入的一个较小正值参数;Ni 为与单元i 相邻的单元j 所组成的集合;Hi , j 为敏度过滤的权重系数。
H i , j = e x p - 1 2 d i , j σ 2 (6)
式中:di,j 为单元i 、j 之间的距离;σ 为高斯参数,与过滤半径R 有关。
1.5 基于新灰度单元抑制算子的更新准则
基于优化准则法求解式(3)时,需要一种启发式更新准则。根据文献[18 ],常用的更新准则为:
x i n e w = m a x 0 , x i - m , x i B i η ≤ m a x 0 , x i - m m i n 1 , x i + m , x i B i η ≥ m i n 1 , x i - m x i B i η , 其 他 (7)
式中:m 和η 分别为正值移动步长和阻尼系数;xi 与Bi η 的结合物xi Bi η 可以视为简单的灰度单元抑制算子,用于更新xi 。
B i = - ∂ c ∂ x i λ ∂ V ∂ x i (8)
为了有效抑制灰度单元的产生,受Groenwold等[21 ] 提出的灰度单元抑制算子的启发,在原灰度单元抑制算子(xi Bi η ) q 的基础上,引入正值比例参数β ,设计了一种新的灰度单元抑制算子,为:
Γ x i B i η = 1 - β x i B i η q + β x i B i η (9)
为了提高拓扑优化方法的收敛速度,在综合考虑每次迭代计算所得的结构柔顺度的基础上,设计了新的参数q 的更新方案,具体为:
q = q 0 , c > c * ξ q + Δ q , c ≤ c * (10)
式中:q 0 和Δq 分别为q 的初值和增量;ξ 和c * 分别为强度因子和结构柔顺度的临界值。
须说明的是,通过测试发现,c * 值设置合理与否会对优化结构的柔顺度、拓扑构型及拓扑优化方法的收敛速度产生较大影响。
根据式(7)和式(9),基于新灰度单元抑制算子的更新准则为:
x i n e w = m a x 0 , x i - m , Γ x i B i η ≤ m a x 0 , x i - m m i n 1 , x i + m , Γ x i B i η ≥ m i n 1 , x i - m Γ x i B i η , 其 他 (11)
1.6 评价指标
选择结构的柔顺度c 、程序循环执行次数I 和灰度单元比例指数G 作为评价拓扑优化方法的指标。其中:c 用于反比表征获取柔顺度小的优化结构的能力; I 用于反比表征收敛速度或收敛性能;G 用于反比表征抑制灰度单元产生的能力。
G = n * n × 100 % (12)
式中:n * 为满足约束条件(xi < 0.01或xi > 0.99)的单元数量。
2 算例分析与讨论
采用2个经典的数值算例,分别测试改进SIMP法、基于高斯权重函数的敏度过滤法、新灰度单元抑制算子和新拓扑优化方法对其的优化效果,并与采用文献[14 ,8 ]提出的拓扑优化方法的优化结果进行对比。设置与文献[14 ,18 ]相同的基本参数:E 0 =1 MPa,E min = 1×10-9 MPa,材料泊松比υ = 0.3,p = 3,m = 0.2,η = 1×10-3 。采用双线性矩形单元对结构进行网格划分。另外,须特别说明的是,本文设计的拓扑优化方法的程序是在文献[18 ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上。
2.1 算例1
MBB (Messerschmitt-Bolkow-Blohm)梁结构如图1 所示。其几何尺寸为280 mm ×50 mm,顶端中部承受大小为2 N的外载荷 F 。对MBB梁结构进行拓扑优化,来测试采用改进SIMP法、基于高斯权重函数的敏度过滤法、新灰度单元抑制算子和新拓扑优化方法的优化效果。根据对称性,采用具有1/2设计域的MBB梁为优化对象。对梁结构进行网格划分,可以得到单元总数量为140×50个的有限元模型。设置R =3.1,V f =0.4。
图1
图1
MBB梁结构
Fig.1
MBB beam structure
不同SIMP法下MBB梁结构的拓扑构型和优化结果分别如图2 和表1 所示;在本文SIMP法下,采用不同敏度过滤法得到的结构拓扑构型和优化结果分别如图3 和表2 所示;在本文SIMP法和敏度过滤法下,采用不同灰度单元抑制算子得到的结构拓扑构型和优化结果分别如图4 和表3 所示;在不同拓扑优化方法下结构拓扑构型和优化结果分别如图5 和表4 所示。
图2
图2
不同SIMP法下MBB梁结构的拓扑构型
Fig.2
Topological configuration of MBB beam structure under different SIMP methods
图3
图3
不同敏度过滤法下MBB梁结构的拓扑构型
Fig.3
Topological configuration of MBB beam structure under different sensitivity filtration methods
图4
图4
不同灰度单元抑制算子下MBB梁结构的拓扑构型
Fig.4
Topological configuration of MBB beam structure under different gray-scale unit suppression operators
图5
图5
不同拓扑优化方法下MBB梁结构的拓扑构型
Fig.5
Topology configuration of MBB beam structure under different topology optimization methods
由图2 可知,采用不同SIMP法得到的MBB梁结构的拓扑构型存在一定差别:采用本文SIMP法所得优化结构在“区域A”和“区域B”处的颜色比较浅。意味着采用文献[18 ]的SIMP法所得结构在这些区域的xi 值比较大。
由表1 可知,相比于文献[18 ]的SIMP法,采用本文SIMP法所得结构的c 和I 分别减少了2.09 %和22.89 %,G 增加了13.80 %。说明采用本文SIMP法可以使拓扑优化方法的收敛速度变得更快,而且可以增强获取柔顺度小的优化结构的能力,但也会削弱抑制灰度单元产生的能力。
由图3 可知,图3 (a)中“区域A”和“区域B”处的形状不同于图3 (b),且在“区域A”处的区别更为明显,说明使用不同的敏度过滤法会得到不同的结构拓扑构型。
由表2 可知,相比于文献[18 ]的敏度过滤法,采用本文敏度过滤法所得的c 和G 分别减少了1.30 %和10.25 %,I 则相同。说明采用本文的敏度过滤法可以增强获取柔顺度小的优化结构和抑制灰度单元产生的能力。
由图4 可知,采用本文灰度单元抑制算子所得优化结构的拓扑构型更为清晰。
由表3 可知,相比于文献[18 ]的灰度单元抑制算子,采用本文灰度单元抑制算子所得的c 、G 和I 分别减少了4.46 %、85.71 %和45.31 %。说明采用本文灰度单元抑制算子不仅可以增强获取柔顺度小的优化结构和抑制灰度单元产生的能力,而且可以提高收敛速度。
由图5 可知,采用文献[14 ,18 ]方法所得结构的拓扑构型存在更多的灰度单元,而采用本文方法所得结构的轮廓非常清晰,且从宏观上看不到灰度单元,所得结构比较好。
由表4 可知,相比于文献[14 ,18 ],采用本文拓扑优化方法所得的c 分别减少了7.14 %和7.68 %,G 分别减少了83.97 %和85.41 %,I 分别减少了51.39 %和57.83 %。说明本文方法在收敛速度、获取柔顺度小和拓扑构型好的优化结构及抑制灰度单元产生等方面具有明显优势。
2.2 算例2
悬臂梁结构如图6 所示。其几何尺寸为160 mm ×70 mm,梁的左端被固定,右端中部承受大小为1 N的竖直向下的外载荷 F 。网格划分后结构的单元总数量为160 × 70 个。设置R = 3.5,V f = 0.4。采用3种拓扑优化方法所得悬臂梁结构的拓扑构型和优化结果分别如图7 和表5 所示。
图6
图6
悬臂梁结构
Fig.6
Cantilever beam structure
图7
图7
不同拓扑优化方法下悬臂梁结构的拓扑构型
Fig.7
Topological configuration of cantilever beam structure under different topology optimization methods
由图7 可知,采用不同拓扑优化方法所得结构的拓扑构型存在较明显差异。相比于文献[14 ,18 ],采用本文方法所得结构的轮廓更为清晰。由表5 可知,采用本文方法所得的c 、I 和G 分别为104.98 N·mm、38次和4.73 %,较文献[14 ]方法分别减少了8.63 %、67.80 %和80.90 %,较文献[18 ]方法分别减少了8.90 %、77.25 %和82.32 %。说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快。
3 结 论
为了实现使连续体结构的体积约束和柔顺度最小的拓扑优化及解决采用经典变密度法引起的数值不稳定问题,在兼顾改进拓扑优化方法的前提下,通过引入改进的SIMP法、基于高斯权重函数的敏度过滤法和新灰度单元抑制算子,并结合优化准则法,提出了新的拓扑优化方法。通过优化算例可知,新的拓扑优化方法能够有效解决连续体结构的拓扑优化问题,具有收敛速度较快、可更好地获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等优势。
参考文献
View Option
[1]
王辉 .不确定性连续体结构的拓扑优化研究
[D].成都 :西南交通大学 ,2021 :1 -154 .
[本文引用: 1]
WANG Hui . Research on topology optimization of uncertain continuum structure
[D]. Chengdu : Southwest Jiaotong University , 2021 : 1 -154 .
[本文引用: 1]
[2]
徐浩然 ,贺福强 ,李赟 ,等 .飞机起落架的拓扑与自由曲面形状优化
[J].组合机床与自动化加工技术 ,2021 (4 ):134 -138 .
[本文引用: 1]
XU Hao-ran , HE Fu-qiang , LI Yun , et al . Topology and free-form surface shape optimization of aircraft landing gear
[J]. Modular Machine Tool & Automatic Manufacturing Technique , 2021 (4 ): 134 -138 .
[本文引用: 1]
[3]
HU Song , HUANG Yong , ZHANG Lu . Optimized topology design of gypsum embedded in filler cast-in-place concrete hollow floor based on ANSYS
[J]. Guizhou Science , 2015 , 33 (1 ): 25 -29 .
DOI:10.3969/j.issn.1003-6563.2015.01.005
[4]
JIAO H Y , LI F , JIANG Z Y , et al . Periodic topology optimization of a stacker crane
[J]. IEEE Access , 2019 , 7 : 186553 -186562 . doi:10.1109/access.2019.2960327
[5]
SHEN W , OHSAKI M . Geometry and topology optimization of plane frames for compliance minimization using force density method for geometry model
[J]. Engineering with Computers , 2021 , 37 (3 ): 2029 -2046 .
[本文引用: 1]
[9]
WANG H , CHENG W M , DU R , et al . Improved proportional topology optimization algorithm for solving minimum compliance problem
[J]. Structural and Multidisciplinary Optimization , 2020 , 62 (2 ): 475 -493 . doi:10.1007/s00158-020-02504-8
[本文引用: 3]
[10]
SIGMUND O . Morphology-based black and white filters for topology optimization
[J]. Structural and Multidisciplinary Optimization , 2007 , 33 (4/5 ): 401 -424 . doi:10.1007/s00158-006-0087-x
[本文引用: 2]
[11]
匡兵 ,刘娟 ,段君伟 ,等 .基于改进灵敏度过滤策略的SIMP方法
[J].计算力学学报 ,2017 ,34 (1 ):81 -87 .
KUANG Bing , LIU Juan , DUAN Jun-wei , et al . SIMP method based on modified sensitivity filtering strategy
[J]. Chinese Journal of Computational Mechanics , 2017 , 34 (1 ): 81 -87 .
[12]
WEI T , YANG X , CHEN C , et al . Grey filter functions for suppression of grey-scale elements
[J]. Engineering Optimization , 2019 , 51 (2 ): 317 -331 . doi:10.1080/0305215x.2018.1454441
[13]
[本文引用: 1]
LIAN Rui-chao , JING Shi-kai , YANG Hai-cheng , et al . Sensitivity filtering method considering partition blending weights
[J]. Journal of Computer-Aided Design & Computer Graphics , 2019 , 31 (5 ): 842 -850 .
DOI:10.3724/sp.j.1089.2019.17336
[本文引用: 1]
[14]
BIYIKLI E , TO A C . Proportional topology optimization: A new non-sensitivity method for solving stress constrained and minimum compliance problems and its implementation in MATLAB
[J]. PLoS One , 2015 , 10 (12 ): e0145041 . doi:10.1371/journal.pone. 0145041
[本文引用: 10]
[15]
张国锋 ,徐雷 ,李大双 ,等 .连续体结构拓扑优化敏度过滤研究
[J].组合机床与自动化加工技术 ,2021 (6 ):29 -32 .
[本文引用: 1]
ZHANG Guo-feng , XU Lei , LI Da-shuang , et al . Research on sensitivity filtering of continuum topology optimization
[J]. Modular Machine Tool & Automatic Manufacturing Technique , 2021 (6 ): 29 -32 .
[本文引用: 1]
[16]
[本文引用: 1]
GAO Xiang , WANG Lin-jun , DU Yi-xian . An improved topology optimization method for suppressing gray elements
[J]. Journal of Computer-Aided Design & Computer Graphics , 2020 , 32 (12 ): 2003 -2012 .
DOI:10.3724/sp.j.1089.2020.18208
[本文引用: 1]
[17]
廉睿超 ,敬石开 ,何志军 ,等 .拓扑优化变密度法的灰度单元分层双重惩罚方法
[J].计算机辅助设计与图形学学报 ,2020 ,32 (8 ):1349 -1356 ,1227 .
[本文引用: 1]
LIAN Rui-chao , JING Shi-kai , HE Zhi-jun , et al . A hierarchical double penalty method of gray-scale elements for SIMP in topology optimization
[J]. Journal of Computer-Aided Design & Computer Graphics , 2020 , 32 (8 ): 1349 -1356 , 1227 .
[本文引用: 1]
[18]
ANDREASSEN E , CLAUSEN A , SCHEVENELS M , et al . Efficient topology optimization in MATLAB using 88 lines of code
[J]. Structural and Multidisciplinary Optimization , 2011 , 43 (1 ): 1 -16 . doi:10.1007/s00158-010-0594-7
[本文引用: 21]
[19]
ZHU H J , ZHANG W H . Integrated layout design of supports and structures
[J]. Computer Methods in Applied Mechanics and Engineering , 2010 , 199 (9/12 ): 557 -569 . doi:10.1016/j.cma.2009.10.011
[本文引用: 2]
[20]
CHENG W , WANG H , ZHANG M , et al . Improved proportional topology optimization algorithm for minimum volume problem with stress constraints
[J]. Engineering Computations , 2021 , 38 (1 ): 392 -412 . doi:10.1108/ec-12-2019-0560
[本文引用: 2]
[21]
GROENWOLD A A , ETMAN L F P . A simple heuristic for gray-scale suppression in optimality criterion-based topology optimization
[J]. Structural and Multidisciplinary Optimization , 2009 , 39 (2 ): 217 -225 . doi:10.1007/s00158-008-0337-1
[本文引用: 1]
不确定性连续体结构的拓扑优化研究
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
不确定性连续体结构的拓扑优化研究
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
飞机起落架的拓扑与自由曲面形状优化
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
飞机起落架的拓扑与自由曲面形状优化
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
基于ANSYS的石膏空腔模无梁楼盖拓扑优化设计
0
2015
基于ANSYS的石膏空腔模无梁楼盖拓扑优化设计
0
2015
Periodic topology optimization of a stacker crane
0
2019
Geometry and topology optimization of plane frames for compliance minimization using force density method for geometry model
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
时速400公里高速列车底架拓扑优化
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
时速400公里高速列车底架拓扑优化
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
档位互换机构壳体拓扑优化
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
档位互换机构壳体拓扑优化
1
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
基于变密度法的送杆机构的运送支架拓扑优化设计
2
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
... 采用2个经典的数值算例,分别测试改进SIMP法、基于高斯权重函数的敏度过滤法、新灰度单元抑制算子和新拓扑优化方法对其的优化效果,并与采用文献[14 ,8 ]提出的拓扑优化方法的优化结果进行对比.设置与文献[14 ,18 ]相同的基本参数:E 0 =1 MPa,E min = 1×10-9 MPa,材料泊松比υ = 0.3,p = 3,m = 0.2,η = 1×10-3 .采用双线性矩形单元对结构进行网格划分.另外,须特别说明的是,本文设计的拓扑优化方法的程序是在文献[18 ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上. ...
基于变密度法的送杆机构的运送支架拓扑优化设计
2
2021
... 拓扑优化是实现结构轻量化的一种有效的设计方法,蕴含如现代力学、现代数学、计算机科学等多学科知识.通过拓扑优化可产生材料分布合理、性能各异、质量最小的结构设计方案,故其已成为结构优化设计的热点研究方向[1 ] .于20世纪90年代提出的备受关注的变密度法作为连续体结构拓扑优化设计的一种主流方法,因自身所具备的技术优势而在航空航天、建筑、机械等领域的结构优化设计中得到广泛应用[2 -5 ] .如:为了实现高速列车底架的轻量化设计,陈秉智等[6 ] 采用基于OptiStruct软件的变密度法对高速列车的底架进行拓扑优化设计,获得了可同时满足强度、刚度及模态频率等设计要求且应力分布得到改善的优质轻质底架;为了得到轻质和性能完善的汽车用挡位互换机构,王铭昭等[7 ] 也采用基于OptiStruct软件的变密度法对挡位互换机构的壳体进行了拓扑优化设计;周围等[8 ] 借助基于ANSYS Workbench软件的变密度法研究了岩心钻机送杆机构运送支架的拓扑优化设计问题,得到了新型轻质运送支架. ...
... 采用2个经典的数值算例,分别测试改进SIMP法、基于高斯权重函数的敏度过滤法、新灰度单元抑制算子和新拓扑优化方法对其的优化效果,并与采用文献[14 ,8 ]提出的拓扑优化方法的优化结果进行对比.设置与文献[14 ,18 ]相同的基本参数:E 0 =1 MPa,E min = 1×10-9 MPa,材料泊松比υ = 0.3,p = 3,m = 0.2,η = 1×10-3 .采用双线性矩形单元对结构进行网格划分.另外,须特别说明的是,本文设计的拓扑优化方法的程序是在文献[18 ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上. ...
Improved proportional topology optimization algorithm for solving minimum compliance problem
3
2020
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
... 建立优化模型是开展结构拓扑优化设计的前提.在式(2) 的基础上,参考文献[9 ,18 ],则使结构柔顺度最小的拓扑优化模型为: ...
... 参考文献[9 ],可得: ...
Morphology-based black and white filters for topology optimization
2
2007
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
... 根据文献[10 ,20 ],得到: ...
基于改进灵敏度过滤策略的SIMP方法
0
2017
基于改进灵敏度过滤策略的SIMP方法
0
2017
Grey filter functions for suppression of grey-scale elements
0
2019
考虑分区混合权重的敏度过滤方法
1
2019
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
考虑分区混合权重的敏度过滤方法
1
2019
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
Proportional topology optimization: A new non-sensitivity method for solving stress constrained and minimum compliance problems and its implementation in MATLAB
10
2015
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
... 采用2个经典的数值算例,分别测试改进SIMP法、基于高斯权重函数的敏度过滤法、新灰度单元抑制算子和新拓扑优化方法对其的优化效果,并与采用文献[14 ,8 ]提出的拓扑优化方法的优化结果进行对比.设置与文献[14 ,18 ]相同的基本参数:E 0 =1 MPa,E min = 1×10-9 MPa,材料泊松比υ = 0.3,p = 3,m = 0.2,η = 1×10-3 .采用双线性矩形单元对结构进行网格划分.另外,须特别说明的是,本文设计的拓扑优化方法的程序是在文献[18 ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上. ...
... ]提出的拓扑优化方法的优化结果进行对比.设置与文献[14 ,18 ]相同的基本参数:E 0 =1 MPa,E min = 1×10-9 MPa,材料泊松比υ = 0.3,p = 3,m = 0.2,η = 1×10-3 .采用双线性矩形单元对结构进行网格划分.另外,须特别说明的是,本文设计的拓扑优化方法的程序是在文献[18 ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上. ...
... Optimization results of MBB beam structure under different topology optimization methods
Table 4 拓扑优化方法 c /(N·mm)I /次G / %文献[14 ] 方法 211.76 72 24.27 文献[18 ] 方法 213.00 83 26.66 本文方法 196.65 35 3.89
由图2 可知,采用不同SIMP法得到的MBB梁结构的拓扑构型存在一定差别:采用本文SIMP法所得优化结构在“区域A”和“区域B”处的颜色比较浅.意味着采用文献[18 ]的SIMP法所得结构在这些区域的xi 值比较大. ...
... 由图5 可知,采用文献[14 ,18 ]方法所得结构的拓扑构型存在更多的灰度单元,而采用本文方法所得结构的轮廓非常清晰,且从宏观上看不到灰度单元,所得结构比较好. ...
... 由表4 可知,相比于文献[14 ,18 ],采用本文拓扑优化方法所得的c 分别减少了7.14 %和7.68 %,G 分别减少了83.97 %和85.41 %,I 分别减少了51.39 %和57.83 %.说明本文方法在收敛速度、获取柔顺度小和拓扑构型好的优化结构及抑制灰度单元产生等方面具有明显优势. ...
... Optimization results of cantilever beam structure under different topology optimization methods
Table 5 拓扑优化方法 c /(N·mm)I /次G / %文献[14 ] 方法 114.90 118 24.77 文献[18 ] 方法 115.23 167 26.75 本文方法 104.98 38 4.73
由图7 可知,采用不同拓扑优化方法所得结构的拓扑构型存在较明显差异.相比于文献[14 ,18 ],采用本文方法所得结构的轮廓更为清晰.由表5 可知,采用本文方法所得的c 、I 和G 分别为104.98 N·mm、38次和4.73 %,较文献[14 ]方法分别减少了8.63 %、67.80 %和80.90 %,较文献[18 ]方法分别减少了8.90 %、77.25 %和82.32 %.说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快. ...
... 由图7 可知,采用不同拓扑优化方法所得结构的拓扑构型存在较明显差异.相比于文献[14 ,18 ],采用本文方法所得结构的轮廓更为清晰.由表5 可知,采用本文方法所得的c 、I 和G 分别为104.98 N·mm、38次和4.73 %,较文献[14 ]方法分别减少了8.63 %、67.80 %和80.90 %,较文献[18 ]方法分别减少了8.90 %、77.25 %和82.32 %.说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快. ...
... 分别为104.98 N·mm、38次和4.73 %,较文献[14 ]方法分别减少了8.63 %、67.80 %和80.90 %,较文献[18 ]方法分别减少了8.90 %、77.25 %和82.32 %.说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快. ...
... ]方法分别减少了8.90 %、77.25 %和82.32 %.说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快. ...
连续体结构拓扑优化敏度过滤研究
1
2021
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
连续体结构拓扑优化敏度过滤研究
1
2021
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
改进的抑制灰度单元的拓扑优化方法
1
2020
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
改进的抑制灰度单元的拓扑优化方法
1
2020
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
拓扑优化变密度法的灰度单元分层双重惩罚方法
1
2020
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
拓扑优化变密度法的灰度单元分层双重惩罚方法
1
2020
... 尽管变密度法有诸多优点,但所得优化结构普遍存在如棋盘格、网格依赖及灰度单元等数值不稳定问题[9 ] .为了解决此类问题,学者们提出了如密度过滤法、敏度过滤法、灰度单元抑制算子等多种方法[10 -13 ] .Biyikli等[14 ] 采用一种简单有效的改进密度过滤法解决了采用变密度法引起的结构棋盘格问题.张国锋等[15 ] 为了消除采用变密度法引起的棋盘格、网格依赖和边界扩散等问题,在保留原卷积因子的基础上,通过引入新卷积因子,提出了一种改进敏度过滤方法.高翔等[16 ] 为了解决采用变密度法引起的灰度单元和现有Heaviside函数须频繁调用密度矩阵的问题,结合优化准则法和Heaviside函数,设计了基于改进Heaviside函数的拓扑优化方法.廉睿超等[17 ] 从变密度法的材料插值方案出发,提出了一种灰度单元分层双重惩罚方法,在有效抑制灰度单元产生的同时提高了拓扑优化的收敛速度.但是,目前所提出的策略或方法尚不能很好地兼顾解决多类型的数值不稳定问题与改进拓扑优化方法. ...
Efficient topology optimization in MATLAB using 88 lines of code
21
2011
... 采用能在0与1之间连续取值的基于SIMP法的材料插值方案可以有效降低结构优化设计的难度.参考Andreassen等[18 ] 的研究成果,目前应用较为普遍的材料插值方案为: ...
... 建立优化模型是开展结构拓扑优化设计的前提.在式(2) 的基础上,参考文献[9 ,18 ],则使结构柔顺度最小的拓扑优化模型为: ...
... 采用经典优化准则法求解式(3) 所示的优化模型时,须设计目标函数和约束关于拓扑设计变量xi 的灵敏度.根据文献[18 ],基于式(2) 和式(3) ,求解c 和V 关于xi 的灵敏度: ...
... 基于优化准则法求解式(3) 时,需要一种启发式更新准则.根据文献[18 ],常用的更新准则为: ...
... 采用2个经典的数值算例,分别测试改进SIMP法、基于高斯权重函数的敏度过滤法、新灰度单元抑制算子和新拓扑优化方法对其的优化效果,并与采用文献[14 ,8 ]提出的拓扑优化方法的优化结果进行对比.设置与文献[14 ,18 ]相同的基本参数:E 0 =1 MPa,E min = 1×10-9 MPa,材料泊松比υ = 0.3,p = 3,m = 0.2,η = 1×10-3 .采用双线性矩形单元对结构进行网格划分.另外,须特别说明的是,本文设计的拓扑优化方法的程序是在文献[18 ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上. ...
... .采用双线性矩形单元对结构进行网格划分.另外,须特别说明的是,本文设计的拓扑优化方法的程序是在文献[18 ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上. ...
... ]所提出方法代码的基础上编写的,本文方法与文献[18 ]方法的区别在于材料插值方案、敏度过滤法和灰度单元抑制算子这几点上. ...
... Optimization results of MBB beam structure under different SIMP methods
Table 1 SIMP法 c /(N·mm)I /次G / %文献[18 ]方法 213.00 83 26.66 本文方法 208.54 64 30.34
图3 不同敏度过滤法下MBB梁结构的拓扑构型 Topological configuration of MBB beam structure under different sensitivity filtration methods Fig.3 ![]()
表2 不同敏度过滤法下MBB 梁结构优化结果 ...
... Optimization results of MBB beam structure under different sensitivity filtration methods
Table 2 敏度过滤法 c /(N·mm)I /次G / %文献[18 ] 方法 208.54 64 30.34 本文方法 205.83 64 27.23
图4 不同灰度单元抑制算子下MBB梁结构的拓扑构型 Topological configuration of MBB beam structure under different gray-scale unit suppression operators Fig.4 ![]()
表3 不同灰度单元抑制算子下MBB 梁结构优化结果 ...
... Optimization results of MBB beam structure under different gray-scale suppression operators
Table 3 灰度单元抑制算子 c /(N·mm)I /次G / %文献[18 ] 方法 205.83 64 27.23 本文方法 196.65 35 3.89
图5 不同拓扑优化方法下MBB梁结构的拓扑构型 Topology configuration of MBB beam structure under different topology optimization methods Fig.5 ![]()
表4 不同拓扑优化方法下MBB 梁结构优化结果 ...
... Optimization results of MBB beam structure under different topology optimization methods
Table 4 拓扑优化方法 c /(N·mm)I /次G / %文献[14 ] 方法 211.76 72 24.27 文献[18 ] 方法 213.00 83 26.66 本文方法 196.65 35 3.89
由图2 可知,采用不同SIMP法得到的MBB梁结构的拓扑构型存在一定差别:采用本文SIMP法所得优化结构在“区域A”和“区域B”处的颜色比较浅.意味着采用文献[18 ]的SIMP法所得结构在这些区域的xi 值比较大. ...
... 由图2 可知,采用不同SIMP法得到的MBB梁结构的拓扑构型存在一定差别:采用本文SIMP法所得优化结构在“区域A”和“区域B”处的颜色比较浅.意味着采用文献[18 ]的SIMP法所得结构在这些区域的xi 值比较大. ...
... 由表1 可知,相比于文献[18 ]的SIMP法,采用本文SIMP法所得结构的c 和I 分别减少了2.09 %和22.89 %,G 增加了13.80 %.说明采用本文SIMP法可以使拓扑优化方法的收敛速度变得更快,而且可以增强获取柔顺度小的优化结构的能力,但也会削弱抑制灰度单元产生的能力. ...
... 由表2 可知,相比于文献[18 ]的敏度过滤法,采用本文敏度过滤法所得的c 和G 分别减少了1.30 %和10.25 %,I 则相同.说明采用本文的敏度过滤法可以增强获取柔顺度小的优化结构和抑制灰度单元产生的能力. ...
... 由表3 可知,相比于文献[18 ]的灰度单元抑制算子,采用本文灰度单元抑制算子所得的c 、G 和I 分别减少了4.46 %、85.71 %和45.31 %.说明采用本文灰度单元抑制算子不仅可以增强获取柔顺度小的优化结构和抑制灰度单元产生的能力,而且可以提高收敛速度. ...
... 由图5 可知,采用文献[14 ,18 ]方法所得结构的拓扑构型存在更多的灰度单元,而采用本文方法所得结构的轮廓非常清晰,且从宏观上看不到灰度单元,所得结构比较好. ...
... 由表4 可知,相比于文献[14 ,18 ],采用本文拓扑优化方法所得的c 分别减少了7.14 %和7.68 %,G 分别减少了83.97 %和85.41 %,I 分别减少了51.39 %和57.83 %.说明本文方法在收敛速度、获取柔顺度小和拓扑构型好的优化结构及抑制灰度单元产生等方面具有明显优势. ...
... Optimization results of cantilever beam structure under different topology optimization methods
Table 5 拓扑优化方法 c /(N·mm)I /次G / %文献[14 ] 方法 114.90 118 24.77 文献[18 ] 方法 115.23 167 26.75 本文方法 104.98 38 4.73
由图7 可知,采用不同拓扑优化方法所得结构的拓扑构型存在较明显差异.相比于文献[14 ,18 ],采用本文方法所得结构的轮廓更为清晰.由表5 可知,采用本文方法所得的c 、I 和G 分别为104.98 N·mm、38次和4.73 %,较文献[14 ]方法分别减少了8.63 %、67.80 %和80.90 %,较文献[18 ]方法分别减少了8.90 %、77.25 %和82.32 %.说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快. ...
... 由图7 可知,采用不同拓扑优化方法所得结构的拓扑构型存在较明显差异.相比于文献[14 ,18 ],采用本文方法所得结构的轮廓更为清晰.由表5 可知,采用本文方法所得的c 、I 和G 分别为104.98 N·mm、38次和4.73 %,较文献[14 ]方法分别减少了8.63 %、67.80 %和80.90 %,较文献[18 ]方法分别减少了8.90 %、77.25 %和82.32 %.说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快. ...
... ]方法分别减少了8.63 %、67.80 %和80.90 %,较文献[18 ]方法分别减少了8.90 %、77.25 %和82.32 %.说明本文方法在获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等方面均明显优于文献[14 ,18 ]的方法,此外,收敛速度也明显加快. ...
... ,18 ]的方法,此外,收敛速度也明显加快. ...
Integrated layout design of supports and structures
2
2010
... 通过进一步研究发现,式(1) 所示的插值方案对改进拓扑优化方法的效果有限.为了进一步改进拓扑优化方法,采用Zhu等[19 ] 提出的改进SIMP法构建材料插值方案,具体如下: ...
... 须说明的是,文献[19 ]提出的材料插值方案中并没有E min ,此处加入E min 的目的是规避奇异刚度矩阵的产生和确保xi 的下限值可取0. ...
Improved proportional topology optimization algorithm for minimum volume problem with stress constraints
2
2021
... 为了有效解决采用变密度法引起的数值不稳定问题及改进拓扑优化方法,选用合适的敏度过滤法极为重要.Cheng等[20 ] 发现,采用基于高斯权重函数的敏度过滤法可以得到轻质且拓扑构型好的优化结构.本文采用基于高斯权重函数的敏度过滤法对c 关于xi 的灵敏度进行过滤处理,具体为: ...
... 根据文献[10 ,20 ],得到: ...
A simple heuristic for gray-scale suppression in optimality criterion-based topology optimization
1
2009
... 为了有效抑制灰度单元的产生,受Groenwold等[21 ] 提出的灰度单元抑制算子的启发,在原灰度单元抑制算子(xi Bi η ) q 的基础上,引入正值比例参数β ,设计了一种新的灰度单元抑制算子,为: ...