Next Article in Journal
Prediction of Strength Properties of Concrete Containing Waste Marble Aggregate and Stone Dust—Modeling and Optimization Using RSM
Next Article in Special Issue
Low-Temperature Sintering of Ag Composite Pastes with Different Metal Organic Decomposition Additions
Previous Article in Journal
Microstructures and Mechanical Properties of an AlCoCrNiFe HEA/WC Reinforcing Particle Composite Coating Prepared by Laser Cladding
Previous Article in Special Issue
Wet Etching of Quartz Using a Solution Based on Organic Solvents and Anhydrous Hydrofluoric Acid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Re-Entrant Honeycomb Auxetic Structure with Enhanced Directional Properties

1
College of Electrical & Mechanical Engineering, National University of Sciences and Technology, Rawalpindi 43701, Pakistan
2
Department of Mechanical Engineering, Faculty of Engineering, Islamic University of Madinah, Madinah 42351, Saudi Arabia
3
Department of Mechanical Engineering, Faculty of Engineering Science and Technology, Hamdard University, Karachi 74600, Pakistan
4
Department of Chemical Engineering, Faculty of Engineering, Islamic University of Madinah, Madinah 42351, Saudi Arabia
*
Author to whom correspondence should be addressed.
Materials 2022, 15(22), 8022; https://doi.org/10.3390/ma15228022
Submission received: 28 September 2022 / Revised: 20 October 2022 / Accepted: 8 November 2022 / Published: 14 November 2022
(This article belongs to the Special Issue Smart Materials and Structures: Characteristics and Applications)

Abstract

:
This paper presents a modified re-entrant honeycomb auxetic structure. The structure is constructed by adding an additional horizontal member between the vertical and re-entrant member of the semi-re-entrant honeycomb model to increase the overall compliance of the structure in order to obtain higher values of negative Poisson’s ratio (NPR). An analytical model of the structure is presented, taking into account the bending, shear, and axial deformations. The model is verified using finite element analysis (FEA) and tensile testing. The results of FEA and tensile testing corroborate the results of the presented mathematical model. The structure is also compared to the existing re-entrant honeycomb structure. The newly added strut has shown a direct effect on the directional properties of the overall structure. With an increase in the newly added strut to re-entrant lengths, NPR was significantly enhanced in the x-direction and reduced in the y-direction loadings. The structure shows an improved Young’s modulus compared to solid material in both loading directions, especially for the low values of the new strut and re-entrant lengths ratio. The structure also shows that high NPR can be achieved for low relative density compared to semi re-entrant honeycomb structure.

1. Introduction

One of the most important properties of a material is its Poisson’s ratio. Most materials possess a positive Poisson’s ratio, meaning that when the material is stretched in one direction, its cross-section decreases in the transverse direction. However, a particular class of materials exhibits a negative Poisson’s ratio (NPR), meaning that their cross-section will increase in the transverse direction when a tensile load is applied. Such materials/structures are also commonly known as auxetic or NPR materials/structures. Auxetic materials provide increased shear modulus, fracture toughness, in-plain strain fracture resistance, energy absorption, indentation resistance, and acoustic response [1,2,3,4,5,6].
A study by Donoghue et al. [7] found that auxetic composite laminates require more energy for a crack to propagate than a conventional laminate under tensile loading, indicating that auxetic composites provide better fracture toughness than conventional composites. Another study by Choi and Lakes [8] compared the fracture stiffness of auxetic foams with conventional foams and observed that auxetic foams provide greater fracture toughness than conventional ones even with a smaller Young’s modulus. Yao et al. [9] proposed and tested an anti-pullout auxetic pedicle screw made using the re-entrant honeycomb structure and found that the greater auxeticity of the screw results in a better anti-pullout performance, which they proposed can resolve the loosening and pulling out issues of typical pedicle screws currently being used. Several studies were found in the literature on the behavior of auxetic structures under impact loading. A study by Zhang et al. [4] showed that sandwich panels with re-entrant honeycomb cores absorbed explosive loads better than the honeycomb structure.
There have been numerous studies on different types of auxetic unit cells in the literature, especially in the past decade, due to the ease of manufacturing such structures because of advancements in additive manufacturing technologies. The traditional periodic cellular structure was proposed by Gibson et al. [10] in the form of a 2D re-entrant honeycomb. The structure is made by inverting the direction of inclined members in a typical honeycomb structure. The resultant structure provides NPR in both the x- and y-directions. A few years later, Gibson et al. [10] provided a mathematical model for the 2D re-entrant model to calculate its Poisson’s ratio and Young’s modulus based on the flexural method. Subsequently, Masters et al. [11] improved on this by developing a model for these elastic constants with the inclusion of stretching, hinging, and flexure. The advent of additive manufacturing gave way to the manufacturing of three-dimensional versions of the auxetic structures, which were difficult to manufacture using conventional means; one such structure was proposed by Lang et al. [12], who provided a model for a 3D version of the re-entrant honeycomb model and compared it with FEA and experimental testing results.
The re-entrant honeycomb model is the most widely researched auxetic structure [13,14,15,16,17,18]. Re-entrant is used to describe something directed inwards or that has a negative angle. The idea of an auxetic structure was first proposed by Gibson et al. [10] in 1982 in the form of a 2D re-entrant honeycomb by introducing a negative angle in the inclined member of a typical honeycomb structure. Recently, Al-Rifaie et al. [19] presented their patented idea of using a re-entrant honeycomb structure manufactured using aluminum sheets by bending in different orientations, and the structure showed good damping properties.
Multiple structure variations have been proposed to enhance or modify the mechanical properties of this structure. These structures are made by either adding additional members to the existing structure or combining it with other types of auxetic structures. One such structure was proposed by Sigmund et al. [19], where they used topology optimization techniques to devise a variation of the re-entrant honeycomb model. Fu et al. [20] proposed a structure created by the additional embedding of a rhombic configuration into a re-entrant honeycomb model. The structure results in enhanced buckling strength and in−plane stiffness. Another such structure was proposed by Harkati et al. [21] by adding a chamfer to the edges of the diagonal member. An analytical model of the structure was presented, supported by finite element analysis (FEA) results. Adding the chamfer increases NPR behavior. Bezazi et al. [22] presented a mathematical model for a novel central symmetric structure by adding a chamfer to the existing honeycomb structure. The authors only considered bending deflection in the model, calculated the strain in the x-direction caused by loading in the respective direction, did the same for the y-direction, and divided them both to obtain the Poisson’s ratio, because of the nature of the central symmetric model. This approach may not be applicable for other unit cells to calculate values of Poisson’s ratio, especially in the case of non−isotropic structures. Zhu et al. [23] introduced a modified version of the re-entrant structure by introducing zigzag inclined members in the re-entrant honeycomb structure resulting in an increase in stiffness and auxeticity in one of the principal directions. Wu et al. [24] combined the re-entrant honeycomb and regular honeycomb structure to create a hybrid structure that can provide positive, negative, and zero Poisson’s ratio depending on its configuration. Zhang et al. [25] studied the energy absorption properties of aluminum composite tubes with hexagonal honeycomb and re-entrant honeycomb fillings and found that the one with the auxetic filling exhibits superior ductility compared to its non-auxetic counterpart. Nguyen et al. [6] analyzed a sandwich panel with a re-entrant honeycomb core and graphene nanoplates as the skin, and such a structure can benefit from the energy absorption and lightweight nature of the structure.
A periodic cellular structure is proposed in this paper by adding an extra member to the existing re-entrant honeycomb structure.
An additional strut is introduced in the existing re-entrant honeycomb structure to enhance the auxetic behavior of the structure. A complete analytical model is introduced to cater to in−plane bending, shear, and axial deformation. Variation in the newly introduced strut length is analyzed analytically, using FEA, and experimentally. A detailed analysis of the proposed structure is conducted by varying the unit cell parameters through the analytical equations and FEA.

