Next Article in Journal
Enhanced Adhesion of Synthetic Discs with Micro-Patterned Margins
Next Article in Special Issue
Has Biomimicry in Architecture Arrived in France? Diversity of Challenges and Opportunities for a Paradigm Shift
Previous Article in Journal
Theoretical Design of a Bionic Spatial 3D-Arrayed Multifocal Metalens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy Absorption Characteristics of Polygonal Bio-Inspired Honeycomb Column Thin-Walled Structure under Quasi-Static Uniaxial Compression Loading

1
College of Civil and Architectural Engineering, Heilongjiang Institute of Technology, Harbin 150050, China
2
School of Civil Engineering, Northeast Forestry University, Harbin 150040, China
3
Longjian Road and Bridge Company Limited, Harbin 150009, China
*
Author to whom correspondence should be addressed.
Biomimetics 2022, 7(4), 201; https://doi.org/10.3390/biomimetics7040201
Submission received: 1 November 2022 / Revised: 13 November 2022 / Accepted: 14 November 2022 / Published: 17 November 2022
(This article belongs to the Special Issue Biomimetic Architectural and Urban Design 2.0)

Abstract

:
In this paper, we investigated the internal structure of the beetle elytra, i.e., two different structural forms I and II of the bio-inspired honeycomb column thin-walled structures (BHTS) that give the honeycomb sandwich structure frequently used in construction projects better mechanical properties and lightweight performance. BHTS specimens were fabricated by additive manufacturing selective laser melting (AM-SLM) using AlSi10Mg Al-Si alloy. In order to understand the effect of section angle number on BHTS during loading, quasi-static uniaxial compression tests were carried out and verified by numerical simulation. The experimental results showed that in the quasi-static uniaxial compression tests, the number of section angles greatly improved the energy absorption (EA) index of each BHTS: the average value of initial peak crushing force (PCF) of hexagonal BHTS increased by 108.82% and 43.44%, respectively, compared to triangular and rectangular BHTS. The average value of the mean crushing force (MCF) increased by 74.87% and 45.48%, respectively. The average value of EA increased by 89.02% and 46.64%, respectively. The results indicate that the number of section angles can be used as an effective way to enhance the EA of BHTS. This work can provide a reference for the design of high-efficiency energy absorbers and will be widely used in EA scenarios such as construction, transportation, etc.

1. Introduction

Thin-walled structures offer the advantages of high mechanical properties and light weight [1,2]. In recent years, thin-walled structures have been found to have good crashworthiness during collapsing [3,4]; therefore, they are widely used in civil construction, transportation, aerospace, and other industries for EA and protection devices. It has been shown that the cross-sectional structure is one of the key factors affecting the EA performance of thin-walled structures [5,6]. The crashworthiness of thin-walled structures with different cross-sectional configurations, including circular [7,8], square [9,10], pentagon [11], tapered tubes [12,13], honeycomb [14,15], multi-cell [16,17], and crisscross tubes [18] has been systematically investigated in detail. In these studies, the most severe plastic deformation occurred near the section angle of the tube, through the membrane deformation and bending deformation along the bending hinge line, in which a lot of energy can be consumed. For example, Ali et al. [11] conducted theoretical, numerical, and experimental studies on the dynamic axial crushing of pentagonal and cross-shaped tubes. They found that the energy absorbed by aluminum and steel pentagonal and crossed tubes was 60% and 92% higher than that of rectangular tubes in the same mass, respectively. Wu et al. [19] found that the number of cells in the cross-section can significantly determine the EA capacity. It can be seen that increasing the section angle and the number of cells in the section is an important way to enhance the EA of the structure from the study above.
Important ways to enhance structural EA, such as increasing the section angle and the number of cells within the cross-section, are also reflected in many biological structures. Inspired by the biological inspiration of nature, many scholars have designed some hierarchical structures to enhance EA capacity. Hu et al. [20] developed a bionic honeycomb tubular nested structure (BHTNS) inspired by the micro-architecture of bamboo vascular bundles. The results showed that BHTNS had excellent specific energy absorption (SEA) and provided a reference for the design of efficient energy absorbers. Inspired by the internal structure of ladybugs, Xiang and Du [21] proposed a novel honeycomb structure BHTS by filling the column in a different way. The results showed that the EA characteristics of the bionic structure were better than those of the bionic honeycomb structure (BHS), which is a conventional honeycomb structure without columns. This honeycomb structure is expected to be used as an alternative structure in automobiles to improve crashworthiness. Hao and Du [22] proposed three bionic honeycomb thin-walled structures based on the internal structure of beetle elytra (Figure 1). The results showed that the SEA of each designed BHTS was greater than that of BHS and showed that BHTS represents an effective improvement in the enhanced EA capacity of honeycomb structures and potential applications in the field of protective equipment.
However, the design of the bionic honeycomb structure of the above investigation is complex, and the production often encounters difficult technical problems that have certain manufacturing difficulties using traditional methods. In recent years, additive manufacturing technology (AM) has become an excellent choice for manufacturing cellular structures with precise porosity and pore size due to its advantages of high precision in manufacturing parts and high material utilization [23,24], which stimulate the research of cellular structures in the AM area [25,26,27]. Xia et al. [28] used the AM-SLM technique to fabricate BHTS for uniaxial compression tests. Ngoc San Ha et al. [29] used fused deposition modeling (FDM) and fabricated a new bio-inspired hierarchical circular honeycomb (BHCH) for uniaxial compression tests. Good results were achieved by the experimental structures and numerical simulation, indicating the high feasibility of AM technology in manufacturing complex cellular structures.
In this study, the three BHTS structures proposed by Hao and Du, which mimic the internal hollow column structure of beetle elytra, were made polygon extensions to propose a novel structural design of triangular, rectangular BHTS. Motivated by these promising findings, the effect of the section angle on the EA of the BHTS was discussed through a series of low-velocity compression experiments and numerical simulations based on the internal structural features of the beetle elytra in previous work [28]. Section 2 briefly describes the new multi-deformation expansion design, the loading scheme with BHTS, and the corresponding numerical simulations. Then, Section 3 presents the results of quasi-static uniaxial compression experiments to obtain the EA behavior of BHTS under axial crushing conditions, and the effect of section angle on the EA performance of BHTS is discussed and verified by experimental results. Finally, some useful conclusions are drawn in Section 4. The results of this study can provide some references for the design of similar structures.

