CN107885929A - 非均匀类蜂窝结构设计方法 - Google Patents
非均匀类蜂窝结构设计方法 Download PDFInfo
- 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
- honeycomb
- cell
- mrow
- wall thickness
- mtd
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 210000004027 cell Anatomy 0.000 claims abstract description 42
- 238000005457 optimization Methods 0.000 claims abstract description 29
- 210000002421 cell wall Anatomy 0.000 claims abstract description 23
- 239000000463 material Substances 0.000 claims abstract description 21
- 238000000265 homogenisation Methods 0.000 claims abstract description 6
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 230000035945 sensitivity Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000013507 mapping Methods 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 210000003850 cellular structure Anatomy 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012679 convergent method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
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>&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>&OverBar;</mo>
</mover>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>t</mi>
<mi>min</mi>
</msub>
<mo><</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>&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)求得的灵敏度对各胞元壁厚进行优化，得到设计结果；根据得到的胞元壁厚信息，进行模型重构，获得最终的非均匀类蜂窝结构。
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 (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109766656A (zh) * | 2019-01-25 | 2019-05-17 | 北京航空航天大学 | 一种基于拓扑优化的梯度点阵结构设计方法 |
| CN111127491A (zh) * | 2019-12-31 | 2020-05-08 | 华中科技大学 | 一种针对蜂窝结构拓扑优化的多变量水平分割方法及设备 |
| CN111125818A (zh) * | 2019-12-12 | 2020-05-08 | 河南大学 | 基于仿生子结构拓扑优化的新型十字交叉节点及优化方法 |
| CN111583243A (zh) * | 2020-05-09 | 2020-08-25 | 中南大学 | 蜂窝规整度检测的相邻点胞元重构方法 |
| CN113962071A (zh) * | 2021-10-12 | 2022-01-21 | 武汉理工大学 | 非均匀中空立方体蜂窝结构及其设计方法 |
| CN118016215A (zh) * | 2024-02-21 | 2024-05-10 | 中南大学 | 一种非均匀壁厚蜂窝材料平台应力计算方法 |
| CN119851835A (zh) * | 2025-03-20 | 2025-04-18 | 中汽研汽车检验中心（天津）有限公司 | 一种蜂窝壁障仿真建模胞元参数确定方法 |
Citations (3)
| 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 | 西北工业大学 | 基于梯度微结构的材料/结构一体化设计方法 |
-
2017
- 2017-11-07 CN CN201711081269.7A patent/CN107885929A/zh active Pending
Patent Citations (4)
| 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)
| Title |
|---|
| PU ZHANG等: "Efficient Design-Optimization", 《JOURNAL OF MANUFACTURING SCIENCE AND ENGINEERING》 * |
| 邢昊等: "拓扑优化密度映射的非均匀蜂窝结构设计方法", 《计算机辅助设计与图形学学报》 * |
Cited By (12)
| 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 | 中南大学 | 蜂窝规整度检测的相邻点胞元重构方法 |
| CN113962071A (zh) * | 2021-10-12 | 2022-01-21 | 武汉理工大学 | 非均匀中空立方体蜂窝结构及其设计方法 |
| CN113962071B (zh) * | 2021-10-12 | 2024-08-13 | 武汉理工大学 | 非均匀中空立方体蜂窝结构及其设计方法 |
| CN118016215A (zh) * | 2024-02-21 | 2024-05-10 | 中南大学 | 一种非均匀壁厚蜂窝材料平台应力计算方法 |
| CN118016215B (zh) * | 2024-02-21 | 2024-10-15 | 中南大学 | 一种非均匀壁厚蜂窝材料平台应力计算方法 |
| CN119851835A (zh) * | 2025-03-20 | 2025-04-18 | 中汽研汽车检验中心（天津）有限公司 | 一种蜂窝壁障仿真建模胞元参数确定方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107885929A (zh) | 非均匀类蜂窝结构设计方法 | |
| Zhang et al. | Topological design of sandwich structures with graded cellular cores by multiscale optimization | |
| CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
| CN107391855B (zh) | 一种面向多种微观结构的材料结构一体化构建方法 | |
| CN104866673B (zh) | 一种轴压加筋柱壳的开口补强方法 | |
| CN107563010B (zh) | 基于形状特征的多尺度结构材料一体化设计方法 | |
| CN103978690A (zh) | 一种面向3d打印的物体内部结构优化方法 | |
| CN111079228B (zh) | 一种基于流场预测的气动外形优化方法 | |
| CN104679956B (zh) | 考虑动态特性的高速压力机底座可靠性稳健设计方法 | |
| CN103294857B (zh) | 以频响位移为优化目标的约束阻尼板的拓扑优化方法 | |
| CN109871574A (zh) | 一种基于代理模型的多尺度拓扑优化方法 | |
| CN107729648A (zh) | 一种基于Shepard插值的曲线纤维复合结构设计瀑布型多级优化方法 | |
| CN105528503A (zh) | 一种基于结构分解的大型构件动态优化设计方法 | |
| CN104091028A (zh) | 一种螺旋油楔轴承的多目标优化设计方法 | |
| CN101811257B (zh) | 机床床身结构优化设计方法 | |
| CN116629079B (zh) | 混合有限元空间构造及求解线弹性力学问题的方法及装置 | |
| CN104036095A (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
| CN104268317B (zh) | 机械零件圆角结构形状优化方法 | |
| CN112182806B (zh) | 一种力流引导的介观结构设计方法 | |
| CN114297779A (zh) | 一种飞行器气动力系数快速预测方法、系统及设备 | |
| Hartmann et al. | Enhanced battery pack for electric vehicle: Noise reduction and increased stiffness | |
| CN102629291B (zh) | 带孔结构和带组件结构的分析设计方法 | |
| CN107220413A (zh) | 基于梯度微结构的材料/结构一体化设计方法 | |
| CN118520607A (zh) | 一种考虑流固换热性能和刚度的曲面胞元结构及其优化设计方法 | |
| CN105138766B (zh) | 基于模糊聚类的高超声速气动热降阶模型的加点方法 |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180406 |