Effective Complex Properties for Three-Phase Elastic Fiber-Reinforced Composites with Different Unit Cells

: The development of micromechanical models to predict the effective properties of multiphase composites is important for the design and optimization of new materials, as well as to improve our understanding about the structure–properties relationship. In this work, the two-scale asymptotic homogenization method (AHM) is implemented to calculate the out-of-plane effective complex-value properties of periodic three-phase elastic ﬁber-reinforced composites (FRCs) with parallelogram unit cells. Matrix and inclusions materials have complex-valued properties. Closed analytical expressions for the local problems and the out-of-plane shear effective coefﬁcients are given. The solution of the homogenized local problems is found using potential theory. Numerical results are reported and comparisons with data reported in the literature are shown. Good agreements are obtained. In addition, the effects of ﬁber volume fractions and spatial ﬁber distribution on the complex effective elastic properties are analyzed. An analysis of the shear effective properties enhancement is also studied for three-phase FRCs. Author Contributions: Conceptualization, F.J.S., R.R.-R.; methodology, Y.E.-A., H.C.-M.; software and validation, R.G.-D., Y.E.-A.; formal analysis and investigation, F.J.S.,


Introduction
Multiphase elastic fiber-reinforced composites (FRCs) are still important in applications because their yields exceed those of their constituents and they offer very interesting properties compared to more conventional materials. Improved levels of functionality can be achieved through the manipulation of physical properties by means of structure optimization [1][2][3]. Therefore, the effective properties prediction for FRCs by means of micromechanical models and numerical approaches is a useful tool for technological innovation [4][5][6][7][8][9][10].
Periodic multiphase elastic FRCs have found applications related to transport problems (conductivity, shear elasticity, dielectric constant, thermal expansion, and others). In this sense, different micromechanical and experimental models have been developed to analyze elastic FRCs. For example, the elastic effective properties of two-phase elastic FRCs with periodic square [11] and hexagonal [12] cells were found by applying the asymptotic homogenization method (AHM). Jiang et al. [13] calculated the effective elastic moduli of FRCs with cylindrical inclusions under longitudinal shear by means of the Eshelby equivalent inclusion concept [14], combining the results of a doubly quasiperiodic Riemann boundary value problem [15]. Artioli et al. [16] applied the AHM to determine the effective longitudinal shear properties of elastic FRCs with radially graded fibers, assuming imperfect interface conditions. Shu and Stanciulescu [17] proposed an analytical approach using multiscale homogenization to characterize FRCs with imperfect interphase through the shear-lag model. Dhimole et al. [18] implemented a multiscale modeling based on homogenization method to predict the accurate mechanical behavior of 3D four-directional braided composites.
In addition, Bisegna and Caselli [19] investigated the effective complex conductivity for a periodic FRC with interfacial impedance and hexagonal symmetry by AHM. Godin computed the complex effective dielectric properties for two-phase FRCs with circular inclusions [20] and for periodic tubular structures [21]. These tubular structures were modeled as a three-phase FRC. Analytic bounds on the complex dielectric constant were reported by the authors of [22]. A correspondence between orthotropic complex-value dielectric media and non-orthotropic elastic media was developed by the authors of [23] through an affine transformation. Mackay and Lakhatia [24] reported the gain and loss enhancement for particulate homogenized composite materials whose active constituents have complex values. Guild et al. [25] analyzed the enhancement of homogenized dielectric effective properties for acoustic waves using multiscale sonic crystals. Luong et al. [26] estimated the complex shear modulus using the least mean square/algebraic helmholtz inversion (LMS/AHI) algorithm for 1D heterogeneous tissue. On the other hand, fibermatrix interaction region has also been studied for multiphase FRCs. In this case, a thin mesophase is added between the fiber and the matrix in a three-phase FRC. This contact zone is commonly defined as imperfect contact region, see, for instance, [17,[27][28][29][30].
The AHM has proved to be advantageous in the description of the multiscale mechanics of composite materials. Many studies have dealt with the theoretical bases of the AHM [31][32][33][34][35]. In general, the AHM makes it possible to obtain an effective characterization of the heterogeneous system or phenomenon under study by encoding the information available at the microscale into the so-called effective coefficients. In particular, multiscale AHM take advantage of the information available at the smaller scales of a given heterogeneous medium to predict the effective properties at its larger scales, see, for instance, [36][37][38][39]. This, in turn, dramatically reduces the computational complexity of the resulting boundary problems. However, the main disadvantage of AHM is that the analytical solution of the local problems has been derived for only a few composite structures [40,41].
The main aim of this work is the estimation by AHM of the effective elastic complexvalues properties for periodic three-phase elastic FRC with complex-valued constituent properties and a parallelogram cell. The out-of-plane case for three-phase composite is considered. Both matrix and fibers have isotropic or transversely isotropic properties, and they are in welded contact. The mathematical statement is presented, and the associated local problems are derived. Simple closed formulas are provided for the shear effective coefficients of three-phase FRCs. Validations of the present model with results reported in literature are shown. The AHM accuracy and convergence is analyzed based on the truncations of the infinite system from the local problems solutions. Also, the effect of volume fraction and spatial fiber distribution on the complex effective elastic properties is analyzed. An example of longitudinal shear enhancement is considered as a function of reinforcement volume fractions for three-phase FRCs with complex-value constituents. Good agreements are obtained.
The novelty of this work is to calculate by AHM the out-of-plane effective properties of periodic three-phase elastic FRC with parallelogram unit cell whose constituents have complex-valued properties. It generalizes earlier works in which the same method has been applied to two-and three-phase FRCs with real-values constituents and square, hexagonal, and parallelogram unit cells, see, for instance, [11,12,40,42]. This generalization allows the study of the shear effective properties enhancement for three-phase FRCs.

Mathematical Formulation for the Elastic Heterogeneous Media
A heterogeneous periodic three-phase elastic FRC Ω ⊂ R 3 (fiber/mesophase/matrix) with a doubly periodic microstructure is studied, which consisted of a parallelepiped array of two concentric and unidirectional cylinders within a homogeneous matrix (Figure 1a). The fibers are infinitely long in the Ox 3 −direction and periodically distributed. At the microstructural level, the composite cross-section (periodic unit cell Y) is defined by a matrix with two concentric circles of different radii located at the parallelogram center, see Figure 1c. The constituent elastic properties, belonging to a crystal symmetry class of 6mm, are assumed to be complex numbers. In addition, as a unidirectional FRC, the composite microstructure is considered to remain constant along the reinforcement's direction (i.e., perpendicular to Oy 1 y 2 −plane of cross-section).

Mathematical Formulation for the Elastic Heterogeneous Media
A heterogeneous periodic three-phase elastic FRC 3 Ω ⊂  (fiber/mesophase/matrix) with a doubly periodic microstructure is studied, which consisted of a parallelepiped array of two concentric and unidirectional cylinders within a homogeneous matrix ( Figure  1a). The fibers are infinitely long in the 3 Ox − direction and periodically distributed. At the microstructural level, the composite cross-section (periodic unit cell Y ) is defined by a matrix with two concentric circles of different radii located at the parallelogram center, see Figure 1c. The constituent elastic properties, belonging to a crystal symmetry class of 6mm, are assumed to be complex numbers. In addition, as a unidirectional FRC, the composite microstructure is considered to remain constant along the reinforcement's direction (i.e., perpendicular to 1 2 Oy y − plane of cross-section). The periodic parallelogram-like unit cell Y is characterized by a principal angle θ and the periods 1 ω and 2 ω in the 1 2 O y y − plane. Also, it is satisfied that In the out-of-plane mechanical deformation state, the mechanical displacement ( ) of a media Ω satisfies that 1 1 2 the non-null stresses involved in this problem are 13 1 2 ( , ) x x σ and 23 1 2 ( , ) x x σ . Thus, the static governing equation for an elastic FRC Ω is defined by the partial differential equations system: where the absence of body forces is considered. Here, ε are Y-periodic and rapidly oscillating coefficients which denotes the elastic stiffness modulus, and , 1,2 = η β . The Equation (1) subject to the prescribed boundary conditions: The periodic parallelogram-like unit cell Y is characterized by a principal angle θ and the periods ω 1 and ω 2 in the Oy 1 y 2 −plane. Also, it is satisfied that Y = 3 ∪ γ=1 S γ and S i ∩ S j = ∅, i = j, where S 1 is the region occupied by matrix phase defined by a parallelogram Σ with a central circular hole of radius R 1 , contour Γ 1 , and volume V 1 . S 2 represents the mesophase of volume V 2 and is surrounded by two circular interfaces Γ 1 and Γ 2 of radius R 1 and R 2 (R 1 > R 2 ), respectively, and S 3 is the central fiber with circular boundary Γ 2 of radius R 2 and volume V 3 . The volume V of the cell satisfies that In the out-of-plane mechanical deformation state, the mechanical displacement and the non-null stresses involved in this problem are σ 13 (x 1 , x 2 ) and σ 23 (x 1 , x 2 ). Thus, the static governing equation for an elastic FRC Ω is defined by the partial differential equations system: where the absence of body forces is considered. Here, C η3β3 (x/ε) are Y-periodic and rapidly oscillating coefficients which denotes the elastic stiffness modulus, and η, β = 1, 2.
The Equation (1) subject to the prescribed boundary conditions: represent the out-of-plane classical elliptic boundary value problem associated with the linear elasticity theory for heterogeneous structures. Here, t 0 (x) is the traction vector on ∂Ω 1 , g 1 (x) is an infinitely differentiable function on the external boundary ∂Ω 2 , and n j is the component of the outward unit normal vector on ∂Ω 1 . The boundary of the composite is partitioned in such a way that ∂Ω = ∂Ω 1 ∪ ∂Ω 2 and ∂Ω 1 ∩ ∂Ω 2 = ∅. In Equations (1)- (3), σ 3j and u 3 are the out-of-plane mechanical stresses and displacements. In addition to Equations (1)-(3), perfect contact conditions at the circular interfaces Γ s (s = 1, 2) are assumed, i.e., means jump of f across Γ s , see for instance [11].

Method of Solution: Local Problems and Effective Coefficients
The two-scale AHM [32,33,43] is used to solve the elliptic boundary value problem (Equations (1)- (4)). The solution is found by means of a two-scale asymptotic representation of u 3 in powers of the small geometrical parameter ε by the ansatz: where the macroscale or fast variable x and the microscale or slow variable y are related by x = εy, and the parameter ε = l/L is the ratio between the periodic unit cell length (l) and the macroscopic dimension of the composite (L). In Equation (5), the second term u 3 (1) (x, y) is a periodic function of y, which represents a correction of the term which is the solution of the local problems. It is possible to obtain an asymptotic solution of the problem (Equations (1)-(4)) when ε → 0 . More details and the rigorous theoretical background of AHM have been described in classical works [31,33,43] and are omitted here.
The out-of-plane local problems on Y, denoted as α3 L (α = 1, 2) for a three-phase elastic FRC (see Figure 1), is stated as follows: where δ 1α and δ 2α are the Kronecker's delta functions related to the 13 L and 23 L local problems, respectively. To guarantee the uniqueness of the local problems solutions, the local functions should satisfy the null average condition α3 Once the α3 L out-of-plane local problems (Equations (6)-(8)) are solved, the corresponding effective elastic coefficients can be calculated by the formulas: Notice that the out-of-plane effective elastic coefficients (Equations (9) and (10)) depend on the local functions 13 N(y) and 23 N(y) relative to the 13 L and 23 L local problems, respectively. Then, 13 N(y) and 23 N(y) need to be found. Therefore, an analytical solution of Equations (6)-(8) is determined.
On the other hand, the homogenized elastic problem equivalent to the boundary value problem (Equations (1)- (3)) is defined by the equation system subject to the homogenized boundary conditions where u (0) 3 is the solution of the homogenized problem, and C * 3α3β are the out-of-plane effective elastic coefficients defined in Equations (9) and (10).

Analytical Solution of the Local Problems
By means the potential methods of a complex variable theory, the solution α3 N of the α3 L local problem (Equations (6)-(8)) is calculated. Here, the doubly periodic Weierstrass' elliptic functions are used to obtain an analytical solution, i.e., the double periodic solution α3 N is found in terms of Laurent and powers expansions as a function of z = y 1 + iy 2 , see, for instance, [44,45], as follows where ,m 2 + n 2 = 0, k + p > 2, and k, p = 1, 3, 5, · · · . By and In Equations (13)-(15), the coefficients α3 a 0 , α3 a p , α3 b p , α3 b −p , and α3 c p are complex and undetermined numbers. They need to be determined in order to find the α3 L local problems solution and the out-of-plane effective elastic coefficients (Equations (9) and (10)). Here, it can be highlighted that the summation symbol with superscript ∑ o means that the sum only runs over odd integers, and the symbols Re and Im represent the real and imaginary parts of complex numbers, respectively. Details of Laurent and powers expansions and its relationship with the double periodic elliptic Weierstrass function ℘(ω 1 , ω 2 ; z) of periods ω 1 and ω 2 , and related expressions can be found in [45,46].
From the double periodicity condition of α3 N, it is satisfied that: see, for instance, [45]. Then, it can be proved that α3 a 0 is linked with α3 a 1 by the equation: where H 1 = δ 1 w 2 − δ 2 w 1 /(w 1 w 2 − w 2 w 1 ), H 2 = (δ 1 w 2 − δ 2 w 1 )/(w 1 w 2 − w 2 w 1 ) and α3 a 1 is the complex conjugate of α3 a 1 . In addition, δ α = 2ζ(w α /2) is the quasiperiodic condition and ζ(z) is the Weierstrass quasiperiodic Zeta function defined by ζ(z) means that the summation does not include the point (0, 0), see [45]. Finally, replacing Equations (13)-(15) into the interface conditions (Equations (7) and (8)) and after some mathematical manipulations, we obtain the normal infinity system of linear equations [47] with unknown complex constants α3 a p , in compact form: where In Equation (18), W kp = pk −1 η kp , V 2 and V 3 are the volume fractions of the mesophase and the central fiber, and the symbol δ 1p is the Kronecker's delta function. Also, the system solution of each local problem α3 L (α = 1, 2) depends on the elastic constituent properties, the phase volume fractions and the fibers spatial distribution within the matrix. Details of the system construction can be found in [40] and are omitted here. Once the unknown constants α3 a p (p = 1, 3, 5, · · · ) are calculated, the local problem solution and the effective coefficients can be determined. Details on the system solution is given in Appendix A.
The non-null out-of-plane effective properties for three-phase elastic FRC are listed as follows: where C 1313 = C (1) 1313 V 3 is the Voigt's average, A 21 = R 2 1 (k 1 + k 2 ) +R 2 2 (k 1 − k 2 ) /2R 1 R 2 , and the constant is defined by Simple closed-form formulas for the effective properties equivalent to Equations (19) and (20) are given in Appendix B.

Numerical Results
The accuracy of the AHM model is determined through comparisons with other results reported in the literature for three-phase elastic FRC (fiber/mesophase/matrix) with complex-values constituent properties. In addition, the effect of the volume fraction of inclusions and the spatial fiber distribution on the complex effective elastic properties is investigated. Also, an example of shear enhancement is reported as a function of reinforcement volume fractions.
Limit cases for the present model can be determined when Equations (19) and (20) and system (18) are reduced to those that represent a two-and three-phase FRC with parallelogram-like unit cell. In these cases, isotropic or transversely isotropic constituents are considered with real-values elastic properties, as reported by the authors of [11,40,42].
The real and imaginary parts of the complex effective elastic properties are given for some biological applications, for example, the behaviors of biological tissues, skeletal muscle, sclera, and other ones in which material properties depends on time [48][49][50][51]. This behavior can be analyzed in viscoelastic materials in which the shear wave speed is connected to the enhancement and lessening of the shear modulus. Besides, in the context of transport properties, the mathematical statement for shear linear responses is identical to conductivity, dielectric permittivity, and so on in equivalent media [52].
In the literature, to the authors' knowledge, the longitudinal shear homogenization problem has not been reported with complex-values coefficients. A comparison with Godin [21] is possible because the governing equation for both models has the same mathematical formula, although they model different physical magnitudes. Both models find an infinity equation system to solve the problem. However, there are differences between the model implementations. The Godin model goes directly to a representative element of analysis and the governing equation is directly solved. The effective properties are proposed as a linear relation between averaged physical magnitudes. The AHM described herein in previous sections is based on a procedure that ends up in the solutions of local problems, the boundary conditions, and the effective properties.
In Table 1, the real (Re) and imaginary (Im) parts of complex effective property C * 1313 are illustrated for two different three-phase FRCs as a function of normalized radius h = R 1 /l with a square and hexagonal unit cell, respectively. Here, h ≤ 1/2, R 1 is the outer interface radius, and l represents the minimum distance between the centers of the fibers. In addition, an analysis of the relative error is shown: Error = [(AHM − Godin)/AHM] × 100% [21]. Hence, the numerical values of the effective property C * 1313 is compared with the complex effective dielectric constant ε * provided by the authors of [21]. The numerical calculations are carried out considering that the matrix, the mesophase, and the fiber have 1313 = 2 − 4i for the FRC with the hexagonal array. In both three-phase FRCs, the concentric fibers radius relation is R 2 2 /R 2 1 = 0.81. As can be seen, good agreements are achieved by AHM for only a system truncation (Equation (18)) to a finite order N 0 = 2. Only a very slight discrepancy is observed near the close-packing condition, although the relative error is below 0.27% in every case. Table 1. Real and imaginary parts of the complex effective property C * 1313 as a function of normalized radius for two three-phase FRCs with a square and hexagonal unit cell, respectively.
Godin [21] Error (%) AHM Godin [21] Error (%) Godin [21] Error (%) AHM Godin [21] Error (%) In Table 2, the variations of the real and imaginary parts of the complex effective shears C * 1313 , C * 1323 , and C * 2323 in terms of system truncate orders N 0 (N 0 = 1, 2, 3, 5, 7, 9, 11) are presented for a three-phase elastic FRC with a parallelogram unit cell of θ = 75 • . In addition, two different normalized radii h = 0.4 and h = 0.499 (close to percolation value) are also considered. For the analysis, the material properties of matrix, mesophase, and fiber have isotropic complex properties, such as C The relation between the concentric fiber's radius is R 2 2 /R 2 1 = 0.81. Notice that the AHM convergence is achieved quickly when low values of h ≤ 0.4 are assumed, i.e., only smaller values of N 0 ≤ 3 are needed. Therefore, truncations of higher order N 0 must be considered for high values of h, as well as for higher contrast among the matrix mesophase and fiber properties, in order to obtain a better accuracy. For example, N 0 ≥ 9 is required to obtain the effective properties values with at least five accuracy digits when h = 0.499. Similar conclusions are achieved when the same analysis is developed for a three-phase elastic FRC with a parallelogram unit cell of θ = 75 • and complex-values constituent properties C (1) Table 2. Real and imaginary parts of the complex effective shears C * 1313 , C * 1323 , and C * 2323 obtained by AHM in term of order system N 0 for a three-phase FRC with a parallelogram unit cell of θ = 75 • and for two different normalized radii h = 0.4 and h = 0.499.

N 0
Effective complex shears C * 1313 , C * 1323 and C * 2323 (in GPa) for a three-phase FRC with parallelogram unit cell of θ = 75 • and h = 0.4. Effective complex shears C * 1313 , C * 1323 , and C * 2323 (in GPa) for a three-phase FRC with a parallelogram unit cell of θ = 75 • and h = 0.499.  Table 3 shows the real and imaginary parts of the overall out-of-plane shear properties C * 1313 , C * 1323 , and C * 2323 for two three-phase FRCs with different parallelogram-like unit cells. The numerical values are computed considering four different parallelogram-like unit cells (i.e., parallelogram cells characterized by a principal angle θ equal to 45 • , 60 • , 75 • , and 90 • ), a system order truncation N 0 = 10, and h = 0.38 (value of the normalized radius near to the percolation point volume 0.38268, for 45 • ). In addition, the composite structure-property relationship is also analyzed. From Table 3, it is noticed that, when the periodic unit cells are characterized by parallelograms with θ different to (60 • ) and (90 • ), the composites belong to monoclinic symmetric class, i.e., 13 non-null effective elastic constants are attained. However, in the out-of-plane case, only C * 1313 = C * 2323 and C * 1323 = C * 2313 = 0 are remained. In the case of the periodic hexagonal (60 • ) and square (90 • ) unit cells, the composite behavior is transversely isotropic, i.e., C * 1313 = C * 2323 and C * 1323 = C * 2313 = 0. These results have also been satisfied in elastic FRCs with real effective properties, see, for instance, [40]. In addition, it can be concluded that a decrease (increase) in the real (imaginary) part of the complex effective shears C * 1313 and C * 2323 resulted as the angle of the periodic unit cell increased. Higher values for the real and imaginary parts of C * 1313 and C * 2323 are obtained when θ = 30 • and the normalized radius is the same. Besides, for the composite with complex-value constituents C (1) 1313 = 80 − 2i, and C (3) 1313 = 2 − 4i, the real and imaginary parts of C * 1323 are negative for 30 • ≤ θ < 60 • and positive for 60 • < θ < 90 • . In a composite with complex-value constituents C  Table 3. Real and imaginary parts of the complex effective shears C * 1313 , C * 1323 , and C * 2323 for two three-phase FRCs with different parallelogram-like unit cells and a normalized radius h = 0.38.
θ Real and imaginary parts of the effective complex shears C * 1313 , C * 1323 and C * 2323 (in GPa) of a three-phase FRC with constituent properties C (1) Real and imaginary parts of the effective complex shears C * 1313 , C * 1323 , and C * 2323 (in GPa) of a three-phase FRC with constituent properties C (1) In Figure 2, an analysis of enhancement of the real part of the shear effective property C * 1313 is illustrated for a three-phase FRC with periodic hexagonal unit cell. Notice that, for this type of unit cell, C * 2323 = C * 1313 and C * 1323 = C * 2313 = 0. Here, the enhancement is studied as function of reduced mesophase filling fraction V 2 , i.e., the enhancement of C * 1313 i obtained for four configurations of the annular inclusion (mesophase and fiber inclusions) with fixed volume fraction V 2 + V 3 . For this analysis, two different three-phase FRCs are considered with complex-values constituents C  (Figure 2b). These properties are considered assuming that the mesophase properties attain an interval of realistic physical properties as a combination of the matrix and fiber phases. From Figure 2, it is important to note that the real part of C * 1313 increase as V 2 + V 3 increase when 0 < V 3 < 1, and that higher values of C * 1313 are always obtained in comparison when V 3 = 0 (two-phase FRCsolid black line). In the figures, the red dashed dot and the blue dashed lines represent the real part values of the matrix and mesophase properties, respectively. In addition, it can be concluded that the imaginary part values of both three-phase FRCs vary monotonically between the values of the matrix and mesophase phases. This picture, in which the threephase FRCs are considered with distinct constituents, is different than Figure 2 reported by the authors of [7]. In this reference, the enhancement was analyzed for a three-phase FRC with the same values of matrix and fiber and different mesophase, as an annular inclusion. Besides, it should be noted that the volume fraction interval where enhancement appeared is much larger. In one case, it is the whole interval.

Conclusions
In this work, the effective shear properties of periodic three-phase fiber-reinforced composites with complex-values constituent properties and parallelogram unit cells are calculated by AHM. Easy-to-handle formulas and fast numerical implementation are derived for all shear effective properties. We conclude that: (i) The fiber spatial distribution, represented as parallelogram-like unit cell, is capable of describing three class of symmetry point group: tetragonal 4 mm (square unit cell), hexagonal 6 mm (hexagonal unit cell), and monoclinic 2 (other parallelogram unit cells) structures.
(ii) The enhancement in the shear effective property * 1313 C is more remarkable for threephase FRC than two-phase FRC. (iii) The volume fraction interval where enhancement appeared is larger for a three-phase FRC than the interval for the two-phase FRC.
(iv) The presence of negative values for the real and imaginary parts of * 1323 C appears for some parallelogram unit cells. (v) The manipulation of the mesophase can be used as a way to enhance the real and imaginary parts of the shear elastic properties. (vi) The numerical results prove that the AHM is an accurate and efficient approach for the study of FRC with a mesophase and for different spatial fiber distributions in a matrix.

Conclusions
In this work, the effective shear properties of periodic three-phase fiber-reinforced composites with complex-values constituent properties and parallelogram unit cells are calculated by AHM. Easy-to-handle formulas and fast numerical implementation are derived for all shear effective properties. We conclude that: (i) The fiber spatial distribution, represented as parallelogram-like unit cell, is capable of describing three class of symmetry point group: tetragonal 4 mm (square unit cell), hexagonal 6 mm (hexagonal unit cell), and monoclinic 2 (other parallelogram unit cells) structures. (ii) The enhancement in the shear effective property C * 1313 is more remarkable for threephase FRC than two-phase FRC. (iii) The volume fraction interval where enhancement appeared is larger for a three-phase FRC than the interval for the two-phase FRC. (iv) The presence of negative values for the real and imaginary parts of C * 1323 appears for some parallelogram unit cells.
(v) The manipulation of the mesophase can be used as a way to enhance the real and imaginary parts of the shear elastic properties. (vi) The numerical results prove that the AHM is an accurate and efficient approach for the study of FRC with a mesophase and for different spatial fiber distributions in a matrix.