# Effect of Fiber Geometry and Representative Volume Element on Elastic and Thermal Properties of Unidirectional Fiber-Reinforced Composites

Authors: Siva Bhaskara Rao Devireddy; Sandhyarani Biswas
Journal: Journal of Composites
 (2014)

---

## Abstract

The aim of present work is focused on the evaluation of elastic and thermal properties of unidirectional fiber-reinforced polymer composites with different volume fractions of fiber up to 0.7 using micromechanical approach. Two ways for calculating the material properties, that is, analytical and numerical approaches, were presented. In numerical approach, finite element analysis was used to evaluate the elastic modulus and thermal conductivity of composite from the constituent material properties. The finite element model based on three-dimensional micromechanical representative volume element (RVE) with a square and hexagonal packing geometry was implemented by using finite element code ANSYS. Circular cross section of fiber and square cross section of fiber were considered to develop RVE. The periodic boundary conditions are applied to the RVE to calculate elastic modulus of composite. The steady state heat transfer simulations were performed in thermal analysis to calculate thermal conductivity of composite. In analytical approach, the elastic modulus is calculated by rule of mixture, Halpin-Tsai model, and periodic microstructure. Thermal conductivity is calculated analytically by using rule of mixture, the Chawla model, and the Hashin model. The material properties obtained using finite element techniques were compared with different analytical methods and good agreement was achieved. The results are affected by a number of parameters such as volume fraction of the fibers, geometry of fiber, and RVE.


---

## Body

## 1. Introduction

There has been a considerable increase in the use of fiber composite materials in various industries like aerospace, automotive, infrastructures, and sporting goods due to their specific properties like strength, stiffness, toughness, high corrosion resistance, high wear resistance, high chemical resistance, and reduced cost. These materials can take advantage of different properties of their constituents, microstructure, and interaction between constituents in order to improve the mechanical behavior of parts made from them. The mechanics of fiber-reinforced composites are complex due to their anisotropic and heterogeneous characteristics. The evaluation of effective mechanical and thermal properties of composite materials is of paramount importance in engineering design and application. Generally, two approaches are considered in obtaining the global properties of composites: (a) macromechanical analysis and (b) micromechanical analysis. In macromechanical analysis the composite material is considered as a homogeneous orthotropic continuum. In micromechanical analysis the study of composite material is at the fiber and matrix level. Typically the unit cell technique combined with the known material properties of fiber and matrix is used to determine the overall behavior of the composite [1]. A number of methods have been developed to predict and to simulate the mechanical and thermal behavior of composites. Basic analytical approaches have been reported [2–4] to predict the composite materials properties, for example, strength, stiffness, and thermal conductivity. Prediction of boron and aluminium composite properties from a representative volume element (RVE) with square and hexagonal geometry has been reported [5]. Micromechanical analysis of unidirectional fiber-reinforced composites with square and hexagonal unit cells has been reported [6, 7] to evaluate the effective material properties.Patnaik et al. [8] studied the micromechanical and thermal characteristics of glass-fiber-reinforced polymer composites. The experimental results were in good agreement with finite element model based on representative area element approach. Melro et al. [9] predict the inelastic deformation and fracture of randomly distributed unidirectional fiber-reinforced polymer composite materials with different RVEs. Several researchers made their significant contribution in studying thermal characteristics of fiber composites. Springer and Tsai [10] studied the composite thermal conductivities of unidirectional composites and expressions are obtained for predicting these conductivities in the directions along the filaments and normal to them. Islam and Pramila [11] predict the effective transverse thermal conductivity of fiber-reinforced composites by using finite element method. Square and circular cross section fibers were used for perfect bonding at fiber-matrix interface and with interfacial barrier by using four different sets of thermal boundary conditions. Al-Sulaiman et al. [12] predict the thermal conductivity of the constituents of fiber-reinforced composite laminates using three empirical formulas. Grove [13] computed transverse thermal conductivity in continuous unidirectional fiber composite materials using finite element and spatial statistical techniques for a range of fiber volume fractions up to 0.5. Lu [14] used boundary collocation scheme for calculation of transverse effective thermal conductivity of 2-dimensional periodic arrays of long circular and square cylinders with square array and long circular cylinders with hexagonal array for a complete range of fiber volume fractions.Although a great deal of work has already been done on fiber-reinforced polymer composites with circular cross-section of fiber, square cross section of fiber with different RVE models using finite element analysis is hardly been reported. To this end, the objective of the present work is developing a three-dimensional micromechanical RVE with a square and hexagonal packing geometry with circular and square fiber cross sections. A numerical homogenization technique based on the finite element analysis was used to evaluate the elastic modulus and thermal conductivity of composite. The finite element results are compared with the analytical methods. The aim is to demonstrate applicability of homogenization technique by using finite element method to predict material characteristics in advance.


