Dynamic response of a sandwich panel cored by butterfly-shaped honeycombs with negative Poisson’s ratio to low-velocity impact
-
摘要: 为了研究负泊松比蝴蝶形蜂窝夹芯板在低速冲击下的动力学响应,采用质量-弹簧模型获得了冲击器与蜂窝夹芯板之间的接触力,同时基于哈密顿原理和一阶剪切变形理论推导了负泊松比蝴蝶形蜂窝夹芯板的运动方程,采用Navier法和Duhamel积分对蜂窝板的振动位移进行了理论解析求解。在理论验证方面,蜂窝夹芯板前5阶固有频率的数值模拟结果与理论模型计算结果的最大相对误差为6.52%,蜂窝夹芯板中心最大横向位移的数值模拟结果与理论模型计算结果的最大相对误差为6.84%,理论模型求解的接触力与文献得到的接触力的最大相对误差为8%,验证了理论模型的有效性。结果表明,随着球形冲击器冲击速度的递增,蜂窝夹芯板的最大横向位移呈现递增的规律。而在相同冲击载荷下,蜂窝夹芯板的抗冲击特性随着胞元壁厚的增大而增强,随着胞元角度的增大而减弱;随着负泊松比蝴蝶形蜂窝夹芯板长宽比以及夹芯层与顶部蒙皮层的高度比的增大,蜂窝夹芯板的横向位移减小,冲击器与蜂窝夹芯板之间的接触力增大。当蜂窝夹芯板的宽长比从1∶1变化到1∶2时,蜂窝夹芯板最大横向位移减小6.1%;当顶部蒙皮层与蜂窝芯层的高度比从1∶6变化到1∶14时,蜂窝夹芯板的最大横向位移减小5.4%,这表明蜂窝夹芯板的抗冲击性能增强,吸能效果明显。Abstract: In order to study the dynamic response of a sandwich panel cored by butterfly-shaped honeycomb with negative Poisson’s ratio to low-velocity impact, a mass-spring (MS) model is applied to obtain the contact force between the spherical impactor and the honeycomb sandwich panel. Meanwhile, based on the Hamilton’s principle and the first-order shear deformation theory, the equation of motion for the butterfly-shaped honeycomb sandwich panel with negative Poisson’s ratio is derived. Besides, the Navier method and Duhamel’s integral are used to solve the vibration displacement of the honeycomb sandwich panel. To validate the theoretical model, the results are compared with the results of ABAQUS’ numerical simulation or published literature. It is shown that the maximum relative error between the numerical modeling results of the first five order natural frequencies and the results of theoretical model calculated in this paper is 6.52%, the maximum relative error between the numerical modeling results of the honeycomb sandwich panel under low-velocity impact and the calculated results of the theoretical model in this paper is 6.84%, and the maximum relative error of the contact force between the theoretical model in this paper and the published studies is 8%, thus verifying the validity of the theoretical model. The results show that the maximum lateral displacement of the honeycomb sandwich panel increases with the increasing velocity of the spherical impactor. Under the same impact load, the impact resistance of the honeycomb sandwich panel increases with the increase of the wall thickness of the unit cell, and decreases with the increase of the unit cell angle. The impact resistance of the honeycomb sandwich panel increases by 3.7% when the thickness of the unit cell wall changes from 1 mm to 3 mm. The lateral displacement of the butterfly-shaped honeycomb sandwich panel decreases while the contact force between the impactor and the honeycomb sandwich panel increases with the increase of the length-width ratio and the height ratio. When the width-length ratio of the honeycomb sandwich panel changes from 1∶1 to 1∶2, the maximum lateral displacement of the honeycomb sandwich panel decreases by 6.1%, and when the height ratio of the top skin layer to the honeycomb core layer changes from 1∶6 to 1∶14, the maximum lateral displacement of the honeycomb sandwich panel decreases by 5.4%, which indicates that the impact resistance of the honeycomb sandwich panel is enhanced and the energy absorption effect is obvious.
-
表 1 网格尺寸对一阶固有频率计算结果的影响
Table 1. Influences of mesh size on the calculated results of the first-order natural frequency
网格尺寸/mm 频率/Hz 15.0 719.26 12.5 699.18 10.0 696.14 7.5 692.37 5.0 691.25 表 2 蝴蝶形蜂窝夹芯板的前5阶固有频率
Table 2. The first five order natural frequencies of the butterfly-shaped honeycomb sandwich panel
阶数 频率/Hz 理论模型解 有限元解 1 646.17 691.25 2 1597.40 163.75 3 1648.40 1706.70 4 2570.50 2629.30 5 3181.90 3200.80 表 3 低速冲击下蜂窝夹芯板中心网格尺寸对中心最大横向位移计算结果的影响
Table 3. Influence of center grid size on the calculated maximum lateral displacement of the sandwich panel center under low-velocity impact
网格尺寸/mm 中心最大横向位移/mm 10.0 0.1928 5.0 0.2129 2.5 0.2218 1.0 0.2516 0.5 0.2556 表 4 不同冲击速度下蝴蝶形蜂窝夹芯板的横向位移
Table 4. Lateral displacement of the butterfly-shaped honeycomb sandwich panel at different impact velocities
冲击速度/(m·s−1) 中心最大横向位移/mm 本文模型 有限元模拟 6 0.177 0.190 8 0.247 0.252 10 0.320 0.311 12 0.395 0.374 表 5 矩形板和冲击器的材料参数和几何尺寸[26-27]
Table 5. Material parameters and geometrical sizes of the homogeneous panel and impactor[26-27]
器件 密度/(kg·m−3) 弹性模量/GPa 泊松比 半径/mm 长度/mm 宽度/mm 厚度/mm 冲击器 7971.8 200 0.3 10 矩形板 7971.8 200 0.3 200 200 8 -
[1] PRAWOTO Y. Seeing auxetic materials from the mechanics point of view: a structural review on the negative Poisson’s ratio [J]. Computational Materials Science, 2012, 58: 140–153. DOI: 10.1016/j.commatsci.2012.02.012. [2] BIRMAN V, KARDOMATEAS G A. Review of current trends in research and applications of sandwich structures [J]. Composites Part B: Engineering, 2018, 142: 221–240. DOI: 10.1016/j.compositesb.2018.01.027. [3] VIJAYASIMHA REDDY B G, SHARMA K V, YELLA REDDY T. Deformation and impact energy absorption of cellular sandwich panels [J]. Materials and Design, 2014, 61: 217–227. DOI: 10.1016/j.matdes.2014.04.047. [4] 王守财, 关志东, 黎增山, 等. 金属蒙皮飞机冲击威胁研究 [J]. 振动与冲击, 2018, 37(4): 13–18. DOI: 10.13465/j.cnki.jvs.2018.4.003.WANG S C, GUAN Z D, LI Z S, et al. A study on the impact threat of metal skin aircrafts [J]. Journal of Vibration and Shock, 2018, 37(4): 13–18. DOI: 10.13465/j.cnki.jvs.2018.4.003. [5] ZHANG J X, YUAN H, LI J F, et al. Dynamic response of multilayer curved aluminum honeycomb sandwich beams under low-velocity impact [J]. Thin-Walled Structures, 2022, 177: 109446. DOI: 10.1016/j.tws.2022.109446. [6] 蔺晓红, 张涛, 张小波, 等. 碳纤维增强铝合金板的抗冲击性能 [J]. 爆炸与冲击, 2013, 33(3): 303–310. DOI: 10.11883/1001-1455(2013)03-0303-08.LIN X H, ZHANG T, ZHANG X B, et al. Impact resistances of carbon fiber-reinforced aluminum laminates [J]. Explosion and Shock Waves, 2013, 33(3): 303–310. DOI: 10.11883/1001-1455(2013)03-0303-08. [7] 谢素超, 井坤坤, 冯哲骏, 等. 铝蜂窝夹层板低速冲击响应及损伤模式的参数化影响 [J]. 复合材料学报, 2022, 40(5): 3046–3060. DOI: 10.13801/j.cnki.fhclxb.20220706.002.XIE S C, JING K K, FENG Z J, et al. Parametric effects of low-velocity impact response and damage mode of aluminum honeycomb sandwich panels [J]. Acta Materiae Compositae Sinica, 2022, 40(5): 3046–3060. DOI: 10.13801/j.cnki.fhclxb.20220706.002. [8] ZHANG D H, JIANG D, FEI Q G, et al. Experimental and numerical investigation on indentation and energy absorption of a honeycomb sandwich panel under low-velocity impact [J]. Finite Elements in Analysis and Design, 2016, 117/118: 21–30. DOI: 10.1016/j.finel.2016.04.003. [9] LIU P F, LI X K, LI Z B. Finite element analysis of dynamic mechanical responses of aluminum honeycomb sandwich structures under low-velocity impact [J]. Journal of Failure Analysis and Prevention, 2017, 17(6): 1202–1207. DOI: 10.1007/s11668-017-0358-4. [10] PALOMBA G, EPASTO G, CRUPI V, et al. Single and double-layer honeycomb sandwich panels under impact loading [J]. International Journal of Impact Engineering, 2018, 121: 77–90. DOI: 10.1016/j.ijimpeng.2018.07.013. [11] ZHANG J X, ZHU Y Q, LI K K, et al. Dynamic response of sandwich plates with GLARE face-sheets and honeycomb core under metal foam projectile impact: experimental and numerical investigations [J]. International Journal of Impact Engineering, 2022, 164: 104201. DOI: 10.1016/j.ijimpeng.2022.104201. [12] HAZIZAN M A, CANTWELL W J. The low velocity impact response of an aluminium honeycomb sandwich structure [J]. Composites Part B: Engineering, 2003, 34(8): 679–687. DOI: 10.1016/S1359-8368(03)00089-1. [13] ZHANG J X, YE Y, QIN Q H, et al. Low-velocity impact of sandwich beams with fibre-metal laminate face-sheets [J]. Composites Science and Technology, 2018, 168: 152–159. DOI: 10.1016/j.compscitech.2018.09.018. [14] ZHANG J X, QIN Q H, ZHANG J T, et al. Low-velocity impact on square sandwich plates with fibre-metal laminate face-sheets: analytical and numerical research [J]. Composite Structures, 2021, 259: 113461. DOI: 10.1016/j.compstruct.2020.113461. [15] SUN M Q, WOWK D, MECHEFSKE C, et al. An analytical study of the plasticity of sandwich honeycomb panels subjected to low-velocity impact [J]. Composites Part B: Engineering, 2019, 168: 121–128. DOI: 10.1016/j.compositesb.2018.12.071. [16] WU Z, MA R, LI Y Z, et al. An advanced five-unknown higher-order theory for free vibration of composite and sandwich plates [J]. Chinese Journal of Aeronautics, 2021, 34(9): 104–118. DOI: 10.1016/J.CJA.2021.03.002. [17] 于靖军, 谢岩, 裴旭. 负泊松比超材料研究进展 [J]. 机械工程学报, 2018, 54(13): 1–14. DOI: 10.3901/JME.2018.13.001.YU J J, XIE Y, PEI X. State-of-art of metamaterials with negative Poisson’s ratio [J]. Journal of Mechanical Engineering, 2018, 54(13): 1–14. DOI: 10.3901/JME.2018.13.001. [18] REN X, DAS R, TRAN P, et al. Auxetic metamaterials and structures: a review [J]. Smart Materials and Structures, 2018, 27(2): 023001. DOI: 10.1088/1361-665X/aaa61c. [19] 关淮桐, 田瑞兰, 张子文. 仿蝴蝶形蜂窝结构夹层板的振动特性研究 [J]. 哈尔滨工程大学学报, 2022, 43(9): 1383–1390. DOI: 10.11990/jheu.202206045.GUAN H T, TIAN R L, ZHANG Z W. Study of the vibrational characteristics of butterfly-shaped honeycomb sandwich panel [J]. Journal of Harbin Engineering University, 2022, 43(9): 1383–1390. DOI: 10.11990/jheu.202206045. [20] THAI H T, CHOI D H. A simple first-order shear deformation theory for laminated composite plates [J]. Composite Structures, 2013, 106: 754–763. DOI: 10.1016/j.compstruct.2013.06.013. [21] ZHANG Z W, TIAN R L, ZHANG X L, et al. A novel butterfly-shaped auxetic structure with negative Poisson’s ratio and enhanced stiffness [J]. Journal of Materials Science, 2021, 56(25): 14139–14156. DOI: 10.1007/S10853-021-06141-4. [22] ZHAO H X, MA S J. High-order Hamilton’s principle and the Hamilton’s principle of high-order Lagrangian function [J]. Communications in Theoretical Physics, 2008, 49(2): 297–302. DOI: 10.1088/0253-6102/49/2/08. [23] AMROUCHE C, REJAIBA A. L-p-theory for Stokes and Navier-Stokes equations with Navier boundary condition [J]. Journal of Differential Equations, 2014, 256(4): 1515–1547. DOI: 10.1016/j.jde.2013.11.005. [24] ROZANSKI M, SIKORA B, SMUDA A, et al. On theoretical and practical aspects of Duhamel’s integral [J]. Archives of Control Sciences, 2021, 31(4): 815–847. DOI: 10.24425/acs.2021.139732. [25] TAN T M, SUN C T. Use of statical indentation laws in the impact analysis of laminated composite plates [J]. Journal of Applied Mechanics, 1985, 52(1): 6–12. DOI: 10.1115/1.3169029. [26] YANG F L, WANG Y Q, LIU Y F. Low-velocity impact response of axially moving functionally graded graphene platelet reinforced metal foam plates [J]. Aerospace Science and Technology, 2022, 123: 107496. DOI: 10.1016/j.ast.2022.107496. [27] WU H Y T, CHANG F K. Transient dynamic analysis of laminated composite plates subjected to transverse impact [J]. Computers and Structures, 1989, 31(3): 453–466. DOI: 10.1016/0045-7949(89)90393-3.