Open access peer-reviewed chapter

FEM Analysis of Mechanical and Structural Properties of Long Fiber-Reinforced Composites

By Michal Petrů and Ondřej Novák

Submitted: May 23rd 2017Reviewed: October 23rd 2017Published: December 20th 2017

DOI: 10.5772/intechopen.71881

Downloaded: 690

Abstract

This chapter deals with studies of the mechanical properties of samples from long fiber-reinforced composite structures that would contribute to the optimization of the developed constructions made of them. First, the basic issues of composite structures reinforced with long fibers (carbon or glass) and generally of composites with the specification of parameters that would lead to the optimization of mechanical properties with respect to the theoretical strength are presented. Further, the possibilities and methods of measurements of composite reinforced with carbon and glass fibers are described. This is followed by the introduction of analytical models for the description of the transversal isotropic composite, where these mathematical relations allow the determination of unknown elastic constants and they are also important for the verification of numerical models. Finally, it is comprehensively outlined the problems of creating a numerical model of advanced composite fibrous structure for determining the mechanical properties, both through the description of the continuum, and complex numerical model with a structural configuration enabling approach to allow closer interaction among fibers and matrix. Compared to the averaged values obtained from experimental samples, numerical simulations show a similar trend of stress on strain, with results obtained from simulations.

Keywords

  • FEM
  • composite structures
  • testing
  • mechanical properties
  • nonlinear properties

1. Introduction

Studies and analyses of mechanical properties of long fiber-reinforced composites provide important information for future lightweight constructions. First of all, it is important to approach the issues and specifics of long fiber-reinforced composite structures to increase the strength and toughness of the resulting structure. The long fiber-reinforced composite structure is typically formed from two dominant components: carrier fiber reinforcement and a matrix. Ideal arrangement of the final composite (fiber-matrix connection), due to synergy, the high specific properties (high strength, stiffness, and toughness) can be achieved, where none of input components reached. It is that the optimal synergistic effect is characterized by a known “illogical” rule 2 + 3 = 7, which characterizes the sum of the properties of the individual input components (fibers + matrix) achieves a higher value of the specific properties of the newly created structure. In general, the highest specific properties can be achieved if the fibers are stressed up to the strength limit σMfF fmaxwith stress transferred with matrix. The matrix transforms the stress into the fibers, and it also has a significant effect on the bonding with the fibers. Thus, the matrix is the binder component of the composite, creating the final geometry of the composite and at the same time protecting the fibers from wear and damage, which would lead to loss of stability and strength of the resulting composite. The description of the properties of composite structures reinforced with long fibers due to their potential and specific characteristics was given by the authors namely Agarwal et al. [1], Guedes [2], Gay and Gambelin [3], Reifsnider [4], Teply and Reddy [5], Berthelot [6], Gibson [7], or Soden et al. [8]. The authors agree that long fiber-reinforced composite structures are unique materials whose mechanical properties cannot be generally described in an analytical or experimental manner. Theories also differ in mathematical relationships derived for unidirectional composite structures, let alone complete synthesis of mechanical properties for geometrically complex structures of frames with multidirectional fibrous arrangement. This is due to the fact that their properties vary significantly with the type of fibers and a matrix (e.g., physical and mechanical properties, surface treatment, chemical compositions, binding agents, density, thermal expansion, etc.) because only a slight change forms various combinations with significantly different properties in mechanical behavior.1 Generally, long fiber-reinforced composite structures can be considered as inhomogeneous and heterogeneous structures with anisotropic properties in terms of physical and mechanical behavior. Their heterogeneity is manifested by a large number of combinations of different variants of the resulting structural materials suitable or unsuitable for the specific design requirements and load.2 If the strength of the composite has to be maximized, the specific surface of the fiber-matrix interface must be high and free of defects. The selection of fiber reinforcement is possible to use a wide range of fibers, whereas their offer is developed and expanded. For structural applications such as frames for machine parts and equipment may be used virtually any organic natural fibers (e.g., coconut, cotton, cellulose fibers, etc.) from a variety of polycrystalline ceramic materials, polymeric fibers, glass, or carbon fibers. The production technology of these fibers is well described by Bareš [9]. Carbon fibers are industrially manufactured with a diameter of 5–12 μm by various methods such as carbonization of organic fibers or pyrolysis. It is generally known that carbon can exist in nature in three forms: diamond, graphite, and glassy (amorphous). Carbon fibers can be considered as fibers produced at 800–1600°C, and graphite fibers are produced at >2200°C. However, only fibers obtained from the crystalline form of carbon arranged in a certain direction (production under tension) have a high elastic modulus and other specific design parameters such as a lower density, higher surface area, lower thermal conductivity, higher electrical resistance, and so forth compared to graphite fibers. Glass fibers with diameter of 3.5–20 μm are produced by fast drawing from the melts (the speed reaches up to 400 m min−1). The spinning speed is also influenced by the viscosity (50–100 Pa s), the melting temperature, and, of course, the chemical composition of the glass. The matrix, which affects the properties and usability of the resulting composite, has been epoxy used for both the testing sample and design of the developed composite construction. The composite production may result in imperfect bonding of the matrix fibers (e.g., low wettability of the fiber reinforcement in the matrix, bubble formation, etc.), which leads in mechanical defects in the composite structure, which often over grow into critical defects with a significant reduction in strength (Figure 1). The resulting strength of the composite structure affects mechanical properties of the selected fiber reinforcement and matrix, which are characterized by mechanical parameters, for example, the elastic modulus, Poisson number, or other parameters such as the creep and fracture properties of the individual components. In terms of strength, a significant role (if not most) plays the interfaces among the fibers and the matrix, which is shown in Figure 2. This is due to the fact that the characteristic properties of the interface create a mechanism that apparently causes the synergistic effect that provides their unique mechanical properties to composite structures. Although a number of theories have been compiled, the synergistic mechanism of the phase interface is not yet clear.