## 2. Materials and Methods

In this present investigation, unidirectional glass fiber as reinforcement phase and epoxy as matrix phase for the composite material were considered. The fiber and matrix materials are considered as isotropic and homogeneous. The properties of the constituent materials are as shown in Table1. In a real unidirectional fiber-reinforced composite, the fibers are arranged randomly and it is difficult to model random fiber arrangement.Table 1
Mechanical and thermal properties of the constituent materials [8, 9].
Properties
Glass fiber
Epoxy matrix
Density (g/cm3)
2.5
1.15
Young’s modulus (GPa)
73
3.76
Poisson’s ratio
0.2
0.39
Shear modulus (GPa)
30.42
1.28
Thermal conductivity (W/mk)
1.3
0.363For this analysis, circular and square cross section fiber composite material is considered. The schematic diagram of the unidirectional fiber composite where the fibers are arranged in the square and hexagonal array is shown in Figure1. By varying the volume fraction of fiber from 0.1 to 0.7 the elastic and thermal properties of composite material are determined.Arrangement of fibers in (a) square array with circular fibers, (b) hexagonal array with circular fibers, (c) square array with square fibers, and (d) hexagonal array with square fibers.
(a)
(b)
(c)
(d)


### 2.1. Constitutive Equations for Fiber Composite Material

The most general form of the anisotropic constitutive equations for homogeneous and elastic composite materials is given by Hook’s law as shown in (1) [15]. Consider
(1)σ1σ2σ3σ4σ5σ6≡σxxσyyσzzτyzτzxτxy=C11C12C13C14C15C16C12C22C23C24C25C26C13C23C33C34C35C36C14C24C34C44C45C46C15C25C35C45C55C56C16C26C36C46C56C66εxxεyyεzzγyzγzxγxy≡ε1ε2ε3ε4ε5ε6,
where σij and τij are normal and shear components of stress, respectively, εij and γij are the normal and shear components of strain, respectively, and Cij is the symmetric stiffness matrix with 21 independent, elastic constants. According to their behaviour, composites may be characterized as generally anisotropic, monoclinic, orthotropic, and transversely isotropic. In present work, transversely isotropic characteristics have been considered for the fiber-reinforced composite. A transversely isotropic material is to be a material whose effective properties are isotropic in one of its planes and the stiffness tensor is represented in
(2)σ1σ2σ3σ4σ5σ6=C11C12C12000C12C22C23000C12C23C2200000012C22-C23000000C66000000C66ε1ε2ε3ε4ε5ε6.
Once the components of the transversely isotropic stiffness tensor C are known, the elastic properties of homogenized material can be computed by (3) [16]. Consider
(3)E1=C11-2C122C22+C23,E2=C11C22+C23-2C122C22-C23C11C22-C122,ν12=C12C22+C23,G23=12C22-C23,
where E1, E2, ν12, and G23 are longitudinal modulus, transverse modulus in plane Poisson’s ratio, and in plane shear modulus, respectively.


### 2.2. Generation of RVE