2. Materials and Methods

2.1. Specimen Design and Loading Method

2.1.1. Tensile Test

Al-Si alloys have low density, high specific strength, good casting performance, and many excellent mechanical properties and are widely used in aerospace and automotive applications [30,31]. In order to determine the mechanical properties of BHTS material (AM-SLM AlSi10Mg), take the mechanical tests of porous structures fabricated through FDM by Ngoc San Ha et al. [29] for reference to test. The constitutive relation of nylon material was determined by a uniaxial tensile test, and numerical simulation was carried out using the obtained test data. The results showed that the finite element model (FEM) could effectively reflect the EA characteristics. According to the test concepts of the above scholars, the uniaxial tensile test was carried out. Five tensile dog-bone rod specimens were designed according to the ASTM E8/E8M standard [31], 200 mm long, 20 mm wide at both ends, 12.5 mm wide in the middle, with a circular arc diameter of 63.5 mm, and a 2 mm diameter hole in the middle (Figure 2a). By using SLM technology with high part-making accuracy and material utilization, each specimen was produced by an industrial-grade iSLM-280 professional metal printer. The tensile test was carried out on a Sanso Vertical 50 kN UTM5504GD universal testing machine with a controlled tensile speed of 0.2 mm/s (Figure 2b). Each sample was cleaned with ethanol before loading to remove dust and other potential contamination on its surface.

2.1.2. Specimen Parameters and Numbering

The printed thin-walled structures models (Figure 3) were available in three configurations: triangular (T), rectangular (R), and hexagonal (H). The thickness T c of the manufactured thin-walled structure was 0.5 mm uniformly. The height H c of all specimens was 20 mm. The diameter D c of the circular thin-walled column structure was 0.5 mm. The side length L c of all polygons was 6 mm. The design of BHTS was based on the thin-walled honeycomb structure with columns of the BHS, and BHTS–S (BHS) was the control group. The BHTS–I and BHTS–II hollow columns were added to the center of the honeycomb wall and the wall connection, respectively. For simplicity, the manufactured thin-walled structures with triangular cross-sections were named as BHTS–S–T, BHTS–Ⅰ–T, BHTS–Ⅱ–T. By analogy, rectangles and hexagons were denoted as BHTS–S–R, BHTS–Ⅰ–R, BHTS–Ⅱ–R; BHTS–S–H, BHTS–Ⅰ–H, BHTS–Ⅱ–H, respectively, as shown in Table 1.
Figure 4 shows the design models made by AM technology. The compression load was applied with a constant speed of 0.02 mm/s by using the Sanso Vertical 50 kN UTM5504GD universal testing machine. All nine groups of specimens were crushed to 75% of their original length at the end of the compression tests, and each group contained two relatively close valid data.

2.2. Finite Element Models

2.2.1. Modeling Process

The explicit finite element commercial software LS-DYNA [32] was used for the out-of-plane crush loading of BHTS. The AlSi10Mg material was simulated by the Cowper–Symonds isotropic hardening models. Considering the strain rate effect [33,34], the model expression is:
σ Y = 1 + ε ˙ C 1 P σ 0 + β E P ε P e f f
In Equation (1), σ 0 is the initial yield stress; ε ˙ is the strain rate; β is the hardening coefficient; C and P are Cowper–Symonds strain rate parameters; ε P e f f is the effective plastic strain and, E P is the plastic hardening modulus.
In order to simulate the tensile process of AlSi10Mg tensile specimens, a FEM for low-velocity tensile analysis was established, as shown in Figure 5. Using coarse grids at the fixture end (tensile end and fixed end; mesh size was 5 mm), and finer grids were used in the fracture zone (mesh size was 1.25 mm) and transition zone (mesh size was 2 mm) to cope with larger plastic deformation in the finite element simulation (Figure 5a). The free translational degrees of the nodes on the contact surface between the fixed-end dog-bone rod and the fixture in three directions were constrained. The parameters in the FEM were adjusted according to the data collected from the tensile tests to be close to the test values and feedback (Figure 5b). It can be seen from the tests that the yield stress of the material fluctuated slightly at the level of 220 MPa, and other material characteristics from the experiments are shown in Table 2.

2.2.2. Uniaxial Compression Test Models

In the BHTS uniaxial compression tests, taking the loading process of BHTS–S–T as an example (Figure 6), the FEM consists of three parts: the fixed support rigid plate, the upper impact rigid plate, and the BHTS in between. The material of BHTS is *MAT_PLASTIC_KINEMATIC. The upper impact plate and the lower support plate are defined as rigid bodies by *MAT_RIGID. The rigid plate element was meshed by a constant stress body element, and the side length of the mesh was 1 mm. In the calculation, *DEFINE_CURVE was used to set the motion path, and *BOUNDARY_PRESCRIBED_MOTION_SET was used to define a time-varying linear displacement for the upper impact plate to ensure its uniform motion. The upper rigid plate constrained all displacements except the impact direction, and the lower rigid plate constrained all displacements. The complete integral Belytschko-Tsay [35,36] membrane element was used in the honeycomb aluminum tube. The honeycomb aluminum tube adopted a quadrilateral grid with a side length of 0.5 mm. Honeycomb structure model self-contact used single-sided automatic contact algorithm *CONTACT_AUTOMATIC_SINGLE_SURFACE honeycomb structure and support plate surface contact type selected point surface automatic contact algorithm *CONTACT_AUTOMATIC_NODES_TO_SURFACE [37], and the same algorithm was used for the contact type between the honeycomb structure and the impact plate.

2.3. Image Acquisition System

In order to capture the complete deformation process of BHTS in each test for a more comprehensive understanding of the EA characteristics of BHTS under quasi-static compression, an industrial-grade global shutter high-speed camera MV-SUA133GC-T with a maximum frame rate of 245, equipped with an industrial-grade HD lens MV-LD-8-4M-G with a focal length of 8 mm was used to record the deformation process of the complete structure (Figure 7).

3. Results and Discussion

3.1. Experimental Phenomena

Two panels of close valid data were used for each sample. The classification explanation of each panel corresponding to each group is shown in Table 3. Figure 8 shows the loading process of BHTS under quasi-static uniaxial compression. The quasi-static out-of-plane uniaxial compression behavior of BHTS was analyzed. With the increase of compression displacement, the BHTS–S–T1 specimen was in the elastic deformation stage, and the middle began to expand. At the same time, the lateral deformation of the specimen showed regularities. As the displacement continued to increase, the specimen was in the elastic nonlinear deformation stage, and the lateral deformation increased greatly. Meanwhile, local bending buckling occurred on both sides of the specimen ribs. As displacement continued to increase, the specimen began to show more obvious shear failure followed by collapsing. The lower middle and upper of the specimen reached the failure strain and gradually crushed into the folded compact stage. In the folded and broken cycle, the specimen was finally compressed into a small broken piece. With the increase of compression displacement, the ribs and the cell walls on both sides of BHTS–I–T1 buckled at the same time and local bending buckling occurred, accompanied by a more obvious torsional response, then the local torsional-flexural buckling led to the collapse of the specimen. The ribs of BHTS–II–T1 buckled on both sides at the same time. After buckling, both sides of the ribs collapsed after obvious symmetrical shear deformation. BHTS–S–R1 had obvious local bending buckling at the top during compression. As the displacement increased, the top collapsed after shear failure. In BHTS–I–R1, with the loading of displacement, local torsional buckling occurred earlier, followed by local torsional–flexural buckling, and then collapsed after shear failure. BHTS–Ⅱ–R1 showed local buckling at the top with the increase of displacement and collapsed after a more obvious shear failure with the increase of displacement. BHTS–S–H1 collapsed after local bending buckling at the bottom with increased displacement, followed by single-celled shear failure. With the increase of displacement, local torsional–flexural buckling occurred in the middle and bottom of BHTS–I–H1. As the displacement continued to increase, single-celled and continuous shear failure occurred and then collapsed. With the increase of displacement, the middle and bottom of the BHTS–Ⅱ–H1 local bending buckling occurred, and with the continued increase of displacement, a single-celled with a continuous form of shear failure collapsed in the BHTS–Ⅱ–H1 later.

3.2. Typical Impact Force-Displacement Response Curves

The force-displacement curves are used to characterize the mechanical and EA properties of the thin-walled structure, as shown in Figure 9. It can be seen that under dynamic impact loading conditions, the crushing force and EA have four different stages: elastic stage, plastic deformation stage, fracture stage, and densification stage. Taking BHTS–S–T1 as an example, it had a linear elastic deformation at 0–0.32 mm during loading. After exceeding the proportional limitation, the sample had obvious nonlinear elastic deformation at 0.32–1.86 mm due to the inevitable lateral displacement during loading. Linear elastic deformation and nonlinear elastic deformation are collectively referred to as the elastic deformation stage, where the crushing force increases rapidly to reach the initial PCF. The curve has no hardening stage; when the loading displacement reached 1.86 mm, local bending buckling and shear failure occurred and then collapsed. This stage is called plastic deformation and fracture, in which the crushing force decreases sharply. Entering the folding phase, the initial folding occurred after the displacement loading at 2.07 mm, and the specimen folded up and became densified into the progressive failure mode of strength failure and buckling coexisting. The strength curve rose at densification, decreased after buckling and strength failure, and repeatedly cycled until all compression was completed, which is called the densification stage. The progressive collapse mode phase after structural buckling was accompanied by plastic deformation and fracture. At this stage, the crushing force fluctuated periodically, and the crushing displacement was longer than the elastic and yielding stages. In addition, this stage achieved most of the EA. Finally, with the saturation of the fragment cracked, the structure was affected by the densification, and the crushing force of the force-displacement curve was on the rise.
Figure 10a,b compares the average PCF with MCF and the finite element simulation results of BHTS–S–T, BHTS–I–T, and BHTS–II–T specimens. It can be seen that the average PCF of all specimens increased with the increase of the section angle. Among the PCF of the three BHTS structures (S/I/II) in BHTS–H was 93.90%, 95.40%, 137.16%, and 10.42%, 50.84%, 69.07% larger than those of BHTS–T and BHTS–R, respectively. In the MCF analysis, except that BHTS–S–R was slightly smaller than BHTS–S–T, the average MCF of all specimens also increased with the increase of the number of section angles. Among the three BHTS structures (S/I/II), the PCF distribution of BHTS–Specimens–H was 43.98%, 71.78%, 108.87%, and 73.74%, 42.27%, 20.45% larger than that of BHTS–Specimens–T and BHTS–Specimens–R.
In addition, the experimental results were not much different from the numerical simulation results. Under quasi-static compression, the experimental mean values of PCF and MCF of triangle, rectangle, and hexagon were compared with the numerical results, and the experimental errors of PCF were 8.17%, 1.66%, and 5.12%, respectively. The test errors of MCF were 9.94%, 6.62%, and 5.58%, respectively, and the error mainly came from the inevitable lateral displacement during the loading process. It can be seen that the numerical simulation results are in good agreement with the experimental data, and the EA behavior of BHTS under impact loading was well simulated. It can be concluded that the present numerical models, including the adopted material models and the corresponding numerical algorithms, can effectively analyze the mechanical behavior of BHTS under quasi-static uniaxial compressive loading.

3.3. Energy Absorption Indexes

In this study, all crushing parameters were calculated directly from the force-displacement curve. PCF [38] is defined as the first peak load. MCF [39] was calculated as the average crushing force over a displacement range of 3 mm to 15 mm. The formula for calculating MCF is:
M C F = 1 b a a b F d δ
In Equation (2), a is 3 mm and b is 15 mm, wherein, d δ is the compression displacement increment at the beginning of the densification state, and F is the compression load.
EA capacity is an important index to characterize the performance of porous materials. The mechanical properties of porous materials require compression shock absorption and EA. In order to study the EA performance of different configurations, two typical indicators were used as evaluation criteria: crushing force efficiency (CFE) [38,40] and SEA [41]. The EA indexes of different types of BHTS under different low-strain rates were evaluated. CFE is defined as the ratio of MCF to PCF, as shown in Equation (3):
C F E = M C F P C F
EA [42] is the energy dissipated during the crushing process of the specimens. According to Equation (3), the area under the force-displacement curve is determined and calculated based on the displacement from 0 to 15 mm.
E A = 0 b F d δ
Due to the high strength-to-weight ratio of cellular structures, it is important to consider the energy absorbed per unit mass of the material. This is also known as SEA, which could be gotten by the dissipated E A divided by the specimen mass ( m ), as shown in Equation (5).
S E A = E A m
Generally speaking, the higher the SEA, the better the EA capacity of the BHTS is. In this paper, the commercial mathematical software Matlab [43] was used to compile the calculation program to summarize the data after calculation. Table 4 summarizes the PCF, MCF, CFE, EA, and SEA of all nine groups of samples. To facilitate comparison, we calculated the mean values of MCF, CFE, EA, and SEA for repeated samples.

3.4. Effect of Section Angle Number on EA in the Crushing Process

The EA and SEA diagrams for each structure are shown in Figure 11. It can be seen that the EA of the designed BHTS rose with the increase of the section angle. In the BHTS–S–Specimens structure, the EA of BHTS–S–H was 76.80% and 64.79% larger than that of BHTS–S–T and BHTS–S–R, respectively. In the BHTS–Ⅰ–Specimens structure, the EA of BHTS–Ⅰ–H was 73.14% and 43.66% larger than that of BHTS–Ⅰ–T and BHTS–Ⅰ–R, respectively; In the BHTS–II–Specimens structure, the EA of BHTS–Ⅱ–H was 117.12% and 31.48% larger than that of BHTS–Ⅱ–T and BHTS–Ⅱ–R, respectively.
The SEA of the designed BHTS did not show a significant correlation with the section angle. The SEA of BHTS–S–H was −17.65% and 3.73% larger than that of BHTS–S–T and BHTS–S–R, respectively. The SEA of BHTS–S–H was −20.93% and −9.45% larger than that of BHTS–S–T and BHTS–S–R, respectively. The SEA of BHTS–S–H was −4.38% and −20.63% larger than that of BHTS–S–T and BHTS–S–R, respectively.
From the above analysis, it can be seen that the EA of the designed BHTS grew linearly with the increase in the number of section angles. In contrast, the SEA of BHTS did not have a significant correlation with the section angle, and there were two possible reasons for this phenomenon. (1) The honeycomb structure is convergent only after a certain number of cellular pores is reached [44]; in this paper, the number of four-cell elements was less according to the design that just meets the Chinese specification. (2) The size of the designed BHTS was small due to size effect [45]; more cell elements and larger sizes should be designed during the subsequent studies to find more accurate patterns.

3.5. Effect of Hollow Column Position on EA

The EA and SEA diagrams for each structure are shown in Figure 12. It can be seen that the EA of the designed BHTS shows a relatively obvious increasing trend with the installation position of the hollow column. In the BHTS–Specimens–T structure, the EA of BHTS–II–T was 77.32% and −23.85% larger than that of BHTS–S–T and BHTS–I–T, respectively. In the structure of BHTS–Specimens-R, the EA of BHTS–S–H was 99.17% and 4.33% larger than that of BHTS–S–R and BHTS–Ⅰ–R, respectively. In the structure of BHTS–I, the EA of BHTS–S–H was 99.17% and 4.33% larger than that of BHTS–S–T and BHTS–S–R, respectively.
The SEA of the designed BHTS also showed a significant upward relationship with the increase of the hollow column. In the BHTS–S–Specimens structure, the EA of BHTS–S–H was larger than that of BHTS–S–T and BHTS–S–R; 40.11% and 16.26%, respectively. The EA of BHTS–Ⅰ–H was 54.11% and 15.52% larger than those of BHTS–Ⅰ–T and BHTS–Ⅰ–R, respectively; The EA of BHTS–II–H was 34.53% and 1.26% larger than that of BHTS–II–T and BHTS–II–R, respectively.
From the analysis above, it is known that the EA of the designed BHTS increases with the increase of the hollow column. The EA of BHTS–I–Specimens and BHTS–II–Specimens was larger than that of BHTS–S–Specimens, and the EA of BHTS–II–Specimens was slightly larger than that of BHTS–I–Specimens. In the SEA analysis, the SEA of triangle BHTS–I–T was slightly smaller than that of BHTS–S–T. The SEAs of BHTS–I–Specimens and BHTS–II–Specimens in rectangle and hexagon were larger than that of BHTS–S–Specimens, and the SEA of BHTS–II–Specimens was slightly larger than that of BHTS–Ⅰ–Specimens. This indicates that the increase of cross-sectional columns and the distribution position can raise the EA and SEA of the honeycomb structure.

