Journal of Applied Mathematics
Volume 2018 (2018), Article ID 3196569, 9 pages
https://doi.org/10.1155/2018/3196569

Research Article

A Theoretical Consideration on the Estimation of Interphase Poisson’s Ratio for Fibrous Polymeric Composites

J. Venetis¹ and E. Sideridis²

¹School of Applied Mathematics and Physical Sciences NTU Athens, Section of Mechanics, 5 Heroes of Polytechnion Avenue, 15773 Athens, Greece
²School of Applied Mathematics and Physical Sciences NTUA, Section of Mechanics, 5 Heroes of Polytechnion Avenue, 15773 Athens, Greece

Correspondence should be addressed to J. Venetis; johnvenetis4@gmail.com

Received 3 July 2018; Revised 23 August 2018; Accepted 5 September 2018; Published 1 October 2018

Academic Editor: Xin-Lin Gao

Copyright © 2018 J. Venetis and E. Sideridis. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An analytical approach on the evaluation of interphase Poisson’s ratio for fibrous composites, consisting of polymeric matrix and unidirectional continuous fibers, is performed. The simulation of the microstructure of the composite was carried out by means of a modified form of Hashin-Rosen cylinder assemblage model. Next, by the use of this three-phase model the authors impose some limitations to the polynomial variation laws which are commonly adopted to approximate the thermomechanical properties of the interphase layer of this type of polymeric composites and then propose an nth-degree polynomial function to approximate the Poisson’s ratio of this layer.

1. Introduction

The rigorous description of a composite system consisting of a matrix in which continuous or short fibers have been dispersed is not an easy task. Indeed, a great number of parameters, geometrical, topological, mechanical, etc. are necessary, the majority of which vary in a stochastic manner or are practically unknown. Theoretical treatments usually attempt to exploit as far as possible readily available information, which, in most cases, consists of the mechanical properties of matrix and fiber and the volume fraction of the latter, whilst appropriate assumptions cover missing data. On the other hand it is well known that in lightweight structures, large specific stiffness, and strength are aiming properties for a material. By combining fibers with an appropriate polymer and controlling the production procedure it is possible to manufacture composites featuring large specific properties. Moreover, the constituents are usually cheap and easily processed, e.g., by injection molding. Hence, fibrous composites have many industrial applications where high performance per weight at a reasonable price is required. One of the most useful forms of composites for the construction of high-performance structural elements is the type of panels made from aligned fibers containing polymerized matrix. Evidently, their mechanical properties depend on the related properties and volume fractions of the constituent materials, the fiber length or aspect ratio, the degree of the alignment, the adhesion between fibers, and matrix and last but not least they are affected by the fabrication techniques [1]. As it was initially stated, the investigation of the elastic properties of unidirectional fibrous composites reinforced with long or short fibers in terms of the related properties of constituents constitutes a very difficult problem of applied mechanics. In a fundamental investigation by Hashin and Rosen [2] bounds and explicit expressions to estimate the effective elastic moduli of unidirectional fibrous composites were derived via a rigorous variational method. On the other hand, amongst the analytical models presented in the literature, some of them take into consideration the existence of an intermediate phase developed during the preparation of the polymeric composite. Evidently this new phase plays an important role in the overall thermomechanical behavior of the composite.

In the meanwhile, the existence of a boundary interphase in polymeric composites (fibrous and pariculate) was shown experimentally by Lipatov [3] who measured its thickness by means of Differential Scanning Calorimetry (DSC) experiments. Also, in this valuable investigation, it was empirically stated that the size of these heat capacity jumps for unfilled and filled polymers is expressed in terms of the thickness of this phase. Yet, the estimation of its physical properties was a remaining problem.

To this end, in a preliminary model performed by Papanicolaou et al. [4] and Theocaris et al. [5] this intermediate phase was initially assumed to be a homogeneous and isotropic material whereas in a better approximation [6] a more complex model was introduced, according to which the fiber was surrounded by a series of successive cylinders, each one of them having a different elastic modulus in a step-function variation with respect to the polar radius.

In addition, problems of obtaining closed-form solutions to many cases of elastic inclusions embedded in an elastic matrix were encountered by Muskhelishvili [7] by the use of conformal mapping techniques. In addition, numerical solution techniques such as finite difference and finite elements have been extensively used to this end [8, 9]. Another remarkable consideration of the variable modulus interphase is the so-called unfolding model, which was carried out on the basis that the boundary interphase layer constitutes a transition zone between fibers with high moduli and matrix with rather low stiffness [10]. An extension of the above fundamental model to more complex composite systems (fibrous and particulate) was exhibited by Theocaris in [11].

In addition, a remarkable study performed by Swain et al [12] has thrown light on the effect of the boundary interphase on the properties and performance of composite materials and laminates.

In the past years, there is a lot of recent research work carried out for the evaluation of elastic and thermal properties of fibrous composites and for the investigation of the influence of various parameters on these properties. A detailed study on the effect of interphase layers on local thermal displacements in fibrous composites was presented in [13]. Also, in [14], an investigation was made by the aid of strength of materials and elasticity approach, to derive explicit expressions for longitudinal, transverse, and shear moduli along with longitudinal and transverse Poisson’s ratio using the concept of boundary interphase between fiber and matrix, in the framework of a three-phase modified form of Hashin-Rosen cylinder assemblage model that was presented in [2]. Here, to evaluate the elastic properties of interphase a unified mode of variation was considered. In particular, an nth degree polynomial variation was initially adopted in terms of the radius of the introduced three-phase coaxial model and finally for facility reasons a quadratic function was taken into account. The above approach of interphase unknown properties was continued and extended to estimate the thermal conductivity of the same class of fibrous composites in [15]. In this work, to cover the whole spectrum of interphase thermal conductivity, five different laws of variation were used of which only two constitute polynomial forms (linear and parabolic). Nevertheless, it can be said that a three-layer cylindrical model to simulate the structure of unidirectional fibrous composites reinforced with continuous fibers may arise based on the theory of self-consistent models and adapting this approach to a three-layered cylinder, delineating the representative volume element for the fibrous composite [11]. Here one should elucidate that a significant variation of the self-consistent model is the three-phase model first introduced by Kerner [16], according to which the inclusion is developed by a matrix annulus which in turn is embedded in an infinite medium with the unknown macroscopic properties of the composite. Further, in [17, 18] the adoption of linear variation laws for interphase elastic constants yielded theoretical predictions of the composite moduli in a reasonable accordance with experimental values. Yet, to approach the interphase modulus and Poisson’s ratio by a quadratic (parabolic) law with respect to the radius of the adopted model of embedded cylinders is many times more convenient since many times the linear variation law cannot alleviate the fact that the transition of the elastic constants from the matrix to fiber is carried out by “jumps” in their characteristic properties [17, 19]. However, a significant problem that remains is the following: letting an interphase property be approached by a quadratic polynomial in the general form , there are three coefficients needed to be found. Unfortunately, in the context of the modified form of Hashin-Rosen model introduced in [14] there are only two boundary conditions available. Specifically, at and at . Another condition may arise by requiring the first derivative to vanish at a critical point over the interval which corresponds to a given local extremum (maximum or minimum) [14]. In particular it was assumed that the critical point coincides with one of the endpoints of this interval and thus . Nevertheless, should this critical point be outside this interval or in the interval such a condition cannot hold [20] and therefore the coefficients are unable to be calculated. To overcome this unfortunate situation, in the present work an alternative approach on the estimation of interphase elastic constants is proposed.

2. Analysis

Let us simulate the microstructure of a unidirectional fibrous composite by means of a coaxial three-phase cylinder unit cell, the cross-sectional area of which can be seen in Figure 1.

Figure 1: Coaxial three-phase cylinder unit cell.

Here the radii denote the bounds of each separate zone, i.e., fiber, interphase, and matrix, respectively. The fiber zone begins at the zero value of the radius of the three-phase cylindrical model and ends at . Next, the interphase zone, starts marginally at and ends at . Finally, the matrix zone starts marginally at and finishes at .

Evidently, the domain of definition of any single-valued continuous function that can be selected to approach the Poisson’s ratio of this intermediate phase between fiber and matrix is the interval .

The above modified form of Hashin-Rosen cylinder assemblage model was initially used in [14]. In the sense of the above-mentioned work the elastic constants of the interphase, i.e., stiffness, shear modulus, and Poisson’s ratio can be generally approached as nth degree polynomials with respect to the polar radius of the above coaxial three-phasemodel of embedded cylinders.

Suggestively, let us focus on Poisson’s ratio of the interphase zone and express it with respect to radius of the above model, as follows:

Hence, (1) denotes the mode of variation for Poisson’s ratio of this intermediate phase the thickness of which is evidently and can be experimentally measured as we shall mention later on.

For facility reasons and to be in accordance with inequality one may consider a polymeric fibrous composite the constituents of which are epoxy resin and e-glass fibers, i.e., .

Moreover, at the interfaces amongst the three different phases of the overall material the following boundary conditions hold [14]:

Here on the basis of inequalities and one may pinpoint that should a linear variation law be adopted to approach the Poisson’s ratio of this layer; the single valued function of the form is strictly increasing. Yet, as stated above such an approach constitutes an oversimplified mode of variation.

In continuing, let us consider the following general quadratic function to approximate the interphase Poisson’s ratio:

with being arbitrarily selected real constants

Now, the mode of variation for the Poisson’s ratio of the boundary interphase is reduced to a second-degree polynomial function in terms of polar radius of the cylindrical three-phase model.

Nonetheless, the remaining problem is the determination of the three coefficients given that only two boundary conditions are generally available, i.e., (2) and (3).

Now, without violating the generality, let us set

where are arbitrary real numbers

Thus (4) becomes

Here one may observe that the latter relationship is completely synonymous to (4) because the totally arbitrary selection of the constant quantities a, b, c from the set of real numbers prohibits any correlation amongst the initial coefficients A, B, C appearing in (4).

Next, by differentiating the above relation with respect to polar radius , one finds

and therefore

Equation (8) can be combined with (6) to yield

Evidently, (9) constitutes a first-order separable ordinary differential equation of the form or and can be solved by trivial techniques [21].

Also, one observes that should the term in the left-hand side of (9) be nonzero the terms agree in sign.

Now, to find the critical value of the variable where the single-valued continuous function has local extremum (minimum or maximum) the quantity should be set equal to zero. Thus it implies that

Besides, at this critical point, (9) yields

Next, the application of the boundary conditions, expressed by (2) and (3), to (9) in association with (11) yields

and

Here, one may conclude that should the terms in the left-hand sides of (12) and (13) be strictly positive the terms and agree in sign.

In this framework, one infers that should the derivatives and be non-zero, the terms agree in sign.

In other words, since , if and , the local extremum does not belong to the interval .

However, given a polymeric composite consisting of matrix and filler (particles or fibers), the values of interphase properties can be neither greater nor less than those of the constituents’ properties [4, 11].

Thus if the critical point yields a value of interphase Poisson’s ratio outside the interval the three coefficients appearing in (4) or in (6) should be calculated only on the basis of the boundary conditions designated by (2), (3) and the validity of any condition concerning the vanishing of over the interval cannot be required.

In addition, by focusing on (4) and taking into account the basic concepts of Polynomial theory [22], one may obtain the critical point of by the following relation which holds identically:

Evidently the fraction on the right hand side of the above identity signifies the unique local expremum of the quadratic variation law of Poisson’s ratio. Also one may point out that (14) is synonymous to the vanishing of the first derivative of since the determination of the critical point can be carried out by setting .

In this context, in order to specify a third condition to complete (2),(3) towards the calculation of the three coefficients or one should know beforehand either the value of in terms of the endpoints or the value of the local extremum in terms of the values corresponding at the endpoints, i.e., with . In [14] it was stated that should the value denote a local minimum it coincides with whereas should denote a local maximum it coincides with . Yet it is the authors’ current opinion on this issue that this assumption is oversimplified and rather unrealistic, although it has yielded theoretical predictions for the composite elastic properties that were found to be in good agreement with several reliable theoretical formulae and experimental data.

Hence should one select a second degree parabola to approach the interphase Poisson’s ratio expressed either by (4) or (6), except the two boundary conditions concerning the interface with fiber and matrix, respectively, another condition is necessary referring either to the critical point of the quadratic function or the value of interphase Poisson’s ratio yielded by this critical point. Evidently, the fact that (14) holds identically prohibits the parallel use of the two above additional conditions designated by (10) and (11), which in general can be presented by the following two expressions:

Here, the notation denotes a binary operation over the set of real numbers combining the elements to produce univocally another real number.

Similarly, the notation designates a binary operation over the set of real numbers combining the elements to produce univocally another real number.

Moreover since andthese binary operations should not be necessarily commutative.

In this context, should the critical value yielded by the operation lie on the interval it is necessary to require the corresponding value of Poisson’s ratio to be strictly greater than and strictly less than .

