Research Article  Open Access
Mozafar Shokri Rad, Zaini Ahmad, Amran Alias, "Computational Approach in Formulating Mechanical Characteristics of 3D Star Honeycomb Auxetic Structure", Advances in Materials Science and Engineering, vol. 2015, Article ID 650769, 11 pages, 2015. https://doi.org/10.1155/2015/650769
Computational Approach in Formulating Mechanical Characteristics of 3D Star Honeycomb Auxetic Structure
Abstract
Auxetic materials exhibit a unique characteristic due to the altered microstructure. Different structures have been used to model these materials. This paper treats a development of finite element model 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.
1. 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â€“6], chiral structures [7â€“9], rotating rigid/semirigid units [10â€“14], angleply laminates [15, 16], hard molecules [17â€“20], microporous polymers [21â€“23], and liquid crystalline polymer [24â€“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]. Manmade 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 polygonalshaped 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 starshaped 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.
2. Definition of the Auxetic Model
Figure 1 shows beamlike 2D and 3D star honeycomb structures formed by eight (8) elastic beams with the same length and crosssectional 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:
(a)
(b)
3. 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 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 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.
4. 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 , , , and are geometrical parameters of the cell shown in Figure 1. and are elastic modulus and density of the base material, respectively. , , , , , , , , and 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 , , and .
For the validity of a set of (4), each function must have twentyseven (27) coefficients. To obtain these coefficients, values of Poissonâ€™s ratios, elastic modulus, and density ratio for twentyseven (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 to 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 function ( to ) from the following equation: Similarly, the following matrix form can be written to calculate coefficients of function ( to ): 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).
5. 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.

(a)
(b)
(c)
(d)
(e)
(f)
(g)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
The primary outcomes of the finite element analyses could be remarked as follows.(i)From to , the structure exhibits auxetic behaviour in all directions.(ii)From to , is negative and other Poissonâ€™s ratios are positive.(iii)From to , the auxeticity decreases as the angle increases.(iv)In terms of elastic modulus and density, and decrease with the increasing of .(v)Dimensionless geometrical parameters used in this research are , , and . Taking into consideration the values of these parameters used for formulation, using the obtained formulas for designing starshaped structures in which domain of dimensionless geometrical parameters is , , and 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 onward, the structure does not exhibit auxetic behaviour.For further details, the result of formulation used to obtain mathematical formula of mechanical properties is shown in Table 3.