CN107885929A - 非均匀类蜂窝结构设计方法 - Google Patents

非均匀类蜂窝结构设计方法 Download PDF

Info

Publication number
CN107885929A
CN107885929A CN201711081269.7A CN201711081269A CN107885929A CN 107885929 A CN107885929 A CN 107885929A CN 201711081269 A CN201711081269 A CN 201711081269A CN 107885929 A CN107885929 A CN 107885929A
Authority
CN
China
Prior art keywords
cell element
mrow
wall thickness
mtd
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.)
Pending
Application number
CN201711081269.7A
Other languages
English (en)
Inventor
王创
朱继宏
张卫红
李韶英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201711081269.7A priority Critical patent/CN107885929A/zh
Publication of CN107885929A publication Critical patent/CN107885929A/zh
Pending legal-status Critical Current

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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Vibration Dampers (AREA)

Abstract

本发明公开了一种非均匀类蜂窝结构设计方法，用于解决现有非均匀类蜂窝结构设计方法设计效率低的技术问题。技术方案是首先建立类蜂窝胞元的拓扑形式，选取胞元壁厚作为设计变量，采用均匀化理论预测类蜂窝胞元的材料等效性能，建立胞元壁厚与材料弹性常数之间的函数关系。然后在优化过程中引入材料用量约束，以结构整体刚度最大为目标，进行结构优化得到设计结果。该方法能够直接获得各蜂窝胞元的壁厚，避免了采用密度映射函数计算所有拓扑微单元对蜂窝胞元的影响，简化了计算流程，提高了设计效率，最终实现类蜂窝结构的非均匀设计。

Description

