Computational Approach in Formulating Mechanical Characteristics of 3 D Star Honeycomb Auxetic Structure

Auxetic materials exhibit a unique characteristic due to the altered microstructure. Different structures have been used to model thesematerials.This paper treats a development of finite elementmodel and theoretical formulation of 3D star honeycomb structure of these materials. Various shape parameters of the structural cell were evaluated with respect to the basic mechanical properties of the cell. Finite element and analytical approach for various geometrical parameters were numerically used to formulate the characteristics of the material. The study aims at quantifying mechanical properties for any domain in which auxetic material is of interest for variations in geometrical parameters. It is evident that mechanical properties of the material could be controlled by changing the base wall angle of the configuration. The primary outcome of the study is a design guideline for the use of 3D star honeycomb auxetic cellular structure in structural applications.


Introduction
For decades, several geometrical structures representing auxetic behaviour have been introduced, fabricated, and tested particularly to obtain their mechanical properties.These geometrical structures are beneficial since they make contribution to researchers to comprehend on how auxetic behaviour can be obtained and how auxetic materials can be fabricated.Computational procedure is the main contribution in predicting their properties.
Reentrant structures which have been investigated by various researchers [1][2][3][4][5][6], chiral structures [7][8][9], rotating rigid/semi-rigid units [10][11][12][13][14], angle-ply laminates [15,16], hard molecules [17][18][19][20], microporous polymers [21][22][23], and liquid crystalline polymer [24][25][26] are among the most important classes of such auxetic structures.Star honeycomb constituent of reentrant structural group has increasingly received attention compared to other structures due to its exceptional ability to model auxetic materials.Within the same group, other structures which are different from star honeycomb structures are double arrow head structures, structurally hexagonal honeycomb, structures formed from lozenge grids, structure formed from square grids, and structure formed from sinusoidal ligaments [27].Man-made auxetic materials were first designed by Theocaris et al. in 1997 [28].In their study, the design of material was based on structure of arrays with polygonal-shaped inclusions and reentrant corners.Based on their research, the moduli of the inclusion material should be lower than the moduli of the material of structure in order to make composite with desired negative Poisson's ratio.Using numerical analysis, they proved that the shape of the star-shaped microinclusions is the main factor causing auxetic behaviour.
In the present study, a simple 3D star honeycomb cellular structure was modelled by adapting a 2D cellular structure.Subsequently, the influence of structure shape on basic properties has also been evaluated.In particular, the basic mechanical properties discussed in this study are elastic modulus, density ratio, and Poisson's ratio which were examined for different values of shape parameters, namely, the wall length, the cell wall thicknesses, and the cellular structure angle.Consisting of eight elastic beams, the 2D cellular structure has a symmetrical configuration.Analytical formulation of mechanical properties for the structure has eventually been developed.This investigation highlights finite element discretization and the developed formulation of 3D star honeycomb structures for auxetic materials that have not been proposed previously.From this study, it is evident that mechanical properties of auxetic material are controllable, thus facilitating the fabrication technique used in preparing samples in the laboratory.

Definition of the Auxetic Model
Figure 1 shows beam-like 2D and 3D star honeycomb structures formed by eight (8) elastic beams with the same length and cross-sectional area.By assembling such single unit of structural cell, an element of auxetic material made of 3D star honeycomb structure is produced as depicted in Figure 1(b).This structural configuration is deemed as cellular structure of a material.From the mechanics point of view, mechanical properties of the cellular material depend on both geometrical parameters of the cells and mechanical properties of the base material [3].This study has attempted to formulate some of the important mechanical properties of the auxetic material.These parameters include Poisson's ratio, elastic modulus, and density ratio.Based on symmetrical properties of the 3D cellular structures, it can be assumed that uniaxial loading of the cell along  direction leads to the same deflection along both  and  directions.Also, uniaxial loading along  and  directions results in the same deflections in the lateral directions.As a result, the following set of equations can be written: The following equation should be satisfied to obtain a starshaped structure as shown in Figure 1:

Finite Element Modelling
An elastic cellular structure has numerically been developed using finite element technique to determine the basic properties of the star honeycomb structures.ABAQUS finite element code has been used to perform these analyses.An element of the auxetic material shown in Figure 1 was taken into consideration.The element consists of 27 cellular structures.Element used is the linear tetrahedral, type C3D4.A total number of elements and nodes are 2,916,604 and 746,044, respectively.The mesh size was determined from the convergence study.It is evident that finite element analysis of this element leads to more accurate and reliable results compared to the analysis of single cellular structure.Furthermore, the element was loaded along  with the imposed normal deflections.The resulting stresses as well as normal and transverse strains could then be obtained.Meanwhile, the Young's modulus was calculated by dividing the averaged uniaxial stresses by the uniaxial strains while the Poisson's ratios were obtained by dividing the averaged transverse strain by the uniaxial imposed ones.Several finite element models were consecutively developed for different values of the length, cross section, and angles.The Poisson's ratios and elastic modulus including ]  , ]  , and  1 were computed using the same method.Figure 2 shows finite element modelling of Model number 1.All elements were then loaded with imposed normal deflections in  direction.The Poisson's ratios and elastic modulus including ]  , ]  , and  2 were computed using the same method.Figure 3 shows finite element modelling of the same model for loading along  direction.Due to the symmetry of the element, loading along  and  direction provides the same results.As such, the loading along  direction has not been performed in the present study.