4. Conclusions

Due to the rapid development of materials science, the requirements for lightweight and high mechanical properties of building materials have grown sharply. Bio-inspired honeycomb column thin-walled structures (BHTS) is a new bionic structure proposed in recent years, with excellent mechanical and energy absorption properties. In order to study the effect of section angle number on the energy absorption of the proposed BHTS and ensure the accurate fabrication of complex honeycomb structures, two BHTS manufactured using AM-SLM technology were proposed. Through the combination of experiment and numerical simulation, the initial peak crushing force, mean crushing force, energy absorption, and specific energy absorption of the proposed structure were compared. The conclusions of this study are as follows:
  • The bionic structure of the beetle elytra greatly enhances the energy absorption performance of the honeycomb structure. The design and installation location of BHTS hollow columns have a significant influence on the mechanical behavior of honeycomb structures.
  • The initial peak crushing force, mean crushing force, and energy absorption of BHTS were greatly improved with the increase of section angle number (triangle, rectangle, and hexagon). The hexagonal structure has the best energy absorption indexes among the three structures, as described above.
  • The finite element simulation results were in good agreement with the experimental results. This research method could provide an important reference for material designing, manufacturing, and modeling based on BHTS.
  • More cell elements and larger sizes should be added to the subsequent design to make the test results obtained by convergence more accurate.
In general, it is proven that the number of section angles has a great influence on the energy absorption capacity of BHTS manufactured by AM-SLM AlSi10Mg. Considering that increase in the section angles has a significant enhancement effect on energy absorption significantly, BHTS has great potential in various engineering applications.

Author Contributions

Conceptualization, H.X.; methodology, S.W. and H.X.; validation, S.W.; formal analysis, S.W.; investigation, S.W. and H.X.; resources, Y.L.; data curation, S.W.; writing—original draft preparation, S.W.; writing—review and editing, H.X.; visualization, S.W.; supervision, H.X.; project administration, H.X.; funding acquisition, H.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Fundamental Research Funds for the Central universities] grant number [2572019AB16] and the APC was funded by [H.X.].