For simplicity reasons, most micromechanical models assume a periodic arrangement of fibers for which a RVE or unit cell can be isolated. The RVE has the same elastic constants and fiber volume fraction as the composite. The periodic fiber sequences commonly used are the square array and the hexagonal array. For a square packing RVE as shown in Figures2(a) and 2(c) the maximum theoretically achievable fiber volume fraction is 78.54%. For square RVE the diameter of fiber is calculated by
(4)Vf=a1π/4df2a1a2a3,
where Vf is volume fraction of fiber; a1, a2, and a3 are the length of square RVE; and df is the diameter of fiber. For the hexagonal packing RVE as shown in Figures 2(b) and 2(d) the maximum theoretically achievable fiber volume fraction is 90.69%. Obviously, with a hexagonal packing geometry a composite can be made more compact than with a square packing geometry. For hexagonal RVE the diameter of fiber is calculated by
(5)Vf=2a1π/4df2a1a2a3,
where a3=a2tan⁡(60°) and a2=4a1.RVE (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)


## 2.1. Constitutive Equations for Fiber Composite Material

The most general form of the anisotropic constitutive equations for homogeneous and elastic composite materials is given by Hook’s law as shown in (1) [15]. Consider
(1)σ1σ2σ3σ4σ5σ6≡σxxσyyσzzτyzτzxτxy=C11C12C13C14C15C16C12C22C23C24C25C26C13C23C33C34C35C36C14C24C34C44C45C46C15C25C35C45C55C56C16C26C36C46C56C66εxxεyyεzzγyzγzxγxy≡ε1ε2ε3ε4ε5ε6,
where σij and τij are normal and shear components of stress, respectively, εij and γij are the normal and shear components of strain, respectively, and Cij is the symmetric stiffness matrix with 21 independent, elastic constants. According to their behaviour, composites may be characterized as generally anisotropic, monoclinic, orthotropic, and transversely isotropic. In present work, transversely isotropic characteristics have been considered for the fiber-reinforced composite. A transversely isotropic material is to be a material whose effective properties are isotropic in one of its planes and the stiffness tensor is represented in
(2)σ1σ2σ3σ4σ5σ6=C11C12C12000C12C22C23000C12C23C2200000012C22-C23000000C66000000C66ε1ε2ε3ε4ε5ε6.
Once the components of the transversely isotropic stiffness tensor C are known, the elastic properties of homogenized material can be computed by (3) [16]. Consider
(3)E1=C11-2C122C22+C23,E2=C11C22+C23-2C122C22-C23C11C22-C122,ν12=C12C22+C23,G23=12C22-C23,
where E1, E2, ν12, and G23 are longitudinal modulus, transverse modulus in plane Poisson’s ratio, and in plane shear modulus, respectively.


## 2.2. Generation of RVE

For simplicity reasons, most micromechanical models assume a periodic arrangement of fibers for which a RVE or unit cell can be isolated. The RVE has the same elastic constants and fiber volume fraction as the composite. The periodic fiber sequences commonly used are the square array and the hexagonal array. For a square packing RVE as shown in Figures2(a) and 2(c) the maximum theoretically achievable fiber volume fraction is 78.54%. For square RVE the diameter of fiber is calculated by
(4)Vf=a1π/4df2a1a2a3,
where Vf is volume fraction of fiber; a1, a2, and a3 are the length of square RVE; and df is the diameter of fiber. For the hexagonal packing RVE as shown in Figures 2(b) and 2(d) the maximum theoretically achievable fiber volume fraction is 90.69%. Obviously, with a hexagonal packing geometry a composite can be made more compact than with a square packing geometry. For hexagonal RVE the diameter of fiber is calculated by
(5)Vf=2a1π/4df2a1a2a3,
where a3=a2tan⁡(60°) and a2=4a1.RVE (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)


## 3. Finite Element Modeling

