Recent decades, textile fabrics used as reinforcements for manufacturing composites are spreading in various fields, such as aerospace, automotive, construction and basic facilities industries     . Woven fabrics, characterized with the special structure form, are widely used as reinforcement of textile composites in structural applications by far   . Generally woven fabrics consist of two sets of yarn components, known as warp yarns and weft yarns  . The overlap structure in the cross section provides with relatively good dimensional stability  .
Roving cloth, as one of the woven fabrics, is characterized with the cross section of the warp fiber bundles and weft fiber bundles, which are usually used as reinforcement for laminate molding   . Roving cloth has the advantage that it can be used for relatively complex shapes or curvatures. However, due to the crimp of the warp and weft fiber bundles, there will be a tendency for the load direction fiber bundles to stretch   , forming a compressive stress against the transverse fiber bundles within the cross-section area     . Generally, final fracture of composites is considered to be caused by the accumulation of these microcracks or fractures   . Thus, initial fracture, which represents the initiation of micro fracture, is considered to be significant in determining the fracture mechanism. GFRP (glass fiber reinforced plastic) materials, characterized by high strength, low cost, low density, corrosion resistance, are also adopted in the form of glass woven fabric as reinforcements. Two kinds of glass roving clothes were used in this study, and effort was made to find out the fracture mechanism, especially the initial fracture behavior.
In previous works, many researchers focused on the mechanical properties as well as the fracture mechanism of glass woven cloth reinforced composites. Demircan, O. carried out the investigation of the fracture process and mechanisms of glass roving and glass cloth composites  . García, I. et al. tested the cross-ply glass-reinforced polyester composites in a direction parallel to one of the directions of reinforcement. The results showed that the crack spacing measurements are against applied stress for specimens with differing ply thicknesses  . Boccardi, S. et al. conducted a research that focused on the behavior of inter-laminar graded interface strength (IGIS) under impact. Results highlight the role played by the stacking sequence in the IGIS laminate  . Jaafer and Muslem  conducted tests on DCB woven roving and CSM interface layers and observed fibre bridging but did not calculate the associated strain energy release rates. Ono, K., Fujii, Y., and Wada, A. investigated the possibility of non-destructive examination with ultrasonic wave testing for mechanical damage of glass fiber reinforced plastics  .
Because of the complicated microstructure of woven fabric, it is of significance to understand of mechanical mechanism of textile fabric reinforced composite materials. There are references   showing that initial fracture occurred in woven fabric reinforced composites has been observed in transverse fiber bundles, which has a significant effect on the mechanical properties of composites. However it is difficult to observe the cracks during the test neither to measure accurately, as a result, simulations carried out by modeling or FEM analysis are often adopted as an effective way to get a knowledge of the initial fracture process  . Researchers try to understand the woven fabric reinforced composites by numerical analysis method     . Zako, M., and Uetsuji, Y.  investigated the damage behavior of woven fabric reinforced FRP, and simulated by finite element analysis using an anisotropic damage model based on damage mechanics. Xu, Z. et al. chose CF/GF hybrid woven fabric as the reinforcement and the initial fracture properties was discussed. Then models of the woven fabric was built by and initial behavior of fracture progress in resin hybrid composite materials were predicted   . Tabiei A. and Jiang Y.  have tried to develop a simple model considering the two-dimensional extent of woven fabric and to have an interface with nonlinear finite element codes. Ichihashi and Hamada et al.  proposed an FEA model for evaluating mechanical properties of single ply of woven fabric structure, focusing on the initial fracture behavior of the woven fabric reinforced composite.
In this paper, investigation on the initial fracture behavior was carried out on roving glass woven fabric reinforced composites manufactured by the hand lay-up method. Two kinds of roving glass woven fabrics of different FAW (Fabric Area Weight) were adopted in this study. Tensile tests combined with AE (acoustic emission) measurements   were conducted and the initial fracture behavior was observed. These two kinds of composites molded by fabrics with different crimp ratio showed difference in tensile properties as well as different fracture mechanism. The results indicated that the crimp ratio decided by the yarn density and fiber bundle tex, is playing an important role in the fracture mechanism.
In the off-axis cases of 5˚/10˚/80˚/85˚, the fracture mechanism and the shear stress were compared and discussed  . Models of the woven fabric were built according to the geometric dimensions of the fabric, and the fracture mechanism was simulated using the Marc software. By combining the experimental results and finite element analyses, the initial behavior of fracture progression was discussed and predicted.
2. Materials and Experimental Methods
2.1. Materials and Specimens
In this paper, two kinds of roving glass woven fabrics of different FAW (Fabric Area Weight), 570 g/m2 and 800 g/m2, were adopted. Unsaturated polyester resin (150 HRBQTNA, Showa Denko K.K.) was adopted as the matrix. The glass roving clothes manufactured by Hokuriku Fiberglass Co., Ltd are in plain structure as showing in Figure 1. The one with lighter FAW is named as Type A (570 g/m2), and the heavier one is named as Type B (800 g/m2). Detailed information of the glass fiber bundle (Nitto Boseki Co., Ltd.), such as linear density and yarn density, is listed in Table 1. It can be seen that fiber bundles of same linear density were used while the ends per inch for the warp direction and the picks per inch in the weft direction are different. And Type A is of relatively higher yarn density. In this study, the warp direction is defined as 0 degree and the weft direction as 90 degree.
Table 1. Information of the two types of roving glass cloth.
Figure 1. Photographs of the two types of roving glass cloth.
The glass roving cloth reinforced composites were manufactured by hand lay-up molding method for only 1 ply for the purpose of reveal the crack propagation and fracture mechanism. After the fabrication, the composite board was cut to the size of 200 mm × 20 mm (length × width) according to ASTM D3039 according to 6 orientations 0˚/5˚/10˚/80˚/85˚/90˚. The thickness of Type A is approximately 0.60 mm, while Type B is approximately 0.65 mm, a little thicker than Type A. Cross-sectional observation was carried out and the crimp percentage for both the 0 and 90 degree directions of Type A and Type B was calculated, and is showed in Table 1. Crimps are formed because of the weaving process. It can be seen from Figure 1 and Table 1 that Type A is of relatively higher yarn density, and the crimp ratio of Type A is slightly higher than that of Type B. It means for the same unit length, the fiber bundle length of Type A is longer than that of Type B.
2.2. Experimental Methods
Mechanical investigation of the tensile properties was carried out on specimens in 6 degrees 0˚/5˚/10˚/80˚/85˚/90˚. The tensile tests were carried out on an Instron universal testing machine at a speed of 1 mm/min and the test room temperature was 22˚C according to the ASTM D3039. (Specimen number N = 3). For 0˚ and 90˚, which are also defined as on-axis cases, the fracture progress during the tensile test has been observed by take video. During the tensile tests, an AE (acoustic emission) device was used in order to detect when the initial fracture occurred. A video was also shot in order to understand the fracture process. In the photographs collected from the video at different strain stages, three main periods can be identified to depict the fracture process. Combined with the AE data and the video footage, the initial fracture behavior has been discussed. For 5˚/10˚/80˚/85˚, which are defined as off-axis cases, tensile modulus, tensile strength were also tested and discussed. Besides, the difference of deferent degrees, and the comparison between Type A and B has been discussed and summarized.
3. Experimental Results and Discussion
3.1. Tensile Results of On-Axis Cases
The stress-stain curves are illustrated in Figure 2, and details of the results of the tensile tests are shown in Figure 3, and Table 2. It can be known that the tensile modulus of Type B is similar with Type A, and the tensile modulus of the 90 degree direction is higher than that of the 0 degree direction.
On the other hand, the tensile strength of Type B turned out to be much higher than that of Type A for both the 0 and 90 degree directions. The fiber bundle cross section area ratio of Type A and Type B was calculated and this value could be referred as an index similarly with Vf (volume fraction of reinforcing fibers). The result showed the value of Type A is slightly higher than that
Figure 2. Stress-strain curves of two types of GFRP.
Figure 3. Tensile modulus and tensile strength of two types of GFRP.
Table 2. Tensile results of two types of GFRP.
of Type B, which is revealing that Type A might get higher tensile strength. However, it does not agree with the tensile strength result.
When it comes to the elongation, it can be known from Table 2 that the elongation rate for Type B is higher than Type A. This would be taken as one of the reasons for the high tensile strength of Type B. But the reason why Type B achieved higher elongation should be found out firstly because the crimp ratio of Type B is lower than that of Type A which indicated Type A should have got higher elongation.
Otherwise, the initial fracture results detected by acoustic emission showing in Table 2 revealed that for both Type A and Type B in both 0 degree and 90 degree directions, the initial fracture stress values are similar.
3.2. Initial Fracture Process and Observation
In order to decide where exactly the initial fracture happened, specimens were under tensile testing and stretched to 1.3 times the initial fracture stress and stopped. Then observation was carried out on the longitudinal cross section by optical microscopy. From the photographs shown in Figure 4, it can be clearly seen that transverse cracks first occurred in the transverse fiber bundles, which can be regarded as contacted to the initial fracture stress.
The fracture processes are illustrated by still images from the video in Figure 5. It can be obviously observed that transverse cracks appeared within the transverse fiber bundles at the very early stage. And with the transverse crack increased, another type of cracks in wave shape began to appear between transverse fiber bundles. The schemas of the transverse crack and the wave shape crack are illustrated in Figure 6(a) and Figure 6(b). The wave shape cracks formed because of the existence of the crimp. With the loading going higher, the fiber bundles along the load direction were stretched and had a tendency to become straight, which gave a compressive stress to the transverse fiber bundles. The resin rich area between the transverse fiber bundles was compressed and wave shape cracks formed along the iso-stress surface, which could be the direct reason leading to the final fiber bundle breaking. The wave shape cracks showing in Figure 6(b) are characteristic cracks appearing in woven fabric reinforced composited because of the woven structure.
There were several period before the final fracture can be summarized during the fracture period for both Type A and Type B involved with transverse cracks
Figure 4. Cracks observed in transverse fiber bundles as initial fracture.
Figure 5. Fracture process during tensile test. (a) Type A; (b) Type B.
1) Firstly, for both Type A and Type B, transverse cracks occurred firstly in the transverse fiber bundles, and then the number of transverse cracks increases.
Figure 6. Schemas of two types of crack modes. (a) Transverse cracks; (b) Wave shape cracks.
Figure 7. Schemas of fracture crack propagation. (a) Type A; (b) Type B.
During this period, the interface between the glass fibers and the resin is subjected to most of the load. White line of cracks can be observed in the cross section area of the warp and weft fiber bundles.
2) During this period, the interface between the glass fibers and the resin is still the main carrier of most of the load. White cracks increased gradually until the cracks went through the transverse fiber bundles. It can be seen that transverse cracks within Type B are much finer than those in Type A.
3) Because of the tendency of the fiber bundle in the crimp to be stretched, the cross section of the warp fiber bundle and weft fiber bundle is compressed and wave shape cracks began to show up both for Type A and Type B. During this period, the interface between the glass fiber and the resin still carried the most of the load but the longitudinal fiber bundles are beginning to take most of the load. More cracks in wave shapes appeared between the transvers fiber bundles were observed in Type A and then increased from Figure 7(a). While for Type B, very few wave shape cracks were also observed, and more fine transverse cracks increased.
4) Fourthly, the longitudinal fiber bundles become carrying the most of the load and finally leading to the fracture of the specimen. It was obviously observed that much more wave shape cracks in Type A and propagated in transverse directions. However for Type B, wave shape cracks did not increase much, and plenty of fine transverse cracks were observed within the transverse fiber bundles.
When comparing Type A with Type B, the difference was observed from the second period. In the case of Type B, more and finer cracks were observed, while for Type A, more wave shape cracks propagated during the third and fourth period. AE results showed in Figure 8. It can been seen that the accumulated AE conus of Type A increased faster than that of Type B and kept higher above Type B. Comparing observation results with AE results, it can be considered that the wave shape cracks resulted in the main fracture and finally leading to the fiber bundle breaking along the load direction. Even the crimp ratio of Type A is
Figure 8. Results of the accumulated AE couns for Type A and Type B.
higher than that of Type B, the wave shape cracks led to fiber bundle breaking and specimen of Type A came to final fracture before it achieved a high elongation. On the contrast, the higher fiber bundle tex and lower crimp ratio of Type B makes it more tolerant with transverse cracks, and it is considered the specimen was stretched smoothly finally with a higher elongation than Type A.
Higher crimp ratio might lead a tendency to fierce wave shape cracks and result in a low tensile strength. But it is also predictable that the fabric reinforcement with a very low crimp ratio will also not achieve good tensile properties because of the low fiber volume fraction. It is considered that for woven fabric reinforced composites, there is a critical crimp ratio which contributes to better tensile properties.
3.3. Tensile Results of Off-Axis Cases
Tensile test was also carried out on specimens of 5˚/10˚/80˚/85˚and the tensile results were summarized in Table 3.
The summary of tensile modulus and strength of all 6 degrees of 0˚/5˚/10˚/80˚/85˚/90˚ are illustrated in Figure 9 and Figure 10 separately. Comparing to Figure 7 and Figure 8, difference between on-axis and off axis
Figure 9. Tensile modulus summary of all 6 degrees.
Figure 10. Tensile strength summary of all 6 degrees.
Table 3. Tensile results of off-axis cases for two kinds of GFRP.
can be seen obviously, that specimens of off-axis have much lower tensile modulus and strength. Even in the case of 5 degree and 85 degree, tensile modulus and tensile strength decreased dramatically comparing with 0 degree and 90 degree, which can be considered as a result of the shear stress.
3.4. Fracture Mechanism Observation
The photographs of fractured specimens of Type A and Type B after the tensile test were shown in Figure 11 and Figure 12 separately. It can be seen from the photograph that the fracture of 5 degree and 85 degree seemed to be fiercer than that of 10 degree and 80 degree. The fracture schema is illustrated in Figure 13. For 5 degree and 85 degree, cracks propagated along fiber bundles and ultimate fiber fracture happened with the fiber bundle split up. It seems no big difference between Type A and Type B. While it also can be observed from the schema that, in the case of 10 degree and 80 degree, especially for Type B, it can be observed that less fiber was pulled out and split, which revealed that the fiber bundle in the loading direction did not take most of the loading. And the shear stress played an important role during the tensile test, leading to the ultimate fracture of the specimen, which can be considered as the reason of the low strength in the case of the 10degree and 80 degree. Especially for Type B, there are very few fiber bundles pulled out and split in the case of the 10degree and 80 degree, and contacted with the tensile strength result which showing a decrease of more than 45% comparing with 5 degree and 85 degree.
4. Modeling and Simulation
4.1. FEM Model Building
To build up the model in software Marc, precise geometric parameters are necessary. Similar work has been done previously   . The specimens before test were cut and buried in resin. After polishing, optical observation was carried out to obtain the geometric parameters. Based on the results of the geometric and mechanical property parameters shown in Table 4 and Table 5, models of both the 0 degree and 90 degree directions were built up first and the interfacial strength was calculated. Models of the 0 degree direction for both Type A and Type B are illustrated in Figure 14 as examples.
Figure 11. Specimens after fracture for off-axis case of Type A.
Figure 12. Specimens after fracture for off-axis case of Type B.
Figure 13. Ultimate fracture schema of off-axis cases. (a) Type A; (b) Type B.
Table 4. Geometric parameters of two types of GFRP.
Table 5. Mechanical parameters for FEM models.
Figure 14. FEM models built by MARC in 0 degree direction. (a) Type A; (b) Type B.
After the model was built, the mechanical properties of the weft and warp yarns were identified according to Table 4. For the mechanical parameters, Vf, Young’s modulus, and the Poisson ratio of a glass fiber bundle and the unsaturated polyester resin are required. Since the elastic modulus of a fiber bundle in the matrix is different from its catalogue elastic modulus value, it needs to be recalculated in order to obtain a precise value in practical case. Here, the Vf of the fiber bundle is calculated by the number and the area of the filaments and the cross-sectional area of the fiber bundle by Equation (1) as follows:
where A is the area of the fiber bundle cross-section; N is the number of filaments in one fibber bundle; and S is the cross-sectional area of the fiber bundle.
4.2. Tensile Modulus Calculation
Before obtaining the interfacial strength, the tensile modulus was first calculated in order to confirm the validity of woven fabric composite model, and the interfacial strength value as well. The simulation was carried out within the elastic tensile area. Unit displacement of 1mm was applied to the models in 0 degree or 90 degree. Deformation in the loading direction was observed. To calculate the simulation tensile modulus, the total stress of all the nodes along the applied force direction was summed up first, and then dividing by the corresponding cross-sectional area.
4.3. Interfacial Strength Calculation
As stated above, initial fracture occurred within the transverse fiber bundle of the specimen with the load increasing during the experiment. And the accumulation of these transverse cracks leads to the propagation into resin. In order to obtain the interfacial strength in the transverse fiber bundle, the principal stress of the interfacial elements was calculated by dividing the axial force of the interface element by the cross-sectional area. Combining this with the initial fracture stress from the experiment and the initial fracture stress from the model, the strength of the glass fiber and carbon fiber bundles in the transverse direction was calculated by Equation (2).
where σinterfacial is the simulating value of interfacial strength in a transverse fiber bundle from modeling; σinitial is experimental value of the initial fracture stress from experiments; is the simulating value of initial fracture stress obtained from model; is the simulating value of interfacial strength obtained from on axis models.
4.4. Simulation Results
The analysis results of the modulus by FEM are compared with the results from the experiments as shown in Table 6. It can be noted that the calculated tensile modulus deviation is below 5%, which confirms the validity of the calculation of the interfacial strength. The calculation results of the interfacial strength are also shown in Table 6. For both the 0 degree and 90 degree directions for Type A and Type B, the values of interfacial strength are similar, both approximately 86
Table 6. Simulation results of two types of GFRP.
MPa. In the case of Type B, the interfacial value is a little higher than Type A with a deviation rate below 5%. In addition, transverse elements with the maximum value obtain from the model can be regarded as the location where initial fracture occurred. As shown in Figure 15, the initial fracture occurred within the cross section of the warp and weft fiber bundles, which is in good agreement with the transverse crack from optical results.
As it is comes to transverse cracks and wave shape cracks, the difference observed between these two crack modes was considered possible to be simulated by FE models. However the wave shape cracks happened beyond the elastic area. In order to simulate these two crack modes, 3D models are planned to be built in future work, and further boundary conditions will be applied.
In this paper, tensile properties, especially the initial fracture properties of two types of roving glass cloth reinforced composites were tested. Tensile test was conducted and tensile properties were discussed on 6 degrees 0˚/5˚/10˚/80˚/85˚/90˚. The tensile results of 0 degree and 90 degree showed that the Type B has higher tensile modulus and strength than Type A. As it is referred to the initial fracture behavior 4 stages have been observed, according to the main load carrier changed from the interface of glass fiber and the resin to the longitudinal fiber bundles. More cracks in wave shape were observed as the crack propagating, while more fine cracks in the transverse fiber bundles were observed in the case of Type B. The important role of crimp ratio was made clear, which can be considered of significance to obtain better tensile property. Among the cases of 5˚/10˚/80˚/85˚, dramatic decrease in tensile modulus and strength was made clear. And it seem that less fiber bundles pulled out from the specimen during the fracture in the cases of 10˚ and 80˚, revealing that less longitudinal fiber bundles carried the load.
According to the geometric parameters from the specimens, a FEM model of the woven structure was built, and the tensile modulus and the interfacial properties of the transverse fiber bundle were simulated and compared. The results show that the interfacial strength values of Type B were slightly higher than Type A, and similarly with a deviation rate below 5%. And although the FAW values are different of Type A and Type B, the simulation calculating results
Figure 15. Initial fracture simulation results.
showed that interfacial strength in the glass transverse fiber bundles is almost the same with a low deviation. Where the initial fracture happens was also predicted and proved agreed with the experiment results well. In future work, 3D models in multiscale are supposed to be built and further simulation of two different crack modes are also planned to be carried out.
 Xu, F., Fan, W., Zhang, Y., Jia, Z., Qiu, Y. and Hui, D. (2016) Modification of Tensile, Wear and Interfacial Properties of Kevlar Fibers under Cryogenic Treatment. Composites Part B: Engineering, 116, 398-405.
 Xie, J., Yao, L., Xu, F., Li, Y., Shan, Z., Hui, D., et al. (2014) Fabrication and Characterization of Three-Dimensional PMR Polyimide Composites Reinforced with Woven Basalt Fabric. Composites Part B Engineering, 66, 268-275.
 Akmar, A.I., Lahmer, T., Bordas, S.P.A., Beex, L.A.A. and Rabczuk, T. (2014) Uncertainty Quantification of Dry Woven Fabrics: A Sensitivity Analysis on Material Properties. Composite Structures, 116, 1-17.
 Sockalingam, S., Chowdhury, S.C., Gillespie Jr, J.W. and Keefe, M. (2017) Recent Advances in Modeling and Experiments of Kevlar Ballistic Fibrils, Fibers, Yarns and Flexible Woven Textile Fabrics—A Review. Textile Research Journal, 87, 984-1010.
 Mostafa, N.H., Ismarrubie, Z.N., Sapuan, S.M. and Sultan, M.T.H. (2016) Effect of fabric Biaxial Prestress on the Fatigue of Woven E-Glass/Polyester Composites. Materials & Design, 92, 579-589.
 Li, L., Lomov, S.V., Yan, X. and Carvelli, V. (2014) Cluster Analysis of acoustic Emission Signals for 2D and 3D Woven Glass/Epoxy Composites. Composite Structures, 116, 286-299.
 Yu, B., Blanc, R., Soutis, C. and Withers, P.J. (2016) Evolution of Damage during the Fatigue of 3D Woven Glass-Fibre Reinforced Composites Subjected to Tension-Tension Loading Observed by Time-Lapse X-Ray Tomography. Composites Part A: Applied Science and Manufacturing, 82, 279-290.
 Nilakantan, G. and Gillespie, J.W. (2013) Yarn Pull-Out Behavior of Plain Woven Kevlar Fabrics: Effect of Yarn Sizing, Pullout Rate, and Fabric Pre-Tension. Composite Structures, 101, 215-224.
 Penava, Z., Simic-Penava, D. and Knezic, Z. (2014) Determination of the Elastic Constants of Plain Woven Fabrics by a Tensile Test in Various Directions. Fibres & Textiles in Eastern Europe, 22, 57-63.
 García, I.G., Mantic, V., Blázquez, A. and París, F. (2014) Transverse Crack Onset and Growth in Cross-Ply [0/90] s Laminates under Tension. Application of a Coupled Stress and Energy Criterion. International Journal of Solids and Structures, 51, 3844-3856.
 Boccardi, S., Meola, C., Carlomagno, G.M., Sorrentino, L., Simeoli, G. and Russo, P. (2016) Effects of Interface Strength Gradation on Impact Damage Mechanisms in Polypropylene/Woven Glass Fabric Composites. Composites Part B: Engineering, 90, 179-187.
 Almuhammadi, K., Alfano, M., Yang, Y. and Lubineau, G. (2014) Analysis of Interlaminar Fracture Toughness and Damage Mechanisms in Composite Laminates Reinforced with Sprayed Multi-Walled Carbon Nanotubes. Materials & Design, 53, 921-927.
 Ono, K., Fujii, Y. and Wada, A. (2015) Investigation of Non-Destructive Examination for Mechanical Damage of FRP. In: ASME 2015 International Mechanical Engineering Congress and Exposition, American Society of Mechanical Engineers, New York, V009T12A029.
 Dixit, A. and Mali, H.S. (2013) Modeling Techniques for Predicting the Mechanical Properties of Woven-Fabric Textile Composites: A Review. Mechanics of Composite Materials, 49, 1-20.
 Zako, M., Uetsuji, Y. and Kurashiki, T. (2003) Finite Element Analysis of Damaged Woven Fabric Composite Materials. Composites Science and Technology, 63, 507-516.
 Xu, Z., Ichikawa, D. and Yang, Y. (2015) Initial Fracture Behavior of Intra CF/GF Woven Fabric Composites. In: ASME 2015 International Mechanical Engineering Congress and Exposition, American Society of Mechanical Engineers, New York, NV02BT02A013.
 Xu, Z., Nakai, A., Yang, Y. and Hiroyuki, H. (2017) A Study on the Initial Fracture Behavior of CF/GF Intra-Hybrid Woven Fabric Reinforced Composites. Open Journal of Composite Materials, 8, 11.
 Bussiba, A., Kupiec, M., Ifergane, S., Piat, R. and Bohlke, T. (2008) Damage Evolution and Fracture Events Sequence in Various Composites by Acoustic Emission Technique. Composites Science and Technology, 68, 1144-1155.
 HaghiKashani, M., Hosseini, A., Sassani, F., Ko, F.K. and Milani, A.S. (2017) The Role of Intra-Yarn Shear in Integrated Multi-Scale Deformation Analyses of Woven Fabrics: A Critical Review. Critical Reviews in Solid State and Materials Sciences, 43, 1-20.