Simulation analysis on impact resistance of aluminum foam sandwich structures using peridynamics
-
摘要: 针对某光学舱所采用的泡沫铝夹层防护结构在破片冲击下的抗冲击性能问题，采用Monte-Carlo方法创建了泡沫铝结构的二维细观模型，在常规态型近场动力学理论中引入了Mises屈服准则和线性各向同性强化模型，建立了近场动力学塑性本构的数值计算框架。基于近场动力学计算程序模拟了低速冲击作用下泡沫铝夹层结构的塑性变形以及有机玻璃背板的裂纹扩展形态，分析了泡沫铝芯材孔隙率对该夹层结构抗冲击性能和损伤模式的影响规律。结果表明：泡沫铝夹层结构良好的塑性变形能力是其发挥缓冲与防护作用的主要因素，并且在一定范围内，泡沫铝芯材孔隙率越高，则夹层结构具有更好的抗冲击性能；当泡沫铝孔隙率从0.4提升到0.7时，泡沫铝对冲击物的动能吸收率从90%提高到99%；模拟结果与实验结果具有较好的一致性，验证了模拟结果的准确性和分析结论的有效性。通过数值模拟，预测了有机玻璃背板的裂纹扩展形态，发现提高泡沫铝的孔隙率能获得更好的防护效果。Abstract: Under impact, aluminum foam undergoes significant plastic deformation, and the kinetic energy of the impactor is dissipated in the process, thereby protecting the structure from damage. The failure modes of aluminum foam sandwich structures under impact are complex, involving plastic deformation, panel failure, and cracking of the bonding interface. Traditional numerical simulation methods are difficult to solve these discontinuous problems. Peridynamics is a non-local numerical method that describes the mechanical behavior of materials by solving spatial integral equations. It has unique advantages in solving crack propagation, material failure, progressive damage of composite materials, and multi-scale problems. Although the basic bond-based peridynamic theory cannot describe plasticity, the ordinary state-based peridynamic method decouples distortion and dilation and can easily simulate the plastic deformation of materials. Therefore, based on ordinary state-based peridynamics, the Mises yield criterion and the linear isotropic hardening model were introduced to study the factors affecting the impact resistance of aluminum foam sandwich structures. Two-dimensional mesoscopic models of aluminum foam sandwich structure were established by the Monte-Carlo method and impact was simulated using the peridynamic method. The influence of the porosity of aluminum foam on the impact resistance and damage mode of the sandwich structure was analyzed. The results show that the good plastic deformation ability of aluminum foam sandwich structure is the main factor for its buffering and protection, and within a certain range, the higher the porosity of aluminum foam core, the better impact resistance of the sandwich structure. When the porosity of aluminum foam increases from 0.4 to 0.7, the kinetic energy absorption rate of aluminum foam to the impactor increases from 90% to 99%. The simulation results are in good agreement with the experimental results, which verifies the accuracy of the simulation results and the effectiveness of the analysis conclusions. The numerical simulation predicts the crack propagation morphology of the plexiglass backplate, and the results show that improving the porosity of aluminum foam can obtain a better protection effect.
-
泡沫铝是由铝或者铝合金基体与内部孔洞组成的新型结构与功能一体化材料，内部存在大量孔洞。与传统的合金材料相比，泡沫铝具有更低的体积密度、更强的比强度和比刚度，并且具有良好的缓冲吸能、阻尼减振、电磁屏蔽和隔音降噪等性能。因此，在航空航天、武器装备、建筑结构和交通运输等领域，泡沫铝具有广阔的应用前景。在实际使用过程中，泡沫铝夹层结构（aluminum foam sandwich，AFS）是最常用的结构形式，理论模型[1-3]表明这类夹芯结构具有优良的抗冲击性能。泡沫铝的孔隙率一般为30%～90%，孔径为0.3～7.0 mm。孔隙率或孔径以及孔洞分布形式的改变，对AFS的学性能有显著影响，研究AFS是力学及材料学科的热点问题之一。郭亚周等[4]采用一级空气炮开展冲击实验，研究了各种形状弹体冲击作用下不同泡沫铝夹层结构的动力性能；张博一等[5]采用一级轻气炮驱动泡沫铝子弹，研究了泡沫铝子弹冲击下泡沫铝夹层结构的动态响应。空气炮是研究AFS抗冲击性能最常用的实验手段，但是实验成本较高，难以广泛应用，因此数值模拟成为目前的重要方法。邓旭辉等[6]采用数值方法，分析了双层泡沫铝夹芯板的抗冲击性能；张永康等[7]采用LS-DYNA进行数值计算，对不同构型泡沫铝夹层结构进行了冰雹冲击分析，研究了中面板位置对缓冲吸能性能的影响；夏志成等[8]设计了落锤实验对泡沫铝夹芯板进行研究，并应用数值模拟软件还原夹芯板冲击过程，计算了泡沫铝结构的吸能效果；Cui 等[9]采用工业CT扫描技术，重建了泡沫铝夹芯结构三维细观有限元模型，分析了不同构型弹道冲击下泡沫铝夹芯结构的失效模式，并在实验和数值观测结果的基础上，建立了预测弹丸残余速度的理论模型。在防爆结构领域，泡沫铝也引起相关学者的注意。苏兴亚等[10]采用弹道冲击摆系统，开展了爆炸载荷下分层梯度泡沫铝夹芯板的变形、失效模式和抗冲击性能的实验研究；王涛等[11]开展了不同爆炸载荷作用下泡沫铝夹芯板变形的实验研究，获得了泡沫铝夹芯结构背板中心点挠度与爆炸冲量之间的关系。
现有的有限元方法能有效计算泡沫铝夹层结构的力学性能，预测其准静态及动态响应，但是传统的有限元方法是从连续介质力学出发，用偏微分方程来描述材料的力学行为，在处理泡沫铝夹层结构界面开裂、芯材溃裂、背板材料裂纹扩展等非连续问题时存在瓶颈。近场动力学（peridynamic，PD）理论[12-14]是一种通过求解空间积分方程来描述材料力学行为的非局部数值方法，因它在解决裂纹扩展、材料破坏、复合材料渐进损伤以及多尺度问题上独特的优势，越来越多的学者将PD方法引入到各种领域解决传统方法难以处理的问题。杨娜娜等[15]采用近场动力学，对破片冲击下复合材料的损伤进行了模拟，分析了冲击速度、复合材料铺层方式等对结构损伤的影响；刘宁等[16]采用短程斥力模型，描述了碰撞过程，模拟了冲击压缩条件下单裂纹圆孔板动态破坏行为；熊伟鹏等[17]采用近场动力学方法，模拟了冰球碰撞过程，并与实验对比取得良好的效果。基于键的近场动力学方法不支持塑性材料本构模型，多用于脆性材料破坏问题。而常规态型近场动力学方法解耦了几何形状变形和体积变形，可以实现材料塑性变形[18]，因此适用于泡沫铝这类具有良好塑性的材料。目前，我国在近场动力学塑性模型方面的研究还较罕见。
本文中，为了明确孔隙率对泡沫铝夹层结构的抗冲击性能的影响，采用Monte-Carlo方法生成泡沫铝的随机模型，基于常规态型近场动力学理论，引入Mises屈服准则和线性各向同性强化模型来描述铝的塑性行为，对泡沫铝夹层结构进行冲击计算，分析孔隙率分别为0.7、0.6、0.5和0.4的泡沫铝夹层结构的缓冲效能以及破坏形态，以期为此类问题提供借鉴。
1. 常规态型近场动力学
1.1 运动方程
如图1所示，质点x与其半径为δ的邻域Hx内其他质点存在相互作用。当材料发生变形时，质点x与其邻域内的另一质点x'分别移动到y和y'，质点x'对质点x产生力密度t，力密度矢量t的方向与质点之间的连线重合，其大小与质点x'的变形以及质点之间的相对位移相关。同理，质点x'也受到其邻域Hx'内其他质点的影响，因此质点x也对质点x'产生力密度矢量t'，其方向与t平行，大小与质点x'的变形以及质点之间的相对位移相关。质点x的运动状态是在与邻域Hx内所有其他质点相互影响下产生，在邻域内对力密度进行积分，可以得到关于质点x的拉格朗日方程[12]：
\rho({\boldsymbol{x}}) \ddot{{\boldsymbol{u}}}({\boldsymbol{x}}, t)=\int_{H_{x}}\left({\boldsymbol{t}}\left({\boldsymbol{u}}^{\prime}-{\boldsymbol{u}}, {\boldsymbol{x}}^{\prime}-{\boldsymbol{x}}, t\right)-{\boldsymbol{t}}^{\prime}\left({\boldsymbol{u}}-{\boldsymbol{u}}^{\prime}, {\boldsymbol{x}}-{\boldsymbol{x}}^{\prime}, t\right)\right) {\rm d} H+{\boldsymbol{b}}({\boldsymbol{x}}, {{t}}) (1) 近场动力学的运动方程不涉及求解空间坐标的偏微分方程，取而代之的是通过求解材料内部质点影响范围内的空间积分来描述材料宏观变形。常用的数值计算方法是将材料离散成空间中均匀分布的一系列质点，质点坐标代表了其所在区域的体积，于是式(1)可写成如下离散形式[13]：
\rho\left({\boldsymbol{x}}_{k}\right) \ddot{{\boldsymbol{u}}}\left({\boldsymbol{x}}_{k}, t\right)=\sum_{j=1}^{N}\left[{\boldsymbol{t}}_{k, j}\left({\boldsymbol{u}}_{j}-{\boldsymbol{u}}_{k}, {\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right)-{\boldsymbol{t}}_{j,k}\left({\boldsymbol{u}}_{k}-{\boldsymbol{u}}_{j}, {\boldsymbol{x}}_{k}-{\boldsymbol{x}}_{j}, t\right)\right]V_{j}+{\boldsymbol{b}}{({\boldsymbol{x}}_k}, {{t}}) (2) 式中：xk和xj为质点坐标，并且xj在xk的影响范围内；N为质点xk影响范围内的质点总数；ρ(xk)为质点的密度；uk为质点xk的位移向量；
\ddot{{\boldsymbol{u}}}\left({\boldsymbol{x}}_{k}, t\right) 为质点xk的加速度矢量；b(xk, t)为作用在质点xk上的外体力密度矢量；Vj为质点xj的体积。1.2 弹脆性材料本构模型
质点间相互作用的力密度可以通过近场动力学应变能密度导出[12-14]：
t_{k, j}=2\left(a_{k}-a_{\mu}\right) \frac{\delta d}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} \varLambda_{k, j} \theta_{k}+2 \delta b s_{k, j} (3) 式中：θk和sk, j分别为质点xk的体积应变和质点xk、xj之间的伸长率。θk和sk, j分别定义为：
\theta_{k}=d \sum_{j=1}^{N} \delta s_{{k}, j} \varLambda_{{k}, j} V_{j} (4) s_{k, j}=\frac{\left|{\boldsymbol{y}}_{j}-{\boldsymbol{y}}_{k}\right|-\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} (5) 参数Λk, j定义为：
\varLambda_{k, j}=\frac{{\boldsymbol{y}}_{j}-{\boldsymbol{y}}_{k}}{\left|{\boldsymbol{y}}_{j}-{\boldsymbol{y}}_{k}\right|} \cdot \frac{{\boldsymbol{x}}_{j} -{\boldsymbol{x}}_{k}} {\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} (6) 通过考虑简单加载情况下的应变能密度，可以将式(3)中的近场动力学参数aκ、aμ、b和d与传统连续介质力学中的体积模量κ和剪切模量μ建立联系，推导出[14, 18]：
\left\{ \begin{array}{l} a_{\kappa}=\dfrac{1}{2} \kappa,\;\;\; a_{\mu}=\dfrac{5 \mu}{6}, \;\;\;b=\dfrac{15 \mu}{2 {\text{π}} \delta^{4}},\;\;\;\;\; d=\dfrac{9}{4 {\text{π}} \delta^{4}} \;\;\;\;\;\;\;\text { 3D } \\ a_{\kappa}=\kappa, \;\;\;a_{\mu}=2 \mu, \;\;\;b=\dfrac{6 \mu}{{\text{π}} h \delta^{5}}, \;\;\;\;\;d=\dfrac{2}{{\text{π}} h \delta^{3}}\;\;\;\;\;\;\;\;\;\;\;\;\mathrm{2D} \\ a_{\kappa}=0, \;\;\;a_{\mu}=0, \;\;\;b=\dfrac{E}{2 \delta^{3} A}, \;\;\;\;\; d=\dfrac{1}{2 \delta^{2} A} \;\;\;\;\;\;\;\;\;\;\;\;\text { 1D } \end{array}\right. (7) 近场动力学中，材料损伤失效是通过截断物质点之间的相互作用描述的，并且截断之后不可恢复，因此，在力密度和体积应变中引入一个与时间相关的标量函数
\varphi 进行修正，其定义为：\varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right)=\left\{ \begin{array}{lll} 1\;\;\;\;\;s_{k, j}{\text{＜}}s_{\text{c}} \; \;({完好}) \\ 0\;\;\;\;\;s_{k, j}{\text{＞}}s_{\text{c}} \;\; {(断裂 }) \end{array}\right. (8) 式中：sc为临界伸长率。当质点之间的伸长率超过临界伸长率，则质点间的相互作用永久性破坏，这意味着这2个质点之间的微势能因裂纹面产生而消失。于是，得到添加损伤修正的力密度和体积应变：
t_{k, j}=2\left(a_{\kappa}-a_{\mu}\right) \frac{\delta d}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} \varLambda_{k, j} \theta_{k}+2 \delta b \varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right) s_{k, j} (9) \theta_{k}=d \delta \sum_{j=1}^{N} \varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right) s_{k, j} \varLambda_{k, j} V_{j} (10) 对穿过裂纹面的微势能在区域内进行积分，与材料断裂的临界能量释放率建立联系，可以得到临界伸长率的表达式为[12-14]：
s_{\text{c}}=\left\{ \begin{array}{l} \sqrt{\dfrac{G_{\text{c}}}{\left[3\mu+\left(\dfrac{3}{4}\right)^4\left(\kappa-\dfrac{5\mu}{3}\right)\right]\delta}}\;\;\;\;\;\;\;\;\;\;{\rm{3D}}\\ \sqrt{\dfrac{G_{\text{c}}}{\left[\dfrac{6}{{\text{π}}}+\dfrac{16}{9{\text{π}}^2}(\kappa-2\mu) \right]\delta}}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\rm{2D}} \end{array}\right. (11) 实际上，对于复杂的材料目前仍无法简单地计算其临界伸长率，此时便需要采用实验标定。通过统计质点影响范围内截断的相互作用占初始相互作用总数的比例，可表征局部的损伤程度：
\phi\left({\boldsymbol{x}}_{k}, t\right)=1-\frac{\displaystyle\sum_{j=1}^{N} \varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right) V_{j}}{\displaystyle\sum_{j=1}^{N} V_{j}} (12) 1.3 近场动力学塑性理论
塑性变形过程与加载路径相关，并且是不可逆的，通常需要采用增量的方式描述载荷与变形之间的关系。将材料变形分解成可恢复的弹性部分和不可恢复的塑性部分，于是可以将质点之间的伸长率增量Δsk, j分解为：
\Delta s_{k, j}=\Delta s_{k, j}^{{\rm{e}}}+\Delta s_{k, j}^{\rm p} (13) 式中：
\Delta s_{k, j}^{\rm e} 和\Delta s_{k, j}^{\rm p} 分别为质点之间的弹性伸长率增量和塑性伸长率增量。当采用步进的方式加载时，第n+1个加载步中质点xk和xj之间的总伸长率为：s_{k, j}^{(n+1)}=s_{k, j}^{(n)}+\Delta s_{k, j}^{(n+1)} =\left(s^{{\rm{e}}(n)}_{k,j}+\Delta s^{{\rm{e}}(n+1)}_{k,j}\right)+\left(\Delta s^{{\rm{p}}(n+1)}_{k,j}\right) (14) 式中：
s_{k, j}^{\rm e} 、s_{k_, j}^{\rm p} 和s_{k, j} 分别为质点之间的弹性伸长率、塑性伸长率和总伸长率。假设在第n+1个加载步中，材料仅存在弹性变形增量，塑性变形增量\Delta s_{k, j}^{{\rm{p}}(n+1)}=0 ，那么弹性伸长率的试算值为：s_{k, j({\rm{trial}})}^{{\rm{e}}(n+1)}=s_{k, j}^{(n+1)}-s_{k, j}^{{\rm p}(n)} (15) 此时，根据式(3)，可以得到力密度的试算值为：
t^{(n+1)}_{k,j({\rm{trial}})}=2(a_{\kappa} -a_{\mu})\frac{\delta d}{|{\boldsymbol{x}}_j-{\boldsymbol{x}}_k|}\varLambda_{k,j}\theta_k \left(s^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})}\right)+2\delta bs^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})} (16) 根据Mises屈服准则，当材料内一点的应力偏量第二不变量达到某一值时，该点便开始进入塑性状态，当考虑线性各向同性强化时，可以建立屈服函数[18]：
F_{k}^{(n+1)}=b \sum_{j=1}^{N} \frac{\delta}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|}\left(s_{k, j}^{{\rm{e}}(n+1)}\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|\right)^{2} V_{j}-a_{\mu}\left(\theta_{k}^{{\rm{e}}(n+1)}\right)^{2}-G\left(\overline s_{k}^{{\rm{p}} (n+1)}\right) (17) 硬化函数定义为：
G\left(\bar{s}_{k}^{\rm p}\right)=\frac{\left(\sigma_{\rm y}+E_{\rm t} \bar{s}_{k}^{\rm p}\right)^{2}}{6 \mu} (18) 式中：
\bar{s}_{k}^{\rm p} 和Et分别为材料的等效塑性应变和切线模量。将弹性伸长率的试算值
s^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})} 代入屈服函数，当F_{k}^{(n+1)} {\text{≤}} 0 时，该加载步位于屈服面内，塑性伸长率增量\Delta s_{k, j}^{{\rm p}(n+1)}=0 ，弹性伸长率的试算值s^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})} 即为真实发生的弹性伸长率。否则，加载状态位于屈服面上方，\Delta s_{k, j}^{\mathrm{p} (n+1)} \neq 0 ，需要确定真实发生的弹性伸长率以使F_{k}^{(n+1)} 退回到屈服面内：\left\{\begin{array}{l}s^{{\rm e}(n+1)}_{k,j}=s^{{\rm e}(n+1)}_{k,j({\rm{trial}})}-\Delta s^{{\rm p}(n+1)}_{k,j}\\ \Delta s^{{\rm p}(n+1)}_{k,j}=C_k\left(\dfrac{1}{2\delta b}|{\boldsymbol{x}}_j-{\boldsymbol{x}}_k|t^{(n+1)}_{k,j({\rm{trial}})}+\dfrac{d}{b}(a_{\kappa}-a_{\mu})\varLambda_{k,j} \theta_k\left( s^{{\rm e}(n+1)}_{k,j({\rm{trial}})} \right) \right) \end{array} \right. (19) 式中：Ck为一个大于0的系数，使得将上式中的弹性伸长率代入屈服函数后
F_{k}^{(n+1)} \to 0 。Ck可以采用牛顿迭代法确定，最终得到此加载步材料各点的塑性伸长率和弹性伸长率分别为：\left\{\begin{array}{l}s^{{\rm p}(n+1)}_{k,j}=s^{{\rm p}(n+1)}_{k,j}-\Delta s^{{\rm p}(n+1)}_{k,j}\\ s^{{\rm e}(n+1)}_{k,j}=s^{(n+1)}_{k,j}- s^{{\rm p}(n+1)}_{k,j} \end{array} \right. (20) 2. 泡沫铝夹层结构模型
如图2所示，光学舱防护结构采用的是泡沫铝夹层复合材料，由上层铝蒙皮面板、中间泡沫铝夹层以及底层有机玻璃背板组成。泡沫铝的孔隙率P有0.7、0.6、0.5和0.4等4种，所选用泡沫铝主孔径为1～6 mm。泡沫铝夹层结构水平放置在刚性平台上，上方有一个直径为25 mm的刚性小球以20 m/s的速度向下冲击。为了便于建模并提高计算效率，将模型简化成二维结构，简化成二维的刚性小球质量为3.5×10−3 kg。采用图3所示的Monte-Carlo方法随机生成泡沫铝模型，分3级进行胞孔投放，第1级孔径为4～6 mm，第2级孔径为2～4 mm，第3级孔径为1～2 mm，先投放大孔径后投放小孔径，以达到提高投放效率并尽可能获得更高的孔隙率的目的。将模型离散成均匀分布的球形粒子，为确保数值计算精度，需保证最小胞孔孔径大于近场范围，确定粒子直径Δx=0.2 mm，近场半径δ=3Δx。最终建立的4种不同孔隙率泡沫铝夹层复合结构的近场动力学离散模型如图4所示。
图 2 破片冲击泡沫铝夹层结构模型Figure 2. Schematic diagram of an aluminum foam sandwich impacted by a fragment
图 3 Monte-Carlo方法生成泡沫铝算法流程Figure 3. Algorithm flow chart of the Monte-Carlo method to generate aluminum foam
图 4 泡沫铝夹层结构的近场动力学离散模型Figure 4. Discrete models of aluminum foam sandwiches for peridynamics泡沫铝夹层结构相邻层之间用环氧树脂进行粘接，计算中不考虑环氧树脂的质量，仅通过环氧树脂的力学性能定义相邻材料界面之间的近场动力学参数，并对材料界面及边界处的参数进行修正[12-14, 18]。实验中，以20 m/s速度冲击时，铝板蒙皮及泡沫铝夹层主要表现为塑性变形，并且未观察到明显破裂，因此在这种工况下可以采用近场动力学塑性模型描述铝材料。有机玻璃和环氧树脂表现为脆性破坏，因此用近场动力学弹脆性模型描述。计算中采用的各材料力学参数如表1所示。
表 1 材料参数Table 1. Material parameters材料 ρ/(kg∙m−3) E/GPa ν σy/MPa Et/MPa sc 铝 2700 69.5 0.33 127 586 有机玻璃 1190 3.6 0.4 0.047 环氧树脂 3 0.37 0.047 3. 数值计算结果及分析
3.1 泡沫铝夹层结构变形形态
图5显示了泡沫铝夹层结构的变形过程，云图颜色表示各点竖直方向位移uy的绝对值：在冲击作用下，先是正对冲击物的部位铝面板发生局部凹陷，凹陷深度不断增大，内部泡沫铝逐渐发生大范围变形，随着泡沫铝变形范围进一步增大，底部有机玻璃背板在上层泡沫铝带动下发生弯曲，结构呈现左右两端向上翘起的整体变形。分析可知，这是由于泡沫铝内部的大量孔洞能大幅削减应力从冲击部位向周围的扩散范围，因此在冲击初期局部变形集中在结构与冲击物的接触部位，随着泡沫铝的孔洞被压实，结构内部应力积累到一定程度才导致发生整体变形。这是泡沫铝能起到缓冲作用的重要原因之一。
此外，影响夹层结构抗冲击性能的另一个主要因素是泡沫铝芯材塑性变形耗能[8]。冲击过程中，随着上层的泡沫铝塑性变形范围增大，冲击物的动能逐渐被耗散，最终使得作用在背板上的冲击力极大缓解，从而达到保护底部背板的目的。由此可认为，冲击作用下泡沫铝芯材的塑性变形程度与范围越大，对背板的保护作用越强，则泡沫铝夹层结构抗冲击性能更好。图6为具有不同孔隙率芯材的夹层结构在冲击结束后的塑性变形形态，可以看出，结构上越靠近冲击物的部位塑性变形程度越大，塑性变形呈以冲击凹坑为中心的放射状分布，并且孔隙率越高，冲击形成的凹坑越深，塑性变形范围越大。
3.2 冲击物的速度与加速度分析
通过分析冲击物在冲击过程中的速度变化以及加速度变化过程，可以了解冲击过程中动能传递规律，间接了解泡沫铝芯材的缓冲作用。图7是冲击物冲击不同孔隙率泡沫铝夹层结构时的速度v变化曲线，孔隙率越高，速度变化曲线越平缓。此外，根据表2的分析结果，孔隙率越高，冲击物的反弹速度越低，冲击结束时冲击物具有更低的残余动能，这说明高孔隙率的泡沫铝通过塑性变形消耗了更多冲击能量。当孔隙率P从0.4提升到0.7时，泡沫铝对冲击物的动能吸收率从90%提高到99%。图8是冲击物冲击不同孔隙率泡沫铝夹层结构时的加速度a变化曲线，可见孔隙率越低，加速度变化曲线越陡峭、峰值越高，冲击碰撞过程越偏向刚性。通过以上分析，可以进一步明确，泡沫铝芯材孔隙率对夹层结构的抗冲击性能有显著影响，在一定范围内，孔隙率越高，缓冲效果越好，抗冲击性能越强。需要注意的是，上述规律在孔隙率进一步提高时将不再适用，虽然提高孔隙率能增强缓冲能力，但过大的孔隙率会造成泡沫铝芯材承载能力显著下降，无法充分发挥其防护作用。
表 2 冲击过程中冲击物的主要运动参数Table 2. Main motion parameters of impact object in the process of impact孔隙率 反弹速度/(m∙s−1) 残余动能/mJ 加速度峰值/(km∙s−2) 动能吸收率/% 0.7 0.68 8 47 99 0.6 2.91 15 53 97 0.5 4.72 39 66 94 0.4 6.08 65 85 90 采用图9所示的落锤实验系统获取冲击物加速度峰值，并与模拟计算结果对比，图10显示模拟计算值稍高于实验值，但是与孔隙率之间关系的变化规律是一致的。出现这一情况的主要原因是计算中将问题简化为二维模型导致的差异，此外，模拟计算中未考虑泡沫铝和铝蒙皮面板的微观裂纹扩展情况也是误差的重要来源。尽管存在这一偏差，模拟计算值与实验结果总体量级上吻合良好，并且表现出相同的规律，证明了近场动力学模拟分析结果与结论的有效性。
图 9 采用落锤实验获取冲击物加速度峰值Figure 9. The peak impact acceleration of the impactor obtained by a drop-weight experiment3.3 有机玻璃背板的裂纹范围
某光学舱采用泡沫铝作为防护结构是为了保护底层光学材料不被破坏，考察泡沫铝夹层结构抗冲击性能最直观的方式便是观察底部有机玻璃背板的破坏形态。近场动力学中采用临界伸长率预测脆性材料断裂问题的方法已被证明是有效的，因此，基于此方法计算出有机玻璃背板的裂纹起始与终止形态如图11所示。
如图11所示，冲击过程中，有机玻璃背板最先发生破坏的位置是中间正对冲击的部位，并逐渐向两侧扩展。随着泡沫铝孔隙率的降低，有机玻璃背板的裂纹扩展范围明显增大，说明高孔隙率的泡沫铝比低孔隙率的泡沫铝有更好的防护效果，进一步验证高孔隙率的泡沫铝夹层结构具有更好的抗冲击性能。此外，上下粘接界面也出现了一定损伤，结合3.1节可知，这是由于在冲击过程中上层铝面板首先发生塑性变形，随着铝面板下凹，铝面板和泡沫铝之间的粘接界面势必会因变形不协调而产生横向剪切力，这是导致上层粘接界面开裂的主要原因，随着泡沫铝芯材变形增加，引起结构发生整体弯曲变形，进而使泡沫铝和有机玻璃背板之间的粘接界面产生剪切和张拉应力，导致粘接界面出现损伤。从图11可以看出不同孔隙率的泡沫铝粘接界面损伤程度基本一致，说明孔隙率对粘接界面开裂程度的影响较小。
4. 结 论
在常规态型近场动力学理论中引入了Mises屈服准则和线性各向同性强化模型，建立了近场动力学塑性本构的数值计算框架。采用随机方法创建了泡沫铝结构的二维细观模型，并基于近场动力学计算程序模拟了冲击作用下不同孔隙率泡沫铝的塑性变形以及有机玻璃背板的裂纹扩展形态，验证了近场动力学解决泡沫铝夹层结构抗冲击问题的可行性，为相关领域问题提供了新的思路。探讨了在20 m/s的低速冲击作用下，孔隙率对泡沫铝夹层结构抗冲击性能的影响规律，得到以下结论。
(1)泡沫铝夹层结构良好的塑性变形能力是其发挥缓冲与防护作用的主要因素。一方面，冲击作用下泡沫铝结构先表现出局部变形，随着泡沫铝内部孔洞被压实，结构逐渐表现出整体弯曲变形，在这一变形过程中，泡沫铝内部孔洞能大幅削减应力从冲击部位向周围的扩散范围，起到缓冲作用；另一方面，泡沫铝芯材发生大范围塑性变形，冲击物的动能在这一过程中被耗散，使作用在背板上的冲击力被缓解。
(2)在一定范围内，泡沫铝芯材孔隙率越高，则夹层结构具有越好的抗冲击性能。孔隙率从0.4升高到0.7时，泡沫铝对冲击物的动能吸收率从90%提高到99%，冲击物加速度峰值从85 km/s2降低到47 km/s2。模拟结果与实验结果具有较好的一致性，验证了模拟结果的准确性和分析结论的有效性。
(3)近场动力学方法能有效预测有机玻璃背板的裂纹扩展以及粘接界面开裂过程，模拟结果表明，泡沫铝芯材孔隙率越高，有机玻璃背板的裂纹扩展范围越小，进一步说明本文所研究的某光学舱防护结构应尽量选择更高孔隙率的泡沫铝。
-
图 2 破片冲击泡沫铝夹层结构模型
Figure 2. Schematic diagram of an aluminum foam sandwich impacted by a fragment
图 3 Monte-Carlo方法生成泡沫铝算法流程
Figure 3. Algorithm flow chart of the Monte-Carlo method to generate aluminum foam
图 4 泡沫铝夹层结构的近场动力学离散模型
Figure 4. Discrete models of aluminum foam sandwiches for peridynamics
图 9 采用落锤实验获取冲击物加速度峰值
Figure 9. The peak impact acceleration of the impactor obtained by a drop-weight experiment
表 1 材料参数
Table 1. Material parameters
材料 ρ/(kg∙m−3) E/GPa ν σy/MPa Et/MPa sc 铝 2700 69.5 0.33 127 586 有机玻璃 1190 3.6 0.4 0.047 环氧树脂 3 0.37 0.047 
下载: 导出CSV
表 2 冲击过程中冲击物的主要运动参数
Table 2. Main motion parameters of impact object in the process of impact
孔隙率 反弹速度/(m∙s−1) 残余动能/mJ 加速度峰值/(km∙s−2) 动能吸收率/% 0.7 0.68 8 47 99 0.6 2.91 15 53 97 0.5 4.72 39 66 94 0.4 6.08 65 85 90
下载: 导出CSV
-
[1] ZHANG J X, QIN Q H, WANG T J. Compressive strengths and dynamic response of corrugated metal sandwich plates with unfilled and foam-filled sinusoidal plate cores [J]. Acta Mechanica, 2013, 224(4): 759–775. DOI: 10.1007/s00707-012-0770-5. [2] ZHANG J X, ZHOU R F, WANG M S, et al. Dynamic response of double-layer rectangular sandwich plates with metal foam cores subjected to blast loading [J]. International Journal of Impact Engineering, 2018, 122: 265–275. DOI: 10.1016/j.ijimpeng.2018.08.016. [3] ZHANG J X, YE Y, QIN Q H. On dynamic response of rectangular sandwich plates with fibre-metal laminate face-sheets under blast loading [J]. Thin-Walled Structures, 2019, 144: 106288. DOI: 10.1016/j.tws.2019.106288. [4] 郭亚周, 刘小川, 何思渊, 等. 不同弹形撞击下泡沫铝夹芯结构动力学性能研究 [J]. 兵工学报, 2019, 40(10): 2032–2041. DOI: 10.3969/j.issn.1000-1093.2019.10.008.GUO Y Z, LIU X C, HE S Y, et al. Research on dynamic properties of aluminum foam sandwich structure impacted by projectiles with different shapes [J]. Acta Armamentarii, 2019, 40(10): 2032–2041. DOI: 10.3969/j.issn.1000-1093.2019.10.008. [5] 张博一, 赵威, 王理, 等. 泡沫铝子弹高速撞击下铝基复合泡沫夹层板的动态响应 [J]. 爆炸与冲击, 2017, 37(4): 600–610. DOI: 10.11883/1001-1455(2017)04-0600-11.ZHANG B Y, ZHAO W, WANG L, et al. Dynamic response of aluminum matrix syntactic foams sandwich panel subjected to foamed aluminum projectile impact loading [J]. Explosion and Shock Waves, 2017, 37(4): 600–610. DOI: 10.11883/1001-1455(2017)04-0600-11. [6] 邓旭辉, 李亚斌. 双层泡沫铝夹芯板抗冲击性能数值研究 [J]. 铁道科学与工程学报, 2019, 16(10): 2603–2611. DOI: 10.19713/j.cnki.43-1423/u.2019.10.029.DENG X H, LI Y B. Numerical study on impact resistance of double-layer foam aluminum sandwich panels [J]. Journal of Railway Science and Engineering, 2019, 16(10): 2603–2611. DOI: 10.19713/j.cnki.43-1423/u.2019.10.029. [7] 张永康, 李玉龙, 汤忠斌, 等. 冰雹撞击下泡沫铝夹芯板的动态响应 [J]. 爆炸与冲击, 2018, 38(2): 373–380. DOI: 10.11883/bzycj-2016-0232.ZHANG Y K, LI Y L, TANG Z B, et al. Dynamic response of aluminum-foam-based sandwich panels under hailstone impact [J]. Explosion and Shock Waves, 2018, 38(2): 373–380. DOI: 10.11883/bzycj-2016-0232. [8] 夏志成, 张建亮, 周竞洋, 等. 泡沫铝夹芯板抗冲击性能分析 [J]. 工程力学, 2017, 34(10): 207–216. DOI: 10.6052/j.issn.1000-4750.2016.06.0494.XIA Z C, ZHANG J L, ZHOU J Y, et al. Analysis on impact resistance of aluminum foam sandwich panels [J]. Engineering Mechanics, 2017, 34(10): 207–216. DOI: 10.6052/j.issn.1000-4750.2016.06.0494. [9] CUI T N, ZHANG J H, LI K K, et al. Ballistic limit of sandwich plates with a metal foam core [J]. Journal of Applied Mechanics, 2022, 89(2): 021006. DOI: 10.1115/1.4052835. [10] 苏兴亚, 敬霖, 赵隆茂. 爆炸载荷下分层梯度泡沫铝夹芯板的失效模式与抗冲击性能 [J]. 爆炸与冲击, 2019, 39(6): 063103. DOI: 10.11883/bzycj-2018-0198.SU X Y, JING L, ZHAO L M. Failure modes and shock resistance of sandwich panels with layered-gradient aluminum foam cores under air-blast loading [J]. Explosion and Shock Waves, 2019, 39(6): 063103. DOI: 10.11883/bzycj-2018-0198. [11] 王涛, 余文力, 秦庆华, 等. 爆炸载荷下泡沫铝夹芯板变形与破坏模式的实验研究 [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. [12] SILLING S A. Reformulation of elasticity theory for discontinuities and long-range forces [J]. Journal of the Mechanics and Physics of Solids, 2000, 48(1): 175–209. DOI: 10.1016/S0022-5096(99)00029-0. [13] SILLING S A, ASKARI E. A meshfree method based on the peridynamic model of solid mechanics [J]. Computers & Structures, 2005, 83(17/18): 1526–1535. DOI: 10.1016/j.compstruc.2004.11.026. [14] SILLING S A, EPTON M, WECKNER O, et al. Peridynamic states and constitutive modeling [J]. Journal of Elasticity, 2007, 88(2): 151–184. DOI: 10.1007/s10659-007-9125-1. [15] 杨娜娜, 赵天佑, 陈志鹏, 等. 破片冲击作用下舰船复合材料结构损伤的近场动力学模拟 [J]. 爆炸与冲击, 2020, 40(2): 023302. DOI: 10.11883/bzycj-2019-0019.YANG N N, ZHAO T Y, CHEN Z P, et al. Peridynamic simulation of damage of ship composite structure under fragments impact [J]. Explosion and Shock Waves, 2020, 40(2): 023302. DOI: 10.11883/bzycj-2019-0019. [16] 刘宁, 胡梦凡, 周飞. 基于键基近场动力学理论的单裂纹圆孔板冲击破坏研究 [J]. 工程力学, 2020, 37(12): 9–17. DOI: 10.6052/j.issn.1000-4750.2020.02.0076.LIU N, HU M F, ZHOU F. The impacted damage study of a single cleavage drilled compression specimen based on bond-based peridynamics [J]. Engineering Mechanics, 2020, 37(12): 9–17. DOI: 10.6052/j.issn.1000-4750.2020.02.0076. [17] 熊伟鹏, 王超, 傅江妍, 等. 冰球冲击试验的近场动力学方法数值模拟 [J]. 振动与冲击, 2020, 39(7): 148–155. DOI: 10.13465/j.cnki.jvs.2020.07.021.XIONG W P, WANG C, FU J Y, et al. Numerical simulation of ice sphere impact test by peridynamics method [J]. Journal of Vibration and Shock, 2020, 39(7): 148–155. DOI: 10.13465/j.cnki.jvs.2020.07.021. [18] MADENCI E, OTERKUS S. Ordinary state-based peridynamics for plastic deformation according to von Mises yield criteria with isotropic hardening [J]. Journal of the Mechanics and Physics of Solids, 2016, 86: 192–219. DOI: 10.1016/j.jmps.2015.09.016. 期刊类型引用(1)
1. 陈洋，王肇喜，翟师慧，盛鹏，王者蓝，朱明亮. 3D打印点阵夹芯结构冲击损伤的近场动力学模拟. 爆炸与冲击. 2024(03): 146-160 . 
本站查看其他类型引用(6)
-
下载:
计量
- 文章访问数: 719
- HTML全文浏览量: 193
- PDF下载量: 255
- 被引次数: 7