CN107220413A - 基于梯度微结构的材料/结构一体化设计方法 - Google Patents
基于梯度微结构的材料/结构一体化设计方法 Download PDFInfo
- Publication number
- CN107220413A CN107220413A CN201710308664.8A CN201710308664A CN107220413A CN 107220413 A CN107220413 A CN 107220413A CN 201710308664 A CN201710308664 A CN 201710308664A CN 107220413 A CN107220413 A CN 107220413A
- Authority
- CN
- China
- Prior art keywords
- design
- mrow
- gradient
- structural
- micro
- 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/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Micromachines (AREA)
Abstract
本发明公开了一种基于梯度微结构的材料/结构一体化设计方法,用于解决现有材料/结构一体化设计方法设计的结构刚度差的技术问题。技术方案是首先给定了梯度微结构的拓扑形式,并采用均匀化方法预测微结构的材料等效性能。引入控制微结构形状的梯度参数和结构设计域上的宏观伪密度两组独立的设计变量,通过SIMP插值模型建立设计变量与材料弹性常数之间的联系。在优化的过程中引入材料用量约束,以结构整体刚度为目标函数,进行结构拓扑优化得到设计结果。优化设计结果显示,在相同的40%材料使用用量情况下,结构柔顺度由背景技术的2296.8J降低到1786.6J,降低了22.2%,大幅提高了结构的刚度。
Description
技术领域
本发明涉及一种材料/结构一体化设计方法,特别涉及一种基于梯度微结构的材料/结构一体化设计方法。
背景技术
文献“刘岭,阎军,程耿东.考虑均一微结构的结构/材料两级协同优化.计算力学学报,2008,25(1):29-34.”中提出了考虑均一微结构的结构/材料两级协同优化方法。该方法在两个尺度上独立定义了单元密度作为设计变量,分别引入了SIMP(Solid IsotropicMaterial with Penalization)和PAMP(Porous Anisotropic Material withPenalization)方法对密度进行惩罚,并且采用了周长约束控制微结构拓扑的复杂度,借助均匀化方法建立了结构和材料之间的联系,将宏/微观两个尺度上的设计纳入到一个优化模型中,以宏观结构最大刚度为目标,实现了结构与材料的协同优化设计。
文献假设了材料微观结构的宏观均一性,将两个尺度上的设计集成到一个优化问题中进行求解。这种方法虽然考虑到了两个尺度之间的相互影响,但由于非常强调材料微结构的单一性,在一定程度上削弱了优化解的优越性。
随着科学技术的快速发展,传统基于均一微结构材料的结构设计已经远不能满足工程应用的需要,消除结构设计和材料设计之间的隔阂,充分发挥材料潜力的一体化设计逐渐成为了研究的热点。为了达到结构宏观性能最优与轻质化要求,理想方案应在不同部位使用大小、形状不同的多孔微结构形式,以便最大限度挖掘材料的潜力。
发明内容
为了克服现有材料/结构一体化设计方法设计的结构刚度差的不足,本发明提供一种基于梯度微结构的材料/结构一体化设计方法。该方法首先给定了梯度微结构的拓扑形式,并采用均匀化方法预测微结构的材料等效性能。为了实现材料与结构的协同设计,引入控制微结构形状的梯度参数和结构设计域上的宏观伪密度两组独立的设计变量,通过SIMP插值模型建立设计变量与材料弹性常数之间的联系。在优化的过程中引入材料用量约束,以结构整体刚度为目标函数,进行结构拓扑优化得到设计结果。该方法将宏观结构布局和微结构设计纳入到一个优化模型中,考虑到两个尺度之间的相互影响,实现了材料与结构最佳匹配和一体化设计。
本发明解决其技术问题所采用的技术方案:一种基于梯度微结构的材料/结构一体化设计方法,其特点是包括以下步骤:
步骤一、给定梯度材料微结构的拓扑形式,采用均匀化方法预测材料的等效性能DH,并建立材料微结构的梯度参数θ与等效性能DH之间的函数关系:DH(θ)=f(θ)×DB,其中f(θ)为多项式插值函数。
步骤二、建立拓扑优化模型,并对模型施加约束和边界载荷。
步骤三、定义悬臂梁为拓扑优化的设计域Ω,并将设计域Ω离散为n个有限单元。xi为单元对应的伪密度,pi是单元对应的材料微结构的梯度参数,Ui为单元位移向量,Ki为单元刚度矩阵,F为载荷向量,U为整体位移向量,K为结构总刚度矩阵,C为结构柔顺度函数。
步骤四、定义优化问题。优化目标函数为结构柔顺度函数最小,约束条件为材料使用量小于
s.t.KU=F
0<xi≤1,i=1,,…n
pmin<pi≤pmax,i=1,,…n
步骤五、有限元分析计算结构的位移响应U。根据U计算结构柔顺度函数对于设计域内单元的伪密度xi和材料微结构的梯度参数pi的灵敏度。
步骤六、根据上述求得的灵敏度进行优化,选取梯度优化算法,优化迭代得到设计结果。
本发明的有益效果是:首先给定了梯度微结构的拓扑形式,并采用均匀化方法预测微结构的材料等效性能。为了实现材料与结构的协同设计,引入控制微结构形状的梯度参数和结构设计域上的宏观伪密度两组独立的设计变量,通过SIMP插值模型建立设计变量与材料弹性常数之间的联系。在优化的过程中引入材料用量约束,以结构整体刚度为目标函数,进行结构拓扑优化得到设计结果。该方法将宏观结构布局和微结构设计纳入到一个优化模型中,考虑到两个尺度之间的相互影响,实现了材料与结构最佳匹配和一体化设计。优化设计结果在相同的40%材料使用用量情况下,采用背景技术方法,即使用均一微结构的情况下,结构柔顺度函数为2296.8J,而本发明方法的结构柔顺度降低到1786.6J,降低了22.2%,大幅提高了结构的刚度。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明方法中悬臂梁结构及尺寸的示意图。
图2是本发明方法中梯度微结构单胞的示意图。
图3是背景技术方法悬臂梁结构拓扑优化设计结果。
图4是本发明方法悬臂梁结构拓扑优化设计结果。
具体实施方式
参照图1-4。本发明基于梯度微结构的材料/结构一体化设计方法具体步骤如下:
步骤一、给定梯度材料微结构的拓扑形式,采用均匀化方法预测材料的等效性能DH,并建立材料微结构的梯度参数θ与等效性能DH之间的函数关系:DH(θ)=f(θ)×DB,f(θ)为多项式插值函数。
步骤二、建立拓扑优化模型,悬臂梁长度800mm,高度300mm。悬臂梁左端固定,右下角施加集中载荷F=30N,方向竖直向下。
步骤三、定义悬臂梁为拓扑优化的设计域Ω,并将设计域Ω离散为600个有限单元。假设微结构基体材料为均质材料,杨氏模量E=1000,泊松比μ=0.3。xi为单元对应的伪密度,pi是单元对应的材料微结构的梯度参数,Ui为单元位移向量,Ki为单元刚度矩阵,F为载荷向量,U为整体位移向量,K为结构总刚度矩阵,C为结构柔顺度函数。
步骤四、定义优化问题。优化目标函数为结构柔顺度函数最小,约束条件为材料使用量小于40%:
s.t.KU=F
V(X)-0.4≤0
0<xi≤1,i=1,,…600
0°<pi≤45°,i=1,,…600
步骤五、用有限元分析软件Ansys计算结构模型的位移响应U。根据U计算结构柔顺度函数对于设计域内单元的伪密度xi和材料微结构的梯度参数pi的灵敏度。
步骤六、根据上述求得的灵敏度,选取梯度优化算法GCMMA(Globally ConvergentMethod of Moving Asymptotes)进行优化迭代,最终取得优化结果。
采用基于梯度微结构的材料/结构一体化设计方法能够有效的降低结构柔顺度,实现材料和结构的最优化匹配。优化结果对比如表1所示。从优化结果(图3、图4)对比可以看出,图3(背景技术方法)中悬臂梁桁架结构较少,刚度较弱,结构柔顺度函数为2296.8J;图4(本发明方法)中桁架结构分布均匀,不同部位微结构构型分布较为合理,结构柔顺度函数降低到1786.6J,相对于基于单一微结构的设计结果,大幅提高了结构的刚度。
表1
本发明所采用的方法很好地解决了宏观结构布局与梯度微结构的一体化设计。与传统基于单一微结构的优化结果对比下,本发明方法的优化结果性能更好。
Claims (1)
1.一种基于梯度微结构的材料/结构一体化设计方法,其特征在于包括以下步骤:
步骤一、给定梯度材料微结构的拓扑形式,采用均匀化方法预测材料的等效性能DH,并建立材料微结构的梯度参数θ与等效性能DH之间的函数关系:DH(θ)=f(θ)×DB,其中f(θ)为多项式插值函数;
步骤二、建立拓扑优化模型,并对模型施加约束和边界载荷;
步骤三、定义悬臂梁为拓扑优化的设计域Ω,并将设计域Ω离散为n个有限单元;xi为单元对应的伪密度,pi是单元对应的材料微结构的梯度参数,Ui为单元位移向量,Ki为单元刚度矩阵,F为载荷向量,U为整体位移向量,K为结构总刚度矩阵,C为结构柔顺度函数;
步骤四、定义优化问题;优化目标函数为结构柔顺度函数最小,约束条件为材料使用量小于
find X=(x1,x2,,…xn);P=(p1,p2,,…pn)
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>,</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
s.t.KU=F
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>V</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
步骤五、有限元分析计算结构的位移响应U;根据U计算结构柔顺度函数对于设计域内单元的伪密度xi和材料微结构的梯度参数pi的灵敏度;
步骤六、根据上述求得的灵敏度进行优化,选取梯度优化算法,优化迭代得到设计结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710308664.8A CN107220413A (zh) | 2017-05-04 | 2017-05-04 | 基于梯度微结构的材料/结构一体化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710308664.8A CN107220413A (zh) | 2017-05-04 | 2017-05-04 | 基于梯度微结构的材料/结构一体化设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107220413A true CN107220413A (zh) | 2017-09-29 |
Family
ID=59943675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710308664.8A Pending CN107220413A (zh) | 2017-05-04 | 2017-05-04 | 基于梯度微结构的材料/结构一体化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107220413A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563010A (zh) * | 2017-08-08 | 2018-01-09 | 西北工业大学 | 基于形状特征的多尺度结构材料一体化设计方法 |
CN107885929A (zh) * | 2017-11-07 | 2018-04-06 | 西北工业大学 | 非均匀类蜂窝结构设计方法 |
CN108491651A (zh) * | 2018-03-29 | 2018-09-04 | 福建工程学院 | 一种孔隙材料结构的设计方法 |
CN108595808A (zh) * | 2018-04-16 | 2018-09-28 | 西北工业大学 | 一种含椭圆形孔洞材料的结构宏微观协同优化方法 |
CN109871574A (zh) * | 2018-12-28 | 2019-06-11 | 华中科技大学 | 一种基于代理模型的多尺度拓扑优化方法 |
CN110717208A (zh) * | 2019-09-20 | 2020-01-21 | 华中科技大学 | 一种基于连续梯度微结构的多尺度频率响应拓扑优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412987A (zh) * | 2013-07-29 | 2013-11-27 | 西北工业大学 | 基于ks函数的多组件结构系统布局优化设计方法 |
CN103425831A (zh) * | 2013-08-06 | 2013-12-04 | 西北工业大学 | 基于多弧段曲线的压力载荷下的结构拓扑-形状联合优化方法 |
CN103440378A (zh) * | 2013-08-27 | 2013-12-11 | 西北工业大学 | 基于应力约束的机翼翼梁结构拓扑优化方法 |
CN104268317A (zh) * | 2014-09-12 | 2015-01-07 | 西北工业大学 | 机械零件圆角结构形状优化方法 |
-
2017
- 2017-05-04 CN CN201710308664.8A patent/CN107220413A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412987A (zh) * | 2013-07-29 | 2013-11-27 | 西北工业大学 | 基于ks函数的多组件结构系统布局优化设计方法 |
CN103425831A (zh) * | 2013-08-06 | 2013-12-04 | 西北工业大学 | 基于多弧段曲线的压力载荷下的结构拓扑-形状联合优化方法 |
CN103440378A (zh) * | 2013-08-27 | 2013-12-11 | 西北工业大学 | 基于应力约束的机翼翼梁结构拓扑优化方法 |
CN104268317A (zh) * | 2014-09-12 | 2015-01-07 | 西北工业大学 | 机械零件圆角结构形状优化方法 |
Non-Patent Citations (1)
Title |
---|
刘岭等: "考虑均一微结构的结构/材料两级协同优化", 《计算力学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563010A (zh) * | 2017-08-08 | 2018-01-09 | 西北工业大学 | 基于形状特征的多尺度结构材料一体化设计方法 |
CN107563010B (zh) * | 2017-08-08 | 2020-09-25 | 西北工业大学 | 基于形状特征的多尺度结构材料一体化设计方法 |
CN107885929A (zh) * | 2017-11-07 | 2018-04-06 | 西北工业大学 | 非均匀类蜂窝结构设计方法 |
CN108491651A (zh) * | 2018-03-29 | 2018-09-04 | 福建工程学院 | 一种孔隙材料结构的设计方法 |
CN108595808A (zh) * | 2018-04-16 | 2018-09-28 | 西北工业大学 | 一种含椭圆形孔洞材料的结构宏微观协同优化方法 |
CN109871574A (zh) * | 2018-12-28 | 2019-06-11 | 华中科技大学 | 一种基于代理模型的多尺度拓扑优化方法 |
CN110717208A (zh) * | 2019-09-20 | 2020-01-21 | 华中科技大学 | 一种基于连续梯度微结构的多尺度频率响应拓扑优化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107220413A (zh) | 基于梯度微结构的材料/结构一体化设计方法 | |
Hasanien | Design optimization of PID controller in automatic voltage regulator system using Taguchi combined genetic algorithm method | |
Yu | Controlling lorenz chaos | |
CN107145657B (zh) | 针对反应堆中子扩散方程的非均匀几何变分节块方法 | |
CN107563010B (zh) | 基于形状特征的多尺度结构材料一体化设计方法 | |
CN103577654B (zh) | 一种大型汽轮发电机定子线棒的有限元精确建模方法 | |
Thunberg et al. | A lifting method for analyzing distributed synchronization on the unit sphere | |
CN105313336A (zh) | 一种薄壳体3d打印优化方法 | |
Zhou et al. | Topology optimization of easy-removal support structures for additive manufacturing | |
Liu et al. | Generating support structures for additive manufacturing with continuum topology optimization methods | |
CN109657378A (zh) | 一种含变尺寸单胞的非均质层级结构拓扑优化方法 | |
CN110210151B (zh) | 基于b样条的点阵结构参数化隐式建模与优化方法 | |
CN111125963A (zh) | 基于拉格朗日积分点有限元的数值仿真系统及方法 | |
CN104268317B (zh) | 机械零件圆角结构形状优化方法 | |
CN108897962A (zh) | 一种基于非线性无网格法的柔性机构拓扑优化设计方法 | |
JP2011133930A (ja) | 形状最適化プログラム、方法及び装置 | |
CN103294861A (zh) | 多弧段曲线形状优化设计方法 | |
CN109815518A (zh) | 基于转动惯量约束的飞行器舵面设计方法 | |
CN105589996A (zh) | 一种以定转子最小对称单元为几何模型的电机有限元热分析方法 | |
Xie et al. | Suboptimal output consensus for a group of weakly nonminimum phase linear systems | |
CN108595808A (zh) | 一种含椭圆形孔洞材料的结构宏微观协同优化方法 | |
CN111127491B (zh) | 一种针对蜂窝结构拓扑优化的多变量水平分割方法及设备 | |
CN109783965B (zh) | 一种结构网格自动分块加密方法 | |
JP2017084081A (ja) | 磁性体シミュレーション装置、マイクロ磁化算出方法及びプログラム | |
Antonietti et al. | High order discontinuous Galerkin methods on surfaces |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170929 |
|
WD01 | Invention patent application deemed withdrawn after publication |