Development of Analytical Formulation
Using analytical approach, Poisson's ratios, elastic modulus, and density ratios of all models may also be predicted.From mathematics and mechanics point of views, there is relationship between the geometrical parameters of the cell, mechanical properties of the base material, and mechanical properties of the star honeycomb structure including Poisson's ratio, elastic modulus, and density ratio.The relationships can be expressed as the following functions: where  1 ,  2 , , and  are geometrical parameters of the cell shown in Figure 1. and   are elastic modulus and density of the base material, respectively.]  , ]  , ]  , ]  , ]  , ]  ,  1 ,  2 , and  3 are Poisson's ratios and elastic modulus of the material.Using dimensionless parameters, the functions of Poisson's ratios and elastic modulus of the material could be predicted as follows: where  = ,  =  2 / 1 , and  = / 1 .For the validity of a set of (4), each function must have twenty-seven (27) coefficients.To obtain these coefficients, values of Poisson's ratios, elastic modulus, and density ratio for twenty-seven (27) models with different geometrical parameters are required.Finite element analysis has been carried out for the models.Table 1 shows the values of , , and  used for inclusion into finite element analysis.
By using a set of (4) for all models, the matrix form of ]  function is written as follows: where (]  )  is ]  of Model number .Using the matrix form in (5), coefficients  000 to  222 can then be determined.Similarly, the coefficients of ]  function can be calculated from the following matrix form: in which (]  )  refers to ]  for Model number .Obtaining ]  from finite element technique for all models also gives the coefficients of ]  function using the following matrix form: where (]  )  refers to ]  for Model number .
Also, the coefficients of ]  function can be determined from the following matrix form: Elastic modulus of all models along  direction gives the coefficients of ( 1 /) function ( 000 to  222 ) from the following equation: Similarly, the following matrix form can be written to calculate coefficients of  2 / function ( 000 to  222 ): Advances in Materials Science and Engineering  Density ratio coefficients can also be calculated from the following matrix form: where   is density of Model number  and   is the density of the base material.
As illustrated, Poisson's ratio and elastic modulus of all defined models have been determined by using finite element approach.In terms of density ratio of the models, the formula of density for auxetic material was used to calculate density ratio.In particular, dependency of material density on geometrical parameters is presented as follows: Using dimensionless geometrical parameters, (12) can be simplified as follows: After calculating the mechanical properties for all models, the coefficients   to   were then determined using ( 5) to (11).

Numerical Analysis of Auxetic Star Honeycomb Structure
This present study has focused on Poisson's ratio, Young's modulus, and the density of a 3D element of auxetic material made of star honeycomb structure using numerical method in conjunction with analytical approach.Details of the results obtained for the models consisting of cells in 27 cases are also presented in Table 1.Some of the numerical results are in details shown in Table 2 and Figures 2, 3, 4, and 5.The primary outcomes of the finite element analyses could be remarked as follows.
(i) From  = 45 ∘ to 70 ∘ , the structure exhibits auxetic behaviour in all directions.
(iv) In terms of elastic modulus and density,  and  decrease with the increasing of .
(v) Dimensionless geometrical parameters used in this research are  = ,  =  2 / 1 , and  = / 1 .Taking into consideration the values of these parameters used for formulation, using the obtained formulas for designing star-shaped structures in which domain of dimensionless geometrical parameters is 0.8202 ≤  ≤ 1.5184, 0.8 ≤  ≤ 0.9, and 0.02 ≤  ≤ 0.06 leads to reasonable results.However, these formulas can only be used to estimate the mechanical properties of star honeycomb structures with dimensionless geometric parameters close to the values of the parameters used in the formulation.
(vi) As illustrated, the aim of this study is to introduce a method for formulating beam like star honeycomb structures.Although the result of that is reasonable for limited domain of geometrical parameters, it offers a method for formulating mechanical properties for any desired domain in which the auxetic material is about to be designed and fabricated.
(vii) Although the shape of the defined structure seems to exhibit auxetic behaviour in all directions, it has interestingly found that it does not exhibit auxetic behaviour in some cases depending on the values of geometrical parameters.For instance, from  = 70 ∘ onward, the structure does not exhibit auxetic behaviour.

Advances in Materials Science and Engineering
For further details, the result of formulation used to obtain mathematical formula of mechanical properties is shown in Table 3.

Conclusion
A symmetrical configuration for the star honeycomb cellular structure of an auxetic material has been proposed and subsequently its analytical principle has also been formulated.The topology allows more degrees of freedom to the designer in controlling the mechanical properties of auxetic material by varying the wall lengths, cell wall thickness, and angle.Manipulation of stiffness properties including Poisson's ratio and fracture toughness can be performed by modification of the cellular structure's structural parameters.The overall mechanical properties can be controlled by modification of the base wall angle of the configurations, with the consequent changes, in particular, of the in-plane Poisson's ratio.This research provides a guideline for applications to 3D star honeycomb auxetic structures and also provides a basis for experimental additive manufacturing processes such as electron beam melting or selective laser sintering.

Figure 1 :
Figure 1: (a) Star honeycomb structure in 2D and 3D.(b) An element of auxetic material made of 3D star honeycomb structures.

Figure 2 :
Figure 2: Von Misses stress of Model number 1 in the case of uniaxial imposed displacement along  direction.

Figure 3 :
Figure 3: Von Misses stress of Model number 1 in the case of uniaxial imposed displacement along  direction.

Figure 4 : 7 −Figure 5 :
Figure 4: Mechanical properties versus  in the case of fixing  and  at 0.8 and 0.02, respectively.

Table 1 :
Values of , , and  used in analysed models.

Table 2 :
Finite element results of Poisson's ratios, elastic modulus, and density ratio of models.