非均匀类蜂窝结构设计方法
技术领域
本发明涉及一种类蜂窝结构设计方法，特别涉及一种非均匀类蜂窝结构设计方法。
背景技术
类蜂窝夹层结构具有重量轻、比刚度大、比强度高和抗冲击等众多优点，被广泛应用于航空航天、船舶、汽车、桥梁建筑等领域。已有研究表明，实现蜂窝结构各胞元的壁厚非均匀变化，可使结构的材料分布更加合理，进而提升结构的承载性能。
文献1“Zhang P,Toman J,Yu Y,et al,Efficient Design-Optimization ofVariable-Density Hexagonal Cellular Structure by Additive Manufacturing:Theory and Validation.Journal of Manufacturing Science & Engineering,2015,137(2):021004.”提出一种基于拓扑优化密度映射的变密度蜂窝结构设计方法。该方法结合拓扑优化与增材制造技术设计了不同直径圆孔阵列构成的类蜂窝结构，实现了类似不均匀壁厚变化的蜂窝设计。
文献2“邢昊,敬石开,张贺,等.拓扑优化密度映射的非均匀蜂窝结构设计方法.计算机辅助设计与图形学学报,2017,29(4):734-741.”将拓扑微单元的密度信息映射到蜂窝胞元的壁厚上，并采用响应面与协同优化方法对蜂窝胞元尺寸参数和密度映射权重参数进行优化，获得了壁厚非均匀的蜂窝结构。
文献1和2都是基于拓扑优化密度信息，建立与蜂窝胞元尺寸参数之间的联系，以实现蜂窝结构的非均匀设计。这种方法首先需要进行拓扑优化设计，得到单元的密度信息，而单元密度与蜂窝胞元尺寸参数之间并不具有简单的对应关系，需要采用与距离相关的密度映射函数计算所有拓扑微单元对蜂窝胞元的影响。这一过程需要进行大量的数值计算，因此该方法的设计效率不高。
发明内容
为了克服现有非均匀类蜂窝结构设计方法设计效率低的不足，本发明提供一种非均匀类蜂窝结构设计方法。该方法首先建立类蜂窝胞元的拓扑形式，选取胞元壁厚作为设计变量，采用均匀化理论预测类蜂窝胞元的材料等效性能，建立胞元壁厚与材料弹性常数之间的函数关系。然后在优化过程中引入材料用量约束，以结构整体刚度最大为目标，进行结构优化得到设计结果。该方法能够直接获得各蜂窝胞元的壁厚，避免了采用密度映射函数计算所有拓扑微单元对蜂窝胞元的影响，可以简化计算流程，提高设计效率，最终实现类蜂窝结构的非均匀设计。
本发明解决其技术问题所采用的技术方案：一种非均匀类蜂窝结构设计方法，其特点是包括以下步骤：
(a)给定类蜂窝胞元的拓扑构型，选取胞元壁厚t作为设计变量。通过改变胞元壁厚，得到一系列壁厚不同的类蜂窝胞元。采用均匀化理论分别计算类蜂窝胞元的材料等效性能，建立胞元壁厚t与材料弹性常数之间的函数关系。
(b)建立拓扑优化模型，并对模型施加约束和边界载荷。
(c)定义拓扑优化的设计域Ω，将设计域Ω离散为n个有限单元。
(d)建立优化问题的数学模型。优化目标为结构整体柔顺度函数最小，约束条件为材料使用量小于
find T＝(t1，t2，...，tn)
式中，Ui为单元位移向量，Ki为单元刚度矩阵，F为载荷向量，U为整体位移向量，K为结构总刚度矩阵，C为结构柔顺度函数，ti是单元i对应类蜂窝胞元的壁厚，tmin和tmax分别指类蜂窝胞元壁厚的最小值和最大值。
(e)有限元分析计算结构的整体位移向量U。根据U计算结构柔顺度函数C对于单元i对应类蜂窝胞元壁厚ti的灵敏度。
(f)选取梯度优化算法，根据步骤(e)求得的灵敏度对各胞元壁厚进行优化，得到设计结果。根据得到的胞元壁厚信息，进行模型重构，获得最终的非均匀类蜂窝结构。
本发明的有益效果是：该方法首先建立类蜂窝胞元的拓扑形式，选取胞元壁厚作为设计变量，采用均匀化理论预测类蜂窝胞元的材料等效性能，建立胞元壁厚与材料弹性常数之间的函数关系。然后在优化过程中引入材料用量约束，以结构整体刚度最大为目标，进行结构优化得到设计结果。本发明直接将类蜂窝胞元壁厚作为设计变量，以结构整体刚度最大为设计目标，建立拓扑优化数学模型，通过灵敏度分析，求得目标函数对于设计变量的灵敏度，采用梯度优化算法进行优化设计，实用性强。在实施例中，基于本发明所提方法获得的设计结果，结构柔顺度函数从136.6J降低到90.3J，证明该方法能够显著提升类蜂窝结构的刚度。与背景文献中所提方法相比，避免了采用密度映射函数计算所有拓扑微单元对蜂窝胞元的影响，简化了计算流程，提高了设计效率，最终实现类蜂窝结构的非均匀设计。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明非均匀类蜂窝结构设计方法中类蜂窝胞元结构示意图。
图2是本发明方法实施例中的工况及尺寸示意图。
图3是基于本发明方法的非均匀类蜂窝结构设计结果图。
具体实施方式
参照图1-3。本发明非均匀类蜂窝结构设计方法，以MBB梁结构为拓扑优化对象，考虑结构整体柔顺度函数最小为例说明本发明。具体步骤如下：
(a)给定类蜂窝胞元的拓扑构型，选取胞元壁厚t作为设计变量。分别使胞元壁厚t＝0.5mm、0.6mm、0.7mm、0.8mm、0.9mm、1.0mm、1.1mm及1.2mm，得到8类壁厚不同的类蜂窝胞元。类蜂窝胞元材料的杨氏模量E＝1000，泊松比μ＝0.3。采用均匀化理论分别计算类蜂窝胞元的材料等效性能，建立胞元壁厚t与材料弹性常数之间的函数关系。
(b)建立拓扑优化模型。MBB梁模型长320mm，宽80mm，左下角点完全固定，右下角点约束y轴方向的自由度，结构上端面中部区域施加均布压力，载荷密度8N/mm，方向竖直向下。
(c)定义MBB梁为拓扑优化的设计域Ω，将设计域Ω离散为256个有限单元。
(d)建立优化问题的数学模型。优化目标为结构整体柔顺度函数最小，约束条件为材料使用量小于10240mm2：
find T＝(t1,t2,...,t256)
式中，Ui为单元位移向量，Ki为单元刚度矩阵，F为载荷向量，U为整体位移向量，K为结构总刚度矩阵，C为结构柔顺度函数，ti是单元i对应类蜂窝胞元的壁厚，0.5和1.2分别指类蜂窝胞元壁厚的最小值和最大值。
(e)使用有限元分析软件Ansys计算结构模型的整体位移向量U。根据U计算结构柔顺度函数C对于单元i对应类蜂窝胞元壁厚ti的灵敏度。
(f)根据上述求得的灵敏度，选取梯度优化算法GCMMA(Globally ConvergentMethod of Moving Asymptotes)，优化迭代得到设计结果。根据得到的胞元壁厚信息ti，进行模型重构，获得最终的非均匀类蜂窝结构。
表1
由图3优化设计结果可以看出，采用本发明方法，可以获得壁厚非均匀的类蜂窝结构，在优化迭代过程中，非均匀类蜂窝结构的柔顺度函数从136.6J降低到90.3J，结构刚度得到大幅度提升，材料分布更加合理。与背景技术中的方法相比，本发明方法直接以各蜂窝胞元的壁厚作为设计变量，避免了繁琐的密度映射过程，简化了计算流程，提升了结构的设计效率。因此，本发明所采用的方法很好地解决了非均匀类蜂窝结构设计问题。