2. Unit Cell Design

The proposed unit cell is shown in Figure 1. Horizontal members have been added to both sides between the re-entrant part and the vertical members of the re-entrant honeycomb structure. The newly added members (shown in red) will increase the compliance of the structure, which should result in increased deflections, especially in the y-direction, which should increase the NPR of the structure. The length of these members is the same and is denoted by a. The structure can be defined entirely using six perimeters which are: the length of the inclined member L, the length of the vertical member H, the length of the newly introduced member a, the thickness of the struts t, the depth of the structure in out-of-plane direction b, and the angle made by the inclined member with the horizontal axis θ.
The total length of the structure in the x- and y-direction can be calculated as follows:
l x = 2 2 a + L c o s θ
l y = 2 H L s i n θ
Here l x is the length in the x-direction and l y is the length in the y-direction. If we are to fit the structure in a cuboid, the volume V s that the structure will occupy can be calculated as:
V s = b l x l y = 4 b 2 a + L c o s θ H L c o s θ
The actual volume V a occupied by the structure can be calculated by adding the individual volumes of every member, which comes out to be:
V a = 2 b t H + 4 a + 2 L
If we consider the structure as it is or model it as a cuboid with the original structure’s properties, the model’s mass m will remain the same. The relative density is defined as the ratio of the cuboid model’s density with the model’s actual density. Hence the relative density of the structure can be calculated as:
ρ s ρ a = m V s × V a m = t H + 4 a + 2 L 2 2 a + L c o s θ H L s i n θ
The parameters can be normalized by taking the ratios with a fixed parameter. This makes an analysis of the structure much more straightforward. For our case, the fixed parameter is the length of the inclined member L. After inserting the normalized parameters, Equation (5) becomes.
ρ s ρ a = t L H L + 4 a L + 2 2 2 a L + c o s θ H L s i n θ
Because we are adding an additional member to an existing structure, making the value of a equal to zero yields the relative density of the re-entrant honeycomb model.

3. Analytical Model

Simple mechanics of materials techniques can be used to calculate the elastic parameters of the structure. By looking at the structure in Figure 1, it can be seen that the structure is symmetric along the X and Y axes. Hence the elastic parameters can be calculated using a quarter model since the strain of the structure will remain the same in both instances because of the differences in lengths. The Poisson’s ratio and Young’s modulus in the x-direction are denoted by v x and E x , respectively, and for the y-direction by v y and E y . The shear modulus is denoted by G . It is assumed that the structure is located in the middle of the lattice to avoid boundary effects. Three types of deformations are considered for the model; these include flexural/bending deformations, shear deformations, and axial deformations. The structure is first loaded in both principal axes, one at a time, to calculate the elastic parameters in the respective direction, then a shear load is applied to calculate the shear modulus of the structure.

3.1. X-Direction

