Next Article in Journal
Machine Learning Method to Explore the Correlation between Fly Ash Content and Chloride Resistance
Previous Article in Journal
Electro-Optical Comb Envelope Engineering Based on Mode Crossing
Previous Article in Special Issue
The Effect of Microballoon Volume Fraction on the Elastic and Viscoelastic Properties of Hollow Microballoon-Filled Epoxy Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanical Behavior of Bio-Inspired Honeycomb–Core Composite Sandwich Structures to Low-Velocity Dynamic Loading

1
Guangzhou Metro Design & Research Institute Co., Ltd., Guangzhou 510010, China
2
School of Mechanical Engineering, Jiangsu University of Science and Technology, Zhenjiang 212000, China
*
Author to whom correspondence should be addressed.
Materials 2024, 17(5), 1191; https://doi.org/10.3390/ma17051191
Submission received: 5 January 2024 / Revised: 29 January 2024 / Accepted: 7 February 2024 / Published: 4 March 2024
(This article belongs to the Special Issue Lightweight and High-Strength Sandwich Panel)

Abstract

:
In order to improve the impact resistance of sandwich panels under low-velocity impact, the lotus leaf vein is selected as a biological prototype to design a bio-inspired honeycomb (BIH) sandwich panel. ABAQUS is used to establish and effectively verify the finite element (FE) model of the BIH sandwich panel. To systematically compare and study the mechanical properties of BIH and conventional hexagonal honeycomb sandwich panels under low-velocity impact, the maximum displacement of face-sheets, the deformation mode, the plastic energy consumption and the dynamic response curve of the impact end are presented. At the same time, the performance differences between them are revealed from the perspective of an energy absorption mechanism. Furthermore, the influence of the circumscribed circle diameter ratio of the BIH trunk to branch (γ), the thickness ratio of the trunk to branch (K) and the impact angle (θ) on impact resistance is studied. Finally, the BIH sandwich panel is further optimized by using the response surface method. It can be concluded that, compared to conventional hexagonal honeycomb sandwich panels, the addition of walls in the BIH sandwich panel reduces the maximum deformation of the rear face-sheet by 10.29% and increases plastic energy consumption by 8.02%. Properly adjusting the structural parameters can effectively enhance the impact resistance of the BIH sandwich panel.

1. Introduction

As a typical lightweight thin-walled structure, honeycomb sandwich panels have been widely used as protective equipment due to their outstanding impact resistance and meeting the requirements of lightweight design [1,2,3]. They also have a wide range of applications in the fields of aviation [4], automobiles [5] and personal protective equipment [6]. Therefore, the honeycomb sandwich structure has always been a hot topic for scholars [7,8,9,10,11,12]. Rejab [13] et al. focused on studying the mechanical properties of corrugated sandwich panels, which exhibit stronger impact resistance. Sun [14] et al. introduced periodic grids to enhance the hexagonal honeycomb core sandwich structure and conducted plane compression experiments and numerical studies. Research has shown that the mechanical properties of sandwich panels were closely related to their core structure. Lan [15] et al. comprehensively analyzed the dynamic response of hexagonal honeycomb cores and cylindrical sandwich panels with auxetic concave honeycomb cores. The cylindrical sandwich panel with auxetic concave honeycomb cores was proven to have a better anti-explosive performance than those with hexagonal honeycomb cores. McShane et al. [16] found that the three special-shaped sandwich panels have better impact resistance.
Meanwhile, the honeycomb sandwich structure is very vulnerable to low-speed impacts from gravel, hail, etc. Therefore, scholars have paid great attention to the mechanical properties of honeycomb sandwich panels under low-speed impact. Sun et al. [17] analyzed the low-velocity impact characteristics of sandwich panels with different structural parameters and carried out relevant experimental research. Zhang et al. [18,19] combined experimental and numerical methods to derive an analytical model of energy balance based on the dynamic response of a spherical indenter impacting the sandwich panel at low speed. The research of Liu et al. [20] showed that sandwich panels with round tube-filled hexagonal honeycomb cores had better impact resistance than conventional hexagonal honeycomb core sandwich panels. Reyes and Børvik [21] conducted quasi-static and dynamic experiments on the impact resistance of polymer foam core sandwich panels. Zhou [22] proposed an auxetic recessed honeycomb filled with new foam concrete and found that, with an increase in foam concrete density, the structure changed from a compression failure mode to a shear failure mode. Dhanarasu et al. [23] systematically analyzed the mechanical response of sandwich panels with CFRP face-sheets. Dai and Akatay [24,25] focused on the effects of repeated impacts on sandwich panels. Zhang and Al [26,27] proposed tube-reinforced honeycomb and comprehensively analyzed its mechanical properties.
Numerous studies have shown that the performance of sandwich plates is closely related to its cell structure [28,29,30,31,32,33]. With the improvement in energy absorption capacity requirements, how to design a honeycomb sandwich structure with better impact resistance has been a challenge for scholars. Researchers have found that applying biological structural features to the design of thin-walled structures can effectively improve their crashworthiness [34,35,36,37,38]. Song et al. [39] designed a new type of bionic tube with grooves and studied its crashworthiness under lateral impact. They found that the energy absorption efficiency of the bionic tube was better than that of the circular tube. Based on the structural characteristics of bamboo, Zou et al. [40] proposed a bionic tube and solved its numerical examples under axial/transverse impact. The results showed that the bionic design improved the specific energy absorption of the tube. Palombini et al. [41] mechanically explored the special geometry of a single vascular bundle in bamboo, and the new design had a noticeable improvement in its strength and energy absorption under axial/transverse loads. It is widely known that the leaf veins not only carry nutrients, but also support the weight and external load of the whole blade structure to ensure the overall stiffness and plane ductility of the blade. He et al. [42] introduced this leaf-vein-branched characteristic into the honeycomb in previous studies, but they only focused on the dynamic response under in-plane dynamic compression.
Inspired by the lotus leaf vein biological structure, we designed a bio-inspired honeycomb (BIH) core sandwich structure. Based on ABAQUS 6.14 software, the dynamic response process of the BIH sandwich panel under low-speed dynamic impact was systematically investigated. Firstly, the reasonableness and accuracy of the simulation method used in this paper were verified by comparing them with previous experimental results. Secondly, the dynamic responses of the BIH sandwich panel and conventional hexagonal honeycomb sandwich panels under low-velocity impact were systematically compared, and the performance difference between them was analyzed from the perspective of the energy absorption mechanism. Then, the coupling influence law of the structural parameters of the BIH sandwich panel and its mechanical properties was further analyzed by parametric study. Finally, the BIH sandwich structure was optimized using the response surface method.