Figure 1.

Example: Low strength of the fiber composite structure due to the poor joint of the fiber with matrix.

Figure 2.

Cross section of composite sample with long fibers (upper), detail: fiber-matrix (bottom left), detail: phase interface (bottom right).

2. Measurement of mechanical properties of composite samples with long fiber reinforcing

The determination of unknown parameters of composite materials has to be performed by experimental measurements. These parameters represent input data for numerical simulations. For a complete description of the properties, it is important to make measurements on both the fiber reinforcement (tow) and the matrix as well as on the resulting long fiber-reinforced composite structures (matrix-bonded fibers). Measurement of the mechanical properties of the samples is carried out according to standard laboratory tests, which are divided according to the time course of the applied load. Tests can be divided into static and dynamic. It can thus perform the tensile test at a constant or cyclic loading of the sample, three-point bending strength, and Charpy impact test, as shown in Figure 3. The samples may be formed in the “dog bone” shape or, optionally, in the form of a rectangle of defined length L, width b, and thickness h,3 whereas they can be used for short- or long-term test.

Figure 3.

Determination of mechanical properties of composite samples: (a) tensile test, (b) three-point bending test, and (c) Charpy impact test.

The characteristic physical properties of samples of long fiber-reinforced composites are influenced by weight and volume ratios of individual input components (fiber reinforcement and matrices) that ultimately affect design parameters (mechanical properties and weight of the structure). The mass and volume amounts of the fibers and the matrix in the composite structure sample can be defined according to the following relationships (1–5).

mc=mf+mmE1
Mf=mfmc,Mm=mmmcE2
Mf=1MmE3
Vf=vfvc,Vm=vmvcE4
Vf=1VmE5

where mcis the total weight of composite, mf,mmis the weight of fibers and matrix, Mf,Mmis the weight amount of fibers and matrix, Vf,Vmis the volume amount of fibers and matrix, vcis the total volume of composite, and vf,vmis the volume of fibers and matrix. Volume amount of fibers Vfand matrix Vmcan be also expressed with the help of the fiber density ρfand matrix density ρm, which is applied in Eq. (6). The total density ρc=mc/vccan then be expressed as the sum of components, that is, the reinforcement and matrix (7). The characteristic thickness of the composite structure can be expressed according to the relationship (8).

Vf=Mf/ρfMf/ρf+Mm/ρm,Mf=VfρfVfρf+VmρmE6
ρc=ρfVf+ρmVmE7
h=mf1ρf+1ρm1MfMfE8

