CN110929358B - 一种参数化功能梯度立方点阵结构设计方法及材料 - Google Patents

一种参数化功能梯度立方点阵结构设计方法及材料 Download PDF

Info

Publication number
CN110929358B
CN110929358B CN202010069303.4A CN202010069303A CN110929358B CN 110929358 B CN110929358 B CN 110929358B CN 202010069303 A CN202010069303 A CN 202010069303A CN 110929358 B CN110929358 B CN 110929358B
Authority
CN
China
Prior art keywords
gradient
lattice
cell
layer
continuum
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.)
Active
Application number
CN202010069303.4A
Other languages
English (en)
Other versions
CN110929358A (zh
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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202010069303.4A priority Critical patent/CN110929358B/zh
Publication of CN110929358A publication Critical patent/CN110929358A/zh
Application granted granted Critical
Publication of CN110929358B publication Critical patent/CN110929358B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种参数化功能梯度立方点阵结构设计方法及材料，包括：1、获取连续体几何结构所需承载的工况载荷参数；2、对连续体几何结构进行冲击性能仿真分析，得到应力分布规律，并根据应力分布规律建立密度梯度排布函数；3、对连续体几何结构进行分层，根据密度梯度排布函数对每层的点阵胞元杆件尺寸进行梯度惩罚计算，得到各层点阵胞元杆件尺寸；4、根据各层点阵胞元杆件尺寸，进行点阵胞元的生成和循环装配，得到功能梯度点阵结构材料模型。本发明通过密度函数使点阵胞元尺寸和相对密度沿高度方向呈梯度连续变化，材料分布趋向于应力较大的区域，结构变形和应力分布趋于均匀化，与均匀点阵结构相比，材料分布更加合理。

Description

一种参数化功能梯度立方点阵结构设计方法及材料
技术领域
本发明属于结构轻量化和轻质功能材料技术领域，尤其涉及一种参数化功能梯度立方点阵结构设计方法及材料。
背景技术
在军用无人装备和航空航天领域，结构轻量化是减轻设备重量、提高机动性、降低成本的重要技术途径。点阵结构材料因其高孔隙率、低相对密度、周期性排布等结构特点，具有轻质化、高比刚度和高比强度等性能优势，在产品设备结构轻量化设计方面具有重要的应用价值。在众多点阵结构材料中，立方点阵结构由于其包络空间为标准立方体，胞元结构形状规则，便于进行点阵排布和映射求解，具备良好的可设计性和适应性。
采用点阵结构材料进行结构轻量化设计，连续材料变成轻质多孔材料，材料微观结构发生变化，力学特性随之改变。在实际应用中，设备运行过程往往存在一些动载荷因素，如振动、冲击等。零件结构受动载荷作用时，主要通过结构和材料的变形发挥承载和缓冲作用，点阵结构材料的变形量、抗冲击强度和能量吸收率等性能与点阵胞元类型、胞元尺寸、密度分布、材料属性等因素密切相关。研究发现，均匀点阵结构材料在承受动载荷作用时，尤其是高速或高频动载荷，容易产生应力集中，导致材料表面或深层的变形、压溃、断裂等问题，严重影响了零件的精度和可靠性。
发明内容
本发明要解决的技术问题是针对现有点阵结构材料在振动冲击等动载荷作用下，导致材料表面或深层的变形、压溃、断裂等问题，严重影响了零件的精度和可靠性的问题，提出了一种减小应力集中，使结构变形和能量耗散趋于均匀化，提高材料抗振动、抗冲击性能的参数化功能梯度立方点阵结构设计方法及材料。
为解决该问题，本发明所采用的技术方案是：
一种参数化功能梯度立方点阵结构设计方法，包括以下步骤：
步骤1：获取连续体几何结构所需承载的工况载荷参数，以及连续体几何结构的点阵胞元类型、单个胞元的主要设计参数P(a,t)，a表示胞元包络立方体边长尺寸，t表示胞元杆件截面边长尺寸；
步骤2：根据所述工况载荷参数，对所述连续体几何结构进行冲击性能仿真分析，得到冲击载荷作用下连续体几何结构内部应力分布规律，并根据应力分布规律建立密度梯度排布函数F(h)，h表示几何结构在应力梯度变化方向上的高度；
步骤3：对连续体几何结构进行分层，根据所述密度梯度排布函数对每层的点阵胞元杆件尺寸进行梯度惩罚计算，得到各层点阵胞元杆件尺寸；
步骤4：根据所述各层点阵胞元杆件尺寸，使用三维软件建模程序进行点阵胞元的生成和循环装配，得到功能梯度点阵结构材料模型。
进一步地，还包括步骤5：对所述功能梯度点阵结构模型进行冲击仿真分析，判断所述功能梯度点阵结构模型是否满足设计要求，若满足，则根据该点阵结构模型得到点阵结构材料，若不满足，则进行结构优化设计直到满足设计要求。
进一步地，对连续体几何结构进行分层是指沿连续体几何结构的梯度变化方向，将点阵胞元排布层的一层或多层分成一个连续体几何结构层。
进一步地，步骤3中根据所述密度梯度排布函数对点阵胞元尺寸进行梯度惩罚计算的方法是指：第k层点阵点阵胞元的杆件尺寸
Figure 628568DEST_PATH_IMAGE001
，
Figure 514484DEST_PATH_IMAGE002
表示初始层胞元杆件尺寸，
Figure 637161DEST_PATH_IMAGE003
表示几何结构沿应力梯度变化方向第
Figure 862606DEST_PATH_IMAGE004
层点阵对应的梯度密度排布函数值。
进一步地，步骤4中所述进行点阵胞元生成和循环装配的方法是：
步骤4.1：给初始点阵胞元尺寸赋值(a,t 0 )，生成初始立方点阵胞元元件；
步骤4.2：设定装配初始原点p、装配基准，并和点阵胞元元件基准相匹配，
步骤4.3：根据步骤3中得到的各层点阵胞元尺寸，生成点阵胞元元件，从第一层第一个点阵胞元元件开始装配，直到完成所有层上点阵胞元的装配，生成梯度点阵结构模型。
进一步地，步骤2中所述根据应力分布规律建立密度梯度排布函数F(h)方法是：
步骤2.1：根据冲击性能仿真分析得到冲击载荷作用下连续体几何结构内部应力分布规律，确定出应力分布函数σ(h)；
步骤2.2：根据所述应力分布函数σ(h)建立初始密度梯度排布函数
Figure 45326DEST_PATH_IMAGE005
：
Figure 621801DEST_PATH_IMAGE006
步骤2.3：根据连续体几何结构的体积、相对密度以及减重比，确定在相同减重比的情况下，初始密度梯度排布函数
Figure 282589DEST_PATH_IMAGE007
的放大系数
Figure 642768DEST_PATH_IMAGE008
：
Figure 793127DEST_PATH_IMAGE009
步骤2.4：输出密度梯度排布函数
Figure 325739DEST_PATH_IMAGE010
：
Figure 55798DEST_PATH_IMAGE011
其中
Figure 52573DEST_PATH_IMAGE012
，、
Figure 577095DEST_PATH_IMAGE013
、
Figure 597004DEST_PATH_IMAGE014
为最小应力、最大应力和第
Figure 865174DEST_PATH_IMAGE015
层可变应力可变应力；
Figure 716455DEST_PATH_IMAGE016
、
Figure 411879DEST_PATH_IMAGE017
、
Figure 184663DEST_PATH_IMAGE018
为胞元最小杆件尺寸、 最大或初始杆件尺寸和可变尺寸；
Figure 990945DEST_PATH_IMAGE019
为相对密度，
Figure 965242DEST_PATH_IMAGE020
为点阵胞元实体杆件体积，
Figure 831566DEST_PATH_IMAGE021
为包络体积计算得出；
Figure 91646DEST_PATH_IMAGE022
为减重比，
Figure 436040DEST_PATH_IMAGE023
为连续体几何结构分布层数。
进一步地，步骤5中所述结构优化设计通过改变梯度密度排布函数F(h)来实现。
进一步地，本发明还提供了一种参数化功能梯度立方点阵结构材料，该材料使用上述设计方法得到的点阵结构材料模型制作而成。
与现有技术相比，本发明所取得的有益效果是：
本发明根据冲击载荷作用仿真结果总结应力分布规律，建立梯度密度排布函数，对点阵胞元杆件尺寸进行梯度惩罚计算，使点阵胞元尺寸和相对密度沿高度方向呈梯度连续变化，材料分布趋向于应力较大的区域，结构变形和应力分布趋于均匀化，与均匀点阵结构相比，材料分布更加合理。通过本发明所提设计方法获得的功能梯度立方点阵结构材料，在实现结构轻量化的同时减小动载荷作用下均匀点阵结构的应力集中问题，使结构变形和能量耗散趋于均匀化，材料分布更加合理，提高结构的缓冲效果和吸能特性。
附图说明
图1为本发明的设计流程图；
图2本发明具体实施例的流程图；
图3为点阵胞元尺寸梯度曲线，（a）为应力分布曲线，（b）为梯度排布函数曲线，（c）为胞元杆件尺寸曲线。
图4为立方点阵胞元参数化设计与建模示意图，（a）为立方点阵胞元结构，（b）为点阵胞元排布模型，（c）为点阵胞元排布模型；
图5为功能梯度立方点阵结构材料示意图；（a）为长方体几何结构平压受力分析，（b）为均匀立方点阵结构模型，（c）为梯度立方点阵结构模型；
图6为两种点阵结构的冲击仿真结果，（a）为均匀立方点阵结构冲击仿真应力分布云图，（b）为梯度立方点阵结构冲击仿真应力分布云图。
具体实施方式
图1至图6示出了本发明一种参数化功能梯度立方点阵结构设计方法的具体实施例，如图1和图2所示，设计方法包括以下步骤：
步骤1：获取连续体几何结构所需承载的工况载荷参数，以及连续体几何结构的点阵胞元类型、单个胞元的主要设计参数
Figure 465176DEST_PATH_IMAGE024
，
Figure 767981DEST_PATH_IMAGE025
表示胞元包络立方体边长尺寸，
Figure 46516DEST_PATH_IMAGE026
表示胞元杆件截面边长尺寸；如图5(a)所示，规则长方体结构底端固定，上端承受冲击载荷(
Figure 460180DEST_PATH_IMAGE027
，
Figure 78243DEST_PATH_IMAGE028
)，设计域长
Figure 551950DEST_PATH_IMAGE029
，宽
Figure 317780DEST_PATH_IMAGE030
，高度
Figure 269556DEST_PATH_IMAGE031
，材料为铝材，弹性模量E=70GPa，给定结构减重要求60％左右，通过功能梯度点阵结构设计使材料分布更加合理，改善结构的冲击动力学性能。
步骤2：根据所述工况载荷参数，对所述连续体几何结构进行冲击性能仿真分析，得到冲击载荷作用下连续体几何结构内部应力分布规律，并根据应力分布规律建立密度梯度排布函数F(h)，h表示几何结构在应力梯度变化方向上的高度；
步骤2.1：根据冲击性能仿真分析得到冲击载荷作用下连续体几何结构内部应力分布规律，确定出应力分布函数σ(h)；
本实施例中采用有限元软件ANSYS进行均匀点阵结构冲击仿真分析，均匀点阵结构设计参数为设计域长
Figure 273284DEST_PATH_IMAGE032
，宽
Figure 652313DEST_PATH_IMAGE033
，高度
Figure 908369DEST_PATH_IMAGE034
，胞元尺寸
Figure 663836DEST_PATH_IMAGE035
，
Figure 990912DEST_PATH_IMAGE036
，
Figure 806421DEST_PATH_IMAGE037
，均匀点阵结构材料模型如图5（b）所示，获得内部应力分布规律，如图6(a)所示，最大应力集中在底端区域，应力沿高度方向呈梯度递减，提取应力分布参数绘制
Figure 750106DEST_PATH_IMAGE038
曲线如图3(a)所示，该曲线近似为一条曲线，按照二阶关系式进行拟合，确定应力分布函数
Figure 106001DEST_PATH_IMAGE039
表达式，得到的应力分布函数为：
Figure 553163DEST_PATH_IMAGE040
；应力分布函数可以是一阶、二阶或者高阶，本实施例中根据应力分布函数曲线选择合适的阶数进行拟合。
步骤2.2：根据所述应力分布函数
Figure 273994DEST_PATH_IMAGE041
建立初始密度梯度排布函数
Figure 501713DEST_PATH_IMAGE042
：
Figure 333403DEST_PATH_IMAGE043
步骤2.3：根据连续体几何结构的体积、相对密度以及减重比，确定在相同减重比的情况下，初始密度梯度排布函数
Figure 166230DEST_PATH_IMAGE044
的放大系数λ：
Figure 854700DEST_PATH_IMAGE045
其中，
Figure 775907DEST_PATH_IMAGE046
、
Figure 411288DEST_PATH_IMAGE013
、
Figure 629780DEST_PATH_IMAGE047
为最小应力、最大应力和第
Figure 692414DEST_PATH_IMAGE048
层可变应力；
Figure 629146DEST_PATH_IMAGE049
、
Figure 333797DEST_PATH_IMAGE050
、
Figure 610057DEST_PATH_IMAGE051
为胞元最小杆件尺寸、 最大或初始杆件尺寸和可变尺寸；
Figure 905909DEST_PATH_IMAGE052
为相对密度，
Figure 533200DEST_PATH_IMAGE053
为点阵胞元实体杆件体积，
Figure 775962DEST_PATH_IMAGE054
为包络体积计算得出；
Figure 906729DEST_PATH_IMAGE055
为减重比，
Figure 842324DEST_PATH_IMAGE056
为连续体几何结构分布层数。
本实施例中，选择点阵结构类型为体心立方点阵结构，如图4(a)所示，那么公式中的一些参数计算方式如下：胞元设计尺寸为点阵胞元边长a和杆件截面正六边形边长尺寸t。其中，胞元边长a根据设计域确定，取
Figure 956911DEST_PATH_IMAGE057
，
Figure 463417DEST_PATH_IMAGE058
，
Figure 979849DEST_PATH_IMAGE059
、
Figure 820766DEST_PATH_IMAGE060
为根据增材制造加工工艺确定的最小、最大杆件尺寸，根据现有研究结果，取
Figure 422648DEST_PATH_IMAGE061
、
Figure 803951DEST_PATH_IMAGE062
，即
Figure 174890DEST_PATH_IMAGE063
。对于体心立方点阵结构，参数计算如下：
直杆与斜杆截面边长关系：
Figure 186708DEST_PATH_IMAGE064
胞元实体杆件体积：
Figure 10307DEST_PATH_IMAGE065
胞元包络体积：
Figure 398563DEST_PATH_IMAGE066
相对密度：
Figure 624008DEST_PATH_IMAGE067
减重比：
Figure 806728DEST_PATH_IMAGE068
根据以上计算公式，由
Figure 383203DEST_PATH_IMAGE069
，计算相对密度
Figure 309571DEST_PATH_IMAGE070
，取
Figure 389522DEST_PATH_IMAGE071
，
Figure 8722DEST_PATH_IMAGE072
，则
Figure 809844DEST_PATH_IMAGE073
。
步骤2.4：输出密度梯度排布函数
Figure 539902DEST_PATH_IMAGE074
：
Figure 474360DEST_PATH_IMAGE075
本实施例中，均匀点阵结构仿真应力曲线如图3（a）所示，根据上述公式计算得到
Figure 264462DEST_PATH_IMAGE076
，
Figure 284370DEST_PATH_IMAGE077
曲线如图3(b)所示，密度梯度函数
Figure 552541DEST_PATH_IMAGE078
：
Figure 607084DEST_PATH_IMAGE079
。
所建立的密度梯度排布函数具有一定的代表性，进行一次或一组仿真分析，可以表征同一类振动或冲击作用应力分布特征，与基于拓扑优化的梯度结构设计方法相比，减少部分重复性工作，在设计方法上具有一定优势。
步骤3：对连续体几何结构进行分层，根据所述密度梯度排布函数对每层的点阵胞元杆件尺寸进行梯度惩罚计算，得到各层点阵胞元杆件尺寸，如图3（c）所示；
对连续体几何结构进行分层是指沿连续体几何结构的梯度变化方向，将点阵胞元排布层的一层或多层分成一个连续体几何结构层。本实施例中是根据点阵胞元的排布层的一层进行分层的。根据连续体给定连续体几何尺寸
Figure 568087DEST_PATH_IMAGE080
和立方点阵胞元包络尺寸
Figure 872029DEST_PATH_IMAGE081
确定点阵胞元排布层数和数量，则x方向胞元排布数目:
Figure 943891DEST_PATH_IMAGE082
；y方向胞元排布数目:
Figure 852941DEST_PATH_IMAGE083
；z方向胞元排布层数:
Figure 719266DEST_PATH_IMAGE084
；
Figure 979346DEST_PATH_IMAGE085
、
Figure 589319DEST_PATH_IMAGE086
、
Figure 618454DEST_PATH_IMAGE087
分别表示连续体几何结构的长、宽、高；
本实施例中，连续体几何结构尺寸为20mm×20mm×60mm，立方点阵胞元包络尺寸为5mm×5mm×5mm，假定使用的立方点阵结构为体心立方点阵结构，如图4(a)所示，胞元建模参考基准如图4(b)所示，点阵胞元排布模型如图4(c)所示，初始胞元设计尺寸为点阵胞元边长
Figure 455348DEST_PATH_IMAGE088
和杆件截面正六边形边长尺寸
Figure 202724DEST_PATH_IMAGE089
；x方向胞元排布数目:
Figure 616388DEST_PATH_IMAGE090
；y方向胞元排布数目:
Figure 500030DEST_PATH_IMAGE091
；z方向胞元排布数目:
Figure 973737DEST_PATH_IMAGE092
；
根据密度梯度排布函数对每层的点阵胞元杆件尺寸进行梯度惩罚计算，梯度惩罚计算得到各层点阵胞元杆件尺寸的方法是：第k层点阵点阵胞元的杆件尺寸
Figure 942830DEST_PATH_IMAGE093
，
Figure 894605DEST_PATH_IMAGE094
表示初始层胞元杆件尺寸，
Figure 632754DEST_PATH_IMAGE095
表示几何结构沿应力梯度变化方向第k层点阵对应的高度上梯度密度排布函数值。
本实施例中所述的梯度立方点阵结构是指点阵胞元杆件尺寸和相对密度沿高度方向呈梯度连续变化，每个胞元包络空间边长尺寸
Figure 277362DEST_PATH_IMAGE096
相同，胞元杆件截面边长尺寸
Figure 733751DEST_PATH_IMAGE097
逐层梯度递减，整体结构呈现梯度分层排布特征。本实施例中，立方点阵结构还可以为边立方、面心立方、OctetFrame立方等，也可以是其它类似的具有包络空间、胞元结构形状规则的立方点阵结构，在不改变点阵胞元杆件连接方式的前提下，杆件截面形状可以是正方形、圆形、正六边形等形状，便于在设计域内进行点阵周期排布和密度映射求解，得到较为规整的三维实体结构。
步骤4：根据所述各层点阵胞元杆件尺寸，使用三维软件建模程序进行点阵胞元的生成和循环装配，得到功能梯度点阵结构材料模型。进行点阵胞元循环装配的方法是：
步骤4.1：给初始点阵胞元尺寸赋值
Figure 223639DEST_PATH_IMAGE098
，生成初始立方点阵胞元元件；
步骤4.2：设定装配初始原点p、装配基准，并和点阵胞元元件基准相匹配，
步骤4.3：根据步骤3中得到的各层点阵胞元尺寸，生成点阵胞元元件，从第一层第一个点阵胞元元件开始装配，直到完成所有层上点阵胞元的装配，生成梯度点阵结构模型，如图5（c）所示。
本实施例中，首先进行第1层点阵胞元元件重生和装配；取
Figure 81873DEST_PATH_IMAGE099
，
Figure 897382DEST_PATH_IMAGE100
，完成第一层第一个点阵胞元元件重生和装配，根据
Figure 841068DEST_PATH_IMAGE101
，进行第一层其它点阵胞元元件重生和装配；然后进行第k+1层其它点阵胞元元件重生和装配，取
Figure 134646DEST_PATH_IMAGE102
，
Figure 647054DEST_PATH_IMAGE103
，保持
Figure 836727DEST_PATH_IMAGE104
不变，根据
Figure 267708DEST_PATH_IMAGE105
，重新生成点阵胞元元件，完成第
Figure 161715DEST_PATH_IMAGE106
层第一个点阵胞元元件装配，根据
Figure 463383DEST_PATH_IMAGE107
，完成对第
Figure 620695DEST_PATH_IMAGE108
层其它点阵胞元元件装配；重复该步骤，直到层数
Figure 538973DEST_PATH_IMAGE109
，完成所有点阵胞元元件重生和装配，生成梯度点阵结构模型，
Figure 971091DEST_PATH_IMAGE110
表示点阵胞元在z方向的胞元层数。
本实施例中，三维建模软件可以是Creo，也可以是Solidworks、UG等其它具有相同建模功能的三维设计软件，进行点阵胞元循环装配可以通过Creo用户自定义辅助建模程序来实现，也可以基于Solidworks、UG等其它建模软件来实现，如果不考虑工作量和建模效率，也可以采用软件现有的“阵列”、“装配”等菜单工具人工装配实现。采用参数化建模方法，基于Creo用户自定义辅助建模程序进行点阵胞元重生和循环装配，以获得功能梯度点阵结构模型，便于进行结构改进和参数更新，二次建模操作简单，提高了建模效率，可适用于较复杂的大尺度几何零件设计。
步骤5：对所述功能梯度点阵结构材料模型进行冲击仿真分析，判断所述功能梯度点阵结构模型是否满足设计要求，若满足，则根据该点阵结构模型得到点阵结构材料，若不满足，则进行结构优化设计直到满足设计要求。本实施例中得到结构优化设计是指通过改变梯度密度排布函数F(h)来实现。本梯度类型可以是沿某一方向的正梯度、负梯度或正负混合梯度。本发明所提设计方法采用理论设计-仿真分析-实验验证相结合的方法，仿真和实验相互验证，不断进行结构优化设计，设计方法较为完善。
对功能梯度立方点阵结构模型进行相同冲击载荷作用仿真分析，获得应力分布云图如图6（b）所示，与均匀点阵结构应力云图6（a）相比，减小了冲击载荷作用的最大应力，
Figure 658424DEST_PATH_IMAGE111
，最大应力降低20％左右，同时，改变了最大应力集中在底端区域的分布状态，使结构变形和应力分布趋于均匀化。通过对比分析，证明了功能梯度立方点阵结构的材料分布更加合理，提高了结构的抗冲击和缓冲吸能特性。
通过本设计方法得到的功能梯度立方点阵材料，其空间结构复杂，但新型加工方法增材制造技术（3D打印）的日渐成熟保证了该复杂点阵结构的可制造性。通过本设计方法最终获得的功能梯度立方点阵结构材料，在实现结构轻量化的同时减小动载荷作用下均匀点阵结构的应力集中问题，使结构变形和能量耗散趋于均匀化，材料分布更加合理，提高结构的缓冲效果和吸能特性。
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。

Claims (7)

1.一种参数化功能梯度立方点阵结构设计方法，其特征在于：包括以下步骤：
步骤1：获取连续体几何结构所需承载的工况载荷参数，以及连续体几何结构的点阵胞元类型、单个胞元的主要设计参数
Figure 583224DEST_PATH_IMAGE001
，
Figure 698948DEST_PATH_IMAGE002
表示胞元包络立方体边长尺寸，
Figure 866624DEST_PATH_IMAGE003
表示胞元杆件截面边长尺寸；
步骤2：根据所述工况载荷参数，对所述连续体几何结构进行冲击性能仿真分析，得到冲击载荷作用下连续体几何结构内部应力分布规律，并根据应力分布规律建立密度梯度排布函数
Figure 409601DEST_PATH_IMAGE004
，
Figure 233200DEST_PATH_IMAGE005
表示几何结构在应力梯度变化方向上的高度；
步骤3：对连续体几何结构进行分层，根据所述密度梯度排布函数对每层的点阵胞元杆件尺寸进行梯度惩罚计算，得到各层点阵胞元杆件尺寸；
所述梯度惩罚计算的方法是指：第
Figure 887035DEST_PATH_IMAGE006
层点阵点阵胞元的杆件尺寸
Figure 174797DEST_PATH_IMAGE007
，
Figure 91938DEST_PATH_IMAGE008
表示初始层胞元杆件尺寸，
Figure 402833DEST_PATH_IMAGE009
表示几何结构沿应力梯度变化方向第
Figure 125939DEST_PATH_IMAGE006
层点阵对应的高度上梯度密度排布函数值；
步骤4：根据所述各层点阵胞元杆件尺寸，使用三维软件建模程序进行点阵胞元的生成和循环装配，得到功能梯度点阵结构材料模型。
2.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法，其特征在于：还包括步骤5：对所述功能梯度点阵结构材料模型进行冲击仿真分析，判断所述功能梯度点阵结构模型是否满足设计要求，若满足，则根据该点阵结构模型得到点阵结构材料，若不满足，则进行结构优化设计直到满足设计要求。
3.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法，其特征在于：对连续体几何结构进行分层是指沿连续体几何结构的梯度变化方向，将点阵胞元排布层的一层或多层分成一个连续体几何结构层。
4.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法，其特征在于：步骤4中所述进行点阵胞元循环装配的方法是：
步骤4.1：给初始点阵胞元尺寸赋值(a,t 0 )，生成初始立方点阵胞元元件；
步骤4.2：设定装配初始原点p、装配基准，并和点阵胞元元件基准相匹配，
步骤4.3：根据步骤3中得到的各层点阵胞元尺寸，生成点阵胞元元件，从第一层第一个点阵胞元元件开始装配，直到完成所有层上点阵胞元的装配，生成梯度点阵结构模型。
5.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法，其特征在于：步骤2中所述根据应力分布规律建立密度梯度排布函数F(h)方法是：
步骤2.1：根据冲击性能仿真分析得到冲击载荷作用下连续体几何结构内部应力分布规律，确定出应力分布函数
Figure 674732DEST_PATH_IMAGE010
；
步骤2.2：根据所述应力分布函数
Figure 825090DEST_PATH_IMAGE011
建立初始密度梯度排布函数
Figure 357703DEST_PATH_IMAGE012
：
Figure 87761DEST_PATH_IMAGE013
步骤2.3：根据连续体几何结构的体积、相对密度以及减重比，确定在相同减重比的情况下，初始密度梯度排布函数
Figure 84536DEST_PATH_IMAGE014
的放大系数
Figure 609058DEST_PATH_IMAGE015
：
Figure 628967DEST_PATH_IMAGE016
步骤2.4：输出密度梯度排布函数
Figure 162717DEST_PATH_IMAGE017
：
Figure 13998DEST_PATH_IMAGE018
其中，
Figure 709422DEST_PATH_IMAGE019
、
Figure 234204DEST_PATH_IMAGE020
、
Figure 306066DEST_PATH_IMAGE021
为最小应力、最大应力和第
Figure 215116DEST_PATH_IMAGE022
层可变应力；
Figure 347020DEST_PATH_IMAGE023
、
Figure 872679DEST_PATH_IMAGE024
、
Figure 748231DEST_PATH_IMAGE025
为胞元最小杆件尺寸、最大或初始杆件尺寸和可变尺寸；
Figure 511788DEST_PATH_IMAGE026
为相对密度，
Figure 611331DEST_PATH_IMAGE027
为点阵胞元实体杆件体积，
Figure 93128DEST_PATH_IMAGE028
为包络体积计算得出；
Figure 241212DEST_PATH_IMAGE029
为减重比，
Figure 124855DEST_PATH_IMAGE030
为连续体几何结构分布层数。
6.根据权利要求2所述的一种参数化功能梯度立方点阵结构设计方法，其特征在于：步骤5中所述结构优化设计通过改变梯度密度排布函数F(h)来实现。
7.一种参数化功能梯度立方点阵结构材料，其特征在于：所述材料使用权利要求1至6中任一项权利要求所述方法制作而成。
CN202010069303.4A 2020-01-21 2020-01-21 一种参数化功能梯度立方点阵结构设计方法及材料 Active CN110929358B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010069303.4A CN110929358B (zh) 2020-01-21 2020-01-21 一种参数化功能梯度立方点阵结构设计方法及材料

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010069303.4A CN110929358B (zh) 2020-01-21 2020-01-21 一种参数化功能梯度立方点阵结构设计方法及材料

Publications (2)

Publication Number Publication Date
CN110929358A CN110929358A (zh) 2020-03-27
CN110929358B true CN110929358B (zh) 2020-05-08

Family

ID=69854422

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010069303.4A Active CN110929358B (zh) 2020-01-21 2020-01-21 一种参数化功能梯度立方点阵结构设计方法及材料

Country Status (1)

Country Link
CN (1) CN110929358B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111680441B (zh) * 2020-06-08 2023-03-28 南京理工大学 一种适用于热力工况的梯度点阵夹芯板结构
CN111967089B (zh) * 2020-08-04 2021-08-03 南京理工大学 一种基于渐变胞元尺寸的异型三维点阵结构隔热性能优化设计方法
CN112356523B (zh) 2020-08-29 2021-12-07 南京航空航天大学 基于可编程刚度的手性胞元构建的梯度点阵吸能结构及其3d打印方法
CN112249509B (zh) * 2020-09-01 2022-08-02 哈尔滨工业大学（深圳） 吸能结构及吸能缓冲装置
CN112268092B (zh) * 2020-09-30 2022-05-27 重庆大学 一种点阵流变智能结构及其制作方法
CN112287491B (zh) * 2020-12-28 2021-03-12 中国人民解放军国防科技大学 一种复合型点阵材料及其设计方法
CN112800553B (zh) * 2021-01-04 2024-01-30 中国科学院力学研究所 多级可控渐进吸能点阵结构
CN113722942B (zh) * 2021-07-08 2022-07-12 北京理工大学 一种考虑3d打印几何缺陷的有限元计算模型
CN114055080A (zh) * 2021-11-12 2022-02-18 吉林大学 一种基于梯度多胞点阵结构的加工方法
CN114379488B (zh) * 2022-01-24 2024-03-08 中南大学 一种仿生梯度多级管状结构
CN114560068A (zh) * 2022-01-26 2022-05-31 重庆大学 一种基于梯度点阵结构的轻质电动舵机抗冲击底座及外壳
CN114636360B (zh) * 2022-03-23 2023-10-27 中国人民解放军海军工程大学 五模冲击隐身复合点阵环状结构及其参数优化方法
CN115958203B (zh) * 2022-10-19 2023-06-20 中国机械总院集团沈阳铸造研究所有限公司 一种具有减振特性的变密度点阵金属
CN116127655B (zh) * 2023-04-17 2023-07-04 之江实验室 一种缓冲组件的制造方法、装置、存储介质及电子设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105912803A (zh) * 2016-04-28 2016-08-31 华南理工大学 一种基于增材制造的产品轻量化设计方法
CN106650148A (zh) * 2016-12-30 2017-05-10 北京航空航天大学 一种位移和应力混合约束下的连续体结构非概率可靠性拓扑优化方法
CN109766656A (zh) * 2019-01-25 2019-05-17 北京航空航天大学 一种基于拓扑优化的梯度点阵结构设计方法
CN109858133A (zh) * 2019-01-25 2019-06-07 北京航空航天大学 一种基于应力映射的点阵结构设计与优化方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9793829B2 (en) * 2013-09-25 2017-10-17 Prime Photonics, Lc Magneto-thermoelectric generator for energy harvesting
CN106694884B (zh) * 2016-12-29 2020-02-21 西安铂力特增材技术股份有限公司 一种具有梯度功能性的镂空点阵夹层及其制造方法
CN107742014B (zh) * 2017-09-29 2020-10-30 北京空间飞行器总体设计部 基于增材制造的相变储能装置点阵夹层结构的设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105912803A (zh) * 2016-04-28 2016-08-31 华南理工大学 一种基于增材制造的产品轻量化设计方法
CN106650148A (zh) * 2016-12-30 2017-05-10 北京航空航天大学 一种位移和应力混合约束下的连续体结构非概率可靠性拓扑优化方法
CN109766656A (zh) * 2019-01-25 2019-05-17 北京航空航天大学 一种基于拓扑优化的梯度点阵结构设计方法
CN109858133A (zh) * 2019-01-25 2019-06-07 北京航空航天大学 一种基于应力映射的点阵结构设计与优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Optimal Design and Modeling of 3D Variable-Density Lattice Structures;Xin Jin 等;《2017 8th International Conference on Mechanical and Aerospace Engineering》;20170830;全文 *
Optimal Design and Modeling of Variable-density Triangular Honeycomb Structures;Xin Jin 等;《2017 8th International Conference on Mechanical and Intelligent Manufacturing Technologies》;20170830;全文 *
基于三维模型的产品模块化配置方法;张萌 等;《计算机集成制造系统》;20121130;全文 *

Also Published As

Publication number Publication date
CN110929358A (zh) 2020-03-27

Similar Documents

Publication Publication Date Title
CN110929358B (zh) 一种参数化功能梯度立方点阵结构设计方法及材料
CN112287491B (zh) 一种复合型点阵材料及其设计方法
CN112214856A (zh) 面向整体结构的精密机床刚度优化设计方法
CN106960096A (zh) 三维负刚度蜂窝结构及其制备方法
CN110837690B (zh) 蜂窝结构非线性本构关系的建立方法、介质和设备
CN111639399B (zh) 考虑吸能刚度的点阵筛选和相对强度点阵填充方法
CN102184298A (zh) 金属体积塑性成形中有限元分析刚度矩阵存储与生成方法
Zhou et al. Design optimization of a NPR structure based on HAM optimization method
CN111382536B (zh) 一种约束阻尼结构的协同优化设计方法
CN112069623A (zh) 一种航空发动机风扇叶片点阵填充方法
CN109766638B (zh) 一种燃油喷嘴及其轻量化设计方法
CN112549670A (zh) 一种基于3d打印的变刚度夹芯复合材料结构及其成型方法
Zhang et al. Crashworthiness optimization of a vertex fractal hexagonal structure
CN112149260B (zh) 一种三维抗冲击负泊松比结构的设计方法
CN116522499A (zh) 一种基于薄壁–点阵填充结构车身前端结构设计方法
Gu et al. Novel gradient design and simulation of Voronoi structures
CN111639391B (zh) 一种作业车车臂截面参数选用方法
CN111723500B (zh) 基于米字形单元孪晶型微桁架结构的复合吸能结构及其3d打印方法
CN114771797A (zh) 粘弹性材料填充曲杆化点阵结构舵机减振底座及制造方法
CN102774042A (zh) 一种蜂窝加强板及其制作方法
CN111723501B (zh) 基于正方形单元多晶型微桁架结构的复合吸能结构及其3d打印方法
Cheng et al. Finite element analysis of the clamping stress of photovoltaic silicon brick in diamond wire sawing
CN218762647U (zh) 一种轻质量高频率的多口型加筋板
CN111985120B (zh) 基于米字形单元多晶型微桁架结构的复合吸能结构及其3d打印方法
Dokainish et al. Vibrations of orthotropic parallelogramic plates with variable thickness

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