On the other hand, should the value of Poisson’s ratio be strictly less than or strictly greater than it is necessary to require the value to be outside the interval .

Obviously, similar analyses on the basis of the same reasoning could take place for the other elastic or thermal properties of interphase region, i.e., stiffness, shear modulus, thermal expansion coefficient etc.

Next, in the sense of [14] let us introduce a general third-degree polynomial to approximate the interphase Poisson’s ratio given as

Evidently, (17) signifies a more complicated mode of variation for interphase Poisson’s ratio when compared with (4).

Yet, it is known from Calculus [20] that if this function will have no local extrema over its domain of definition and therefore, except (2) and (3), no additional conditions can be set regarding the first derivative of .

Hence, the consideration of a general third-degree polynomial to approach the interphase Poisson’s ratio or any other elastic property of this phase can be made on the premise that the quantity is strictly positive.

Thus we can write out

Provided that the critical points are given as

Thus in a quite analogous reasoning with that resulting in (15), (16) one may require the validity of the following additional conditions:

or

Apparently the above disjunction is exclusive. Now, should even one of the critical values yielded by the operations and , respectively, lie on the interval it is necessary to require the corresponding value of Poisson’s ratio to be strictly greater than and concurrently strictly less that .

On the other hand, should even one the values of Poisson’s ratio be strictly less than or strictly greater than it is necessary to require the corresponding critical point to lie outside the interval .

Moreover, on the basis of (17), it can be proved that the following relationship holds:

For a thorough presentation of the mathematical derivations resulting in (21c) let us refer to Appendix Section.

Nevertheless, since (15), (16), (20a), (20b), (21a), (21b), (21c) have a rather theoretical character let us give a more specific example of the necessity of these conditions by proposing the following single-valued polynomial representation to approximate interphase Poisson’s ratio

with and being arbitrary real numbers

An application of the boundary conditions expressed by (2) and (3) at the interfaces with fiber and matrix, respectively, yields

and therefore

In the sequel, let us calculate the first derivative of with respect to polar radius

and therefore

Here denotes the arithmetic mean of .

Equation (28) can be combined with (25) to yield

At one obtains

By focusing on the above representation one may point out that the critical points at which vanishes cannot coincide with the endpoints neither with their arithmetic mean.

On the other hand, in an analogous manner with that resulting in (9) let us try to relate in an explicit manner the interphase Poisson’s ratio to its first derivative.

In this context, (27) yields

Now, given that the coefficient has been already calculated, the following relation links the critical points , where vanishes, with the coefficient

Equation (32) can be combined with (22) to yield

Solving for one obtains

Equation (34) can be combined with (25) and (26) to yield explicitly the set of the local extrema of interphase Poisson’s ratio in terms of the polynomial degree n. Thus we can write out

Here one observes that the set of local extrema of interphase Poisson’s ratio has been evaluated without the knowledge of the coefficient which ought to be calculated by (30) or by (32).

Evidently, these two equations are equivalent and cannot be considered as a system.

Right here is the necessity of the theoretical representations designated by (15), (16), (20a), (20b) and (21a), (21b), (21c).

In [14] it was assumed in a rather subjective and very simplified mode that whilst .

Apparently, such a consideration is rather unrealistic. However, a reasonable conjecture for the calculation of coefficient is to focus on (30) or (35) which yields the Poisson’s ratio in a closed form and try to search for the critical points by means of other types of means as discussed in [23] since it has been proved that mean-theory is very useful both from theoretical point of view and for practical (engineering) purposes. Yet, such an approach cannot include the case when or and then one may examine the interval . Alternatively, given that most materials (except those ones with auxetic core) have Poisson’s ratio values ranging between 0 and 0.5, one may carry out a parametrical study considering an upper bound of interphase Poisson’s and then substitute it on (35) in order to solve it for .

3. Discussion

The boundary interphase zone was assumed as a natural phase which is developed in reality between fiber and polymer matrix. In this context, it can be said that this intermediate phase is neither an artificial one, e.g., by the immersion of the fibers in an agent, nor a pseudophase being contrived to simulate the microstructure of the composite. Hence it is not possible to know beforehand or to determine the interphase properties, a fact that renders necessary to make assumptions about them. In order to approximate the mode of variation of the variable elastic properties of interphase layer, such as Poisson’s ratio, an nth degree polynomial function was considered by one of the authors in [14]. This function for n = 2 yielded a parabolic law. Also, to assume that at , i.e. at the interface between matrix and interphase, is realistic indeed. Moreover, the assumption that at , i.e., at the interface between fiber and interphase seems reasonable too. However, the remaining problem is to find the coefficients of such a polynomial function given that only the above two boundary conditions are available. The vanishing of the first derivative of interphase Poisson’s ratio (or stiffness, shear modulus, etc.) at the endpoints of the interval that was adopted in [14] seems to be a rather superficial and oversimplified consideration though it yielded realistic theoretical predictions for the elastic constants of the composite.

In the current investigation the authors made an endeavor to shed some light on this difficult task by setting some limitations to the polynomial variation laws of second and third degree in the form of some theoretical formulae expressed by (15), (16), (20a), (20b), (21a), (21b), (21c) and in the sequel they propose an nth degree polynomial function to approximate the variation of interphase Poisson’s ratio.

In this framework, this fundamental elastic property can be estimated provided of course an experimental estimation of the interphase thickness and given that is known beforehand.

To illustrate the physical meaning of , one may emphasize that in reality any polymeric composite (fibrous or particulate) necessarily consists of three distinct phases (matrix, filler, and interphase) as it is stated in [3, 11]

Lipatov [3] has shown that if calorimetric measurements are performed in the neighborhood of the glass transition zone of the composite, energy jumps are observed. These jumps are too sensitive to the amount of filler added to the matrix and can be used to evaluate the boundary layers developed around the inclusions. Apparently, as the filler volume fraction is increased, the proportion of macromolecules characterized by a reduced mobility is also increased. This is equivalent to an augmentation of the interphase concentration by volume and evidently it is in consensus with the conclusion of [3] that the extent of interphase expressed by its thickness is the reason of the variation of the amplitudes of heat capacity jumps appearing at the glass transition zones of the matrix material and the composite with various fiber contents. The size of heat capacity jumps for unfilled and filled materials is directly related to by the next empirical relationship:

where the coefficient is evaluated as

Here the numerator and the denominator of the fraction in the right-hand side of (37) are the sudden changes of the heat capacity for the reinforced polymer and pure one, respectively.

In this context, by using Lipatov’s theory [3] interrelating the abrupt jumps in the specific heat of a composite at its respective glass transition temperatures with the values of the extent of this boundary layer, the interphase thickness can be calculated in an accurate manner and thus the proposed analytical method for the theoretical determination of Poisson’s ratio of this unknown intermediate phase is further enhanced.

4. Conclusions

This paper, as a continuation into a previous research work carried out by one of the authors, had two main objectives:

The first aim was to set some restrictions to the use of a single-valued polynomial function of nth degree to approach the interphase elastic properties for a general class of fibrous composites.

These restrictions concern the set of critical points and the corresponding local extrema (minima and /or peak points).

A modified form of the well-known Hashin-Rosen cylinder assemblage model was adopted to simulate the microstructure of the overall material.

Suggestively, the authors focused on interphase Poisson’s ratio and examined in a rigorous theoretical manner its variations over the interval defined by the interphase thickness , as obtained from a second and third degree parabola.

Next, the authors introduced an nth degree polynomial function with respect to the radius of the coaxial three-phase cylinder model in order to approach the interphase Poisson’s ratio. This function renders very clear the necessity of the indicated restrictions for the set of critical points.

In this context, this fundamental elastic property of the interphase layer can be accurately estimated on the premise that an experimental measurement of the interphase thickness is carried out.

Appendix

In this section, let as perform in detail the mathematical derivations resulting in (21c) which links the maximum and minimum value of interphase Poisson’s ratio when the latter is approached by a third-degree polynomial with respect to the radius of the coaxial three-phase cylindrical model that was adopted to simulate the composite structure.

Evidently, (17) can be formulated as

or equivalently

Here, for facility reasons we have put

Now, by the use of the trivial substitution it implies that

where is a single-valued continuous function.

In addition, note that .

Thus we can write out

and therefore

Equation (A.6) can be combined with (A.3) to yield

Now, for facility reasons, let us set

Thus one may deduce that

The first derivative with respect to auxiliary variable y is given as

Evidently, the necessary and sufficient condition in order for the quantity to have real roots is the coefficients and to disagree in sign. In this framework, one also infers that these real roots should disagree in sign too. Also, since , it implies that the polynomial has a unique local minimum and a unique local maximum.

Let be the real roots of (A.10):

Then without violating the generality we can write out

and

