Measures describing curvilinear short fiber distributions

In this article, we discuss measures for fibers having a curvilinear shape. This is the case, for example, for man-made cellulose fibers having a weak stiffness. The fibers are bent during the injection molding process of short fiber reinforced plastics. For this purpose, μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu $$\end{document}-CT data can be evaluated and several measures can be introduced defining the geometrical orientation of the fibers. These measures are the length, a mean curvature, and the mean torsion. Furthermore, a mean orientation of a fiber and a mean deviation to a straight line can be defined. Additionally, to these measures, which are based on a continuous interpolation of given data points, discretized quantities only considering the data points are compared. Finally, the distributions of these measures at real μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu $$\end{document}-CT data are provided.

(a) Fig. 1 Geometrical description of curves can be provided. There are two possibilities: first, the simplest approach is to draw on polygon approximations using the coordinates describing the fiber, which is easy to be implemented. Second, continuous formulations can be applied implying further measures such as curvature, and torsion. The latter terms will be explained in detail in Sect. 2. We draw on a spline concept to interpolate the coordinates so that each fiber obtains a continuous representation, and, accordingly, the possibility to apply concepts in differential geometry.
The article is structured as follows: since we are interested to introduce measures of curvilinear fibers, a summary of the description of the local behavior of curves is provided. In a second step, global mean values are proposed. These mean values are the orientation, curvature, torsion, and a mean deviation from a straight line. Afterwards, one special interpolation method, which is based on classical, natural splines, is drawn on to evaluate the proposed measures. Afterwards, a discretized, simplified version is offered. Both schemes are compared using real μ-CT data of a man-made cellulose fiber reinforced polypropylene specimen. In this article, the distributions of the measures are discussed.

Measures of curvilinear curves
First, we recap some basics in the description of curves, and, subsequently, some measures of curvilinear fibers are proposed characterizing their orientation, curvature and torsion.

Properties of curves
A curve in the three-dimensional case is described by the position vector r(ξ ), r ∈ E 3 , which depends on a parameter ξ , ξ ∈ [ξ 1 , ξ 2 ]. The tangent vector is given by see, for example, [3,11]. The arc length represents the length of the covered path from point ξ 1 until point ξ on the curve, whereas denotes the length of the curve between two points ξ 1 and ξ 2 (ξ 2 > ξ 1 ), see Fig. 1a.
defines the rate of the arc length representing the norm of the tangent vector or the square root of the metric coefficient concerned. Since s, ξ = 0 holds, there exists a functionξ(s) having the derivativê In this sense, we will need the mapping between the parameter ξ and the arc length s. In this case, the curve can also be related to the arc lengthr (s(ξ )) = r(ξ ), (6) so that the tangent vector (1) reads g ξ = r, ξ =r, s s, ξ =r, s g ξ (7) leading to the obvious property e t is called the unit tangent vector. Using the property we can see that the vector e t , s =r, ss is orthogonal to e t , see Eq. (8). This property is chosen to define the unit normal vector The norm ofr, ss represents the curvature of the curve κ(s) := r, ss = e t , s .
The two vectors e t and e n are chosen to define the binormal unit vector see Fig. 1a, i.e. the natural basis or trihedron of the curve. Using the Frenet equations, see [11], e t , s = κe n , e n , s = −κe n + τ e b , e b , s = −τ e n (13) motivates an additional local measure-apart the curvature-namely the torsion τ . Eq. (13) 1 is already provided by Eq. (10) using definition (11), i.e. κ represents the (local) deviation of the curve from a straight line. τ is a measure for the (local) deviation of a curve from a flat course, which can be seen at Eq. (13) 2,3 , where the curve remains in the (e t , e n )-plane. The torsion itself can be calculated by see Eq. (13) 2,3 . Its sign represents the rotation (right-or left handed). All these expressions look quite simple if an arc length representation is given. However, later on only the parametric representation r(ξ ) is provided, which requires a different representation. Here, we follow the presentation in [2]. For the subsequent calculations, d dξ r, ξ n = n r, ξ n−2 (r, ξ ·r, ξξ ) is required. Then, the product and chain rule leads to e t , s =r, ss = d ds r, ξ r, ξ −1 = 1 r, ξ 4 r, ξ 2 r, ξξ −(r, ξ ·r, ξξ )r, ξ , see Eqs. (5) and (7). This expression is necessary to determine the curvature in Eq. (13) 1 by calculating the cross product with e t e t , s ×e t = κe n × e t = −κe b , i.e. using the norm of both sides, we obtain after some calculations κ = e t , s ×e t = 1 r, ξ 3 r, ξ ×r, ξξ .

