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

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

Info

Publication number
CN110929358A
CN110929358A CN202010069303.4A CN202010069303A CN110929358A CN 110929358 A CN110929358 A CN 110929358A CN 202010069303 A CN202010069303 A CN 202010069303A CN 110929358 A CN110929358 A CN 110929358A
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.)
Granted
Application number
CN202010069303.4A
Other languages
English (en)
Other versions
CN110929358B (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 (8)

1.一种参数化功能梯度立方点阵结构设计方法,其特征在于:包括以下步骤:
步骤1:获取连续体几何结构所需承载的工况载荷参数,以及连续体几何结构的点阵胞元类型、单个胞元的主要设计参数
Figure 477617DEST_PATH_IMAGE001
Figure 27547DEST_PATH_IMAGE002
表示胞元包络立方体边长尺寸,
Figure 985881DEST_PATH_IMAGE003
表示胞元杆件截面边长尺寸;
步骤2:根据所述工况载荷参数,对所述连续体几何结构进行冲击性能仿真分析,得到冲击载荷作用下连续体几何结构内部应力分布规律,并根据应力分布规律建立密度梯度排布函数
Figure 545038DEST_PATH_IMAGE004
Figure 257780DEST_PATH_IMAGE005
表示几何结构在应力梯度变化方向上的高度;
步骤3:对连续体几何结构进行分层,根据所述密度梯度排布函数对每层的点阵胞元杆件尺寸进行梯度惩罚计算,得到各层点阵胞元杆件尺寸;
步骤4:根据所述各层点阵胞元杆件尺寸,使用三维软件建模程序进行点阵胞元的生成和循环装配,得到功能梯度点阵结构材料模型。
2.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法,其特征在于:还包括步骤5:对所述功能梯度点阵结构材料模型进行冲击仿真分析,判断所述功能梯度点阵结构模型是否满足设计要求,若满足,则根据该点阵结构模型得到点阵结构材料,若不满足,则进行结构优化设计直到满足设计要求。
3.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法,其特征在于:对连续体几何结构进行分层是指沿连续体几何结构的梯度变化方向,将点阵胞元排布层的一层或多层分成一个连续体几何结构层。
4.根据权利要求3所述的一种参数化功能梯度立方点阵结构设计方法,其特征在于:步骤3中根据所述密度梯度排布函数对点阵胞元尺寸进行梯度惩罚计算的方法是指:第
Figure 244190DEST_PATH_IMAGE006
层点阵点阵胞元的杆件尺寸
Figure 409592DEST_PATH_IMAGE007
Figure 506861DEST_PATH_IMAGE008
表示初始层胞元杆件尺寸,
Figure 339688DEST_PATH_IMAGE009
表示几何结构沿应力梯度变化方向第
Figure 231421DEST_PATH_IMAGE010
层点阵对应的高度上梯度密度排布函数值。
5.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法,其特征在于:步骤4中所述进行点阵胞元循环装配的方法是:
步骤4.1:给初始点阵胞元尺寸赋值(a,t 0 ),生成初始立方点阵胞元元件;
步骤4.2:设定装配初始原点p、装配基准,并和点阵胞元元件基准相匹配,
步骤4.3:根据步骤3中得到的各层点阵胞元尺寸,生成点阵胞元元件,从第一层第一个点阵胞元元件开始装配,直到完成所有层上点阵胞元的装配,生成梯度点阵结构模型。
6.根据权利要求1所述的一种参数化功能梯度立方点阵结构设计方法,其特征在于:步骤2中所述根据应力分布规律建立密度梯度排布函数F(h)方法是:
步骤2.1:根据冲击性能仿真分析得到冲击载荷作用下连续体几何结构内部应力分布规律,确定出应力分布函数σ(h);
步骤2.2:根据所述应力分布函数σ(h)建立初始密度梯度排布函数
Figure 884119DEST_PATH_IMAGE011
Figure 785079DEST_PATH_IMAGE012
步骤2.3:根据连续体几何结构的体积、相对密度以及减重比,确定在相同减重比的情况下,初始密度梯度排布函数
Figure 269150DEST_PATH_IMAGE013
的放大系数
Figure 597363DEST_PATH_IMAGE014
Figure 737357DEST_PATH_IMAGE015
步骤2.4:输出密度梯度排布函数
Figure 176429DEST_PATH_IMAGE016
Figure 718269DEST_PATH_IMAGE017
其中,
Figure 220313DEST_PATH_IMAGE018
Figure 847603DEST_PATH_IMAGE019
Figure 90366DEST_PATH_IMAGE020
为最小应力、最大应力和第
Figure 221133DEST_PATH_IMAGE021
层可变应力;
Figure 156728DEST_PATH_IMAGE022
Figure 271314DEST_PATH_IMAGE023
Figure 52188DEST_PATH_IMAGE024
为胞元最小杆件尺寸、最大或初始杆件尺寸和可变尺寸;
Figure 303041DEST_PATH_IMAGE025
为相对密度,
Figure 143958DEST_PATH_IMAGE026
为点阵胞元实体杆件体积,
Figure 542579DEST_PATH_IMAGE027
为包络体积计算得出;
Figure 861564DEST_PATH_IMAGE028
为减重比,
Figure 232503DEST_PATH_IMAGE029
为连续体几何结构分布层数。
7.根据权利要求2所述的一种参数化功能梯度立方点阵结构设计方法,其特征在于:步骤5中所述结构优化设计通过改变梯度密度排布函数F(h)来实现。
8.一种参数化功能梯度立方点阵结构材料,其特征在于:所述材料使用权利要求1至7中任一项权利要求所述方法制作而成。
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 true CN110929358A (zh) 2020-03-27
CN110929358B 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)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111680441A (zh) * 2020-06-08 2020-09-18 南京理工大学 一种适用于热力工况的梯度点阵夹芯板结构
CN111967089A (zh) * 2020-08-04 2020-11-20 南京理工大学 一种基于渐变胞元尺寸的异型三维点阵结构隔热性能优化设计方法
CN112249509A (zh) * 2020-09-01 2021-01-22 哈尔滨工业大学(深圳) 吸能结构及吸能缓冲装置
CN112268092A (zh) * 2020-09-30 2021-01-26 重庆大学 一种点阵流变智能结构及其制作方法
CN112287491A (zh) * 2020-12-28 2021-01-29 中国人民解放军国防科技大学 一种复合型点阵材料及其设计方法
CN112800553A (zh) * 2021-01-04 2021-05-14 中国科学院力学研究所 多级可控渐进吸能点阵结构
CN113722942A (zh) * 2021-07-08 2021-11-30 北京理工大学 一种考虑3d打印几何缺陷的有限元计算模型
CN114055080A (zh) * 2021-11-12 2022-02-18 吉林大学 一种基于梯度多胞点阵结构的加工方法
WO2022042580A1 (zh) * 2020-08-29 2022-03-03 南京航空航天大学 梯度点阵吸能结构及其可编程刚度的手性胞元、3d打印方法
CN114239145A (zh) * 2021-12-16 2022-03-25 北京理工大学 一种变形稳定的高强度轻质防护点阵结构材料
CN114386117A (zh) * 2021-12-31 2022-04-22 西北工业大学 环向梯度泡沫材料的三维细观模型的建模方法及其系统
CN114379488A (zh) * 2022-01-24 2022-04-22 中南大学 一种仿生梯度多级管状结构
CN114560068A (zh) * 2022-01-26 2022-05-31 重庆大学 一种基于梯度点阵结构的轻质电动舵机抗冲击底座及外壳
CN114636360A (zh) * 2022-03-23 2022-06-17 中国人民解放军海军工程大学 五模冲击隐身复合点阵环状结构及其参数优化方法
CN115958203A (zh) * 2022-10-19 2023-04-14 沈阳铸造研究所有限公司 一种具有减振特性的变密度点阵金属及其增材制造方法
CN116127655A (zh) * 2023-04-17 2023-05-16 之江实验室 一种缓冲组件的制造方法、装置、存储介质及电子设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150083196A1 (en) * 2013-09-25 2015-03-26 Prime Photonics, Lc Magneto-thermoelectric generator for energy harvesting
CN105912803A (zh) * 2016-04-28 2016-08-31 华南理工大学 一种基于增材制造的产品轻量化设计方法
CN106650148A (zh) * 2016-12-30 2017-05-10 北京航空航天大学 一种位移和应力混合约束下的连续体结构非概率可靠性拓扑优化方法
CN106694884A (zh) * 2016-12-29 2017-05-24 西安铂力特激光成形技术有限公司 一种具有梯度功能性的镂空点阵夹层及其制造方法
CN107742014A (zh) * 2017-09-29 2018-02-27 北京空间飞行器总体设计部 基于增材制造的相变储能装置点阵夹层结构的设计方法
CN109766656A (zh) * 2019-01-25 2019-05-17 北京航空航天大学 一种基于拓扑优化的梯度点阵结构设计方法
CN109858133A (zh) * 2019-01-25 2019-06-07 北京航空航天大学 一种基于应力映射的点阵结构设计与优化方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150083196A1 (en) * 2013-09-25 2015-03-26 Prime Photonics, Lc Magneto-thermoelectric generator for energy harvesting
CN105912803A (zh) * 2016-04-28 2016-08-31 华南理工大学 一种基于增材制造的产品轻量化设计方法
CN106694884A (zh) * 2016-12-29 2017-05-24 西安铂力特激光成形技术有限公司 一种具有梯度功能性的镂空点阵夹层及其制造方法
CN106650148A (zh) * 2016-12-30 2017-05-10 北京航空航天大学 一种位移和应力混合约束下的连续体结构非概率可靠性拓扑优化方法
CN107742014A (zh) * 2017-09-29 2018-02-27 北京空间飞行器总体设计部 基于增材制造的相变储能装置点阵夹层结构的设计方法
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
XIN JIN 等: "Optimal Design and Modeling of 3D Variable-Density Lattice Structures", 《2017 8TH INTERNATIONAL CONFERENCE ON MECHANICAL AND AEROSPACE ENGINEERING》 *
XIN JIN 等: "Optimal Design and Modeling of Variable-density Triangular Honeycomb Structures", 《2017 8TH INTERNATIONAL CONFERENCE ON MECHANICAL AND INTELLIGENT MANUFACTURING TECHNOLOGIES》 *
张萌 等: "基于三维模型的产品模块化配置方法", 《计算机集成制造系统》 *

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111680441B (zh) * 2020-06-08 2023-03-28 南京理工大学 一种适用于热力工况的梯度点阵夹芯板结构
CN111680441A (zh) * 2020-06-08 2020-09-18 南京理工大学 一种适用于热力工况的梯度点阵夹芯板结构
CN111967089A (zh) * 2020-08-04 2020-11-20 南京理工大学 一种基于渐变胞元尺寸的异型三维点阵结构隔热性能优化设计方法
WO2022042580A1 (zh) * 2020-08-29 2022-03-03 南京航空航天大学 梯度点阵吸能结构及其可编程刚度的手性胞元、3d打印方法
US11731389B2 (en) 2020-08-29 2023-08-22 Nanjing University Of Aeronautics And Astronautics Graded lattice energy-absorbing structure, chiral cell thereof having programmable stiffness, and 3D printing method
CN112249509A (zh) * 2020-09-01 2021-01-22 哈尔滨工业大学(深圳) 吸能结构及吸能缓冲装置
CN112249509B (zh) * 2020-09-01 2022-08-02 哈尔滨工业大学(深圳) 吸能结构及吸能缓冲装置
CN112268092A (zh) * 2020-09-30 2021-01-26 重庆大学 一种点阵流变智能结构及其制作方法
CN112287491A (zh) * 2020-12-28 2021-01-29 中国人民解放军国防科技大学 一种复合型点阵材料及其设计方法
CN112287491B (zh) * 2020-12-28 2021-03-12 中国人民解放军国防科技大学 一种复合型点阵材料及其设计方法
CN112800553A (zh) * 2021-01-04 2021-05-14 中国科学院力学研究所 多级可控渐进吸能点阵结构
CN112800553B (zh) * 2021-01-04 2024-01-30 中国科学院力学研究所 多级可控渐进吸能点阵结构
CN113722942A (zh) * 2021-07-08 2021-11-30 北京理工大学 一种考虑3d打印几何缺陷的有限元计算模型
CN114055080A (zh) * 2021-11-12 2022-02-18 吉林大学 一种基于梯度多胞点阵结构的加工方法
CN114055080B (zh) * 2021-11-12 2024-06-21 吉林大学 一种基于梯度多胞点阵结构的加工方法
CN114239145A (zh) * 2021-12-16 2022-03-25 北京理工大学 一种变形稳定的高强度轻质防护点阵结构材料
CN114239145B (zh) * 2021-12-16 2024-08-02 北京理工大学 一种变形稳定的高强度轻质防护点阵结构材料
CN114386117A (zh) * 2021-12-31 2022-04-22 西北工业大学 环向梯度泡沫材料的三维细观模型的建模方法及其系统
CN114386117B (zh) * 2021-12-31 2024-09-06 西北工业大学 环向梯度泡沫材料的三维细观模型的建模方法及其系统
CN114379488A (zh) * 2022-01-24 2022-04-22 中南大学 一种仿生梯度多级管状结构
CN114379488B (zh) * 2022-01-24 2024-03-08 中南大学 一种仿生梯度多级管状结构
CN114560068A (zh) * 2022-01-26 2022-05-31 重庆大学 一种基于梯度点阵结构的轻质电动舵机抗冲击底座及外壳
CN114636360B (zh) * 2022-03-23 2023-10-27 中国人民解放军海军工程大学 五模冲击隐身复合点阵环状结构及其参数优化方法
CN114636360A (zh) * 2022-03-23 2022-06-17 中国人民解放军海军工程大学 五模冲击隐身复合点阵环状结构及其参数优化方法
CN115958203A (zh) * 2022-10-19 2023-04-14 沈阳铸造研究所有限公司 一种具有减振特性的变密度点阵金属及其增材制造方法
CN116127655A (zh) * 2023-04-17 2023-05-16 之江实验室 一种缓冲组件的制造方法、装置、存储介质及电子设备
WO2024078165A1 (zh) * 2023-04-17 2024-04-18 之江实验室 一种缓冲组件的制造方法、装置、存储介质及电子设备

Also Published As

Publication number Publication date
CN110929358B (zh) 2020-05-08

Similar Documents

Publication Publication Date Title
CN110929358B (zh) 一种参数化功能梯度立方点阵结构设计方法及材料
CN112287491B (zh) 一种复合型点阵材料及其设计方法
CN109359321B (zh) 一种拓扑优化与点阵结构结合的轻量化连接铰优化方法
CN109918712B (zh) 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法
CN112140647A (zh) 一种具有负泊松比特性的抗冲击、高吸能的点阵夹芯结构
CN112492864B (zh) 一种可控电磁屏蔽构件及其制备方法
CN111709085A (zh) 一种约束阻尼薄板结构拓扑优化设计方法
CN110837690B (zh) 蜂窝结构非线性本构关系的建立方法、介质和设备
CN112069623A (zh) 一种航空发动机风扇叶片点阵填充方法
Zhou et al. Design optimization of a NPR structure based on HAM optimization method
CN111382536B (zh) 一种约束阻尼结构的协同优化设计方法
CN112283277A (zh) 结构承载与振动控制一体化的复合点阵结构及其制备方法
CN116522499A (zh) 一种基于薄壁–点阵填充结构车身前端结构设计方法
CN112149260B (zh) 一种三维抗冲击负泊松比结构的设计方法
Roland et al. Numerical simulation and comparison of a real Al–Si alloy with virtually generated alloys
McConaha Graded lattice structure density optimization for additive manufacturing
CN111639391B (zh) 一种作业车车臂截面参数选用方法
CN111703067B (zh) 基于正方形单元孪晶型微桁架结构的复合吸能结构及其3d打印方法
CN113987676A (zh) 一种约束阻尼结构的材料结构一体化优化方法
CN114117635A (zh) 一种多维多刚度轻质吸能结构的设计方法
CN111723500B (zh) 基于米字形单元孪晶型微桁架结构的复合吸能结构及其3d打印方法
CN201112434Y (zh) 大容量锂离子动力电池壳体筋板结构
Zhao et al. A Parametric Design Method of Functionally Graded Lattice Structure Based on the Stress Distribution
CN111985120B (zh) 基于米字形单元多晶型微桁架结构的复合吸能结构及其3d打印方法
CN218762647U (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
GR01 Patent grant
GR01 Patent grant