The free body diagram (FBD) of the structure under loading in the x-direction is given in Figure 2.
The bending deflection of the re-entrant member consists of the bending caused by the force P and the deflection caused as a result of the moment induced in the horizontal member, and it can be calculated as [26]:
δ f X = P L 3 s i n θ 12 E s I 6 a L + 1
The bending deflection caused due to moment M on member 2–3 is:
δ f a X = P L a 2 s i n θ 2 E s I
The shear deflection δ s i x caused in the inclined member can be calculated by [10]:
δ s i X = P L 3 s i n θ 12 E s I 2.4 + 1.5 v s t L 2
The axial deformation of the inclined member can be calculated as follows [27]:
δ a i X = P L c o s θ E s b t
Similarly, the deformations for the horizontal members can be calculated by [27]:
δ a h X = 2 P a E s b t
The total deflection ( δ x x ) of the structure in the x-direction equals:
δ x x = δ f x s i n θ + δ s i x s i n θ + δ a i x c o s θ + δ a h x
After simplifying, the equation becomes
δ x x = P L 3 s i n θ 2 12 E s I 1 + 6 a L + 2.4 + 1.5 v s + 1 t a n θ 2 + 1 s i n θ 2 a L t L 2
The total deflection ( δ x y ) of the structure in the y-direction equals:
δ x y = δ f x c o s θ + δ s i x c o s θ + δ f a X
δ x y = P L 3 s i n θ c o s θ 12 E s I 1 + 6 a L + 2.4 + 1.5 v s t L 2 + 1 c o s θ a L 2
Strain is defined as the ratio of the deflections of a material to its un-deformed length. Hence the strain ϵ x x in the x-direction becomes [28]:
ϵ x x = δ x x 2 a + L c o s θ
ϵ x x = P L 2 s i n θ 2 12 E s I 2 a L + c o s θ 1 + 6 a L + 2.4 + 1.5 v s + 1 t a n θ 2 + 1 s i n θ 2 a L t L 2
Similarly, the strain ϵ x y in the y-direction [28]:
ϵ x y = δ x y H L s i n θ
ϵ x y = P L 2 s i n θ c o s θ 12 E s I H L s i n θ 1 + 6 a L + 2.4 + 1.5 v s t L 2 + 1 c o s θ a L 2
The Poisson’s ratio is the ratio of lateral direction strain to the strain in the direction of the applied force. Hence the Poisson’s ratio for loading in the x-direction equals [28]:
v x = ϵ x y ϵ x x
v x = 2 a L + c o s θ c o s θ H L s i n θ s i n θ 1 + 6 a L + 2.4 + 1.5 v s t L 2 + 1 c o s θ a L 2 1 + 6 a L + 2.4 + 1.5 v s + 1 t a n θ 2 + 1 s i n θ 2 a L t L 2
The Young’s modulus for loading in the x-direction is defined as [28]:
E x = σ x x ε y y
E x = E s t L 3 2 a L + c o s θ H L s i n θ s i n 2 θ 1 1 + 6 a L + 2.4 + 1.5 v s + 1 t a n θ 2 + 1 s i n θ 2 a L t L 2

3.2. Y-Direction

The FBD of a quarter structure under loading in the y-direction is given in Figure 3.
The bending deflection caused in the re-entrant member can be calculated as [28]
δ r Y = W L 3 c o s θ 12 E s I
The deflection caused as a result of moment My2-3 equals [28]
δ r m Y = W a L 2 c o s θ 2 E s I
The deflection of the horizontal member 4–5 is
δ h m Y = W a 3 2 E s I + W a 2 L c o s θ 2 E s I
The shear deflection of the structure can be calculated as [10]:
δ s r Y = W L 3 c o s θ 12 E s I 2.4 + 1.5 v s t L 2
Moreover, the axial deformation of the re-entrant member can be calculated by [28]
δ a r Y = W L s i n θ E s b t
Similarly, the axial deformation in the vertical members can be calculated by [28]
δ a v Y = W H E s b t
The total deflection of the horizontal members can be calculated as follows [28]:
δ h Y = W a 3 6 E s I
The total deformation in the y-direction can be calculated by
δ y y = δ r Y c o s θ + δ s r Y c o s θ + δ a r Y s i n θ + δ r m Y c o s θ + δ h m Y + δ a v Y
After simplification, we get
δ y y = P L 3 c o s 2 θ 12 E s I 1 + 6 a L + 6 a L 2 1 c o s θ + 8 a L 3 1 c o s 2 θ + 2.4 + 1.5 v s + 1 tan 2 θ + H L 1 c o s 2 θ t L 2
Similarly, for deflections in the x-direction, we get
δ y x = P L 3 s i n θ c o s θ 12 E s I 1 + 6 a L + 2 a L 2 1 c o s θ + 2 a L 3 1 s i n θ c o s θ + 2.4 + 1.5 v s t L 2
The strain ϵ y y in the y-direction equals
ϵ y y = δ y y H L s i n θ
ϵ y y = F 2 L 3 c o s θ 2 12 E s I H L s i n θ 1 + 6 a L + 6 a L 2 1 c o s θ + 8 a L 3 1 c o s 2 θ + 2.4 + 1.5 v s + 1 tan 2 θ + H L 1 c o s 2 θ t L 2
The strain ϵ y x in the x-direction is
ϵ y x = δ y x 2 a + L c o s θ
ϵ y x = F 2 L 3 c o s θ s i n θ 12 E s I 2 a + L c o s θ 1 + 6 a L + 2 a L 2 1 c o s θ + 2 a L 3 1 s i n θ c o s θ + 3.4 + 1.5 v s t L 2
Using the ratio of the two strains, the Poisson’s ratio v y for loading in the y-direction in normalized form can be calculated as follows:
v y = ϵ y x ϵ x x
v y = H L s i n θ s i n θ 2 a L + c o s θ c o s θ 1 + 6 a L + 2 a L 2 1 c o s θ + 2 a L 3 1 s i n θ c o s θ + 3.4 + 1.5 v s t L 2 1 + 6 a L + 6 a L 2 1 c o s θ + 8 a L 3 1 c o s 2 θ + 2.4 + 1.5 v s + 1 tan 2 θ + H L 1 c o s 2 θ t L 2
Young’s modulus E y can be calculated using the stress σ y y and strain ϵ y y in the y-direction:
E y = σ y y ϵ y y
E y = E s   t L 3 H L s i n θ 2 a L + c o s θ c o s 2 θ 1 1 + 6 a L + 6 a L 2 1 c o s θ + 8 a L 3 1 c o s 2 θ + 2.4 + 1.5 v s + 1 tan 2 θ + H L 1 c o s 2 θ t L 2