2. Geometric Model of the BIH Sandwich Structure

Inspired by the biological structure of nature, the lotus leaf vein was chosen as a biological prototype. Figure 1 shows the design process and geometric configuration of the bio-inspired honeycomb core sandwich structure. The geometric structure of the lotus leaf vein in Figure 1a is characterized by the main trunk branch, and the cross-sectional size of the trunks is larger than that of the branches. Figure 1b shows the BIH core sandwich structure and its cell geometry is shown in Figure 1c. The dotted circles, as seen in the picture, are the circumscribed circles of the main trunks and branches. tt and tb are the wall thicknesses of the trunks and branches, respectively, and Φb and Φt are the circumscribed circle diameter of the branches and trunks, respectively. The BIH structure parameter, K, is the ratio of tt to tb, and γ is the ratio of Φt to Φb, as shown in Equations (1) and (2). The branch circumscribed circle diameter Φb is 12 mm, and the sum of tt and tb is 0.15 mm, as shown in Equations (3) and (4). The range of values for γ and K are 0.1 to 0.7 and 0 to 10, respectively. The structure size and wall thickness of the BIH change with parameters K and γ, respectively. Figure 1d depicts the cross-sectional evolution of BIH cells, and Table 1 gives the specific parameters of the BIH sandwich panel.
K = t t t b   0 < K < 10
γ = Φ t Φ b   0.1 γ 0.7
t t + t b = 0.15   m m
Φ b = 12   m m

3. Finite Element Method of BIH Sandwich Structure

3.1. FE Model

ABAQUS is used to simulate the mechanical properties of the BIH sandwich panel under low-speed impact, and the numerical simulation model is depicted in Figure 2. The material model of the BIH sandwich structure adopts a bilinear elastoplastic model, the panel material is AL5083 and the honeycomb core material is AL3003 [18]. The specific material parameters of AL5083 and AL3003 are shown in Table 2. A binding constraint is adopted between the honeycomb core and the panel; the trunk and branches are connected together by means of co-nodes. The diameter and mass of the spherical impactor are 40 mm and 0.2 kg, respectively. The mass and velocity are set at the center of the sphere and the impact direction is downward and perpendicular to the sandwich panel. The impact point is at the center of the sandwich panel. When the spherical impactor impacts the BIH interlayer structure, the speed of the spherical impactor is 35 m/s [20]. The spherical impact block is set as a rigid body and meshed by 4-node 3D bilinear rigid elements. The 4-node doubly curved shell element is used to mesh the face-sheets and the honeycomb core. Meanwhile, the contact between the spherical impactor and the sandwich panel is simulated by a surface-to-surface contact. A general contact is set for the rest of the contacts to avoid interpenetration of the honeycomb cores. The static and dynamic friction coefficients between each structure are both set to 0.3 [20]. The three sides of the face-sheets are fixed [18], as reflected in the red line in Figure 2.
The grid size of the finite element model has an inseparable influence on the calculation cost and accuracy. After weighing the calculation cost and accuracy, a gradient grid is used to divide the finite element model into several areas. According to Figure 3, the light blue area is meshed by a grid with a size of 0.5 mm × 0.5 mm, the grid size of the orange area is 1.0 mm × 1.0 mm and the other areas are 2.0 mm × 2.0 mm.

3.2. Validity Verification

The experimental results of Zhang [18] were used to verify the validity of the FE model. The detailed information of drop-weight tests was reported in [18]. In order to be consistent with the test conditions, the weight of the impactor was 0.3 kg with a radius of 20 mm and an impact velocity of 5.05 m/s. Other relevant parameters and settings, such as the structural dimensions, the material properties and the boundary conditions, are consistent with those in Section 3.1. The numerical simulation results were compared with the experimental results, as shown in Figure 4. It can be seen from the map that the experimental and numerical deformation modes are in good agreement. The errors of the dent depth on the rear face-sheet and the energy absorption are 3.24% and 4.67%, respectively. Therefore, the FE method and the material model are accurate enough to be used to simulate the dynamic response of the BIH sandwich panels.