Gay and Hoa [10] reported that winding of the fibers onto-shaped geometry may achieve the volume fraction of fibers in the composite maximally 55–80% of the total volume of the composite structure. Ideally, these values can be increased by precisely placing fiber tows side by side. A limit state of volume fiber fraction, that is, 100%, cannot be achieved due to the necessity of the presence of the matrix. Also by perfectly precise laying of fiber strands, the strands will always have a certain fill value that will never be equal to 1 in the geometric configuration. Perfectly precise laying of fiber tows does not provide 100% of volume filling due to fiber cross section. However, it should be noted that the optimum ratio of fiber reinforcement is in the range of the synergistic effect, that is, Vf0.4÷0.65.

The influence of selected physical parameters on the geometric parameter hof some tested samples from long fiber-reinforced composite structures is shown in Table 1. These parameters can then be used to establish numerical models. Other input parameters that are required for the numerical model have to be obtained by measuring the test samples. Other input parameters that are necessary for the numerical model are obtained by testing composite samples.

Name of fibersmf[g m−2]mm[g m−2]Mf[%]MV[%]Vf[%]Vm[%]ρf[g cm−3]ρm[g cm−3]h[mm]
GF 1600 tex/PUR Huntsman560600485230702.451.11.2
CF prepreg HEXPLY-M10R15091.96623852481.81.20.22
CF 24K/PUR Huntsman213747227815851.81.11.2

Table 1.

Examples of physical and geometrical parameters of the selected samples of composite structures reinforced with long fibers.

Note: mf∗,mmrepresent area weight, mcis the total weight of composite.

3. Analytical models for the study of mechanical properties of long fiber-reinforced composites

Numerical modeling of the mechanical properties of the composite is a very difficult problem because there are many unknown parameters that come into model simulations, which are discussed in this chapter. Therefore, some parameters need to be properly verified with analytical models. It is assumed that though mechanical properties of the sample are formed from uniformly spaced transversely isotropic structure, its theoretical description is difficult, as shown in Figure 4.

Figure 4.

Model of idealized transverse isotropic fiber composite structure.

The model of the transverse isotropic fiber composite structure can be defined by six independent elastic constants through the constitutive Eq. (9). The mechanical properties, such as composite structures, are also affected by the volume of fibers Vfand matrix Vm.

ε11ε22ε33γ12γ23γ13=1/E11ν12/E11ν12/E11000ν12/E111/E22ν23/E22000ν12/E11ν23/E221/E220000001/G230000001/G120000001/G12σ11σ22σ33τ12τ23τ13E9

where σii,εiiare the principal stresses and strains in the transversal isotropic composite in individual directions of the coordinate system x1,x2,x3, whereas σ11>σ22=σ33, ε11ε22=ε33, and τ12,τ23=τ13are the shear stresses in the given planes, γ12,γ23=γ13expresses the shear to individual planes, E11,E22=E33expresses the longitudinal and transverse modulus of elasticity, G12,G23is the shear modulus in the plane of the principal load direction and in a plane perpendicular to the principal load direction, ν12is the Poisson ratio in the principal direction of the load, and ν23is the Poisson ratio in a plane perpendicular to the principal load direction.

For the corresponding model, the interconnection of individual components A, B, C must be included (see Figure 2) to create a multiphase system approaching the behavior of composite structures. Therefore, the problem of modeling a composite can be treated as a continuum (a solid model without a geometric arrangement of individual components) or by creating a completely new model with structural parameters, that is, the individual components will be included in the structured unit. The problem of analytical modeling of mechanical properties of general fiber structures through a structural unit is described, among others, by Wyk for the study of interfiber contacts [11] and by Neckář [12]. However, the description of the mechanical properties of the fibrous composite structure is more difficult and has not yet been properly described. This is probably due to the fact that knowledge of the deformation mechanism and damage process is more important for understanding the mechanical properties than the knowledge of the absolute strength that cannot be determined with sufficient precision. This is due to the fact that it is not possible to comprehensively construct a general energy theory (to derive empirical relationships for deformation work) based on statistical characteristics, as can be done with very good accuracy for other anisotropic structures (Petrů et al. [13, 14]). The problem is that the individual components composing the composite structure cannot be reliably quantified even with homogeneous isotropic materials (matrices, glass fibers), let alone anisotropic structures such as carbon fibers (the theoretical value presented in the data sheet is different than value determined experimentally). Therefore, the main problems are related to the complexity of the description and modeling of deformation and the consequent character of the stress (stress concentration under loading). This is mainly due to technological influences in composite production (influence of temperature, humidity, and initial microporosity) that cannot be predicted for model simulations, and it is also relatively difficult to experimentally identify these parameters.4

