Identiﬁcation of the Physical and Mechanical Properties of Moroccan Sisal Yarns Used as Reinforcements for Composite Materials

: This work aims to investigate the physical and mechanical properties of sisal ﬁber and yarn of Moroccan origin. The cellulosic and non-cellulosic constituents of the Moroccan sisal ﬁber were identiﬁed by FTIR spectroscopy. The thermal properties were studied by thermogravimetric analysis. The hydrophilicity of the ﬁber was evaluated by the contact angle. The results show that the sisal ﬁber has a low thermal stability. The mechanical properties of the ﬁber analyzed by the Impregnated Fiber Bundle Test (IFBT) method show that the porosity of the impregnated yarns and the twist angle of the yarns inﬂuence the elastic modulus of the sisal ﬁber. The physical and mechanical properties of the manufactured sisal yarns were also characterized and analyzed. The obtained results reveal an interesting potential to use the Moroccan sisal ﬁber in development of bio-sourced composite materials. The thermogravimetric analysis of sisal ﬁber shows a low thermal stability. The presence of the hydroxyl groups (OH) in the structure of the ﬁber favors the hydrophilic behavior of the ﬁber, with a contact angle of 53.06 . The moisture absorption rate of the ﬁber is 13.2%. The mechanical properties of the ﬁber analyzed by the IFBT method shows that the porosity of the impregnated yarns and the twist angle of the yarns inﬂuence the elastic modulus of sisal ﬁber. The physical and mechanical properties of sisal yarns show a variability in the values of diameter and a linear density that depends on the properties of the twisted ﬁbers used in the spinning process. This investigation on the properties of sisal ﬁber of Moroccan origin proves that this ﬁber can be an alternative ﬁber to synthetic ﬁbers because of its recyclability, low eco-friendly and its abundance in several regions, which are highly demanded points in the strategy of sustainable development. The limitations of the ﬁber such as moisture absorption and low thermal stability can be optimized for the development of reinforcements of composite materials.


