Research Article | Open Access
Jiangu Qian, Jean-Pierre Bardet, Maosong Huang, "Spectral Classification of Non-Coaxiality for Two-Dimensional Incremental Stress-Strain Response", Mathematical Problems in Engineering, vol. 2010, Article ID 963043, 20 pages, 2010. https://doi.org/10.1155/2010/963043
Spectral Classification of Non-Coaxiality for Two-Dimensional Incremental Stress-Strain Response
The present study examines the non-coaxial aspects of incremental material behavior, and attempts to classify the incremental non-coaxiality that relates stress and strain increments. In the solid mechanics literature, non-coaxiality (NC) refers usually to incremental strains and stress states having different principal directions. Departing from conventional non-coaxiality, the analysis investigates the incremental non-coaxiality (INC) of linearized rate-type solids. This study uses the concept of deviatoric second-order work for examining the relations between stability and incremental non-coaxiality. Based on a spectral analysis of the constitutive compliance matrix, it proposes three classifications for distinguishing various degrees of incremental non-coaxiality and stability. These classifications determine the conditions for the existence of incremental coaxiality (i.e., colinearity of stress and strain increments), stability, instability, and stable-instable transition (i.e., positive, negative, or zero second-order deviatoric work). The study illustrates these classifications in the cases of generic elastic and elastoplastic constitutive models. The analysis pertains to two-dimensional cases. Additional research is required to extend the analysis from two to three dimensions.
The anisotropic and non-coaxial behaviors of geomaterials are challenging to model using the conventional flow theory of plasticity, which assumes that the strain increments and principal stress have identical direction, that is, are coaxial [1–3]. For instance, the associative flow rule of plasticity, which assumes that strain increments are normal to the yield surface, disagrees with many experimental evidences [4, 5] and micromechanical observations [6, 7], which show non-coaxiality, that is, different principal directions for stress states and strain increments.
As early as Hill , several theories have been proposed to introduce noncoaxiality. For instance, one approach added tangent plasticity to classical coaxial models [1, 9]. Another approach defined strain increment in terms of stress states and material anisotropy, which may lead to an anisotropic hardening law [10–12]. Others have described noncoaxial behaviors using double-shear models [13–15]. Micromechanical studies have related non-coaxiality to the anisotropic fabric resulting from the arrangement of material particles and associated voids [16, 17].
Most studies of non-coaxiality focused on the case of non-proportional loading, and assumed coaxiality under proportional loading . However, experiments showed that anisotropic solids, such as geomaterials, may exhibit considerable degrees of noncoaxiality even under proportional loading [5, 19, 20]. This results in two distinct mechanisms relevant to non-coaxiality. According to theoretical [21, 22] and experimental studies [5, 20], non-coaxiality is largely induced by tangent yield effects under non-proportional loading. For proportional loading, experiments showed that non-coaxiality tends to result from induced or inherent material anisotropy [5, 19]. To distinguish the two different mechanisms, this analysis expands the definition of conventional non-coaxiality to incremental non-coaxiality, which is hereafter referred to as the difference between principal directions of strain and stress increments.
Figure 1 illustrates the relationship among non-coaxiality, incremental non-coaxiality and stress path. It relates these three concepts using stress state , stress increment , and strain increment . The principal direction angles of , and are defined as and , respectively (Figure 2). Like in the solid mechanics literature , Figure 1 defines coaxiality (C, ) and non-coaxiality (NC, ) in terms of and . Proportional (P, ) loading corresponds to proportional to , and non-proportional (NP, ) loading refers to changes in principal direction of , when moves along nonlinear stress paths. Incremental coaxiality (IC, ) and incremental non-coaxiality (INC,) are defined in terms of and . Therefore, NC can be generated in three different ways through INC + P (Figure 2(a)), INC + NP (Figure 2(b)) and IC + NP (Figure 2(c)).
(a) Incremental non-coaxiality with proportional loading ( 𝜃 𝑑 𝜖 ≠ 𝜃 𝑑 𝜎 = 𝜃 𝜎 )
(b) Incremental non-coaxiality with non-proportional loading ( 𝜃 𝑑 𝜖 ≠ 𝜃 𝑑 𝜎 ≠ 𝜃 𝜎 )
(c) Incremental coaxiality with non-proportional loading ( 𝜃 𝑑 𝜖 = 𝜃 𝑑 𝜎 ≠ 𝜃 𝜎 )
The incremental non-coaxiality (INC) and non-proportionality (NP) of loading can be defined, respectively, as
Therefore, non-coaxiality (NC), incremental non-coaxiality (INC), and non-proportionality (NP) of loading are related by where can be determined from the prescribed loading path and is a type-rate material property that needs to be measured from incremental stress-strain response. Clearly, in the particular case of proportional loading (Figure 2(a)).
For stress-controlled problems, the above considerations restrict our attention to incremental non-coaxiality (INC). The present analysis focuses on two aspects: (1) for which materials incremental non-coaxiality (INC) may become negligible or large and (2) in which direction materials may preserve incremental coaxiality (IC). The analysis is presented in two dimensions, but can be extended to three-dimensions with some additional effort.
2. Lame’s Representation of Incremental Non-Coaxiality
Material behaviors are modeled using rate-type linearized solids, that is, where is the incremental strain, is the incremental stress, and is the incremental compliance tensor that represents some underlying anisotropy state. The indices are in two-dimensions. can be written in matrix form as follows:
The physical meaning of incremental non-coaxiality can be illustrated using Lame’s ellipse (Figure 3), which represents a two-dimensional second-order tensor using an inclined ellipse. is uniquely expressed in terms of principal stress values and () and the orientation angle of principal direction. Similarly, the strain increment is defined using the principal values of strain increment and () and the orientation angle of principal direction. The incremental response is coaxial when and non-coaxial when .
(a) Incremental stress ellipse
(b) Incremental strain ellipse
In Lame’s representation, (2.1) is restated as follows: where the nonlinear functions , , and depend on and .
Without loss of generality, (2.3) may be simplified assuming and (see the appendix) where functions and are analogous to , and , and denotes the ratio of incremental principal stresses.
Mapping uniquely a stress increment onto a Lame’s ellipse requires assigning a positive sign to principal values. Otherwise, different stress increments would correspond to the same Lame’s ellipses. For instance, assuming , Table 1 shows five tensors that have the same Lame’s ellipse but different values of and .
3. Classifications of Incremental Non-Coaxiality
Three different classifications are proposed for incremental non-coaxiality, based on the values of coefficients that indicate some kind of fabric anisotropy. These classifications addresses the following two issues: (1) relation of stability and INC, and (2) number of directions that preserve incremental coaxiality, that is, stress and strain increments have identical or opposite directions.
3.1. Classification 1: Energy-Based Classification for INC
The first classification intends to distinguish between materials with strong or weak INC, noting that the degree of INC is intrinsically related to material stability. Since the strain increment has the same principal direction as its deviatoric part, the deviatoric constitutive equation is considered where , and .
The eigenvalues of can be analytically calculated by solving a quadratic equation, which is simpler than solving a cubic equation for the eigenvalues of . This is because is rank deficient and always possess one zero eigenvalue, leaving only two eigenvalues to determine.
The deviatoric second-order work is where is the deviatoric part of , and denote the magnitude of and , respectively. One obtains:
The trivial cases and can be ignored without loss of generality. In general, there are three independent deviation angles between the principal axes of stress and strain increments in three-dimensional case, for which given by (3.3) does not represent a certain angular orientation of incremental non-coaxiality. However, is an energy-based parameter that can be extended to higher-dimensional problems.
Assuming , can be explicitly expressed as which implies that and have the same principal directions.
In the particular case , that is, which corresponds to pressure increments producing no deviatoric strains. With , equation (3.3) becomes where isAccording to (3.6a)-(3.6b), is independent from , and and have identical sign.
Equation (3.1) can be equivalently expressed in matrix form as follows: or
Complicated spectral analyses for nonsymmetric matrices are avoided by noting that any nonsymmetric matrix A satisfies for any vector x. Equation (3.3) can be rewritten where (i.e., ) is the symmetric part of , whose equivalent matrix form is
The eigenvalues of are denoted by , and , and the corresponding orthogonal eigentensors are denoted by , and . Any stress increments can be expressed as where , , and are the components in the eigendirections.
Equation (3.8) becomes which shows that depends on the signs of , and .
When , the eigenvalues of are where Parameters and solely control the values of eigenvalues. Invoking that , Figure 4 shows four distinct domains:(I)stable; ; and ; ,(II)unstable; ; and ; ,(III)partially stable; ; and ; .
Point : Stable-Unstable Transition; ; and ; .
When , the above spectral approach does not fully apply however. This case is similar to partially stable INC (Domain III in Figure 4).
3.2. Classification 2—Spectral Classification for IC
The second classification relates to the material directions that preserve stable incremental coaxiality (SIC) or produce unstable incremental coaxiality (UIC). The analysis looks for the principal directions (if any) of anisotropic materials, which may also represent the optimal orientation of anisotropic solids, see, for example, .
are assumed to have three eigenvalues , , and and three corresponding eigendirections , , and . The eigenvalues of are where Consider that a given stress increment where the coefficients a, b, and c are variables comparable to , and in (3.10). In the case (i.e., )Thus, Equation (3.16) implies that are coaxial with when or .
As shown in (3.13), (3.16) and Figure 5, parameters A and B control the values of eigenvalues and , which determine the stability and degree of IC. The stability of IC (i.e., SIC or UIC) depends on the signs of and . The degree of SIC (or UIC) is indicated with the number n of directions along which the incremental stress and deviator strain have identical (or opposite) direction.
3.3. Classification 3—Couped Spectral Classification of INC
The third classification builds upon the two previous classifications. Hereafter, we focus on the case because creates incremental non-coaxiality (INC) in all directions without possible stable or unstable incremental coaxiality. implies that
Figures 4 and 5 were combined into Figure 6 to produce the third classification, which gives information on stability as well as the number of directions for stable or unstable incremental coaxiality. Table 2 is a tabular representation of classifications in Figure 6. This classification will be illustrated using examples in each classification domains.
|SIC = : number of directions with identical direction of incremental stress and deviator strain, (SIC alone means that incremental stress and deviator strain have identical direction). UIC = : number of directions with opposite direction of incremental stress and deviator strain, (UIC alone means that all incremental stress and deviator strain have opposite direction). |
Stable-Unstable Transition: ; stable: ; unstable: ; partially stable: .
When , which results in (i.e., and ), points O, , and a+ coincide, and Domains I, II and III remain only in Figure 6.
Incremental non-coaxiality is modeled for compliance matrices (2.2) that do not necessarily have a major symmetry, that is, . As shown in Table 3, all components of are dimensionless. With respect to distinct eigenvalue cases, examples are given to demonstrate how the classifications distinguish various degrees of incremental non-coaxiality and stability. Tables 3(a)–3(e) and Figures 7–11 contain 17 examples that apply to anisotropic elasticity as well as elastoplasticity.
|(a) Model parameters for stable INC ()|
|(b) Model parameters for instable INC ()|
|(c) Model parameters for stable-instable INC transition ()|
|(d) Model parameters for partially stable INC ()|
|(e) Model parameters for partially stable INC ()|
(a) 4 SIC
(b) 2 SIC
(d) No SIC
(a) 4 UIC
(b) 2 UIC
(d) No UIC
(a) 4 SIC
(b) 2 SIC
(c) 4 UIC
(d) 2 UIC
(e) 2 SIC and 2 UIC
(f) No SIC/UIC
(a) No SIC/UIC (real 𝜆 2 , 𝜆 3 )
(b) No SIC/UIC (complex 𝜆 2 , 𝜆 3 )
Figure 7 and Table 3(a) show four cases of stable incremental non-coaxiality when . Figures 7(a) and 7(b) show that there are four and two stable IC (SIC) (i.e., ). Figure 7(a) falls within Domain I (Figure 6) while Figure 7(b) is on Boundary (Figure 6). Figure 7(c) shows a particular case of stable incremental coaxiality in all directions, which is available only if , corresponding to Boundary ) in Figure 6. Finally, Figure 7(d) corresponds to complex and and , which falls within Domain VI (Figure 6). In this particular case, INC is stable in all directions. As illustrated in Figure 6, and are two cases that cannot coexist with .
Figure 8 and Table 3(b) show four cases of unstable INC, which corresponds to . In contrast to Figure 7, unstable IC is possible along a few directions. Figures 8(a) and 8(b) show four and two unstable IC (UIC), respectively, which correspond to , and fall within Domain II ( and and along boundary ( and ) in Figure 6. Figure 8(c) shows the particular case of unstable IC in all directions, which is available only if , corresponding to boundary () in Figure 6. Finally, Figure 8(d) shows that there is unstable INC in all directions, which corresponds to Domain VII in Figure 6.
Figure 10 and Table 3(d) show six cases of partially stable INC, which corresponds to . The case can be obtained for any values of and . The classifications correspond to Domains III, IV, V, VIII and boundaries and . Stable IC (SIC) can be obtained in 4 directions in Figure 10(a), 2 in Figure 10(b), and 2 in Figure 10(e). Unstable IC (UIC) can be obtained in 4 directions in Figure 10(c), 2 in Figure 10(d) and 2 in Figure 10(e). Figure 10(f) shows that partially stable INC is obtained for all directions, which corresponds to Domain VIII in Figure 6.
As shown in Figures 7–10, INC is independent of , which results from . On the other hand, Figure 11 presents the variation of INC in terms and , when . The intensity of shaded contours indicates the stability of INC. Table 3(e) lists the corresponding material coefficients. In this case, INC cannot be classified using spectral classification, because it depends not only on the principal direction but also on the ratio of principal stress increments, . Similar results are founds when and are real and complex.
In the above examples, the compliance matrices are clearly related to classical anisotropic materials. For instance, isotropic elasticity always creates stable IC (), while orthogonally anisotropy may create stable INC ( and ) presented in Figure 7(a), or partially stable INC () as shown in Figure 11(a). On the other hand, the coefficients of compliance matrix for elastoplastic materials may change and generate several distinct eigenvalue cases with the development of work hardening (softening).
In addition, as shown in Figures 7 and 8, pure stable (or unstable) INC that requires the condition of , may be reproduced by Tresca-type or Mises-type models, which are generally used to simulate incompressible material behaviors. However, pressure sensitive materials () create partially stable INC only, as shown in Figure 11.
The present study has presented a spectral approach to explore the non-coaxial aspects of incremental material behavior and has classified the non-coaxiality between stress and strain increments. It expands the definition of non-coaxiality (NC), which usually refers to the difference between the principal directions of incremental strain and stress state in the solid mechanics literature. The analysis has investigated the incremental non-coaxiality (INC) induced by incremental stress-strain relations, for example, linearized rate-type solids. Based on the concept of deviatoric second-order work, this study has examined the relations between incremental non-coaxiality and stability. It has proposed three classifications that are based on eigenvalues of the constitutive compliance matrix. These classifications distinguish various degrees of incremental non-coaxiality and stability. They determine the conditions for which stress and strain increments are collinear and result into stability, instability, and stable-instable transition. The classifications have been illustrated using examples of constitutive matrices that are relevant to elastic as well as elastoplastic constitutive modeling. The present classifications are useful for examining the incremental stability of incompressible materials and determining anisotropic state and optimal material orientation.
The incremental stress-strain relationship given by (2.1) can be rewritten as follows: The principal direction of strain increment can be defined as Substitution of (A.1), (A.2) and (A.3) into (A.4) leads to Invoking that stress is related to its principal values by With the introduction of (A6), (A.7) and (A.8) into (A.5), thus (A.5) can be restated as where , , and denote , , and , respectively. Recall that in (A.9). Based on (A.9), we can concluded (2.4).
Analogous to the above analysis, it can be checked that it is also true for (2.5).
The financial support provided by National Natural Science Foundation of China (Grant nos. 10972159 and 50825803) and Natural Science Foundation of Shanghai, China (Grant no. 08ZR1420100).
- J. W. Rudnicki and J. R. Rice, “Conditions for the localization of deformation in pressure-sensitive dilatant materials,” Journal of the Mechanics and Physics of Solids, vol. 23, no. 6, pp. 371–394, 1975.
- E. Papamichos and I. Vardoulakis, “Shear band formation in sand according to non-coaxial plasticity model,” Géotechnique, vol. 45, no. 4, pp. 649–661, 1995.
- J. G. Qian, J. Yang, and M. S. Huang, “Three-dimensional noncoaxial plasticity modeling of shear band formation in geomaterials,” Journal of Engineering Mechanics, vol. 134, no. 4, pp. 322–329, 2008.
- K. H. Roscoe, “The influence of strain in soil mechanics,” Géotechnique, vol. 20, pp. 129–170, 1970.
- M. Gutierrez, K. Ishihara, and I. Towhata, “Flow theory for sand during rotation of principal stress direction,” Soils and Foundations, vol. 31, no. 4, pp. 121–132, 1991.
- A. Drescher and G. de Josselin de Jong, “Photoelastic verification of a mechanical model for the flow of a granular material,” Journal of the Mechanics and Physics of Solids, vol. 20, no. 5, pp. 337–340, 1972.
- M. Oda, J. Konishi, and S. Nemat-Nasser, “Experimental micromechanical evaluation of strength of granular materials: effects of particle rolling,” Mechanics of Materials, vol. 1, no. 4, pp. 269–283, 1982.
- R. Hill, The Mathematical Theory of Plasticity, The Clarendon Press, Oxford, UK, 1950.
- I. Vardoulakis and B. Graf, “Calibration of constitutive models for granular materials using data from biaxial experiments,” Géotechnique, vol. 35, pp. 299–317, 1985.
- Z. Mroz, “An attempt to describe the behavior of metals under cyclic loads using a more general work hardening model,” Acta Mecanica, vol. 7, pp. 199–212, 1969.
- W. Wu, “Rational approach to anisotropy of sand,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 22, no. 11, pp. 921–940, 1998.
- X. S. Li and Y. F. Dafalias, “A constitutive, framework for anisotropic sand including non-proportional loading,” Géotechnique, vol. 54, no. 1, pp. 41–55, 2004.
- A. J. M. Spencer, “A theory of the kinematics of ideal soils under plane strain conditions,” Journal of the Mechanics and Physics of Solids, vol. 12, pp. 337–351, 1964.
- G. de Josselin de Jong, “The double sliding, free rotating model for granular assemblies,” Géotechnique, vol. 21, pp. 155–162, 1971.
- J. A. M. Teunissen, “On double shearing in frictional materials,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 31, no. 1, pp. 23–51, 2007.
- M. Oda, “Inherent and induced anisotropy in plasticity theory of granular soils,” Mechanics of Materials, vol. 16, no. 1-2, pp. 35–45, 1993.
- S. Nemat-Nasser, “A micromechanically-based constitutive model for frictional deformation of granular materials,” Journal of the Mechanics and Physics of Solids, vol. 48, no. 6-7, pp. 1541–1563, 2000.
- H. S. Yu, “Non-coaxial theories of plasticity for granular materials,” in Proceedings of the 12th International Association for Computer Methods and Advances in Geomechanics (IACMAG '08), pp. 361–378, Goa, India, 2008.
- R. K. S. Wong and J. R. F. Arthur, “Induced and inherent anisotropy in sand,” Géotechnique, vol. 35, no. 4, pp. 471–481, 1985.
- K. Miura, S. Miura, and S. Toki, “Deformation behavior of anisotropic dense sand under principal stress axes rotation,” Soils and Foundations, vol. 26, no. 1, pp. 36–52, 1986.
- H. S. Yu and X. Yuan, “On a class of non-coaxial plasticity models for granular soils,” Proceedings of the Royal Society A, vol. 462, no. 2067, pp. 725–748, 2006.
- S. Tsutsumi and K. Hashiguchi, “General non-proportional loading behavior of soils,” International Journal of Plasticity, vol. 21, no. 10, pp. 1941–1969, 2005.
- G. Mandl and R. F. Luque, “Fully developed plastic flow of granular materials,” Géotechnique, vol. 20, pp. 277–307, 1970.
- A. N. Norris, “Optimal orientation of anisotropic solids,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 59, no. 1, pp. 29–53, 2006.
Copyright © 2010 Jiangu Qian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.