Over time, there have been widespread analytical relationships to form the approach to obtain all elastic constants that can be used by these models, which are given as follows:

  • phenomenological models.

  • semiempirical models.

  • homogenized models.

3.1. Phenomenological models

In the past, phenomenological models have been created as the primary mathematical derivation of the mechanical properties of transversally isotropic fibrous composite structures but can be used well today. Such models include the Voigt and Reuss models. These are models using the mixing rule (mixing of the individual input components, i.e., fibers and matrices), while the Voigt model is very well usable for determining the elastic constants E11,ν12defined by relationships (10 and 11) and is characterized by isotropic strain. The Reuss model is usable for determination E22, G12defined by relationships (12 and 13) and unlike the Voigt model is characterized by isotropic stress.

dσ11dε11=Vfdσfdεf+VmdσmdεmE11=VfE11f+VmEmE10
ν12=Vfν12f+VmνmE11
E22=E22fEmEmVf+E22fEmE12
G12=G12fGmGmVf+G12fEmE13

where E11f,E22fare the longitudinal and transverse modulus of the fiber elasticity, G12fis a fiber shear module, and ν12fis the Poisson ratio in the plane of the principal direction of the load of the fiber.

3.2. Semiempirical models

Semiempirical models were created later than phenomenological models, and based on the new information and knowledge, they are still being updated. Their development led in particular to the further expansion of the Voigt and Reuss models because these models have been modified by correction factors to specify the resulting elastic constants for the given types of input components. This category includes models that are implemented in certain modifications in finite element softwares such as the Halpin-Tsai model or the Chamis model.

  • Modified model according to the mixing rule

Modified model according to the mixing rule is derived from Voigt [16] and Reuss [17], and for elastic constants, E11,ν12is defined according to Eqs. (10 and 11). Modification occurs with constants E22,G12,because the resulting difference between the results obtained by the measurements and the relationships (12–13) is usually noticeable. Therefore, it was necessary to make a correction for E22,G12according to Eqs. (14 and 15).

1E22=ζfVfE22f+ζmVmEmE14
1G12=VfG12f+ζ'VmGmVf+ζ'VmE15

where ζf,ζmare correction factors,5 according to Younes et al. [18].

  • Halpin-Tsai model

This is a model that is implemented in a number of numerical programs by using finite element method (FEM). This model is developed as a semiempirical model [19] with correction of E22,G12.Its semiempirical derivation (16–17) using correction factors ζ,ξhas a very good agreement with experiments.

E22=Em1+ξζVf1ζVfE16
G12=Gm1+ξζVf1ζVfE17

where ζcorrection factor, for which it is valid ζ=Mf/Mm1Mf/Mm+ξ, ξis constant, which is for E22equal to 1, and for G12is equal to 2, M=Eor Gin the case of an expression E22and G12according to Eqs. (1213).

  • Chamis model

This is another semiempirical model [20], which was unlike previous models developed not only for independent elastic constants E11,E22,G12,ν12but also for G23.The determination of E11,ν12is again based on Voigt and Reusse according to Eqs. (1011). The Chamis model for calculating other elastic constants introduces a square root of the volume of fiber Vf,which has in Eqs. (1820), the meaning of fiber incompressibility, which is in line with principle of mass conservation.

E22=Em1Vf1Em/E22fE18
G12=Gm1Vf1Gm/G12fE19
G23=Gm1Vf1Gm/G23fE20

where G23fis the shear modulus of the fiber elasticity in a plane perpendicular to the principle direction of loading.

3.3. Homogenized models

Homogenized models are generalized models that can be used to determine very accurate values of elastic constants for developed composite structures reinforced by longitudinally laid fibers. Such models include, for example, the Mori-Tanaka model [21], a consistent model created by Hill [22] or the Bridging model. Their applicability compared to phenomenological or semiempirical models largely limits the more difficult determination of all constants entering to homogenized models. An example is the Eshelby toughness tensor that can be used for inclusion, which is introduced in both the Mori-Tanaka model and the consistent model. In view of this, from homogenized models, the Brindling model can be used to determine the elastic constants.

  • Bridging model