Introduction
Composite materials are high-performance materials used in various fields of application such as automotive, aeronautics, construction. Despite the advantages of conventional composites such as mechanical strength, corrosion resistance and design flexibility, the limitations of composites are related to the nature of their petroleum source (nonbiodegradability, manufacturing process with high energy consumption . . . ) [1][2][3].
The use of alternative fibers derived from natural resources to synthetic fibers in the manufacture of composite materials reduces the economic and ecological footprint of these materials and meets environmental regulations on industrial activity [4][5][6]. Natural fibers have become more attractive to the composite materials industry, which has increased the world production of these fibers (an increase of 2 million tons in ten years 2008-2018). The automotive sector remains the largest market for natural fiber reinforced composites [5,6]. The use of natural fibers in the automotive industry not only concerns the environmental aspect but also the economic aspect by reducing production costs and resin consumption [6]. The low density of natural fibers compared to the density of glass fibers allows the development of lightweight vehicles, which causes reduction in CO 2 emissions. Different works have focused on a comparative study of composite materials based on synthetic and natural fibers [7][8][9][10]. However, while the use of natural fibers as a reinforcement has shown important advantages, major problems have been identified that are related to the nature of the fiber structure, such as moisture absorption and variability of mechanical performance and odor release due to the emission of oxygenated Volatile Fibers 2021, 9,13 2 of 16 organic compounds (VOC's). These problems can be solved, respectively, by chemical treatment of the fibers and by the incorporation of non-porous inorganic fillers to trap gaseous substitutes [11][12][13][14].
In general, natural fibers can be either of plant origin such as Flax [15], bamboo [16], abaca [17] or animal origin such as wool [18] and silk [19]. According to the location of the fibers in the plant, the fibers are classified as: leaf fiber (such as banana fiber [20], pineapple fiber [21]), bast fiber (such as kenaf fiber [22], ramie fiber [23]) and seed fiber (such as cotton fiber [24]). The different types of plant fibers allow some variability in mechanical performance and it also depends on the growing conditions of the plant.
The performance of natural fibers has been analyzed in different studies in order to describe their structure and identify their specific properties. Murali et al. [25] have studied the physical properties of the natural fibers (bamboo, palm tree date, vakka) by determining the length, diameter and density of the fibers. The results show variability in diameter for all the three fibers. The calculated densities of the three fibers are lower than the density of the glass fibers [25]. The chemical composition of vegetable fiber has been studied by infrared spectroscopy in several works. The rigidity of vegetable fiber depends on the cellulosic component of the fiber [26].
After the analysis of the chemical composition of the fiber, the hydrophilicity of the fiber proves to be an interesting property to study by measuring the contact angle and surface energy to identify the appropriate matrix for each fiber to obtain a good adhesion between the reinforcement and the matrix. Singh et al. [27] show that jute fiber shows hydrophilic behavior with a contact angle less than 90 • and moisture absorption rate of 14%.
Vegetable yarns made by spinning twisted fibers have been studied in the literature. Corbin et al. [28] have compared the properties of flax yarns and rovings. The results obtained show that the twisting of the yarns improves the tenacity of the yarns at breakage, as well as the irregularity of the yarn properties which is due to the spinning process. Kim and Netravali [29] reveal that the breakage of hemp yarns with a tensile strength of 450 MPa and a linear density of 40.74 Tex propagates around the first broken fibers [29].
The mechanical performance of natural fibers is an important factor in the replacement of synthetic fibers. The tensile properties of natural fibers are tested by two methods either with the unit of the tensile test of the fiber or with the back-calculation of the properties of fibers from the results of the impregnated bundles or yarns using the rule of mixtures [30,31]. Lansiaux et al. [32] have studied the mechanical properties of flax fibers by the Impregnated Fiber Bundle Test (IFBT) method by varying the twist angle of the flax yarns. The obtained results are similar for a twist angle between 0 and 100 tpm with the modified rule of mixtures taking into account the porosity and orientation of the fibers. Among the most efficient natural fibers, sisal fiber is a fiber extracted from the leaves of the agave sisalana plant (Table 1). Thanks to the ease of cultivation of the agave sisalana plant, the production of this fiber represents almost half of the production of textile natural fibers [33,34]. The sisal fibers are widely cultivated in the tropical and subtropical regions of north and south America and African countries [34][35][36][37][38]. Belaadi et al. [39,40] studied the properties of sisal fiber of Algerian origin. The results prove that the tensile strength and the Young modulus of the Algerian fibers are superior to the properties of fibers of Indian origin [40]. The growing conditions of each region allow for variability in the performance of the sisal fiber (Table 2). The sisal fiber cultivated in several African countries (precisely in Morocco) is exploited only in traditional medicine. The performance of Moroccan sisal fiber was investigated in a previous study by showing the important mechanical resistance of the Moroccan sisal fiber [41]. This work aims to show the potential of sisal fiber of Moroccan origin for the reinforcement of composite materials at the fiber and yarn scale. This multi-scale approach (fiber/yarn) allows the understanding of the complexity of 2D and 3D fibrous reinforcements. The characterization of sisal yarns is the last step of this study in order to study the smallest entity of the reinforcement. The important factors related to the fiber to assess the adhesion between the matrix and the Moroccan fiber (sisal) were analyzed (hydrophilicity, moisture absorption) according to standards.

Materials
Moroccan-Sisal fibers were extracted mechanically from the leaves of the Agave Sisalana plant. After the fiber extraction process, the fibers were washing with distilled water (2 h) and dried in an oven at 40 • C for 24 h [41].
Sisal yarn (Figure 1) was purchased from Moroccan company Sonajute (reference product (300/1 (m/Kg)). It is a simple twisted yarn with a twist level of 80 tpm and a linear density of 3300 Tex. The sisal fibers are widely cultivated in the tropical and subtropical regions of north and south America and African countries [34][35][36][37][38]. Belaadi et al. [39,40] studied the properties of sisal fiber of Algerian origin. The results prove that the tensile strength and the Young modulus of the Algerian fibers are superior to the properties of fibers of Indian origin [40]. The growing conditions of each region allow for variability in the performance of the sisal fiber ( Table 2). The sisal fiber cultivated in several African countries (precisely in Morocco) is exploited only in traditional medicine. The performance of Moroccan sisal fiber was investigated in a previous study by showing the important mechanical resistance of the Moroccan sisal fiber [41]. This work aims to show the potential of sisal fiber of Moroccan origin for the reinforcement of composite materials at the fiber and yarn scale. This multi-scale approach (fiber/yarn) allows the understanding of the complexity of 2D and 3D fibrous reinforcements. The characterization of sisal yarns is the last step of this study in order to study the smallest entity of the reinforcement. The important factors related to the fiber to assess the adhesion between the matrix and the Moroccan fiber (sisal) were analyzed (hydrophilicity, moisture absorption) according to standards.

Materials
Moroccan-Sisal fibers were extracted mechanically from the leaves of the Agave Sisalana plant. After the fiber extraction process, the fibers were washing with distilled water (2 h) and dried in an oven at 40 °C for 24 h [41].
Sisal yarn (Figure 1) was purchased from Moroccan company Sonajute (reference product (300/1 (m/Kg)). It is a simple twisted yarn with a twist level of 80 tpm and a linear density of 3300 Tex.

FTIR Analysis
Fourier Transform Infrared analysis was conducted to identify organic and inorganic components of the provided sisal fibers. A Thermo Scientific Smart OMNI-Sampler spectrometer was employed with a wavenumber range from 4000 to 400 cm −1 . The thermal behavior of the sisal fiber was investigated by thermogravimetric analysis operated on TGA-Q500 of TA Instruments (ISO 11358). A sample weighing 5.7 mg of the sisal fibers was heated from 20 • C to 600 • C with a heating rate of 10 • C/min in the nitrogen atmosphere.

Contact Angle Measurement
The wettability of sisal fibers was evaluated by calculating the contact angle θ according to tensiometric method by using a GBX Digidrop (Dublin, Ireland). The fiber was immersed in distilled water with a surface tension of 73.13 (mN/m) at room temperature. The meniscus weight (m) and the wetting perimeter (p) were measured. Then the contact angle was deduced by Equation (1) (Wilhelmy [45]).
where F is the capillary Force (N), g is the acceleration of gravity (m.s −2 ), Y L (mN/m) is the surface tension of distilled water, P is the wetting perimeter and θ is the contact angle.

Moisture Absorption Analysis
The moisture absorption analysis consisted of calculating the rate of moisture absorption under the hygrometric conditions indicated in the NF G 08-001-4. A sample weighing 100 g of sisal fibers was dried at 60 • C for 12 h, until a constant mass was reached. Then, the sample was placed in conditioned room (T = 22 • C, RH = 62%) and weighted every 15 min until the mass stabilized again. The moisture content in sisal fiber at time t was expressed by the following Equation (2): where M a (%) is the moisture content of sisal fibers, M o is the initial weight of samples and M t is the weight of samples at time t [46]. Fick's law was used to model the moisture diffusion in a composite material by a diffusion coefficient given in Equation (3). This law can be employed to model the moisture diffusion in the inner structure of sisal fiber since it is considered as composite material regarding its composition, cellulosic and non-cellulosic components (lignin, hemicellulose, pectin): where D is the Diffusion coefficient, b is the diameter of sisal fiber, M m is the rate of moisture absorption and M 1 and M 2 are the moisture content at t 1 and t 2 , respectively [47,48].

Fiber Geometry
The morphology of the fiber depends on the length and diameter of the fiber. The diameter was measured using an optical microscope (Olympus BX51, Tokyo, Japan). Twenty-five samples were analyzed through images in the longitudinal direction of the fiber. The average value of the fiber diameter was calculated. The lengths of twenty-five fibers were measured using a special device called a "maillemeter", where one end of the fiber is grasped by a clamp and the second end is attached to a mobile carriage which exerts a force to obtain the real value of the length of the fiber (without any crimp). A longitudinal view of sisal fiber was observed by an electron-scanning microscope (JEOL JSM 6380LA, Jeol Ltd., Tokyo, Japan).

Impregnated Fiber Bundle Test
The mechanical properties of sisal fibers were determined by the Impregnated Fiber Bundle Test (IFBT). The method consists of making composite samples from sisal yarns and identifying the tensile properties of these samples. Fiber-scale properties are deduced, by inverse computation, thanks to the rule of mixtures [32,33].
Bundles of the sisal fiber were impregnated with an epoxy resin (SR8200/SD7204 sicomin) in a rectangular mold (250 mm × 10 mm). The composite material was pressed using a heating press at 60 • C for 4 h to promote resin curing. Different pressure values were applied in order to investigate the pressure effect on the porosity content in the impregnated samples. Thickness of the composite sample was controlled to 2 mm.
Tensile test was performed on the cured bundle with a speed of 1 mm/min using the Instron 5985 tensile machine with a maximum load of 250 KN. Stiffness of the sisal fibers was determined thanks to the modified rule of mixtures [49,50], as given in Equation (4). In comparison with the classic rule of mixtures, the modified formulae take into account the porosity content in composite samples in addition to the two associated factors, respectively, fiber orientation and fiber length. (4) where: • E f , E c and E m are the stiffness of fibers, composites and the matrix, respectively. • V f , V m and V p are the volume fraction of fibers, the matrix and porosity, respectively. • n is the porosity efficiency exponent. If n = 0, Equation (4) does not consider the level of porosity and will be denoted here as rule of mixtures R1. n = 2 corresponds to the general value of the composites based on natural fibers (denoted here, rule of mixtures R2). • η 0 is the fiber orientation efficiency factor. • η 1 is the fiber length efficiency factor (η 1 =1 [30,48,49]).
The fiber orientation efficiency factor can be expressed either as a function of linear density of yarns and density of sisal fiber, given by Equation (5) (Rule of mixtures R3) or as a function of the twist angle of the sisal yarns [4][5][6], given by Equation (6) (Rule of mixture R4).
From these models, stiffness of the sisal fibers was computed by inverse computation from the properties of the elaborated composites in two ranges of deformation: between 0 and 0.1% for the first modulus, denoted E1f and between 0.3 and 0.5% for the second modulus, denoted E2f.

Diameter and Linear Density Measurement
The yarns were conserved in lab conditions (21 • C, 65% for 24 h). The diameter of sisal yarns was measured using the mini optical microscope. Twenty-five samples of yarns with a length of 40 mm were visualized and the average value of their diameters was calculated. The linear density of sisal yarns was measured according to ISO 1889. Furthermore, the angle of the fiber located on the outer surface of the yarn was measured with the aid of captured pictures taken on the optical microscope and using ImageJ software.

Tensile Analysis
The tensile behavior of sisal yarns was investigated according to NF ISO 2062 by a tensile tester (MTS SYSTEM) with a maximal load 10KN. Twenty specimens of yarns were tested with a gauge length of 100 mm at room temperature.

FTIR spectroscopy of Sisal Fibers
Different spectral bands describing the chemical composition of the sisal fiber are shown in Figure 2. The first band between 3200 and 3500 cm −1 is attributed to an O-H, group showing the presence of hemicellulose and cellulose. The peak identified between 2900 and 3100 cm −1 corresponds to a C-H group bound to cellulose and hemicellulose. These two bands represent both the macromolecular interactions of cellulose and hemicellulose as well as the presence of water in the fibers, both of which promote the binding of water molecules to these characteristic hydroxyl groups. The C-H Stretching vibration region at the band level around 1735 cm −1 is due to the hemi cellulosic component. Shear vibration of the OH-bond illustrates the presence of free water in the fiber structure. Pectin constituting the fiber appears at the vibration of the carboxyl bond between 1430 and 1420 cm −1 . The carboxyl functions of pectins which make up the middle lamella, are also impacted by absorption of water. Stretching of the C-O bonds at the 1225 cm −1 peak identifies the presence of lignin. The cellulose is shown by stretching the C-O-C bond at wavelength 1023 cm −1 which is the characteristic of the polysaccharide nature of the fiber. This signature is also impacted by the increase in humidity [51][52][53]. The infrared spectrum of the fiber shows the ability of the fiber to absorb moisture due to its chemical composition [54].

Diameter and Linear Density Measurement
The yarns were conserved in lab conditions (21 °C, 65% for 24 h). The diameter of sisal yarns was measured using the mini optical microscope. Twenty-five samples of yarns with a length of 40 mm were visualized and the average value of their diameters was calculated. The linear density of sisal yarns was measured according to ISO 1889. Furthermore, the angle of the fiber located on the outer surface of the yarn was measured with the aid of captured pictures taken on the optical microscope and using ImageJ software.

Tensile Analysis
The tensile behavior of sisal yarns was investigated according to NF ISO 2062 by a tensile tester (MTS SYSTEM) with a maximal load 10KN. Twenty specimens of yarns were tested with a gauge length of 100 mm at room temperature.

FTIR spectroscopy of Sisal Fibers
Different spectral bands describing the chemical composition of the sisal fiber are shown in Figure 2. The first band between 3200 and 3500 cm −1 is attributed to an O-H, group showing the presence of hemicellulose and cellulose. The peak identified between 2900 and 3100 cm −1 corresponds to a C-H group bound to cellulose and hemicellulose. These two bands represent both the macromolecular interactions of cellulose and hemicellulose as well as the presence of water in the fibers, both of which promote the binding of water molecules to these characteristic hydroxyl groups. The C-H Stretching vibration region at the band level around 1735 cm −1 is due to the hemi cellulosic component. Shear vibration of the OH-bond illustrates the presence of free water in the fiber structure. Pectin constituting the fiber appears at the vibration of the carboxyl bond between 1430 and 1420 cm −1 . The carboxyl functions of pectins which make up the middle lamella, are also impacted by absorption of water. Stretching of the C-O bonds at the 1225 cm −1 peak identifies the presence of lignin. The cellulose is shown by stretching the C-O-C bond at wavelength 1023 cm −1 which is the characteristic of the polysaccharide nature of the fiber. This signature is also impacted by the increase in humidity [51][52][53]. The infrared spectrum of the fiber shows the ability of the fiber to absorb moisture due to its chemical composition [54].

Thermal Behavior of Sisal Fibers
The thermal behavior of the sisal fiber is illustrated in Figure 3. The first phase is dehydration of the fiber. The dehydration phase takes place between 50 • C and 150 • C while the dehydration is linked to the chemical structure of the fiber components (cellulose, hemicellulose). The components of the fiber (cellulose, hemicellulose) occur between 150 • C and 200 • C. The weight loss at the dehydration phase is 9.67% of the sample weight. The second phase involves the degradation of hemicellulose in view of its branch-rich amorphous structure which facilitates the degradation of its constituents. The decomposition temperature is between 201 • C and 320 • C [55,56]. The weight loss at the decomposition of hemicellulose is 20.9%. The third phase is related to the degradation of the cellulose. It is due to the cleavage of the glycosidic linkage of the cellulosic component in the fiber [57,58]. The last decomposition phase occurs between 390 • C and 475 • C is attributed to free radicals at oxidative atmosphere. The sisal fiber residue at 600 • C is 10.56%. The temperature range of the thermal degradation of the fiber components (cellulose, hemicellulose and lignin) corresponds to the same range of the decomposition of thermoplastic and thermosetting resins [59]. The sisal fiber is resistant to a temperature of 220 • C with a residual mass of 3%. The sensitivity to temperature of the sisal fiber is an important factor in the processing of composite materials based on this fiber.

Thermal Behavior of Sisal Fibers
The thermal behavior of the sisal fiber is illustrated in Figure 3. The first phase is dehydration of the fiber. The dehydration phase takes place between 50 °C and 150 °C while the dehydration is linked to the chemical structure of the fiber components (cellulose, hemicellulose). The components of the fiber (cellulose, hemicellulose) occur between 150 °C and 200 °C. The weight loss at the dehydration phase is 9.67% of the sample weight. The second phase involves the degradation of hemicellulose in view of its branch-rich amorphous structure which facilitates the degradation of its constituents. The decomposition temperature is between 201 °C and 320 °C [55,56]. The weight loss at the decomposition of hemicellulose is 20.9%. The third phase is related to the degradation of the cellulose. It is due to the cleavage of the glycosidic linkage of the cellulosic component in the fiber [57,58]. The last decomposition phase occurs between 390 °C and 475 °C is attributed to free radicals at oxidative atmosphere. The sisal fiber residue at 600 °C is 10.56%. The temperature range of the thermal degradation of the fiber components (cellulose, hemicellulose and lignin) corresponds to the same range of the decomposition of thermoplastic and thermosetting resins [59]. The sisal fiber is resistant to a temperature of 220 °C with a residual mass of 3%. The sensitivity to temperature of the sisal fiber is an important factor in the processing of composite materials based on this fiber.

Contact Angle of Sisal Fibers
Hydrophobicity of the sisal fiber is determined by measuring the contact angle with distilled water. Hydrophobic material is characterized by a T contact angle exceeding 90°, whereas this angle is less than 90° for a hydrophilic material. The contact angle on the surface of the tested sisal fiber is 53.06°; Table 3 lists the obtained result in comparison with Flax and Jute fibers. In general, the natural fibers in the presence of lignocellulosic constituents promote water absorption by hydroxyl groups. Otherwise, this observed hydrophilicity property of the sisal fiber limits its use as reinforcement with a non-polar hydrophobic matrix due to poor interfacial adhesion [60]. The hydrophilicity of the sisal fiber shows the capacity of the fiber to bind water molecules. This factor is very important for the choice of the right matrix for the structure of the fiber.

Contact Angle of Sisal Fibers
Hydrophobicity of the sisal fiber is determined by measuring the contact angle with distilled water. Hydrophobic material is characterized by a T contact angle exceeding 90 • , whereas this angle is less than 90 • for a hydrophilic material. The contact angle on the surface of the tested sisal fiber is 53.06 • ; Table 3 lists the obtained result in comparison with Flax and Jute fibers. In general, the natural fibers in the presence of lignocellulosic constituents promote water absorption by hydroxyl groups. Otherwise, this observed hydrophilicity property of the sisal fiber limits its use as reinforcement with a non-polar hydrophobic matrix due to poor interfacial adhesion [60]. The hydrophilicity of the sisal fiber shows the capacity of the fiber to bind water molecules. This factor is very important for the choice of the right matrix for the structure of the fiber.

Moisture Absorption
Evolution of moisture absorption in the sisal fiber is illustrated in Figure 4. The hydrophilicity and the hygroscopic property of sisal fiber favor moisture absorption with a rate of 13.6% by mass. Sisal fiber shows more attractive moisture behavior than ramie fiber with a recovery rate of 6% [63]. The moisture absorption by the fiber passes through two phases: the first phase corresponds to the absorption of moisture due to capillary action because of the porous structure of the fiber; this phase is shown by the linear part of the absorption curve ( Figure 4). The second phase is characterized by the concave part of the curve explained by the change in the internal structure of the fiber due to the fiber swelling. The absorption of moisture in the sisal fiber contributes to the insertion of water molecules between the molecules of the fiber constituents, which affects the interactions between the cellulose fibrils and the pectin matrix and thus the modulus [64]. absorption curve (Figure 4). The second phase is characterized by the concave part of the curve explained by the change in the internal structure of the fiber due to the fiber swelling. The absorption of moisture in the sisal fiber contributes to the insertion of water molecules between the molecules of the fiber constituents, which affects the interactions between the cellulose fibrils and the pectin matrix and thus the modulus [64].
The diffusion coefficient defined b they Fick law is determined for the moisture absorption in sisal fiber. The fiber is modelled as a cylindrical shaped material according to the Equation (3). The diffusion coefficient of the sisal fiber is 1.77 × 10 −4 mm 2 /s. The diffusion coefficients of the jute and flax fibers tested by Célino et al. are listed in Table 4 [46]. The difference between the moisture absorption capacity of sisal fiber and jute fiber is due to the higher hydrophobicity of sisal fiber (contact angle 53.06°) compared to the hydrophobicity of jute fiber (contact angle 63.9° [62]). However, the Moisture absorption capacity of natural fibers is a parameter affecting the long-term mechanical performance of composites.    Table 4 [46]. The difference between the moisture absorption capacity of sisal fiber and jute fiber is due to the higher hydrophobicity of sisal fiber (contact angle 53.06 • ) compared to the hydrophobicity of jute fiber (contact angle 63.9 • [62]). However, the Moisture absorption capacity of natural fibers is a parameter affecting the long-term mechanical performance of composites.

Structural Fiber Geometry
The length and diameter of sisal fiber were measured to determine the fiber consistency and fitness. The length of sisal fibers (Table 5) shows a variability that depends on the morphology and the extraction method of the fiber. The average value of sisal fiber was 74.5 ± 30.23 mm. The length of the sisal fiber is small compared to the length of the ramie fiber, which is between 50 and 110 mm [64]. The diameter of the sisal fiber is around 121.6-411.0 µm (Figure 5a). The sisal fiber shows an irregular diameter like other vegetable fibers, unlike the uniform shape of synthetic fibers such as glass fiber. The structure of the fiber (Figure 5b) visualized by the scanning electron microscope shows the rough surface consisting of micro fibrils. The nature of the surface of the sisal fiber is due to the presence of impurities (oil, pectin, wax). The surface of the raw fiber can be an obstacle for the good adhesion between the fiber and the matrix. This problem can be solved by the chemical treatments of the fiber (alkaline) [65].  ble fibers, unlike the uniform shape of synthetic fibers such as glass fiber. The structure of the fiber (Figure 5b) visualized by the scanning electron microscope shows the rough surface consisting of micro fibrils. The nature of the surface of the sisal fiber is due to the presence of impurities (oil, pectin, wax). The surface of the raw fiber can be an obstacle for the good adhesion between the fiber and the matrix. This problem can be solved by the chemical treatments of the fiber (alkaline) [65].

Tensile Properties of Sisal Fibers
The tensile properties of sisal fibers were evaluated using the IFBT Test. The fiber volume fraction of the prepared impregnated samples is about 43.62% as an average with a consideration of the fiber density to 1.45 g/cm 3 . The porosity content in these samples was also assessed and low porosity was obtained with about 4.12% as an average. Porosities can be observed on a micrograph made on the cross section of the sample, Figure 6, and located inter the fiber-bundles (yarns). In addition, porosities of smaller size are observed inside the bundles and at the matric-bundle interface. However, the micrograph shows penetration of the resin inside the fiber bundles. fiber testing of sisal fibers [39,[66][67][68].
The tensile modulus of sisal fiber from different origins ( Table 6) shows that the range of values of the elastic modulus of the sisal fiber with Indian origin is similar to the range of the elastic modulus from Moroccan origin analyzed in the current study. The elastic modulus of the flax fiber has been studied by Lansiaux et al. [30] using IFBT method and back-calculated from the rule of mixtures R3 using twisted flax yarns with a twist level between 0 and 100 tpm, the value obtained is around 30 GPa. Figure 6. Cross sectional images of specimens of the IFBT method. Figure 6. Cross sectional images of specimens of the IFBT method.
The elastic modulus of the sisal fibers is back-calculated from the three modified rules of mixtures (Equations (4)-(6)) in addition to the standard rule of mixtures. The elastic modulus of used epoxy resin is Eav = 3 GPa. Figure 7 shows that a significant variation in the modulus values of E1f and E2f by the different rules of mixtures. The elastic modulus back-calculated by the modified rule of mixtures R4 taking into account the porosity content and the fiber orientation gives the upper value of the modulus E1f and E2f. The back-calculated values of the modulus E1f and E2f from different equations are between 9 and 18 GPa. These results are in accordance with other works studying the elementary single fiber testing of sisal fibers [39,[66][67][68].
served inside the bundles and at the matric-bundle interface. However, the micrograph shows penetration of the resin inside the fiber bundles.
The elastic modulus of the sisal fibers is back-calculated from the three modified rules of mixtures (Equations (4)-(6)) in addition to the standard rule of mixtures. The elastic modulus of used epoxy resin is Eav = 3 GPa. Figure 7 shows that a significant variation in the modulus values of E1f and E2f by the different rules of mixtures. The elastic modulus back-calculated by the modified rule of mixtures R4 taking into account the porosity content and the fiber orientation gives the upper value of the modulus E1f and E2f. The backcalculated values of the modulus E1f and E2f from different equations are between 9 and 18 GPa. These results are in accordance with other works studying the elementary single fiber testing of sisal fibers [39,[66][67][68].
The tensile modulus of sisal fiber from different origins (Table 6) shows that the range of values of the elastic modulus of the sisal fiber with Indian origin is similar to the range of the elastic modulus from Moroccan origin analyzed in the current study. The elastic modulus of the flax fiber has been studied by Lansiaux et al. [30] using IFBT method and back-calculated from the rule of mixtures R3 using twisted flax yarns with a twist level between 0 and 100 tpm, the value obtained is around 30 GPa.  The tensile modulus of sisal fiber from different origins ( Table 6) shows that the range of values of the elastic modulus of the sisal fiber with Indian origin is similar to the range of the elastic modulus from Moroccan origin analyzed in the current study. The elastic modulus of the flax fiber has been studied by Lansiaux et al. [30] using IFBT method and back-calculated from the rule of mixtures R3 using twisted flax yarns with a twist level between 0 and 100 tpm, the value obtained is around 30 GPa.

Diameter and Linear Density of Sisal Yarns
The sisal yarn is made of a bundle of z-twisted fibers with a twist level about 80 tpm, Figure 8. Mostly, the yarns are z-twisted because the spinning machines are made for righthanded manipulators [69]. The average value of the sisal yarns' diameter is 2.01 ± 0.9 mm ( Table 7). The average value of the linear density of the sisal yarns is around 3.3 ± 0.7 kTex. The variability of the linear density along the yarn is due to the irregular value of twisted sisal fibers diameter. This variability is due to the variability of the climatic conditions and the cultivation conditions of the plant. The sisal yarns with Algerian origin were studied by Belaadi et al. [70] and the obtained results show that the average diameter value is 2.04 ± 0.26 mm.

Diameter and Linear Density of Sisal Yarns
The sisal yarn is made of a bundle of z-twisted fibers with a twist level about 80 tpm, Figure 8. Mostly, the yarns are z-twisted because the spinning machines are made for right-handed manipulators [69]. The average value of the sisal yarns' diameter is 2.01 ± 0.9 mm ( Table 7). The average value of the linear density of the sisal yarns is around 3.3 ± 0.7 kTex. The variability of the linear density along the yarn is due to the irregular value of twisted sisal fibers diameter. This variability is due to the variability of the climatic conditions and the cultivation conditions of the plant. The sisal yarns with Algerian origin were studied by Belaadi et al. [70] and the obtained results show that the average diameter value is 2.04 ± 0.26 mm.

Yarn
Average Diameter (mm)

Average Linear Density (Ktex) Twist Level (Tpm)
Sisal Yarn 2.01 ± 0.9 3.3 ± 0.7 80 The twist angle of the sisal yarns was calculated theoretically based on the fiber orientation factor which can be expressed as a function of the twist angle on the outer surface of the yarn (Equation (6)) [49,50].
is evaluated by two methods with the assumption of circular cross-section shape of the twisted yarn:  The twist angle of the sisal yarns was calculated theoretically based on the fiber orientation factor which can be expressed as a function of the twist angle α on the outer surface of the yarn (Equation (6)) [49,50].
α is evaluated by two methods with the assumption of circular cross-section shape of the twisted yarn: Method 1: the section of the yarn is not considered to be full of fibers (Equation (5)) (see Section 2.2.6.). The fibers are packed in the yarn section can be expected to behave according to the semi-empirical model proposed by Pan [49,50], and given by Equation (5).
Thus, the fiber twist angle on the outer surface of the yarn is computed by Equation (7). α = tan −1 T 10 −3 4π kTex ρφ (7) Method 2: the section of the yarn is considered full of fibers so the fiber packing factor in the yarn section is equal to one and the fiber twist angle is computed by Equation (8): where: • ρ is the density of Moroccan sisal fibers (1.45 kg/m 3 ). • T is the twist level in tpm (turns per meters).
The twist angle of the sisal yarns was calculated based on the two configurations of the yarn cross-section. The twist angle of the yarns obtained with a packing factor is 26.92 • , while the twist angle without a packing factor is 23.16 • . From the elaborated samples of the IFBT method (Section 2.2.6.), the packing factor was verified by an Equation by measuring the yarn cross-sectional area in Figure 7.
The Yarn Packing Factor: where: • ρ is the density of Moroccan sisal fibers (1.45 kg/m 3 ).

•
KTex is the linear density of yarns (Ktex). • A fil is the yarn cross-sectional area (m 2 ).
The yarn packing factor is 70.2% with a standard deviation of 5.63%. With this packing factor verification, the yarn twist angle is 23.65 ± 3.65 • . In order to evaluate the accuracy of the different methods, the twist angle was measured from the yarn pictures using the ImageJ software (Figure 9), the value obtained is 23.14 • with a standard deviation of 2.75 • . The obtained results (Table 8) show that the twist angle calculated by both methods is included in the range of variation of the measured twist angle. Method 1: the section of the yarn is not considered to be full of fibers (Equation (5)) (see Section 2.2.6.). The fibers are packed in the yarn section can be expected to behave according to the semi-empirical model proposed by Pan [49,50], and given by Equation (5).
Thus, the fiber twist angle on the outer surface of the yarn is computed by Equation (7).

= (7)
Method 2: the section of the yarn is considered full of fibers so the fiber packing factor in the yarn section is equal to one and the fiber twist angle is computed by Equation (8): where: • ρ is the density of Moroccan sisal fibers (1.45 kg/m 3 ).

•
T is the twist level in tpm (turns per meters).
The twist angle of the sisal yarns was calculated based on the two configurations of the yarn cross-section. The twist angle of the yarns obtained with a packing factor is 26.92°, while the twist angle without a packing factor is 23.16°. From the elaborated samples of the IFBT method (Section 2.2.6.), the packing factor was verified by an Equation by measuring the yarn cross-sectional area in Figure 7.
The Yarn Packing Factor: where: • ρ is the density of Moroccan sisal fibers (1.45 kg/m 3 ).
The yarn packing factor is 70.2% with a standard deviation of 5.63%. With this packing factor verification, the yarn twist angle is 23.65 ± 3.65°. In order to evaluate the accuracy of the different methods, the twist angle was measured from the yarn pictures using the ImageJ software (Figure 9), the value obtained is 23.14° with a standard deviation of 2.75°. The obtained results (Table 8) show that the twist angle calculated by both methods is included in the range of variation of the measured twist angle.

Tensile Behavior of Sisal Yarns
The tensile properties of sisal yarns were studied under quasi-static loading. Figure 10 shows load versus strain of sisal yarns. The curves can be divided into two zones, the first zone is non-linear between 1 and 4% of the deformation value, which can be explained by the realignment of the twisted fibers during the tensile test, which creates a space between the fibers in the yarn structure [71]. The second zone is quasi-linear between 4 and 7% of the deformation before yarn breakage. The observed behavior of the sisal yarns of Moroccan origin studied in this work is similar to the behavior of sisal yarns of Algerian origin [70] and flax yarns [71].

Tensile Behavior of Sisal Yarns
The tensile properties of sisal yarns were studied under quasi-static loading. Figure  10 shows load versus strain of sisal yarns. The curves can be divided into two zones, the first zone is non-linear between 1 and 4% of the deformation value, which can be explained by the realignment of the twisted fibers during the tensile test, which creates a space between the fibers in the yarn structure [71]. The second zone is quasi-linear between 4 and 7% of the deformation before yarn breakage. The observed behavior of the sisal yarns of Moroccan origin studied in this work is similar to the behavior of sisal yarns of Algerian origin [70] and flax yarns [71].
The average value of the maximum load is 672.36 ± 236 N. The dispersion of the results of this test is due to the variability of the properties of the sisal fibers. The average value of the tenacity of sisal yarns is 20.24 cN/Tex (±3.25) ( Table 9). The tenacity of sisal yarns of Indian origin is of the order of 18.38 cN/Tex with a maximum strength of 514.7 N [72]. The difference in the mechanical properties of the sisal yarns according to the origin depends on the climatic and growing conditions of the agave sisalana plant in different regions.

Conclusions
This work aims to analyze the physical and mechanical properties of sisal fibers and yarns of Moroccan origin. The chemical composition of the fiber contains cellulosic ele- The average value of the maximum load is 672.36 ± 236 N. The dispersion of the results of this test is due to the variability of the properties of the sisal fibers. The average value of the tenacity of sisal yarns is 20.24 cN/Tex (±3.25) ( Table 9). The tenacity of sisal yarns of Indian origin is of the order of 18.38 cN/Tex with a maximum strength of 514.7 N [72]. The difference in the mechanical properties of the sisal yarns according to the origin depends on the climatic and growing conditions of the agave sisalana plant in different regions.

Conclusions
This work aims to analyze the physical and mechanical properties of sisal fibers and yarns of Moroccan origin. The chemical composition of the fiber contains cellulosic elements that promote a mechanical performance of the reinforcement of the composite material. The thermogravimetric analysis of sisal fiber shows a low thermal stability. The presence of the hydroxyl groups (OH) in the structure of the fiber favors the hydrophilic behavior of the fiber, with a contact angle of 53.06 • . The moisture absorption rate of the fiber is 13.2%. The mechanical properties of the fiber analyzed by the IFBT method shows that the porosity of the impregnated yarns and the twist angle of the yarns influence the elastic modulus of sisal fiber. The physical and mechanical properties of sisal yarns show a variability in the values of diameter and a linear density that depends on the properties of the twisted fibers used in the spinning process. This investigation on the properties of sisal fiber of Moroccan origin proves that this fiber can be an alternative fiber to synthetic fibers because of its recyclability, low density, eco-friendly and its abundance in several regions, which are highly demanded points in the strategy of sustainable development. The limitations of the fiber such as moisture absorption and low thermal stability can be optimized for the development of reinforcements of composite materials.