In order to evaluate the effective properties of composite, the finite element software package ANSYS is used. The program is written in APDL (ANSYS Programming Design Language), which is delivered by the software and it makes the handling much more comfortable. For simplification, there are many assumptions considered for the present analysis such as fibers which are arranged in a particular pattern (square and hexagonal) in a matrix. The composite is free of voids and other irregularities, all fibers are uniformly distributed in the matrix and perfectly aligned, and the interface between the fiber and matrix is perfectly debonded. In the study of the micromechanics of fiber-reinforced materials, it is convenient to use an orthogonal coordinate system that has one axis aligned with the fiber direction. The axis 1 is aligned with the fiber direction, the axis 2 is in the plane of the RVE and is perpendicular to the fibers, and the axis 3 is perpendicular to the plane of the RVE and is also perpendicular to the fibers as shown in Figure2. Dimensions considered for the analysis are a1=1.0×10-5 m, a2=1.0×10-5 m, and a3=1.0×10-5 m for square RVE. For hexagonal RVE a1=1.0×10-5 m, a2 and a3 are calculated by using (5). The radius of fibers corresponds to volume fractions ranging from 0.1 to 0.7. Three-dimensional structural solid element SOLID186 is used to determine elastic properties and is defined by 20 nodes having three degrees of freedom at each node. They are translations in the nodal 1, 2, and 3 directions. For thermal conductivity a three-dimensional quadratic brick element SOLID90 is used for discretization of the constituents and is defined by 20 nodes with a single degree of freedom (temperature) at each node. The meshed model of square and hexagonal RVE at 0.4 of fiber volume fraction is shown in Figure 3.Meshed model of (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)


### 3.1. Boundary Conditions for Evaluation of Elastic Properties

Composite materials can be represented as a periodic array of the RVEs. Therefore, the periodic boundary conditions must be applied to the RVE models. This implies that each RVE in the composite has the same deformation mode and there is no separation or overlap between the neighboring RVEs after deformation [17, 18]. The resumed boundary conditions applied are given in Table 2. Note that U, V, and W are the displacements along 1, 2, and 3 directions, respectively, applied on the AEDH, BFCG, ABCD, EFGH, DHGC, and AEFB faces as shown in Figure 2. After applying boundary conditions and the displacement constant, the corresponding engineering constants are calculated as follows in terms of corresponding stresses and strains shown in (6):
(6)σij¯=1V∫VσijdV,εij¯=1V∫VεijdV,
where σ¯ and ε¯ are the average stresses and average strains and V is the volume of the RVE. The elastic properties can be calculated by using the constitutive equations of the material properties as the ratio of corresponding average stresses and average strains as shown in (3). Figures 4 and 5 show the counter of stress and strain in square and hexagonal RVE at 0.4 of volume fraction.Table 2
Boundary conditions along the 1, 2, and 3 directions of the RVE.
Load
Constraint
1-direction faces
2-direction faces
3-direction faces
AEDH
BFCG
ABCD
EFGH
DHGC
AEFB
Load 11
U
0
a
1
V
0
0
W
0
0
Load 22
U
0
0
V
0
a
2
W
0
0
Load 33
U
0
0
V
0
0
W
0
a
3Counter of stress in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)Counter of strain in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)


### 3.2. Boundary Conditions for Evaluation of Thermal Conductivity

The steady state heat transfer simulations are performed by using finite element analysis to predict thermal conductivity of composite along the longitudinal and transverse direction. The thermal boundary conditions considered in the present analysis are shown in Figure2. One wall (BFGC) is kept isothermal at elevated temperature, while the corresponding wall (AEHD) is subjected 100 Kelvin to maintain temperature difference for the calculation of longitudinal thermal conductivity. All other surfaces are subjected to insulation boundary conditions. One wall (ABCD) is kept isothermal at elevated temperature, while the corresponding wall (EFGH) is subjected 100 Kelvin to maintain temperature difference for the calculation of transverse thermal conductivity. All other surfaces are subjected to insulation boundary conditions. Using these temperature gradients the heat flux was obtained from ANSYS software.The effective thermal conductivity is established from the fundamental heat conduction law, found by Fourier’s, which states that the heat flux is proportional to the temperature gradient [19]:
(7)q=-KdTdx,
where dT/dx is temperature gradient between two isothermal surfaces, q is the heat flux W/m2, and k is the thermal conductivity. The temperature distribution in square and hexagonal RVE along longitudinal and transverse direction at 0.4 of volume fraction of fiber is shown in Figures 6 and 7.Longitudinal temperature distribution in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)Transverse temperature distribution in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)


