CN111680441B - 一种适用于热力工况的梯度点阵夹芯板结构 - Google Patents
一种适用于热力工况的梯度点阵夹芯板结构 Download PDFInfo
- Publication number
- CN111680441B CN111680441B CN202010510048.2A CN202010510048A CN111680441B CN 111680441 B CN111680441 B CN 111680441B CN 202010510048 A CN202010510048 A CN 202010510048A CN 111680441 B CN111680441 B CN 111680441B
- Authority
- CN
- China
- Prior art keywords
- gradient
- lattice
- value
- stress
- node
- 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
Links
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Laminated Bodies (AREA)
Abstract
本发明提出了一种适用于热力工况的梯度点阵夹芯板结构,包括两层金属蒙皮以及设置在两层金属蒙皮中间的点阵芯子,所述点阵芯子的表面由三周期极小曲面函数表示。本发明增强了结构的整体承力、隔热性能,可广泛用于航空、航天等的热防护结构上。
Description
技术领域
本发明属于三维点阵夹芯结构优化设计领域,具体为一种适用于热力工况的梯度点阵夹芯板结构。
背景技术
航天器在大气层中高速飞行时,表面会与空气摩擦而急剧升温,其温度可达几百甚至上千度,在这样一个高温环境中,航天器除本身承受的外力载荷以外,还需要承受温升带来的热应力。因此,除提高***隔热性能以外,航天器结构本身必须具备良好的热力耦合性能。
传统的夹芯板通常采用固定杆径的桁架点阵、固定壁厚的蜂窝结构或固定孔隙率的泡沫多孔结构作为芯子层,芯子层相对密度处处相等,这类结构在设计时很难兼顾承力、隔热性能,在实际力热载荷耦合作用下往往会出现应力、温度分布不均的情况,使得材料利用率大大降低,且可能因局部失效导致整体破坏。
发明内容
本发明的目的在于提出了一种适用于热力工况的梯度点阵夹芯板结构。
实现本发明目的的技术解决方案为:一种适用于热力工况的梯度点阵夹芯板结构,包括两层金属蒙皮以及设置在两层金属蒙皮中间的点阵芯子,所述点阵芯子的表面由三周期极小曲面函数表示。
优选地,所述点阵芯子的曲面函数表达式具体为:
F(x,y,z)=cos(2πx/L)+cos(2πy/L)+cos(2πz/L)-t
式中,x,y,z表示笛卡尔直角坐标系坐标,L表示点阵结构周期大小或胞元尺寸大小,t为梯度参数。
优选地,所述点阵芯子的相对密度为:
式中,t为梯度参数。
优选地,梯度参数的确定方法为:
步骤1:设定点阵芯子的梯度参数t变化范围为[a,b],其中-1<a<b<1;;
步骤2:设初始点阵芯子梯度参数t初始值均为0,构建均匀密度点阵夹芯板模型并导入有限元软件ABAQUS;
步骤3:设定均质点阵夹芯板约束及载荷并进行热力耦合仿真计算,获得温度以及冯·米塞斯等效应力;
步骤4:确定点阵芯子温度最大值tempmax和最小值tempmin及其对应的结点;
步骤5:建立每个结点处梯度参数与温度的线性模型:
t1i=k1tempi+b1
式中t1i表示第i个结点处的梯度参数值,tempi表示第i个结点处的温度值,k1、b1为待定系数;
令温度最高的结点对应梯度参数t1的最大值b,令温度最低的结点对应梯度参数t1的最小值a,代入梯度参数与温度的线性模型,获得待定系数k1、b1的取值;
步骤6:确定点阵芯子冯·米塞斯等效应力最大值stressmax和最小值stressmin及其对应的结点;
步骤7:建立每个结点处梯度参数与应力的线性模型:
t2i=k2stressi+b2
式中,t2i表示第i个结点处的梯度参数值,stressi表示第i个结点处的冯·米塞斯等效应力值,k2、b2为待定系数;
令冯·米塞斯等效应力值最高的结点对应参数t2的最大值b,令冯·米塞斯等效应力值最低的结点对应参数t2的最小值a,代入梯度参数与应力的线性模型,获得待定系数k2、b2的取值;
步骤8:根据梯度参数与温度的线性模型以及梯度参数与应力的线性模型确定每个结点处的梯度参数为:tcouple=f(t1,t2)=t1-t2。
优选地,设定的约束及载荷为:底面固定,顶面受压强,底面对流换热,顶面受温度载荷。
优选地,所述金属蒙皮长宽均为200mm,厚度为5mm。
优选地,所述点阵芯子周期大小L为10mm。
本发明与现有技术相比,其显著优点为:本发明结合有限元仿真,将结构实际工作下的温度、应力分布映射到点阵的相对密度分布上,得到变密度的点阵夹芯结构,增强了结构的整体承力、隔热性能,可广泛用于航空、航天等的热防护结构上。
下面结合附图对本发明做进一步详细的描述。
附图说明
图1是梯度Primitive点阵夹芯板结构示意图。
图2是梯度Primitive点阵芯子俯视图。
图3是梯度Primitive点阵芯子正视图。
图4是点阵夹芯板载荷示意图。
图5是均质Primitive点阵夹芯板热力耦合仿真位移云图。
图6是均质Primitive点阵夹芯板热力耦合仿真温度云图。
图7是梯度Primitive点阵夹芯板热力耦合仿真位移云图。
图8是梯度Primitive点阵夹芯板热力耦合仿真温度云图。
具体实施方式
如图1~4所示,一种适用于热力工况的梯度点阵夹芯板结构,包括两层金属蒙皮1以及设置在两层金属蒙皮1中间的点阵芯子2,所述点阵芯子2为一类由三周期极小曲面衍生而来的点阵结构,曲面类型为Schwartz Primitive,函数表达式为:
F(x,y,z)=cos(2πx/L)+cos(2πy/L)+cos(2πz/L)-t
点阵相对密度越大,承力性能越好,隔热性能越差;点阵相对密度越小,承力性能越差,隔热性能越好。因此,参数t越小,则承力性能越好,隔热性能越差;参数t越大,则承力性能越差,隔热性能越好。对于具体结构来说,应力高的区域需要使参数t减小,应力低的区域需要使参数t增大;温度高的区域需要使参数t增大,温度低的区域需要使参数t减小。
本发明中点阵夹芯板整体长宽各200mm,高40mm,其中金属蒙皮长宽均为200mm,厚度为5mm。点阵芯子周期大小L为10mm,参数t的具体设计过程如下:
步骤1:设定Primitive点阵的梯度参数t变化范围为[a,b],其中-1<a<b<1;
步骤2:设原始点阵夹芯板芯子部分参数t初始值均为0,将生成的均匀密度点阵夹芯板模型导入有限元软件ABAQUS;
步骤3:利用有限元软件ABAQUS为均质点阵夹芯板设定约束及载荷并进行热力耦合仿真计算,输出温度以及冯·米塞斯等效应力。设定的约束及载荷如图3所示,底面固定,顶面受压强,底面对流换热,顶面受温度载荷;
步骤4:读取计算结果文件中点阵芯子部分所有结点的温度数据temp={temp1,temp2,…tempN},其中N为结点总数,并找出最大值tempmax和最小值tempmin;
步骤5:建立每个结点处梯度参数与温度的线性模型。设t1={t11,t12,…t1N},其中N为结点总数,令每个结点处t1的值与该点温度数据temp值线性相关,具体计算方法如下:
设t1与temp函数关系为:
t1i=k1tempi+b1
式中t1i表示第i个结点处的t1值,tempi表示第i个结点处的温度值,k1、b1为待定系数。令温度最高的结点对应参数t1的最大值b,令温度最低的结点对应参数t1的最小值a,即:
解得
由此得t1与temp函数关系为:
由该方法计算出的t1表示仅针对温度分布进行优化的参数;
步骤6:读取计算结果文件中所有结点的冯·米塞斯等效应力数据stress={stress1,stress2,…stressN},其中N为结点总数,并找出最大值stressmax和最小值stressmin;
步骤7:建立每个结点处梯度参数与应力的线性模型。设t2={t21,t22,…t2N},其中N为结点总数,令每个结点处t2的值与该点冯·米塞斯等效应力数据stress值线性相关,具体计算方法如下:
设t2与stress函数关系为:
t2i=k2stressi+b2
式中t2i表示第i个结点处的t2值,stressi表示第i个结点处的冯·米塞斯等效应力值,k2、b2为待定系数。令冯·米塞斯等效应力值最高的结点对应参数t2的最大值b,令冯·米塞斯等效应力值最低的结点对应参数t2的最小值a,即:
解得
由此得t2与stress函数关系为:
由该方法计算出的t2表示仅针对冯·米塞斯等效应力分布进行优化的参数;
步骤8:设tcouple为最终优化结果对应的梯度参数,令tcouple=f(t1,t2)=t1-t2;
根据计算出的梯度参数tcouple生成对应的点阵芯子模型,加上蒙皮之后形成完整的点阵夹芯板;
实施例
一种适用于热力工况的梯度点阵夹芯板结构的具体实施和验证如下:
根据上述步骤求出优化后的Primitive点阵梯度参数t分布,并生成对应的点阵芯子stl文件。
将点阵芯子与上下蒙皮进行布尔并运算,合并为整体,导入3-matic软件划分网格。
导入abaqus软件进行热力耦合仿真计算。
本实施例设计一组实例进行热力耦合仿真验证,包含体积近似相等的均质Primitive点阵夹芯板与变密度Primitive点阵夹芯板,胞元数量均为20×20×3,单胞尺寸10mm,给定参数t范围为[-0.3,0.3]。图5和图6分别为均质Primitive点阵夹芯板位移云图和温度云图,图7和图8分别为梯度Primitive点阵夹芯板位移云图和温度云图。表1为均质点阵夹芯板和梯度点阵夹芯板的顶面最大变形和底面平均温度,在体积近似相等的情况下,点阵夹芯板顶面最大变形降低了15.46%,底面平均温度降低了0.28%。在力热耦合的复杂工况下,变密度点阵夹芯板在保证足够隔热性能的同时,有效降低了面板最大变形,具有实际应用价值。
表1
均质点阵夹芯板 | 变密度点阵夹芯板 | |
顶面最大变形(mm) | 0.841 | 0.711 |
底面平均温度(℃) | 49.466 | 49.328 |
体积(mm<sup>3</sup>) | 964572.239 | 964519.928 |
Claims (4)
1.一种适用于热力工况的梯度点阵夹芯板结构,其特征在于,包括两层金属蒙皮(1)以及设置在两层金属蒙皮(1)中间的点阵芯子(2),所述点阵芯子(2)的表面由三周期极小曲面函数表示;
所述点阵芯子的曲面函数表达式具体为:
F(x,y,z)=cos(2πx/L)+cos(2πy/L)+cos(2πz/L)-t
式中,x,y,z表示笛卡尔直角坐标系坐标,L表示点阵结构周期大小或胞元尺寸大小,t为梯度参数;
所述点阵芯子(2)的相对密度为:
式中,t为梯度参数;
梯度参数的确定方法为:
步骤1:设定点阵芯子(2)的梯度参数t变化范围为[a,b],其中-1<a<b<1;
步骤2:设初始点阵芯子(2)梯度参数t初始值均为0,构建均匀密度点阵夹芯板模型并导入有限元软件ABAQUS;
步骤3:设定均质点阵夹芯板约束及载荷并进行热力耦合仿真计算,获得温度以及冯·米塞斯等效应力;
步骤4:确定点阵芯子温度最大值tempmax和最小值tempmin及其对应的结点;
步骤5:建立每个结点处梯度参数与温度的线性模型:
t1i=k1tempi+b1
式中t1i表示第i个结点处的梯度参数值,tempi表示第i个结点处的温度值,k1、b1为待定系数;
令温度最高的结点对应梯度参数t1的最大值b,令温度最低的结点对应梯度参数t1的最小值a,代入梯度参数与温度的线性模型,获得待定系数k1、b1的取值;
步骤6:确定点阵芯子冯·米塞斯等效应力最大值stressmax和最小值stressmin及其对应的结点;
步骤7:建立每个结点处梯度参数与应力的线性模型:
t2i=k2stressi+b2
式中,t2i表示第i个结点处的梯度参数值,stressi表示第i个结点处的冯·米塞斯等效应力值,k2、b2为待定系数;
令冯·米塞斯等效应力值最高的结点对应参数t2的最大值b,令冯·米塞斯等效应力值最低的结点对应参数t2的最小值a,代入梯度参数与应力的线性模型,获得待定系数k2、b2的取值;
步骤8:根据梯度参数与温度的线性模型以及梯度参数与应力的线性模型确定每个结点处的梯度参数为:tcouple=f(t1,t2)=t1-t2。
2.根据权利要求1所述的适用于热力工况的梯度点阵夹芯板结构,其特征在于,设定的约束及载荷为:底面固定,顶面受压强,底面对流换热,顶面受温度载荷。
3.根据权利要求1所述的适用于热力工况的梯度点阵夹芯板结构,其特征在于,所述金属蒙皮长宽均为200mm,厚度为5mm。
4.根据权利要求1所述的适用于热力工况的梯度点阵夹芯板结构,其特征在于,所述点阵芯子周期大小L为10mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010510048.2A CN111680441B (zh) | 2020-06-08 | 2020-06-08 | 一种适用于热力工况的梯度点阵夹芯板结构 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010510048.2A CN111680441B (zh) | 2020-06-08 | 2020-06-08 | 一种适用于热力工况的梯度点阵夹芯板结构 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111680441A CN111680441A (zh) | 2020-09-18 |
CN111680441B true CN111680441B (zh) | 2023-03-28 |
Family
ID=72435441
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010510048.2A Active CN111680441B (zh) | 2020-06-08 | 2020-06-08 | 一种适用于热力工况的梯度点阵夹芯板结构 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111680441B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112687253A (zh) * | 2020-12-31 | 2021-04-20 | 广州大学 | 一种基于三重周期极小曲面的轻质隔声结构及其制备方法 |
CN114239369B (zh) * | 2022-02-25 | 2022-07-15 | 成都飞机工业(集团)有限责任公司 | 一种蜂窝夹芯结构成型过程的受力分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110929358A (zh) * | 2020-01-21 | 2020-03-27 | 中国人民解放军国防科技大学 | 一种参数化功能梯度立方点阵结构设计方法及材料 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2943174B1 (fr) * | 2009-03-12 | 2011-04-15 | Soitec Silicon On Insulator | Adaptation du parametre de maille d'une couche de materiau contraint |
-
2020
- 2020-06-08 CN CN202010510048.2A patent/CN111680441B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110929358A (zh) * | 2020-01-21 | 2020-03-27 | 中国人民解放军国防科技大学 | 一种参数化功能梯度立方点阵结构设计方法及材料 |
Also Published As
Publication number | Publication date |
---|---|
CN111680441A (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109145427B (zh) | 一种基于三周期极小曲面的多孔结构设计与优化方法 | |
CN108009336B (zh) | 一种微桁架结构承载和热防护结构的多目标优化方法 | |
CN111680441B (zh) | 一种适用于热力工况的梯度点阵夹芯板结构 | |
CN111709096B (zh) | 一种强化自然对流换热的异型翅片结构设计方法 | |
CN107103138B (zh) | 一种激光喷丸变刚度轻量化方法 | |
US20220129595A1 (en) | Design and optimization method of porous structure for 3d heat dissipation based on triply periodic minimal surface (tpms) | |
Niknam et al. | Architected multi-directional functionally graded cellular plates | |
CN107885945A (zh) | 面向3d打印的多尺度孔洞结构轻量化建模方法 | |
CN107563010B (zh) | 基于形状特征的多尺度结构材料一体化设计方法 | |
CN111967089B (zh) | 一种基于渐变胞元尺寸的异型三维点阵结构隔热性能优化设计方法 | |
CN102184298B (zh) | 金属体积塑性成形中有限元分析刚度矩阵存储与生成方法 | |
CN109960878B (zh) | 一种面向高超声速飞行器总体的主/被动热防护***耦合设计方法 | |
CN114756934B (zh) | 一种三维多尺度超材料结构优化设计方法 | |
CN115618503A (zh) | 一种舵翼类结构增材工艺仿真及工艺优化方法 | |
CN112728297B (zh) | 一种基于増材制造提高隔热效果的装置 | |
CN113326582B (zh) | 一种基于应力分布的变密度点阵结构及其设计方法 | |
Sun et al. | Design of a bionic aviation material based on the microstructure of beetle’s elytra | |
CN117637076A (zh) | 一种tpms梯度多孔复合材料制备方法 | |
CN110489912B (zh) | 一种太阳能跨季节土壤蓄热分层切片数值模拟的方法 | |
CN105426640A (zh) | 基于罚函数的多组件结构***布局优化设计方法 | |
CN112836399A (zh) | 一种基于有限元算法的非线性接触热阻热分析求解方法 | |
CN109299499B (zh) | 考虑修正因子的多步骤结构优化设计方法及飞行器 | |
CN106909747A (zh) | 一种热对流扩散***中模糊参数隶属度函数识别方法 | |
CN112115563B (zh) | 一种热压罐成型框架式模具的整体拓扑优化设计方法 | |
CN115935730B (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 |