and therefore

In continuing, the above relation can be combined with (A.8) to yield

Finally, since one may deduce that

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

T. W. Clyne and P. J. Withers, An Introduction to composites, Cambridge University Press, Cambridge, 1993.

Z. Hashin and B. W. Rosen, “The elastic moduli of fiber-reinforced materials,” Journal of Applied Mechanics, vol. 31, no. 2, pp. 223–232, 1964.

Y. Lipatov, Physical Chemistry of Filled Polymers, Khimiya Moscow, 1977.

G. C. Papanicolaou, S. A. Paipetis, and P. S. Theocaris, “The concept of boundary interphase in composite mechanics,” Colloid and Polymer Science, vol. 256, no. 7, pp. 625–630, 1978.

P. S. Theocaris and G. C. Papanicolaou, “The effect of the boundary interphase on the thermomechanical behaviour of composites reinforced with short fibers,” Fibre Science and Technology, vol. 12, no. 6, pp. 421–433, 1979.

G. C. Papanicolaou, P. S. Theocaris, and G. D. Spathis, “Adhesion efficiency between phases in fiber-reinforced polymers by means of the concept of boundary interphase,” Colloid and Polymer Science, vol. 258, no. 11, pp. 1231–1237, 1980.

N. I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticity, Noordhoff, Groningen, The Netherlands, 1953.

K. H. Murray, Finite Element Analysis of a Composite Material Interface, Old Dominion University Research Foundation, 1973.

S. N. Atlur, R. H. Gallagher, and O. C. Zienkiewicz, Hybrid and mixed finite element methods, Wiley Publishing House, 1983.

P. S. Theocaris, “The unfolding model for the representation of the mesophase layer in composites,” Journal of Applied Polymer Science, vol. 30, no. 2, pp. 621–645, 1985.

P. S. Theocaris, The Mesophase Concept in Composites, Polymers - Properties and Applications, G. Henrici-Olivé and S. Olivé, Eds., vol. 11, Springer Berlin Heidelberg, Berlin, Heidelberg, 1987.

R. E. Swain, K. L. Reifsnider, K. Jayaraman, and M. el-Zein, “Interface/Interphase Concepts in Composite Material Systems,” Journal of Thermoplastic Composite Materials, vol. 3, no. 1, pp. 13–23, 1990.

N. R. Sottos, R. L. McCullough, and W. R. Scott, “The influence of interphase regions on local thermal displacements in composites,” Composites Science and Technology, vol. 44, no. 4, pp. 319–332, 1992.

E. Sideridis, G. A. Papadopoulos, and E. Kyriazi, “Strength of materials and elasticity approach to stiffness of fibrous composites using the concept of interphase,” Journal of Applied Polymer Science, vol. 95, no. 6, pp. 1578–1588, 2005.

J. Venetis and E. Sideridis, “Thermal conductivity coefficients of unidirectional fiber composites defined by the concept of interphase,” The Journal of Adhesion, vol. 91, no. 4, pp. 262–291, 2015.

E. H. Kerner, “The elastic and Thermo-elastic properties of composite media,” Proceedings of the Physical Society, Section B, vol. 69, no. 8, pp. 808–813, 1956.

J. Venetis and E. Sideridis, “Elastic constants of fibrous polymer composite materials reinforced with transversely isotropic fibers,” AIP Advances, vol. 5, no. 3, 2015.

J. Venetis and E. Sideridis, “Upper bounds of transverse modulus for unidirectional fibrous composites by strength of materials approach,” WSEAS Transactions on Applied and Theoretical Mechanics, vol. 13, pp. 103–116, 2018.

E. Sideridis, E. E. Theotokoglou, and I. Giannopoulos, “Analytical and computational study of the moduli of fiber-reinforced composites and comparison with experiments,” Composite Interfaces, vol. 22, no. 7, pp. 563–578, 2015.

J. Stewart, Calculus, Cengage Learning, USA, Boston, 2015.

P. Hartman, Ordinary differential equations, vol. 38 of Classics in Applied Mathematics, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2002.

P. Borwein and T. Erdélyi, Polynomials and Polynomial Inequalities, vol. 161 of Graduate Texts in Mathematics, Springer, New York, NY, USA, 1995.

M. Raïssouli and M. Chergui, “On Some Inequalities Involving Three or More Means,” Abstract and Applied Analysis, vol. 2016, Article ID 1249604, 8 pages, 2016.