A Theoretical Consideration on the Estimation of Interphase Poisson's Ratio for Fibrous Polymeric Composites.
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 . 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  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  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.  and Theocaris et al.  this intermediate phase was initially assumed to be a homogeneous and isotropic material whereas in a better approximation  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 stepfunction 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  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 . An extension of the above fundamental model to more complex composite systems (fibrous and particulate) was exhibited by Theocaris in .
In addition, a remarkable study performed by Swain et al  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 . Also, in , 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 . 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 . 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 selfconsistent models and adapting this approach to a threelayered cylinder, delineating the representative volume element for the fibrous composite . Here one should elucidate that a significant variation of the self-consistent model is the three-phase model first introduced by Kerner , 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 Mi(r) be approached by a quadratic polynomial in the general form Ar2 + Br + C, there are three coefficients A;B;C needed to be found. Unfortunately, in the context of the modified form of Hashin-Rosen model introduced in  there are only two boundary conditions available. Specifically, at r = r(: Mj(r) = Mf and at r = r; : M;(r) = Mm.
Another condition may arise by requiring the first derivative dMt(r)/dr to vanish at a critical point over the interval [rf, rj which corresponds to a given local extremum (maximum or minimum) . In particular it was assumed that the critical point coincides with one of the endpoints of this interval and thus dMj(r()/dr = 0 V dMt(rt)/dr = 0. Nevertheless, should this critical point be outside this interval or in the interval (rf, r;) such a condition cannot hold  and therefore the coefficients A; B; C 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.
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.
Here the radii [r.sub.f]; [r.sub.f]; [r.sub.m] 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 r = [r.sub.j]. Next, the interphase zone, starts marginally at r [greater than or equal to] [r.sub.j] and ends at r = [r.sub.i]. Finally, the matrix zone starts marginally at r [greater than or equal to] [r.sub.i] and finishes at r = [r.sub.m].
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 [[r.sub.f], [r.sub.i].
The above modified form of Hashin-Rosen cylinder assemblage model was initially used in . 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 r of the above coaxial threephasemodel of embedded cylinders.
Suggestively, let us focus on Poisson's ratio of the interphase zone and express it with respect to radius r of the above model, as follows:
[v.sub.i] (r) = [Ar.sup.n] + [Br.sup.n-1] + [Cr.sup.n-]2 + ***. (1)
with [r.sub.f] [less than or equal to] r [less than or equal to] [r.sub.i] (1a)
and [v.sub.f] [less than or equal to] [v.sub.i] (r) [less than or equal to] [v.sub.m] (1b)
Hence, (1) denotes the mode of variation for Poisson's ratio of this intermediate phase the thickness of which is evidently [r.sub.i]-[r.sub.j] 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., [v.sub.f] = 0.2; [v.sub.m] = 0.35.
Moreover, at the interfaces amongst the three different phases of the overall material the following boundary conditions hold :
At r = [r.sub.f]: [v.sub.i] (r) = [v.sub.f] (2)
At r = [r.sub.i]: [v.sub.i] (r) = [v.sub.m] (3)
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 [v.sub.i](r) = Ar + B 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:
[v.sub.i] (r) =[ Ar.sup.2] + Br + C (4)
with A; B; C 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 r of the cylindrical three-phase model.
Nonetheless, the remaining problem is the determination of the three coefficients A; B; C given that only two boundary conditions are generally available, i.e., (2) and (3).
Now, without violating the generality, let us set
[mathematical expression not reproducible] (5)
where a, b, c are arbitrary real numbers
Thus (4) becomes
[v.sub.i] (r) = a/4 [r.sup.2] + a x br + (a x [b.sup.2] + c) (6)
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 r, one finds
d/dr [v.sub.i] (r) = a (r/2 + b) (7)
[mathematical expression not reproducible] (8)
Equation (8) can be combined with (6) to yield
[mathematical expression not reproducible] (9)
Evidently, (9) constitutes a first-order separable ordinary differential equation of the form N(y)y' = Q(x) or N(y)dy = Q(x)dx and can be solved by trivial techniques .
Also, one observes that should the term in the left-hand side of (9) be nonzero the terms a; ([v.sub.i](r) - c) agree in sign.
Now, to find the critical value of the variable r where the single-valued continuous function v;(r) has local extremum (minimum or maximum) the quantity (d/dr)[v.sub.i](r) should be set equal to zero. Thus it implies that
[mathematical expression not reproducible] (10)
Besides, at this critical point, (9) yields
[mathematical expression not reproducible] (11)
Next, the application of the boundary conditions, expressed by (2) and (3), to (9) in association with (11) yields
[mathematical expression not reproducible] (12)
[mathematical expression not reproducible] (13)
Here, one may conclude that should the terms in the left-hand sides of (12) and (13) be strictly positive the terms a; ([v.sub.f] - [v.sub.i]; ([r.sub.cr])) and a; ([v.sub.m] - [v.sub.i]([r.sub.cr])) agree in sign.
In this framework, one infers that should the derivatives (d/dr)[v.sub.i]([r.sub.f]) and (d/dr)[v.sub.i]([r.sub.m]) be non-zero, the terms ([v.sub.f] - [v.sub.i]]; ([r.sub.cr]); ([v.sub.m] - [v.sub.i]([r.sub.cr]) agree in sign.
In other words, since [r.sub.cr] = -2b, if [r.sub.f] [not equal to] -2b and rm [not equal to] -2b, the local extremum [v.sub.i] ([r.sub.c]r) does not belong to the interval [[v.sub.f], [v.sub.m]].
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 [r.sub.cr] yields a value of interphase Poisson's ratio outside the interval [[v.sub.f], [v.sub.m]] the three coefficients A;B;C appearing in (4) or a;b;c 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 ((d/dr)[v.sub.i](r)) over the interval [[r.sub.f], [r.sub.i] cannot be required.
In addition, by focusing on (4) and taking into account the basic concepts of Polynomial theory , one may obtain the critical point of [v.sub.i] (r) by the following relation which holds identically:
[v.sub.i] (B-2A) [equivalent to] 4AC - [B.sup.2]/4A (14)
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 [v.sub.i] since the determination of the critical point can be carried out by setting (d/dr)[v.sub.i](r) = 0.
In this context, in order to specify a third condition to complete (2),(3) towards the calculation of the three coefficients A, B, C or a, b, c one should know beforehand either the value of [r.sub.cr] in terms of the endpoints [r.sub.f], [r.sub.i] or the value of the local extremum in terms of the values corresponding at the endpoints, i.e., [v.sub.f], [v.sub.m] with [v.sub.f] < [v.sub.m]. In  it was stated that should the value (4AC - [B.sup.2])/4A denote a local minimum it coincides with [v.sub.f] whereas should denote a local maximum it coincides with [v.sub.m]. 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:
[r.sub.cr] = [r.sub.f] * [r.sub.i] (15)
Here, the notation * denotes a binary operation over the set of real numbers combining the elements [r.sub.f], [r.sub.i] to produce univocally another real number.
[v.sub.i] ([r.sub.cr]) = [v.sub.f] x [v.sub.m] (16)
Similarly, the notation x designates a binary operation over the set of real numbers combining the elements [v.sub.f], [v.sub.m] to produce univocally another real number.
Moreover since [r.sub.f] < [r.sub.i] and [v.sub.f] < [v.sub.m] these binary operations should not be necessarily commutative.
In this context, should the critical value [r.sub.cr] yielded by the operation [r.sub.f] * [v.sub.f] x [v.sub.m] be on the interval [[r.sub.f], [r.sub.i] it is necessary to require the corresponding value of Poisson's ratio [v.sub.i] ([r.sub.cr]) = [v.sub.i] ([r.sub.f] * [r.sub.i]) to be strictly greater than [v.sub.f] and strictly less than [v.sub.m].
On the other hand, should the value of Poisson's ratio [v.sub.i] ([r.sub.cr]) = [v.sub.f] x [v.sub.m] be strictly less than [v.sub.f] or strictly greater than [v.sub.m] it is necessary to require the value [r.sub.cr] to be outside the interval [[r.sub.f], [r.sub.i]].
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  let us introduce a general thirddegree polynomial to approximate the interphase Poisson's ratio given as
[v.sub.i] (r) = [Ar.sup.3] + [Br.sup.2] + Cr + D (17)
Evidently, (17) signifies a more complicated mode of variation for interphase Poisson's ratio when compared with (4).
Yet, it is known from Calculus  that if [B.sup.2] - 3AC [less than or equal to] 0 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 [v.sub.i] (r).
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 [B.sup.2] - 3AC is strictly positive.
Thus we can write out
d/[d.sub.r] [v.sub.i] (r) = 3[Ar.sup.2] + 2f?r + C (18)
Provided that [B.sup.2] - 3AC > 0 the critical points are given as
[r.sub.cr,1] = -Br - [square root of [B.sup.2] - 3AC]/A; (19a)
[r.sub.cr,2] = -Br + [square root of [B.sup.2] - 3AC]/A; (19a)
Thus in a quite analogous reasoning with that resulting in (15), (16) one may require the validity of the following additional conditions:
[r.sub.cr,1] = [r.sub.f] * [r.sub.i]; (20a)
[r.sub.cr,2] = [r.sub.f] <> [r.sub.i] (20b)
[v.sub.i] ([r.sub.cr,1]) = [v.sub.f] x [v.sub.m]; (21a)
[v.sub.i] ([r.sub.cr,2]) = [v.sub.f] o [v.sub.m] (21b)
Apparently the above disjunction is exclusive. Now, should even one of the critical values [r.sub.cr,1]; [r.sub.cr,1] yielded by the operations [r.sub.f] * [r.sub.i] and [r.sub.cr,2] = rf <> [r.sub.i], respectively, lie on the interval [[r.sub.f], [r.sub.i]] it is necessary to require the corresponding value of Poisson's ratio to be strictly greater than [v.sub.f] and concurrently strictly less that [v.sub.m].
On the other hand, should even one the values of Poisson's ratio [mathematical expression not reproducible] be strictly less than [v.sub.f] or strictly greater than [v.sub.m] it is necessary to require the corresponding critical point to lie outside the interval [[r.sub.f], [r.sub.i]].
Moreover, on the basis of (17), it can be proved that the following relationship holds:
[mathematical expression not reproducible] (21c)
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
[v.sub.i] (r) = A ([v.sub.m] - [v.sub.f]) [(r - [r.sub.f]).sup.n] [([r.sub.i] - r).sup.n] + Br + C (22)
with n [member of] [N.sup.*] and A, B, C being arbitrary real numbers
An application of the boundary conditions expressed by (2) and (3) at the interfaces with fiber and matrix, respectively, yields
[Br.sub.f] + C = [v.sub.i] ([r.sub.f]) = [v.sub.f] (23)
[Br.sub.i] + C = [v.sub.i]([r.sub.i]) = [v.sub.m] (24)
B = [v.sub.m] - vf/[r.sub.i] - [r.sub.f] (25)
C = [v.sub.f] - [v.sub.m] - [v.sub.f]/[r.sub.i] - [r.sub.f] [r.sub.f] (26)
In the sequel, let us calculate the first derivative of [v.sub.i](r) with respect to polar radius
[mathematical expression not reproducible] (27)
[mathematical expression not reproducible] (28)
Here [M.sub.Ar] ([r.sub.f], [r.sub.i]) denotes the arithmetic mean of [r.sub.f], [r.sub.i].
Equation (28) can be combined with (25) to yield
[mathematical expression not reproducible] (29)
At r = [r.sub.c,r] one obtains
[mathematical expression not reproducible] (30)
By focusing on the above representation one may point out that the critical points [r.sub.cr] at which [dv.sub.i](r)/dr vanishes cannot coincide with the endpoints [r.sub.f], [r.sub.i] 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
[mathematical expression not reproducible] (31)
Now, given that the coefficient B has been already calculated, the following relation links the critical points [r.sub.cr], where [dv.sub.i](r)/dr vanishes, with the coefficient A
[mathematical expression not reproducible] (32)
Equation (32) can be combined with (22) to yield
[mathematical expression not reproducible] (33)
Solving for [v.sub.i] ([r.sub.cr]) one obtains
[mathematical expression not reproducible] (34)
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
[mathematical expression not reproducible] (35)
Here one observes that the set of local extrema of interphase Poisson's ratio has been evaluated without the knowledge of the coefficient A 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  it was assumed in a rather subjective and very simplified mode that [v.sub.i]([r.sub.cr]) = [v.sub.f] v [v.sub.i]([r.sub.cr]) = [v.sub.m] whilst [r.sub.cr] = [r.sub.i] v [r.sub.cr] = [r.sub.f].
Apparently, such a consideration is rather unrealistic. However, a reasonable conjecture for the calculation of coefficient A 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 r = [r.sub.cr] by means of other types of means as discussed in  since it has been proved that meantheory is very useful both from theoretical point of view and for practical (engineering) purposes. Yet, such an approach cannot include the case when [r.sub.cr] < [r.sub.f] or [r.sub.cr] > [r.sub.i] and then one may examine the interval [0, [r.sub.m]]. 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 [r.sub.cr].
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 . This function for n = 2 yielded a parabolic law. Also, to assume that at r = [r.sub.i], i.e. at the interface between matrix and interphase, [v.sub.i]([r.sub.i]) = [v.sub.m] is realistic indeed. Moreover, the assumption that at r = [r.sub.f], i.e., at the interface between fiber and interphase [v.sub.i](rf) = [v.sub.f] 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 [[r.sub.f], [r.sub.i]] that was adopted in  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 [DELTA]r = [r.sub.i] - [r.sub.f] and given that [r.sub.f] is known beforehand.
To illustrate the physical meaning of [DELTA]r, 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  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  that the extent of interphase expressed by its thickness [DELTA]r 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 [DELTA]r by the next empirical relationship:
([r.sub.f] + [DELTA]r).sup.2]/[r.sup.2.sub.f] - 1 = [lambda][U.sub.f]/1 - [U.sub.f] (36)
where the coefficient [lambda] is evaluated as
[lambda] = 1 - [DELTA][C.sup.f.sub.p]/[DELTA] [C.sup.0.sub.p] (37)
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  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.
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 [DELTA]r, 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 threephase 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.
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
[v.sub.i] (r) = A([r.sup.3] + B/A [r.sup.2] + C/A r + D/A) (A.l)
[v.sub.i] (r) = A([r.sup.3] + [B.sub.i] [r.sup.2] + [C.sub.1] r + [D.sub.1]) (A.2)
Here, for facility reasons we have put
[mathematical expression not reproducible] (A.3)
Now, by the use of the trivial substitution r = y - [B.sub.1]/3 it implies that
[mathematical expression not reproducible] (A.4)
where h : R [right arrow] R is a single-valued continuous function.
In addition, note that y [member of] [[r.sub.f] - [B.sub.1]/3, [r.sub.i] - [B.sub.1]/3j.
Thus we can write out
[mathematical expression not reproducible] (A.5)
[mathematical expression not reproducible] (A.6)
Equation (A.6) can be combined with (A.3) to yield
[mathematical expression not reproducible] (A.7)
Now, for facility reasons, let us set
[mathematical expression not reproducible] (A.8)
Thus one may deduce that
h (y) = a[y.sup.3] + py + q (A.9)
The first derivative with respect to auxiliary variable y is given as
dh (y)/dy = 3a[y.sup.2] + p (A.10)
Evidently, the necessary and sufficient condition in order for the quantity dh(y)/dy to have real roots is the coefficients a and p to disagree in sign. In this framework, one also infers that these real roots should disagree in sign too. Also, since [d.sup.2] h(y)/[dy.sup.2] = 6ay, it implies that the polynomial h(y) has a unique local minimum and a unique local maximum.
Let [y.sub.1]; [y.sub.2] be the real roots of (A.10):
Evidently, [mathematical expression not reproducible] (A.11)
Then without violating the generality we can write out
[mathematical expression not reproducible] (A.12)
[mathematical expression not reproducible] (A.13)
max h (y) x min h (y)
[mathematical expression not reproducible] (A.14)
max h (y) * min h (y)
[mathematical expression not reproducible]
In continuing, the above relation can be combined with (A.8) to yield
max h (y) * min h (y)
[mathematical expression not reproducible] (A.15)
Finally, since [v.sub.i] (r) [equivalent to] h(y) one may deduce that
[mathematical expression not reproducible] (A.16)
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.
 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-Olive and S. Olive, Eds., vol. 11, Springer Berlin Heidelberg, Berlin, Heidelberg, 1987.
 R. E. Swain, K. L. Reifsnider, K. Jayaraman, and M. elZein, "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. Erdelyi, Polynomials and Polynomial Inequalities, vol. 161 of Graduate Texts in Mathematics, Springer, New York, NY, USA, 1995.
 M. Raissouli and M. Chergui, "On Some Inequalities Involving Three or More Means," Abstract and Applied Analysis, vol. 2016, Article ID 1249604, 8 pages, 2016.
J. Venetis [ID] (1) and E. Sideridis [ID] (2)
(1) School of Applied Mathematics and Physical Sciences NTU Athens, Section of Mechanics, 5 Heroes of Polytechnion Avenue, 15773 Athens, Greece
(2) 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; firstname.lastname@example.org
Received 3 July 2018; Revised 23 August 2018; Accepted 5 September 2018; Published 1 October 2018
Academic Editor: Xin-Lin Gao
Caption: Figure 1: Coaxial three-phase cylinder unit cell.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Venetis, J.; Sideridis, E.|
|Publication:||Journal of Applied Mathematics|
|Date:||Jan 1, 2018|
|Previous Article:||Bridging the Gap between Economic Modelling and Simulation: A Simple Dynamic Aggregate Demand-Aggregate Supply Model with Matlab.|
|Next Article:||Numerical Solution to Coupled Burgers' Equations by Gaussian-Based Hermite Collocation Scheme.|