CN105956292A - 一种进化水平集结构拓扑优化方法 - Google Patents
一种进化水平集结构拓扑优化方法 Download PDFInfo
- Publication number
- CN105956292A CN105956292A CN201610300744.4A CN201610300744A CN105956292A CN 105956292 A CN105956292 A CN 105956292A CN 201610300744 A CN201610300744 A CN 201610300744A CN 105956292 A CN105956292 A CN 105956292A
- Authority
- CN
- China
- Prior art keywords
- optimization
- topology
- level set
- topology optimization
- optimization method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Geometry (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Civil Engineering (AREA)
- Evolutionary Computation (AREA)
- Structural Engineering (AREA)
- Architecture (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种进化水平集结构拓扑优化方法,其特征在于:所述的一种进化水平集结构拓扑优化方法综合了进化结构优化算法和原来结构拓扑优化中所用的水平集方法的优点。所提算法在优化迭代过程中,能更好地选择相邻区域的节点并自由地增加新孔,从而降低变形能。该方法特别适用于那些预先无法准确确定孔的数目和位置的复杂结构的优化中,提高拓扑优化方法在寻优过程中的搜索能力。
Description
技术领域
本发明涉及拓扑优化技术领域,特别是涉及一种进化水平集结构拓扑优化方法。
背景技术
进化算法是一种智能计算方法,它的思想来源于达尔文《进化论》中的进化与自然选择思想,现已广泛应用于各种工程领域,如计算机科学、人工智能、复杂系统、机械设计等。
拓扑优化通常是指层次优化或一般的形状优化,它可以让工程师在设计之初,得到优化了的结构拓扑或新的结构拓扑,近而完成结构尺寸和形状的设计。但是拓扑优化在使用过程中也存在着缺陷,之一就是其数学模型推导复杂,求解困难。
结构拓扑优化领域使用了许多基于生物思维的方法,进化算法就是其中的一种流行且稳健的方法。这些基于人类思维的算法的主要特点就是可并行计算并能够实现全局优化搜索。
为了得到较好的最终的结构拓扑,Osher和Senthia提出了水平集的概念,结果证明作为数值化的工具,水平集方法是非常有效的。由于水平集方法能够得到边界平滑的最终拓扑且没有中间密度等诸多优点,它吸引了很多工程人员和研究者的广泛兴趣。
在使用水平集方法进行拓扑优化过程中,仅仅通过孔的移动或合并就可以实现拓扑优化,而不能产生新的孔。对一些工程问题来说,基于水平集的拓扑优化方法存在明显的缺点。要想克服这些缺点,一种方法就是在初始化假想的设计域时包含尽可能多的孔,以得到尽可能多的拓扑形式。为使优化结果良好,在优化之前初始假想的结构形式时应遵循下面的两条基本原则:
(1)孔的数量足够多,以包含所有可能的拓扑;
(2)孔的布置应该合理。
使用本发明这种方法,可以在优化过程中在不同的位置插入新孔,以改变最终的结构拓扑。
发明内容
为了克服上述现有技术的不足,本发明提供了一种进化水平集结构拓扑优化方法。
本发明所采用的技术方案是:本发明提出了一种改进了的LSM算法,该方法将ESO方法中可任意插入孔的技术引入到LSM方法中,克服了LSM方法的缺点。
本发明可以在优化过程中在不同的位置插入新孔,以改变最终的结构拓扑。从原始结构的方面来看,该方法扩大了LSM的搜索空间。
附图说明
图1第i个节点应变能的邻域
图2悬臂梁
图3(a)迭代步为6时的优化拓扑
图3(b)迭代步为16时的优化拓扑
图3(c)迭代步为27时的优化拓扑
图3(d)迭代步为38时的优化拓扑
图3(e)迭代步为72时的优化拓扑
图3悬臂梁迭代过程
图4 MBB梁
图5(a)迭代步为11时的优化拓扑
图5(b)迭代步为22时的优化拓扑
图5(c)迭代步为33时的优化拓扑
图5(d)迭代步为44时的优化拓扑
图5(e)迭代步为55时的优化拓扑
图5(f)迭代步为68时的优化拓扑
图5 MBB梁迭代过程
具体实施方式
下面结合附图对本发明进一步说明。
在水平集方法中,结构的边界被定义为零水平集,它能够表示复杂表面形状,如在表面上形成孔,把孔分开或者合成以改变表面边界。零水平集由目标函数如变形能、应力、或奇异值来决定,通过结构外边界的拼合和移动得到最优的结构形式。
为了克服拓扑结果依赖初始结构的问题,需要新的法则,使其能够在优化到一定迭代步时,在合适的位置插入新的孔。
本发明建立在单元应变能的基础上,如图1所示,计算公式如下:
本发明的方法的实现步骤如下:
(1)初始化。初始化结构的假想拓扑,根据距离函数表示外部边界;
(2)求解结构的平衡方程。使用有限元分析法计算被移去的区域,并通过线性适应系统连接移去的区域;
(3)计算所选节点的适应度。该值就是节点邻域αi的应变能;
(4)插孔。在有材料的区域内,根据上面计算的数值,移去应变能较低的单元。一般来说移去单元的比例为影响体积约束的单元的2%-3%;
(5)结构拓扑的进化。求解水平集方程以更新嵌入里面的函数。这一点与水平集算法相同;
(6)一致性验证。如果满足体积约束,则迭代结束,否则,重复(2)-(6)步,直到满足。
在第(3)步和第(4)步中,研究人员可以任意控制孔的位置,此外,孔的数量和迭代步数也可以分别通过编程实现。对于不同的优化问题,相应的参数应该适当地进行调整。
悬臂梁问题是拓扑优化时的常用问题。我们也以此为例。如图2所示,一根长64mm,高40mm的悬臂梁,在自由端的中部集中承载80N。该问题的目标函数是在材料体积的约束下,结构的应变能最小。示例中材料的弹性模数和泊松比分别为200GPa和0.3。δ(x)和H(x)的近似值分别用参数α=10-9和Δ=1.0来表示。体积比限制在25%。在64x 40平面上,使用4节点等参数单元划分网格,将该问题作为平面应力问题来处理。要想得到初始化时没有孔的拓扑的优化结果,传统的LSM算法不容易做到这一点。根据进化的结构优化方法理论,本文提出了在优化过程中可以任意插入孔的方法。此结构拓扑的进化拓扑表示在图3(a)至图3(c)中。可以得出所提方法是有效的。
为了进一步说明本发明的效率,应用本发明提出的一种进化水平集结构拓扑优化方法。对图4所示MBB梁进行了拓扑优化。此结构拓扑的进化拓扑表示在图5(a)至图5(f)中。
Claims (3)
1.一种进化水平集结构拓扑优化方法,其特征在于:所述的一种进化水平集结构拓扑优化方法综合了进化结构优化算法和原来结构拓扑优化中所用的水平集方法的优点。本发明可以在优化过程中在不同的位置插入新孔,以改变最终的结构拓扑。从原始结构的方面来看,该方法扩大了拓扑优化的搜索空间。
2.根据权利要求1所述的一种进化水平集结构拓扑优化方法,其特征在于:所述的优点是在水平集方法基础上,结合了进化结构优化算法,克服了原有的方法在结构优化过程中不能产生新孔及结果拓扑极大地依赖原始结构拓扑的缺点。
3.根据权利要求1和2所述的一种进化水平集结构拓扑优化方法,其特征在于:所述算法在优化迭代过程中,能更好地选择相邻区域的节点并自由地增加新孔,从而降低变形能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610300744.4A CN105956292A (zh) | 2016-05-05 | 2016-05-05 | 一种进化水平集结构拓扑优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610300744.4A CN105956292A (zh) | 2016-05-05 | 2016-05-05 | 一种进化水平集结构拓扑优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105956292A true CN105956292A (zh) | 2016-09-21 |
Family
ID=56914562
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610300744.4A Pending CN105956292A (zh) | 2016-05-05 | 2016-05-05 | 一种进化水平集结构拓扑优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105956292A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106570258A (zh) * | 2016-11-03 | 2017-04-19 | 华中科技大学 | 一种基于水平集拓扑优化的局部模态识别方法 |
CN106600696A (zh) * | 2016-12-01 | 2017-04-26 | 华中科技大学 | 一种结构及其多类型边界的水平集隐式表达方法 |
CN107025340A (zh) * | 2017-03-30 | 2017-08-08 | 华中科技大学 | 一种适用于增材制造的自支撑网状结构拓扑优化设计方法 |
CN107315872A (zh) * | 2017-06-23 | 2017-11-03 | 华中科技大学 | 一种高效的结构频率响应拓扑优化方法 |
CN109190233A (zh) * | 2018-08-24 | 2019-01-11 | 华南理工大学 | 一种结构拓扑优化方法 |
-
2016
- 2016-05-05 CN CN201610300744.4A patent/CN105956292A/zh active Pending
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106570258A (zh) * | 2016-11-03 | 2017-04-19 | 华中科技大学 | 一种基于水平集拓扑优化的局部模态识别方法 |
CN106570258B (zh) * | 2016-11-03 | 2019-06-28 | 华中科技大学 | 一种基于水平集拓扑优化的局部模态识别方法 |
CN106600696A (zh) * | 2016-12-01 | 2017-04-26 | 华中科技大学 | 一种结构及其多类型边界的水平集隐式表达方法 |
CN107025340A (zh) * | 2017-03-30 | 2017-08-08 | 华中科技大学 | 一种适用于增材制造的自支撑网状结构拓扑优化设计方法 |
CN107315872A (zh) * | 2017-06-23 | 2017-11-03 | 华中科技大学 | 一种高效的结构频率响应拓扑优化方法 |
CN109190233A (zh) * | 2018-08-24 | 2019-01-11 | 华南理工大学 | 一种结构拓扑优化方法 |
CN109190233B (zh) * | 2018-08-24 | 2020-11-24 | 华南理工大学 | 一种结构拓扑优化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105956292A (zh) | 一种进化水平集结构拓扑优化方法 | |
Li et al. | Stiffness design of machine tool structures by a biologically inspired topology optimization method | |
He | Lagrange crisis and generalized variational principle for 3D unsteady flow | |
CN100583164C (zh) | 一种树状形体的立体分解和分级骨架提取方法 | |
CN102262699A (zh) | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 | |
Liu et al. | Topology optimization of stiffened plate/shell structures based on adaptive morphogenesis algorithm | |
Alejandrino et al. | A machine learning approach of lattice infill pattern for increasing material efficiency in additive manufacturing processes | |
Shobeiri | The topology optimization design for cracked structures | |
CN101783016A (zh) | 一种基于形状分析的树冠外形提取方法 | |
Ahlquist et al. | Development of a digital framework for the computation of complex material and morphological behavior of biological and technological systems | |
CN109408939A (zh) | 一种兼顾应力和位移约束的薄板结构加强筋分布优化的改进方法 | |
CN104036550A (zh) | 基于形状语义的建筑立面激光雷达点云解译与重建的方法 | |
CN104463972A (zh) | 观音手状文物骨架线提取方法 | |
CN104484527B (zh) | 一种离散结构拓扑优化过程中均布载荷自动动态修改方法 | |
CN104462414A (zh) | 一种基于拓扑结构的流程图相似性方法 | |
CN106445882A (zh) | 一种由vof函数快速构造符号距离函数的改进clsvof方法 | |
Yalçın et al. | High-order detached eddy simulation of unsteady flow around NREL S826 airfoil | |
Zheng et al. | Algorithm composition of Chinese folk music based on swarm intelligence | |
CN105513136A (zh) | 一种基于水平集中心聚类的3d角色模型骨架提取方法 | |
Chen et al. | Cellular automata based level set method for image segmentation | |
Bermudez et al. | Supracellular measurement of spatially varying mechanical heterogeneities in live monolayers | |
CN104504283A (zh) | 一种确定随机分布碳纳米管网络刚度阈值的方法 | |
Li et al. | An intelligent computational approach for design optimization of stiffened engineering structures | |
Che | Bim-based artificial engineering integration method for building engineering database | |
Lukashevich | Computational solution of contact problems with unilateral constraints and friction by step-by-step analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Jia Haipeng Inventor after: Ma Wenrui Inventor after: Jiang Chundong Inventor after: Liu Congyi Inventor before: Jia Haipeng Inventor before: Liu Congyi |
|
CB03 | Change of inventor or designer information | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160921 |
|
RJ01 | Rejection of invention patent application after publication |