Fabric Tensor Characterization of Tensor-Valued Directional Data: Solution, Accuracy, and Symmetrization
Fabric tensor has proved to be an effective tool statistically characterizing directional data in a smooth and frame-indifferent form. Directional data arising from microscopic physics and mechanics can be summed up as tensor-valued orientation distribution functions (ODFs). Two characterizations of the tensor-valued 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 tensor-valued ODFs are presented in terms of orthogonal irreducible decompositions in both two- and three-dimensional (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 tensor-valued ODF and fabric tensor in continuum damage mechanics is presented.
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 frame-indifferent form, or a tensorial form which is invariant to coordinate transformations; see, for example, Kanatani , Advani and Tucker , and Yang et al. .
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 stereo-angle along a certain direction. The classically defined ODFs can be abstracted as scalar-valued directional functions, because from a vector-space viewpoint, their values on each orientation are nothing but normalized scalars. Therefore, the classical ODFs are referred to as scalar-valued ODFs in this paper, following Yang et al. , in regard to the generality of tensor-valued ODFs.
The scalar-valued 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 , 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 , who characterized the scalar-valued 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.  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 tensor-valued 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 vector-valued ODF, which cannot be fully described by a Cauchy-like stress tensor of order two. In fabric-tensor-based micromechanics of granular media; see, for example, Rahmoun et al.  and Li and Yu , the contact forces and contact vectors should both be taken for vector-valued ODFs. In the microplane theory of Bažant and Gambarova  and Carol and Bažant , the microplane elasticity tensor is essentially a second-order tensor-valued ODF that defines a one-to-one linear map between the stress and strain vector-valued ODFs. The fabric tensor algebra of vector-valued ODFs is preliminarily studied by Yang et al. , who derived the explicit expressions of symmetric fabric tensors up to the sixth order. Li and Yu  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  suggested the test of uniformity to evaluate the fitness. Li and Yu  simply defined the proximity index for both scalar- and vector-valued ODFs. One can intuitively estimate the accuracy of the zeroth-, second-, and fourth-order fabric tensor approximations to randomly distribute scalar-valued 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 scalar-valued ODFs.
In this paper, we complete the fabric tensor algebra for the largest generality of tensor-valued ODFs. A strict but asymmetric characterization is presented in complete parallel to that of the scalar-valued ODFs by Kanatani . The relationship of fabric tensors of different orders  holds for the tensor-valued 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 tensor-valued ODFs. Analytic solutions of the two types of fabric tensors characterizing centrosymmetric and anticentrosymmetric tensor-valued ODFs in 2D and 3D are presented in the forms of orthogonal irreducible decompositions by Zheng and Zou . 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 fourth-order 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 second-order identity tensor and means the th power tensor product of : The alternating product is actually a modified inner product in which the left-hand 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 , 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 , which states that any symmetric tensor in either 2D or 3D can be orthogonally decomposed into deviator-tensor-induced irreducible bases in the following unified forms:where and are two clusters of deviator tensors. The deviatorization of any deviator-tensor-induced 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 Tensor-Valued ODFs
A tensor-valued ODF is a tensor function dependent on the predefined orientation . A generic tensor-valued 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 tensor-valued ODF . Let be approximated by another tensor-valued 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 tensor-valued ODF if the following condition of least-square 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 tensor-valued ODFs is the intrinsic extension of that of scalar-valued ODFs by Kanatani , in that the components of the tensor-valued ODF are handled as independent scalar-valued 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  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 left-hand 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 ,
For the anticentrosymmetric tensor-valued 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 tensor-valued 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 tensor-valued 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 .
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 scalar-valued ODFs is sufficient to evaluate the fitness of the asymmetric characterization, because each component of the tensor-valued ODFs is characterized in terms of independent scalar-valued ODFs. Nevertheless, vector-valued ODFs are investigated in the next section to calibrate the accuracy of the symmetric characterization.
The original scalar-valued 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 fourth-order (). 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.
4. Symmetric Characterization of Tensor-Valued ODFs
For quite a few non-scalar-valued ODFs in practice, their values on all orientations are objectively symmetric, that is, . Such circumstances include all the vector-valued ODFs and many second-order tensor-valued 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 least-square 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),
For vector-valued ODFs (), the closed-form expression of the coefficient can be shown as However, for generic tensor-valued ODFs of higher orders, the closed-form 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 tensor-valued 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
4.2. Nonconvergence of the Symmetric Characterization
The symmetric characterization is proved to be nonconvergent simply by a counter example. Assume an anticentrosymmetric vector-valued ODF that is precisely determined by its second-order symmetrized fabric tensor,
As demanded by the convergence condition, also must be precisely determined by its fourth-order symmetrized fabric tensor . Applying (4.19) to (4.18), one obtains the fourth-order 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 fourth-order fabric tensor of the second kind ,
Despite of such non-convergence, the symmetric characterization is still available for most irregular tensor-valued 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 vector-valued and second-order tensor-valued ODFs, because the symmetric characterization is identical to the asymmetric one for scalar-valued ODFs. The original vector-valued 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 second-order tensor-valued 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 vector-valued 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.
5. Application: Microplane-Based Damage Effective Stress
In this section, we show an application of tensor-valued 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  and Rabotnov  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 cross-sectional area subjected to a uniaxial stress . The uniaxial damage effective stress is defined by Kachanov  and Rabotnov : where is the scalar damage variable and the damaged cross-sectional area.
The multiaxial generalization of the Kachanov-Rabotnov 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  into five types. All these models start from a subjective premise that the damage effective stress is a Cauchy-like 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 vector-valued ODF owing to its microplane-based 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 second-order 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 scalar-valued ODF , Obviously, because the term is nontensorial, the second-order 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 second-order damage effective stress tensor is so controversial; see, for example, the discussion by Zheng and Betten , Kachanov , and Lemaitre et al. .
In fact, the second-order tensor can be identified with the second-order fabric tensor of the vector-valued 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 fourth-order symmetrized fabric tensors of in 3D are given bywhere and are the symmetrized fabric tensors of the first kind,
However, the definition of the damage effective stress ODF in (5.3) is quite narrow, which only deals with coaxial damage effects; that is, damage only affects the net area of each microplane without changing its orientation. In general, noncoaxial damage effects should be described by a second-order tensor-valued ODF , which connects the damage effective stress ODF to the Cauchy stress ODF on each microplane with full geometry equivalence. In noncoaxial cases, it can be shown that if the second-order fabric tensor is applied to characterize , the widely applied anisotropic model based on the fourth-order damage effect tensor [16, 29, 32, 33, 35, 36] in Voigt symmetry can be realized. The noncoaxial cases as well as microplane-based physics equivalence hypotheses and damage evolution laws will be elaborated in another paper, which may help to clarify some uncertainties over the basic concepts of continuum damage mechanics.
The fabric tensor approach for the largest generality of tensor-valued orientation distribution functions (ODFs) is presented, including the analytical solution and accuracy analysis of the asymmetric and symmetric characterizations. The asymmetric characterization of tensor-valued ODFs is in complete parallel to that of scalar-valued ODFs by Kanatani , where the components of tensor-valued ODFs are handled as independent scalar-valued ODFs. The asymmetric characterization proves to be convergent  and featured by the relationship of fabric tensors of different orders . The symmetric characterization does not possess such convergence property and relationship, but is still an available solution for most irregular tensor-valued ODFs whose fabric tensors are prescribed in full symmetry. Accuracy analyses are performed on normally distributed random ODFs to evaluate the approximation quality of the two characterizations. The results indicate that the fitness is dominated by the dispersion degree of the original ODFs rather than the orders of the fabric tensors. The accuracy and convergence rate of the symmetric characterization are lower than those of the asymmetric one, but the differences keep in an acceptable range for lower approximation orders.
This work is supported by the National Natural Science Foundation of China under Projects 50925931 and 11172150 and subsidized by the special funds for major state basic research projects with Grant no. 2009CB724604.
H. J. Bunge, “Fabric analysis by orientation distribution functions,” Tectonophysics, vol. 78, no. 1–4, pp. 1–21, 1981.View at: Google Scholar
J. H. Bunge, Texture Analysis in Materials Science: Mathematical Methods, Butterworths, London, UK, 1982.
M. Oda, “Fabric tensors for discontinuous geological materials,” Soils and Foundations, vol. 22, no. 4, pp. 96–108, 1982.View at: Google Scholar
V. A. Lubarda and D. Krajcinovic, “Damage tensors and the crack density distribution,” International Journal of Solids and Structures, vol. 30, no. 20, pp. 2859–2877, 1993.View at: Google Scholar
P. D. Wu and E. Van Der Giessen, “On improved network models for rubber elasticity and their applications to orientation hardening in glassy polymers,” Journal of the Mechanics and Physics of Solids, vol. 41, no. 3, pp. 427–456, 1993.View at: Google Scholar
Z. P. Bažant and P. G. Gambarova, “Crack shear in concrete: crack band microplane model,” Journal of Structural Engineering, vol. 110, no. 9, pp. 2015–2035, 1984.View at: Google Scholar
I. Carol and Z. P. Bažant, “Damage and plasticity in microplane theory,” International Journal of Solids and Structures, vol. 34, no. 29, pp. 3807–3835, 1997.View at: Google Scholar
L. M. Kachanov, “Time of the Rupture Process under Creep Conditions,” Izvestiya Sibirskogo Otdeleniya,. Akademii Nauk SSSR., Seriya Tekhnicreskikh Nauk, vol. 8, pp. 26–31, 1958.View at: Google Scholar
I. N. Rabotnov, “On the equations of state for creep,” Progress in Applied Mechanics, the Prager Anniversary Volume, vol. 8, pp. 307–315, 1963.View at: Google Scholar
Q. S. Zheng and J. Betten, “On damage effective stress and equivalence hypothesis,” International Journal of Damage Mechanics, vol. 5, no. 3, pp. 219–240, 1996.View at: Google Scholar
G. Z. Voyiadjis and T. Park, “Anisotropic damage effect tensors for the symmetrization of the effective stress tensor,” Journal of Applied Mechanics, Transactions ASME, vol. 64, no. 1, pp. 106–110, 1997.View at: Google Scholar
S. Murakami, “Mechanical modeling of material damage,” Journal of Applied Mechanics, Transactions ASME, vol. 55, no. 2, pp. 280–286, 1988.View at: Google Scholar
J. L. Chaboche, “Continuous damage mechanics—a tool to describe phenomena before crack initiation,” Nuclear Engineering and Design, vol. 64, no. 2, pp. 233–247, 1981.View at: Google Scholar
D. Krajcinovic, “Continuous damage mechanics,” Applied Mechanics Review, vol. 37, pp. 1–6, 1984.View at: Google Scholar