Next Article in Journal
Evaluation of a Hydrophobic Coating Agent Based on Cellulose Nanofiber and Alkyl Ketone Dimer
Next Article in Special Issue
The Development and Study of a New Silylated Polyurethane-Based Flexible Adhesive—Part 2: Joint Testing and Numerical Modelling
Previous Article in Journal
Development of Pore Pressure in Cementitious Materials under Low Thermal Effects: Evidence from Optimization of Pore Structure by Incorporation of Fly Ash
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy Equivalence Based Estimation of Hybrid Composites Mechanical Properties

Department of Applied Mechanics and Biomechanics, Faculty of Mechanical Engineering, Tadeusz Kosciuszko Cracow University of Technology, 37 Jana Pawla II Av., 31-864 Cracow, Poland
*
Author to whom correspondence should be addressed.
Materials 2023, 16(12), 4215; https://doi.org/10.3390/ma16124215
Submission received: 28 March 2023 / Revised: 1 June 2023 / Accepted: 2 June 2023 / Published: 6 June 2023
(This article belongs to the Special Issue Research on Material Durability and Mechanical Properties)

Abstract

:
Hybrid composites, usually combining natural and synthetic reinforcing filaments, have gained a lot of attention due to their better properties than traditional two-component materials. For structural applications of hybrid composites, there is a need to precisely determine their mechanical properties on the basis of the mechanical properties, volume fractions, and geometrical distributions of constituent materials. The most common methods, such as the rule of mixture, are inaccurate. More advanced methods, giving better results in the case of classic composites, are difficult to apply in the case of several types of reinforcement. In the present research, a new estimation method is considered, which is simple and accurate. The approach is based on the definition of two configurations: the real, heterogeneous, multi-phase hybrid composite configuration, and the fictitious, quasi-homogeneous one, in which the inclusions are “smeared out” over a representative volume. A hypothesis of the internal strain energy equivalence between the two configurations is formulated. The effect of reinforcing inclusions on the mechanical properties of a matrix material is expressed by functions of constituent properties, their volume fractions, and geometrical distribution. The analytical formulas are derived for an isotropic case of a hybrid composite reinforced with randomly distributed particles. The validation of the proposed approach is performed by comparing the estimated hybrid composite properties with the results of other methods, and with experimental data available in the literature. It is shown that a very good agreement is obtained between experimentally measured hybrid composite properties and their predictions resulting from the proposed estimation method. The estimation errors are much lower than the errors of other methods.

1. Introduction