4. Materials and Methods

4.1. Finite Element Analysis (FEA)

For initial analysis, two commercial software were used (Solidworks simulation suite and Abaqus) to carry out the FEA of the structure. There was no variation in results for both software. For ease of modeling and the faster ability to perform multiple design studies, the Solidworks simulation suite was chosen for the FEA. A linear static study was performed to calculate the deflections of the structure. The shell-type element was used. A mesh size of 0.5 mm was used for all the studies conducted. For the FEA of such cellular structures, to avoid edge and size effects during deformation, it is recommended to use a minimum of four unit cells [12]. In this study, a 5 × 5-unit cell array was considered, and the deflections were measured at the center cell of the lattice. The center unit cell and boundary conditions used in the FEA are shown in Figure 4. The center unit cell was drawn in different colors and mentioned points to measure the deflection. The difference of these points divided by the length of the cell in the respective direction gave the strain in that direction. The ratio of the two strains provides us with Poisson’s ratio of the structure. The structure was set to be fixed in both directions at one end, and a load was applied at the other end, replicating the boundary conditions of a tensile test.

4.2. Samples Preparation and Experimental Procedures

A total of 10 test samples were manufactured from 3.75 mm acrylic sheets. Thickness greater than that resulted in melting of the edges and sharp corners of the structure with the laser cutting. Half of the samples were for testing the properties in the x-direction, while the other half were for the y-direction. The parameters of the samples are shown in Table 1. All the samples were laser cut from the same sheet of acrylic to achieve the same material properties for each sample. In addition, five dog-bone-shaped samples were manufactured from the same sheet using the dimensions provided by the standard ISO 527-1 used for calculating the material properties of plastics. The laser-cut samples and the hinged support are shown in Figure 5.
In order to study the effect of the horizontal member, all the variables were fixed, except the length of the newly introduced strut a . The length is varied with 0.8 mm increments resulting in a / L values that range from 0.2 mm to 0.5 mm with 0.1 mm increments. Each sample was connected to a rectangular portion 120 mm in length and 20 mm in width. This portion contains six holes for 8 mm bolts to pass through. The samples with a = 0.8 mm ( a / L = 0.1) showed clear manufacturing defects as it was too close to the edges of the connecting struts and difficult to manufacture with the existing setup. Thus, they were excluded, and experimental results were discarded for those samples.
These 2D structures had to be tested with a tensile load, as the compressive load would result in buckling and out-of-plane deformation. A customized mild steel mount was designed to apply a uniform tension load to the structure, as shown in Figure 6. The mount was designed as a two-piece assembly that can hold the variable thickness of the sample for future work.
Tensile testing of the samples was performed using an Autograph AGS-J tabletop Universal Testing Machine (UTS). Tensile testing of the dog-bone samples was carried out at different extension rates according to ISO 527-1 to determine the Young’s Modulus of the material. The results of these tests can be seen in Figure 7. The strain rate is sensitive for polymers since the molecular mobility of the polymer chains becomes stiffer with the increment in strain rate [29]. The deflection rate for the modified re-entrant structure specimen was set at 10 mm/min, as this value shows less sensitivity to materials behavior.
The black color specimens were sprayed with a white dotted pattern to determine the deflections through image processing during the tensile tests. The experiment was filmed at 1080p using a DSLR camera at 30 frames per second. The deflection of the structure was analyzed by tracking the movement of each individual dot using a program, GOM Correlate [30].

5. Results and Discussion

The force vs. displacement results obtained for the structure under loading in the x- and y-directions can be seen in Figure 8. As we increase the values of a/L, the structure becomes more compliant at the expense of failure force.
In the next section, a detailed comparison of the values of Poisson’s ratios in both principal directions is given. However, looking at the test data from Figure 8 and Figure 9, we can see that the area under the graph increases by increasing the values of a/L, suggesting that the structure can absorb more energy compared to the re-entrant honeycomb model. Figure 10 compares the deflections in both the x- and y-directions obtained by FEA and testing results when overlapped with the image processing software. The behavior of the structure was almost the same, exhibiting the auxetic nature of the structure, and the relative deflection pattern in both cases was similar.

5.1. Effect on NPR with Different UC Parameters

5.1.1. Loading in the X-Direction

The comparison between the values of NPR obtained using the mathematical model, FEA, and the experimental testing is shown in Table 2 and Figure 11. The NPR obtained by all three methods exhibited close values and followed the same trend with approximately constant error margins. It is interesting to note that the modified structure’s auxetic behavior increased with the added strut value. When the a/L ratio was increased from 0.2 to 0.5, the tensile test showed that the NPR enhanced by 127%, which is a significant improvement. The additional strut provides greater flexibility to the re-entrant strut to enhance the Poisson’s ratio in the x-direction.
Once the FEA model was validated with the experimental results, a more detailed comparison between the mathematical model and FEA was conducted by varying the other parameters of the unit cell. In Figure 12, the comparison is performed by keeping the angle θ constant at 60° while the value of H/L is varied. Even in the extended range, the trends of both calculations followed the same pattern. Furthermore, it can be observed that by reducing the value of H/L, the value of NPR increased. Reducing the value of H/L means a reduction in the mass or relative density of the structure, which can be favorable in certain scenarios. Setting the value of a/L equal to zero results in the structure becoming a standard re-entrant honeycomb structure. Comparing the proposed structure with the re-entrant honeycomb model, we can see a significant increase in auxeticity for varying values of H/L.
On the other hand, in Figure 13, H/L is kept constant, and the angle θ is varied for an extended range of a/L. The results suggested that decreasing the re-entrant angle resulted in an enhancement of NPR. The mathematical model started deviating from the FEA when approaching higher re-entrant angles, but the trend remained the same, and the variation remained low. Previously significant gains in NPR can be seen compared to the re-entrant honeycomb model across the board.

5.1.2. Loading in the Y-Direction

The comparison between the values of NPR obtained using the mathematical model, FEA, and the experimental testing is shown in Table 3 and Figure 14. A sharp contrast can be observed here when compared to the NPR values when loading was in the x-direction. The increment in a/L showed a counter impact on NPR. Interestingly, the absolute values of NPR when the load was applied in the y-direction were high. During the tensile test, when a/L was 0.2, the NPR when loading in the x- and y-directions was −0.51 and −1.41, respectively. Therefore, high NPR can be achieved if the structure is oriented to apply a unidirectional load in the y-direction. However, when the a / L value was 0.4, the NPR value was approximately the same during tensile testing (NPRs when loading in the x- and y-directions were −0.89 and −0.9, respectively). This was the limiting ratio for a/L; after that, when the a/L equal to 0.5 was analyzed during tensile testing, the NPR was higher for the x-direction rather than in y-direction loading (NPRs when loading in the x- and y-directions were −1.16 and −0.72, respectively). These results are intriguing, as the orthotropic unit cell can exhibit isotropic behavior at certain unit cell parameters. Here, in this case, it was θ constant at 60 degrees and a/L ≃ 0.4.
Figure 15 and Figure 16 give a detailed comparison of NPR between the mathematical model and FEA for different H/L and re-entrant angles. The results gave an insight into the effects of these variables on NPR. As observed before, the trends exhibited by the mathematical model corroborate the trend obtained from FEA. In both cases, it was observed that NPR decreased as the length of the horizontal member increased. The NPR lost auxeticity rapidly for smaller values of a/L but started leveling out at the extreme values of a/L.
Contrary to x-direction loading, here it can be observed that the NPR increased with the value of H/L (see the arrow direction in Figure 15). Similarly, it can be noticed that NPR increases as the re-entrant angle decreases (see the arrow direction in Figure 16). Comparing the structure with the re-entrant honeycomb model, we can see that, in this specific case, better NPR can be achieved using the re-entrant honeycomb model.

5.1.3. Comparison of NPR

The combined effect on NPR calculated using the mathematical model for the variation of H/L, a/L, and re-entrant angle for both loading directions is shown in Figure 17. In general, when the loading was applied in the x-direction, the NPR was improved with low values of H/L, re-entrant angle, and higher values of a/L. This was because the added strut a is in the x-direction, and the lower value was favorable for the re-entrant strut to improve NPR. In contrast to the x-direction loading, the NPR showed mixed behavior when the loading was applied in the y-direction. NPR improved with either low values of H/L and high values of a/L or high values of re-entrant angle and low value of a/L.

5.2. Young’s Modulus of the Structure

The ratios of Young’s modulus for respective loading direction and the solid material of the structure, calculated by the mathematical model, are shown in Figure 18 (θ = 60°) and Figure 19 (H/L is 2.5). There was a significant improvement in NPR in both loading directions. The unit cell improved the strength of the material by adding the strut and making it stiffer. For x-direction loading, increasing the length of additional strut a or ratio a/L did not show a considerable variation in enhancement in the strength. However, lower values of a/L showed a stiffer behavior of the structure. The additional strut a is parallel along the x-direction, and thus the bending effect was low when the load was applied in the same direction. As a result, the variation in Ex/Es was also low.
On the other hand, the strut is perpendicular to the y-direction and can be subjected to extensive bending during the tensile load. Therefore, significant variation was observed with the increasing a/L ratio or varying re-entrant angle. However, the strength in the y-direction is significantly higher than in the x-direction. In fact, the structure showed stiffness in the y-direction for smaller values of a/L, which diminished with the a/L ratio being approximately more than 1.25. Nonetheless, a tradeoff for the strength can be made in both directions depending on the application of the structure.

5.3. Effect of Relative Density on NPR

The relative density and NPR of the structure for different values of a/L during loading in both directions (H/L is 2.5 and θ = 60°), calculated from the mathematical model, is shown in Figure 20. The relative density is the ratio of the volumes of the structure to the total volume it occupies if it is solid. It can be seen that with the increasing length of additional strut a, the relative density decreases while the NPR significantly increases in x-direction loading. There is a reverse effect on y-direction loading, but the variation is not high. It can be concluded that for a given thickness of the structure inside a fixed volume, higher NPR is possible while decreasing the overall mass of the structure for loading in the x-direction only. For comparison, the relative density and NPR for both directional loading for a standard re-entrant honeycomb structure are also shown in Figure 20. It can be observed that added strut highly improved the relative density and significantly improved NPR in x-direction loading.

6. Conclusions

