Research Article  Open Access
Fabric Tensor Characterization of TensorValued Directional Data: Solution, Accuracy, and Symmetrization
Abstract
Fabric tensor has proved to be an effective tool statistically characterizing directional data in a smooth and frameindifferent form. Directional data arising from microscopic physics and mechanics can be summed up as tensorvalued orientation distribution functions (ODFs). Two characterizations of the tensorvalued ODFs are proposed, using the asymmetric and symmetric fabric tensors respectively. The later proves to be nonconvergent and less accurate but still an available solution for where fabric tensors are required in full symmetry. Analytic solutions of the two types of fabric tensors characterizing centrosymmetric and anticentrosymmetric tensorvalued ODFs are presented in terms of orthogonal irreducible decompositions in both two and threedimensional (2D and 3D) spaces. Accuracy analysis is performed on normally distributed random ODFs to evaluate the approximation quality of the two characterizations, where fabric tensors of higher orders are employed. It is shown that the fitness is dominated by the dispersion degree of the original ODFs rather than the orders of fabric tensors. One application of tensorvalued ODF and fabric tensor in continuum damage mechanics is presented.
1. Introduction
Many quantities arising from microscopic physics and mechanics are directional dependent, and as is often the case, are randomly and discontinuously distributed along the orientations. Statistical characterization of such directional data is essential [1–4]. With regard to physical problems, such characterization must take a frameindifferent form, or a tensorial form which is invariant to coordinate transformations; see, for example, Kanatani [5], Advani and Tucker [6], and Yang et al. [7].
Directional data can be generalized as orientation distribution functions (ODFs). In its classical definition [3–5, 8, 9], an ODF is actually a probability density function of a continuous random orientational variable, which practically represents the fraction of the total number of directional elements per unit stereoangle along a certain direction. The classically defined ODFs can be abstracted as scalarvalued directional functions, because from a vectorspace viewpoint, their values on each orientation are nothing but normalized scalars. Therefore, the classical ODFs are referred to as scalarvalued ODFs in this paper, following Yang et al. [7], in regard to the generality of tensorvalued ODFs.
The scalarvalued ODFs can be statistically characterized by their fabric tensors [5, 6, 9–21]. Such tensors are labeled fabric because they were first and have been extensively introduced to describe the fabric (or texture) of structured solid materials, such as granular [5, 10–12] and microcracked [13–16] media, polycrystals and polymers [17, 18], and bones [19], and prove to be an effective tool bridging the multiscale behaviors of the materials. Recent applications of fabric tensors extend beyond material studies to wherever directional variables are focused, such as mixed fluids [6, 20] and medical imaging [9, 21], whether being called fabric tensors or not. The analytic solutions of fabric tensors of higher orders are obtained by Kanatani [5], who characterized the scalarvalued ODFs in expansions of three kinds of fabric tensors and addressed their equivalence to the Fourier series expansion in 2D and to the spherical harmonics expansion in 3D. Yang et al. [22] showed that fabric tensors of different orders are not independent: fabric tensors of lower orders can be precisely determined by those of higher orders.
The concept of the tensorvalued ODFs is relatively new, but their objectivity in physics and mechanics is visible. Yang et al. [7, 23] proved that the damage effective stress in damage mechanics is inherently a vectorvalued ODF, which cannot be fully described by a Cauchylike stress tensor of order two. In fabrictensorbased micromechanics of granular media; see, for example, Rahmoun et al. [11] and Li and Yu [12], the contact forces and contact vectors should both be taken for vectorvalued ODFs. In the microplane theory of Bažant and Gambarova [24] and Carol and Bažant [25], the microplane elasticity tensor is essentially a secondorder tensorvalued ODF that defines a onetoone linear map between the stress and strain vectorvalued ODFs. The fabric tensor algebra of vectorvalued ODFs is preliminarily studied by Yang et al. [7], who derived the explicit expressions of symmetric fabric tensors up to the sixth order. Li and Yu [12] provided a stepwise repetitive method to determine asymmetric fabric tensors of higher orders. However, general analytic solution is not found in previous studies.
The approximation quality of the fabric tensor characterization is not so frequently stressed in previous studies, perhaps because there are few alternative ways to characterize directional data in a tensorial form. Kanatani [5] suggested the test of uniformity to evaluate the fitness. Li and Yu [12] simply defined the proximity index for both scalar and vectorvalued ODFs. One can intuitively estimate the accuracy of the zeroth, second, and fourthorder fabric tensor approximations to randomly distribute scalarvalued ODFs in many of the referred studies [5, 6, 9, 11, 13, 14, 16]. Nevertheless, systematic accuracy analysis has not been reported yet, even for scalarvalued ODFs.
In this paper, we complete the fabric tensor algebra for the largest generality of tensorvalued ODFs. A strict but asymmetric characterization is presented in complete parallel to that of the scalarvalued ODFs by Kanatani [5]. The relationship of fabric tensors of different orders [22] holds for the tensorvalued ODFs in this case. However, the asymmetric fabric tensors may be inconvenient for practical use because of potential inconsistency with physical hypotheses and observations. Therefore, a nonstandard symmetric characterization is proposed, which is proved to be nonconvergent but still available for most irregular tensorvalued ODFs. Analytic solutions of the two types of fabric tensors characterizing centrosymmetric and anticentrosymmetric tensorvalued ODFs in 2D and 3D are presented in the forms of orthogonal irreducible decompositions by Zheng and Zou [26]. Accuracy analysis is performed on normally distributed random ODFs to evaluate the approximation quality of the two characterizations, where fabric tensors of higher orders are employed.
2. Notation and Operation of Tensor and Orientation
In this paper, we mainly use the abstract notation of tensors to simplify the expressions. For this purpose, quite a few notations and operators must be predefined. The abstract form of a generic tensor is written in bold and nonitalic, with its order shown in the subscript, for example, a fourthorder tensor (Cartesian components in ). In particular, scalars are written like and vectors like .
2.1. Tensor Symmetrization
In the component notation, the symmetrization of a generic tensor on the indices is denoted by , which is defined as the sum of the components of indexed by the permutations of the overbared subscripts divided by If two symmetrizations are required at the same time (luckily, two is enough for our use), double overbars must be used. Here are two examples:
In the abstract notation, the symmetrization is simply represented by the overbared orders, for example, , , and . The operator means sequenced addition of orders, which removes ambiguities when the orders are operated. For example, precisely indicates that the second and third orders are symmetrized, while may be mistaken for or .
Full symmetrization of a tensor can be realized by defining the alternating product operation signified by , where is the secondorder identity tensor and means the th power tensor product of : The alternating product is actually a modified inner product in which the lefthand operand tensor participates on alternate orders. Equation (2.2) makes itself explicit in the component form
2.2. Tensor Deviatorization and Orthogonal Irreducible Decomposition
In this subsection, we only consider full symmetric tensors. Successive contraction of a full symmetric tensor for times is denoted by
A symmetric tensor is termed a deviator tensor if it is traceless, . The deviatorization operation of a symmetric tensor means to segregate its deviator part, as signified by the operator on its order,where is the th order identity tensor, . The coefficients are determined in such a way that and any contraction of (2.6a) and (2.6b) vanishes. The results of are provided by Kanatani [5], where is the order of , whether even or odd. For examples, the deviatorizations of the 3D symmetric tensors of orders from 1 to 4 are shown in the component forms
The orthogonal irreducible decomposition theorem of symmetric tensors of higher orders is proposed by Zheng and Zou [26], which states that any symmetric tensor in either 2D or 3D can be orthogonally decomposed into deviatortensorinduced irreducible bases in the following unified forms:where and are two clusters of deviator tensors. The deviatorization of any deviatortensorinduced irreducible basis vanishes except for the last term where , namely, () and .
2.3. Definition and Integral of Orientation
Any orientation about a point can be defined by a unit vector starting at that point and ending on a unit circle denoted in 2D or a unit sphere denoted in 3D, where and are the stereo angles defining the orientation .
The contour integral in 2D and 3D are, respectively, defined as Noting that and are the circumference of and the spherical area of , respectively, so the integral represents an average value of the distribution over all orientations.
The identity tensors are related to the integral of orientations by where and the coefficient is given by Noting that , the successive contraction of for times takes the form For example, in 3D, , , and .
3. Asymmetric Characterization of TensorValued ODFs
A tensorvalued ODF is a tensor function dependent on the predefined orientation . A generic tensorvalued ODF of order is denoted by , which may be randomly and discontinuously distributed along the orientations. The fabric tensor characterization first requires that be decomposed into the centrosymmetric and anticentrosymmetric parts namely, , where
We first consider the centrosymmetric tensorvalued ODF . Let be approximated by another tensorvalued ODF , which is precisely determined by a tensor (even orders must be added due to the centrosymmetry of ), for example, . Because is full symmetric, it is quite visible that .
The tensor is defined as the fabric tensor of the second kind (or simply referred to as fabric tensor where no ambiguity arises) of the tensorvalued ODF if the following condition of leastsquare minimization of error is imposed: Applied with (3.2), (3.3) yields which can be further recast as where is defined as the fabric tensor of the first kind of , One contracting property of is favorable, .
Such asymmetric characterization of tensorvalued ODFs is the intrinsic extension of that of scalarvalued ODFs by Kanatani [5], in that the components of the tensorvalued ODF are handled as independent scalarvalued ODFs, whose fabric tensor expansion proves to be equivalent to the Fourier series expansion in 2D and to the spherical harmonics expansion in 3D. In other words, the asymmetric characterization is convergent; namely, the approximation error is a nonincreasing function of the used order and tends to zero when tends to infinity.
3.1. Analytic Solution of
Evidently, can be approximately determined by observed values of on finite numbers of orientations, so what left for us to do is to find the explicit expression of with respect to , or the inversed form of (3.5).
Resorting to (2.9a), the orthogonal irreducible decomposition of takes the following form: where () are defined by [5] as the fabric tensors of the third kind of , which are, however, barely perceptible in practical use.
Applying (3.7) to (3.5), and noting that any contraction of vanishes, one obtains where the coefficient is determined by (if one substitutes (3.7) to (3.5) and observes the result, he may realize that is all about a ball and box puzzle. Assume that distinct balls, colored red and blue, are paired and placed into sequenced boxes, each box with a pair. And is actually the probability of that boxes are filled with a pair in different colors)
Successively contracting (3.8) for times on the symmetrized orders, taking in the property , where is the product of the coefficients arising from each contraction,
Taking the deviatorization of (3.10) on the remaining symmetric orders, noting that any term in the lefthand side vanishes except the last one, one finally acquires the fabric tensors of the third kind of , The coefficient is fully simplified by taking in (2.13), which turns out to be independent of ,
The analytic solution of fabric tensors of the second kind of is finally obtained by applying (3.12) and (2.6a) to (3.7),
For the anticentrosymmetric tensorvalued ODF , odd orders must be appended to the fabric tensors of the second kind,
Accordingly, the fabric tensors of the first kind are defined as
The analytic solution of with respect to can be acquired completely parallel to that of , which takes the following form: where is calculated just by replacing with in (3.13).
For the generic tensorvalued ODF , its centrosymmetric part and anticentrosymmetric part are characterized independently, which means that the order of is unnecessarily one larger than that of (anyway, we take this assumption in the accuracy analysis). For quick reference, the ultimately solved coefficients of the fabric tensors of the second kind up to the sixth and seventh approximation orders are listed in Table 1.

3.2. Relationship between Fabric Tensors of Different Orders
Assume that two fabric tensors and , where , are used to characterize the centrosymmetric tensorvalued ODF In view of (3.4), the fabric tensors of the first kind defined by (3.6) take the following form:The following relation between and is evident: where the second equality is achieved by contracting (3.19b) for times.
Equation (3.20) indicates that once a fabric tensor of the original ODF is obtained, those of lower orders can be precisely determined by the approximate ODF induced from , needless to use the original ODF again. Simply put, fabric tensors of lower orders can be precisely determined by those of higher orders. Such dependence implies that fabric tensors of different orders cannot be treated as independent state variables in thermodynamics [22].
3.3. Accuracy Analysis
The approximation quality of the fabric tensor characterization is quantified by the squared error ratio defined as where the superscript determines the approximation orders applied for both the centrosymmetric and anticentrosymmetric parts of ,
Accuracy analysis performed on scalarvalued ODFs is sufficient to evaluate the fitness of the asymmetric characterization, because each component of the tensorvalued ODFs is characterized in terms of independent scalarvalued ODFs. Nevertheless, vectorvalued ODFs are investigated in the next section to calibrate the accuracy of the symmetric characterization.
The original scalarvalued ODF is assigned on each direction with a random real number sampling from the normal distribution with mean and standard deviation , denoted . In the accuracy analysis, is fixed at 1.0 and ranges from 0.1 to 0.5 with increment 0.1. The order parameter ranges from 0 to 5, namely, as high as the tenth, and eleventh fabric tensors are involved. The squared error ratios in 2D and 3D are shown in Figures 1 and 2, respectively.
As evidenced by the figures, the increase of will to a considerable extent reduce the accuracy when the approximation order is at a lower level, for example, the mostly used second and fourthorder (). The loss of accuracy can be remedied by increasing the order to a higher level such as the tenth order (), which, however, cannot be effective due to upsurging computational expense. Therefore, the approximation quality of the fabric tensor characterization is dominated by the dispersion degree of the original ODFs rather than the orders of fabric tensors. The original and approximate ODFs related to distribution in 2D and 3D are plotted in Figures 3 and 4, respectively, from which the fitness can be intuitively estimated.
(a)
(b)
(c)
(d)
4. Symmetric Characterization of TensorValued ODFs
For quite a few nonscalarvalued ODFs in practice, their values on all orientations are objectively symmetric, that is, . Such circumstances include all the vectorvalued ODFs and many secondorder tensorvalued ODFs in physics, such as the elasticity ODF and the conservative diffusion gradient ODF . The symmetry of the fabric tensors of such ODFs is obvious, and . On many occasions, however, these fabric tensors are supposed to be fully symmetrized to keep consistency with physical hypotheses or observations. Therefore, the symmetric characterization is proposed, where the fabric tensors are fully symmetrized.
Consider again that the centrosymmetric part of denoted , which is characterized by the symmetrized fabric tensor of the second kind , The condition of leastsquare minimization of error goes as follows:
Noting that is related to by the alternating product defined in (2.2), (4.2) can be recast as which can be further reduced to where is termed the symmetrized fabric tensor of the first kind, Note that does not have the contracting property, that is, , which means that the relationship between fabric tensors of different orders does hold for the symmetric characterization. Equation (4.4) looks somewhat confusing, so we just show its component form for and ,
4.1. Analytic Solution of
Likewise, our task is to find the analytic solution of with respect to , or the inversed form of (4.4).
Known from (4.4), is also a full symmetric tensor (just look at (4.6)). To avail of the orthogonal irreducible decomposition, the parity of must be discussed. We only consider the even case in the derivation. When is an even integer, let , and the decomposition of takes the following form: where () are accordingly termed the symmetrized fabric tensors of the third kind, of which any contractions vanish.
Applying (4.7) to (4.4), and after some manipulations, one obtains Note that each term in the summation takes the following recursive pattern (it is actually another more complicated ball and box puzzle),
Thus, (4.8) can be reduced to where the coefficient is given by and is a recursive function drawn from (4.9), compatible for any parities of the parameters,
The solving of (4.10) is through successive contraction and deviatorization, identical to that of (3.8), which yields the symmetrized fabric tensors of the third kind,
For vectorvalued ODFs (), the closedform expression of the coefficient can be shown as However, for generic tensorvalued ODFs of higher orders, the closedform expression of is not found because of greater recursion depth.
The analytic solution of the symmetrized fabric tensors of the second kind is finally obtained, considering the parity of the approximation order,
The anticentrosymmetric tensorvalued ODF is characterized by where the symmetrized fabric tensors of the second kind are solved aswith the symmetrized fabric tensors of the first kind defined as
The ultimately solved coefficients of the symmetrized fabric tensors of the second kind of lower orders for vectorvalued ODFs and secondorder tensorvalued ODFs are listed in Tables 2 and 3.
4.2. Nonconvergence of the Symmetric Characterization
The symmetric characterization is proved to be nonconvergent simply by a counter example. Assume an anticentrosymmetric vectorvalued ODF that is precisely determined by its secondorder symmetrized fabric tensor,
As demanded by the convergence condition, also must be precisely determined by its fourthorder symmetrized fabric tensor . Applying (4.19) to (4.18), one obtains the fourthorder fabric tensor of the first kind in 3D, which is then substituted to (4.17b) (or referring to the second row of the right half of Table 2) to determine the fourthorder fabric tensor of the second kind ,
Finally, applying (4.21) to (4.16), one can find the convergence condition undermined because
Despite of such nonconvergence, the symmetric characterization is still available for most irregular tensorvalued ODFs. On one hand, the approximation condition and accuracy are definitely determined for a given order, while the accuracy may well rise with the used order owing to more independent components in the fabric tensors, as illustrated in the next subsection. On the other hand, if full symmetry is prescribed on the fabric tensors, the symmetric characterization is of the highest accuracy, particularly much more accurate than the ungrounded direct symmetrization of the asymmetric fabric tensors.
4.3. Accuracy Analysis
Accuracy analyses are performed on vectorvalued and secondorder tensorvalued ODFs, because the symmetric characterization is identical to the asymmetric one for scalarvalued ODFs. The original vectorvalued ODF is generated in the following way: the magnitude samples from , with ranging from 0.1 to 0.5, and the stereo angles and from and , respectively. The secondorder tensorvalued ODF is formed by . Both the asymmetric and symmetric characterizations are investigated and compared.
The results indicate that the symmetric characterization still follows the general rule of the asymmetric one; that is to say, the fitness majorly decreases with the dispersion degree of the original data () and improves only when the approximation order reaches very high. So we only display and discuss the results for . The squared error ratios of the two characterizations in 2D and 3D are shown in Figures 5 and 6, respectively.
Known from Figures 5 and 6, the accuracy and convergence rate of the symmetric characterization are lower than those of the asymmetric one, because of fewer independent components of cause. But the difference does not reach an unacceptable degree, especially for lower orders. Therefore, the symmetric characterization is by all means an available solution for where fabric tensors are required in full symmetry; after all, it is the most accurate way to obtain full symmetric fabric tensors. The original and approximate vectorvalued ODFs related to in 2D and 3D are plotted in Figures 7 and 8, respectively, where the magnitude and the cosine of the deviation angle serve as the radius and color functions, respectively.
(a)
(b)
(c)
(d)
5. Application: MicroplaneBased Damage Effective Stress
In this section, we show an application of tensorvalued ODF and fabric tensor in continuum damage mechanics, which helps to clarify the multiaxial generalization of the damage effective stress in the classical models. Continuum damage mechanics deals with the degradation of solid materials due to small defects (voids, microcracks, and dislocations) using continuity modeling. Kachanov [27] and Rabotnov [28] raised the concept of damage effective stress, which constitutes the core of continuum damage mechanics.
We start from the original concept of uniaxial damage. Consider a smooth tensile bar of nominal crosssectional area subjected to a uniaxial stress . The uniaxial damage effective stress is defined by Kachanov [27] and Rabotnov [28]: where is the scalar damage variable and the damaged crosssectional area.
The multiaxial generalization of the KachanovRabotnov damage effective stress is not a straightforward task, especially for anisotropic damage effects. Quite a few models have been developed, as summarized by Zheng and Betten [29] into five types. All these models start from a subjective premise that the damage effective stress is a Cauchylike stress tensor of order two, which is critically reviewed by Yang et al. [7, 23], who proposed that damage effective stress should be generalized as a vectorvalued ODF owing to its microplanebased nature.
A microplane is a plane of any orientation cutting the material at a given point, which is represented by the unit normal of that orientation, . The stress state on a microplane is described by the Cauchy stress ODF , which is precisely determined by the secondorder Cauchy stress tensor ,
Equation (5.1) can be naturally extended onto each microplane to define the damage effective stress ODF , with anisotropic damage effects measured by a scalarvalued ODF , Obviously, because the term is nontensorial, the secondorder tensor cannot fully describe the damage effective stress ODF but is only an approximate estimation. This is the essential reason why the definition of a secondorder damage effective stress tensor is so controversial; see, for example, the discussion by Zheng and Betten [29], Kachanov [30], and Lemaitre et al. [31].
In fact, the secondorder tensor can be identified with the secondorder fabric tensor of the vectorvalued ODF , and of course, the order can reach higher than two. Because only the symmetric part of accounts for constitutive equations [32–34], the symmetric characterization should be applied. For example, the second and fourthorder symmetrized fabric tensors of in 3D are given bywhere and are the symmetrized fabric tensors of the first kind,