Institutional Review Board Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baroutaji, A.; Sajjia, M.; Olabi, A.G. On the crashworthiness performance of thin-walled energy absorbers: Recent advances and future developments. Thin-Walled Struct. 2017, 118, 137–163. [Google Scholar] [CrossRef] [Green Version]
  2. Cwi, A.; Ce, B. A review of the crashworthiness performance of energy absorbing composite structure within the context of materials, manufacturing and maintenance for sustainability. Compos. Struct. 2020, 257, 113081. [Google Scholar]
  3. Zhang, Y.; He, N.; Song, X.; Chen, T.; Chen, H. On impacting mechanical behaviors of side fractal structures. Thin-Walled Struct. 2020, 146, 106490. [Google Scholar] [CrossRef]
  4. Gao, Q.; Liao, W.H.; Huang, C. Theoretical predictions of dynamic responses of cylindrical sandwich filled with auxetic structures under impact loading. Aerosp. Sci. Technol. 2020, 107, 106270. [Google Scholar] [CrossRef]
  5. Alkoles, O.M.S.; Mahdi, E.; Hamouda, A.M.S.; Sahari, B.B. Ellipticity Ratio Effects in the Energy Absorption of Axially Crushed Composite Tubes. Appl. Compos. Mater. 2003, 10, 339–363. [Google Scholar] [CrossRef]
  6. Czaplicki, M.J.; Robertson, R.E.; Thornton, P.H. Comparison of bevel and tulip triggered pultruded tubes for energy absorption. Compos. Sci. Technol. 1991, 40, 31–46. [Google Scholar] [CrossRef] [Green Version]
  7. Reuter, C.; Sauerland, K.-H.; Tröster, T. Experimental and numerical crushing analysis of circular CFRP tubes under axial impact loading. Compos. Struct. 2017, 174, 33–44. [Google Scholar] [CrossRef]
  8. Gao, Q.; Wang, L.; Wang, Y.; Wang, C. Crushing analysis and multiobjective crashworthiness optimization of foam-filled ellipse tubes under oblique impact loading. Thin-Walled Struct. 2016, 100, 105–112. [Google Scholar] [CrossRef]
  9. Zhang, X.W.; Su, H.; Yu, T.X. Energy Absorption of An Axially Crushed Square Tube With A Buckling Initiator. Int. J. Impact Eng. 2009, 36, 402–417. [Google Scholar] [CrossRef]
  10. Xing, B.F.; Hu, D.Y.; Sun, Y.X.; Yang, J.L.; Yu, T.X. Effects of hinges and deployment angle on the energy absorption characteristics of a single cell in a deployable energy absorber. Thin-Walled Struct. 2015, 94, 107–119. [Google Scholar] [CrossRef]
  11. Ali, M.; Ohioma, E.; Kraft, F.; Alam, K. Theoretical, numerical, and experimental study of dynamic axial crushing of thin walled pentagon and cross-shape tubes. Thin-Walled Struct. 2015, 94, 253–272. [Google Scholar] [CrossRef]
  12. Song, X.; Sun, G.; Li, G.; Gao, W.; Li, Q. Crashworthiness optimization of foam-filled tapered thin-walled structure using multiple surrogate models. Struct. Multidiscip. Optim. 2013, 47, 221–231. [Google Scholar] [CrossRef]
  13. Li, G.; Xu, F.; Sun, G.; Li, Q. A comparative study on thin-walled structures with functionally graded thickness (FGT) and tapered tubes withstanding oblique impact loading. Int. J. Impact Eng. 2015, 77, 68–83. [Google Scholar] [CrossRef]
  14. Wang, Z.; Liu, J.; Hui, D. Mechanical behaviors of inclined cell honeycomb structure subjected to compression. Compos. Part B Eng. 2016, 110, 307–314. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Liu, T.; Tizani, W. Experimental and numerical analysis of dynamic compressive response of Nomex honeycombs. Compos. Part B 2018, 148, 27–39. [Google Scholar] [CrossRef]
  16. Guangyong, S.; Tangying, L.; Jianguang, F.; Steve, P.G.; Qing, L. Configurational optimization of multi-cell topologies for multiple oblique loads. Struct. Multidiscip. Optim. 2018, 57, 469–488. [Google Scholar]
  17. Wang, Z.; Liu, J.; Yao, S. On folding mechanics of multi-cell thin-walled square tubes. Compos. Part B 2018, 132, 17–27. [Google Scholar] [CrossRef]
  18. Sun, G.; Pang, T.; Fang, J.; Li, G.; Li, Q. Parameterization of Criss-Cross Configurations for Multiobjective Crashworthiness Optimization. Int. J. Mech. Sci. 2017, 124–125, 145–157. [Google Scholar] [CrossRef]
  19. Wu, S.; Zheng, G.; Sun, G.; Liu, Q.; Li, G.; Li, Q. On design of multi-cell thin-wall structures for crashworthiness. Int. J. Impact Eng. 2016, 88, 102–117. [Google Scholar] [CrossRef]
  20. Hu, D.; Wang, Y.; Song, B.; Dang, L.; Zhang, Z. Energy-absorption characteristics of a bionic honeycomb tubular nested structure inspired by bamboo under axial crushing. Compos. Part B 2019, 162, 21–32. [Google Scholar] [CrossRef]
  21. Xiang, J.; Du, J. Energy absorption characteristics of bio-inspired honeycomb structure under axial impact loading. Mater. Sci. Eng. A 2017, 696, 283–289. [Google Scholar] [CrossRef]
  22. Hao, P.; Du, J. Energy absorption characteristics of bio-inspired honeycomb column thin-walled structure under impact loading. J. Mech. Behav. Biomed. Mater. 2018, 79, 301–308. [Google Scholar] [CrossRef] [PubMed]
  23. Gibson, I.; Rosen, D.W.; Stucker, B. Development of Additive Manufacturing Technology; Springer: New York, NY, USA, 2010. [Google Scholar]
  24. Dizon, J.; Espera, A.H.; Chen, Q.; Advincula, R.C. Mechanical characterization of 3D-printed polymers. Addit. Manuf. 2018, 20, 44–67. [Google Scholar] [CrossRef]
  25. Ahangar, P.; Cooke, M.E.; Weber, M.H.; Rosenzweig, D.H. Current Biomedical Applications of 3D Printing and Additive Manufacturing. Appl. Sci. 2019, 9, 1713. [Google Scholar] [CrossRef] [Green Version]
  26. Genovese, K.; Leeflang, S.; Zadpoor, A.A. Microscopic full-field three-dimensional strain measurement during the mechanical testing of additively manufactured porous biomaterials. J. Mech. Behav. Biomed. Mater. 2017, 69, 327–341. [Google Scholar] [CrossRef]
  27. Javaid, M.; Haleem, A. Additive manufacturing applications in medical cases: A literature based review. Alex. J. Med. 2017, 54, 411–422. [Google Scholar] [CrossRef] [Green Version]
  28. Xia, H.; Sun, Q.; Liu, Y. Energy Absorption Characteristics of Bio-Inspired Honeycomb Column Thin-Walled Structure under Low Strain Rate Uniaxial Compression Loading. Energies 2022, 15, 6957. [Google Scholar] [CrossRef]
  29. Ngoc, S.H.; Pham, T.P.; Tran, T.T.; Hong, H.; Guoxing, L. Mechanical properties and energy absorption of bio-inspired hierarchical circular honeycomb. Compos. Part B Eng. 2022, 236, 109818. [Google Scholar]
  30. Tiryakioglu, M. Si particle size and aspect ratio distributions in an Al-7%Si-0.6%Mg alloy during solution treatment. Mater. Sci. Eng. A 2008, 473, 1–6. [Google Scholar] [CrossRef]
  31. Laurencon, M.; Resseguier, T.D.; Loison, D.; Baillargeat, J.; Ngnekou, J.; Nadot, Y. Effects of additive manufacturing on the dynamic response of AlSi10Mg to laser shock loading. Mater. Sci. Eng. 2019, 748, 407–417. [Google Scholar] [CrossRef]
  32. Yin, H.; Huang, X.; Scarpa, F.; Wen, G.; Chen, Y.; Zhang, C. In-plane crashworthiness of bio-inspired hierarchical honeycombs. Compos. Struct. 2018, 192, 516–527. [Google Scholar] [CrossRef] [Green Version]
  33. Wan, Y.; Zhu, L.; Fang, H.; Liu, W.; Mao, Y. Experimental testing and numerical simulations of ship impact on axially loaded reinforced concrete piers. Int. J. Impact Eng. 2018, 125, 246–262. [Google Scholar] [CrossRef]
  34. Wang, H.; Yang, B.; Chen, K.; Elchalakani, M. Parametric analysis and simplified approach for steel-framed subassemblies with reverse channel connection under falling-debris impact. Eng. Struct. 2020, 225, 111263. [Google Scholar] [CrossRef]
  35. Gunes, R.; Arslan, K. Development of numerical realistic model for predicting low-velocity impact response of aluminium honeycomb sandwich structures. J. Sandw. Struct. Mater. 2016, 18, 1. [Google Scholar] [CrossRef]
  36. Xie, S.; Jing, K.; Zhou, H.; Liu, X. Mechanical properties of Nomex honeycomb sandwich panels under dynamic impact. Compos. Struct. 2020, 235, 111814. [Google Scholar] [CrossRef]
  37. Available online: http://lsdyna.ru/wp-content/uploads/2018/10/LS-DYNA-KEYWORD-USERS-MANUAL-VOL1.pdf (accessed on 31 October 2022).
  38. Ahmed, N.; Pu, X.; Zafar, N. Dynamic axial crushing of bitubular tubes with curvy polygonal inner-tube sections. Int. J. Comput. Mater. Sci. Eng. 2017, 06, 1750024. [Google Scholar] [CrossRef]
  39. Tran, T.N.; Ton, T. Lateral crushing behavior and theoretical prediction of thin-walled rectangular and square tubes. Compos. Struct. 2016, 154, 374–384. [Google Scholar] [CrossRef]
  40. Jhd, A.; Yjw, C.; Fnj, A.; Hlfb, C. Crushing behaviors of buckling-induced metallic meta-lattice structures. Def. Technol. 2021, 18, 10. [Google Scholar]
  41. Ryzińska, G.; David, M.; Prusty, G.; Tarasiuk, J.; Wroński, S. Effect of fibre architecture on the specific energy absorption in carbon epoxy composite tubes under progressive crushing. Compos. Struct. 2019, 227, 111292. [Google Scholar] [CrossRef]
  42. Hamid, W.L.H.W.A.; Aminanda, Y.; Dawood, M.S.I.S. Experimental Investigation on the Energy Absorption Capability of Foam-Filled Nomex Honeycomb Structure. Appl. Mech. Mater. 2013, 2622, 393. [Google Scholar]
  43. Philipp, B. CircStat: A MATLAB toolbox for circular statistics. J. Stat. Softw. 2009, 31, 1–21. [Google Scholar] [CrossRef] [Green Version]
  44. Sun, D.; Zhang, W.; Wei, Y. Mean out-of-plane dynamic plateau stresses of hexagonal honeycomb cores under impact loadings. Compos. Struct. 2010, 92, 2609–2621. [Google Scholar]
  45. Jin, L.; Yu, W.; Li, D.; Du, X. Numerical and theoretical investigation on the size effect of concrete compressive strength considering the maximum aggregate size. Int. J. Mech. Sci. 2021, 192, 106130. [Google Scholar] [CrossRef]