Claims (1)

1.一种非均匀类蜂窝结构设计方法，其特征在于包括以下步骤：
(a)给定类蜂窝胞元的拓扑构型，选取胞元壁厚t作为设计变量；通过改变胞元壁厚，得到一系列壁厚不同的类蜂窝胞元；采用均匀化理论分别计算类蜂窝胞元的材料等效性能，建立胞元壁厚t与材料弹性常数之间的函数关系；
(b)建立拓扑优化模型，并对模型施加约束和边界载荷；
(c)定义拓扑优化的设计域Ω，将设计域Ω离散为n个有限单元；
(d)建立优化问题的数学模型；优化目标为结构整体柔顺度函数最小，约束条件为材料使用量小于
find T＝(t1,t2,...,tn)
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mi>min</mi> </mtd> <mtd> <mrow> <mi>C</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <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>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> </mrow> </mtd> <mtd> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>K</mi> <mi>U</mi> <mo>=</mo> <mi>F</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>V</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>-</mo> <mover> <mi>V</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;le;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>t</mi> <mi>min</mi> </msub> <mo>&lt;</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>t</mi> <mi>max</mi> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>n</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mtd> </mtr> </mtable> </mfenced>
式中，Ui为单元位移向量，Ki为单元刚度矩阵，F为载荷向量，U为整体位移向量，K为结构总刚度矩阵，C为结构柔顺度函数，ti是单元i对应类蜂窝胞元的壁厚，tmin和tmax分别指类蜂窝胞元壁厚的最小值和最大值；
(e)有限元分析计算结构的整体位移向量U；根据U计算结构柔顺度函数C对于单元i对应类蜂窝胞元壁厚ti的灵敏度；
(f)选取梯度优化算法，根据步骤(e)求得的灵敏度对各胞元壁厚进行优化，得到设计结果；根据得到的胞元壁厚信息，进行模型重构，获得最终的非均匀类蜂窝结构。
CN201711081269.7A 2017-11-07 2017-11-07 非均匀类蜂窝结构设计方法 Pending CN107885929A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711081269.7A CN107885929A (zh) 2017-11-07 2017-11-07 非均匀类蜂窝结构设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711081269.7A CN107885929A (zh) 2017-11-07 2017-11-07 非均匀类蜂窝结构设计方法

Publications (1)

Publication Number Publication Date
CN107885929A true CN107885929A (zh) 2018-04-06

Family

ID=61779091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711081269.7A Pending CN107885929A (zh) 2017-11-07 2017-11-07 非均匀类蜂窝结构设计方法

Country Status (1)