Hybrid composites are materials that consist of a matrix and at least two different types of reinforcement. They can be manufactured from a combination of a variety of organic or inorganic materials such as metals, ceramics, polymers, etc., to obtain the desired mechanical, thermal, electrical, or other properties. Virtually unlimited possibilities of combining several organic or inorganic constituents allow the designing of a huge number of new products with high precision, as well as modifying conventional composite materials.
Hybrid composites usually exhibit properties that neither organic nor inorganic constituent materials can reach separately; they may be flexible and at the same time be characterized by high mechanical strength and thermal stability. Moreover, the resulting properties of the combination of several reinforcing materials may exhibit a synergistic effect [1]. For example, composites often have better fatigue properties than metals [2]. Although initial fatigue damage may appear earlier in composites than in metals, further damage development may slow down due to arresting microcracks in the internal structure of a composite, while in metals cracks usually rapidly develop [3].
Increasing demand for eco-friendly materials encourages implementing agricultural waste as natural fibers, increasingly used as an addition or replacement to traditional synthetic fibers, which have harmful effects on the environment. For this reason, the hybrid use of natural and synthetic reinforcing materials has gained a considerable attention from researchers [4,5]. Generally, synthetic fibers exhibit better mechanical properties, while natural fibers, despite their weaker mechanical properties, may have a positive effect e.g., on reducing the weight of the composite, which is a feature desirable especially for lightweight structures. One of the sectors for which weight reduction is important is the automotive industry. For example, the combination of kevlar and date palm fibers gives the composite better enhancement in tensile strength and tensile modulus. Such dependence has been observed by Muthalagu et al. [6] who studied the mechanical properties of kevlar fiber/date palm fiber/epoxy resin hybrid composite in the context of the application for passenger car bumper beams. Composites with different weight percentages of reinforcements were fabricated by hand lay-up method and tested in a tensile test. A bumper model was also performed by CATIA and subjected to a finite element analysis in ANSYS. Ganesarajan et al. [7] provide a comprehensive study on the development of hybrid composites containing engineered poly-saccharide to deliver optimum mechanical properties suited for body interior and under-the-hood applications for passenger vehicles or light-duty trucks. The proposed solutions exhibit progress in light-weighting, manufacturing cycle time, and overall material sustainability. In paper [8] natural rubber toughened silica/kenaf fiber/epoxy resin composites are investigated for possible automotive applications. The addition of liquid rubber MG30 produced rubber particles in the composite helped to increase the overall mechanical properties. Mansor et al. [9] used the Analytical Hierarchy Process (AHP) method to select the most proper natural fiber for hybridization with glass fiber reinforced polymer composites for the design of a passenger vehicle brake lever. Of the 13 fiber types tested, kenaf fiber turned out to be the best option which fulfills performance requirements and the design objectives. Dattatreya et al. [10] studied the mechanical properties of hybrid composites made of waste natural fibers or fillers and epoxy resin. The research included composites with jute fiber, coconut coir, sugarcane, and neem wood powder. Due to its high tensile properties, jute fiber was chosen as the base material and the other three fillers of hybrid composites had different weight percentages. The composite composed of 15% neem wood waste, 5% coconut fiber, and 5% sugarcane bagasse achieved the best results in terms of flexural, tensile, and impact strength. Which makes these composites potentially useful in the automotive industry.
Considering other fields, Bagheri et al. [11] investigated carbon fiber/flax fiber/epoxy resin hybrid composite for bone implants. Compared to commonly used metal plates for long bone fracture treatment, the mechanical properties of the hybrid composite were closer to human cortical bone. The presented results suggest that this material may be considered for the use in long bone fracture fixation. Scutaru and Baba [12] studied the impact properties of carbon fiber/hemp fiber/polyester resin laminates. They concluded a high stiffness of investigated hybrid composite. Additionally, carbon fiber/hemp fiber composite laminates have good thermal properties. Due to this, they are used e.g., for insulators in spaceships, airplanes, and nuclear reactors. Cicala et al. [13] tested the possibility of using glass fiber/natural fiber/epoxy resin hybrid composites for applications in the piping industry. Hemp, kenaf, and flax fiber were used. Manufactured hybrid composites allowed for a reduction in cost by 20% and a reduction in weight by 23% compared to the current commercial solutions. Hybrid composites are also used on boat components such as propellers and masts. Research on one of the first prototypes of a glass fiber/carbon fiber hybrid composite mast revealed improvements in ballistic performance as well as stiffness. A weight reduction of 20–50% compared to an aluminum shaft of the same size has also been reported [14].
Panthapulakkal and Sain [15] investigated water absorption and mechanical and thermal properties of glass fiber/hemp fiber/polypropylene hybrid composites. Results show that the addition of glass fibers to hemp fiber/polypropylene composite improved all studied material properties and that such a hybrid material can be successfully used for structural applications where thermal resistance and high stiffness are desired.
Following the expanding application of hybrid composites, it is necessary to extend the use of classical averaging methods for two-component composite properties to hybrid composites with a larger number (in general: n) of reinforcing materials. The simplest classical models for the estimation of the average elastic moduli of a two-component composite are the Voigt [16] and Reuss [17] mixture rules based on volume fractions of component materials, and are proved (cf Hill [18]) to constitute respectively the upper and lower bounds of the real effective stiffness values. However, neither the Reuss nor the Voigt assumptions are exact. For the Voigt approach, the tractions and phase boundaries may not be in equilibrium, while for the implied Reuss strains the matrix and inclusions may not remain bonded. Additionally, the mentioned bounds are usually distant from each other, therefore there is a wide range of possible predicted values of the effective elastic properties. Hashin and Shtrikman [19] used the concept of polarization and the principle of minimum potential energy to provide narrower bounds for isotropic composites. Another popular homogenization method named the Self-Consistent Scheme (SCS) [20] assumes that the single particle is embedded in an effective medium of unknown properties. This method provides the implicit expressions for the effective shear and bulk moduli. A different concept is due to Nazarenko et al. [21], who considered the effective properties of composites with the interphases characterized by the surface-varying properties. Many classical theories and homogenization methods account both for the topology of reinforcement and the material symmetry of a constitutive model describing a composite (cf. for ex. Sun and Vaidya [22], Gan et al. [23], Liu et al. [24], Würkner et al. [25], Selvadurai and Nikopour [26] and others). Generally, it is assumed that the distribution of component materials has a periodic manner, therefore it is possible to assume the existence of the so-called representative unit cell (RUC).
Predicting the effective properties of hybrid composites becomes more complicated with a growing number of components, and so far, there has been little theoretical modeling of property relationships in these materials. In [27] the classical lamination approach to predict the theoretical elastic modulus of a natural fiber-reinforced hybrid composite is used. It is shown that the theoretical values of elastic and flexural moduli were not accurately estimated by the hybrid mixture method. The application of the Mori–Tanaka method [28] or self-consistent method to the solution of this problem gives overall elastic moduli tensors of such composites that do not exhibit the necessary symmetry [29]. The hybrid graphene/carbon nanotube (CNT)/elastomer nanocomposite is considered in [30] as a pressure sensor under bending. Both graphene and CNTs are assumed to be randomly distributed in the composite and treated as ellipsoidal inclusions in an isotropic matrix. The effective elastic modulus was estimated by the use of the Mori–Tanaka method and compared to the experimental data.
In the majority of research papers, multi-step homogenization is applied. Hasanzadeh et al. [31] propose a multi-procedure micromechanics approach based on the Mori–Tanaka model to obtain the properties of piezoelectric hybrid composites containing carbon nanotubes. In the first step, the elastic properties of a nanocomposite consisting of randomly distributed CNTs in the polymer matrix are modeled. In the second step, considering the nanocomposite as a matrix and piezoelectric fiber as reinforcement, the elastic properties of CNT–piezoelectric fiber-reinforced hybrid composites are predicted. Bazan et al. [32] proposed a two-step homogenization procedure to predict the effective Young modulus of the three-phase short-fiber-reinforced composite. In the first step, the estimation of the effective elastic modulus is performed using the Halpin-Tsai [33] approximation, applied to the matrix and one of the reinforcing components. In the second step, the property estimated in the first step is used as the “matrix” modulus and homogenized with the second reinforcing phase. In paper [34] a multiscale finite element technique and a sequential methodology are used to estimate the mechanical properties of hybrid composites made of a polymer matrix, reinforced with carbon fibers and silica microparticles. The results show a very good agreement with the experimental data.
Although the classical Voigt’s rule of mixture was shown to result in significant discrepancies with the experimental results, especially in the case of natural fibers [35], it can still be found in use due to its simplicity [36,37].
The approach to predict the elastic properties of a hybrid composite presented in this research utilizes a concept of a fictitious quasi-homogeneous continuum [38], in which the inclusions are “smeared out” within a representative volume element (RVE). The hypothesis of equivalence between the fictitious quasi-homogeneous continuum and the real composite material configuration is then postulated. The approach is analogous to the formalism of the continuum damage mechanics (CDM), which deals with materials containing microdefects that degrade the material properties. The subject of the current approach is a material with inclusions that improve its properties, which is the opposite effect; nevertheless, it can be properly reflected with the use of the same thermodynamic formalism.
The first attempts to extend the CDM-based approach to a general case of thermodynamically based constitutive modeling were performed by Egner and Ryś [39,40,41]. The idea to apply such reasoning to the composite material properties estimation was started in the paper by Wiśniewska et al. [42], where the basic concept of mapping between a real composite and a fictitious quasi-homogeneous configuration was formulated and used to estimate elastic properties of isotropic two-component composite materials. The approach was furtherly developed by Wiśniewska and Egner in [38]. Different forms of an inclusion-effect tensor were discussed in terms of its influence on the stress and strain states in a fictitious quasi-homogeneous configuration. Additionally, the concentration factors were derived and compared with other homogenization methods, and an optimization problem of an FGM shaft was analyzed with the use of the proposed estimation method. The extension of the method for inelastic material properties of two-component isotropic composites was presented in [43].
In the present paper, this concept is further developed and applied for the first time to hybrid composite materials. A significant advantage of the method used in the presented research in comparison with the above-mentioned classical approaches is the simplicity of calculations and the possibility to predict with significantly better accuracy the effective properties of composite materials. What is more, there is no need to incorporate any multi-step procedures.

2. Theoretical Formulation

2.1. Basic Ideas