The modified re-entrant honeycomb auxetic structure with enhanced directional properties is proposed and analyzed through a mathematical model, FEA, and experimental test. The modification was made by adding a strut connecting the vertical and re-entrant members of the re-entrant honeycomb unit cell. A detailed analytical model has been derived that considers bending, shear, and axial deformations to calculate Poisson’s ratio and other properties. Tensile tests were performed on laser−cut samples to validate the mathematical and FEA models. All three methods showed close agreement with the proposed structure’s negative Poisson’s ratio (NPR). The structure was analyzed by applying the loading in two principal directions on the doubly symmetrical unit cell. When the loading was applied in the x-direction, the NPR was improved with lower values of H/L and re-entrant angle, and higher values of a/L. On the other hand, when the loading is applied in the y-direction, the NPR showed mixed behavior. NPR improved with a combination of either low values of H/L and high values of a/L or high values of re-entrant angle and low values of a/L.
For a range of values of ‘a/L’, the relative density of the structure decreases with the increasing value of ‘a’. A lower relative density suggested that the structure utilized the occupied volume more effectively, meaning that it could provide a better NPR-to-mass ratio. To relate the relative density to NPR, the value of ‘a’ should be increased to obtain better relative density and an increase in Poisson’s ratio in the x-direction at the cost of a decrease in Poisson’s ratio in the y-direction. To improve its properties, the structure can be further modified by increasing the thickness of struts or adding cross-link members, making it a hybrid structure at the cost of losing NPR.
The proposed structure could have many applications in areas that require a greater NPR. For example, filters, energy absorption, lightweight structures, etc. It would also be interesting to analyze the bi−directional loading on the proposed auxetic cellular structure to analyze the combined loading effect.

Author Contributions

Conceptualization, S.Z.K. and E.R.I.M.; Formal analysis, F.M.; Funding acquisition, Y.H.A.; Investigation, A.A.Z. and H.A.; Methodology, S.Z.K., A.A.Z. and Y.H.A.; Project administration, S.Z.K. and A.A.Z.; Resources, S.Z.K. and H.A.; Software, F.M.; Validation, S.Z.K. and Y.H.A.; Writing—original draft, F.M. and S.Z.K.; Writing—review & editing, S.Z.K., E.R.I.M. and H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Islamic University of Madinah, Tamayuz−2 grant number 703.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

There is no external data available with this work. All data, models, and code generated or used during the study appear in the submitted article.

Acknowledgments

The authors would like to express their appreciation for the support provided by the Scientific Research Deanship, Islamic University of Madinah, through the research program Tamayuz-2 grant number 703.

Conflicts of Interest

The authors have no conflict of interest with any entity for producing and working on this paper.