4. Comparative Analysis of Two Honeycomb Core Sandwich Structures

4.1. Dynamic Response Analysis

In order to compare the performance difference between the BIH and the conventional hexagonal honeycomb sandwich panel, the cell wall thickness of the honeycomb core is adjusted to ensure the same mass. The sizes of both honeycomb cores are 150   m m × 150   m m × 16   m m . The other structural parameters are shown in Table 3. Figure 5 shows the dynamic response curve at the impact end. The dynamic response of the impact end of these two structures shows a trend of first increasing and then decreasing. When the spherical block impacts the two sandwich structures, it goes through three stages, namely, an elastic stage (I), a plastic deformation stage (II) and a rebound stage (III). In the elastic phase, the BIH sandwich structure has the same impact level as the conventional hexagonal honeycomb sandwich panel, which shows that the initial structural rigidity of the BIH sandwich panel is similar to that of the conventional structure. In the plastic deformation stage, the impact load level of the BIH sandwich structure is higher than that of the conventional hexagonal honeycomb sandwich structure.
Figure 6 exhibits the maximum displacement of the face-sheet of the BIH sandwich structure and the conventional hexagonal honeycomb sandwich structure. The maximum displacement of the front face-sheet of the BIH sandwich structure is 16.35 mm, which is 3.22% higher than the 15.84 mm conventional sandwich structure displacement. However, the maximum displacement of the rear panel decreased by 10.29% compared to the conventional sandwich structure. This means that the BIH sandwich cores are easier to compress and deform than the conventional hexagonal honeycomb sandwich cores, which is beneficial for the absorption of more impact kinetic energy. Therefore, the BIH sandwich structure shows better impact resistance.

4.2. Energy Absorption and Its Enhancement Mechanism

Figure 7 exhibits the plastic energy consumption of each part of the BIH sandwich structure and the conventional hexagonal honeycomb sandwich structure. The plastic energy consumptions of the front face-sheet and the honeycomb core of the BIH sandwich panel are higher than those of the conventional hexagonal honeycomb sandwich panel. The plastic energy consumption of the rear face-sheet of the BIH sandwich panel is lower than that of the conventional panel. Among them, the plastic energy consumption of the BIH is 8.02% higher than that of the hexagonal honeycomb. The damage to the rear face-sheet may be reduced considering that the front face-sheet and the honeycomb core of the BIH sandwich panel absorb most of the impact energy.
The deformation modes are given in Figure 8. The interaction between honeycomb cells mainly occurs near the local recessed area, and the buckling of the cells occurs in the fold area. Evidently, the interaction enhancement effect between the BIH cells is more obvious than that of the conventional hexagonal honeycomb, and the folding wavelength of the BIH fold area is smaller than that of the hexagonal honeycomb. Generally, the more interactions between honeycomb cells, the better the impact energy dissipation. The shorter folded wavelengths are beneficial to the formation of more folds, which causes more energy to be absorbed. This analysis explains why the BIH sandwich structure has better impact resistance than the conventional hexagonal honeycomb sandwich structure.

5. Parameter Study

In order to reveal the coupling effects of the structural parameters and their impact resistance, the influence of the circumscribed circle diameter ratio of the BIH trunks and branches (γ), the thickness ratio of trunks to branches (K) and the impact angle (θ) were investigated.

5.1. Effect of Ratio γ

Four different γ values were designed as follows: 0.1, 0.3, 0.5, 0.7. At the same time, the parameter K remained unchanged (K = 1). Figure 9 shows the maximum displacement and the plastic energy dissipation of the BIH sandwich structure with different γ ratios. The maximum displacements of front and rear face-sheets both increase and then decrease as γ increases. In particular, when γ is equal to 0.3, the front and rear face-sheets have the largest displacement. Figure 9b describes the plastic energy dissipation of the BIH sandwich structure. As γ increases, the plastic energy dissipation of the whole structure and the honeycomb core decreases and then increases, while the plastic energy dissipation law of the structural face-sheets is exactly the opposite to the previous two. It should be noted that when γ is 0.3, the plastic energy dissipation values of the overall structure and the honeycomb core are the lowest, while the plastic energy dissipation value of the face-sheet is the highest. This means that the honeycomb core has important impacts on its total plastic dissipated energy, which may further determine the impact resistance performance of this structure.
Therefore, Figure 10 exhibits the deformation modes of the BIH sandwich structure at 0.3 ms and 0.9 ms under different γ values. When the BIH sandwich structure is impacted by a solid sphere at low speed, the central area is in a depression state and the cell walls in this region buckle plastically. The contour line of cell buckling is enlarged inside the red dashed box. When γ is equal to 0.7, the folding wavelength is the shortest, followed by γ equal to 0.1. Interestingly, the structures with γ = 0.3 and γ = 0.5 have similar folding wavelengths. According to the plastic collapse of the cells in the orange dotted box, it can be clearly seen that the cell walls of the sandwich structure, with γ = 0.1, γ = 0.5 and γ = 0.7, gradually fold in an orderly manner, which is different from the irregular folding in the sandwich structure with γ = 0.3. Ordered cell folding means that more impact energy is effectively absorbed. Therefore, when designing the BIH sandwich energy-absorbing structure, it should be noted that γ equal to 0.3 is not conducive to impact resistance.

5.2. Effect of Ratio K