## 3.1. Boundary Conditions for Evaluation of Elastic Properties

Composite materials can be represented as a periodic array of the RVEs. Therefore, the periodic boundary conditions must be applied to the RVE models. This implies that each RVE in the composite has the same deformation mode and there is no separation or overlap between the neighboring RVEs after deformation [17, 18]. The resumed boundary conditions applied are given in Table 2. Note that U, V, and W are the displacements along 1, 2, and 3 directions, respectively, applied on the AEDH, BFCG, ABCD, EFGH, DHGC, and AEFB faces as shown in Figure 2. After applying boundary conditions and the displacement constant, the corresponding engineering constants are calculated as follows in terms of corresponding stresses and strains shown in (6):
(6)σij¯=1V∫VσijdV,εij¯=1V∫VεijdV,
where σ¯ and ε¯ are the average stresses and average strains and V is the volume of the RVE. The elastic properties can be calculated by using the constitutive equations of the material properties as the ratio of corresponding average stresses and average strains as shown in (3). Figures 4 and 5 show the counter of stress and strain in square and hexagonal RVE at 0.4 of volume fraction.Table 2
Boundary conditions along the 1, 2, and 3 directions of the RVE.
Load
Constraint
1-direction faces
2-direction faces
3-direction faces
AEDH
BFCG
ABCD
EFGH
DHGC
AEFB
Load 11
U
0
a
1
V
0
0
W
0
0
Load 22
U
0
0
V
0
a
2
W
0
0
Load 33
U
0
0
V
0
0
W
0
a
3Counter of stress in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)Counter of strain in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)


## 3.2. Boundary Conditions for Evaluation of Thermal Conductivity

The steady state heat transfer simulations are performed by using finite element analysis to predict thermal conductivity of composite along the longitudinal and transverse direction. The thermal boundary conditions considered in the present analysis are shown in Figure2. One wall (BFGC) is kept isothermal at elevated temperature, while the corresponding wall (AEHD) is subjected 100 Kelvin to maintain temperature difference for the calculation of longitudinal thermal conductivity. All other surfaces are subjected to insulation boundary conditions. One wall (ABCD) is kept isothermal at elevated temperature, while the corresponding wall (EFGH) is subjected 100 Kelvin to maintain temperature difference for the calculation of transverse thermal conductivity. All other surfaces are subjected to insulation boundary conditions. Using these temperature gradients the heat flux was obtained from ANSYS software.The effective thermal conductivity is established from the fundamental heat conduction law, found by Fourier’s, which states that the heat flux is proportional to the temperature gradient [19]:
(7)q=-KdTdx,
where dT/dx is temperature gradient between two isothermal surfaces, q is the heat flux W/m2, and k is the thermal conductivity. The temperature distribution in square and hexagonal RVE along longitudinal and transverse direction at 0.4 of volume fraction of fiber is shown in Figures 6 and 7.Longitudinal temperature distribution in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)Transverse temperature distribution in (a) square RVE with circular fiber, (b) hexagonal RVE with circular fiber, (c) square RVE with square fiber, and (d) hexagonal RVE with square fiber.
(a)
(b)
(c)
(d)


## 4. Analytical Methods

In order to validate the finite element results, three well-existing analytical methods such as rule of mixture, semiempirical model, and periodic microstructure models have been used for the current study.


### 4.1. Rule of Mixture Model