References

  1. Yang, W.; Li, Z.-M.; Shi, W.; Xie, B.-H.; Yang, M.-B. Review on auxetic materials. J. Mater. Sci. 2004, 39, 3269–3279. [Google Scholar] [CrossRef]
  2. Ren, X.; Das, R.; Tran, P.; Ngo, T.D.; Xie, Y.M. Auxetic metamaterials and structures: A review. Smart Mater. Struct. 2018, 27, 23001. [Google Scholar] [CrossRef]
  3. Bronder, S.; Herter, F.; Röhrig, A.; Bähre, D.; Jung, A. Design Study for Multifunctional 3D Re-entrant Auxetics. Adv. Eng. Mater. 2022, 24, 2100816. [Google Scholar] [CrossRef]
  4. Zhang, J.; Zhu, X.; Yang, X.; Zhang, W. Transient nonlinear responses of an auxetic honeycomb sandwich plate under impact loads. Int. J. Impact Eng. 2019, 134, 103383. [Google Scholar] [CrossRef]
  5. Kelkar, P.; Kim, H.S.; Cho, K.-H.; Kwak, J.Y.; Kang, C.-Y.; Song, H.-C. Cellular Auxetic Structures for Mechanical Metamaterials: A Review. Sensors 2020, 20, 3132. [Google Scholar] [CrossRef]
  6. Nguyen, N.V.; Nguyen-Xuan, H.; Nguyen, T.N.; Kang, J.; Lee, J. A comprehensive analysis of auxetic honeycomb sandwich plates with graphene nanoplatelets reinforcement. Compos. Struct. 2021, 259, 113213. [Google Scholar] [CrossRef]
  7. Donoghue, J.P.; Alderson, K.L.; Evans, K.E. The fracture toughness of composite laminates with a negative Poisson’s ratio. Phys. Status Solidi (b) 2009, 246, 2011–2017. [Google Scholar] [CrossRef]
  8. Choi, J.B.; Lakes, R.S. Fracture toughness of re-entrant foam materials with a negative Poisson’s ratio: Experiment and analysis. Int. J. Fract. 1996, 80, 73–83. [Google Scholar] [CrossRef]
  9. Yao, Y.; Yuan, H.; Huang, H.; Liu, J.; Wang, L.; Fan, Y. Biomechanical design and analysis of auxetic pedicle screw to resist loosening. Comput. Biol. Med. 2021, 133, 104386. [Google Scholar] [CrossRef]
  10. Cellular Solids by Lorna J. Gibson. Available online: https://www.cambridge.org/core/books/cellular-solids/BC25789552BAA8E3CAD5E1D105612AB5 (accessed on 19 September 2019).
  11. Masters, I.; Evans, K. Models for the elastic deformation of honeycombs. Compos. Struct. 1996, 35, 403–422. [Google Scholar] [CrossRef]
  12. Yang, L.; Harrysson, O.; West, H.; Cormier, D. Mechanical properties of 3D re-entrant honeycomb auxetic structures realized via additive manufacturing. Int. J. Solids Struct. 2015, 69–70, 475–490. [Google Scholar] [CrossRef]
  13. Khan, S.Z.; Masood, S.; Cottam, R. Mechanical properties in tensile loading of H13 re-entrant honeycomb auxetic structure manufactured by direct metal deposition. MATEC Web Conf. 2015, 34, 01004. [Google Scholar] [CrossRef]
  14. Berinskii, I. Elastic networks to model auxetic properties of cellular materials. Int. J. Mech. Sci. 2016, 115–116, 481–488. [Google Scholar] [CrossRef]
  15. Khan, K.A.; Al-Mansoor, S.; Khan, S.Z. Piezoelectric Metamaterial with Negative and Zero Poisson’s Ratios. J. Eng. Mech. 2019, 145. [Google Scholar] [CrossRef] [Green Version]
  16. Shen, L.; Wang, Z.; Wang, X.; Wei, K. Negative Poisson’s ratio and effective Young’s modulus of a vertex-based hierarchical re-entrant honeycomb structure. Int. J. Mech. Sci. 2021, 206, 106611. [Google Scholar] [CrossRef]
  17. Khan, S.Z.; Mustahsan, F.; Mahmoud, E.R.; Masood, S. A novel modified re-entrant honeycomb structure to enhance the auxetic behavior: Analytical and numerical study by FEA. Mater. Today: Proc. 2021, 39, 1041–1045. [Google Scholar] [CrossRef]
  18. Al-Rifaie, H.; Novak, N.; Vesenjak, M.; Ren, Z.; Sumelka, W. Fabrication and Mechanical Testing of the Uniaxial Graded Auxetic Damper. Materials 2022, 15, 387. [Google Scholar] [CrossRef] [PubMed]
  19. Sigmund, O.; Torquato, S. Design of smart composite materials using topology optimization. Smart Mater. Struct. 1999, 8, 365–379. [Google Scholar] [CrossRef]
  20. Fu, M.-H.; Chen, Y.; Hu, L.-L. A novel auxetic honeycomb with enhanced in-plane stiffness and buckling strength. Compos. Struct. 2017, 160, 574–585. [Google Scholar] [CrossRef]
  21. Harkati, E.; Daoudi, N.; Bezazi, A.; Haddad, A.; Scarpa, F. In-plane elasticity of a multi re-entrant auxetic honeycomb. Compos. Struct. 2017, 180, 130–139. [Google Scholar] [CrossRef]
  22. Bezazi, A.; Scarpa, F.; Remillat, C. A novel centresymmetric honeycomb composite structure. Compos. Struct. 2005, 71, 356–364. [Google Scholar] [CrossRef]
  23. Zhu, Y.; Luo, Y.; Gao, D.; Yu, C.; Ren, X.; Zhang, C. In-plane elastic properties of a novel re-entrant auxetic honeycomb with zigzag inclined ligaments. Eng. Struct. 2022, 268, 114788. [Google Scholar] [CrossRef]
  24. Wu, W.; Liu, P.; Kang, Z. A novel mechanical metamaterial with simultaneous stretching- and compression-expanding property. Mater. Des. 2021, 208, 109930. [Google Scholar] [CrossRef]
  25. Zhang, X.G.; Jiang, W.; Zhang, Y.; Luo, C.; Zhang, X.Y.; Han, D.; Hao, J.; Teng, X.C.; Xie, Y.M.; Ren, X. Energy absorption properties of composite tubes with hexagonal and re-entrant honeycomb fillers. Constr. Build. Mater. 2022, 356, 129298. [Google Scholar] [CrossRef]
  26. Timoshenko, S.; Goodier, J.N. Theory of Elasticity, 3rd ed.; McGraw-Hill Education: Auckland, New Zealand, 1970. [Google Scholar]
  27. Elementary Continuum Mechanics for Everyone. Available online: https://link.springer.com/book/10.1007/978-94-007-5766-0 (accessed on 18 October 2022).
  28. Ghavami, P. Mechanics of Materials; Springer International Publishing: Cham, Switzerland, 2015. [Google Scholar] [CrossRef]
  29. Hussein, M. Effects of strain rate and temperature on the mechanical behavior of carbon black reinforced elastomers based on butyl rubber and high molecular weight polyethylene. Results Phys. 2018, 9, 511–517. [Google Scholar] [CrossRef]
  30. GOM Correlate Pro: Strains and Displacements from Motion Pictures. Available online: https://www.gom.com/en/products/gom-suite/gom-correlate-pro (accessed on 6 March 2022).
Figure 1. Unit cell of the proposed structure.
Figure 1. Unit cell of the proposed structure.
Materials 15 08022 g001
Figure 2. Free body diagram for loading in the x-direction.
Figure 2. Free body diagram for loading in the x-direction.
Materials 15 08022 g002
Figure 3. Free body diagram for loading in the y-direction.
Figure 3. Free body diagram for loading in the y-direction.
Materials 15 08022 g003
Figure 4. Boundary conditions for the FEA with points for measuring deflections.
Figure 4. Boundary conditions for the FEA with points for measuring deflections.
Materials 15 08022 g004
Figure 5. Laser-cut acrylic samples for testing loads in both the x and y-directions.
Figure 5. Laser-cut acrylic samples for testing loads in both the x and y-directions.
Materials 15 08022 g005
Figure 6. (a) Custom mount for holding the structure and (b) tension testing setup.
Figure 6. (a) Custom mount for holding the structure and (b) tension testing setup.
Materials 15 08022 g006
Figure 7. Load vs. displacement graph for acrylic dog-bone samples.
Figure 7. Load vs. displacement graph for acrylic dog-bone samples.
Materials 15 08022 g007
Figure 8. Force vs. displacement plot for loading in the x-direction.
Figure 8. Force vs. displacement plot for loading in the x-direction.
Materials 15 08022 g008
Figure 9. Force vs. displacement plot for loading in the y-direction.
Figure 9. Force vs. displacement plot for loading in the y-direction.
Materials 15 08022 g009
Figure 10. Deflection pattern between the tensile test and FEA for loading in (a,b) the x-direction and (c,d) y-direction.
Figure 10. Deflection pattern between the tensile test and FEA for loading in (a,b) the x-direction and (c,d) y-direction.
Materials 15 08022 g010
Figure 11. Poisson’s ratio obtained via a mathematical model, FEA, and tensile testing for loading in the x-direction.
Figure 11. Poisson’s ratio obtained via a mathematical model, FEA, and tensile testing for loading in the x-direction.
Materials 15 08022 g011
Figure 12. Poisson’s ratio obtained from the FEA and mathematical model for varying H/L for x-direction loading.
Figure 12. Poisson’s ratio obtained from the FEA and mathematical model for varying H/L for x-direction loading.
Materials 15 08022 g012
Figure 13. Poisson’s ratio obtained from the mathematical model and FEA for varying θ for x-direction loading.
Figure 13. Poisson’s ratio obtained from the mathematical model and FEA for varying θ for x-direction loading.
Materials 15 08022 g013
Figure 14. Poisson’s ratio obtained via the mathematical model, FEA, and tensile testing for loading in the y-direction.
Figure 14. Poisson’s ratio obtained via the mathematical model, FEA, and tensile testing for loading in the y-direction.
Materials 15 08022 g014
Figure 15. Comparison of Poisson’s ratio obtained from the FEA and mathematical model for varying H/L for y-direction loading.
Figure 15. Comparison of Poisson’s ratio obtained from the FEA and mathematical model for varying H/L for y-direction loading.
Materials 15 08022 g015
Figure 16. Comparison of Poisson’s ratio obtained from the FEA and mathematical model for varying angle θ for y-direction loading.
Figure 16. Comparison of Poisson’s ratio obtained from the FEA and mathematical model for varying angle θ for y-direction loading.
Materials 15 08022 g016
Figure 17. NPR for different values of H/L, a/L, and re-entrant angle for (a,b) x-direction loading and (c,d) y-direction loading.
Figure 17. NPR for different values of H/L, a/L, and re-entrant angle for (a,b) x-direction loading and (c,d) y-direction loading.
Materials 15 08022 g017
Figure 18. Young’s modulus in the x-direction for different values of H/L and re-entrant angle.
Figure 18. Young’s modulus in the x-direction for different values of H/L and re-entrant angle.
Materials 15 08022 g018
Figure 19. Young’s modulus in they-direction for different values of H/L and re-entrant angle (inset shows a portion of the extreme values).
Figure 19. Young’s modulus in they-direction for different values of H/L and re-entrant angle (inset shows a portion of the extreme values).
Materials 15 08022 g019
Figure 20. Relative density and NPR for loading in both directions.
Figure 20. Relative density and NPR for loading in both directions.
Materials 15 08022 g020
Table 1. Dimensions of the test samples.
Table 1. Dimensions of the test samples.
L H t b a θ H / L t / L a / L Young’s Modulus
E s
mmmmmmmmmmDegreesMPa
8201.753.751.6
2.4
3.2
4.0
602.250.2180.2
0.3
0.4
0.5
1050
Table 2. Comparison of Poisson’s ratio obtained from the mathematical model, FEA, and tensile testing for loading in the x-direction.
Table 2. Comparison of Poisson’s ratio obtained from the mathematical model, FEA, and tensile testing for loading in the x-direction.
a/LMathematical ModelFEATensile Test
0.2−0.45−0.46−0.51
0.3−0.61−0.62−0.67
0.4−0.82−0.83−0.89
0.5−1.06−1.08−1.16
Table 3. Poisson’s ratio obtained from the mathematical model, FEA, and tensile testing for loading in the y-direction.
Table 3. Poisson’s ratio obtained from the mathematical model, FEA, and tensile testing for loading in the y-direction.
a/LMathematical ModelFEATensile Test
0.2−1.36−1.46−1.41
0.3−1.1−1.18−1.12
0.4−0.88−0.94−0.9
0.5−0.7−0.75−0.72
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mustahsan, F.; Khan, S.Z.; Zaidi, A.A.; Alahmadi, Y.H.; Mahmoud, E.R.I.; Almohamadi, H. Re-Entrant Honeycomb Auxetic Structure with Enhanced Directional Properties. Materials 2022, 15, 8022. https://doi.org/10.3390/ma15228022

AMA Style

Mustahsan F, Khan SZ, Zaidi AA, Alahmadi YH, Mahmoud ERI, Almohamadi H. Re-Entrant Honeycomb Auxetic Structure with Enhanced Directional Properties. Materials. 2022; 15(22):8022. https://doi.org/10.3390/ma15228022

Chicago/Turabian Style

Mustahsan, Farrukh, Sohaib Z. Khan, Asad A. Zaidi, Yaser H. Alahmadi, Essam R. I. Mahmoud, and Hamad Almohamadi. 2022. "Re-Entrant Honeycomb Auxetic Structure with Enhanced Directional Properties" Materials 15, no. 22: 8022. https://doi.org/10.3390/ma15228022

APA Style

Mustahsan, F., Khan, S. Z., Zaidi, A. A., Alahmadi, Y. H., Mahmoud, E. R. I., & Almohamadi, H. (2022). Re-Entrant Honeycomb Auxetic Structure with Enhanced Directional Properties. Materials, 15(22), 8022. https://doi.org/10.3390/ma15228022

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