CN102789533B - 基于密度阈值的结构拓扑优化设计灵敏度过滤方法 - Google Patents
基于密度阈值的结构拓扑优化设计灵敏度过滤方法 Download PDFInfo
- Publication number
- CN102789533B CN102789533B CN201210268937.8A CN201210268937A CN102789533B CN 102789533 B CN102789533 B CN 102789533B CN 201210268937 A CN201210268937 A CN 201210268937A CN 102789533 B CN102789533 B CN 102789533B
- Authority
- CN
- China
- Prior art keywords
- density
- topology optimization
- sensitivity
- density threshold
- design
- 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.)
- Active
Links
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于密度阈值的结构拓扑优化设计灵敏度过滤方法,用于解决现有的结构拓扑优化设计灵敏度过滤方法中拓扑优化设计结构可制造性差的技术问题。技术方案是首先建立拓扑优化模型,再通过有限元计算,得到单元灵敏度和单元伪密度,计算模型单元伪密度xi值介于给定范围内的单元数n,再计算中间单元百分比C,若优化迭代第k步的百分比C小于密度阈值T,则在第k+1中禁用灵敏度过滤直到迭代结束,采用优化准则法进行优化设计,得到二维简支梁的拓扑优化结果,拓扑优化设计结构可制造性变强。由于采用密度阈值法,结构拓扑优化设计的所需的迭代步数由背景技术的112步收敛减少为41步收敛,减少了71步。
Description
技术领域
本发明涉及一种结构拓扑优化设计灵敏度过滤方法,特别是涉及一种基于密度阈值的结构拓扑优化设计灵敏度过滤方法。
背景技术
参照图1、图2。在航空航天、汽车制造等领域,为了满足结构轻量化设计需求,拓扑优化设计技术广泛应用于其零件结构设计中。由于连续结构进行有限元离散处理后数值解的非唯一性,不带附加约束的拓扑优化设计模型是一个病态问题,会出现棋盘格现象。导致最终拓扑优化设计结构可制造性降低,失去了拓扑优化设计技术的工程意义。
文献1“99 line topology optimization code written in Matlab.Sigmund,Structuraland Multidisciplinary Optimization,2001,vol.21,pp.120-127”公开了一种基于灵敏度过滤技术的结构拓扑优化设计方法。通过有限元计算,得到设计目标值和单元灵敏度,对单元灵敏度进行灵敏度过滤,避免了在优化迭代结果出现棋盘格现象。
文献2“Morphology-based black and white filters for topology optimization.Sigmund,Structural and Multidisciplinary Optimization,2007,vol.33,pp.401-424”公开了采用文献1中公开的灵敏度过滤方法优化简支梁的优化结果,该方法避免了优化迭代结果出现棋盘格现象,但优化迭代步数为112步,仍较多。
文献公开的方法虽然通过灵敏度过滤方法避免了拓扑优化结果出现棋盘格现象,但是该方法增加了优化迭代步数,降低了拓扑优化的计算效率。
发明内容
为了克服现有的结构拓扑优化设计灵敏度过滤方法中拓扑优化设计结构可制造性差的不足,本发明提供一种基于密度阈值的结构拓扑优化设计灵敏度过滤方法。该方法采用密度阈值法,通过计算迭代步中有限元模型单元伪密度数值介于给定范围内的单元数目百分比,记为C,给定密度阈值T;当迭代步中的密度值C小于密度阈值T时,则在当前迭代步之后禁用灵敏度过滤。该方法在避免拓扑优化结果出现棋盘格现象,拓扑优化设计结构可制造性变强。同时,可以减少拓扑优化迭代步数,使结构拓扑优化设计灵敏度过滤方法优化迭代过程变简单。
本发明解决其技术问题所采用的技术方案是:一种基于密度阈值的结构拓扑优化设计灵敏度过滤方法,其特点是包括以下步骤:
(a)建立拓扑优化模型,设计域为整个二维简支梁,定义拓扑优化的设计域、目标函数和约束条件,设定密度过滤下限ra,上限rb,密度阈值T。
(b)通过有限元计算,得到单元灵敏度和单元伪密度。
(c)计算模型单元伪密度xi值介于给定范围内的单元数n,即:
ra≤xi≤rb(i=1,2,…,N) (1)
计算中间单元百分比C:
式中,N是有限元单元总数。
(d)若优化迭代第k步的百分比C小于密度阈值T,则在第k+l中禁用灵敏度过滤直到迭代结束。
(e)采用优化准则法进行优化设计,得到二维简支梁的拓扑优化结果。
所述二维简支梁结构的体积分数小于0.5。
所述密度过滤下限ra=0.1。
所述密度过滤上限rb=0.8。
所述密度阈值T=0.6。
本发明的有益效果是:由于采用密度阈值法,通过计算迭代步中有限元模型单元伪密度数值介于给定范围内的单元数目百分比,记为C,给定密度阈值T;当迭代步中的密度值C小于密度阈值T时,则在当前迭代步之后禁用灵敏度过滤。该方法在避免拓扑优化结果出现棋盘格现象,拓扑优化设计结构可制造性变强。同时,减少了拓扑优化迭代步数,使结构拓扑优化设计灵敏度过滤方法优化迭代过程变简单。实施例中结构拓扑优化设计的所需的迭代步数由背景技术的112步收敛减少为41步收敛,减少了71步。
下面结合附图和实施例对本发明作详细说明。
附图说明
图1是背景技术方法的简支梁拓扑优化结果示意图。
图2是二维简支梁结构示意图。
具体实施方式
参照图2,以简支梁在对称压力载荷下的拓扑优化设计为例说明本发明。二维简支梁杨氏模量E=1GPa,泊松比μ=0.3,二维简支梁尺寸为长240mm,高40mm,载荷为在简支梁的上边中间部位施加集中力F=1kN,方向向下。方法步骤如下:
(a)建立拓扑优化模型,设计域为整个简支梁,定义拓扑优化的设计域、目标函数和约束条件,设计目标为二维简支梁结构整体的刚度最大,约束条件为结构体积分数小于0.5,设定密度过滤下限ra=0.1,密度过滤上限rb=0.8,密度阈值T=0.6。
(b)在Matlab软件中进行有限元计算,求得单元的灵敏度和单元伪密度。
(c)计算模型单元伪密度xi值介于给定范围内的单元数n,即:
ra≤xi≤rb(i=1,2,…,N) (1)
计算中间单元百分比C:
其中N有限元单元总数。
(d)判断优化迭代第k步的百分比C,若小于密度阈值,则在第k+1步中禁用灵敏度过滤直到迭代结束。
(e)采用优化准则法进行优化设计,得到二维简支梁拓扑优化结果。
采用基于密度阈值的结构拓扑优化设计灵敏度过滤法与文献2中的方法优化迭代步数的对比如表1所示。
表1
Claims (2)
1.一种基于密度阈值的结构拓扑优化设计灵敏度过滤方法,其特征在于包括以下步骤:
(a)建立拓扑优化模型,设计域为整个二维简支梁,定义拓扑优化的设计域、目标函数和约束条件,设定密度过滤下限ra=0.1,密度过滤上限rb=0.8,密度阈值T=0.6;
(b)通过有限元计算,得到单元灵敏度和单元伪密度;
(c)计算模型单元伪密度xi值介于给定范围内的单元数n,即:
ra≤xi≤rb(i=1,2,...,N) (1)
计算中间单元百分比C:
(d)若优化迭代第k步的百分比C小于密度阈值T,则在第k+1中禁用灵敏度过滤直到迭代结束;
(e)采用优化准则法进行优化设计,得到二维简支梁的拓扑优化结果。
2.根据权利要求1所述的基于密度阈值的结构拓扑优化设计灵敏度过滤方法,其特征在于:所述二维简支梁结构的体积分数小于0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210268937.8A CN102789533B (zh) | 2012-07-31 | 2012-07-31 | 基于密度阈值的结构拓扑优化设计灵敏度过滤方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210268937.8A CN102789533B (zh) | 2012-07-31 | 2012-07-31 | 基于密度阈值的结构拓扑优化设计灵敏度过滤方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102789533A CN102789533A (zh) | 2012-11-21 |
CN102789533B true CN102789533B (zh) | 2014-07-16 |
Family
ID=47154936
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210268937.8A Active CN102789533B (zh) | 2012-07-31 | 2012-07-31 | 基于密度阈值的结构拓扑优化设计灵敏度过滤方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102789533B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020361B (zh) * | 2012-12-17 | 2015-08-26 | 华南理工大学 | 一种柔顺机构无棋盘格拓扑图提取方法 |
CN103440378B (zh) * | 2013-08-27 | 2016-06-08 | 西北工业大学 | 基于应力约束的机翼翼梁结构拓扑优化方法 |
CN105092479B (zh) * | 2014-05-04 | 2021-12-21 | 睿励科学仪器(上海)有限公司 | 一种用于优化测量系统的灵敏度的方法 |
CN105653779B (zh) * | 2015-12-28 | 2019-02-12 | 大连理工大学 | 基于温度比拟连通性的可制造性约束拓扑优化方法 |
CN106372347B (zh) * | 2016-09-08 | 2019-05-07 | 厦门大学嘉庚学院 | 改进双向渐进法的等效静载荷法动态响应拓扑优化方法 |
US11914934B2 (en) | 2018-09-14 | 2024-02-27 | Siemens Industry Software Inc. | Active region adaptations for design domains in topology optimizations |
CN110069864B (zh) * | 2019-04-24 | 2023-07-21 | 南京航空航天大学 | 一种结合变密度法的改进双向渐进结构拓扑优化方法 |
CN110348102B (zh) * | 2019-07-04 | 2022-11-25 | 广州大学 | 基于反正切的动态进化率beso拓扑优化方法及其应用 |
CN113191077B (zh) * | 2021-04-25 | 2022-12-09 | 西安交通大学 | 一种基于连续纤维复材3d打印的变纤维含量拓扑优化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504678A (zh) * | 2009-03-09 | 2009-08-12 | 西北工业大学 | 薄壁加筋结构加筋布局优化设计方法 |
CN101814103A (zh) * | 2010-04-01 | 2010-08-25 | 西北工业大学 | 基于超单元的多组件布局建模与结构优化设计方法 |
CN101980222A (zh) * | 2010-11-03 | 2011-02-23 | 北京航空航天大学 | 一种基于拓扑优化的仿蜘蛛网夹层筋板复合结构设计方法 |
-
2012
- 2012-07-31 CN CN201210268937.8A patent/CN102789533B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504678A (zh) * | 2009-03-09 | 2009-08-12 | 西北工业大学 | 薄壁加筋结构加筋布局优化设计方法 |
CN101814103A (zh) * | 2010-04-01 | 2010-08-25 | 西北工业大学 | 基于超单元的多组件布局建模与结构优化设计方法 |
CN101980222A (zh) * | 2010-11-03 | 2011-02-23 | 北京航空航天大学 | 一种基于拓扑优化的仿蜘蛛网夹层筋板复合结构设计方法 |
Non-Patent Citations (4)
Title |
---|
Bi-directional evolutionary topology optimization using element replaceable method;J.H.Zhu等;《Comput Mech》;20071231;全文 * |
J.H.Zhu等.Bi-directional evolutionary topology optimization using element replaceable method.《Comput Mech》.2007,全文. |
Jihong Zhu等.Simultaneous design of components layout and supporting structures using coupled shape and topology optimization technique.《Struct Multidisc Optim》.2008,全文. |
Simultaneous design of components layout and supporting structures using coupled shape and topology optimization technique;Jihong Zhu等;《Struct Multidisc Optim》;20081231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN102789533A (zh) | 2012-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102789533B (zh) | 基于密度阈值的结构拓扑优化设计灵敏度过滤方法 | |
CN102663695B (zh) | 基于小波变换的dr图像去噪方法及系统 | |
CN102879824B (zh) | 一种基于迭代收缩的快速稀疏Radon变换方法 | |
CN103017802B (zh) | 基于小波变换的布里渊谱去噪方法 | |
CN103743435A (zh) | 一种多传感器数据融合方法 | |
CN105527649A (zh) | 一种多域多次分离的高效采集多震源混合数据分离方法 | |
CN105354362A (zh) | 模拟二维水流运动的三次样条多尺度有限元方法 | |
CN102819679B (zh) | 基于周长阈值的结构拓扑优化设计灵敏度过滤方法 | |
CN103941692B (zh) | 基于驻留概率信息方法的随机切换系统h∞滤波信息处理方法 | |
CN104700361A (zh) | 基于边缘检测的图像插值方法及系统 | |
CN105319593A (zh) | 基于曲波变换和奇异值分解的联合去噪方法 | |
CN104142908A (zh) | 电波传播抛物方程分步傅里叶变换解的上边界处理方法 | |
CN109800500B (zh) | 存储介质、尾矿库渗流稳定性的计算方法及其装置 | |
CN103869376A (zh) | 一种基于深度变化的三维重力位场正则化下延方法及应用 | |
CN103268595B (zh) | 基于侧抑制系数的红外图像增强处理方法 | |
Kang et al. | Parameter estimation of variable-parameter nonlinear Muskingum model using excel solver | |
CN105664726A (zh) | 净水机的冲洗方法及系统 | |
CN104657952A (zh) | 一种结合时空信息的遥感逐日积雪产品去云方法 | |
CN102104560A (zh) | 信道估计方法及装置 | |
CN106871891B (zh) | 一种基于累加累减混合微分器的信号去噪方法 | |
CN105160436A (zh) | 一种基于n邻近Lipschitz支撑面的广义增强群体全局优化方法 | |
CN105093221B (zh) | 一种小波域干涉合成孔径雷达相位滤波方法及装置 | |
CN103295233A (zh) | 合成孔径雷达影像配准中最优匹配窗口的选择方法 | |
CN105068374B (zh) | 光学临近修正中的二维图形快速识别方法 | |
CN102955895B (zh) | 一种基于面积的曲线亚像素配准方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |