CN116401726A - 一种基于曲面密度分布的梯度极小曲面结构的设计方法 - Google Patents
一种基于曲面密度分布的梯度极小曲面结构的设计方法 Download PDFInfo
- Publication number
- CN116401726A CN116401726A CN202310671920.5A CN202310671920A CN116401726A CN 116401726 A CN116401726 A CN 116401726A CN 202310671920 A CN202310671920 A CN 202310671920A CN 116401726 A CN116401726 A CN 116401726A
- Authority
- CN
- China
- Prior art keywords
- curved surface
- gradient
- minimum
- equation
- surface structure
- 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.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 32
- 238000013461 design Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000010521 absorption reaction Methods 0.000 claims abstract description 24
- 239000000654 additive Substances 0.000 claims abstract description 7
- 230000000996 additive effect Effects 0.000 claims abstract description 7
- 238000004519 manufacturing process Methods 0.000 claims abstract description 7
- 238000005516 engineering process Methods 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims abstract description 6
- 238000012216 screening Methods 0.000 claims abstract description 5
- 238000006073 displacement reaction Methods 0.000 claims abstract description 4
- 238000012669 compression test Methods 0.000 claims description 4
- 230000001276 controlling effect Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000005457 optimization Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000010146 3D printing Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- 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
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
Abstract
本发明属于极小曲面结构设计技术领域,涉及一种基于曲面密度分布的梯度极小曲面结构的设计方法,包括以下步骤:建立结构密度分布的二维曲面;极小曲面结构模型建立;梯度结构模型建立;通过建模软件得到极小曲面结构的三维模型,通过增材制造技术制备梯度结构试样;通过试验获得的结构应力‑应变相应曲线和力‑位移曲线计算梯度结构的平台应力和能量吸收;评估不同类型的梯度极小曲面结构在力学响应和能量吸收性能方面的表现,筛选得到最优的结构类型,实现基于曲面密度梯度极小曲面结构的优化设计。本发明采用上述一种基于曲面密度分布的梯度极小曲面结构的设计方法,提升了极小曲面结构的力学响应稳定性和能量吸收性能。
Description
技术领域
本发明涉及极小曲面结构设计技术领域,特别是涉及一种基于曲面密度分布的梯度极小曲面结构的设计方法。
背景技术
极小曲面结构作为一种新型的多孔结构,由于其稳定的力学响应、优异的能量吸收性能及稳定的变形行为而受到研究人员的广泛关注。根据极小曲面建模时所采用数学方程的不同,极小曲面结构可以有多种不同的结构形式,常见的结构类型有Gyroid、Diamond、IWP等。对于Gyroid形式的极小曲面结构,其在轴向压缩载荷下的变形模式主要为“X”型的且结构整体会出现类似于鼓状变形。因此,结构在进行承载时其力学响应会出现一定的软化现象,而这种软化现象会在一定程度上降低其能量吸收性能。已有的研究结果表明梯度极小曲面结构在能量吸收方面有极大的提升,但是其持续的硬化行为并不利于对防护对象的保护。故减弱极小曲面的软化行为提升结构的能量吸收性能具有十分重要的意义。
一般的,极小曲面结构的优化设计主要有以下几种:(1)功能梯度设计,通过改变结构在一个方向上的壁厚/结构单胞数量来建立梯度极小曲面结构;(2)胞元混杂设计,通过几种不同的极小曲面结构之间相互组合进行设计;(3)多级结构设计,基于多级结构的设计理念将结构单胞以极小曲面的结构进行替换;(4)模型参数优化,通过改变在模型构建时方程的参数建立不同参数构型的结构模型。
以上方法并未充分解决结构在力学响应稳定和能量吸收提升方面存在的问题。虽然在一定程度上提升了结构的能量吸收性能,但是对于结构承载、受保护目标防护能力的提升等方面并未得到解决。
发明内容
本发明的目的是提供一种基于曲面密度分布的梯度极小曲面结构的设计方法,解决了Gyroid极小曲面结构在力学响应和能量吸收方面存在的问题,通过研究几种基于不同二维曲面类型的相对密度变化形式,进一步提升极小曲面结构的力学响应稳定性和能量吸收性能。
为实现上述目的,本发明提供了一种基于曲面密度分布的梯度极小曲面结构的设计方法,包括以下步骤:
S1、建立结构密度分布的二维曲面
基于二维曲面方程构建二维曲面,通过改变曲面函数中的不同参数值调整曲面中各个点的值大小;
S2、极小曲面结构模型建立
根据Gyroid极小曲面结构的建模方程建立极小曲面结构,通过改变方程控制参数来调整结构的壁厚,从而实现结构相对密度的变化;
S3、梯度结构模型建立
基于S1、S2中的二维曲面方程和极小曲面结构方程,在极小曲面方程中引入二维曲面方程,使结构在平面内的相对密度按照二维曲面值分布并得到梯度极小曲面结构的建模方程;
S4、在S3中得到的梯度结构方程基础上得到极小曲面结构的三维模型,通过增材制造技术制备梯度结构试样;
S5、对S4中获得的试样进行准静态压缩试验,通过试验获得的结构应力-应变相应曲线和力-位移曲线计算梯度结构的平台应力和能量吸收;
S6、评估不同类型的梯度极小曲面结构在力学响应和能量吸收性能方面的表现,筛选得到目标结构类型,实现基于曲面密度梯度极小曲面结构的优化设计。
优选的,步骤S1中不同形式的二维曲面包括二次曲面、三角函数曲面和马鞍面;对应的二维曲面方程分别为:
优选的,步骤S2中,Gyroid极小曲面结构的建模方程如式(1):
优选的,步骤S4中,构建的极小曲面结构三维模型的方程如式(2):
本发明的有益效果:
1)本发明公开的一种基于曲面密度分布的梯度极小曲面结构设计方法,利用二维曲面的分布规律特征将极小曲面结构的相对密度按照二维曲面的特征分布;通过增材制造技术制备结构试样并进行压缩试验,通过对比传统均匀结构筛选评估得到力学性能最优的结构。弥补了传统极小曲面结构出现应变软化的缺点,同时使结构的变形更加稳定。考虑3D打印几何缺陷的有限元计算模型,基于3D打印构件的真实几何信息创建,能够反映打印构件的主要几何特征。
2)本发明公开的一种基于曲面密度分布的梯度极小曲面结构设计方法,能够有效提升结构的平台应力强度。
3)本发明公开的一种基于曲面密度分布的梯度极小曲面结构设计方法,有效提升了结构变形稳定性和能量吸收性能。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1是本发明的三种类型的二维曲面图;
图2是本发明的均匀Gyroid结构和三种梯度结构的3D模型及密度分布规律图;
图3是本发明的三种梯度极小曲面结构图;
图4是本发明的三种梯度极小曲面和均匀Gyroid极小曲面结构的力学响应和平台应力曲线图;
图5为本发明中三种梯度极小曲面和均匀Gyroid极小曲面结构的能量吸收柱状图;
图6为本发明中基于马鞍面梯度结构变形过程的示意图。
具体实施方式
下面结合实施例,对本发明进一步描述。除非另外定义,本发明使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本发明提到的上述特征或具体实例提到的特征可以任意组合,这些具体实施例仅用于说明本发明而不用于限制本发明的范围。
实施例
S1、建立结构密度分布的二维曲面
基于不同类型的二维曲面方程,采用数学软件构建不同形式的二维曲面。
图1是本发明的三种类型的二维曲面图,其中图1中的(a)为二次曲面,图1中的(b)为三角函数曲面,图1中的(c)为马鞍面,如图所示,不同形式的二维曲面包括二次曲面、三角函数曲面和马鞍面。对应的二维曲面方程分别为:
通过改变曲面函数中的不同参数值调整曲面中各个点的值大小。
S2、极小曲面结构模型建立
根据Gyroid极小曲面结构的建模方程建立极小曲面结构。Gyroid极小曲面结构的建模方程如式(1):
S3、梯度结构模型建立
基于S1、S2中的二维曲面方程和极小曲面结构方程,在极小曲面方程中引入二维曲面方程,使结构在平面内的相对密度按照二维曲面值分布并得到梯度极小曲面结构的建模方程。
S4、在S3中得到的梯度结构方程基础上通过建模软件得到极小曲面结构的三维模型,通过增材制造技术制备梯度结构试样;
构建的极小曲面结构三维模型的方程如式(2):
本发明中梯度极小曲面结构是使极小曲面结构的相对密度在平面内按照二维曲面的变化规律进行分布,将表达式/>用二维曲面方程进行表示,使极小曲面结构在/>平面内相对密度的变化分布规律与二维曲面中/>值的大小分布相同。
图2是本发明的均匀Gyroid结构和三种梯度结构的3D模型及密度分布规律图,图2中的(a)为均匀Gyroid结构,图2中的(b)为基于二次曲面的相对密度分布规律,图2中的(c)为基于三角函数曲面的相对密度分布规律,图2中的(d)为基于马鞍面的相对密度分布规律,如图2中的(b)、(c)、(d)所示分别为三种梯度结构相对密度的分布规律。其中,极小曲面结构在平面内每一个点位置的相对密度与曲面值的大小一一对应。通过调控曲面方程中a,b,和c的值来实现极小曲面结构在/>平面内相对密度大小的变化区间,采用数学建模软件编写程序生成三维几何模型。
图3是本发明的三种梯度极小曲面结构图,其中图3中的(a)为基于二次曲面的梯度结构图,图3中的(b)为基于三角曲面的梯度结构图,图3中的(c)为基于二次曲面的梯度结构图,如图所示,图3中的(a)、(b)和(c)分别为基于二次曲面、三角函数曲面和马鞍面三种曲面生成得到的梯度极小曲面结构。
完成结构设计并建立三维梯度结构模型之后,采用FDM增材制造技术制备梯度Groid极小曲面结构,结构的基体材料采用增强的PLA材料。增材制造得到的试样的各项参数如表1所示。采用准静态实验机开展试验,几种结构进行准静态压缩时的加载速度均为1.5mm/min。
表1
试样编号 | 实际尺寸/mm | 质量/g | 加载速度 |
JY-0-1# | 25.06×25.02×25.46 | 4.95 | 1.5mm/min |
JY-0-2# | 25.10×24.98×25.54 | 4.93 | 1.5mm/min |
TD-1-1# | 25.10×25.08×25.08 | 4.74 | 1.5mm/min |
TD-1-2# | 25.04×25.06×25.08 | 4.73 | 1.5mm/min |
TD-2-1# | 25.10×25.16×25.02 | 4.89 | 1.5mm/min |
TD-2-2# | 25.06×25.12×25.04 | 4.90 | 1.5mm/min |
TD-3-1# | 25.00×25.10×25.10 | 5.10 | 1.5mm/min |
TD-3-2# | 25.00×25.08×25.00 | 5.10 | 1.5mm/min |
S5、对S4中获得的试样进行准静态压缩试验,通过试验获得的结构应力-应变相应曲线和力-位移曲线计算梯度结构的平台应力和能量吸收。
平台应力可以通过公式(3)得到:
能量吸收可以通过公式(5)计算得到:
S6、评估不同类型的梯度极小曲面结构在力学响应和能量吸收性能方面的表现,筛选得到最优的结构类型,实现基于曲面密度梯度极小曲面结构的优化设计。
S6.1吸能指标:根据实验结果,计算所有结构的平台应力和总能量吸收。
S6.2变形模式特性:通过对比不同梯度形式结构和均匀结构的变形模式,筛选变形更加稳定的梯度结构类型。
S6.3优化结果:根据研究结果选出变形及响应更急稳定的梯度结构,实现基于曲面梯度极小曲面结构的性能优化设计。
图4是本发明的三种梯度极小曲面和均匀Gyroid极小曲面结构的力学响应和平台应力曲线图,其中图4中的(a)为三种梯度极小曲面和均匀Gyroid极小曲面结构的应力-应变曲线,图4中的(b)为三种梯度极小曲面和均匀Gyroid极小曲面结构的平台应力,如图4所示,三种梯度结构的应力-应变曲线中的软化现象明显降低,同时基于马鞍面的梯度结构在曲线平台部分高出均匀结构约18.32%。
图5为本发明中三种梯度极小曲面和均匀Gyroid极小曲面结构的能量吸收柱状图,从图5中可以看出,基于曲面密度分布的梯度极小曲面结构在相同相对密度情况下,能量吸收性能明显提升,实现了结构力学性能的大幅优化。基于马鞍面的梯度结构相较于均匀结构能量吸收提升了19.21%,表明基于马鞍面的梯度结构可以在提升结构平台应力的同时使结构的承载力曲线更加平稳,更有利于吸收能量并减少受保护目标受到高峰值载荷的冲击。
图6为本发明中基于马鞍面梯度结构变形过程的示意图,从图6中的基于马鞍面的梯度结构作为力学性能最优异的结构在受到压缩载荷时的变形过程,可以看出结构的变形非常稳定并未出现明显的局部变形,这也充分说明了基于曲面分布的结构梯度特性设计实现了对极小曲面结构变形过程的调控,减小了结构承载力曲线中的软化现象,使结构的整体力学性能得到了提升。
最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。
Claims (4)
1.一种基于曲面密度分布的梯度极小曲面结构的设计方法,其特征在于,包括以下步骤:
S1、建立结构密度分布的二维曲面
基于二维曲面方程构建二维曲面,通过改变曲面函数中的不同参数值调整曲面中各个点的值大小;
S2、极小曲面结构模型建立
根据Gyroid极小曲面结构的建模方程建立极小曲面结构,通过改变方程的控制参数来调整结构的壁厚,实现结构相对密度的变化;
S3、梯度结构模型建立
基于S1、S2中的二维曲面方程和Gyroid极小曲面结构方程,在极小曲面结构方程中引入二维曲面方程,使极小曲面结构在平面内的相对密度按照二维曲面值分布并得到梯度极小曲面结构的建模方程;
S4、在S3中得到的梯度结构方程基础得到极小曲面结构的三维模型,通过增材制造技术制备梯度结构试样;
S5、对S4中获得的试样进行准静态压缩试验,通过试验获得的结构应力-应变相应曲线和力-位移曲线计算梯度结构的平台应力和能量吸收;
S6、评估不同类型的梯度极小曲面结构在力学响应和能量吸收性能方面的表现,筛选得到目标结构类型,实现基于曲面密度梯度极小曲面结构的优化设计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310671920.5A CN116401726B (zh) | 2023-06-08 | 2023-06-08 | 一种基于曲面密度分布的梯度极小曲面结构的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310671920.5A CN116401726B (zh) | 2023-06-08 | 2023-06-08 | 一种基于曲面密度分布的梯度极小曲面结构的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116401726A true CN116401726A (zh) | 2023-07-07 |
CN116401726B CN116401726B (zh) | 2023-08-11 |
Family
ID=87010900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310671920.5A Active CN116401726B (zh) | 2023-06-08 | 2023-06-08 | 一种基于曲面密度分布的梯度极小曲面结构的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116401726B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859693A (zh) * | 2020-07-28 | 2020-10-30 | 大连理工大学 | 一种高效的多孔结构表示和优化方法 |
CN112395746A (zh) * | 2020-11-03 | 2021-02-23 | 山东大学 | 微结构族等效材料性质的计算方法、微结构、系统及介质 |
CN113158526A (zh) * | 2021-05-07 | 2021-07-23 | 杭州电子科技大学 | 一种基于tpms可控孔隙结构设计方法 |
US20220129595A1 (en) * | 2020-06-28 | 2022-04-28 | Dalian University Of Technology | Design and optimization method of porous structure for 3d heat dissipation based on triply periodic minimal surface (tpms) |
CN116059012A (zh) * | 2023-02-01 | 2023-05-05 | 武汉理工大学 | 一种多维空间梯度孔型极小曲面骨植入体材料及其设计方法 |
-
2023
- 2023-06-08 CN CN202310671920.5A patent/CN116401726B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220129595A1 (en) * | 2020-06-28 | 2022-04-28 | Dalian University Of Technology | Design and optimization method of porous structure for 3d heat dissipation based on triply periodic minimal surface (tpms) |
CN111859693A (zh) * | 2020-07-28 | 2020-10-30 | 大连理工大学 | 一种高效的多孔结构表示和优化方法 |
CN112395746A (zh) * | 2020-11-03 | 2021-02-23 | 山东大学 | 微结构族等效材料性质的计算方法、微结构、系统及介质 |
CN113158526A (zh) * | 2021-05-07 | 2021-07-23 | 杭州电子科技大学 | 一种基于tpms可控孔隙结构设计方法 |
CN116059012A (zh) * | 2023-02-01 | 2023-05-05 | 武汉理工大学 | 一种多维空间梯度孔型极小曲面骨植入体材料及其设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116401726B (zh) | 2023-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11238194B2 (en) | Structural design method of product | |
CN110008512B (zh) | 一种考虑承载特性的负泊松比点阵结构拓扑优化方法 | |
Zhao et al. | Design and mechanical performances of a novel functionally graded sheet-based lattice structure | |
Liang et al. | Parametric design strategy of a novel self-similar hierarchical honeycomb for multi-stage energy absorption demand | |
CN101691012B (zh) | 箱形支撑结构内部加筋板分布优化设计方法 | |
CN110929358A (zh) | 一种参数化功能梯度立方点阵结构设计方法及材料 | |
CN111709085B (zh) | 一种约束阻尼薄板结构拓扑优化设计方法 | |
Feng et al. | Stiffness optimization design for TPMS architected cellular materials | |
CN110795873A (zh) | 一种考虑尺寸控制的跨尺度拓扑优化方法 | |
CN112182929A (zh) | 一种考虑尺寸控制的多孔材料跨尺度可靠性拓扑优化方法 | |
CN112287491B (zh) | 一种复合型点阵材料及其设计方法 | |
CN114636360B (zh) | 五模冲击隐身复合点阵环状结构及其参数优化方法 | |
CN113158273A (zh) | 一种恒定孔隙尺寸的极小曲面连续梯度多孔结构生成方法 | |
CN110210079A (zh) | 一种面向整机动态特性的机床支承件质量匹配方法 | |
CN114756934B (zh) | 一种三维多尺度超材料结构优化设计方法 | |
CN116522499A (zh) | 一种基于薄壁–点阵填充结构车身前端结构设计方法 | |
CN112446163A (zh) | 基于参数化水平集的能量有限元拓扑优化方法 | |
Bai et al. | Influence of relative density distribution rules on the mechanical compression responses of additive manufactured Ti6Al4V graded lattice structures | |
CN115203997A (zh) | 一种基于多变量设计的点阵-实体复合结构拓扑优化方法 | |
CN116401726B (zh) | 一种基于曲面密度分布的梯度极小曲面结构的设计方法 | |
CN111737908B (zh) | 一种基于动载荷静力等效的蒙皮桁条结构快速动态优化设计方法 | |
Zhang et al. | Protective performance of hybrid triply periodic minimal surface lattice structure | |
CN110704912B (zh) | 一种应力约束下的桥架托臂结构拓扑优化方法 | |
Shi et al. | Mechanical Response and Energy Absorption of 316L Stainless Steel Body‐Centered‐Cubic Functionally Graded Lattice Structures under Quasistatic Compressive Test | |
CN111659892B (zh) | 基于对角形单元沉淀型微桁架结构的复合吸能结构及其3d打印方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |