Failure mechanism and dynamic response of a composite lattice structure under intense explosion loadings
-
摘要: 基于碳纤维增强复合材料面板与金属芯层，设计出金字塔型复合点阵夹芯结构。利用地面爆炸冲击实验，研究复合点阵结构在强爆炸载荷作用下的损伤机理和失效模式。基于材料的细观损伤机理，构建复合材料面板的三维渐进损伤模型和金属芯层的Johnson-Cook损伤模型，并结合有限元方法发展了复合点阵结构的爆炸冲击响应预报模型。开展了不同载荷工况下结构的动态响应特性分析，结合实验测试结果分析了结构抗爆性能的主要影响因素。研究表明，在近距离强爆炸载荷作用下，复合点阵结构整体构型基本保持完好，仅局部出现失效现象，主要失效形式为边缘区域面芯脱粘和局部芯层杆件断裂，但结构整体上仍具有较好的承载能力。探讨了考虑多种载荷条件和结构参数相关变量的毁伤函数，给出了结构的可行设计域。研究结果可为装备关键部件轻量化/抗爆设计提供参考。Abstract: In order to comply with the requirements of explosive shock wave protection for the new generation of equipment structures, it is necessary to design a lightweight, high energy absorption ratio structure and further systematically understand its dynamic responses under explosion loadings. A composite lattice sandwich structure with pyramidal truss core was designed, which consisting of carbon fiber reinforced composite panels and metal cores. The explosion experiments were carried out. The failure mechanism and damage mode of this composite lattice structure under intense explosion shock loadings were analyzed. Based on the failure mechanism in mesoscale of the material, both the three-dimensional progressive damage model of the composite panels and the Johnson-Cook damage model of the metal cores were constructed. By combining with the finite element method, a numerical model for predicting explosion shock response of the composite lattice structure was developed. Both the bonding properties between layers of the composite panels, and the performances of the adhesive layer between the panels and the cores were considered in the numerical model. The initial damage criterion based on strain description was established, and the damage dynamic evolution equations corresponding to different damage modes were given. A damage variable was introduced to characterize the attenuation degree of stiffness properties of the damaged elements. Furthermore, the stress of damaged elements could be obtained. The dynamic responses of this structure under different loadings were analyzed using the developed model. The main mechanisms influencing the explosion protection properties of the composite lattice structure were discussed based on both simulated and experimental results. It is revealed that the local failures occur when the composite lattice structure is exposed and close to explosion loadings. The main failure modes are the debonding between the composite panels and the pyramidal truss cores in the edge area, and the fracture of local struts. However, the overall configuration of this composite lattice sandwich structure is basically intact and it still has a good carrying capacities. The damage function considering various variables of load conditions and structural parameters was discussed. The feasible domain for this structure design was given. These research results can provide theoretical basis and technical support for the designing and safety evaluation of lightweight, explosion protection structure of key equipment components.
-
新一代舰船、战机、军车等装备系统在服役过程中通常会遭受强爆炸冲击载荷[1-3]。轻质点阵夹芯结构具有高结构效率、高强度、多功能、设计性强的优点，且具有良好的抗冲击吸能特性，是新一代装备系统的理想承载结构方案[4-6]。强爆炸冲击载荷下点阵夹芯结构会发生复杂的渐进失效过程[7-8]，阐明在爆炸冲击载荷作用下点阵夹芯结构的损伤演化规律和动态响应特性，可为装备承载/防护结构设计、安全性评价提供必要的理论基础和科学依据，具有十分重要的意义。
根据芯子构型的不同，点阵夹芯结构通常包括金字塔型、四面体型和 Kagome 型等[9]。近年来，国内外学者针对爆炸冲击载荷下点阵夹芯结构的吸能特性已开展了系统深入的研究工作，研究表明，与传统结构相比，点阵夹芯结构的抗爆吸能特性优异、结构效率高，且充分利用芯层的空隙可开展储能、吸声、降噪、热控等多种功能设计[4-6]。Evans等[10]、Li等[11]系统地研究了随机芯层夹芯结构（如金属泡沫夹芯结构等）与周期芯层夹芯结构（如点阵夹芯结构）的吸能、减振特性，认为通过合理设计和优化，周期芯层夹芯结构可实现力学性能和其他性能的最大化，具有良好的应用潜力。Qi等[12]、王涛等[13]利用爆炸冲击实验详细研究了泡沫铝夹芯结构的变形破坏过程，总结了泡沫铝夹芯结构各部件的失效模式，构建了爆炸冲量与夹芯结构最大变形量之间的关系。Yungwirth等[14]、Wadley等[15]研究了多层金字塔点阵夹芯结构在水中爆炸冲击载荷作用下的动态响应规律，与同等密度的实体板相比，多层金字塔点阵夹芯结构的透射冲击波强度降低了28%。Dharmasena等[16]研究了球形沙子爆炸冲击载荷下金字塔点阵夹芯结构的动态结构响应，相同载荷下金字塔点阵结构背面板的最大变形量明显小于相同质量的实体板，进一步证实了点阵夹芯结构的良好抗爆吸能能力。亓昌等[17]研究了芯子密度、排布方式、结构参数等对金字塔点阵夹芯结构抗爆、抗侵彻性能的影响规律，确定了结构吸能特性的关键影响参数。
点阵夹芯结构在爆炸冲击载荷下的毁伤特性与载荷工况、结构参数、材料参数等多种物理参数密切相关，涉及微细观材料损伤、宏观结构失效等多尺度特征，构建点阵夹芯结构的毁伤数据库，准确地描述结构的易损特性，可以为国防工业中的防护设计与毁伤评估提供有益参考。国内外学者针对承载结构毁伤效应评估方法及毁伤函数构建方法开展了深入的研究。彭航[18]研究了典型建筑物爆炸载荷作用下的整体毁伤效应，引入跨中损伤范围比例系数、残余挠度两个无量纲参数，用于表征结构的破坏程度，拟合获得了载荷条件、材料参数相关的建筑物结构毁伤程度工程计算公式。冯晓伟等[19]、Han等[20]详细地研究了爆炸冲击作用下飞机机翼结构的毁伤效应评估方法，确定了机翼结构发生失效时的超压、作用时间，建立了基于冲击波超压、冲量的机翼结构临界失效准则。胡榕等[21]研究了机场跑道在内爆载荷作用下的毁伤效应，利用爆炸实验分析了弹坑形态及裂纹破坏特征，采用弹坑半径、环向裂纹半径、最大爆腔半径、弹坑深度等关键破坏特征参量，定量地描述了跑道内爆毁伤场形态，建立了不同等级跑道下装药量、装药埋深相关的毁伤函数。
复合材料化是实现装备结构轻量化的重要手段，也是新一代装备结构的发展趋势。为了进一步提高点阵夹芯结构的比刚度和比强度，通常采用复合材料面板作为点阵夹芯结构的上、下面板，芯层仍采用金属材料，从而组成复合材料/金属混杂点阵夹芯结构。复合点阵夹芯结构在满足高承载、抗冲击性能的前提下，具有更高的结构效率，因此，可满足新一代装备的承载、防护需求。在爆炸冲击载荷作用下，复合点阵夹芯结构的失效模式呈现出多样性特点，包括欧拉屈曲、面板塌陷、起皱、分层、面芯脱粘以及芯层杆件屈曲、断裂、压溃等多种可能的失效模式[22-24]。复合点阵夹芯结构在空间上有材料与结构双重多尺度特征，爆炸冲击载荷作用又增加了时间多尺度特征。点阵夹芯结构爆炸冲击载荷下的损伤机理尚不十分清楚，失效模型、毁伤函数等方面的理论与数值方法研究尚不完善。
本文中，以碳纤维增强复合材料面板和金属芯层组成的复合金字塔型点阵夹芯结构为研究对象，开展爆炸冲击实验，揭示复合点阵夹芯结构在强爆炸冲击载荷下的损伤演化规律；基于复合材料三维渐进损伤准则和金属材料Johnson-Cook失效判据，建立复合点阵夹芯结构的爆炸冲击响应模型，分析典型爆炸载荷下复合点阵结构的吸能特性、损伤机理及失效模式，给出复合点阵夹芯结构的毁伤函数；研究结果以期为装备关键部件轻量化/抗爆设计提供参考。
1. 复合点阵夹芯结构设计与制备
复合材料具有高比刚度和比强度的特点，能够在满足承载、吸能要求的前提下减轻结构质量；金属点阵芯子具有高强度、高承载特性，能够抵抗面板传递进来的冲击波载荷，且可通过塑性变形等方式吸收掉部分冲击能量；金字塔型点阵芯子兼具承载和其他多功能设计，制造工艺相对简单，是目前应用较为广泛的一种芯子构型。本文中，选用碳纤维增强环氧复合材料作为面板母材料，选用铝合金作为金字塔型芯子母材料，开展金字塔型复合点阵夹芯结构设计。碳/环氧复合材料预浸料是由上海力硕复合材料科技有限公司提供，单层厚度为0.1 mm。
复合点阵夹芯结构如图1(a)所示，主要参数为：碳/环氧复合材料面板尺寸为300 mm×300 mm，上、下面板厚度均为2 mm；面板由20层单向碳/环氧复合材料预浸料层合而成，铺层角度为[0, 90, 45, −45]5，如图1(b)所示；金字塔型芯子代表性单胞的设计参数如图1(c)所示，具体数值在表1中给出。复合点阵夹芯结构的制备工艺为：首先，利用模压成型工艺制备出复合材料面板，采用激光切割、线切割等方法制备出芯子杆件；然后，利用嵌锁组装工艺，将32条长度为186 mm的芯子杆件制成金字塔芯层；最后，采用高强度环氧胶膜（将复合材料面板和金字塔芯层粘接，经固化后，最终制备形成复合点阵夹芯结构，如图2(a)所示。
图 1 复合点阵夹芯结构示意图及设计参数Figure 1. Schematic diagram and design parameters of a composite lattice sandwich structure表 1 金字塔型芯子代表性单胞主要设计参数Table 1. Main design parameters of representative structure cell of pyramidal truss coreh/mm p/mm t1/mm s/mm t2/mm θ/ (°) ha/mm hb/mm 15 35.75 1.59 8 1.5 45 3 1.5 复合点阵夹芯结构的面密度为7.38 kg/m2，相对密度为2.40%，与金属点阵夹芯结构相比，减重约37.9%。碳/环氧复合材料和铝合金的基本材料性能参数在表2中给出。
表 2 碳/环氧复合材料的基本材料性能参数Table 2. Mechanical properties of carbon/epoxy composite materials密度/(kg·m−3) 拉压模量/GPa 泊松比 剪切模量/GPa 拉伸强度/MPa 压缩强度/MPa 面内剪切强度/MPa 层间剪切强度/MPa 1600 123 (0°) 0.41 (0°) 4.8 1400 (0°) 850 (0°) 60 16 8.3 (90°) 0.26 (90°) 18 (90°) 96 (90°) 注：0°是指力作用在纤维方向；90°是指力作用在垂直纤维方向。 2. 爆炸冲击实验
复合点阵夹芯结构的爆炸冲击实验在西安近代化学研究所爆炸实验塔内开展（图3）。实验过程中，将高能炸药悬挂于空中，距离地面高度约为1.5 m；将复合点阵夹芯结构样件固定在支架上（图2(b)），并调整支架高度使样件中心点与炸药高度相同，以确保炸药引爆后产生的爆炸冲击波垂直作用于结构。设置支架与炸药之间的水平距离，可以研究不同爆炸距离下复合点阵结构的动态冲击响应特性。为了测量到达样件表面的爆炸冲击波压力，在样件同等水平距离、同等高度位置处放置压力传感器，记录冲击波脉冲变化规律。
爆炸力学中，通常采用比例距离衡量到达物体表面的爆炸冲击载荷强弱，比例距离
ˉR 可以定义为：ˉR=L3√q (1) 式中：q为炸药的TNT当量，kg；L为爆炸距离，m。根据式（1）可知，在相同大气条件下，当爆炸距离与爆炸威力立方根的比值相等，即比例距离相同时，则到达物体表面的冲击波阵面超压相同。分别设置了900 g TNT当量−1.2 m爆距以及1 000 g TNT当量−1.86 m爆距2种载荷工况（表3），根据式（1），2种工况的比例距离分别为1.24和1.86 m/kg1/3。
表 3 爆炸实验载荷工况Table 3. Loading conditions of explosion experiments工况 TNT当量/g 爆炸距离/m 比例距离/(m·kg−1/3) 1 900 1.20 1.24 2 1000 1.86 1.86 3 1500 1.20 1.05 4 3000 1.20 0.83 图4给出了实验过程中压力传感器测量得到的爆炸源引爆后的冲击波超压曲线。可以看出，冲击波脉冲的整体作用时间较短，迅速达到峰值后，在1～2 ms内几乎衰减为零。为了估算不同载荷条件下的冲击波超压数值，国内外科研人员开展了大量的行之有效的研究工作[25]，其中CONWEP爆炸分析经验模型具有较高的精度，是常用的爆炸载荷估算模型。本文中采用CONWEP爆炸分析模型获得的两种工况下的冲击波超压数值也在图4中给出。可以看出，工况1下，实验测得的冲击波超压峰值为2.11 MPa，估算值为2.6 MPa，二者之间的偏差为23.2%；工况2下，冲击波超压的实验测量值为0.68 MPa，估算值为0.74 MPa，二者之间的偏差仅为8.8%。误差原因可能是：一方面，实验测量值受到场地布置尺寸偏差、爆轰产物干扰等因素影响；另一方面，经验模型忽略了空气与结构之间的流固耦合作用，会导致冲击波作用在结构表面上的反射超压偏大。需要指出的是，爆炸过程中，爆炸载荷与结构的作用过程较为复杂，除了产生强冲击波超压外，冲击波入射、透射、反射还会对结构产生一定的振动冲击效应。本文中重点研究冲击波超压作用下复合点阵结构的动态响应，忽略振动冲击影响。
3. 基于渐进损伤机理的复合点阵结构动态响应预报模型
3.1 复合材料面板的渐进损伤模型
复合点阵夹芯结构的上、下面板采用碳/环氧复合材料，由于复合材料各向异性的特点，复合点阵夹芯结构的失效过程较为复杂。在爆炸冲击载荷作用下，结构上、下面板产生的面外应力较为明显，平面应力状态假设不再适用，必须考虑复合材料面板厚度方向上的应力。
首先构建复合材料面板单个铺层的细观几何模型，然后定义铺层的材料主方向，赋予相应的材料性能参数。对于单个铺层，假设纤维方向为方向1，垂直于纤维方向以及面外厚度方向分别为方向2、3，则可以认为单个铺层材料为横观各向同性材料，其中平面2-3为横观各向同性面。在弹性变形条件下，单个铺层材料的应力-应变关系为：
(εe11εe22εe33γe23γe13γe12)=(1E11−ν12E11−ν13E11−1E12−ν23E111E331G231G131G12)(σ11σ22σ33σ23σ13σ12) (2) 式中：上标e表示各向异性材料的应力-应变关系在线弹性范围内。
为了描述复合材料结构的渐进损伤行为，Hashin将复合材料结构的损伤归因于基体损伤、纤维损伤以及两种组分相的耦合损伤，同时考虑了拉伸、压缩等力学载荷作用下复合材料损伤机理的多样性，建立了横观各向同性复合材料结构的损伤初始准则。借鉴Hashin损伤初始准则：
{Afσ11+Bfσ211+1S2L(σ212+σ213)=1纤维模式Am(σ22+σ33)+Bm(σ22+σ33)2+1S2T(σ223−σ22σ33)+1S2L(σ212+σ213)=1基体模式 (3) 式中：下标f、m分别代表纤维和基体，
A 、B 为未知系数，SL 、ST 分别为材料轴向和横向剪切强度。面板铺层材料出现损伤后，其刚度发生衰减，单元应力也逐渐降低，直至单元发生失效。由于应变在材料渐进损伤过程中保持连续状态，可以采用应变作为变量描述复合材料面板的失效行为，进而建立基于应变描述的损伤初始准则：
\left\{ \begin{array}{lll}{f}_{\rm{ft}}=\left(\dfrac{{\varepsilon }_{11}}{{\varepsilon}'_{11,{\rm{t}}}}\right)^{2}+\alpha \left(\dfrac{\gamma_{12}^2+\gamma_{13}^2}{\gamma'_{12}}^2\right)& \alpha =0,1 ，\varepsilon _{11} {\text{＞}} 0 & {\text{纤维拉伸损伤}}\\ f_{\rm{fc}}=\left(\dfrac{\varepsilon _{11}}{\varepsilon '_{11,{\rm{c}}}}\right)^{2} & \varepsilon _{11}{\text{＜}} 0 & {\text{纤维压缩损伤}}\\ f_{\rm{mt}}=\left(\dfrac{\varepsilon_{22}+\varepsilon_{33}}{\varepsilon '_{22,{\rm{t}}}}\right)^2+\dfrac{\gamma_{23}^2-\varepsilon_{22}\varepsilon _{33}}{\gamma '_{23}^2}+\dfrac{\gamma_{12}^2+\gamma _{13}^2}{{\gamma'_{12}}^2}& \varepsilon _{22}+\varepsilon _{33} {\text{＞}} 0 & {\text{基体拉伸损伤}}\\ f_{\rm{mc}}=\left[{\left(\dfrac{{\varepsilon }'_{22,{\rm{c}}}}{2{\bar{\gamma }}_{23}}\right)}^{2}-1\right]\left(\dfrac{{\varepsilon }_{22}+{\varepsilon }_{33}}{{\varepsilon }'_{22,{\rm{c}}}}\right)+{\left(\dfrac{{\varepsilon }_{22}+{\varepsilon }_{33}}{2{\gamma }'_{23}}\right)}^{2}+\dfrac{{\gamma }_{23}^{2}-{\varepsilon }_{22}{\varepsilon }_{33}}{{\left({\gamma }'_{23}\right)}^{2}}+\dfrac{{\gamma }_{12}^{2}+{\gamma }_{13}^{2}}{{\left({\gamma }'_{12}\right)}^{2}}& {\varepsilon }_{22}+{\varepsilon }_{33}{\text{＜}} 0 & {\text{基体压缩损伤}}\\ {f}_{ij}={\left(\dfrac{{\gamma }_{ij}}{{\gamma }'_{ij}}\right)}^{2}& i,j=1,2,3； i\ne j & {\text{剪切损伤}}\end{array} \right. (4) 式中：
ε′ 和γ′ 分别表示材料拉压失效应变和剪切失效应变，下标t、c分别表示拉伸和压缩性能。为了模拟复合材料面板的渐进损伤过程，首先，构建损伤初始准则，根据损伤初始准则来判断材料是否发生损伤，即当达到损伤初始准则时，表明材料开始出现损伤，此时材料发生刚度衰减。其次，建立损伤动态演化方程，复合材料结构内部的损伤状态与损伤的扩展速率有关，针对5种损伤模式，通过引入损伤状态变量
ω(t+Δt) 来描述当前时刻材料的损伤状态，则损伤扩展速率即为损伤状态变量的导数˙ω(t+Δt) 。结合损伤初始准则和损伤动态演化方程，构建复合材料面板的三维渐进损伤模型，进而描述结构的渐进损伤过程。不同损伤模式对应的损伤动态演化方程为：{˙ωft(t+Δt)=Ω0,f+Ω1,fω1t(t)[(ε11(t+Δt)ε′11,t(1−ω1t(t)))2−1]纤维拉伸˙ωfc(t+Δt)=Ω0,f+Ω1,fω1c(t)[(ε11(t+Δt)ε′11,c(1−ω1c(t)))2−1]纤维压缩˙ωmt(t+Δt)=Ω0,m+Ω1,mω2t(t)[(ε22(t+Δt)ε′22,t(1−ω2t(t)))2−1]基体拉伸˙ωmc(t+Δt)=Ω0,m+Ω1,mω2c(t)[(ε22(t+Δt)ε′22,c(1−ω2c(t)))2−1]基体压缩˙ωij(t+Δt)=Ω0,s+Ω1,sωij(t)[(εj(t+Δt)ε′ij(1−ωij(t)))2−1]i,j=1,2,3；i≠j剪切损伤 (5) 式中：
Ω0 、Ω1 分别为损伤形成速率和损伤扩展速率，下标s表示剪切性能。爆炸冲击波载荷的整体作用时间较短，其时间步长
Δt 通常处于亚微秒量级，则可采用线性方程对损伤演化规律进行近似，即可得到相邻两个时刻损伤状态变量之间的关联关系：ωi(t+Δt)=ωi(t)+˙ωi(t+Δt)Δtωi(t)∈[0,1] (6) 引入损伤变量
di ，可以计算复合材料结构发生损伤后，其损伤单元刚度性能的衰减程度，从而实现损伤过程的模拟。损伤变量di 的表达式为：{d1=1−(1−ωft)(1−ωfc)d2=1−(1−ωmt)(1−ωmc)d12=ω12d23=ω23d13=ω13 (7) 则损伤单元材料各方向上的刚度参数为：
Eij,d=Eij,0(1−dij)，Gij,d=Gij,0(1−dij) (8) 式中：
Eij,d 、Gij,d 代表损伤发生后单元材料的模量，Eij,0 、Gij,0 为原始材料的模量。损伤单元的应力为：
σd(t+Δt)=Cd(t+Δt)ε(t+Δt) (9) 式中：
Cd 为损伤刚度矩阵。3.2 金属芯层的Johnson-Cook损伤模型
爆炸冲击载荷作用下，点阵夹芯结构的金属芯子通过塑性变形、裂纹扩展等行为吸收冲击能量。研究表明，Johnson-Cook模型(简称JC模型)可以很好地描述金属材料的应变率硬化与塑性硬化过程，其表达式为：
ˉσ=(A+Bˉεnp)[1+Cln(˙ˉεp˙ε0)] (10) 式中：
ˉσ 为等效应力，A为屈服强度，B、n为塑性硬化参数，ˉεp 为等效塑性应变，C为材料应变率强化参数，˙ˉεp 为等效塑性应变率，˙ε0 为参考塑性应变率。在高应变率假设下，金属芯子单元材料的损伤状态变量可以由等效塑性应变描述：
ω=∑ˉεp/ˉεp,F (11) 式中：
ˉεp,F 为断裂应变。3.3 复合点阵夹芯结构的动态冲击响应模型
将上述复合材料面板的三维渐进损伤模型和金属芯子的Johnson-Cook损伤模型编写为VUMAT子程序，利用有限元方法，建立复合点阵夹芯结构的爆炸冲击响应分析模型，如图5所示。为了考虑层间开裂与面芯脱粘失效，采用双线性粘接关系，在面板各铺层之间、面板与芯子之间分别赋予粘接接触属性。在显示动力学分析中，考虑爆炸冲击载荷的脉冲宽度，计算总时间设置为2 ms。为保证计算准确度，面板与芯子均采用C3D8R实体单元。将结构的上面板上表面设置为迎爆面，上、下面板的4条边均设置为固支边界条件。
图 5 复合点阵结构的爆炸冲击响应有限元分析模型Figure 5. The finite element model for predicting explosion shock response of composite lattice structure利用爆炸冲击响应有限元模型，对复合点阵结构开展显示动力学分析，可以获得爆炸冲击波载荷作用下复合点阵夹芯结构的变形、损伤、失效规律，进而可分析结构的吸能机理。
4. 强爆炸载荷下复合点阵结构动态响应规律分析
图6给出了工况1条件下复合点阵夹芯结构爆炸冲击实验后的样件照片及局部失效情况。从图中可以看出，结构上面板内侧与芯子粘接处出现了脱粘现象，上面板与芯层完全分离（图6(a)~(b)）。部分芯层杆件发生屈曲，部分杆件在连接结点处发生断裂（图6(c)），这些破坏现象主要出现在结构与支架的固定边缘附近。此外，迎爆面边缘固支处出现一些局部裂缝（图6(d)）。同时可以发现，在爆炸冲击载荷作用下，结构样件下面板保持完好，未出现明显的损伤，也未出现塑性变形，这反映了轻质复合点阵结构能够抵抗一定强度的爆炸冲击载荷，可以保护结构背面的器件不受冲击破坏。
图 6 复合点阵结构爆炸冲击实验后的样件照片及局部失效情况Figure 6. Photographs and local failures of composite lattice structure specimens after explosion shock experiments由工况1条件下复合点阵结构上面板承受的冲击波超压脉冲曲线（图4）可以看出，炸药引爆后，经过约0.75 ms爆炸冲击波到达结构上表面，0.78 ms时上面板承受的超压达到峰值2.6 MPa。随后冲击波超压迅速减小，脉冲有效作用宽度约为0.8 ms。利用本文中建立的爆炸冲击响应有限元模型，开展复合点阵结构工况1载荷条件下的显示动力学分析，可以模拟结构的动态响应过程。图7给出了复合点阵结构工况1载荷下的位移场。从图7可以看出，由于炸药空中自由爆炸冲击波为球面波，爆心恰好到达结构上面板的中心，因此，上面板中心位置处的变形幅度最大，1 ms时刻上面板的最大位移为6.42 mm。
图 7 复合点阵夹芯结构爆炸冲击载荷下的位移变化云图Figure 7. Strain response of composite lattice sandwich structure under explosion shock loadings图8给出了结构上、下面板中心点处的位移-时间变化曲线，可以看出，当爆炸冲击波到达上面板后，上面板瞬时产生较大位移，在爆炸开始后的1.125 ms位移达到最大值−8.38 mm。随后上面板开始回弹，在1.5 ms时越过初始位置，继续回弹至2.98 mm处。回弹位移相比最大位移下降了64.47%，表明结构在变形过程中发生了能量耗散。与上面板相比，下面板结构响应滞后约0.05 ms，最大位移为−7.34 mm。上、下面板最大位移之差即为结构金字塔芯层的最大压缩量，为1.04 mm。从图中还可以看出，上面板的回弹量（2.98 mm）大于下面板回弹量（1.76 mm），且回弹方向一致，即上、下面板之间的距离大于初始距离，这表明结构芯层与面板之间发生了脱粘失效。在爆炸实验后的样件中同样出现了典型的面芯脱粘失效模式，这充分验证了本文中爆炸冲击响应有限元模型的正确性。
图 8 复合点阵夹芯结构爆炸冲击载荷下的位移变化曲线Figure 8. Strain response curves of composite lattice sandwich structure under explosion shock loadings工况1爆炸载荷作用下复合点阵夹芯结构的应力变化规律，如图9所示。从图中可以看出，复合点阵夹芯结构上、下面板的边缘区域出现较大应力，最大应力约为600 MPa，这可能是由于边缘位置施加了固支边界条件；而上、下面板其余区域的应力值均处于较低水平。金字塔芯层由于使用较小密度的结构，承受了来自面板较大的爆炸载荷，其整体应力水平略高于面板大部分区域，最大应力为483.3 MPa，未超过铝合金材料的强度极限（图9(b)）。芯层总体上保持较好的完整性，整体受力情况比较均匀，这反映了结构设计的合理性。需要指出的是，在航空、航天、船舶领域的实际工程问题中，承载结构除了承受高频冲击载荷外，还需抵抗低频振动载荷。而在防护结构轻量化设计时，高频冲击防护与低频振动防护往往是相互制约、相互矛盾的。比如：高频冲击防护结构设计时，可以采用吸能特性优异、比密度小的多孔夹层结构，然而这不利于减小冲击振动的振幅；另一方面，低频振动防护结构设计时，可以通过提高材料的密度和刚度实现，但这对轻量化、吸能是不利的。因此，在两种载荷都很突出的工程问题中，需要综合考虑、权衡设计，确定出最优结构形式。
图 9 复合点阵夹芯结构爆炸冲击载荷下的应力变化规律Figure 9. Stress response of composite lattice sandwich structure under explosion shock loadings图10 给出了爆炸载荷作用下复合点阵夹芯结构上、下面板单个铺层的应力云图。在上、下面板厚度方向上每隔0.4 mm选取一个铺层，分析其应力分布规律。从图中可以看出，由于面板各铺层纤维铺设方向不同，载荷到达每个铺层的时间也不相同，因此，各铺层的应力分布规律都不相同。同时，靠近迎爆面铺层的应力水平整体高于远离迎爆面的铺层。各铺层面内的应力情况并不均匀，面板边界区域的应力较大，这是由于数值计算时在面板四周施加了固支边界条件。因边界区域为非有效考查区域，评价结构毁伤效果时，可以忽略边界处的应力情况。比较图10(a)和(b)，在1.25 ms时刻，下面板各铺层的应力明显高于上面板各铺层应力，这可能是因为上面板通过变形、面芯脱粘等方式吸收了部分冲击能量。由于上、下面板采用相同的铺层角度设计方案，因此对应铺层处的应力分布较为相似。
图 10 复合点阵夹芯结构上、下面板单个铺层的应力变化规律Figure 10. Cloud maps of stress distribution between layers of composite lattice structure panels爆炸冲击载荷作用下复合点阵夹芯结构内部的能量转化关系由图11给出。从图中可以看出，随着冲击过程的进行，爆炸冲击波作用在结构上的能量逐渐转化为结构内能和结构动能。冲击波到达结构面板后，冲击能量首先转化为结构的动能，随后在结构塑性变形、黏性耗散、摩擦耗散等效应作用下，冲击能量迅速转化为内能。1.11 ms时刻，爆炸冲击波作用在结构上的总能量达到最大值401.8 J。此时结构吸收的内能达到393.2 J，能量吸收率为97.86%。此后冲击能量不再增加，结构动能在回弹震荡中逐渐耗散为内能。
图 11 复合点阵夹芯结构爆炸冲击过程中的能量转化关系Figure 11. Energy absorption curves of composite lattice sandwich structure exposed to explosion loadings在航空航天领域，为了满足结构轻量化与高抗爆设计要求，通常采用面比吸能（areal specific energy absorption, ASEA）作为结构吸能效率的评价指标。本文中，同样采用面比吸能评价复合点阵夹芯结构的吸能效率，面比吸能
EA 是指单位面密度结构所能吸收的能量，其表达式为：EA=Enma (12) 式中：ma为结构单位面积的质量，En为结构吸收的爆炸冲击总能量。结合复合点阵结构的吸能曲线（图10）和主要结构参数，利用式（12），可以得到本文中复合点阵夹芯结构的单位面积质量为7.38 kg/m2，工况1条件下的面比吸能为54.4 J·m2/kg。
5. 复合点阵结构抗冲击性能影响参数分析与毁伤函数探讨
图12给出了复合点阵夹芯结构爆炸冲击载荷（工况2）下的失效模式。从爆炸冲击实验后的样件照片（图12(a)）中，发现了明显的面芯脱粘失效模式。同时，在有限元模拟结果（图12(b)）中，也发现了类似的面芯脱粘现象。这是由于复合点阵夹芯结构的面板与芯层之间采用胶膜粘接，由于芯层与面板之间的粘接接触面积较小，粘接强度较低，在强动态载荷冲击下，粘接界面容易发生破坏，进而导致面芯之间出现脱粘失效。模型预报结果与实验结果基本吻合，这也进一步验证了数值模型的正确性。此外，爆炸实验与数值模拟结果均证明，面芯脱粘是复合点阵夹芯结构在爆炸冲击载荷下的主要失效模式，如何提高面板与芯层之间的界面粘接强度是提高轻质点阵夹芯结构抗爆能力的关键。
图 12 复合点阵夹芯结构在爆炸冲击载荷下的失效模式Figure 12. Failure modes of composite lattice sandwich structure under explosion loadings为了研究爆距、药量等载荷条件以及面板厚度、芯子构型等结构参数对复合点阵夹芯结构抗冲击性能的影响规律，利用爆炸冲击响应有限元模型，开展了不同载荷条件和结构参数下复合点阵结构的动态响应特性分析（图13），获得了一些有益的结果。图13(a)给出了1 500 g TNT当量爆炸冲击载荷（工况3）下，金字塔点阵芯层的位移分布规律。图中黑色部分表示金字塔点阵芯层的初始位置，彩色云图代表爆炸冲击后的结构变形。通过比较两者之间的相对位移可以看出，金字塔点阵芯层除沿爆炸冲击方向发生较大的弯曲变形外，各结点在水平方向上也出现较大程度的滑移，这表明结构已出现大范围面芯脱粘失效，这与其他工况条件下获得的失效模式比较类似。
图 13 复合点阵夹芯结构不同载荷条件下的结构响应Figure 13. Structural responses of composite lattice sandwich structures under different loadings图13(b)给出了3 000 g TNT当量爆炸冲击载荷（工况4）作用下不同时刻金字塔芯层的位移变化规律。可以看出，爆心下方金字塔点阵芯层杆件承受的载荷超出了铝合金材料的断裂极限，出现大范围断裂失效。芯层阵列中心处的单胞杆件最先发生断裂，然后向外蔓延。而边缘处的芯层杆件发生局部扭曲、倒伏。图13(c)给出了工况4条件下1.05 ms时刻芯层中心局部区域的应力场，从图中可以看出，金字塔单胞杆件的未搭接部位为薄弱区，爆心对应区域的杆件发生断裂；而杆件相互嵌锁的十字形区域由于相互之间的加固作用，并未出现明显的断裂现象。
比较4种工况条件下复合点阵夹芯结构的结构响应云图可以看出，随着药量的增加、爆距的减小，即爆炸载荷量的增加，结构的破坏情况更加严重，除了出现典型的面芯脱粘失效模式外，芯层中心区域也会出现杆件断裂、压溃等破坏模式。此外，随着爆距的减小，即迎爆面越靠近炸药中心，上面板表面也出现了局部裂纹、铺层撕裂等破坏模式。
为了更准确地表达装备防护结构的易损特性，工程上往往需要建立结构的毁伤数据库。然而，综上可知，复合点阵结构的失效模式形式多样，且与载荷条件、结构参数密切相关，描述其毁伤效果并不容易。本文中尝试提出毁伤变量，基于毁伤变量抽象出复合点阵结构在冲击载荷下的毁伤函数。假定毁伤变量为
Rf ，则可以将毁伤函数写成关于各变量的多元函数：Rf=f(q,L,N,Score,⋯) (13) 式中：
N 为面板铺层的总层数，Score为芯子构型。爆炸冲击载荷作用下，复合点阵夹芯结构的上、下面板均会发生不同程度的变形。在装备设计过程中，应该尽可能减小面板变形，以避免对内部人员或元器件产生影响。爆炸冲击载荷下，结构面板的最大挠度在一定程度上反映了结构抵抗变形的能力。因此，选取面板最大挠度X作为描述结构毁伤程度的变量Rf。根据数值模型获得的复合点阵结构不同工况条件下的最大挠度数据，可以拟合出结构面板最大挠度X与炸药TNT当量q、爆炸距离L之间的二元函数关系式
X=f(q,L) ，并将其作为复合点阵夹芯结构爆炸冲击下的毁伤函数。复合点阵夹芯结构不同爆炸载荷下的毁伤情况由表4给出。表 4 复合点阵夹芯结构不同工况条件下的毁伤情况Table 4. Damage of composite lattice structures under different loading conditions芯子构型 爆炸距离L/mm TNT当量q/g 最大挠度X/mm 金字塔芯子 800 900 16.58 1 200 900 8.38 1 500 12.14 3 000 26.08 采用回归模型，将最大挠度
X 表示为关于TNT当量q 、爆炸距离L 的二元函数X=f(q,L) ，结合表4中的数据，可以得到：X=26.41−0.022L+0.0086q (14) 从式(14)可以看出，各项的系数代表该项对结构毁伤程度的影响权重，因此，可以看出金字塔型复合点阵夹芯结构受爆距的影响较大。根据式(14)可以绘制出复合点阵夹芯结构爆炸冲击下的二元毁伤函数图像，如图14所示。经计算，二元毁伤函数拟合的相关系数为0.982，这说明拟合效果较好。根据图14，一方面，从打击与破坏角度，可以获得复合点阵结构多种变量相关的毁伤区域；另一方面，从防护与承载角度，可以进一步得到结构的设计与优化空间。
图 14 复合点阵结构的二元毁伤函数图像Figure 14. Damage function related to explosive weight and explosion distance for the composite lattice structure6. 结 论
本文中采用爆炸冲击实验与基于复合材料三维渐进损伤准则的动态响应数值预报模型相结合的方法，系统地分析了金字塔型复合点阵夹芯结构在不同爆炸冲击载荷条件下的动态响应规律和吸能特性，获得的主要结论如下。
(1)碳纤维增强复合材料面板与金属芯层组合成的复合点阵夹芯结构具有轻质、高强、高抗爆、优异吸能等特点，在近距离强爆炸载荷作用下，结构整体上基本保持完好。初步探索了复合点阵结构多种影响参数相关的毁伤函数，给出了结构的可行设计域，可根据工程需求，进一步开展优化设计。
(2)基于有限元方法的复合点阵夹芯结构爆炸冲击响应数值预报模型具有较高的准确度，模型获得的结构失效模式与爆炸冲击实验结果基本吻合，利用数值模型可以模拟爆炸载荷在结构内部的传递过程和能量转化规律。
(3)在强爆炸冲击载荷作用下，复合点阵夹芯结构出现了局部失效现象，失效模式随药量、爆距变化呈现出多样化特征，其中面芯脱粘失效为主要破坏模式，复合材料面板与金属芯层之间的粘接强度是影响结构抗爆性能的关键因素，如何提高面芯粘接性能是复合点阵结构的重要研究方向。
-
图 1 复合点阵夹芯结构示意图及设计参数
Figure 1. Schematic diagram and design parameters of a composite lattice sandwich structure
图 5 复合点阵结构的爆炸冲击响应有限元分析模型
Figure 5. The finite element model for predicting explosion shock response of composite lattice structure
图 6 复合点阵结构爆炸冲击实验后的样件照片及局部失效情况
Figure 6. Photographs and local failures of composite lattice structure specimens after explosion shock experiments
图 7 复合点阵夹芯结构爆炸冲击载荷下的位移变化云图
Figure 7. Strain response of composite lattice sandwich structure under explosion shock loadings
图 8 复合点阵夹芯结构爆炸冲击载荷下的位移变化曲线
Figure 8. Strain response curves of composite lattice sandwich structure under explosion shock loadings
图 9 复合点阵夹芯结构爆炸冲击载荷下的应力变化规律
Figure 9. Stress response of composite lattice sandwich structure under explosion shock loadings
图 10 复合点阵夹芯结构上、下面板单个铺层的应力变化规律
Figure 10. Cloud maps of stress distribution between layers of composite lattice structure panels
图 11 复合点阵夹芯结构爆炸冲击过程中的能量转化关系
Figure 11. Energy absorption curves of composite lattice sandwich structure exposed to explosion loadings
图 12 复合点阵夹芯结构在爆炸冲击载荷下的失效模式
Figure 12. Failure modes of composite lattice sandwich structure under explosion loadings
图 13 复合点阵夹芯结构不同载荷条件下的结构响应
Figure 13. Structural responses of composite lattice sandwich structures under different loadings
图 14 复合点阵结构的二元毁伤函数图像
Figure 14. Damage function related to explosive weight and explosion distance for the composite lattice structure
表 1 金字塔型芯子代表性单胞主要设计参数
Table 1. Main design parameters of representative structure cell of pyramidal truss core
h/mm p/mm t1/mm s/mm t2/mm θ/ (°) ha/mm hb/mm 15 35.75 1.59 8 1.5 45 3 1.5 
下载: 导出CSV
表 2 碳/环氧复合材料的基本材料性能参数
Table 2. Mechanical properties of carbon/epoxy composite materials
密度/(kg·m−3) 拉压模量/GPa 泊松比 剪切模量/GPa 拉伸强度/MPa 压缩强度/MPa 面内剪切强度/MPa 层间剪切强度/MPa 1600 123 (0°) 0.41 (0°) 4.8 1400 (0°) 850 (0°) 60 16 8.3 (90°) 0.26 (90°) 18 (90°) 96 (90°) 注：0°是指力作用在纤维方向；90°是指力作用在垂直纤维方向。
下载: 导出CSV
表 3 爆炸实验载荷工况
Table 3. Loading conditions of explosion experiments
工况 TNT当量/g 爆炸距离/m 比例距离/(m·kg−1/3) 1 900 1.20 1.24 2 1000 1.86 1.86 3 1500 1.20 1.05 4 3000 1.20 0.83
下载: 导出CSV
表 4 复合点阵夹芯结构不同工况条件下的毁伤情况
Table 4. Damage of composite lattice structures under different loading conditions
芯子构型 爆炸距离L/mm TNT当量q/g 最大挠度X/mm 金字塔芯子 800 900 16.58 1 200 900 8.38 1 500 12.14 3 000 26.08
下载: 导出CSV
-
[1] 金键, 朱锡, 侯海量, 等. 大型舰船在水下接触爆炸下的毁伤与防护研究综述 [J]. 爆炸与冲击, 2020, 40(11): 111401. DOI: 10.11883/bzycj-2020-0105.JIN J, ZHU X, HOU H L, et al. Review on damage and protection of large ships under underwater contact explosion [J]. Explosion and Shock Waves, 2020, 40(11): 111401. DOI: 10.11883/bzycj-2020-0105. [2] 张军, 黄含军, 王军评, 等. 炸药驱动式爆炸管的载荷计算 [J]. 装备环境工程, 2021, 18(5): 21–27. DOI: 10.7643/issn.1672-9242.2021.05.004.ZHANG J, HUANG H J, WANG J P, et al. Load calculation of explosive driven explosive tube [J]. Equipment Environment Engineering, 2021, 18(5): 21–27. DOI: 10.7643/issn.1672-9242.2021.05.004. [3] 马福临, 杨娜娜, 赵天佑, 等. 冲击波-破片群联合作用下舰船复合材料结构近场动力学损伤模拟 [J]. 爆炸与冲击, 2022, 42(3): 033304. DOI: 10.11883/bzycj-2021-0080.MA F L, YANG N N, ZHAO T Y, et al. Near field dynamic damage simulation of ship composite structure under the combined action of shock wave and fragment group [J]. Explosion and Shock Waves, 2022, 42(3): 033304. DOI: 10.11883/bzycj-2021-0080. [4] 王志鹏, 李海波, 韦冰峰, 等. 嵌锁式CFRP方形蜂窝夹芯梁低速冲击响应及失效机理[J] [J]. 爆炸与冲击, 2022, 42(7): 073102. DOI: 10.11883/bzycj-2021-0525.WANG Z P, LI H B, WEI B F, et al. Low-velocity impact response and failure mechanism of CFRP sandwich beams with a square honeycomb core fabricated by the interlocking method [J]. Explosion and Shock Waves, 2022, 42(7): 073102. DOI: 10.11883/bzycj-2021-0525. [5] 陈东, 吴永鹏, 李忠盛, 等. 轻质高强多功能点阵夹层结构研究进展 [J]. 装备环境工程, 2020, 17(4): 77–84. DOI: 10.7643/issn.1672-9242.2020.04.013.CHEN D, WU Y P, LI Z S, et al. Research progress of light, high strength and multi-function lattice sandwich structures [J]. Equipment Environment Engineering, 2020, 17(4): 77–84. DOI: 10.7643/issn.1672-9242.2020.04.013. [6] 孙晓旺, 陶晓晓, 王显会, 等. 负泊松比蜂窝材料抗爆炸特性及优化设计研究 [J]. 爆炸与冲击, 2020, 40(9): 095101. DOI: 10.11883/bzycj-2020-0011.SUN X W, TAO X X, WANG X H, et al. Study on anti-explosion characteristics and optimum design of honeycomb material with negative Poisson’s ratio [J]. Explosion and Shock Waves, 2020, 40(9): 095101. DOI: 10.11883/bzycj-2020-0011. [7] 吴艳青, 刘彦, 黄风雷, 等. 爆炸力学理论及应用[M]. 1版. 北京: 北京理工大学出版社, 2021: 1–20. [8] 程帅, 刘文祥, 童念雪, 等. 爆炸载荷下飞机典型加筋结构毁伤特性 [J]. 爆炸与冲击, 2021, 41(1): 013302. DOI: 10.11883/bzycj-2020-0077.CHENG S, LIU W X, TONG N X, et al. Damage characteristics of typical stiffened aircraft structures under explosive loads [J]. Explosion and Shock Waves, 2021, 41(1): 013302. DOI: 10.11883/bzycj-2020-0077. [9] WU Q Q, GAO Y, XIONG J. Quasi-static mechanical properties of composite lattice sandwich structures with enhanced face panels [J]. European Journal of Mechanics: A/Solids, 2023, 97: 104808. DOI: 10.1016/j.euromechsol.2022.104808. [10] EVANS A G, HUTCHINSON J W, FLECK N A, et al. The topological design of multifunctional cellular metals [J]. Progress in Materials Science, 2001, 46(3/4): 309–327. DOI: 10.1016/S0079-6425(00)00016-5. [11] LI X, WANG Z H, ZHU F, et al. Response of aluminium corrugated sandwich panels under air blast loadings: experiment and numerical simulation [J]. International Journal of Impact Engineering, 2014, 65: 79–88. DOI: 10.1016/j.ijimpeng.2013.11.002. [12] QI C, YANG S, YANG L J, et al. Blast resistance and multi-objective optimization of aluminum foam-cored sandwich panels [J]. Composite Structures, 2013, 105: 45–57. DOI: 10.1016/j.compstruct.2013.04.043. [13] 王涛, 余文力, 秦庆华, 等. 爆炸载荷下泡沫铝夹芯板变形与破坏模式的实验研究 [J]. 兵工学报, 2016, 37(8): 1456–1463. DOI: 10.3969/j.issn.1000-1093.2016.08.017.WANG T, YU W L, QIN Q H, et al. Experimental investigation into deformation and damage patterns of sandwich plates with aluminum foam core subjected to blast loading [J]. Acta Armamentarii, 2016, 37(8): 1456–1463. DOI: 10.3969/j.issn.1000-1093.2016.08.017. [14] YUNGWIRTH C J, WADLEY H N G, O’CONNOR J H, et al. Impact response of sandwich plates with a pyramidal lattice core [J]. International Journal of Impact Engineering, 2008, 35(8): 920–936. DOI: 10.1016/j.ijimpeng.2007.07.001. [15] WADLEY H, DHARMASENA K, CHEN Y, et al. Compressive response of multilayered pyramidal lattices during underwater shock loading [J]. International Journal of Impact Engineering, 2008, 35(9): 1102–1114. DOI: 10.1016/j.ijimpeng.2007.06.009. [16] DHARMASENA K P, WADLEY H N G, LIU T, et al. The dynamic response of edge clamped plates loaded by spherically expanding sand shells [J]. International Journal of Impact Engineering, 2013, 62: 182–195. DOI: 10.1016/j.ijimpeng.2013.06.012. [17] 亓昌, 郝鹏程, 舒剑, 等. 金字塔型点阵材料夹芯板抗爆性能仿真与优化 [J]. 振动与冲击, 2019, 38(16): 245–252. DOI: 10.13465/j.cnki.jvs.2019.16.035.QI C, HAO P C, SHU J, et al. Simulation and optimization for blast-resistant performances of pyramidal lattice cored sandwich panels [J]. Journal of Vibration and Shock, 2019, 38(16): 245–252. DOI: 10.13465/j.cnki.jvs.2019.16.035. [18] 彭航. 侵爆战斗部对典型建筑物内爆毁伤效应与评估方法研究[D]. 南京: 南京理工大学, 2021. DOI: 10.27241/d.cnki.gnjgu.2021.000774.PENG H. Study on damage effect and evaluation method of penetration warhead on typical building internal explosion [D]. Nanjing, Jiangsu, China: Nanjing University of Science and Technology, 2021. DOI: 10.27241/d.cnki.gnjgu.2021.000774. [19] 冯晓伟, 卢永刚, 李永泽. 飞机目标在爆炸冲击波作用下的毁伤效应评估方法 [J]. 高压物理学报, 2019, 33(4): 045101. DOI: 10.11858/gywlxb.20180687.FENG X W, LU YG, LI Y Z. Damage effect evaluation method of aircraft target under explosion shock wave [J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 045101. DOI: 10.11858/gywlxb.20180687. [20] HAN L, HAN Q, GE Y X, et al. Vulnerability assessment of combat aircraft to blast loading [J]. Proceedings of the Institution of Mechanical Engineers, 2019, 233(2): 604–615. DOI: 10.1177/0954410017740162. [21] 胡榕,姜春兰,卢广照,等. 机场跑道内爆炸毁伤效应及工程化函数模型 [J]. 兵工学报, 2023, 44(4): 929–939. DOI: 10.12382/bgxb.2022.1220.HU R, JIANG C L, LU G Z, et al. Explosion damage effect and engineering function model in airport runway [J]. Acta Armamentarii, 2023, 44(4): 929–939. DOI: 10.12382/bgxb.2022.1220. [22] WEI X Y, XIONG J, WANG J, et al. New advances in fiber-reinforced composite honeycomb materials [J]. Science China Technological Sciences, 2020, 63(8): 1348–1370. DOI: 10.1007/s11431-020-1650-9. [23] LI Z B, GAO Y, WEI X Y, et al. Fabrication and failure mechanisms of all-composite honeycomb sandwich cylinder under the axial compression [J]. Composites Part A: Applied Science and Manufacturing, 2022, 161: 107075. DOI: 10.1016/j.compositesa.2022.107075. [24] ZHAO H N, FANG H Y, ZHAO X H. Experimental and numerical investigation on dynamic behaviors of glass fiber reinforced polymer plates under explosion loadings [J]. International Journal of Impact Engineering, 2023, 171: 104362. DOI: 10.1016/j.ijimpeng.2022.104362. [25] 李伟, 王鹏, 李佳, 等. 高速冲击下多胞材料细观结构中波传播规律 [J]. 装备环境工程, 2021, 18(5): 87–93. DOI: 10.7643/issn.1672-9242.2021.05.013.LI W, WANG P, LI J, et al. Wave propagation in mesostructure of multicellular materials under high speed impact [J]. Equipment Environment Engineering, 2021, 18(5): 87–93. DOI: 10.7643/issn.1672-9242.2021.05.013. 期刊类型引用(4)
1. 张昊，段晟昱，符瑞，彭小洋，雷红帅. 异形点阵夹芯结构底板抗爆炸冲击设计与优化研究. 应用数学和力学. 2025(05): 582-590 . 
百度学术2. 陈洋，王肇喜，翟师慧，盛鹏，王者蓝，朱明亮. 3D打印点阵夹芯结构冲击损伤的近场动力学模拟. 爆炸与冲击. 2024(03): 146-160 . 
本站查看3. 黄斌，范启东，梁剑烽. 薄壁仿生护具微冲击响应智能预测方法仿真. 计算机仿真. 2024(10): 438-442 . 百度学术4. 陈清华，吕可，王德俊，王建刚，王建业，冯鹏. 爆炸载荷下实验舱功能梯度防爆结构的性能. 中国机械工程. 2023(21): 2568-2576 . 百度学术其他类型引用(3)
-
下载:
百度学术
计量
- 文章访问数: 463
- HTML全文浏览量: 145
- PDF下载量: 206
- 被引次数: 7