This is a model that is developed to predict the stiffness and strength of transverse isotropic fiber composites. The elastic properties are for the elastic modulus E11,ν12the same as for Voigt and Reusse models (10–11). Elastic constants E22,G12,G23can be expressed using the Bridging model (21–23).

E22=Vf+Vma11Vf+Vma22Vf+Vma11VfS11f+Vma22S22m+VfVmS21mS21fa12E21
G12=Vf+Vma66G12fGmVfGm+Vma66G12fE22
G23=1/2Vf+Vma44VfS22fS23f+Vma44S22mS23mE23

where aii,Siif,mare the matrix components, which relate to fiber and matrix ratios in the composite structure as reported by Huang [23, 24], where a11=Em/E22f, a22=a44=0,35+1035Em/E22f,  a66=0,3+10305Em/(1+νm/G12f,   a12=S12fS12ma11a22/S11fS11m,S11f,m=1/E11f,m,S22f,m=1/E22f,mS12f,m=S21f,m=ν12f,m/E11f,m,S23f,m=S32f,m=ν23f,m/E22f,m.

4. Numerical models for the study of mechanical properties of long fiber-reinforced composites

Measurement and analytical models of long fiber-reinforced composite structures designed to study mechanical properties are generally able to provide only limited information. This is due to the fact that the measurements are limited by the possibilities of positioning of the sensors and also by the fact that some properties cannot be measured well (e.g., the distribution of the main stress and deformation in the composite structure). The knowledge of the distribution of the main stresses and deformations in the structure is important for assessing how the structure is changed and under which stress. In this case, the corresponding model simulation using numeric methods represents a significant support for the development. Very suitable is to build model simulation in finite element method (FEM), but other numerical methods, such as discrete element (DEM), boundary element (BEM) or finite volume method (FVM) method, are also available. The mechanical loading of composite causes many different processes in the inner structure that varies with the actual deformation. Therefore, it is necessary to simplify or neglect some characteristic features in modeling of such structures. A major problem of mechanical properties modeling of composite structures is in particular the description of the principal stresses in short time Δt=ti+1ti. The solution of problem of composite with boundary conditions under tensile loading lies not only in the specification of the correct boundary conditions and material properties but also in the design of the proposed finite element mesh. The FEM programs are currently very sophisticated and allow the solution of a continuous problem transform into a final solution where the corresponding geometric simple subareas (finite elements) can be designed in the preprocessor. Let Ω3is the continuous area of the three-dimensional space in which the problem is solved. Its borders will be denoted Γ, where Γis Lipschitz border and let the approximation of the selected base functions is derived above each finite element of length lelement,because any continuous function can be represented by a linear combination of algebraic polynomials converging to a continuous solution, that is, limlelement01.Thus, the FEM method can be understood as a special type of variation method by using the mathematical description of the problem solution. The current commercial software and FEM programs (e.g., Ansys, Abaqus, Permas, LS-Dyna, Marc, PAM CRASH) allow to assemble and subsequently solve a series of problems with nonlinear materials not only with elastic but also plastic behavior corresponding to the properties of the long fiber-reinforced composite.

4.1. FEM simulation of mechanical properties of long fiber-reinforced composite

Model simulations in FEM were performed for different combinations of reinforcement arrangements of long fiber-reinforced composites, which are important for comparison with experiments and analytical relationships. This gives the material properties for numerical simulation of the strength characteristics of whole frames.

This chapter describes the creation of two numerical models and their comparison:

  1. continuum model

  2. extended continuous model with structural unit

The simulations were performed for a complete assessment of the mechanical properties σ11,σ22,ε11,ε22,γ12,γ23and elastic constants E11,E22,G12,G23,ν12,ν23,whereas also information explaining the shape changes of the samples observed especially during tensile stress. Model simulations were performed in the following steps:

  • creating two model simulations of the long fiber-reinforced composite,

  • creating the corresponding mesh of finite elements of the computational model in the preprocessor,

  • defining the corresponding initial and boundary conditions,

  • assembling a material model of the long fiber-reinforced composite,

  • the evaluation and comparison of model simulation results in postprocessor.

4.1.1. Assembling of continuum of the long fiber-reinforced composite model

The FEM model was created in the concept of coherent continuum consisting of a surface geometry corresponding to the test sample with length L = 100 mm, width b = 20 mm, and thickness h = 1.7 mm. The finite element mesh of the numerical model was created from SHELL elements (2D elements) with a constant element size of 2 mm. The boundary conditions affecting the magnitude of the displacements and stress can be defined in two ways, that is, the boundary conditions of the first and second types. The first way is to determine the displacement and stress distributions if force conditions are known, that is, volume forces, surface forces, and nodal loads. The other way is to determine the displacement and tension distribution if the geometric conditions are known, that is, the size of node displacement, the deformations, and so on. Both ways can be also combined. It is mixed boundary conditions, as shown by Li [25]. Boundary and initial conditions for the model were made by the boundary conditions of the second type. One side of the sample was fixed against the displacement and rotation of nodes Ui=Ri=0in all directions (layout), and the opposite side of the edge of the sample was fixed identically; only in X direction, the movement was allowed, whereas Ux=1mmthat corresponds to the deformation εx=1%. The strain rate was 2 mm.min−1. The boundary conditions are shown in Figure 5.

Figure 5.

Continuum FEM model of the composite reinforced with long fibers with defined boundary and initial conditions.

4.1.2. Assembling of extended continuous composite model of long fiber-reinforced composite

The second numerical FEM model, which was created in the concept of structure unit, is formed from three components: fiber matrix—the interfacial interface, where the microscopic dimensions of such a model are closer to the more real composite. Such a model can be created from a structural unit with the 1,2,,nfibers, wherein the volume geometric configuration (e.g., structural unit is a cube, cuboid, sphere) can affect the volumetric quantity of fibers and a matrix Vf,Vm,as shown by Neckář [12]. The change in volume ratio Vfcan be given on the example of the structural unit of the cuboid, which is shown in the section in Figure 6. The structural unit consists of six fibers represented by circles with the same spacing mi, which are bounded by a matrix (rectangle). By changing of the spacing can then be increased or decreased volume ratio of fibers Vf.The finite element mesh of the numerical model was created from a combination of following elements: BRICK elements (3D elements) with a designed element size of 0.0002 mm defined for fibers and matrix (Figure 7).

Figure 6.

The influence of fiber spacing in the structural unit on V f fiber volume ratio.

Figure 7.

Structural FEM model of a composite reinforced with long fibers.

It will be assumed that E11=σ11ε11, E22=σ22ε22, G12=τ12γ12, G23=τ23γ23.

The problem lies in joining of fibers with the matrix because the interconnections form an interphase. The structural FEM model assembling presents a problem of the determination of appropriate boundary conditions, which is important in terms of accuracy and model verification. Incorrect design may result in concentrators and singularities of stress. The boundary conditions are created by the second type (geometric boundary conditions) as follows: the perimeter surfaces of the model perpendicular to the plane of the stretching direction have defined symmetry conditions on one side (symmetry in axis yand z) and on the opposite side, the boundary conditions are not prescribed. On surfaces in the plane of the stretch direction, that is, in the direction of the X axis, the displacements and rotations were not allowed Ui=Ri=0in all directions. On the opposite surface of the specimen, the condition was the same, only displacement in the stretching direction was allowed. The displacement was defined constantly to the maximum strain 1%, that is, Ux=k.εx=1%, k.=const., with strain rate 2 mm min−1. Boundary conditions are shown in Figure 8 and Table 2. The material properties applied in both FEM models (I. Continuum Model and the II. Continuum Model with the Structure Unit) are based on the generally known values reported by fiber and matrix manufacturers. The fiber and epoxy matrix parameters are listed in Table 3. The results of both numerical simulations have exhibited approximately the same stress at the defined strain εi=1%under tensile load in applied direction for a given fiber reinforcement (carbon or glass). The resulting dependence of force on the displacement of the samples obtained from the models showed an approximately linear course, both for carbon and glass fiber-reinforced composite. Figure 9 shows the tensile test for volume ratio Vf=0.3, where carbon fiber-reinforced composite with the epoxy matrix exhibits approximately 2.2 times higher force response than the glass fiber-reinforced composite with epoxy matrix.

Figure 8.

The boundary conditions of the structural FE model of the composite structure reinforced with long fibers.

Planes in axis xPlanes in axis yPlanes in axis z
+xx+yy+zz
E11,μ12Ux=k.εx=1%Ui,Ri=0Uy,Rz,Rx=0Uz,Rx,Ry=0
E22,μ23Ux,Ry,Rz=0Uy=k.εy=1%Ui,Ri=0Uz,Rx,Ry=0
G12Uy=k.εy=1%Uz=0Uy=0Uz=0Ux=0
Uz=0
Ux=0
Uz=0
Uz=0Uz=0
G23Ux=0Ux=0Ux=0
Uz=0
Ux=0
Uz=0
Uy=k.εy=1%
Ux=0
Ux=0
Uy=0

Table 2.

FEM model boundary conditions for obtaining elastic constants.

MaterialDensity
[kg m−3]
Modulus of elasticity
[GPa]
Shear module
[GPa]
Poisson’s ratio [−]Tensile strength
[GPa]
Elongation
[%]
E11f,mE22f,mG12f,mG23f,mν12f,mν23f,m
Carbon fibers1750 ± 15023015245.40.2790.492.3 ± 1.21.9 ± 0.6
Glass fibers2370 ± 23072.472.428.728.70.220.221.06 ± 0.654.8 ± 0.7
Epoxy matrix1150 ± 3703.5733.5731.311.310.3450.3450.067 ± 0.0333.6

Table 3.

Material and mechanical properties of individual constituent of composite.

Figure 9.

Comparison of experiments and FEM models: Dependence of applied force on sample strain.

The obtained results shown in Figures 10 and 11 can be stated that the continuous model (FE model I) has an approximately steady monotonic course manifested not only in continuous distribution of deformation (Figure 12 above) but also in uniform distribution of the principal stress σI=σ11acting in the load direction (Figure 12 center).

Figure 10.

Dependence of modulus G23 (left above), ν 12 (left below), ν 23 (right) on Vf of transversally isotropic composites with epoxy matrix and carbon fibers and glass fibers.

Figure 11.

Dependence of modulus E11, E22, G12 on Vf of transversally isotropic composites with epoxy matrix and carbon fibers (left) and glass fibers (right).

Figure 12.

Distribution of deformation (above) and principal stress in the loading direction (center) of the FEM model I of long fiber-reinforced composite with epoxy matrix. The resulting distribution of axial strain (below left) and the principal stress acting in the loading direction (below right) and the FEM model II of composite reinforced with long carbon fiber with epoxy matrix.

Due to the simplicity of the FEM model I, it appears to be very suitable for determining the mechanical properties of composite structures and their optimization. However, such a model does not provide information about strain and stress between the fibers and the matrix, let alone the interphase. The continuous model with the structural unit (FEM model II) is significantly more complex, and for some elastic constants, its resultant course is not monotone (G23,ν23); in other words, the resulting dependency is not stable and may not be accurate enough but more complex in terms of results. FEM model II allows to approximate the layout distribution of the structure unit in the loading direction (Z axis) as shown in Figure 12 (left) and also the principal stress distribution. The principal stress can be determined in isosurfaces or in sequential sections (Figure 12 right), which allow to analyze the stress distribution between the fibers and the matrix including the interphase. By comparing the maximum values of the stress of 189.1 and 190.9 MPa at same strain ε11=1%and with the volume ratio Vf=0.5can be stated that the models have significant agreement. This is affected not only by the same material parameters and boundary conditions but also appropriately selected types of elements of the finite element mesh as discussed earlier. FEM model II of the composite structure reinforced with longitudinal fibers with the epoxy matrix is more complex, and the time of the calculation is larger than the FEM model I. However, it must be added that the FEM model II shows valuable scientific knowledge of the approximate distribution of the maximum stress between the fibers and the matrix, which is the largest in the interphase (Figure 13 left). This confirms the theoretical assumption of the system mechanism (fiber-matrix interphase), where the highest stress transmits newly created component, that is, the interphase, which causes the synergic effect of the resultant construction of the composite structure. Figure 13 also shows the information that FEM model II shows a nonuniform maximum stress in the composite structure (unlike the FEM model I) and also shows that maximum stress is concentrated only on two fibers (instead of six in the structural unit) of the structural unit. It will reduce the maximum synergic effect that theoretically in the composite structure can occur. The distribution of the interphase in the numerical model and in the real measurement is shown in Figure 13. It is necessary to add that from the analyses carried out by measurements on real samples was in all cases evident that the identification of the interphase is very complicated. Due to the unidentifiable technological process, interphase (third component) could not be created. Also, it is problem to identify and measure the interphase that is important both for verifying of numerical models and for a statistical evaluation how many fibers are involved in synergistic effect.

Figure 13.

Distribution of principal stress in direction of loading with maximal stress in interphase (left), real sample with visible interphase (right).

5. Conclusion

In this chapter, analyses and numerical simulations of mechanical properties of samples from composite structures were presented. Several studies and experiments have been carried out on samples reinforced with carbon and glass fibers, and mechanical properties allow them to form structural reinforcements of composite materials. Further, analytical models with mathematical relationships (e.g., Voigt, Reuss or Chamis model) allow to determine the unknown elastic constants E11,E22,G12,G23,ν12,ν23of the resulting composite structure. This is followed by a more extensive description of the creation of a numerical model of a composite fiber structure pattern for determining mechanical properties, both through the description of a general continuum and a more complex numerical model with a structural arrangement to allow closer interaction of the fiber and the matrix. From the numerical models, the stress and strain distribution can be determined over a given time interval under chosen packing density Vfas well as the elastic constants. The course of elastic constants has to be compared in some cases only with analytical models because unknown constants cannot be appropriately measured. In summary, the I. continuous model is more user-friendly for numerical simulation and that is suitable for describing the principal stresses, but it does not allow to analyze and study the composite on microlevel. Thus, it does not allow the distribution of the stress between the fiber and the matrix (interphase). This can be done with more complex II. extended continuous model with a structural unit. The results of numerical models establish valuable knowledge and information, including the determination of elastic constants for a particular specific composite design. These results can be used for modeling of large samples and complicated geometries to optimize the design solution.

Notes

  • This can be mentioned in the example given by Bareš [9]. A simple combination of three homogeneous isotropic light metals to form ternary cast iron is obtained 82,160 possible variants of alloys, if more six metals are combined, more than 300 million different alloy variants could be obtained. (The composite structure reinforced with long fibers has a similar behavior, where the change of the matrix, directional arrangement, and type of fiber significantly affect resulting mechanical properties because it leads to qualitatively different structure [10]).
  • Transmission of static stress applied to the composite and transferred with fibers is required excellent consistency of fibers and matrix; on the contrary, a dynamic impact requires energy absorption by the crack propagation along fibers.
  • Note: geometrical dimensions h , b , L can be also smaller, but it can lead to a problem with clamping of the sample to the jaws of dynamometer.
  • In the advanced model, simulations can be assembled material models with any parameters, including the statistical parameters, which describe technological production factors, for example, through the theory of random fields as defined by Bittnar and Šejnoha [15]. The problem lies in the identification of the effects and the subsequent statistical evaluation.
  • Correction factors ζ f , ζ m can be express as ζ f = E 11 f V f + 1 − ν 12 f ν 21 f E m + ν m ν 21 f E 11 f V m E 11 f V f + E m V m , ζ m = E m V m + 1 − ν m ν m E 11 f − 1 − ν m ν 12 f E m V f E 11 f V f + E m V m , ζ ' is a variable function 0 < ζ ' < 1 , whereas preferred is ζ ' ≈ 0.5 − 0.6 .

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Michal Petrů and Ondřej Novák (December 20th 2017). FEM Analysis of Mechanical and Structural Properties of Long Fiber-Reinforced Composites, Finite Element Method - Simulation, Numerical Analysis and Solution Techniques, Răzvan Păcurar, IntechOpen, DOI: 10.5772/intechopen.71881. Available from:

chapter statistics

690total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Finite Element Dynamic Analysis on Residual Stress Distribution of Titanium Alloy and Titanium Matrix Composite after Shot Peening Treatment

By Lechun Xie, Zhou Wang, Chengxi Wang, Yan Wen, Liqiang Wang, Chuanhai Jiang, Weijie Lu, Lai-Chang Zhang and Lin Hua

Related Book

First chapter

Modeling the Physical Phenomena Involved by Laser Beam – Substance Interaction

By Marian Pearsica, Stefan Nedelcu, Cristian-George Constantinescu, Constantin Strimbu, Marius Benta and Catalin Mihai

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us