[发明专利]非平面波在多胞材料中传播的波阵面识别方法

说明书

本发明公开了一种非平面波在多胞材料中传播的波阵面识别方法，基于数值仿真软件生成的网格节点位移场，采用基于最小二乘法离散局部变形梯度的局部应变场计算方法，采用连续介质力学计算方法对Green应变张量转化为工程应变张量，通过八面体等效应变计算方法，确定非平面波阵面位置，从而解决多胞材料内部非平面波的识别问题，为泡沫宏观结构呈现复杂结构或者载荷非单一方向等一般应用环境下的非平面波波阵面识别提供了较准确的计算方法，同时为非平面波阵面在胞孔材料中的传播规律认识提供了有力工具和奠定了研究基础。该发明可应用于航空航天、军事防护、汽车制造、产品包装等抗爆抗冲击防护设计的基础科学研究和工程设计。

技术领域

本发明涉及多胞材料冲击波波阵面识别技术领域，尤其涉及一种在爆炸或冲击条件下非平面波在多胞材料中传播的波阵面识别方法。

背景技术

多胞材料作为冲击防护类缓冲材料，广泛应用于航空航天、军事防护、汽车制造、产品包装等，冲击波在多胞材料中的传播规律对于抗爆抗冲击防护设计不仅具有基础科学研究意义，而且对于军事和民用防护也有重要指导意义，而冲击波波阵面在多胞材料这类非连续介质中的识别技术是该领域研究的重要课题。

目前在多胞材料冲击波波阵面识别技术领域，分为两条路线。第一条路线是实验手段，如附图1所示，通过高速摄影拍摄多胞材料试件在动态压缩下的高速图像，然后通过DIC数字相关技术进行分析，得到试件表面的应变分布场，根据应变突跃位置，确定波阵面位置；第二条路线是数值仿真+理论计算，基于实验手段无法获取试件内部应变分布场，采用多胞材料细观数值模拟方法，获取冲击条件下的单元节点处的位移场数据，再采用基于最小二乘法离散局部变形梯度的局部应变场计算方法，得到节点处的Green应变张量，并通过连续介质力学转化为冲击方向的应变，从而识别应变突跃位置，确定波阵面位置。

以上多胞材料冲击波波阵面识别技术领域的两条路线中，第一条的实验路线受限于高速摄影只能拍摄试件外表面动态变形，对内部变形及冲击波波阵面的形态无法获取；第二条的数值仿真+理论计算路线只能识别多胞材料中的平面冲击波波阵面，由于只考虑了正应变分量而未涉及剪应变分量，因此非平面波无法识别或界面识别不准确。基于此，本发明的提出是为了解决多胞材料内部非平面波的识别问题。

发明内容

本发明的目的就在于为了解决上述问题而提供一种非平面波在多胞材料中传播的波阵面识别方法。

为了实现上述目的，本公开提供一种非平面波在多胞材料中传播的波阵面识别方法，包括：

通过细观数值仿真方法，建立多胞材料的数值仿真模型，并开展爆炸或冲击条件下的动态数值模拟；

通过局部应变场计算方法，得到多胞材料的局部应变场；

采用连续介质力学计算方法，将Green应变张量转化为工程应变张量；

通过八面体等效应变计算方法，得到材料各节点所组成的八面体等效应变场，并通过软件确定非平面波阵面位置。

可选地，多胞材料构建有限元模型，在指定的材料宏观区域Π₀内按照均匀概率分布散布N个成核点，并且控制成核点间的距离δ不小于给定的距离δ_min；

通过MATLAB软件生成胞棱、胞面几何信息，并对胞面进行网格划分，在有限元软件ABAQUS中得到多胞材料的数值仿真模型。

可选地，对于规则型细观结构的多胞材料，利用有限元方法周期性构建有限元模型。

可选地，对于随机型细观结构的多胞材料，利用3D-Voronoi随机细观结构构建有限元模型。

可选地，在多胞材料中，采用基于最小二乘法的离散局部变形梯度来描述多胞性带来的离散位移场，从而得到多胞材料的局部应变场。