The well-known models that have been proposed and used to evaluate the properties of unidirectional composites are the Voigt and Reuss models which are also known as the rule of mixture model and the inverse rule of mixture model. The rule of mixture is the mathematical expressions which give the property of the composite in terms of the properties, quantity, and arrangement of its constituents. The longitudinal property of the composite is calculated by using rule of mixture as shown in(8)ϕ1=ϕfVf+ϕmVm,
where ϕ1 is the property of composite in longitudinal direction. ϕf and ϕm are the material properties of fiber and epoxy matrix. Vf and Vm are the volume fractions of fiber and epoxy matrix. The transverse modulus of composite is calculated by using rule of mixture as shown in
(9)1ϕ2=Vfϕf+Vmϕm,
where ϕ2 is the property of composite in transverse direction.


### 4.2. Halpin-Tsai Model

The most useful semiempirical models are those by Halpin and Tsai. Halpin and Tsai developed their models as simple equations by curve fitting to results that are based on elasticity. The material properties calculated in longitudinal direction by using Halpin-Tsai model are the same as those in rule of mixtures. For the transverse direction, the material property of composite is calculated by(10)ϕ2ϕm=1+ζηVf1-ηVf,
where η=((ϕf/ϕm)-1)/((ϕf/ϕm)+ζ).Halpin and Tsai found that the valueξ=2 gave an excellent fit to the finite difference elasticity solution of Adams and Doner [20] for the transverse modulus of a square array of circular fibers.


### 4.3. Periodic Microstructure

If the composite has periodic microstructure, then Fourier series can be used to estimate all the components of the stiffness tensor of the composite. Explicit formulas for a composite reinforced by long circular cylindrical fibers were discussed in [21].


### 4.4. Thermal Conductivity Models

Thermal conductivity obtained by finite element analysis, comparing with three well-existing analytical methods such as rule of mixture, the Chawla model, and the Hashin model [22], has been used for the current study.


## 4.1. Rule of Mixture Model

The well-known models that have been proposed and used to evaluate the properties of unidirectional composites are the Voigt and Reuss models which are also known as the rule of mixture model and the inverse rule of mixture model. The rule of mixture is the mathematical expressions which give the property of the composite in terms of the properties, quantity, and arrangement of its constituents. The longitudinal property of the composite is calculated by using rule of mixture as shown in(8)ϕ1=ϕfVf+ϕmVm,
where ϕ1 is the property of composite in longitudinal direction. ϕf and ϕm are the material properties of fiber and epoxy matrix. Vf and Vm are the volume fractions of fiber and epoxy matrix. The transverse modulus of composite is calculated by using rule of mixture as shown in
(9)1ϕ2=Vfϕf+Vmϕm,
where ϕ2 is the property of composite in transverse direction.


## 4.2. Halpin-Tsai Model

The most useful semiempirical models are those by Halpin and Tsai. Halpin and Tsai developed their models as simple equations by curve fitting to results that are based on elasticity. The material properties calculated in longitudinal direction by using Halpin-Tsai model are the same as those in rule of mixtures. For the transverse direction, the material property of composite is calculated by(10)ϕ2ϕm=1+ζηVf1-ηVf,
where η=((ϕf/ϕm)-1)/((ϕf/ϕm)+ζ).Halpin and Tsai found that the valueξ=2 gave an excellent fit to the finite difference elasticity solution of Adams and Doner [20] for the transverse modulus of a square array of circular fibers.


## 4.3. Periodic Microstructure

If the composite has periodic microstructure, then Fourier series can be used to estimate all the components of the stiffness tensor of the composite. Explicit formulas for a composite reinforced by long circular cylindrical fibers were discussed in [21].


## 4.4. Thermal Conductivity Models

Thermal conductivity obtained by finite element analysis, comparing with three well-existing analytical methods such as rule of mixture, the Chawla model, and the Hashin model [22], has been used for the current study.


## 5. Results and Discussion