Measures defining a straight line (fiber)
Short glass or carbon fibers can be modeled as straight lines. For a straight line,r(s) = r 1 + se t , we have κ = 0, and τ is defined to be zero (e n is in this case arbitrary, but orthonormal to e t ), i.e. only the length L and the orientation e t are chosen for characterizing purposes. The orientation is given by Eq. (8),r, s = e t , see Fig. 2. Thus, we need the starting and the end point of a fiber, r 1 and r 2 , respectively, The orientation vector e t can be expressed in a spherical coordinate system by the two angles ϕ and ϑ, see, for example, [12]. Frequently, the two angles are chosen for characterizing the orientation of fibers in composites, see, for example, [1,17]. Thus, straight fibers are simple to be evaluated. Since we are interested in scalar measures representing some "mean-value" of measures characterizing curvilinear fibers, we introduce two approaches. First, given data points representing a spatial curve are interpolated to obtain a continuous representation (piece-wise C ∞ , but globally C 2 continuous). Based on these approximations, the measures (mean curvature, torsion, direction, and deviation from straight line) are calculated. Afterwards, these measures are provided without the interpolation concept, i.e. only the data points are evaluated directly.

Mean values using continuous interpolation
Before defining the mean values of interest, a fiber-given by n d data points d k , k = 1, . . . , n d -should be represented by an interpolating function, As a first approach, piece-wise polynomials of the third order are chosen, where the sampling points ξ i are estimated using the Chord-approach with ξ 1 = 0. The index i defines the interval, i = 1, . . . , n d − 1, and k = 1, 2, 3 the coordinate direction.
In other words, we havex k (ξ ) =x ki (ξ ) in interval i, ξ i ≤ ξ < ξ i+1 . Obviously, we obtain for ξ = ξ i the parameters of each coordinate, k = 1, 2, 3 (starting point implies r(ξ 1 ) = d 1 ). Using this approach, the first three derivatives can be determinedx In each interval i, four unknowns are given, i.e. we have 4 × (n d − 1) unknowns. Thus, we need the same number of equations. With Eq. (31) and the end point r(ξ n d ) = d n d , there are n d equations. Furthermore, we require that the following conditions are fulfilled at the point ξ i+1 i = 1, . . . , n d − 2, leading to 3 × (n d − 2) equations. Thus, two additional equations are required. Here, we draw on natural splines, i.e. we assume curvature free starting and ending pointŝ The index ki + 1 in Eq. (33) defines the kth component (or its derivatives), see Eq. (28), in interval i + 1. To avoid a comma separation, which is chosen for derivatives, a comma-free notation is chosen. This holds for kn d − 1 in Eq. (34) as well, i.e. the evaluation of the kth function defined in interval n d − 1. The numerical implementation to solve the linear system of equations is provided in [5], see [4] as well. In other words, the given data points are exactly fulfilled. The curves are continuous of second order. In conclusion, the parameters a ki , b ki , c ki ,d ki , i = 1, n d − 1, k = 1, 2, 3, are assumed to be known in the following. This implies that in each interval the position vector and its derivatives In the following, we define the measures determining fibers. First, the length of a fiber has to be computed using definition (3) The integrals can, for example, be evaluated using Simpson's rule, see, for instance, [16], with Δξ = (ξ i+1 − ξ i )/n S , n S is the number of subintervals between two data points. Then, we have to calculated at the integration points ξ (l) = ξ i + (l − 1)Δξ , In the following, all integrals are evaluated using this integration scheme. (Of course, other schemes or adaptive quadrature rules are possible.) Next, the mean curvature the mean torsion An additional quantity is the deviation of the curved line to a straight line. The straight line is defined by see, for instance, [10], so that the mean value is defined by In Particularly, the length might be of interest to estimate-together with a mean diameter-the volume of the fiber, or the orientation vector n = r M , ξ / r M , ξ in Eq. (40) is drawn on to determine an orientation tensor N = n ⊗ n, see [15]. However, new measures such as mean curvature or mean deviation to a straight line give insight into the flow behavior and the weakness of natural fibers relative to the matrix material.