Country Link
CN (1) CN107885929A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109766656A (zh) * 2019-01-25 2019-05-17 北京航空航天大学 一种基于拓扑优化的梯度点阵结构设计方法
CN111125818A (zh) * 2019-12-12 2020-05-08 河南大学 基于仿生子结构拓扑优化的新型十字交叉节点及优化方法
CN111127491A (zh) * 2019-12-31 2020-05-08 华中科技大学 一种针对蜂窝结构拓扑优化的多变量水平分割方法及设备
CN111583243A (zh) * 2020-05-09 2020-08-25 中南大学 蜂窝规整度检测的相邻点胞元重构方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103425830A (zh) * 2013-08-06 2013-12-04 西北工业大学 基于多点位移协调约束的结构拓扑优化方法
CN106678221A (zh) * 2016-07-09 2017-05-17 李萌 一种提高结构强度的蜂窝结构及设计方法
CN107220413A (zh) * 2017-05-04 2017-09-29 西北工业大学 基于梯度微结构的材料/结构一体化设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103425830A (zh) * 2013-08-06 2013-12-04 西北工业大学 基于多点位移协调约束的结构拓扑优化方法
CN103425830B (zh) * 2013-08-06 2017-02-08 西北工业大学 基于多点位移协调约束的结构拓扑优化方法
CN106678221A (zh) * 2016-07-09 2017-05-17 李萌 一种提高结构强度的蜂窝结构及设计方法
CN107220413A (zh) * 2017-05-04 2017-09-29 西北工业大学 基于梯度微结构的材料/结构一体化设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
PU ZHANG等: "Efficient Design-Optimization", 《JOURNAL OF MANUFACTURING SCIENCE AND ENGINEERING》 *
邢昊等: "拓扑优化密度映射的非均匀蜂窝结构设计方法", 《计算机辅助设计与图形学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109766656A (zh) * 2019-01-25 2019-05-17 北京航空航天大学 一种基于拓扑优化的梯度点阵结构设计方法
CN109766656B (zh) * 2019-01-25 2021-03-09 北京航空航天大学 一种基于拓扑优化的梯度点阵结构设计方法
CN111125818A (zh) * 2019-12-12 2020-05-08 河南大学 基于仿生子结构拓扑优化的新型十字交叉节点及优化方法
CN111127491A (zh) * 2019-12-31 2020-05-08 华中科技大学 一种针对蜂窝结构拓扑优化的多变量水平分割方法及设备
CN111127491B (zh) * 2019-12-31 2022-09-27 华中科技大学 一种针对蜂窝结构拓扑优化的多变量水平分割方法及设备
CN111583243A (zh) * 2020-05-09 2020-08-25 中南大学 蜂窝规整度检测的相邻点胞元重构方法
CN111583243B (zh) * 2020-05-09 2021-03-30 中南大学 蜂窝规整度检测的相邻点胞元重构方法

Similar Documents

Publication Publication Date Title
CN107885929A (zh) 非均匀类蜂窝结构设计方法
CN107391855B (zh) 一种面向多种微观结构的材料结构一体化构建方法
CN105512385B (zh) 采用序列幂函数插值方法实现结构多材料拓扑优化的方法
Long et al. A meshless local Petrov-Galerkin method for solving the bending problem of a thin plate
CN106096158B (zh) 一种柔性铰链的拓扑优化设计方法
CN107563010B (zh) 基于形状特征的多尺度结构材料一体化设计方法
CN110941924B (zh) 一种多组件系统集成一体化的多尺度拓扑优化设计方法
CN103399996B (zh) 快速伺服刀架用柔性机构多目标拓扑优化设计方法
CN104765922A (zh) 基于保形约束的悬臂梁结构拓扑优化设计方法
JP2018537784A (ja) 異なる厚さの構造セグメントについて短縮された長さの境界を使用する形態最適化
CN116629079B (zh) 混合有限元空间构造及求解线弹性力学问题的方法及装置
CN108897935A (zh) 面向增材制造的点阵结构动力学响应优化设计方法
CN109670207A (zh) 一种面向多种多孔材料结构的动力学一体化设计方法
CN109657378A (zh) 一种含变尺寸单胞的非均质层级结构拓扑优化方法
CN113345536A (zh) 一种基于极限各向异性点阵材料的结构拓扑优化方法
CN102629291B (zh) 带孔结构和带组件结构的分析设计方法
CN114999591A (zh) 一种多构型点阵结构的拓扑优化方法
CN104535040B (zh) 用于叶片的有限元单元划分方法和叶片的检测方法
CN103294861A (zh) 多弧段曲线形状优化设计方法
CN103065015A (zh) 一种基于内力路径几何形态的承载结构低碳节材设计方法
CN109190173A (zh) 一种构造蜂窝夹层结构模型的方法
CN110717208B (zh) 一种基于连续梯度微结构的多尺度频率响应拓扑优化方法
Koltakov et al. Moving grid method for numerical simulation of stratified flows
CN103217906B (zh) 基于约束方程的固体重量压力载荷下的拓扑优化设计方法
Smith et al. Application specific mesh partition improvement

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: 20180406

WD01 Invention patent application deemed withdrawn after publication