### 5.1. Effect of Volume Fraction on Elastic Properties

The longitudinal modulus is the response of composite during the application of load parallel to the fiber direction. It can be defined as the ratio of longitudinal stress to the longitudinal strain. Figure8 shows the effect of fiber content on the longitudinal modulus of composites using rule of mixtures, Halpin-Tsai model, periodic microstructure, and finite element analysis with square and hexagonal RVE with circular and square fiber geometry. It can be observed from the figure that the longitudinal modulus increases with the increase in volume fraction of fiber and there is a good agreement between finite element results and analytical methods. This is because the stiffness of the composite increases with increase in volume fraction of fiber.Figure 8
Longitudinal modulus validation with different volume fraction of fiber.The transverse modulus is the response of composite during the application of load perpendicular to the fiber direction. It can be defined as the ratio of transverse stress to the transverse strain. The effect of fiber volume fraction on transverse modulus of composites using finite element analysis and three analytical methods is graphically shown in Figure9. As expected, it is clear from the figure that the transverse modulus increases with increase in fiber volume fraction. The transverse modulus evaluated by finite element analysis with hexagonal RVE is more close to the Halpin-Tsai model and periodic microstructure as compared to the results obtained from rule of mixture and finite element analysis with square RVE.Figure 9
Transverse modulus validation with different volume fraction of fiber.Figure10 shows the effect of fiber volume fraction on the in-plane Poisson’s ratio of composite. It is evident from the figure that the major Poisson’s ratio decreases with increase in the volume fraction of fiber due to increase in material resistance. The finite element results are in good agreement with analytical methods.Figure 10
In-plane Poisson’s ratio validation with different volume fraction of fiber.In-plane shear modulus of composite is the ratio of shear stress to the shear strain in longitudinal direction. Figure11 shows the effect of fiber volume fraction on the in-plane shear modulus of composite. It is clear from the figure that the shear modulus increases with increases in fiber volume fraction. Also, it can be observed that there is a good agreement between results obtained from finite element analysis and hexagonal RVE with Halpin-Tsai model and periodic microstructure as compared to rule of mixture and finite element analysis with square RVE.Figure 11
In-plane shear modulus validation with different volume fraction of fiber.


### 5.2. Effect of Volume Fraction on Thermal Conductivity

Longitudinal thermal conductivity of composite is the property of a material to conduct heat in parallel to the direction of the fibers. Figure12 shows the effect of fiber content on the longitudinal thermal conductivity using rule of mixture, the Chawla model, the Hashin model, and finite element analysis. Finite element results match well the analytical methods and increases in linear manner with the increase in volume fraction since the fiber thermal conductivity dominates the matrix thermal conductivity.Figure 12
Longitudinal thermal conductivity validation with different volume fraction of fiber.Transverse thermal conductivity of composite is the property of a material to conduct heat, in the direction perpendicular to the fibers. Figure13 shows the effect of fiber content on the transverse thermal conductivity using analytical methods and finite element analysis.Figure 13
Transverse thermal conductivity validation with different volume fraction of fiber.It is clear from the figure that the transverse thermal conductivity increases with increases in fiber volume fraction. The transverse thermal conductivity evaluated by finite element analysis with hexagonal RVE is more close to the Hashin model as compared to the results obtained from rule of mixture and the Chawla model.


## 5.1. Effect of Volume Fraction on Elastic Properties