The main assumption adopted in the considered homogenization method states that the effective properties of an isotropic composite can be described by volume averaging of the actual properties of its components (matrix and reinforcing inclusions) over a domain called the representative volume element (RVE) which maps a finite volume of the real heterogeneous hybrid composite material into a material point of the fictitious pseudo-homogeneous configuration. This effective quasi-homogeneous configuration method is based on the following assumptions:
  • Each inclusion within the RVE matrix is subjected to the same stress field derived from the external tractions applied at the boundary of this element.
  • The effect of all other inclusions within the RVE matrix on the observed inclusion is measured through the change of effective composite stiffness or compliance.
  • For the isotropic short fiber-reinforced composites the exact spatial correlation of the inclusions within the RVE has negligible influence on the effective composite properties defined within the element.
A transition between the microscale and the effective composite properties on the macroscale requires a proper selection of the RVE size, bounded by two constraints. On one hand, it must be large enough to include a sufficient number of reinforcing inclusions that justify homogenization within the RVE. However, on the other hand, the size of the RVE must be small enough for the stress state to be considered homogeneous. The existence of the RVE, which allows the material to be considered statistically homogeneous within the volume element, is the condition for the presented approach, in which there are no scale parameters involved. The true distribution of the inclusions within the RVE matrix, described by the (scalar or tensorial) parameter ξ , results in the change of the macroscopic constitutive properties of a composite, which may be estimated on the basis of the mechanical equivalence hypothesis. From among several physical quantities (strain, stress, elastic strain energy, etc.), the equivalence of the total strain energy between the real composite configuration and the fictitious quasi-homogeneous configuration is here assumed, which allows determining not only elastic mechanical properties but also inelastic parameters [43].
The total strain energy equivalence (TEE) hypothesis can be described as follows [38,39,42,43,44]:
At any time, to an RVE in its real (multiphase, hybrid) configuration, described by the set of state variable pairs V α (stresses, strains), we associate a homogeneous (monophase) equivalent fictive configuration, the state of which is described by the effective state variables V ˜ α (effective stresses, effective strains)—in such a manner that the total internal strain energy U defined over the two (real and fictive) configurations is the same:
U V α , ξ = U V ˜ α , 0  
Accordingly, the definitions of effective stresses σ ˜ i j and effective elastic strains ε ˜ i j e are related to the two above mentioned configurations: real (R) and fictitious (F). Assuming a linear relation between two second-order symmetric tensors, the relation between the variables V α and V ˜ α may be expressed in the most general form with the use of the fourth-order tensor N i j k l e , called here the inclusion-effect tensor. This tensor allows mapping from heterogeneous (R) to quasi-homogeneous (F) configuration (cf [45,46,47], see Figure 1):
σ ˜ i j = N i j k l e ξ 1 σ k l ,   ε ˜ i j e = N i j k l e ξ T ε k l e
In the simplest case of an isotropic composite with randomly distributed and oriented short reinforcing fibers, the inclusion-effect tensor was proposed in [42] in the following isotropic two-parameter form:
N i j k l e ξ = f 1 ξ δ i j δ k l + f 2 ξ δ i k δ j l + δ i l δ j k ,
where f 1 ξ , f 2 ξ are scalar functions of the parameter ξ , defined in the isotropic case as the reinforcement volume fraction d V I in the total volume d V R V E of the RVE:
ξ = d V I d V R V E
By the use of the total strain energy equivalence hypothesis, the elastic stiffness tensor E i j k l ξ of a composite material results from the equivalence between elastic strain energies:
U e = 1 2 σ i j ε i j e = 1 2 σ ˜ i j ε ˜ i j e ,
Taking into account Equations (2) and (5) and the Hooke law, the elastic stiffness tensor E i j k l ξ takes the following form, dependent on the matrix stiffness tensor E i j k l M ξ and the inclusions volume fraction ξ [38,42,43]:
E i j k l ξ = N i j p q ξ E p q r s M N r s k l ξ ,
The above equation can be transformed with the use of Equation (3) into the following form:
E i j k l ξ = λ ξ δ i j δ k l + μ ξ δ i k δ j l + δ i l δ j k ,
where λ ξ and μ ξ are the effective Lamé constants dependent on the matrix constants λ M , μ M , and the volume fraction of inclusions:
λ ξ = λ M 9 f 1 2 + 12 f 1 f 2 + 4 f 2 2 + 2 μ M 3 f 1 2 + 4 f 1 f 2 ,
μ ξ = 4 μ M f 2 2 .
The details of the above transformations can be found in [42], and will not be repeated here.
The elastic properties estimation based on the hypothesis of total strain energy equivalence allows obtaining a very good agreement with the experimental data, while the calculations are relatively simple. Figure 2 shows several examples of composite material properties estimations compared with the experimental data. In this graph, solid lines represent total strain energy equivalence-based (TEE) estimations, while dots stand for the experimental results of the Young modulus for carbon short-fiber reinforced polyacetal (POM-CF), glass short-fiber reinforced polyamide (PA-GF) [48], and hydroxyapatite reinforced polyethylene (PE-HAp) [49,50].
There are some small discrepancies for the POM-CF composite, while for the other two materials, the predictions are excellent. It should be pointed out that the method proposed here, like all other homogenization methods, is based on simplifying assumptions and therefore approximate (for example, the matrix/inclusion interphase interaction is disregarded). Nevertheless, the errors in predicted values are very small in comparison with other classical methods. As it can be seen in Table 1, for more than 60% of results the difference between experimental and predicted Young modulus values, related to the experimental value (error), is the smallest for the TEE method proposed in the present research (the smallest errors are highlighted and underlined).
The above-described concept will be further extended in the next section for the case of hybrid composites.

2.2. Effective Elastic Properties of Hybrid Composite Material

In a general case of a hybrid composite, the volume fraction of reinforcement (4) can be decomposed as follows:
ξ = d V I 1 + + d V I n d V R V E = i = 1 n ξ i .
In engineering practice, the most common are hybrid composites made of two types of reinforcement I 1 , I 2 and one matrix material M . For such a three-phase isotropic composite with randomly oriented and distributed inclusions, the total fraction ξ of the reinforcement volume d V I in the total volume d V R V E of the RVE takes the form:
ξ = d V I d V R V E = d V I 1 + d V I 2 d V R V E = ξ 1 + ξ 2 .
In the case of a three-phase hybrid composite, the unknown functions f 1 ξ 1 , ξ 2 and f 2 ξ 1 , ξ 2 appearing in Equation (3) are therefore the functions of two variables, namely the volume fractions of both reinforcing materials ξ 1 and ξ 2 . These functions must satisfy the boundary conditions at the two characteristic points:
(1)
ξ = ξ 1 = ξ 2 = 0 when a homogeneous pure matrix material is considered, therefore the stiffness tensor (7) becomes the stiffness tensor of the matrix material:
ξ 1 = ξ 2 = ξ = 0 E i j k l 0 , 0 = E i j k l M λ 0 , 0 = λ M μ 0 , 0 = μ M f 1 0 , 0 = 0 f 2 0 , 0 = 1 2 ,
(2)
ξ i = 1 when the volume fraction of one of the reinforcement materials is equal to 1 (a homogeneous pure inclusion material is then considered). This second condition involves two cases: ξ 1 = 1 ξ 2 = 0 , and ξ 1 = 0 ξ 2 = 1 :
ξ i = 1 ξ = 1 ξ 1 = 1 ξ 2 = 0 o r ξ 1 = 0 ξ 2 = 1 E i j k l 1 , 0 = E i j k l I 1 o r E i j k l 0 , 1 = E i j k l I 2 λ 1 , 0 = λ I 1 μ 1 , 0 = μ I 1 λ 0 , 1 = λ I 2 μ 0 , 1 = μ I 2 .
In the above expression ( λ I 1 , μ I 1 ) and ( λ I 2 , μ I 2 ) denote pairs of Lamé constants for the first and second inclusion material, respectively.
Taking into account expressions (7)–(9), the second boundary condition, (13), results in the following solutions for functions f 1 ξ 1 , ξ 2 and f 2 ξ 1 , ξ 2 :
f 1 ξ i = 1 , 0 = 1 3 μ I i μ M K I i K M ,   f 2 ξ i = 1 , 0 = 1 2 μ I i μ M ,   i = 1 , 2 .
where K M , K I 1 , and K I 2 denote respectively the bulk moduli for matrix, inclusion-1, and inclusion-2 materials.
The first boundary condition therefore gives a single point solution:
f 1 0 , 0 = 0 f 2 0 , 0 = 1 2 .
whereas two points result from the second boundary condition:
f 1 1 , 0 = 1 3 μ I 1 μ M K I 1 K M f 2 1 , 0 = 1 2 μ I 1 μ M ,
and
f 1 0 , 1 = 1 3 μ I 2 μ M K I 2 K M f 2 0 , 1 = 1 2 μ I 2 μ M .
The simplest linear approximation of the functions f 1 ξ 1 , ξ 2 and f 2 ξ 1 , ξ 2 between the boundary points (15), (16), and (17) has been here proposed:
f 1 ξ 1 , ξ 2 = 1 3 μ I 1 μ M K I 1 K M ξ 1 1 3 μ I 2 μ M K I 2 K M ξ 2 ,
f 2 ξ 1 , ξ 2 = 1 2 μ I 1 μ M 1 ξ 1 + 1 2 μ I 2 μ M 1 ξ 2 + 1 2 .
The linear approximation between the boundary points together with the total strain energy equivalence hypothesis result in the second-degree polynomial description of the effective Lamé constants λ ξ 1 , ξ 2 and μ ξ 1 , ξ 2 for the three-component hybrid composite. Taking into account (8)–(9) and (18)–(19) the following relations are finally obtained:
λ ξ 1 , ξ 2 = λ M C 1 ξ 1 + C 2 ξ 2 2 2 C 1 ξ 1 + C 2 ξ 2 D 1 ξ 1 + D 2 ξ 2 + 1 + D 1 ξ 1 + D 2 ξ 2 + 1 2 + 2 3 μ M C 1 ξ 1 + C 2 ξ 2 2 2 C 1 ξ 1 + C 2 ξ 2 D 1 ξ 1 + D 2 ξ 2 + 1 ,
μ ξ 1 , ξ 2 = μ M D 1 ξ 1 + D 2 ξ 2 + 1 2 .
where
C i = μ I i μ M K I i K M ,                 i = 1 , 2     D i = μ I i μ M 1 ,                               i = 1 , 2
The approach may be furtherly developed by replacing in (18)–(19) the linear functions with more advanced interpolations.

3. Results

3.1. Parametric Studies

It should be underlined that the proposed method does not involve any model parameters that need to be identified experimentally. Parametric studies in this case concern therefore only the analysis of the model response (prediction of effective composite properties) for different volume fractions of the fillers. To perform such parametric studies the hybrid composite material described in [32] was selected, which is composed of a polyamide matrix (PA) with the addition of aramid fibers (AF) and basalt fibers (BF). Aramid fibers are a kind of synthetic organic fibers with high fatigue and creep resistance. Basalt fibers are inorganic reinforcement characterized by high strength and stiffness. They have good chemical and high temperature resistance. The advantage is also their non-toxicity and environmental friendliness. Aramid fibers are here denoted as reinforcement material number 1 ( I 1 ), while basalt fibers are named as filler number 2 ( I 2 ). The Lamé constants and the bulk moduli for all the composite components are presented in Table 2.
Figure 3a shows the results assuming a constant volume fraction of the aramid fibers ξ 1 and a variable volume fraction of the basalt fibers ξ 2 . Dashed lines indicate the limit cases when the volume fraction of one of the composite components (either matrix or one of the filaments) is zero, therefore a two-component material is then considered: polyamide/basalt at the lower curve, basalt/aramid at the upper line, and polyamide/aramid at the vertical dashed line. Also, at three characteristic points pure homogeneous material properties are properly reflected: ξ 1 = ξ 2 = 0 (no fillers, pure polyamide matrix), ξ 1 = 0 and ξ 2 = 1 (no aramid fibers, no polyamide matrix, pure basalt), ξ 1 = 1 and ξ 2 = 0 (no basalt fibres, no polyamid matrix, pure aramid). Analogical simulations were performed for the constant volume fraction of the second reinforcement (basalt), while the content of the first one (aramid) is variable (Figure 3b). Similar conclusions may be drawn as for Figure 3a.
The obtained results were next compared with the Voigt and Reuss methods. Figure 4 shows the distribution of the effective Young modulus assuming a constant volume fraction of the aramid reinforcement ξ 1 and a variable volume fraction of the basalt reinforcement ξ 2 . All the total strain energy equivalence-based predictions (TEE) are correctly placed between the lower (Reuss) and upper (Voigt) bounds.

3.2. Validation

Model validation was carried out for hybrid composites described in the available literature. The results for the (PA-AF-BF) composite described in the previous Section 3.1 will be discussed first. Table 3 contains the Young modulus and Poisson ratio for the components of this composite (see also Table 2).
The space distribution of the Young modulus E ξ 1 , ξ 2 as a function of two variables is shown in Figure 5. The 3D surface predicted by the proposed method (TEE) is shown juxtaposed with the Voigt and Reuss surfaces, and the available experimental data. In the enlarged fragment of the drawing, a very good agreement can be seen between the TEE predictions and the experimental data. Similar results are presented in Figure 6, where the experimental data are placed along with the lines created by the intersection of the surface E ξ 1 , ξ 2 with planes ξ 2 = c o n s t (variable volume fraction of aramid, constant volume fraction of basalt, Figure 6a) and ξ 1 = c o n s t (variable volume fraction of basalt, constant volume fraction of aramid, Figure 6b).
The same validation was performed with the use of another hybrid composite material available in the literature, which is a polylactide (PLA) reinforced with carbon (CF) and basalt (BF) fibers [51]—see Figure 6c,d. In this case, the hybridization was aimed at producing a composite with comparable strength and stiffness as two-component PLA-BF composites, but with lower density, thanks to the synergistic combination of carbon and basalt fibers. Table 4 summarizes the material data for the components.
The validation results are summarized in Table 5. The best predictions (exhibiting the smallest error) are highlighted and underlined. In all the cases they appear to be the TEE method estimations.

4. Conclusions

The use of hybrid composites for semi-structural applications is increasing in many branches of industry. Hybrid composites are often more advantageous than classical two-component composites due to their environmental friendliness, weight reduction, and lower production costs. Structural applications of hybrid composites require knowledge of their mechanical properties, which result from the properties of constituent materials as well as their volume fractions and geometrical configuration. The commonly used estimation methods (for example the rule of mixtures) usually give results with a large error. For this reason, there is a need to develop a homogenization technique that would be accurate and preferably conceptually simple, not requiring complicated calculations.
The paper describes a new method of hybrid composite properties estimation. The approach is based on the mechanical equivalence between a real, heterogeneous material, composed of a matrix and reinforced inclusions, and a fictitious pseudo-homogeneous configuration. The fictitious pseudo-homogeneous material is mechanically equivalent to the real composite in the sense of the total internal strain energy being the same for both materials.
The proposed approach to properties estimation exhibits several important advantages in comparison with the most common homogenization methods. First of all, it has a physical, energy-based background. Additionally, the analysis is performed at the macro-scale, with the use of a mathematical formalism that does not require any assumptions about the exact size, type and distribution of reinforcing inclusions. The method can be used for an arbitrary number of composite components without any further simplifications or the need for multi-step calculations. The predictions are very close to the experimental values, which confirms that the physical assumptions of the approach are correct.
One of the most important features of the proposed estimation method is its simplicity: there is no need to numerically homogenize the composite components. What is more, the estimates are given by the smooth analytical functions of the component material properties, while there is no need to identify any additional parameters. Such functions are convenient for example for optimization problems.
In the present paper, the method was validated for isotropic short-fiber reinforced hybrid composites in the elastic range of response. However, it is possible to extend the approach to anisotropic materials. In that case instead of a scalar volume fraction parameter ξ , a tensorial quantity should be utilized, which can reflect not only the fraction of volume but also the geometrical configuration of inclusions. The proposed approach can also be applied to inelastic composite parameters estimations due to the contribution of both elastic and inelastic strains to the total strain energy.
Finally, it should be pointed out that the proposed approach does not take into account some of the conditions that affect the effective properties of hybrid composites. For example, the influence of the reinforcing particle size is here disregarded, as well as the (surface) varying properties of interphases.

Author Contributions

Conceptualization, methodology, writing—original draft preparation, A.J. and H.E.; formal analysis, writing—review and editing, H.E.; validation, visualization, A.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Science Centre of Poland through the Grant No. 2018/31/N/ST8/01052.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Haris, N.I.N.; Hassan, M.Z.; Ilyas, R.A.; Suhot, M.A.; Sapuan, S.M.; Dolah, R.; Mohammad, R.; Asyraf, M.R.M. Dynamic mechanical properties of natural fiber reinforced hybrid polymer composites: A review. J. Mater. Res. Technol. 2022, 19, 167–182. [Google Scholar] [CrossRef]
  2. Jamir, M.R.M.; Majid, M.S.A.; Khasri, A. Natural Lightweight Hybrid Composites for Aircraft Structural Applications. In Sustainable Composites for Aerospace Applications; Jawaid, M., Thariq, M., Eds.; Woodhead Publishing Series in Composites Science and Engineering; Woodhead Publishing: Cambridge, UK, 2018; pp. 155–170. [Google Scholar]
  3. Shahzad, A. Investigation into Fatigue Strength of Natural/Synthetic Fiber-Based Composite Materials. In Mechanical and Physical Testing of Biocomposites, Fibre-Reinforced Composites and Hybrid Composites; Jawaid, M., Thariq, M., Saba, N., Eds.; Woodhead Publishing Series in Composites Science and Engineering; Woodhead Publishing: Cambridge, UK, 2019; pp. 215–239. [Google Scholar]
  4. Dong, C. Review of natural fibre-reinforced hybrid composites. J. Reinforc. Plast. Compos. 2017, 37, 331–348. [Google Scholar] [CrossRef] [Green Version]
  5. Saleem, A.; Medina, L.; Skrifvars, M.; Berglin, L. Hybrid Polymer Composites of Bio-Based Bast Fibers with Glass, Carbon and Basalt Fibers for Automotive Applications—A Review. Molecules 2020, 25, 4933. [Google Scholar] [CrossRef] [PubMed]
  6. Muthalagu, R.; Murugesan, J.; Sathees Kumar, S.; Sridhar Babu, B. Tensile attributes and material analysis of kevlar and date palm fibers reinforced epoxy composites for automotive bumper applications. Mater. Today Proc. 2021, 46, 433–438. [Google Scholar] [CrossRef]
  7. Ganesarajan, D.; Simon, L.; Tamrakar, S.; Kiziltas, A.; Mielewski, D.; Behabtu, N.; Lenges, C. Hybrid composites with engineered polysaccharides for automotive lightweight. Compos. Part C Open Access 2022, 7, 100222. [Google Scholar] [CrossRef]
  8. Jaafar, C.N.A.; Zainol, I.; Ishak, N.S.; Ilyas, R.A.; Sapuan, S.M. Effects of the liquid natural rubber (LNR) on mechanical properties and microstructure of epoxy/silica/kenaf hybrid composite for potential automotive applications. J. Mater. Res. Technol. 2021, 12, 1026–1038. [Google Scholar] [CrossRef]
  9. Mansor, M.R.; Sapuan, S.M.; Zainudin, E.S.; Nuraini, A.A.; Hambali, A. Hybrid natural and glass fibers reinforced polymer composites material selection using Analytical Hierarchy Process for automotive brake lever design. Mater. Des. 2013, 51, 484–492. [Google Scholar] [CrossRef]
  10. Dattatreya, K.; Sathees Kumar, S.; Prasad, V.V.S.H.; Pati, P.R. Mechanical properties of waste natural fibers/fillers reinforced epoxy hybrid composites for automotive applications. Mater. Today Proc. 2023. [Google Scholar] [CrossRef]
  11. Bagheri, Z.S.; Sawi, I.E.; Schemitsch, E.H.; Zdero, R.; Bougherara, H. Biomechanical properties of an advanced new carbon/flax/epoxy composite material for bone plate applications. J. Mech. Behav. Biomed. Mater. 2013, 20, 398–406. [Google Scholar] [CrossRef]
  12. Scutaru, M.L.; Baba, M. Investigation of the Mechanical Properties of Hybrid Carbon-Hemp Laminated Composites Used as Thermal Insulation for Different Industrial Applications. Adv. Mech. Eng. 2014, 6, 829426. [Google Scholar] [CrossRef] [Green Version]
  13. Cicala, G.; Cristaldi, G.; Recca, G.; Ziegmann, G.; El-Sabbagh, A.; Dickert, M. Properties and performances of various hybrid glass/natural fibre composites for curved pipes. Mater. Des. 2009, 30, 2538–2542. [Google Scholar] [CrossRef]
  14. Critchfield, M.O.; Judy, T.D.; Kurzweil, A.D. Low-cost design and fabrication of composite ship structures. Mar. Struct. 1994, 7, 475–494. [Google Scholar] [CrossRef]
  15. Panthapulakkal, S.; Sain, M. Injection-molded short hemp fiber/glass fiber-reinforced polypropylene hybrid composites—Mechanical, water absorption and thermal properties. J. Appl. Polym. Sci. 2007, 103, 2432–2441. [Google Scholar] [CrossRef]
  16. Voigt, W. Über die Beziehungen zwischen den beiden Elastizitäskonstanten isotroper Körper. Wied. Ann. 1889, 38, 573–587. [Google Scholar] [CrossRef] [Green Version]
  17. Reuss, A. Berechnung der Fliessgrenze von Mischkristallen auf Grund der Plastizitätsbedingung für Einkristalle. Z. Angew. Math. Mech. 1929, 9, 49–58. [Google Scholar] [CrossRef]
  18. Hill, R. A theory of the yielding and plastic flow of anisotropic materials. Proc. R. Soc. Lond. A 1948, 193, 281–297. [Google Scholar]
  19. Hashin, Z.; Shtrikman, S. A variational approach to the theory of the elastic behaviour of polycrystals. J. Mech. Phys. Solid 1962, 10, 343–352. [Google Scholar] [CrossRef]
  20. Miloh, T.; Benveniste, Y. A Generalized Self-Consistent Method for the Effective Conductivity of Composites With Ellipsoidal Inclusions and Cracked Bodies. J. Appl. Phys. 1988, 63, 789–796. [Google Scholar] [CrossRef]
  21. Nazarenko, L.; Stolarski, H.; Altenbach, H. Effective properties of particulate composites with surface-varying interphases. Compos. B Eng. 2018, 149, 268–284. [Google Scholar] [CrossRef]
  22. Sun, C.T.; Vaidya, R.S. Prediction of composite properties from a representative volume element. Compos. Sci. Technol. 1996, 56, 171–179. [Google Scholar] [CrossRef]
  23. Gan, H.; Orozco, C.E.; Herakovich, C.T. A strain-compatible method for micromechanical analysis of multi-phase composites. Int. J. Solid. Struct. 2000, 37, 5097–5122. [Google Scholar] [CrossRef]
  24. Liu, B.; Feng, X.; Zhang, S.-M. The effective Young's modulus of composites beyond the Voigt estimation due to Poisson effect. Compos. Sci. Technol. 2009, 69, 2198–2204. [Google Scholar] [CrossRef]
  25. Würkner, M.; Berger, H.; Gabbert, U. On numerical evaluation of effective material properties for composite structures with rhombic fiber arrangements. Int. J. Eng. Sci. 2011, 49, 322–332. [Google Scholar] [CrossRef]
  26. Selvadurai, A.P.S.; Nikopour, H. Transverse elasticity of a unidirectionally reinforced composite with an irregular fibre arrangement: Experiments, theory and computations. Compos. Struct. 2012, 94, 1973–1981. [Google Scholar] [CrossRef]
  27. Arjun, K.P.; de Barros, S.; Budhe, S. Theoretical determination of elastic and flexural modulus for inter-ply and intra-ply hybrid composite material. Compos. Struct. 2022, 281, 114971. [Google Scholar] [CrossRef]
  28. Mori, T.; Tanaka, K. Average stress in matrix and average elastic energy of materials with misfitting inclusions. Acta Metal. 1973, 21, 571–574. [Google Scholar] [CrossRef]
  29. Kanaun, S.K.; Jeulin, D. Elastic properties of hybrid composites by the effective field approach. J. Mech. Phys. Solids 2001, 49, 2339–2367. [Google Scholar] [CrossRef]
  30. Du, H.; Mazzeo, A.D.; Shan, J.W.; Xia, X.; Weng, G.J. Electrical response, elastic property, and pressure sensing under bending of hybrid graphene/CNT/elastomer nanocomposites. Compos. Struct. 2023, 311, 116838. [Google Scholar] [CrossRef]
  31. Hasanzadeh, M.; Ansari, R.; Hassanzadeh-Aghdam, M.K. Evaluation of effective properties of piezoelectric hybrid composites containing carbon nanotubes. Mech. Mater. 2019, 129, 63–79. [Google Scholar] [CrossRef]
  32. Bazan, P.; Nosal, P.; Wierzbicka-Miernik, A.; Kuciel, S. A novel hybrid composites based on biopolyamide 10.10 with basalt/aramid fibers: Mechanical and thermal investigation. Compos. Part B Eng. 2021, 223, 109125. [Google Scholar] [CrossRef]
  33. Halpin, J.C.; Tsai, S.W. Environmental factors in composite materials design. US Air Force Materials Lab. Rep. AFML Tech. Rep. 1969, 67423, 749–767. Available online: https://www.semanticscholar.org/paper/Effects-of-Environmental-Factors-on-Composite-Halpin/1a8502fad082b55755965085de1b6dad8edcbca2 (accessed on 13 February 2023).
  34. Ribeiro, F.J.V.; Filho, S.L.M.R.; Silveira, E.M.; Panzera, T.H.; Scarpa, F.; Tonatto, M.L.P. A sequential multiscale technique to evaluate the mechanical behaviour of hybrid composites containing carbon fibre and silica microparticles. Compos. Struct. 2023, 314, 116977. [Google Scholar] [CrossRef]
  35. Facca, A.G.; Kortschot, M.T.; Yan, N. Predicting the elastic modulus of natural fibre reinforced thermoplastics. Compos. Part A 2006, 37, 1660–1671. [Google Scholar] [CrossRef]
  36. Alkandary, A.; Netravali, A.N. Hybrid green composites using rice straw and jute fabric as reinforcement for soy protein-based resin. Compos. Part B Eng. 2023, 256, 110626. [Google Scholar] [CrossRef]
  37. Wang, A.; Liu, X.; Yue, Q.; Xian, G. Tensile properties hybrid effect of unidirectional flax/carbon fiber hybrid reinforced polymer composites. J. Mater. Res. Technol. 2023, 24, 1373–1389. [Google Scholar] [CrossRef]
  38. Wiśniewska, A.; Egner, H. Optimization of Functionally Graded Structural Members by Means of New Effective Properties Estimation Method. Materials 2019, 12, 3139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Egner, H.; Ryś, M. Total energy equivalence in constitutive modelling of multidissipative materials. Int. J. Damage Mech. 2017, 3, 417–446. [Google Scholar] [CrossRef]
  40. Ryś, M.; Egner, H. Energy equivalence based constitutive model of austenitic stainless steel at cryogenic temperatures. Int. J. Solids Struct. 2019, 164, 52–65. [Google Scholar] [CrossRef]
  41. Egner, H.; Egner, W. Classification of Constitutive Equations for Dissipative Materials-General Review. In Mechanics of Anisotropic Materials; Skrzypek, J.J., Ganczarski, A.W., Eds.; Springer International: Cham, Switzerland, 2015; pp. 247–294. [Google Scholar]
  42. Wiśniewska, A.; Hernik, S.; Liber-Kneć, A.; Egner, H. Effective properties of composite material based on total strain energy equivalence. Compos. Part B Eng. 2019, 166, 213–220. [Google Scholar] [CrossRef]
  43. Wiśniewska, A.; Hernik, S.; Egner, H. Effective Properties of Composite Material Based on Total Strain Energy Equivalence. In Plasticity, Damage and Fracture in Advanced Materials. Advanced Structured Materials; Altenbach, H., Brünig, M., Kowalewski, Z., Eds.; Springer: Cham, Switzerland, 2020; Volume 121, pp. 201–213. [Google Scholar]
  44. Saanouni, K.; Forster, C.; Ben Hatira, F. On the anelastic flow with damage. Int. J. Damage Mech. 1994, 3, 140–169. [Google Scholar] [CrossRef]
  45. Murakami, S. Continuum Damage Mechanics: A Continuum Mechanics Approach to the Analysis of Damage and Fracture; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  46. Skrzypek, J.J.; Ganczarski, A.W.; Rustichelli, F.; Egner, H. Advanced Materials and Structures for Extreme Operating Conditions; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  47. Egner, H. On the full coupling between thermo-plasticity and thermo-damage in thermodynamic modeling of dissipative materials. Int. J. Solids Struct. 2012, 49, 279–288. [Google Scholar] [CrossRef] [Green Version]
  48. Koszkul, J. Kompozyty poliamidu 6 z włóknem szklanym. Composites 2001, 1, 159–162. [Google Scholar]
  49. Charrière, E.; Terrazzoni, S.; Pittet, C.; Mordasini, P.; Dutoit, M.; Lemaître, J.; Zysset, P. Mechanical characterization of brushite and hydroxyapatite cements. Biomaterials 2001, 22, 2937–2945. [Google Scholar] [CrossRef] [PubMed]
  50. Bonfield, W. Hydroxyapatite-Reinforced Polyethylene as an Analogous Material for Bone Replacement. Ann. N. Y. Acad. Sci. 1988, 523, 173–177. [Google Scholar] [CrossRef] [PubMed]
  51. Kuciel, S.; Kornobis, M. Kompozyty hybrydowe na osnowie polilaktydu z włóknami węglowymi i bazaltowymi. Przetw. Tworz. 2017, 4, 344–350. [Google Scholar]
  52. Lai, J. Moldflow Material Testing Report. MAT2238, NatureWorks PLA; Moldflow: Melbourne, Australia, 2007. [Google Scholar]
  53. Bino, D.; Retnam, B.S.J.; Ramachandran, M.; Sivapragash, M. Analysis of mechanical properties of glass and carbon fiber reinforced polymer material. Int. J. Appl. Eng. 2015, 10, 10387–10391. [Google Scholar]
Figure 1. Schematic representation of energy equivalence idea between real (R) and fictitious (F) configurations.
Figure 1. Schematic representation of energy equivalence idea between real (R) and fictitious (F) configurations.
Materials 16 04215 g001
Figure 2. Comparison between the theoretical estimations of composite elastic properties and experimental data. Data adapted from [48,49,50].
Figure 2. Comparison between the theoretical estimations of composite elastic properties and experimental data. Data adapted from [48,49,50].
Materials 16 04215 g002
Figure 3. Distribution of effective Young modulus depending on the volume fraction of (a) basalt fibers ξ 2 with a constant content of aramid fibers ξ 1 ; (b) aramid fibers ξ 1 with a constant content of basalt fibers ξ 2 .
Figure 3. Distribution of effective Young modulus depending on the volume fraction of (a) basalt fibers ξ 2 with a constant content of aramid fibers ξ 1 ; (b) aramid fibers ξ 1 with a constant content of basalt fibers ξ 2 .
Materials 16 04215 g003
Figure 4. Comparison of effective Young modulus predictions depending on the volume fraction of basalt fibers ξ 2 with constant content of aramid fibers ξ 1 for total strain energy equivalence (TEE), Reuss (R), and Voigt (V) methods.
Figure 4. Comparison of effective Young modulus predictions depending on the volume fraction of basalt fibers ξ 2 with constant content of aramid fibers ξ 1 for total strain energy equivalence (TEE), Reuss (R), and Voigt (V) methods.
Materials 16 04215 g004
Figure 5. Distribution of effective Young modulus depending on the volume fraction of both types of fibers for the TEE, Reuss (R), and Voigt (V) methods.
Figure 5. Distribution of effective Young modulus depending on the volume fraction of both types of fibers for the TEE, Reuss (R), and Voigt (V) methods.
Materials 16 04215 g005
Figure 6. Comparison TEE, Reuss (R), and Voigt (V) methods with experimental results for: (a) polyamide reinforced with aramid fibers ξ 1 . and basalt fibers ξ 2 , variable volume fraction of aramid and constant volume fraction of basalt; (b) polyamide reinforced with aramid fibers ξ 1 and basalt fibers ξ 2 , variable volume fraction of basalt and constant volume fraction of aramid; (c) polylactide reinforced with carbon fibers ξ 1 and basalt fibers ξ 2 , variable volume fraction of carbon and constant volume fraction of basalt; (d) polylactide reinforced with carbon fibers ξ 1 and basalt fibers ξ 2 , variable volume fraction of basalt and constant volume fraction of carbon.
Figure 6. Comparison TEE, Reuss (R), and Voigt (V) methods with experimental results for: (a) polyamide reinforced with aramid fibers ξ 1 . and basalt fibers ξ 2 , variable volume fraction of aramid and constant volume fraction of basalt; (b) polyamide reinforced with aramid fibers ξ 1 and basalt fibers ξ 2 , variable volume fraction of basalt and constant volume fraction of aramid; (c) polylactide reinforced with carbon fibers ξ 1 and basalt fibers ξ 2 , variable volume fraction of carbon and constant volume fraction of basalt; (d) polylactide reinforced with carbon fibers ξ 1 and basalt fibers ξ 2 , variable volume fraction of basalt and constant volume fraction of carbon.
Materials 16 04215 g006
Table 1. Comparison of estimation errors for data from Figure 2.
Table 1. Comparison of estimation errors for data from Figure 2.
Composite ξ (-) Experimental   Young   Modulus   E (GPa)Estimation Method Error (%)
Energy Equivalence Based Method (TEE)Voigt (V)Reuss (R)Hashin- Shtrikman Upper (HSU)Hashin- Shtrikman Lower (HSL)Mori-Tanaka (MT)
POM-CF0.003.420.000.000.000.000.000.00
0.057.0310.31107.8248.8332.6446.1746.17
0.109.644.39167.5360.6560.4556.5756.57
0.2016.5821.83190.3974.3173.2069.0369.03
0.3024.0040.94193.7779.7781.1773.5873.58
PA-GF0.054.0413.1469.7816.8228.2612.5912.59
0.075.082.0363.0732.4917.8427.7027.70
0.137.776.8960.8152.978.9246.8346.83
0.168.572.2270.1255.9213.7048.8748.87
0.2311.632.3166.9264.7310.6456.7256.72
0.3113.0617.0990.7365.1827.9154.7054.70
PE-HAp0.001.300.000.000.000.000.000.00
0.101.4049.99110.342.2559.0911.7311.73
0.202.0053.18124.3820.3760.585.985.98
0.252.5044.16109.5132.5049.2617.4717.47
0.303.0039.4199.2840.1842.4424.4624.46
0.353.7029.6681.3948.2230.7932.6332.63
0.404.4023.8469.0753.2923.4537.5537.55
0.455.904.0038.3362.442.6148.5448.54
0.507.7010.9315.3368.7812.8856.2856.28
Table 2. Material properties of composite components.
Table 2. Material properties of composite components.
Material Bulk   Modulus   K
(GPa)
Lam é   Constant   λ
(GPa)
Lam é   Constant   μ
(GPa)
PA2.361.950.62
AF49.1234.8721.37
BF61.8138.2635.32
Table 3. Material properties of matrix polyamide, aramid fibers, and basalt fibers (data taken from [32]).
Table 3. Material properties of matrix polyamide, aramid fibers, and basalt fibers (data taken from [32]).
Material Young   Modulus   E
(GPa)
Poisson   Ratio   ν
(-)
PA1.700.38
AF56.000.31
BF89.000.26
Table 4. Material properties of matrix polylactide, carbon fibers, and basalt fibers (data taken from [32,51,52,53]).
Table 4. Material properties of matrix polylactide, carbon fibers, and basalt fibers (data taken from [32,51,52,53]).
Material Young   Modulus   E
(GPa)
Poisson   Ratio   ν
(-)
PLA3.760.36
CF242.000.32
BF89.000.26
Table 5. Comparison of numerical values and estimation errors for data from Figure 5 and Figure 6.
Table 5. Comparison of numerical values and estimation errors for data from Figure 5 and Figure 6.
Composite ξ 1 ( ) ξ 2   ( ) Experimental   E (GPa) Estimation Method
TEEVoigt (V)Reuss (R)
Value   E (GPa) Error (%) Value   E (GPa) Error (%) Value   E (GPa) Error (%)
PA-AF-BF0.0360.0202.892.870.615.4187.191.8037.77
0.0550.0303.313.598.617.33121.401.8543.99
0.0740.0404.324.422.229.29115.071.9155.69
PLA-CF-BF0.0520.0358.848.493.9319.12116.274.1153.52
0.0880.06112.1312.936.6229.85146.124.4063.77
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jerzyńska, A.; Egner, H. Energy Equivalence Based Estimation of Hybrid Composites Mechanical Properties. Materials 2023, 16, 4215. https://doi.org/10.3390/ma16124215

AMA Style

Jerzyńska A, Egner H. Energy Equivalence Based Estimation of Hybrid Composites Mechanical Properties. Materials. 2023; 16(12):4215. https://doi.org/10.3390/ma16124215

Chicago/Turabian Style

Jerzyńska, Anna, and Halina Egner. 2023. "Energy Equivalence Based Estimation of Hybrid Composites Mechanical Properties" Materials 16, no. 12: 4215. https://doi.org/10.3390/ma16124215

APA Style

Jerzyńska, A., & Egner, H. (2023). Energy Equivalence Based Estimation of Hybrid Composites Mechanical Properties. Materials, 16(12), 4215. https://doi.org/10.3390/ma16124215

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop