Changing Coordinate Systems and Mohr's Circle

Lecture 5: Changing Coordinate Systems and Mohr’s Circle Lecture 2 explained that temperature is a zeroth-order tensor, force is a first-order tensor, and stress is a second-order tensor. The order of a tensor is called its rank and is defined by its law of transformation under a change of coordinates. This lecture explains the transformation laws for first- and second-order tensors, and uses these laws to derive a convenient representation of stress called Mohr’s circle. 1. Transforming Tensors into Different Coordinate Systems Changing coordinate systems In order to transform a tensor into a different coordinate system, one must first understand how to transform the coordinate system itself. Consider the following figure: x3 x3' φ11 x1 P φ23 φ12 x2' x2 x1' Figure 5.1 If xˆi ' and x̂ j represent unit vectors that are the axes of two coordinate systems with the same origin, they are related by the equation xˆ i ' = α ij xˆ j where αij is the cosine of the angle between the primed axis xˆi ' and the unprimed axis x̂ j . For example, α12 is the cosine of the angle between xˆ1 ' and x̂ 2 . αji represents a 9component matrix called the transformation matrix. Unlike the stress tensor, it is not symmetric (αij ≠ αji). ⎡α 11 α 12 α 13 ⎤ α ij = ⎢⎢α 21 α 22 α 23 ⎥⎥ ⎢⎣α 31 α 32 α 33 ⎥⎦ In matix equations, the transformation law is written ⎡ x1 ' ⎤ ⎡α 11 α 12 α 13 ⎤ ⎡ x1 ⎤ ⎥⎢ ⎥ ⎢ x '⎥ = ⎢α ⎢ 2 ⎥ ⎢ 21 α 22 α 23 ⎥ ⎢ x 2 ⎥ ⎢⎣ x3 '⎥⎦ ⎢⎣α 31 α 32 α 33 ⎥⎦ ⎢⎣ x3 ⎥⎦ The inverse transformation law is written xˆ i = α ji xˆ j ' Example of coordinate transformations Consider the following transformation of coordinates: x2 x2' x1' 30o x3, x3' Figure 5.2 The transformation matrix is ⎡ 30° 60° 0⎤ α ij = ⎢⎢120° 30° 0⎥⎥ ⎢⎣ 0 0 0⎥⎦ x1 The explicit transformation equations are xˆ1 ' = cos 30° xˆ1 + cos 60° xˆ 2 xˆ 2 ' = cos120° xˆ1 + cos 30° xˆ 2 Since xˆi ' and x̂ j are both unit length, these equations are easy to verify from the picture. First-order tensors First-order tensors or vectors have two components in 2D coordinates and three components in 3D coordinates. They transform according to the same laws as coordinate axes because coordinate axes are themselves vectors. If uj is a vector in the x̂ j coordinate system and ui’ is a vector in the xˆi ' coordinate system, then the following equations describe their transformation: u i ' = α ij u j u i = α ji u j ' Note that αij is positive if the angle is measured counterclockwise from xˆi ' to x̂ j . It is negative if the angle is measured clockwise. Second-order tensors The transformation law for second-order tensors like stress and strain is more complicated than the transformation law for first-order tensors. It may be derived as follows: a. Begin with the vector transformation of traction Ti to Tk’: Tk ' = α ki Ti b. Rewrite Tk’ and Ti using Cauchy’s formulas: Tk ' = σ ' kl nl ' and Ti = σ ij n j Substitute Cauchy’s formulas into the original transformation equation: σ ' kl nl ' = α kiσ ij n j c. Transform the normal vector nj to nl’ and substitute into the previous equation: n j = α lj nl ' σ ' kl nl ' = α kiσ ij α lj nl ' d. Cancel the nl’ term on each side and group the αs: σ ' kl = α kiα jl σ ij Note that changing the position of the last α term changes the order of its subscripts. In vector notation, the equation is σ ' = α σα T where the double underbars denote second-rank tensors and the superscript T denotes the transpose of matrix α. 2. Mohr’s Circle Motivation Lecture II explained that an object resting on a slope will slide down when the shear traction on the slope is greater than or equal to the product of the normal traction and the coefficient of friction. τ = f sσ n On a shallow slope, σn is large and the object will not slide. On a steep slow, τ is large and the object will slide. For any plane with normal n̂ , we can calculate if the plane will fail if the stress tensor σij at the interface between the object and the slope is known. Calculate σn and , τ as follows: Vector and Tensor Notation r T = σ nˆ r σ n = T ⋅ nˆ r r τ = T − T ⋅ nˆ Summation Notation Ti = σ ij n j σ n = Ti ni τ = Ti − Ti ni This method is straightforward but cumbersome. A different approach involves rotating the coordinate system such that x1’ is along n̂ . In this case σn and τ are much easier to derive: σ n = σ '11 τ = σ '12 Deriving Mohr’s Circle Mohr’s circle may be derived in two or three dimensions. This lecture explains the derivation in two dimensions because it is more straightforward and the results are easier to graph and understand. The derivation assumes that x1, x2, and x3 are principle directions. Consider the following figure in which the xi coordinate system is rotated clockwise about the x3 axis to xi’: x1 θ x1' φ11 x2' x2 φ22 φ21 Figure 5.3 The rotation matrix is: ⎡α 11 α 12 α 13 ⎤ ⎡cos θ α ij = ⎢⎢α 21 α 22 α 23 ⎥⎥ = ⎢⎢ sin θ ⎢⎣α 31 α 32 α 33 ⎥⎦ ⎢⎣ 0 − sin θ cos θ 0 0⎤ 0⎥⎥ 1⎥⎦ The stress tensor σ in the xi coordinate system is transformed to σ ' in the xi’ coordinate system by the following equation: σ ' = α σα T ⎡cos θ σ ' = ⎢⎢ sin θ ⎢⎣ 0 − sin θ 0⎤ ⎡σ 11 0 0⎥⎥ ⎢⎢ 0 σ 22 1⎥⎦ ⎢⎣ 0 0 cos θ 0 ⎡σ 11 cos θ σ ' = ⎢⎢σ 11 sin θ ⎢⎣ 0 − σ 22 sin θ σ 22 cos θ 0 0 ⎤ ⎡ cos θ 0 ⎥⎥ ⎢⎢− sin θ σ 33 ⎥⎦ ⎢⎣ 0 0 ⎤ ⎡ cos θ 0 ⎥⎥ ⎢⎢− sin θ σ 33 ⎥⎦ ⎢⎣ 0 sin θ cos θ 0 sin θ cos θ 0 0⎤ 0⎥⎥ 1⎥⎦ 0⎤ 0⎥⎥ 1⎥⎦ ⎡σ 11 cos 2 θ + σ 22 sin 2 θ (σ 11 − σ 22 ) sin θ cos θ 0 ⎤ ⎢ ⎥ 2 2 σ ' = ⎢ (σ 11 − σ 22 ) sin θ cos θ σ 11 sin θ + σ 22 cos θ 0 ⎥ ⎢ 0 0 σ 33 ⎥⎦ ⎣ Use the double-angle identities for sine and cosine to simplify the expressions for the normal stress σ11’ and the shear stress σ12’ become in the new coordinate system: sin 2θ = 2 sin θ cos θ cos 2θ = 2 cos 2 θ − 1 = 1 − 2 sin 2 θ σ 11 ' = σ n = σ 11 + σ 22 + σ 11 − σ 22 2 2 σ − σ 22 σ 12 ' = τ = 11 sin 2θ 2 cos 2θ The expression for the normal stress and shear stress can be shown graphically in shear space: τ Allowable tractions 2θ σ2 σ1 τ Figure 5.4 This figure is called Mohr’s circle. σn Interpreting Mohr’s circle Mohr’s circle plots in stress space. Admonton’s law may also be plotted in stress space as a line with slope fs. When this line and Mohr’s circle intersect, the criterion for failure across a plane is met. Consider a common experiment in rock mechanics in which scientists apply a uniaxial stress σ2 to a cylindrical sample confined by a uniform stress σ1. x2 x1 x3 σ2 σ1 σ1 Figure 5.5 Admonton’s law and the Mohr’s circle that represents the state of stress are plotted below. The Mohr’s circle plots on the negative σn axis because these notes follow the convention that compressive stresses are negative. σ1 plots farther to the right because of the convention that σ1 is greater than σ2, and σ2 is greater than σ3. The line that represents Admonton’s law crosses the τ axis at σo, the shear strength of the rock. τ |τ| =|fsσn|+|σ0| σ0 σ2 σ1 σn Figure 5.6 As the uniaxial stress σ2 is increased, the Mohr’s circle becomes larger. When the circle intersects the line of Admonton’s law, the rock breaks. The angle 2θ at which the circle intersects the line is twice the angle between σ1 and the normal vector to the failure plane. τ σ2 2θ φ = tan-1fs σ1 Figure 5.7 σn The above picture shows that the coefficient of friction determines the slope of Admonton’s law, and the slope of Admonton’s law determines 2θ. Consequently, predicting the failure plane of a rock only requires knowing the direction of the most compressive stress σ2, the direction of the least compressive stress σ1, and the coefficient of friction fs. The table below lists the angle between σ1 and the normal vector to the failure plane for different values of fs. fs 0.0 0.6 1.0 lim f s →∞ θ 45° 30° 23° 0° Since most rocks have a coefficient of friction of about 0.6, the normal vector to the failure plane is typically 30°from the direction of the least compressive stress. Another way of saying this is that the failure plane is 30°from the direction of the most compressive stress. Faulting Faults are large-scale failure planes. Since failure planes are normal to the plane containing the least compressive and the most compressive stresses and are typically at 30°from the direction of the most compressive stress, different styles of faulting can be used to infer the directions of σ1, σ2, and σ3. The following diagram shows the deviatoric stresses associated with thrust faults, normal faults, and strike-slip faults. The pictures assume that x2 is vertical and that x1 and x3 lie in the plane of the surface of the Earth. Figure 5.8