Figure 1. BHTS and its application: (a) the beetle Allomyrina dichotoma; (b) micromorphology of the honeycomb structure with columns [22] Copyright©5430520399371.
Figure 1. BHTS and its application: (a) the beetle Allomyrina dichotoma; (b) micromorphology of the honeycomb structure with columns [22] Copyright©5430520399371.
Biomimetics 07 00201 g001
Figure 2. Specimen size and loading: (a) specimen size; (b) loading process.
Figure 2. Specimen size and loading: (a) specimen size; (b) loading process.
Biomimetics 07 00201 g002
Figure 3. BHTS models.
Figure 3. BHTS models.
Biomimetics 07 00201 g003
Figure 4. BHTS loading.
Figure 4. BHTS loading.
Biomimetics 07 00201 g004
Figure 5. FEM: (a) boundary setting of the FEM; (b) comparison of the results.
Figure 5. FEM: (a) boundary setting of the FEM; (b) comparison of the results.
Biomimetics 07 00201 g005
Figure 6. Schematic diagram of the compression process: (a) Compression start; (b) Slight compression; (c) Severe compression.
Figure 6. Schematic diagram of the compression process: (a) Compression start; (b) Slight compression; (c) Severe compression.
Biomimetics 07 00201 g006aBiomimetics 07 00201 g006b
Figure 7. Site layout diagram.
Figure 7. Site layout diagram.
Biomimetics 07 00201 g007
Figure 8. The crushing process of BHTS specimens under quasi-static uniaxial compression loading.
Figure 8. The crushing process of BHTS specimens under quasi-static uniaxial compression loading.
Biomimetics 07 00201 g008
Figure 9. Force-displacement diagram of each structure under quasi-static uniaxial compression loading: (ac) BHTS–S; (df) BHTS–Ⅰ; (gi) BHTS–Ⅱ.
Figure 9. Force-displacement diagram of each structure under quasi-static uniaxial compression loading: (ac) BHTS–S; (df) BHTS–Ⅰ; (gi) BHTS–Ⅱ.
Biomimetics 07 00201 g009aBiomimetics 07 00201 g009b
Figure 10. Statistics of PCF and MCF corresponding to different section angles: (a) PCF; (b) MCF.
Figure 10. Statistics of PCF and MCF corresponding to different section angles: (a) PCF; (b) MCF.
Biomimetics 07 00201 g010
Figure 11. EA evaluation indexes of BHTS at different section angle number: (a,c,e) EA; (b,d,f) SEA.
Figure 11. EA evaluation indexes of BHTS at different section angle number: (a,c,e) EA; (b,d,f) SEA.
Biomimetics 07 00201 g011
Figure 12. EA evaluation indexes of BHTS at different hollow column positions: (a,c,e) EA; (b,d,f) SEA.
Figure 12. EA evaluation indexes of BHTS at different hollow column positions: (a,c,e) EA; (b,d,f) SEA.
Biomimetics 07 00201 g012
Table 1. BHTS parameters of the study.
Table 1. BHTS parameters of the study.
Specimen LabelCell Side Length
L c   ( mm )
Cell Wall Thickness T c (mm) Column Diameter
D c (mm)
Column Height
H c (mm)
BHTS–S–T6.00.5/20.0
BHTS–Ⅰ–T6.00.50.520.0
BHTS–Ⅱ–T6.00.50.520.0
BHTS–S–R6.00.5/20.0
BHTS–Ⅰ–R6.00.50.520.0
BHTS–Ⅱ–R6.00.50.520.0
BHTS–S–H6.00.5/20.0
BHTS–Ⅰ–H6.00.50.520.0
BHTS–Ⅱ–H6.00.50.520.0
Table 2. Properties of Si-Al alloy AlSi10Mg.
Table 2. Properties of Si-Al alloy AlSi10Mg.
DensityYoung’s ModulusPoisson’s RatioInitial Yield Strength
2670 kg/m369 ± 5 GPa0.3220 ± 10 MPa
Table 3. BHTS test set number.
Table 3. BHTS test set number.
TypeGroup Number|Panel
Specimen1BHTS–S–T14BHTS–S–R17BHTS–S–H1
FEMFEM–BHTS–S–TFEM-BHTS–S–RFEM-BHTS–S–H
Specimen2BHTS–Ⅰ–T15BHTS–Ⅰ–R18BHTS–Ⅰ–H1
FEMFEM-BHTS–Ⅰ–TFEM-BHTS–Ⅰ–RFEM-BHTS–Ⅰ–H
Specimen3BHTS–Ⅱ–T16BHTS–Ⅱ–R19BHTS–Ⅱ–H1
FEMFEM–BHTS–Ⅱ–TFEM-BHTS–Ⅱ–RFEM-BHTS–Ⅱ–H
Table 4. Evaluation index values of mechanical properties of BHTS under different strain rates.
Table 4. Evaluation index values of mechanical properties of BHTS under different strain rates.
GroupPanelSpecimenMassPCFMCFCFEEASEAAverage Value
Number PCFMCFCFEEASEA
(g)(kN)(kN)(%)(J)(J/g)(kN)(kN)(%)(J)(J/g)
11BHTS–S–T1.381.662.21 133.29 25.2818.321.642.16131.8724.9617.90
21.411.612.10 130.44 24.6517.48
21BHTS–Ⅰ–T1.787.912.43 30.64 44.4124.957.822.4130.8444.2625.08
21.757.722.39 31.04 44.1025.20
31BHTS–II–T1.605.042.08 41.19 34.4321.525.142.0339.5433.7021.00
21.615.231.98 37.88 32.9720.48
41BHTS–S–R1.892.861.79 62.50 27.0114.292.881.7961.9626.7814.21
21.882.891.78 61.41 26.5514.12
51BHTS–Ⅰ–R2.439.962.91 29.30 52.5021.6010.132.9128.7653.3421.90
22.4410.302.91 28.25 54.1722.20
61BHTS–II–R2.227.303.58 48.93 55.9325.197.213.5248.7255.6525.30
22.187.123.45 48.50 55.3725.40
71BHTS–S–H3.023.203.1397.81 44.2914.673.183.1197.9644.1314.74
22.973.153.0998.1043.9614.80
81BHTS–Ⅰ–H3.8515.304.19 27.39 77.8920.2315.284.1427.1076.6319.83
23.8815.264.09 26.81 75.3719.43
91BHTS–II–H3.6312.134.19 34.50 71.8419.7912.194.2434.7873.1720.08
23.6612.244.29 35.06 74.5020.36
Note: The number of the panel represents the first or second valid test data.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, S.; Xia, H.; Liu, Y. Energy Absorption Characteristics of Polygonal Bio-Inspired Honeycomb Column Thin-Walled Structure under Quasi-Static Uniaxial Compression Loading. Biomimetics 2022, 7, 201. https://doi.org/10.3390/biomimetics7040201

AMA Style

Wang S, Xia H, Liu Y. Energy Absorption Characteristics of Polygonal Bio-Inspired Honeycomb Column Thin-Walled Structure under Quasi-Static Uniaxial Compression Loading. Biomimetics. 2022; 7(4):201. https://doi.org/10.3390/biomimetics7040201

Chicago/Turabian Style

Wang, Shijie, Hongxiang Xia, and Yancheng Liu. 2022. "Energy Absorption Characteristics of Polygonal Bio-Inspired Honeycomb Column Thin-Walled Structure under Quasi-Static Uniaxial Compression Loading" Biomimetics 7, no. 4: 201. https://doi.org/10.3390/biomimetics7040201

APA Style

Wang, S., Xia, H., & Liu, Y. (2022). Energy Absorption Characteristics of Polygonal Bio-Inspired Honeycomb Column Thin-Walled Structure under Quasi-Static Uniaxial Compression Loading. Biomimetics, 7(4), 201. https://doi.org/10.3390/biomimetics7040201

Article Metrics

Back to TopTop