The longitudinal modulus is the response of composite during the application of load parallel to the fiber direction. It can be defined as the ratio of longitudinal stress to the longitudinal strain. Figure8 shows the effect of fiber content on the longitudinal modulus of composites using rule of mixtures, Halpin-Tsai model, periodic microstructure, and finite element analysis with square and hexagonal RVE with circular and square fiber geometry. It can be observed from the figure that the longitudinal modulus increases with the increase in volume fraction of fiber and there is a good agreement between finite element results and analytical methods. This is because the stiffness of the composite increases with increase in volume fraction of fiber.Figure 8
Longitudinal modulus validation with different volume fraction of fiber.The transverse modulus is the response of composite during the application of load perpendicular to the fiber direction. It can be defined as the ratio of transverse stress to the transverse strain. The effect of fiber volume fraction on transverse modulus of composites using finite element analysis and three analytical methods is graphically shown in Figure9. As expected, it is clear from the figure that the transverse modulus increases with increase in fiber volume fraction. The transverse modulus evaluated by finite element analysis with hexagonal RVE is more close to the Halpin-Tsai model and periodic microstructure as compared to the results obtained from rule of mixture and finite element analysis with square RVE.Figure 9
Transverse modulus validation with different volume fraction of fiber.Figure10 shows the effect of fiber volume fraction on the in-plane Poisson’s ratio of composite. It is evident from the figure that the major Poisson’s ratio decreases with increase in the volume fraction of fiber due to increase in material resistance. The finite element results are in good agreement with analytical methods.Figure 10
In-plane Poisson’s ratio validation with different volume fraction of fiber.In-plane shear modulus of composite is the ratio of shear stress to the shear strain in longitudinal direction. Figure11 shows the effect of fiber volume fraction on the in-plane shear modulus of composite. It is clear from the figure that the shear modulus increases with increases in fiber volume fraction. Also, it can be observed that there is a good agreement between results obtained from finite element analysis and hexagonal RVE with Halpin-Tsai model and periodic microstructure as compared to rule of mixture and finite element analysis with square RVE.Figure 11
In-plane shear modulus validation with different volume fraction of fiber.


## 5.2. Effect of Volume Fraction on Thermal Conductivity

Longitudinal thermal conductivity of composite is the property of a material to conduct heat in parallel to the direction of the fibers. Figure12 shows the effect of fiber content on the longitudinal thermal conductivity using rule of mixture, the Chawla model, the Hashin model, and finite element analysis. Finite element results match well the analytical methods and increases in linear manner with the increase in volume fraction since the fiber thermal conductivity dominates the matrix thermal conductivity.Figure 12
Longitudinal thermal conductivity validation with different volume fraction of fiber.Transverse thermal conductivity of composite is the property of a material to conduct heat, in the direction perpendicular to the fibers. Figure13 shows the effect of fiber content on the transverse thermal conductivity using analytical methods and finite element analysis.Figure 13
Transverse thermal conductivity validation with different volume fraction of fiber.It is clear from the figure that the transverse thermal conductivity increases with increases in fiber volume fraction. The transverse thermal conductivity evaluated by finite element analysis with hexagonal RVE is more close to the Hashin model as compared to the results obtained from rule of mixture and the Chawla model.


## 6. Conclusions

In this work, the micromechanical approach has been studied to evaluate the elastic and thermal properties of unidirectional fiber-reinforced polymer composites with different volume fractions of fiber up to 0.7. The material properties obtained by finite element techniques were compared with different analytical methods. The following conclusions can be drawn.(1)
The finite element model based on three-dimensional RVE with a square and hexagonal packing geometry was successfully implemented by using finite element code ANSYS to calculate elastic and thermal properties.(2)
It has been observed that longitudinal modulus and in-plane Poisson’s ratio predicted by the finite element analysis agree well with all the existing analytical predictions. Transverse modulus and in-plane shear modulus predicted by the finite element analysis with hexagonal RVE is more close to the periodic microstructure model as compared to finite element analysis with square RVE and other analytical methods.(3)
It has been observed that longitudinal thermal conductivity predicted by the finite element analysis agrees well with all analytical methods. Transverse thermal conductivity predicted by the finite element analysis with hexagonal array is more close to the Hashin model as compared to rule of mixture and Chawla model.(4)
In finite element analysis the elastic properties and thermal conductivity of fiber-reinforced composites are affected by a number of parameters such as cross section fiber, fiber geometry of RVE, and volume fraction of the fibers.



---

*Source: 629175-2014-11-18.xml*