To analyze the influence of the K ratio on the dynamic response of the BIH sandwich panel, seven different K values were designed as follows: 0.4, 0.6, 0.8, 1, 3, 5, 7. At the same time, the parameter γ remained unchanged (γ = 0.5). Figure 11 shows the maximum displacement and the plastic energy dissipation of the BIH sandwich structure with different values of K. When K ≤ 1, the plastic energy consumption of the face-sheets and the honeycomb core of the BIH sandwich structure changes very little with the change in K. However, when K > 1, they all have significant differences with the change in K. The maximum displacement of both the front and rear face-sheets decreases and then increases with the increase in K. In particular, the rear face-sheet of the BIH sandwich structure has a minimum depth of depression when K = 0.6. This means that when K = 0.6, the BIH sandwich structure has the best impact resistance, and its deformation mode is shown in the red dashed box of the picture.

5.3. Effect of Impact Angle θ

Considering that sandwich panels may be subjected to impact loads from different directions, an analysis of the impact angle θ is of great significance. Four different θ values were designed as follows: 0°, 10°, 20°, 30°. At the same time, K and γ were set to 1 and 0.5, respectively. Figure 12 depicts the maximum displacement and the plastic energy dissipation of the BIH sandwich structure at different impact angles. The maximum displacement of the front and rear face-sheets gradually decreases as θ increases. In Figure 12b, the plastic energy dissipation of the front and rear face-sheets of the BIH sandwich panel gradually decreases with the increase in θ. However, the plastic dissipated energy of the BIH core is the opposite. The impact from different angles will cause varying degrees of damage to the BIH sandwich panel.
Figure 13 gives the deformation mode of the BIH sandwich structure for different values of θ. As the impact angle increases, the number of wrinkles at the local buckling of the sandwich structure increases. In addition, the interaction between the cells of the BIH gradually increases. Therefore, it can be concluded that the energy absorption characteristics of the structure are closely related to the impact angle.

6. Optimization Design

6.1. Determination of Optimization Problem

In order to obtain the BIH sandwich panel with the best impact resistance, the structural parameters of the BIH sandwich structure were further optimized. In the process of design optimization, the weight efficiency control of the BIH sandwich structure cannot be ignored. The dent depth and the mass of the BIH sandwich structure are represented by U and M, respectively. The response surface method [43,44,45] is used to optimize the design of the BIH sandwich panel. The expression of the optimization problem can be described as follows:
M i n i m i z e   [ U ( K , γ ) ] 0.4 K 1.7 0.1 γ 0.7   [ M ( K , γ ) ] 100   g

6.2. Response Surface Method

6.2.1. The Process of Design Optimization

The optimization design process using the response surface method is shown in Figure 14. The optimization problem is first determined. Then, the design of experiments is carried out according to the optimization problem to obtain the corresponding design sample points. The corresponding value of the objective function is obtained through finite element simulation. Subsequently, a response surface model is established through regression analysis based on the design sample points and their response values, and the optimal solution is then obtained from the response surface model. Finally, the optimal solution is tested. If the error between the optimal result obtained from the prediction model and the numerical simulation result is within the allowable range, the optimal solution is considered acceptable. If the error between them is large, the design domain needs to be adjusted, re-optimized and analyzed until the error is within the allowable range.

6.2.2. Design of Experiments and Metamodel

Typical design experiments include the Taguchi orthogonal experiment method [46,47], the central composite design, the full factor design and the Latin hypercube. The full factor design was adopted to generate 20 design sampling points (four levels for K, five levels for γ), and the response values of the selected sample points were obtained by simulation. Polynomial regression analysis was used to deal with the design optimization problems [48,49]. The least squares method was used to fit U and M. Figure 15 is the response surface and the predictive models of U and M are as follows:
U = 3.698 0.430 K + 2.957 γ 1.485 K γ + 0.817 K 2 1.984 γ 2
M = 130.298 32.632 K 41.814 γ + 21.091 K γ + 6.007 K 2 + 15.25 γ 2
Figure 16 compares the predicted value of the response surface model with the actual value. It can be seen that all the square dots fluctuate around the red diagonal line with a slope of 1. The predicted values of U and M are not significantly different from the actual values, which further verifies the accuracy of the response surface model.

6.3. Results and Discussion

In order to obtain the optimal BIH sandwich structure, Figure 17a shows a partial solution of the response surface model during the optimization process. The blue dot is the optimal design point to obtain the smallest U value, and its detailed cross-sectional parameters are given in Figure 17b. Interestingly, when the mass is less than 97.286 g, the dent depth of the rear face-sheet increases as the mass increases. However, when the mass is greater than 97.286 g, the dent depth of the rear face-sheet decreases as the mass increases. This shows that when the mass is greater than 97.286 g, reducing the mass of the BIH sandwich structure is conducive to the enhancement of its impact resistance. This is of great benefit in engineering applications.
The optimal solution of the response surface model is U = 4.135 mm, mass M = 97.286 g, and the corresponding ratios are K = 0.9 and γ = 0.7. The BIH sandwich structure with the same structural parameters must be simulated to further verify the credibility of the optimal result. The simulation value of U is 4.246 mm, and the error between prediction and simulation is 2.61% (Table 4), which is within an acceptable range. Therefore, the optimal solution obtained by the response surface model is highly reliable. The maximum rear face-sheet displacement U was used as the evaluation index for the sandwich panel’s impact resistance. It can be concluded that the impact resistance of the optimized sandwich panel increased by 5.16% compared to the un-optimized sandwich panel.

7. Conclusions

Inspired by the vein structure of a lotus leaf, a bio-inspired honeycomb (BIH) core sandwich structure was proposed to further enhance the impact resistance of conventional honeycomb sandwich structures. ABAQUS was used to establish and effectively verify the numerical model of the BIH sandwich structure. Then, the dynamic responses of the BIH sandwich panel and the conventional hexagonal honeycomb sandwich panel under low-velocity impact were compared and analyzed, followed by a detailed analysis of the coupling effect of structural parameters and impact resistance. Finally, the BIH sandwich panel was optimized using the response surface method.
Considering that the BIH is more prone to plastic deformation and there is more contact between cells, it may be concluded that the BIH has better impact resistance than conventional hexagonal honeycomb sandwich panels. Meanwhile, the impact resistance of the sandwich panel decreased and then increased as the circumscribed circle diameter ratio of the BIH trunk to branch (γ) increased; the sandwich panel with γ = 0.3 showed the worst impact resistance. However, the thickness ratio of trunk to branch (K) showed the opposite trend. The impact resistance increased and then decreased as K increased and the sandwich panel with K = 0.6 showed the best performance. As for the influence of the impact angle (θ), the impact resistance of the sandwich panel gradually increased as θ increased. The coupling effect between the structural parameters and the energy absorption capacity can provide a reference for design optimization and engineering applications.
Based on the sample points obtained by the full-factor design of the experiments and their response values, the prediction models of U and M were obtained by polynomial regression analysis. The optimization results were further tested, which verified the reliability of the optimization results obtained from the response surface model. These results can provide valuable suggestions for the research and design of new honeycomb sandwich structures.

Author Contributions

Conceptualization, Q.H.; Methodology, T.T. and Q.H.; Software, L.L. and J.G.; Data curation, Y.W.; Writing—original draft, L.L.; Writing—review & editing, T.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 51705215), the Chinese Postdoctoral Science Foundation (2022M712932) and the Natural Science Fundamental Research Project of Jiangsu Universities (No. 22KJA460003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

Author Tao Tao was employed by the company Guangzhou Metro Design & Research Institute Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Yahaya, M.A.; Ruan Lu, D.G.; Dargusch, M.S. Response of aluminium honeycomb sandwich panels subjected to foam projectile impact—An experimental study. Int. J. Impact Eng. 2015, 75, 100–109. [Google Scholar] [CrossRef]
  2. Wowk, D.; Reyno, T.; Yeung, R.; Marsden, C. An experimental and numerical investigation of core damage size in honeycomb sandwich panels subject to low-velocity impact. Compos. Struct. 2020, 254, 112739. [Google Scholar] [CrossRef]
  3. Alavi Nia, A.; Mokari, S.; Zakizadeh, M.; Kazemi, M. Experimental and numerical investigations of the effect of cellular wired core on the ballistic resistance of sandwich structures. Aerosp. Sci. Technol. 2017, 70, 445–452. [Google Scholar] [CrossRef]
  4. Vaziri, A.; Hutchinson, J.W. Metal sandwich plates subject to intense air shocks. Int. J. Solids Struct. 2007, 44, 2021–2035. [Google Scholar] [CrossRef]
  5. Palomba, G.; Epasto, G.; Crupi, V.; Guglielmino, E. Single and Double-Layer Honeycomb Sandwich Panels under Impact Loading. Int. J. Impact Eng. 2018, 121, 77–90. [Google Scholar] [CrossRef]
  6. Foster, L.; Peketi, P.; Allen, T.; Senior, T.; Duncan, O.; Alderson, A. Application of auxetic foam in sports helmets. Appl. Sci. 2018, 8, 354. [Google Scholar] [CrossRef]
  7. Dharmasenaa, K.P.; Wadleya HN, G.; Xue, Z.Y.; Hutchinson, J.W. Mechanical response of metallic honeycomb sandwich panel structures to high-intensity dynamic loading. Int. J. Impact Eng. 2008, 35, 1063–1074. [Google Scholar] [CrossRef]
  8. Bahei-El-Din, Y.A.; Dvorak, G.J. Enhancement of blast resistance of sandwich plates. Compos. Part B 2008, 39, 120–127. [Google Scholar] [CrossRef]
  9. Van Quyen, N.; Van Thanh, N.; Quan, T.Q.; Duc, N.D. Nonlinear forced vibration of sandwich cylindrical panel with negative Poisson’s ratio auxetic honeycombs core and CNTRC face sheets. Thin-Walled Struct. 2021, 162, 107571. [Google Scholar] [CrossRef]
  10. Abbadi, A.; Tixier, C.; Gilgert, J.; Azari, Z. Experimental study on the fatigue behavior of honeycomb sandwich panels with artificial defects. Compos. Struct. 2015, 120, 394–405. [Google Scholar] [CrossRef]
  11. Sezgin, F.E.; Tanoğlu, M.; Eğilmez, O.Ö.; Dönmez, C. Mechanical Behavior of Polypropylene-based Honeycomb-Core Composite Sandwich Structures. J. Reinf. Plast. Compos. 2010, 29, 1569–1579. [Google Scholar] [CrossRef]
  12. Baştürk, S.; Tanoğlu, M.; Çankaya, M.; Eğilmez, O. Dynamic behavior predictions of fiber-metal laminate/aluminum foam sandwiches under various explosive weights. J. Sandw. Struct. Mater. 2016, 18, 321–342. [Google Scholar] [CrossRef]
  13. Rejab MR, M.; Cantwell, W.J. The mechanical behavior of corrugated-core sandwich panels. Compos. Part B 2013, 47, 267–277. [Google Scholar] [CrossRef]
  14. Sun, Z.; Shi, S.S.; Guo, X.; Hu, X.Z.; Chen, H.R. On compressive properties of composite sandwich structures with grid reinforced honeycomb core. Compos. Part B-Eng. 2016, 94, 245–252. [Google Scholar] [CrossRef]
  15. Lan, X.K.; Feng, S.S.; Qi, H.; Zhou, T. A comparative study of blast resistance of cylindrical sandwich panels with aluminum foam and auxetic honeycomb cores. Aerosp. Sci. Technol. 2019, 87, 37–47. [Google Scholar] [CrossRef]
  16. Mcshane, G.; Deshpande, V.; Fleck, N. The underwater blast resistance of metallic sandwich beams with prismatic lattice core. J. Appl. Mech. 2007, 72, 352–367. [Google Scholar] [CrossRef]
  17. Sun, G.Y.; Huo, X.T.; Wang, H.X.; Hazell, P.J.; Li, Q. On the structural parameters of honeycomb-core sandwich panels against low-velocity impact. Compos. Part B 2021, 216, 108881. [Google Scholar] [CrossRef]
  18. Zhang, D.H.; Fei, Q.G.; Zhang, P.W. Drop-weight impact behavior of honeycomb sandwich panels under a spherical impactor. Compos. Struct. 2017, 168, 633–645. [Google Scholar] [CrossRef]
  19. Zhang, D.; Jiang, D.; Fei, Q.; Wu, S. Experimental and numerical investigation on indentation and energy absorption of a honeycomb sandwich panel under low velocity impact. Finite Elem. Anal. Des. 2016, 117–118, 21–30. [Google Scholar] [CrossRef]
  20. Liu, J.F.; Chen, W.S.; Hao, H.; Wang, Z.G. Numerical study of low-speed impact response of sandwich panel with tube filled honeycomb core. Compos. Struct. 2019, 220, 736–748. [Google Scholar] [CrossRef]
  21. Reyes, A.; Børvik, T. Low velocity impact on crash components with steel skins and polymer foam cores. Int. J. Impact Eng. 2019, 132, 1879–3509. [Google Scholar] [CrossRef]
  22. Zhou, H.Y.; Jia, K.C.; Wang, X.J.; Xiong, M.-X.; Wang, Y. Experimental and numerical investigation of low velocity impact response of foam concrete filled auxetic honeycombs. Thin-Walled Struct. 2020, 154, 106898. [Google Scholar] [CrossRef]
  23. Dhanarasu, M.; Vm, S. Structural response of sandwich structures with CFRP face sheets under quasi-static indentation and high velocity impact: An experimental and numerical study. Proc. Inst. Mech. Eng. Part L J. Mater. Des. Appl. 2023, 237, 737–752. [Google Scholar] [CrossRef]
  24. Sabah SH, A.; Kueh AB, H.; Al-Fasih, M.Y. Bio-inspired vs. conventional sandwich beams: A low-velocity repeated impact behavior exploration. Constr. Build. Mater. 2018, 169, 193–204. [Google Scholar] [CrossRef]
  25. Akatay, A.; Bora, M.; Çoban, O.; Fidan, S.; Tuna, V. The influence of low velocity repeated impacts on residual compressive properties of honeycomb sandwich structures. Compos. Struct. 2015, 125, 425–433. [Google Scholar] [CrossRef]
  26. Zhang, Y.; Yan, L.; Zhang, C.; Guo, S. Low-velocity impact response of tube-reinforced honeycomb sandwich structure. Thin-Walled Struct. 2021, 158, 107188. [Google Scholar] [CrossRef]
  27. Al Antali, A.; Umer, R.; Zhou, J.; Cantwell, W.J. The energy-absorbing properties of composite tube-reinforced aluminum honeycomb. Compos. Struct. 2017, 176, 630–639. [Google Scholar] [CrossRef]
  28. Hundley, J.M.; Clough, E.C.; Jacobsen, A.J. The low velocity impact response of sandwich panels with lattice core reinforcement. Int. J. Impact Eng. 2015, 84, 64–77. [Google Scholar] [CrossRef]
  29. Baran, T.; Öztürk, M. In-plane elasticity of a strengthened re-entrant honeycomb cell. Eur. J. Mech.-A/Solids 2020, 83, 104037. [Google Scholar] [CrossRef]
  30. Hassan, M.Z.; Cantwell, W.J. The influence of core properties on the perforation resistance of sandwich structures—An experimental study. Compos. Part B-Eng. 2012, 43, 3231–3238. [Google Scholar] [CrossRef]
  31. Tatlıer, M.S.; Öztürk, M.; Baran, T. Linear and non-linear in-plane behaviour of a modified re-entrant core cell. Eng. Struct. 2021, 234, 111984. [Google Scholar] [CrossRef]
  32. San Ha, N.; Pham, T.M.; Tran, T.T.; Hao, H.; Lu, G. Mechanical properties and energy absorption of bio-inspired hierarchical circular honeycomb. Compos. Part B Eng. 2022, 236, 109818. [Google Scholar]
  33. Al-Dhaheri, M.; Khan, K.A.; Umer, R.; van Liempt, F.; Cantwell, W.J. Process-induced Deformation in U-shaped Honeycomb Aerospace Composite Structures. Compos. Struct. 2020, 248, 112503. [Google Scholar] [CrossRef]
  34. Ha, N.S.; Jin, T.L.; Goo, N.S.; Park, H.C. Anisotropy and non-homogeneity of an Allomyrina Dichotoma beetle hind wing membrane. Bioinspir. Biomim. 2011, 6, 046003. [Google Scholar] [CrossRef] [PubMed]
  35. Ha, N.S.; Le, V.T.; Goo, N.S. Investigation of punch resistance of the Allomyrira dichtoloma beetle forewing. J. Bionic Eng. 2018, 15, 57–68. [Google Scholar] [CrossRef]
  36. San Ha, N.; Lu, G.; Shu, D.W.; Yu, T.X. Mechanical properties and energy absorption characteristics of tropical fruit durian (Durio zibethinus). J. Mech. Behav. Biomed. Mater. 2020, 104, 103603. [Google Scholar]
  37. San Ha, N.; Pham, T.M.; Hao, H.; Lu, G. Energy absorption characteristics of bio-inspired hierarchical multi-cell square tubes under axial crushing. Int. J. Mech. Sci. 2021, 201, 106464. [Google Scholar]
  38. San Ha, N.; Pham, T.M.; Chen, W.; Hao, H.; Lu, G. Crashworthiness analysis of bio-inspired fractal tree-like multi-cell circular tubes under axial crushing. Thin-Walled Struct. 2021, 169, 108315. [Google Scholar]
  39. Song, J.F.; Xu, S.C.; Liu, S.F.; Huang, H.; Zou, M. Design and numerical study on bionic columns with grooves under lateral impact. Thin-Walled Struct. 2020, 148, 106546. [Google Scholar] [CrossRef]
  40. Zou, M.; Xu, S.C.; Wei, C.G.; Wang, H.; Liu, Z. A bionic method for the crashworthiness design of thin-walled structures inspired by bamboo. Thin-Walled Struct. 2016, 101, 222–230. [Google Scholar] [CrossRef]
  41. Palombini, F.L.; Araujo Mariathb, J.E.; Oliveiraa, B.F. Bionic design of thin-walled structure based on the geometry of the vascular bundles of bamboo. Thin-Walled Struct. 2020, 155, 106936. [Google Scholar] [CrossRef]
  42. He, Q.; Wang, Y.; Gu, H.; Feng, J.; Zhou, H. Dynamic crushing analysis of a circular honeycomb with leaf vein branched characteristic. Mech. Mater. 2021, 153, 103566. [Google Scholar] [CrossRef]
  43. Kazemi, M.; Serpoush, J. Energy absorption parameters of multi-cell thin-walled structure with various thicknesses under lateral loading. Mater. Des. Appl. 2021, 235, 513–526. [Google Scholar] [CrossRef]
  44. Li, L.H.; Guo, C.; Chen, Y.T.; Chen, Y.H. Optimization design of lightweight structure inspired by glass sponges (Porifera, Hexacinellida) and its mechanical properties. Bioinspir. Biomim. 2020, 15, 036006. [Google Scholar] [CrossRef]
  45. Liu, X.; Liang, R.; Hu, Y.Z. Body optimization approach of sedan structure for improving small overlap impact rating. Int. J. Crashworthiness 2019, 26, 182–190. [Google Scholar] [CrossRef]
  46. Rathod, S.; Tiwari, G. Structural analysis of cylindrical composite structures for high velocity applications: A parametric optimization study. Thin-Walled Struct. 2023, 192, 111178. [Google Scholar]
  47. Benyounis, K.Y.; Olabi, A.G.; Hashmi, M.S.J. Ulti-response optimization of CO2 laser-welding process of austenitic stainless steel. Opt. Laser Technol. 2008, 40, 76–87. [Google Scholar] [CrossRef]
  48. Johnson, W.R.; Zhu, X.-K.; Sindelar, R.; Wiersma, B. A parametric finite element study for determining burst strength of thin andthick-walled pressure vessels. Constr. Build. Mater. 2023, 204, 104968. [Google Scholar]
  49. Rajput, A.; Sunny, M.R.; Sarkar, A. Optimization of honeycomb parameters of sandwich composites for energy and specific energy absorption using particle swarm optimization. Mar. Struct. 2023, 92, 103498. [Google Scholar] [CrossRef]
Figure 1. The design process and geometric configuration of the bio-inspired honeycomb sandwich structure: (a) geometry of lotus leaf veins; (b) BIH core sandwich structure; (c) cell geometry; (d) cross-sectional evolution of BIH cells.
Figure 1. The design process and geometric configuration of the bio-inspired honeycomb sandwich structure: (a) geometry of lotus leaf veins; (b) BIH core sandwich structure; (c) cell geometry; (d) cross-sectional evolution of BIH cells.
Materials 17 01191 g001
Figure 2. FE model of BIH under low-velocity impact.
Figure 2. FE model of BIH under low-velocity impact.
Materials 17 01191 g002
Figure 3. Meshing of the finite element model.
Figure 3. Meshing of the finite element model.
Materials 17 01191 g003
Figure 4. Comparison of experimental [18] and numerical results: (a) deformation mode; (b) dent depth and energy absorption.
Figure 4. Comparison of experimental [18] and numerical results: (a) deformation mode; (b) dent depth and energy absorption.
Materials 17 01191 g004
Figure 5. Comparison of time history curves of impact load of sandwich structure.
Figure 5. Comparison of time history curves of impact load of sandwich structure.
Materials 17 01191 g005
Figure 6. Comparison of maximum displacement of sandwich structure face-sheets: (a) conventional hexagonal honeycomb sandwich structure; (b) BIH sandwich structure.
Figure 6. Comparison of maximum displacement of sandwich structure face-sheets: (a) conventional hexagonal honeycomb sandwich structure; (b) BIH sandwich structure.
Materials 17 01191 g006
Figure 7. Comparison of plastic energy consumption of the hexagonal honeycomb and BIH sandwich structures.
Figure 7. Comparison of plastic energy consumption of the hexagonal honeycomb and BIH sandwich structures.
Materials 17 01191 g007
Figure 8. Deformation mode of sandwich structures with partially enlarged areas.
Figure 8. Deformation mode of sandwich structures with partially enlarged areas.
Materials 17 01191 g008
Figure 9. Maximum displacement and plastic energy dissipation of BIH sandwich structure with different γ ratios: (a) maximum displacement; (b) plastic energy dissipation.
Figure 9. Maximum displacement and plastic energy dissipation of BIH sandwich structure with different γ ratios: (a) maximum displacement; (b) plastic energy dissipation.
Materials 17 01191 g009
Figure 10. Deformation mode of BIH sandwich structure under different γ ratios.
Figure 10. Deformation mode of BIH sandwich structure under different γ ratios.
Materials 17 01191 g010
Figure 11. The maximum displacement and plastic energy dissipation of BIH sandwich structure under different values of ratio K.
Figure 11. The maximum displacement and plastic energy dissipation of BIH sandwich structure under different values of ratio K.
Materials 17 01191 g011
Figure 12. Maximum displacement and plastic energy dissipation of BIH sandwich structure under different values of θ: (a) maximum displacement; (b) plastic energy dissipation.
Figure 12. Maximum displacement and plastic energy dissipation of BIH sandwich structure under different values of θ: (a) maximum displacement; (b) plastic energy dissipation.
Materials 17 01191 g012
Figure 13. Deformation mode of BIH sandwich structure under different values of θ.
Figure 13. Deformation mode of BIH sandwich structure under different values of θ.
Materials 17 01191 g013
Figure 14. Optimization process of response surface method.
Figure 14. Optimization process of response surface method.
Materials 17 01191 g014
Figure 15. (a) The response surface of U; (b) the response surface of M.
Figure 15. (a) The response surface of U; (b) the response surface of M.
Materials 17 01191 g015
Figure 16. Comparison of predicted and actual results: (a) dent depth, U; (b) structural mass, M.
Figure 16. Comparison of predicted and actual results: (a) dent depth, U; (b) structural mass, M.
Materials 17 01191 g016
Figure 17. (a) The solutions of the response surface model; (b) optimal BIH cell cross section.
Figure 17. (a) The solutions of the response surface model; (b) optimal BIH cell cross section.
Materials 17 01191 g017
Table 1. Structural parameters of the BIH sandwich panel.
Table 1. Structural parameters of the BIH sandwich panel.
Face-Sheets (AL5083)BIH Core (AL3003)
DimensionThicknessKγ
150   m m × 150   m m × 16   m m 0.5   m m 0.4; 0.6; 0.8; 1; 3; 5; 70.1; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7
Table 2. Material parameters [18].
Table 2. Material parameters [18].
ρ   ( k g / m 3 ) E   ( G P a ) v Y i e l d   S t r e s s   σ y  (MPa) T a n g e n t   M o d u l u s   E t  (MPa)
AL5083 (face-sheets)2700720.33280933
AL3003 (BIH core)2700700.33185720
Table 3. Structural parameters of the honeycomb core.
Table 3. Structural parameters of the honeycomb core.
Hexagonal Honeycomb CoreBIH Core
Cell wall lengthWall thicknessΦbtbKγ
6 mm0.193 mm12 mm0.075 mm10.5
Table 4. Comparison of optimal solution and numerical result.
Table 4. Comparison of optimal solution and numerical result.
Optimal SolutionFEMError
U (mm)4.1354.2462.61%
M (g)97.286980.729%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tao, T.; Li, L.; He, Q.; Wang, Y.; Guo, J. Mechanical Behavior of Bio-Inspired Honeycomb–Core Composite Sandwich Structures to Low-Velocity Dynamic Loading. Materials 2024, 17, 1191. https://doi.org/10.3390/ma17051191

AMA Style

Tao T, Li L, He Q, Wang Y, Guo J. Mechanical Behavior of Bio-Inspired Honeycomb–Core Composite Sandwich Structures to Low-Velocity Dynamic Loading. Materials. 2024; 17(5):1191. https://doi.org/10.3390/ma17051191

Chicago/Turabian Style

Tao, Tao, Lizheng Li, Qiang He, Yonghui Wang, and Junlan Guo. 2024. "Mechanical Behavior of Bio-Inspired Honeycomb–Core Composite Sandwich Structures to Low-Velocity Dynamic Loading" Materials 17, no. 5: 1191. https://doi.org/10.3390/ma17051191

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop