CN109783950A - 增材制造中连通结构的拓扑优化设计方法 - Google Patents
增材制造中连通结构的拓扑优化设计方法 Download PDFInfo
- Publication number
- CN109783950A CN109783950A CN201910056048.7A CN201910056048A CN109783950A CN 109783950 A CN109783950 A CN 109783950A CN 201910056048 A CN201910056048 A CN 201910056048A CN 109783950 A CN109783950 A CN 109783950A
- Authority
- CN
- China
- Prior art keywords
- design
- function
- porous characteristic
- connectivity
- present
- 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
Landscapes
- Image Generation (AREA)
Abstract
本发明公开了一种增材制造中连通结构的拓扑优化设计方法，用于解决现有连通结构的拓扑优化方法实用性差的技术问题。技术方案是在实体设计区域内布置一定数目封闭的孔洞特征，通过孔洞的移动、变形、融合、缩小和膨胀等行为来驱动结构的拓扑布局演变。此外，每个孔洞特征的中心点被控制在设计域之外，使孔洞特征不能完全进入到设计域之内形成封闭的孔洞，来达到结构连通的目的。相比背景技术的设计方法，本发明无需进行温度场分析，计算量少。此外，设计变量和单元无关，数目少、易收敛。本发明设计的结构边界光滑，内部不存在中间密度单元，有利于工程设计人员进行模型重构等后处理过程，实用性好。
Description
技术领域
本发明涉及一种连通结构的拓扑优化方法，特别涉及一种增材制造中连通结构的拓扑优化设计方法。
背景技术
增材制造过程中，对于内部含封闭孔洞的结构，由于无法去除支撑材料或者未熔融的金属粉末，需要二次修正或者结构分区制造，这大大增加了制造的工艺难度，提高了成本。因此在结构拓扑优化设计时考虑结构的连通性，即在寻找最优结构的同时满足结构内部不存在封闭孔洞的要求，具有非常重要的工程应用价值。
文献“Q.Li,W.Chen,S.Liu,L.Tong,Structural topology optimizationconsidering connectivity constraint,Structural and MultidisciplinaryOptimization,54(2016)971-984”公开了一种连通结构构型的拓扑优化方法。文献提出的虚拟温度场理论是以密度法为基础，空单元中存在高导热且自发热的材料，实体单元中则填充隔热材料。通过这种方法，结构封闭孔洞处产生的热量不能耗散，保持较高的温度。反之，对于连通结构，空单元产生的热量直接传递到散热边界，使整体结构温度降低。所以，文献将连通性条件转化为结构的最大温度小于有限值的约束，并应用于结构拓扑优化设计中，实现了最优拓扑结构的连通性控制。
文献所述方法虽然可以得到连通结构，但是为了实现结构的连通性，除了力学分析之外还需进行一次温度场分析来获得结构的最大温度，增加了计算量。此外，设计变量数目多、收敛慢也是不可忽视的问题。中间密度单元的存在使设计的结果不明确，且边界成锯齿状，这大大地增加了重构的工作量，难以应用于工程实际。
发明内容
为了克服现有连通结构的拓扑优化方法实用性差的不足，本发明提供一种增材制造中连通结构的拓扑优化设计方法。该方法在实体设计区域内布置一定数目封闭的孔洞特征，通过孔洞的移动、变形、融合、缩小和膨胀等行为来驱动结构的拓扑布局演变。此外，每个孔洞特征的中心点被控制在设计域之外，使孔洞特征不能完全进入到设计域之内形成封闭的孔洞，来达到结构连通的目的。相比背景技术的设计方法，本发明无需进行温度场分析，计算量少。此外，设计变量和单元无关，数目少、易收敛。本发明设计的结构边界光滑，内部不存在中间密度单元，有利于工程设计人员进行模型重构等后处理过程，实用性好。
本发明解决其技术问题所采用的技术方案：一种增材制造中连通结构的拓扑优化设计方法，其特点是包括以下步骤：
步骤一、实体区域初始分布m个孔洞特征，并且孔洞特征的中心点全部位于实体区域的边界上。
步骤二、利用KS函数得到整体结构的水平集函数Φ。
式中，φi为第i个孔洞特征的水平集函数。ln和e为自然对数和自然指数计算，p为KS函数的参数，且p>0，代表布尔并操作，负号代表特征是孔洞特征。
步骤三、采用固定网格将实体区域离散，同时定义载荷和边界条件。
步骤四、定义拓扑优化问题为：
Min J＝FU
式中J，V和表示结构的柔顺度，总体积和最大体积约束。K、F和U分别表示结构的总体刚度矩阵，总体载荷向量和位移向量。Ω表示实体区域，H是指Heaviside函数，用来筛选参与计算的积分点，x是积分点的坐标。是设计变量向量，包括了优化问题中所有的设计变量,dt表示第t个设计变量，它的上限和下限分别是d t和
步骤五、对上面建立的模型进行一次有限元分析，分别对目标函数和约束函数进行灵敏度分析。选取优化算法GCMMA进行优化设计，得到最优化结果。
本发明的有益效果是：该方法在实体设计区域内布置一定数目封闭的孔洞特征，通过孔洞的移动、变形、融合、缩小和膨胀等行为来驱动结构的拓扑布局演变。此外，每个孔洞特征的中心点被控制在设计域之外，使孔洞特征不能完全进入到设计域之内形成封闭的孔洞，来达到结构连通的目的。相比背景技术的设计方法，本发明无需进行温度场分析，计算量少。此外，设计变量和单元无关，数目少、易收敛。本发明设计的结构边界光滑，内部不存在中间密度单元，有利于工程设计人员进行模型重构等后处理过程，实用性好。
本发明方法经过实施例150步迭代后得到设计结果，收敛速度快。初始结构的体积为30262.40，柔顺度为219.07，优化得到的结构柔顺度为238.37，体积为9600，体积减少了68.28％，柔顺度仅增加了8.81％。在不考虑连通性的情况下，即不限制孔洞特征中心点的位置，用相同初始结构优化得到的结果内部存在一个封闭的孔洞，结构柔顺度为232.50，体积同样为9600，所以实现结构连通性仅牺牲了1.25％的柔顺度。本实施例的结果结构边界清晰，且不存在中间密度。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明方法实施例中模型几何尺寸与边界条件示意图。
图2是本发明实施例方法中孔洞特征初始布局的示意图。
图3是本发明方法实施例的设计结果图。
图4是本发明方法实施例设计结果的截面图。
图5是本发明方法实施例不考虑结构连通性的设计结果图。
图6是本发明方法实施例不考虑结构连通性设计结果的截面图。
具体实施方式
参照图1-6。
以简支六面体梁结构为例说明本发明。六面体梁结构尺寸长度为40，宽40，高20，底面存在厚度为2的非设计域。材料杨氏模量为1，泊松比为0.3。设计六面体梁中的材料分布使得其柔顺度最小，总材料用量体积分数最大为30％。具体步骤为：
步骤一、定义初始结构。在六面体实体区域上初始分布13个B样条封闭曲面孔洞特征，其中5个孔洞特征的中心点位于实体区域的上表面，其余8个孔洞特征的中心点分布在实体区域的4个侧表面。每个孔洞特征有120个控制点，给定中心点到每个控制点距离的初始值，使每个孔洞特征的初始形状为球形。
步骤二、水平集函数的构建。根据每个B样条封闭曲面孔洞特征的中心点坐标以及中心点到B样条曲面控制点的距离，得到所有孔洞特征的水平集函数，利用KS函数将所有孔洞特征进行布尔并操作，得到整体结构的水平集函数，计算公式为：
式中，ln和e为自然对数和自然指数计算，φi为第i个孔洞特征的水平集函数。Φ为整体结构的水平集函数，p为KS函数的参数，p＝7，代表布尔并操作。负号表示特征为孔洞特征。
步骤三、采用0.8×0.8×0.8的正六面体网格将实体区域离散。定义边界条件为：六面体梁底面4个角点的自由度全部固定，在底面的中点处施加竖直向下的载荷，载荷大小为10。
步骤四、定义拓扑优化问题：
Min J＝FU
式中J，V和表示结构的柔顺度，总体积和最大体积约束。K、F和U分别表示结构的总体刚度矩阵，总体载荷向量和位移向量。Ω表示实体区域，H是Heaviside函数，用来筛选参与计算的积分点，x是积分点坐标。是设计变量向量,dt表示第t个设计变量，t＝1,2,…,13×(120+3)，它的上限和下限分别是d t和本实施例的设计变量包括中心点坐标，中心点到每个控制点的距离两种类型。根据每个孔洞特征的位置确定中心点坐标设计变量的上下限，限制孔洞特征的中心点在优化过程中始终位于实体区域之外。距离设计变量的范围为[0,20]。
步骤五、有限元分析和优化求解。对模型进行有限元分析以及对目标函数和约束函数的灵敏度的求解都是通过Matlab的编程实现，再通过结构优化平台Boss-Quattro中选取优化算法GCMMA对该实施例进行优化设计，得到最优的结果。
本发明方法经过实施例150步迭代后得到设计结果，收敛速度快。初始结构的体积为30262.40，柔顺度为219.07，优化得到的结构柔顺度为238.37，体积为9600，体积减少了68.28％，柔顺度仅增加了8.81％。在不考虑连通性的情况下，即不限制孔洞特征中心点的位置，用相同初始结构优化得到的结果内部存在一个封闭的孔洞，结构柔顺度为232.50，体积同样为9600，所以实现结构连通性仅牺牲了1.25％的柔顺度。本实施例的结果结构边界清晰，且不存在中间密度，工程实用性强。
Claims (1)
1.一种增材制造中连通结构的拓扑优化设计方法，其特征在于包括以下步骤：
步骤一、实体区域初始分布m个孔洞特征，并且孔洞特征的中心点全部位于实体区域的边界上；
步骤二、利用KS函数得到整体结构的水平集函数Φ；
式中，vi为第i个孔洞特征的水平集函数；ln和e为自然对数和自然指数计算，p为KS函数的参数，且p>0，代表布尔并操作，负号代表特征是孔洞特征；
步骤三、采用固定网格将实体区域离散，同时定义载荷和边界条件；
步骤四、定义拓扑优化问题为：
Min J＝FU
式中J，V和表示结构的柔顺度，总体积和最大体积约束；K、F和U分别表示结构的总体刚度矩阵，总体载荷向量和位移向量；Ω表示实体区域，H是指Heaviside函数，用来筛选参与计算的积分点，x是积分点的坐标；是设计变量向量，包括了优化问题中所有的设计变量,dt表示第t个设计变量，它的上限和下限分别是d t和
步骤五、对上面建立的模型进行一次有限元分析，分别对目标函数和约束函数进行灵敏度分析；选取优化算法GCMMA进行优化设计，得到最优化结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910056048.7A CN109783950A (zh) | 2019-01-22 | 2019-01-22 | 增材制造中连通结构的拓扑优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910056048.7A CN109783950A (zh) | 2019-01-22 | 2019-01-22 | 增材制造中连通结构的拓扑优化设计方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109783950A true CN109783950A (zh) | 2019-05-21 |
Family
ID=66501761
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910056048.7A Pending CN109783950A (zh) | 2019-01-22 | 2019-01-22 | 增材制造中连通结构的拓扑优化设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109783950A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111319268A (zh) * | 2020-02-20 | 2020-06-23 | 西北工业大学 | 一种考虑增材制造打印方向的自支撑结构优化设计方法 |
| CN112233242A (zh) * | 2020-10-09 | 2021-01-15 | 西北工业大学 | 一种三维自支撑结构的拓扑优化设计方法 |
| CN115004188A (zh) * | 2019-08-28 | 2022-09-02 | 西门子工业软件有限公司 | 用于处理经拓扑优化的几何结构的***和方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107341316A (zh) * | 2017-07-13 | 2017-11-10 | 西北工业大学 | 设计相关压力载荷作用下的结构形状‑拓扑联合优化方法 |
| CN107391824A (zh) * | 2017-07-11 | 2017-11-24 | 西北工业大学 | 增材制造中自支撑结构的拓扑优化设计方法 |
-
2019
- 2019-01-22 CN CN201910056048.7A patent/CN109783950A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107391824A (zh) * | 2017-07-11 | 2017-11-24 | 西北工业大学 | 增材制造中自支撑结构的拓扑优化设计方法 |
| CN107341316A (zh) * | 2017-07-13 | 2017-11-10 | 西北工业大学 | 设计相关压力载荷作用下的结构形状‑拓扑联合优化方法 |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115004188A (zh) * | 2019-08-28 | 2022-09-02 | 西门子工业软件有限公司 | 用于处理经拓扑优化的几何结构的***和方法 |
| CN111319268A (zh) * | 2020-02-20 | 2020-06-23 | 西北工业大学 | 一种考虑增材制造打印方向的自支撑结构优化设计方法 |
| CN111319268B (zh) * | 2020-02-20 | 2021-12-28 | 西北工业大学 | 一种考虑增材制造打印方向的自支撑结构优化设计方法 |
| CN112233242A (zh) * | 2020-10-09 | 2021-01-15 | 西北工业大学 | 一种三维自支撑结构的拓扑优化设计方法 |
| CN112233242B (zh) * | 2020-10-09 | 2022-08-05 | 西北工业大学 | 一种三维自支撑结构的拓扑优化设计方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10613496B2 (en) | Support structure constrained topology optimization for additive manufacturing | |
| US20210073428A1 (en) | Structure topology optimization method based on material-field reduced series expansion | |
| CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
| CN111737835B (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
| CN109783950A (zh) | 增材制造中连通结构的拓扑优化设计方法 | |
| CN107563010B (zh) | 基于形状特征的多尺度结构材料一体化设计方法 | |
| CN107391824A (zh) | 增材制造中自支撑结构的拓扑优化设计方法 | |
| WO2022000132A1 (zh) | 基于三周期极小曲面的三维多孔散热结构的设计与优化方法 | |
| CN105313336B (zh) | 一种薄壳体3d打印优化方法 | |
| CN105787226A (zh) | 四边有限元网格模型的参数化模型重建 | |
| CN109543207B (zh) | 考虑可变分型线实现双模铸造件多组件设计的方法 | |
| CN104036095B (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
| CN111428397B (zh) | 一种考虑增材制造结构自支撑约束的拓扑优化设计方法 | |
| CN109255185A (zh) | 一种基于城市地表地下管网的一二维水动力学耦合分析方法 | |
| CN111768502A (zh) | 一种基于gpu加速技术的非结构网格二维洪水模拟*** | |
| CN110414127B (zh) | 一种面向增材制造的支撑体积约束拓扑优化方法 | |
| CN109255164B (zh) | 一种基于空间拓扑的地表二维地下管网水动力学耦合方法 | |
| Feng et al. | An improved two-level support structure for extrusion-based additive manufacturing | |
| CN106373184A (zh) | 一种三维打印模型摆放所需支撑量快速估算方法 | |
| CN110210151B (zh) | 基于b样条的点阵结构参数化隐式建模与优化方法 | |
| CN111523270A (zh) | 一种改进的连续体结构拓扑优化后处理方法 | |
| CN109657378A (zh) | 一种含变尺寸单胞的非均质层级结构拓扑优化方法 | |
| CN106844963B (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
| Zhu et al. | Design of lightweight tree-shaped internal support structures for 3D printed shell models | |
| Meng et al. | Shape–thickness–topology coupled optimization of free-form shells |
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: 20190521 |
|
| WD01 | Invention patent application deemed withdrawn after publication |