Mean values using simplified approach
Apart from the more precise definitions using interpolated points in Sect. 2.3.1, an alternative, simplified access drawing on rough approximations of the derivatives by numerical differentiation and the integrals applying a trapezoidal rule can be provided.
Again the chord method (30) is applied using In a first step, we compute at all data points i = 1, . . . , n d the first three derivatives using the approximations in "Appendix A". The tangent vector reads, according to Eq. (65), with The tangent vectors at the first and the last point can be estimated using Eq. (67), The second derivative reads according to Eq. (70) Center

Example
In the following, we apply the two concepts to real μ-CT data of a man-made cellulose fiber reinforced polypropylene specimen. The data are generated using a 3D X-ray microscope Xradia 520 Versa from Zeiss. The fibers were detected with the program Avizo (version 9.4.0 of the company FEI). The considered volume element of the sample has a size of approximately 208.5 µm × 629.6 µm × 629.6 µm. This produces the fibers shown in Fig. 4. Here, we are not interested in discussing whether it makes sense considering fibers which are cut at the edge surfaces of the cube-like region or not. We are interested in how the applied schemes work with some data.
To obtain more insight in the fiber distribution, we draw on the concepts of Sects. 2.3.1 and 2.3.2. Figures 5  and 6 compile the results. The length computation shows very similar results of both methods, the interpolation concept, and the simplified scheme, see Fig. 5a and b.
The mean curvature in Fig. 5c and d, however, shows slightly different results, since the simplified method draws on the second derivative, and derivatives using differential formulas roughen the results. This is the essential drawback in the mean torsion computation shown in Fig. 5e-f, where the third derivative of the curve vector is evaluated. Since there are some fibers with "kinks", totally unrealistic values are computed using the simplified version to determine the torsion of a fiber. The estimation of the mean fiber orientation, represented by the angles ϑ M and ϕ M , is appropriate for both schemes, see Fig. 6a-d, which holds for the mean distance of the fibers relative to the center line (42) in Fig. 6e-f.

Conclusions
In this paper, the quantitative properties of man-made cellulose fibers are proposed. Since these fibers represent curvilinear functions in the matrix material, some characteristic measures are of particular interest. We propose several measures describing curvilinear lines: length, mean curvature, mean torsion, mean direction, and a mean distance to a straight line. To describe the curvilinear fibers, the discrete data points are interpolated by a first concept. This interpolation is based on splines to obtain a continuous representation. In a second scheme, only the discrete data points are evaluated. There, the required first three derivatives are based on differential derivative formulas leading to a much simpler evaluation process. However, it turns out that the third derivative leads to very inadequate results in some situations so that the measure of mean torsion, i.e. a measure of the deviation of a curve from a plane, is not suitable. All other measures lead-for the distribution of a real μ-CT scan data-to similar results. Thus, both methods are a promising tool for evaluating the fiber measures, except for torsion, where only the more sophisticated evaluation procedure is necessary. The third derivative can be computed by inserting the second derivative into Eq. (65) instead of the function values (quadratic interpolation) i = 3, . . . , n d − 2; thus, five points are required at least to estimate the third derivative (the values f i+1 requires y i+2 and f i−1 requires y i−2 ) if the differentials (70) would be used. Since we have estimated the second derivatives at the first and last point, see Eq. (71), formula (72) can be evaluated also in the range i = 2, . . . , n d − 1. The third-order derivatives at the "starting" and "ending" points are estimated by (linear interpolation) Therefore, at each point ξ i , i = 1, . . . , n d , a derivative can be estimated.